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.

Nitrous oxide (N2O) is an important greenhouse gas and ozone-depleting substance that has anthropogenic as well as natural marine and terrestrial sources. The tropospheric N2O concentrations have varied substantially in the past in concert with changing climate on glacial-interglacial and millennial timescales. It is not well understood, however, how N2O emissions from marine and terrestrial sources change in response to varying environmental conditions. The distinct isotopic compositions of marine and terrestrial N2O sources can help disentangle the relative changes in marine and terrestrial N2O emissions during past climate variations. Here we present N2O concentration and isotopic data for the last deglaciation, from 16,000 to 10,000 years before present, retrieved from air bubbles trapped in polar ice at Taylor Glacier, Antarctica. With the help of our data and a box model of the N2O cycle, we find a 30 per cent increase in total N2O emissions from the late glacial to the interglacial, with terrestrial and marine emissions contributing equally to the overall increase and generally evolving in parallel over the last deglaciation, even though there is no a priori connection between the drivers of the two sources. However, we find that terrestrial emissions dominated on centennial timescales, consistent with a state-of-the-art dynamic global vegetation and land surface process model that suggests that during the last deglaciation emission changes were strongly influenced by temperature and precipitation patterns over land surfaces. The results improve our understanding of the drivers of natural N2O emissions and are consistent with the idea that natural N2O emissions will probably increase in response to anthropogenic warming.

Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO2 remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO2 fertilization, land use, wildfire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO2 dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO2 and δ13CO2 changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO2 dynamics from 8 ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO2 changes after 122 ka BP. This failure to simulate late-Eemian CO2 dynamics could be a result of the imposed forcings such as prescribed CaCO3 accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO2 dynamics - shallow water CaCO3 accumulation, peat and permafrost carbon dynamics - are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO2 dynamics.