Search

The Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership (SNPP) satellite regularly performs on-orbit radiometric calibration of its reflective solar bands (RSBs), primarily through observations of an onboard sunlit solar diffuser (SD). The on-orbit change of the SD bidirectional reflectance distribution function (BRDF) value, quantified by a numerical factor called the H-factor, is determined by the onboard SD stability monitor. Our previous study showed that the H-factor is solar angle- and view direction-dependent. In this study, we determine the dependence of the H-factor on the detector SD view footprint location. We fit an empirical model to the NASA Collection 1 SNPP VIIRS Level 1B (L1B) spectral reflectance difference across the detectors in an RSB over uniform Earth scenes of the Libya 4 desert and deep convective clouds (DCCs). We apply the model predicted SD-positional-dependent H-factor to calibrate the RSBs. Under this new calibration scheme, the original unreal striping is removed from the homogeneous Libya 4 desert and the DCC images, as well as the original unreal striping from the Dunhuang desert image. The SD-positional-dependent H-factor has been used to calculate the SNPP VIIRS RSB radiometric correction factor for the NASA Collection 2.0 SNPP VIIRS L1B products.

Calibration of Terra and Aqua MODIS reflective solar bands (RSB) has evolved significantly since the launch of the first MODIS instrument on the Terra satellite more than 21 years ago. In NASA’s current Collection 6 and 6.1 Level 1B products (C6/C6.1 L1B), the RSB calibration algorithm continues to rely primarily on the onboard solar diffuser to calibrate the instrument gain. Lunar observations are used to track on-orbit changes in the response versus scan angle (RVS), and data from pseudo-invariant desert sites are used to apply adjustments to the gain and RVS calibration for select bands. The resulting reflectance products have in general shown a very stable performance. In recent years, some performance degradation has been noted for a few bands and algorithm changes have been tested to further improve the calibration accuracy for the upcoming Collection 7 (C7) L1B reprocess. In this paper, we present the MODIS RSB calibration improvements that will be included in C7. Major improvements include: applying polarization correction to the desert data before using it to generate RVS for Terra bands 8, 9, 3, and 10; using ocean scene data and an interband calibration approach to correct for long-term drift of Terra bands 11 and 12; applying an updated crosstalk correction to Terra SWIR bands over the entire mission; and using data from deep convective clouds in Terra SWIR band calibration, including the addition of time-dependent RVS for bands 5 and 26. All other minor calibration changes are also covered. Overall, the reflectance differences at nadir between C6.1 and C7 are within a few percent, though the differences increase in some cases at large scan angles. The Terra visible (3, 8-12) and SWIR bands (5-7, 26) have the most significant improvements. For all other Terra bands and all Aqua bands, the C7-C6.1 differences are mostly within 1%.

The design of the Moderate-Resolution Imaging Spectroradiometer (MODIS) instrument was driven by the scientific community’s desire to have near-daily global coverage at moderate resolution (~1km) with comprehensive spectral coverage from visible to long-wave infrared wavelengths. Since their launches in 1999 and 2002, respectively, the Terra and Aqua MODIS instruments have made continuous global observations and generated numerous data products to help users worldwide with their studies of the Earth’s system and its short- and long-term changes. The 20 reflective solar bands (RSBs) with wavelengths from 0.41 to $2.2~\mu \text{m}$ collect data at three nadir spatial resolutions: 250 m, 500 m, and 1 km. The solar diffuser (SD) coupled with the SD stability monitor (SDSM) provides a reflectance-based calibration on-orbit. In addition, lunar observations and response trends from pseudoinvariant desert sites are used to characterize the response versus scan-angle changes on-orbit. This paper provides a brief overview of MODIS RSB calibration algorithms, as implemented in the latest Level 1B version 6.1, operational activities, on-orbit performance, remaining challenges, and potential improvements. Results from the SD and SDSM measurements show a wavelength and mirror-side-dependent degradation in RSB responses, with the largest degradation at the shortest wavelengths, particularly for Terra MODIS. Aqua MODIS has experienced far less degradation of its optics and on-board calibrators compared with Terra MODIS, resulting in an overall better performance. With the exception of Aqua band 6, there have been no new noisy or inoperable detectors in the RSB of either instrument during postlaunch operations. As the instruments age and continue to endure the space environment, the detectors and the optical systems degrade. The challenges associated with incorporating these on-orbit changes to ensure a production of high-quality calibrated L1B data products are also discussed in this paper.

The Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the NOAA-20 (N20) satellite was launched on November 18, 2017. The N20 VIIRS reflective solar bands (RSBs) are calibrated on-orbit using a solar diffuser. An accurate on-orbit calibration is crucial to the high-quality downstream products facilitating atmosphere, ocean and land applications. In this study, the stability of the Level 1B (L1B) reflectance product is investigated using measurements over deep convective clouds (DCCs) for M-bands M1-M5, M7-M11, and I-bands I1-I3. The methodologies developed previously for Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) sensors and Suomi National Polar-orbiting Partnership (SNPP) VIIRS are extended and applied to the N20 RSB to derive DCC-based trends. The Collection 2 L1B data produced by NASA Land Science Investigator-led Processing Systems (SIPS) is used to evaluate the performance of the N20 VIIRS RSB calibration. At nadir, the reflectance trends for M1, M5, M8-M11, and I3 are insignificant compared to their corresponding variations (STDs) except for bands M2-M4, M7, and I1-I2, whose trends are larger than or equivalent to their STDs. The reflectance is relatively stable compared to their STDs for all the study RSBs at six aggregation zones across the entire scan angle range. Also discussed in this paper are the detector-to-detector differences and half-angle mirror side differences using the DCCs. Future applications using DCCs, which include an intercomparison with SNPP VIIRS, are also discussed.

An accurate on-orbit characterization of the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors onboard the Terra and Aqua satellites is essential to satisfy the requirements from the scientific community for high-quality operational and research products. In this study, NASA’s C6.1 L1B data are utilized to assess the stability of the calibrated top of atmosphere reflectance retrieved from individual detectors within a spectral band over deep convective clouds (DCCs) and desert sites. The results from these invariant Earth targets show detector-to-detector (D2D) differences in the L1B reflectance products for select MODIS reflective solar bands. For Aqua MODIS, the D2D differences have no significant change over time for the bands studied. However, the D2D differences exhibit an increase in magnitude over time for Terra MODIS bands 1, 3-7, and 26, which can induce noticeable striping and higher uncertainties in the downstream data products. The D2D stability assessment results over DCCs are consistent with those observed over desert sites. In addition, the view-angle dependence of the D2D differences is investigated by separating the frames along scan into five zones. The D2D assessment is performed for each of the five zones, and the results will further benefit any improvements in the calibration of a detector-dependent response versus scan angle that can be considered for future MODIS Level 1B Collections.

The Moderate-Resolution Imaging Spectroradiometer (MODIS) is the key instrument of the NASA’s Earth Observing System (EOS) Terra and Aqua missions, launched in December 1999 and May 2002, respectively. The Visible Infrared Imaging Radiometer Suite (VIIRS) expands the MODIS legacy, launched onboard the Suomi National Polar-orbiting Partnership (SNPP) satellite in October 2011 and NOAA20 satellite in November 2017, respectively. The MODIS and VIIRS sensors have a similar design with spectrally matched reflective solar spectral bands (RSB). Their on-board calibration components include a solar diffuser and a solar diffuser stability monitor for RSB, a V-grooved blackbody for the thermal emissive bands (TEB), and a space view as a background reference. This study evaluates the calibration consistency of the Terra and Aqua MODIS RSB based on the current Collection 6.1 L1B data products generated by the NASA MODAPS, which have used consistent calibration coefficient look up tables (LUT) for the entire data-record. In the case of SNPP and NOAA20 VIIRS, the latest L1B data produced by NASA Land SIPS are used. Several independent vicarious approaches are used to examine the stability and consistency of the at-sensor reflectance among MODIS and VIIRS instruments. Vicarious approaches include observations from simultaneous nadir overpasses (SNO), the Libya-4 desert and Dome C snow sites, as well as deep convective clouds (DCC). Impact of existing band spectral difference on the reflectance is corrected using hyperspectral observations provided by Europe Space Agency’s SCIAMACHY sensor. Results of this study provide comprehensive assessments of calibration performance, radiometric agreement and associated uncertainties.

The Moderate Resolution Imaging Spectroradiometer (MODIS) instruments on-board the Aqua and Terra spacecraft have provided valuable science data for the last 18 and 20 years, respectively. Each instrument is equipped with 36 spectral bands, 20 of which are reflective solar bands (RSBs). These bands cover a wavelength range of 0.4 - 2.2 μm and are calibrated on-orbit using several on-board calibrators (OBCs), such as a solar diffuser (SD) and a solar diffuser stability monitor (SDSM), along with regularly-scheduled lunar observations through the space view (SV) port. The gain (1/m1) and response-verses-scan angle (RVS) are updated on a near-monthly basis and act as the primary look-up-tables (LUTs) for the RSB calibration. A set of separate uncertainty LUTs for each of the RSBs are also delivered regularly and incorporated into the Level 1B (L1B) product to generate a pixel-level Uncertainty Index (UI). In addition to the gain, RVS and uncertainty, there are several other LUTs associated with the reflective bands that are either updated less frequently or remain static. The accuracy of both the forward-predicted and historical RSB LUTs, which are derived by the MODIS Characterization Support Team (MCST), is important in maintaining the quality and accuracy of the L1B and science products. To ensure a timely and accurate LUT update, MCST has established a comprehensive set of procedures. This paper provides an overview of the calibration process, along with the current LUT delivery process for the RSBs in Collection 6 (C6) and Collection 6.1 (C6.1). Improvements to be implemented in future collections are also discussed.

Terrestrial evapotranspiration (ET) for each plant functional type (PFT) is a key variable for linking the energy, water and carbon cycles of the atmosphere, hydrosphere and biosphere. Process-based algorithms have been widely used to estimate global terrestrial ET, yet each ET individual algorithm has exhibited large uncertainties. In this study, the support vector machine (SVM) method was introduced to improve global terrestrial ET estimation by integrating three process-based ET algorithms: MOD16, PT-JPL and SEMI-PM. At 200 FLUXNET flux tower sites, we evaluated the performance of the SVM method and others, including the Bayesian model averaging (BMA) method and the general regression neural networks (GRNNs) method together with three process-based ET algorithms. We found that the SVM method was superior to all other methods we evaluated. The validation results showed that compared with the individual algorithms, the SVM method driven by tower-specific (Modern Era Retrospective Analysis for Research and Applications, MERRA) meteorological data reduced the root mean square error (RMSE) by approximately 0.20 (0.15) mm/day for most forest sites and 0.30 (0.20) mm/day for most crop and grass sites and improved the squared correlation coefficient (R2) by approximately 0.10 (0.08) (95% confidence) for most flux tower sites. The water balance of basins and the global terrestrial ET calculation analysis also demonstrated that the regional and global estimates of the SVM-merged ET were reliable. The SVM method provides a powerful tool for improving global ET estimation to characterize the long-term spatiotemporal variations of the global terrestrial water budget. (Resume d'auteur)

For satellite sensor, the striping observed in images is typically associated with the relative multiple detector gain difference derived from the calibration. A method using deep convective cloud (DCC) measurements to assess the difference among detectors after calibration is proposed and demonstrated for select reflective solar bands (RSBs) of the Moderate Resolution Imaging Spectroradiometer (MODIS). Each detector of MODIS RSB is calibrated independently using a solar diffuser (SD). Although the SD is expected to accurately characterize detector response, the uncertainties associated with the SD degradation and characterization result in inadequacies in the estimation of each detector's gain. This work takes advantage of the DCC technique to assess detector uniformity and scan mirror side difference for RSB. The detector differences for Terra MODIS Collection 6 are less than 1% for bands 1, 3-5, and 18 and up to 2% for bands 6, 19, and 26. The largest difference is up to 4% for band 7. Most Aqua bands have detector differences less than 0.5% except bands 19 and 26 with up to 1.5%. Normally, large differences occur for edge detectors. The long-term trending shows seasonal oscillations in detector differences for some bands, which are correlated with the instrument temperature. The detector uniformities were evaluated for both unaggregated and aggregated detectors for MODIS band 1 and bands 3-7, and their consistencies are verified. The assessment results were validated by applying a direct correction to reflectance images. These assessments can lead to improvements to the calibration algorithm and therefore a reduction in striping observed in the calibrated imagery.

Terra MODIS has provided continuous global observations for science research and applications for more than 18 years. The MODIS Thermal emissive bands (TEB) radiometric calibration uses a quadratic function for instrument response. The calibration coefficients are updated using the response of an on-board blackbody (BB) in quarterly warm-up and cool-down (WUCD) events. As instrument degradation and electronic crosstalk of long-wave infrared (LWIR) bands 27 to 30 developed substantial issues, accurate calibration is crucial for a high-quality L1B product. The on-board BB WUCD temperature ranges from 270 K to 315 K and the derived nonlinear response has a relatively large uncertainty for the offset, especially for these LWIR bands, which affects the measurements of low brightness temperature (BT) scenes. In this study, the TEB radiometric calibration impact on the L1B product is assessed using selected cold targets and the measurements during regular lunar rolls. The cold targets include Antarctic Dome Concordia (Dome-C) and deep convective clouds (DCC) for the calibration assessment, focusing on bands 27 to 30. Dome-C area is covered with uniformly-distributed permanent snow, and the atmospheric effect is small and relatively constant. Usually the DCC is treated as an invariant earth target to evaluate the reflective solar band calibration. The DCC can also be treated as a stable target to assess the performance of TEB calibration. During a scheduled lunar observation event with a spacecraft roll maneuver to view the moon through the space view port, the instrument cavity provides a stable reference for calibration assessment. The long-term trending of BT measurements and the relative difference between scan mirror sides and detectors are used for the assessment of the calibration consistency and stability. The comparison of L1B products over the selected targets before and after the calibration coefficients update can be used to assess the impact of a calibration look-up table (LUT) update. This assessment is beneficial for future calibration algorithm and LUT update procedure improvements for enhancing the L1B product quality.

The Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP) satellite has operated successfully since its launch in October, 2011. The calibration using regular solar diffuser (SD) observations for the 14 reflective solar bands facilitates the generation of high-quality calibrated products, known as sensor data records (SDR). As SNPP VIIRS nears six years on-orbit, monitoring the on-orbit calibration performance using multiple techniques is vital. One such technique, using the deep convective clouds (DCC), has been widely used to monitor the on-orbit calibration performance of earth observing sensors, such as the MODIS instrument on Terra and Aqua spacecrafts. In this study, the DCC technique is utilized to evaluate the stability of SNPP VIIRS SDR reflectance product for 10 moderate resolution bands (M-bands, M1–M5 and M7– M11) and three imaging bands (I-bands, I1–I3). An empirical Bidirectional Reflectance Distribution Function (BRDF) correction for the long-term reflectance measurements over DCC is formulated and implemented. The fluctuation in the long-term trending is analyzed to evaluate the uncertainties. The BRDF-corrected reflectances over DCC are used to evaluate the stability of the response versus scan-angle (RVS) and its effects in the SDR products. The RVS effects are analyzed based on the difference in mode reflectances over DCC among six different frame zones. Results indicate that the RVS effects should be monitored for an improved VIIRS RSB calibration in the future.

Accurate estimation of the terrestrial latent heat flux (LE) for each plant functional type (PFT) at high spatial and temporal scales remains a major challenge. We developed a satellite-based hybrid algorithm to determine the Priestley–Taylor (PT) parameter for estimating global terrestrial LE across multiple biomes. The hybrid algorithm combines a simple empirical equation with physically based ecophysiological constraints to obtain the sum of the weighted ecophysiological constraints (f(e)) from satellite-based normalized difference vegetation index (NDVI) and ground-measured air temperature (Ta), relative humidity (RH), vapor pressure deficit (VPD) and LE for 2000 to 2009 provided by 240 globally distributed FLUXNET eddy covariance (ECOR) tower sites. Cross-validation analysis indicated that the optimization at a PFT level performed well with a RMSE of less than 0.15 and a R2 between 0.61 and 0.88 for estimated monthly f(e). Cross-validation analysis also revealed good performance of the hybrid-based PT method in estimating seasonal variability with a RMSE of the monthly LE varying from 4.3 W/m2 (for 6 deciduous needleleaf forest sites) to 18.1 W/m2 (for 34 crop sites) and with a R2 of more than 0.67. The algorithm's performance was also good for predicting among-site and inter-annual variability with a R2 of more than 0.78 and 0.70, respectively. We implemented the global terrestrial LE estimation from 2003 to 2005 for a spatial resolution of 0.05°by recalibrating the coefficients of the hybrid algorithm using Modern Era Retrospective Analysis for Research and Applications (MERRA) meteorological data, Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI product and ground-measured LE. This simple but accurate hybrid algorithm provides an alternative method for mapping global terrestrial LE, with a performance generally improved as compared to other satellite algorithms that are not calibrated with tower. The calibrated f(e) differs for different PFTs, and all driving forces of the algorithm can be acquired from satellite and meteorological observations.

This study compares six evapotranspiration ET products for Canada’s landmass, namely, eddy covariance EC measurements; surface water budget ET; remote sensing ET from MODIS; and land surface model (LSM) ET from the Community Land Model (CLM), the Ecological Assimilation of Land and Climate Observations (EALCO) model, and the Variable Infiltration Capacity model (VIC). The ET climatology over the Canadian landmass is characterized and the advantages and limitations of the datasets are discussed. The EC measurements have limited spatial coverage, making it difficult for model validations at the national scale. Water budget ET has the largest uncertainty because of data quality issues with precipitation in mountainous regions and in the north. MODIS ET shows relatively large uncertainty in cold seasons and sparsely vegetated regions. The LSM products cover the entire landmass and exhibit small differences in ET among them. Annual ET from the LSMs ranges from small negative values to over 600 mm across the landmass, with a countrywide average of 256 ± 15 mm. Seasonally, the countrywide average monthly ET varies from a low of about 3 mm in four winter months (November–February) to 67 ± 7 mm in July. The ET uncertainty is scale dependent. Larger regions tend to have smaller uncertainties because of the offset of positive and negative biases within the region. More observation networks and better quality controls are critical to improving ET estimates. Future techniques should also consider a hybrid approach that integrates strengths of the various ET products to help reduce uncertainties in ET estimation.

MODIS reflective solar bands are calibrated on-orbit using a solar diffuser and near-monthly lunar observations. To monitor the performance and effectiveness of the on-orbit calibrations, pseudo-invariant targets such as deep convective clouds (DCCs), Libya-4, and Dome-C are used to track the long-term stability of MODIS Level 1B product. However, the current MODIS operational DCC technique (DCCT) simply uses the criteria set for the 0.65-µm band. We optimize several critical DCCT parameters including the 11-µm IR-band Brightness Temperature (BT11) threshold for DCC identification, DCC core size and uniformity to help locate DCCs at convection centers, data collection time interval, and probability distribution function (PDF) bin increment for each channel. The mode reflectances corresponding to the PDF peaks are utilized as the DCC reflectances. Results show that the BT11 threshold and time interval are most critical for the Short Wave Infrared (SWIR) bands. The Bidirectional Reflectance Distribution Function model is most effective in reducing the DCC anisotropy for the visible channels. The uniformity filters and PDF bin size have minimal impacts on the visible channels and a larger impact on the SWIR bands. The newly optimized DCCT will be used for future evaluation of MODIS on-orbit calibration by MODIS Characterization Support Team.

An accurate on-orbit characterization of the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership satellite is essential to satisfy the requirements from NOAA, NASA, and the general research community for high-quality operational and research products. NASA’s Land Science Investigator-led Processing System (LSIPS) sensor data records (SDRs) are utilized to assess the stability of the calibrated top of atmosphere reflectance over the deep convective clouds (DCCs) and over the Libya-4 desert site. The results from DCCs and desert show detector-to-detector (D2D) reflectance differences in the VIIRS reflective solar bands, 10 moderate-resolution bands (M-bands, M1 to M5, M7 to M11), and three imagery resolution bands (I-bands, I1-I3). More importantly, the D2D differences exhibit an increase in magnitude over time for bands M1 to M5 and I1 to I2, which induces noticeable striping and higher uncertainties in the downstream data products. The largest D2D difference is observed in bands M1 and M2, with magnitudes greater than 1.8% in trends among the 16 detectors and standard deviations less than 0.12% over time. The D2D stability assessment results over DCCs are consistent with those observed over the Libya-4 desert site. To improve the quality of the calibrated SDR reflectance data, NASA’s VIIRS characterization support team has improved the calibration algorithm to incorporate corrections based on these DCC measurements to mitigate the striping (detector differences) observed in the current version of the SDR. These improvements are planned to be included in the next mission reprocessing of the LSIPS land products.

Evapotranspiration (ET) is a key component of terrestrial ecosystems because it links the hydrological, energy, and carbon cycles. Several satellite-based ET models have been developed for extrapolating local observations to regional and global scales, but recent studies have shown large model uncertainties in ET simulations. In this study, we compared eight ET models, including five empirical and three process-based models, with the objective of providing a reference for choosing and improving methods. The results showed that the eight models explained between 61 and 80% of the variability in ET at 23 eddy covariance towers in China and adjacent regions. The mean annual ET for all of China varied from 535 to 852 mm yr(-1) among the models. The interannual variability of yearly ET varied significantly between models during 1982-2009 because of different model structures and the dominant environmental factors employed. Our evaluation results showed that the parameters of the empirical methods may have different combination because the environmental factors of ET are not independent. Although the three process-based models showed high model performance across the validation Sites, there were substantial differences among them in the temporal and spatial patterns of ET, the dominant environment factors and the energy partitioning schemes. The disagreement among current ET models highlights the need for further improvements and validation, which can be achieved by investigating model structures and examining the ET component estimates and the critical model parameters. (C) 2013 Elsevier Inc. All rights reserved.

Regional drought and flooding from extreme climatic events are increasing in frequency and severity, with significant adverse ecosocial impacts. Detecting and monitoring drought at regional to global scales remains challenging, despite the availability of various drought indices and widespread availability of potentially synergistic global satellite observational records. The authors have developed a method to generate a near-real-time remotely sensed drought severity index (DSI) to monitor and detect drought globally at 1-km spatial resolution and regular 8-day, monthly, and annual frequencies. The new DSI integrates and exploits information from current operational satellite-based terrestrial evapo-transpiration (ET) and vegetation greenness index [normalized difference vegetation index (NDVI)] products, which are sensitive to vegetation water stress. Specifically, this approach determines the annual DSI departure from its normal (2000–11) using the remotely sensed ratio of ET to potential ET (PET) and NDVI. The DSI results were derived globally and captured documented major regional droughts over the last decade, including severe events in Europe (2003), the Amazon (2005 and 2010), and Russia (2010). The DSI corresponded favorably (correlation coefficient r = 0.43) with the precipitation-based Palmer drought severity index (PDSI), while both indices captured similar wetting and drying patterns. The DSI was also correlated with satellite-based vegetation net primary production (NPP) records, indicating that the combined use of these products may be useful for assessing water supply and ecosystem interactions, including drought impacts on crop yields and forest productivity. The remotely sensed global terrestrial DSI enhances capabilities for nearreal-time drought monitoring to assist decision makers in regional drought assessment and mitigation efforts, and without many of the constraints of more traditional drought monitoring methods.

The MODerate Resolution Imaging Spectroradiometer (MODIS) sensors onboard Terra and Aqua satellites are calibrated on-orbit with a solar diffuser (SD) for the reflective solar bands (RSB). The MODIS sensors are operating beyond their designed lifetime and hence present a major challenge to maintain the calibration accuracy. The degradation of the onboard SD is tracked by a solar diffuser stability monitor (SDSM) over a wavelength range from 0.41 to 0.94 μm. Therefore, any degradation of the SD beyond 0.94 μm cannot be captured by the SDSM. The uncharacterized degradation at wavelengths beyond this limit could adversely affect the Level 1B (L1B) product. To reduce the calibration uncertainties caused by the SD degradation, invariant Earth-scene targets are used to monitor and calibrate the MODIS L1B product. The use of deep convective clouds (DCCs) is one such method and particularly significant for the short-wave infrared (SWIR) bands in assessing their long-term calibration stability. In this study, we use the DCC technique to assess the performance of the Terra and Aqua MODIS Collection-6 L1B for RSB 1 3-7 , and 26, with spectral coverage from 0.47 to 2.13 μm. Results show relatively stable trends in Terra and Aqua MODIS reflectance for most bands. Careful attention needs to be paid to Aqua band 1, Terra bands 3 and 26 as their trends are larger than 1% during the study time period. We check the feasibility of using the DCC technique to assess the stability in MODIS bands 17-19. The assessment test on response versus scan angle (RVS) calibration shows substantial trend difference for Aqua band 1between different angles of incidence (AOIs). The DCC technique can be used to improve the RVS calibration in the future.

MODIS global evapotranspiration (ET) products by Mu et al. [Mu, Q., Heinsch, F. A., Zhao, M., Running, S. W. (2007). Development of a global evapotranspiration algorithm based on MODIS and global meteorology data. Remote Sensing of Environment, 111, 519–536. doi: 10.1016/j.rse.2007.04.015] are the first regular 1-km 2 land surface ET dataset for the 109.03 Million km 2 global vegetated land areas at an 8-day interval. In this study, we have further improved the ET algorithm in Mu et al. (2007a, hereafter called old algorithm) by 1) simplifying the calculation of vegetation cover fraction; 2) calculating ET as the sum of daytime and nighttime components; 3) adding soil heat flux calculation; 4) improving estimates of stomatal conductance, aerodynamic resistance and boundary layer resistance; 5) separating dry canopy surface from the wet; and 6) dividing soil surface into saturated wet surface and moist surface. We compared the improved algorithm with the old one both globally and locally at 46 eddy flux towers. The global annual total ET over the vegetated land surface is 62.8 × 10 3 km 3 , agrees very well with other reported estimates of 65.5 × 10 3 km 3 over the terrestrial land surface, which is much higher than 45.8 × 10 3 km 3 estimated with the old algorithm. For ET evaluation at eddy flux towers, the improved algorithm reduces mean absolute bias (MAE) of daily ET from 0.39 mm day −1 to 0.33 mm day −1 driven by tower meteorological data, and from 0.40 mm day −1 to 0.31 mm day −1 driven by GMAO data, a global meteorological reanalysis dataset. MAE values by the improved ET algorithm are 24.6% and 24.1% of the ET measured from towers, within the range (10–30%) of the reported uncertainties in ET measurements, implying an enhanced accuracy of the improved algorithm. Compared to the old algorithm, the improved algorithm increases the skill score with tower-driven ET estimates from 0.50 to 0.55, and from 0.46 to 0.53 with GMAO-driven ET. Based on these results, the improved ET algorithm has a better performance in generating global ET data products, providing critical information on global terrestrial water and energy cycles and environmental changes.

The reflective solar bands (RSB) of the Visible and Infrared Imaging Radiometer Suite (VIIRS) on Suomi National Polar-orbiting Partnership satellite are calibrated using solar diffuser and lunar observations. To evaluate the performance of the VIIRS RSB calibration, the deep convective clouds (DCC) are utilized as invariant targets to assess the stability of NASA’s VIIRS sensor data record (SDR) reflectance product for 10 moderate resolution bands (M-bands, M1 to M5 and M7 to M11) and three imagery resolution bands (I-bands, I1 to I3). All the frames along with each scan are divided into six aggregation zones for the scan angle. An empirical bidirectional reflectance distribution function (BRDF) correction is performed to SDR reflectances over DCC to reduce the scattering anisotropy in view and solar angles. The BRDF-corrected reflectances are used to assess the stability and to explore the residual RVS effects for selected bands. Results of this study show that after more than 6 years of on-orbit operation, the use of the prelaunch-based RVS still meets the radiometric requirement for the VIIRS RSB. Examination of the small but noticeable RVS residuals indicates they occur in a few shortest wavelength bands based on a time-dependent quadratic fit of the reflectance trends obtained at all scan angle zones. As the mission continues, continuing monitoring of on-orbit RVS stability and its effects on VIIRS SDR reflectance products are necessary with the DCC technique in the future VIIRS calibrations to ensure the quality of the RSB SDR products.

Model-data fusion offers considerable promise in remote sensing for improved state and parameter estimation particularly when applied to multi-sensor image products. This paper demonstrates the application of a ‘multiple constraints’ model-data fusion (MCMDF) scheme to integrating AMSR-E soil moisture content (SMC) and MODIS land surface temperature (LST) data products with a coupled biophysical model of surface moisture and energy budgets for savannas of northern Australia. The focus in this paper is on the methods, difficulties and error sources encountered in developing an MCMDF scheme and enhancements for future schemes. An important aspect of the MCMDF approach emphasized here is the identification of inconsistencies between model and data, and among data sets. The MCMDF scheme was able to identify that an inconsistency existed between AMSR-E SMC and LST data when combined with the coupled SEB-MRT model. For the example presented, an optimal fit to both remote sensing data sets together resulted in an 84% increase in predicted SMC and 0.06% increase for LST relative to the fit to each data set separately. That is the model predicted on average cooler LST's (∼ 1.7 K) and wetter SMC values (∼ 0.04 g cm − 3 ) than the satellite image products. In this instance we found that the AMSR-E SMC data on their own were poor constraints on the model. Incorporating LST data via the MCMDF scheme ameliorated deficiencies in the SMC data and resulted in enhanced characterization of the land surface soil moisture and energy balance based on comparison with the MODIS evapotranspiration (ET) product of Mu et al. [Mu, Q., Heinsch, F.A, Zhao, M. and Running, S.W. (in press), Development of a global evapotranspiration algorithm based on MODIS and global meteorology data, Remote Sensing of Environment .]. Canopy conductance, g C , and latent heat flux, λE , from the MODIS ET product were in good agreement with RMSEs for g C = 0.5 mm s − 1 and for λE = 18 W m − 2 , respectively. Differences were attributable to a greater canopy-to-air vapor pressure gradient in the MCMDF approach obtained from a more realistic partitioning of soil surface and canopy temperatures.

The objective of this research is to develop a global remote sensing evapotranspiration (ET) algorithm based on Cleugh et al.'s [Cleugh, H.A., R. Leuning, Q. Mu, S.W. Running (2007) Regional evaporation estimates from flux tower and MODIS satellite data. Remote Sensing of Environment 106, page 285–304- 2007 (doi: 10.1016/j.rse.2006.07.007).] Penman–Monteith based ET (RS-PM). Our algorithm considers both the surface energy partitioning process and environmental constraints on ET. We use ground-based meteorological observations and remote sensing data from the MODerate Resolution Imaging Spectroradiometer (MODIS) to estimate global ET by (1) adding vapor pressure deficit and minimum air temperature constraints on stomatal conductance; (2) using leaf area index as a scalar for estimating canopy conductance; (3) replacing the Normalized Difference Vegetation Index with the Enhanced Vegetation Index thereby also changing the equation for calculation of the vegetation cover fraction (FC); and (4) adding a calculation of soil evaporation to the previously proposed RS-PM method. We evaluate our algorithm using ET observations at 19 AmeriFlux eddy covariance flux towers. We calculated ET with both our Revised RS-PM algorithm and the RS-PM algorithm using Global Modeling and Assimilation Office (GMAO v. 4.0.0) meteorological data and compared the resulting ET estimates with observations. Results indicate that our Revised RS-PM algorithm substantially reduces the root mean square error (RMSE) of the 8-day latent heat flux (LE) averaged over the 19 towers from 64.6 W/m2 (RS-PM algorithm) to 27.3 W/m2 (Revised RS-PM) with tower meteorological data, and from 71.9 W/m2 to 29.5 W/m2 with GMAO meteorological data. The average LE bias of the tower-driven LE estimates to the LE observations changed from 39.9 W/m2 to − 5.8 W/m2 and from 48.2 W/m2 to − 1.3 W/m2 driven by GMAO data. The correlation coefficients increased slightly from 0.70 to 0.76 with the use of tower meteorological data. We then apply our Revised RS-PM algorithm to the globe using 0.05° MODIS remote sensing data and reanalysis meteorological data to obtain the annual global ET (MODIS ET) for 2001. As expected, the spatial pattern of the MODIS ET agrees well with that of the MODIS global terrestrial gross and net primary production (MOD17 GPP/NPP), with the highest ET over tropical forests and the lowest ET values in dry areas with short growing seasons. This MODIS ET product provides critical information on the regional and global water cycle and resulting environment changes.

Two models were evaluated for their ability to estimate land surface evaporation at 16-day intervals using MODIS remote sensing data and surface meteorology as inputs. The first was the aerodynamic resistance–surface energy balance model, and the second was the Penman–Monteith (P–M) equation, where the required surface conductance is estimated from remotely-sensed leaf area index. The models were tested using 3 years of evaporation and meteorological measurements from two contrasting Australian ecosystems, a cool temperate, evergreen Eucalyptus forest and a wet/dry, tropical savanna. The aerodynamic resistance–surface energy balance approach failed because small errors in the radiative surface temperature translate into large errors in sensible heat, and hence into estimates of evaporation. The P–M model adequately estimated the magnitude and seasonal variation in evaporation in both ecosystems (RMSE = 27 W m − 2 , R 2 = 0.74), demonstrating the validity of the proposed surface conductance algorithm. This, and the ability to constrain evaporation estimates via the energy balance, demonstrates the superiority of the P–M equation over the surface temperature-based model. There was no degradation in the performance of the P–M model when gridded meteorological data at coarser spatial (0.05°) and temporal (daily) resolution were substituted for locally-measured inputs. The P–M approach was used to generate a monthly evaporation climatology for Australia from 2001 to 2004 to demonstrate the potential of this approach for monitoring land surface evaporation and constructing monthly water budgets from 1-km to continental spatial scales.

Both MODIS and VIIRS use a solar diffuser (SD) to calibrate their reflective solar bands (RSB), covering wavelengths from 0.41 to 2.3 μm. On-orbit changes of the SD bi-directional reflectance factor (BRF) are tracked by an on-board solar diffuser stability monitor (SDSM). The current SDSM design only covers the spectral range from 0.41 to 0.93 μm. In general, the SD degradation is strongly wavelength-dependent with larger degradation occurring at shorter wavelengths, and the degradation in the SWIR region is expected to be extremely small. As each mission continues, however, the impact due to SD degradation at SWIR needs to be carefully examined and the correction if necessary should be applied in order to maintain the calibration quality. For Terra MODIS, alternative approaches have been developed and used to estimate the SD degradation for its band 5 at 1.24 μm and a time-dependent correction has already been applied to the current level 1B (L1B) collection 6 (C6). In this paper, we present different methodologies that can be used to examine the SD degradation and their applications for both Terra and Aqua MODIS and S-NPP VIIRS SWIR calibration. These methodologies include but not limited to the use of lunar observations, Pseudo Invariant Calibration Sites (PICS), and deep convective clouds (DCC). A brief description of relative approaches and their use is also provided in this paper.