Clay minerals are the most reactive and important inorganic components in soils, but soil mineralogy classifies as a minor topic in soil sciences. Revisiting soil mineralogy has been gradually required. Clay minerals in soils are more complex and less well crystallized than those in sedimentary rocks, and thus, they display more complicated X-ray diffraction (XRD) patterns. Traditional characterization methods such as XRD are usually expensive and time-consuming, and they are therefore inappropriate for large datasets, whereas visible and near-infrared reflectance spectroscopy (VNIR) is a quick, cost-efficient, and nondestructive technique for analyzing soil mineralogic properties of large datasets. The main objectives of this review are to bring readers up to date with information and understanding of VNIR as it relates to soil mineralogy and attracts more attention from a wide variety of readers to revisit soil mineralogy. We begin our review with a description of fundamentals of VNIR. We then review common methods to process soil VNIR spectra and summary spectral features of soil minerals with particular attention to those <2 μm fractions. We further critically review applications of chemometric methods and related model building in spectroscopic soil mineral studies. We then compare spectral measurement with multivariate calibration methods, and we suggest that they both produce excellent results depending on the situation. Finally, we suggest a few avenues of future research, including the development of theoretical calibrations of VNIR more suitable for various soil samples worldwide, better elucidation of clay mineral-soil organic carbon (SOC) interactions, and building the concept of integrated soil mapping through combined information (e.g., mineral composition, soil organic matter-SOM, SOC, pH, and moisture).
Special Funding for Soil MineralogyCUG170106National Natural Science Foundation of China4177203241472041NSFC for Young Scholars4140203641602037NSF of Hubei for Young Scholars2016CFB183China Postdoctoral Science Foundation2015M582301China Scholarship Council2017064100172017064100061. Introduction
Soils are open, complex, and dynamic systems as well as fundamental natural environments for animals, plants, microorganisms, and human interaction [1]. Mineral composition is the most fundamental property of a soil, and soil minerals account generally for half the soil volume [2]. According to Churchman [3], clay minerals in the soil context are “secondary inorganic compounds of <2 μm size” including Fe, Al, and Mn oxides (hydroxides and oxyhydroxides), as well as noncrystalline phases. Importantly, they are the most reactive and important inorganic components in soils, and they occur commonly in close association with the most reactive organic matter [4, 5]. Clays influence soil function through both their bulk properties and their associations with their huge outer/inner surfaces (e.g., cation exchange capacity [6]). The effort involved in comprehensive understanding of the nature of soil minerals is of particular importance as they may help us explain and predict how different soil types function [7].
However, soil mineralogy (mainly clay mineralogy) is still a minor topic in soil sciences. This may be due partly to the unjustified assumption that a given soil mineral will have the same properties as those of its better-crystallized counterpart that formed in a more “geologic” context (e.g., sedimentary kaolinite will have the same properties as pedogenic kaolinite) [4]. Revisiting soil mineralogy has been gradually important, for instance, in terms of the manner by which soil minerals are defined and investigated [8].
The most commonly used method to characterize soil minerals is XRD, which is fundamentally qualitative. Since soil clay minerals are generally more complex and less well crystallized than those of geological environments [9–11], they display more complicated XRD patterns [12, 13]. Despite quantitative improvements of XRD [14], mineral characterization is usually expensive and time-consuming [2]. Some chemical extraction procedures can be useful in the analysis of Fe oxides. However, this is expensive, time-consuming, and can complicate our scientific interpretation of the soil by changing the chemical equilibrium between soil solution and solid phases in soil specimens [15, 16]. Thus, these conventional analyses are not appropriate for larger scale soil studies, and we must use an alternative method to target and characterize soil minerals.
Visible and near-infrared reflectance spectroscopy (VNIR, 350–2500 nm), that is, the study of light of the visible and near-infrared reflected from material surfaces, is a quick, cost-efficient, and nondestructive technique in soil sciences [17, 18]. This technique has been greatly developed in soil sciences in the past several decades and has seen apparent exponential growth over the past 20 years [19]. VNIR has been of increasing interest for the analyses of soil parameters including soil organic carbon, pH, bulk texture, elemental concentration, and cation exchange capacity [20, 21]. In soil mineralogy, VNIR can be used to characterize various soil mineralogic properties such as clay mineral composition, clay content, and mineral weathering/alteration degree, although quartz and feldspar have weak/nonexistent absorption in the VNIR range [22–24]. In this paper, we aim to bring readers up to date with VNIR as it relates to soil mineralogy and we seek to attract more attention from readers to revisit soil mineralogy.
2. Fundamentals of VNIR
The VNIR part of the electromagnetic spectrum includes both the visible (350–780 nm) and near-infrared (780–2500 nm) ranges, which overlaps with the optical radiation range (100–1000 nm; Figure 1). Sometimes, the 350–1000 wavelength range is referred as VNIR (visible-near-infrared), and the 1000–2500 range is referred as the SWIR (short-wave infrared) in remote sensing literature [25]. The human eyes and brain can process spectral information from the visible region and see color, while modern spectroscopy can observe precise details over a much broader wavelength range.
Optical spectrum between 102 to 106 nm.
2.1. Absorption, Scattering, and Emission
When photons enter a solid, liquid, or gaseous material, they will either be absorbed, reflected from its surface, or pass through it [26]. The reflective process is defined as scattering, and the scattered photons can be detected and measured. Photons can also be detected when they are emitted from a surface with a temperature above absolute zero [25]. Three general physical processes (i.e., electronic transitions, vibrational transitions, and rotational transitions) result in the absorption bands in the spectra of materials. The absorption bands in the VNIR range are derived from both the electronic and vibrational transitions [27, 28].
2.2. Causes of Absorptions in the VNIR Region2.2.1. Electronic Transitions
Discrete ions and atoms have independent energy states. A photon is emitted from an atom when one of its electrons moves to a lower energy state. When an atom absorbs a photon of a given wavelength, its electrons move from a relatively low electron state to a higher one [25]. These electron processes occur because of their high energy and mobility. The electronic processes are mainly caused by (1) crystal-field effects. Since iron is a very common transition element in minerals, a common electronic process revealed in the visible region is due to unfilled d-orbitals of Fe-oxide minerals [24, 29]. Electron energy levels are influenced by many factors, including the valence state of the atom (e.g., Fe2+ and Fe3+), the type of ligands, the asymmetry of the location it occupies, the distance between the metal ion and the ligand, and the deformation degree of the site [28]. (2) Charge transfer: it is dominated by mineralogy, and it is hundred times more powerful than the crystal-field effects. It is the main reason of the red color of hydroxides and Fe oxides. Moreover, the conduction bands and color centers can also be causes of the electronic transitions in some minerals [25].
2.2.2. Vibrational Transitions
The bonds in a crystal lattice or molecule vibrate like springs. The molecule’s mass and the strength of each molecular bond dominate their vibration frequency [25]. The absorption bands in the VNIR range are observed as a consequence of molecular vibrations [30]. Soil minerals (e.g., phyllosilicate and carbonate minerals), in particular, have unique absorption features in the VNIR region due to overtones and vibrational combinations related to the stretching and bending of the molecular bonds such as O-H, C-H, C-C, and N-H [31].
The raw spectra are usually preprocessed through various approaches to accentuate features and remove signal noise [32]. The processed soil spectra facilitate mineral identification, and the accuracy of soil mineral prediction is greatly improved through the use of various preprocessing methods [33]. The following preprocessing methods for spectra have been used in previous soil mineralogic studies.
3.1.1. Continuum Removal Approaches
The continuum removal approach aims to remove background noise and isolate particular absorption features for identification and analysis [34]. The continuum is usually determined using local maxima to generate a hull of boundary points (Figure 2(a)) [22]. All the boundary points are fitted by straight-line segments, and then, the continuum removal is calculated by removing the original reflectance intensities from corresponding intensities of the continuum (Figure 2(b)) [23]. Continuum removal analysis is a particularly robust tool for detecting and predicting iron oxides and phyllosilicate minerals. Thus, it is feasible to substitute a statistical method to apply to soil mineralogy studies [10, 20, 22, 24].
(a) Continuum removal of the spectrum corresponding to a red earth soil sample [24]; (b) a continuum-removed spectrum and some spectral parameters (i.e., W, D, and F).
Absorption bands in the VNIR region can be described by geometrical parameters derived from the continuum removal curve (Figure 2(b)). Four parameters are directly displayed in Figure 2(b), which include position (P), width (W), depth (D), and full width at half maximum (FWHM, abbreviated to “F”). The parameter asymmetry (AS) is calculated as follows:(1)AS=FleftFright,where Fleft represents the left width at half maximum, and Fright represents the right width at half maximum [20].
3.1.2. Smoothing Techniques
Smoothing techniques are used to extract the maximum amount of information from each spectrum possibly by minimizing the influence of background noise [32]. Commonly used smoothing techniques include the Savitzky–Golay transform (SG [35]), Norris smoothing filter (NG [36]), and averaging spectra [37]. SG smoothing eliminates the influences of ground interference noise and baseline float, thus enhancing the signal-to-noise ratio. NG smoothing removes the effects of particle-size variation when the soil samples vary in texture, moisture, and grain size [32].
3.1.3. Derivative Algorithms
Derivative algorithms can rapidly identify characteristic positions of spectral minimum, maximum, and inflection point values [32]. Additionally, the effect of variation in optical setup and sample grinding is eliminated after derivative transformation [38]. Because the spectral noise tends to amplify with derivative transform, a smoothing technique is often used before the derivative algorithm [37]. The spectral curve after the first derivative, for example, is better at discriminating goethite and hematite and estimating their abundance, with two peaks at 435 and 535 nm for goethite and a single absorption at ∼570 nm for hematite (Figure 3) [39].
Raw spectra for goethite (a) and hematite (b) and the first-derivative spectra of goethite (c) and hematite (d) (modified from Zheng et al. [39]).
3.2. Spectral Features of Soil Minerals3.2.1. Fe-Oxide Minerals
Fe-oxide minerals are known to be pedogenic indicators for investigating soil temperature and moisture regimes, which are directly related to pedogenic climate evolution [24, 40]. Fe-oxide minerals are the main active components in the VNIR region (350–1000 nm) since most electron transitions are caused by various kinds of iron oxides [41, 42]. The most common Fe-oxide minerals in soils are goethite (α-FeOOH) and hematite (α-Fe2O3), which can track climate change [43, 44]. Goethite and hematite exhibit diagnostic spectral features in the VNIR region, and the absorption bands are generally broad and smooth (Figure 3). A strong absorption band near 920 nm indicates the presence of goethite (Figure 3(a)), and four absorption bands at 420, 480, 600, and 1700 nm can be used to map its distribution [39]. Hematite is dominated by three absorption bands at 520, 650, and 880 nm [45]. Both goethite and hematite have an absorption band at around 500 nm (480 for goethite and 520 for hematite, resp.); the band for goethite (at 480 nm) is narrow with intense reflectance, while the band for hematite (at 520 nm) is wide with low reflectance (Figure 3(a)). The absorptions in the VNIR region cause the vivid colors of Fe oxides, for example, yellow goethite and red hematite [37]. For a spectral curve representing a sample soil mixture, the width of the absorption band at ∼870 nm (W870) is higher when the soil sample contains more Fe oxides [46]. The concave shape of the 800–1000 nm range indicates the crystallinity of the Fe-oxide minerals. When a soil sample is composed of well-crystallized minerals, the corresponding spectrum reveals a symmetric and deeper feature in this range [47].
3.2.2. Clay Minerals
Clay minerals are frequently used as climatic indicators since their nature is directly influenced by the temperature and amount of precipitation at the site during pedogenesis [9, 48]. As climate conditions shift from cool/dry to warm/moist, the dominant clay minerals go from chlorite/illite → vermiculite → montmorillonite → kaolinite [24, 49]. The dominant clays in soils show diagnostic absorptions in the SWIR domain [39]. These absorption bands are caused by vibrational transitions and commonly display sharp and narrow features (Figure 4). The diagnostic bands are mainly focused on ∼1400 nm (overtones caused by OH), ∼1900 nm (overtones caused by molecular water), and ∼2200 nm (combination tones caused by Al-OH [50, 51]). Additionally, some weak absorption bands in the 2300–2500 nm region are related to the presence of Fe-OH and/or Mg-OH in the clay minerals [24].
Raw and continuum-removed spectra of five common clay minerals.
The spectral characteristics of some clay minerals are showed in Figure 4 and Table 1. Chlorites are a group of clay minerals containing specific octahedral cations such as Fe, Mg, and Al [52]. Their reflectance spectra exhibit a weak absorption band at approximately 1400 nm and triple absorption features near 2300 nm. The bands at 2250 and 2350 nm are related to Fe-OH and Mg-OH, respectively [53]. Illite is characterized by three prominent absorptions at ∼1400, ∼1900, and ∼2200 nm. Two secondary diagnostic Al-OH absorption peaks close to 2344 and 2445 nm are modified by Fe and Mg tschermak cation exchange [24, 31]. Vermiculite has two broad absorptions at 1400 and 1900 nm and two weak absorptions near 2200 and 2300 nm [39]. Montmorillonite has three strong and sharp absorption bands at ∼1400, ∼1900, and ∼2200 nm, which are similar to but generally stronger than illite.
Spectral features of some predominant soil minerals.
Absorption band
Origin
Diagnostic feature
Description
Assigned soil minerals
∼465 nm
Electronic transition
6A1g → 4T1g,4Eg
Width
Narrow
Goethite
Wide
Hematite
∼650 nm
6A1g → 4T2g
Position
Left → right
Goethite → hematite
∼900 nm
6A1g → 4T1g
Position
Left → right
Hematite → goethite
∼1400 nm
Molecular vibration
OH
Double absorption band (1395 nm and 1415 nm)
Kaolinite
Asymmetry
>1
Kaolinite
<1
Montmorillonite
∼1900 nm
H2O
Position
<1910 nm
Montmorillonite/vermiculite
>1910 nm
Kaolinite/illite
Depth
Intense → shallow
Montmorillonite → illite → kaolinite
Width
Narrow
Montmorillonite
Wide
Vermiculite/illite
∼2200 nm
Al-OH
Double absorption band (2163 nm and 2208 nm)
Kaolinite
Depth
Intense
Kaolinite
Shallow
Vermiculite/montmorillonite/illite
Asymmetry
>1
Kaolinite
<1
Montmorillonite
>2300 nm
Fe-OH and Mg-OH
Additional absorptions (2340 nm and 2445 nm)
Illite
Triple absorption features between 2330 and 2360 nm
Chlorite
Additionally, the combination bands produced by the vibrations of absorbed water cause weak shoulders near 1468 nm and 1970 nm for montmorillonite spectra [37]. Kaolinite is featured by two spectral doublets: one is near 1400 nm (1390 and 1410 nm), and the other is near 2200 nm (2160 and 2210 nm).
3.2.3. Carbonates
In soils, carbonates are leached from the surface with time and accumulate in the subsoil at a certain depth [54]. The presence of carbonate is widely used as a basic soil characteristic to describe soil types and quantify soil erosion [22]. Carbonates are characterized by several absorptions in the VNIR domain, caused by overtones and combinations of fundamental vibrations of the CO32− ion (Figure 5) [31, 37]. A strong absorption band at ∼2350 nm and three weak absorption bands at ∼1900, ∼2000, and ∼2160 nm were reported by Hunt and Salisbury [55] for carbonates in the NIR region, with the ∼2350 nm absorption showing obvious double-band structures (Figure 5).
Absorptions for calcite in the VNIR domain, caused by overtones and combinations of fundamental vibrations of the CO32− ion.
3.3. Prediction from the Continuum Removal Spectra
As discussed in Section 3.1.1, several geometrical features of the absorption bands can be extracted through the continuum removal method. Those parameters (e.g., P, D, and AS) from the continuum removal spectra are key to characterizing and predicting mineral compositions in soils. Viscarra Rossel et al. [23] quantitatively estimated the mineral composition by using the continuum removal method. Compositions of soil minerals such as kaolinite, illite, Al-smectite, goethite, and hematite are considered in this study, and the parameter D is selected for prediction. The spectroscopic predictions are generally in consistence with those interpreted by XRD analysis. According to Dufrechou et al. [20], the parameter D at ∼1400, ∼1900, and ∼2200 nm was strongly affected by the amounts of kaolinite, illite, and montmorillonite in soil mixtures. Additionally, the estimation of montmorillonite abundance shows reliability when compared with XRD results. Five parameters (P, D, W, F, and AS) were used in the work by Zhao et al. [24] for assessing the utility of the continuum removal method. We compared these parameters with the results from both XRD and DRS analyses and found that some of the parameters are good at mineral content prediction. Furthermore, some parameters (e.g., AS at ∼2200 nm) are confirmed as reliable proxies for soil weathering and paleoclimate reconstruction.
4. Chemometric Methods
VNIR spectra of soil mixtures are commonly weak and nonspecific due to (1) low concentration of particular soil minerals, (2) scatter effects caused by soil structure, (3) overlapping absorptions of soil attributes, and (4) influences of specific constituents such as quartz [37]. All of these factors pose a challenge for VNIR analyses. Therefore, useful information needs to be mathematically extracted from the spectra and correlated with soil attributes [45]. The development of VNIR in soil studies would have been impossible without the parallel application of chemometric methods [56].
Building a predictive soil mineral abundance model (i.e., multivariate calibration) is an important first step in chemometric analysis. Overall, we should understand the data and the objective of the modeling prior to building a model. Then, the spectral dataset is preprocessed and subdivided. Finally, we can proceed to build, evaluate, and select models [57].
4.1. Prior to Model Building
The first step in any model building process for the study of spectral pedology is to understand the characteristics of the dataset. We need to consider three main concepts in understanding the dataset process [57]: (1) understanding the distribution of the responses (i.e., outcomes): the responses are either numerical or categorical. In the model building process for soil mineral analysis, the outcomes (e.g., contents of clay/Fe-oxide minerals) are described numerically. Understanding the characteristics of responses provides better ways for partitioning the data into calibration and validation sets; (2) understanding the nature of the predictors: the predictors in the spectral dataset are numerical, since they are usually the spectral signals between 350 and 2500 nm. In fact, these predictors are highly related, leading to numerically redundant information. Different predictors are suitable for different kinds of models. For example, partial least squares can be used for correlated predictors, while recursive partitioning can manage missing predictor information [58]; (3) the relationship between the amount of the predictor set (P) and the sample set (N): when building a model for a soil mineral study, the dataset commonly has far fewer samples (N<200) than predictors (P>2000). Therefore, a model that can handle dataset where N<P is preferred.
After understanding the dataset, a preprocessing procedure is often used for improving the performance of the model [47, 59]. For a model used in a soil mineral study, the data transformations for multiple predictors contain the following methods: (1) data reduction: principal component analysis (PCA [60]) is a commonly used data reduction technique. In this technique, the number of datasets is largely reduced by seeking principle components (PC)—linear combinations of the predictors that capture the greatest possible variance; (2) removing predictors: in some cases, removing predictors prior to modeling has potential advantages. For example, Adeline et al. [59] showed that performances of the predictive models were globally stable and accurate when the spectral resolution decreased from 3 nm to 60 nm. Additionally, for a model based on a spectral signal dataset, the spectra were transformed to apparent absorbance: A=log1/R prior to developing a regression model, and the spectral preprocessing methods discussed in Section 3 also have potential for model performance improvement [18, 23].
4.2. Candidate Models
Once we fully understand the dataset, the next step is to setup several candidate models. The most commonly used type of model in soil mineral analysis is a regression model, which is defined as a model that predicts numerical outcomes [57]. Establishing a regression model related to the soil VNIR spectral data is the basic role of chemometric analysis [61]. The regression models are subdivided into linear and nonlinear regressions. Linear regressions are the dominant calibration methods for spectral pedology and include partial least squares regression (PLSR [62]) and principal component regression (PCR [63]). The nonlinear data are managed by data mining techniques, namely, multivariate adaptive regression splines (MARS [38]), neural networks (NN [64]), and regression tree analysis (RTA [65]).
4.2.1. Linear Regression Models
Both PLSR and PCR can deal with predictors that are highly collinear and are effective in situations where the number of predictors is far beyond the number of available samples [37]. Furthermore, PLSR and PCR are closely related and share similar prediction errors in most situations [61]. Regardless, the PLSR algorithm is usually preferred in spectral pedology analysis because (1) it maximizes covariance between response variables and predictors so that the model is more interpretable, and (2) it is a faster algorithm [45].
PLSR has been widely and successfully used in predicting the mineralogic compositions of weathering levels of soils. Viscarra Rossel et al. [2] accurately predicted the concentrations of kaolinite, illite, and smectite (R2=0.94,0.96, and 0.92, resp.) in mineral mixtures, although the prediction for Fe oxides was biased against measurement. Summers et al. [66] and Ostovari et al. [67] showed that the PLSR method is good at predicting CaCO3 content with R2 values of 0.69 and 0.71 for soil samples from Australia and Iran, respectively. The total clay content and free iron in soils were also proven to be predictable attributes by the PLSR model [68, 69].
4.2.2. Nonlinear Regression Models
The use of models that are inherently nonlinear in nature (i.e., data mining techniques) has gained increasing attention in recent years [37, 61]. A more detailed description of the nonlinear models is available in Kuhn and Johnson [57]. Previous studies have suggested that nonlinear regression models or the combination of nonlinear and linear models may provide better predictions for soil properties. Mouazen et al. [70] showed that a combined PLSR-NN model was better at predicting soil properties than a PLSR model. Viscarra Rossel and Behrens [45] proposed that the combined FSVIP-ANN and FSMARS-ANN models were the best models for predicting clay content, pH, and soil organic carbon (SOC) when both the parsimony and accuracy of the model were taken into consideration. Mulder et al. [71] determined the mineral composition of a soil by coupling an RTA model with exponential Gaussian optimization results. The abundances of kaolinite and calcite were predicted with acceptable RMSE values (<0.1) in both laboratory and field samples.
4.3. Model Evaluation
Two techniques are commonly used to test the prediction performance of the model [33, 56]. In the first, the soil spectral configuration database is randomly divided into a calibration dataset and a validation dataset [72]. The calibration dataset (usually ∼2/3 of the complete database) is used to derive the model, while the validation dataset (commonly contains 1/3 of the complete database) is set aside to exclusively validate the derived model. This process is used to obtain realistic estimates of prediction accuracy. The second is a procedure called cross validation. It uses the “leave-group-out” method (namely, repeated random subsampling validation method [73]) and was adopted to verify the predictive capability for the calibration dataset. A calibration dataset containing X samples is built from the total database N (N≥X+1). The soil property value of the other N−X samples for validation is predicted. The prediction of relative soil mineral abundance is obtained by repeating the cross-validation process [74].
ParLeS version 3.0 is usually used for multivariate calibration performance [75]. The bias and accuracy of the prediction models are assessed by adjusting the coefficient of determination (R2) between observed and predicted values, the mean error (ME)(2)ME=∑i=1NYi′−Yi,and the root mean-square error (RMSE)(3)RMSE=1N∑i=0nYi′−Yi2,where N is the number of the dataset, Yi is the observed value, and Yi′ represents the predicted value [18, 23].
We compromise between model parsimony and model accuracy to find the most satisfactory model [76]. The Akaike information criterion (AIC) is suggested for best-performing algorithm selection [45]:(4)AIC=2p+n×lnRMSE,where p is the number of factors, and n is the number of samples used in the prediction. The best model will have the minimum AIC value.
4.4. Feature Selection
Feature selection is mainly applied to remove redundant and/or noninformative predictors from the model [57] and may improve model accuracy. Some models such as PLSR, MARS, and RT will provide a feature selection procedure by default.
The variable importance of the projection (VIP) and b-coefficient scores obtained by the PLSR model help us measure the statistical significance of predictors and select the most important ones [77]. The VIP score of the kth predictor is calculated as follows:(5)VIPk=K⋅∑a=1aoptwak2⋅Ra2,where K represents the total number of the predictor variables, aopt is an optimal number of latent variables selected by the PLSR model, wak is the loading weight for the kth latent variable, and Ra2 represents the adjust coefficient of determination of the ath latent variable in the PLSR model [78]. A predictor (such as wavelength) is selected and considered to be very important if (1) the VIP exceeds the threshold value of one (Chong and Jun [77]) and (2) its b-coefficient is higher than the b-coefficient based on all spectral bands [23].
According to Gomez et al. [22], the important spectral bands selected by the PLSR model are related to the presence of clay minerals such as kaolinite and illite. Additionally, surrogate spectral features selected by VIP and b-coefficient approach contain enough information to satisfactorily estimate the studied soil attributes. According to Viscarra Rossel and Behrens [45], a combined NN and feature selection model (FSvip-ANN) is the best method to predict clay content and pH and produce smaller RSME and AIC values.
5. Comparison between Spectral Measurement and Multivariate Calibration
The relative abundance of minerals in a soil sample can be predicted either by spectral analyses (e.g., continuum removal) or chemometric methods (e.g., PLSR and NN) [20, 24, 45, 47]. Although both types of methods correlate the spectral signal with information about the soil minerals, they differ in many ways, including their focused spectral bands, complexity, and how they are applied.
5.1. Focused Spectral Bands
Spectral analyses focus on specific absorption bands representative of the corresponding soil minerals, while the multivariable regression algorithms commonly use the signals from the whole 350–2500 nm region. In some cases, the 350–400 nm and 2450–2500 nm ranges with low instrumental signal-to-noise ratios are removed [59, 79]. Therefore, a multivariable regression model deals with over 1000 spectral bands—many more than the number of focused bands in a continuum removal study. Moreover, several geometric parameters can be extracted from each band in a spectral measurement, including P, W, D, F, and AS, whereas only the information of depth for each band can be gleaned from a chemometric study. Note that some algorithms intrinsically provide a feature selection method (e.g., SMLR and PLSR), and it has been shown that the most important features selected by a regression model are the ones that we should pay the most attention to in a spectral measurement study [22].
5.2. Complexity
Theoretically, multivariate calibration is very complicated because it involves a larger number of algorithms and because different algorithms have the potential to be combined into better predictive models, depending on the situation [45, 70, 79]. However, in practice, multivariate modeling and prediction is not that complicated. Thanks to the development of executable and fast running software such as ParLeS and Unscrambler [75, 80], the difficult calculation process can be done much more easily. On the other hand, spectral measurement studies cost more time because we must (1) identify a soil mineral based on the spectral features, (2) extract parameters from the bands, and (3) relate those parameters with the information about the soil mineral.
5.3. Application Preference
The geometric features of the spectra are more suitable for monitoring the molecular structural changes of soil minerals, since the variations of the absorption bands are caused by electron transitions (e.g., Fe2+ to Fe3+) and molecular vibration (e.g., Al-OH versus H2O). Thus, spectral measurement is widely and successfully applied to (1) measure mineral physicochemistry that is sensitive to changes in metamorphic grade [53, 81], (2) map and monitor mineral erosion, deposition, and weathering of minerals [24, 51], and (3) explore water and potential life on extraterrestrial objects [10, 21]. Chemometric methods are more often used in monitoring overall soil properties, since almost all of the signals in the VNIR domain are involved in the modeling process. Several soil attributes are successfully determined by an appropriate multivariate calibration technique, including soil clay [23, 69], organic matter [32, 67], and nitrogen content [82, 83].
Table 2 is a review of some soil mineralogic attributes predicted by VNIR spectroscopy using either chemometric analysis or spectral-based measurement. In this summary, most of the studies used soil samples for analysis, and many of them are among diverse soil types (Table 2). The predictions of the soil properties are still good when there is great range of soil types (e.g., 22, 45, 84, and 85). A single mineral (e.g., kaolinite and goethite) is more precisely predicted when mineral mixtures are used in the measurements [2, 86]. The studies in Table 2 include both data collected in the laboratory and data based on field soil sensing. In the lab, the sample pretreatment and illumination conditions can be controlled to eliminate the influences of the moisture and the grain size of the soil sample [18]. While in the field, the VNIR spectroscopy may be affected by many potential problems such as variable distances between the sensor and the soil, the smearing of soil surfaces, the size of the soil aggregates, and the amount of moisture [87]. These potential problems may reduce the prediction accuracy of field-based analysis [22, 85]. However, the field-based VNIR spectroscopy is more attractive because it (1) enables the potential analysis of soil properties with promising results in previous studies [87] and (2) reduces the cost of the measurement by simplifying the sample preparation. Based on the results of the studies, PLSR is proved to be the most robust soil mineralogic analysis method amongst all of the multivariate calibrations (Table 2; 0.43<R2<0.96). The CR-based model is good at predicting clay mineral concentration (e.g., 20, 47, and 88; R2>0.79). In some cases, the nonlinear models (NN and MARS) exhibit better estimation in predicting soil mineralogy than the PLSR model (e.g., 45 and 89). In general, when a soil mineral is investigated by spectroscopy, the PLSR and the CR-based models are the most promising methods to provide estimates of mineral abundance.
A review of VNIR spectroscopy for the prediction of soil attributes using feature-based and chemometric methods in previous studies.
Soil mineralogic attributea
Soil typeb
Location
Samplec
ncal/nvald
Methode
R2f
RMSEg
Source
ClayT
Alfisols and entisols
Romania
Lab/sieved (2)-air-dried
210/90
SVR
0.27
84.4
[90]
ClayT
Not mentioned
France
Lab/sieved (2)-air-dried
99/49
PLSR
0.76
33.6
[59]
ClayT
Diverse
Romania
Lab/sieved (2)-air-dried
210/90
SVR
0.46
85.3
[91]
ClayT
Vertisols and alfisols
India
Lab/sieved (2)-air-dried
175/58
PLSR
0.80
2.2
[92]
ClayT
Diverse
Denmark
Lab/sieved (2)-air-dried
480CV
PLSR
0.94
17
[84]
ClayT
Entisol and inceptisols
Italy
Lab/sieved (2)-oven-dried
175/60
PLSR
0.80
0.17
[93]
ClayT
Diverse
Australia
Lab/sieved (2)-air-dried
1104CV
DWT-ANN
0.88
64.2
[45]
ClayT
Diverse
Australia
Lab/sieved (2)-air-dried
1104CV
PLSR
0.83
77.7
[45]
ClayT
Entisol
Turkey
Lab/sieved (2)-air-dried
502CV
PLSR
0.84
38.20
[89]
ClayT
Entisol
Turkey
Lab/sieved (2)-air-dried
502CV
MARS
0.86
35.10
[89]
ClayT
Diverse
France
Lab/sieved (2)-oven-dried
52CV
PLSR
0.85
31.2
[22]
ClayT
Diverse
France
Lab/sieved (2)-oven-dried
52CV
CR
0.73
44
[22]
ClayT
Diverse
France
Field/fresh
52CV
PLSR
0.64
49.6
[22]
ClayT
Diverse
France
Field/fresh
52CV
CR
0.58
82
[22]
ClayT
Diverse
USA
Lab/sieved (2)-air-dried
72CV
PLSR
0.92
41
[85]
ClayT
Diverse
USA
Field/fresh
72CV
PLSR
0.83
61
[85]
ClayK
Diverse
Australia
Lab/sieved (2)-air-dried
102CV
PLSR
0.95
111.5
[94]
ClayK
Not mentioned
China
Lab/sieved (2)-air-dried
20CV
CR-MLR
0.80
10.41
[88]
ClayK
Diverse topsoil
Australia
Lab/sieved (2)-air-dried
4606CV
Model trees
0.52
0.8
[95]
ClayK
Diverse subsoil
Australia
Lab/sieved (2)-air-dried
2492CV
Model trees
0.46
1.1
[95]
ClayK
Mineral-organic mixes
Australia
Lab/sieved (0.2)-oven-dried
8CV
PLSR
0.94
0.36
[2]
ClayI
Diverse
Australia
Lab/sieved (2)-air-dried
90CV
PLSR
0.96
102.1
[94]
ClayI
Not mentioned
China
Lab/sieved (2)-air-dried
20CV
CR-MLR
0.79
28.14
[88]
ClayI
Diverse topsoil
Australia
Lab/sieved (2)-air-dried
4606CV
Model trees
0.41
1.2
[95]
ClayI
Diverse subsoil
Australia
Lab/sieved (2)-air-dried
2492CV
Model trees
0.40
1.5
[95]
ClayI
Mineral-organic mixes
Australia
Lab/sieved (0.2)-oven-dried
8CV
PLSR
0.96
0.34
[2]
ClayS
Diverse
Australia
Lab/sieved (2)-air-dried
98CV
PLSR
0.94
118.7
[94]
ClayS
Not mentioned
China
Lab/sieved (2)-air-dried
20CV
CR-MLR
0.84
21.35
[88]
ClayS
Diverse
France
Lab/sieved (2)-air-dried
63CV
CR
0.83
158
[20]
ClayS
Diverse topsoil
Australia
Lab/sieved (2)-air-dried
4606CV
Model trees
0.61
0.8
[95]
ClayS
Diverse subsoil
Australia
Lab/sieved (2)-air-dried
2492CV
Model trees
0.44
1.2
[95]
ClayS
Mineral-organic mixes
Australia
Lab/sieved (0.2)-oven-dried
8CV
PLSR
0.92
0.34
[2]
ClayS
Not mentioned
USA
Lab/sieved (2)-air-dried
178CV
PLSR
0.83
1
[96]
Clay2:1
Argiudoll, hapludolls, and eutrudox
Brazil
Lab/sieved (2)-oven-dried
29CV
CR
0.80
—
[47]
FeT
Not mentioned
France
Lab/sieved (2)-air-dried
99/49
PLSR
0.84
24h
[59]
FeT
Vertisols and alfisols
India
Lab/sieved (2)-air-dried
175/58
PLSR
0.78
0.15h
[92]
FeT
Diverse
Moravia
Lab/sieved (2)-air-dried
97CV
MLR
0.37
12.76h
[97]
FeT
Cambisols and luvisols
South Africa
Lab/sieved (2)-air-dried
123/40
CR-MLR
0.21
9.3
[98]
FeT
Cambisols and luvisols
South Africa
Field/fresh
94/31
CR-MLR
0.23
16.3
[98]
FeT
Cambisols and luvisols
South Africa
Lab/sieved (2)-air-dried
123/40
PLSR
0.43
9.8
[98]
FeT
Cambisols and luvisols
South Africa
Field/fresh
94/31
PLSR
0.43
8.0
[98]
FeH
Typic hapludox and typic eutrudox
Brazil
Lab/sieved (2)-air-dried
206CV
PLSR
0.77
10.7
[99]
FeG
Fe-oxide mixes
USA
Lab/sieved (2)-oven-dried
6CV
CR
0.95
25.3
[86]
C
Entisols, inceptisols, and aridisols
Iran
Lab/sieved (2)-air-dried
30/10
PLSR
0.71
52.4
[67]
C
Not mentioned
France
Lab/sieved (2)-air-dried
99/49
PLSR
0.89
42.1
[59]
C
Diverse topsoil
Iran
Lab/sieved (2)-air-dried
251CV
PLSR
0.58
117
[100]
C
Diverse topsoil
Iran
Lab/sieved (2)-air-dried
251CV
CR
0.35
149
[100]
C
Entisol
Turkey
Lab/sieved (2)-air-dried
502CV
PLSR
0.65
6.80
[89]
C
Entisol
Turkey
Lab/sieved (2)-air-dried
502CV
MARS
0.71
6.20
[89]
C
Diverse
France
Lab/sieved (2)-oven-dried
52CV
PLSR
0.94
38.34
[22]
C
Diverse
France
Lab/sieved (2)-oven-dried
52CV
CR
0.92
52
[22]
C
Diverse
France
Field/fresh
52CV
PLSR
0.77
76.67
[22]
C
Diverse
France
Field/fresh
52CV
CR
0.47
132
[22]
aSoil attributes: clayT, total clay content; clay2:1,2:1 clay mineral content; ClayK, kaolinite content; ClayI, illite content; ClayS, smectite content; FeT, total iron content; FeH, hematite content; FeG, goethite content; C, carbonate content. bsoil type: soil type is not given in the study (not mentioned); the soil samples in the study include more than three soil types (diverse). csample treatment: oven-dried, sieved (<2 mm), and measured under laboratory conditions (lab/sieved (2)-oven-dried); oven-dried, sieved (<0.2 mm), and measured under laboratory conditions (lab/sieved (0.2)-oven-dried); air-dried, sieved (<2 mm), and measured under laboratory conditions (lab/sieved (2)-air-dried); fresh (wet and unprocessed) and measured under field conditions (field/fresh). dncal/nval shows the number of samples used in the spectral calibration and validation, respectively. CV indicates that the validation was conducted independently using a statistical cross-validation technique. epartial least squares regression (PLSR); multivariate adaptive regression splines (MARS); support vector regression (SVR); the discrete wavelet transform-artificial neural networks (DWT-ANN); continuum removal (CR); feature-based multiple linear regression (CR-MLR). fcoefficient of determination. groot mean-square error (g·kg−1). —, RMSE is not given in the study. hmg·kg−1.
6. Conclusions and Future Research Directions
Clay minerals in soils are more complex and less well crystallized than those in sedimentary rocks. Traditional characterization methods such as XRD are usually expensive and time-consuming, whereas VNIR is a quick, cost-efficient, and nondestructive technique for analyzing the soil mineralogic properties of large datasets. The major strength of soil mineralogy studies is that there is a direct relationship between soil minerals and their spectra, since the diagnostic absorption bands of soil minerals lie within the VNIR region. Therefore, the nature of soil mineralogy can be approached through both spectral measurement and multivariate calibration. The spectral measurement is focused on geometric information extracted from several bands (e.g., 350–400, ∼1900, ∼2200, and 2450–2500 nm) that relate to soil minerals. The parameters derived from the continuum removal method are mainly used for mineral identification and prediction. In a multivariate calibration analysis, the dataset contains the entire VNIR domain. The most robust model for soil mineral estimation is selected after understanding the data, data preprocessing, candidate model building, and performance assessment.
Firstly, VNIR has been greatly developed in soil sciences over the past several decades. However, no definitive results on theoretical calculations have yet been found because most soil studies occur on a regional scale so their results are only regionally representative. Thus, it is essential to further develop the theoretical calibrations of VNIR that are more suitable for soil samples worldwide, despite difficulties due to high soil variability across the globe.
Secondly, more field analyses are required for obtaining full potential of VNIR. The in situ data collection in the field is one of the advantages compared with conventional techniques. The heterogeneity of the technical and environmental factors (e.g., soil moisture, soil surface condition, and biological residue) will directly influence the characteristics of the absorption bands, causing increased uncertainty of the spectral measurements. Nevertheless, multivariate calibration models for field data show good or even better mineral prediction than laboratory data. There has been a lack of more systematic studies on the various effects of field sample data and variations in mineralogy, moisture, organic matter, and their interactions. Therefore, future work should focus on these types of studies rather than laboratory spectra.
Thirdly, VNIR may have the potential to help us investigate interactions between soil clay minerals and SOC. Mechanisms of SOC stabilization have attracted increasing interest due to their potential to influence the global carbon cycle. It is widely suggested that soil clay minerals play a central role in capturing and permanently sequestering atmospheric CO2. Both clay content and clay mineral type exert important influences on the carbon sequestration. Because VNIR is capable of characterizing most of carbon- and hydroxyl-related properties, it should allow us to study clay-SOC interactions when combined with the other common or state-of-the-art techniques.
Finally, integrated soil mapping is needed in future large-scale soil analysis. The VNIR spectrum contains integrative information (e.g., mineral composition, SOM, SOC, pH, and moisture) of the soil attributes that reflect the nature of a soil system. Thus, we could use VNIR to map soils. More collaborative and strategic spectral studies are needed to better understand the complete nature of soil [101, 102]. Some global or national spectral libraries [103, 104] have been established to build collaborative networks for soil spectroscopy, but more spectral libraries will facilitate the wider use of VNIR and make global-scale soil monitoring possible.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the Special Funding for Soil Mineralogy (CUG170106), NSF of China (41772032 and 41472041), NSFC for Young Scholars (41402036 and 41602037), NSF of Hubei for Young Scholars (2016CFB183), and Postdoctoral Science Foundation of China (2015M582301). Thanks to Jiacheng Liu and Feng Cheng for their valuable suggestions and Yeqing Liu for his help with sample analysis. Qian Fang and Lulu Zhao acknowledge the China Scholarship Council (CSC) for financial support (201706410017 for Qian Fang and 201706410006 for Lulu Zhao).
RiedelF.DenkM.MüllerI.BarthN.GläßerC.Prediction of soil parameters using the spectral range between 350 and 15,000 nm: a case study based on the permanent soil monitoring program in Saxony, Germany201831518819810.1016/j.geoderma.2017.11.0272-s2.0-85037527235Viscarra RosselR. A.McGlynnandR. N.McBratneyA. B.Determining the composition of mineral-organic mixes using UV–vis–NIR diffuse reflectance spectroscopy20061371708210.1016/j.geoderma.2006.07.0042-s2.0-33845196221ChurchmanG. J.Is the geological concept of clay minerals appropriate for soil science? A literature-based and philosophical analysis20103515–1892794010.1016/j.pce.2010.05.0092-s2.0-78649633498ChurchmanG. J.LoweD. J.HuangP. M.LiY.SumnerM. E.Alteration, formation, and occurrence of minerals in soils201212ndBoca Raton, FL, USACRC PressSinghM.SarkarB.SarkarS.Stabilization of soil organic carbon as influenced by clay mineralogy2018New York, NY, USAElsevier3384ArthurE.Rapid estimation of cation exchange capacity from soil water content201768336537310.1111/ejss.124182-s2.0-85017375178AbdelkaderB.AhmedA.MosterfaB.IsamS.Laboratory study of geotextiles performance on reinforced sandy soil20172761060107010.1007/s12583-015-0621-02-s2.0-85002674519ChurchmanG. J.Game changer in soil science. Functional role of clay minerals in soil201818119910310.1002/jpln.2017006052-s2.0-85039786775HongH.ChengF.YinK.ChurchmanG. J.WangC.Three-component mixed-layer illite/smectite/kaolinite (I/S/K) minerals in hydromorphic soils, south China20151008-91883189110.2138/am-2015-51702-s2.0-84939495644FangQ.ChurchmanG. J.HongH.New insights into microbial smectite illitization in the Permo-Triassic boundary K-bentonites, South China20171409611110.1016/j.clay.2017.01.0292-s2.0-85012177371FangQ.HongH.ZhaoL.Tectonic uplift-influenced monsoonal changes promoted hominin occupation of the Luonan Basin: Insights from a loess-paleosol sequence, eastern Qinling Mountains, central China201716931232910.1016/j.quascirev.2017.05.0252-s2.0-85020941320HongH.ChurchmanG. J.YinK.LiR.LiZ.Randomly interstratified illite–vermiculite from weathering of illite in red earth sediments in Xuancheng, southeastern China2014214-215424910.1016/j.geoderma.2013.10.0042-s2.0-84888205039YinK.HongH.ChurchmanG. J.LiZ.FangQ.Mixed-layer illite-vermiculite as a paleoclimatic indicator in the Pleistocene red soil sediments in Jiujiang, southern China2018In press10.1016/j.palaeo.2017.06.0342-s2.0-85023642445SrodonJ.DritsV. A.McCartyD. K.HsiehJ. C.EberlD. D.Quantitative X-ray diffraction analysis of clay-bearing rocks from random preparations200149651452810.1346/ccmn.2001.0490604HuangX.JiangH.LiY.The role of poorly crystalline iron oxides in the stability of soil aggregate-associated organic carbon in a rice–wheat cropping system201627911010.1016/j.geoderma.2016.05.0112-s2.0-84971221518CowardE. K.ThompsonA. T.PlanteA. F.Iron-mediated mineralogical control of organic matter accumulation in tropical soils201730620621610.1016/j.geoderma.2017.07.0262-s2.0-85026499028BrownD. J.ShepherdK. D.WalshM. G.Dewayne MaysM.ReinschT. G.Global soil characterization with VNIR diffuse reflectance spectroscopy20061323-427329010.1016/j.geoderma.2005.04.0252-s2.0-33646251126NocitaM.StevensA.NoonC.van WesemaelB.Prediction of soil organic carbon for different levels of soil moisture using Vis-NIR spectroscopy2013199374210.1016/j.geoderma.2012.07.0202-s2.0-84875638315Bellon-MaurelV.McBratneyA.Near-infrared (NIR) and mid-infrared (MIR) spectroscopic techniques for assessing the amount of carbon stock in soils–critical review and research perspectives20114371398141010.1016/j.soilbio.2011.02.0192-s2.0-79955961271DufrechouG.GrandjeanG.BourguignonA.Geometrical analysis of laboratory soil spectra in the short-wave infrared domain: clay composition and estimation of the swelling potential2015243-2449210710.1016/j.geoderma.2014.12.0142-s2.0-84920369992GoudgeT. A.RussellJ. M.MustardJ. F.HeadJ. W.BijaksanaS.A 40,000 yr record of clay mineralogy at Lake Towuti, Indonesia: paleoclimate reconstruction from reflectance spectroscopy and perspectives on paleolakes on Mars20171297-880681910.1130/b31569.12-s2.0-85026724974GomezC.LagacherieP.CouloumaG.Continuum removal versus PLSR method for clay and calcium carbonate content estimation from laboratory and airborne hyperspectral measurements2008148214114810.1016/j.geoderma.2008.09.0162-s2.0-56249104807Viscarra RosselR. A.CattleS. R.OrtegaA.FouadY.In situ measurements of soil colour, mineral composition and clay content by vis–NIR spectroscopy2009150325326610.1016/j.geoderma.2009.01.0252-s2.0-64549156055ZhaoL.HongH.LiuJ.Assessing the utility of visible-to-shortwave infrared reflectance spectroscopy for analysis of soil weathering intensity and paleoclimate reconstruction2018In press10.1016/j.palaeo.2017.07.0072-s2.0-85026428274ClarkR. N.Spectroscopy of rocks and minerals, and principles of spectroscopy19993New York, NY, USAJohn Wiley and Sons358ZudeM.Optical monitoring of fresh and processed agricultural crops2008Boca Raton, FL, USACRC PressFarmerV. C.Infra-red spectroscopy in mineral chemistry1975Berlin, GermanySpringer357388BurnsR. G.19932ndCambridge, UKCambridge University PressTappertM.RivardB.GilesD.TappertR.MaugerA.Automated drill core logging using visible and near-infrared reflectance spectroscopy2011106228929610.2113/econgeo.106.2.2892-s2.0-79955140461SieslerH. W.OzakiY.KawataS.HeiseH. M.20081stWeinheim, GermanyJohn Wiley and SonsClarkR. N.KingT.KlejwaM.SwayzeG. A.VergoN.High spectral resolution reflectance spectroscopy of minerals199095B8126531268010.1029/jb095ib08p12653TianY.ZhangJ.YaoX.CaoW.ZhuY.Laboratory assessment of three quantitative methods for estimating the organic matter content of soils in China based on visible/near-infrared reflectance spectra2013202-20316117010.1016/j.geoderma.2013.03.0182-s2.0-84876739026CambuleA. H.RossiterD. G.StoorvogelJ. J.SmalingE. M. A.Building a near infrared spectral library for soil organic carbon estimation in the Limpopo national park, Mozambique2012183-184414810.1016/j.geoderma.2012.03.0112-s2.0-84860551001ClarkR. N.RoushT. L.Reflectance spectroscopy: quantitative analysis techniques for remote sensing applications198489B76329634010.1029/jb089ib07p06329SavitzkyA.GolayM. J.Smoothing and differentiation of data by simplified least squares procedures19643681627163910.1021/ac60214a0472-s2.0-0042553279NorrisK. H.WilliamsP. C.Optimization of mathematical treatments of raw near-infrared signal in the measurement of protein in hard red spring wheat. I. influence of particle size1984612158165BoS.RosselR. A. V.MouazenA. M.WetterlindJ.Chapter five–visible and near infrared spectroscopy in soil science201010710716321510.1016/s0065-2113(10)07005-72-s2.0-77954634711ShepherdK. D.WalshM. G.Development of reflectance spectral libraries for characterization of soil properties200266398899810.2136/sssaj2002.0988ZhengG.JiaoC.ZhouS.ShangG.Analysis of soil chronosequence studies using reflectance spectroscopy20163781881190110.1080/01431161.2016.11637512-s2.0-84963857251HuP.JiangZ.LiuQ.Estimating the concentration of aluminum-substituted hematite and goethite using diffuse reflectance spectrometry and rock magnetism: feasibility and limitations201612164180419410.1002/2015jb0126352-s2.0-84977502404CornellR. M.SchwertmannU.20032ndWeinheim, GermanyWiley-VCH VerlagMurphyR. J.MonteiroS. T.Mapping the distribution of ferric iron minerals on a vertical mine face using derivative analysis of hyperspectral imagery (430-970 nm)201375293910.1016/j.isprsjprs.2012.09.0142-s2.0-84869886422TorrentJ.LiuQ.BloemendalJ.BarronV.Magnetic enhancement and iron oxides in the upper luochuan Loess-paleosol sequence, Chinese Loess plateau20077151570157810.2136/sssaj2006.03282-s2.0-36749032646ZhaoL.HongH.FangQ.Monsoonal climate evolution in southern China since 1.2 Ma: new constraints from Fe-oxide records in red earth sediments from the Shengli section, Chengdu Basin201747311510.1016/j.palaeo.2017.02.0272-s2.0-85013798043Viscarra RosselR. A.BehrensT.Using data mining to model and interpret soil diffuse reflectance spectra20101581465410.1016/j.geoderma.2009.12.0252-s2.0-77953621099BaumgardnerM. F.SilvaL. F.BiehlL. L.StonerE. R.Reflectance properties of soils19863814410.1016/s0065-2113(08)60672-02-s2.0-77954747792SheldonN. D.TaborN. J.Quantitative paleoenvironmental and paleoclimatic reconstruction using paleosols2009951-215210.1016/j.earscirev.2009.03.0042-s2.0-65549107813DemattêJ. A. M.Da Silva TerraF.Spectral pedology: a new perspective on evaluation of soils along pedogenetic alterations2014217-21819020010.1016/j.geoderma.2013.11.0122-s2.0-84891132113NesbittH. W.FedoC. M.YoungG. M.Quartz and feldspar stability, steady and non-steady-state weathering, and petrogenesis of siliciclastic sands and muds1997105217319210.1086/515908Da CruzR. S.Dias FernandesC. M.Nobre VillasR. N.A study of the hydrothermal alteration in Paleoproterozoic volcanic centers, Sao Felix do Xingu region, Amazonian Craton, Brazil, using short-wave infrared spectroscopy201530432433510.1016/j.jvolgeores.2015.09.0052-s2.0-84942626862CudahyT.CaccettaM.ThomasM.Satellite-derived mineral mapping and monitoring of weathering, deposition and erosion201662370211210.1038/srep237022-s2.0-84963649734AbwenyM. S.van RuitenbeekF. J. A.de SmethB.Short-Wavelength Infrared (SWIR) spectroscopy of low-grade metamorphic volcanic rocks of the Pilbara Craton201611712413410.1016/j.jafrearsci.2016.01.0242-s2.0-84957581680DoublierM. P.RoacheT.PotelS.Short-wavelength infrared spectroscopy: a new petrological tool in low-grade to very low-grade pelites201038111031103410.1130/g31272.12-s2.0-78249274888BadíaD.MartíC.CasanovaJ.A quaternary soil chronosequence study on the terraces of the Alcanadre river (semiarid Ebro Basin, NE Spain)2015241-24215816710.1016/j.geoderma.2014.11.0172-s2.0-84913580887HuntG. R.SalisburyJ. W.Visible and near infrared spectra of minerals and rocks. II. carbonates197122330HeiseandH. M.WinzenR.SieslerH. W.OzakiY.KawataS.HeiseH. M.Chemometrics in near-infrared spectroscopy2002Weinheim, GermanyWiley-VCH125162KuhnM.JohnsonK.2013New York, NY, USASpringerHothornT.HornikK.ZeileisA.Unbiased recursive partitioning: a conditional inference framework200615365167410.1198/106186006x1339332-s2.0-33749677657AdelineK. R. M.GomezC.GorrettaN.RogerJ. M.Predictive ability of soil properties to spectral degradation from laboratory Vis-NIR spectroscopy data201728814315310.1016/j.geoderma.2016.11.0102-s2.0-84996555053AbdiH.WilliamsL. J.Principal component analysis20102443345910.1002/wics.1012-s2.0-77957553895ShaoX.BianX.LiuJ.ZhangM.CaiW.Multivariate calibration methods in near infrared spectroscopic analysis20102111662166610.1039/c0ay00421a2-s2.0-78650313708WoldS.MartensH.WoldH.RuheA.KagstromB.The multivariate calibration problem in chemistry solved by the PLS method1982Heidelberg, GermanySpringer VerlagChangC.LairdD. A.MausbachM. J.HurburghC. R.Near-infrared reflectance spectroscopy–principal components regression analyses of soil properties200165248049010.2136/sssaj2001.652480xFidencioP. H.PoppiR. J.de AndradeJ. C.Determination of organic matter in soils using radial basis function networks and near infrared spectroscopy2002453112513410.1016/s0003-2670(01)01506-92-s2.0-0037127460BreimanL.FriedmanJ.StoneC. J.OlshenR. A.1984Belmont, CA, USACRC PressSummersD.LewisM.OstendorfB.ChittleboroughD.Visible near-infrared reflectance spectroscopy as a predictive indicator of soil properties201111112313110.1016/j.ecolind.2009.05.0012-s2.0-79954451216OstovariY.Ghorbani-DashtakiS.BahramiH.Towards prediction of soil erodibility, SOM and CaCO3 using laboratory Vis-NIR spectra: a case study in a semi-arid region of Iran201831410211210.1016/j.geoderma.2017.11.0142-s2.0-85033598504GomezC.LagacherieP.CouloumaG.Regional predictions of eight common soil properties and their spatial structures from hyperspectral Vis–NIR data2012189-19017618510.1016/j.geoderma.2012.05.0232-s2.0-84865060367NouriM.GomezC.GorrettaN.RogerJ. M.Clay content mapping from airborne hyperspectral Vis-NIR data by transferring a laboratory regression model2017298546610.1016/j.geoderma.2017.03.0112-s2.0-85016139261MouazenA. M.KuangB.De BaerdemaekerJ.RamonH.Comparison among principal component, partial least squares and back propagation neural network analyses for accuracy of measurement of selected soil properties with visible and near infrared spectroscopy20101581-2233110.1016/j.geoderma.2010.03.0012-s2.0-77953620835MulderV. L.PlötzeM.de BruinS.Quantifying mineral abundances of complex mixtures by coupling spectral deconvolution of SWIR spectra (2.1–2.4 μm) and regression tree analysis2013207-20827929010.1016/j.geoderma.2013.05.0112-s2.0-84879558525VasquesG. M.DemattêJ. A. M.Viscarra RosselR. A.Ramírez-LópezL.TerraF. S.Soil classification using visible/near-infrared diffuse reflectance spectra from multiple depths2014223–225737810.1016/j.geoderma.2014.01.0192-s2.0-84897111897VapnikV.2013New York, NY, USASpringer-VerlagGomezC.Le BissonnaisY.AnnabiM.BahriH.RaclotD.Laboratory vis–NIR spectroscopy as an alternative method for estimating the soil aggregate stability indexes of Mediterranean soils2013209-210869710.1016/j.geoderma.2013.06.0022-s2.0-84880023002Viscarra RosselR. A.ParLeS: software for chemometric analysis of spectroscopic data2008901728310.1016/j.chemolab.2007.06.0062-s2.0-36749102560AkaikeH.Information theory as an extension of the maximum likelihood PrincipleProceedings of the Second International Symposium on Information Theory1973Akadémia Kiadó, Budapest, HungaryChongI.JunC.Performance of some variable selection methods when multicollinearity is present2005781-210311210.1016/j.chemolab.2004.12.0112-s2.0-21244436700WoldS.SjöströmM.ErikssonL.PLS-regression: a basic tool of chemometrics200158210913010.1016/s0169-7439(01)00155-12-s2.0-0035965476ShiZ.JiW.Viscarra RosselR. A.ChenS.ZhouY.Prediction of soil organic matter using a spatially constrained local partial least squares regression and the Chinese vis-NIR spectral library201566467968710.1111/ejss.122722-s2.0-84934437782BaoN.WuL.YeB.YangK.ZhouW.Assessing soil organic matter of reclaimed soil from a large surface coal mine using a field spectroradiometer in laboratory2017288475510.1016/j.geoderma.2016.10.0332-s2.0-84995560653CarrinoT. A.CrostaA. P.Bemfica ToledoC. L.SilvaA. M.Unveiling the hydrothermal mineralogy of the Chapi Chiara gold prospect, Peru, through reflectance spectroscopy, geochemical and petrographic data201567299315MorellosA.PantaziX.MoshouD.Machine learning based prediction of soil total nitrogen, organic carbon and moisture content by using VIS-NIR spectroscopy201615210411610.1016/j.biosystemseng.2016.04.0182-s2.0-84971620502JiW.ZhouS.JingyiH.ShuoL.In situ measurement of some soil properties in paddy soil using visible and near-infrared spectroscopy201498e015978510.1371/journal.pone.01597852-s2.0-84978967337HermansenC.KnadelM.MoldrupP.GreveM. H.GislumR.de JongeL. W.Visible–near-infrared spectroscopy can predict the clay/organic carbon and mineral fines/organic carbon ratios2016806148610.2136/sssaj2016.05.01592-s2.0-85007137743WaiserT. H.MorganC. L. S.BrownD. J.HallmarkC. T.In situ characterization of soil clay content with visible near-infrared diffuse reflectance spectroscopy200771238939610.2136/sssaj2006.02112-s2.0-33947411704MurphyR. J.SchneiderS.MonteiroS. T.Consistency of measurements of wavelength position from hyperspectral imagery: use of the ferric iron crystal field absorption at ∼900 nm as an indicator of mineralogy20145252843285710.1109/tgrs.2013.22666722-s2.0-84896314555XuS.ZhaoY.WangM.ShiX.Comparison of multivariate methods for estimating selected soil properties from intact soil cores of paddy fields by Vis–NIR spectroscopy2018310294310.1016/j.geoderma.2017.09.0132-s2.0-85029405976ChenJ.ChenJ.WangQ.ZhangY.DingH.HuangZ.Retrieval of soil dispersion using hyperspectral remote sensing201644456357210.1007/s12524-015-0530-92-s2.0-84957955102Volkan BilgiliA.van EsH. M.AkbasF.DurakA.HivelyW. D.Visible-near infrared reflectance spectroscopy for assessment of soil properties in a semi-arid area of Turkey201074222923810.1016/j.jaridenv.2009.08.0112-s2.0-71549129818RajA.ChakrabortyS.DudaB. M.Soil mapping via diffuse reflectance spectroscopy based on variable indicators: an ordered predictor selection approach201831414615910.1016/j.geoderma.2017.10.0432-s2.0-85034087380DudaB. M.WeindorfD. C.ChakrabortyS.Soil characterization across catenas via advanced proximal sensors2017298789110.1016/j.geoderma.2017.03.0172-s2.0-85016395370SarathjithM. C.DasB. S.WaniS. P.SahrawatK. L.Variable indicators for optimum wavelength selection in diffuse reflectance spectroscopy of soils20162671910.1016/j.geoderma.2015.12.0312-s2.0-84953221075ConfortiM.FroioR.MatteucciG.ButtafuocoG.Visible and near infrared spectroscopy for predicting texture in forest soil: an application in southern Italy20158333934710.3832/ifor1221-0072-s2.0-84929322638RobinsonN. J.KitchingM.The 3D distribution of phyllosilicate clay minerals in western Victoria201628415217710.1016/j.geoderma.2016.07.0062-s2.0-84988009611Viscarra RosselR. A.Fine-resolution multiscale mapping of clay minerals in Australian soils measured with near infrared spectra2011116F410.1029/2011jf0019772-s2.0-82355164119GoetzA.ChabrillatS.LuZ.Field reflectance spectrometry for detection of swelling clays at construction sites20015314315510.1002/fact.10152-s2.0-0035596698VašátR.KodešováR.BorůvkaL.KlementA.JakšíkO.GholizadehA.Consideration of peak parameters derived from continuum-removed spectra to predict extractable nutrients in soils with visible and near-infrared diffuse reflectance spectroscopy (VNIR-DRS)2014232-23420821810.1016/j.geoderma.2014.05.0122-s2.0-84901660529BayerA.BachmannM.MüllerA.KaufmannH.A comparison of feature-based MLR and PLS regression techniques for the prediction of three soil constituents in a degraded South African ecosystem201220122097125210.1155/2012/9712522-s2.0-84873862061CamargoL. A.Marques JuniorJ.BarronV.Ferracciu AlleoniL. R.BarbosaR. S.PereiraG. T.Mapping of clay, iron oxide and adsorbed phosphate in Oxisols using diffuse reflectance spectroscopy201525112413210.1016/j.geoderma.2015.03.0272-s2.0-84926433418KhayamimF.WetterlindJ.KhademiH.RobertsonA. H. J.CanoA. F.StenbergB.Using visible and near infrared spectroscopy to estimate carbonates and gypsum in soils in Arid and Subhumid regions of Isfahan, Iran201523315516510.1255/jnirs.11572-s2.0-84940391377FangQ.HongH.ZhaoL.ChengF.YinK.WangC.Climatic implication of authigenic minerals formed during pedogenic weathering processes201843375376910.3799/dqkx.2018.905HongH.FangQ.WangC.ZhaoL.GongN.Constraints of parent magma on altered clay minerals: a case study on the ashes near the Permian-Triassic Boundary in Xinmin section, Guizhou Province2017422161172Viscarra RosselR.The Soil Spectroscopy Group and the development of a global soil spectral library2009204141510.1255/nirn.1131ShiZ.WangQ.PengJ.Development of a national VNIR soil-spectral library for soil classification and prediction of organic matter concentrations20145771671168010.1007/s11430-013-4808-x2-s2.0-84903128662