Mean global ocean temperatures during the last glacial transition

Deeper and stronger North Atlantic Gyre during the Last Glacial Maximum

Subtropical gyre (STG) depth and strength are controlled by wind stress curl and surface buoyancy forcing1,2. Modern hydrographic data reveal that the STG extends to a depth of about 1 km in the Northwest Atlantic, with its maximum depth defined by the base of the subtropical thermocline. Despite the likelihood of greater wind stress curl and surface buoyancy loss during the Last Glacial Maximum (LGM)3, previous work suggests minimal change in the depth of the glacial STG4. Here we show a sharp glacial water mass boundary between 33° N and 36° N extending down to between 2.0 and 2.5 km-approximately 1 km deeper than today. Our findings arise from benthic foraminiferal δ18O profiles from sediment cores in two depth transects at Cape Hatteras (36-39° N) and Blake Outer Ridge (29-34° N) in the Northwest Atlantic. This result suggests that the STG, including the Gulf Stream, was deeper and stronger during the LGM than at present, which we attribute to increased glacial wind stress curl, as supported by climate model simulations, as well as greater glacial production of denser subtropical mode waters (STMWs). Our data suggest (1) that subtropical waters probably contributed to the geochemical signature of what is conventionally identified as Glacial North Atlantic Intermediate Water (GNAIW)5-7 and (2) the STG helped sustain continued buoyancy loss, water mass conversion and northwards meridional heat transport (MHT) in the glacial North Atlantic.

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.

Global oceanic oxygenation controlled by the Southern Ocean through the last deglaciation

Ocean dissolved oxygen (DO) can provide insights on how the marine carbon cycle affects global climate change. However, the net global DO change and the controlling mechanisms remain uncertain through the last deglaciation. Here, we present a globally integrated DO reconstruction using thallium isotopes, corroborating lower global DO during the Last Glacial Maximum [19 to 23 thousand years before the present (ka B.P.)] relative to the Holocene. During the deglaciation, we reveal reoxygenation in the Heinrich Stadial 1 (~14.7 to 18 ka B.P.) and the Younger Dryas (11.7 to 12.9 ka B.P.), with deoxygenation during the Bølling-Allerød (12.9 to 14.7 ka B.P.). The deglacial DO changes were decoupled from North Atlantic Deep Water formation rates and imply that Southern Ocean ventilation controlled ocean oxygen. The coherence between global DO and atmospheric CO2 on millennial timescales highlights the Southern Ocean's role in deglacial atmospheric CO2 rise.

Enhanced subglacial discharge from Antarctica during meltwater pulse 1A

Subglacial discharge from the Antarctic Ice Sheet (AIS) likely played a crucial role in the loss of the ice sheet and the subsequent rise in sea level during the last deglaciation. However, no direct proxy is currently available to document subglacial discharge from the AIS, which leaves significant gaps in our understanding of the complex interactions between subglacial discharge and ice-sheet stability. Here we present deep-sea coral 234U/238U records from the Drake Passage in the Southern Ocean to track subglacial discharge from the AIS. Our findings reveal distinctively higher seawater 234U/238U values from 15,400 to 14,000 years ago, corresponding to the period of the highest iceberg-rafted debris flux and the occurrence of the meltwater pulse 1A event. This correlation suggests a causal link between enhanced subglacial discharge, synchronous retreat of the AIS, and the rapid rise in sea levels. The enhanced subglacial discharge and subsequent AIS retreat appear to have been preconditioned by a stronger and warmer Circumpolar Deep Water, thus underscoring the critical role of oceanic heat in driving major ice-sheet retreat.

Tropical mountain ice core δ18O: A Goldilocks indicator for global temperature change

Very high tropical alpine ice cores provide a distinct paleoclimate record for climate changes in the middle and upper troposphere. However, the climatic interpretation of a key proxy, the stable water oxygen isotopic ratio in ice cores (δ18Oice), remains an outstanding problem. Here, combining proxy records with climate models, modern satellite measurements, and radiative-convective equilibrium theory, we show that the tropical δ18Oice is an indicator of the temperature of the middle and upper troposphere, with a glacial cooling of -7.35° ± 1.1°C (66% CI). Moreover, it severs as a "Goldilocks-type" indicator of global mean surface temperature change, providing the first estimate of glacial stage cooling that is independent of marine proxies as -5.9° ± 1.2°C. Combined with all estimations available gives the maximum likelihood estimate of glacial cooling as -5.85° ± 0.51°C.

Widespread global disparities between modelled and observed mid-depth ocean currents

The mid-depth ocean circulation is critically linked to actual changes in the long-term global climate system. However, in the past few decades, predictions based on ocean circulation models highlight the lack of data, knowledge, and long-term implications in climate change assessment. Here, using 842,421 observations produced by Argo floats from 2001-2020, and Lagrangian simulations, we show that only 3.8% of the mid-depth oceans, including part of the equatorial Pacific Ocean and the Antarctic Circumpolar Current, can be regarded as accurately modelled, while other regions exhibit significant underestimations in mean current velocity. Knowledge of ocean circulation is generally more complete in the low-latitude oceans but is especially poor in high latitude regions. Accordingly, we propose improvements in forecasting, model representation of stochasticity, and enhancement of observations of ocean currents. The study demonstrates that knowledge and model representations of global circulation are substantially compromised by inaccuracies of significant magnitude and direction, with important implications for modelled predictions of currents, temperature, carbon dioxide sequestration, and sea-level rise trends.

Multi-proxy constraints on Atlantic circulation dynamics since the last ice age

Uncertainties persist in the understanding of the Atlantic meridional overturning circulation and its response to external perturbations such as freshwater or radiative forcing. Abrupt reduction of the Atlantic circulation is considered a climate tipping point that may have been crossed when Earth's climate was propelled out of the last ice age. However, the evolution of the circulation since the Last Glacial Maximum (22-18 thousand years ago) remains insufficiently constrained due to model and proxy limitations. Here we leverage information from both a compilation of proxy records that track various aspects of the circulation and climate model simulations to constrain the Atlantic circulation over the past 20,000 years. We find a coherent picture of a shallow and weak Atlantic overturning circulation during the Last Glacial Maximum that reconciles apparently conflicting proxy evidence. Model-data comparison of the last deglaciation-starting from this new, multiple constrained glacial state-indicates a muted response during Heinrich Stadial 1 and that water mass geometry did not fully adjust to the strong reduction in overturning circulation during the comparably short Younger Dryas period. This demonstrates that the relationship between freshwater forcing and Atlantic overturning strength is strongly dependent on the climatic and oceanic background state.

Revisiting the Holocene global temperature conundrum

Recent global temperature reconstructions for the current interglacial period (the Holocene, beginning 11,700 years ago) have generated contrasting trends. This Review examines evidence from indicators and drivers of global change, as inferred from proxy records and simulated by climate models, to evaluate whether anthropogenic global warming was preceded by a long-term warming trend or by global cooling. Multimillennial-scale cooling before industrialization requires extra climate forcing and major climate feedbacks that are not well represented in most climate models at present. Conversely, global warming before industrialization challenges proxy-based reconstructions of past climate. The resolution of this conundrum has implications for contextualizing post-industrial warming and for understanding climate sensitivity to several forcings and their attendant feedbacks, including greenhouse gases. From a large variety of available evidence, we find support for a relatively mild millennial-scale global thermal maximum during the mid-Holocene, but more research is needed to firmly resolve the conundrum and to advance our understanding of slow-moving climate variability.

Evidence for late-glacial oceanic carbon redistribution and discharge from the Pacific Southern Ocean

Southern Ocean deep-water circulation plays a vital role in the global carbon cycle. On geological time scales, upwelling along the Chilean margin likely contributed to the deglacial atmospheric carbon dioxide rise, but little quantitative evidence exists of carbon storage. Here, we develop an X-ray Micro-Computer-Tomography method to assess foraminiferal test dissolution as proxy for paleo-carbonate ion concentrations ([CO₃^2-]). Our subantarctic Southeast Pacific sediment core depth transect shows significant deep-water [CO₃^2-] variations during the Last Glacial Maximum and Deglaciation (10-22 ka BP). We provide evidence for an increase in [CO₃^2-] during the early-deglacial period (15-19 ka BP) in Lower Circumpolar Deepwater. The export of such low-carbon deep-water from the Pacific to the Atlantic contributed to significantly lowered carbon storage within the Southern Ocean, highlighting the importance of a dynamic Pacific-Southern Ocean deep-water reconfiguration for shaping late-glacial oceanic carbon storage, and subsequent deglacial oceanic-atmospheric CO₂ transfer.

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 change¹. 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.

Thermal coupling of the Indo-Pacific warm pool and Southern Ocean over the past 30,000 years

The role of the tropical Pacific Ocean and its linkages to the southern hemisphere during the last deglacial warming remain highly controversial. Here we explore the evolution of Pacific horizontal and vertical thermal gradients over the past 30 kyr by compiling 340 sea surface and 7 subsurface temperature records, as well as one new ocean heat content record. Our records reveal that La Niña-like conditions dominated during the deglaciation as a result of the more intense warming in the western Pacific warm pool. Both the subsurface temperature and ocean heat content in the warm pool rose earlier than the sea surface temperature, and in phase with South Pacific subsurface temperature and orbital precession, implying that heat exchange between the tropical upper water column and the extratropical Southern Ocean facilitated faster warming in the western Pacific. Our study underscores the key role of the thermal coupling between the warm pool and the Southern Ocean and its relevance for future global warming.

The mechanism of the last deglacial global warming is key for future climate. Here, the authors shed light on the pivotal role of the thermal coupling between the western Pacific warm pool and the Southern Ocean.

Multiple carbon cycle mechanisms associated with the glaciation of Marine Isotope Stage 4

Here we use high-precision carbon isotope data (δ¹³C-CO₂) to show atmospheric CO₂ during Marine Isotope Stage 4 (MIS 4, ~70.5-59 ka) was controlled by a succession of millennial-scale processes. Enriched δ¹³C-CO₂ during peak glaciation suggests increased ocean carbon storage. Variations in δ¹³C-CO₂ in early MIS 4 suggest multiple processes were active during CO₂ drawdown, potentially including decreased land carbon and decreased Southern Ocean air-sea gas exchange superposed on increased ocean carbon storage. CO₂ remained low during MIS 4 while δ¹³C-CO₂ fluctuations suggest changes in Southern Ocean and North Atlantic air-sea gas exchange. A 7 ppm increase in CO₂ at the onset of Dansgaard-Oeschger event 19 (72.1 ka) and 27 ppm increase in CO₂ during late MIS 4 (Heinrich Stadial 6, ~63.5-60 ka) involved additions of isotopically light carbon to the atmosphere. The terrestrial biosphere and Southern Ocean air-sea gas exchange are possible sources, with the latter event also involving decreased ocean carbon storage.

Summary for general audience: We used carbon stable isotope data from an Antarctic ice core to evaluate which mechanisms caused changes in atmospheric carbon dioxide 74-59 thousand years ago, including a ~40 ppm decrease at the beginning of the last ice age.

Evolution of tropical land temperature across the last glacial termination

The tropical West Pacific hosts the warmest part of the surface ocean and has a considerable impact on the global climate system. Reconstructions of past temperature in this region can elucidate climate connections between the tropics and poles and the sensitivity of tropical temperature to greenhouse forcing. However, existing data are equivocal and reliable information from terrestrial archives is particularly sparse. Here we constrain the magnitude and timing of land temperature change in the tropical West Pacific across the last deglaciation using an exceptionally precise paleothermometer applied to a well-dated stalagmite from Northern Borneo. We show that the cave temperature increased by 4.4 ± 0.3 °C (2 SEM) from the Last Glacial Maximum to the Holocene, amounting to 3.6 ± 0.3 °C (2 SEM) when correcting for sea-level induced cave altitude change. The warming closely follows atmospheric CO₂ and Southern Hemisphere warming. This contrasts with hydroclimate, as reflected by drip water δ¹⁸O, which responds to Northern Hemisphere cooling events in the form of prominent drying, while temperature was rising. Our results thus show a close response of tropical temperature to greenhouse forcing, independent of shifts in the tropical circulation patterns.

Ancient drip water in a Borneo stalagmite reveals a strong land temperature rise across the last glacial termination in close correspondence with atmospheric CO2, and an intriguing decoupling between tropical temperature and hydroclimate.

Abrupt climate changes in the last two deglaciations simulated with different Northern ice sheet discharge and insolation

There were significant differences between the last two deglaciations, particularly in Atlantic Meridional Overturning Circulation (AMOC) and Antarctic warming in the deglaciations and the following interglacials. Here, we present transient simulations of deglaciation using a coupled atmosphere–ocean general circulation model for the last two deglaciations focusing on the impact of ice sheet discharge on climate changes associated with the AMOC in the first part, and the sensitivity studies using a Northern Hemisphere ice sheet model in the second part. We show that a set of abrupt climate changes of the last deglaciation, including Bolling–Allerod warming, the Younger Dryas, and onset of the Holocene were simulated with gradual changes of both ice sheet discharge and radiative forcing. On the other hand, penultimate deglaciation, with the abrupt climate change only at the beginning of the last interglacial was simulated when the ice sheet discharge was greater than in the last deglaciation by a factor of 1.5. The results, together with Northern Hemisphere ice sheet model experiments suggest the importance of the transient climate and AMOC responses to the different orbital forcing conditions of the last two deglaciations, through the mechanisms of mass loss of the Northern Hemisphere ice sheet and meltwater influx to the ocean.

Globally resolved surface temperatures since the Last Glacial Maximum

Climate changes across the past 24,000 years provide key insights into Earth system responses to external forcing. Climate model simulations^1,2 and proxy data^3-8 have independently allowed for study of this crucial interval; however, they have at times yielded disparate conclusions. Here, we leverage both types of information using paleoclimate data assimilation^9,10 to produce the first proxy-constrained, full-field reanalysis of surface temperature change spanning the Last Glacial Maximum to present at 200-year resolution. We demonstrate that temperature variability across the past 24 thousand years was linked to two primary climatic mechanisms: radiative forcing from ice sheets and greenhouse gases; and a superposition of changes in the ocean overturning circulation and seasonal insolation. In contrast with previous proxy-based reconstructions^6,7 our results show that global mean temperature has slightly but steadily warmed, by ~0.5 °C, since the early Holocene (around 9 thousand years ago). When compared with recent temperature changes¹¹, our reanalysis indicates that both the rate and magnitude of modern warming are unusual relative to the changes of the past 24 thousand years.

Sea level and deep-sea temperature reconstructions suggest quasi-stable states and critical transitions over the past 40 million years

Sea level and deep-sea temperature variations are key indicators of global climate changes. For continuous records over millions of years, deep-sea carbonate microfossil-based δ¹⁸O (δ_c) records are indispensable because they reflect changes in both deep-sea temperature and seawater δ¹⁸O (δ_w); the latter are related to ice volume and, thus, to sea level changes. Deep-sea temperature is usually resolved using elemental ratios in the same benthic microfossil shells used for δ_c, with linear scaling of residual δ_w to sea level changes. Uncertainties are large and the linear-scaling assumption remains untested. Here, we present a new process-based approach to assess relationships between changes in sea level, mean ice sheet δ¹⁸O, and both deep-sea δ_w and temperature and find distinct nonlinearity between sea level and δ_w changes. Application to δ_c records over the past 40 million years suggests that Earth's climate system has complex dynamical behavior, with threshold-like adjustments (critical transitions) that separate quasi-stable deep-sea temperature and ice-volume states.

Widespread six degrees Celsius cooling on land during the Last Glacial Maximum

The magnitude of global cooling during the Last Glacial Maximum (LGM, the coldest multimillennial interval of the last glacial period) is an important constraint for evaluating estimates of Earth's climate sensitivity^1,2. Reliable LGM temperatures come from high-latitude ice cores^3,4, but substantial disagreement exists between proxy records in the low latitudes^1,5-8, where quantitative low-elevation records on land are scarce. Filling this data gap, noble gases in ancient groundwater record past land surface temperatures through a direct physical relationship that is rooted in their temperature-dependent solubility in water^9,10. Dissolved noble gases are suitable tracers of LGM temperature because of their complete insensitivity to biological and chemical processes and the ubiquity of LGM-aged groundwater around the globe^11,12. However, although several individual noble gas studies have found substantial tropical LGM cooling^13-16, they have used different methodologies and provide limited spatial coverage. Here we use noble gases in groundwater to show that the low-altitude, low-to-mid-latitude land surface (45 degrees south to 35 degrees north) cooled by 5.8 ± 0.6 degrees Celsius (mean ± 95% confidence interval) during the LGM. Our analysis includes four decades of groundwater noble gas data from six continents, along with new records from the tropics, all of which were interpreted using the same physical framework. Our land-based result broadly supports a recent reconstruction based on marine proxy data assimilation¹ that suggested greater climate sensitivity than previous estimates^5-7.

Combining Modern and Paleoceanographic Perspectives on Ocean Heat Uptake

Monitoring Earth's energy imbalance requires monitoring changes in the heat content of the ocean. Recent observational estimates indicate that ocean heat uptake is accelerating in the twenty-first century. Examination of estimates of ocean heat uptake over the industrial era, the Common Era of the last 2,000 years, and the period since the Last Glacial Maximum, 20,000 years ago, permits a wide perspective on modern-day warming rates. In addition, this longer-term focus illustrates how the dynamics of the deep ocean and the cryosphere were active in the past and are still active today. The large climatic shifts that started with the melting of the great ice sheets have involved significant ocean heat uptake that was sustained over centuries and millennia, and modern-ocean heat content changes are small by comparison.

Timing and structure of the Younger Dryas event and its underlying climate dynamics

The Younger Dryas (YD), arguably the most widely studied millennial-scale extreme climate event, was characterized by diverse hydroclimate shifts globally and severe cooling at high northern latitudes that abruptly punctuated the warming trend from the last glacial to the present interglacial. To date, a precise understanding of its trigger, propagation, and termination remains elusive. Here, we present speleothem oxygen-isotope data that, in concert with other proxy records, allow us to quantify the timing of the YD onset and termination at an unprecedented subcentennial temporal precision across the North Atlantic, Asian Monsoon-Westerlies, and South American Monsoon regions. Our analysis suggests that the onsets of YD in the North Atlantic (12,870 ± 30 B.P.) and the Asian Monsoon-Westerlies region are essentially synchronous within a few decades and lead the onset in Antarctica, implying a north-to-south climate signal propagation via both atmospheric (decadal-time scale) and oceanic (centennial-time scale) processes, similar to the Dansgaard-Oeschger events during the last glacial period. In contrast, the YD termination may have started first in Antarctica at ∼11,900 B.P., or perhaps even earlier in the western tropical Pacific, followed by the North Atlantic between ∼11,700 ± 40 and 11,610 ± 40 B.P. These observations suggest that the initial YD termination might have originated in the Southern Hemisphere and/or the tropical Pacific, indicating a Southern Hemisphere/tropics to North Atlantic-Asian Monsoon-Westerlies directionality of climatic recovery.

The Biological Pump During the Last Glacial Maximum

Much of the global cooling during ice ages arose from changes in ocean carbon storage that lowered atmospheric CO₂. A slew of mechanisms, both physical and biological, have been proposed as key drivers of these changes. Here we discuss the current understanding of these mechanisms with a focus on how they altered the theoretically defined soft-tissue and biological disequilibrium carbon storage at the peak of the last ice age. Observations and models indicate a role for Antarctic sea ice through its influence on ocean circulation patterns, but other mechanisms, including changes in biological processes, must have been important as well, and may have been coordinated through links with global air temperature. Further research is required to better quantify the contributions of the various mechanisms, and there remains great potential to use the Last Glacial Maximum and the ensuing global warming as natural experiments from which to learn about climate-driven changes in the marine ecosystem.

Two-million-year-old snapshots of atmospheric gases from Antarctic ice

Over the past eight hundred thousand years, glacial-interglacial cycles oscillated with a period of one hundred thousand years ('100k world'¹). Ice core and ocean sediment data have shown that atmospheric carbon dioxide, Antarctic temperature, deep ocean temperature, and global ice volume correlated strongly with each other in the 100k world^2-6. Between about 2.8 and 1.2 million years ago, glacial cycles were smaller in magnitude and shorter in duration ('40k world'⁷). Proxy data from deep-sea sediments suggest that the variability of atmospheric carbon dioxide in the 40k world was also lower than in the 100k world^8-10, but we do not have direct observations of atmospheric greenhouse gases from this period. Here we report the recovery of stratigraphically discontinuous ice more than two million years old from the Allan Hills Blue Ice Area, East Antarctica. Concentrations of carbon dioxide and methane in ice core samples older than two million years have been altered by respiration, but some younger samples are pristine. The recovered ice cores extend direct observations of atmospheric carbon dioxide, methane, and Antarctic temperature (based on the deuterium/hydrogen isotope ratio δD_ice, a proxy for regional temperature) into the 40k world. All climate properties before eight hundred thousand years ago fall within the envelope of observations from continuous deep Antarctic ice cores that characterize the 100k world. However, the lowest measured carbon dioxide and methane concentrations and Antarctic temperature in the 40k world are well above glacial values from the past eight hundred thousand years. Our results confirm that the amplitudes of glacial-interglacial variations in atmospheric greenhouse gases and Antarctic climate were reduced in the 40k world, and that the transition from the 40k to the 100k world was accompanied by a decline in minimum carbon dioxide concentrations during glacial maxima.

Rabbits and Bergmann’s rule: how cold was Portugal during the last glaciation?

Osteometric data from > 450 modern wild rabbits, mostly from Portugal, Spain and France, show an inverse correlation between their size and the temperature of the environment, in accordance with Bergmann’s rule. Similar measurements made on some 1660 rabbit bones from 14 Portuguese late Pleistocene and Holocene archaeological sites indicate that rabbits became considerably smaller at the Pleistocene–Holocene boundary. Thus, rabbit size varies or varied with temperature both today and in the past. A direct temperature–size relationship was assumed, and the regression of modern rabbit bone size on temperature was then used to calibrate the temperature equivalent for the change of size of rabbit bones in the past. The result indicates a Last Glacial Maximum to present-day difference, Δt°, of 7 or 8 °C. An alternative interpretation that does not assume a direct temperature–size relationship would indicate that the environment in Portugal 15 000–30 000 years ago was similar to that in northern France today.

Earth's radiative imbalance from the Last Glacial Maximum to the present

The energy imbalance at the top of the atmosphere determines the temporal evolution of the global climate, and vice versa changes in the climate system can alter the planetary energy fluxes. This interplay is fundamental to our understanding of Earth's heat budget and the climate system. However, even today, the direct measurement of global radiative fluxes is difficult, such that most assessments are based on changes in the total energy content of the climate system. We apply the same approach to estimate the long-term evolution of Earth's radiative imbalance in the past. New measurements of noble gas-derived mean ocean temperature from the European Project for Ice Coring in Antarctica Dome C ice core covering the last 40,000 y, combined with recent results from the West Antarctic Ice Sheet Divide ice core and the sea-level record, allow us to quantitatively reconstruct the history of the climate system energy budget. The temporal derivative of this quantity must be equal to the planetary radiative imbalance. During the deglaciation, a positive imbalance of typically +0.2 W⋅m^-2 is maintained for ∼10,000 y, however, with two distinct peaks that reach up to 0.4 W⋅m^-2 during times of substantially reduced Atlantic Meridional Overturning Circulation. We conclude that these peaks are related to net changes in ocean heat uptake, likely due to rapid changes in North Atlantic deep-water formation and their impact on the global radiative balance, while changes in cloud coverage, albeit uncertain, may also factor into the picture.

Air-sea disequilibrium enhances ocean carbon storage during glacial periods

The prevailing hypothesis for lower atmospheric carbon dioxide (CO₂) concentrations during glacial periods is an increased efficiency of the ocean's biological pump. However, tests of this and other hypotheses have been hampered by the difficulty to accurately quantify ocean carbon components. Here, we use an observationally constrained earth system model to precisely quantify these components and the role that different processes play in simulated glacial-interglacial CO₂ variations. We find that air-sea disequilibrium greatly amplifies the effects of cooler temperatures and iron fertilization on glacial ocean carbon storage even as the efficiency of the soft-tissue biological pump decreases. These two processes, which have previously been regarded as minor, explain most of our simulated glacial CO₂ drawdown, while ocean circulation and sea ice extent, hitherto considered dominant, emerge as relatively small contributors.

Central Europe temperature constrained by speleothem fluid inclusion water isotopes over the past 14,000 years

The reasons for the early Holocene temperature discrepancy between northern hemispheric model simulations and paleoclimate reconstructions-known as the Holocene temperature conundrum-remain unclear. Using hydrogen isotopes of fluid inclusion water extracted from stalagmites from the Milandre Cave in Switzerland, we established a mid-latitude European mean annual temperature reconstruction for the past 14,000 years. Our Milandre Cave fluid inclusion temperature record (MC-FIT) resembles Greenland and Mediterranean sea surface temperature trends but differs from recent reconstructions obtained from biogenic proxies and climate models. The water isotopes are further synchronized with tropical precipitation records, stressing the Northern Hemisphere signature. Our results support the existence of a European Holocene Thermal Maximum and data-model temperature discrepancies. Moreover, data-data comparison reveals a significant latitudinal temperature gradient within Europe. Last, the MC-FIT record suggests that seasonal biases in the proxies are not the primary cause of the Holocene temperature conundrum.

Highly sensitive temperature sensor based on a polymer-infiltrated Mach-Zehnder interferometer created in graded index fiber

A highly sensitive temperature sensor is proposed and demonstrated based on a UV-curable polymer-infiltrated Mach-Zehnder interferometer (MZI) created in a graded index fiber (GIF). The device was constructed by splicing a half-pitch GIF between two single-mode fibers and creating an inner air cavity in one lateral side of the GIF core by means of femtosecond laser micromachining. The air cavity and the residual GIF core functioned as two interference arms of the MZI. Moreover, the GIF was used as a miniature in-fiber collimator to reduce insertion loss of the air cavity. Experimental results show such an MZI device has a high refractive index (RI) sensitivity of 24611.54 nm/RIU (RI=1.545-1.565). Subsequently, thermo-sensitive polymer liquid was infiltrated into the air cavity, then cured with UV illumination, and annealed at 50°C for 12 h. The infiltrated MZI exhibits a high temperature sensitivity of -13.27  nm/°C. In addition, this MZI also has excellent thermal stability and repeatability, compact structure, low insertion loss, and high fringe visibility. As such, the proposed MZI could be developed for high-accuracy temperature measurements in many areas such as biomedical or oceanographic applications.

Sedimentary record from Patagonia, southern Chile supports cosmic-impact triggering of biomass burning, climate change, and megafaunal extinctions at 12.8 ka

The Younger Dryas (YD) impact hypothesis posits that fragments of a large, disintegrating asteroid/comet struck North America, South America, Europe, and western Asia ~12,800 years ago. Multiple airbursts/impacts produced the YD boundary layer (YDB), depositing peak concentrations of platinum, high-temperature spherules, meltglass, and nanodiamonds, forming an isochronous datum at >50 sites across ~50 million km² of Earth's surface. This proposed event triggered extensive biomass burning, brief impact winter, YD climate change, and contributed to extinctions of late Pleistocene megafauna. In the most extensive investigation south of the equator, we report on a ~12,800-year-old sequence at Pilauco, Chile (~40°S), that exhibits peak YD boundary concentrations of platinum, gold, high-temperature iron- and chromium-rich spherules, and native iron particles rarely found in nature. A major peak in charcoal abundance marks an intense biomass-burning episode, synchronous with dramatic changes in vegetation, including a high-disturbance regime, seasonality in precipitation, and warmer conditions. This is anti-phased with northern-hemispheric cooling at the YD onset, whose rapidity suggests atmospheric linkage. The sudden disappearance of megafaunal remains and dung fungi in the YDB layer at Pilauco correlates with megafaunal extinctions across the Americas. The Pilauco record appears consistent with YDB impact evidence found at sites on four continents.

Using Noble Gases to Assess the Ocean's Carbon Pumps

Natural mechanisms in the ocean, both physical and biological, concentrate carbon in the deep ocean, resulting in lower atmospheric carbon dioxide. The signals of these carbon pumps overlap to create the observed carbon distribution in the ocean, making the individual impact of each pump difficult to disentangle. Noble gases have the potential to directly quantify the physical carbon solubility pump and to indirectly improve estimates of the biological organic carbon pump. Noble gases are biologically inert, can be precisely measured, and span a range of physical properties. We present dissolved neon, argon, and krypton data spanning the Atlantic, Southern, Pacific, and Arctic Oceans. Comparisons between deep-ocean observations and models of varying complexity enable the rates of processes that control the carbon solubility pump to be quantified and thus provide an important metric for ocean model skill. Noble gases also provide a powerful means of assessing air-sea gas exchange parameterizations.

CO2 storage and release in the deep Southern Ocean on millennial to centennial timescales

The cause of changes in atmospheric carbon dioxide (CO₂) during the recent ice ages is yet to be fully explained. Most mechanisms for glacial-interglacial CO₂ change have centred on carbon exchange with the deep ocean, owing to its large size and relatively rapid exchange with the atmosphere¹. The Southern Ocean is thought to have a key role in this exchange, as much of the deep ocean is ventilated to the atmosphere in this region². However, it is difficult to reconstruct changes in deep Southern Ocean carbon storage, so few direct tests of this hypothesis have been carried out. Here we present deep-sea coral boron isotope data that track the pH-and thus the CO₂ chemistry-of the deep Southern Ocean over the past forty thousand years. At sites closest to the Antarctic continental margin, and most influenced by the deep southern waters that form the ocean's lower overturning cell, we find a close relationship between ocean pH and atmospheric CO₂: during intervals of low CO₂, ocean pH is low, reflecting enhanced ocean carbon storage; and during intervals of rising CO₂, ocean pH rises, reflecting loss of carbon from the ocean to the atmosphere. Correspondingly, at shallower sites we find rapid (millennial- to centennial-scale) decreases in pH during abrupt increases in CO₂, reflecting the rapid transfer of carbon from the deep ocean to the upper ocean and atmosphere. Our findings confirm the importance of the deep Southern Ocean in ice-age CO₂ change, and show that deep-ocean CO₂ release can occur as a dynamic feedback to rapid climate change on centennial timescales.

New methods for measuring atmospheric heavy noble gas isotope and elemental ratios in ice core samples

RATIONALE

The global ocean constitutes the largest heat buffer in the global climate system, but little is known about its past changes. The isotopic and elemental ratios of heavy noble gases (krypton and xenon), together with argon and nitrogen in trapped air from ice cores, can be used to reconstruct past mean ocean temperatures (MOTs). Here we introduce two successively developed methods to measure these parameters with a sufficient precision to provide new constraints on past changes in MOT.

METHODS

The air from an 800-g ice sample - containing roughly 80 mL STP air - is extracted and processed to be analyzed on two independent dual-inlet isotope ratio mass spectrometers. The primary isotope ratios (δ¹⁵ N, δ⁴⁰ Ar and δ⁸⁶ Kr values) are obtained with precisions in the range of 1 per meg (0.001‰) per mass unit. The three elemental ratio values δKr/N₂ , δXe/N₂ and δXe/Kr are obtained using sequential (non-simultaneous) peak-jumping, reaching precisions in the range of 0.1-0.3‰.

RESULTS

The latest version of the method achieves a 30% to 50% better precision on the elemental ratios and a twofold better sample throughput than the previous one. The method development uncovered an unexpected source of artefactual gas fractionation in a closed system that is caused by adiabatic cooling and warming of gases (termed adiabatic fractionation) - a potential source of measurement artifacts in other methods.

CONCLUSIONS

The precisions of the three elemental ratios δKr/N₂ , δXe/N₂ and δXe/Kr - which all contain the same MOT information - suggest smaller uncertainties for reconstructed MOTs (±0.3-0.1°C) than previous studies have attained. Due to different sensitivities of the noble gases to changes in MOT, δXe/N₂ provides the best constraints on the MOT under the given precisions followed by δXe/Kr, and δKr/N₂ ; however, using all of them helps to detect methodological artifacts and issues with ice quality.