2.1. MODIS
For both MODIS instruments, a linear relationship or algorithm between the incident radiance (
L) and detector response (
dn*) is applied for the RSB calibration and retrieval,
where
m1 is the calibration coefficient (inversely proportional to the gain of a given detector) derived with reference to the SD bi-directional reflectance factor (BRF),
Esun is the solar spectral irradiance at an Earth–Sun distance of 1 astronomical unit (AU) and integrated over the RSR for each detector,
dn∗ is the detector digital response corrected for instrument background and temperature effects, and
RVS is the response versus scan angle, which accounts for the instrument gain variations as a function of the angle of incidence (AOI) of light relative to the scan mirror. Since the SD is used primarily for MODIS RSB calibration, the RSB RVS is conveniently normalized at the AOI of its SD view, i.e.,
RVSSD = 1. For the EV observations, the MODIS RSB L1B primary data product is the top-of-atmosphere (TOA) reflectance factor,
, where
θEV is the solar zenith angle of the EV pixel. The EV radiance,
LEV, and the reflectance factor,
, can be easily converted to each other by multiplying or dividing a factor of
, with
dES being the Earth–Sun distance (normalized at 1 AU) at the time of sensor observation. The solar spectral irradiance used in the MODIS RSB is a combination of Thuillier et al. (1998; 0.4–0.8 μm), Neckel and Labs (1994; 0.8–1.1 μm), and Smith and Gottlieb (above 1.1 μm) [
9]. MODIS L1B calibration algorithms produce both radiance and reflectance data products for the RSB.
The on-orbit calibration coefficient m1 and RVS in Equation (1) change with time and are thus updated regularly. For both MODIS instruments, the m1 and the RVS are currently derived by using the SD calibration, lunar calibration, and EV observations over select desert sites at multiple AOIs.
The MODIS SD is a flat and near-rectangular panel made of Spectralon with a near-Lambertian reflectance profile. It is located inside the instrument cavity. The SD provides diffusely reflected sunlight that can be used for the RSB calibration. The SD panel can be illuminated by the sun when the instrument passes the Earth terminator from the nighttime side to the daytime side. Only the responses to the fully illuminated SD are used to compute the calibration coefficients.
Figure 1a shows a schematic of the MODIS scan operation that enables data to be collected each scan from its on-board calibrators and the EV. The left side of
Figure 1b shows the AOI of each of the OBC in relation to the EV data. During each SD calibration, the solar radiance diffusely reflected from the SD can be accurately calculated by
where
ρSD is the SD BRF derived from prelaunch measurements,
θSD is the solar zenith angle relative to the SD, Δ
SD is the SD on-orbit degradation, and
τSDS is the SD screen (SDS) transmission function, which is also referred to as the vignetting function (VF). During sensor nominal operations, the SDS can be commanded to an open or a closed position, thus providing two different levels of the intensity for the sunlight illuminated on the SD surface. Placed in front of the SDS is an aperture door that is only opened during nominally scheduled SD and SDSM calibration events. For Terra MODIS, however, the SD door has been fixed in the open position with the SDS in the closed position since 2 July 2003, resulting from an anomaly related to its SD door and/or SDS operation. The SD BRF was measured prelaunch, and its relative profile was validated on-orbit using measurements made during spacecraft yaw maneuvers, which were performed early in the mission for both Terra and Aqua MODIS. The SD on-orbit degradation, Δ
SD, is tracked by the onboard SDSM. For MODIS, the SDS VF was not fully characterized prelaunch, and it was derived on-orbit from measurements during yaw maneuvers made with and without the SDS in place. When the SDS is placed in the open position during an SD calibration event,
τSDS in Equation (2) becomes a constant of 1. Otherwise, it varies with the solar illumination angle relative to the SDS. In Equation (2), an assumption that the SD degrades uniformly with respect to incident and outgoing directions has been applied, meaning that the degradation measured at the output angle of the SDSM will be the same as in the direction of the RSBs. Therefore, the SD on-orbit BRF can be expressed as the product of its prelaunch BRF,
ρSD, and its on-orbit degradation, Δ
SD [
47]. However, evidence of non-uniformity in the SD degradation has led to the use of EV data to supplement the calibration of the MODIS RSB, which will be discussed further below.
By applying Equation (1) to the SD view and substituting
LSD in Equation (2) to Equation (1), we can derive the calibration coefficient by
The calibration coefficient is calculated for each band, detector, and mirror side for the 1-km RSB as well as for each sub-frame for the 500-m (2) and 250-m (4) resolution RSBs.
As previously mentioned, the SD on-orbit degradation, Δ
SD, in Equation (3) is tracked by the SDSM, which functions as a ratioing radiometer that views the SD, the Sun through its Sun-view port, and an internal dark scene, alternately. The MODIS SDSM has nine detectors, and each detector tracks the SD degradation at a discrete wavelength. The center wavelengths of the SDSM detectors cover a spectral range from 412 nm to 936 nm. After corrections applied for the view geometry effects, the ratios of the background subtracted digital count for the SD view to the background subtracted digital count of the Sun view provide the trends of the SD on-orbit degradation. The response of the SDSM detectors is assumed to be linear. It was also assumed that the SD degradation between its prelaunch characterization and its first on-orbit measurement was negligible since the SD exposure to the environment was minimal. Consequently, the SDSM ratios normalized to its first on-orbit SD measurements are used to track the SD on-orbit degradation at the wavelengths of its detectors. A linear interpolation approach is applied to obtain the SD degradation at any wavelength in the range from 0.412 µm to 0.936 µm from the measured SD degradations at the nine center wavelengths of the SDSM detectors [
48,
49].
The lunar surface has a long history of use as a calibration target for the reflective bands of a number of satellite instruments, owing to its stable reflective surface and lack of atmosphere [
50]. Since the lunar surface is not smooth, only the integrated lunar irradiance is used in the MODIS lunar calibration methodology. Using Equation (1), the measured lunar radiance from individual detectors can be easily calculated and their corresponding integrated lunar irradiance (
I) can be expressed by
where
N is number of the scans used in the computation, each of which fully covers the lunar surface, and
ω = 1/(705
Stot)
2 is the solid angle (steradians) of each pixel.
Stot is the number of sub-frames of each band (or detector) and 705 km is the nominal orbital altitude. The value of 1/
Stot corresponds to the HSR at the nadir in km. The summation is made over detectors, frames, and select scans. In this analysis, the lunar irradiance is calculated using the MODIS C7 LUTs prepared for the upcoming L1B reprocessing. Instead of using only scans that fully cover the lunar surface, the lunar irradiance can also be calculated using the measurements from all scans in a lunar observation event. Using all scans requires the creation of lunar images corresponding to each detector in order to be able to use images where the Moon is only partially on the FPA. These detector images are created by extracting single lines from each scan corresponding to the detector as the Moon moves across the FPA. This approach, however, requires a correction for the oversampling effects. Oversampling occurs because the Moon does not move on the FPA by 1 pixel per scan, causing some parts of the lunar surface to be observed more than once in a detector image (higher resolution bands can be undersampled under certain observation geometries). The all-scan approach has a relatively large uncertainty due to corrections applied for the oversampling effect, but it can help examine calibration differences among individual detectors. In this analysis we focus on the methodology that uses the scans with full coverage of the lunar surface as described by Equation (4).
Figure 2a,b show examples of the lunar images acquired by Aqua MODIS bands 1 and 8 during the scheduled lunar observation on 24 January 2021. Also shown in
Figure 2 are the lunar images acquired by the SNPP VIIRS on the same day for its bands I1 (c) and M1 (d).
By comparing the integrated lunar irradiance predicted by the ROLO model (
IROLO) with that measured by the MODIS using Equation (4), the band-averaged calibration coefficient,
, can be computed by
where
are the calibration coefficients derived from the first SD on-orbit calibration and the angled brackets indicate an average over all detectors of the band. To calibrate the RSB using the Moon, a reference for the lunar irradiance is required. In this analysis, the lunar irradiance for each calibration in Equation (5) was provided by the Robotic Lunar Observatory (ROLO) model prediction, developed by the USGS [
31]. Since the absolute uncertainty of the current ROLO lunar model is larger than the MODIS calibration specification of 2%, the MODIS lunar calibration is only used to track the RSB on-orbit changes. In this case, the constant terms in Equation (5) can be omitted in routine data processing [
28].
Among the RSBs, bands 13–16 partially saturate when they observe the Moon. This saturation occurs at the inner part of the illuminated lunar surface (with unsaturated pixels near the disk edge and terminator), which is typically at the highest radiance levels. To correct for saturation, a ratio approach is applied to replace the saturated pixels using band 18 as a reference band [
51]. To obtain the ratio, the saturated band data is plotted versus the spatially co-registered reference band data at the pixel level and fit to a linear equation for all unsaturated pixels, where the slope represents the ratio between the two bands. The saturated data can then be replaced by multiplying the reference band data by the ratio at the location of the saturated pixels. For SWIR bands, there are strong crosstalk contaminations among themselves and from mid-wave infrared bands as well as the large out-of-band (OOB) RSR contributions at the wavelength of 5.3 µm. These contaminations need to be mitigated before the calculation of the lunar irradiance using Equation (4). Accurate mitigation of these effects is still a challenging issue and needs more effort [
52]. In this analysis, MODIS lunar calibration is mainly focused on the VIS and NIR bands.
For MODIS RSBs, the calibration coefficients,
m1, and the RVS are needed to produce the L1B products as shown in Equation (1). Due to the non-uniformity of the SD degradation with respect to the incident and outgoing directions, the SD degradation measured at the SDSM view direction may deviate from that at the RSB view direction, resulting in a long-term bias in the calibration coefficients derived from the SD, especially for short wavelength RSB that have experienced more significant degradation on-orbit. As a result, EV response trends from pseudo-invariant desert sites at the SD AOI are used to correct the long-term drifts in SD-based calibration coefficients for the short wavelength bands. Combination of the SD calibration results and EV response trends at the same AOI help produce the calibration coefficients with both long-term accuracy and short-term stability [
53].
MODIS RSBs view the SV, through which the Moon is also observed, and the SD at different AOIs to the scan mirror, with the SV at 11.25° and the SD at 50.25°. The trending differences of the two calibration results provide the information that is directly related to on-orbit changes in the RVS.
Figure 3 shows the SD and lunar gain trending for MODIS bands 1 and 8. For both Terra and Aqua MODIS, the shortest wavelengths have experienced the most gain and RVS changes. To date, the band 8 (412 nm) gains have changed (decreased) up to 40% for Terra MODIS and more than 45% for Aqua MODIS based on their SD and lunar calibrations. In comparison, the NIR band 1 (646 nm) shows a gain change of less than 20%. The temporal divergence between the SD and lunar gain measurements is a result of the evolution of the on-orbit RVS. Accurate characterization of the on-orbit RVS is extremely important for the MODIS RSB on-orbit calibration, especially for the short wavelength bands.
Initially, the RSB time-dependent RVS was derived by using the lunar and SD calibration differences with an approximation that the RVS on-orbit change for a given RSB is a linear function of the AOI. As each mission continues to operate beyond its designed lifetime, this approximation no longer meets the L1B calibration accuracy requirements, especially at short wavelengths. As a result, the EV response trends at multiple AOIs have been used together with on-orbit SD and lunar measurements to track on-orbit changes in the RVS for a few select bands, starting from L1B Collection 6 (C6) for both Terra and Aqua MODIS. In C6, we note that the lunar results are not used for bands 1 and 2 EV time-dependent RVS derivation due to the disagreement of lunar measurements with the EV response trending from the desert sites [
53].
2.2. VIIRS
Similar to MODIS, a simple smooth function is applied to establish the relationship between the incident radiance and detector digital response for the VIIRS RSB. For all SNPP VIIRS RSB and N20 VIIRS VIS and NIR bands, a quadratic approximation is applied, while for N20 VIIRS SWIR a third order polynomial is used due to a significant nonlinearity effect for these bands [
54]. The relationship between the incident radiance and instrument response for the VIIRS RSBs can be written as
where
ci (
i = 0, 1, 2, 3), are the prelaunch measured calibration coefficients of the polynomial,
F, called F-factor, which is the ratio of the on-orbit coefficients of the polynomial at the time of the measurement to the prelaunch coefficients, assuming that the coefficients of the polynomial change proportionally with each other on-orbit,
dn is the background subtracted instrument response, and
RVS is the response versus scan angle of the half-angle mirror (HAM). The calibration coefficients,
c0,
c1,
c2, and
c3, are instrument and electronic temperature-dependent. Both
F and the
RVS in Equation (6) may, in principle, change temporally on-orbit. However, there is no evidence to suggest that the
RVS has a noticeable on-orbit change for either VIIRS instrument. As a result, only the
F-factors have been updated regularly on-orbit by using SD and lunar calibrations on an as-needed basis. For the VIIRS EV, the TOA radiance is the primary L1B product, which can be easily converted to its TOA reflectance factor. The VIIRS L1B products are also referred to as the sensor data records (SDR).
VIIRS has the same type of SD as MODIS. The radiance of the sunlight diffusely reflected from the VIIRS SD can also be calculated by Equation (2). The VIIRS SD BRF, ρSD, and the SDS transmittance, τSDS, were measured prelaunch and refined on orbit by measurements made during yaw maneuvers. The VIIRS SD degradation is tracked by the on-board SDSM at eight discrete wavelengths, compared to MODIS at nine different wavelengths of the same spectral range. The VIIRS SD port has a permanently fixed attenuation screen, but it does not have a dedicated door cover as MODIS does. This means that the SD is illuminated by the Sun every orbit. The assumption applied to the MODIS SD calibration that the SD degrades uniformly with respect to incident and outgoing directions is also initially applied to the VIIRS SD calibration. However, there is evidence of non-uniform degradation based on the differences between the lunar and SD trends as will be discussed below.
Comparing the predicted solar radiance (
) with that measured using Equation (6) via detector response (
to the SD, the F-factors for the VIIRS RSB calibration can be calculated by
where
is the RVS at the AOI of the SD. The VIIRS RSB calibration is performed for each scan using the detector’s average response to the SD averaged over the scans in a selected illumination angle range, known as the “sweet spot”. The
RSR(λ,
t) in Equation (7) is time-dependent for SNPP VIIRS RSBs as a result of the wavelength-dependent degradation in SNPP rotating telescope assembly (RTA) optics and large OOB RSR contributions [
52,
55]. The SD calibration is performed for each orbit and the F-factor is derived for each RSB detector, HAM side, and gain stage for the dual gain bands. Compared to Equation (4), designed to derive the MODIS RSB reflectance calibration coefficient, Equation (7) is used to compute the radiance calibration coefficients for the VIIRS RSB. In addition to the SD bi-directional reflectance function and solar attenuation screen transmission, the sensor’s solar spectral irradiance is also needed to determine the predicted radiance reflected off the SD (
) used to compute the F in Equation (7). The SNPP uses Kurucz spectra from MODTRAN 4.3 while N20 uses the Thuillier spectra [
56].
Applying Equation (6) to the SV lunar observations, the integrated lunar irradiance measured by a VIIRS RSB can be calculated using
where
ω = (
S/824)
2 is the solid angle (steradians) of each pixel of the band.
S is 0.375 km for an I-band and 0.75 km for an M-band and 824 km is the nominal orbital altitude. The VIIRS SV has the same AOI as its SD as shown in
Figure 1b, at which the RSB RVS is normalized; thus, the
in Equation (8) is equal to 1. Similar to the lunar calibration for MODIS, the VIIRS lunar calibrations use only the
N scans in which the full disk of the Moon can be observed during each scheduled lunar observation. For all regularly scheduled lunar observations of both SNPP and N20 VIIRS, the gain stages of all dual-gain bands are fixed at high gain. Examples of SNPP lunar images for bands I1 and M1 are also shown in
Figure 2. For VIIRS lunar observations, only band M7 in N20 has shown any signs of saturation, and even then, for only a few pixels. To correct this, the same approach as used for MODIS bands 13–16 is employed, this time with band M5 as a reference. Similar to the MODIS configuration, a sector rotation is applied to collect lunar data in the EV data sector. However, the VIIRS EV has three different aggregation regions. SNPP and N20 VIIRS lunar data are collected in different aggregation regions and special attention should be paid to the summation over pixels along the scan direction in Equation (8). In this analysis, prelaunch
RVS and C2 F-factor LUTs are applied in Equation (8) for SNPP and C2.1 LUTs for N20 VIIRS.
Same as for the MODIS RSBs, the impact of the detector difference on the calibration coefficients derived from a lunar observation for a VIIRS RSB can be assumed to be negligible. Then the detector-averaged relative
F factor can be derived from each of the scheduled lunar observations using
The predicted lunar irradiance for each lunar observation event, I, is provided by the ROLO model. To distinguish the F-factor derived from the SD/SDSM calibration, a superscript “Moon” is added to the F. Similar to MODIS, the VIIRS lunar calibration is only used to track its RSB on-orbit changes and some of the constant parameters in Equation (8) are dropped in Equation (9). The VIIRS lunar calibration coefficients are scaled by normalizing the lunar F-factors derived from Equation (8) to the corresponding F-factors derived from the SD/SDSM calibration at the time of instrument launch to obtain absolute values of the lunar calibration coefficients.
VIIRS RSBs view the SD and the SV at same AOI of the HAM and thus the SD calibration and lunar calibration should provide identical on-orbit changes for the RSBs if both the SD and lunar calibration results are accurate. It is known that the SD degrades non-uniformly with respect to incident and outgoing direction [
45]. Thus, the SD degradation from the SDSM view direction when applied to the RSB view direction may result in long-term biases in the F-factors derived from the SD calibration, especially for short wavelength bands, as also confirmed by EV measurements. A comparison between the two sets of F-factors can identify the long-term biases in the SD F-factors and can be used to obtain the SD degradation differences between the two view directions. Combining the SD and lunar calibration results provides the RSB F-factors with both long-term accuracy and short-term stability.
Similar to the MODIS gains shown in
Figure 3, the VIIRS gains are shown in
Figure 4 for the VIS band M1 (412 nm) and NIR band I1 (645 nm). Unlike the RVS-caused separation between the lunar and SD gains in MODIS, the separation between the two sources is a result of the inadequacy in SDSM to accurately characterize the non-uniform degradation in the SD, which manifests as a divergence with the lunar gain. The lunar data can be used as a method for correcting the SD degradation trends so that the SD and lunar trends agree.