Global and regional temperature change over the past 4.5 million years.

Much of our understanding of Cenozoic climate is based on the record of δ18O measured in benthic foraminifera. However, this measurement reflects a combined signal of global temperature and sea level, thus preventing a clear understanding of the interactions and feedbacks of the climate system in causing global temperature change. Our new reconstruction of temperature change over the past 4.5 million years includes two phases of long-term cooling, with the second phase of accelerated cooling during the Middle Pleistocene Transition (1.5 to 0.9 million years ago) being accompanied by a transition from dominant 41,000-year low-amplitude periodicity to dominant 100,000-year high-amplitude periodicity. Changes in the rates of long-term cooling and variability are consistent with changes in the carbon cycle driven initially by geologic processes, followed by additional changes in the Southern Ocean carbon cycle.

A well-oxygenated eastern tropical Pacific during the warm Miocene.

The oxygen content of the oceans is susceptible to climate change and has declined in recent decades1, with the largest effect in oxygen-deficient zones (ODZs)2, that is, mid-depth ocean regions with oxygen concentrations <5 µmol kg-1 (ref. 3). Earth-system-model simulations of climate warming predict that ODZs will expand until at least 2100. The response on timescales of hundreds to thousands of years, however, remains uncertain3-5. Here we investigate changes in the response of ocean oxygenation during the warmer-than-present Miocene Climatic Optimum (MCO; 17.0-14.8 million years ago (Ma)). Our planktic foraminifera I/Ca and δ15N data, palaeoceanographic proxies sensitive to ODZ extent and intensity, indicate that dissolved-oxygen concentrations in the eastern tropical Pacific (ETP) exceeded 100 µmol kg-1 during the MCO. Paired Mg/Ca-derived temperature data suggest that an ODZ developed in response to an increased west-to-east temperature gradient and shoaling of the ETP thermocline. Our records align with model simulations of data from recent decades to centuries6,7, suggesting that weaker equatorial Pacific trade winds during warm periods may lead to decreased upwelling in the ETP, causing equatorial productivity and subsurface oxygen demand to be less concentrated in the east. These findings shed light on how warm-climate states such as during the MCO may affect ocean oxygenation. If the MCO is considered as a possible analogue for future warming, our findings seem to support models suggesting that the recent deoxygenation trend and expansion of the ETP ODZ may eventually reverse3,4.

Accelerated mafic weathering in Southeast Asia linked to late Neogene cooling.

Arc-continent collision in Southeast Asia during the Neogene may have driven global cooling through chemical weathering of freshly exposed ophiolites resulting in atmospheric CO2 removal. Yet, little is known about the cause-and-effect relationships between erosion and the long-term evolution of tectonics and climate in this region. Here, we present an 8-million-year record of seawater chemistry and sediment provenance from the eastern Indian Ocean, near the outflow of Indonesian Throughflow waters. Using geochemical analyses of foraminiferal shells and grain size-specific detrital fractions, we show that erosion and chemical weathering of ophiolitic rocks markedly increased after 4 million years (Ma), coincident with widespread island emergence and gradual strengthening of Pacific zonal sea-surface temperature gradients. Together with supportive evidence for enhanced mafic weathering at that time from re-analysis of the seawater 87Sr/86Sr curve, this finding suggests that island uplift and hydroclimate change in the western Pacific contributed to maintaining high atmospheric CO2 consumption throughout the late Neogene.

North Atlantic cooling triggered a zonal mode over the Indian Ocean during Heinrich Stadial 1.

Abrupt changes in the Atlantic meridional overturning circulation (AMOC) are thought to affect tropical hydroclimate through adjustment of the latitudinal position of the intertropical convergence zone (ITCZ). Heinrich Stadial 1 (HS1) involves the largest AMOC reduction in recent geological time; however, over the tropical Indian Ocean (IO), proxy records suggest zonal anomalies featuring intense, widespread drought in tropical East Africa versus generally wet but heterogeneous conditions in the Maritime Continent. Here, we synthesize proxy data and an isotope-enabled transient deglacial simulation and show that the southward ITCZ shift over the eastern IO during HS1 strengthens IO Walker circulation, triggering an east-west precipitation dipole across the basin. This dipole reverses the zonal precipitation anomalies caused by the exposed Sunda and Sahul shelves due to glacial lower sea level. Our study illustrates how zonal modes of atmosphere-ocean circulation can amplify or reverse global climate anomalies, highlighting their importance for future climate change.

Warm pool ocean heat content regulates ocean-continent moisture transport.

The Indo-Pacific Warm Pool (IPWP) exerts a dominant role in global climate by releasing huge amounts of water vapour and latent heat to the atmosphere and modulating upper ocean heat content (OHC), which has been implicated in modern climate change1. The long-term variations of IPWP OHC and their effect on monsoonal hydroclimate are, however, not fully explored. Here, by combining geochemical proxies and transient climate simulations, we show that changes of IPWP upper (0-200 m) OHC over the past 360,000 years exhibit dominant precession and weaker obliquity cycles and follow changes in meridional insolation gradients, and that only 30%-40% of the deglacial increases are related to changes in ice volume. On the precessional band, higher upper OHC correlates with oxygen isotope enrichments in IPWP surface water and concomitant depletion in East Asian precipitation as recorded in Chinese speleothems. Using an isotope-enabled air-sea coupled model, we suggest that on precessional timescales, variations in IPWP upper OHC, more than surface temperature, act to amplify the ocean-continent hydrological cycle via the convergence of moisture and latent heat. From an energetic viewpoint, the coupling of upper OHC and monsoon variations, both coordinated by insolation changes on orbital timescales, is critical for regulating the global hydroclimate.

Cyclic evolution of phytoplankton forced by changes in tropical seasonality.

Although the role of Earth's orbital variations in driving global climate cycles has long been recognized, their effect on evolution is hitherto unknown. The fossil remains of coccolithophores, a key calcifying phytoplankton group, enable a detailed assessment of the effect of cyclic orbital-scale climate changes on evolution because of their abundance in marine sediments and the preservation of their morphological adaptation to the changing environment1,2. Evolutionary genetic analyses have linked broad changes in Pleistocene fossil coccolith morphology to species radiation events3. Here, using high-resolution coccolith data, we show that during the last 2.8 million years the morphological evolution of coccolithophores was forced by Earth's orbital eccentricity with rhythms of around 100,000 years and 405,000 years-a distinct spectral signature to that of coeval global climate cycles4. Simulations with an Earth System Model5 coupled with an ocean biogeochemical model6 show a strong eccentricity modulation of the seasonal cycle, which we suggest directly affects the diversity of ecological niches that occur over the annual cycle in the tropical ocean. Reduced seasonality in surface ocean conditions favours species with mid-size coccoliths, increasing coccolith carbonate export and burial; whereas enhanced seasonality favours a larger range of coccolith sizes and reduced carbonate export. We posit that eccentricity pacing of phytoplankton evolution contributed to the strong 405,000-year cyclicity that is seen in global carbon cycle records.

Remote and local drivers of Pleistocene South Asian summer monsoon precipitation: A test for future predictions.

South Asian precipitation amount and extreme variability are predicted to increase due to thermodynamic effects of increased 21st-century greenhouse gases, accompanied by an increased supply of moisture from the southern hemisphere Indian Ocean. We reconstructed South Asian summer monsoon precipitation and runoff into the Bay of Bengal to assess the extent to which these factors also operated in the Pleistocene, a time of large-scale natural changes in carbon dioxide and ice volume. South Asian precipitation and runoff are strongly coherent with, and lag, atmospheric carbon dioxide changes at Earth's orbital eccentricity, obliquity, and precession bands and are closely tied to cross-equatorial wind strength at the precession band. We find that the projected monsoon response to ongoing, rapid high-latitude ice melt and rising carbon dioxide levels is fully consistent with dynamics of the past 0.9 million years.

Seasonal origin of the thermal maxima at the Holocene and the last interglacial.

Proxy reconstructions from marine sediment cores indicate peak temperatures in the first half of the last and current interglacial periods (the thermal maxima of the Holocene epoch, 10,000 to 6,000 years ago, and the last interglacial period, 128,000 to 123,000 years ago) that arguably exceed modern warmth1-3. By contrast, climate models simulate monotonic warming throughout both periods4-7. This substantial model-data discrepancy undermines confidence in both proxy reconstructions and climate models, and inhibits a mechanistic understanding of recent climate change. Here we show that previous global reconstructions of temperature in the Holocene1-3 and the last interglacial period8 reflect the evolution of seasonal, rather than annual, temperatures and we develop a method of transforming them to mean annual temperatures. We further demonstrate that global mean annual sea surface temperatures have been steadily increasing since the start of the Holocene (about 12,000 years ago), first in response to retreating ice sheets (12 to 6.5 thousand years ago), and then as a result of rising greenhouse gas concentrations (0.25 ± 0.21 degrees Celsius over the past 6,500 years or so). However, mean annual temperatures during the last interglacial period were stable and warmer than estimates of temperatures during the Holocene, and we attribute this to the near-constant greenhouse gas levels and the reduced extent of ice sheets. We therefore argue that the climate of the Holocene differed from that of the last interglacial period in two ways: first, larger remnant glacial ice sheets acted to cool the early Holocene, and second, rising greenhouse gas levels in the late Holocene warmed the planet. Furthermore, our reconstructions demonstrate that the modern global temperature has exceeded annual levels over the past 12,000 years and probably approaches the warmth of the last interglacial period (128,000 to 115,000 years ago).

Pacific warm pool subsurface heat sequestration modulated Walker circulation and ENSO activity during the Holocene.

Dynamics driving the El Niño-Southern Oscillation (ENSO) over longer-than-interannual time scales are poorly understood. Here, we compile thermocline temperature records of the Indo-Pacific warm pool over the past 25,000 years, which reveal a major warming in the Early Holocene and a secondary warming in the Middle Holocene. We suggest that the first thermocline warming corresponds to heat transport of southern Pacific shallow overturning circulation driven by June (austral winter) insolation maximum. The second thermocline warming follows equatorial September insolation maximum, which may have caused a steeper west-east upper-ocean thermal gradient and an intensified Walker circulation in the equatorial Pacific. We propose that the warm pool thermocline warming ultimately reduced the interannual ENSO activity in the Early to Middle Holocene. Thus, a substantially increased oceanic heat content of the warm pool, acting as a negative feedback for ENSO in the past, may play its role in the ongoing global warming.

Interglacial instability of North Atlantic Deep Water ventilation.

Disrupting North Atlantic Deep Water (NADW) ventilation is a key concern in climate projections. We use (sub)centennially resolved bottom water δ13C records that span the interglacials of the last 0.5 million years to assess the frequency of and the climatic backgrounds capable of triggering large NADW reductions. Episodes of reduced NADW in the deep Atlantic, similar in magnitude to glacial events, have been relatively common and occasionally long-lasting features of interglacials. NADW reductions were triggered across the range of recent interglacial climate backgrounds, which demonstrates that catastrophic freshwater outburst floods were not a prerequisite for large perturbations. Our results argue that large NADW disruptions are more easily achieved than previously appreciated and that they occurred in past climate conditions similar to those we may soon face.

A low climate threshold for south Greenland Ice Sheet demise during the Late Pleistocene.

The Greenland Ice Sheet (GIS) has been losing mass at an accelerating rate over the recent decades. Models suggest a possible temperature threshold between 0.8 and 3.2 °C, beyond which GIS decline becomes irreversible. The duration of warmth above a given threshold is also a critical determinant for GIS survival, underlining the role of ocean warming, as its inertia prolongs warmth and triggers longer-term feedbacks. The exact point at which these feedbacks are triggered remains equivocal. Late Pleistocene interglacials provide potential case examples for constraining the past response of the GIS to a range of climate states, including conditions warmer than present. However, little is known about the magnitude and duration of warming near Greenland during these periods. Using high-resolution multiproxy surface ocean climate records off southern Greenland, we show that the previous 4 interglacials over the last ∼450 ka all reached warmer than present climate conditions and exceeded the modeled temperature threshold for GIS collapse but by different magnitudes and durations. Complete deglaciation of the southern GIS in Marine Isotope Stage 11c (MIS 11c; 394.7 to 424.2 ka) occurred under climates only slightly warmer than present (∼0.5 ± 1.6 °C), placing the temperature threshold for major GIS retreat in the lower end of model estimates and within projections for this century.

Asynchronous Antarctic and Greenland ice-volume contributions to the last interglacial sea-level highstand.

The last interglacial (LIG; ~130 to ~118 thousand years ago, ka) was the last time global sea level rose well above the present level. Greenland Ice Sheet (GrIS) contributions were insufficient to explain the highstand, so that substantial Antarctic Ice Sheet (AIS) reduction is implied. However, the nature and drivers of GrIS and AIS reductions remain enigmatic, even though they may be critical for understanding future sea-level rise. Here we complement existing records with new data, and reveal that the LIG contained an AIS-derived highstand from ~129.5 to ~125 ka, a lowstand centred on 125-124 ka, and joint AIS + GrIS contributions from ~123.5 to ~118 ka. Moreover, a dual substructure within the first highstand suggests temporal variability in the AIS contributions. Implied rates of sea-level rise are high (up to several meters per century; m c-1), and lend credibility to high rates inferred by ice modelling under certain ice-shelf instability parameterisations.

Capturing the global signature of surface ocean acidification during the Palaeocene-Eocene Thermal Maximum.

Geologically abrupt carbon perturbations such as the Palaeocene-Eocene Thermal Maximum (PETM, approx. 56 Ma) are the closest geological points of comparison to current anthropogenic carbon emissions. Associated with the rapid carbon release during this event are profound environmental changes in the oceans including warming, deoxygenation and acidification. To evaluate the global extent of surface ocean acidification during the PETM, we present a compilation of new and published surface ocean carbonate chemistry and pH reconstructions from various palaeoceanographic settings. We use boron to calcium ratios (B/Ca) and boron isotopes (δ11B) in surface- and thermocline-dwelling planktonic foraminifera to reconstruct ocean carbonate chemistry and pH. Our records exhibit a B/Ca reduction of 30-40% and a δ11B decline of 1.0-1.2‰ coeval with the carbon isotope excursion. The tight coupling between boron proxies and carbon isotope records is consistent with the interpretation that oceanic absorption of the carbon released at the onset of the PETM resulted in widespread surface ocean acidification. The remarkable similarity among records from different ocean regions suggests that the degree of ocean carbonate change was globally near uniform. We attribute the global extent of surface ocean acidification to elevated atmospheric carbon dioxide levels during the main phase of the PETM.This article is part of a discussion meeting issue 'Hyperthermals: rapid and extreme global warming in our geological past'.

Paleoceanography. Antarctic role in Northern Hemisphere glaciation.

Earth's climate underwent a major transition from the warmth of the late Pliocene, when global surface temperatures were ~2° to 3°C higher than today, to extensive Northern Hemisphere glaciation (NHG) ~2.73 million years ago (Ma). We show that North Pacific deep waters were substantially colder (4°C) and probably fresher than the North Atlantic Deep Water before the intensification of NHG. At ~2.73 Ma, the Atlantic-Pacific temperature gradient was reduced to <1°C, suggesting the initiation of stronger heat transfer from the North Atlantic to the deep Pacific. We posit that increased glaciation of Antarctica, deduced from the 21 ± 10-meter sea-level fall from 3.15 to 2.75 Ma, and the development of a strong polar halocline fundamentally altered deep ocean circulation, which enhanced interhemispheric heat and salt transport, thereby contributing to NHG.

The Holocene temperature conundrum.

A recent temperature reconstruction of global annual temperature shows Early Holocene warmth followed by a cooling trend through the Middle to Late Holocene [Marcott SA, et al., 2013, Science 339(6124):1198-1201]. This global cooling is puzzling because it is opposite from the expected and simulated global warming trend due to the retreating ice sheets and rising atmospheric greenhouse gases. Our critical reexamination of this contradiction between the reconstructed cooling and the simulated warming points to potentially significant biases in both the seasonality of the proxy reconstruction and the climate sensitivity of current climate models.

Rapid reductions in North Atlantic Deep Water during the peak of the last interglacial period.

Deep ocean circulation has been considered relatively stable during interglacial periods, yet little is known about its behavior on submillennial time scales. Using a subcentennially resolved epibenthic foraminiferal δ(13)C record, we show that the influence of North Atlantic Deep Water (NADW) was strong at the onset of the last interglacial period and was then interrupted by several prominent centennial-scale reductions. These NADW transients occurred during periods of increased ice rafting and southward expansions of polar water influence, suggesting that a buoyancy threshold for convective instability was triggered by freshwater and circum-Arctic cryosphere changes. The deep Atlantic chemical changes were similar in magnitude to those associated with glaciations, implying that the canonical view of a relatively stable interglacial circulation may not hold for conditions warmer and fresher than at present.

Pacific Ocean heat content during the past 10,000 years.

Observed increases in ocean heat content (OHC) and temperature are robust indicators of global warming during the past several decades. We used high-resolution proxy records from sediment cores to extend these observations in the Pacific 10,000 years beyond the instrumental record. We show that water masses linked to North Pacific and Antarctic intermediate waters were warmer by 2.1 ± 0.4°C and 1.5 ± 0.4°C, respectively, during the middle Holocene Thermal Maximum than over the past century. Both water masses were ~0.9°C warmer during the Medieval Warm period than during the Little Ice Age and ~0.65° warmer than in recent decades. Although documented changes in global surface temperatures during the Holocene and Common era are relatively small, the concomitant changes in OHC are large.

Aragonite precipitation by "proto-polyps" in coral cell cultures.

The mechanisms of coral calcification at the molecular, cellular and tissue levels are poorly understood. In this study, we examine calcium carbonate precipitation using novel coral tissue cultures that aggregate to form "proto-polyps". Our goal is to establish an experimental system in which calcification is facilitated at the cellular level, while simultaneously allowing in vitro manipulations of the calcifying fluid. This novel coral culturing technique enables us to study the mechanisms of biomineralization and their implications for geochemical proxies. Viable cell cultures of the hermatypic, zooxanthellate coral, Stylophora pistillata, have been maintained for 6 to 8 weeks. Using an enriched seawater medium with aragonite saturation state similar to open ocean surface waters (Ω(arag)~4), the primary cell cultures assemble into "proto-polyps" which form an extracellular organic matrix (ECM) and precipitate aragonite crystals. These extracellular aragonite crystals, about 10 µm in length, are formed on the external face of the proto-polyps and are identified by their distinctive elongated crystallography and X-ray diffraction pattern. The precipitation of aragonite is independent of photosynthesis by the zooxanthellae, and does not occur in control experiments lacking coral cells or when the coral cells are poisoned with sodium azide. Our results demonstrate that proto-polyps, aggregated from primary coral tissue culture, function (from a biomineralization perspective) similarly to whole corals. This approach provides a novel tool for investigating the biophysical mechanism of calcification in these organisms.