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 simulations1,2 and proxy data3-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 assimilation9,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 reconstructions6,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 changes11, our reanalysis indicates that both the rate and magnitude of modern warming are unusual relative to the changes of the past 24 thousand years.

Industrial-era decline in subarctic Atlantic productivity.

Marine phytoplankton have a crucial role in the modulation of marine-based food webs1, fishery yields2 and the global drawdown of atmospheric carbon dioxide3. However, owing to sparse measurements before satellite monitoring in the twenty-first century, the long-term response of planktonic stocks to climate forcing is unknown. Here, using a continuous, multi-century record of subarctic Atlantic marine productivity, we show that a marked 10 ± 7% decline in net primary productivity has occurred across this highly productive ocean basin over the past two centuries. We support this conclusion by the application of a marine-productivity proxy, established using the signal of the planktonic-derived aerosol methanesulfonic acid, which is commonly identified across an array of Greenlandic ice cores. Using contemporaneous satellite-era observations, we demonstrate the use of this signal as a robust and high-resolution proxy for past variations in spatially integrated marine productivity. We show that the initiation of declining subarctic Atlantic productivity broadly coincides with the onset of Arctic surface warming4, and that productivity strongly covaries with regional sea-surface temperatures and basin-wide gyre circulation strength over recent decades. Taken together, our results suggest that the decline in industrial-era productivity may be evidence of the predicted5 collapse of northern Atlantic planktonic stocks in response to a weakened Atlantic Meridional Overturning Circulation6-8. Continued weakening of this Atlantic Meridional Overturning Circulation, as projected for the twenty-first century9,10, may therefore result in further productivity declines across this globally relevant region.

Nonlinear rise in Greenland runoff in response to post-industrial Arctic warming.

The Greenland ice sheet (GrIS) is a growing contributor to global sea-level rise1, with recent ice mass loss dominated by surface meltwater runoff2,3. Satellite observations reveal positive trends in GrIS surface melt extent4, but melt variability, intensity and runoff remain uncertain before the satellite era. Here we present the first continuous, multi-century and observationally constrained record of GrIS surface melt intensity and runoff, revealing that the magnitude of recent GrIS melting is exceptional over at least the last 350 years. We develop this record through stratigraphic analysis of central west Greenland ice cores, and demonstrate that measurements of refrozen melt layers in percolation zone ice cores can be used to quantifiably, and reproducibly, reconstruct past melt rates. We show significant (P < 0.01) and spatially extensive correlations between these ice-core-derived melt records and modelled melt rates5,6 and satellite-derived melt duration4 across Greenland more broadly, enabling the reconstruction of past ice-sheet-scale surface melt intensity and runoff. We find that the initiation of increases in GrIS melting closely follow the onset of industrial-era Arctic warming in the mid-1800s, but that the magnitude of GrIS melting has only recently emerged beyond the range of natural variability. Owing to a nonlinear response of surface melting to increasing summer air temperatures, continued atmospheric warming will lead to rapid increases in GrIS runoff and sea-level contributions.

North Atlantic jet stream projections in the context of the past 1,250 years.

Reconstruction of the North Atlantic jet stream (NAJ) presents a critical, albeit largely unconstrained, paleoclimatic target. Models suggest northward migration and changing variance of the NAJ under 21st-century warming scenarios, but assessing the significance of such projections is hindered by a lack of long-term observations. Here, we incorporate insights from an ensemble of last-millennium water isotope-enabled climate model simulations and a wide array of mean annual water isotope ([Formula: see text]O) and annually accumulated snowfall records from Greenland ice cores to reconstruct North Atlantic zonal-mean zonal winds back to the 8th century CE. Using this reconstruction we provide preobservational constraints on both annual mean NAJ position and intensity to show that late 20th- and early 21st-century NAJ variations were likely not unique relative to natural variability. Rather, insights from our 1,250 year reconstruction highlight the overwhelming role of natural variability in thus far masking the response of midlatitude atmospheric dynamics to anthropogenic forcing, consistent with recent large-ensemble transient modeling experiments. This masking is not projected to persist under high greenhouse gas emissions scenarios, however, with model projected annual mean NAJ position emerging as distinct from the range of reconstructed natural variability by as early as 2060 CE.

Calibration of Arctic ice core bromine enrichment records for past sea ice reconstructions.

Bromine in ice cores has been proposed as a qualitative sea ice proxy to produce sea ice reconstructions for the polar regions. Here we report the first statistical validation of this proxy with satellite sea ice observations by combining bromine enrichment (with respect to seawater, Brenr) records from three Greenlandic ice cores (SIGMA-A, NU and RECAP) with satellite sea ice imagery, over three decades. We find that during the 1984-2016 satellite-era, ice core Brenr values are significantly correlated with first-year sea ice formed in the Baffin Bay and Labrador Sea supporting that the gas-phase bromine enrichment processes, preferentially occurring over the sea ice surface, are the main driver for the Brenr signal in ice cores. Moreover, in assessing Brenr's capability to record historical sea ice variability, we compare 20th-century Arctic Sea ice historical and proxy records with our reconstructions, based on an autoregressive-moving-average (ARMA) model, finding overall good agreement. While further enhancements are warranted, including site-specific calibrations and a comprehensive investigation into bromine transport-related concerns, this study presents a new method to quantitatively reconstruct past seasonal sea ice variability through bromine enrichment in ice cores.

Last Glacial Maximum pattern effects reduce climate sensitivity estimates.

Here, we show that the Last Glacial Maximum (LGM) provides a stronger constraint on equilibrium climate sensitivity (ECS), the global warming from increasing greenhouse gases, after accounting for temperature patterns. Feedbacks governing ECS depend on spatial patterns of surface temperature ("pattern effects"); hence, using the LGM to constrain future warming requires quantifying how temperature patterns produce different feedbacks during LGM cooling versus modern-day warming. Combining data assimilation reconstructions with atmospheric models, we show that the climate is more sensitive to LGM forcing because ice sheets amplify extratropical cooling where feedbacks are destabilizing. Accounting for LGM pattern effects yields a median modern-day ECS of 2.4°C, 66% range 1.7° to 3.5°C (1.4° to 5.0°C, 5 to 95%), from LGM evidence alone. Combining the LGM with other lines of evidence, the best estimate becomes 2.9°C, 66% range 2.4° to 3.5°C (2.1° to 4.1°C, 5 to 95%), substantially narrowing uncertainty compared to recent assessments.

Reconciling ice core CO2 and land-use change following New World-Old World contact.

Ice core records of carbon dioxide (CO2) throughout the last 2000 years provide context for the unprecedented anthropogenic rise in atmospheric CO2 and insights into global carbon cycle dynamics. Yet the atmospheric history of CO2 remains uncertain in some time intervals. Here we present measurements of CO2 and methane (CH4) in the Skytrain ice core from 1450 to 1700 CE. Results suggest a sudden decrease in CO2 around 1610 CE in one widely used record may be an artefact of a small number of anomalously low values. Our analysis supports a more gradual decrease in CO2 of 0.5 ppm per decade from 1516 to 1670 CE, with an inferred land carbon sink of 2.6 PgC per decade. This corroborates modelled scenarios of large-scale reorganisation of land use in the Americas following New World-Old World contact, whereas a rapid decrease in CO2 at 1610 CE is incompatible with even the most extreme land-use change scenarios.