Using Imaging Spectroscopy to study soil properties

Imaging Spectroscopy (IS) is a promising tool for studying soil properties in large spatial domains. Going from point to image spectrometry is not only a journey from micro to macro scales, but also a long stage where problems such as dealing with data having a low signal-to-noise level, contamination of the atmosphere, large data sets, the BRDF effect and more are often encountered. In this paper we provide an up-to-date overview of some of the case studies that have used IS technology for soil science applications. Besides a brief discussion on the advantages and disadvantages of IS for studying soils, the following cases are comprehensively discussed: soil degradation (salinity, erosion, and deposition), soil mapping and classification, soil genesis and formation, soil contamination, soil water content, and soil swelling. We review these case studies and suggest that the 15 data be provided to the end-users as real reflectance and not as raw data and with better signal-to-noise ratios than presently exist. This is because converting the raw data into reflectance is a complicated stage that requires experience, knowledge, and specific infrastructures not available to many users, whereas quantitative spectral models require good quality data. These limitations serve as a barrier that impedes potential end-users, inhibiting researchers from trying this technique for their needs. The paper ends with a general call to the soil science audience to extend the utilization of the IS technique, and it provides some ideas on how to propel this technology forward to enable its widespread adoption in order to achieve a breakthrough in the field of soil science and remote sensing. (C) 2009 Elsevier Inc. All rights reserved.

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and interannually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R(2) = 0.80 and R(2) = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages. (C) 2010 Elsevier Inc. All rights reserved.

Spectral behavior of some modal soil profiles from São Paulo State, Brazil

Remote sensing has a high potential for environmental evaluation. However, a necessity exists for a better understanding of the relations between the soil attributes and spectral data. The objective of this work was to analyze the spectral behavior of some soil profiles from the region of Piracicaba, São Paulo State, using a laboratory spectroradiometer (400 to 2500 nm). The relations between the reflected electromagnetic energy and the soil physical, chemical and mineralogical attributes were analyzed, verifying the spectral variations of soil samples in depth along the profiles with their classification and discrimination. Sandy soil reflected more, presenting a spectral curve with an ascendant form, opposite to clayey soils. The 1900 nm band discriminated soil with 2:1 mineralogy from the 1:1 and oxidic soils. It was possible to detect the presence of kaolinite, gibbsite, hematite and goethite in the soils through the descriptive aspects of curves, absorption features and reflectance intensity. A relation exists between the weathering stage and spectral data. The evaluation of the superficial and subsuperficial horizon samples allowed characterizing and discriminating the analytical variability of the profile, helping to soil distinguishing and classification.

Using numerical classification of profiles based on Vis-NIR spectra to distinguish soils from the Piracicaba Region, Brazil

Considering that information from soil reflectance spectra is underutilized in soil classification, this paper aimed to evaluate the relationship of soil physical, chemical properties and their spectra, to identify spectral patterns for soil classes, evaluate the use of numerical classification of profiles combined with spectral data for soil classification. We studied 20 soil profiles from the municipality of Piracicaba, State of São Paulo, Brazil, which were morphologically described and classified up to the 3rd category level of the Brazilian Soil Classification System (SiBCS). Subsequently, soil samples were collected from pedogenetic horizons and subjected to soil particle size and chemical analyses. Their Vis-NIR spectra were measured, followed by principal component analysis. Pearson's linear correlation coefficients were determined among the four principal components and the following soil properties: pH, organic matter, P, K, Ca, Mg, Al, CEC, base saturation, and Al saturation. We also carried out interpretation of the first three principal components and their relationships with soil classes defined by SiBCS. In addition, numerical classification of the profiles based on the OSACA algorithm was performed using spectral data as a basis. We determined the Normalized Mutual Information (NMI) and Uncertainty Coefficient (U). These coefficients represent the similarity between the numerical classification and the soil classes from SiBCS. Pearson's correlation coefficients were significant for the principal components when compared to sand, clay, Al content and soil color. Visual analysis of the principal component scores showed differences in the spectral behavior of the soil classes, mainly among Argissolos and the others soils. The NMI and U similarity coefficients showed values of 0.74 and 0.64, respectively, suggesting good similarity between the numerical and SiBCS classes. For example, numerical classification correctly distinguished Argissolos from Latossolos and Nitossolos. However, this mathematical technique was not able to distinguish Latossolos from Nitossolos Vermelho férricos, but the Cambissolos were well differentiated from other soil classes. The numerical technique proved to be effective and applicable to the soil classification process.

Soil density evaluated by spectral reflectance as an evidence of compaction effects

Soil compaction, reflected by high bulk density, is an environmental degradation process and new technologies are being developed for its detection. Despite the proven efficiency of remote sensing, it has not been widely used for soil density. Our objective was to evaluate the density of two soils: a Typic Quartzpisament (TQ) and a Rhodic Paleudalf (RP), using spectral reflectance obtained by a laboratory spectroradiometer between 450 and 2500 nm. Undisturbed samples were taken at two depths (0-20 and 60-80 cm), and were artificially compacted. Spectral data, obtained before and after compaction, were compared for both wet and dried compacted samples. Results demonstrated that soil density was greater in RP than in TQ at both depths due to its clayey texture. Spectral data detected high density (compacted) from low density (non-compacted) clayey soils under both wet and dry conditions. The detection of density in sandy soils by spectral reflectance was not possible. The intensity of spectral reflectance of high soil bulk density (compacted) samples was higher than for low density (non-compacted) soils due to changes in soil structure and porosity. Dry samples with high bulk density showed differences in the spectral intensity, but not in the absorption features. Wet samples in equal condition had statistically higher reflectance intensity than that of the low soil bulk density (non-compacted), and absorption differences at 1920 nm, which was due to the altered position of the water molecules. Soil line and spectral reflectance used together could detect soil bulk density variations for the clay soil. This technique could assist in the detection of high soil density in the laboratory by providing new soil information.

The combined use of reflectance, emissivity and elevation Aster/Terra data for tropical soil studies

Reflectance, emissivity and elevation data of the sensor ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer)/Terra were used to characterize soil composition variations according to the toposequence position. Normalized data of SWIR (shortwave infrared) reflectance and TIR (thermal infrared) emissivity, coupled to a soil-fraction image from a spectral mixture model, were evaluated to separate bare soils from nonphotosynthetic vegetation. Regression relationships of some soil properties with reflectance and emissivity data were then applied on the exposed soil pixels. The resulting estimated values were plotted on the ASTER-derived digital elevation model. Results showed that the SWIR bands 5 and 6 and the TIR bands 10 and 14 measured the clay mineral absorption band and the quartz emissivity feature, respectively. These bands improved also the discrimination between nonphotosynthetic vegetation and soils. Despite the differences in pixel size and field sampling size, some soil properties were correlated with reflectance (R² of 0.65 for Al2O3 in band 6; 0.61 for Fe2O3 in band 3) and emissivity (R² of 0.65 for total sand fraction in the 10/14 band ratio). The combined use of reflectance, emissivity and elevation data revealed variations in soil composition with topography in specific parts of the landscape. From higher to lower slope positions, a general decrease in Al2O3 and increase in total sand fraction was observed, due to the prevalence of Rhodic Acrustox at the top and its gradual transition to Typic Acrustox at the bottom.

The relationship between diffuse spectral reflectance of the soil and its cation exchange capacity is scale-dependent

Diffuse reflectance spectroscopy (DRS) is increasingly being used to predict numerous soil physical, chemical and biochemical properties. However, soil properties and processes vary at different scales and, as a result, relationships between soil properties often depend on scale. In this paper we report on how the relationship between one such property, cation exchange capacity (CEC), and the DRS of the soil depends on spatial scale. We show this by means of a nested analysis of covariance of soils sampled on a balanced nested design in a 16 km × 16 km area in eastern England. We used principal components analysis on the DRS to obtain a reduced number of variables while retaining key variation. The first principal component accounted for 99.8% of the total variance, the second for 0.14%. Nested analysis of the variation in the CEC and the two principal components showed that the substantial variance components are at the > 2000-m scale. This is probably the result of differences in soil composition due to parent material. We then developed a model to predict CEC from the DRS and used partial least squares (PLS) regression do to so. Leave-one-out cross-validation results suggested a reasonable predictive capability (R2 = 0.71 and RMSE = 0.048 molc kg− 1). However, the results from the independent validation were not as good, with R2 = 0.27, RMSE = 0.056 molc kg− 1 and an overall correlation of 0.52. This would indicate that DRS may not be useful for predictions of CEC. When we applied the analysis of covariance between predicted and observed we found significant scale-dependent correlations at scales of 50 and 500 m (0.82 and 0.73 respectively). DRS measurements can therefore be useful to predict CEC if predictions are required, for example, at the field scale (50 m). This study illustrates that the relationship between DRS and soil properties is scale-dependent and that this scale dependency has important consequences for prediction of soil properties from DRS data

Alteration of soil properties through a weathering sequence as evaluated by spectral reflectance

In the study of physical, chemical, and mineralogical data related to the weathering of soils and the quantification of their properties, remote sensing constitutes an important technique that, in addition to conventional analyses, can contribute to soil survey. The objectives of this research were to characterize and differentiate soils developed from basaltic rocks that occur in the Parana state, Brazil and to quantify soil properties based on their spectral reflectance. These observations were used to verify the relationship between the soils and reflectance with regard to weathering, organic matter (OM), and forms of Fe. From the least to the most weathered soil, we used a Typic Argiudoll (Reddish Brunizem), Rhodudalf (Terra Roxa Estruturada), and Rhodic Hapludox (Very Dark Red Latosol). The spectral reflectances between 400 and 2500 nm were obtained in the laboratory from soil samples collected at two depth increments, 0- to 20- and 40- to 60-cm, using an Infra Red Intelligent Spectroradiometer (IRIS). Correlation, regression, and discriminant estimates were used in analyzing the soil and spectral data. Results of this study indicated that soils could be separated at the soil-type level based on reflectance intensity in various absorption bands. Soil collected in the 40- to 60-cm depth appeared to have higher reflectance intensities than those from the 0- to 20-cm depth. Removal of OM from soil samples promoted higher reflectance intensity in the entire spectrum. Amorphous and crystalline Fe influenced reflectance differently. Weathering of basaltic soils was correlated with alterations in the reflectance intensities and absorption features of the spectral curves. Multivariate analysis demonstrated that this technique was efficient in the estimation of clay, silt, kaolinite, crystalline Fe, amorphous Fe, and Mg through the use of reflected energy of the soils.

Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Relationships between soil nitrate desorbed from anion-exchange membranes, canopy reflectance and nitrate leaching losses from cool-season lawn turf.

Nutrient leaching studies are expensive and require expertise in water collection and analyses. Less expensive or easier methods that estimate leaching losses would be desirable. The objective of this study was to determine if anion-exchange membranes (AEMs) and reflectance meters could predict nitrate (NO3-N) leaching losses from a cool-season lawn turf. A two-year field study used an established 90% Kentucky bluegrass (Poa pratensis L.)-10% creeping red fescue (Festuca rubra L.) turf that received 0 to 98 kg N ha-1 month-1, from May through November. Soil monolith lysimeters collected leachate that was analyzed for NO3-N concentration. Soil NO3-N was estimated with AEMs. Spectral reflectance measurements of the turf were obtained with chlorophyll and chroma meters. No significant (p > 0.05) increase in percolate flow-weighted NO3-N concentration (FWC) or mass loss occurred when AEM desorbed soil NO3-N was below 0.84 µg cm-2 d-1. A linear increase in FWC and mass loss (p < 0.0001) occurred, however, when AEM soil NO3-N was above this value. The maximum contaminant level (MCL) for drinking water (10 mg L-1 NO3-N) was reached with an AEM soil NO3-N value of 1.6 µg cm-2 d-1. Maximum meter readings were obtained when AEM soil NO3 N reached or exceeded 2.3 µg cm-2 d-1. As chlorophyll index and hue angle (greenness) increased, there was an increased probability of exceeding the NO3-N MCL. These data suggest that AEMs and reflectance meters can serve as tools to predict NO3-N leaching losses from cool-season lawn turf, and to provide objective guides for N fertilization.

Near infrared spectroscopy for soil bulk density assessment

P>Soil bulk density values are needed to convert organic carbon content to mass of organic carbon per unit area. However, field sampling and measurement of soil bulk density are labour-intensive, costly and tedious. Near-infrared reflectance spectroscopy (NIRS) is a physically non-destructive, rapid, reproducible and low-cost method that characterizes materials according to their reflectance in the near-infrared spectral region. The aim of this paper was to investigate the ability of NIRS to predict soil bulk density and to compare its performance with published pedotransfer functions. The study was carried out on a dataset of 1184 soil samples originating from a reforestation area in the Brazilian Amazon basin, and conventional soil bulk density values were obtained with metallic ""core cylinders"". The results indicate that the modified partial least squares regression used on spectral data is an alternative method for soil bulk density predictions to the published pedotransfer functions tested in this study. The NIRS method presented the closest-to-zero accuracy error (-0.002 g cm-3) and the lowest prediction error (0.13 g cm-3) and the coefficient of variation of the validation sets ranged from 8.1 to 8.9% of the mean reference values. Nevertheless, further research is required to assess the limits and specificities of the NIRS method, but it may have advantages for soil bulk density predictions, especially in environments such as the Amazon forest.

Spatial Modeling of a Soil Fertility Index using Visible-Near-Infrared Spectra and Terrain Attributes

Our objective was to develop a methodology to predict soil fertility using visible near-infrared (vis-NIR) diffuse reflectance spectra and terrain attributes derived from a digital elevation model (DEM). Specifically, our aims were to: (i) assemble a minimum data set to develop a soil fertility index for sugarcane (Sarcharum officinarum L.) (SFI-SC) for biofuel production in tropical soils; (ii) construct a model to predict the SFI-SC using soil vis-NIR spectra and terrain attributes; and (iii) produce a soil fertility map for our study area and assess it by comparing it with a green vegetation index (GVI). The study area was 185 ha located in sao Paulo State, Brazil. In total, 184 soil samples were collected and analyzed for a range of soil chemical and physical properties. Their vis-NIR spectra were collected from 400 to 2500 nm. The Shuttle Radar Topographic Mission 3-arcsec (90-m resolution) DEM of the area was used to derive 17 terrain attributes. A minimum data set of soil properties was selected to develop the SFI-SC. The SFI-SC consisted of three classes: Class 1, the highly fertile soils; Class 2, the fertile soils; and Class 3, the least fertile soils. It was derived heuristically with conditionals and using expert knowledge. The index was modeled with the spectra and terrain data using cross-validated decision trees. The cross-validation of the model correctly predicted Class 1 in 75% of cases, Class 2 in 61%, and Class 3 in 65%. A fertility map was derived for the study area and compared with a map of the GVI. Our approach offers a methodology that incorporates expert knowledge to derive the SFI-SC and uses a versatile spectro-spatial methodology that may be implemented for rapid and accurate determination of soil fertility and better exploration of areas suitable for production.

The Reflectance Signature of Canopy Components: Implications for the Interpretation of Remotely Sensed Images

The reflectance signatures of plantation pine canopy and understorey components were measured using a spectro-radiometer. The aim was to establish whether differences observed in the reflectance signature of stressed and unstressed pine needles were consistent with observed differences in the reflectance of multispectral Landsat Thematic Mapper (TM) images of healthy and stressed forest. Because overall scene reflectance includes the contribution of each scene component, needle reflectance may not be representative of canopy reflectance. In this investigation, a limited dataset of reflectance signatures from stressed and unstressed needles confirmed the negative relationship between pine needle health and reflectance which was observed in visible red wavelengths. However, the reflectance contribution from bushes, pine needle litter and bare soil tended to reinforce this relationship suggesting that in this instance, overall scene reflectance is comprised of the proportional reflectance of each scene component. In near infrared wavelengths, differences between healthy and stressed needle reflectance suggested a strong positive relationship between reflectance and tree health. For Landsat TM images, previous research had only observed a weak positive relationship between stand health and near infrared reflectance in these pine canopies. This suggests that for multispectral Landsat TM images, reflectance of near infrared light from pine canopies may be affected by other factors which may include the scattering of light within canopies. These results are seen as promising for the use of hyperspectral images to detect stand health, provided that pixel reflectance is not influenced by other scene components.

Chemical characterization and infrared spectroscopy of soil organic matter from two southern brazilian soils

Soil organic matter from the surface horizon of two Brazilian soils (a Latosol and a Chernosol), in bulk samples (in situ SOM) and in HF-treated samples (SOM), was characterized by elemental analyses, diffuse reflectance (DRIFT) and transmission Fourier transform infrared spectroscopy (T-FTIR). Humic acids (HA), fulvic acids (FA) and humin (HU) isolated from the SOM were characterized additionally by ultraviolet-visible spectroscopy (UV-VIS). After sample oxidation and alkaline treatment, the DRIFT technique proved to be more informative for the detection of "in situ SOM" and of residual organic matter than T-FTIR. The higher hydrophobicity index (HI) and H/C ratio obtained in the Chernosol samples indicate a stronger aliphatic character of the organic matter in this soil than the Latosol. In the latter, a pronounced HI decrease was observed after the removal of humic substances (HS). The weaker aliphatic character, the higher O/C ratio, and the T-FTIR spectrum obtained for the HU fraction in the Latosol suggest the occurrence of surface coordination of carboxylate ions. The Chernosol HU fraction was also oxygenated to a relatively high extent, but presented a stronger hydrophobic character in comparison with the Latosol HU. These differences in the chemical and functional group composition suggest a higher organic matter protection in the Latosol. After the HF treatment, decreases in the FA proportion and the A350/A550 ratio were observed. A possible loss of FA and condensation of organic molecules due to the highly acid medium should not be neglected.