Uncertainties in CERES Top-of-Atmosphere Fluxes Caused by Changes in Accompanying Imager

: The Clouds and the Earth’s Radiant Energy System (CERES) project provides observations of Earth’s radiation budget using measurements from CERES instruments on board the Terra, Aqua, Suomi National Polar-orbiting Partnership (S-NPP), and NOAA-20 satellites. The CERES top-of-atmosphere (TOA) ﬂuxes are produced by converting radiance measurements using empirical angular distribution models, which are functions of cloud properties that are retrieved from imagers ﬂying with the CERES instruments. As the objective is to create a long-term climate data record, not only calibration consistency of the six CERES instruments needs to be maintained for the entire time period, it is also important to maintain the consistency of other input data sets used to produce this climate data record. In this paper, we address aspects that could potentially affect the CERES TOA ﬂux data quality. Discontinuities in imager calibration can affect cloud retrieval which can lead to erroneous ﬂux trends. When imposing an artiﬁcial 0.6 per decade decreasing trend to cloud optical depth, which is similar to the trend difference between CERES Edition 2 and Edition 4 cloud retrievals, the decadal SW ﬂux trend changed from − 0.35 ± 0.18 Wm − 2 to 0.61 ± 0.18 Wm − 2 . This indicates that a 13% change in cloud optical depth results in about 1% change in the SW ﬂux. Furthermore, different CERES instruments provide valid ﬂuxes at different viewing zenith angle ranges, and including ﬂuxes derived at the most oblique angels unique to S-NPP ( > 66 ◦ ) can lead to differences of 0.8 Wm − 2 and 0.3 Wm − 2 in global monthly mean instantaneous SW ﬂux and LW ﬂux. To ensure continuity, the viewing zenith angle ranges common to all CERES instruments (<66 ◦ ) are used to produce the long-term Earth’s radiation budget climate data record. The consistency of cloud properties retrieved from different imagers also needs to be maintained to ensure the TOA ﬂux consistency.


Introduction
Earth's radiation budget (ERB) at the top-of-atmosphere (TOA) is one of the most critical variables that shows the combined effects of many processes within the Earth-atmosphere system. Currently, the Clouds and the Earth's Radiant Energy System (CERES) instruments [1,2] have been providing continuous global TOA reflected shortwave (SW) radiation and outgoing longwave (LW) radiation (OLR) since 2000. Seven CERES instruments have flown on five different satellites thus far. Table 1 listed the CERES flight mode (FM) on different satellites, imagers that they fly alongside with, their orbit altitude, and their nominal nadir footprint size. The CERES instrument consists of a three-channel broadband scanning radiometer [1]. The scanning radiometer measures radiances in shortwave (SW, 0.3-5 µm) channel, total (0.3-200 µm) channel, and window (WN, 8-12 µm) channel for FM1-FM5 and longwave (LW, 5-35 µm) channel for FM6. For FM1-FM5, LW radiance is derived as the difference between total and SW channels; for FM6, the LW channel measurement offers an independent check to the LW radiance derived by differencing total and SW channels. These measured radiances (I) at a given sun-Earth-satellite geometry are converted to outgoing reflected solar and emitted thermal TOA radiative fluxes (F) as: where θ 0 is the solar zenith angle, θ is the CERES viewing zenith angle, φ is the relative azimuth angle between CERES and the solar plane, and R j (θ 0 , θ, φ) is the anisotropic factor for scene type j. Here, scene type is a combination of variables (e.g., surface type, cloud fraction, cloud optical depth, cloud phase, aerosol optical depth, precipitable water, lapse rate, etc.) that are used to group the data to develop distinct angular distribution models (ADMs). The scene type classifications are made possible by the imagers that are flying alongside the CERES instruments. Note the SW ADMs are developed as a function of θ 0 , θ, φ for each scene type, whereas the LW ADMs are a weak function of θ 0 and φ and are developed only as a function of θ [3,4].
To create a long-term climate data record of the ERB using CERES observations, consistency in instrument calibration and algorithm implementation must be maintained throughout the entire time series. Loeb et al. [2] detailed how the gains and spectral response functions are determined for CERES instruments to ensure radiance consistency over time and across different instruments. The CERES measured radiances are converted to fluxes using scene-type dependent anisotropic factors (Equation (1)). These anisotropic factors are determined using the CERES rotating azimuth plane data collected at the beginning of the Terra and Aqua missions and the cloud properties retrieved from MODIS spectral measurements. It is evident that any drift in CERES measured radiances will result in drift in fluxes, but how will drift in imager spectral radiances affect the flux inversion by changing the anisotropic factor selection is unknown. In this paper, we will quantify the impact of imager radiance drift on flux inversion. Furthermore, the maximum viewing zenith angle (VZA) of S-NPP CERES footprints with VIIRS cloud retrievals is about 71.4 • , whereas it is only about 66 • for CERES on Terra and Aqua. The impact of the VZA ranges on cloud properties and fluxes is discussed in Section 4. As shown in Table 1, the imagers on Aqua and S-NPP are different, and VIIRS does not have any infrared channels sensitive to water vapor and carbon dioxide. Therefore the MODIS cloud retrieval algorithms need to be modified to apply to S-NPP observations. The impact on cloud property retrieved from MODIS and VIIRS is detailed in Section 5. Discussions and conclusions are presented in Section 6.

CERES ADMs
To facilitate the construction of ADMs, there are pairs of identical CERES instruments on both Terra and Aqua (see Table 1). At the beginning of these missions one of the instruments on each satellite was always placed in a rotating azimuth plane (RAP) scan mode, while the other one was placed in cross-track mode to provide spatial coverage. When in RAP mode, the instrument scans in elevation as it rotates in azimuth, thus acquiring radiance measurements from a wide range of viewing combinations. There are about 60 months of RAP data collected on Terra and about 32 months of RAP data collected on Aqua.
The spatial resolution of the CERES instruments on Terra and Aqua is about 20 km at nadir, and they are flying with MODIS instruments. The spatial resolution of CERES instruments on S-NPP and NOAA-20 is about 24 km at nadir, and they are flying with VIIRS instruments. The higher-resolution imagers provide cloud retrievals for every CERES field of view (FOV). The cloud algorithms developed by the CERES cloud working group retrieve cloud fraction ( f , in %), cloud optical depth, cloud phase (liquid = 1, and ice = 2), cloud top and effective temperature/pressure (among other variables) based on imager pixel-level measurements [5][6][7]. These pixel-level cloud properties are spatially and temporally matched with the CERES FOV. If the cloud top pressure within a CERES FOV is significantly different, up to two non-overlapping cloud layers are defined [8]. For each cloud layer, pixel-level cloud properties are weighted by the CERES point spread function (PSF, [9]) to provide layer-mean cloud properties. Similarly, spectral radiances from MODIS/VIIRS observations are averaged over the CERES footprints weighted by the CERES PSF. Layer-mean cloud properties are then weighted by the layer cloud fractions to derive the cloud properties over the CERES FOV. These cloud properties are used for scene type classification. Cloud detection and cloud property retrievals are assessed through comparison with CALIPSO and CloudSat [7,10]. Instantaneous comparison of nadir daytime cloud detection agrees with CALIPSO in 90-99% of the cases over snow-free areas, and in 89-94% of the cases over snow-covered areas. With a comparison using only single-phase CALIPSO pixels, the cloud top phase from MODIS retrievals agree with CALIPSO 95.8% and 92.0% of the time over snow-free and snow/ice-covered areas during the daytime. Comparisons of other cloud properties can be found in [10].
Besides cloud properties, surface types, and the snow and ice percentage cover are also used for scene type classification. The International Geosphere Biosphere Program ( [11], IGBP) global land cover data set and a combination of the National Snow and Ice Data Center (NSIDC) microwave snow/ice map and the National Environmental Satellite, Data and Information Service (NESDIS) snow/ice map are utilized for this purpose. The spatial resolution of IGBP is 10 min and up to eight surface types are kept for each CERES FOV. The spatial resolutions for both snow/ice maps are also 10 min. These two maps are complementary to each other. NESDIS uses imager data to identify snow and sea ice and provide snow and sea ice information near the coast, whereas NSIDC does not provide microwave retrievals within 50 km of the coast.
The Edition 4 CERES ADMs [4] are constructed for various scene types. Over a given scene type (j), the CERES measured radiances are sorted into discrete angular bins. Averaged radiances (Î) in all angular bins are calculated and all radiances in the upwelling directions are integrated to provide the ADM flux (F). The ADM anisotropic factors (R) for scene type j are then calculated as: Over cloudy ocean and land scenes, the SW ADMs were constructed using a sigmoidal function that relates CERES radiances and cloud properties retrieved from MODIS [4]. Under cloudy conditions, the magnitude of CERES radiance in a given angular bin is most sensitive to cloud fraction ( f , in %), cloud optical depth (τ), and effective cloud phase (ECP), which is a PSF-weighted average of the cloud phase derived from imager pixel data (1 for liquid clouds, and 2 for ice clouds). The analytical fits were derived as a function of ln( f τ) separately for liquid (1.00 < ECP < 1.01), mixed (1.01 ≤ ECP ≤ 1.75), and ice (1.75 < ECP ≤ 2.00) clouds. For each CERES radiance measurement (I o ), the scene type is determined based on the cloud properties retrieved from the imager and other ancillary datasets to select the appropriate anisotropic factor to convert the radiance to flux using Equation (1).

CERES Energy Balanced and Filled (EBAF) Product
CERES Level-3 data products provide monthly TOA flux at 1 • latitude by 1 • longitude spatial resolution, and it is the only global dataset that can be used to study the variations of ERB over a range of time-space scales. However, the absolute accuracy requirement necessary to quantify the Earth's energy imbalance (EEI) is much higher than any current instruments can deliver. The EEI ranges between 0.5 and 1 Wm −2 [12]. This is a small residual of TOA flux of the order of 340 Wm −2 , about 0.15% of the total incoming and outgoing radiation at the TOA. Given the current uncertainty in instrument calibration and algorithm implementation, the net imbalance is about 4.3 Wm −2 using the most recent CERES Edition 4 data product [13]. The CERES EBAF dataset uses an objective constrainment algorithm to adjust SW and LW TOA fluxes within their ranges of uncertainty to remove the inconsistency between average global net TOA flux and heat storage in the Earth-atmosphere system [13,14]. Additionally, because of the relatively coarse spatial resolution of the CERES instruments (20-km at nadir), the standard CERES Level-3 data products have many spatial gaps in monthly mean clear-sky TOA flux maps due to the absence of cloud-free areas occurring at the CERES footprint scale in some 1 • ×1 • regions. In EBAF, this problem is mitigated by inferring clear-sky fluxes from both CERES and MODIS measurements to produce a new clear-sky TOA flux climatology for every 1 • ×1 • grid cell every month.

CERES TOA SW Flux Trend
The CERES Ed4.0 EBAF data set represents the best effort in producing climate data record to achieve consistency throughout the entire time series for assessing climate variability (e.g., [15][16][17]), and for evaluating climate models (e.g., [18][19][20]). When the TOA reflected SW flux trend from EBAF Ed4.0 is compared with that from EBAF Ed2.8, we note that the decreasing trend in Ed4.0 is about four times larger than that from Ed2.8 ( Figure 1). On the contrary, the trends of the daytime TOA LW flux from these two editions are very similar. To understand if this difference is caused by the CERES instrument calibration difference between these two editions, we compared the TOA flux anomaly produced using the same cloud retrievals and ADMs (those used to produce Ed2.8 EBAF) but different calibration coefficients used in Ed2.8 and Ed4.0 EBAF. The resulting anomalies are nearly identical with the decadal trends agree within their 95% uncertainty range (not shown), indicating that the CERES calibration differences are not the cause for the differences seen in Figure 1.   [21]. The abrupt change in 2009 was caused by the solar diffuser degradation, which was about 1.5% for Terra-MODIS [22]. These degradations in the Terra MODIS visible channels were not corrected in Collection 4 and Collection 5 data that were used for CERES Ed2 cloud retrievals, but have been corrected in Collection 6. Therefore the Ed4 cloud retrievals are free from this artificial trend.
MODIS cloud optical depths are used to select anisotropic factors in converting CERES radiances to radiative fluxes, the different optical depth trends shown in Figure 2 may have contributed to the different SW TOA flux trends as well. To quantify how artificial cloud optical depth trend impacts the TOA SW flux trend, we impose an Ed2-like cloud optical depth trend to the Ed4 retrievals. The global mean cloud optical depth trends differ by 0.6 per decade between these two editions. To force the Ed4 cloud optical depth to have a trend similar to that of Ed2, cloud optical depths of all cloudy footprints in Ed4 are decreased by 0.005 per month and are set to be 0.03 when negative cloud optical depths are encountered. The CERES TOA fluxes inverted by using ADMs selected based on Ed4 cloud properties are compared with the TOA fluxes inverted by using the same ADMs and cloud properties except that the cloud optical depths are imposed with an Ed2-like trend. All simulations are based on the CERES monthly mean Single Scanner Footprint TOA/Surface Fluxes and Cloud (SSF) data processing for daytime.   Figure 4 shows the SW anisotropic factors for liquid clouds over ocean solar zenith angles of 30 • and 50 • for a typical side scattering angle that CERES collects measurements when in cross track mode. The CERES SW flux is determined by accounting the anisotropy of the radiance field, and the anisotropic factors are selected based upon the cloud properties. When cloud optical depths decrease (while the cloud fraction stays more or less the same), the anisotropic factors also decrease for observations taken at most viewing zenith angles less than ∼55 • but increase for more oblique views. The decreasing cloud optical depth thus causes an overall increase in TOA SW flux, as the number of footprints with viewing zenith angles less than ∼55 • far outweigh those with viewing zenith angles greater than ∼55 • . This simulation shows that the TOA SW flux trend in EBAF Ed2.8 was indeed 'flattened' by the trend in cloud optical depth caused by the MODIS calibration issues. Note the trend difference shown in Figure 3a is greater than the trend difference between EBAF 2.8 and EBAF 4.0 (Figure 1a). The possible explanations for the larger TOA SW flux trend change in the imposed cloud optical depth trend simulation are: (1) the imposed trend is continuous, whereas the Ed2 cloud optical depth shows step changes; (2) the impact of the MODIS calibration changes on cloud optical depth is different for optically thin and thick clouds; (3) ADMs are more sensitive to optical depth changes for optically thin than thick clouds.  For the TOA LW flux, imposing a 0.005 per month decreasing trend to the Ed4 cloud optical depth has little impact on the trend. This is because the LW ADMs for cloudy conditions are constructed as a function of "pseudoradiance", which characterizes scenes in terms of cloud fraction, precipitable water, surface skin temperature, surface-cloud temperature difference, surface infrared emissivity, and cloud infrared emissivity [4]. The pseudoradiance has a weak dependence on visible cloud optical depth via cloud infrared emissivity and thus the imposed trend on cloud optical depth has little effect on the inverted LW flux.
This result demonstrates that to produce the ERB climate data record, not only it is paramount to maintain the CERES calibration stability, it is also important to maintain the imager calibration stability. Spurious changes in imager radiances can impact the retrieved cloud properties and thus lead to an erroneous TOA SW flux trend.

Viewing Zenith Angle Impact on Imager Retrieved Cloud Property and CERES TOA Flux
Since the MODIS imager can only scan to a maximum viewing zenith angle (VZA) of ∼66 • , this means that only CERES footprints with VZA < 66 • are retained in the SSF data when CERES is in the cross track scan mode. The VIIRS imager on S-NPP and NOAA-20, however, scans to a maximum VZA of ∼71 • and the accompanying CERES instruments thus provide fluxes for all footprints with VZA < 72 • . Cloud amount from satellite retrievals increases with the imager viewing zenith angle, and the increment is not uniform across different cloud regimes, with a more pronounced increase observed for broken cloud fields [23]. As demonstrated in Section 3, changes in cloud properties can affect the scene identification used to select the ADMs and thus affect the inverted fluxes. Figure 5 shows the multi-year (2014-2018) mean cloud fraction from S-NPP VIIRS and the cloud fraction difference between all footprints and footprints with VZA < 66 • . The differences are positive everywhere and the global mean cloud fractions derived using all footprints is about 0.5% greater than only using footprints with VZA < 66 • . Note the sharp line around 65 • N is because north of this latitude the number of observations taken at VZA > 66 • is very low. The largest differences (up to 2%) are observed over tropical oceans and Antarctica where multi-year mean cloud fractions are around 50%. Regions where cloud fractions exceed 80% and desert areas with very low cloud amount show smaller differences. The VZA dependency is also smaller off the west coast of Africa, where a special cloud mask is implemented over this dust laden region [7], than adjacent regions. Many studies [23][24][25] have noted that satellite retrieved cloud fraction increases as the VZA increases and attributed this to viewing perspective of the sensor, increase in path length, and in pixel size. Although the increase of cloud fraction at the very oblique angles (>66 • ) is obvious, the magnitude is relatively small when compared to the cloud fraction difference between near-nadir (VZA < 10 • ) views and oblique views (VZA > 60 • ). Global mean cloud fraction retrieved from near-nadir view footprints is about 8% less than using oblique view footprints (not shown). The dependency of cloud fraction on VZA also propagates into other cloud properties. Figure 6 shows the multi-year mean daytime liquid and ice cloud optical depth from S-NPP VIIRS retrievals and the difference between all footprints and footprints with VZA < 66 • . For liquid clouds, optical depth differences are mostly negative with the largest differences seen over northern high latitude land regions. Positive differences are observed over desert and high elevation regions, and over the Southern Ocean and Arctic Ocean. For ice clouds, the difference shows a large land and sea contrast. Negative differences of up to 0.4 are observed over the land regions, whereas mostly positive differences are observed over ocean except negative differences over high latitude oceanic regions. Many factors contribute to the VZA dependency of cloud optical depth retrievals. The dependency is more pronounced under oblique solar zenith angles. Studies also show that in the backward scattering directions cloud optical depth decreases with VZA whereas it increases with VZA in the forward scattering directions, and that the cloud optical depth dependency on VZA tends to be larger in the backward scattering directions than in the forward scattering directions (e.g., [26][27][28]). Thus the overall biases depend on the sampling differences in the forward and backward directions. The biases also depend on the cloud types and cloud spatial distribution pattern such that the broken clouds are more susceptible to these VZA dependency biases than uniformly-distributed clouds.
Given that cloud properties depend on VZA and that CERES empirical ADMs are functions of these cloud properties, different ranges of VZAs on different CERES instruments can affect the continuity of the ERB climate data record. Figure 7 shows the monthly mean instantaneous daytime SW flux difference averaged using all footprints and using footprints with VZA < 66 • only from CERES S-NPP observations in 2014. For all four months, the global monthly mean daytime SW flux would be about 0.7 to 0.9 Wm −2 smaller if the VZA range is not restricted. Regionally, the flux can differ by up to 10 Wm −2 . Large differences are observed around the ITCZ and around 60 • N and 60 • S. Figures 8 and 9 show the daytime and nighttime monthly mean LW flux differences averaged using all footprints and also using footprints with VZA < 66 • only from CERES S-NPP observation. For daytime, the LW flux would be about 0.1-0.3 Wm −2 smaller if all VZAs are included, and the nighttime LW flux would be about 0.1 Wm −2 smaller. Regionally, the LW flux differences of up to 4 Wm −2 are seen over the tropical oceans. Similar to SW flux, large differences are observed near the ITCZ and over the western Pacific Ocean. However, it is problematic to isolate the contributing factors that cause the flux differences as they depend on cloud type, cloud amount, cloud optical depth, and phase, as well as the relative number of observations taken with VZA > 66 degrees and ≤66 degrees.
To ensure continuity of the CERES flux data record on the global and regional scales, we decided to only include footprints with VZA < 66 • for CERES S-NPP and NOAA-20 Level 3 and Level 4 product to ensure continuity of ERB climate data record. However, footprints with VZA > 66 • are kept in the Level 2 product for completeness and for process studies.

Cloud Property Consistency between MODIS and VIIRS
Both Aqua and S-NPP are in sun-synchronous orbits, but fly at different altitudes (see Table 1). The S-NPP flies directly over Aqua in about every 64 h. The overlap provides a great opportunity to exam the consistency of cloud property retrieved from MODIS and VIIRS. These two imagers have many similar channels, but the most striking difference is the lack of water vapor (at 6.7 µm) and carbon dioxide (15 µm) channels on VIIRS. The lack of these channels makes the application of MODIS-based cloud retrieval algorithm to VIIRS problematic, and algorithm changes have to be made to accommodate these differences. Here we use matched footprints on Aqua and S-NPP to assess the cloud fraction and cloud optical depth difference.
The matching criteria used are that both the latitude and longitude differences are less than 0.05 • , the solar zenith angle and viewing zenith angle differences are less than 2 • , and the relative azimuth angle difference is less than 5 • . To minimize the impact of footprint size difference on the comparison, we only use footprints with viewing zenith angles less than 10 • . The cloud retrieval algorithms are different for polar and non-polar regions, therefore the comparisons are done separately for polar and non-polar regions and only daytime comparisons are presented here. We used Aqua Edition 4 and NPP Edition 1 datasets for this comparison. Figure 10 shows the probability distribution of the cloud fraction and cloud optical depth differences between S-NPP and Aqua separately for regions between 60 • S-60 • N and for regions poleward of 60 • . The number of matched footprints over the non-polar regions is 32,222, and most of them are over ocean; whereas the number of matched footprints over polar regions is only 5880. As expected, the polar regions have more and thicker clouds than the non-polar regions. The cloud fraction from S-NPP on average is smaller than that from Aqua by about 1%, with the root mean square (RMS) error around 7% for both the non-polar and polar regions. The cloud optical depth from S-NPP and Aqua agree fairly well with each other over non-polar regions, with the bias close to 0 and the RMS error of 2.3. However, the optical depth difference over polar regions is much larger. The cloud optical depth from S-NPP is about 0.4 greater than that from Aqua on average and the RMS error is over 9.
As mentioned early, the CERES ADMs are functions of cloud fraction and cloud optical depth. It is important to maintain cloud property consistency in order to achieve flux consistency. Algorithm development is underway to reconcile the cloud property differences.

Discussion and Conclusions
Producing a multi-decade-long data record that is suitable for detecting and quantifying the climate variability and changes is a very challenging task. Climate data records that are decades long generally involve multiple satellites and consistency of the measurements across multiple platforms is key to maintain the stability of the measurements. Additionally, for data records that use other input data sets, it is also important to maintain the consistency of the input data sets.
CERES TOA flux data set is produced by using the radiances measured by the CERES instruments, which are then converted to fluxes by using the scene-type dependent ADMs. These ADMs are functions of cloud properties that are retrieved from imagers that are flying alongside the CERES instruments. In this paper, we illustrated that a spurious discontinuity in imager calibration introduced an erroneous trend in cloud optical depth retrieval and thus affected the TOA SW flux trend. When the Ed4 cloud optical depth was decreased by 0.005 per month to mimic the Ed2 cloud optical depth trend, the resulting TOA SW flux trend changed from −0.35 ± 0.18 Wm −2 per decade to 0.61 ± 0.18 Wm −2 per decade, with all other inputs and algorithms being identical. The global daytime mean cloud optical depth is approximately 4.5, thus 0.6 per decade trend represents a 13% change of the global mean value. The global mean SW flux is on the order of 100 Wm −2 , thus 0.96 Wm −2 per decade trend change represents a 1% change of the global mean value. This means that a 13% change in cloud optical depth results in about 1% change in the SW flux.
CERES instruments on different satellites provide valid fluxes at different viewing zenith angle ranges. CERES instruments on Terra and Aqua only provide fluxes for VZA less than 66 • whereas CERES instruments on S-NPP and NOAA-20 provide fluxes for VZA less than 72 • . Satellite retrieved cloud fraction generally increases as the VZA increases because of viewing perspective of the sensor, increase in path length, and in pixel size. Global mean cloud fractions retrieved only using near-nadir view footprints can be 8% less than the cloud fraction retrieved only using oblique-view footprints. The global mean cloud fraction difference is 0.5% with and without the footprints with VZA > 66 • . This VZA dependency also propagates to other cloud properties, and the retrieved cloud optical depth generally decreases as VZA increases. However, there are regional exceptions and the dependency of liquid cloud optical depth on VZA differs from that of ice clouds. Different range of VZAs on different CERES instruments can affect the continuity of the ERB climate data record, as the cloud properties depend on VZA and that CERES empirical ADMs are functions of these cloud properties. The global monthly mean instantaneous SW flux difference averaged using all footprints and using only footprints with VZA < 66 • from CERES S-NPP observations are about 0.7 to 0.9 Wm −2 . Regionally, the flux can differ by up to 10 Wm −2 with large differences around ITCZ and near 60 • N and 60 • S. LW flux differences due to different VZA ranges are smaller comparing with the SW flux difference, with daytime and nighttime monthly mean LW flux difference less than 0.3 Wm −2 . Regionally, the LW flux differences up to 4 Wm −2 are seen over the tropical oceans near the ITCZ and over the western Pacific Ocean. These flux differences caused by using different VZA ranges can result in a discontinuity in CERES TOA flux. To ensure continuity, the viewing zenith angle ranges common to all CERES instruments are used to produce the long-term Earth's radiation budget climate data record.
Using collocated CERES footprints on S-NPP and Aqua, cloud fraction and cloud optical depth retrieved from MODIS and VIIRS are compared against each other. The cloud fractions from these two imagers agree reasonable well for both polar and non-polar regions, with the bias around 1% and the RMS errors around 7%. The cloud optical depths show good agreement over the non-polar regions, but VIIRS produces higher cloud optical depth than MODIS over the polar regions. Algorithm development is underway to minimize the cloud property differences between MODIS and VIIRS, which includes developing algorithm that uses channels that are common to both MODIS and VIIRS and placing them on the same radiometric scale.
Creating CERES TOA flux climate data records needs to ensure consistency not only in the radiance measurements across all CERES instruments over the last two decades and beyond, but also in the input data sets used. Producing CERES TOA flux represents the highest degree of data fusion of CERES, MODIS, VIIRS, and many operational geostationary satellites, together with input data sets of snow/ice map, aerosol assimilation model, and meteorological fields from GEOS assimilation system. To ensure the highest fidelity of the TOA SW flux for detecting and quantifying climate variability, consistency in satellite measurements as well as in input data used must be maintained.
Author Contributions: W.S., L.L., H.W. and Z.A.E. analyzed the data; W.S. wrote the manuscript; All authors provided edits to the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.