Description: Paleo-climate records and geodynamic modelling indicate the existence of complex interactions between glacial sea level changes, volcanic degassing and atmospheric CO2, which may have modulated the climate system's descent into the last ice age. Between ∼85 and 70 kyr ago, during an interval of decreasing axial tilt, the orbital component in global temperature records gradually declined, while atmospheric CO2, instead of continuing its long-term correlation with Antarctic temperature, remained relatively stable. Here, based on novel global geodynamic models and the joint interpretation of paleo-proxy data as well as biogeochemical simulations, we show that a sea level fall in this interval caused enhanced pressure-release melting in the uppermost mantle, which may have induced a surge in magma and CO2 fluxes from mid-ocean ridges and oceanic hotspot volcanoes. Our results reveal a hitherto unrecognized negative feedback between glaciation and atmospheric CO2 predominantly controlled by marine volcanism on multi-millennial timescales of ∼5,000-15,000 years.

Description: The mobilization of glacial permafrost carbon during the last glacial–interglacial transition has been suggested by indirect evidence to be an additional and significant source of greenhouse gases to the atmosphere, especially at times of rapid sea-level rise. Here we present the first direct evidence for the release of ancient carbon from degrading permafrost in East Asia during the last 17 kyrs, using biomarkers and radiocarbon dating of terrigenous material found in two sediment cores from the Okhotsk Sea. Upscaling our results to the whole Arctic shelf area, we show by carbon cycle simulations that deglacial permafrost-carbon release through sea-level rise likely contributed significantly to the changes in atmospheric CO2 around 14.6 and 11.5 kyrs BP.

Description: Significance: Cold and dry glacial-state climate conditions persisted in the Southern Hemisphere until approximately 17.7 ka, when paleoclimate records show a largely unexplained sharp, nearly synchronous acceleration in deglaciation. Detailed measurements in Antarctic ice cores document exactly at that time a unique, ∼192-y series of massive halogen-rich volcanic eruptions geochemically attributed to Mount Takahe in West Antarctica. Rather than a coincidence, we postulate that halogen-catalyzed stratospheric ozone depletion over Antarctica triggered large-scale atmospheric circulation and hydroclimate changes similar to the modern Antarctic ozone hole, explaining the synchronicity and abruptness of accelerated Southern Hemisphere deglaciation. Abstract: Glacial-state greenhouse gas concentrations and Southern Hemisphere climate conditions persisted until ∼17.7 ka, when a nearly synchronous acceleration in deglaciation was recorded in paleoclimate proxies in large parts of the Southern Hemisphere, with many changes ascribed to a sudden poleward shift in the Southern Hemisphere westerlies and subsequent climate impacts. We used high-resolution chemical measurements in the West Antarctic Ice Sheet Divide, Byrd, and other ice cores to document a unique, ∼192-y series of halogen-rich volcanic eruptions exactly at the start of accelerated deglaciation, with tephra identifying the nearby Mount Takahe volcano as the source. Extensive fallout from these massive eruptions has been found 〉2,800 km from Mount Takahe. Sulfur isotope anomalies and marked decreases in ice core bromine consistent with increased surface UV radiation indicate that the eruptions led to stratospheric ozone depletion. Rather than a highly improbable coincidence, circulation and climate changes extending from the Antarctic Peninsula to the subtropics—similar to those associated with modern stratospheric ozone depletion over Antarctica—plausibly link the Mount Takahe eruptions to the onset of accelerated Southern Hemisphere deglaciation ∼17.7 ka.

Description: We reanalyze existing paleodata of global mean surface temperature ΔTg and radiative forcing ΔR of CO2 and land ice albedo for the last 800,000 years to show that a state‐dependency in paleoclimate sensitivity S, as previously suggested, is only found if ΔTg is based on reconstructions, and not when ΔTg is based on model simulations. Furthermore, during times of decreasing obliquity (periods of land‐ice sheet growth and sea level fall) the multi‐millennial component of reconstructed ΔTg diverges from CO2, while in simulations both variables vary more synchronously, suggesting that the differences during these times are due to relatively low rates of simulated land ice growth and associated cooling. To produce a reconstruction‐based extrapolation of S for the future we exclude intervals with strong ΔTg‐CO2 divergence and find that S is less state‐dependent, or even constant (state‐independent), yielding a mean equilibrium warming of 2–4 K for a doubling of CO2.

Description: Radiocarbon (14C) ages cannot provide absolutely dated chronologies for archaeological or paleoenvironmental studies directly but must be converted to calendar age equivalents using a calibration curve compensating for fluctuations in atmospheric 14C concentration. Although calibration curves are constructed from independently dated archives, they invariably require revision as new data become available and our understanding of the Earth system improves. In this volume the international 14C calibration curves for both the Northern and Southern Hemispheres, as well as for the ocean surface layer, have been updated to include a wealth of new data and extended to 55,000 cal BP. Based on tree rings, IntCal20 now extends as a fully atmospheric record to ca. 13,900 cal BP. For the older part of the timescale, IntCal20 comprises statistically integrated evidence from floating tree-ring chronologies, lacustrine and marine sediments, speleothems, and corals. We utilized improved evaluation of the timescales and location variable 14C offsets from the atmosphere (reservoir age, dead carbon fraction) for each dataset. New statistical methods have refined the structure of the calibration curves while maintaining a robust treatment of uncertainties in the 14C ages, the calendar ages and other corrections. The inclusion of modeled marine reservoir ages derived from a three-dimensional ocean circulation model has allowed us to apply more appropriate reservoir corrections to the marine 14C data rather than the previous use of constant regional offsets from the atmosphere. Here we provide an overview of the new and revised datasets and the associated methods used for the construction of the IntCal20 curve and explore potential regional offsets for tree-ring data. We discuss the main differences with respect to the previous calibration curve, IntCal13, and some of the implications for archaeology and geosciences ranging from the recent past to the time of the extinction of the Neanderthals.

Description: The concentration of radiocarbon (14C) differs between ocean and atmosphere. Radiocarbon determinations from samples which obtained their 14C in the marine environment therefore need a marine-specific calibration curve and cannot be calibrated directly against the atmospheric-based IntCal20 curve. This paper presents Marine20, an update to the internationally agreed marine radiocarbon age calibration curve that provides a non-polar global-average marine record of radiocarbon from 0–55 cal kBP and serves as a baseline for regional oceanic variation. Marine20 is intended for calibration of marine radiocarbon samples from non-polar regions; it is not suitable for calibration in polar regions where variability in sea ice extent, ocean upwelling and air-sea gas exchange may have caused larger changes to concentrations of marine radiocarbon. The Marine20 curve is based upon 500 simulations with an ocean/atmosphere/biosphere box-model of the global carbon cycle that has been forced by posterior realizations of our Northern Hemispheric atmospheric IntCal20 14C curve and reconstructed changes in CO2 obtained from ice core data. These forcings enable us to incorporate carbon cycle dynamics and temporal changes in the atmospheric 14C level. The box-model simulations of the global-average marine radiocarbon reservoir age are similar to those of a more complex three-dimensional ocean general circulation model. However, simplicity and speed of the box model allow us to use a Monte Carlo approach to rigorously propagate the uncertainty in both the historic concentration of atmospheric 14C and other key parameters of the carbon cycle through to our final Marine20 calibration curve. This robust propagation of uncertainty is fundamental to providing reliable precision for the radiocarbon age calibration of marine based samples. We make a first step towards deconvolving the contributions of different processes to the total uncertainty; discuss the main differences of Marine20 from the previous age calibration curve Marine13; and identify the limitations of our approach together with key areas for further work. The updated values for ΔR, the regional marine radiocarbon reservoir age corrections required to calibrate against Marine20, can be found at the data base http://calib.org/marine/.

Description: The importance of volcanic CO2 release, continental weathering, and coral reef growth on the global carbon cycle has been highlighted by several different studies. Based on these independent approaches, we here revisit the last 800 kyr with the box model BICYCLE, which has been extended to be able to address these solid Earth contributions to the carbon cycle. We show that the volcanic outgassing of CO2 as a function of sea level change from mid‐ocean ridges and hot spot island volcanoes cannot be the generic process that leads during phases of falling obliquity to a sea level‐CO2 decoupling as has been suggested before. The combined contribution from continental and marine volcanism, if both lagging sea level change by 4 kyr, might have added up to 13 ppm to the glacial/interglacial CO2 rise. Coral reef growth as suggested by an independent model is during glacial terminations about an order of magnitude too high to be reconciled with meaningful carbon cycle dynamics. Global riverine input of bicarbonate caused by silicate and carbonate weathering is suggested to have been stable over Termination I. However, if weathering fluxes are changed by up to 50% in sensitivity experiments, the corresponding bicarbonate input might contribute less than 20 ppm to the deglacial atmospheric CO2 rise. The overall agreement of results with the new process‐based sediment module and the previously applied time‐delayed response function to mimic carbonate compensation gives confidence in the results obtained in previous applications of the BICYCLE model without solid Earth processes.

Description: Numerical models are important tools for understanding the processes and feedbacks in the Earth system, including those involving changes in atmospheric CO2 (CO2,atm) concentrations. Here, we compile 55 published model studies (consisting of 778 individual simulations) that assess the impact of six forcing mechanisms on millennial-scale CO2,atm variations: changes in freshwater supply to the North Atlantic and Southern Ocean, the strength and position of the southern-hemisphere westerlies, Antarctic sea ice extent, and aeolian dust fluxes. We generally find agreement on the direction of simulated CO2,atm change across simulations, but the amplitude of change is inconsistent, primarily due to the different complexities of the model representation of Earth system processes. When freshwater is added to the North Atlantic, a reduced Atlantic Meridional Overturning Circulation (AMOC) is generally accompanied by an increase in Southern Ocean- and Pacific overturning, reduced Antarctic sea ice extent, spatially varying export production, and changes in carbon storage in the Atlantic (rising), in other ocean basins (generally decreasing) and on land (more varied). Positive or negative CO2,atm changes are simulated during AMOC minima due to a spatially and temporally varying dominance of individual terrestrial and oceanic drivers (and compensating effects between them) across the different models. In contrast, AMOC recoveries are often accompanied by rising CO2,atm levels, which are mostly driven by ocean carbon release (albeit from different regions). The magnitude of simulated CO2,atm rise broadly scales with the duration of the AMOC perturbation (i.e., the stadial length). When freshwater is added to the Southern Ocean, reduced deep-ocean ventilation drives a CO2,atm drop via reduced carbon release from the Southern Ocean. Although the impacts of shifted southern-hemisphere westerlies are inconsistent across model simulations, their intensification raises CO2,atm via enhanced Southern Ocean Ekman pumping. Increased supply of aeolian dust to the ocean, and thus iron fertilisation of marine productivity, consistently lowers modelled CO2,atm concentrations via more efficient nutrient utilisation. The magnitude of CO2,atm change in response to dust flux variations, however, largely depends on the complexity of models' marine ecosystem and iron cycle. This especially applies to simulations forced by Antarctic sea ice changes, in which the direction of simulated CO2,atm change varies greatly across model hierarchies. Our compilation highlights that no single (forcing) mechanism can explain observed past millennial-scale CO2,atm variability, and identifies important future needs in coupled carbon cycle-climate modelling to better understand the mechanisms governing CO2,atm changes in the past.

Description: During the last deglaciation (18–8 kyr BP), shelf flooding and warming presumably led to a large-scale decomposition of permafrost soils in the mid-to-high latitudes of the Northern Hemisphere. Microbial degradation of old organic matter released from the decomposing permafrost potentially contributed to the deglacial rise in atmospheric CO2 and also to the declining atmospheric radiocarbon contents (Δ14C). The significance of permafrost for the atmospheric carbon pool is not well understood as the timing of the carbon activation is poorly constrained by proxy data. Here, we trace the mobilization of organic matter from permafrost in the Pacific sector of Beringia over the last 22 kyr using mass-accumulation rates and radiocarbon signatures of terrigenous biomarkers in four sediment cores from the Bering Sea and the Northwest Pacific. We find that pronounced reworking and thus the vulnerability of old organic carbon to remineralization commenced during the early deglaciation (~16.8 kyr BP) when meltwater runoff in the Yukon River intensified riverbank erosion of permafrost soils and fluvial discharge. Regional deglaciation in Alaska additionally mobilized significant fractions of fossil, petrogenic organic matter at this time. Permafrost decomposition across Beringia's Pacific sector occurred in two major pulses that match the Bølling-Allerød and Preboreal warm spells and rapidly initiated within centuries. The carbon mobilization likely resulted from massive shelf flooding during meltwater pulses 1A (~14.6 kyr BP) and 1B (~11.5 kyr BP) followed by permafrost thaw in the hinterland. Our findings emphasize that coastal erosion was a major control to rapidly mobilize permafrost carbon along Beringia's Pacific coast at ~14.6 and ~11.5 kyr BP implying that shelf flooding in Beringia may partly explain the centennial-scale rises in atmospheric CO2 at these times. Around 16.5 kyr BP, the mobilization of old terrigenous organic matter caused by meltwater-floods may have additionally contributed to increasing CO2 levels.

Description: Continuous records of the atmospheric greenhouse gases (GHGs) CO2, CH4, and N2O are necessary input data for transient climate simulations, and their associated radiative forcing represents important components in analyses of climate sensitivity and feedbacks. Since the available data from ice cores are discontinuous and partly ambiguous, a well-documented decision process during data compilation followed by some interpolating post-processing is necessary to obtain those desired time series. Here, we document our best possible data compilation of published ice core records and recent measurements on firn air and atmospheric samples spanning the interval from the penultimate glacial maximum ( ∼  156 kyr BP) to the beginning of the year 2016 CE. We use the most recent age scales for the ice core data and apply a smoothing spline method to translate the discrete and irregularly spaced data points into continuous time series. These splines are then used to compute the radiative forcing for each GHG using well-established, simple formulations. We compile only a Southern Hemisphere record of CH4 and discuss how much larger a Northern Hemisphere or global CH4 record might have been due to its interpolar difference. The uncertainties of the individual data points are considered in the spline procedure. Based on the given data resolution, time-dependent cutoff periods of the spline, defining the degree of smoothing, are prescribed, ranging from 5000 years for the less resolved older parts of the records to 4 years for the densely sampled recent years. The computed splines seamlessly describe the GHG evolution on orbital and millennial timescales for glacial and glacial–interglacial variations and on centennial and decadal timescales for anthropogenic times. Data connected with this paper, including raw data and final splines, are available at doi:10.1594/PANGAEA.871273.