skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: Closure of Earth’s Global Seasonal Cycle of Energy Storage

The global seasonal cycle of energy in Earth’s climate system is quantified using observations and reanalyses. After removing long-term trends, net energy entering and exiting the climate system at the top of the atmosphere (TOA) should agree with the sum of energy entering and exiting the ocean, atmosphere, land, and ice over the course of an average year. Achieving such a balanced budget with observations has been challenging. Disagreements have been attributed previously to sparse observations in the high-latitude oceans. However, limiting the local vertical integration of new global ocean heat content estimates to the depth to which seasonal heat energy is stored, rather than integrating to 2000 m everywhere as done previously, allows closure of the global seasonal energy budget within statistical uncertainties. The seasonal cycle of energy storage is largest in the ocean, peaking in April because ocean area is largest in the Southern Hemisphere and the ocean’s thermal inertia causes a lag with respect to the austral summer solstice. Seasonal cycles in energy storage in the atmosphere and land are smaller, but peak in July and September, respectively, because there is more land in the Northern Hemisphere, and the land has more thermal inertia than the atmosphere. Global seasonal energy storage by ice is small, so the atmosphere and land partially offset ocean energy storage in the global integral, with their sum matching time-integrated net global TOA energy fluxes over the seasonal cycle within uncertainties, and both peaking in April.

more » « less
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Surveys in Geophysics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The thermal component of oceanic eddy available potential energy (EPE) generation due to air‐sea interaction is proportional to the product of anomalous sea surface temperature (SST) and net air‐sea heat flux (SHF). In this study we assess EPE generation and its timescale and space‐scale dependence from observations and a high‐resolution coupled climate model. A dichotomy exists in the literature with respect to the sign of this term, that is, whether it is a source or a sink of EPE. We resolve this dichotomy by partitioning the SST and net heat flux into climatological mean, climatological seasonal cycle, and remaining transient contributions, thereby separating the mesoscale eddy variability from the forced seasonal cycle. In this decomposition the mesoscale air‐sea SST‐SHF feedbacks act as a 0.1 TW global sink of EPE. In regions of the ocean with a large seasonal cycle, for example, midlatitudes of the Northern Hemisphere, the EPE generation by the forced seasonal cycle exceeds the mesoscale variability sink, such that the global generation by seasonal plus eddy variability acts as a 0.8 TW source. EPE destruction is largest in the midlatitude western boundary currents due to mesoscale air‐sea interaction and in the tropical Pacific where SST variability is due mainly to the El Niño–Southern Oscillation. The EPE sink in western boundary currents is spatially aligned with SST gradients and offset to the poleward side of currents, while the mean and seasonal generation are aligned with the warm core of the current. By successively smoothing the data in space and time we find that half of the EPE sink is confined to timescales less than annual and length scales less than 2°, within the oceanic mesoscale band.

    more » « less
  2. Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land-use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO2 products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2022, EFOS increased by 0.9 % relative to 2021, with fossil emissions at 9.9±0.5 Gt C yr−1 (10.2±0.5 Gt C yr−1 when the cement carbonation sink is not included), and ELUC was 1.2±0.7 Gt C yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1±0.8 Gt C yr−1 (40.7±3.2 Gt CO2 yr−1). Also, for 2022, GATM was 4.6±0.2 Gt C yr−1 (2.18±0.1 ppm yr−1; ppm denotes parts per million), SOCEAN was 2.8±0.4 Gt C yr−1, and SLAND was 3.8±0.8 Gt C yr−1, with a BIM of −0.1 Gt C yr−1 (i.e. total estimated sources marginally too low or sinks marginally too high). The global atmospheric CO2 concentration averaged over 2022 reached 417.1±0.1 ppm. Preliminary data for 2023 suggest an increase in EFOS relative to 2022 of +1.1 % (0.0 % to 2.1 %) globally and atmospheric CO2 concentration reaching 419.3 ppm, 51 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959–2022, with a near-zero overall budget imbalance, although discrepancies of up to around 1 Gt C yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living-data update documents changes in methods and data sets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at (Friedlingstein et al., 2023).

    more » « less
  3. Abstract

    In situ observation networks and reanalyses products of the state of the atmosphere and upper ocean show well-defined, large-scale patterns of coupled climate variability on time scales ranging from seasons to several decades. We summarize these phenomena and their physics, which have been revealed by analysis of observations, by experimentation with uncoupled and coupled atmosphere and ocean models with a hierarchy of complexity, and by theoretical developments. We start with a discussion of the seasonal cycle in the equatorial tropical Pacific and Atlantic Oceans, which are clearly affected by coupling between the atmosphere and the ocean. We then discuss the tropical phenomena that only exist because of the coupling between the atmosphere and the ocean: the Pacific and Atlantic meridional modes, the El Niño–Southern Oscillation (ENSO) in the Pacific, and a phenomenon analogous to ENSO in the Atlantic. For ENSO, we further discuss the sources of irregularity and asymmetry between warm and cold phases of ENSO, and the response of ENSO to forcing. Fundamental to variability on all time scales in the midlatitudes of the Northern Hemisphere are preferred patterns of uncoupled atmospheric variability that exist independent of any changes in the state of the ocean, land, or distribution of sea ice. These patterns include the North Atlantic Oscillation (NAO), the North Pacific Oscillation (NPO), and the Pacific–North American (PNA) pattern; they are most active in wintertime, with a temporal spectrum that is nearly white. Stochastic variability in the NPO, PNA, and NAO force the ocean on days to interannual times scales by way of turbulent heat exchange and Ekman transport, and on decadal and longer time scales by way of wind stress forcing. The PNA is partially responsible for the Pacific decadal oscillation; the NAO is responsible for an analogous phenomenon in the North Atlantic subpolar gyre. In models, stochastic forcing by the NAO also gives rise to variability in the strength of the Atlantic meridional overturning circulation (AMOC) that is partially responsible for multidecadal anomalies in the North Atlantic climate known as the Atlantic multidecadal oscillation (AMO); observations do not yet exist to adequately determine the physics of the AMO. We review the progress that has been made in the past 50 years in understanding each of these phenomena and the implications for short-term (seasonal-to-interannual) climate forecasts. We end with a brief discussion of advances of things that are on the horizon, under the rug, and over the rainbow.

    more » « less
  4. Abstract. Accurate assessment of anthropogenic carbon dioxide(CO2) emissions and their redistribution among the atmosphere,ocean, and terrestrial biosphere – the “global carbon budget” – isimportant to better understand the global carbon cycle, support thedevelopment of climate policies, and project future climate change. Here wedescribe data sets and methodology to quantify the five major components ofthe global carbon budget and their uncertainties. Fossil CO2emissions (EFF) are based on energy statistics and cementproduction data, while emissions from land use and land-use change (ELUC),mainly deforestation, are based on land use and land-use change data andbookkeeping models. Atmospheric CO2 concentration is measureddirectly and its growth rate (GATM) is computed from the annualchanges in concentration. The ocean CO2 sink (SOCEAN)and terrestrial CO2 sink (SLAND) are estimated withglobal process models constrained by observations. The resulting carbonbudget imbalance (BIM), the difference between the estimatedtotal emissions and the estimated changes in the atmosphere, ocean, andterrestrial biosphere, is a measure of imperfect data and understanding ofthe contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2008–2017), EFF was9.4±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.7±0.02 GtC yr−1,SOCEAN 2.4±0.5 GtC yr−1, and SLAND 3.2±0.8 GtC yr−1, with a budget imbalance BIM of0.5 GtC yr−1 indicating overestimated emissions and/or underestimatedsinks. For the year 2017 alone, the growth in EFF was about 1.6 %and emissions increased to 9.9±0.5 GtC yr−1. Also for 2017,ELUC was 1.4±0.7 GtC yr−1, GATM was 4.6±0.2 GtC yr−1, SOCEAN was 2.5±0.5 GtC yr−1, and SLAND was 3.8±0.8 GtC yr−1,with a BIM of 0.3 GtC. The global atmosphericCO2 concentration reached 405.0±0.1 ppm averaged over 2017.For 2018, preliminary data for the first 6–9 months indicate a renewedgrowth in EFF of +2.7 % (range of 1.8 % to 3.7 %) basedon national emission projections for China, the US, the EU, and India andprojections of gross domestic product corrected for recent changes in thecarbon intensity of the economy for the rest of the world. The analysispresented here shows that the mean and trend in the five components of theglobal carbon budget are consistently estimated over the period of 1959–2017,but discrepancies of up to 1 GtC yr−1 persist for the representationof semi-decadal variability in CO2 fluxes. A detailed comparisonamong individual estimates and the introduction of a broad range ofobservations show (1) no consensus in the mean and trend in land-use changeemissions, (2) a persistent low agreement among the different methods onthe magnitude of the land CO2 flux in the northern extra-tropics,and (3) an apparent underestimation of the CO2 variability by oceanmodels, originating outside the tropics. This living data update documentschanges in the methods and data sets used in this new global carbon budgetand the progress in understanding the global carbon cycle compared withprevious publications of this data set (Le Quéré et al., 2018, 2016,2015a, b, 2014, 2013). All results presented here can be downloaded from 
    more » « less
  5. Abstract

    In the tropics, the absorbed solar radiation is larger than the outgoing longwave radiation, while the opposite is true at high latitudes. This basic fact implies a poleward energy transport (PET) in both hemispheres, which is accomplished by the atmosphere and oceans. The magnitude of PET is determined by the top of atmosphere gradient in the net radiation flux, and small changes to this quantity must change the total PET in the absence of changes in heat uptake. We analyze a large ensemble of 50 historic climate simulations from the CESM LENS2 project and find a significant PET anomaly in the latter half of the twentieth century. The temporal evolution of this anomaly—with a rapid increase after 1950, a peak near 1975, and a rapid decrease in the 1990s—mirrors the historic trend of sulfur dioxide (SO2, a significant aerosol predecessor) emissions from Europe and North America. This anomaly also appears in an analysis of the PET calculated from ERA5 reanalyses and from the CESM2 Single Forcing Large Ensemble. Consistent with previous studies, we find that historic SO2emissions from Europe and North America brightened clouds, which reflected additional solar radiation back to space in the midlatitudes: this shortwave anomaly increased the meridional gradient in the net TOA radiation flux and induced an anomalous northward energy transport. Finally, our results suggest that cryosphere processes become an additional important factor in setting the PET anomaly during the first years of the twenty-first century by contributing to the difference in absorbed solar radiation between hemispheres alongside cloud radiative effects.

    significance statement

    In this study, we analyze a large group of climate model simulations from 1850 to 2014 and find that this historical pollution changed the way that heat was transported from the tropics to Earth’s poles. We also find this change in heat transport when we analyzed an atmospheric reanalysis, which is a historical dataset that combines many meteorological observations into a best estimate of the past climate state. This extra reflection of sunlight from polluted clouds cooled the Northern Hemisphere, and we hypothesize that this cooling caused more heat transport out of the tropics. Last, we find that similar pollution emitted from China and India in more recent decades has not led to a change in Earth’s heat transport because of counteracting changes in snow and ice in the Northern Hemisphere.

    more » « less