The algorithm to produce the Clouds and the Earth’s Radiant Energy System (CERES) Edition 4.0 (Ed4) Energy Balanced and Filled (EBAF)-surface data product is explained. The algorithm forces computed top-of-atmosphere (TOA) irradiances to match with Ed4 EBAF-TOA irradiances by adjusting surface, cloud, and atmospheric properties. Surface irradiances are subsequently adjusted using radiative kernels. The adjustment process is composed of two parts: bias correction and Lagrange multiplier. The bias in temperature and specific humidity between 200 and 500 hPa used for the irradiance computation is corrected based on observations by Atmospheric Infrared Sounder (AIRS). Similarly, the bias in the cloud fraction is corrected based on observations by Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) and CloudSat. Remaining errors in surface, cloud, and atmospheric properties are corrected in the Lagrange multiplier process. Ed4 global annual mean (January 2005 through December 2014) surface net shortwave (SW) and longwave (LW) irradiances increase by 1.3 W m−2 and decrease by 0.2 W m−2, respectively, compared to EBAF Edition 2.8 (Ed2.8) counterparts (the previous version), resulting in an increase in net SW + LW surface irradiance of 1.1 W m−2. The uncertainty in surface irradiances over ocean, land, and polar regions at various spatial scales are estimated. The uncertainties in all-sky global annual mean upward and downward shortwave irradiance are 3 and 4 W m−2, respectively, and the uncertainties in upward and downward longwave irradiance are 3 and 6 W m−2, respectively. With an assumption of all errors being independent, the uncertainty in the global annual mean surface LW + SW net irradiance is 8 W m−2.

Earth's Energy Imbalance (EEI) is a relatively small (presently ∼0.3%) difference between global mean solar radiation absorbed and thermal infrared radiation emitted to space. EEI is set by natural and anthropogenic climate forcings and the climate system's response to those forcings. It is also influenced by internal variations within the climate system. Most of EEI warms the ocean; the remainder heats the land, melts ice, and warms the atmosphere. We show that independent satellite and in situ observations each yield statistically indistinguishable decadal increases in EEI from mid‐2005 to mid‐2019 of 0.50 ± 0.47 W m⁻² decade⁻¹ (5%–95% confidence interval). This trend is primarily due to an increase in absorbed solar radiation associated with decreased reflection by clouds and sea‐ice and a decrease in outgoing longwave radiation (OLR) due to increases in trace gases and water vapor. These changes combined exceed a positive trend in OLR due to increasing global mean temperatures.

This study examines changes in Earth’s energy budget during and after the global warming “pause” (or “hiatus”) using observations from the Clouds and the Earth’s Radiant Energy System. We find a marked 0.83 ± 0.41 Wm⁻² reduction in global mean reflected shortwave (SW) top-of-atmosphere (TOA) flux during the three years following the hiatus that results in an increase in net energy into the climate system. A partial radiative perturbation analysis reveals that decreases in low cloud cover are the primary driver of the decrease in SW TOA flux. The regional distribution of the SW TOA flux changes associated with the decreases in low cloud cover closely matches that of sea-surface temperature warming, which shows a pattern typical of the positive phase of the Pacific Decadal Oscillation. Large reductions in clear-sky SW TOA flux are also found over much of the Pacific and Atlantic Oceans in the northern hemisphere. These are associated with a reduction in aerosol optical depth consistent with stricter pollution controls in China and North America. A simple energy budget framework is used to show that TOA radiation (particularly in the SW) likely played a dominant role in driving the marked increase in temperature tendency during the post-hiatus period.

Both to reconcile the large range in satellite-based estimates of the aerosol direct radiative effect (DRE) and to optimize the design of future observing systems, this study builds a framework for assessing aerosol DRE uncertainty. Shortwave aerosol DRE radiative kernels (Jacobians) were derived using the MERRA-2 reanalysis data. These radiative kernels give the differential response of the aerosol DRE to perturbations in the aerosol extinction coefficient, aerosol single-scattering albedo, aerosol asymmetry factor, surface albedo, cloud fraction, and cloud optical depth. This comprehensive set of kernels provides a convenient way to consistently and accurately assess the aerosol DRE uncertainties that result from observational or model-based uncertainties. The aerosol DRE kernels were used to test the effect of simplifying the full vertical profile of aerosol scattering properties into column-integrated quantities. This analysis showed that, although the clear-sky aerosol DRE can be had fairly accurately, more significant errors occur for the all-sky DRE. The sensitivity in determining the broadband spectral dependencies of the aerosol scattering properties directly from a limited set of wavelengths was quantified. These spectral dependencies can be reasonably constrained using column-integrated aerosol scattering properties in the midvisible and near-infrared wavelengths. Separating the aerosol DRE and its kernels by scene type shows that accurate aerosol properties in the clear sky are the most crucial component of the global aerosol DRE. In cloudy skies, determining aerosol properties in the presence of optically thin cloud is more radiatively important than doing so when optically thick cloud is present.

Satellite observations from the Clouds and the Earth’s Radiant Energy System show that Earth’s energy imbalance has doubled from 0.5 ± 0.2 Wm⁻² during the first 10 years of this century to 1.0 ± 0.2 Wm⁻² during the past decade. The increase is the result of a 0.9 ± 0.3 Wm⁻² increase absorbed solar radiation (ASR) that is partially offset by a 0.4 ± 0.25 Wm⁻² increase in outgoing longwave radiation (OLR). Despite marked differences in ASR and OLR trends during the hiatus (2000–2010), transition-to-El Niño (2010–2016) and post-El Niño (2016–2022) periods, trends in net top-of-atmosphere flux (NET) remain within 0.1 Wm⁻² per decade of one another, implying a steady acceleration of climate warming. Northern and southern hemisphere trends in NET are consistent to 0.06 ± 0.31 Wm⁻² per decade due to a compensation between weak ASR and OLR hemispheric trend differences of opposite sign. We find that large decreases in stratocumulus and middle clouds over the sub-tropics and decreases in low and middle clouds at mid-latitudes are the primary reasons for increasing ASR trends in the northern hemisphere (NH). These changes are especially large over the eastern and northern Pacific Ocean, and coincide with large increases in sea-surface temperature (SST). The decrease in cloud fraction and higher SSTs over the NH sub-tropics lead to a significant increase in OLR from cloud-free regions, which partially compensate for the NH ASR increase. Decreases in middle cloud reflection and a weaker reduction in low-cloud reflection account for the increase in ASR in the southern hemisphere, while OLR changes are weak. Changes in cloud cover in response to SST increases imply a feedback to climate change yet a contribution from radiative forcing or internal variability cannot be ruled out.