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.     

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.     

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.     

Estimation of soil moisture using optical/thermal infrared remote sensing in the Canadian Prairies

A new approach to estimate soil moisture (SM) based on evaporative fraction (EF) retrieved from optical/thermal infrared MODIS data is presented for Canadian Prairies in parts of Saskatchewan and Alberta. An EF model using the remotely sensed land surface temperature (Ts)/vegetation index concept was modified by incorporating North American Regional Reanalysis (NAAR) Ta data and used for SM estimation. Two different combinations of temperature and vegetation fraction using the difference between Ts from MODIS Aqua and Terra images and Ta from NARR data (Ts−Ta Aqua-day and Ts−Ta Terra-day, respectively) were proposed and the results were compared with those obtained from a previously improved model (ΔTs Aqua-DayNight) as a reference. For the estimation of SM from EF, two empirical models were tested and discussed to find the most appropriate model for converting MODIS-derived EF data to SM values. Estimated SM values were then correlated with in situ SM measurements and their relationships were statistically analyzed. Results indicated statistically significant correlations between SM estimated from all three EF estimation approaches and field measured SM values (R² = 0.42–0.77, p values < 0.04) exhibiting the possibility to estimate SM from remotely sensed EF models. The proposed Ts−Ta MODIS Aqua-day and Terra-day approaches resulted in better estimations of SM (on average higher R² values and similar RMSEs) as compared with the ΔTs reference approach indicating that the concept of incorporating NARR Ta data into Ts/Vegetation index model improved soil moisture estimation accuracy based on evaporative fraction. The accuracies of the predictions were found to be considerably better for intermediate SM values (from 12 to 22 vol/vol%) with square errors averaging below 11 (vol/vol%)². This indicates that the model needs further improvements to account for extreme soil moisture conditions. The findings of this research can be potentially used to downscale SM estimations obtained from passive microwave remote sensing techniques.     

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.     

Forecast of wheat yield throughout the agricultural season using optical and radar satellite images

The aim of this study is to estimate the capabilities of forecasting the yield of wheat using an artificial neural network combined with multi-temporal satellite data acquired at high spatial resolution throughout the agricultural season in the optical and/or microwave domains. Reflectance (acquired by Formosat-2, and Spot 4–5 in the green, red, and near infrared wavelength) and multi-configuration backscattering coefficients (acquired by TerraSAR-X and Radarsat-2 in the X- and C-bands, at co- (abbreviated HH and VV) and cross-polarization states (abbreviated HV and VH)) constitute the input variable of the artificial neural networks, which are trained and validated on the successively acquired images, providing yield forecast in near real-time conditions. The study is based on data collected over 32 fields of wheat distributed over a study area located in southwestern France, near Toulouse. Among the tested sensor configurations, several satellite data appear useful for the yield forecasting throughout the agricultural season (showing coefficient of determination (R²) larger than 0.60 and a root mean square error (RMSE) lower than 9.1 quintals by hectare (q ha⁻¹)): C_VH, C_HV, or the combined used of X_HH and C_HH, C_HH and C_HV, or green reflectance and C_HH. Nevertheless, the best accurate forecast (R² = 0.76 and RMSE = 7.0 q ha⁻¹) is obtained longtime before the harvest (on day 98, during the elongation of stems) using the combination of co- and cross-polarized backscattering coefficients acquired in the C-band (C_VV and C_VH). These results highlight the high interest of using synthetic aperture radar (SAR) data instead of optical ones to early forecast the yield before the harvest of wheat.     

Modelling high arctic percent vegetation cover using field digital images and high resolution satellite data

In this study, digital images collected at a study site in the Canadian High Arctic were processed and classified to examine the spatial-temporal patterns of percent vegetation cover (PVC). To obtain the PVC of different plant functional groups (i.e., forbs, graminoids/sedges and mosses), field near infrared-green-blue (NGB) digital images were classified using an object-based image analysis (OBIA) approach. The PVC analyses comparing different vegetation types confirmed: (i) the polar semi-desert exhibited the lowest PVC with a large proportion of bare soil/rock cover; (ii) the mesic tundra cover consisted of approximately 60% mosses; and (iii) the wet sedge consisted almost exclusively of graminoids and sedges. As expected, the PVC and green normalized difference vegetation index (GNDVI; (R_NIR − R_Green)/(R_NIR + R_Green)), derived from field NGB digital images, increased during the summer growing season for each vegetation type: i.e., ∼5% (0.01) for polar semi-desert; ∼10% (0.04) for mesic tundra; and ∼12% (0.03) for wet sedge respectively. PVC derived from field images was found to be strongly correlated with WorldView-2 derived normalized difference spectral indices (NDSI; (R_x − R_y)/(R_x + R_y)), where R_x is the reflectance of the red edge (724.1 nm) or near infrared (832.9 nm and 949.3 nm) bands; R_y is the reflectance of the yellow (607.7 nm) or red (658.8 nm) bands with R²’s ranging from 0.74 to 0.81. NDSIs that incorporated the yellow band (607.7 nm) performed slightly better than the NDSIs without, indicating that this band may be more useful for investigating Arctic vegetation that often includes large proportions of senescent vegetation throughout the growing season.     

Developing water quality retrieval models with in situ hyperspectral data in Poyang Lake,China

Applying remote sensing techniques to develop the retrieval models and further to obtain the spatiotemporal information of water quality parameters is necessary for understanding, managing, and protecting lake ecosystems. This study aimed to calibrate and validate the retrieval models for estimating the concentrations of chlorophyll a (C_CHL), suspended particulate matter (C_SPM), and dissolved organic carbon (C_DOC) with the in situ hyperspectral measurements in Poyang Lake, China in 2010 and 2011. The model calibration and validation results indicated that: (1) for C_CHL retrieval, significantly strong and moderate correlations existed between the measured and estimated values (with the correlation coefficient r = 0.92 and r = 0.76) using the exponential model and the three-band model, respectively, with biased estimation observed for the exponential model; (2) for retrieving C_SPM, there was a strong correlation between the measured and estimated values (r = 0.95) using the exponential model; and (3) no significant correlation between measured and estimated C_DOC values was found with our developed models. More work is needed to allow the water quality of Poyang Lake to be accurately and steadily estimated, especially for C_CHL and C_DOC.     

Combining UAV-based plant height from crop surface models,visible, and near infrared vegetation indices for biomass monitoring in barley

In this study we combined selected vegetation indices (VIs) and plant height information to estimate biomass in a summer barley experiment. The VIs were calculated from ground-based hyperspectral data and unmanned aerial vehicle (UAV)-based red green blue (RGB) imaging. In addition, the plant height information was obtained from UAV-based multi-temporal crop surface models (CSMs). The test site is a summer barley experiment comprising 18 cultivars and two nitrogen treatments located in Western Germany. We calculated five VIs from hyperspectral data. The normalised ratio index (NRI)-based index GnyLi (Gnyp et al., 2014) showed the highest correlation (R² = 0.83) with dry biomass. In addition, we calculated three visible band VIs: the green red vegetation index (GRVI), the modified GRVI (MGRVI) and the red green blue VI (RGBVI), where the MGRVI and the RGBVI are newly developed VI. We found that the visible band VIs have potential for biomass prediction prior to heading stage. A robust estimate for biomass was obtained from the plant height models (R² = 0.80–0.82). In a cross validation test, we compared plant height, selected VIs and their combination with plant height information. Combining VIs and plant height information by using multiple linear regression or multiple non-linear regression models performed better than the VIs alone. The visible band GRVI and the newly developed RGBVI are promising but need further investigation. However, the relationship between plant height and biomass produced the most robust results. In summary, the results indicate that plant height is competitive with VIs for biomass estimation in summer barley. Moreover, visible band VIs might be a useful addition to biomass estimation. The main limitation is that the visible band VIs work for early growing stages only.     

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.     

Deriving leaf mass per area (LMA) from foliar reflectance across a variety of plant species using continuous wavelet analysis

Leaf mass per area (LMA), the ratio of leaf dry mass to leaf area, is a trait of central importance to the understanding of plant light capture and carbon gain. It can be estimated from leaf reflectance spectroscopy in the infrared region, by making use of information about the absorption features of dry matter. This study reports on the application of continuous wavelet analysis (CWA) to the estimation of LMA across a wide range of plant species. We compiled a large database of leaf reflectance spectra acquired within the framework of three independent measurement campaigns (ANGERS, LOPEX and PANAMA) and generated a simulated database using the PROSPECT leaf optical properties model. CWA was applied to the measured and simulated databases to extract wavelet features that correlate with LMA. These features were assessed in terms of predictive capability and robustness while transferring predictive models from the simulated database to the measured database. The assessment was also conducted with two existing spectral indices, namely the Normalized Dry Matter Index (NDMI) and the Normalized Difference index for LMA (ND_LMA).Five common wavelet features were determined from the two databases, which showed significant correlations with LMA (R²: 0.51–0.82, p < 0.0001). The best robustness (R² = 0.74, RMSE = 18.97 g/m² and Bias = 0.12 g/m²) was obtained using a combination of two low-scale features (1639 nm, scale 4) and (2133 nm, scale 5), the first being predominantly important. The transferability of the wavelet-based predictive model to the whole measured database was either better than or comparable to those based on spectral indices. Additionally, only the wavelet-based model showed consistent predictive capabilities among the three measured data sets. In comparison, the models based on spectral indices were sensitive to site-specific data sets. Integrating the ND_LMA spectral index and the two robust wavelet features improved the LMA prediction. One of the bands used by this spectral index, 1368 nm, was located in a strong atmospheric water absorption region and replacing it with the next available band (1340 nm) led to lower predictive accuracies. However, the two wavelet features were not affected by data quality in the atmospheric absorption regions and therefore showed potential for canopy-level investigations. The wavelet approach provides a different perspective into spectral responses to LMA variation than the traditional spectral indices and holds greater promise for implementation with airborne or spaceborne imaging spectroscopy data for mapping canopy foliar dry biomass.     

A sun–crown–sensor model and adapted C-correction logic for topographic correction of high resolution forest imagery

Canopy shadowing mediated by topography is an important source of radiometric distortion on remote sensing images of rugged terrain. Topographic correction based on the sun–canopy–sensor (SCS) model significantly improved over those based on the sun–terrain–sensor (STS) model for surfaces with high forest canopy cover, because the SCS model considers and preserves the geotropic nature of trees. The SCS model accounts for sub-pixel canopy shadowing effects and normalizes the sunlit canopy area within a pixel. However, it does not account for mutual shadowing between neighboring pixels. Pixel-to-pixel shadowing is especially apparent for fine resolution satellite images in which individual tree crowns are resolved. This paper proposes a new topographic correction model: the sun–crown–sensor (SCnS) model based on high-resolution satellite imagery (IKONOS) and high-precision LiDAR digital elevation model. An improvement on the C-correction logic with a radiance partitioning method to address the effects of diffuse irradiance is also introduced (SCnS + C). In addition, we incorporate a weighting variable, based on pixel shadow fraction, on the direct and diffuse radiance portions to enhance the retrieval of at-sensor radiance and reflectance of highly shadowed tree pixels and form another variety of SCnS model (SCnS + W). Model evaluation with IKONOS test data showed that the new SCnS model outperformed the STS and SCS models in quantifying the correlation between terrain-regulated illumination factor and at-sensor radiance. Our adapted C-correction logic based on the sun–crown–sensor geometry and radiance partitioning better represented the general additive effects of diffuse radiation than C parameters derived from the STS or SCS models. The weighting factor W_t also significantly enhanced correction results by reducing within-class standard deviation and balancing the mean pixel radiance between sunlit and shaded slopes. We analyzed these improvements with model comparison on the red and near infrared bands. The advantages of SCnS + C and SCnS + W on both bands are expected to facilitate forest classification and change detection applications.     

The brightness temperature adjusted dust index: An improved approach to detect dust storms using MODIS imagery

Moderate Resolution Imaging Spectroradiometer (MODIS) imagery provides a good data source for timely and accurate monitoring of dust storms. However, effective MODIS-based dust indices are inadequate. In this study, we proposed an improved brightness temperature adjusted dust index (BADI) by integrating the brightness temperatures of three thermal infrared MODIS bands: band₂₀ (3.66–3.84 μm), band₃₁ (10.78–11.28 μm) and band₃₂ (11.77–12.27 μm). We used the BADI to monitor several representative dust storms over the Northeast Asia between 2000 and 2011. When compared to commonly used MODIS-based dust indices, such as the brightness temperature difference index in band₃₂ and band₃₁ (BTD_32-31) and the normalized difference dust index (NDDI), the BADI captured the spatial extent and density of dust storms more accurately. The BADI detected dust storm extent with an overall accuracy >90%, which was 7% and 29% higher than the results derived from BTD_32-31 and NDDI, respectively. The BADI also demonstrated good agreement with the density indicator of MODIS Deep Blue Aerosol Optical Depth (R² = 0.59, P < 0.01). We suggest that the BADI is an effective tool to monitor large-scale dust storms.     

Assessing the potential of MODIS/Terra version 5 images to improve near shore lake bathymetric surveys

Estimation of storage volume of inland water bodies is made using bathymetric surveys. Such surveys are often terminated far inward from shores due to operational limitations. The gap in depth data for the shore area is often overcome by extrapolating the established volume-area relationship beyond the survey boundary. This introduces large errors in volume estimation. In this study, the potential of the Moderate Resolution Imaging Spectroradiometer MODIS version-5 images to improve an existing bathymetric model is assessed over Lake Tana, Ethiopia. MODIS images are used to generate a new bathymetric model of the unsurveyed shore. The new model is validated through its ability to reproduce measured lake levels. The result is compared with lake level estimates made from an existing bathymetric survey. Results showed that MODIS improved bathymetric model reproduced lake water level with a RMSE of 0.14 m as compared to 0.82 m using the existing bathymetric model. The Mean Absolute Error (MAE) has also improved from 0.78 m to 0.13 m. These results suggest that despite their unavailability on cloudy days, MODIS images provide significant advantages for improving bathymetric models.     

Empirical models for estimating the suspended sediment concentration in Amazonian white water rivers using Landsat 5/TM

Suspended sediment yield is a very important environmental indicator within Amazonian fluvial systems, especially for rivers dominated by inorganic particles, referred to as white water rivers. For vast portions of Amazonian rivers, suspended sediment concentration (SSC) is measured infrequently or not at all. However, remote sensing techniques have been used to estimate water quality parameters worldwide, from which data for suspended matter is the most successfully retrieved. This paper presents empirical models for SSC retrieval in Amazonian white water rivers using reflectance data derived from Landsat 5/TM. The models use multiple regression for both the entire dataset (global model, N = 504) and for five segmented datasets (regional models) defined by general geological features of drainage basins. The models use VNIR bands, band ratios, and the SWIR band 5 as input. For the global model, the adjusted R² is 0.76, while the adjusted R² values for regional models vary from 0.77 to 0.89, all significant (p-value < 0.0001). The regional models are subject to the leave-one-out cross validation technique, which presents robust results. The findings show that both the average error of estimation and the standard deviation increase as the SSC range increases. Regional models were more accurate when compared with the global model, suggesting changes in optical proprieties of water sampled at different sampling stations. Results confirm the potential for the estimation of SSC from Landsat/TM historical series data for the 1980s and 1990s, for which the in situ database is scarce. Such estimates supplement the SSC temporal series, providing a more comprehensive SSC temporal series which may show environmental dynamics yet unknown.     

Retrieval of leaf area index in alpine wetlands using a two-layer canopy reflectance model

In this paper, we focused on the retrieval of the LAI in an alpine wetland located in western part of China in late August and early July 2011. A two-layer canopy reflectance model (ACRM) was used to establish the relationships between the LAI and the reflectance of near-infrared (NIR) and red (RED) wavebands. The reflectance data were derived from Landsat TM L1T product and the Terra and Aqua MODIS 16-day and 8-day composite reflectance products (MOD/MYD09) at 250 m resolution. Due to the lack of the information about some major input parameters for ACRM, which are sensitive to model outputs in the reflectance of NIR and RED wavebands, the inverse problem was ill-posed. To overcome this problem, a method of increasing the sensitivity of the LAI while reducing the influence of other model free parameters based on the study of free parameters’ sensitivity to the ACRM outputs and the region’s features was studied. The area of interest was divided into two parts using the approximately statistic normalized difference vegetation index (NDVI) value around 0.5. One part was sparse vegetation (0.1 < NDVI < 0.5), which is more sensitive to soil background effects and less sensitive to the canopy biophysical and biochemical variables. The other part was dense vegetation (0.5 ≤ NDVI < 1.0), which is less sensitive to soil background effects and more sensitive to plant canopies and leaf parameters. Then, the relationships of ρ_nir–LAI and ρ_red–LAI were established using a look-up table algorithm for the two parts. Furthermore, a regularization technique for fast pixel-wise retrieval was introduced to reduce the elements of LUT sets while maintaining a relatively high accuracy. The results were very promising compared to the field measured LAI values that the correlation (R²) of the measured LAI values and retrieved LAI values reached 0.95, and the root-mean-square deviation (RMSD) was 0.33 for late August, 2011, while the R² reached 0.82 and RMSD was 0.25 for early July 2011.     

Brown and green LAI mapping through spectral indices

When crops senescence, leaves remain until they fall off or are harvested. Hence, leaf area index (LAI) stays high even when chlorophyll content degrades to zero. Current LAI approaches from remote sensing techniques are not optimized for estimating LAI of senescent vegetation. In this paper a two-step approach has been proposed to realize simultaneous LAI mapping over green and senescent croplands. The first step separates green from brown LAI by means of a newly proposed index, ‘Green Brown Vegetation Index (GBVI)’. This index exploits two shortwave infrared (SWIR) spectral bands centred at 2100 and 2000 nm, which fall right in the dry matter absorption regions, thereby providing positive values for senescent vegetation and negative for green vegetation. The second step involves applying linear regression functions based on optimized vegetation indices to estimate green and brown LAI estimation respectively. While the green LAI index uses a band in the red and a band in the red-edge, the brown LAI index uses bands located in the same spectral region as GBVI, i.e. an absorption band located in the region of maximum absorption of cellulose and lignin at 2154 nm, and a reference band at 1635 nm where the absorption of both water and dry matter is low. The two-step approach was applied to a HyMap image acquired over an agroecosystem at the agricultural site Barrax, Spain.     

Hyperspectral-based predictive modelling of grapevine water status in the Portuguese Douro wine region

In this study, hyperspectral reflectance (HySR) data derived from a handheld spectroradiometer were used to assess the water status of three grapevine cultivars in two sub-regions of Douro wine region during two consecutive years. A large set of potential predictors derived from the HySR data were considered for modelling/predicting the predawn leaf water potential (Ψ_pd) through different statistical and machine learning techniques. Three HySR vegetation indices were selected as final predictors for the computation of the models and the in-season time trend was removed from data by using a time predictor. The vegetation indices selected were the Normalized Reflectance Index for the wavelengths 554 nm and 561 nm (NRI_554;561), the water index (WI) for the wavelengths 900 nm and 970 nm, and the D1 index which is associated with the rate of reflectance increase in the wavelengths of 706 nm and 730 nm. These vegetation indices covered the green, red edge and the near infrared domains of the electromagnetic spectrum. A large set of state-of-the-art analysis and statistical and machine-learning modelling techniques were tested. Predictive modelling techniques based on generalized boosted model (GBM), bagged multivariate adaptive regression splines (B-MARS), generalized additive model (GAM), and Bayesian regularized neural networks (BRNN) showed the best performance for predicting Ψ_pd, with an average determination coefficient (R²) ranging between 0.78 and 0.80 and RMSE varying between 0.11 and 0.12 MPa. When cultivar Touriga Nacional was used for training the models and the cultivars Touriga Franca and Tinta Barroca for testing (independent validation), the models performance was good, particularly for GBM (R² = 0.85; RMSE = 0.09 MPa). Additionally, the comparison of Ψ_pd observed and predicted showed an equitable dispersion of data from the various cultivars. The results achieved show a good potential of these predictive models based on vegetation indices to support irrigation scheduling in vineyard.     

Spectral resampling based on user-defined inter-band correlation filter: C3 and C4 grass species classification

In this paper, a user-defined inter-band correlation filter function was used to resample hyperspectral data and thereby mitigate the problem of multicollinearity in classification analysis. The proposed resampling technique convolves the spectral dependence information between a chosen band-centre and its shorter and longer wavelength neighbours. Weighting threshold of inter-band correlation (WTC, Pearson's r) was calculated, whereby r = 1 at the band-centre. Various WTC (r = 0.99, r = 0.95 and r = 0.90) were assessed, and bands with coefficients beyond a chosen threshold were assigned r = 0. The resultant data were used in the random forest analysis to classify in situ C₃ and C₄ grass canopy reflectance. The respective WTC datasets yielded improved classification accuracies (kappa = 0.82, 0.79 and 0.76) with less correlated wavebands when compared to resampled Hyperion bands (kappa = 0.76). Overall, the results obtained from this study suggested that resampling of hyperspectral data should account for the spectral dependence information to improve overall classification accuracy as well as reducing the problem of multicollinearity.     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.