Nature Communications

Nat Commun. 9: 4954

Major intensification of Atlantic overturning circulation at the onset of Paleogene greenhouse warmth

http://orcid.org/0000-0002-4076-1248S. J. Batenburg

sietske.batenburg@earth.ox.ac.uk

¹², S. Voigt¹, O. Friedrich³, A. H. Osborne⁴, http://orcid.org/0000-0003-4341-2366A. Bornemann⁵, T. Klein¹, L. Pérez-Díaz⁶, M. Frank⁴

1. 0000 0004 1936 9721grid.7839.5Institut für Geowissenschaften, Goethe-Universität Frankfurt, Altenhöferallee 1, Frankfurt am Main, 60438 Germany

2. 0000 0004 1936 8948grid.4991.5Department of Earth Sciences, University of Oxford, South Parks Road, Oxford, OX1 3AN UK

3. 0000 0001 2190 4373grid.7700.0Institut für Geowissenschaften, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 234-236, 69120 Heidelberg, Germany

4. 0000 0000 9056 9663grid.15649.3fGEOMAR Helmholtz-Zentrum für Ozeanforschung Kiel, Wischhofstr. 1-3, Kiel, 24148 Germany

5. 0000 0001 2155 4756grid.15606.34Bundesanstalt für Geowissenschaften und Rohstoffe, Stilleweg 2, 30655 Hannover, Germany

6. 0000 0001 2161 2573grid.4464.2Department of Earth Sciences, Royal Holloway, University of London, Egham, TW20 0EX UK

© The Author(s) 2018

DOI: 10.1038/s41467-018-07457-7

Published online: 23 November 2018

Abstract

During the Late Cretaceous and early Cenozoic the Earth experienced prolonged climatic cooling most likely caused by decreasing volcanic activity and atmospheric CO₂ levels. However, the causes and mechanisms of subsequent major global warming culminating in the late Paleocene to Eocene greenhouse climate remain enigmatic. We present deep and intermediate water Nd-isotope records from the North and South Atlantic to decipher the control of the opening Atlantic Ocean on ocean circulation and its linkages to the evolution of global climate. The marked convergence of Nd-isotope signatures 59 million years ago indicates a major intensification of deep-water exchange between the North and South Atlantic, which coincided with the turning point of deep-water temperatures towards early Paleogene warming. We propose that this intensification of Atlantic overturning circulation in concert with increased atmospheric CO₂ from continental rifting marked a climatic tipping point contributing to a more efficient distribution of heat over the planet.

Abstract

Ocean circulation likely played a major role in past episodes of global warming. Here the authors show that converging neodymium isotope values indicate intensification of the Atlantic overturning circulation at 59 million years ago, which may have contributed to a more efficient distribution of heat over the planet.

Introduction

The Earth underwent long-term climatic cooling between the peak-greenhouse intervals of the mid-Cretaceous and the Eocene¹–⁵. Globally averaged deep-water temperatures gradually declined by almost 10 °C from 72 to 59 Ma, as estimated from benthic foraminiferal oxygen-isotope data³,⁶. This cooling has been ascribed to decreasing atmospheric CO₂ levels⁷–⁹ through global reduction of volcanism and sea-floor spreading rates¹⁰ combined with changes in ocean circulation patterns³. In contrast, there is no comprehensive model explaining how the greenhouse conditions of the Eocene were established and what the roles of atmospheric CO₂ and ocean circulation were in promoting global warming. Mechanisms proposed so far have solely focussed on increased atmospheric CO₂ levels either induced by carbon cycle changes⁶, rates of continental rifting¹¹, or by enhanced volcanism of the North Atlantic igneous province¹²,¹³. The role of changes in overturning circulation caused by the opening of the Atlantic Ocean and related changes in oceanic heat transport has, however, not been addressed yet.

While circum-equatorial flow, which had dominated circulation in the proto-North Atlantic earlier in the Cretaceous, gradually declined¹⁴, the ongoing opening and deepening of the Atlantic basin¹⁵,¹⁶ led to increased North-South connectivity, although the timing of the establishment of a deep-water connection remains debated¹⁷–²². Enhanced latitudinal water-mass exchange likely promoted the distribution of heat across the planet via the thermohaline conveyor and resulted in reduced temperature contrasts between the equator and the poles. To distinguish tectonic constraints on circulation from climatically driven changes, the role of subsiding submarine barriers has to be assessed. We determine the timing of the establishment of a persistent deep-water connection between the North and South Atlantic by combining deep-water neodymium (Nd) isotope and temperature records.

Assessing the role of ocean circulation on Earth’s climate in the latest Cretaceous and early Paleogene requires tight constraints on the modes and locations of deep-water formation and the extent of mixing of different deep-water masses. Information on past water mass mixing and exchange can be derived from Nd-isotope signatures (¹⁴³Nd/¹⁴⁴Nd, expressed as εNd(t)) of authigenic, seawater-derived sedimentary archives such as ferromanganese coatings of sediment particles or fish debris, which have been demonstrated to incorporate the Nd-isotope composition of ambient deep waters²³. Deep-water masses mainly acquire their Nd-isotope signatures from continental contributions via rivers and dust inputs in their source areas²³, as well as through exchange processes with ocean margin sediments²⁴. These characteristic Nd-isotope compositions of deep-water masses are then conservatively advected and mixed over large distances in the open ocean given that the average Nd residence time of 400–2000 years is similar to the global ocean mixing time²³,²⁵. Analysis of the Nd-isotope composition of authigenic sedimentary archives thus allows the reconstruction of changes in deep-water mixing over time.

Existing Late Cretaceous Atlantic seawater εNd(t) signatures display a large spread in values (Fig. 1) that led to the suggestion that different mechanisms and locations of deep-water formation operated simultaneously²¹,²⁶–³⁰. There are indications that intermediate and deep-water exchange commenced as early as at 90 Ma in the Late Turonian²⁰, although the deep Atlantic Ocean potentially operated as a number of sub-basins with limited connectivity until the Maastrichtian²². The large variability in Cretaceous εNd(t) values has so far been interpreted to reflect different modes and locations of deep-water formation either in the southern high latitudes²¹,²²,²⁹,³⁰ or in the North Atlantic²⁶–²⁸, or local deep-water formation in relatively shallow sub-basins separated at depth²². The spread in εNd(t) values was likely further enhanced by local boundary exchange and weathering inputs into the relatively small Atlantic basins. Since the Cretaceous Atlantic Ocean was limited in depth and width, its contact area with the margins was large compared to its volume. Regional processes such as boundary exchange thus had a profound effect on water-column chemistry, as is the case for modern near-shore settings²⁴,³¹ or restricted sub-basins³². The Cretaceous Nd-isotope records from Demerara Rise and the Cape Verde Basin exemplify this effect by local weathering inputs of highly unradiogenic Nd from the old cratons of South America and Africa²⁷,²⁸ (Fig. 1). Despite the potential influence of continental inputs near ocean margins, several open-ocean sites in the Late Cretaceous North and South Atlantic show parallel Nd-isotope trends. These parallel trends have been interpreted to reflect the formation and northward flow of a southern-sourced deep-water mass, “Southern Component Water”³⁰, although the behaviour of individual εNd(t) records is highly variable on a time scale of millions of years and patterns of change are dissimilar between localities.

Fig. 1View larger versionFig. 1. Nd-isotope data across the latest Cretaceous – early Paleogene. In this compilation of published Atlantic εNd(t) records and new Nd-isotope data from Site 369 (this study) only records with four or more data points over the time-interval 72–50 Ma are included. The εNd(t) data of Demerara Rise (Sites 1258, 1260, 1261), Maud Rise, Blake Nose and Site 525 have been smoothed by LOWESS regression (f = 0.3) with the dashed lines representing the 97.5% confidence interval. The εNd(t) data are flanked by the geomagnetic polarity time scale with chrons C31-C24. PETM: Paleocene/Eocene thermal maximum, K/Pg: Cretaceous/Paleogene Boundary

From the Paleocene–Eocene Thermal Maximum (PETM) at 56 Ma onwards, most open-ocean εNd(t) signatures from the North and South Atlantic were within a narrow range of -8 to -10²⁸,³³,³⁴, indicating common water masses at bathyal and abyssal depths. There is, however, a lack of data for the Paleocene, which limits our understanding of when and to what extent deep waters exchanged and when the Atlantic started to play a key role in hemispheric oceanic heat exchange. A compilation of existing εNd(t) records for the period of time from 72 to 50 Ma (Fig. 1) shows that Nd-isotope data for the Paleocene are only available from a limited number of sites and, with the exception of Demerara Rise²⁶,²⁸, are of limited resolution (less than one sample per two million years). The Paleocene, however, marks the time when the Atlantic significantly widened and deepened, which potentially paved the way for similar-to-modern ocean overturning processes¹⁷. Here we fill this gap and present new Paleocene intermediate- and deep-water Nd-isotope records from the North and South Atlantic Ocean. Five ocean drilling sites were selected from paleo-water depths between 500 and 4500 m (Supplementary Table 1/Fig. 2) to obtain seawater Nd-isotope records covering the critical time span from the end-Cretaceous to the early Paleogene.

Fig. 2View larger versionFig. 2. Paleogeographic setting. a Paleogeographic reconstruction at 60 Ma. Plate tectonic configuration after ref. ⁶³ and global bathymetry from ref. ³⁹. Circles indicate Atlantic Ocean drilling sites for which Nd-isotope data are available (at least four points in the time-interval 72–50 Ma), DR: Demerara Rise. Diamonds indicate locations for which new data are presented here. Numbers indicate DSDP/ODP/IODP Sites. b Three-dimensional plot of the reconstructed bathymetry of the South Atlantic with the Rio Grande Rise – Walvis Ridge barrier at 60 Ma from ref. ⁴⁰. c Contour plot of the bathymetry in panel B, with black numbers indicating depths in meters

Results

Seawater origin of Nd-isotope signatures

Seawater Nd-isotope signatures were obtained by leaching ferromanganese coatings of bulk sediments that are considered a reliable archive if sufficiently weak leaching procedures are applied³⁵. The εNd(t) variability of the detrital material was also determined for selected samples in this study (details in “methods” section), to evaluate the potential influence of local weathering inputs. The εNd(t) signatures of the detrital fractions and the leached ferromanganese oxide coatings show similar long-term trends at Sites 516 and U1403 and parts of the records at Sites 1267 and 525. Despite following parallel trends, most detrital εNd(t) values are significantly offset from those of the coatings supporting the validity of the seawater εNd(t) signatures extracted from the coatings at the offshore locations of our studied sites as faithful recorders of past water mass mixing (Fig. 3). The Nd-isotope composition of the water-masses themselves may have been influenced to some extent by local factors such as boundary exchange processes that mainly occur when deep-water circulation is slow and/or the sites were located in small or partly isolated basins with high detrital input³¹,³²,³⁶. In addition, the dissolved seawater Nd-isotope signature may have been incorporated into the hydrogenous component of pelagic clays²⁰,³⁷, which may partly explain the similarity in the long-term evolution of the detrital and leached εNd(t) values.

Fig. 3View larger versionFig. 3. New Nd-isotope results. Neodymium isotope values of ferromanganese coatings (coloured symbols) and detrital fractions (open symbols) measured in this study and previously for Site 525²². Error bars indicate 2 s.d. external reproducibility

Parallel trends and convergence of Nd-isotope values

Our new seawater Nd-isotope records from the North and South Atlantic (Fig. 3 and Supplementary Tables 2 to 6) display a wide range of values (−2 to −11) in the Maastrichtian interval (72.1–66 Ma) with parallel trends that converge to a common value of -8 to -9 at 59 Ma (Fig. 4). Our North Atlantic record from Site U1403 ends at 58 Ma, but εNd(t) values between −9.2 and 8 around 57 Ma at northern Site 549³⁸ corroborate our findings (Fig. 1).

Fig. 4View larger versionFig. 4. The convergence of Nd-isotope values. Maastrichtian to early Eocene Nd-isotope records compared to atmospheric CO₂ levels and deep-sea temperature evolution. Open symbols depict new εNd(t) data for North Atlantic Site U1403 (Newfoundland margin, dark blue) and South Atlantic Sites 516 (Rio Grande Rise, red), 1267 (northern flank of the Walvis Ridge, green) and 525 (top of the Walvis Ridge, yellow). Filled symbols depict selected high-resolution εNd(t) data from Demerara Rise in the North Atlantic²⁶,²⁸ (purple), Sites 527³³ and 1262³⁴ on the flank of the Walvis Ridge (green), Site 1263³⁴ and Site 525²² on the Walvis Ridge (yellow), and Maud Rise in the Southern Ocean²⁰ (dark blue). The εNd(t) data have been smoothed by LOWESS regression (f = 0.3) with the shaded areas representing the 97.5% confidence interval. The εNd(t) data are flanked by the geomagnetic polarity time scale on the left and on the right by reconstructed atmospheric CO₂ levels¹¹, benthic oxygen isotope data from the Pacific⁴¹ and the Atlantic¹³, and the suggested evolution of the Atlantic Ocean. The Paleogene age model of the Sites rests on a combination of magnetic and nannofossil stratigraphy plotted against the GTS 2012⁶⁴. C31-C24: paleomagnetic chrons, δ¹⁸Ob: oxygen isotope ratio from benthic foraminifera, VPDB: Vienna Pee Dee Belemnite standard, PETM: Paleocene/Eocene Thermal Maximum, DR: Demerara Rise, K/Pg: Cretaceous/Paleogene boundary, MR: Maud Rise

Sites 525, 1267 and 516 in the South Atlantic, and Site U1403 in the North Atlantic show a trend of decreasing εNd(t) from approximately 70 to 63 Ma, with lowest values reached in the first half of the Paleocene. This decrease may reflect the reduction in active volcanism and exposed volcanic terrains in and around the Atlantic Ocean²⁰. Nd-isotope values at Site 525 were positively offset from εNd(t) signatures at comparably shallow Site 516 on the Rio Grande Rise and nearby deeper Site 1267 at the base of the north-western slope of the Walvis Ridge until the end of the Cretaceous. This positive offset was most likely caused by the weathering influx of volcanic material from the partially subaerially exposed Walvis Ridge in the latest Cretaceous¹⁵,³⁹. The offset decreased as the ridge and Site 525 subsided.

From approximately 64 Ma onwards, average εNd(t) values display an increasing trend until 60–59 Ma. We assign this trend to the enhanced volumetric flow of deep and intermediate water masses in the opening South Atlantic Basin which likely led to a decrease of the influence of local inputs and boundary change effects. In addition, the observed trend coincides with a first phase of magmatic activity of the North Atlantic Igneous Province from 62 to 61 Ma¹³, which may have supplied radiogenic Nd, and ongoing deepening of the study sites that may have reduced unradiogenic weathering inputs from nearby continents.

From 59 Ma onwards, the Nd-isotope signatures at all newly studied sites, as well as Demerara Rise²⁶, decrease together and our εNd(t) results fall within a narrow range of -7 to -9.5 for the period 58.5–56.5 Ma. This convergence may reflect increasing admixture of southern-sourced deep water, which would have carried a εNd(t) signature similar to that at Maud Rise of approximately -9 in the Paleocene (Fig. 4)²⁰.

Opening of the Atlantic Ocean and climatic implications

Recent paleobathymetric reconstructions show that deep oceanic basins in the Atlantic Ocean, like the Cape and the Angola basins, were constricted until the end of the Cretaceous¹⁵. Deeper structures, such as the Vema and Hunter channels flanking Rio Grande Rise only allowed intermediate-water exchange at depths shallower than 2500 m. In the Paleocene, the South Atlantic deepened and widened, with the western portion of the Rio Grande Rise having subsided below 2500 m water depth at 60 Ma and the Argentine and Brazil basins reaching depths of over 5500 m in the early Eocene¹⁵,⁴⁰ (Fig. 2). The close correspondence in εNd(t) signatures at 59 Ma suggests a common deep-water signature (εNd(t) −9 to −8) in the South and North Atlantic (Fig. 4). We interpret the converging trend of Nd-isotope signatures to reflect an increasingly efficient deep-ocean circulation in the Atlantic Ocean with the dominant deep-water masses most likely originating in the high southern latitudes. Such a southern origin of deep water is consistent with recent modelling results¹⁷ suggesting locations of major deep-water formation in the Southern Ocean, potentially supplemented by a minor source of deep water formed offshore North America. At the same time, our data show that the sub-basins of the deep Atlantic became fully connected by subsidence of the Rio Grande Rise near 59 Ma, accompanied by the widening and deepening of the equatorial gateway¹⁷,²⁰. The improved connectivity and the increased volumetric exchange of water masses in the Atlantic Ocean at 59 Ma allowed modern-like open-ocean processes and water-mass mixing to be established, which decreased the sensitivity of the Nd-isotope composition of seawater to local effects such as terrigenous and coastal sedimentary inputs. The convergence of Nd-isotope signatures across the entire Atlantic Ocean spanning paleo-waterdepths of 500 to 4500 m, further suggests that between 62 and 59 Ma, both local tectonic restrictions as well as the vertical stratification of the deep Atlantic Ocean decreased and a global mode of thermohaline circulation was initiated.

The close correspondence in Nd-isotope values among sites at 59 Ma coincided with the onset of the mid-Paleocene global climate warming as evident from benthic foraminiferal oxygen isotopes¹³,⁴¹ (Fig. 4). Based on a recent reconstruction of continental rift length histories¹¹ in comparison to the long-term evolution of atmospheric pCO₂⁸, the underlying cause of this warming may lie in the increased cumulative length of incipient continental rifts. Despite a reconstructed gradual increase in pCO₂ levels during the end of the Cretaceous and earliest Paleocene⁸,¹¹ (Fig. 4), as well as an initial magmatic phase of the North Atlantic Igneous Province from 62 to 61 Ma¹³, the long-term increasing trend in bottom-water temperatures did not start until 59 Ma¹³,⁴¹, when pCO₂ started to increase at a higher rate⁸,¹¹ (Fig. 4).

Discussion

We hypothesize that the strengthened Atlantic overturning circulation suggested by our data enhanced oceanic poleward heat transport thereby contributing to global climate warming culminating in the peak greenhouse conditions of the Eocene. Global warming may itself have enhanced vertical mixing through increased occurrence of storms and cyclones⁴² that enabled more efficient overturning circulation in the Atlantic Ocean. Both the deepening of the Rio Grande Rise and enhanced mixing associated with global warming would have increased the capacity of the overturning circulation in the Atlantic Ocean to transport heat. These interpretations of our new Nd-isotope data are consistent with observed changes in Late Cretaceous to early Paleogene Nd-isotope records from the Pacific Ocean⁴² and Earth system modelling results, which indicate that vigorous ocean circulation and strong vertical mixing resulted in increased oceanic heat transport and reduced equator–pole temperature gradients⁴²,⁴³. Higher oceanic heat transport efficiency likely also set the stage for the occurrence of brief hyperthermals which were frequently superimposed on the overall temperature rise of the Eocene hothouse⁴¹. Together with increasing atmospheric CO₂ levels⁸,¹¹, the changing paleogeography of the Atlantic Ocean may have contributed to the boundary conditions that pushed the Earth’s climate into a greenhouse state.

Methods

Extraction of Nd isotopes

For Nd-isotope analyses of past seawater extracted from ferromanganese oxide coatings, bulk sediment samples consisting mainly of nannofossil oozes and chalks were dried and homogenised in an agate mortar. To extract the authigenic, seawater-derived Nd-isotope signature, ~2.5 g of powder was treated following the procedure described in ref. ⁴⁴, omitting the carbonate removal step⁴⁵. Powdered samples were rinsed three times with de-ionized (MQ) water, after which 10 ml of MQ was added and 10 ml of a 0.05 M hydroxylamine hydrochloride/15% acetic acid solution, buffered with NaOH to a pH of 4. Samples were placed on a shaker for 1 h and centrifuged. The supernatant containing the seawater Nd-isotope signature of the ferromanganese oxide coatings was pipetted off and dried down. For determining the detrital εNd signature, selected samples underwent an additional 12 h leaching step with 20 ml of the hydroxylamine solution (above), after which samples were rinsed with MQ three times and ~50 mg of dried sample was dissolved in a mixture of aqua regia and HF. As preparatory steps for column chemistry, all samples were refluxed in concentrated HNO₃ at 80 °C overnight, centrifuged, and 80% of the supernatant was dried down. Twice, 0.5 ml of 1 M HCl was added and the sample was dried down, after which the samples were redissolved in 0.5 ml 1 M HCl. Samples were passed through cation-exchange columns with 0.8 ml AG50W-X12 resin (mesh size 200‒400 μm), using standard procedures, to separate Sr and the Rare Earth Elements (REEs), as well as removing most of the Ba⁴⁶. A second set of columns with 2 ml Ln-Spec resin (mesh size 50‒100 μm) was used to separate Nd from the other REEs and remaining Ba⁴⁷.

Analytical procedure

Neodymium isotope ratios were measured on a Nu Instruments Multiple Collector Inductively Coupled Plasma Mass Spectrometer (MC-ICPMS). The majority of samples were measured at GEOMAR Kiel, Germany, and a subset of samples at the department of Earth Sciences of Oxford University, UK (Supplementary Tables 2, 5 and 6). Measured ¹⁴³Nd/¹⁴⁴Nd results were mass-bias corrected to a ¹⁴⁶Nd/¹⁴⁴Nd ratio of 0.7219 and were normalized to the accepted ¹⁴³Nd/¹⁴⁴Nd value of 0.512115 for the JNdi-1 standard⁴⁸, which was measured after every third sample.

The results were decay-corrected for the time of deposition by (¹⁴³Nd/¹⁴⁴Nd)sample(t) = (¹⁴³Nd/¹⁴⁴Nd)sample(0) – [(¹⁴⁷Sm/¹⁴⁴Nd)sample(0) * (℮ʎt – 1)] where t is time, the decay constant ʎ is 6.54 × 10⁻¹², and using an average ¹⁴⁷Sm/¹⁴⁴Nd ratio of 0.124²². Nd-isotope ratios are reported as εNd(t) values with respect to the Chondritic Uniform Reservoir (CHUR), which are calculated as εNd(t) = [(¹⁴³Nd/¹⁴⁴Nd)sample(t) / (¹⁴³Nd/¹⁴⁴Nd)CHUR(t) − 1] × 10⁴ using a (¹⁴³Nd/¹⁴⁴Nd)CHUR(0) value of 0.512638, and a (¹⁴⁷Sm/¹⁴⁴Nd)CHUR(0) of 0.1966⁴⁹. External reproducibility (2σ) of the measurements was between 0.15 and 0.54 εNd units and procedural Nd blanks were ≤ 30 pg Nd and thus negligible.

Age models

Age models for the individual sites were generated by an integrated approach of magneto- and biostratigraphy and if available astrochronology. All datum levels are tied to the Geological Timescale GTS2012. Ages of polarity chrons are from ref. ⁵⁰, and of calcareous nannofossils (NP zonation) from ref. ⁵¹ as compiled in ref. ⁵². In detail the following data are used and summarized in Supplementary Table 7. Tie points for Site 516 are defined by magneto- and calcareous nannofossil stratigraphy given in ref. ⁵³. Tie points for Site 525 are defined by polarity chrons in the Maastrichtian⁵⁴ and by calcareous nannofossils in the Paleocene⁵⁵. Tie points for Site 1267 are derived from precession cycle counting for the upper Paleocene (until 58.2 Ma ago)⁵⁶ and polarity chrons for the lower to middle Paleocene and the Maastrichtian⁵⁷. Ages of neodymium isotope data from Site 1262³⁴ and Site 527³³ were converted to GTS 2012. Tie points for Site 369 follow the age model of the Shipboard Scientific Party⁵⁸. Tie points for Site U1403 are defined by first occurrences (FO) of calcareous nannofossils for the Paleocene⁵⁹ with an adjustment for the FO of Lithoptychius spp. at 227 m depth rCCSF (corresponding to the first radiation of fasciculithids according to refs. ⁶⁰,]⁶¹) and by astronomical tuning of 405 kyr cycles and carbon isotope stratigraphy⁶².

Supplementary Material

Open In Web Browser

Footnotes

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

Electronic supplementary material

Supplementary Information accompanies this paper at 10.1038/s41467-018-07457-7.

Acknowledgements

We would like to thank the International Ocean Discovery Program (IODP) for providing samples. Authors were funded by the German Research Foundation (DFG) under grant numbers DFG VO 687/14, FR2544/8, and FR1198/11, BO2505/8 and EH89/20. We would like to thank Jutta Heinze and Chris Siebert at GEOMAR, Kiel and Alan Hsieh at Oxford University, UK for smooth operation of the laboratory and the mass spectrometers.

Author contributions

S.J.B., S.V., O.F. and M.F. developed the project. S.J.B., A.H.O. and T.K. performed Nd-isotope analyses. L.P-D. provided bathymetric information. All authors contributed to the writing of the manuscript.

Data availability

The authors declare that all the data generated during this study are available within the manuscript and its supplementary information file.

Notes

Competing interests

The authors declare no competing interests.

Articles from Nature Communications are provided here courtesy of Nature Publishing Group

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Clarke LJ, Jenkyns HC, authors. New oxygen isotope evidence for long-term Cretaceous climatic change in the Southern Hemisphere. Geology. 1999;27:699–702. [CrossRef]

Cramer BS, Toggweiler JR, Wright JD, Katz ME, Miller KG, authors. Ocean overturning since the Late Cretaceous: Inferences from a new benthic foraminiferal isotope compilation. Paleoceanography. 2009;24:PA4216 [CrossRef]

Friedrich O, Norris RD, Erbacher J, authors. Evolution of middle to Late Cretaceous oceans--A 55 m.y. record of Earth’s temperature and carbon cycle. Geology. 2012;40:107–110. [CrossRef]

Huber BT, MacLeod KG, Watkins DK, Coffin MF, authors. The rise and fall of the Cretaceous hot greenhouse climate. Glob. Planet. Change. 2018;167:1–23. [CrossRef]

Huber BT, Hodell DA, Hamilton CP, authors. Middle–Late Cretaceous climate of the southern high latitudes: Stable isotopic evidence for minimal equator-to-pole thermal gradients. Geol. Soc. Am. Bull. 1995;107:1164–1191. [CrossRef]

Zachos JC, Dickens GR, Zeebe RE, authors. An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature. 2008;451:279–283. [CrossRef] [PubMed]

Anagnostou E, et al., authors. Changing atmospheric CO2 concentration was the primary driver of early Cenozoic climate. Nature. 2016;533:380–384. [CrossRef] [PubMed]

Foster GL, Royer DL, Lunt DJ, authors. Future climate forcing potentially without precedent in the last 420 million years. Nat. Commun. 2017;8:14845 [CrossRef] [PubMed]

Pucéat E, et al., authors. Fish tooth δ18O revising Late Cretaceous meridional upper ocean water temperature gradients. Geology. 2007;35:107–110. [CrossRef]

Larson RL, author. Latest pulse of Earth: evidence for a mid-Cretaceous superplume. Geology. 1991;19:547–550. [CrossRef]

Brune S, Williams SE, Müller RD, authors. Potential links between continental rifting, CO2 degassing and climate change through time. Nat. Geosci. 2017;10:941–946. [CrossRef]

Gutjahr M, et al., authors. Very large release of mostly volcanic carbon during the Palaeocene–Eocene thermal maximum. Nature. 2017;548:573–577. [CrossRef] [PubMed]

Westerhold T, Röhl U, Donner B, McCarren HK, Zachos JC, authors. A complete high-resolution Paleocene benthic stable isotope record for the central Pacific (ODP Site 1209). Paleoceanography. 2011;26:PA2216 [CrossRef]

Pucéat E, Lecuyer C, Reisberg L, authors. Neodymium isotope evolution of NW Tethyan upper ocean waters throughout the Cretaceous. Earth Planet. Sci. Lett. 2005;236:705–720. [CrossRef]

Pérez-Díaz L, Eagles G, authors. South Atlantic paleobathymetry since early Cretaceous. Sci. Rep. 2017;7:11819 [CrossRef] [PubMed]

Sewall JO, et al., authors. Climate model boundary conditions for four Cretaceous time slices. Clim. Past. 2007;3:647–657. [CrossRef]

Donnadieu Y, Pucéat E, Moiroud M, Guillocheau F, Deconinck JF, authors. A better-ventilated ocean triggered by Late Cretaceous changes in continental configuration. Nat. Commun. 2016;7:10316 [CrossRef] [PubMed]

Frank TD, Arthur MA, authors. Tectonic forcings of Maastrichtian ocean‐climate evolution. Paleoceanography. 1999;14:103–117. [CrossRef]

Friedrich O, Erbacher J, authors. Benthic foraminiferal assemblages from Demerara Rise (ODP Leg 207, western tropical Atlantic): possible evidence for a progressive opening of the Equatorial Atlantic Gateway. Cretac. Res. 2006;27:377–397. [CrossRef]

Moiroud Mathieu, Pucéat Emmanuelle, Donnadieu Yannick, Bayon Germain, Guiraud Michel, Voigt Silke, Deconinck Jean-François, Monna Fabrice, authors. Evolution of neodymium isotopic signature of seawater during the Late Cretaceous: Implications for intermediate and deep circulation. Gondwana Research. 2016;36:503–522. [CrossRef]

Robinson SA, Murphy DP, Vance D, Thomas DJ, authors. Formation of “Southern Component Water” in the Late Cretaceous: evidence from Nd-isotopes. Geology. 2010;38:871–874. [CrossRef]

Voigt S, et al., authors. Tectonically restricted deep-ocean circulation at the end of the Cretaceous greenhouse. Earth Planet. Sci. Lett. 2013;369–370:169–177. [CrossRef]

Frank, M. Radiogenic isotopes: tracers of past ocean circulation and erosional input. Rev. Geophys. 40, 10.1029/2000RG000094 (2002).

Lacan F, Jeandel C, authors. Neodymium isotopes as a new tool for quantifying exchange fluxes at the continent–ocean interface. Earth Planet. Sci. Lett. 2005;232:245–257. [CrossRef]

Rempfer J, Stocker TF, Joos F, Dutay JC, Siddall M, authors. Modelling Nd-isotopes with a coarse resolution ocean circulation model: Sensitivities to model parameters and source/sink distributions. Geochim. Cosmochim. Acta. 2011;75:5927–5950. [CrossRef]

MacLeod KG, Isaza Londoño C, Martin EE, Jiménez Berrocoso Aacute, Basak C, authors. Changes in North Atlantic circulation at the end of the Cretaceous greenhouse interval. Nat. Geosci. 2011;4:779–782. [CrossRef]

MacLeod KG, Martin EE, Blair SW, authors. Nd isotopic excursion across Cretaceous ocean anoxic event 2 (Cenomanian-Turonian) in the tropical North. Atl. Geol. 2008;36:811–814. [CrossRef]

Martin E.E.; MacLeod K.G.; Jiménez Berrocoso A.; Bourbon E., authors. Water mass circulation on Demerara Rise during the Late Cretaceous based on Nd isotopes. Earth and Planetary Science Letters. 2012;327-328:111–120. [CrossRef]

Murphy DP, Thomas DJ, authors. The evolution of Late Cretaceous deep-ocean circulation in the Atlantic basins: Neodymium isotope evidence from South Atlantic drill sites for tectonic controls. Geochem. Geophys. Geosystems. 2013;14:5323–5340. [CrossRef]

Robinson SA, Vance D, authors. Widespread and synchronous change in deep-ocean circulation in the North and South Atlantic during the Late Cretaceous. Paleoceanography. 2012;27:PA1102 [CrossRef]

Grasse P, et al., authors. Short-term variability of dissolved rare earth elements and neodymium isotopes in the entire water column of the Panama Basin. Earth Planet. Sci. Lett. 2017;475:242–253. [CrossRef]

Osborne AH, Haley BA, Hathorne EC, Flögel S, Frank M, authors. Neodymium isotopes and concentrations in Caribbean seawater: Tracing water mass mixing and continental input in a semi-enclosed ocean basin. Earth Planet. Sci. Lett. 2014;406:174–186. [CrossRef]

Abbott AN, Haley BA, Tripati AK, Frank M, authors. Constraints on ocean circulation at the Paleocene–Eocene Thermal Maximum from neodymium isotopes. Clim. Discuss. 2015;11:2557–2583. [CrossRef]

Via RK, Thomas DJ, authors. Evolution of Atlantic thermohaline circulation: Early Oligocene onset of deep-water production in the North. Atl. Geol. 2006;34:441–444. [CrossRef]

Blaser P, et al., authors. Extracting foraminiferal seawater Nd isotope signatures from bulk deep sea sediment by chemical leaching. Chem. Geol. 2016;439:189–204. [CrossRef]

Yu Z, et al., authors. Seasonal variations in dissolved neodymium isotope composition in the Bay of Bengal. Earth Planet. Sci. Lett. 2017;479:310–321. [CrossRef]

Goldstein SL, O’Nions RK, authors. Nd and Sr isotopic relationships in pelagic clays and ferromanganese deposits. Nature. 1981;292:324–327. [CrossRef]

Thomas DJ, Bralower TJ, Jones CE, authors. Neodymium isotopic reconstruction of late Paleocene–early Eocene thermohaline circulation. Earth Planet. Sci. Lett. 2003;209:309–322. [CrossRef]

Müller RD, Sdrolias M, Gaina C, Steinberger B, Heine C, authors. Long-term sea-level fluctuations driven by ocean basin dynamics. Science. 2008;319:1357–1362. [CrossRef] [PubMed]

Pérez-Díaz L, Eagles G, authors. Estimating palaeobathymetry with quantified uncertainties: a workflow illustrated with South Atlantic data. Geol. Soc. Lond. Spec. Publ. 2018;SP476:1

Littler K, Röhl U, Westerhold T, Zachos JC, authors. A high-resolution benthic stable-isotope record for the South Atlantic: Implications for orbital-scale changes in Late Paleocene–Early Eocene climate and carbon cycling. Earth Planet. Sci. Lett. 2014;401:18–30. [CrossRef]

Thomas DJ, Korty R, Huber M, Schubert JA, Haines B, authors. Nd isotopic structure of the Pacific Ocean 70-30 Ma and numerical evidence for vigorous ocean circulation and ocean heat transport in a greenhouse world. Paleoceanography. 2014;29:454–469. [CrossRef]

Sijp WP, England MH, authors. The effect of low ancient greenhouse climate temperature gradients on the ocean’s overturning circulation. Clim. Past. 2016;12:543–552. [CrossRef]

Gutjahr M, et al., authors. Reliable extraction of a deepwater trace metal isotope signal from Fe–Mn oxyhydroxide coatings of marine sediments. Chem. Geol. 2007;242:351–370. [CrossRef]

Stumpf R, Frank M, Schönfeld J, Haley BA, authors. Late quaternary variability of Mediterranean outflow water from radiogenic Nd and Pb isotopes. Quat. Sci. Rev. 2010;29:2462–2472. [CrossRef]

Barrat J, et al., authors. Determination of rare earth elements in sixteen silicate reference samples by ICP_MS after Tm addition and ion exchange separation. Geostand. Newsl. 1996;20:133–139. [CrossRef]

Le Fèvre B, Pin C, authors. A straightforward separation scheme for concomitant Lu–Hf and Sm–Nd isotope ratio and isotope dilution analysis. Anal. Chim. Acta. 2005;543:209–221. [CrossRef]

Tanaka T, et al., authors. JNdi-1: a neodymium isotopic reference in consistency with LaJolla neodymium. Chem. Geol. 2000;168:279–281. [CrossRef]

Jacobsen SB, Wasserburg GJ, authors. Sm-Nd isotopic evolution of chondrites. Earth Planet. Sci. Lett. 1980;50:139–155. [CrossRef]

Ogg J.G., author. Geomagnetic Polarity Time Scale. The Geologic Time Scale. 2012. p. 85–113

Berggren WA, Kent DV, Swisher CC, Aubry MP, authors. A Cenozoic geochronology and chronostratigraphy. Soc. Sed. Geol. Spec. Publ. 1995;54:129–212

Anthonissen D. Erik, Ogg James G., authors. Cenozoic and Cretaceous Biochronology of Planktonic Foraminifera and Calcareous Nannofossils. The Geologic Time Scale. 2012. p. 1083–1127

Berggren WA, et al., authors. Magnetobiostratigraphy of deep sea drilling project leg 72, sites 515-518, Rio Grande rise (South Atlantic). Init. Rep. Deep Sea Drill. Proj. 1983;72:938–948

Chave AD, author. Lower paleocene-upper cretaceous magnetostratigraphy, sites 525, 527, 528, and 529, deep sea drilling project leg 74. Init. Rep. Deep Sea Drill. Proj. 1984;74:525–531

Manivit H, author. Paleogene and Upper Cretaceous calcareous nannofossils from Deep Sea Drilling Project Leg 74. Init. Rep. Deep Sea Drill. Proj. 1984;74:475–499

Westerhold T, et al., authors. On the duration of magnetochrons C24r and C25n and the timing of early Eocene global warming events: Implications from the Ocean Drilling Program Leg 208 Walvis Ridge depth transect. Paleoceanography. 2007;22:PA2201 [CrossRef]

Zachos JC, et al., authors. Leg 208 summary. Proc. Ocean Drill. Prog Init. Rep. 2004;208:1–112

Lancelot, Y. Seibold, E., etal. Volume 41. Init. Rep. DeepSeaDrill. Proj. 41, (U.S. Government Printing Office, Washington, 1977).

Norris, R. D., Wilson, P. A., Blum, P. & Expedition 342 Scientists. Site 1403. In: Paleogene Newfoundland Sediment Drifts and MDHDS Test. Proc. Integrat. OceanDrill. Prog. 342, (2014).

Monechi S, Reale V, Bernaola G, Balestra B, authors. The Danian/Selandian boundary at Site 1262 (South Atlantic) and in the Tethyan region: biomagnetostratigraphy, evolutionary trends in fasciculiths and environmental effects of the Latest Danian Event. Mar. Micro. 2013;98:28–40. [CrossRef]

Dinarès-Turell J, Westerhold T, Pujalte V, Röhl U, Kroon D, authors. Astronomical calibration of the Danian stage (Early Paleocene) revisited: Settling chronologies of sedimentary records across the Atlantic and Pacific Oceans. Earth Planet. Sci. Lett. 2014;405:119–131. [CrossRef]

Batenburg SJ, et al., authors. Late Maastrichtian carbon isotope stratigraphy and cyclostratigraphy of the Newfoundland Margin (Site U1403, IODP Leg 342). Newsl. Stratigr. 2017;10:1127

Hay William W.; DeConto Robert M.; Wold Christopher N.; Wilson Kevin M.; Voigt Silke, Schulz Michael, Wold Adrienne Rossby, Dullo Wolf-Christian, Ronov Alexander B.; Balukhovsky Alexander N.; Söding Emanuel, authors. Alternative global Cretaceous paleogeography. Special Paper 332: Evolution of the Cretaceous Ocean-Climate System. 1999. p. 1–47

Gradstein, F., Ogg, J., Schmitz, M. & Ogg, G. (Eds) The Geologic Time Scale 2012 (Elsevier, Amsterdam, 2012).

[Back]

Fig. 1.

Nd-isotope data across the latest Cretaceous – early Paleogene. In this compilation of published Atlantic εNd(t) records and new Nd-isotope data from Site 369 (this study) only records with four or more data points over the time-interval 72–50 Ma are included. The εNd(t) data of Demerara Rise (Sites 1258, 1260, 1261), Maud Rise, Blake Nose and Site 525 have been smoothed by LOWESS regression (f = 0.3) with the dashed lines representing the 97.5% confidence interval. The εNd(t) data are flanked by the geomagnetic polarity time scale with chrons C31-C24. PETM: Paleocene/Eocene thermal maximum, K/Pg: Cretaceous/Paleogene Boundary

[Back]

[Back]

Fig. 2.

Paleogeographic setting. a Paleogeographic reconstruction at 60 Ma. Plate tectonic configuration after ref. ⁶³ and global bathymetry from ref. ³⁹. Circles indicate Atlantic Ocean drilling sites for which Nd-isotope data are available (at least four points in the time-interval 72–50 Ma), DR: Demerara Rise. Diamonds indicate locations for which new data are presented here. Numbers indicate DSDP/ODP/IODP Sites. b Three-dimensional plot of the reconstructed bathymetry of the South Atlantic with the Rio Grande Rise – Walvis Ridge barrier at 60 Ma from ref. ⁴⁰. c Contour plot of the bathymetry in panel B, with black numbers indicating depths in meters

[Back]

[Back]

Fig. 3.

New Nd-isotope results. Neodymium isotope values of ferromanganese coatings (coloured symbols) and detrital fractions (open symbols) measured in this study and previously for Site 525²². Error bars indicate 2 s.d. external reproducibility

[Back]

[Back]

Fig. 4.

The convergence of Nd-isotope values. Maastrichtian to early Eocene Nd-isotope records compared to atmospheric CO₂ levels and deep-sea temperature evolution. Open symbols depict new εNd(t) data for North Atlantic Site U1403 (Newfoundland margin, dark blue) and South Atlantic Sites 516 (Rio Grande Rise, red), 1267 (northern flank of the Walvis Ridge, green) and 525 (top of the Walvis Ridge, yellow). Filled symbols depict selected high-resolution εNd(t) data from Demerara Rise in the North Atlantic²⁶,²⁸ (purple), Sites 527³³ and 1262³⁴ on the flank of the Walvis Ridge (green), Site 1263³⁴ and Site 525²² on the Walvis Ridge (yellow), and Maud Rise in the Southern Ocean²⁰ (dark blue). The εNd(t) data have been smoothed by LOWESS regression (f = 0.3) with the shaded areas representing the 97.5% confidence interval. The εNd(t) data are flanked by the geomagnetic polarity time scale on the left and on the right by reconstructed atmospheric CO₂ levels¹¹, benthic oxygen isotope data from the Pacific⁴¹ and the Atlantic¹³, and the suggested evolution of the Atlantic Ocean. The Paleogene age model of the Sites rests on a combination of magnetic and nannofossil stratigraphy plotted against the GTS 2012⁶⁴. C31-C24: paleomagnetic chrons, δ¹⁸Ob: oxygen isotope ratio from benthic foraminifera, VPDB: Vienna Pee Dee Belemnite standard, PETM: Paleocene/Eocene Thermal Maximum, DR: Demerara Rise, K/Pg: Cretaceous/Paleogene boundary, MR: Maud Rise

[Back]