Validation of CERES Flight Model 5 in-orbit calibrations using lunar observations Janet L. Daniels*a, George L. Smitha, Susan Thomasa, Kory J. Priestleyb aScience Systems and Applications, Inc., Hampton, Virginia, U. S.; bNASA Langley Research Center, Hampton, Virginia, U. S.

ABSTRACT Scientific studies require radiation fluxes over the to be accurate within 1% for shortwave fluxes and 0.5% for outgoing longwave fluxes. The validation of in-orbit instrument performance requires both stability in calibration source and also calibration corrections to compensate for instrument changes. The offers an external source whose signal variance is predictable and non-degrading. CERES detectors register the signal output from the entire face of the Moon. Lunar observations performed by CERES Flight Models (FM) 1 through 4 have been successful in assisting validation of radiances to the required accuracy. CERES Flight Model 5 (FM-5) is on Suomi/NPP spacecraft, orbiting since October 2011. This paper uses lunar measurements to validate detector output of FM-5. These measurements are adjusted to remove orbital effects due to variations in distance between Moon and Sun, distance between the satellite and the Moon and lunar phase angle. These effects create a total variation in lunar irradiance of 20% in the total channel and 8% in the shortwave channel. The change in orientation of the Moon as seen by the detector is called and causes variations of about 1% of the irradiance. A consistent dataset spanning at least 2 years in length is required to remove variations due to libration. The major uncertainties remaining in the measurements are assumed to be due to changes of the spectral responses of the channels due to degradation of optical surfaces in orbit. The results demonstrate that lunar observations can be used to validate FM-5 measurements. Keywords: Clouds and Earth Radiant Energy System, Earth radiation budget, lunar measurements, radiometry, remote measurements, calibration, validation.

1. INTRODUCTION

The understanding of inter-annual climate variations requires measurements of the Earth’s radiation budget over a period of years. The Clouds and Earth Radiant Energy System (CERES) Project was created for this need1. CERES Flight Models FM-1 and -2 have been operating on the Terra spacecraft starting in February 2000 and on Aqua FM-3 and -4 since June 2002. To assure continuity of the Earth Radiation Budget Climate Data Record, CERES FM-5 was placed into orbit on the Suomi/NPP in October 20112.

The radiation fluxes over the Earth should be accurate within 1% for shortwave fluxes and 0.5% for outgoing longwave fluxes. Maintaining this accuracy demands frequent in-flight calibration3,4. The basis of in-flight calibration is the Internal Black Body System and the Shortwave Internal Calibration System5. These are augmented by using Deep Convective clouds as a vicarious calibration target6.

A thorough validation protocol7 is implemented to assure the resulting accuracy. The first set of validation methods is comparing radiances from the different CERES instruments. The FM-1 and -2 instruments can be compared directly as they are both on the Terra spacecraft. Likewise the FM-3 and -4 on Aqua can be compared. FM-1 and FM-3 can be compared near 70oN where the orbits cross8. Similarly FM-5 measurements are compared with FM-1 radiances9. Because the Aqua and Suomi/NPP orbits are so close, there are ways to compare FM-5 and FM-3 frequently10. The Tropical Mean Method11 provides a check between the channels. Daniels et al.12 developed a method of validation using lunar observations near by the Clouds and the Earth’s Radiant Energy System (CERES) and applied the method to the four CERES instruments on the Terra and Aqua spacecraft. The present paper uses lunar observations to validate the FM-5 measurements from 2011 to 2019. First the CERES method and its application to each of the three channels are described, then the lunar observation used to compute the changes of each of the FM-5 channels. The results are compared to those obtained from FM-1 through -4. *[email protected]; phone 1 757 864-2778; fax 1 757 864-6319; nasa.gov

2. METHODOLOGY

The validation of in-orbit instrument performance requires both stability in calibration source and also a method to correct calibrations to compensate for instrument changes. Unlike internal calibration systems, the Moon is an external source whose signal variance is predictable and non-degrading. Figure 1 shows the viewing geometry during lunar observations. The observations are made near full Moon so as to have near the brightest and hottest condition to get a good dynamic range. At this time the Moon is near the Sun-Earth line. As the Suomi/NPP spacecraft comes over Antarctica, the instrument elevation gimbal is fixed and FM-5 is rotated slowly in azimuth to bring the Moon through the field of view several times during the viewing opportunity. Thus, there are several measurements for each full Moon. Because the instrument is not scanning in elevation angle and the rotation in azimuth is slow compared to the detector response time, the dynamic point response function does not apply. The static point response function is constant over the field of view except near the edges, where the optical blur (0.16o) causes it to decrease.

Figure 1: Suomi/NPP orbital location during lunar observations. Figure 2 shows the image of the Moon at mean Earth-Moon distance within the CERES field of view. Because of the variation of albedo over its surface, as the Moon rotates relative to the instrument the net albedo changes. These create changes in the measurements. The radius of the Moon’s image varies inversely with the Satellite-to-Moon distance due to the ellipticity of the Moon’s orbit, so the area of the image and the measurement decrease as the inverse square of the distance.

Figure 2: Image of Moon in CERES field of view. The insolation on the Moon varies because of the variation of the distance between Sun and Moon. The measurement also varies due to the change of distance between the instrument and the Moon. The lunar phase angle, i.e. the lunar central angle between the lunar sub-solar point and the sub-satellite point, determines how well the part of the Moon which is viewed is sun-lit. Daniels et al.12 demonstrated that when these factors are taken into account, there is about a 1% variation remaining, due to lunar libration. They showed that these changes of irradiance for the three channels averaged out well below the accuracy requirements over a two year period. This procedure is be used for validating the in-flight calibration of FM-5. 3. ANALYSIS

Edition 1-CV data are used in this paper. For each channel, lunar measurements were normalized. Figure 3a shows the normalized measurements by the total channel for lunar observations from February 2012 through June 2016. The satellite- to-Moon distances for the measurements are also shown, as a solid line. The total channel measurements are adjusted to the mean satellite-to-Moon distance using the inverse square and is shown in figure 3b. There is still a strong cycle present. The distance from Moon-to-Sun is also plotted and the two curves are quite similar and are out of phase 180o. The inverse square relation is used to adjust the measurements to a distance of one astronautical unit and the result is plotted in figure 3c. The time variation is nearly periodic with a range of ±0.02. Much of this variation is due to the lunar libration.

Figure 4 shows the shortwave channel measurements analyzed in the same manner. The scatter of the measurements for each are greater than for the total channel. The albedo of the Moon is about 0.1, so the variations of albedo over the surface of the Moon cause a large relative change in the solar radiation. Because the total channel is measuring both the reflected and the emitted longwave fluxes, these variations of albedo make very little difference to that channel. When the satellite to Moon distance and the Moon to Sun distance are taken into account, the variations due to lunar librations remain in figure 4c. These are at the ±0.02 level, as for the total channel. The longwave channel measurements are shown in figure 5, which looks very much like those of the total channel in figure 3. As noted earlier, 0.9 of the lunar radiation is emitted radiation.

a)

b)

c)

Figure 3: Total channel radiances: (a) Variation of normalized radiances and satellite to Moon distance (solid line), (b) Variation of normalized radiances adjusted for satellite to distance and Moon to Sun distance (solid line), (c) Radiance after adjustments for satellite-to-Moon and Moon-to-Sun distances.

a)

b)

c)

Figure 4: Shortwave channel radiances: (a) Variation of normalized radiances and satellite to Moon distance (solid line), (b) Variation of normalized radiances adjusted for satellite to distance and Moon to Sun distance (solid line), (c) Radiance after adjustments for satellite-to-Moon and Moon-to-Sun distances.

a)

b)

c)

Figure 5: Window channel radiances: (a) Variation of normalized radiances and satellite to Moon distance (solid line), (b) Variation of normalized radiances adjusted for satellite to distance and Moon to Sun distance (solid line), (c) Radiance after adjustments for satellite-to-Moon and Moon-to-Sun distances. The effect of lunar phase angle (LPA) on the measurements is determined empirically. Figure 6 shows the measurements after they have been adjusted to a mean satellite-to-Moon distance and Moon-to-Sun distance of 1 A.U. as a function of lunar phase angle LPA. Measurements for all five CERES flight models are included. Note that the ordinate for the shortwave channel is greater than for the total and longwave channels. A second order fit has been made for each channel for LPA between 5o and 12o. For LPA less than 5o the irradiance increases sharply for the shortwave channel due to the opposition effect13. If the Moon is in opposition to the Sun from the viewpoint of the observer, there are no shadows cast by the terrain and the surface is fully sun-lit. The albedo of the Moon is about 0.1, so while a 0.01 change of albedo is significant in the shortwave channel, it is small compared to the absorption of 0.9, which governs the emitted radiation, of which the 8 to 12 window channel is a part. The increase of reflected shortwave flux is measured by the total channel, causing an increase for LPA less than 5o but not as much as for the shortwave channel. For use of the Moon for validating the calibrations, only LPA between 5o and 12o are used.

Figure 6: Normalized CERES detector responses dependence on Lunar Phase Angle. The measurements as adjusted in figures 3c, 4c and 5c are now adjusted to eliminate the LPA effect and are shown by figure 7. The top panels of figure 7 show the measurements adjusted for distances over the time period. The second-order fit for effect of LPA is divided out, giving the bottom panels. Much of the variability is removed from the data, leaving primarily an oscillation which is attributed to lunar librations. Before the phase angle effect is removed, the total channel had a variation of about ±0.15. With the LPA taken into account, the total channel has an oscillation of about 0.02 amplitude remaining. The shortwave channel has a variation of ±0.2 before the LPA is taken into account, and afterwards shows a cycle with amplitude 0.05. The window channel has variations of ±0.15 before the LPA is taken into account and ±0.03 afterwards. These lunar observations were taken by FM-5 in the same manner and processed in the same way as those by FM-1 through FM-4. The results are consistent and demonstrate that the method of using lunar observations developed by Daniels et al.12 can be used with FM-5 to validate in-flight calibrations through its operational lifetime.

Figure 7: Normalized measurements of lunar observations from all three CERES channels adjusted for satellite to Moon and Moon to Sun distances without accounting for lunar phase angle LPA (top panels) and with accounting for lunar phase angle LPA (bottom panel). The final cyclic perturbation to be removed is due to lunar librations. Adjusted data for each channel are plotted with respect to sub-satellite lunar latitude and longitude. The plots shown in figure 8 show the repeatability based on lunar orientation to the satellite at the time of observation.

Figure 8: Normalized measurements of lunar observations from all three CERES channels, with lunar phase angle effects removed, plotted with respect to the sub-satellite latitude and longitude. To remove this final effect, a second-order linear fit is calculated and applied, described in more detail in Daniels et al.12. Figure 9 shows the measurements for all lunar observations for each of the three channels of FM-5 over the data period. The y-axis shows percent change per channel and the x-axis is elapsed days. The average for each month is denoted by an asterisk symbol. The trend lines have slopes of 0.019% per year for the total channel, 0.058% for the shortwave channel and 0.016% for the window channel.

Fig. 9: Response history of total, shortwave and window channels adjusted for spacecraft-to-Moon, Moon-to-Sun distance, lunar phase angle and lunar librations effects.

4. CONCLUSIONS

CERES Flight Model 5 measurements of the Moon’s irradiance have been used to validate the in-flight calibration results of that instrument. The measurements from the three channels of FM-5 are adjusted for orbital geometry and lunar phase angle and libration effects, and the results are consistent with those from FM-1 through -4, demonstrating that FM-5 can be validated by use of lunar observations over its operational lifetime. It is demonstrated that over the data period each of the three channel varies by less than ±0.1%.

ACKNOWLEDGEMENTS

The authors thank the Science Directorate of Langley Research Center and the Science Mission Directorate of the Earth Science Division of NASA for the support of the CERES Project.

REFERENCES

[1] Barkstrom, B. R.: “Earth radiation budget measurements, Pre-ERBE, ERBE, and CERES,” Proc. SPIE, v 1299, 52-60, 1990. [2] Priestley, K. J., S. Thomas and G. L. Smith: 'CERES Flight Model 5 On NPP: Post-Launch Performance And Initial SDR Validation Results,' IGARSS, 2012. [3] Priestley, K. J. et al., “Radiometric Performance of the CERES Earth radiation budget climate data record sensors on the EOS Aqua and Terra spacecraft through April 2007,” J. Atmos. Oceanic Tech., vol. 28, pp. 3-21, 2011. [4] Thomas, S. et al., “Performance Assessment of the Clouds and Earth’s Radiant Energy System (CERES) Instruments on Terra and Aqua Spacecraft,” Proc. SPIE vol. 8866, 2013. [5] Smith, G. L., K. J. Priestley and N. G. Loeb: “Clouds and Earth Radiant Energy System: from Design to Data,” IEEE Trans. Geosci. and Rem. Sens., 52, pp. 1729-1738, 2014. [6] Kratz, D. P., K. J. Priestley, and R. N. Green, “Establishing the relationship between the CERES window channel and total channel measured radiances for conditions involving deep convective clouds at night,” J. Geophys. Res., vol. 107, no. D15, pp. ACL 5-1– ACL 5-11, 2002. [7] Priestley, K. L., G. L. Smith, S. Thomas and G. Matthews, “Validation Protocol for Climate Quality CERES Measurements,’ Proc. SPIE, 2007. [8] Szewczyk, Z. Peter, G. L. Smith and K. J. Priestley: “Validation of Clouds and Earth Radiant Energy System instruments aboard the Terra and Aqua satellites”, J. Geophys. Res._Atmos., 110(D2), 2010. [9] Szewczyk, Z. Peter, G. L. Smith and K. J. Priestley: “Comparison of unfiltered radiances measured by CERES instruments aboard the Suomi/NPP, Terra and Aqua satellites,” Proc. SPIE, 2014. [10] Szewczyk, Z. Peter, G. L. Smith and K. J. Priestley, “Comparison of unfiltered radiances measured by CERES instruments aboard the Suomi/NPP, Terra and Aqua satellites,” Proc. SPIE, 2014. [11] Smith, G. L. Susan Thomas1, Kory J. Priestley2, Dale Walikainen, “Tropical Mean Fluxes: A Tool for Calibration and Validation of CERES Radiometers,” Transactions of Geophysics and Remote Sensing, vol. 54, pp: 5135 - 5142, 2016. [12] Daniels, J. L., G. L. Smith, K. J. Priestley and S. Thomas: Using Lunar Observations to Validate In-Flight Calibrations of Cloud and Earth’s Radiant Energy System Instruments,” Trans. Geophys. And Rem. Sens., vol. 53, 5110-5116, 2015. [13] Pohn, H. A., H. W. Radin and R. L. Wildey, “The Moon’s photometric function near zero phase angle from Apollo 8 photography., Astrophys. J., vol. 157, pp. L193-L195, 1969.