# Burial-induced oxygen-isotope re-equilibration of fossil foraminifera explains ocean paleotemperature paradoxes

Samuel Bernard, D. Daval, P. Ackerer, S. Pont, A. Meibom

► To cite this version:

Samuel Bernard, D. Daval, P. Ackerer, S. Pont, A. Meibom. Burial-induced oxygen-isotope re-equilibration of fossil foraminifera explains ocean paleotemperature paradoxes. Nature Communications, 2017, 8, pp.1134. 10.1038/s41467-017-01225-9. hal-01633900

HAL Id: hal-01633900

https://hal.sorbonne-universite.fr/hal-01633900

Submitted on 13 Nov 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

DOI: 10.1038/s41467-017-01225-9

OPEN

# Burial-induced oxygen-isotope re-equilibration of fossil foraminifera explains ocean paleotemperature paradoxes

S. Bernard<fnref n="1" />, D. Daval², P. Ackerer², S. Pont<fnref n="1" /> & A. Meibom³,⁴

Oxygen-isotope compositions of fossilised planktonic and benthic foraminifera tests are used as proxies for surface- and deep-ocean paleotemperatures, providing a continuous benthic record for the past 115 Ma. However, visually imperceptible processes can alter these proxies during sediment burial. Here, we investigate the diffusion-controlled re-equilibration process with experiments exposing foraminifera tests to elevated pressures and temperatures in isotopically heavy artificial seawater (H₂¹⁸O), followed by scanning electron microscopy and quantitative NanoSIMS imaging: oxygen-isotope compositions changed heterogeneously at submicrometer length scales without any observable modifications of the test ultrastructures. In parallel, numerical modelling of diffusion during burial shows that oxygen-isotope re-equilibration of fossil foraminifera tests can cause significant overestimations of ocean paleotemperatures on a time scale of 10⁷ years under natural conditions. Our results suggest that the late Cretaceous and Paleogene deep-ocean and high-latitude surface-ocean temperatures were significantly lower than is generally accepted, thereby explaining the paradox of the low equator-to-pole surface-ocean thermal gradient inferred for these periods.

<footnote n="1">¹IMPMC, Sorbonne Universités, CNRS UMR 7590, MNHN, UPMC, IRD UMR 206, 61 Rue Buffon, 75005 Paris, France. ²LHyGeS, CNRS UMR 7517, Université de Strasbourg/EOST, 1 Rue Blessig, 67084 Strasbourg, France. ³Laboratory for Biological Geochemistry, School of Architecture, Civil and Environmental Engineering, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland. ⁴Center for Advanced Surface Analysis, Institute of Earth Sciences, University of Lausanne, 1015 Lausanne, Switzerland. Correspondence and requests for materials should be addressed to S.B. (email: sbernard@mnhn.fr)</footnote>
Precise knowledge of the past ocean temperature is essential for understanding hydrosphere evolution and for placing anthropogenic global climate change in a geologic perspective. Since the early 1950s, the O isotope compositions of fossil foraminifera tests have been an important tool in paleoclimate research¹,². Field and laboratory studies have demonstrated that the ¹⁸O/¹⁶O ratio of living foraminifera calcite tests is a function of both the temperature and the O isotope composition of the seawater, the latter of which varies with global ice volume, pH and salinity³⁻⁵. The most widely accepted interpretation of the continuous benthic foraminifera O isotope record for the past 115 Myr is that the Cretaceous deep ocean was very warm and continuously cooled by ~15 °C during the late Cretaceous and the Paleogene⁶⁻⁹.

The corresponding planktonic O isotope record was initially interpreted to indicate relatively warm high-latitude sea-surface temperatures and relatively cold tropical sea-surface temperatures during the late Cretaceous and the Paleogene¹⁰, giving rise to the cool tropics paradox¹¹. This paradox was related to secondary calcite precipitation that distorts paleotemperature reconstructions¹²⁻¹⁶. In fact, tests showing no visible alterations, also known as glassy foraminifera, yielded significantly higher temperatures for the Cretaceous and Paleogene low-latitude surface oceans¹⁷,¹⁸. Still, the currently accepted interpretation of the planktonic and benthic foraminifera records is that both the vertical thermal gradient in the tropical ocean and the equator-to-pole surface-ocean temperature gradient were much less steep during these periods than those in the present ocean⁶⁻¹⁰,¹⁷,¹⁸. However, such thermal gradients cannot be reconciled with the most recent climate and ocean circulation models¹⁹,²⁰.

The O isotope composition of fossil foraminifera tests can be modified during burial without any visible structural changes. Despite the early warning of Urey et al.²¹ and recent evidence that carbonates can undergo isotope re-equilibration at low temperatures²², the potential bias of paleotemperature estimates resulting from such visually imperceptible processes has never been quantified. Isotope re-equilibration of foraminifera tests with sediment pore water likely occurs over millions of years

during sediment burial at relatively low temperatures (~20–30° C) and pressures (200–500 bars). However, such a re-equilibration cannot be visualised in natural specimens, even when using high-resolution analytical techniques, such as the NanoSIMS ion microprobe. In fact, although the NanoSIMS can quantify isotope compositions at a lateral resolution of ~100 nm²³, it does not have the required analytical precision to reveal the resulting permil-level variations.

Here, we conducted experiments with water strongly enriched in ¹⁸O at high temperatures to visualise the burial-induced isotope re-equilibration of foraminifera tests using quantitative NanoSIMS imaging. The results demonstrate that the O isotope compositions of foraminifera tests can change without any observable modifications of their ultrastructures. Numerical simulations suggest that the late Cretaceous and Paleogene deep-ocean and high-latitude surface-ocean temperatures were significantly lower than is generally accepted.

## Results

**Isotope re-equilibration experiments.** Sealed gold capsules, each containing 160 µg of cleaned foraminifera tests (i.e. ~12 specimens of *Globigerina bulloides* (bulk δ¹⁸O = 1.35‰ Vienna Pee Dee Belemnite (VPDB) (±0.05‰; 2σ)) from modern sediments of the Lion Gulf, France) and 100 µL of an artificial H₂¹⁸O seawater solution with 0.55 mol L⁻¹ of NaCl and 0.003 mol L⁻¹ of NaHCO₃, were submitted to a temperature of 300 °C under a confining pressure of 200 bars for 3 months. Control experiments were performed under identical conditions with chemically similar but isotopically normal artificial seawater. The foraminifera tests were ultrasonically cleaned in pure ethanol before these experiments to remove the clays attached to their surfaces and were cleaned after the experiments to remove as much of the ¹⁸O-pure artificial seawater adhering to their surfaces as possible. The fine-scale morphology of these tests was documented before and after the experiments with scanning electron microscopy (SEM). Finally, the tests were embedded in epoxy resin and polished for quantitative NanoSIMS isotope imaging.

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Foraminifera tests prior to isotope re-equilibration experiments. **a–e** SEM images showing the overall morphologies and ultrastructures of the foraminifera tests (*G. bulloides*) after their ultrasonic cleaning in pure ethanol. **f** Backscattered SEM image of a polished section of a test embedded in epoxy. **g, h** Energy dispersive X-ray spectroscopy (EDXS) maps showing the spatial distributions of carbonates (Ca appears in yellow) and clays (Si appears in pink). **i** NanoSIMS map displaying the isotopically homogeneous nature of the carbonates and clays (¹⁸O/¹⁶O = 0.002). The red rectangles in **f, g** indicate the locations of **i**. The red rectangle in **i** indicates the location of **h**. Scale bars are 100 µm (**a, b, f, g**), 50 µm (**c**), 20 µm (**d, i**) and 10 µm (**e, h**)</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Foraminifera tests after isotope re-equilibration experiments. **a–e** SEM images showing the overall morphologies and ultrastructures of the foraminifera tests (*G. bulloides*) after their ultrasonic cleaning in pure ethanol. **f** Backscattered SEM image of a polished section of a test embedded in epoxy. **g, h** EDXS maps showing the spatial distributions of carbonates (Ca appears in yellow) and clays (Si appears in pink). **i, j** NanoSIMS maps showing the ¹⁸O/¹⁶O ratio distributions in carbonates and clays. Carbonates exhibit a ¹⁸O/¹⁶O ratio ranging from 0.1 to 0.8, while clays exhibit a ¹⁸O/¹⁶O ratio ranging from 0.9 to greater than 1.0. Note that the isotope exchange occurred very heterogeneously, leading to areas that are more or less enriched in ¹⁸O. The red rectangles in **f, g** indicate the locations of **i, j**. The red rectangle in **i** indicates the location of **h**. Scale bars are 100 µm (**a, b, f, g**), 50 µm (**c**), 20 µm (**d, i**) and 10 µm (**e, h, j**).</figcaption></figure>

The morphology of the foraminifera tests that were submitted to these experiments did not discernibly differ from the starting materials, even at submicrometer length scales (Figs. 1 and 2). In contrast, the NanoSIMS imaging revealed that the average ¹⁸O/¹⁶O ratio of the tests immersed within the ¹⁸O-pure artificial seawater, initially 0.002 (Fig. 1), reached up to 0.15 during the experiments (Fig. 2), which is equivalent to the replacement of ~15 vol% of the initial biogenic calcite by pure CaC¹⁸O¹⁸O¹⁸O. Note that areas extremely enriched in ¹⁸O were observed (Fig. 2), showing that the isotope exchange did not occur homogeneously within the biogenic calcite matrix. These results unambiguously demonstrate a process in which the foraminifera tests undergo modifications to their bulk O isotope composition without any appreciable modification of their fine scale morphology, raising serious questions about the reliability of fossil foraminifera tests (including the so-called glassy specimens) as proxies of seawater temperatures in the geological past.

**Isotope re-equilibration processes.** Mechanistically, the observed O isotope exchange can occur either through coupled dissolution and reprecipitation at mineral-fluid interfaces (replacement that leads to the formation of secondary calcites while retaining the original shape of the primary calcite structures)²⁴,²⁵, or through (much slower) solid-state grain boundary and volume diffusion (movement and transport of oxygen atoms within the biogenic calcite matrix without dissolution)²⁶⁻²⁸. Coupled dissolution and reprecipitation at mineral-fluid interfaces requires the aqueous phase and the dissolving solids to remain in contact to propagate a replacement front²⁴,²⁵. This replacement process is a volume deficit reaction that generates porous secondary products²⁴,²⁵. Because there is no volume change in the replacement of ¹⁸O-poor calcite by ¹⁸O-rich calcite, this process is not considered here. In the following, we thus investigate the O isotope re-equilibration of foraminifera tests during sediment burial under the conservative assumption that this re-equilibration occurs exclusively through solid-state grain boundary and volume
diffusion, thereby exploring the slowest process that can alter the original foraminifera isotope signal.

Once incorporated into sediments, fossil foraminifera tests experience increasing temperatures with increasing sediment burial depths. Because O isotope fractionation between calcite

and water decreases with increasing temperature (see ‘Methods’), the burial-induced increase of the sediment temperature establishes isotopic disequilibrium between the tests and the surrounding pore water. This disequilibrium constitutes the driving force for diffusion, which works continuously towards re-establishing isotopic equilibrium between the pore water and tests.

Electron and atomic force microscopy of the foraminifera tests show that they consist of 50–250 nm calcite domains<sup>29,30</sup> embedded in an organic matrix that is highly prone to degradation at shallow burial depths<sup>31</sup>. This, in combination with the formation of cracks of various length scales, allows most of the submicrometric calcite domains of fossil foraminifera tests to be directly in contact with pore water. Subsequent solid-state grain boundary and volume diffusion occurs along the boundaries and into individual nanoscale calcite domains. In natural settings, at temperatures characterising a typical sedimentary stack, solid-state grain boundary and volume diffusion in calcite is a slow process that occurs over very small distances (<tens of nm). However, given the small sizes of the calcite domains<sup>29,30</sup>, diffusion can significantly impact the bulk O isotope composition of the foraminifera tests over geologic timescales. Because full re-equilibration cannot be achieved in natural settings on relevant timescales, fossil foraminifera tests remain isotopically heterogeneous at the nanoscale, as is the case in the present experiments (Fig. 2).

**Isotope re-equilibration impact.** Using numerical modelling, we quantified changes in fossil foraminifera test O isotope composition in response to burial-induced O isotope re-equilibration through solid-state diffusion assuming a constant chemical composition of the pore water. Here, the sediment pore water was assumed to have the same O isotope composition as the seawater in which foraminifera formed<sup>32</sup>, i.e. a δ<sup>18</sup>O of −1‰ Vienna Standard Mean Ocean Water (VSMOW), consistent with the absence of polar ice caps<sup>2</sup>. Realistic values were chosen for the burial rate, the geothermal gradient and the activation energy. The burial rate (i.e. the sedimentation rate) was adopted from a compilation of the representative DSDP, ODP and IODP sites (http://deepseadrilling.org/; Fig. 3a). Geothermal gradients from 40 to 60 °C km<sup>-1</sup>, a range consistent with reported observations<sup>31</sup>, were tested. Note that in these scenarios, foraminifera are never exposed to temperatures exceeding 25–30 °C (Fig. 3a). The activation energy for oxygen diffusion was varied between 85 and 95 kJ mol<sup>-1</sup>, consistent with values proposed in the literature<sup>26,27</sup> (see ‘Methods’).

Figure 3 shows the results of the simulations that assume benthic foraminifera tests that formed in equilibrium with a

**Fig. 3** Assessment of the impact of burial-induced isotope re-equilibration on the benthic foraminifera O isotope record. **a** A five-point moving average of the burial depths at the DSDP-ODP-IODP sites from which the most comprehensive O isotope compilation is principally derived<sup>9</sup>. The red line indicates the depth profile used for the numerical simulations. The corresponding sediment temperatures are indicated in purple. **b, c** The benthic foraminifera O isotope compilation from Friedrich et al.<sup>9</sup> (grey) is shown with the results from numerical simulations run with a diffusion activation energy of 90 kJ mol<sup>-1</sup> and a range of geothermal gradients in **b**, and with a geothermal gradient of 50 °C km<sup>-1</sup> and a range of diffusion activation energies in **c**. These curves represent the present-day O isotope compositions of benthic foraminifera tests that underwent diffusion-controlled O isotope re-equilibration during burial, assuming an ice-free ocean (δ<sup>18</sup>O = −1‰ VSMOW) and a temperature of 3.5 °C at the water-sediment interface

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Assessment of the impact of burial-induced O isotope re-equilibration on the planktonic foraminifera O isotope records. **a** Planktonic foraminifera O isotope data from Pearson et al.¹⁷ (Tanzania, Paleolatitude ~20° S) and from Wilson and Norris ¹⁸ (ODP Site 1052, Paleolatitude ~20° N) are shown together with the results from numerical simulations run with a geothermal gradient of 50 °C km⁻¹ and a diffusion activation energy of 90 kJ mol⁻¹, assuming an ice-free ocean (δ¹⁸O = -1‰ VSMOW) and a temperature of 3.5 °C at the water-sediment interface. The curves represent the present-day O isotope compositions of fossil planktonic foraminifera tests formed in equilibrium with seawaters at -2, 2, 7, 13, 20, 25, 27 and 27.5 °C. These temperatures roughly correspond to those of surface waters in the modern ocean at 70, 60, 50, 40, 30, 20, 10 and 0° latitude¹⁰. **b** Latitudinal profiles of the O isotope compositions of the planktonic foraminifera for different simulated burial durations: a latitude gradient similar to the modern one (yellow curve) will flatten substantially after 45 Myr (red curve) and 100 Myr (brown curve) in response to the burial-induced O isotope re-equilibration. Also shown are the Eocene (45 Ma) tropical planktonic foraminifera O isotope data from Pearson et al.¹⁷ (red) and the high-latitude planktonic foraminifera O isotope data from Zachos et al.¹⁰ (grey). The temperature conversion is from Anderson and Arthur ³ and assumes an ice-free ocean with a seawater δ¹⁸O = -1 ‰ (VSMOW)</figcaption></figure>

seawater at 3.5 °C. Simulations show that burial-induced O isotope re-equilibration through diffusion can cause a significant shift in the O isotope compositions of fossil benthic foraminifera tests on a time scale on the order of 10⁷ years (>0.1‰, i.e. larger than twice the typical standard deviation of δ¹⁸O measurements) (Fig. 3b, c). Worse, foraminifera tests that formed 100 Myr ago in equilibrium with a seawater at 3.5 °C would today have experienced isotope re-equilibration equivalent to a decrease in their bulk δ¹⁸O of about 3‰ (Fig. 3b, c). The bulk O isotope compositions of these re-equilibrated foraminifera tests could be interpreted to indicate a seawater about 15 °C warmer. Simulations thus show that burial-induced isotope re-equilibration of foraminifera tests leads to substantial paleotemperature overestimations. Our results suggest that instead of indicating a global cooling of the deep ocean during the late Cretaceous and Paleogene, the benthic foraminifera O isotope record principally reflects burial-induced isotope re-equilibration through solid-state diffusion. In other words, the late Cretaceous and Paleogene deep oceans were likely much colder than is currently thought.

Because the sediment stack heats up, the diffusive isotope re-equilibration impacts the O isotope composition of the fossil foraminifera tests that formed in cold waters (i.e. benthic species and high-latitude planktonic species) more than those that formed in warm waters (i.e. tropical planktonic species), such as

is shown by the slopes of the curves in Fig. 4a. The paleotemperatures interpreted from the cold-water foraminifera species are therefore more prone to be overestimated than the paleotemperatures from tropical planktonic species. Numerical simulations show that burial-induced O isotope re-equilibration significantly flattens the inferred paleo-latitudinal temperature gradient on a time scale of 45 Myr (Fig. 4b). Corrected for isotope re-equilibration, a steeper temperature gradient between low and high-latitude surface-ocean waters is re-established for the late Cretaceous and the Paleogene, i.e. a gradient similar to the modern one. This resolves the paradox of the late Cretaceous and the Paleogene low equator-to-pole surface-ocean thermal gradient and is consistent with climate and ocean circulation models¹⁹,²⁰.

In conclusion, accounting for the diffusion-controlled burial-induced O isotope re-equilibration of fossilised benthic foraminifera removes the requirement for a strong, continuous and global cooling of the deep-ocean (on the order of 15 °C) during the late Cretaceous and the Paleogene. Furthermore, the present study suggests that the vertical and latitudinal temperature gradients of the late Cretaceous and Paleogene oceans were likely not very different from the current ones. Importantly, because O isotope re-equilibration through solid-state diffusion is a slow process, it likely had little impact on recent (<10 Ma) high frequency signals, such as the glacial to interglacial fluctuations³³
(driven by oscillations of Earth's orbit and mainly related to fluctuations of the seawater O isotope composition). However, these processes have potentially attenuated the relative amplitudes of older, transient signals, such as the Eocene Oligocene transition or the Palaeocene Eocene thermal maximum⁷⁻⁹.

## Methods

**Isotope re-equilibration experiments.** A synthetic pure ¹⁸O seawater solution was prepared by adding NaCl and NaHCO₃ to pure H₂¹⁸O until reaching 0.55 and 0.003 mol L⁻¹, respectively. Sub-modern planktonic foraminifera (G. bulloides) from modern sediments of the Gulf of Lion, France, with bulk δ¹⁸O ~ 1.35 ± 0.05‰ (VPDB) were rinsed three times in ethanol using an ultrasonic bath. Gold capsules, each containing 160 µg of cleaned foraminifera tests (~12 specimens) and 100 µL of this synthetic seawater solution, were sealed using a Lampert PUK 4 welding machine and placed in Parr autoclaves for 82 days at T = 300 °C and P = 200 bars. Note that with a ¹⁸O/(¹⁶O+¹⁸O) ratio of the foraminifera test close to 0 (because of the rarity of ¹⁸O) when the ratio of the solution is close to 1, the mass balance is such that the water δ¹⁸O can be considered constant during the experiments, regardless of the magnitude of the isotope exchange during the experiments.

Speciation calculations were performed by running the geochemical code CHESS³⁴ using the thermodynamic database of the EQ3/6 code³⁵. Activity coefficients for aqueous species were calculated using the Davies equation³⁶; the electrical balance was achieved using the H⁺ concentration. At ambient temperatures (293 K), the solution is slightly undersaturated with respect to calcite (within uncertainty bounds). At this temperature, assuming that the dissolution of calcite is thermodynamically controlled (i.e. neglecting kinetic barriers), equilibrium would be reached after the dissolution of <2% of the tests. At the experimental temperature (573 K), because of the retrograde solubility of calcite, saturation with respect to calcite is reached for a dissolution progress corresponding to only 0.1% of the tests. Therefore, we can assume that the NanoSIMS images reflect an isotope re-equilibration very close to chemical equilibrium, with little or no contributions from secondary crystallisation from the bulk solution.

Because the carbonate 'building blocks' of planktonic and benthic foraminifera are the same at microscales²⁹,³⁰, we confidently use the results of experiments conducted on planktonic foraminifera as representing the results for both planktonic and benthic foraminifera.

**Characterisation techniques.** SEM observations were performed on foraminifera tests deposited on aluminium stubs and with 5-nm thick coatings of gold using a SEM-FEG Ultra 55 Zeiss (IMPMC, Paris, France) microscope operating at a 2-kV accelerating voltage and a working distance of 2 mm for imaging with secondary electrons and at a 15-kV accelerating voltage and a working distance of 7.5 mm for imaging with backscattered electrons and EDXS mapping.

Isotope maps were produced on triplicate samples with the Cameca NanoSIMS 50 (IMPMC, Paris, France). A thorough technical explanation of the NanoSIMS instrument has been provided in a recent review²³. Briefly, the NanoSIMS ion microprobe is a secondary ion mass spectrometry instrument characterised by an extremely high spatial resolution, a high sensitivity, a high mass-resolving power and a multicollection capability³⁴. Using a focused primary beam of ¹³³Cs⁺ ions, secondary ions were sputtered from the sample surface, typically to a depth of ~100 nm. ¹⁸O⁻ and ¹⁶O⁻ ions from the sample were simultaneously detected (multicollection mode) by electron multipliers at a mass-resolving power of ~9000 (M/ΔM). At this mass-resolving power, the measured secondary ions were resolved from the potential interferences. Images were obtained from a presputtered surface area (the same surface area previously imaged using SEM) by rastering the primary beam across the sample surface. The primary beam was focused to a spot size of ~150 nm, and the pixel size was adjusted so that it was smaller than the size of the primary beam. An electron gun supplied the electrons to the sputtered surface during the analysis to compensate for positive charge deposition from the primary beam and to minimise the surface charging effects. Imaging data were processed using custom-made software (LIMAGE, L. Nittler, Carnegie Institution of Washington).

**Modelling the experiments.** In the present experiments, the driving force for the isotope re-equilibration through solid-state diffusion is the large ¹⁸O excess of the fluid. The most generic expression describing the temperature dependence of the diffusion coefficient can be written as follows:

$$D = D_0 \exp \left( - \frac{E_a}{RT} \right), \quad (1)$$

where D₀ is the intrinsic diffusion constant, Eₐ is the activation energy of oxygen diffusion in calcite, R is the ideal gas constant and T is the absolute temperature. Although the bulk diffusion of oxygen in calcite is the sum of the contributions of both the grain boundary and volume diffusion, it is commonly assumed that the volume diffusion is slow and can be ignored³⁷,³⁸.

Applied to the present case, Eq. (1) becomes:

$$D_{\text{foram}} = D_{0,\text{foram}} \exp \left( - \frac{E_{\text{aforam}}}{RT_{\text{xp}}} \right), \quad (2)$$

where Tₓp is the absolute temperature of the experiment and Eₐₗoram is the activation energy of oxygen solid-state diffusion in foraminifera calcite. No consensus exists in the literature on the value of this last parameter (see below).

In the present experiments, the initial and boundary conditions can be expressed as follows:

$$C(x, t = 0) = 0, \quad \forall x > 0 \quad (3)$$

$$C(x = 0, t) = C_0 = 1 \quad (4)$$

with C standing for the ¹⁸O concentration (i.e. the molar ¹⁸O/(¹⁶O+¹⁸O) ratio) in foraminifera calcite, t is time and x is the length (positive distance from the fluid/ solid interface into the solid).

The diffusion equation for C reads as follows (Fick's second law):

$$\frac{\partial C}{\partial t} = D_{0,\text{foram}} \exp \left( - \frac{E_{\text{aforam}}}{RT_{\text{xp}}} \right) \frac{\partial^2 C}{\partial x^2}, \quad (5)$$

According to Crank³⁹, Because the calcite fraction quantitatively affected by the diffusion process is generally small compared to the total calcite size/volume, the analytical solution for the diffusion into an infinite one-dimensional medium can be applied to calculate the ¹⁸O concentration along the profile after the experimental duration tₓp:

$$C(x) = C_0 \left( 1 - \operatorname{erf} \left( x / \left( 2 \sqrt{D_{0,\text{foram}} \exp \left( - \frac{E_{\text{aforam}}}{RT_{\text{xp}}} \right) t_{\text{xp}}} \right) \right) \right), \quad (6)$$

Because the error function is defined as follows:

$$\operatorname{erf}(y) = \frac{2}{\sqrt{\pi}} \int_{0}^{y} \exp(-\tau^2) d\tau, \quad (7a)$$

with

$$y = \frac{x}{2\sqrt{D_{0,\text{foram}} \exp\left(-\frac{E_{\text{aforam}}}{RT_{\text{xp}}}\right) t_{\text{xp}}}}, \quad (7b)$$

its derivative can be written as follows:

$$\frac{\partial}{\partial y}(\operatorname{erf}(y)) = \frac{2}{\sqrt{\pi}} \exp(-y^2). \quad (8)$$

The flux of matter F entering the calcite domain at the interface x = 0 can be expressed as follows:

$$F(0, t) = -D_{\text{foram}} \frac{\partial C}{\partial x}, \quad (9)$$

with (according to Eqs. (6) and (8)):

$$\frac{\partial C}{\partial x} = \frac{\partial C \partial y}{\partial y \partial x} = \frac{-C_0 \exp(-y^2)}{\sqrt{\pi} D_{\text{foram}} t}, \quad (10)$$

leading to the following expression for the flux F at the interface x = 0:

$$F(0, t) = \frac{C_0 \sqrt{D_{\text{foram}}}}{\sqrt{\pi t}}, \quad (11)$$

The cumulative flux (Fₐ), which represents the amount of ¹⁸O that has entered the calcite domain during the experimental duration tₓp, can be expressed as follows:

$$F_{\text{c}}(0, t_{\text{xp}}) = \int_{0}^{t_{\text{xp}}} F(0, \tau) d\tau = C_{0} \frac{\sqrt{D_{\text{foram}}}}{\sqrt{\pi}} \int_{0}^{t_{\text{xp}}} \frac{1}{\sqrt{\tau}} d\tau = 2C_{0} \frac{\sqrt{D_{\text{foram}} t_{\text{xp}}}}{\sqrt{\pi}}, \quad (12)$$

From a physical standpoint, this result illustrates that the integral of any diffusion profile that satisfies the abovementioned boundary conditions is equivalent to that of a step function of the height C₀ with a length d of:

$$d = \frac{2}{\sqrt{\pi}} \sqrt{D_{\text{foram}} t_{\text{xp}}}, \quad (13)$$

Because the initial ¹⁸O/(¹⁶O+¹⁸O) value of the foraminifera is close to 0, the ¹⁸O/(¹⁶O+¹⁸O) value of the foraminifera calcite crystals at the end of the experiment (Q) is equivalent to the ratio of the re-equilibrated volume to the initial volume. Assuming that the foraminifera calcite crystals have a simple spherical geometry with a radius r₀ and that d << r₀ (as suggested by NanoSIMS imaging),
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 5** Determination of the $Ea_{foram}$ and $D_{0,foram}$ values for the oxygen solid-state diffusion in foraminifera calcites. **a** Foraminifera calcite grain sizes (2\*r₀ values) obtained using Eq. (15) for the different x values (Eq. (16)). **b** $D_{0,foram}$ and $Ea_{foram}$ values obtained using Eq. (16) for the different x values</figcaption></figure>

Q can be written as follows:

$$Q = \frac{\frac{4}{3}\pi r_0^3 - \frac{4}{3}\pi(r_0 - d)^3}{\frac{4}{3}\pi r_0^3}. \quad (14)$$

Combining Eqs. (2), (13) and (14), the mean radius of the calcite grains ($r_0$) can thus be expressed as follows:

$$r_0 = \frac{d}{1 - (1 - Q)^{1/3}} = \frac{2}{\sqrt{\pi}} \frac{\sqrt{D_{\text{foram}} t_{\text{xp}}}}{1 - (1 - Q)^{1/3}} = \frac{2}{\sqrt{\pi}} \frac{\sqrt{D_{0,\text{foram}} \exp\left(-\frac{Ea_{\text{foram}}}{RT_{\text{xp}}}\right) t_{\text{xp}}}}{1 - (1 - Q)^{1/3}}. \quad (15)$$

Given the actual calcite grain size of foraminifera (50–250 nm)²⁹,³⁰, i.e. $r_0$ values of 25–125 nm, the quantity of ¹⁸O that diffused within the foraminifera tests during the experiments can be used to estimate the activation energy and the intrinsic diffusion coefficient of oxygen solid-state diffusion in foraminifera calcite ($Ea_{foram}$ and $D_{0,foram}$). A large range of several tens of kJ mol⁻¹ exists in the literature for the activation energy of oxygen solid-state diffusion in calcite aggregates. Farver and

Yund²⁷ estimated Ea values of 127 ± 17 kJ mol⁻¹ ($D_{0, Farver} = 7.6 \cdot 10^{-9} \text{ m}^2 \text{ s}^{-1}$), while Anderson²⁶ suggested that diffusion within calcites can even proceed with an activation energy as low as 70 kJ mol⁻¹ ($D_{0, Anderson} = 4.6 \cdot 10^{-20} \text{ m}^2 \text{ s}^{-1}$) at low temperatures. Applied to the present experiments, Eq. (15) yields $r_0$ values of 7 nm (i.e. a calcite grain size of ~14 nm, assuming spherical grains) with the values reported by Anderson²⁶, while an $r_0$ value of 8.25 µm (i.e. a calcite grain size of ~16.5 µm) is obtained when using the values reported by Farver and Yund²⁷. This discrepancy from the actual grain size of foraminifera is easily explained by (1) the large uncertainty (typically several tens of kJ mol⁻¹) in the experimental determination of the Ea values and (2) the large dispersion of foraminifera grain sizes and effective grain boundaries, which require intermediate values of $Ea_{foram}$ and $D_{0,foram}$ computed as weighted average defined by:

$$\left\{ \begin{aligned} Ea_{\text{foram}} &= xEa_{\text{Farver}} + (1-x)Ea_{\text{Anderson}} \\ \log(D_{0,\text{foram}}) &= x\log(D_{0,\text{Farver}}) + (1-x)\log(D_{0,\text{Anderson}}) \end{aligned} \right., \quad (16)$$

with x varying between 0 and 1.

Figure 5 shows that x has to be between 0.18 and 0.41 to yield foraminifera calcite grain size values between 50 and 250 nm (i.e. $r_0$ values of 25–125 nm). This range corresponds to $Ea_{foram}$ values ranging from 81 to 94 kJ mol⁻¹ and $D_{0,foram}$ values ranging from $4.8 \times 10^{-18}$ to $1.83 \times 10^{-15} \text{ m}^2 \text{ s}^{-1}$.

**Modelling natural settings.** We developed a numerical model to estimate the present-day δ¹⁸O of fossil tests of benthic foraminifera assumed to have formed during the Cretaceous and the Cenozoic in equilibrium with a seawater/sediment interface at 3.5 °C and that underwent partial isotope re-equilibration via solid-state diffusion during sediment burial. This model was also used to estimate the present-day δ¹⁸O values of fossil tests of planktonic foraminifera assumed to have formed during the Cretaceous and the Cenozoic at different latitudes, i.e. at equilibrium with seawater at different temperatures, and that underwent a similar burial history.

Because the volumes impacted by re-equilibration may not be negligible compared to the size of the foraminifera calcite building blocks, it is preferable to solve the diffusion equation with a spherical geometry rather than with a simpler 1D semi-infinite medium. Because the spherical diffusion equation expressed below does not have a simple analytical solution, the diffusion equation was numerically solved using the finite difference with an implicit time discretisation.

The driving force for isotope re-equilibration through solid-state diffusion over geological times within sediments is the temperature-dependent O isotope fractionation between foraminifera test calcite and seawater at thermodynamic equilibrium ($\Delta_{cc-sw}(\theta)$), where $\theta$ stands for the temperature in °C. A number of formulations of $\Delta_{cc-sw}(\theta)$ can be found in the literature². The most widely used equation is the one reported by Anderson and Arthur³), which relates $\Delta_{cc-sw}$ and $\theta$ as follows:

$$\theta = 16.0 - 4.14(\Delta_{cc-sw}) + 0.13(\Delta_{cc-sw})^2, \quad (17)$$

with

$$\Delta_{cc-sw}(\theta) = \delta^{18}\text{O}_{\text{cc}} - \delta^{18}\text{O}_{\text{sw}}, \quad (18)$$

where $\delta^{18}\text{O}_{\text{cc}}$ is the O isotope composition of the test calcite (in ‰, relative to VPDB) and $\delta^{18}\text{O}_{\text{sw}}$ is the O isotope composition of seawater (in ‰, relative to VSMOW). This equation allows the $\delta^{18}\text{O}_{\text{cc}}$ to be written as a function of $\delta^{18}\text{O}_{\text{sw}}$ and T (in K):

$$\delta^{18}\text{O}_{\text{cc}} = \frac{4.14 - \sqrt{(4.14)^2 - 4*0.13*(16 - T + 273.15)}}{2*0.13} + \delta^{18}\text{O}_{\text{sw}}. \quad (19)$$

Assuming spherical calcite crystals, the mathematical model describing diffusion in spherical coordinates is:

$$\frac{\partial C}{\partial t} = D_{0,\text{foram}} \exp\left(-\frac{Ea_{\text{foram}}}{RT}\right) \frac{1}{r^2} \frac{\partial}{\partial r} \left(r^2 \frac{\partial}{\partial r}\right), \quad (20)$$

with the following initial and boundary conditions in natural settings:

$$\delta^{18}\text{O}_{\text{cc}}(r, t=0) = \Delta_{cc-sw}(T(t=0)) + \delta^{18}\text{O}_{\text{sw}}(t=0), \quad \forall r \in [0; r_0], \quad (21)$$

$$\delta^{18}\text{O}_{\text{cc}}(r = r_0, t) = \Delta_{cc-sw}(T(t)) + \delta^{18}\text{O}_{\text{sw}}(t), \quad (22)$$

where r represents the distance to the solid centre of the spherical foraminifera calcite crystals and $r_0$ represents their radii. Equation (20) was numerically solved using the implicit finite differences method⁴⁰ with M cells numbered from the solid/fluid interface to the solid centre. The numerical discretisation leads to the following mass balance equation:

$$\begin{aligned} & -\frac{S_{i,j+1}D(z)}{dr} \delta^{18}\text{O}_{cc-i}^{n+1} + \left( \frac{V_i}{\Delta t} + \frac{S_{i,j+1}D(z)}{dr} + \frac{S_{i,j+1}D(z)}{dr} \right) \delta^{18}\text{O}_{cc_i}^{n+1} \\ & - \frac{S_{i,j+1}D(z)}{dr} \delta^{18}\text{O}_{cc_i}^{n+1} = \frac{V_i}{\Delta t} \delta^{18}\text{O}_{cc_i}^{n} \end{aligned} \quad (23)$$
where $S_{i,j-1} = 4\pi(r - dr/2)^2$ (respectively $S_{i,j+1}$) is the surface area between the ith and (i-1)th (respectively the ith and (i+1)th) spherical shells of calcite with thicknesses of dr, located at a distance r-dr from the centre. $\delta^{18}\text{O}_{cc_{i-1}}^{n+1}$, $\delta^{18}\text{O}_{cc_{i+1}}^{n+1}$ and $\delta^{18}\text{O}_{cc_{i+1}}^{n}$ represent the $^{18}$O concentrations within the (i-1)th, ith and (i+1)th shells at time step n+1; and $\delta^{18}\text{O}_{cc_i}^n$ is the $^{18}$O concentration within the ith shell at the time step n. $V_i$ is the volume of the ith shell defined as $\frac{4}{3}\pi(r + dr/2)^3 - (r - dr/2)^3$, $\Delta t$ is the time interval and $D(z)$ is the diffusion coefficient at a given burial depth z.

Applying the appropriate boundary conditions permits the calculation of the mass balance for the cell of the mesh in contact with the pore water:

$$V_1 \frac{\delta^{18}\text{O}_{cc_{i-1}}^{n+1} - \delta^{18}\text{O}_{cc_{i-1}}^n}{\Delta t} = S_0 D(z) \frac{\delta^{18}\text{O}_{cc_{i-1}}^{n+1} - \delta^{18}\text{O}_{cc_{i-1}}^{n+1}}{\frac{dr}{\Delta t}} + S_{1,2} D(z) \frac{\delta^{18}\text{O}_{cc_{i+1}}^{n+1} - \delta^{18}\text{O}_{cc_{i+1}}^{n+1}}{\frac{dr}{\Delta t}}, \quad (24)$$

with $S_0$ representing the surface area of the calcite crystals of radius $r_0$.

Similarly, the mass balance can be calculated for the Mth cell corresponding to the centre of a calcite sphere by supposing that the inner interface is impermeable, which sets a no flux boundary condition:

$$V_M \frac{\delta^{18}\text{O}_{cc_M}^{n+1} - \delta^{18}\text{O}_{cc_M}^n}{\Delta t} = S_{M,M-1} D(z) \frac{\delta^{18}\text{O}_{cc_{M-1}}^{n+1} - \delta^{18}\text{O}_{cc_M}^{n+1}}{dr}. \quad (25)$$

Finally, the overall present-day isotope composition of the foraminifera ($\delta^{18}\text{O}_{\text{overall}}(t_{\text{age}})$) is calculated via the cumulative isotope composition of calcite on each cell in the domain $(t_{\text{age}}$ being the time needed to reach present day):

$$\delta^{18}\text{O}_{\text{overall}}(t_{\text{age}}) = \frac{\sum_{i=1}^{M} V_i \delta^{18}\text{O}_{cc_i}(t_{\text{age}})}{(4/3)\pi r_0^3}. \quad (26)$$

**Input parameters of the simulations.** The main parameters of the present simulations were the seawater temperature, the O isotope composition of seawater, the sediment burial rate, the geothermal gradient and the activation of oxygen diffusion in foraminifera calcite (as demonstrated below, the present simulations do not depend on the intrinsic diffusion coefficient).

Here, benthic foraminifera were assumed to have grown in equilibrium with the contemporary deep seawater at a constant temperature of 3.5 °C, and the planktonic foraminifera were assumed to have grown in equilibrium with the surface seawaters at -2, 2, 7, 13, 20, 25, 27 and 27.5 °C, roughly corresponding to the surface waters of the modern ocean at 70, 60, 50, 40, 30, 20, 10 and 0° latitude¹⁰. In other words, we assumed that the ocean temperatures remained constant during the late Cretaceous and the entire Cenozoic.

The O isotope composition of seawater has varied through geologic times, primarily as a consequence of ice sheet growth and decay⁶⁻⁹. The calculation of the influence of the global ice volume on the average seawater $\delta^{18}\text{O}$ depends on how much continental ice was on the planet at a given time. Zachos et al.⁷ suggested that the growth of the ice sheets on Antarctica and in the Northern Hemisphere have been responsible for a total decline of ~2.3‰ of the seawater $\delta^{18}\text{O}$ values (1.2‰ for Antarctic ice sheets and 1.1‰ for Northern Hemisphere ice sheets). However, a value of -1.2‰ for an ice-free ocean is commonly encountered in the literature². Models of present-day ice sheet growth have indicated that an ice-free ocean would have an average $\delta^{18}\text{O}$ value between -0.89 and -1.1‰ (VSMOW)⁴¹,⁴². Assuming a more or less ice-free ocean over the late Cretaceous and the entire Cenozoic, the $\delta^{18}\text{O}$ of the seawater was fixed to -1‰ (VSMOW) for the present simulations. The sediment pore water in which the foraminifera have been fossilised was assumed to be isotopically similar to the seawater in which they lived³² (as supported by the $\delta^{18}\text{O}$ close to -1‰ (VSMOW) of the pore water in the old sediments).

To properly estimate the burial rate and consider the inevitable sediment compaction, the burial depth of the sediments from the DSDP/ODP/IODP sites used in the Friedrich et al.⁹ compilation of the deep sea benthic $\delta^{18}\text{O}$ values have been plotted against their ages, excluding sites that have undergone significant erosion. A sliding five-point average of these data is reported in Fig. 3a with a reasonable fit with $z_{\text{max}} = 500$ m, $\alpha = 10$ and $\beta = 1.5$:

$$z = z_{\text{max}} \left(1 - \exp\left(-\frac{t}{\alpha}\right)\right)^{\beta}. \quad (27)$$

The geothermal gradients in oceanic sediments have been estimated to range between 30 and 70 °C km⁻¹³². Because the temperature experienced by the fossil foraminifera directly controls the extent of solid-state diffusion, the present model was computed using different geothermal gradients, ranging from 40 to 60 °C km⁻¹. The temperature experienced by the sediments $\theta$ (in °C) is the product of the burial depth $z$ (in metres) and the geothermal gradient $V\theta$ (in °C m⁻¹):

$$\theta = \nabla\theta * z \left(1 - \exp\left(-\frac{t}{\alpha}\right)\right)^{\beta}. \quad (28)$$

As detailed above, $E_{\text{foram}}$ values ranging from 81 to 94 kJ mol⁻¹ had to be considered. The present numerical simulations were thus conducted using different activation energies ranging from 85 to 95 kJ mol⁻¹.

**Non-dependence on the intrinsic diffusion coefficient.** The present numerical simulations do not depend on the intrinsic diffusion coefficient ($D_{0,\text{foram}}$), as demonstrated here. If calcite crystals are approximated by simple rods, the Q ratio corresponds to:

$$Q = \frac{d}{r_0} = \frac{2}{r_0 \sqrt{\pi}} \sqrt{D_{0,\text{foram}} t_{\text{xp}}}. \quad (29)$$

where $r_0$ is the radius of the calcite domain. Combining Eq. (29) with Eqs. (2) and (13) yields the following:

$$r_0 = \frac{2}{Q\sqrt{\pi}} \sqrt{D_{0,\text{foram}} \exp\left(-\frac{E_{\text{foram}}}{RT_{\text{xp}}}\right) t_{\text{xp}}}. \quad (30)$$

As emphasised here, the integral of any diffusion profile can be approximated via a step function with length d and height $C_0$. Therefore, at any time, the instantaneous isotope composition of the calcites ($\delta^{18}\text{O}_{\text{cc}}(t)$) can be expressed using the following mass balance:

$$\delta^{18}\text{O}_{\text{cc}}(t) = \frac{\delta^{18}\text{O}_{\text{cc}}(r=r_0, t)d + \delta^{18}\text{O}_{\text{cc}}(r, t=0)(r_0-d)}{r_0}, \quad (31)$$

which yields, after rearrangement,

$$\delta^{18}\text{O}_{\text{cc}}(t) = \frac{d}{r_0} \left( \delta^{18}\text{O}_{\text{cc}}(r_0, t) - \delta^{18}\text{O}_{\text{cc}}(r, 0) \right) + \delta^{18}\text{O}_{\text{cc}}(r, 0). \quad (32)$$

Therefore, after a given duration $t_{\text{age}}$, the isotope composition of the foram ($\delta^{18}\text{O}_{\text{cc}}(t_{\text{age}})$) as previously defined is as follows:

$$\delta^{18}\text{O}_{\text{cc}}^{\text{overall}}(t_{\text{age}}) = \frac{1}{t_{\text{age}}} \frac{1}{r_0} \int_{0}^{t_{\text{age}}} d(\tau) (\delta^{18}\text{O}_{\text{cc}}(r_0, \tau) - \delta^{18}\text{O}_{\text{cc}}(r, 0)) d\tau + \delta^{18}\text{O}_{\text{cc}}(r, 0). \quad (33)$$

Replacing $r_0$ and $d(\tau)$ with the values taken from Eqs. (2), (13) and Eq. (30) yields the following:

$$\delta^{18}\text{O}_{\text{cc}}^{\text{overall}}(t_{\text{age}}) = \frac{1}{t_{\text{age}}} \frac{1}{\frac{2}{Q\sqrt{\pi}} \sqrt{D_{0,\text{foram}} \exp\left(-\frac{E_{\text{foram}}}{RT_{\text{xp}}}\right) t_{\text{xp}}}} \int_{0}^{t_{\text{age}}} \frac{2}{\sqrt{\pi}} \sqrt{D_{0,\text{foram}} \exp\left(-\frac{E_{\text{foram}}}{RT_{\text{xp}}}\right)} (\delta^{18}\text{O}_{\text{cc}}(r_0, \tau) - \delta^{18}\text{O}_{\text{cc}}(r, 0)) d\tau + \delta^{18}\text{O}_{\text{cc}}(r, 0), \quad (34)$$

which yields, after rearrangement,

$$\delta^{18}\text{O}_{\text{cc}}^{\text{overall}}(t_{\text{age}}) = \frac{1}{t_{\text{age}}} \frac{1}{\frac{1}{Q\sqrt{\pi}} \sqrt{\exp\left(-\frac{E_{\text{foram}}}{RT_{\text{xp}}}\right) t_{\text{xp}}}} \int_{0}^{t_{\text{age}}} \sqrt{\exp\left(-\frac{E_{\text{foram}}}{RT(\tau)}\right)} (\delta^{18}\text{O}_{\text{cc}}(r_0, \tau) - \delta^{18}\text{O}_{\text{cc}}(r, 0)) d\tau + \delta^{18}\text{O}_{\text{cc}}(r, 0). \quad (35)$$

Thus, $\delta^{18}\text{O}_{\text{cc}}^{\text{overall}}(t_{\text{age}})$ does not depend on $D_{0,\text{foram}}$.

This result holds true in the case of spherical calcites, at least when the characteristic length of diffusion is negligible compared to the calcite grain size, which is the case for most of the present simulations. In fact, rewriting Eq. (13) and Eq. (15) yields the following:

$$d = \frac{2}{\sqrt{\pi}} \sqrt{D_{0,\text{foram}} \exp\left(-\frac{E_{\text{foram}}}{RT(t)}\right) t} = \sqrt{D_{0,\text{foram}} f(t)} \quad (36)$$

$$\text{and} \quad r_0 = \frac{2}{\sqrt{\pi}} \sqrt{\frac{D_{0,\text{foram}} \exp\left(-\frac{E_{\text{foram}}}{RT_{\text{xp}}}\right) t_{\text{xp}}}{1 - (1 - Q)^{1/3}}} = \sqrt{D_{0,\text{foram}} C t_{\text{exp}}}, \quad (37)$$

in which all the constant parameters of Eq. (15) but $\sqrt{D_{0,\text{foram}}}$ are included in $C_{\text{texp}}$.

If the characteristic length of diffusion is negligible compared to the calcite grain size, the instantaneous isotope composition of calcite ($\delta^{18}\text{O}_{\text{cc}}(t)$) can be written as follows:

$$\begin{aligned} \delta^{18}\text{O}_{\text{cc}}(t) &= \frac{\delta^{18}\text{O}_{\text{cc}}(r_0, t) V_{\text{shell}} + \delta^{18}\text{O}_{\text{cc}}(r, 0) V_{\text{int}}}{V_0} \\ &= \frac{\delta^{18}\text{O}_{\text{cc}}(r_0, t) 4\pi(r_0 - d)^3 d + \delta^{18}\text{O}_{\text{cc}}(r, 0) \frac{d}{3} (r_0 - d)^3}{\frac{d}{3} r_0^3}, \end{aligned} \quad (38)$$

where $V_{\text{shell}}$ stands for a fully re-equilibrated external shell of calcite crystals, $V_{\text{int}}$ is the internal volume unaffected by the re-equilibration and $V_0$ is the initial volume.
Replacing *d* and *r*₀ by their values yields the following:

$$ \delta^{18}\text{O}_{\text{r}}(t) = \frac{\delta^{18}\text{O}_{\text{r}}(r_0, t)4\pi(\sqrt{D_{\text{benth}}C\text{te}_{\text{exp}}} - \sqrt{D_{\text{benth}}f(t)})^2\sqrt{D_{\text{benth}}f(t)} + \delta^{18}\text{O}_{\text{r}}(r, 0)\sqrt{(\sqrt{D_{\text{benth}}C\text{te}_{\text{exp}}} - \sqrt{D_{\text{benth}}f(t)})^3}}{\frac{3}{8}(\sqrt{D_{\text{benth}}C\text{te}_{\text{exp}}})} \quad (39) $$

which yields, after rearrangement,

$$ \delta^{18}\text{O}_{\text{cc}}(t) = \frac{\left( \delta^{18}\text{O}_{\text{cc}}(r_0, t)(\text{Cte}_{\text{exp}} - f(t))^2 f(t) + \frac{1}{3}\delta^{18}\text{O}_{\text{cc}}(r, 0)(\text{Cte}_{\text{exp}} - f(t))^3 \right)}{\frac{1}{3}\text{Cte}_{\text{exp}}^3} \quad (40) $$

After a given duration $t_{\text{age}}$, the isotope composition of a foraminifera ($\delta^{18}\text{O}_{\text{overall}}(t_{\text{age}})$) is written as follows:

$$ \delta^{18}\text{O}_{\text{cc}}(t_{\text{age}}) = \frac{\frac{2}{t_{\text{age}}} \int_{0}^{t_{\text{age}}} \left( \delta^{18}\text{O}_{\text{cc}}(r_0, \tau)(\text{Cte}_{\text{exp}} - f(t))^2 f(t) + \frac{1}{3}\delta^{18}\text{O}_{\text{cc}}(r, 0)(\text{Cte}_{\text{exp}} - f(t))^3 \right) d\tau}{\frac{1}{3}\delta^{18}\text{O}_{\text{cc}}(r, 0)(\text{Cte}_{\text{exp}} - f(t))^3} \quad (41) $$

Thus, $\delta^{18}\text{O}_{\text{overall}}(t_{\text{age}})$ does not depend on $D_{0,\text{foram}}$ in the case of spherical calcites.

**Data availability.** Original data and numerical codes are available from the corresponding author upon request.

Received: 17 January 2017 Accepted: 25 August 2017

Published online: 26 October 2017

## References

1. Urey, H. C. Oxygen isotopes in nature and in the laboratory. *Science* **5**, 489–496 (1948).

2. Pearson, P. N. Oxygen isotopes in foraminifera: overview and historical review. *Paleontol. Soc. Pap.* **18**, 1–38 (2012).

3. Anderson, T. F. & Arthur, M. A. Stable isotopes of oxygen and carbon and their application to sedimentologic and paleoenvironmental problems. In Arthur, M. A., Anderson, T.F., Kaplan, I.R., Veizer, J., Land, L.S. (eds.) *Stable Isotopes in Sedimentary Geology, SEPM Short Course*, Vol. 10, Ch. 1, 1–15 (1983).

4. Chappell, J. & Shackleton, N. J. Oxygen isotopes and sea level. *Nature* **324**, 137–140 (1986).

5. Spero, H. J., Bijma, J., Lea, D. W. & Bemis, B. E. Effect of seawater carbonate concentration on foraminiferal carbon and oxygen isotopes. *Nature* **390**, 497–500 (1997).

6. Savin, S. M. The history of the Earth's surface temperature during the past 100 million years. *Annu. Rev. Earth Planet. Sci.* **5**, 319–355 (1977).

7. Zachos, J., Pagani, M., Sloan, L., Thomas, E. & Billups, K. Trends, rhythms, and aberrations in global climate 65 Ma to present. *Science* **292**, 686–693 (2001).

8. Cramer, B. S., Toggweiler, J. R., Wright, J. D., Katz, M. W. & Miller, K. G. Ocean overturning since the Late Cretaceous: inferences from a new benthic foraminiferal isotope compilation. *Paleoceanography* **24**, PA4216 (2009).

9. Friedrich, O., Norris, R. D. & Erbacher, J. Evolution of middle to Late Cretaceous oceans – a 55 My record of Earth’s temperature and carbon cycle. *Geology* **40**, 107–110 (2012).

10. Zachos, J. C., Stott, L. D. & Lohmann, K. C. Evolution of early Cenozoic marine temperatures. *Paleoceanography* **9**, 353–387 (1994).

11. D'Hondt, S. & Arthur, M. A. Late cretaceous oceans and the cool tropic paradox. *Science* **271**, 1838–1841 (1998).

12. Killingley, J. S. Effects of diagenetic recrystallization on $^{18}$O/$^{16}$O values of deep-sea sediments. *Nature* **301**, 594–597 (1983).

13. Schrag, D. P. Effects of diagenesis on the isotopic record of late Paleogene tropical sea surface temperatures. *Chem. Geol.* **161**, 215–224 (1999).

14. Kozdon, R. et al. In situ $\delta^{18}$O and Mg/Ca analyses of diagenetic and planktic foraminiferal calcite preserved in a deep-sea record of the Paleocene–Eocene thermal maximum. *Paleoceanography* **28**, 517–528 (2013).

15. Edgar, K. M., Anagnostou, E., Pearson, P. N. & Foster, G. L. Assessing the impact of diagenesis on $\delta^{11}$B, $\delta^{13}$C, $\delta^{18}$O, Sr/Ca and B/Ca values in fossil planktic foraminiferal calcite. *Geochim. Cosmochim. Acta* **166**, 189–209 (2015).

16. Lohmann, G. P. A model for variation in the chemistry of planktonic foraminifera due to secondary calcification and selective dissolution. *Paleoceanography* **10**, 445–457 (1995).

17. Pearson, P. N. et al. Warm tropical sea surface temperatures in the Late Cretaceous and Eocene epochs. *Nature* **413**, 481–487 (2001).

18. Wilson, P. A. & Norris, R. D. Warm tropical ocean surface and global anoxia during the mid-Cretaceous period. *Nature* **412**, 425–429 (2001).

19. Huber, M. & Sloan, L. C. Heat transport, deep waters, and thermal gradients: coupled simulation of an Eocene greenhouse climate. *Geophys. Res. Lett.* **28**, 3481–3484 (2001).

20. Valdes, P. Built for stability. *Nat. Geosci.* **4**, 414–416 (2011).

21. Urey, H. C., Lowenstam, A., Epstein, S. & McKinney, C. R. Measurement of paleotemperatures and temperatures of the Upper Cretaceous of England, Denmark, and the Southeastern United States. *Bull. Geol. Soc. Am.* **62**, 399–416 (1951).

22. Mavromatis, V. et al. The continuous re-equilibration of carbon isotope compositions of hydrous Mg carbonates in the presence of cyanobacteria. *Chem. Geol.* **404**, 41–51 (2015).

23. Hoppe, P., Cohen, S. & Meibom, A. NanoSIMS: technical aspects and applications in cosmochemistry and biological geochemistry. *Geostand. Geoanal. Res.* **37**, 111–134 (2013).

24. Ruiz-Agudo, E., Putnis, C. V. & Putnis, A. Coupled dissolution and precipitation at mineral-fluid interfaces. *Chem. Geol.* **383**, 132–146 (2014).

25. Putnis, A. Why mineral interfaces matter. *Science* **343**, 1441–1442 (2014).

26. Anderson, T. F. Self-diffusion of carbon and oxygen in calcite by isotope exchange with carbon dioxide. *J. Geophys. Res.* **74**, 3918–3932 (1969).

27. Farver, J. R. & Yund, R. A. Oxygen grain boundary diffusion in natural and hot-pressed calcite aggregates. *Earth Planet. Sci. Lett.* **161**, 189–200 (1998).

28. Stipp, S. L. S., Konnerup-Madsen, J., Franzreb, K., Kulik, A. & Mathieu, H. J. Spontaneous movement of ions through calcite at standard temperature and pressure. *Nature* **396**, 356–359 (1998).

29. Cuif, J. P., Dauphin, Y. & Sorauf, J. E. *Biominerals and Fossils Through Time* (Cambridge University Press Editions, 2011).

30. Iwasaki, S. et al. Observation of the dissolution process of *Globigerina bulloides* tests (planktic foraminifera) by X-ray microcomputed tomography. *Paleoceanography* **30**, 317–331 (2015).

31. Malinverno, A. & Martinez, E. A. The effect of temperature on organic carbon degradation in marine sediments. *Sci. Rep.* **5**, 17861 (2015).

32. Schrag, D. P., Hampt, G. & Murray, D. W. Pore fluid constraints on the temperature and oxygen isotopic composition of the glacial ocean. *Science* **272**, 1930–1932 (1996).

33. Herbert, T. D., Peterson, L. C., Lawrence, K. T. & Liu, Z. Tropical ocean temperatures over the past 3.5 million years. *Science* **328**, 1350–1354 (2010).

34. van der Lee, J. & De Windt, L. CHESS Tutorial and Cookbook. Updated for Version 3.0. *Manual Nr. LMH/RD/02/13, Paris* **116** (2002).

35. Wolery, T. EQ3/6: Software package for geochemical modeling of aqueous systems: Package overview and installation guide (version 7.0). *Report UCRL-MA-210662*. Lawrence Livermore National Laboratory, Livermore, CA, USA (1992).

36. Samson, E. G., Lemaire, J. & Beaudoin, J. J. Modeling chemical activity effects in strong ionic solutions. *Comput. Mater. Sci.* **15**, 285–294 (1999).

37. Cole, D. R. & Chakraborty, S. Rates and mechanisms of isotopic exchange. *Rev. Mineral. Geochem.* **43**, 83–223 (2001).

38. Watson, E. B. & Baxter, E. F. Diffusion in solid-Earth systems. *Earth Planet. Sci. Lett.* **253**, 307–327 (2007).

39. Crank, J. *The Mathematics of Diffusion* 2nd edn (Oxford University Press, 1975).

40. LeVeque, R. J. Finite difference methods for ordinary and partial differential equations: steady-state and time-dependent problems. *SIAM Publ.* **98**, 328 (2007).

41. Lhomme, N., Clarke, G. K. C. & Ritz, C. Global budget of water isotopes inferred from polar ice sheets. *Geophys. Res. Lett.* **32**, L205002 (2005).

42. Cramer, B. S., Miller, K. G., Barrett, P. J. & Wright, J. D. Late Cretaceous–Neogene trends in deep ocean temperature and continental ice volume: reconciling records of benthic foraminiferal geochemistry ($\delta^{18}$O and Mg/Ca) with sea level history. *J. Geophys. Res.* **116**, C12023 (2011).

## Acknowledgements

This work was supported by the European Research Council Advanced Grant 246749 BIOCARB to A.M. The SEM facility of the IMPMC is supported by Region Ile de France grant SESAME NI-07-593/R, INSU-CNRS, INP-CNRS, UPMC-Paris 6, and by the Agence Nationale de la Recherche (ANR) Grant No. ANR-07-BLAN-0124-01. The National NanoSIMS Facility at the MNHN is supported by the CNRS, the Région Ile de France, the Ministère délégué à l'enseignement supérieur et à la recherche, and the MNHN. Special thanks go to O. Lengléré for stimulative discussions and to five anonymous reviewers for their constructive comments.

## Author contributions

S.B., D.D. and A.M. designed the experimental study; S.B. and D.D. performed the isotope re-equilibration experiments; S.B. and S.P. performed the SEM observations; S.B. and A.M. performed the NanoSIMS imaging; D.D. and P.A. developed analytical and
numerical models with critical inputs from S.B.; S.B., D.D. and A.M. wrote the manuscript.

**Additional information**

**Supplementary Information** accompanies this paper at https://doi.org/10.1038/s41467-017-01225-9.

**Competing interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017