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

Robust acceleration of Earth system heating observed over the past six decades

Global heating of the Earth system is unequivocal. However, detecting an acceleration of Earth heating has remained elusive to date, despite suggestive evidence of a potential increase in heating rates. In this study, we demonstrate that since 1960, the warming of the world ocean has accelerated at a relatively consistent pace of 0.15 ± 0.05 (W/m2)/decade, while the land, cryosphere, and atmosphere have exhibited an accelerated pace of 0.013 ± 0.003 (W/m2)/decade. This has led to a substantial increase in ocean warming, with a magnitude of 0.91 ± 0.80 W/m2 between the decades 1960-1970 and 2010-2020, which overlies substantial decadal-scale variability in ocean warming of up to 0.6 W/m2. Our findings withstand a wide range of sensitivity analyses and are consistent across different observation-based datasets. The long-term acceleration of Earth warming aligns qualitatively with the rise in CO2 concentrations and the decline in aerosol concentration during the same period, but further investigations are necessary to properly attribute these changes.

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.

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.

A global long-term ocean surface daily/0.05° net radiation product from 1983–2020

The all-wave net radiation (Rn) on the ocean surface characterizes the available radiative energy balance and is important to understand the Earth’s climate system. Considering the shortcomings of available ocean surface Rn datasets (e.g., coarse spatial resolutions, discrepancy in accuracy, inconsistency, and short duration), a new long-term global daily Rn product at a spatial resolution of 0.05° from 1983 to 2020, as part of the Global High Resolution Ocean Surface Energy (GHOSE) products suite, was generated in this study by fusing several existing datasets including satellite and reanalysis products based on the comprehensive in situ measurements from 68 globally distributed moored buoy sites. Evaluation against in-situ measurements shows the root mean square difference, mean bias error and correlation coefficient squared of 23.56 Wm−2, 0.88 Wm−2 and 0.878. The global average ocean surface Rn over 1983–2020 is estimated to be 119.71 ± 2.78 Wm−2 with a significant increasing rate of 0.16 Wm−2 per year. GHOSE Rn product can be valuable for oceanic and climatic studies.

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.

A New Empirical Estimation Scheme for Daily Net Radiation at the Ocean Surface

Ocean surface net radiation (Rn) is significant in research on the Earth’s heat balance systems, air–sea interactions, and other applications. However, there have been few studies on Rn until now. Based on radiative and meteorological measurements collected from 66 globally distributed moored buoys, it was found that Rn was dominated by downward shortwave radiation (Rg↓) when the length ratio of daytime (LRD) was greater than 0.4 but dominated by downward longwave radiation (Rl↓) for the other cases (LRD ≤ 0.4). Therefore, an empirical scheme that includes two conditional models named Case 1 (LRD > 0.4) utilizing Rg↓ as a major input and Case 2 (LRD ≤ 0.4) utilizing Rl↓ as a major input for Rn estimation was successfully developed. After validation against in situ Rn, the performance of the empirical scheme was satisfactory with an overall R2 value of 0.972, an RMSE of 9.768 Wm−2, and a bias of −0.092 Wm−2. Specifically, the accuracies of the two conditional models were also very good, with RMSEs of 9.805 and 2.824 Wm−2 and biases of −0.095 and 0.346 Wm−2 for the Case 1 and Case 2 models, respectively. However, due to the limited number of available samples, the performances of these new models were poor in coastal and high-latitude areas, and the models did not work when the LRD was too small (i.e., LRD < 0.3). Overall, the newly developed empirical scheme for Rn estimation has strong potential to be widely used in practical use because of its simple format and high accuracy.

20th century cooling of the deep ocean contributed to delayed acceleration of Earth's energy imbalance

The historical evolution of Earth’s energy imbalance can be quantified by changes in the global ocean heat content. However, historical reconstructions of ocean heat content often neglect a large volume of the deep ocean, due to sparse observations of ocean temperatures below 2000 m. Here, we provide a global reconstruction of historical changes in full-depth ocean heat content based on interpolated subsurface temperature data using an autoregressive artificial neural network, providing estimates of total ocean warming for the period 1946-2019. We find that cooling of the deep ocean and a small heat gain in the upper ocean led to no robust trend in global ocean heat content from 1960-1990, implying a roughly balanced Earth energy budget within −0.16 to 0.06 W m⁻² over most of the latter half of the 20th century. However, the past three decades have seen a rapid acceleration in ocean warming, with the entire ocean warming from top to bottom at a rate of 0.63 ± 0.13 W m⁻². These results suggest a delayed onset of a positive Earth energy imbalance relative to previous estimates, although large uncertainties remain.

Cooling of the global ocean below 2000 m counteracted some of the warming of the shallow ocean over much of the late 20th century. Here the authors show that this trend has shifted to warming, leading the deep ocean to absorb a meaningful fraction of total ocean heat during the 21st century.

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.