An approach for developing Landsat-5 TM-based retrieval models of suspended particulate matter concentration with the assistance of MODIS

It is challenging to develop Landsat-5 TM (TM5) image-based retrieval models for estimating the suspended particulate matter concentration (C_SPM) in water when missing coincident ground C_SPM measurements. This study, with the Poyang Lake in China as a case study, proposed an approach for developing TM5-based C_SPM retrieval models with the assistance of moderate resolution imaging spectroradiometer (MODIS) images. After validation with an independent dataset, a cubic C_SPM retrieval model of 250 m MODIS red band was used to estimate the C_SPM values at 100 sampling points from the MODIS images (MODIS-based C_SPM) captured at three time periods. The MODIS-based C_SPM values at the time period with the largest C_SPM variation were combined with their coincident TM5 image reflectance for TM5-based model calibrations. The linear, quadratic, cubic, power and exponential models of MODIS-based C_SPM against TM5 single bands and their combinations were calibrated, respectively. Four best-fitting TM5-based C_SPM models were selected to retrieve the C_SPM values at 100 sampling points from the TM5 images (TM5-based C_SPM) at the other two time periods, and the coincident MODIS- and TM5-based C_SPM values were compared to assess TM5-based model performances. Model calibration results showed that the cubic and exponential models of TM5 red band (band 3) and red subtracting mid-infrared band (band 5) obtained the best fitting for estimating C_SPM from the TM5 image on 12 August 2005, and they explained 94–97% of the variation of MODIS-based C_SPM values with an estimated standard error of 6.617–8.457 mg/l. Model validations indicated that the exponential model of TM5 red band got the best result for estimating C_SPM from TM5 images when the MODIS-based C_SPM values were assumed as ground truths (correlation coefficient between MODIS- and TM5-based C_SPM values = 0.96, root mean square error = 4.60 mg/l). We concluded that the TM5-based C_SPM retrieval models could be developed with the assistance of MODIS, and the approach proposed in this study will be helpful for other researchers who also want to retrieve C_SPM from TM5 image archive but without coincident ground C_SPM measurements.     

Vegetation chlorophyll estimates in the Amazon from multi-angle MODIS observations and canopy reflectance model

As a preparatory study for future hyperspectral missions that can measure canopy chemistry, we introduce a novel approach to investigate whether multi-angle Moderate Resolution Imaging Spectroradiometer (MODIS) data can be used to generate a preliminary database with long-term estimates of chlorophyll. MODIS monthly chlorophyll estimates between 2000 and 2015, derived from a fully coupled canopy reflectance model (ProSAIL), were inspected for consistency with eddy covariance fluxes, tower-based hyperspectral images and chlorophyll measurements. MODIS chlorophyll estimates from the inverse model showed strong seasonal variations across two flux-tower sites in central and eastern Amazon. Marked increases in chlorophyll concentrations were observed during the early dry season. Remotely sensed chlorophyll concentrations were correlated to field measurements (r² = 0.73 and r² = 0.98) but the data deviated from the 1:1 line with root mean square errors (RMSE) ranging from 0.355 μg cm⁻² (Tapajós tower) to 0.470 μg cm⁻² (Manaus tower). The chlorophyll estimates were consistent with flux tower measurements of photosynthetically active radiation (PAR) and net ecosystem productivity (NEP). We also applied ProSAIL to mono-angle hyperspectral observations from a camera installed on a tower to scale modeled chlorophyll pigments to MODIS observations (r² = 0.73). Chlorophyll pigment concentrations (Chl_A+B) were correlated to changes in the amount of young and mature leaf area per month (0.59 ≤  r² ≤  0.64). Increases in MODIS observed Chl_A+B were preceded by increased PAR during the dry season (0.61 ≤ r² ≤  0.62) and followed by changes in net carbon uptake. We conclude that, at these two sites, changes in LAI, coupled with changes in leaf chlorophyll, are comparable with seasonality of plant productivity. Our results allowed the preliminary development of a 15-year time series of chlorophyll estimates over the Amazon to support canopy chemistry studies using future hyperspectral sensors.     

Testing estimation of water surface in Italian rice district from MODIS satellite data

Recent changes in rice crop management within Northern Italy rice district led to a reduction of seeding in flooding condition, which may have an impact on reservoir water management and on the animal and plant communities that depend on the flooded paddies. Therefore, monitoring and quantifying the spatial and temporal variability of water presence in paddy fields is becoming important. In this study we present a method to estimate dynamics of presence of standing water (i.e. fraction of flooded area) in rice fields using MODIS data. First, we produced high resolution water presence maps from Landsat by thresholding the Normalised Difference Flood Index (NDFI) made: we made it by comparing five Landsat 8 images with field-obtained information about rice field status and water presence. Using these data we developed an empirical model to estimate the flooding fraction of each MODIS cell. Finally we validated the MODIS-based flooding maps with both Landsat and ground information. Results showed a good predictability of water surface from Landsat (OA = 92%) and a robust usability of MODIS data to predict water fraction (R² = 0.73, EF = 0.57, RMSE = 0.13 at 1 × 1 km resolution). Analysis showed that the predictive ability of the model decreases with the greening up of rice, so we used NDVI to automatically discriminate estimations for inaccurate cells in order to provide the water maps with a reliability flag. Results demonstrate that it is possible to monitor water dynamics in rice paddies using moderate resolution multispectral satellite data. The achievement is a proof of concept for the analysis of MODIS archives to investigate irrigation dynamics in the last 15 years to retrieve information for ecological and hydrological studies.     

Evaluating MODIS snow products for modelling snowmelt runoff: Case study of the Rio Grande headwaters

Snow-covered area (SCA) is a key variable in the Snowmelt-Runoff Model (SRM) and in other models for simulating discharge from snowmelt. Landsat Thematic Mapper (TM), Enhanced Thematic Mapper (ETM + ) or Operational Land Imager (OLI) provide remotely sensed data at an appropriate spatial resolution for mapping SCA in small headwater basins, but the temporal resolution of the data is low and may not always provide sufficient cloud-free dates. The coarser spatial resolution Moderate Resolution Imaging Spectroradiometer (MODIS) offers better temporal resolution and in cloudy years, MODIS data offer the best alternative for mapping snow cover when finer spatial resolution data are unavailable. However, MODIS’ coarse spatial resolution (500 m) can obscure fine spatial patterning in snow cover and some MODIS products are not sensitive to end-of-season snow cover. In this study, we aimed to test MODIS snow products for use in simulating snowmelt runoff from smaller headwater basins by a) comparing maps of TM and MODIS-based SCA and b) determining how SRM streamflow simulations are changed by the different estimates of seasonal snow depletion. We compared gridded MODIS snow products (Collection 5 MOD10A1 fractional and binary SCA; SCA derived from Collection 6 MOD10A1 Normalised Difference Snow Index (NDSI) Snow Cover), and the MODIS Snow Covered-Area and Grain size retrieval (MODSCAG) canopy-corrected fractional SCA (SCA_MG), with reference SCA maps (SCA_REF) generated from binary classification of TM imagery. SCA_MG showed strong agreement with SCA_REF; excluding true negatives (where both methods agreed no snow was present) the median percent difference between SCA_REF and SCA_MG ranged between −2.4% and 4.7%. We simulated runoff for each of the four study years using SRM populated with and calibrated for snow depletion curves derived from SCA_REF. We then substituted in each of the MODIS-derived depletion curves. With efficiency coefficients ranging between 0.73 and 0.93, SRM simulation results from the SCA_MG runs yielded the best results of all the MODIS products and only slightly underestimated discharge volume (between 7 and 11% of measured annual discharge). SRM simulations that used SCA derived from Collection 6 NDSI Snow Cover also yielded promising results, with efficiency coefficients ranging between 0.73 and 0.91.In conclusion, we recommend that when simulating snowmelt runoff from small basins (<4000 km²) with SRM, we recommend that users select either canopy-corrected MODSCAG or create their own site-specific products from the Collection 6 MOD10A1 NDSI.     

Upscaling plot-scale soil respiration in winter wheat and summer maize rotation croplands in Julu County,North China

Soil respiration (R_s) data from 45 plots were used to estimate the spatial patterns of R_s during the peak growing seasons of winter wheat and summer maize in Julu County, North China, by combining satellite remote sensing data, field-measured data, and a support vector regression (SVR) model. The observed R_s values were well reproduced by the model at the plot scale, with a root-mean-square error (RMSE) of 0.31 μmol CO₂ m⁻² s⁻¹ and a coefficient of determination (R²) of 0.73. No significant difference was detected between the prediction accuracy of the SVR model for winter wheat and summer maize. With forcing from satellite remote sensing data and gridded soil property data, we used the SVR model to predict the spatial distributions of R_s during the peak growing seasons of winter wheat and summer maize rotation croplands in Julu County. The SVR model captured the spatial variations of R_s at the county scale. The satellite-derived enhanced vegetation index was found to be the most important input used to predict R_s. Removal of this variable caused an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.42 μmol CO₂ m⁻² s⁻¹. Soil properties such as soil organic carbon (SOC) content and soil bulk density (SBD) were the second most important factors. Their removal led to an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.37 μmol CO₂ m⁻² s⁻¹. The SVR model performed better than multiple regression in predicting spatial variations of R_s in winter wheat and summer maize rotation croplands, as shown by the comparison of the R² and RMSE values of the two algorithms. The spatial patterns of R_s are better captured using the SVR model than performing multiple regression, particularly for the relatively high and relatively low R_s values at the center and northeast study areas. Therefore, SVR shows promise for predicting spatial variations of R_s values on the basis of remotely sensed data and gridded soil property data at the county scale.     

Comparison of MODIS-based models for retrieving suspended particulate matter concentrations in Poyang Lake,China

Suspended particulate matter (SPM) is a key parameter describing water quality, and developing the retrieval model of SPM concentration (C_SPM) is fundamental for obtaining the spatiotemporal information of C_SPM and further for understanding, managing and protecting aquatic ecosystems. This study aimed to compare moderate resolution imaging spectroradiometer (MODIS)-based C_SPM retrieval models in order to find the optimal model for improving the C_SPM estimation in Poyang Lake. The C_SPM measurements on 27 September 2007 and their coincident MODIS Terra image were used to calibrate retrieval models with the least-squares technique. The C_SPM measurements on 31 August 2012 and the MODIS Terra image on 30 August 2012 were applied to validate the calibrated models, and the correlation coefficient (r) between the measured and estimated C_SPM values, the root mean square error (RMSE) and relative root mean square error (RRMSE) of estimation as well as the model bias evaluation result were compared to determine the optimal model for estimating the C_SPM values of Poyang Lake from MODIS images. Model calibration showed that, after two samples were removed, the exponential models of blue, green and red band, the linear model of infrared band, the cubic model of red band as well as the exponential model of red minus infrared band explained about 92%, 88%, 90%, 89%, 90% and 76% of the variation of C_SPM, respectively; while model validation indicated that, after removing two samples, the exponential models of blue and green band got biased C_SPM estimations, the agreement between the measured and estimated C_SPM values was not very high (r = <0.8) for the models with single red and infrared band, and the exponential model of red minus infrared band got the best result among all calibrated models (r = 0.87, RMSE = 22.1 mg/l, RRMSE = 52.8%). We concluded that the exponential model of red minus infrared band obtained stable C_SPM estimation and was the optimal model for C_SPM estimation in this study, and more independent datasets should be obtained to further validate our finding for improving the C_SPM estimation in Poyang Lake.     

Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data

A sufficient number of satellite acquisitions in a growing season are essential for deriving agronomic indicators, such as green leaf area index (GLAI), to be assimilated into crop models for crop productivity estimation. However, for most high resolution orbital optical satellites, it is often difficult to obtain images frequently due to their long revisit cycles and unfavorable weather conditions. Data fusion algorithms, such as the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), have been developed to generate synthetic data with high spatial and temporal resolution to address this issue. In this study, we evaluated the approach of assimilating GLAI into the Simple Algorithm for Yield Estimation model (SAFY) for winter wheat biomass estimation. GLAI was estimated using the two-band Enhanced Vegetation Index (EVI2) derived from data acquired by the Operational Land Imager (OLI) onboard the Landsat-8 and a fusion dataset generated by blending the Moderate-Resolution Imaging Spectroradiometer (MODIS) data and the OLI data using the STARFM and ESTARFM models. The fusion dataset had the temporal resolution of the MODIS data and the spatial resolution of the OLI data. Key parameters of the SAFY model were optimised through assimilation of the estimated GLAI into the crop model using the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm. A good agreement was achieved between the estimated and field measured biomass by assimilating the GLAI derived from the OLI data (GLAI_L) alone (R² = 0.77 and RMSE = 231 g m⁻²). Assimilation of GLAI derived from the fusion dataset (GLAI_F) resulted in a R² of 0.71 and RMSE of 193 g m⁻² while assimilating the combination of GLAI_L and GLAI_F led to further improvements (R² = 0.76 and RMSE = 176 g m⁻²). Our results demonstrated the potential of using the fusion algorithms to improve crop growth monitoring and crop productivity estimation when the number of high resolution remote sensing data acquisitions is limited.     

Inter- and intra-annual variations of clumping index derived from the MODIS BRDF product

Clumping index quantifies the level of foliage aggregation, relative to a random distribution, and is a key structural parameter of plant canopies and is widely used in ecological and meteorological models. In this study, the inter- and intra-annual variations in clumping index values, derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) BRDF product, are investigated at six forest sites, including conifer forests, a mixed deciduous forest and an oak-savanna system. We find that the clumping index displays large seasonal variation, particularly for the deciduous sites, with the magnitude in clumping index values at each site comparable on an intra-annual basis, and the seasonality of clumping index well captured after noise removal. For broadleaved and mixed forest sites, minimum clumping index values are usually found during the season when leaf area index is at its maximum. The magnitude of MODIS clumping index is validated by ground data collected from 17 sites. Validation shows that the MODIS clumping index can explain 75% of variance in measured values (bias = 0.03 and rmse = 0.08), although with a narrower amplitude in variation. This study suggests that the MODIS BRDF product has the potential to produce good seasonal trajectories of clumping index values, but with an improved estimation of background reflectance.     

Intercomparison and validation of MODIS and GLASS leaf area index (LAI) products over mountain areas: A case study in southwestern China

The validation study of leaf area index (LAI) products over rugged surfaces not only gives additional insights into data quality of LAI products, but deepens understanding of uncertainties regarding land surface process models depended on LAI data over complex terrain. This study evaluated the performance of MODIS and GLASS LAI products using the intercomparison and direct validation methods over southwestern China. The spatio-temporal consistencies, such as the spatial distributions of LAI products and their statistical relationship as a function of topographic indices, time, and vegetation types, respectively, were investigated through intercomparison between MODIS and GLASS products during the period 2011–2013. The accuracies and change ranges of these two products were evaluated against available LAI reference maps over 10 sampling regions which standed for typical vegetation types and topographic gradients in southwestern China.The results show that GLASS LAI exhibits higher percentage of good quality data (i.e. successful retrievals) and smoother temporal profiles than MODIS LAI. The percentage of successful retrievals for MODIS and GLASS is vulnerable to topographic indices, especially to relief amplitude. Besides, the two products do not capture seasonal dynamics of crop, especially in spring over heterogeneously hilly regions. The yearly mean LAI differences between MODIS and GLASS are within ±0.5 for 64.70% of the total retrieval pixels over southwestern China. The spatial distribution of mean differences and temporal profiles of these two products are inclined to be dominated by vegetation types other than topographic indices. The spatial and temporal consistency of these two products is good over most area of grasses/cereal crops; however, it is poor for evergreen broadleaf forest. MODIS presents more reliable change range of LAI than GLASS through comparison with fine resolution reference maps over most of sampling regions. The accuracies of direct validation are obtained for GLASS LAI (r = 0.35, RMSE = 1.72, mean bias = −0.71) and MODIS LAI (r = 0.49, RMSE = 1.75, mean bias = −0.67). GLASS performs similarly to MODIS, but may be marginally inferior to MODIS based on our direct validation results. The validation experience demonstrates the necessity and importance of topographic consideration for LAI estimation over mountain areas. Considerable attention will be paid to the improvements of surface reflectance, retrieval algorithm and land cover types so as to enhance the quality of LAI products in topographically complex terrain.     

Assessment of leaf carotenoids content with a new carotenoid index: Development and validation on experimental and model data

Leaf carotenoids content (LCar) is an important indicator of plant physiological status. Accurate estimation of LCar provides valuable insight into early detection of stress in vegetation. With spectroscopy techniques, a semi-empirical approach based on spectral indices was extensively used for carotenoids content estimation. However, established spectral indices for carotenoids that generally rely on limited measured data, might lack predictive accuracy for carotenoids estimation in various species and at different growth stages. In this study, we propose a new carotenoid index (CARI) for LCar assessment based on a large synthetic dataset simulated from the leaf radiative transfer model PROSPECT-5, and evaluate its capability with both simulated data from PROSPECT-5 and 4SAIL and extensive experimental datasets: the ANGERS dataset and experimental data acquired in field experiments in China in 2004. Results show that CARI was the index most linearly correlated with carotenoids content at the leaf level using a synthetic dataset (R² = 0.943, RMSE = 1.196 μg/cm²), compared with published spectral indices. Cross-validation results with CARI using ANGERS data achieved quite an accurate estimation (R² = 0.545, RMSE = 3.413 μg/cm²), though the RBRI performed as the best index (R² = 0.727, RMSE = 2.640 μg/cm²). CARI also showed good accuracy (R² = 0.639, RMSE = 1.520 μg/cm²) for LCar assessment with leaf level field survey data, though PRI performed better (R² = 0.710, RMSE = 1.369 μg/cm²). Whereas RBRI, PRI and other assessed spectral indices showed a good performance for a given dataset, overall their estimation accuracy was not consistent across all datasets used in this study. Conversely CARI was more robust showing good results in all datasets. Further assessment of LCar with simulated and measured canopy reflectance data indicated that CARI might not be very sensitive to LCar changes at low leaf area index (LAI) value, and in these conditions soil moisture influenced the LCar retrieval accuracy.     

Soil respiration mapped by exclusively use of MODIS data for forest landscapes of Saskatchewan,Canada

Soil respiration (R_s) is of great importance to the global carbon balance. Remote sensing of R_s is challenging because of (1) the lack of long-term R_s data for model development and (2) limited knowledge of using satellite-based products to estimate R_s. Using 8-years (2002–2009) of continuous R_s measurements with nonsteady-state automated chamber systems at a Canadian boreal black spruce stand (SK-OBS), we found that R_s was strongly correlated with the product of the normalized difference vegetation index (NDVI) and the nighttime land surface temperature (LSTn) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery. The coefficients of the linear regression equation of this correlation between R_s and NDVI × LSTn could be further calibrated using the MODIS leaf area index (LAI) product, resulting in an algorithm that is driven solely by remote sensing observations. Modeled R_s closely tracked the seasonal patterns of measured R_s and explained 74–92% of the variance in R_s with a root mean square error (RMSE) less than 1.0 g C/m²/d. Further validation of the model from SK-OBS site at another two independent sites (SK-OA and SK-OJP, old aspen and old jack pine, respectively) showed that the algorithm can produce good estimates of R_s with an overall R² of 0.78 (p < 0.001) for data of these two sites. Consequently, we mapped R_s of forest landscapes of Saskatchewan using entirely MODIS observations for 2003 and spatial and temporal patterns of R_s were well modeled. These results point to a strong relationship between the soil respiratory process and canopy photosynthesis as indicated from the greenness index (i.e., NDVI), thereby implying the potential of remote sensing data for detecting variations in R_s. A combination of both biological and environmental variables estimated from remote sensing in this analysis may be valuable in future investigations of spatial and temporal characteristics of R_s.     

Retrieval of land surface temperature from the Kalpana-1 VHRR data using a single-channel algorithm and its validation over western India

Indian geostationary satellite Kalpana-1 (K1) offers a potential to capture the diurnal cycle of land surface temperature (LST) through thermal infrared channel (10.5–12.5 μm) observations of the Very High Resolution Radiometer (VHRR) sensor. A study was carried out to retrieve LST by adapting a generalized single-channel (SC) algorithm (Jiménez-Muñoz and Sobrino, 2003) for the VHRR sensor over India. The basis of SC algorithm depends on the concept of Atmospheric Functions (AFs) that are dependent on transmissivity, upwelling and downwelling radiances of the atmosphere. In the present study AFs were computed for the VHRR sensor through the MODTRAN simulations based upon varying atmospheric and surface inputs. The AFs were fitted with the atmospheric columnar water vapour content and a set of coefficients was derived for LST retrieval. The K1-LST derived with the SC algorithm was validated with (a) in situ measurements at two sites located in western parts of India and (b) the MODIS LST products. Comparison of K1-LST with the in situ measurements demonstrated that SC algorithm was successful in capturing the prominent diurnal variations of 283–332 K in the LST at desert and agriculture experimental sites with a rmse of 1.6 K and 2.7 K, respectively. Inter comparison of K1-LST and MODIS LST showed a reasonable agreement between these two retrievals up to LST of 300 K, however a cold bias up to 7.9 K was observed in MODIS LST for higher LST values (310–330 K) over the hot desert region.     

Dry-season vegetation mass and cover fraction from SWIR1.6 and SWIR2.1 band ratio: Ground-radiometer and MODIS data in the Sahel

The potential of the short-wave infrared (SWIR) bands to detect dry-season vegetation mass and cover fraction is investigated with ground radiometry and MODIS data, confronted to vegetation data collected in rangeland and cropland sites in the Sahel (Senegal, Niger, Mali). The ratio of the 1.6 and 2.1 μm bands (called STI) acquired with a ground radiometer proved well suited for grassland mass estimation up to 2500 kg/ha with a linear relation (r² = 0.89). A curvilinear regression is accurate for masses ranging up to 3500 kg/ha. STI proved also well suited to retrieve vegetation cover fraction in crop fields, fallows and rangelands. Such dry-season monitoring, with either ground or satellite data, has important applications for forage, erosion risk and fire risk assessment in semi-arid areas.     

Retrieval of leaf water content spanning the visible to thermal infrared spectra

The objective of this study was to investigate the entire spectra (from visible to the thermal infrared; 0.390–14.0 μm) to retrieve leaf water content in a consistent manner. Narrow-band spectral indices (calculated from all possible two band combinations) and a partial least square regression (PLSR) were used to assess the strength of each spectral region. The coefficient of determination (R²) and root mean square error (RMSE) were used to report the prediction accuracy of spectral indices and PLSR models. In the visible-near infrared and shortwave infrared (VNIR–SWIR), the most accurate spectral index yielded R² of 0.89 and RMSE of 7.60%, whereas in the mid infrared (MIR) the highest R² was 0.93 and RMSE of 5.97%. Leaf water content was poorly predicted using two-band indices developed from the thermal infrared (R² = 0.33). The most accurate PLSR model resulted from MIR reflectance spectra (R² = 0.96, RMSE = 4.74% and RMSE cross validation RMSE_CV = 6.17%) followed by VNIR–SWIR reflectance spectra (R² = 0.91, RMSE = 6.90% and RMSE_CV = 7.32%). Using thermal infrared (TIR) spectra, the PLSR model yielded a moderate retrieval accuracy (R² = 0.67, RMSE = 13.27% and RMSE_CV = 16.39%). This study demonstrated that the mid infrared (MIR) and shortwave infrared (SWIR) domains were the most sensitive spectral region for the retrieval of leaf water content.     

Asian Dust Storm Monitoring Combining Terra and Aqua MODIS SRB Measurements

Sand and dust storms (SDSs), which present environmental risks and affect the regional climate, have been worsened in the East Asian regions over the last decade. Monitoring SDS from space using satellite remote sensing (RS) has become one of the most important issues in this field. At present, satellite RS of SDS is limited to using true-color images or aerosol optical thickness (AOT), or a new algorithm called "Deep Blue". Using current existing approaches makes it difficult to identify SDS from clouds. The authors have detected SDS by combining Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) solar reflectance band (SRB) measurements. Based on the dust spectral characteristic, this letter proposes a normalized difference dust index (NDDI) using MODIS reflectance measurements and applies it to the Asian SDS cases. The simple NDDI index is found to be able to identify SDS and clouds easily. The results suggest that NDDI could be used to detect SDS over bright surfaces where the MODIS AOT product is not available     

A visible band index for remote sensing leaf chlorophyll content at the canopy scale

Leaf chlorophyll content is an important variable for agricultural remote sensing because of its close relationship to leaf nitrogen content. The triangular greenness index (TGI) was developed based on the area of a triangle surrounding the spectral features of chlorophyll with points at (670 nm, R₆₇₀), (550 nm, R₅₅₀), and (480 nm, R₄₈₀), where R_λ is the spectral reflectance at wavelengths of 670, 550 and 480, respectively. The equation is TGI = −0.5[(670 − 480)(R₆₇₀ − R₅₅₀) − (670 − 550)(R₆₇₀ − R₄₈₀)]. In 1999, investigators funded by NASA's Earth Observations Commercialization and Applications Program collaborated on a nitrogen fertilization experiment with irrigated maize in Nebraska. Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and Landsat 5 Thematic Mapper (TM) data were acquired along with leaf chlorophyll meter and other data on three dates in July during late vegetative growth and early reproductive growth. TGI was consistently correlated with plot-averaged chlorophyll-meter values at the spectral resolutions of AVIRIS, Landsat TM, and digital cameras. Simulations using the Scattering by Arbitrarily Inclined Leaves (SAIL) canopy model indicate an interaction among TGI, leaf area index (LAI) and soil type at low crop LAI, whereas at high LAI and canopy closure, TGI was only affected by leaf chlorophyll content. Therefore, TGI may be the best spectral index to detect crop nitrogen requirements with low-cost digital cameras mounted on low-altitude airborne platforms.     

Spatial analysis of human-induced vegetation productivity decline over eastern Africa using a decade (2001–2011) of medium resolution MODIS time-series data

Climate variation and land transformations related to exploitative land uses are among the main drivers of vegetation productivity decline and ongoing land degradation in East Africa. We combined analysis of vegetation trends and cumulative rain use efficiency differences (CRD), calculated from 250-m MODIS NDVI time-series data, to map vegetation productivity loss over eastern Africa between 2001 and 2011. The CRD index values were furthermore used to discern areas of particular severe vegetation productivity loss over the observation period. Monthly 25-km Tropical Rainfall Measuring Mission (TRMM) data metrics were used to mask areas of rainfall declines not related to human-induced land productivity loss. To provide insights on the productivity decline, we linked the MODIS-based vegetation productivity map to land transformation processes using very high resolution (VHR) imagery in Google Earth (GE) and a Landsat-based land-cover change map. In total, 3.8 million ha experienced significant vegetation loss over the monitoring period. An overall agreement of 68% was found between the rainfall-corrected MODIS productivity decline map and all reference pixels discernable from GE and the Landsat map. The CRD index showed a good potential to discern areas with ‘severe’ vegetation productivity losses under high land-use intensities.     

Impact of sub-pixel heterogeneity on modelled brightness temperature for an agricultural region

Knowledge of sub-pixel heterogeneity, particularly at the passive microwave scale, can improve the brightness temperature (and ultimately the soil moisture) estimation. However, the impact of surface heterogeneity (in terms of soil moisture, soil temperature and vegetation water content) on brightness temperature in an agricultural setting is relatively unknown. The Soil Moisture Active Passive Validation Experiment 2012 (SMAPVEX12) provided an opportunity to evaluate sub-pixel heterogeneity at the scale of a Soil Moisture Ocean Salinity (SMOS) or the Soil Moisture Active Passive (SMAP) radiometer footprint using field measured data. The first objective of this study was to determine if accounting for surface heterogeneity reduced the error between estimated brightness temperature (T_b) and T_b measured by SMOS. It was found that when accounting for variation in surface soil moisture, temperature and vegetation water content within the pixel footprint, the error between the modelled T_b and the measured T_b was less than if a homogeneous pixel were modelled. The correlation between the surface parameters and the error associated with not accounting for surface heterogeneity were investigated. It was found that there was low to moderate correlation between the error and the coefficient of variance associated with the measured soil moisture, soil temperature and vegetation volumetric water content during the field campaign. However, it was found that the correlations changed depending on the stage of vegetation growth and the amount of time following a precipitation event. At the start of the field campaign (following a precipitation event), there was strong correlation between the error and all three surface parameters (r ≥ 0.75). Following a precipitation event close to the middle of the field campaign (during which there was rapid growth in vegetation), there was strong correlation between the error and the variability in vegetation water content (r = 0.89), moderate correlation with soil moisture (r = 0.61) and low correlation with soil temperature (r = 0.26).     

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.