S-NPP VIIRS Thermal Emissive Bands 10-Year On-Orbit Calibration and Performance

The Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership Program (S-NPP) satellite, launched in late 2011, has reached the decade landmark under successful operations. VIIRS has 22 spectral bands, 7 of which are thermal emissive bands (TEB) that cover the 3.70 to 11.84 μm wavelength range. Over the years, VIIRS TEB observations have been used to generate several data products (e.g., surface/cloud/atmospheric temperatures, cloud top altitude, and water vapor properties). The VIIRS TEB calibration uses a quadratic algorithm and is referenced to an on-board blackbody with temperature measurements traceable to the National Institute of Standards and Technology standard. This manuscript provides an overview of the VIIRS instrument operations and TEB calibration activities and algorithms used in the level 1B data and describes the TEB on-orbit performance for S-NPP VIIRS. The 10-year on-orbit performance of the S-NPP VIIRS TEB has generally been stable, and the degradations in the S-NPP TEB detector responses are minor after a decade in orbit. The noise characterization performance repeatedly meets the design requirements for all TEB detectors as well. On-orbit changes in the TEB response-versus-scan-angle, based on pitch maneuver observations, have been demonstrated to be extremely small. Moreover, multiple time series over select ground targets have shown that the sensor’s on-orbit performance is quite stable.


Introduction
Launched on 28 October 2011, the Suomi National Polar-orbiting Partnership (S-NPP) satellite's primary sensor, the Visible Infrared Imaging Radiometer Suite (VIIRS), has operated successfully for a decade. Designed with substantial Moderate Resolution Imaging Spectroradiometer (MODIS) heritage, VIIRS has 22 spectral bands that span from the visible to infrared wavelengths (0.4-12.5 µm). These spectra are intended to support a cohort of environmental data records (EDR) that assist users in the land, ocean, and atmospheric science communities [1][2][3]. VIIRS bridges the gap between the data collected by the Earth Observing System (EOS) satellites and the future generation of Earth observatories-the Joint Polar Satellite System (JPSS), which carries (NOAA-20) or will carry instruments like the ones onboard S-NPP. The S-NPP spacecraft operates in a near Sun-synchronous polar orbit (nominal altitude~828 km) and at a 98 • inclination angle (relative to the Equator; crossing time is 1:30 p.m.) [4,5]. With an Earth-view (EV) scan-angle range of about ±56 degrees, the VIIRS sensor can make continuous global observations twice daily. The VIIRS sensor intends to perpetuate, advance, and safeguard the datasets-whose calibration is traceable to ground-based references-provided by the MODIS instruments (effectively operating onboard the EOS Terra and Aqua spacecraft). As such, VIIRS was built to succeed and refine the MODIS sensor design. The S-NPP VIIRS instrument was built and tested prelaunch by its vendor, Raytheon Santa Barbara Remote Sensing, while support on the calibration data analyses was provided by the VIIRS Characterization Support Team (VCST) at the VIIRS has five M-thermal emissive bands (TEB) and two I-TEB. Located in the S/MWIR FPA are bands I4, M12, and M13, while bands I5, M14, M15, and M16 are situated on the LWIR FPA. Band M16 has two separate (and time delay integrated) arrays of 16 detectors. Only band M13 is dual gain-its low gain is primarily used for fire detection. Table 1 shows a summary of key design specifications for the VIIRS TEB. The VIIRS TEB VIIRS has five M-thermal emissive bands (TEB) and two I-TEB. Located in the S/MWIR FPA are bands I4, M12, and M13, while bands I5, M14, M15, and M16 are situated on the LWIR FPA. Band M16 has two separate (and time delay integrated) arrays of 16 detectors. Only band M13 is dual gain-its low gain is primarily used for fire detection. Table 1 shows a summary of key design specifications for the VIIRS TEB. The VIIRS TEB are calibrated by referencing an onboard blackbody (BB). Key TEB performance parameters (i.e., detector responses (gains), dynamic ranges, noise equivalent differential temperatures (NEdT), and relative spectral responses) were carefully evaluated during prelaunch calibration and characterization using accurate and traceable sources and test equipment [8,9]. Postlaunch, the quality of the VIIRS TEB data and their downstream Remote Sens. 2021, 13, 3917 3 of 19 products relies vigorously on devoted efforts to regularly monitor instrument operations and onboard calibrator (OBC) functions. Moreover, it is equally important to evaluate and update the calibration parameters (when necessary) and employ new calibration strategies to address and correct on-orbit changes in the sensor characteristics [10,11]. Both the VIIRS TEB prelaunch and on-orbit calibrations are traceable to the standards maintained by the National Institute of Standards and Technology (NIST) [12]. This paper provides a comprehensive assessment of one decade of S-NPP VIIRS TEB on-orbit calibration and performance, which are critical to the SDR and EDR data quality. It serves as an update to previous reports [13,14]. Overall, the TEB on-orbit performance, in terms of detector responses, noise characteristics, and vicarious calibration, has remained extremely stable throughout the entire mission. In Section 2, an overview of the S-NPP VIIRS on-orbit operation and calibration activities, along with any changes to its calibration strategies over the last decade, is provided. Section 3 describes the VIIRS TEB calibration algorithm. The sensor's TEB on-orbit performance using responses derived from its OBC BB and select ground targets will be discussed in Section 4. To conclude, Section 5 provides a brief summary of this manuscript.

Instrument On-Orbit Operation and Calibration Activities
Following its launch in late October 2011, but prior to opening the nadir aperture door, S-NPP VIIRS underwent a sequence of sensor and OBC functional tests. Constructed by the VIS/NIR bands, the "first light" images were originated on 21 November 2011. Yet, observations from the S/MWIR and LWIR channels became technically reliable after the opening of the cryo-cooler door (18 January 2012). Afterwards, two more days were required for the S/MWIR and LWIR FPAs to reach operational temperatures. Important calibration activities executed during the sensor's early thorough calibration and validation stages involved routine operations and special calibration maneuvers. Table 2 itemizes some of the main events for the S-NPP VIIRS on-orbit operation and TEB calibration. Only the first event of each instrument anomaly and TEB routine calibration activity is registered in Table 2.
The onboard BB is the primary TEB calibration target. On a scan-by-scan basis, the BB temperature is derived based on measurements from six thermistors embedded on the backside of the BB substrate. The OBC BB uniformity-regarded as the standard deviation of the temperatures measured by the six thermistors-the requirement is 30 mK (under temperature-controlled or unpowered conditions) [12]. The S-NPP VIIRS BB temperature is nominally controlled at 292.7 K. Since the S-NPP launch (starting from February 2012) up until mid-2018, VIIRS BB warm-up and cool-down (WUCD) activities used to be scheduled on a quarterly basis. Afterwards, the BB WUCD operation frequency was changed to once annually. This was done to minimize the marginal thermal effects during each BB WUCD, which have an undesirable impact on data quality. Throughout a WUCD operation, the BB temperature sequences over discrete temperatures varying from instrument ambient temperature (~267 K) to 315 K. The process is performed recurrently to characterize and monitor changes in the TEB detector nonlinear responses while at the same time allowing for the examination of the detector noise characterization for the TEB over an ample temperature range. As of 31 December 2021, a total of 28 S-NPP VIIRS BB WUCD events have been completed thus far. Figure 2a shows a typical S-NPP VIIRS BB temperature profile during an entire WUCD cycle; the process takes roughly 2 days to complete. In Figure 2b, the BB uniformity for the same event is illustrated. It demonstrates how the BB operates within its design uniformity requirement (30 mK; red line) throughout the CD and most of the WU process. This behavior is consistent for all WUCD operations. The onboard BB is the primary TEB calibration target. On a scan-by-scan basis, the BB temperature is derived based on measurements from six thermistors embedded on the backside of the BB substrate. The OBC BB uniformity-regarded as the standard deviation of the temperatures measured by the six thermistors-the requirement is 30 mK (under temperature-controlled or unpowered conditions) [12]. The S-NPP VIIRS BB temperature is nominally controlled at 292.7 K. Since the S-NPP launch (starting from February 2012) up until mid-2018, VIIRS BB warm-up and cool-down (WUCD) activities used to be scheduled on a quarterly basis. Afterwards, the BB WUCD operation frequency was changed to once annually. This was done to minimize the marginal thermal effects during each BB WUCD, which have an undesirable impact on data quality. Throughout a WUCD operation, the BB temperature sequences over discrete temperatures varying from instrument ambient temperature (~267 K) to 315 K. The process is performed recurrently to characterize and monitor changes in the TEB detector nonlinear responses while at the same time allowing for the examination of the detector noise characterization for the TEB over an ample temperature range. As of 31 December 2021, a total of 28 S-NPP VIIRS BB WUCD events have been completed thus far. Figure 2a shows a typical S-NPP VIIRS BB temperature profile during an entire WUCD cycle; the process takes roughly 2 days to complete. In Figure 2b, the BB uniformity for the same event is illustrated. It demonstrates how the BB operates within its design uniformity requirement (30 mK; red line) throughout the CD and most of the WU process. This behavior is consistent for all WUCD operations.  In order to evaluate the BB short-term performance during normal operations, Figure 3a illustrates the average BB temperature over three orbits (approximately 5 h) on 1 February 2021; it varies with the orbital cycle at about an amplitude of approximately 27 mK (during daytime). This pattern was recorded during early mission [10,13] and remains in present time. Moreover, it was reported that Thermistors 3 and 6 (farthest from the EV port) drive these average BB temperature variations (with the possible cause attributed to the BB heating on one side only due to Earth illumination). Nonetheless, the OBC BB still operates within its design stability and uniformity requirements (see Figure 3b; daytime (σ ≈ 27 mK) and nighttime orbit values (σ ≈ 6 mK) are within the 30 mK requirement). Figure 4a illustrates the long-term average BB temperature; its mission-long time series demonstrates that it has been stable (within 25 mK), excluding WUCD processes and instrument anomalies. Early mission discontinuities (~15 mK) were due to the application of two different BB settings (driven by a WUCD event and an instrument anomaly) [13]. Subsequent operational processing was updated to avoid this artifact and the BB temperature has remained steady-within 10 mK-since the update.
(during daytime). This pattern was recorded during early mission [10,13] and rem present time. Moreover, it was reported that Thermistors 3 and 6 (farthest from t port) drive these average BB temperature variations (with the possible cause attribu the BB heating on one side only due to Earth illumination). Nonetheless, the OBC B operates within its design stability and uniformity requirements (see Figure 3b; d (σ ≈ 27 mK) and nighttime orbit values (σ ≈ 6 mK) are within the 30 mK require Figure 4a illustrates the long-term average BB temperature; its mission-long time demonstrates that it has been stable (within 25 mK), excluding WUCD processes a strument anomalies. Early mission discontinuities (~15 mK) were due to the appl of two different BB settings (driven by a WUCD event and an instrument anomaly Subsequent operational processing was updated to avoid this artifact and the BB te ature has remained steady-within 10 mK-since the update.  Also used in the thermal model and to inform about the instrument's health, 4b-d displays the mission-long time series of the other crucial telemetry points (s the cold FPAs temperatures). The non-BB temperatures used in the thermal mod shown in Figure 4b. All temperatures show yearly variations no larger than 3 K; th maximums correspond to the Earth at perihelion, while some other out-of-behavi ues are due to single board computer lock-up/petulant mode events. The electroni instrument temperatures, shown in Figure 4c, exhibited an increase of about 2.5 K ruary 2021; it varies with the orbital cycle at about an amplitude of approximately 27 mK (during daytime). This pattern was recorded during early mission [10,13] and remains in present time. Moreover, it was reported that Thermistors 3 and 6 (farthest from the EV port) drive these average BB temperature variations (with the possible cause attributed to the BB heating on one side only due to Earth illumination). Nonetheless, the OBC BB still operates within its design stability and uniformity requirements (see Figure 3b; daytime (σ ≈ 27 mK) and nighttime orbit values (σ ≈ 6 mK) are within the 30 mK requirement). Figure 4a illustrates the long-term average BB temperature; its mission-long time series demonstrates that it has been stable (within 25 mK), excluding WUCD processes and instrument anomalies. Early mission discontinuities (~15 mK) were due to the application of two different BB settings (driven by a WUCD event and an instrument anomaly) [13]. Subsequent operational processing was updated to avoid this artifact and the BB temperature has remained steady-within 10 mK-since the update.    Figure 4b. All temperatures show yearly variations no larger than 3 K; the local maximums correspond to the Earth at perihelion, while some other out-of-behavior values are due to single board computer lock-up/petulant mode events. The electronics and instrument temperatures, shown in Figure 4c, exhibited an increase of about 2.5 K over a decade. Figure 4d displays the cold FPAs temperatures. The LWIR FPA temperature has been rather stable (79.94 K) since it reached its operational temperature. The S/MWIR FPA temperature has risen approximately 40 mK (from~79.83 K to~79.87 K) over 10 years. Because the S/MWIR FPA temperature is passively controlled, as opposed to the LWIR FPA temperature, it shows more fluctuation [4]. Under nominal operation conditions, the impact due to such temperature variations on the calibrated L1B product is minimal, since the TEB calibration is updated on a scan-by-scan basis by referencing detector responses to the onboard BB. Table 3 lists all the acronyms used in Section 3.1 to facilitate readability. The VIIRS TEB apply a quadratic polynomial algorithm to retrieve EV scene spectral radiance, L EV , using background-subtracted detector responses (dn EV ) and prelaunch-derived calibration coefficients [4,5,7]. Additionally, these radiances are corrected for on-orbit detector degradation effects, instrument self-emission, and response-versus-scan-angle (RVS) effects,

Calibration Algorithm
where c i are the instrument temperature-dependent calibration coefficients derived during prelaunch characterization; RVS EV and RVS SV are the detector's RVS at the EV and space view (SV) HAM angle-of-incidence (AOI), correspondingly. RVS EV was derived from prelaunch measurements and RVS SV is equal to 1 since the SV is chosen as the normalization point of the RVS function. Moreover, ∆L bg is the instrument background emission radiance difference; ρ RTA is the RTA reflectivity; and L RTA and L H AM are the RTA and HAM emitted radiances, respectively. F is the on-orbit degradation or calibration scaling factor (F-factor), which is derived scan-by-scan from on-orbit BB measurements by comparing the known BB spectral radiance (L BB_thermal_model ) with that retrieved by the sensor (L BB_prelaunch ) using the prelaunch calibration coefficients for each band, HAM-side, and detector, where L BB_prelaunch = ∑ 2 i=0 c i ·dn i BB and dn BB is the background-subtracted BB detector response. Because the VIIRS TEB on-orbit calibration is performed referencing the instrument's onboard BB, L BB_thermal_model is modeled as the radiance difference between the BB and SV paths as: where L BB is the BB radiance and RVS BB is the RVS at the BB HAM AOI. The TEB BB radiance, L BB , is the sum of the emitted-from the thermal sources around the BB (i.e., BB shield, cavity, and telescope originated radiances)-and reflected (off the BB into the optical path) radiances. The BB emission (>99.6%) dominates L BB . The instrument background emission term, ∆L bg , in Equation (2) accounts for the RTA and HAM emissions that do not cancel out in the path differences (due to RVS differences at the two viewing angles) in Equations (1) and (4). Figure 5 shows a simplified diagram of the RTA as it rotates and views the different sectors used for the VIIRS TEB calibration. Using Equation (3), the TEB F-factors are estimated in alternating scans (for every HAM side) for all TEBexcept for the band M13 high gain (HG). The M13 HG F-factors are calculated sequentially every fourth scan (because L BB is too low to trend accurately, the M13 low gain (LG) F-factors are set to 1). A methodology to calibrate the M13 LG state using regular Moon observations, as the maximum surface temperature is around 390 K, from the SV port was proposed by McIntire et al. [15]. The general idea was that because dual gain bands alternate gain states every two scans during the calibration sector view, an opportunity presents itself to compare the M13 LG and HG responses over the intersecting portion of their dynamic ranges. McIntire et al. [15] describe that a 60 K range (about 40 low gain counts) over which to compare the calibration of both gain states is available (M13 HG saturation: 360 K; M13 LG signal becomes detectable at around 300 K). (3), the TEB F-factors are estimated in alternating scans (for every HAM side) for all TE-Bexcept for the band M13 high gain (HG). The M13 HG F-factors are calculated sequentially every fourth scan (because is too low to trend accurately, the M13 low gain (LG) F-factors are set to 1). A methodology to calibrate the M13 LG state using regular Moon observations, as the maximum surface temperature is around 390 K, from the SV port was proposed by McIntire et al. [15]. The general idea was that because dual gain bands alternate gain states every two scans during the calibration sector view, an opportunity presents itself to compare the M13 LG and HG responses over the intersecting portion of their dynamic ranges. McIntire et al. [15] describe that a 60 K range (about 40 low gain counts) over which to compare the calibration of both gain states is available (M13 HG saturation: 360 K; M13 LG signal becomes detectable at around 300 K).

On-Orbit Calibration Coefficients
The BB WUCD operation provides the basis to calibrate radiance over a large temperature range of nearly 50 K-nine temperature planes during warm-up; continual dataset during cool-down-over which the detector's on-orbit offset, linear response, and quadratic terms (c 0on−orbit , c 1on−orbit , and c 2on−orbit ) can be derived with a second-order polynomial fit to the incident detector radiance (L BB_thermal_model ) as it relates to the detector's response (dn BB ) as such: These on-orbit-derived calibration coefficients support the characterization and monitoring of changes in the TEB detector offset and nonlinear responses.

NEdT at Nominal BB and Typical Temperatures
The NEdT describes the instrument's sensitivity and denotes the (per sample) noise in the temperature space. Considered to be a noise estimate at the BB temperature level, the NEdT is derived over the 48-sample background-subtracted BB signal using the standard deviation [4,13]. NEdT can thus be calculated as, where NEdL is the Noise Equivalent radiance difference, which itself is dependent upon the BB detector response's signal-to-noise ratio, and ∂L/∂T is the derivative of the Planck function evaluated at the BB temperature. During a WUCD, because the BB temperature range is quite large, the NEdT can be assessed at each band's T typ , with which the instrument requirements are defined (see Table 1). Under normal operations, the NEdT is monitored at the nominal BB temperature (~292.7 K).

On-Orbit Performance
This section discusses the S-NPP VIIRS TEB on-orbit calibration performance by assessing the detector response trends, noise characterization, and WUCD calibration coefficients. It also discusses RVS stability assessments based on results from a designated spacecraft maneuver. Lastly, it presents results from vicarious calibration techniques to support the sensor's on-orbit calibration performance evaluation. Overall, the S-NPP VIIRS TEB calibration has been quite stable since launch. Note that all the results shown begin on 20 January 2012 (cold FPAs reached operational temperatures).

Detector Gain Trends
As presented in Section 3.1, the VIIRS TEB calibration uses a quadratic algorithm, and the detector gain (1/F-factor) is computed on a scan-by-scan basis using its response to the BB radiance. Using Equation (3), the TEB F-factors during nominal operations were calculated (data collected from two granules during every orbit's solar observation). Afterwards, the TEB gains were trended over on-orbit operations. Figure 6 shows the band-averaged results. All results shown henceforth are for HAM-side A only to avoid redundancy and prevent all figures from being overcrowded, since the HAM-side B results are essentially similar. All values are normalized to 20 January 2012. The degradations of the S-NPP VIIRS TEB have generally been small after a decade on-orbit. Band I5 exhibits the most degradation (around 3.1%), while the other LWIR bands show degradations no larger than 0.3%. Moreover, while bands M15 and M16 cover the same spectra as I5, these do not share the same level of degradation; this is possibly due to different electronics and sampling. All the LWIR bands display an annual cycle wherein local minimums are seen whenever the Earth is at perihelion. Moreover, a small jump early in the mission (mostly noticeable for the LWIR bands) can be attributed to a different BB temperature setting during a short period. These annual gain sequences are correlated to those in the temperatures used in the instrument thermal calibration model; though, since these temperatures are highly interrelated, it is difficult to pinpoint any specific telemetry temperature. All changes in the S-NPP MWIR detector responsivities are smaller than 1% after 10 years of sensor operations. Moreover, the MWIR bands had a minor gain decrease during early mission, it has since flattened. Additional discontinuities can be attributed to single board lock-up/petulant mode events and early mission instrument and spacecraft anomalies (see Table 2).
Remote Sens. 2021, 13, 3917 9 of 19 are highly interrelated, it is difficult to pinpoint any specific telemetry temperature. All changes in the S-NPP MWIR detector responsivities are smaller than 1% after 10 years of sensor operations. Moreover, the MWIR bands had a minor gain decrease during early mission, it has since flattened. Additional discontinuities can be attributed to single board lock-up/petulant mode events and early mission instrument and spacecraft anomalies (see Table 2). In the interest of evaluating the on-orbit short-term performance at the detector level, Figure 7 illustrates the bands I4 (detector 16, Product Order) and M14 (detector 8, Product Order) scan-by-scan and granule-averaged F-factors for one day's worth of data (1 February 2021). The F-factors were normalized to the first data point of the day and are only shown for HAM-side A. These two bands were selected to represent both cold FPAs and band types. Minor orbital variations (~0.1% and ~0.2% (or smaller) on a granule-averaged and scan-by-scan basis, respectively) can be observed for all TEB; all partially attributable to BB temperature variations throughout every orbit [13]. Only HAM-side A is illustrated. Note: Throughout this paper, the detectors are listed/numbered in Product Order. This means that the highest numbered detector (e.g., detector 16 for the M-bands and detector 32 for the I-bands) of a given scan is adjacent to detector 1 in the next scan. This is opposite to the prelaunch numbering convention (also known as Santa Barbara Remote Sensing (SBRS) Order), which is reversed. In the interest of evaluating the on-orbit short-term performance at the detector level, Figure 7 illustrates the bands I4 (detector 16, Product Order) and M14 (detector 8, Product Order) scan-by-scan and granule-averaged F-factors for one day's worth of data (1 February 2021). The F-factors were normalized to the first data point of the day and are only shown for HAM-side A. These two bands were selected to represent both cold FPAs and band types. Minor orbital variations (~0.1% and~0.2% (or smaller) on a granuleaveraged and scan-by-scan basis, respectively) can be observed for all TEB; all partially attributable to BB temperature variations throughout every orbit [13].
are highly interrelated, it is difficult to pinpoint any specific telemetry temperature. All changes in the S-NPP MWIR detector responsivities are smaller than 1% after 10 years of sensor operations. Moreover, the MWIR bands had a minor gain decrease during early mission, it has since flattened. Additional discontinuities can be attributed to single board lock-up/petulant mode events and early mission instrument and spacecraft anomalies (see Table 2). In the interest of evaluating the on-orbit short-term performance at the detector level, Figure 7 illustrates the bands I4 (detector 16, Product Order) and M14 (detector 8, Product Order) scan-by-scan and granule-averaged F-factors for one day's worth of data (1 February 2021). The F-factors were normalized to the first data point of the day and are only shown for HAM-side A. These two bands were selected to represent both cold FPAs and band types. Minor orbital variations (~0.1% and ~0.2% (or smaller) on a granule-averaged and scan-by-scan basis, respectively) can be observed for all TEB; all partially attributable to BB temperature variations throughout every orbit [13]. Only HAM-side A is illustrated. Note: Throughout this paper, the detectors are listed/numbered in Product Order. This means that the highest numbered detector (e.g., detector 16 for the M-bands and detector 32 for the I-bands) of a given scan is adjacent to detector 1 in the next scan. This is opposite to the prelaunch numbering convention (also known as Santa Barbara Remote Sensing (SBRS) Order), which is reversed. Only HAM-side A is illustrated. Note: Throughout this paper, the detectors are listed/numbered in Product Order. This means that the highest numbered detector (e.g., detector 16 for the M-bands and detector 32 for the I-bands) of a given scan is adjacent to detector 1 in the next scan. This is opposite to the prelaunch numbering convention (also known as Santa Barbara Remote Sensing (SBRS) Order), which is reversed.
While the S-NPP VIIRS TEB F-factors have been demonstrated to exhibit little variability and overall good stability during nominal operations, it is important to assess the on-orbit F-factors throughout non-nominal OBC BB temperatures (i.e., during WUCD operations). In Figure 8, the bands I4 (detector 16, Product Order) and M14 (detector 8, Product Order) F-factors, along with the BB temperature, are shown during the WUCD event performed from 15 March to 17 March 2021. All the F-factor values are granuleaveraged and only illustrated for HAM-side A. The M14 F-factors were shifted upward by 2.75% so that these can be in the same range as the I4 F-factors and in the same figure. Similar to Figure 7, these bands were selected to represent both cold FPAs and band types. Small F-factor anomalies and thus TEB calibration biases can be observed during the BB WUCD. F-factor transients at each WU step can be linked to each 5-K BB temperature step/interval during the WUCD process. Moreover, there is a noticeable increase in the F-factor during CD. This behavior is consistent for all WUCD operations. As shown in Figure 8, the band M14 F-factor decreases during WU and vice versa throughout the BB WUCD. Moreover, the WU period is shorter (~15 h) than the CD (24 h), and the latter's change in magnitude is larger as well. Because the F-factor performs as a scaling agent to the calibrated radiances, any flawed variations in F will transfer biases in the EV radiances. Other infrared bands exhibit similar behavior during the WUCD operation, albeit smaller in magnitude as the wavelengths become shorter. As mentioned previously, the 292.5 K nominal temperature is close to peak thermal emission for the long-wave infrared channels spectrally. Thus, any temperature change will have larger effects on the calibration of these bands. Cao et al. and Wang et al. have discussed this in detail [10,16], and Wang and Cao provided the methodology to attenuate the TEB calibration anomalies during the WUCD processes [14]. The NOAA IDPS SDRs have been updated to handle these changes since 25 July 2019. As discussed in detail by Wang and Cao, the F-factor computation will not change during nominal operations, and the updated TEB calibration algorithm supports three WUCD bias correction techniques: the Nominal-F method, the Ltrace method, and the Ltrace-2 method [14]. Moreover, the VIIRS emissive LUT was amended appropriately to arrange the WUCD bias correction technique used and to fit in the required correction coefficients. Their results indicate that implementing this correction efficiently decreases the F-factor anomalies for S-NPP VIIRS. More importantly, the daily-averaged F-factor anomalies subsided from 0.18% to 0.02% and -0.17% to −0.05% for M15 and M13, respectively, during the cool-down portion of the WUCD operation. The F-factor anomalies for the other TEB were also successfully mitigated. Lastly, Wang and Cao also provide a VIIRS-Cross-track Infrared Sounder (CrIS) brightness temperature (BT) bias time series comparison during the WUCD from March 2020 and show that the WUCD bias correction works properly for S-NPP (residual biases of~0.01 K). erations). In Figure 8, the bands I4 (detector 16, Product Order) and M14 (detector 8, Product Order) F-factors, along with the BB temperature, are shown during the WUCD event performed from 15 March to 17 March 2021. All the F-factor values are granule-averaged and only illustrated for HAM-side A. The M14 F-factors were shifted upward by 2.75% so that these can be in the same range as the I4 F-factors and in the same figure. Similar to Figure 7, these bands were selected to represent both cold FPAs and band types. Small Ffactor anomalies and thus TEB calibration biases can be observed during the BB WUCD. F-factor transients at each WU step can be linked to each 5-K BB temperature step/interval during the WUCD process. Moreover, there is a noticeable increase in the F-factor during CD. This behavior is consistent for all WUCD operations. As shown in Figure 8, the band M14 F-factor decreases during WU and vice versa throughout the BB WUCD. Moreover, the WU period is shorter (~15 h) than the CD (24 h), and the latter's change in magnitude is larger as well. Because the F-factor performs as a scaling agent to the calibrated radiances, any flawed variations in F will transfer biases in the EV radiances. Other infrared bands exhibit similar behavior during the WUCD operation, albeit smaller in magnitude as the wavelengths become shorter. As mentioned previously, the 292.5 K nominal temperature is close to peak thermal emission for the long-wave infrared channels spectrally. Thus, any temperature change will have larger effects on the calibration of these bands. Cao et al. and Wang et al. have discussed this in detail [10,16], and Wang and Cao provided the methodology to attenuate the TEB calibration anomalies during the WUCD processes [14]. The NOAA IDPS SDRs have been updated to handle these changes since 25 July 2019. As discussed in detail by Wang and Cao, the F-factor computation will not change during nominal operations, and the updated TEB calibration algorithm supports three WUCD bias correction techniques: the Nominal-F method, the Ltrace method, and the Ltrace-2 method [14]. Moreover, the VIIRS emissive LUT was amended appropriately to arrange the WUCD bias correction technique used and to fit in the required correction coefficients. Their results indicate that implementing this correction efficiently decreases the F-factor anomalies for S-NPP VIIRS. More importantly, the daily-averaged F-factor anomalies subsided from 0.18% to 0.02% and -0.17% to −0.05% for M15 and M13, respectively, during the cool-down portion of the WUCD operation. The F-factor anomalies for the other TEB were also successfully mitigated. Lastly, Wang and Cao also provide a VIIRS-Cross-track Infrared Sounder (CrIS) brightness temperature (BT) bias time series comparison during the WUCD from March 2020 and show that the WUCD bias correction works properly for S-NPP (residual biases of ~0.01 K).

Noise Characterization
As mentioned in Section 3.3, the noise performance for each TEB detector is characterized using the NEdT during nominal operating temperatures (~292.7 K) as well as the NEdT at T typ -derived using measurements from periodic BB WUCD operations. Figure 9 shows the band-averaged normalized NEdT trends (HAM-side A only) during nominal operations for all the S-NPP VIIRS TEB. These trends were plotted with the same granules used for Figure 6. All data points are normalized to the NEdT specification for each band. It can be observed that all TEB are well within the specification for the entirety of the mission. Bands M13, M15, and M16 are at about the halfway mark (~0.40-0.50) from their respective NEdT specifications, while bands I4, I5, and M12 have NEdT values no higher than 15% of their specification. Overall, the noise performance for all bands is quite stable. As discussed previously in Section 3.1, the small jump during early mission can be attributed to a different BB temperature setting. It should be noted that detector 2 in band I5, detector 16 in band M12, and detector 9 in band M16A (all in Product Order) have shown higher NEdT, albeit well within the specification, than other detectors in their respective bands and have been consistently like this throughout the mission. However, while well within their respective requirements, some of the stripe noise from these three bands propagates from the L1B clear-sky BT onto the Level 2 Sea Surface Temperature (SST) products. This striping has been shown to cause up to a 0.2 K SST difference in some cases [17]. Bouali and Ignatov [17] proposed an empirical destriping approach and demonstrated it to be feasible over ocean scenes via quantitative and qualitative analyses using an S-NPP VIIRS 3-day global dataset. Table 4 summarizes the band-averaged NEdT at T typ -derived from BB CD data for all S-NPP VIIRS TEB. The NEdT values are provided for three different years (2012, 2016, and 2021) in the sensor's mission. It also includes the prelaunch NEdT values as well as the NEdT requirement for each band. Similar to the results during nominal operations, the NEdT values for all bands are well within their specification at T typ . Moreover, all bands exhibit on-orbit NEdT slightly lower than their prelaunch-derived noise characterization. Lastly, the NEdT at T typ for the three years selected demonstrates that the noise performance has been quite stable over the instrument's mission as well.
terized using the NEdT during nominal operating temperatures (~292.7 K) as well as the NEdT at Ttyp-derived using measurements from periodic BB WUCD operations. Figure 9 shows the band-averaged normalized NEdT trends (HAM-side A only) during nominal operations for all the S-NPP VIIRS TEB. These trends were plotted with the same granules used for Figure 6. All data points are normalized to the NEdT specification for each band. It can be observed that all TEB are well within the specification for the entirety of the mission. Bands M13, M15, and M16 are at about the halfway mark (~0.40-0.50) from their respective NEdT specifications, while bands I4, I5, and M12 have NEdT values no higher than 15% of their specification. Overall, the noise performance for all bands is quite stable. As discussed previously in Section 3.1, the small jump during early mission can be attributed to a different BB temperature setting. It should be noted that detector 2 in band I5, detector 16 in band M12, and detector 9 in band M16A (all in Product Order) have shown higher NEdT, albeit well within the specification, than other detectors in their respective bands and have been consistently like this throughout the mission. However, while well within their respective requirements, some of the stripe noise from these three bands propagates from the L1B clear-sky BT onto the Level 2 Sea Surface Temperature (SST) products. This striping has been shown to cause up to a 0.2 K SST difference in some cases [17]. Bouali and Ignatov [17] proposed an empirical destriping approach and demonstrated it to be feasible over ocean scenes via quantitative and qualitative analyses using an S-NPP VIIRS 3-day global dataset. Table 4 summarizes the band-averaged NEdT at Ttyp-derived from BB CD data for all S-NPP VIIRS TEB. The NEdT values are provided for three different years (2012, 2016, and 2021) in the sensor's mission. It also includes the prelaunch NEdT values as well as the NEdT requirement for each band. Similar to the results during nominal operations, the NEdT values for all bands are well within their specification at Ttyp. Moreover, all bands exhibit on-orbit NEdT slightly lower than their prelaunch-derived noise characterization. Lastly, the NEdT at Ttyp for the three years selected demonstrates that the noise performance has been quite stable over the instrument's mission as well.

WUCD Calibration Coefficients
As briefly discussed in Section 3.2, the data from the WUCD operation provide a calibration radiance foundation over a broad 50 K temperature range that is independent from the prelaunch calibration source. Moreover, band, HAM-side, and detector-dependent 2nd order polynomial calibration coefficients (c 0on−orbit , c 1on−orbit , and c 2on−orbit ) can be derived using Equation (5). These WUCD-derived coefficients can be used to describe the VIIRS TEB on-orbit calibration nonlinearity over time. Furthermore, and as discussed previously, these can be utilized to lessen calibration biases during the WUCD cycle [14]. Figure 10 illustrates the band-averaged on-orbit-derived calibration coefficients (HAMside A only) for the S-NPP VIIRS TEB over the entire mission. Results are separated into WU and CD processes. Prelaunch coefficients are also included and defined by horizontal black lines. The offset (c 0on−orbit ) and nonlinear (c 2on−orbit ) coefficients derived from each WUCD are steady over the twenty-eight WUCD events thus far. These are also consistent with the LUT values derived prelaunch. There are some differences between the WU and CD calibration coefficients; the WU sequence generates larger nonlinear terms. For the S-NPP MWIR bands, the WUCD-derived linear coefficients (c 1on−orbit ) appear to be trending upward. This behavior is consistent with the slow degradations shown in the F-factors time series (see Figure 6). Band I5 also exhibits similar performance. on-orbit calibration updates. There is a major hurdle when trying to transfer prelaunch nonlinear calibration to on-orbit: three calibration coefficients (Equation (1)) and only one blackbody at nominal temperature. This results in one known parameter with two unknowns during on-orbit calibration. Ultimately, it was decided that an assumption was necessary to make this transfer possible. Thus, the degradation was represented by a single scaling (F) factor. This F-factor is one for the prelaunch test data. Cao et al. [16] proved this assumption to have its setbacks and have since then developed and implemented a methodology to remove F-factor fluctuations during BB WUCD events in the NOAA IDPS SDRs. The NASA SIPS L1B products have yet to implement it. Nonetheless, monitoring these coefficients is essential to track changes in sensor performance and characterization. It is important to mention that the F-factor efficiently adjusts for detector response changes at nominal BB temperatures. Over a broader range of scene temperatures, the F-factor is sensitive to gain changes, while changes in the nonlinear terms would affect the calibration at lower and higher temperatures than the nominal 292.7 K. Hence, these two on-orbit derived nonlinear terms are less accurate than their prelaunch counterparts, which were derived over a larger range of temperatures (190-345 K) using a NIST-calibrated BB source. Thus, any significant discrepancies between the prelaunch-and WUCD-derived coefficients can be evaluated using vicarious calibration results (Section 4.5) to determine if the prelaunch values should be updated. Currently, the prelaunch calibration coefficients, as opposed to the on-orbit-derived, are used to generate the NASA SIPS L1B product. Generally, successful calibration can be achieved this way because gain changes are corrected for using the calibration scaling F-factor (which is a function of detector response, the prelaunch calibration coefficients, instrument self-emission, and scan angle effects) discussed in Section 4.1. This is later confirmed in Section 4.5 through vicarious calibration. The NOAA IDPS SDR use a combination of the prelaunch coefficients (during nominal operations) and the WUCD-derived coefficients (during BB WUCD unsteady states; as discussed by Cao et al. and Wang et al. [10,14]).
As part of the VIIRS TEB calibration, it is quite challenging to accurately account for the instrument's nonlinear response. This is due to the fact that the sensor's nonlinearity cannot be simply determined after launch, and thus has to be transferred from prelaunch with an external blackbody calibration source in a laboratory. This is not uncommon for infrared radiometers. However, this philosophy (prelaunch to on-orbit calibration transfer) is quite different in the VIIRS radiometric calibration algorithm theoretical basis when compared to its heritage approach. As mentioned in Section 3.1, a nonlinear calibration curve is determined using the prelaunch calibration coefficients. In the case of heritage sensors such as MODIS, the linear term is dynamically calculated and applied on-orbit based on onboard BB measurements, while the nonlinear terms are updated frequently (when necessary). Hence, the VIIRS algorithm carries these three terms (c 0 ,c 1 ,c 2 ) from prelaunch to on-orbit assuming that the prelaunch laboratory test using the NIST blackbody calibration source provides the best sensor characterization and the derived prelaunch coefficients are reliable. Therefore, all postlaunch calibration is traceable to prelaunch. If indeed the instrument response remains identical from prelaunch to postlaunch, the on-orbit calibration can be achieved using Equation (1) with the prelaunch calibration coefficients. However, instrument degradation and different operating conditions may call for on-orbit calibration updates. There is a major hurdle when trying to transfer prelaunch nonlinear calibration to on-orbit: three calibration coefficients (Equation (1)) and only one blackbody at nominal temperature. This results in one known parameter with two unknowns during on-orbit calibration. Ultimately, it was decided that an assumption was necessary to make this transfer possible. Thus, the degradation was represented by a single scaling (F) factor. This F-factor is one for the prelaunch test data. Cao et al. [16] proved this assumption to have its setbacks and have since then developed and implemented a methodology to remove F-factor fluctuations during BB WUCD events in the NOAA IDPS SDRs. The NASA SIPS L1B products have yet to implement it. Nonetheless, monitoring these coefficients is essential to track changes in sensor performance and characterization.

RVS Validation during Designated Spacecraft Maneuver
A spacecraft deep space pitch maneuver is a designated postlaunch activity for RVS characterization. The S-NPP VIIRS pitch maneuver was designed and performed to validate its prelaunch TEB RVS characterization. It was completed on 20 February 2012 as part of the VIIRS rigorous calibration and validation stage to provide data for the on-orbit characterization of the TEB RVS [18]. In this maneuver, the spacecraft rotates around its y-axis and VIIRS gets a clear deep space view for a quarter of an orbit. When observing the deep space, the radiance perceived by the detectors comes from the sensor components' thermal emissions (including the HAM; whose emission varies with scan angle), while the EV external radiance is fundamentally zero. As such, the emissive radiation's relative change across the EV scan range (proportional to the change in detector's relative responses) can be examined to determine the TEB RVS. Data used for the TEB RVS characterization include the deep space and OBC BB detector responses. Furthermore, there is an underlying assumption that there are no polarization effects. This would impede HAM self-emission with scan angle variations from being effortlessly translated to RVS. During the prelaunch RVS derivation, an external unpolarized source was used to estimate the RVS [19]. Figure 11 illustrates the results published by Wu et al. [18,20]. These demonstrate good agreement between the on-orbit pitch maneuver-based and prelaunch RVS values-within 0.4% (close to the 0.3% requirement, and only exceeded by bands I4 and M13 at the edge-of-scan angles). It was inferred then that any HAM emission polarization effects are insignificant, and thus the RVS prelaunch characterization is valid.
A spacecraft deep space pitch maneuver is a designated postlaunch activity for RVS characterization. The S-NPP VIIRS pitch maneuver was designed and performed to validate its prelaunch TEB RVS characterization. It was completed on 20 February 2012 as part of the VIIRS rigorous calibration and validation stage to provide data for the on-orbit characterization of the TEB RVS [18]. In this maneuver, the spacecraft rotates around its y-axis and VIIRS gets a clear deep space view for a quarter of an orbit. When observing the deep space, the radiance perceived by the detectors comes from the sensor components' thermal emissions (including the HAM; whose emission varies with scan angle), while the EV external radiance is fundamentally zero. As such, the emissive radiation's relative change across the EV scan range (proportional to the change in detector's relative responses) can be examined to determine the TEB RVS. Data used for the TEB RVS characterization include the deep space and OBC BB detector responses. Furthermore, there is an underlying assumption that there are no polarization effects. This would impede HAM self-emission with scan angle variations from being effortlessly translated to RVS. During the prelaunch RVS derivation, an external unpolarized source was used to estimate the RVS [19]. Figure 11 illustrates the results published by Wu et al. [18,20]. These demonstrate good agreement between the on-orbit pitch maneuver-based and prelaunch RVS values-within 0.4% (close to the 0.3% requirement, and only exceeded by bands I4 and M13 at the edge-of-scan angles). It was inferred then that any HAM emission polarization effects are insignificant, and thus the RVS prelaunch characterization is valid.

Vicarious Calibration
Vicarious calibration practices make use of natural or artificial sites on the Earth's surface to calibrate sensors postlaunch. These are frequently used as an alternative to the common on-orbit calibration methods used to monitor sensor performance. Hence, Earth observations offer valuable information for calibration assessments. In this subsection, the on-orbit calibration performance of the S-NPP VIIRS TEB is evaluated using both an in situ ocean buoy and the Infrared Atmospheric Sounding Interferometer (IASI) instrument as references. The idea is to provide performance assessments via a ground site and through an inter-sensor comparison. Overall, the S-NPP VIIRS TEB show excellent performance over the instrument's mission through vicarious calibration.

VIIRS-In Situ Ocean Buoy Comparison
The in situ ocean buoy used belongs to the National Water Level Observation Network. The site is overseen and preserved by the NOAA's National Data Buoy Center

Vicarious Calibration
Vicarious calibration practices make use of natural or artificial sites on the Earth's surface to calibrate sensors postlaunch. These are frequently used as an alternative to the common on-orbit calibration methods used to monitor sensor performance. Hence, Earth observations offer valuable information for calibration assessments. In this subsection, the on-orbit calibration performance of the S-NPP VIIRS TEB is evaluated using both an in situ ocean buoy and the Infrared Atmospheric Sounding Interferometer (IASI) instrument as references. The idea is to provide performance assessments via a ground site and through an inter-sensor comparison. Overall, the S-NPP VIIRS TEB show excellent performance over the instrument's mission through vicarious calibration.

VIIRS-In Situ Ocean Buoy Comparison
The in situ ocean buoy used belongs to the National Water Level Observation Network. The site is overseen and preserved by the NOAA's National Data Buoy Center (NDBC) and has operated mostly uninterrupted for all the on-orbit S-NPP VIIRS lifetime. It is often referred to as Station SNDP5 (or #1619910), and is situated near the Sand Island, Midway Islands, HI (28.215 N, 177.361 W). The station consists of a 3-m discus buoy with a standard Self-Contained Ocean Observing Payload of sensors. It records SST measurements-at the Mean Lower Low Water level with high-accuracy (±0.002 • C) inductive sensors; these deliver SST data with excellent long-term stability (0.0001 K/y). Established in 2005, the data have been recorded every 6 min. Significant data gaps coinciding with the S-NPP VIIRS mission include 2 months in 2014 and another 2 months in 2019 for a total of 4 months.
S-NPP VIIRS Collection 1 L1B calibrated radiance, geolocation, and cloud mask products were downloaded for the 6-min granules collocated with the ocean buoy to extract the necessary radiance, location, and cloud cover information, respectively. Generally, there are two S-NPP overpasses over the site per day. Only nighttime observations were used to avoid the solar reflectance impacts on the MWIR bands. The methodology em-ployed for the VIIRS-buoy comparison is described in detail by Perez Diaz et al. [21,22]. In simplified terms, using the VIIRS L1B calibrated radiance and geolocation data, a 576 km 2 (24 km × 24 km) region-of-interest aligned around the site is found using the great-circle distance. Radiance data are converted to BT for all the TEB (via the band's spectral response function and Planck's law). All overcast pixels, by way of the VIIRS cloud mask, are sifted and the average BT of the residual pixels is computed. Data from all detectors for each TEB are utilized to compute the band-averaged BT. Lastly, all buoy SST data are temporally matched with the coincident VIIRS overpass. Figure 12 shows the S-NPP VIIRS bands M13 and M15 BT and in situ SST over the ocean buoy site near the Hawaii Islands. These two bands were selected to represent one band from each cold FPA and because their T typ are quite close to the average SST (see Table 1). The VIIRS and buoy temperatures are collocated in space and coincident in time. The brightness temperatures for the VIIRS bands were plotted as monthly averages. The in situ SST shown are nighttime observations only. All time series start on 20 January 2012. Data gaps in the SST observations are due to outages. It can be observed that both VIIRS bands exhibit similar temperatures in terms of magnitude and seasonal oscillations. Moreover, while there is an offset between the SST and VIIRS BT, it seems apparent that their overall trends are quite similar. This offset can be attributed to slight differences between the BT perceived by the VIIRS infrared channels and the water depth of the temperature sensors. To confirm this, Table 5 shows the VIIRS-buoy yearly change rates for each TEB over the sensor's mission. A detailed description and more thorough analyses of the results presented next is discussed by Pérez Díaz and Xiong [22]. Results indicate that all TEB display excellent on-orbit performance through vicarious calibration. All bands exhibit change rates within 0.02 K/y. Band M13 shows the largest temperature trend when compared to the buoy data with a total drift of 0.19 K after a decade in orbit. It is followed by band M12 with a downward drift of 0.13 K over 10 years of operations. These results demonstrate that the S-NPP VIIRS TEB continue to be well-calibrated throughout the mission at typical temperature scenes. In Section 4.5.2, the S-NPP VIIRS TEB on-orbit performance will be evaluated over a broader temperature range using an inter-sensor comparison.  The VIIRS TEB on-orbit performance was also evaluated by using IASI simultaneous nadir observations (SNOs) as reference. IASI, onboard the Meteorological Operational satellites (MetOp)--A and MetOp-C polar-orbiting meteorological satellites (altitude ~817 km), SNO data were used from May 2007 to April 2020, and from July 2019 to April 2021. The IASI spectral range spans from 3.6 μm to 15.5 μm with 8700 spectral samples covering

VIIRS-IASI Comparison
The VIIRS TEB on-orbit performance was also evaluated by using IASI simultaneous nadir observations (SNOs) as reference. IASI, onboard the Meteorological Operational satellites (MetOp)-A and MetOp-C polar-orbiting meteorological satellites (altitude~817 km), SNO data were used from May 2007 to April 2020, and from July 2019 to April 2021. The IASI spectral range spans from 3.6 µm to 15.5 µm with 8700 spectral samples covering nearly the whole VIIRS TEB spectrum. Hyperspectral IASI measurements were integrated to obtain simulated radiances for each VIIRS spectral band. All VIIRS pixels within each IASI SNO pixel (6-km radius) were aggregated to match the IASI pixel footprint. The radiances were then converted to BT values for cross-comparison [23]. MetOp-A shares overpasses with S-NPP at higher latitudes (73.0 • ± 1.1 • ). The number of IASI SNOs coincident with VIIRS retrievals are limited by each satellite's overpass time periodicity. There is a coincident S-NPP-IASI overpass event approximately every 50 days. There are roughly 270 VIIRS pixels within each IASI field-of-view.
In order to assess the calibration consistency over time, the BT difference (∆BT) between S-NPP VIIRS and IASI at different BT levels were trended over time. These three BT tiers were selected to evaluate the sensor performance at different scene temperatures, albeit with the caveat that the typical temperature (see Table 1) for each TEB cannot be reached using SNOs. Figure 13 displays the ∆BT time series for bands M13 and M15 at IASI-simulated M15 BT centered at 220 K, 240 K, and 260 K (described by different colors) with bin sizes of ±5 K. Empty symbols define IASI-A, while filled symbols represent IASI-C. It is observed that the ∆BT trends are stable over time (within 0.4 K; peak-to-peak variation), except for band M13 at 220 K (1.5 K peak-to-peak fluctuation).

Discussion
After close to a decade in-orbit, a solid S-NPP VIIRS TEB calibration methodology foundation has been established. Some challenges and room for improvement remain, with new challenges and anomalies more than likely to emerge as the instrument continues to operate past its designed lifetime (5 years). Moreover, improvements to the higher-level science products will surely make it necessary for the TEB to be accurately calibrated. There are remaining topics that need additional investigation as well as current issues that need resolution.
The S-NPP VIIRS TEB have dynamic ranges defined prelaunch. Oftentimes, when a channel's incident radiance reaches its maximum defined radiance level, the band's response is typically saturated. This is known as digital saturation. However, for some channels, as is the case for bands I4 and M12, the detector response decreases with increasing radiance. This phenomenon is defined as analog saturation or radiance rollover. This was defined during prelaunch testing and confirmed postlaunch. Generally, the concern is that such a band's response may seem valid and be recorded in the L1B when in reality, it is not. This artifact would then trickle down to the higher-level products. Band I4 is the most impacted by this artifact. It is crucial to identify the impacted pixels and flag them appropriately in the L1B product. VCST is currently investigating the analog saturation phenomenon and developing a mitigation strategy for future use.
The TEB calibration uncertainty is essential to evaluate the quality and accuracy of the TEB calibration. A detailed description of the uncertainty contributors to the TEB calibration was presented by McIntire et al. [24]. The preliminary results demonstrated that the leading terms in the MWIR bands EV uncertainties are the relative L BB and relative dn EV uncertainties, while the LWIR channels uncertainties are dominated by the c 0 , RVS, and dn EV relative uncertainties. Supplementary investigation and analyses are ongoing. The VIIRS TEB uncertainty algorithm is under refinement. Ultimately, the L1B product will have uncertainty information for each pixel once implemented.

Conclusions
S-NPP VIIRS has completed nearly 10 years of successful on-orbit operation and provided high-quality, long-term Earth data records to the scientific user community. Its OBC BB continues to function and is routinely operated to provide baseline calibration forand distinctively categorize-the S-NPP VIIRS TEB. The S-NPP VIIRS TEB calibration performance has remained stable throughout the mission. The gains of the S-NPP VIIRS detectors have experienced changes on orbit (within 3%), while the noise characterization performance has remained stable with most TEB detectors continuously meeting their design requirements. Moreover, the S-NPP VIIRS TEB on-orbit-derived calibration coefficients continue to exhibit and confirm this stability. The TEB RVS were validated via an on-orbit deep space maneuver and agree with the prelaunch values (within 0.4%). Lastly, calibration stability assessments using an in situ site and IASI SNOs demonstrate excellent calibration stability for all the S-NPP VIIRS TEB over a broad range of scenes and brightness temperatures. Uninterrupted stable calibration performance for any Earth observing instrument is key for the continuity and longevity of its data records. This presents a challenge towards the development of long-term climate data records. As evidenced in this paper, this requires devoted efforts to track (and possibly correct) on-orbit changes in sensor characteristics. The Joint Polar Satellite System VIIRS instruments intend to extend the VIIRS data records presented here for decades to come.
Author Contributions: Conceptualization, X.X.; methodology, C.L.P.D. and X.X.; formal analysis, C.L.P.D. and Y.L.; investigation, C.L.P.D. and Y.L.; writing-original draft preparation, C.L.P.D. and Y.L.; writing-review and editing, C.L.P.D., X.X and K.C. All authors have read and agreed to the published version of the manuscript.