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.     

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

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

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

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

Considering combined or separated roughness and vegetation effects in soil moisture retrievals

For more than six years, the Soil Moisture and Ocean Salinity (SMOS) mission has provided multi angular and full-polarization brightness temperature (TB) measurements at L-band. Geophysical products such as soil moisture (SM) and vegetation optical depth at nadir (τ_nad) are retrieved by an operational algorithm using TB observations at different angles of incidence and polarizations. However, the quality of the retrievals depends on several surface effects, such as vegetation, soil roughness and texture, etc. In the microwave forward emission model used in the retrievals (L-band Microwave Emission Model, L-MEB), soil roughness is modelled with a semi-empirical equation using four main parameters (Q_r, H_r, N_rp, with p = H or V polarizations). At present, these parameters are calibrated with data provided by airborne studies and in situ measurements made at a local scale that is not necessarily representative of the large SMOS footprints (43 km on average) at global scale. In this study, we evaluate the impact of the calibrated values of N_rp and H_r on the SM and τ_nad retrievals based on SMOS TB measurements (SMOS Level 3 product) over the Soil Climate Analysis Network (SCAN) network located in North America over five years (2011–2015). In this study, Q_r was set equal to zero and we assumed that N_rH = N_rV. The retrievals were performed by varying N_rp from −1 to 2 by steps of 1 and H_r from 0 to 0.6 by steps of 0.1. At satellite scale, the results show that combining vegetation and roughness effects in a single parameter provides the best results in terms of soil moisture retrievals, as evaluated against the in situ SM data. Even though our retrieval approach was very simplified, as we did not account for pixel heterogeneity, the accuracy we obtained in the SM retrievals was almost systematically better than those of the Level 3 product. Improved results were also obtained in terms of optical depth retrievals. These new results may have key consequences in terms of calibration of roughness effects within the algorithms of the SMOS (ESA) and the SMAP (NASA) space missions.     

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.     

Reducing background effects in orchards through spectral vegetation index correction

Satellite remote sensing provides an alternative to time-consuming and labor intensive in situ measurements of biophysical variables in agricultural crops required for precision agriculture applications. In orchards, however, the spatial resolution causes mixtures of canopies and background (i.e. soil, grass and shadow), hampering the estimation of these biophysical variables. Furthermore, variable background mixtures obstruct meaningful comparisons between different orchard blocks, rows or within each row. Current correction methodologies use spectral differences between canopies and background, but struggle with a vegetated orchard floor. This background influence and the lack of a generic solution are addressed in this study.Firstly, the problem was demonstrated in a controlled environment for vegetation indices sensitive to chlorophyll content, water content and leaf area index. Afterwards, traditional background correction methods (i.e. soil-adjusted vegetation indices and signal unmixing) were compared to the proposed vegetation index correction. This correction was based on the mixing degree of each pixel (i.e. tree cover fraction) to rescale the vegetation indices accordingly and was applied to synthetic and WorldView-2 satellite imagery. Through the correction, the effect of background admixture for vegetation indices was reduced, and the estimation of biophysical variables was improved (ΔR² = 0.2–0.31).     

Determination of the factors governing soil erodibility using hyperspectral visible and near-infrared reflectance spectroscopy

Soil erodibility, which is difficult to estimate and upscaling, was determined in this study using multiple spectral models of soil properties (soil organic matter (SOM), water-stable aggregates (WSA) > 0.25 mm, the geometric mean radius (D_g)). Herein, the soil erodibility indicators were calculated, and soil properties were quantitatively analyzed based on laboratory simulation experiments involving two selected contrasting soils. In addition, continuous wavelet transformation was applied to the reflectance spectra (350–2500 nm) of 65 soil samples from the study area. To build the relationship, the soil properties that control erodibility were identified prior to the spectral analysis. In this study, the SOM, D_g and WSA >0.25 mm were selected to represent the most significant soil properties controlling erodibility and describe the erodibility indicator based on a logarithmic regression model as a function of SOM or WSA > 0.25 mm. Five, six and three wavelet features were observed to calibrate the estimated soil properties model, and the best performance was obtained with a combination feature regression model for SOM (R² = 0.86, p < 0.01), D_g (R² = 0.79, p < 0.01) and WSA >0.25 mm (R² = 0.61, p < 0.01), respectively. One part of the wavelet features captured amplitude variations in the broad shape of the reflectance spectra, and another part captured variations in the shape and depth of the soil dry substances. The wavelet features for the validated dataset used to predict the SOM, WSA >0.25 mm and D_g were not significantly different compared with the calibrated dataset. The synthesized spectral models of soil properties, and the formation of a new equation for soil erodibility transformed from the spectral models of soil properties are presented in this study. These results show that a spectral analytical approach can be applied to complex datasets and provide new insights into emerging dynamic variation with erodibility estimation.     

DisPATCh as a tool to evaluate coarse-scale remotely sensed soil moisture using localized in situ measurements: Application to SMOS and AMSR-E data in Southeastern Australia

Validating coarse-scale satellite soil moisture data still represents a big challenge, notably due to the large mismatch existing between the spatial resolution (> 10 km) of microwave radiometers and the representativeness scale (several m) of localized in situ measurements. This study aims to examine the potential of DisPATCh (Disaggregation based on Physical and Theoretical scale Change) for validating SMOS (Soil Moisture and Ocean Salinity) and AMSR-E (Advanced Microwave Scanning Radiometer-Earth observation system) level-3 soil moisture products. The ∽40–50 km resolution SMOS and AMSR-E data are disaggregated at 1 km resolution over the Murrumbidgee catchment in Southeastern Australia during a one year period in 2010–2011, and the satellite products are compared with the in situ measurements of 38 stations distributed within the study area. It is found that disaggregation improves the mean difference, correlation coefficient and slope of the linear regression between satellite and in situ data in 77%, 92% and 94% of cases, respectively. Nevertheless, the downscaling efficiency is lower in winter than during the hotter months when DisPATCh performance is optimal. Consistently, better results are obtained in the semi-arid than in a temperate zone of the catchment. In the semi-arid Yanco region, disaggregation in summer increases the correlation coefficient from 0.63 to 0.78 and from 0.42 to 0.71 for SMOS and AMSR-E in morning overpasses and from 0.37 to 0.63 and from 0.47 to 0.73 for SMOS and AMSR-E in afternoon overpasses, respectively. DisPATCh has strong potential in low vegetated semi-arid areas where it can be used as a tool to evaluate coarse-scale remotely sensed soil moisture by explicitly representing the sub-pixel variability.     

Using multi-polarization C- and L-band synthetic aperture radar to estimate biomass and soil moisture of wheat fields

Biomass and soil moisture are two important parameters for agricultural crop monitoring and yield estimation. In this study, the Water Cloud Model (WCM) was coupled with the Ulaby soil moisture model to estimate both biomass and soil moisture for spring wheat fields in a test site in western Canada. This study exploited both C-band (RADARSAT-2) and L-band (UAVSAR) Synthetic Aperture Radars (SARs) for this purpose. The WCM-Ulaby model was calibrated for three polarizations (HH, VV and HV). Subsequently two of these three polarizations were used as inputs to an inversion procedure, to retrieve either soil moisture or biomass without the need for any ancillary data. The model was calibrated for total canopy biomass, the biomass of only the wheat heads, as well as for different wheat growth stages. This resulted in a calibrated WCM-Ulaby model for each sensor-polarization-phenology-biomass combination. Validation of model retrievals led to promising results. RADARSAT-2 (HH-HV) estimated total wheat biomass with root mean square (RMSE) and mean average (MAE) errors of 78.834 g/m² and 58.438 g/m²; soil moisture with errors of 0.078 m³/m³ (RMSE) and 0.065 m³/m³ (MAE) are reported. During the period of crop ripening, L-band estimates of soil moisture had accuracies of 0.064 m³/m³ (RMSE) and 0.057 m³/m³ (MAE). RADARSAT-2 (VV-HV) produced interesting results for retrieval of the biomass of the wheat heads. In this particular case, the biomass of the heads was estimated with accuracies of 38.757 g/m² (RSME) and 33.152 g/m² (MAE). For wider implementation this model will require additional data to strengthen the model accuracy and confirm estimation performance. Nevertheless this study encourages further research given the importance of wheat as a global commodity, the challenge of cloud cover in optical monitoring and the potential of direct estimation of the weight of heads where wheat production lies.     

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.     

Modeling soil parameters using hyperspectral image reflectance in subtropical coastal wetlands

Developing spectral models of soil properties is an important frontier in remote sensing and soil science. Several studies have focused on modeling soil properties such as total pools of soil organic matter and carbon in bare soils. We extended this effort to model soil parameters in areas densely covered with coastal vegetation. Moreover, we investigated soil properties indicative of soil functions such as nutrient and organic matter turnover and storage. These properties include the partitioning of mineral and organic soil between particulate (>53 μm) and fine size classes, and the partitioning of soil carbon and nitrogen pools between stable and labile fractions. Soil samples were obtained from Avicennia germinans mangrove forest and Juncus roemerianus salt marsh plots on the west coast of central Florida. Spectra corresponding to field plot locations from Hyperion hyperspectral image were extracted and analyzed. The spectral information was regressed against the soil variables to determine the best single bands and optimal band combinations for the simple ratio (SR) and normalized difference index (NDI) indices. The regression analysis yielded levels of correlation for soil variables with R² values ranging from 0.21 to 0.47 for best individual bands, 0.28 to 0.81 for two-band indices, and 0.53 to 0.96 for partial least-squares (PLS) regressions for the Hyperion image data. Spectral models using Hyperion data adequately (RPD > 1.4) predicted particulate organic matter (POM), silt + clay, labile carbon (C), and labile nitrogen (N) (where RPD = ratio of standard deviation to root mean square error of cross-validation [RMSECV]). The SR (0.53 μm, 2.11 μm) model of labile N with R² = 0.81, RMSECV= 0.28, and RPD = 1.94 produced the best results in this study. Our results provide optimism that remote-sensing spectral models can successfully predict soil properties indicative of ecosystem nutrient and organic matter turnover and storage, and do so in areas with dense canopy cover.     

Glacier surface velocity estimation using repeat TerraSAR-X images: Wavelet- vs. correlation-based image matching

For the observation and monitoring of glacier surface velocity (GSV), remote sensing is an increasingly suitable tool thanks to the high temporal and spatial resolution of the data. Radar sensors have the specific advantage over optical sensors of being nearly weather and time-independent.Two image pairs separated by 11 days, acquired with the high-resolution spotlight (HS) and stripmap (SM) modes of the German sensor TerraSAR-X, were used to estimate GSV over Switzerland’s Aletsch Glacier. The SM mode covers larger ground swaths, making it more suitable for glacier-wide observations, while the HS images cover less area but offer the highest-possible spatial resolution, approximately 1 × 1 m on the ground. The images were acquired during the summer to maximise feature visibility by minimal snow cover.GSV estimation was performed using two methods, the comparison of which was a major goal of this study: traditional cross-correlation optimisation and a dense image matching algorithm based on complex wavelet decomposition. Each method was found to have unique advantages and disadvantages, but it was concluded that for GSV monitoring, cross-correlation is probably preferable to the wavelet-based approach. While it generates fewer estimates per unit area, this is not necessarily a critical requirement for all glaciological applications, and the method requires less initial “tuning” (calibration) than the wavelet algorithm, making it a slightly better tool in operational contexts. Also, the use of the highest-resolution spotlight datasets is recommended over stripmap mode images when large-area coverage is less critical. The comparative lack of visible features at the resolution of the stripmap images made reliable GSV estimation difficult, with the exception of several small areas dominated by large crevasses.     

Comparison of UAV and WorldView-2 imagery for mapping leaf area index of mangrove forest

Unmanned Aerial Vehicle (UAV) remote sensing has opened the door to new sources of data to effectively characterize vegetation metrics at very high spatial resolution and at flexible revisit frequencies. Successful estimation of the leaf area index (LAI) in precision agriculture with a UAV image has been reported in several studies. However, in most forests, the challenges associated with the interference from a complex background and a variety of vegetation species have hindered research using UAV images. To the best of our knowledge, very few studies have mapped the forest LAI with a UAV image. In addition, the drawbacks and advantages of estimating the forest LAI with UAV and satellite images at high spatial resolution remain a knowledge gap in existing literature. Therefore, this paper aims to map LAI in a mangrove forest with a complex background and a variety of vegetation species using a UAV image and compare it with a WorldView-2 image (WV2).In this study, three representative NDVIs, average NDVI (AvNDVI), vegetated specific NDVI (VsNDVI), and scaled NDVI (ScNDVI), were acquired with UAV and WV2 to predict the plot level (10 × 10 m) LAI. The results showed that AvNDVI achieved the highest accuracy for WV2 (R² = 0.778, RMSE = 0.424), whereas ScNDVI obtained the optimal accuracy for UAV (R² = 0.817, RMSE = 0.423). In addition, an overall comparison results of the WV2 and UAV derived LAIs indicated that UAV obtained a better accuracy than WV2 in the plots that were covered with homogeneous mangrove species or in the low LAI plots, which was because UAV can effectively eliminate the influence from the background and the vegetation species owing to its high spatial resolution. However, WV2 obtained a slightly higher accuracy than UAV in the plots covered with a variety of mangrove species, which was because the UAV sensor provides a negative spectral response function(SRF) than WV2 in terms of the mangrove LAI estimation.     

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

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

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.     

Determining iron content in Mediterranean soils in partly vegetated areas,using spectral reflectance and imaging spectroscopy

The possibility of quantifying iron content in the topsoil of the slopes of the El Hacho Mountain complex in Southern Spain using imaging spectroscopy is investigated. Laboratory, field and airborne spectrometer (ROSIS) data are acquired, in combination with soil samples, which are analysed for dithionite extractable iron (Fe_d) content. Analysis of the properties of two iron related absorption features present in laboratory spectra demonstrates good relations, especially between the standard deviation (S.D.) of the values in an absorption feature and the Fe_d content (R² = 0.67) as well as the ratio based Redness Index (R² = 0.51). Such derived relations are less strong for the ROSIS data (R² for S.D. = 0.26 and R² for Redness Index = 0.22). The spatial distribution of iron in vegetated areas shows a strong sensitivity of these relations with the presence of vegetation. A combination of both methods shows that the overestimation of the Fe_d content with the one method is (partly) compensated by the underestimation with the other method.     

Estimation of rice grain yield from dual-polarization Radarsat-2 SAR data by integrating a rice canopy scattering model and a genetic algorithm

Fast and accurate estimation of rice yield plays a role in forecasting rice productivity for ensuring regional or national food security. Microwave synthetic aperture radar (SAR) data has been proved to have a great potential for rice monitoring and parameters retrieval. In this study, a rice canopy scattering model (RCSM) was revised and then was applied to simulate the backscatter of rice canopy. The combination of RCSM and genetic algorithm (GA) was proposed for retrieving two important rice parameters relating to grain yield, ear length and ear number density, from a C-band, dual-polarization (HH and HV) Radarsat-2 SAR data. The stability of retrieved results of GA inversion was also evaluated by changing various parameter configurations.Results show that RCSM can effectively simulate backscattering coefficients of rice canopy at HH and HV mode with an error of <1 dB. Reasonable selection of GA’s parameters is essential for stability and efficiency of rice parameter retrieval. Two rice parameters are retrieved by the proposed RCSM-GA technology with better accuracy. The rice ear length are estimated with error of <1.5 cm, and ear number density with error of <23 #/m². Rice grain yields are effectively estimated and mapped by the retrieved ear length and number density via a simple yield regression equation. This study further illustrates the capability of C-band Radarsat-2 SAR data on retrieval of rice ear parameters and the practicability of radar remote sensing technology for operational yield estimation.     

Quantifying post-fire fallen trees using multi-temporal lidar

Massive tree-felling due to root damage is a common fire effect on burnt areas in Scandinavia, but has so far not been analyzed in detail. Here we explore if pre- and post-fire lidar data can be used to estimate the proportion of fallen trees. The study was carried out within a large (14,000 ha) area in central Sweden burnt in August 2014, where we had access to airborne lidar data from both 2011 and 2015. Three data-sets of predictor variables were tested: POST (post-fire lidar metrics), DIF (difference between post- and pre-fire lidar metrics) and combination of those two (POST_DIF). Fractional logistic regression was used to predict the proportion of fallen trees. Training data consisted of 61 plots, where the number of fallen and standing trees was calculated both in the field and with interpretation of drone images. The accuracy of the best model was tested based on 100 randomly selected validation plots with a size of 25 × 25 m.Our results showed that multi-temporal lidar together with field-collected training data can be used for quantifying post-fire tree felling over large areas. Several height-, density- and intensity metrics correlated with the proportion of fallen trees. The best model combined metrics from both datasets (POST_DIF), resulting in a RMSE of 0.11. Results were slightly poorer in the validation plots with RMSE of 0.18 using pixel size of 12.5 m and RMSE of 0.15 using pixel size of 6.25 m. Our model performed least well for stands that had been exposed to high-intensity crown fire. This was likely due to the low amount of echoes from the standing black tree skeletons. Wall-to-wall maps produced with this model can be used for landscape level analysis of fire effects and to explore the relationship between fallen trees and forest structure, soil type, fire intensity or topography.     

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.     

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.