Determining crop acreage estimates for specific winter crops using shape attributes from sequential MODIS imagery

There are increasing societal and plant industry demands for more accurate, objective and near real-time crop production information to meet both economic and food security concerns. The advent of the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite platform has augmented the capability of satellite-based applications to monitor large agricultural areas at acceptable pixel scale, cost and accuracy. Fitting parametric profiles to growing season vegetation index time series reduces the volume of data and provides simple quantitative parameters that relates to crop phenology (sowing date, flowering). In this study, we modelled various Gaussian profiles to time sequential MODIS enhanced vegetation index (EVI) images over winter crops in Queensland, Australia. Three simple Gaussian models were evaluated in their effectiveness to identify and classify various winter crop types and coverage at both pixel and regional scales across Queensland's main agricultural areas. Equal to or greater than 93% classification accuracies were obtained in determining crop acreage estimates at pixel scale for each of the Gaussian modelled approaches. Significant high to moderate correlations (log-linear transformation) were also obtained for determining total winter crop (R² = 0.93) areas as well as specific crop acreage for wheat (R² = 0.86) and barley (R² = 0.83). Conversely, it was much more difficult to predict chickpea acreage (R² ≤ 0.26), mainly due to very large uncertainties in survey data. The quantitative approach utilised here further had additional benefits of characterising crop phenology in terms of length of growing season and providing regression diagnostics of how well the fitted profiles matched the EVI time series. The Gaussian curve models utilised here are novel in application and therefore will enhance the use and adoption of remote sensing technologies in targeted agricultural application. With innate simplicity and accuracies comparable to other more convoluted multi-temporal approaches it is a good candidate in determining total and specific crop acreage estimates in future national and global food security frameworks.     

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.     

Evaluating a satellite-based seasonal evapotranspiration product and identifying its relationship with other satellite-derived products and crop yield: A case study for Ethiopia

Satellite-derived evapotranspiration anomalies and normalized difference vegetation index (NDVI) products from Moderate Resolution Imaging Spectroradiometer (MODIS) data are currently used for African agricultural drought monitoring and food security status assessment. In this study, a process to evaluate satellite-derived evapotranspiration (ETa) products with a geospatial statistical exploratory technique that uses NDVI, satellite-derived rainfall estimate (RFE), and crop yield data has been developed. The main goal of this study was to evaluate the ETa using the NDVI and RFE, and identify a relationship between the ETa and Ethiopia’s cereal crop (i.e., teff, sorghum, corn/maize, barley, and wheat) yields during the main rainy season. Since crop production is one of the main factors affecting food security, the evaluation of remote sensing-based seasonal ETa was done to identify the appropriateness of this tool as a proxy for monitoring vegetation condition in drought vulnerable and food insecure areas to support decision makers. The results of this study showed that the comparison between seasonal ETa and RFE produced strong correlation (R² > 0.99) for all 41 crop growing zones in Ethiopia. The results of the spatial regression analyses of seasonal ETa and NDVI using Ordinary Least Squares and Geographically Weighted Regression showed relatively weak yearly spatial relationships (R² < 0.7) for all cropping zones. However, for each individual crop zones, the correlation between NDVI and ETa ranged between 0.3 and 0.84 for about 44% of the cropping zones. Similarly, for each individual crop zones, the correlation (R²) between the seasonal ETa anomaly and de-trended cereal crop yield was between 0.4 and 0.82 for 76% (31 out of 41) of the crop growing zones. The preliminary results indicated that the ETa products have a good predictive potential for these 31 identified zones in Ethiopia. Decision makers may potentially use ETa products for monitoring cereal crop yields and early warning of food insecurity during drought years for these identified zones.     

Improved maize cultivated area estimation over a large scale combining MODIS–EVI time series data and crop phenological information

The accurate and timely information of crop area is vital for crop production and food security. In this study, the Enhanced Vegetation Index (EVI) data from MODerate resolution Imaging Spectroradiometer (MODIS) integrated crop phenological information was used to estimate the maize cultivated area over a large scale in Northeast China. The fine spatial resolution China’s Environment Satellite (HJ-1 satellite) images and the support vector machine (SVM) algorithm were employed to discriminate distribution of maize in the reference area. The mean MODIS–EVI time series curve of maize was extracted in the reference area by using multiple periods MODIS–EVI data. By analysing the temporal shift of crop calendars from northern to southern parts in Northeast China, the lag value was derived from phenological data of twenty-one agro-meteorological stations; here integrating with the mean MODIS–EVI time series image of maize, a standard MODIS–EVI time series image of maize was obtained in the whole study area. By calculating mean absolute distances (MAD) map between standard MODIS–EVI image and mean MODIS–EVI time series images, and setting appropriate thresholds in three provinces, the maize cultivated area was extracted in Northeast China. The results showed that the overall classification accuracy of maize cultivated area was approximately 79%. At the county level, the MODIS-derived maize cultivated area and statistical data were well correlated (R² = 0.82, RMSE = 283.98) over whole Northeast China. It demonstrated that MODIS–EVI time series data integrated with crop phenological information can be used to improve the extraction accuracy of crop cultivated area over a large scale.     

The assimilation of spectral sensing and the WOFOST model for the dynamic simulation of cadmium accumulation in rice tissues

The accurate detection of heavy metal-induced stress on crop growth is important for food security and agricultural, ecological and environmental protection. Spectral sensing offers an efficient and undamaged observation tool to monitor soil and vegetation contamination. This study proposed a methodology for dynamically estimating the total cadmium (Cd) accumulation in rice tissues by assimilating spectral information into WOFOST (World Food Study) model. Based on the differences among ground hyperspectral data of rice in three experiments fields under different Cd concentration levels, the spectral indices MCARI1, NREP and RH were selected to reflect the rice stress condition and dry matter production of rice. With assimilating these sensitive spectral indices into the WOFOST + PROSPECT + SAIL model to optimize the Cd pollution stress factor f_wi, the dynamic dry matter production processes of rice were adjusted. Based on the relation between dry matter production and Cd accumulation, we dynamically simulating the Cd accumulation in rice tissues. The results showed that the method performed well in dynamically estimating the total amount of Cd accumulation in rice tissues with R² over 85%. This study suggests that the proposed method of integrating the spectral information and the crop growth model could successfully dynamically simulate the Cd accumulation in rice tissues.     

Estimating crop net primary production using national inventory data and MODIS-derived parameters

National estimates of spatially-resolved cropland net primary production (NPP) are needed for diagnostic and prognostic modeling of carbon sources, sinks, and net carbon flux between land and atmosphere. Cropland NPP estimates that correspond with existing cropland cover maps are needed to drive biogeochemical models at the local scale as well as national and continental scales. Existing satellite-based NPP products tend to underestimate NPP on croplands. An Agricultural Inventory-based Light Use Efficiency (AgI-LUE) framework was developed to estimate individual crop biophysical parameters for use in estimating crop-specific NPP over large multi-state regions. The method is documented here and evaluated for corn (Zea mays L.) and soybean (Glycine max L. Merr.) in Iowa and Illinois in 2006 and 2007. The method includes a crop-specific Enhanced Vegetation Index (EVI), shortwave radiation data estimated using the Mountain Climate Simulator (MTCLIM) algorithm, and crop-specific LUE per county. The combined aforementioned variables were used to generate spatially-resolved, crop-specific NPP that corresponds to the Cropland Data Layer (CDL) land cover product. Results from the modeling framework captured the spatial NPP gradient across croplands of Iowa and Illinois, and also represented the difference in NPP between years 2006 and 2007. Average corn and soybean NPP from AgI-LUE was 917 g C m⁻² yr⁻¹ and 409 g C m⁻² yr⁻¹, respectively. This was 2.4 and 1.1 times higher, respectively, for corn and soybean compared to the MOD17A3 NPP product. Site comparisons with flux tower data show AgI-LUE NPP in close agreement with tower-derived NPP, lower than inventory-based NPP, and higher than MOD17A3 NPP. The combination of new inputs and improved datasets enabled the development of spatially explicit and reliable NPP estimates for individual crops over large regional extents.     

Remotely sensed biomass over steep slopes: An evaluation among successional stands of the Atlantic Forest,Brazil

Remotely sensed images have been widely used to model biomass and carbon content on large spatial scales. Nevertheless, modeling biomass using remotely sensed data from steep slopes is still poorly understood. We investigated how topographical features affect biomass estimation using remotely sensed data and how such estimates can be used in the characterization of successional stands in the Atlantic Rainforest in southeastern Brazil. We estimated forest biomass using a modeling approach that included the use of both satellite data (LANDSAT) and topographic features derived from a digital elevation model (TOPODATA). Biomass estimations exhibited low error predictions (Adj. R² = 0.67 and RMSE = 35 Mg/ha) when combining satellite data with a secondary geomorphometric variable, the illumination factor, which is based on hill shading patterns. This improved biomass prediction helped us to determine carbon stock in different forest successional stands. Our results provide an important source of modeling information about large-scale biomass in remaining forests over steep slopes.     

Predicting stem borer density in maize using RapidEye data and generalized linear models

Average maize yield in eastern Africa is 2.03 t ha⁻¹ as compared to global average of 6.06 t ha⁻¹ due to biotic and abiotic constraints. Amongst the biotic production constraints in Africa, stem borers are the most injurious. In eastern Africa, maize yield losses due to stem borers are currently estimated between 12% and 21% of the total production. The objective of the present study was to explore the possibility of RapidEye spectral data to assess stem borer larva densities in maize fields in two study sites in Kenya. RapidEye images were acquired for the Bomet (western Kenya) test site on the 9th of December 2014 and on 27th of January 2015, and for Machakos (eastern Kenya) a RapidEye image was acquired on the 3rd of January 2015. Five RapidEye spectral bands as well as 30 spectral vegetation indices (SVIs) were utilized to predict per field maize stem borer larva densities using generalized linear models (GLMs), assuming Poisson (‘Po’) and negative binomial (‘NB’) distributions. Root mean square error (RMSE) and ratio prediction to deviation (RPD) statistics were used to assess the models performance using a leave-one-out cross-validation approach. The Zero-inflated NB (‘ZINB’) models outperformed the ‘NB’ models and stem borer larva densities could only be predicted during the mid growing season in December and early January in both study sites, respectively (RMSE = 0.69–1.06 and RPD = 8.25–19.57). Overall, all models performed similar when all the 30 SVIs (non-nested) and only the significant (nested) SVIs were used. The models developed could improve decision making regarding controlling maize stem borers within integrated pest management (IPM) interventions.     

Retrieval of wheat leaf area index from AWiFS multispectral data using canopy radiative transfer simulation

Accurate representation of leaf area index (LAI) from high resolution satellite observations is obligatory for various modelling exercises and predicting the precise farm productivity. Present study compared the two retrieval approach based on canopy radiative transfer (CRT) method and empirical method using four vegetation indices (VI) (e.g. NDVI, NDWI, RVI and GNDVI) to estimate the wheat LAI. Reflectance observations available at very high (56 m) spatial resolution from Advanced Wide-Field Sensor (AWiFS) sensor onboard Indian Remote Sensing (IRS) P6, Resourcesat-1 satellite was used in this study. This study was performed over two different wheat growing regions, situated in different agro-climatic settings/environments: Trans-Gangetic Plain Region (TGPR) and Central Plateau and Hill Region (CPHR). Forward simulation of canopy reflectances in four AWiFS bands viz. green (0.52–0.59 μm), red (0.62–0.68 μm), NIR (0.77–0.86 μm) and SWIR (1.55–1.70 μm) were carried out to generate the look up table (LUT) using CRT model PROSAIL from all combinations of canopy intrinsic variables. An inversion technique based on minimization of cost function was used to retrieve LAI from LUT and observed AWiFS surface reflectances. Two consecutive wheat growing seasons (November 2005–March 2006 and November 2006–March 2007) datasets were used in this study. The empirical models were developed from first season data and second growing season data used for validation. Among all the models, LAI-NDVI empirical model showed the least RMSE (root mean square error) of 0.54 and 0.51 in both agro-climatic regions respectively. The comparison of PROSAIL retrieved LAI with in situ measurements of 2006–2007 over the two agro-climatic regions produced substantially less RMSE of 0.34 and 0.41 having more R² of 0.91 and 0.95 for TGPR and CPHR respectively in comparison to empirical models. Moreover, CRT retrieved LAI had less value of errors in all the LAI classes contrary to empirical estimates. The PROSAIL based retrieval has potential for operational implementation to determine the regional crop LAI and can be extendible to other regions after rigorous validation exercise.     

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.     

Satellite-based investigation of flood-affected rice cultivation areas in Chao Phraya River Delta,Thailand

The occurrence of catastrophic floods in Thailand in 2011 caused significant damage to rice agriculture. This study investigated flood-affected rice cultivation areas in the Chao Phraya River Delta (CRD) rice bowl, Thailand using time-series moderate resolution imaging spectroradiometer (MODIS) data. The data were processed for 2008 (normal flood year) and 2011, comprising four main steps: (1) data pre-processing to construct time-series MODIS vegetation indices (VIs), to filter noise from the time-series VIs by the empirical mode decomposition (EMD), and to mask out non-agricultural areas in respect to water-related cropping areas; (2) flood-affected area classification using the unsupervised linear mixture model (ULMM); (3) rice crop classification using the support vector machines (SVM); and (4) accuracy assessment of flood and rice crop mapping results. The comparisons between the flood mapping results and the ground reference data indicated an overall accuracy of 97.9% and Kappa coefficient of 0.62 achieved for 2008, and 95.7% and 0.77 for 2011, respectively. These results were reaffirmed by close agreement (R² > 0.8) between comparisons of the two datasets at the provincial level. The crop mapping results compared with the ground reference data revealed that the overall accuracies and Kappa coefficients obtained for 2008 were 88.5% and 0.82, and for 2011 were 84.1% and 0.76, respectively. A strong correlation was also found between MODIS-derived rice area and rice area statistics at the provincial level (R² > 0.7). Rice crop maps overlaid on the flood-affected area maps showed that approximately 16.8% of the rice cultivation area was affected by floods in 2011 compared to 4.9% in 2008. A majority of the flood-expanded area was observed for the double-cropped rice (10.5%), probably due to flood-induced effects to the autumn–summer and rainy season crops. Information achieved from this study could be useful for agricultural planners to mitigate possible impacts of floods on rice production.     

Agriculture pest and disease risk maps considering MSG satellite data and land surface temperature

Pest risk maps for agricultural use are usually constructed from data obtained from in-situ meteorological weather stations, which are relatively sparsely distributed and are often quite expensive to install and difficult to maintain. This leads to the creation of maps with relatively low spatial resolution, which are very much dependent on interpolation methodologies. Considering that agricultural applications typically require a more detailed scale analysis than has traditionally been available, remote sensing technology can offer better monitoring at increasing spatial and temporal resolutions, thereby, improving pest management results and reducing costs. This article uses ground temperature, or land surface temperature (LST), data distributed by EUMETSAT/LSASAF (with a spatial resolution of 3 × 3 km (nadir resolution) and a revisiting time of 15 min) to generate one of the most commonly used parameters in pest modeling and monitoring: “thermal integral over air temperature (accumulated degree-days)”. The results show a clear association between the accumulated LST values over a threshold and the accumulated values computed from meteorological stations over the same threshold (specific to a particular tomato pest). The results are very promising and enable the production of risk maps for agricultural pests with a degree of spatial and temporal detail that is difficult to achieve using in-situ meteorological stations.     

Detecting forest damage after a low-severity fire using remote sensing at multiple scales

Remote sensing technologies are an ideal platform to examine the extent and impact of fire on the landscape. In this study we assess that capacity of the RapidEye constellation and Landsat (Thematic Mapper and Operational Land Imager to map fine-scale burn attributes for a small, low severity prescribed fire in a dry Western Canadian forest. Estimates of burn severity from field data were collated into a simple burn index and correlated with a selected suite of common spectral vegetation indices. Burn severity classes were then derived to map fire impacts and estimate consumed woody surface fuels (diameter ≥2.6 cm). All correlations between the simple burn index and vegetation indices produced significant results (p < 0.01), but varied substantially in their overall accuracy. Although the Landsat Soil Adjusted Vegetation Index provided the best regression fit (R² = 0.56), results suggested that RapidEye provided much more spatially detailed estimates of tree damage (Soil Adjusted Vegetation Index, R² = 0.51). Consumption estimates of woody surface fuels ranged from 3.38 ± 1.03 Mg ha⁻¹ to 11.73 ± 1.84 Mg ha⁻¹, across four derived severity classes with uncertainties likely a result of changing foliage moisture between the before and after fire images. While not containing spectral information in the short wave infrared, the spatial variability provided by the RapidEye imagery has potential for mapping and monitoring fine scale forest attributes, as well as the potential to resolve fire damage at the individual tree level.     

Effects of the partitioning of diffuse and direct solar radiation on satellite-based modeling of crop gross primary production

Modeling crop gross primary production (GPP) is critical to understanding the carbon dynamics of agro-ecosystems. Satellite-based studies have widely used production efficiency models (PEM) to estimate cropland GPP, wherein light use efficiency (LUE) is a key model parameter. One factor that has not been well considered in many PEMs is that canopy LUE could vary with illumination conditions. This study investigates how the partitioning of diffuse and direct solar radiation influences cropland GPP using both flux tower and satellite data. The field-measured hourly LUE under cloudy conditions was 1.50 and 1.70 times higher than that under near clear-sky conditions for irrigated corn and soybean, respectively. We applied a two-leaf model to simulate the canopy radiative transfer process, where modeled photosynthetically active radiation (PAR) absorbed by canopy agreed with tower measurements (R² = 0.959 and 0.914 for corn and soybean, respectively). Derived canopy LUE became similar after accounting for the impact of light saturation on leaf photosynthetic capacity under varied illumination conditions. The impacts of solar radiation partitioning on satellite-based modeling of crop GPP was examined using vegetation indices (VI) derived from MODIS data. Consistent with the field modeling results, the relationship between daily GPP and PAR × VI under varied illumination conditions showed different patterns in terms of regression slope and intercept. We proposed a function to correct the influences of direct and diffuse radiation partitioning and the explained variance of flux tower GPP increased in all experiments. Our results suggest that the non-linear response of leaf photosynthesis to light absorption contributes to higher canopy LUE on cloudy days than on clear days. We conclude that accounting for the impacts of solar radiation partitioning is necessary for modeling crop GPP on a daily or shorter basis.     

Aboveground biomass mapping in French Guiana by combining remote sensing,forest inventories and environmental data

Mapping forest aboveground biomass (AGB) has become an important task, particularly for the reporting of carbon stocks and changes. AGB can be mapped using synthetic aperture radar data (SAR) or passive optical data. However, these data are insensitive to high AGB levels (>150 Mg/ha, and >300 Mg/ha for P-band), which are commonly found in tropical forests. Studies have mapped the rough variations in AGB by combining optical and environmental data at regional and global scales. Nevertheless, these maps cannot represent local variations in AGB in tropical forests. In this paper, we hypothesize that the problem of misrepresenting local variations in AGB and AGB estimation with good precision occurs because of both methodological limits (signal saturation or dilution bias) and a lack of adequate calibration data in this range of AGB values. We test this hypothesis by developing a calibrated regression model to predict variations in high AGB values (mean >300 Mg/ha) in French Guiana by a methodological approach for spatial extrapolation with data from the optical geoscience laser altimeter system (GLAS), forest inventories, radar, optics, and environmental variables for spatial inter- and extrapolation. Given their higher point count, GLAS data allow a wider coverage of AGB values. We find that the metrics from GLAS footprints are correlated with field AGB estimations (R² = 0.54, RMSE = 48.3 Mg/ha) with no bias for high values. First, predictive models, including remote-sensing, environmental variables and spatial correlation functions, allow us to obtain “wall-to-wall” AGB maps over French Guiana with an RMSE for the in situ AGB estimates of ∼50 Mg/ha and R² = 0.66 at a 1-km grid size. We conclude that a calibrated regression model based on GLAS with dependent environmental data can produce good AGB predictions even for high AGB values if the calibration data fit the AGB range. We also demonstrate that small temporal and spatial mismatches between field data and GLAS footprints are not a problem for regional and global calibrated regression models because field data aim to predict large and deep tendencies in AGB variations from environmental gradients and do not aim to represent high but stochastic and temporally limited variations from forest dynamics. Thus, we advocate including a greater variety of data, even if less precise and shifted, to better represent high AGB values in global models and to improve the fitting of these models for high values.     

An empirical standardized soil moisture index for agricultural drought assessment from remotely sensed data

We propose a simple, spatially invariant and probabilistic year-round Empirical Standardized Soil Moisture Index (ESSMI) that is designed to classify soil moisture anomalies from harmonized multi-satellite surface data into categories of agricultural drought intensity. The ESSMI is computed by fitting a nonparametric empirical probability density function (ePDF) to historical time-series of soil moisture observations and then transforming it into a normal distribution with a mean of zero and standard deviation of one. Negative standard normal values indicate dry soil conditions, whereas positive values indicate wet soil conditions. Drought intensity is defined as the number of negative standard deviations between the observed soil moisture value and the respective normal climatological conditions. To evaluate the performance of the ESSMI, we fitted the ePDF to the Essential Climate Variable Soil Moisture (ECV SM) v02.0 data values collected in the period between January 1981 and December 2010 at South–Central America, and compared the root-mean-square-errors (RMSE) of residuals with those of beta and normal probability density functions (bPDF and nPDF, respectively). Goodness-of-fit results attained with time-series of ECV SM values averaged at monthly, seasonal, half-yearly and yearly timescales suggest that the ePDF provides triggers of agricultural drought onset and intensity that are more accurate and precise than the bPDF and nPDF. Furthermore, by accurately mapping the occurrence of major drought events over the last three decades, the ESSMI proved to be spatio-temporal consistent and the ECV SM data to provide a well calibrated and homogenized soil moisture climatology for the region. Maize, soybean and wheat crop yields in the region are highly correlated (r > 0.82) with cumulative ESSMI values computed during the months of critical crop growing, indicating that the nonparametric index of soil moisture anomalies can be used for agricultural drought assessment.     

Improving the prediction of arsenic contents in agricultural soils by combining the reflectance spectroscopy of soils and rice plants

Visible and near-infrared reflectance spectroscopy provides a beneficial tool for investigating soil heavy metal contamination. This study aimed to investigate mechanisms of soil arsenic prediction using laboratory based soil and leaf spectra, compare the prediction of arsenic content using soil spectra with that using rice plant spectra, and determine whether the combination of both could improve the prediction of soil arsenic content. A total of 100 samples were collected and the reflectance spectra of soils and rice plants were measured using a FieldSpec3 portable spectroradiometer (350–2500 nm). After eliminating spectral outliers, the reflectance spectra were divided into calibration (n = 62) and validation (n = 32) data sets using the Kennard-Stone algorithm. Genetic algorithm (GA) was used to select useful spectral variables for soil arsenic prediction. Thereafter, the GA-selected spectral variables of the soil and leaf spectra were individually and jointly employed to calibrate the partial least squares regression (PLSR) models using the calibration data set. The regression models were validated and compared using independent validation data set. Furthermore, the correlation coefficients of soil arsenic against soil organic matter, leaf arsenic and leaf chlorophyll were calculated, and the important wavelengths for PLSR modeling were extracted. Results showed that arsenic prediction using the leaf spectra (coefficient of determination in validation, R_v² = 0.54; root mean square error in validation, RMSE_v = 12.99 mg kg⁻¹; and residual prediction deviation in validation, RPD_v = 1.35) was slightly better than using the soil spectra (R_v² = 0.42, RMSE_v = 13.35 mg kg⁻¹, and RPD_v = 1.31). However, results also showed that the combinational use of soil and leaf spectra resulted in higher arsenic prediction (R_v² = 0.63, RMSE_v = 11.94 mg kg⁻¹, RPD_v = 1.47) compared with either soil or leaf spectra alone. Soil spectral bands near 480, 600, 670, 810, 1980, 2050 and 2290 nm, leaf spectral bands near 700, 890 and 900 nm in PLSR models were important wavelengths for soil arsenic prediction. Moreover, soil arsenic showed significantly positive correlations with soil organic matter (r = 0.62, p < 0.01) and leaf arsenic (r = 0.77, p < 0.01), and a significantly negative correlation with leaf chlorophyll (r = −0.67, p < 0.01). The results showed that the prediction of arsenic contents using soil and leaf spectra may be based on their relationships with soil organic matter and leaf chlorophyll contents, respectively. Although RPD of 1.47 was below the recommended RPD of >2 for soil analysis, arsenic prediction in agricultural soils can be improved by combining the leaf and soil spectra.     

Mapping land cover gradients through analysis of hyper-temporal NDVI imagery

The green cover of the earth exhibits various spatial gradients that represent gradual changes in space of vegetation density and/or in species composition. To date, land cover mapping methods differentiate at best, mapping units with different cover densities and/or species compositions, but typically fail to express such differences as gradients. Present interpretation techniques still make insufficient use of freely available spatial-temporal Earth Observation (EO) data that allow detection of existing land cover gradients. This study explores the use of hyper-temporal NDVI imagery to detect and delineate land cover gradients analyzing the temporal behavior of NDVI values. MODIS-Terra MVC-images (250 m, 16-day) of Crete, Greece, from February 2000 to July 2009 are used. The analysis approach uses an ISODATA unsupervised classification in combination with a Hierarchical Clustering Analysis (HCA). Clustering of class-specific temporal NDVI profiles through HCA resulted in the identification of gradients in landcover vegetation growth patterns. The detected gradients were arranged in a relational diagram, and mapped. Three groups of NDVI-classes were evaluated by correlating their class-specific annual average NDVI values with the field data (tree, shrub, grass, bare soil, stone, litter fraction covers). Multiple regression analysis showed that within each NDVI group, the fraction cover data were linearly related with the NDVI data, while NDVI groups were significantly different with respect to tree cover (adj. R² = 0.96), shrub cover (adj. R² = 0.83), grass cover (adj. R² = 0.71), bare soil (adj. R² = 0.88), stone cover (adj. R² = 0.83) and litter cover (adj. R² = 0.69) fractions. Similarly, the mean Sorenson dissimilarity values were found high and significant at confidence interval of 95% in all pairs of three NDVI groups. The study demonstrates that hyper-temporal NDVI imagery can successfully detect and map land cover gradients. The results may improve land cover assessment and aid in agricultural and ecological studies.     

A new 500-m resolution map of canopy height for Amazon forest using spaceborne LiDAR and cloud-free MODIS imagery

In the present study, we aimed to map canopy heights in the Brazilian Amazon mainly on the basis of spaceborne LiDAR and cloud-free MODIS imagery with a new method (the Self-Organizing Relationships method) for spatial modeling of the LiDAR footprint. To evaluate the general versatility, we compared the created canopy height map with two different canopy height estimates on the basis of our original field study plots (799 plots located in eight study sites) and a previously developed canopy height map. The compared canopy height estimates were obtained by: (1) a stem diameter at breast height (D) – tree height (H) relationship specific to each site on the basis of our original field study, (2) a previously developed D–H model involving environmental and structural factors as explanatory variables (Feldpausch et al., 2011), and (3) a previously developed canopy height map derived from the spaceborne LiDAR data with different spatial modeling method and explanatory variables (Simard et al., 2011). As a result, our canopy height map successfully detected a spatial distribution pattern in canopy height estimates based on our original field study data (r = 0.845, p = 8.31 × 10⁻³) though our canopy height map showed a poor correlation (r = 0.563, p = 0.146) with the canopy height estimate based on a previously developed model by Feldpausch et al. (2011). We also confirmed that the created canopy height map showed a similar pattern with the previously developed canopy height map by Simard et al. (2011). It was concluded that the use of the spaceborne LiDAR data provides a sufficient accuracy in estimating the canopy height at regional scale.     

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.