Monitoring Land Management Practices Using Vis–NIR Spectroscopy Provides Insights into Predicting Soil Organic Carbon and Limestone Levels in Agricultural Plots

.


Introduction
Agricultural soils often exhibit signs of degradation characterized by weak structure, low organic carbon content, and nutrient scarcity.In Mediterranean semi-arid environments, these conditions are exacerbated by erosion processes and dryness, especially during spring and summer.The frequent tillage practiced in land management has been identified as a contributing factor to soil degradation, ultimately hindering optimal crop production.Recognizing the urgent need for soil protection [1], avoiding intensive tillage, and implementing soil conservation measures such as cover cropping have emerged as crucial Agronomy 2024, 14, 1150 2 of 15 strategies [2,3].By preventing soil erosion and maintaining soil fertility, these measures not only enhance crop yields but also bolster the resilience of agricultural systems against climate variability and extreme weather events common in Mediterranean semi-arid environments.Moreover, they help preserve other ecosystem services provided by healthy soils, as emphasized by the Food and Agriculture Organization (FAO).These circumstances are especially concerning when considering woody crops like olive groves.
If more sustainable land management practices enhance soil conditions, an efficient and cost-effective monitoring of land conditions is the next imperative step in response to changes in land management practices.Classical methods for monitoring land status across extensive areas may be impractical due to economic and timing constraints.In recent decades, numerous authors have explored various methods for assessing land degradation, including soil erosion modeling and remote sensing techniques, particularly in Mediterranean semi-arid environments.Such research endeavors aim to identify indicators that can effectively track the impacts of soil degradation or soil recovery on agricultural land.Near-infrared reflectance spectroscopy (Vis-NIR) has been seen as a cost-effective and accurate tool for soil quality assessment [4].The accurate, widespread adoption of spectroscopy for soil monitoring can be achieved through the establishment of standard protocols for collecting laboratory soil spectra to facilitate the sharing of spectral libraries and sharing scanning existing soil archives, which would minimize the necessity for expensive sampling campaigns.Particularly, by analyzing the Vis-NIR spectra and correlating them with soil properties such as mineral composition and organic carbon content, researchers can develop predictive models and spectral indices to estimate soil attributes non-destructively and rapidly.This approach is widely used in soil science, remote sensing, and precision agriculture for monitoring soil health, assessing land degradation, and optimizing agricultural management practices [5][6][7].
Correlations between the depth of spectral features and soil organic carbon or carbonate concentration have been established [8].High soil organic carbon (SOC) levels lead to a decrease in albedo throughout the visible, shortwave infrared, and near-infrared (Vis-NIR-SWIR) reflectance spectrum [9,10].Nonetheless, it has been observed that a threshold of 2% organic carbon exists, below which the impact of SOC on soil reflectance diminishes significantly [11].
The identification of soil carbonate poses challenges due to its absorption feature, which shifts to longer or shorter wavelengths depending on the presence of impurities [12].Moreover, the depths of these spectral features are also influenced by particle size and porosity [13].
Correlations have been found between soil carbonate concentration and reflectance spectra, discerned through alterations in color and albedo [8,14].
In this research, we investigated the characteristics of agricultural areas in semiarid environments by comparing the soil conditions under different land management practices, which have an effect on the spectral signature of soils.Specifically, we focused on traditional soil-intensive tillage methods, and we compared them to similar soils used for cropping but managed in an alternative way.We examined how various soil characteristics, including soil organic carbon content, bulk density, and mineral composition, affect the visible and infrared response observed in the Vis-NIR spectra.By doing so, we aim to increase knowledge and enrich existing databases for monitoring and tracking alterations in soil compositions over both temporal and spatial scales to enhance the precision of soil data for local and regional inventories and also encompass an examination of the divergent impacts of soil organic carbon and calcium carbonate on soil spectra.This investigation aims to elucidate the nuanced interactions between these variables and their influence on spectral signatures.

Materials and Methods
The study area is located in central Spain (municipalities: Chinchón, Valdelaguna, Perales de Tajuña, Belmonte del Tajo, Arganda del Rey, Torres de la Alameda, and Villarejo de Salvanés).The area has a semi-arid Mediterranean climate, with average annual temperatures of 14.5 • C, 473 mm annual accumulated precipitation, and 1112 mm potential evapotranspiration.The altitude of sampling points was between 550 m and 560 masl.The aforementioned municipalities are dedicated to agricultural use, which is also one of the most important land uses in the country.Of Spain's total area of 50.59 million hectares, agriculture constitutes a significant portion.Of this, 32% is agricultural land, and 5.5%, or 2.77 million hectares, is dedicated specifically to olive orchards, which is the focus of this paper.The remaining areas are classified as forest or shrublands, grazing lands, urban or industrial areas, and inland waters, totaling 33.77 million hectares.
According to the FAO soil classification map for the study region, which pertains to the Community of Madrid [15], Cambisols occupy approximately 34% of the area, Luvisols around 21%, and Leptosols about 19%.
Sixteen crop fields, eight pairs of fields, were selected to study and compare two different soil management types in each field (Table 1).On the one hand, eight of them have management A, or alternative to tillage management, involving typically no-till or vegetative cover.These plots had soil covered by spontaneous vegetation that comprised a diverse array of plant species.This includes grasses such as Avena spp., Lolium spp., and Poa spp., as well as leguminous plants like Trifolium spp.Additionally, it encompasses annual broadleaf weeds such as Papaver rhoeas, Sinapis arvensis, Amaranthus spp., and Diplotaxis spp., along with perennial weeds like Carduus nutans.Other common species found in spontaneous vegetation are Bromus spp.and Chenopodium spp.On the other hand, the rest eight plots with management B, or conventional tillage, were used as a reference for business-as-usual management (Figure 1).Traditional tillage consists of labor passes usually three to four times a year to a depth of 15-20 cm.Samples were gathered from a single geographical area to evaluate the prediction of soil properties across a naturally diverse range, aiming to offer predictions relevant to regional image analysis.
For each sampling point, four composite-samples of soil were obtained randomly to collect 1 kg at four established depths: 0 to 5 cm, 5 to 10 cm, 10 to 20 cm, and, finally, 20 to 30 cm.This resulted in a total of 64 soil samples (8 × 8 × 4).These samples were air-dried and sieved (2 mm mesh) to remove coarser fractions.
were used as a reference for business-as-usual management (Figure 1).Traditional tillage consists of labor passes usually three to four times a year to a depth of 15-20 cm.Samples were gathered from a single geographical area to evaluate the prediction of soil properties across a naturally diverse range, aiming to offer predictions relevant to regional image analysis.For each sampling point, four composite-samples of soil were obtained randomly to collect 1 kg at four established depths: 0 to 5 cm, 5 to 10 cm, 10 to 20 cm, and, finally, 20 to 30 cm.This resulted in a total of 64 soil samples (8 × 8 × 4).These samples were air-dried and sieved (2 mm mesh) to remove coarser fractions.Texture was determined through particle size analyses utilizing the hydrometer method [16].
The soil organic matter (SOM) content was measured using the Walkley-Black method, in which the soil undergoes an oxidation reaction with a standardized potassium dichromate solution.For these measurements, soil samples were previously milled [17].The bulk density (g/cm 3 ) was determined at each sampling point using the cylinder method [18].Differences in SOM were measured in different soils and depths.
The mineral composition of samples was estimated by XR diffraction using SIEMENS D-5000 equipment with a scanning speed of 1 • 2θ/min and Cu-Kα radiation (40 kV, 20 mA).XR analysis was carried out on randomly oriented samples (bulk samples).Powdered whole-rock samples were scanned from 5 • to 31 • 2θ.The mineral intensity factors (MIFs) method was applied to XRD peak intensity ratios normalized to 100% with calibration constants for the quantitative estimation of mineral contents [19][20][21].
The air-dried soil sieved samples' reflectance was determined using an Analytical Spectral Device ASD FieldSpecPro Vis-NIR spectroradiometer in the laboratory, measuring wavelengths from 350 nm to 1100 nm with a spectral resolution of 3 nm.Following the recommendations of the Global Soil Laboratory Network Initiative on Soil spectroscopy (GLOSOLAN-Spec) of the Global Soil Partnership [22], the samples used for this process were sieved and air-dried but not crushed, as the sample reflectance is influenced by their size and would increase the reflectance by brightening the samples, making them appear lighter than in their initial conditions.Spectroradiometry allows for the recording of spectral signatures of various surfaces with high spectral resolution [23].The scanning of the samples was conducted in a dark room, connecting it via fiber optics to an ASD contact probe along with a halogen bulb at a temperature of 2900 K.The setup was positioned vertically on a tripod, with the probe placed 5 cm above the sample, aiming to minimize lighting differences.Black plastic plates with a diameter of 5 cm were filled with soil samples with 1 cm thickness of soil.An average of 10 internal scans were performed to reduce random noise in the entire spectrum and obtain accurate results; noise was present at extremes of the wavelength range; noisy segments of spectra were eliminated prior to analysis, resulting in a final range from 400 nm to 900 nm.A spectrum on the standard panel was used to obtain the white reference.The spectrometer was calibrated every 10 min.
Statistical analyses were performed to find significant differences between different agricultural soil management practices regarding soil organic carbon (SOC) concentration, estimated as 58% of SOM [24], bulk density and mineral composition in the different soil layers, were performed using the IBM ® SPSS ® Statistics 26 software.A factorial analysis of variance (ANOVA) with repeated measures was conducted to identify significant differences between variables.To understand the interaction between depth and management while considering the locations, the locations were grouped according to their soil type to increase the robustness of the ANOVA.The normal distribution of variables was assessed using the Shapiro-Wilk test [25] (see supplemental data).
Quantitative analysis of reflectance spectra was conducted using Partial Least Squares Regression (PLSR) with ParLeS version 3.1 software [26].This procedure addressed collinearity issues arising from numerous highly correlated variables and a limited number of samples.The matrix of 66 samples was randomly split into two sets; the first one consisted of 2/3 of the samples to perform model calibration, and the second group, which consisted of the remaining 1/3 of the samples, was used to validate the model.Data conforming to the soil diffuse reflectance spectra were transformed previously as preprocessing to produce more robust models being standardized to sum 1, and they were also mean-centered prior to the multivariate modeling.
Subsequently, cross-validation, which is intended to assess whether a statistical model can be generalized to an independent dataset, was performed to determine the optimal factors for the regression model [27].The validation accuracy was assessed using the root mean square error (RMSE; Equation ( 1)).
where, ŷi is the predicted value, and y i is the observed value.
To estimate the predictive capability of the model, the value of R 2 should fall within the range (0 < R 2 < 1), with 1 indicating a perfect prediction or goodness of fit [28].The residual predictive deviation (RPD) represents the ratio between the standard deviation of the Y data, denoted as SD, and the RMSE of cross-validated predictions.RPD is used to determine the quality of predictions: the consensus holds that values below 1.5 signify a deficient predictive model when applied to soil property prediction, while values between 1.5 and 2.0 are deemed acceptable, and those exceeding 2.0 are considered satisfactory [29,30].

Results
Soils under study exhibit alkaline characteristics, with a high content of calcium carbonate, which exceeds 20% (Table 2).The pH for all the sets of locations ranged from 7.3 to 8.4; low electrical conductivity was found, ranging from 0.3 to 1 dS/m.From a classification standpoint according to the WRB [31] 2 aims to characterize the soil from 0 to 30 cm depth as a whole.

Soil Organic Matter Comparing Tillage and Alternative Management
There is a wider variation in the median values of SOM in soils under alternative management; this is because alternative management does not experience the homogenization that tillage provides to the soil properties, considering the soil depth (Figure 2).These agricultural soils are situated within an area spanning approximately 60 square kilometers; the highest values of SOM are around 10 g/kg and the lowest around 4 g/kg.There were significant differences (p < 0.05) between traditional tillage (SOM 6.6 ± 1.4 g/kg) and alternative managements (7.2 ± 2.3 g/kg) considering all the different layers (0 to 30 cm depth).The most vulnerable soils, classified as Leptosols, showed a lower concentration for SOC, below average, for both managements.

Soil Organic Matter Comparing Tillage and Alternative Management
There is a wider variation in the median values of SOM in soils under alternative management; this is because alternative management does not experience the homogenization that tillage provides to the soil properties, considering the soil depth (Figure 2).These agricultural soils are situated within an area spanning approximately 60 square kilometers; the highest values of SOM are around 10 g/kg and the lowest around 4 g/kg.There were significant differences (p < 0.05) between traditional tillage (SOM 6.6 ± 1.4 g/kg) and alternative managements (7.2 ± 2.3 g/kg) considering all the different layers (0 to 30 cm depth).The most vulnerable soils, classified as Leptosols, showed a lower concentration for SOC, below average, for both managements.The general pattern of SOM distribution throughout the soil profile is shown in Figure 3, where the average of soil organic carbon in the different layers changes over depth.In this case, as expected, there are no differences between managements in the deepest layers, but there is a difference in the upper layer.Moreover, tillage leads to similar SOM content, which is not statistically different, ranging from 7.1% in the topsoil to 5.9% at 30 cm; this was caused by the abovementioned homogenization of soil materials due to plow.The general pattern of SOM distribution throughout the soil profile is shown in Figure 3, where the average of soil organic carbon in the different layers changes over depth.In this case, as expected, there are no differences between managements in the deepest layers, but there is a difference in the upper layer.Moreover, tillage leads to similar SOM content, which is not statistically different, ranging from 7.1% in the topsoil to 5.9% at 30 cm; this was caused by the abovementioned homogenization of soil materials due to plow.The factorial analysis of variance (Table 3) shows that the three main variables considered alone-management, depth, and type of soil-have a significant effect on SOM.Both management and type of soil are also significant.However, other combinations of variables cannot explain the variance in SOM; however, in Figure 3, we can observe that Calcisols and Cambisols experience an increase in SOM, and there is a tendency to increase SOM up to 30 cm depth in both soil types.The factorial analysis of variance (Table 3) shows that the three main variables considered alone-management, depth, and type of soil-have a significant effect on SOM.Both management and type of soil are also significant.However, other combinations of variables cannot explain the variance in SOM; however, in Figure 3, we can observe that Calcisols and Cambisols experience an increase in SOM, and there is a tendency to increase SOM up to 30 cm depth in both soil types.Taking into account only SOM and limestone, we can find several differences in their relationship in the different soils of study.In Figure 4, we observe that, except for Cambisols, where there is a significant negative correlation between SOM content and calcium carbonate content, in the rest of the samples, this relationship is not clear from different soils.This fact can lead to false predictions of SOM content masked by the reflectance of calcium carbonate.
Agronomy 2024, 14, x FOR PEER REVIEW 8 of 15 Taking into account only SOM and limestone, we can find several differences in their relationship in the different soils of study.In Figure 4, we observe that, except for Cambisols, where there is a significant negative correlation between SOM content and calcium carbonate content, in the rest of the samples, this relationship is not clear from different soils.This fact can lead to false predictions of SOM content masked by the reflectance of calcium carbonate.

Prediction of Soil Characteristics Based on Spectral Data
The variations observed in the SOM for the different managements can also be noticed in the variations of the Vis-NIR spectra (400-900 nm) of samples (Figure 5).The colors and, therefore, the spectral features of soils are related to the magnitude of organic matter and minerals present in soils [32].

Prediction of Soil Characteristics Based on Spectral Data
The variations observed in the SOM for the different managements can also be noticed in the variations of the Vis-NIR spectra (400-900 nm) of samples (Figure 5).The colors and, therefore, the spectral features of soils are related to the magnitude of organic matter and minerals present in soils [32].
In Figure 5, we can observe the different patterns of soil reflectance caused by management in which each soil shows its unique spectral signature; on the right position at each site, tilled soils usually exhibit narrow spectra, indicating similar reflectance for soils coming from different depths.When there is a difference between spectra, soils coming from the deepest layers, shown in grey or orange line colors, usually show higher reflectance.On the left columns, in the charts corresponding to the alternative managements (A), different depths of soils exhibited higher scattering than the corresponding similar soil managed by tillage.We hypothesize that these patterns indicate less SOC, high content of lime, or both.At the same time, the absence of the disturbance of soil intensifies the different reflectance of soil layers.In blue, the samples are in the topsoil (0-5 cm depth); in orange, the layers are 5-10 cm, in gr to 20 cm, and in yellow, 20 to 30 cm.The samples corresponding to the alternative manageme labeled with the letter A (on the left), while those associated with traditional tillage end with the right).Additionally, each sample is presented with scanned soil colors in the Vis range ( with a black frame), color code according to the Munsell scale and the average (in brackets) fo organic carbon (%).In blue, the samples are in the topsoil (0-5 cm depth); in orange, the layers are 5-10 cm, in grey, 10 to 20 cm, and in yellow, 20 to 30 cm.The samples corresponding to the alternative management are labeled with the letter A (on the left), while those associated with traditional tillage end with B (on the right).Additionally, each sample is presented with scanned soil colors in the Vis range (boxes with a black frame), color code according to the Munsell scale and the average (in brackets) for soil organic carbon (%).
The increase of SOC content or decrease of lime content are associated with a significant decrease in the reflectance, usually ranging over the entire wavelength, from 400 to 900 nm, in the study.This is based on the effects of the dark soil produced by the increase of SOC [33].
Exploratory analyses were conducted on all samples to identify outliers before constructing the PLRS model [34]; eliminating outliers can enhance prediction accuracy; therefore, they were excluded [35].The number of outliers eliminated for each variable is provided in Table 4 (N out).The PLSR models were able to explain 85%, 87%, 79%, and 75% of the variance in SOC, limestone, clay, and quartz, respectively.The variations of the rest of the parameters did not contribute to the explanation of the variance within the dataset.For the PLSR calibration set, R 2 values ranged between 0.14 and 0.79.The best predictive models were obtained for SOC (R 2 = 0.79, RMSE = 0.09, RPD = 2.05), and limestone (R 2 = 0.70, RMSE = 11.1,RPD = 1.75), and an acceptable prediction was found for quartz (R 2 = 0.56, RMSE = 10.2,RPD = 1.5).The rest of the variables would yield inaccurate predictions due to their low R 2 , high RMSE, or RPD lower than 1.5.The variability of the sample set, like the colors shown in Figure 5, may affect the accuracy of PLSR calibration models developed for the rest of the minerals: quartz, phyllosilicate, feldspar, and plagioclase, which in this study do not provide significant information.
In PLSR, distinct wavelengths were identified as highly responsive to changes in the concentration of soil organic carbon (SOC) within the range of 850-860 nm.Additionally, separate wavelengths were found to be sensitive to variations in the concentration of lime (720-730 nm).
In order to establish the accuracy of the prediction models of SOC, limestone, or the rest of the variables considered in this study, with 66% of samples, the model was subjected to cross-validation.Samples for the validations were not used in the model.Figure 6 shows the correlations obtained between the values of variables predicted by the Vis-NR models and the corresponding values observed in the laboratory in the set of the rest of the 33% of samples used for the validation of models, which were not included in the set for calibration.
rest of the variables considered in this study, with 66% of samples, the model was subjected to cross-validation.Samples for the validations were not used in the model.Figure 6 shows the correlations obtained between the values of variables predicted by the Vis-NR models and the corresponding values observed in the laboratory in the set of the rest of the 33% of samples used for the validation of models, which were not included in the set for calibration.

Discussion
The different patterns of soil reflectance and unique spectral signature over depth caused by management indicate the vertical distribution of soil organic carbon or limestone and highlight the importance of incorporating vertical dimension into models of SOC dynamics, as demonstrated by other soil research on the distribution of SOC over depth [25,26].These soils, which often exhibit lighter-colored subsurface layers due to higher concentrations of calcium carbonate, are particularly susceptible to erosion induced by tillage practices.As a consequence of erosion, subsurface layers are frequently exposed to the surface.One of the primary objectives of this study is to establish soil characteristics, based on their reflectance properties, as indicators of soil erosion.Reflectance

Discussion
The different patterns of soil reflectance and unique spectral signature over depth caused by management indicate the vertical distribution of soil organic carbon or limestone and highlight the importance of incorporating vertical dimension into models of SOC dynamics, as demonstrated by other soil research on the distribution of SOC over depth [25,26].These soils, which often exhibit lighter-colored subsurface layers due to higher concentrations of calcium carbonate, are particularly susceptible to erosion induced by tillage practices.As a consequence of erosion, subsurface layers are frequently exposed to the surface.One of the primary objectives of this study is to establish soil characteristics, based on their reflectance properties, as indicators of soil erosion.Reflectance measurements are conducted across different layers to capture variations in soil color, as deeper layers tend to be lighter in color compared to uppermost layers.
Soil management alternatives to tillage have demonstrated an increase in SOC, confirming that management practices can significantly improve soil characteristics, even in the least developed soils like Leptosols.
Understanding the benefits of sustainable land management practices extends beyond merely assessing the concentration of SOM; evaluating its quality and stability is paramount.Previous research conducted in this region on cover crops in vineyards revealed contrasting characteristics between soils managed with tillage and those managed with cover crops.Soils subjected to tillage exhibited little but well-humified SOM, with a higher proportion of large, complex, and more humified organic molecules.Conversely, soils managed with cover crops showed a greater presence of fresh SOM [36].While this pool of fresh SOM may not directly contribute to carbon sequestration, it plays a crucial role in enhancing soil aggregation, water infiltration, water holding capacity, erosion reduction, and the promotion of root and microbial activity.Consequently, the next step involves a further examination of the nature of SOM in subsequent projects.
As has been frequently reported in the scientific literature, visible-near-infrared reflectance spectra used with soil samples can be used to predict certain soil properties.In our study, we observed that prediction of soil organic carbon (SOC) levels can be influenced by the presence of calcium carbonate in the soil.This is because calcium carbonate has the tendency to increase the reflectance of soils, regardless of the actual SOC content.As a result, even in cases where the SOC levels are high, the presence of calcium carbonate can artificially elevate the reflectance readings, leading to potentially misleading predictions.This phenomenon poses a challenge in accurately assessing SOC levels using reflectance spectroscopy, as the influence of calcium carbonate must be carefully accounted for to avoid misinterpretation of results.Therefore, it is imperative to perform complementary chemical analyses to directly measure SOC and calcium carbonate content or any relevant minerals, which can then validate and refine our spectral predictions to develop robust calibration models that effectively account for the impact of calcium carbonate on soil reflectance in order to improve the accuracy of SOC predictions.Our study demonstrates the feasibility of predicting soil organic carbon, lime content, and, less accurately, quartz content using reflectance data obtained with a high-resolution laboratory spectrometer.It is important to highlight that the recalibration of Partial Least Squares (PLS) predictive functions would be necessary for other soil types and mineral compositions.
Soil organic carbon and calcium carbonate contents exhibited the most accurate predictions, followed by quartz.The validation R 2 for these analyses exceeded 0.5, and the RPD (Ratio of Performance to Deviation) was satisfactory.We found the effectiveness of specific spectral regions within the 400-900 nm spectrum for predicting carbon content (at 850-860 nm), which is consistent with the spectral regions contributing to the prediction of SOC in other regions [37] and carbonate (720-730 nm), although other authors found absorption bands by 2250 to 2465 [38].This underscores the potential of this methodology as a tool for swift and dependable soil mapping.However, the presence of both soil organic carbon and calcium carbonate can lead to overlapping spectral signatures, making it challenging to accurately differentiate between the two components based solely on remote sensing data.This can result in spectral confounding, where the contributions of each constituent to the overall reflectance spectrum are difficult to discern.Soils with complex compositions may exhibit spectral mixing effects, where the spectral signature observed remotely is a combination of the spectral characteristics of multiple soil constituents.This can complicate the interpretation of remote sensing data and make it challenging to isolate the contributions of individual soil properties or monitor soil health.
In this research, soils with higher levels of organic carbon often exhibit reduced reflectance values in their spectra.The prediction accuracy for SOC has been found to be high in different studies, with R 2 by 0.9 and RPD higher than 3 in the scientific literature.Samples coming from different depths are arguably increasing the heterogeneity of the set of data [39].This type of study, including samples from different depths, shows moderate predictions, like in those found by Morgan et al. [40], R 2 was 0.82 and RPD = 1.5, and in Brown et al. [41], R 2 ranged from 0.68 to 0.86 and RPD from 1.73 to 2.62.
Calcium carbonate is a highly reflective mineral, particularly in the near-infrared region of the electromagnetic spectrum.As a result, soils with elevated calcium carbonate content often display increased reflectance in their spectra.The prediction of calcium carbonate content on soils by PLSR is not so successful: e.g., Summers et [42], Turkey [43], and Iran [44], respectively.In this study, the R 2 was 0.70 for Calcisols in semi-arid environments in central Spain.

Figure 1 .
Figure 1.Example of the election of sampling points among plots having close locations, so similar soil but different land management.(a) sample O-12, (b) sample O-16.A = Alternative to Tillage; B = Traditional tillage.

Figure 1 .
Figure 1.Example of the election of sampling points among plots having close locations, so similar soil but different land management.(a) sample O-12, (b) sample O-16.A = Alternative to Tillage; B = Traditional tillage.

Figure 2 .
Figure 2. Soil organic matter (SOM) in the agricultural plots considered in the study.Values are considering the whole profile from the surface to 30 cm depth.

Figure 2 .
Figure 2. Soil organic matter (SOM) in the agricultural plots considered in the study.Values are considering the whole profile from the surface to 30 cm depth.

Figure 3 .
Figure 3. Soil organic matter (SOM, g/kg) change over depth in alternative (A) and traditional tillage management (B).The effect of management is shown over depth (upper left chart), between soil types (upper right chart), and considering both depth and soil type (down chart).Vertical lines denote 0.95 confidence intervals.

Figure 3 .
Figure 3. Soil organic matter (SOM, g/kg) change over depth in alternative (A) and traditional tillage management (B).The effect of management is shown over depth (upper left chart), between soil types (upper right chart), and considering both depth and soil type (down chart).Vertical lines denote 0.95 confidence intervals.

Figure 4 .
Figure 4. Correlations between soil organic matter content (SOM g/kg) and calcium carbonate (Limestone g/g) in the different groups of soils.

Figure 4 .
Figure 4. Correlations between soil organic matter content (SOM g/kg) and calcium carbonate (Limestone g/g) in the different groups of soils.

Figure 5 .
Figure 5. Vis-NIR spectra of soil samples (from O-3 to O-24 agricultural plots) at different deIn blue, the samples are in the topsoil (0-5 cm depth); in orange, the layers are 5-10 cm, in gr to 20 cm, and in yellow, 20 to 30 cm.The samples corresponding to the alternative manageme labeled with the letter A (on the left), while those associated with traditional tillage end with the right).Additionally, each sample is presented with scanned soil colors in the Vis range ( with a black frame), color code according to the Munsell scale and the average (in brackets) fo organic carbon (%).

Figure 5 .
Figure 5. Vis-NIR spectra of soil samples (from O-3 to O-24 agricultural plots) at different depths.In blue, the samples are in the topsoil (0-5 cm depth); in orange, the layers are 5-10 cm, in grey, 10 to 20 cm, and in yellow, 20 to 30 cm.The samples corresponding to the alternative management are labeled with the letter A (on the left), while those associated with traditional tillage end with B (on the right).Additionally, each sample is presented with scanned soil colors in the Vis range (boxes with a black frame), color code according to the Munsell scale and the average (in brackets) for soil organic carbon (%).

Figure 6 .
Figure 6.Comparison between measured and cross-validation predicted values of soil contents in soil organic carbon, limestone, quartz, clay, phyllosilicates, and K-feldspars.

Figure 6 .
Figure 6.Comparison between measured and cross-validation predicted values of soil contents in soil organic carbon, limestone, quartz, clay, phyllosilicates, and K-feldspars.

Table 1 .
Location and management of sampling points.

Table 3 .
Effects of management (tillage or alternative), depth, and soil type (FAO, WRB) on soil organic matter.Significant differences (p < 0.05) are marked with *.

Table 3 .
Effects of management (tillage or alternative), depth, and soil type (FAO, WRB) on soil organic matter.Significant differences (p < 0.05) are marked with *.

Table 4 .
Performance of Vis-NIR spectroscopy-based predictive models for the soil variables, including statistical parameters obtained for calibration and cross-validation.Var Y is the variance explained.
al., Volkan Bilgili et al., and Khayamim et al. reported R 2 values of 0.69, 0.64, and 0.54 in soils of dry environments in Australia