Evaluating suitability of MODIS-Terra images for reproducing historic sediment concentrations in water bodies: Lake Tana,Ethiopia

Government and NGO funded conservation programs are being implemented in developing countries with the potential benefit of reduced sediment inflow into fresh water lakes. However, these claims are difficult to verify due to limited historical sediment concentration data in lakes and rivers. Remote sensing can potentially aid in monitoring sediment concentration. With almost daily availability over the past ten years and consistent atmospheric correction applied to the images, Moderate Resolution Imaging Spectroradiometer (MODIS) 250 meter images are potential resources capable of monitoring future concentrations and reconstructing historical sediment concentration records. In this paper, site-specific relationships are developed between reflectance in near-infrared (NIR) images and three factors: total suspended solids (TSS), turbidity and Secchi depth for Lake Tana near the mouth of the Gumara River. The first two sampling campaigns on November 27, 2010 and May 13, 2011 are used in calibration. Reflectance in the NIR varies linearly with turbidity (R² = 0.89) and TSS (R² = 0.95). Secchi depth fit best to an exponential relation with R² of 0.74. The relationships are validated using a third sample set collected on November 7, 2011 with RMSE of 11 Nephelometric Turbidity Units (NTU) for Turbidity, 16.5 mg l⁻¹ for TSS and 0.12 meters for Secchi depth. The MAE was 10% for TSS, 14% for turbidity and 0.1% for Secchi depth. Using the relationship for TSS, a 10-year time series of sediment concentration in Lake Tana near the Gumara River was plotted. It was found that after the severe drought of 2002 and 2003 the concentration in the lake increased significantly. The results showed that MODIS images are potential cost effective tools to monitor suspended sediment concentration and obtain a past history of concentration for evaluating the effect of best management practices.     

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.     

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.     

Assessment of satellite ocean color products of MERIS,MODIS and SeaWiFS along the East China Coast (in the Yellow Sea and East China Sea)

The validation of satellite ocean-color products is an important task of ocean-color missions. The uncertainties of these products are poorly quantified in the Yellow Sea (YS) and East China Sea (ECS), which are well known for their optical complexity and turbidity in terms of both oceanic and atmospheric optical properties. The objective of this paper is to evaluate the primary ocean-color products from three major ocean-color satellites, namely the Moderate Resolution Imaging Spectroradiometer (MODIS), Medium Resolution Imaging Spectrometer (MERIS), and Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Through match-up analysis with in situ data, it is found that satellite retrievals of the spectral remote sensing reflectance R_rs(λ) at the blue-green and green bands from MERIS, MODIS and SeaWiFS have the lowest uncertainties with a median of the absolute percentage of difference (APD_m) of 15–27% and root-mean-square-error (RMS) of 0.0021–0.0039 sr⁻¹, whereas the R_rs(λ) uncertainty at 412 nm is the highest (APD_m 47–62%, RMS 0.0027–0.0041 sr⁻¹). The uncertainties of the aerosol optical thickness (AOT) τ_a, diffuse attenuation coefficient for downward irradiance at 490 nm K_d(490), concentrations of suspended particulate sediment concentration (SPM) and Chlorophyll a (Chl-a) were also quantified. It is demonstrated that with appropriate in-water algorithms specifically developed for turbid waters rather than the standard ones adopted in the operational satellite data processing chain, the uncertainties of satellite-derived properties of K_d(490), SPM, and Chl-a may decrease significantly to the level of 20–30%, which is true for the majority of the study area. This validation activity advocates for (1) the improvement of the atmosphere correction algorithms with the regional aerosol optical model, (2) switching to regional in-water algorithms over turbid coastal waters, and (3) continuous support of the dedicated in situ data collection effort for the validation task.     

Estimation of higher chlorophylla concentrations using field spectral measurement and HJ-1A hyperspectral satellite data in Dianshan Lake,China

Based on in situ water sampling and field spectral measurements in Dianshan Lake, a semi-analytical three-band algorithm was used to estimate Chlorophylla (Chla) content in case II waters. The three bands selected to estimate Chla for high concentrations included 653, 691 and 748 nm. An equation, based on the difference in reciprocal reflectance between 653 and 691 nm, multiplied by reflectance at 748 nm as [R_rs⁻¹(653) − R_rs⁻¹ (691)] R_rs(748), explained 85.57% of variance in Chla concentration with a root mean square error (RMSE) of <6.56 mg/m³. In order to test the utility of this model with satellite data, HJ-1A Hyperspectral Imager (HSI) data were analyzed using comparable wavelengths selected from the in situ data [B₆₇⁻¹(656) − B₈₀⁻¹(716)] B₈₇(753). This model accounted for 84.3% of Chla variation, estimating Chla concentrations with an RMSE of <4.23 mg/m³. The results illustrate that, based on the determined wavelengths, the spectrum-based model can achieve a high estimation accuracy and can be applied to hyperspectral satellite imagery especially for higher Chla concentration waters.     

A radiative transfer model-based method for the estimation of grassland aboveground biomass

This paper presents a novel method to derive grassland aboveground biomass (AGB) based on the PROSAILH (PROSPECT + SAILH) radiative transfer model (RTM). Two variables, leaf area index (LAI, m²m⁻², defined as a one-side leaf area per unit of horizontal ground area) and dry matter content (DMC, gcm⁻², defined as the dry matter per leaf area), were retrieved using PROSAILH and reflectance data from Landsat 8 OLI product. The result of LAI × DMC was regarded as the estimated grassland AGB according to their definitions. The well-known ill-posed inversion problem when inverting PROSAILH was alleviated using ecological criteria to constrain the simulation scenario and therefore the number of simulated spectra. A case study of the presented method was applied to a plateau grassland in China to estimate its AGB. The results were compared to those obtained using an exponential regression, a partial least squares regression (PLSR) and an artificial neural networks (ANN). The RTM-based method offered higher accuracy (R² = 0.64 and RMSE = 42.67 gm⁻²) than the exponential regression (R² = 0.48 and RMSE = 41.65 gm⁻²) and the ANN (R² = 0.43 and RMSE = 46.26 gm⁻²). However, the proposed method offered similar performance than PLSR as presented better determination coefficient than PLSR (R² = 0.55) but higher RMSE (RMSE = 37.79 gm⁻²). Although it is still necessary to test these methodologies in other areas, the RTM-based method offers greater robustness and reproducibility to estimate grassland AGB at large scale without the need to collect field measurements and therefore is considered the most promising methodology.     

A new calibration of the effective scattering albedo and soil roughness parameters in the SMOS SM retrieval algorithm

This study focuses on the calibration of the effective vegetation scattering albedo (ω) and surface soil roughness parameters (H_R, and N_Rp, p = H,V) in the Soil Moisture (SM) retrieval from L-band passive microwave observations using the L-band Microwave Emission of the Biosphere (L-MEB) model. In the current Soil Moisture and Ocean Salinity (SMOS) Level 2 (L2), v620, and Level 3 (L3), v300, SM retrieval algorithms, low vegetated areas are parameterized by ω = 0 and H_R = 0.1, whereas values of ω = 0.06 − 0.08 and H_R = 0.3 are used for forests. Several parameterizations of the vegetation and soil roughness parameters (ω, H_R and N_Rp, p = H,V) were tested in this study, treating SMOS SM retrievals as homogeneous over each pixel instead of retrieving SM over a representative fraction of the pixel, as implemented in the operational SMOS L2 and L3 algorithms. Globally-constant values of ω = 0.10, H_R = 0.4 and N_Rp = −1 (p = H,V) were found to yield SM retrievals that compared best with in situ SM data measured at many sites worldwide from the International Soil Moisture Network (ISMN). The calibration was repeated for collections of in situ sites classified in different land cover categories based on the International Geosphere-Biosphere Programme (IGBP) scheme. Depending on the IGBP land cover class, values of ω and H_R varied, respectively, in the range 0.08–0.12 and 0.1–0.5. A validation exercise based on in situ measurements confirmed that using either a global or an IGBP-based calibration, there was an improvement in the accuracy of the SM retrievals compared to the SMOS L3 SM product considering all statistical metrics (R = 0.61, bias = −0.019 m³ m⁻³, ubRMSE = 0.062 m³ m⁻³ for the IGBP-based calibration; against R = 0.54, bias = −0.034 m³ m⁻³ and ubRMSE = 0.070 m³ m⁻³ for the SMOS L3 SM product). This result is a key step in the calibration of the roughness and vegetation parameters in the operational SMOS retrieval algorithm. The approach presented here is the core of a new forthcoming SMOS optimized SM product.     

Towards decadal soil salinity mapping using Landsat time series data

Salinization is one of the major soil problems around the world. However, decadal variation in soil salinization has not yet been extensively reported. This study exploited thirty years (1985–2015) of Landsat sensor data, including Landsat-4/5 TM (Thematic Mapper), Landsat-7 ETM+ (Enhanced Thematic Mapper Plus) and Landsat-8 OLI (Operational Land Imager), for monitoring soil salinity of the Yellow River Delta, China. The data were initially corrected for atmospheric effects, and then matched the spectral bands of EO-1 (Earth Observing One) ALI (Advanced Land Imager). Subsequently, soil salinity maps were derived with a previously developed PLSR (Partial Least Square Regression) model. On intra-annual scale, the retrievals showed that soil salinity increased in February, stabilized in March, and decreased in April. On inter-annual scale, soil salinity decreased within 1985–2000 (−0.74 g kg⁻¹/10a, p < 0.001), and increased within 2000–2015 (0.79 g kg⁻¹/10a, p < 0.001). Our study presents a new perspective for use of multiple Landsat data in soil salinity retrieval, and further the understanding of soil salinization development over the Yellow River Delta.     

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.     

Uncertainty analysis of remote sensing of colored dissolved organic matter: Evaluations and comparisons for three rivers in North America

The uncertainties involved in remote sensing inversion of CDOM (Colored Dissolved Organic Matter) were analyzed in estuarine and coastal regions of three North American rivers: Mississippi, Hudson, and Neponset. Water optical and biogeochemical properties, including CDOM absorption and above-surface spectra, were collected in very high resolution. CDOM’s concentrations (a_g(440), absorption coefficient at 440 nm) were inverted from EO-1 Hyperion images, using a quasi-analytical algorithm for CDOM (QAA-CDOM). Uncertainties are classified to five levels, in which the underwater measurement uncertainty (level 1), image preprocessing uncertainty (level 4) and inverse model uncertainty (level 5) were evaluated. Results indicate that at level 1, in situ CDOM measurement is significant with 0.1 in the unit of QSU and 0.01 in the unit of a_g(440) (m⁻¹). At level 4, surface wave is a potential uncertainty source for high-resolution images in estuarine and coastal regions. The remote sensing reflectance of wavy water is about 10 times of the truth. At level 5, the overall uncertainty of QAA-CDOM inversion is 0.006 m⁻¹, with accuracy R² = 0.77, k = 1.1 and RMSE_log = 0.33 m⁻¹. The correlations between uncertainties and other water properties indicate that the large uncertainty in some rivers, such as the Neponset and Atchafalaya, might be caused by high-concentration chlorophyll or sediments. The relationships among the three level uncertainties show that the level 1 uncertainty generally does not propagate into level 4 and 5, but the large uncertainty at level 4 usually introduce large uncertainty at level 5.     

Evaluation of forest fire on Madeira Island using Sentinel-2A MSI imagery

A forest fire started on August 8th, 2016 in several places on Madeira Island causing damage and casualties. As of August 10th the local media had reported the death of three people, over 200 people injured, over 950 habitants evacuated, and 50 houses damaged. This study presents the preliminary results of the assessment of several spectral indices to evaluate the burn severity of Madeira fires during August 2016. These spectral indices were calculated using the new European satellite Sentinel-2A launched in June 2015. The study confirmed the advantages of several spectral indices such as Normalized Difference Vegetation Index (NDVI), Green Normalized Difference Vegetation Index (GNDVI), Normalized Burn Ratio (NBR) and Normalized Difference Vegetation Index (NDVIreXn) using red-edge spectral bands to assess the post-fire conditions. Results showed high correlation between NDVI, GNDVI, NBR and NDVIre1n spectral indices and the analysis performed by Copernicus Emergency Management Service (EMSR175), considered as the reference truth. Regarding the red-edge spectral indices, the NDVIre1n (using band B5, 705 nm) presented better results compared with B6 (740 nm) and B7 (783 nm) bands. These preliminary results allow us to assume that Sentinel-2 will be a valuable tool for post-fire monitoring. In the future, the two twin Sentinel-2 satellites will offer global coverage of the Madeira Archipelago every five days, therefore allowing the simultaneous study of the evolution of the burnt area and reforestation information with high spatial (up to 10 m) and temporal resolution (5 days).     

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.     

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.     

Modeling biophysical properties of broad-leaved stands in the hyrcanian forests of Iran using fused airborne laser scanner data and ultraCam-D images

Inventories of mixed broad-leaved forests of Iran mainly rely on terrestrial measurements. Due to rapid changes and disturbances and great complexity of the silvicultural systems of these multilayer forests, frequent repetition of conventional ground-based plot surveys is often cost prohibitive. Airborne laser scanning (ALS) and multispectral data offer an alternative or supplement to conventional inventories in the Hyrcanian forests of Iran. In this study, the capability of a combination of ALS and UltraCam-D data to model stand volume, tree density, and basal area using random forest (RF) algorithm was evaluated. Systematic sampling was applied to collect field plot data on a 150 m × 200 m sampling grid within a 1100 ha study area located at 36°38′- 36°42′N and 54°24′–54°25′E. A total of 308 circular plots (0.1 ha) were measured for calculation of stand volume, tree density, and basal area per hectare. For each plot, a set of variables was extracted from both ALS and multispectral data. The RF algorithm was used for modeling of the biophysical properties using ALS and UltraCam-D data separately and combined. The results showed that combining the ALS data and UltraCam-D images provided a slight increase in prediction accuracy compared to separate modeling. The RMSE as percentage of the mean, the mean difference between observed and predicted values, and standard deviation of the differences using a combination of ALS data and UltraCam-D images in an independent validation at 0.1-ha plot level were 31.7%, 1.1%, and 84 m³ ha⁻¹ for stand volume; 27.2%, 0.86%, and 6.5 m² ha⁻¹ for basal area, and 35.8%, −4.6%, and 77.9 n ha⁻¹ for tree density, respectively. Based on the results, we conclude that fusion of ALS and UltraCam-D data may be useful for modeling of stand volume, basal area, and tree density and thus gain insights into structural characteristics in the complex Hyrcanian forests.     

Winter wheat yield forecasting in Ukraine based on Earth observation,meteorological data and biophysical models

Ukraine is one of the most developed agriculture countries and one of the biggest crop producers in the world. Timely and accurate crop yield forecasts for Ukraine at regional level become a key element in providing support to policy makers in food security. In this paper, feasibility and relative efficiency of using moderate resolution satellite data to winter wheat forecasting in Ukraine at oblast level is assessed. Oblast is a sub-national administrative unit that corresponds to the NUTS2 level of the Nomenclature of Territorial Units for Statistics (NUTS) of the European Union. NDVI values were derived from the MODIS sensor at the 250 m spatial resolution. For each oblast NDVI values were averaged for a cropland map (Rainfed croplands class) derived from the ESA GlobCover map, and were used as predictors in the regression models. Using a leave-one-out cross-validation procedure, the best time for making reliable yield forecasts in terms of root mean square error was identified. For most oblasts, NDVI values taken in April–May provided the minimum RMSE value when comparing to the official statistics, thus enabling forecasts 2–3 months prior to harvest. The NDVI-based approach was compared to the following approaches: empirical model based on meteorological observations (with forecasts in April–May that provide minimum RMSE value) and WOFOST crop growth simulation model implemented in the CGMS system (with forecasts in June that provide minimum RMSE value). All three approaches were run to produce winter wheat yield forecasts for independent datasets for 2010 and 2011, i.e. on data that were not used within model calibration process. The most accurate predictions for 2010 were achieved using the CGMS system with the RMSE value of 0.3 t ha⁻¹ in June and 0.4 t ha⁻¹ in April, while performance of three approaches for 2011 was almost the same (0.5–0.6 t ha⁻¹ in April). Both NDVI-based approach and CGMS system overestimated winter wheat yield comparing to official statistics in 2010, and underestimated it in 2011. Therefore, we can conclude that performance of empirical NDVI-based regression model was similar to meteorological and CGMS models when producing winter wheat yield forecasts at oblast level in Ukraine 2–3 months prior to harvest, while providing minimum requirements to input datasets.     

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.     

Estimating soil zinc concentrations using reflectance spectroscopy

Soil contamination by heavy metals has been an increasingly severe threat to nature environment and human health. Efficiently investigation of contamination status is essential to soil protection and remediation. Visible and near-infrared reflectance spectroscopy (VNIRS) has been regarded as an alternative for monitoring soil contamination by heavy metals. Generally, the entire VNIR spectral bands are employed to estimate heavy metal concentration, which lacks interpretability and requires much calculation. In this study, 74 soil samples were collected from Hunan Province, China and their reflectance spectra were used to estimate zinc (Zn) concentration in soil. Organic matter and clay minerals have strong adsorption for Zn in soil. Spectral bands associated with organic matter and clay minerals were used for estimation with genetic algorithm based partial least square regression (GA-PLSR). The entire VNIR spectral bands, the bands associated with organic matter and the bands associated with clay minerals were incorporated as comparisons. Root mean square error of prediction, residual prediction deviation, and coefficient of determination (R²) for the model developed using combined bands of organic matter and clay minerals were 329.65 mg kg⁻¹, 1.96 and 0.73, which is better than 341.88 mg kg⁻¹, 1.89 and 0.71 for the entire VNIR spectral bands, 492.65 mg kg⁻¹, 1.31 and 0.40 for the organic matter, and 430.26 mg kg⁻¹, 1.50 and 0.54 for the clay minerals. Additionally, in consideration of atmospheric water vapor absorption in field spectra measurement, combined bands of organic matter and absorption around 2200 nm were used for estimation and achieved high prediction accuracy with R² reached 0.640. The results indicate huge potential of soil reflectance spectroscopy in estimating Zn concentrations in soil.     

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.     

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.