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.     

Global rain-fed,irrigated, and paddy croplands: A new high resolution map derived from remote sensing,crop inventories and climate data

Irrigation accounts for 70% of global water use by humans and 33–40% of global food production comes from irrigated croplands. Accurate and timely information related to global irrigation is therefore needed to manage increasingly scarce water resources and to improve food security in the face of yield gaps, climate change and extreme events such as droughts, floods, and heat waves. Unfortunately, this information is not available for many regions of the world. This study aims to improve characterization of global rain-fed, irrigated and paddy croplands by integrating information from national and sub-national surveys, remote sensing, and gridded climate data sets. To achieve this goal, we used supervised classification of remote sensing, climate, and agricultural inventory data to generate a global map of irrigated, rain-fed, and paddy croplands. We estimate that 314 million hectares (Mha) worldwide were irrigated circa 2005. This includes 66 Mha of irrigated paddy cropland and 249 Mha of irrigated non-paddy cropland. Additionally, we estimate that 1047 Mha of cropland are managed under rain-fed conditions, including 63 Mha of rain-fed paddy cropland and 985 Mha of rain-fed non-paddy cropland. More generally, our results show that global mapping of irrigated, rain-fed, and paddy croplands is possible by combining information from multiple data sources. However, regions with rapidly changing irrigation or complex mixtures of irrigated and non-irrigated crops present significant challenges and require more and better data to support high quality mapping of irrigation.     

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.     

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.     

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.     

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.     

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.     

Discharge and suspended sediment flux estimated along the mainstream of the Amazon and the Madeira Rivers (from in situ and MODIS Satellite Data)

Water and suspended sediment fluxes are considered during the period 2000–2008 in a region including the full Amazon River from the confluence of the Negro River to Santarém, the end part of the Solimões River, and the lower part of the Madeira River. Three types of data are used: water discharge estimated from field measurements, and suspended sediment obtained from field measurements and derived from MODIS satellite data. A generalized least square method including a propagating term is developed in order to propagate the signal upward and downward the river. The approach is introduced and tested. Several experiments are considered in order, first, to estimate the ability to propagate the signal from stations located before the confluences of Negro and Madeira Rivers to stations located on the Amazon River; second to investigate the possibility to propagate the signal along the Amazon River which dynamics is coupled with floodplains dynamics; and third to produce optimal solutions of water and sediment fluxes. For each experiment, the influence of field and satellite data is compared. The approach is efficient in the upper part of the region of study where the Solimões, the Negro and the Madeira Rivers meet and fails in the lower part of the region where interactions between Amazon River and floodplains play an important role on the fluxes’ dynamics. The optimal experiment includes in situ and satellite data from all the stations available and is used to analyse the recent evolution of suspended sediment flux along the Amazon River and its interaction with the large coupled floodplains. A high accumulation rate is observed during the 2000–2002 period, followed by decreasing rates until 2005 and by increasing values in 2006 and 2007. Our results suggest that floodplains extending along a river reach of 390 km-long between Itacoatiara and Óbidos trap about 15% of the suspended sediment flux passing at Óbidos. The simulated deposition rate is of about 0.3 Mt km⁻¹ yr⁻¹ corresponding to an accretion rate of about 27 mm yr⁻¹.     

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.     

Using remote sensing to monitor the influence of river discharge on watershed outlets and adjacent coral Reefs: Magdalena River and Rosario Islands,Colombia

Worldwide, coral reef ecosystems are being increasingly threatened by sediments loads from river discharges, which in turn are influenced by changing rainfall patterns due to climate change and by growing human activity in their watersheds. In this case study, we explored the applicability of using remote sensing (RS) technology to estimate and monitor the relationship between water quality at the coral reefs around the Rosario Islands, in the Caribbean Sea, and the rainfall patterns in the Magdalena River watershed. From the Moderate Resolution Imaging Spectroradiometer (MODIS), this study used the water surface reflectance product (MOD09GQ) to estimate water surface reflectance as a proxy for sediment concentration and the land cover product (MCD12Q1 V51) to characterize land cover of the watershed. Rainfall was estimated by using the 3B43 V7 product from the Tropical Rainforest Measuring Mission (TRMM). For the first trimester of each year, we investigated the inter-annual temporal variation in water surface reflectance at the Rosario Islands and at the three main mouths of the Magdalena River watershed. No increasing or decreasing trends of water surface reflectance were detected for any of the sites for the study period 2001–2014 (p > 0.05) but significant correlations were detected among the trends of each site at the watershed mouths (r = 0.57–0.90, p < 0.05) and between them and the inter-annual variation in rainfall on the watershed (r = 0.63–0.67, p < 0.05). Those trimesters with above-normal water surface reflectance at the mouths and above-normal rainfall at the watershed coincided with La Niña conditions while the opposite was the case during El Niño conditions. Although, a preliminary analysis of inter-annual land cover trends found only cropland cover in the watershed to be significantly correlated with water surface reflectance at two of the watershed mouths (r = 0.58 and 0.63, p < 0.05), the validation analysis draw only a 40.7% of accuracy in this land cover classification. This requires further analysis to confirm the impact of the cropland on the water quality at the watershed outlets. Spatial analysis with MOD09GQ imagery detected the overpass of river plumes from Barbacoas Bay over the Rosario Islands waters.     

A GPP assimilation model for the southeastern Tibetan Plateau based on CO2 eddy covariance flux tower and remote sensing data

The gross primary production (GPP) at individual CO₂ eddy covariance flux tower sites (GPP_Tower) in Dali (DL), Wenjiang (WJ) and Linzhi (LZ) around the southeastern Tibetan Plateau were determined by the net ecosystem exchange of CO₂ (NEE) and ecosystem respiration (R_e). The satellite remote sensing-VPM model estimates of GPP values (GPP_MODIS) used the satellite-derived 8-day surface reflectance product (MOD09A1), including satellite-derived enhanced vegetation index (EVI) and land surface water index (LSWI). In this paper, we assembled a subset of flux tower data at these three sites to calibrate and test satellite-VPM model estimated GPP_MODIS, and introduced the satellite data and site-level environmental factors to develop four new assimilation models. The new assimilation models’ estimates of GPP values were compared with GPP_MODIS and GPP_Tower, and the final optimum model among the four assimilation models was determined and used to calibrate GPP_MODIS. The results showed that GPP_MODIS had similar temporal variations to the GPP_Tower, but GPP_MODlS were commonly higher in absolute magnitude than GPP_Tower with relative error (RE) about 58.85%. While, the assimilation models’ estimates of GPP values (GPP_MODEL) were much more closer to GPP_Tower with RE approximately 6.98%, indicating that the capacity of the simulation in the new assimilation model was greatly improved, the R² and root mean square error (RMSE) of the new assimilation model were 0.57–4.90% higher and 0.74–2.47 g C m⁻² s⁻¹ lower than those of the GPP_MODIS, respectively. The assimilation model was used to predicted GPP dynamics around the Tibetan Plateau and showed a reliable result compared with other researches. This study demonstrated the potential of the new assimilation model for estimating GPP around the Tibetan Plateau and the performances of site-level biophysical parameters in related to satellite-VPM model GPP.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Validation of the ESA CCI soil moisture product in China

The quality of a newly merged soil moisture product (ECV_SM v0.1) from active and passive microwave sensors has attracted widespread international attention. The performance evaluation of this product will benefit studies on climate, meteorology, agriculture, hydrology, ecology and the environment. In this study, meteorological station data and the Noah soil moisture product were used to validate the ECV_SM product in China. First, some conventional statistical measures, such as correlation coefficients, bias, root mean square difference (RMSD) and mean relative error (MRE), were computed to describe the level of agreement between the meteorological station data and ECV_SM values. The accuracy was moderately high (the correlation was significant at the 0.05 level), although the two datasets differed slightly for various types of land cover. Compared with cropland and urban and built-up areas, the performance of ECV_SM was best in grassland regions. Second, the triple collocation technique was used to assess the random error in the meteorological station data, Noah soil moisture product and ECV_SM product. The mean errors in these three datasets were 0.108, 0.079 and 0.075 m³ m⁻³, respectively, on July 8, 2010 and 0.099, 0.061 and 0.059 m³ m⁻³, respectively, on October 8, 2010. Only two days of data were used for the triple collocation test as a representative, but this cannot precisely indicate that the test results on any other day correspond with the test results on these two days. Additionally, a trend analysis of ECV_SM during 2003–2010 was carried out using the Mann–Kendall trend test.