Generating global products of LAI and FPAR from SNPP-VIIRS data: theoretical background and implementation

Leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) absorbed by vegetation have been successfully generated from the Moderate Resolution Imaging Spectroradiometer (MODIS) data since early 2000. As the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard, the Suomi National Polar-orbiting Partnership (SNPP) has inherited the scientific role of MODIS, and the development of a continuous, consistent, and well-characterized VIIRS LAI/FPAR data set is critical to continue the MODIS time series. In this paper, we build the radiative transfer-based VIIRS-specific lookup tables by achieving minimal difference with the MODIS data set and maximal spatial coverage of retrievals from the main algorithm. The theory of spectral invariants provides the configurable physical parameters, i.e., single scattering albedos (SSAs) that are optimized for VIIRS-specific characteristics. The effort finds a set of smaller red-band SSA and larger near-infraredband SSA for VIIRS compared with the MODIS heritage. The VIIRS LAI/FPAR is evaluated through comparisons with one year of MODIS product in terms of both spatial and temporal patterns. Further validation efforts are still necessary to ensure the product quality. Current results, however, imbue confidence in the VIIRS data set and suggest that the efforts described here meet the goal of achieving the operationally consistent multisensor LAI/FPAR data sets. Moreover, the strategies of parametric adjustment and LAI/FPAR evaluation applied to SNPP-VIIRS can also be employed to the subsequent Joint Polar Satellite System VIIRS or other instruments.Accepted manuscrip

Implications of whole-disc DSCOVR EPIC spectral observations for estimating Earth's spectral reflectivity based on low-earth-orbiting and geostationary observations

Earth’s reflectivity is among the key parameters of climate research. National Aeronautics and Space Administration (NASA)’s Earth Polychromatic Imaging Camera (EPIC) onboard National Oceanic and Atmospheric Administration (NOAA)’s Deep Space Climate Observatory (DSCOVR) spacecraft provides spectral reflectance of the entire sunlit Earth in the near backscattering direction every 65 to 110 min. Unlike EPIC, sensors onboard the Earth Orbiting Satellites (EOS) sample reflectance over swaths at a specific local solar time (LST) or over a fixed area. Such intrinsic sampling limits result in an apparent Earth’s reflectivity. We generated spectral reflectance over sampling areas using EPIC data. The difference between the EPIC and EOS estimates is an uncertainty in Earth’s reflectivity. We developed an Earth Reflector Type Index (ERTI) to discriminate between major Earth atmosphere components: clouds, cloud-free ocean, bare and vegetated land. Temporal variations in Earth’s reflectivity are mostly determined by clouds. The sampling area of EOS sensors may not be sufficient to represent cloud variability, resulting in biased estimates. Taking EPIC reflectivity as a reference, low-earth-orbiting-measurements at the sensor-specific LST tend to overestimate EPIC values by 0.8% to 8%. Biases in geostationary orbiting approximations due to a limited sampling area are between −0.7% and 12%. Analyses of ERTI-based Earth component reflectivity indicate that the disagreement between EPIC and EOS estimates depends on the sampling area, observation time and vary between −10% and 23%.The NASA/GSFC DSCOVR project is funded by NASA Earth Science Division. W. Song, G. Yan, and X. Mu were also supported by the key program of National Natural Science Foundation of China (NSFC; Grant No. 41331171). This research was conducted and completed during a 13-month research stay of the lead author in the Department of Earth and Environment, Boston University as a joint Ph.D. student, which was supported by the Chinese Scholarship Council (201606040098). DSCOVR EPIC L1B data were obtained from the NASA Langley Research Center Atmospheric Science Data Center. The authors would like to thank the editor who handled this paper and the two anonymous reviewers for providing helpful and constructive comments and suggestions that significantly helped us improve the quality of this paper. (NASA Earth Science Division; 41331171 - key program of National Natural Science Foundation of China (NSFC); 201606040098 - Chinese Scholarship Council)Accepted manuscrip

Satellite-Detected Contrasting Responses of Canopy Structure and Leaf Physiology to Drought

Disentangling drought impacts on plant photosynthesis is crucial for projecting future terrestrial carbon dynamics. We examined the separate responses of canopy structure and leaf physiology to an extreme summer drought that occurred in 2011 over Southwest China, where the weather was humid and radiation was the main growth-limiting factor. Canopy structure and leaf physiology were, respectively, represented by near-infrared reflectance of vegetation (NIRv) derived from MODIS data and leaf scale fluorescence yield (Φf) derived from both continuous SIF (CSIF) and global OCO-2 SIF (GOSIF). We detected contrasting responses of canopy structure and leaf physiology to drought with a 14.0% increase in NIRv, compared with 12.6 or 19.3% decreases in Φf from CSIF and GOSIF, respectively. The increase in structure resulted in a slight carbon change, due to water deficit-induced physiological constraints. The net ecosystem effect was a 7.5% (CSIF), 1.2% (GOSIF), and-2.96% (EC-LUE GPP) change in photosynthesis. Our study improves understanding of complex vegetation responses of plant photosynthesis to drought and may contribute to the reconciliation of contrasting observed directions in plant responses to drought in cloudy regions via remote sensing

Sun-Angle Effects on Remote-Sensing Phenology Observed and Modelled Using Himawari-8

Satellite remote sensing of vegetation at regional to global scales is undertaken at considerable variations in solar zenith angle (SZA) across space and time, yet the extent to which these SZA variations matter for the retrieval of phenology remains largely unknown. Here we examined the effect of seasonal and spatial variations in SZA on retrieving vegetation phenology from time series of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) across a study area in southeastern Australia encompassing forest, woodland, and grassland sites. The vegetation indices (VI) data span two years and are from the Advanced Himawari Imager (AHI), which is onboard the Japanese Himawari-8 geostationary satellite. The semi-empirical RossThick-LiSparse-Reciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model was inverted for each spectral band on a daily basis using 10-minute reflectances acquired by H-8 AHI at different sun-view geometries for each site. The inverted RTLSR model was then used to forward calculate surface reflectance at three constant SZAs (20°, 40°, 60°) and one seasonally varying SZA (local solar noon), all normalised to nadir view. Time series of NDVI and EVI adjusted to different SZAs at nadir view were then computed, from which phenological metrics such as start and end of growing season were retrieved. Results showed that NDVI sensitivity to SZA was on average nearly five times greater than EVI sensitivity. VI sensitivity to SZA also varied among sites (biome types) and phenological stages, with NDVI sensitivity being higher during the minimum greenness period than during the peak greenness period. Seasonal SZA variations altered the temporal profiles of both NDVI and EVI, with more pronounced differences in magnitude among NDVI time series normalised to different SZAs. When using VI time series that allowed SZA to vary at local solar noon, the uncertainties in estimating start, peak, end, and length of growing season introduced by local solar noon varying SZA VI time series, were 7.5, 3.7, 6.5, and 11.3 days for NDVI, and 10.4, 11.9, 6.5, and 8.4 days for EVI respectively, compared to VI time series normalised to a constant SZA. Furthermore, the stronger SZA dependency of NDVI compared with EVI, resulted in up to two times higher uncertainty in estimating annual integrated VI, a commonly used remote-sensing proxy for vegetation productivity. Since commonly used satellite products are not generally normalised to a constant sun-angle across space and time, future studies to assess the sun-angle effects on satellite applications in agriculture, ecology, environment, and carbon science are urgently needed. Measurements taken by new-generation geostationary (GEO) satellites offer an important opportunity to refine this assessment at finer temporal scales. In addition, studies are needed to evaluate the suitability of different BRDF models for normalising sun-angle across a broad spectrum of vegetation structure, phenological stages and geographic locations. Only through continuous investigations on how sun-angle variations affect spatiotemporal vegetation dynamics and what is the best strategy to deal with it, can we achieve a more quantitative remote sensing of true signals of vegetation change across the entire globe and through time

An Improved Approach to Monitoring Wheat Stripe Rust with Sun-Induced Chlorophyll Fluorescence

Sun-induced chlorophyll fluorescence (SIF) has shown potential in quantifying plant responses to environmental changes by which abiotic drivers are dominated. However, SIF is a mixed signal influenced by factors such as leaf physiology, canopy structure, and sun-sensor geometry. Whether the physiological information contained in SIF can better quantify crop disease stresses dominated by biological drivers, and clearly explain the physiological variability of stressed crops, has not yet been sufficiently explored. On this basis, we took winter wheat naturally infected with stripe rust as the research object and conducted a study on the responses of physiological signals and reflectivity spectrum signals to crop disease stress dominated by biological drivers, based on in situ canopy-scale and leaf-scale data. Physiological signals include SIF, SIFyield (normalized by absorbed photosynthetically active radiation), fluorescence yield (ΦF) retrieved by NIRvP (non-physiological components of canopy SIF) and relative fluorescence yield (ΦF-r) retrieved by near-infrared radiance of vegetation (NIRvR). Reflectance spectrum signals include normalized difference vegetation index (NDVI) and near-infrared reflectance of vegetation (NIRv). At the canopy scale, six signals reached extremely significant correlations (P < 0.001) with disease severity levels (SL) under comprehensive experimental conditions (SL without dividing the experimental samples) and light disease conditions (SL < 20%). The strongest correlation between NDVI and SL (R = 0.69) was observed under the comprehensive experimental conditions, followed by NIRv (R = 0.56), ΦF-r (R = 0.53) and SIF (R = 0.51), and the response of ΦF (R = 0.45) and SIFyield (R = 0.34) to SL was weak. Under lightly diseased conditions, ΦF-r (R = 0.62) showed the strongest response to disease, followed by SIFyield (R = 0.60), SIF (R = 0.56) and NIRv (R = 0.54). The weakest correlation was observed between ΦF and SL (R = 0.51), which also showed a result approximating NDVI (R = 0.52). In the case of a high level of crop disease severity, NDVI showed advantages in disease monitoring. In the early stage of crop diseases, which we pay more attention to, compared with SIF and reflectivity spectrum signals, ΦF-r estimated by the newly proposed ‘NIRvR approach’ (which uses SIF together with NIRvR (i.e., SIF/ NIRvR) as a substitute for ΦF) showed superior ability to monitor crop physiological stress, and was more sensitive to plant physiological variation. At the leaf scale, the response of SIF to SL was stronger than that of NDVI. These results validate the potential of ΦF-r estimated by the NIRvR approach to monitoring disease stress dominated by biological drivers, thus providing a new research avenue for quantifying crop responses to disease stress

Regional Leaf Area Index Retrieval Based on Remote Sensing: The Role of Radiative Transfer Model Selection

Physically-based approaches for estimating Leaf Area Index (LAI) using remote sensing data rely on radiative transfer (RT) models. Currently, many RT models are freely available, but determining the appropriate RT model for LAI retrieval is still problematic. This study aims to evaluate the necessity of RT model selection for LAI retrieval and to propose a retrieval methodology using different RT models for different vegetation types. Both actual experimental observations and RT model simulations were used to conduct the evaluation. Each of them includes needleleaf forests and croplands, which have contrasting structural attributes. The scattering from arbitrarily inclined leaves (SAIL) model and the four-scale model, which are 1D and 3D RT models, respectively, were used to simulate the synthetic test datasets. The experimental test dataset was established through two field campaigns conducted in the Heihe River Basin. The results show that the realistic representation of canopy structure in RT models is very important for LAI retrieval. If an unsuitable RT model is used, then the root mean squared error (RMSE) will increase from 0.43 to 0.60 in croplands and from 0.52 to 0.63 in forests. In addition, an RT model’s potential to retrieve LAI is limited by the availability of a priori information on RT model parameters. 3D RT models require more a priori information, which makes them have poorer generalization capability than 1D models. Therefore, physically-based retrieval algorithms should embed more than one RT model to account for the availability of a priori information and variations in structural attributes among different vegetation types

Global Land Cover Heterogeneity Characteristics at Moderate Resolution for Mixed Pixel Modeling and Inversion

Spatial heterogeneity is present in the land surface at every scale and is one of the key factors that introduces inherent uncertainty into simulations of land surface processes and parameter retrieval based on remotely sensed data. Because of a lack of understanding of the heterogeneous characteristics of global mixed pixels, few studies have focused on modeling and inversion algorithms in heterogeneous areas. This paper presents a parameterization scheme to describe land cover heterogeneity quantitatively by composition and boundary information based on high-resolution land cover products. Global heterogeneity features at the 1-km scale are extracted from the ‘GlobeLand30’ land cover dataset with a spatial resolution of 30 m. The composition analysis of global mixed pixels shows that only 35% of pixels over the land surface of Earth are covered by a single land cover type, namely, pure pixels, and only 25.8% are located in vegetated areas. Pixels mixed with water are more common than pixels mixed with any other non-vegetation type. The fragmentation analysis of typical biomes based on the boundary length shows that the savanna is the most heterogeneous biome, while the evergreen broadleaf forest is the least heterogeneous. Deciduous needleleaf forests are significantly affected by canopy height differences, while crop and grass biomes are less affected. Lastly, the strengths and limitations of the method and the application of the land cover heterogeneity characteristics extracted in this study are discussed

TCNIRv: Topographically corrected near-infrared reflectance of vegetation for tracking gross primary production over mountainous areas

The near-infrared reflectance of vegetation (NIRv) has been increasingly used as a proxy of gross primary production (GPP) across various temporal scales, ecosystems, and climate conditions. However, topography significantly distorts NIRv and GPP estimations over mountainous areas. We evaluated the topographic effects on NIRv and applied a path length correction (PLC) for improving its performance over mountainous areas. The proposed topographically corrected NIRv (referred to TCNIRv) was evaluated by multiple Landsat-8 operational land imager (OLI) images with concurrent in situ GPP measurements over the Lägeren mountainous forest area. TCNIRv reduced topographic effects in the original NIRv and it was comparable to the normalized difference vegetation index (NDVI) and the green normalized difference vegetation index (GNDVI), which are often deemed to be independent of topographic effects. In addition, TCNIRv better agreed with GPP than the other vegetation indices (VIs): coefficient of determination R2 = 0.90 and root mean square error RMSE = 1.40 g⋅ Cm −2⋅d −1 for TCNIRv compared to R2 = 0.71 and RMSE = 2.47 g⋅ Cm −2 ⋅d −1 for NIRv. The evaluation shows that TCNIRv is a reliable proxy of GPP, and because of its simplicity and physical soundness, it will facilitate vegetation monitoring over complex topography mountainous areas

A Cost-Constrained Sampling Strategy in Support of LAI Product Validation in Mountainous Areas

Increasing attention is being paid on leaf area index (LAI) retrieval in mountainous areas. Mountainous areas present extreme topographic variability, and are characterized by more spatial heterogeneity and inaccessibility compared with flat terrain. It is difficult to collect representative ground-truth measurements, and the validation of LAI in mountainous areas is still problematic. A cost-constrained sampling strategy (CSS) in support of LAI validation was presented in this study. To account for the influence of rugged terrain on implementation cost, a cost-objective function was incorporated to traditional conditioned Latin hypercube (CLH) sampling strategy. A case study in Hailuogou, Sichuan province, China was used to assess the efficiency of CSS. Normalized difference vegetation index (NDVI), land cover type, and slope were selected as auxiliary variables to present the variability of LAI in the study area. Results show that CSS can satisfactorily capture the variability across the site extent, while minimizing field efforts. One appealing feature of CSS is that the compromise between representativeness and implementation cost can be regulated according to actual surface heterogeneity and budget constraints, and this makes CSS flexible. Although the proposed method was only validated for the auxiliary variables rather than the LAI measurements, it serves as a starting point for establishing the locations of field plots and facilitates the preparation of field campaigns in mountainous areas