Two new soil moisture indices based on the NIR-red triangle space of Landsat-8 data

In this study, the NIR-red spectral space of Landsat-8 images, which is manifested by a triangle shape, is deployed for developing two new Soil Moisture (SM) indices. First, ten parameters consisting of six distances and four angles were extracted using the position of a random pixel in this triangle. Then, some correlation assessments were made to derive those parameters that were useful for SM estimation, which were five parameters. To build a soil moisture index, all combinations of these five parameters, which were in total 31 different regression equations, were considered, and the best model was named the Triangle Soil Moisture Index (TSMI). The TSMI consists of three parameters. It showed a RMSE of 0.08 and correlation coefficient (R) of 0.67. Since the TSMI does not consider vegetation interface in SM estimation, the Modified TSMI (MTSMI), which takes into account the fraction of soil cover in each pixel, beside those parameters which were used in the TSMI, was developed (MTSMI: RMSE = 0.07, R = 0.74). The results of the TSMI and MTSMI were compared with each other, and with another soil moisture index (SMMRS introduced by Zhan et al. (2007)). It was concluded that the TSMI and MTSMI provide similar results for bare soil or sparsely vegetated surfaces. However, the MTSMI demonstrated a much better performance in densely vegetated surfaces. The accuracy of both the TSMI and MTSMI were significantly higher than the SMMRS. Moreover, the TSMI and MTSMI were validated by comparison with field measured SM data at five different depths. The results showed that satellite estimated SM by these two indices was more correlated with in situ data at 5 cm soil depth compared to other depths. Also, to show the high applicability of the proposed approach for SM estimation, we selected another set of field SM data collected in Australia. The results proved the effectiveness of the method in different study areas.     

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).     

Satellite surface soil moisture from SMOS and Aquarius: Assessment for applications in agricultural landscapes

Satellite surface soil moisture has become more widely available in the past five years, with several missions designed specifically for soil moisture measurement now available, including the Soil Moisture and Ocean Salinity (SMOS) mission and the Soil Moisture Active/Passive (SMAP) mission. With a wealth of data now available, the challenge is to understand the skill and limitations of the data so they can be used routinely to support monitoring applications and to better understand environmental change. This paper examined two satellite surface soil moisture data sets from the SMOS and Aquarius missions against in situ networks in largely agricultural regions of Canada. The data from both sensors was compared to ground measurements on both an absolute and relative basis. Overall, the root mean squared errors for SMOS were less than 0.10 m³ m⁻³ at most sites, and less where the in situ soil moisture was measured at multiple sites within the radiometer footprint (sites in Saskatchewan, Manitoba and Ontario). At many sites, SMOS overestimates soil moisture shortly after rainfall events compared to the in situ data; however this was not consistent for each site and each time period. SMOS was found to underestimate drying events compared to the in situ data, however this observation was not consistent from site to site. The Aquarius soil moisture data showed higher root mean squared errors in areas where there were more frequent wetting and drying cycles. Overall, both data sets, and SMOS in particular, showed a stable and consistent pattern of capturing surface soil moisture over time.     

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.     

Soil moisture retrieval from airborne L-band passive microwave using high resolution multispectral data

For the soil moisture retrieval from passive microwave sensors, such as ESA’s Soil Moisture and Ocean Salinity (SMOS) and the NASA Soil Moisture Active and Passive (SMAP) mission, a good knowledge about the vegetation characteristics is indispensable. Vegetation cover is a principal factor in the attenuation, scattering and absorption of the microwave emissions from the soil; and has a direct impact on the brightness temperature by way of its canopy emissions. Here, brightness temperatures were measured at three altitudes across the TERENO (Terrestrial Environmental Observatories) Rur catchment site in Germany to achieve a range of spatial resolutions using the airborne Polarimetric L-band Multibeam Radiometer 2 (PLMR2). The L-band Microwave Emission of the Biosphere (L-MEB) model which simulates microwave emissions from the soil–vegetation layer at L-band was used to retrieve surface soil moisture for all resolutions. A Monte Carlo approach was developed to simultaneously estimate soil moisture and the vegetation parameter b’ describing the relationship between the optical thickness τ and the Leaf Area Index (LAI). LAI was retrieved from multispectral RapidEye imagery and the plant specific vegetation parameter b′ was estimated from the lowest flight altitude data for crop, grass, coniferous forest, and deciduous forest. Mean values of b’ were found to be 0.18, 0.07, 0.26 and 0.23, respectively. By assigning the estimated b′ to higher flight altitude data sets, a high accuracy soil moisture retrieval was achieved with a Root Mean Square Difference (RMSD) of 0.035 m³ m⁻³ when compared to ground-based measurements.     

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.     

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.     

Combination of multi-sensor remote sensing data for drought monitoring over Southwest China

Drought is one of the most frequent climate-related disasters occurring in Southwest China, where the occurrence of drought is complex because of the varied landforms, climates and vegetation types. To monitor the comprehensive information of drought from meteorological to vegetation aspects, this paper intended to propose the optimized meteorological drought index (OMDI) and the optimized vegetation drought index (OVDI) from multi-source satellite data to monitor drought in three bio-climate regions of Southwest China. The OMDI and OVDI were integrated with parameters such as precipitation, temperature, soil moisture and vegetation information, which were derived from Tropical Rainfall Measuring Mission (TRMM), Moderate Resolution Imaging Spectroradiometer Land Surface Temperature (MODIS LST), AMSR-E Soil Moisture (AMSR-E SM), the soil moisture product of China Land Soil Moisture Assimilation System (CLSMAS), and MODIS Normalized Difference Vegetation Index (MODIS NDVI), respectively. Different sources of satellite data for one parameter were compared with in situ drought indices in order to select the best data source to derive the OMDI and OVDI. The Constrained Optimization method was adopted to determine the optimal weights of each satellite-based index generating combined drought indices. The result showed that the highest positive correlation and lowest root mean square error (RMSE) between the OMDI and 1-month standardized precipitation evapotranspiration index (SPEI-1) was found in three regions of Southwest China, suggesting that the OMDI was a good index in monitoring meteorological drought; in contrast, the OVDI was best correlated to 3-month SPEI (SPEI-3), and had similar trend with soil relative water content (RWC) in temporal scale, suggesting it a potential indicator of agricultural drought. The spatial patterns of OMDI and OVDI along with the comparisons of SPEI-1 and SPEI-3 for different months in one year or one month in different years showed significantly varied drought locations and areas, demonstrating regional and seasonal fluctuations, and suggesting that drought in Southwest China should be monitored in seasonal and regional level, and more fine distinctions of seasons and regions need to be considered in the future studies of this area.     

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 salinity in Pingluo County of China using QuickBird data and soil reflectance spectra

Soil salinization is a worldwide environmental problem with severe economic and social consequences. In this paper, estimating the soil salinity of Pingluo County, China by a partial least squares regression (PLSR) predictive model was carried out using QuickBird data and soil reflectance spectra. At first, a relationship between the sensitive bands of soil salinity acquired from measured reflectance spectra and the spectral coverage of seven commonly used optical sensors was analyzed. Secondly, the potentiality of QuickBird data in estimating soil salinity by analyzing the correlations between the measured reflectance spectra and reflectance spectra derived from QuickBird data and analyzing the contributions of each band of QuickBird data to soil salinity estimation Finally, a PLSR predictive model of soil salinity was developed using reflectance spectra from QuickBird data and eight spectral indices derived from QuickBird data. The results indicated that the sensitive bands covered several bands of each optical sensor and these sensors can be used for soil salinity estimation. The result of estimation model showed that an accurate prediction of soil salinity can be made based on the PLSR method (R² = 0.992, RMSE = 0.195). The PLSR model's performance was better than that of the stepwise multiple regression (SMR) method. The results also indicated that using spectral indices such as intensity within spectral bands (Int1, Int2), soil salinity indices (SI1, SI2, SI3), the brightness index (BI), the normalized difference vegetation index (NDVI) and the ratio vegetation index (RVI) as independent model variables can help to increase the accuracy of soil salinity mapping. The NDVI and RVI can help to reduce the influences of vegetation cover and soil moisture on prediction accuracy. The method developed in this paper can be applied in other arid and semi-arid areas, such as western China.     

Soil moisture variability over Odra watershed: Comparison between SMOS and GLDAS data

Monitoring of temporal and spatial soil moisture variability is an important issue, both from practical and scientific point of view. It is well known that passive, L-band, radiometric measurements provide best soil moisture estimates. Unfortunately as it was observed during Soil Moisture and Ocean Salinity (SMOS) mission, which was specially dedicated to measure soil moisture, these measurements suffer significant data loss. It is caused mainly by radio frequency interference (RFI) which strongly contaminates Central Europe and even in particularly unfavorable conditions, might prevent these data from being used for regional or watershed scale analysis. Nevertheless, it is highly awaited by researchers to receive statistically significant information on soil moisture over the area of a big watershed. One of such watersheds, the Odra (Oder) river watershed, lies in three European countries – Poland, Germany and the Czech Republic. The area of the Odra river watershed is equal to 118,861 km² making it the second most important river in Poland as well as one of the most significant one in Central Europe.This paper examines the SMOS soil moisture data in the Odra river watershed in the period from 2010 to 2012. This attempt was made to check the possibility of assessing, from the low spatial resolution observations of SMOS, useful information that could be exploited for practical aims in watershed scale, for example, in water storage models even while moderate RFI takes place. Such studies, performed over the area of a large watershed, were recommended by researchers in order to obtain statistically significant results. To meet these expectations, Centre Aval de Traitement des Donnes SMOS (CATDS), 3-days averaged data, together with Global Land Data Assimilation System (GLDAS) National Centers for Environmental Prediction/Oregon State University/Air Force/Hydrologic Research Lab (NOAH) model 0.25 soil moisture values were used for statistical analyses and mutual comparisons.The results obtained using various statistical tools unveil high scientific potential of CATDS SMOS data to study soil moisture over the Odra river watershed. This was also confirmed by reasonable agreement between results derived from CATDS SMOS Ascending and GLDAS data sets. This agreement was achieved mainly by using these data spatially averaged over the whole watershed area, and for observations performed in the period longer than three-day averaging time. Comparisons of separate three-day data in a given pixel position, or at smaller areas would be difficult because of data gaps. Hence, the results of the work suggest that despite of RFI interferences, SMOS observations can provide effective input for analysis of soil moisture at regional scales. Moreover, it was shown that CATDS SMOS soil moisture data are better correlated with rainfall rate than GLDAS ones.     

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.     

A GIHS-based spectral preservation fusion method for remote sensing images using edge restored spectral modulation

High spatial resolution and spectral fidelity are basic standards for evaluating an image fusion algorithm. Numerous fusion methods for remote sensing images have been developed. Some of these methods are based on the intensity–hue–saturation (IHS) transform and the generalized IHS (GIHS), which may cause serious spectral distortion. Spectral distortion in the GIHS is proven to result from changes in saturation during fusion. Therefore, reducing such changes can achieve high spectral fidelity. A GIHS-based spectral preservation fusion method that can theoretically reduce spectral distortion is proposed in this study. The proposed algorithm consists of two steps. The first step is spectral modulation (SM), which uses the Gaussian function to extract spatial details and conduct SM of multispectral (MS) images. This method yields a desirable visual effect without requiring histogram matching between the panchromatic image and the intensity of the MS image. The second step uses the Gaussian convolution function to restore lost edge details during SM. The proposed method is proven effective and shown to provide better results compared with other GIHS-based methods.     

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.     

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.     

闪电河流域农牧交错带微波遥感土壤水分产品评价

空间网格分辨率为9 km的SMAP (Soil Moisture Active and Passive)、0.1D (Degree)的ASCAT (The Advanced Scatterometer)、 25 km的FY-3B以及25 km ESA-CCI (European Space Agency-Climate Change Initiative)是较为广泛应用的卫星遥感土壤水分产品,对数据质量的评价是进一步应用于旱情监测、蒸散发估算等研究的前提。本研究基于2018年9月在闪电河流域内蒙古农牧交错带区域开展的碳、水循环与能量平衡遥感综合试验,采用近似同步的两种尺度观测数据即点尺度地面实测土壤水分数据以及面尺度(1 km×1 km)机载土壤水分数据,利用RMSE (Root Mean Square Error), MAE (Mean Absolute Error),R (Correlation Coefficient),Bias以及ubRMSE (unbiased Root Mean Square Error)等评价指标分别对SMAP, ASCAT, FY-3B, ESA-CCI土壤水分卫星遥感产品进行了评价。本研究利用机载土壤水分数据作为桥梁,实现了从点尺度地面实测土壤水分数据、至面尺度(1 km×1 km)机载土壤水分数据、再至粗格网面尺度(9 km×9 km、0.1 D×0.1 D、25 km×25 km)卫星遥感土壤水分产品的对比分析过程。利用地面观测值对机载观测土壤水分开展评价分析,发现在裸土区域,机载土壤水分数据与地面实测数据较为一致,RMSE, MAE, Bias, ubRMSE以及R值分别为0.033 cm~3/cm~3,0.030 cm~3/cm~3,-0.004 cm~3/cm~3, 0.033 cm~3/cm~3, 0.474。对卫星土壤水分产品的评价结果显示,SMAP的9 km土壤水分卫星产品与地面观测更为一致,其RMSE,MAE,Bias,ub RMSE以及R值分别为0.037 cm~3/cm~3,0.032 cm~3/cm~3,-0.008 cm~3/cm~3, 0.036 cm~3/cm~3, 0.507。SMAP, ASCAT, FY-3B以及ESA-CCI与机载土壤水分数据有更高的相关性,R值分别为0.735, 0.558, 0.558, 0.575。综上,闪电河流域实验区内的4种卫星遥感土壤水分产品中,SMAP产品与地面土壤水分、机载土壤水分数据均较为一致,其次是FY-3B与ESA-CCI。     

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.     

Spectra and vegetation index variations in moss soil crust in different seasons,and in wet and dry conditions

Similar to vascular plants, non-vascular plant mosses have different periods of seasonal growth. There has been little research on the spectral variations of moss soil crust (MSC) over different growth periods. Few studies have paid attention to the difference in spectral characteristics between wet MSC that is photosynthesizing and dry MSC in suspended metabolism. The dissimilarity of MSC spectra in wet and dry conditions during different seasons needs further investigation. In this study, the spectral reflectance of wet MSC, dry MSC and the dominant vascular plant (Artemisia) were characterized in situ during the summer (July) and autumn (September). The variations in the normalized difference vegetation index (NDVI), biological soil crust index (BSCI) and CI (crust index) in different seasons and under different soil moisture conditions were also analyzed. It was found that (1) the spectral characteristics of both wet and dry MSCs varied seasonally; (2) the spectral features of wet MSC appear similar to those of the vascular plant, Artemisia, whether in summer or autumn; (3) both in summer and in autumn, much higher NDVI values were acquired for wet than for dry MSC (0.6 ∼ 0.7 vs. 0.3 ∼ 0.4 units), which may lead to misinterpretation of vegetation dynamics in the presence of MSC and with the variations in rainfall occurring in arid and semi-arid zones; and (4) the BSCI and CI values of wet MSC were close to that of Artemisia in both summer and autumn, indicating that BSCI and CI could barely differentiate between the wet MSC and Artemisia.     

Integrating age in the detection and mapping of incongruous patches in coffee (Coffea arabica) plantations using multi-temporal Landsat 8 NDVI anomalies

The development of cost-effective, reliable and easy to implement crop condition monitoring methods is urgently required for perennial tree crops such as coffee (Coffea arabica), as they are grown over large areas and represent long term and higher levels of investment. These monitoring methods are useful in identifying farm areas that experience poor crop growth, pest infestation, diseases outbreaks and/or to monitor response to management interventions. This study compares field level coffee mean NDVI and LSWI anomalies and age-adjusted coffee mean NDVI and LSWI anomalies in identifying and mapping incongruous patches across perennial coffee plantations. To achieve this objective, we first derived deviation of coffee pixels from the global coffee mean NDVI and LSWI values of nine sequential Landsat 8 OLI image scenes. We then evaluated the influence of coffee age class (young, mature and old) on Landsat-scale NDVI and LSWI values using a one-way ANOVA and since results showed significant differences, we adjusted NDVI and LSWI anomalies for age-class. We then used the cumulative inverse distribution function (α ≤ 0.05) to identify fields and within field areas with excessive deviation of NDVI and LSWI from the global and the age-expected mean for each of the Landsat 8 OLI scene dates spanning three seasons. Results from accuracy assessment indicated that it was possible to separate incongruous and healthy patches using these anomalies and that using NDVI performed better than using LSWI for both global and age-adjusted mean anomalies. Using the age-adjusted anomalies performed better in separating incongruous and healthy patches than using the global mean for both NDVI (Overall accuracy = 80.9% and 68.1% respectively) and for LSWI (Overall accuracy = 68.1% and 48.9% respectively). When applied to other Landsat 8 OLI scenes, the results showed that the proportions of coffee fields that were modelled incongruent decreased with time for the young age category and while it increased for the mature and old age classes with time. We concluded that the method could be useful for the identification of anomalous patches using Landsat scale time series data to monitor large coffee plantations and provide an indication of areas requiring particular field attention.