CERES Data Management Working Group Report
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The present study evaluates the ability of a cloud-resolving model (CRM) to simulate the physical properties of tropical deep convective cloud objects identified from a Clouds and the Earth s Radiant Energy System (CERES) data product. The emphasis of this study is the comparisons among the small-, medium- and large-size categories of cloud objects observed during March 1998 and between the large-size categories of cloud objects observed during March 1998 (strong El Ni o) and March 2000 (weak La Ni a). Results from the CRM simulations are analyzed in a way that is consistent with the CERES retrieval algorithm and they are averaged to match the scale of the CERES satellite footprints. Cloud physical properties are analyzed in terms of their summary histograms for each category. It is found that there is a general agreement in the overall shapes of all cloud physical properties between the simulated and observed distributions. Each cloud physical property produced by the CRM also exhibits different degrees of disagreement with observations over different ranges of the property. The simulated cloud tops are generally too high and cloud top temperatures are too low except for the large-size category of March 1998. The probability densities of the simulated top-of-the-atmosphere (TOA) albedos for all four categories are underestimated for high albedos, while those of cloud optical depth are overestimated at its lowest bin. These disagreements are mainly related to uncertainties in the cloud microphysics parameterization and inputs such as cloud ice effective size to the radiation calculation. Summary histograms of cloud optical depth and TOA albedo from the CRM simulations of the large-size category of cloud objects do not differ significantly between the March 1998 and 2000 periods, consistent with the CERES observations. However, the CRM is unable to reproduce the significant differences in the observed cloud top height while it overestimates the differences in the observed outgoing longwave radiation and cloud top temperature between the two periods. Comparisons between the CRM results and the observations for most parameters in March 1998 consistently show that both the simulations and observations have larger differences between the large- and small-size categories than between the large- and medium-size, or between the medium- and small-size categories. However, the simulated cloud properties do not change as much with size as observed. These disagreements are likely related to the spatial averaging of the forcing data and the mismatch in time and in space between the numerical weather prediction model from which the forcing data are produced and the CERES observed cloud systems.
Monthly mean Global Horizontal Irradiances (GHI) and Diffuse Horizontal Irradiances (DHI) are more widely available than monthly mean Direct Normal Irradiances (DNI) and Global Tilted Irradiances (GTI). Empirical methods have been developed to derive monthly mean DNI and GTI from monthly mean GHI or from GHI and DHI. In this paper, we evaluate two such methods. The first one was the Whitlock Method developed by Charles H. Whitlock (2005) for the NASA POWER project by means of regression of the BSRN data. The method expresses the monthly mean DHI-to-GHI ratio as polynomial functions of monthly mean clearness index, sunset hour angle and noon solar elevation angle on the monthly-average-day. The monthly mean DNI is calculated by dividing the monthly mean GHI-DHI difference, or DirHI, by the cosine of the solar zenith angle at the mid-time between sunrise and solar noon on the monthly-average-day. The second method is the LJCR Method developed by Liu and Jorden (1960) and Collares-Pereira and Rabl (1979), and this method empirically splits monthly mean GHI and DHI into hourly means on the monthly-average-day, and the resulting hourly mean GHI and DHI and their difference, DirHI, can then be used to compute the monthly mean DNI, GTI and the global solar tracker irradiance (GTrI). This method is also used by RETScreen. We recently produced a set of hourly DNI and DHI by bias-correcting the CERES hourly DNI and DHI, and computed hourly GTI and GTrI as well. The data span twenty plus years from March 2000 to near present on a 1 by 1 grid system. The monthly mean CERES GHI and the corrected DHI are used as inputs to the above two methods to compute monthly mean DNI, GTI and GTrI. Through comparisons with the BSRN data, it is found that the Whitlock Method, with slight modification, and the LJCR Method can produce results that are nearly as good as the results derived from the CERES hourly data.
The NASA CERES project has provided the climate quality observed TOA and computed surface fluxes to the scientific community. CERES instruments are onboard the Terra, Aqua, NPP, and NOAA-20 spacecraft. CERES uses MODIS and VIIRS imagers to retrieve cloud properties needed to convert CERES footprint radiance observations into fluxes using empirically derived angular directional models obtained during early CERES record. CERES utilizes geostationary sensors to infer the regional diurnal fluxes in between Terra and Aqua CERES observations. The imagers are also used to radiometrically scale the geostationary sensors (GEO) radiances to the imager calibration reference to ensure that the GEO derived cloud properties and broadband TOA fluxes are consistent in both space and time. Both the imager and GEO retrieved cloud properties are used to compute the surface fluxes. The Aqua-MODIS, SNPP-VIIRS, and NOAA-20 VIIRS afternoon imagers will also need to be radiometrically scaled to the same common calibration reference. Although all three imagers employ onboard solar diffusers, the calibration is not consistent over time due to instrument anomalies and ageing of the onboard calibrators. The CERES imager and GEO calibration group (IGCG) has been tasked to remove the imager channel calibration drifts for the next CERES reprocessing effort. The team will rely primarily on deep convective clouds, desert, and polar ice invariant targets to monitor the imager channel stability. The team uses ray-matched radiance pairs to radiometrically scale the SNPP and NOAA-20 VIIRS sensors to Aqua-MODIS. The radiometric scaling is further validated using geostationary imagers as transfer radiometers. The DSCOVR satellite was launched on February 25, 2015 and orbits around the L1 Lagrange point directly between the Earth and the sun. The EPIC instrument onboard DSCOVR employs a CCD array to image the Earth approximately every 2-hours. The EPIC sensor contains no onboard calibration systems. However, multiple inter-calibration studies have shown that the EPIC imager is very stable in time. This is likely due to the DSCOVR satellite being located about a 1.5M km from the Earth, where very little Earth reflected solar radiation degrades the optics. The excellent radiometric stability of EPIC allows the CERES IGCG to utilize the EPIC observations as a stable reference for monitoring the calibration stability of the three afternoon imagers, as well as to validate the radiometric scaling factors between them. Examples of the imager relative calibration using EPIC before and after radiometric scaling will be shown along with the results from the use of invariant targets to remove imager calibration drifts.
Cloud properties are critical for understanding the Earth’s radiation budget and cloud feedbacks. At NASA Langley Research Center, the Satellite ClOud and Radiative Property retrieval System (SatCORPS) provides real-time and historical analyses of clouds derived from Geostationary satellite (GEOsat) data for weather and climate applications. For the Clouds and the Earth’s Radiant Energy System (CERES) program, the global constellation of GEOsats has been analyzed since 2000 to help characterize and account for the diurnal cycle of clouds and their radiative impacts in the CERES climate data record. Obtaining consistent cloud properties over the GEOsat data record during the CERES era is a major objective but a significant challenge considering the diversity of imaging capabilities deployed during that time. The GEOsat data analysis approach for the current CERES Edition-4 (Ed4) data products was focused on accuracy and consistency with MODIS by employing as much spectral information as possible from each satellite. However, the inconsistent use of spectral information across GEOsats led to marked discontinuities in the spatial and temporal record of cloud properties that had to be accounted for post facto in downstream CERES processing. This paper reports progress in developing a new GEOsat analysis system for the next CERES edition (Ed5) that has potential to improve cross-platform consistency and continuity. In this approach, the spectral channel complement is limited to just 3-channels during daytime, ~0.65 µm (VIS), ~3.9 µm (NIR), and ~10.8 µm (IR), common to nearly all of the satellites in the record. At night, a 2-channel approach is taken with the NIR and IR, and ~6.7 µm bands that includes a machine learning approach for optically thick cloud properties. A tradeoff is the potential for reduced accuracy particularly using data from the more advanced satellites that have more spectral channels (e.g. SEVIRI, AHI and ABI) that are known to help improve thin cirrus detection, cloud-aerosol discrimination and estimates in other difficult conditions that challenge cloud remote sensing. The new continuity approach is applied to one month of global GEOSat data for each year of the CERES record since 2000 and compared with the Ed4 GEO and MODIS cloud property time series in order to evaluate the level of improved consistency in the GEOsat record and to assess the accuracy impacts. Cloud fraction will also be assessed with CALIPSO data. Outstanding issues and challenges will be discussed. The results are expected to guide future work needed to develop a more robust GEOsat cloud data record for CERES.
The NASA Clouds and the Earth's Radiant Energy System (CERES) project provides over 24 years of accurately observed top-of-the-atmosphere (TOA) and surface flux data record for climate monitoring and diagnostic studies. The major level-3 data products include SSF1deg, SYN1deg, FluxByCldTyp, and EBAF data. The CERES SYN1deg product provides 1° gridded observed TOA fluxes and computed surface fluxes based on hourly GEO clouds and derived broadband TOA fluxes. The hourly GEO derived fluxes fill in the temporal gaps between the Terra and Aqua CERES observations. In order to maintain a climate quality dataset the GEO cloud and flux retrievals are consistent over the record. The CERES project is preparing for Edition 5 product providing an opportunity to improve the GEO flux retrieval algorithms. To achieve this, the GEO satellite channel radiances are calibrated against the Aqua-MODIS calibration reference. For SW the GEO channel radiances are then converted to broadband radiances using empirical and theoretical models. The SW broadband radiances are then converted to fluxes using the same CERES angular directional models. For LW the GEO channel radiances are directly converted to LW fluxes utilizing empirical models. The resulting GEO derived broadband fluxes are then normalized to the CERES fluxes regionally by regressing the instantaneous coincident flux pairs. The normalized GEO fluxes are then tied to the CERES instrument calibration, which is very stable over time. This study will attempt to improve the GEO derived LW fluxes based on machine learning technique for the future Edition5 SYN1deg data product.
It has been well known that convection tends to be more intense over land than over ocean and continental convection generally contains wider cores that are protected from entrainment than their oceanic counterparts. How does this difference in convective intensity impact the rest of cloud types and associated cloud radiative effects (CREs)? We select five regions in the tropical “chimney” zones, with two regions over land (Africa and Amazon) and three regions over ocean (eastern and western Pacific, and Atlantic) to understand the differences in the cloud-type mean cloud properties and associated CREs using 20-years high-resolution (2x2 km2) CERES satellite data. The cloud types are based upon the ISCCP classification using the joint cloud-top pressure and cloud optical depth distribution. These five regions are compared to the entire tropics (25°S to 25°N). First, we compare the frequencies of occurrence of cloud types among the regions. Compared to the entire tropics, the “chimney” zones are cloudier, particularly, for cloud types with moderate cloud optical depths in the middle and upper troposphere. Except for the eastern Pacific region, the low-level cloud types are less abundant over land (due to higher boundary layers) and ocean (due to weaker lower-tropospheric subsidence). Over Africa, upper-level clouds are scarcer, due perhaps to weak production of convective anvils related to the drier atmosphere. Second, liquid water path (LWP) and ice water path (IWP) are compared among the regions. Compared to the entire tropics, mean LWPs for moderate cloud optical depths are higher for oceanic regions but lower for Africa while they are slightly lower in the lower troposphere but slightly higher in the middle troposphere for Amazon. Not surprisingly, mean IWPs are higher over land than over ocean, which is directly related to the difference in convective intensity. That is, the explanations to this result are the relatively weak convective intensity over the eastern Pacific and the abundance of anvil clouds over the western Pacific. Both reduce the mean IWPs. The CRE results are still being analyzed because the totally-clear grids (1° x 1°) of the CERES data product should be removed before taking a regional mean of clear-sky fluxes, which cause the wrong signs of CREs for a few cloud types.
The NASA Clouds and the Earth’s Radiant Energy System (CERES) SYN1deg(Ed4.1) provides hourly global horizontal irradiance (GHI), diffuse horizontal irradiance (DHI) and direct horizontal irradiance (DirHI) at 1-degree latitude by 1-degree longitude resolution, and the time span is from March 2000 to near present. While the GHI agrees well with the ground-based Baseline Surface Radiation Network (BSRN) data, the DHI and DirHI are, respectively, positively and negatively biased against the BSRN data, and so is the direct normal irradiance (DNI) derived simply by dividing the DirHI by the cosine of the hourly solar zenith angle (SZA). We found that the biases tend to be well-defined functions of cos(SZA) and the cloud fraction (CLFR). We thus performed a bias-correction on the hourly DHI and DNI in the latitude-cos(SZA)-CLFR phase space, and therefrom, we derived the hourly global tilted irradiance (GTI) at a number of tilt angles using the isotropic diffuse irradiance model, and the results agree well with the GTI derived from the original BSRN near-instantaneous records. Meanwhile, the National Renewable Energy Laboratory (NREL) has produced the National Solar Radiation Database (NSRDB) with finer spatiotemporal resolution using its Physical Solar Model (PSM) and surface-based measured data, model results and satellite-based data as inputs. Now the NSRDB data covers not only the United States and its neighboring regions, but Europe, Africa, Asia and Oceania. The available temporal resolutions are 5, 10, 15, 30 and 60 minutes and the spatial resolutions are 2, 4 and 10 km, depending on the selected region. In this presentation, we compare both the CERES-based and NSRDB GHI, DHI and DNI with their BSRN counterparts and show how finer spatial resolution may give us an advantage in data accuracy and usability.
The NASA Clouds and the Earth's Radiant Energy System (CERES) product provides over 20 years of accurately observed top-of-the-atmosphere (TOA) and surface flux data record for climate monitoring and diagnostic studies. The interaction between clouds and radiation interaction is a key factor that dominate climate feedbacks but is not well understood. To further advance our understanding of the cloud-radiation interaction, a new CERES FluxByCldTyp (FBCT) product has been developed that contains radiative fluxes by cloud-type, which can provide more stringent constraints when validating models and reveal more insight into the interactions between clouds and climate. For CERES partly cloudy and multiple cloud-type footprints, the FBCT product utilizes Moderate Resolution Imaging Spectroradiometer (MODIS) narrow-band (NB) imager channel radiances partitioned by cloud-type within a CERES footprint to estimate the cloud-type broadband fluxes. The MODIS multi-channel derived broadband fluxes were compared with the CERES observed footprint fluxes and were found to be within 1% and 2.5% for LW and SW, respectively, as well as being mostly free of cloud property dependencies. The FBCT all-sky and clear-sky monthly averaged fluxes were found to be consistent with the CERES SSF1deg product. This study takes advantage of recent progress in machine learning (ML) field by applying deep neural network algorithm to improve fluxes based on MODIS NB radiances. The preliminary study shows ML produce are an improvement over the current FBCT Edition 4 NB2BB algorithm. Furthermore, unlike Ed4 NB2BB, the new ML method convert NB radiances directly to broadband fluxes. For future Ed5, new NB radiances are proposed and used by ML to improve fluxes calculation. Prelimary results show significant LW improvement.
The wide field of view instrument on Earth Radiation Budget Experiment (ERBE) Earth Radiation Budget Satellite (ERBS) took measurements from November 1984 through September 1999. It provides top-of-atmosphere (TOA) reflected shortwave and emitted longwave irradiance at a 10 degree by 10 degree spatial resolution. We have recently re-processed ERBS wide field of view data with algorithms similar to those used in the CERES process (Shrestha et al. 2019). We also calibrated ERBS shortwave and longwave irradiance with CERES-derived irradiances. Combining ERBS and CERES data, we analyze the time series of global mean net TOA irradiance from 1985. The change of TOA net radiation over 30 years derived from ERBS and CERES data is larger than the change derived by Allan et al. (2014). The difference in the net irradiance change from two data sets is, however, within the uncertainty of ocean heating rates derived from ocean temperature measurements in 80?s. We also investigate the spatial distribution of changes in reflected shortwave irradiances derived from ERBS and CERES periods. The spatial distribution of reflected shortwave irradiance changes over the tropics and sub-tropical regions derived from our data set is qualitatively consistent with the result of Norris et al. (2016).
The Clouds and Earth’s Radiant Energy System (CERES) project provides satellite-based observations of how Earth’s energy flows are varying in time and space and how clouds and aerosols are affecting the Earth’s radiation budget. Nominally, CERES data products require months of validation and calibration before releasing a climate quality data. The Fast Longwave and Shortwave Radiative Flux (FLASHFlux) data product was developed as to provide key radiative flux data products within a week of observations for scientific, applied science research and educational usage. FLASHFlux achieves this by using simplified calibration process and optimized CERES data fusion production system, an operation meteorology product from Global Modeling and Assimilation Office(GMAO), and surface parameterized radiative flux models. The CERES FLASHFlux provides two data products: 1) an hourly Level 2 Single Scanner Footprint (SSF) data separately for Aqua and Terra observations, and 2) a daily Level 3 Time Interpolated and Spatially Averaged (TISA) 1ox1o gridded data that combines Aqua and Terra observations.The current operational Version 3C is being transitioned to Version4A. We compare the FLASHFlux Version4A fluxes using the specialized calibration to Version3C and to the climate quality CERES Edition4A TOA and surface fluxes over time for Terra and Aqua to determine the uncertainties of the calibration process. In addition, we also compare the impact of the fluxes by using the near real-time GMAO GEOS FP-IT instead of GMAO G5-CERES meteorology data.We also compare both SSF overpass products (for Terra and Aqua) and TISA radiative flux products to a set of surface measurement sites that are globally distributed and assess agreement.Lastly, we introduce the first CERES FF SSF data products from the CERES instrument on board NOAA-20.
The NASA Clouds and the Earth's Radiant Energy System (CERES) product provides over 20 years of accurately observed top-of-the-atmosphere (TOA) and surface flux data record for climate monitoring and diagnostic studies. The interaction between clouds and radiation interaction is a key factor that dominate climate feedbacks but is not well understood. To further advance our understanding of the cloud-radiation interaction, a new CERES FluxByCldTyp (FBCT) product has been developed that contains radiative fluxes by cloud-type, which can provide more stringent constraints when validating models and reveal more insight into the interactions between clouds and climate. For CERES partly cloudy and multiple cloud-type footprints, the FBCT product utilizes Moderate Resolution Imaging Spectroradiometer (MODIS) narrow-band (NB) imager channel radiances partitioned by cloud-type within a CERES footprint to estimate the cloud-type broadband fluxes. The MODIS multi-channel derived broadband fluxes were compared with the CERES observed footprint fluxes and were found to be within 1% and 2.5% for LW and SW, respectively, as well as being mostly free of cloud property dependencies. The FBCT all-sky and clear-sky monthly averaged fluxes were found to be consistent with the CERES SSF1deg product. This study takes advantage of recent progress in machine learning (ML) field by applying deep neural network algorithm to improve fluxes based on MODIS NB radiances. The preliminary study shows ML produce are an improvement over the current FBCT Edition 4 NB2BB algorithm. Furthermore, unlike Ed4 NB2BB, the new ML method convert NB radiances directly to broadband fluxes. For future Ed5, new NB radiances are proposed and used by ML to improve fluxes calculation. Preliminary results show significant LW improvement.
There are 2 CERES scanning radiometer instruments aboard the TERRA spacecraft, one for mapping the solar radiation reflected from the Earth and the outgoing longwave radiation and the other for measuring the anisotropy of the radiation. Each CERES instrument has on-board calibration devices, which have demonstrated that from ground to orbit the broadband total and shortwave sensor responses maintained their ties to the International Temperature Scale of 1990 at precisions approaching radiances have been validated in orbit to +/- 0.3 % (0.3 W/sq m sr). Top of atmosphere fluxes are produced by use of the CERES data alone. By including data from other instruments, surface radiation fluxes and radiant fluxes within the atmosphere and at its top, shortwave and longwave, for both up and down components, are derived. Validation of these data products requires ground and aircraft measurements of fluxes and of cloud properties.
The NASA Clouds and the Earth's Radiant Energy System (CERES) product provides over 20 years of accurately observed top-of-the-atmosphere and surface flux data for climate studies. The interaction between clouds and radiation interaction is a key factor that dominate climate feedbacks but is not well understood. To further advance our understanding of the cloud-radiation interaction, a new CERES FluxByCldTyp (FBCT) product has been developed that contains radiative fluxes by cloud-type, which can provide more stringent constraints when validating models. The FBCT product utilizes Moderate Resolution Imaging Spectroradiometer (MODIS) narrow-band (NB) imager channel radiances partitioned by cloud-type within a CERES footprint to estimate their broadband fluxes. The MODIS multi-channel derived broadband fluxes were compared with the CERES observed footprint fluxes and were found to be within 1% and 2.5% for LW and SW, respectively, as well as being mostly free of cloud property dependencies. The FBCT all-sky and clear-sky monthly averaged fluxes were found to be consistent with the CERES SSF1deg product. This study takes advantage of recent progress in machine learning (ML) field by applying deep neural network algorithm to improve fluxes based on MODIS NB radiances. The preliminary study shows ML produce are an improvement over the current FBCT Edition 4 NB2BB algorithm. Furthermore, unlike Ed4 NB2BB, the new ML method convert NB radiances directly to broadband fluxes. For future Ed5, new NB radiances are proposed and used by ML to improve fluxes calculation. Preliminary results show significant LW improvement.
The NASA Clouds and the Earth's Radiant Energy System (CERES) product provides over 20 years of accurately observed top-of-the-atmosphere (TOA) and surface flux data record for climate monitoring and diagnostic studies. The interaction between clouds and radiation interaction is a key factor that dominate climate feedbacks but is not well understood. To further advance our understanding of the cloud-radiation interaction, a new CERES FluxByCldTyp (FBCT) product has been developed that contains radiative fluxes by cloud-type, which can provide more stringent constraints when validating models and reveal more insight into the interactions between clouds and climate. For CERES partly cloudy and multiple cloud-type footprints, the FBCT product utilizes Moderate Resolution Imaging Spectroradiometer (MODIS) narrow-band (NB) imager channel radiances partitioned by cloud-type within a CERES footprint to estimate the cloud-type broadband fluxes. The MODIS multi-channel derived broadband fluxes were compared with the CERES observed footprint fluxes and were found to be within 1% and 2.5% for LW and SW, respectively, as well as being mostly free of cloud property dependencies. The FBCT all-sky and clear-sky monthly averaged fluxes were found to be consistent with the CERES SSF1deg product. This study takes advantage of recent progress in machine learning (ML) field by applying deep neural network algorithm to improve fluxes based on MODIS NB radiances. The preliminary study shows ML produce are an improvement over the current FBCT Edition 4 NB2BB algorithm. Furthermore, unlike Ed4 NB2BB, the new ML method convert NB radiances directly to broadband fluxes. For future Ed5, new NB radiances are proposed and used by ML to improve fluxes calculation. Prelimary results show significant LW improvement.
Three boundary-layer cloud object types, stratus, stratocumulus and cumulus, that occurred over the Pacific Ocean during January-August 1998, are identified from the CERES (Clouds and the Earth s Radiant Energy System) single scanner footprint (SSF) data from the TRMM (Tropical Rainfall Measuring Mission) satellite. This study emphasizes the differences and similarities in the characteristics of each cloud-object type between the tropical and subtropical regions and among different size categories and among small geographic areas. Both the frequencies of occurrence and statistical distributions of cloud physical properties are analyzed. In terms of frequencies of occurrence, stratocumulus clouds dominate the entire boundary layer cloud population in all regions and among all size categories. Stratus clouds are more prevalent in the subtropics and near the coastal regions, while cumulus clouds are relatively prevalent over open ocean and the equatorial regions, particularly, within the small size categories. The largest size category of stratus cloud objects occurs more frequently in the subtropics than in the tropics and has much larger average size than its cumulus and stratocumulus counterparts. Each of the three cloud object types exhibits small differences in statistical distributions of cloud optical depth, liquid water path, TOA albedo and perhaps cloud-top height, but large differences in those of cloud-top temperature and OLR between the tropics and subtropics. Differences in the sea surface temperature (SST) distributions between the tropics and subtropics influence some of the cloud macrophysical properties, but cloud microphysical properties and albedo for each cloud object type are likely determined by (local) boundary-layer dynamics and structures. Systematic variations of cloud optical depth, TOA albedo, cloud-top height, OLR and SST with cloud object sizes are pronounced for the stratocumulus and stratus types, which are related to systematic variations of the strength of inversion with cloud object sizes, produced by large-scale subsidence. The differences in cloud macrophysical properties over small regions are significantly larger than those of cloud microphysical properties and TOA albedo, suggesting a greater control of (local) large-scale dynamics and other factors on cloud object properties. When the three cloud object types are combined, the relative population among the three types is the most important factor for determining the cloud object properties in a Pacific transect where the transition of boundary-layer cloud types takes place.
Earth’s radiation budget (ERB) is a key indicator of the climate system as it describes the delicate balance between how much of the sun’s energy reaching Earth is absorbed and how much thermal infrared radiation is emitted back to space. It is now well established that there is a positive net imbalance in ERB of 0.5-1 Wm-2. This imbalance, referred to as Earth’s Energy Imbalance (EEI), is shaped by natural and anthropogenic climate forcings and the climate system’s response to the forcing.When averaged over a long period, a positive EEI leads to an increase in surface temperature, sea levelrise, ocean heat content, melting of snow and sea-ice, and an acceleration of the hydrological cycle.EEI is also influenced by internal climate variability at interannual (e.g., El Niño/Southern Oscillation(ENSO)) and longer timescales (e.g., Pacific Decadal Oscillation (PDO)). Here we use CERES data toshow that since 2000, EEI has been increasing at a rate of 0.04+-0.03 Wm-2 per yr, corresponding to 0.8 Wm-1 over the 20-year CERES record. To understand the underlying causes for this increase, we analyze the CERES data in the context of internal climate variability indices and regional pattern changes. We also perform a partial radiative perturbation analysis that decomposes the TOA SW and LW radiation variations in terms of contributions by clouds, aerosols, surface albedo, temperature,water vapor and other gases (e.g., ozone, carbon dioxide, methane, etc.).