Next Article in Journal
Exogenous Glomalin-Related Soil Proteins Differentially Regulate Soil Properties in Trifoliate Orange
Next Article in Special Issue
Estimating Leaf Water Content through Low-Cost LiDAR
Previous Article in Journal
Crop Establishment Methods and Integrated Nutrient Management Improve: Part II. Nutrient Uptake and Use Efficiency and Soil Health in Rice (Oryza sativa L.) Field in the Lower Indo-Gangetic Plain, India
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating Soil Properties and Nutrients by Visible and Infrared Diffuse Reflectance Spectroscopy to Characterize Vineyards

by
José Ramón Rodríguez-Pérez
1,*,
Víctor Marcelo
2,
Dimas Pereira-Obaya
1,
Marta García-Fernández
1 and
Enoc Sanz-Ablanedo
1
1
Grupo de Investigación en Geomática e Ingeniería Cartográfica (GEOINCA), Universidad de León, Avenida de Astorga s/n, 24401 Ponferrada, León, Spain
2
Departamento de Ingeniería y Ciencias Agrarias, Universidad de León, Avenida de Astorga s/n, 24401 Ponferrada, León, Spain
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(10), 1895; https://doi.org/10.3390/agronomy11101895
Submission received: 6 September 2021 / Accepted: 15 September 2021 / Published: 22 September 2021
(This article belongs to the Special Issue Advances in Field Spectroscopy in Agriculture)

Abstract

:
Visible, near, and shortwave infrared (VIS-NIR-SWIR) reflectance spectroscopy, a cost-effective and rapid means of characterizing soils, was used to predict soil sample properties for four vineyards (central and north-western Spain). Sieved and air-dried samples were measured using a portable spectroradiometer (350–2500 nm) and compared for pistol grip (PG) versus contact probe (CP) setups. Raw data processed using standard normal variate (SVN) and detrending transformation (DT) were grouped into four subsets (VIS: 350–700 nm; NIR: 701–1000 nm; SWIR: 1001–2500 nm; and full range: 350–2500 nm) in order to identify the most suitable range for determining soil characteristics. The performance of partial least squares regression (PLSR) models in predicting soil properties from reflectance spectra was evaluated by cross-validation. The four spectral subsets and transformed reflectances for each setup were used as PLSR predictor variables. The best performing PLSR models were obtained for pH, electrical conductivity, and phosphorous (R2 values above 0.92), while models for sand, nitrogen, and potassium showed moderately good performances (R2 values between 0.69 and 0.77). The SWIR subset and SVN + DT processing yielded the best PLSR models for both the PG and CP setups. VIS-NIR-SWIR reflectance spectroscopy shows promise as a technique for characterizing vineyard soils for precision viticulture purposes. Further studies will be carried out to corroborate our findings.

1. Introduction

Knowledge of soil properties and mapping is regarded as key to decision making in precision viticulture, mainly because of a growing interest in more environmentally friendly and sustainable practices [1]. Chemical and physical characteristics are especially important in evaluating soil fertility and understanding soil dynamics [2]. Modern viticulture requires the evaluation of a wide range of soil properties in a timely and cost-effective way. However, conventional methods for laboratory analysis of soils are expensive, time-consuming, and non-environmentally friendly (they require the use of chemical reagents), and need a whole range of sophisticated protocols and equipment [3]. Soil assessment using visible (VIS), near infrared (NIR), and shortwave infrared (SWIR) spectroscopy, although it cannot replace laboratory chemical analysis, is fast, cost-effective, environmentally friendly, non-destructive, reproducible, and repeatable analytical technique [4]. It is also easy to use since samples only require minimal preparation, and, furthermore, it requires no chemicals or reagents and so does not generate chemical waste [5]. A single wavelength spectrum may contain comprehensive information that can predict various soil components [6]. Spectroscopic applications to the soil include NIR, VIS-NIR, and mid-infrared (MIR) analyses comprising Fourier transform infrared (FTIR), FTIR-attenuated total reflection (FTIR-ATR), and Raman spectroscopy [3].
Spectroscopic techniques are physical characterisation methods that involve studying electromagnetic wave interaction with the material under consideration in the ultraviolet, VIS, and infrared (IR) wavelengths [7]. Furthermore, spectroscopy, when coupled with multivariate data analysis, has been shown to be a powerful tool for developing quantitative and classification models in many disciplines, including food technology [8], petroleum engineering [9], and soil science [10], as described by Barra et al. [3]. VIS-NIR spectroscopy is an empirical method based on an analysis of diffuse reflectance radiation in relation to a material’s characteristics and the assumption that the concentration of a given constituent is a linear combination of several absorption features [11].
Ben-Dor [12] described the principles and mechanisms of soil–radiation interactions in relation to quantitative remote sensing of soil properties, noting problematic factors that prevent direct spectral analysis of electromagnetic signals and reviewing studies that describe advances in this quantitative method. The same author previously published research focused on the reflectance spectrum in the VIS-NIR-SWIR regions, together with proposals for practical applications [13]. Stenberg et al. [14] comprehensively reviewed the literature on soil VIS and IR diffuse reflectance spectroscopy (including fundamentals, studied soil properties, conditioning factors, calibrations, field analyses, and practical applications), while Kuang et al. [15] reviewed the sensing concept applied to soil properties (basics and brief theory, factors affecting results, and relationship between sensor output and soil properties).
When electromagnetic radiation is directed to a soil sample, it causes individual molecular bonds to vibrate (they bend or stretch), resulting in a characteristic absorption spectrum [15]. The resulting spectrum has a specific shape dependent on soil composition that can be used for physical and chemical analyses [14]. Soil content in carbon (C), nitrogen (N), water, and clay minerals are properties with direct NIR spectral responses that can be attributed to overtones of OH and overtones and/or combinations of C-H + C-H, C-H + C-C, OH+ minerals, and N-H. Moreover, absorption bands in the VIS range (400–700 nm), due to electron excitation, are related to soil colour [15,16]. Numerous studies have used VIS-NIR spectroscopy in an attempt to predict soil content in total and organic C, total N, clay minerals, and water. Other studies have focused specifically on sand and silt content, pH, electrical conductivity (Ec), total content in N, extractable phosphorous (P), extractable potassium (K), extractable calcium (Ca), extractable iron (Fe), extractable sodium (Na), extractable manganese (Mn), extractable magnesium (Mg), and cation exchange capacity (CEC) [16,17,18,19]. However, results for those studies have been typically modest and also highly variable, as they were based on co-variations in constituents that are spectrally active.
The availability of commercial spectroscopic equipment and software packages for multivariate calibration has led to VIS-NIR spectroscopy becoming widely used for soil characterisation purposes. Standards and protocols for reflectance measurements of soils in the laboratory have been proposed by Pimstein et al. [20] and Ben-Dor et al. [21], while Kuang et al. [15] have reviewed several studies of different VIS-NIR reflectance sensors, including laboratory, non-mobile/field (in situ), and mobile/field (online) sensors.
Diffuse reflectance spectra in soil are non-specific, since scatter effects caused by structure result in overlapping absorption features. Therefore, multivariate techniques are required to extract absorption patterns and to correlate spectra with soil properties. Calibration methods for soil applications include linear regression approaches, such as stepwise multiple linear regression (SMLR), principal component regression (PCR), and partial least squares regression (PLSR), and also data mining techniques, such as neural networks (NN), multivariate adaptive regression splines (MARS), boosted regression trees (BRT), random forests (RF), and support vector machines (SVM), along with their combinations [14].
In agriculture, quantitative and qualitative analyses of soil properties yield accurate information to guide the management of soil fertility and productivity through adjusted fertiliser formulations and recommendations [22]. The rapid development of portable and handheld spectrometers allows analyses to be conducted in situ [23]. As a key factor for site-specific management practices, Angelopoulou et al. [24] recently reviewed laboratory and proximal sensing spectroscopy in the VIS, NIR, and SWIR wavelength regions for soil organic matter estimates. MIR spectroscopy and laser diffraction analysis (LDA) have also been demonstrated to be useful for calculating organic matter and clay content in soils [25].
Spectroscopy has previously been applied to viticulture. For vineyards located in Australia, Cozzolino et al. [23] evaluated use of a portable NIR spectrophotometer in the field to predict soil chemical properties, fitting PLSR models with coefficients of determination (R2) that ranged from 0.69 for P to 0.95 for total N content. Muganu et al. [26] demonstrated the great potential of NIR-acoustic optically tuneable filter (AOTF) spectroscopy in assessing grape quality, noting the influence of soil management practices on vine and grape characteristics. Páscoa et al. [27] developed a method for indirect soil differentiation based on grapevine leaf spectra, demonstrating that leaf spectral information can be used to define soil maps for vineyards. For Northern Portugal, Lopo et al. [28] demonstrated the ability of NIR spectroscopy to discriminate between vineyard soil types, showing that water content is not a significant factor in differentiating between soils.
As reported by Marín-González et al. [19], VIS-NIR spectroscopy can be used to detect soil properties using laboratory, in situ, and online measurements. This technique is effective mainly for assessing primary soil properties with direct spectral responses in the VIS-NIR range, e.g., water, C, N, and clay [15], as well as other soil chemical parameters in the laboratory [17]. Few studies, however, have described evaluations of soil properties without direct spectral responses in the VIS-NIR-SWIR range or have compared different approaches to spectral preprocessing and the use of different accessories. Marín-González et al. [19] evaluated models to estimate soil properties without direct spectral responses in the NIR spectroscopy range (CEC, pH, and extractable Ca and Mg), reporting very good accuracy for pH and moderately good accuracy for CEC and Mg. Munnaf et al. [29] explored accuracy improvements to visible NIR spectroscopy estimates of secondary soil properties (pH and extractable K, Mg, Ca and Na) by laboratory fusion approaches, finding that exclusively online spectrum or hybrid models (50% online scanned spectra and laboratory spectra) significantly improved online prediction accuracies. Note, however, that since those works were based on online spectral measurements obtained by specialist industrial-grade instruments mounted in heavy soil-tilling machinery, and so they are not applicable to multi-year crops such as vineyards.
The objectives of this study were (1) to compare spectral signatures of soils as measured in two setups, using a pistol grip (PG) and fibre optic cable, with light provided by an external illuminator lamp, and using a contact probe (CP), with light provided by an internal halogen bulb; and (2) to assess the ability of linear regression models to calculate soil properties (mainly without direct spectral responses in the VIS-NIR-SWIR range) from preprocessed and non-preprocessed spectral data. Thus, two measurement methods (PG and CP) and two modelling approaches (with and without preprocessing) were applied and compared in order to define a suitable protocol to predict vineyard soil composition by VIS-NIR-SWIR spectroscopy.

2. Materials and Methods

2.1. Study Area and Soil Sampling

Soils were sampled in four different commercial vineyards belonging to three Designations of Origin (DOs): Bierzo (northwest Spain), Ribera del Duero (north-central Spain), and Rueda (northwest-central Spain). A total of 12 soil samples were collected from each of the vineyards, yielding 48 samples in total. Table 1 summarises the main characteristics of the sampled sites, which were very diverse in terms of soil textures, crops, and landscapes.
Soil samples were collected in the 0–0.40 m layer between June and August 2015. Soil cores were air-dried and were sieved (10-mesh) by hand selecting fractions <2 mm before chemical analyses, performed in the Instrumental Techniques Laboratory attached to León University (certified by UNE-EN ISO 9001). The following official analytical measurement methods [31] were used: particle-size distribution of clay, silt, and sand (%) by the pipette method, pH at 1:2.5 soil/water suspension, Ec (dS m−1) at 1:5 soil/water suspension, organic matter (%) by the Walkley–Black method, N (%) by total Kjeldahl nitrogen, P extracted with NaHCO3 0.5 M at pH 8.5 by optical spectrometer UV/VIS analysis (mg kg−1), K and Ca extracted with AcONH4 1N at pH 7 by ICP-AES analysis (cmol kg−1), Mn and Fe extracted with DTPA at pH 7.3 by ICP-AES analysis (mg kg−1), and CEC measured by extraction with ClBa 0.1 M by ICP-AES analysis (cmol kg−1).

2.2. Spectral Reflectance Acquisition

Soil samples were air-dried and spread in black soil cores (20 × 20 cm). Spectral reflectances were recorded at 1 nm intervals from 350 nm to 2500 nm using an ASD FieldSpec 4 Portable Spectroradiometer (Analytical Spectral Devices, Inc., Boulder, CO, USA). Measurements were made, using a 1.5 m fibre optic cable (25° field-of-view), in two ways: (1) PG setup, with two tungsten halogen lamps supporting the fibre optic; and (2) CP setup, with an internal halogen bulb attached by cable.
Data were collected following spectroradiometer manufacturer recommendations [32]. Spectral measurements corresponded to reflectance calculated as the ratio of reflected soil sample energy to reflected energy of a reference calibration panel, consisting of a white reflectance panel providing a diffuse homogeneous mix of full-source energy at nearly 100%. Recalibration was performed after each measurement of five soil samples.

2.2.1. PG Setup Measurements

The geometry parameters of measurements (lamp to soil sample and fibre optic to soil sample distances, and the angle between those two distances) were set to ensure homogenous illumination, with the spot area over the sample surface. To ensure a representative spectrum for each soil sample, four reflectance readings (turning the soil core 90° clockwise before each capture) were calculated, each representing the average of 15 individual measurements.

2.2.2. CP Setup Measurements

The CP accessory with an internal halogen bulb allowed the fibre optic to be attached at a fixed measurement angle of 35°, reducing noise caused by shadows and other errors associated with stray light [33]. The sensed spot had a diameter of 10 mm, so measurements were made five times at five different points of the samples and then averaged.

2.2.3. Preprocessing

The spectral signatures were preprocessed to identify outliers, and the spectra measured for each sample were averaged. To identify the most suitable range to estimate soil properties, wavelengths were grouped into four spectral subsets: VIS (350–700 nm), NIR (701–1000 nm), SWIR (1001–2500 nm), and full range (350–2500 nm).
Standard normal variate (SVN) and detrending transformation (DT) were used for scatter correction following previous studies of soil composition estimation by spectroscopy [23,34]. SVN removes multiplicative interferences of scatter and particle size effects from spectral data by centring and scaling each spectral signature [34]. DT removes nonlinear trends in spectroscopic data by calculating a baseline function as the least squares fit of a polynomial to the sample spectrum [34].

2.2.4. Soil Property Estimation by PLSR

We used PLSR to estimate soil properties (predicted variables) from spectral signatures (predictor variables), given that (as explained above) diffuse reflectance spectra are correlated with soil properties. Since soil spectra show an overlap of weak overtones and combinations of fundamental vibrational bands, multivariate calibration methods were required to quantitatively determine soil properties [35]. PLSR is a generalisation of linear multiple regression that reduces a large number of collinear variables (e.g., reflectance values) to a few non-correlated hidden (latent) variables or factors (see Geladi and Kowalski [36] and Wold et al. [37] for comprehensive descriptions of PLSR).
We fitted several models in order to identify the most suitable procedure. The three reflectance datasets considered were non-preprocessed data and SVN and DT processed data. Additionally, in order to fit simpler and more effective models, an independent model was fitted for each dataset considering the following subsets as independent variables in the PLSR: VIS (350–700 nm), NIR (701–1000 nm), SWIR (1001–2500 nm), and the full range (350–2500 nm).
The resulting models were compared regarding requirements to fit a robust PLSR model: a small number of factors, small errors in leave-one-out cross-validation (CV), and a high R2 [38]. Because of the small number of soil samples, we used the leave-one-out CV procedure to validate the regression models. R2 and root mean square error (RMSE) values for CV were calculated to test the prediction accuracy of each model; also calculated for CV were standard error (SE) values. The ratio of performance to deviation (RPD), i.e., the standard deviation (SD) to SE ratio, was used to test the usability of the calibrated models [38], with an RPD value of 2 or more considered appropriate for soil analysis by spectroscopy [35]. Statistics were calculated according to the following expressions:
R 2 = ( n ( x y ) ( x ) ( y ) [ ( n x 2 ( x ) 2 ) ( n y 2 ( y ) 2 ) ] ) 2
where y is the predicted values, z is the measured values, and n is the number of samples;
RMSE = ( y z ) 2 n
where y is the predicted values, z is the measured values, and n is the number of samples;
SE = ( y z ) 2 n 1
where y is the predicted values, z is the measured values, n is the number of samples; and
RPD = SD SE
where SD reflects the SD values of the measured variable, y is the predicted values, z is the measured values, and n is the number of samples.
The PLSR factors used in the models were selected on the basis of the lowest RMSE and highest R2 [39]. The criterion to choose the optimal number of factors was based on RMSE and the explained variance of the model: another factor was added to the model if the RMSE was reduced by >2% and the explained variance increased. The maximum number of factors ultimately selected was seven.

3. Results

3.1. Soil Reflectance Spectra

Soil spectra were mainly dominated by combinations of fundamental vibrational bands for H-C, H-N, and H-O bonds and by weak overtones, especially from the MIR region [35]. The range of reflectance values for the sampled soils and average spectral signatures for the PG and CP setups are shown in Figure 1.
As was expected, the spectral signatures derived from the PG and CP setups were similar, while the reflectance values were higher for the CP setup due to its greater illumination intensity. Reflectance is influenced by the physical structure of soil [35]; the size, shape, and arrangement of particles and voids affect the length of the light transmitted through a soil sample, thereby influencing spectral signatures [40,41]. All spectral signatures followed the typical shape in each wavelength region, i.e., low values in VIS that rise in NIR and SWIR, while showing water absorbance features at around 1400 nm, 1900 nm, and 2200 nm. The 1400–1900 nm absorption bands dominated for water (O–H bonds), even though the peak at 1400 nm was associated with aliphatic C–H and the peak at 1900 nm was associated with amide N–H [15,42]. The spectral shape at 2200 nm was associated with groups such as phenolic O–H, amide N–H, amine N–H, and aliphatic C–H [2]. In sum, the three major reflectance peaks identified were caused by absorbances of O–H bonds of hygroscopically bound water, clay lattices, and various oxides [43].

3.2. Laboratory Analysis

Table 2 shows basic statistics for the chemical and physical properties of the soil samples. Since the soil dataset reflected four different locations with different chemical and physical soil properties, values were very diverse. In fact, the coefficients of variation (CoV) obtained for P, Ca, and Fe were large. Generally, the variability observed in the soil samples for some chemical and physical properties was considered appropriate for spectroscopic calibrations, while the variability of other properties (clay and organic matter) was not great enough to build robust PLSR models.
An important issue in chemometric calibration is collinearity in the analytical values conditioning the validity of the results [23]. Pearson correlations (not shown) were calculated for all the soil properties; the greatest correlations were observed for silt with sand (r = 0.95) and for pH with Ec (r = 0.94), while the lowest correlation was found for Fe with K.

3.3. PLSR Model Predictions

Only variables with R2 values above 0.60 were considered for PLSR in this research. The reference for the other preprocessing results was PLSR prediction results using the full VIS + NIR + SWIR range (350–2500 nm) and non-preprocessed data, as summarised in Table 3. Broadly speaking, the PLSR calibration results indicated good predictions for pH, Ec, and P, and reasonably good predictions for sand, N, K, and Mn. The CP setup models had higher R2 and lower RMSE values for pH, Ec, P, Ca, and Mn, while the PG setup models had higher R2 and lower RMSE values for sand, N, and K. Regarding the number of factors, CP setup models required fewer factors that PG setup models.
The PLSR predictions explained about 90% of the variance in the laboratory analyses of pH, Ec, and P, with RMSE values of 0.340–0.329, 0.01–0.01 dS m−1, and 6.530–5.770 mg kg−1, respectively, for the PG and CP setups. For sand, N, and K, R2 values were better for the PG setup (0.75, 0.68, and 0.65, respectively) than for the CP setup (0.70, 0.62, and 0.64, respectively). RMSE results were very similar; the most different were predictions for sand: 7.678% for the PG setup and 8.327% for the CP setup. Predictions for Mn ranged from R2 = 0.62 (RMSE = 4.399 mg kg−1) for the PG setup to R2 = 0.66 (RMSE = 4.195 mg kg−1) for the CP setup.
Cozzolino and Morón [2] suggest that calibration models developed for soil composition by spectroscopy can be classified according to RPD as poor (<1.6), acceptable (1.6–2.0), or excellent (>2.0). According to this classification, the fitted PLSR models proved excellent for pH, P, and Ca and acceptable for sand, Ec, N, K, and Mn. Chang et al. [35] suggest that spectroscopic prediction models in the intermediate category could be improved using different calibration strategies. The strategy used in this research was SVN and DT preprocessing to achieve models that reduce errors (RMSE and SE) and number of factors and increase R2.
Table 4 shows PLSR results for preprocessed reflectance. Regarding the SVN transformation, R2 increased for all variables except for P, Ec, and N. PLSR performance improved more for the CP setup models. RMSE decreased except for N and K, while the number of factors was also reduced except for N. Regarding DT preprocessing, R2 did not increase except for K. Results were generally better for the PG setup models. RMSE values were maintained or increased except for P and Mn, while the number of factors was reduced except for N. The reduction in the number of factors was less for PG setup models. For SVN preprocessing, R2 values increased except for Ec and pH, which remained constant, while RMSE values decreased. Results were better for the models based on the CP setup, while the number of factors was also reduced, with the exception of the PG models estimating Mn (+2 factors) and the CP setup models estimating N and P (+1 factor). For DT preprocessing, although not significantly greater, R2 and RMSE values were better for the CP setup than the PG setup. The main improvement was the simplification of the models in reducing the number of factors. Finally, applying SVN + DT, R2 values increased and RMSE values decreased, while the number of factors was reduced, with the exception of N (7 factors).

4. Discussion

4.1. Soil Property Predictions

Cross-validation results for the PLSR models were different for the three particle-size distributions (clay, silt, and sand). For sand, results were satisfactory (R2 = 0.75 and R2 = 0.70 for the PG and CP data, respectively), and also corroborated other published results [17,39,44]. For clay, however, results were quite poor (R2 = 0.53 and R2 = 0.51 for the PG and CP data, respectively), and likewise for silt (R2 = 0.51 and R2 = 0.49 for the PG and CP data, respectively). Those unexpectedly low R2 values may be due to narrow variability in clay content (min = 10 and max = 32) and silt content (min = 14 and max = 42) of the analysed soils.
Soil content in N was estimated by spectroscopy because it is quite sensitive to IR radiation. The R2 values obtained ranged from R2 = 0.68 (RMSE = 0.017%) to R2 = 0.62 (RMSE = 0.018%) for the PG and CP setups, respectively, lower than the values of R2 = 0.80–0.98 cited elsewhere [45] and the R2 = 0.92 (SE = 2.19) obtained by Cozzolino et al. [23]. Our result can be explained by the fact that N estimation by spectroscopy is soil-dependent, due mainly to varying carbonate contents [46]. While MIR-ATR spectroscopy can, in fact, predict nitrate concentration in soil pastes by direct measurement, prediction accuracy is strongly conditioned for water and soil constituents [47,48,49].
Previous research has reported the ability of soil reflectance spectroscopy to accurately determine soil organic matter [24]. However, our results for organic matter were R2 = 0.29 (RMSE = 0.438%) and R2 = 0.27 (RMSE = 0.445%) for the PG and CP setups, respectively. These poor results may be due to the fact that the analysed soils have low organic matter content (0.37–2.40%) and high sand content (24–76%). In fact, spectroscopic predictions of organic matter are poorly accurate in soils with low C content [50] and high sand content [14].
Soil pH is a key factor for agriculture as an important fertility regulator of nutrient solubility and plant root development, biological activity, decomposition, mineralisation, etc. Because pH is a soil property with no direct spectral responses in the NIR spectroscopy range [15], calibrations rarely perform better than an RMSE of one-third or half a pH unit [14]. However, soil pH has been predicted quite successfully in several studies [19,51,52].
The pH prediction performance of our PLSR models was excellent (RPD > 3.3) for both PG and CP reflectances (R2 = 0.92 (RMSE = 0.340) and R2 = 0.92 (RMSE = 0.329), respectively). Our results were similar to those reported by Kuang et al. (RMSE = 0.36; RPD = 2.02) and better than those reported by Sorenson et al. [52] (R2 = 0.68) and by Marín-González et al. [19] (R2 = 0.86; RPD = 2.69). Results for pH predictions may be explained by co-variation to spectrally active soil constituents such as organic matter and clay [35] or by soil mineralogy and carbonate content [15]. Note that pH calibrations tend to vary from one dataset to another because they reflect different scenarios.
Our results for Ca calibrations were reasonably good, at R2 = 0.87 (RMSE = 2.332 cmol kg−1) and R2 = 0.88 (RMSE = 2.141 cmol kg−1) for the PG and CP setups, respectively, while RPD > 2 indicated excellent accuracy. Those results improved on those obtained elsewhere: R2 = 0.75 (RMSE = 4.00 cmol kg−1) by Chang et al. [35]; R2 = 0.72 (RMSE = 4.20 cmol kg−1) by Islam et al. [6]; and R2 = 0.67 (RMSE = 3.89 cmol kg−1) by Soriano-Disla et al. [53]. Similar results to ours were obtained by Shepherd and Walsh [54] (R2 = 0.88), Cozzolino and Morón [2] (R2 = 0.90), and Marín-González et al. [19] (R2 = 0.89; RMSE = 22.05).
As for Ec, the results of the PLSR models indicate accurate predictions for both the PG setup (R2 = 0.89; RMSE = 0.011; RPD = 2.73) and the CP setup (R2 = 0.90; RMSE = 0.011; RPD = 2.73). Our results corroborate those of Lei et al. [55], who fitted a PLSR model with R2 = 0.69 (RMSE = 0.039) for semi-arid grasslands, and Farifteh et al. [56], who fitted models that ranged from R2 = 0.80 (RMSE = 0.070) to R2 = 0.80 (RMSE = 0.060) for sandy and agricultural areas, while Mashimbye et al. [57] reported PLSR models with R2 = 0.65–0.85 (validation) for soil in South Africa.
According to Marín-González et al. [19], K is a difficult property to estimate with NIR spectroscopy. Islam et al. [6] reported poor predictions for extractable K for calibrations (R2 = 0.29) using PCR. Using PLSR, Volkan et al. [18] reported K soil content estimates (depending on the validation set) with R2 values between 0.32 and 0.25 (RMSE 0.21 and 0.22, respectively). Our models predicted extractable K with moderate accuracy for both the PG setup (R2 = 0.65; RMSE = 0.086; RPD = 1.61) and the CP setup (R2 = 0.64; RMSE = 0.088; RPD = 1.59). Those accuracy levels corroborate Zornoza et al. [42], who reported PLSR models with R2 = 0.79 (RMSE = 0.11 g kg−1; RPD = 2.19), but using soil with high K content (mean 0.60 g kg−1).
Our extractable Mn estimates were moderately accurate for both PG (R2 = 0.62; RMSE = 4.399 mg kg−1; RPD = 1.59) and CP (R2 = 0.66; RMSE = 4.195 mg kg−1; RPD = 1.67) setups. There are few references in the literature for K content estimations using VIR-NIR, although Chang et al. [35] built a PCR model resulting in R2 = 0.70 (RMSE = 56.40 mg kg−1; RPD = 1.79).
Since our PG and CP setup models had comparable predictive capacities (R2) and accuracies (RMSE), it can be concluded that the CP setup, more versatile for field measurements, is preferable for soil property estimates by VIS-NIR-SWIR spectroscopy. Note that while Rosero-Blasova et al. [33] reported a PG setup to perform better than a CP setup, they detected no statistically significant differences between the two setups.

4.2. PLSR Model Performance

Figure 2 shows distributions of the weighted regression coefficients over the full spectral range for both the PG and CP setups and the considered soil properties (to highlight differences, the regression coefficients for each soil property are offset by 3.0 units). Evident are several peaks in wavelength bands located in the VIS and NIR regions, attributable to colour, water, organic matter, and clay minerals [16]. Regarding sand, K, P, and Mn, the main peaks in the VIS range are associated with the blue and green regions around 450 nm and 550 nm, respectively, demonstrating that colour contributes similarly to predicting those properties. Mouazen et al. [16] reported a similar distribution of regression coefficients to ours, identifying the spectral range between 1800 nm and 2450 nm as the most active for P and K estimates. As for pH and Ec, these are mainly associated with the blue and green regions, denoting the influence of Fe oxides associated with clay minerals [58]. Predictions of N content are little affected by colour, while Ca predictions are influenced in the red region.
As was expected, regression coefficient distributions were very similar for both PG and CP setups, thereby corroborating measurement and prediction consistency between both. The main difference was reported for N estimates, which can be attributed to PLSR algorithm calculations instead of actual differences in spectral signatures.
Figure 3 shows the R2 values for the PLSR models obtained for the different data subsets, different pre-processing approaches, and PG and CP setups. Using only part of the spectrum (VIS, NIR, or SWIR), general trends for R2 were similar: R2 values with SWIR were the best, followed by R2 values with VIS, then R2 values with NIR. RMSE values were highest with NIR, while the lowest values were obtained with SWIR (see Table 4 above). The main disadvantage of using SWIR subsets was that some models (e.g., those for N and Mn) needed a greater number of factors than the VIS and NIR subsets, which potentially delays computational calculations.

4.3. Data Preprocessing

Table 5 shows the best fitting models for each soil property, each type of preprocessing, each spectral subset, and each setup. Generally, sand, pH, Ec, N, P, and K were best predicted with models using the SWIR subset, Ca with models using the VIS subset, and Mn with models using the full spectrum. The best models for sand, Ec, K, Ca, and Mn were obtained for SVN + DT preprocessing, while models for N and P only required SVN preprocessing, and models for pH required no preprocessing.
No one specific kind of preprocessing ensures the effectiveness of models. Spectral signatures of soils are influenced by chemical composition and structural properties that produce non-linear light scattering effects. Regression model performance depends on the soil dataset, the analysed soil property, and the variability of the data [59], so a specific model needs to be fitted that reflects each scenario. Furthermore, it has been reported that spectral preprocessing has a minor influence on results when PLSR models are used [60].
Stenberg et al. [14] report that SVN combined with DT is one of the more commonly used means of improving PLSR performance, as this approach usually enhances weak soil spectral signals. In our research, while SVN + DT increased R2 and reduced RMSE (see Table 3 and Table 4), improvement depended on the studied soil property, and was not so great probably because the raw reflectance data were quite stable and consistent. Other authors [33,61] report, for VIS-NIR spectroscopy, that preprocessing of spectral samples is data-specific, so no single or combination technique is generally applicable to preprocessing. In fact, different preprocessing methods should be used for different calibration techniques, different datasets, and different soil conditions [59].
Table 5 confirms that the predictive performance of soil property PLSR spectroscopic models varies with different kinds of preprocessing. Furthermore, use of different accessories results in different illumination setups and observation geometries that condition measurement and that consequently may affect the performance of models [21]. Model effectiveness is also probably conditioned by variability in the data [59]. In fact, for properties where standard deviations are greater, more variance is explained and greater accuracy is achieved.
Figure 4 represents Pearson coefficient values reflecting correlations between soil properties and wavelengths. The correlograms grouped by correlation structure are Fe and Mn; N and organic matter; pH and Ec; and sand and K. Analysing the correlograms, within groups, the correlation structure is quite redundant; only Fe and organic matter have direct optical features, while predictions for the remaining properties are based on spurious correlations [62]. Patterns in the groups can be attributed to dominant chemical characteristics (e.g., iron oxides and clay minerals in the group consisting of Fe and Mn) and to the aggregate effect of several optically active minerals [62]. Regression models based on spurious correlations depend on underlying geology and soil parameters and so only are useful for our studied plots.
Interpretation of the regression coefficient curves and correlograms (Figure 3 and Figure 4, respectively) is complicated due to the complexity in overlapping soil constituent absorption patterns. The studied chemical properties do not have direct spectral responses in the considered spectral regions. The prediction of these properties, namely, sand, pH, Ec, N, P, K, Ca, and Mn, can be attributed to locally present co-variations in spectrally active constituents (mainly organic C and clay minerals). Furthermore, correlations of some soil properties with NIR spectroscopy are still unknown and so require further investigation [16]. In fact, Miller [63] acknowledged that it is difficult to identify relevant effects in the NIR spectrum based on chemistry and spectroscopy of samples alone. Therefore, further studies are needed to understand why, in our study and using VIS-NIR-SWIR spectroscopy, properties were estimated with excellent accuracy (pH, Ec, P and Ca) and acceptable accuracy (sand, N, and Mn).
Our results suggest that it is possible to estimate variables such as sand, pH, Ec, N, P, K, Ca, and Mn that are optically non-active chemical properties with featureless spectra, because those elements are bonded to spectrally active soil components, mainly iron oxides, organic matter, and clay minerals, in such a way that the bonds constitute a key predictive mechanism [62]. Similar conclusions have been published by Martínez-Carreras et al. [64] and Wu et al. [65].

5. Conclusions

Vineyard soil parameters were calculated by relating spectral signatures and laboratory analytical determinations using PLSR. Reflectance measurements were made using PG and CP setups. Our findings suggest that proximal soil spectroscopy is a useful technique for soil characterisation and monitoring. The great advantage of the spectroscopic approach is that it is cost-effective and rapid, although prediction accuracy is less than for laboratory analyses. The predictive capacity (R2) and accuracy (RMSE) of the PLSR models depends on setup (PG or CP), preprocessing (SVN and/or DT), spectral subset (VIS, NIR, SWIR, or full spectrum), and individual soil properties. The best predictions, with R2 values above 0.915, were obtained for pH, Ec, and P, while moderately accurate predictions, with R2 values of 0.69 to 0.77, were obtained for sand, N, and K.
In conclusion, PLSR models can be useful for monitoring overall changes in soil properties. Further studies aimed at more effective precision viticulture practices will focus on vineyard soil characterisation using VIS-NIR-SWIR spectroscopy combined with geographical information system (GIS) data.

Author Contributions

Conceptualisation, M.G.-F., V.M. and J.R.R.-P.; methodology, M.G.-F., V.M. and J.R.R.-P.; software, E.S.-A.; formal analysis, V.M., E.S.-A. and J.R.R.-P.; investigation, D.P.-O.; writing—original draft preparation, M.G.-F., D.P.-O. and J.R.R.-P.; writing—review and editing, D.P-O. and J.R.R.-P.; supervision and project administration, J.R.R.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Education Department of the Junta de Castilla y León, grant number LE112G18.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was supported by funding from the Education Department of the Junta de Castilla y León-Spain (grant number LE112G18) under a call for research by recognised research groups attached to public universities starting in 2018 (Order 20 November 2017). Marta García Fernández gratefully acknowledges financial support provided by the Fundación Carolina Rodríguez and the Universidad de León.

Conflicts of Interest

The authors declare no conflict of interest. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Mashalaba, L.; Galleguillos, M.; Seguel, O.; Poblete-Olivares, J. Predicting Spatial Variability of Selected Soil Properties Using Digital Soil Mapping in a Rainfed Vineyard of Central Chile. Geoderma Reg. 2020, 22, e00289. [Google Scholar] [CrossRef]
  2. Cozzolino, D.; Morón, A. The Potential of Near-Infrared Reflectance Spectroscopy to Analyse Soil Chemical and Physical Characteristics. J. Agric. Sci. 2003, 140, 65–71. [Google Scholar] [CrossRef]
  3. Barra, I.; Haefele, S.M.; Sakrabani, R.; Kebede, F. Soil Spectroscopy with the Use of Chemometrics, Machine Learning and Pre-Processing Techniques in Soil Diagnosis: Recent Advances—A Review. TrAC Trends Anal. Chem. 2021, 135, 116166. [Google Scholar] [CrossRef]
  4. Nocita, M.; Stevens, A.; van Wesemael, B.; Aitkenhead, M.; Bachmann, M.; Barthès, B.; Ben Dor, E.; Brown, D.J.; Clairotte, M.; Csorba, A.; et al. Chapter Four—Soil Spectroscopy: An Alternative to Wet Chemistry for Soil Monitoring. In Advances in Agronomy; Sparks, D.L., Ed.; Academic Press: Cambridge, MA, USA, 2015; Volume 132, pp. 139–159. [Google Scholar]
  5. Osborne, B.G. Near-Infrared Spectroscopy in Food Analysis. In Encyclopedia of Analytical Chemistry; American Cancer Society: Atlanta, GA, USA, 2006; ISBN 978-0-470-02731-8. [Google Scholar]
  6. Islam, K.; Singh, B.; McBratney, A. Simultaneous Estimation of Several Soil Properties by Ultra-Violet, Visible, and near-Infrared Reflectance Spectroscopy. Soil Res. 2003, 41, 1101–1114. [Google Scholar] [CrossRef]
  7. Zumdahl, S.S.; Zumdahl, S.A.; DeCoste, D.J. Molecular spectroscopy. In Chemistry; Cengage Learning: Boston, MA, USA; University of Illinois: Boston, MA, USA, 2017; ISBN 978-1-305-95740-4. [Google Scholar]
  8. Kharbach, M.; Marmouzi, I.; Kamal, R.; Yu, H.; Barra, I.; Cherrah, Y.; Alaoui, K.; Heyden, Y.V.; Bouklouze, A. Extra Virgin Argan Oils’ Shelf-Life Monitoring and Prediction Based on Chemical Properties or FTIR Fingerprints and Chemometrics. Food Control 2021, 121, 107607. [Google Scholar] [CrossRef]
  9. Barra, I.; Kharbach, M.; Qannari, E.M.; Hanafi, M.; Cherrah, Y.; Bouklouze, A. Predicting Cetane Number in Diesel Fuels Using FTIR Spectroscopy and PLS Regression. Vib. Spectrosc. 2020, 111, 103157. [Google Scholar] [CrossRef]
  10. Martínez-España, R.; Bueno-Crespo, A.; Soto, J.; Janik, L.J.; Soriano-Disla, J.M. Developing an Intelligent System for the Prediction of Soil Properties with a Portable Mid-Infrared Instrument. Biosyst. Eng. 2019, 177, 101–108. [Google Scholar] [CrossRef]
  11. Ben-Dor, E.; Banin, A. Near-Infrared Analysis as a Rapid Method to Simultaneously Evaluate Several Soil Properties. Soil Sci. Soc. Am. J. 1995, 59, 364–372. [Google Scholar] [CrossRef]
  12. Ben-Dor, E. Quantitative remote sensing of soil properties. In Advances in Agronomy; Academic Press: Cambridge, MA, USA, 2002; Volume 75, pp. 173–243. [Google Scholar]
  13. Ben-Dor, E.; Irons, J.R.; Epema, G.F.; Rencz, A.N. Soil reflectance. In Remote Sensing for the Earth Sciences: Manual of Remote Sensing 3/3; John Wiley & Sons: Hoboken, NJ, USA, 1999; pp. 111–188. [Google Scholar]
  14. Stenberg, B.; Viscarra Rossel, R.A.; Mouazen, A.M.; Wetterlind, J. Visible and Near Infrared Spectroscopy in Soil Science. Adv. Agron. 2010, 107, 163–215. [Google Scholar] [CrossRef] [Green Version]
  15. Kuang, B.; Mahmood, H.S.; Quraishi, Z.; Hoogmoed, W.B.; Mouazen, A.M.; van Henten, E. Sensing Soil Properties in the Laboratory, in Situ, and on-Line: A Review. Adv. Agron. 2012, 114, 155–223. [Google Scholar] [CrossRef]
  16. Mouazen, A.M.; Kuang, B.; De Baerdemaeker, J.; Ramon, H. 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 Spectroscopy. Geoderma 2010, 158, 23–31. [Google Scholar] [CrossRef]
  17. Viscarra-Rossel, R.A.; Walvoort, D.J.J.; McBratney, A.B.; Janik, L.J.; Skjemstad, J.O. Visible, near Infrared, Mid Infrared or Combined Diffuse Reflectance Spectroscopy for Simultaneous Assessment of Various Soil Properties. Geoderma 2006, 131, 59–75. [Google Scholar] [CrossRef]
  18. Bilgili, A.V.; Van Es, H.M.; Akbas, F.; Durak, A.; Hively, W.D. Visible-near Infrared Reflectance Spectroscopy for Assessment of Soil Properties in a Semi-Arid Area of Turkey. J. Arid. Env. 2010, 74, 229–238. [Google Scholar] [CrossRef]
  19. Marín-González, O.; Kuang, B.; Quraishi, M.Z.; Munóz-García, M.Á.; Mouazen, A.M. On-Line Measurement of Soil Properties without Direct Spectral Response in near Infrared Spectral Range. Soil Tillage Res. 2013, 132, 21–29. [Google Scholar] [CrossRef] [Green Version]
  20. Pimstein, A.; Notesco, G.; Ben-Dor, E. Performance of Three Identical Spectrometers in Retrieving Soil Reflectance under Laboratory Conditions. Soil Sci. Soc. Am. J. 2011, 75, 746–759. [Google Scholar] [CrossRef]
  21. Ben-Dor, E.; Ong, C.; Lau, I.C. Reflectance Measurements of Soils in the Laboratory: Standards and Protocols. Geoderma 2015, 245–246, 112–124. [Google Scholar] [CrossRef]
  22. Silveira, M.L.; Kohmann, M.M. Chapter 3—Maintaining soil fertility and health for sustainable pastures. In Management Strategies for Sustainable Cattle Production in Southern Pastures; Rouquette, M., Aiken, G.E., Eds.; Academic Press: Cambridge, MA, USA, 2020; pp. 35–58. ISBN 978-0-12-814474-9. [Google Scholar]
  23. Cozzolino, D.; Cynkar, W.U.; Dambergs, R.G.; Shah, N.; Smith, P. In Situ Measurement of Soil Chemical Composition by Near-Infrared Spectroscopy: A Tool Toward Sustainable Vineyard Management. Commun. Soil Sci. Plant Anal. 2013, 44, 1610–1619. [Google Scholar] [CrossRef]
  24. Angelopoulou, T.; Balafoutis, A.; Zalidis, G.; Bochtis, D. From Laboratory to Proximal Sensing Spectroscopy for Soil Organic Carbon Estimation—A Review. Sustainability 2020, 12, 443. [Google Scholar] [CrossRef] [Green Version]
  25. Thomas, C.L.; Hernandez-Allica, J.; Dunham, S.J.; McGrath, S.P.; Haefele, S.M. A Comparison of Soil Texture Measurements Using Mid-Infrared Spectroscopy (MIRS) and Laser Diffraction Analysis (LDA) in Diverse Soils. Sci. Rep. 2021, 11, 16. [Google Scholar] [CrossRef]
  26. Muganu, M.; Paolocci, M.; Gnisci, D.; Barnaba, F.E.; Bellincontro, A.; Mencarelli, F.; Grosu, I. Effect of Different Soil Management Practices on Grapevine Growth and on Berry Quality Assessed by NIR-AOTF Spectroscopy. Acta Hortic 2013, 978, 125. [Google Scholar] [CrossRef]
  27. Páscoa, R.N.M.J.; Lopo, M.; Teixeira dos Santos, C.A.; Graça, A.R.; Lopes, J.A. Exploratory Study on Vineyards Soil Mapping by Visible/near-Infrared Spectroscopy of Grapevine Leaves. Comput. Electron. Agric. 2016, 127, 15–25. [Google Scholar] [CrossRef]
  28. Lopo, M.; Teixeira dos Santos, C.A.; Páscoa, R.N.M.J.; Graça, A.R.; Lopes, J.A. Near Infrared Spectroscopy as a Tool for Intensive Mapping of Vineyards Soil. Precis. Agric 2018, 19, 445–462. [Google Scholar] [CrossRef]
  29. Munnaf, M.A.; Nawar, S.; Mouazen, A.M. Estimation of Secondary Soil Properties by Fusion of Laboratory and On-Line Measured Vis–NIR Spectra. Remote Sens. 2019, 11, 2819. [Google Scholar] [CrossRef] [Green Version]
  30. IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015 International Soil Classification System Fort Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports No. 106; FAO: Rome, Italy, 2015. [Google Scholar]
  31. Ministerio de Agricultura Pesca y Alimentación–MAPA. Métodos Oficiales de Análisis; Secretaría General Técnica del Ministerio de Agricultura, Pesca y Alimentación: Madrid, Spain, 1994; Volume III, ISBN 978-84-491-0003-1.
  32. ASD Inc. FieldSpect4 User Manual. Available online: https://www.malvernpanalytical.com/en/learn/knowledge-center/user-manuals/fieldspec-4-user-guide (accessed on 25 June 2021).
  33. Rosero-Vlasova, O.A.; Pérez-Cabello, F.; Montorio Llovería, R.; Vlassova, L. Assessment of Laboratory VIS-NIR-SWIR Setups with Different Spectroscopy Accessories for Characterisation of Soils from Wildfire Burns. Biosyst. Eng. 2016, 152, 51–67. [Google Scholar] [CrossRef]
  34. Barnes, R.J.; Dhanoa, M.S.; Lister, S.J. Standard Normal Variate Transformation and De-Trending of Near-Infrared Diffuse Reflectance Spectra. Appl. Spectrosc. 1989, 43, 772–777. [Google Scholar] [CrossRef]
  35. Chang, C.W.; Laird, D.A.; Mausbach, M.J.; Hurburgh, C.R. Near-Infrared Reflectance Spectroscopy–Principal Components Regression Analyses of Soil Properties. Soil Sci. Soc. Am. J. 2001, 65, 480–490. [Google Scholar] [CrossRef] [Green Version]
  36. Geladi, P.; Kowalski, B.R. Partial Least-Squares Regression: A Tutorial. Anal. Chim. Acta 1986, 185, 1–17. [Google Scholar] [CrossRef]
  37. Wold, S.; Sjöström, M.; Eriksson, L. PLS-Regression: A Basic Tool of Chemometrics. Chemom. Intell. Lab. Syst. 2001, 58, 109–130. [Google Scholar] [CrossRef]
  38. Williams, P.; Dardenne, P.; Flinn, P. Tutorial: Items to Be Included in a Report on a near Infrared Spectroscopy Project. J. Near Infrared Spectrosc. JNIRS 2017, 25, 85–90. [Google Scholar] [CrossRef]
  39. Conforti, M.; Froio, R.; Matteucci, G.; Buttafuoco, G. Visible and near Infrared Spectroscopy for Predicting Texture in Forest Soil: An Application in Southern Italy. Iforest—Biogeosciences For. 2015, 8, 339. [Google Scholar] [CrossRef]
  40. Wetzel, D.L. Near-Infrared Reflectance Analysis. Anal. Chem. 1983, 55, 1165A–1176A. [Google Scholar] [CrossRef]
  41. Barthès, B.G.; Brunet, D.; Ferrer, H.; Chotte, J.L.; Feller, C. Determination of Total Carbon and Nitrogen Content in a Range of Tropical Soils Using near Infrared Spectroscopy: Influence of Replication and Sample Grinding and Drying. J. Near Infrared Spectrosc. 2006, 14, 341–348. [Google Scholar] [CrossRef]
  42. Zornoza, R.; Guerrero, C.; Mataix-Solera, J.; Scow, K.M.; Arcenegui, V.; Mataix-Beneyto, J. Near Infrared Spectroscopy for Determination of Various Physical, Chemical and Biochemical Properties in Mediterranean Soils. Soil Biol. Biochem. 2008, 40, 1923–1930. [Google Scholar] [CrossRef] [Green Version]
  43. Martin, P.D.; Malley, D.F.; Manning, G.; Fuller, L. Determination of Soil Organic Carbon and Nitrogen at the Field Level Using Near-Infrared Spectroscopy. Can. J. Soil. Sci. 2002, 82, 413–422. [Google Scholar] [CrossRef]
  44. Curcio, D.; Ciraolo, G.; D’Asaro, F.; Minacapilli, M. Prediction of Soil Texture Distributions Using VNIR-SWIR Reflectance Spectroscopy. Procedia Environ. Sci. 2013, 19, 494–503. [Google Scholar] [CrossRef] [Green Version]
  45. Stevens, A.; van Wesemael, B.; Bartholomeus, H.; Rosillon, D.; Tychon, B.; Ben-Dor, E. Laboratory, Field and Airborne Spectroscopy for Monitoring Organic Carbon Content in Agricultural Soils. Geoderma 2008, 144, 395–404. [Google Scholar] [CrossRef] [Green Version]
  46. Du, C.; Zhou, J. Evaluation of Soil Fertility Using Infrared Spectroscopy—A Review. In Climate Change, Intercropping, Pest Control and Beneficial Microorganisms; Lichtfouse, E., Ed.; Sustainable Agriculture Reviews; Springer: Dordrecht, The Netherlands, 2009; pp. 453–483. ISBN 978-90-481-2716-0. [Google Scholar]
  47. Linker, R. Waveband Selection for Determination of Nitrate in Soil Using Mid-Infrared Attenuated Total Reflectance Spectroscopy. Appl. Spectrosc. 2004, 58, 1277–1281. [Google Scholar] [CrossRef]
  48. Linker, R.; Shmulevich, I.; Kenny, A.; Shaviv, A. Soil Identification and Chemometrics for Direct Determination of Nitrate in Soils Using FTIR-ATR Mid-Infrared Spectroscopy. Chemosphere 2005, 61, 652–658. [Google Scholar] [CrossRef]
  49. Borenstein, A.; Linker, R.; Shmulevich, I.; Shaviv, A. Determination of Soil Nitrate and Water Content Using Attenuated Total Reflectance Spectroscopy. Appl. Spectrosc. 2006, 60, 1267–1272. [Google Scholar] [CrossRef]
  50. Sithole, N.J.; Ncama, K.; Magwaza, L.S. Robust Vis-NIRS Models for Rapid Assessment of Soil Organic Carbon and Nitrogen in Feralsols Haplic Soils from Different Tillage Management Practices. Comput. Electron. Agric. 2018, 153, 295–301. [Google Scholar] [CrossRef]
  51. Kuang, B.; Tekin, Y.; Mouazen, A.M. Comparison between Artificial Neural Network and Partial Least Squares for On-Line Visible and near Infrared Spectroscopy Measurement of Soil Organic Carbon, PH and Clay Content. Soil Tillage Res. 2015, 146, 243–252. [Google Scholar] [CrossRef]
  52. Sorenson, P.T.; Small, C.; Tappert, M.C.; Quideau, S.A.; Drozdowski, B.; Underwood, A.; Janz, A. Monitoring Organic Carbon, Total Nitrogen, and pH for Reclaimed Soils Using Field Reflectance Spectroscopy. Can. J. Soil. Sci. 2017, 97, 241–248. [Google Scholar] [CrossRef]
  53. Soriano-Disla, J.M.; Janik, L.J.; Allen, D.J.; McLaughlin, M.J. Evaluation of the Performance of Portable Visible-Infrared Instruments for the Prediction of Soil Properties. Biosyst. Eng. 2017, 161, 24–36. [Google Scholar] [CrossRef]
  54. Shepherd, K.D.; Walsh, M.G. Development of Reflectance Spectral Libraries for Characterization of Soil Properties. Soil Sci. Soc. Am. J. 2002, 66, 988–998. [Google Scholar] [CrossRef]
  55. Lei, S.; Bao, N.; Liu, S.; Liu, X. Modelling and Predicting of Soil Electrical Conductivity and pH from Semi-Arid Grassland Using VIS-NIR Spectroscopy Technology. In Proceedings of the Computer and Computing Technologies in Agriculture; Li, D., Zhao, C., Eds.; Springer International Publishing: Cham, Switzerland, 2019; pp. 442–453. [Google Scholar]
  56. Farifteh, J.; Van der Meer, F.; Atzberger, C.; Carranza, E.J.M. Quantitative Analysis of Salt-Affected Soil Reflectance Spectra: A Comparison of Two Adaptive Methods (PLSR and ANN). Remote Sens. Environ. 2007, 110, 59–78. [Google Scholar] [CrossRef]
  57. Mashimbye, Z.E.; Cho, M.A.; Nell, J.P.; De Clercq, W.P.; Van niekerk, A.; Turner, D.P. Model-Based Integrated Methods for Quantitative Estimation of Soil Salinity from Hyperspectral Remote Sensing Data: A Case Study of Selected South African Soils. Pedosphere 2012, 22, 640–649. [Google Scholar] [CrossRef]
  58. Sherman, D.M.; Waite, T.D. Electronic Spectra of Fe3 + Oxides and Oxide Hydroxides in the near IR to near UV. Am. Mineral. 1985, 70, 8. [Google Scholar]
  59. Vašát, R.; Kodešová, R.; Klement, A.; Borůvka, L. Simple but Efficient Signal Pre-Processing in Soil Organic Carbon Spectroscopic Estimation. Geoderma 2017, 298, 46–53. [Google Scholar] [CrossRef]
  60. Buddenbaum, H.; Steffens, M. The Effects of Spectral Pretreatments on Chemometric Analyses of Soil Profiles Using Laboratory Imaging Spectroscopy. Appl. Environ. Soil Sci. 2012, 2012, e274903. [Google Scholar] [CrossRef] [Green Version]
  61. Engel, J.; Gerretzen, J.; Szymańska, E.; Jansen, J.J.; Downey, G.; Blanchet, L.; Buydens, L.M.C. Breaking with Trends in Pre-Processing? TrAC Trends Anal. Chem. 2013, 50, 96–106. [Google Scholar] [CrossRef]
  62. Hill, J.; Udelhoven, T.; Vohland, M.; Stevens, A. The Use of Laboratory Spectroscopy and Optical Remote Sensing for Estimating Soil Properties. In Precision Crop Protection—The Challenge and Use of Heterogeneity; Oerke, E.-C., Gerhards, R., Menz, G., Sikora, R.A., Eds.; Springer: Dordrecht, The Netherlands, 2010; pp. 67–85. ISBN 978-90-481-9277-9. [Google Scholar]
  63. Miller, C.E.; Williams, P.C.; Norris, K.H. Chemical principles of near infrared technology. In Near Infrared Technology in the Agricultural and Food Industries; American Association of Cereal Chemist: St. Paul, MN, USA, 2001; pp. 19–39. [Google Scholar]
  64. Martínez-Carreras, N.; Krein, A.; Udelhoven, T.; Gallart, F.; Iffly, J.F.; Hoffmann, L.; Pfister, L.; Walling, D.E. A Rapid Spectral-Reflectance-Based Fingerprinting Approach for Documenting Suspended Sediment Sources during Storm Runoff Events. J. Soils Sediments 2010, 10, 400–413. [Google Scholar] [CrossRef]
  65. Wu, Y.; Chen, J.; Ji, J.; Gong, P.; Liao, Q.; Tian, Q.; Ma, H. A Mechanism Study of Reflectance Spectroscopy for Investigating Heavy Metals in Soils. Soil Sci. Soc. Am. J. 2007, 71, 918–926. [Google Scholar] [CrossRef]
Figure 1. Range of reflectance values for the sampled soils and average spectral signatures for the PG and CP setups.
Figure 1. Range of reflectance values for the sampled soils and average spectral signatures for the PG and CP setups.
Agronomy 11 01895 g001
Figure 2. Weighted regression coefficient distribution over the spectral range obtained for PLSR models for the CP setup (coloured unbroken lines) and PG setup (coloured broken lines). Soil properties analysed for cross-validation are sand content (Sand), pH, electrical conductivity (Ec), total nitrogen content (N), extractable phosphorous (P), extractable potassium (K), extractable calcium (Ca), and extractable manganese (Mn). Black broken lines represent zero correlation, offset by 3.0 units for clarity of presentation.
Figure 2. Weighted regression coefficient distribution over the spectral range obtained for PLSR models for the CP setup (coloured unbroken lines) and PG setup (coloured broken lines). Soil properties analysed for cross-validation are sand content (Sand), pH, electrical conductivity (Ec), total nitrogen content (N), extractable phosphorous (P), extractable potassium (K), extractable calcium (Ca), and extractable manganese (Mn). Black broken lines represent zero correlation, offset by 3.0 units for clarity of presentation.
Agronomy 11 01895 g002
Figure 3. Variation in coefficients of determination for spectral subsets and preprocessing approaches. R2 values obtained for PLSR models for the CP setup (coloured unbroken lines) and PG setup (coloured broken lines). Spectral subsets were VIS (350–700 nm), NIR (701–1000 nm), SWIR (1001–2500 nm), and VIS + NIR + SWIR (350–2500 nm). Spectral preprocessing approaches were standard normal variate (SVN), detrending transformation (DT), and SVN + DT. Soil properties are sand content (Sand), pH, electrical conductivity (Ec), total nitrogen content (N), extractable phosphorous (P), extractable potassium (K), extractable calcium (Ca), and extractable manganese (Mn).
Figure 3. Variation in coefficients of determination for spectral subsets and preprocessing approaches. R2 values obtained for PLSR models for the CP setup (coloured unbroken lines) and PG setup (coloured broken lines). Spectral subsets were VIS (350–700 nm), NIR (701–1000 nm), SWIR (1001–2500 nm), and VIS + NIR + SWIR (350–2500 nm). Spectral preprocessing approaches were standard normal variate (SVN), detrending transformation (DT), and SVN + DT. Soil properties are sand content (Sand), pH, electrical conductivity (Ec), total nitrogen content (N), extractable phosphorous (P), extractable potassium (K), extractable calcium (Ca), and extractable manganese (Mn).
Agronomy 11 01895 g003
Figure 4. Correlograms of soil properties obtained from simple linear correlations for each wavelength for the CP setup (coloured unbroken lines) and the PG setup (coloured broken lines). Soil properties are sand content (Sand), pH, electrical conductivity (Ec), total nitrogen content (N), extractable phosphorous (P), extractable potassium (K), extractable calcium (Ca), extractable manganese (Mn), organic matter (OM), and extractable iron (Fe).
Figure 4. Correlograms of soil properties obtained from simple linear correlations for each wavelength for the CP setup (coloured unbroken lines) and the PG setup (coloured broken lines). Soil properties are sand content (Sand), pH, electrical conductivity (Ec), total nitrogen content (N), extractable phosphorous (P), extractable potassium (K), extractable calcium (Ca), extractable manganese (Mn), organic matter (OM), and extractable iron (Fe).
Agronomy 11 01895 g004
Table 1. Sites sampled for soils.
Table 1. Sites sampled for soils.
MunicipalityDesignation of OriginGrape CultivarLongitudeLatitudeSoil Classification
CacabelosBierzoMencía6.754 W42.626 NDystric Cambisol
CamponarayaBierzoGodello6.692 W42.606 NChromic Cambisol
Valbuena de DueroRibera de DueroTempranillo4.391 W41.631 NLithic Leptosol
MatapozuelosRuedaVerdejo4.765 W41.364 NAlbic Arenosol
Geographic coordinates refer to WGS84. The soil classification system is that of the IUSS Working Group WRB [30].
Table 2. Descriptive statistics for soil sample properties.
Table 2. Descriptive statistics for soil sample properties.
Soil PropertyNMinMaxRangeMedianMeanSDCoV
Silt (%)481456423634.5811.0031.81
Clay (%)481032221918.465.9132.02
Sand (%)482476524646.9615.1032.15
pH 485.338.473.147.637.241.1616.07
Ec (dS m−1)480.020.120.100.080.070.0349.27
Organic matter (%)480.372.402.030.811.020.5251.04
Total N (%)480.050.160.110.080.090.0333.76
P (mg kg−1)355.6558.3852.7316.3926.1220.6779.14
K (cmol kg−1)480.130.800.670.400.380.1437.55
Ca (cmol kg−1)481.5620.9019.349.8010.496.4461.37
Mn (mg kg−1)481.9028.4026.5010.6511.027.0964.35
Fe (mg kg−1)482.56212.41209.858.0028.6041.49145.11
Table 3. CV statistics for PLSR for the full VIS + NIR + SWIR range (350–2500 nm).
Table 3. CV statistics for PLSR for the full VIS + NIR + SWIR range (350–2500 nm).
PG Setup CP Setup
Soil PropertyR2RMSESERPDFactorsR2RMSESERPDFactors
Sand0.757.6787.7591.9570.708.3278.4151.796
pH0.950.3400.3433.3840.920.3290.3343.474
Ec0.890.0110.0112.7330.900.0110.0112.734
N0.680.0170.0171.7660.620.0180.0181.673
P0.906.5306.6193.0970.906.6476.7413.044
K0.650.0860.0871.6160.640.0870.0881.596
Ca0.872.3322.3572.7360.892.1412.1632.986
Mn0.624.3994.4461.5930.664.1954.2391.675
Table 4. CV statistics for PLSR for the full VIS + NIR + SWIR range (350–2500 nm) and spectral transformations.
Table 4. CV statistics for PLSR for the full VIS + NIR + SWIR range (350–2500 nm) and spectral transformations.
PG Setup CP Setup
Soil PropertyR2RMSESERPDFactorsR2RMSESERPDFactors
SVN
Sand0.767.5567.6351.9860.728.0288.1121.865
pH0.960.3170.3203.6340.920.3360.3393.423
Ec0.890.0110.0112.7320.910.0100.0112.733
N0.700.0160.0161.8850.650.0180.0181.674
P0.925.8915.9773.4350.935.7015.7843.545
K0.660.0840.0851.6550.660.0840.0851.655
Ca0.892.1192.1413.0150.911.9601.9813.255
Mn0.654.1154.1541.7150.693.9924.0321.765
DT
Sand0.747.7257.8051.9340.737.9027.9851.897
pH0.910.3560.3593.2330.920.3430.3463.353
Ec0.900.0110.0112.7330.900.0070.0112.733
N0.680.0170.0171.7640.670.0170.0171.763
P0.916.3076.3983.2040.935.7195.7953.534
K0.680.0820.0831.6940.640.0870.0871.444
Ca0.862.4122.4372.6440.902.0432.0643.125
Mn0.684.0754.1181.7240.644.3154.3601.631
SVN + DT
Sand0.777.2807.741.9560.767.5507.621.986
pH0.930.3150.323.6440.920.3420.353.353
Ec0.900.0110.012.7320.910.0110.012.733
N0.710.0160.021.8870.700.0160.021.887
P0.926.0836.173.3240.925.9005.993.424
K0.680.0830.081.6960.630.0880.091.574
Ca0.892.1982.222.9040.912.0052.033.184
Mn0.713.8873.931.8050.693.9834.031.765
Preprocessing: SVN: standard normal variate; DT: detrending transformation; SVN + DT: standard normal variate plus detrending transformation.
Table 5. Best PLSR models for predicting soil properties.
Table 5. Best PLSR models for predicting soil properties.
PG SetupCP Setup
Soil PropertyTransformationSpectral SubsetR2RMSESERPDFactorsR2RMSESERPDFactors
SandSVN + DTSWIR0.777.2767.3532.0560.777.3267.3962.046
pHNoneSWIR0.940.2840.2874.0430.940.2870.294.005
EcSVN + DTSWIR0.920.0100.013.0030.920.0100.013.003
NSVNSWIR0.770.0140.0142.1470.840.0120.0161.887
PSVNSWIR0.925.8165.8913.4850.935.5715.6533.625
KSVN + DTSWIR0.720.0760.0771.8240.710.0780.0791.774
CaSVN + DTVIS0.941.6031.9833.2550.921.8592.2772.835
MnSVN + DTVIS + NIR + SWIR0.713.8873.9281.8050.693.9834.0251.765
Transformations: SVN: standard normal variate; DT: detrending, SVN + DT: standard normal variate plus detrending. Spectral subsets: VIS (350–700 nm), SWIR (1001–2500 nm), VIS + NIR + SWIR (350–2500 nm).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rodríguez-Pérez, J.R.; Marcelo, V.; Pereira-Obaya, D.; García-Fernández, M.; Sanz-Ablanedo, E. Estimating Soil Properties and Nutrients by Visible and Infrared Diffuse Reflectance Spectroscopy to Characterize Vineyards. Agronomy 2021, 11, 1895. https://doi.org/10.3390/agronomy11101895

AMA Style

Rodríguez-Pérez JR, Marcelo V, Pereira-Obaya D, García-Fernández M, Sanz-Ablanedo E. Estimating Soil Properties and Nutrients by Visible and Infrared Diffuse Reflectance Spectroscopy to Characterize Vineyards. Agronomy. 2021; 11(10):1895. https://doi.org/10.3390/agronomy11101895

Chicago/Turabian Style

Rodríguez-Pérez, José Ramón, Víctor Marcelo, Dimas Pereira-Obaya, Marta García-Fernández, and Enoc Sanz-Ablanedo. 2021. "Estimating Soil Properties and Nutrients by Visible and Infrared Diffuse Reflectance Spectroscopy to Characterize Vineyards" Agronomy 11, no. 10: 1895. https://doi.org/10.3390/agronomy11101895

APA Style

Rodríguez-Pérez, J. R., Marcelo, V., Pereira-Obaya, D., García-Fernández, M., & Sanz-Ablanedo, E. (2021). Estimating Soil Properties and Nutrients by Visible and Infrared Diffuse Reflectance Spectroscopy to Characterize Vineyards. Agronomy, 11(10), 1895. https://doi.org/10.3390/agronomy11101895

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop