Next Article in Journal
Sorption Processes in Soils and Sediments
Previous Article in Journal
Short-Term Waterlogging Depresses Early Growth of Sunflower (Helianthus annuus L.) on Saline Soils with a Shallow Water Table in the Coastal Zone of Bangladesh
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ensemble Modeling on Near-Infrared Spectra as Rapid Tool for Assessment of Soil Health Indicators for Sustainable Food Production Systems

1
CGIAR Research Program on Climate Change, Agriculture and Food Security (CCAFS) East Africa, International Livestock Research Institute (ILRI), Nairobi P.O. Box 30709-00100, Kenya
2
Department of Chemistry, School of Pure and Applied Sciences, Kisii University, Kisii P.O. Box 408-40209, Kenya
3
World Agroforestry (ICRAF), United Nations Avenue, Nairobi P.O. Box 30677-00100, Kenya
*
Author to whom correspondence should be addressed.
Soil Syst. 2021, 5(4), 69; https://doi.org/10.3390/soilsystems5040069
Submission received: 7 October 2021 / Revised: 3 November 2021 / Accepted: 6 November 2021 / Published: 12 November 2021

Abstract

:
A novel total ensemble (TE) algorithm was developed and compared with random forest optimization (RFO), gradient boosted machines (GBM), partial least squares (PLS), Cubist and Bayesian additive regression tree (BART) algorithms to predict numerous soil health indicators in soils with diverse climate-smart land uses at different soil depths. The study investigated how land-use practices affect several soil health indicators. Good predictions using the ensemble method were obtained for total carbon (R2 = 0.87; RMSE = 0.39; RPIQ = 1.36 and RPD = 1.51), total nitrogen (R2 = 0.82; RMSE = 0.03; RPIQ = 2.00 and RPD = 1.60), and exchangeable bases, m3. Cu, m3. Fe, m3. B, m3. Mn, exchangeable Na, Ca (R2 > 0.70). The performances of algorithms were in order of TE > Cubist > BART > PLS > GBM > RFO. Soil properties differed significantly among land uses and between soil depths. In Kenya, however, soil pH was not significant, except at depths of 45–100 cm, while the Fe levels in Tanzanian grassland were significantly high at all depths. Ugandan agroforestry had a substantially high concentration of ExCa at 0–15 cm. The total ensemble method showed better predictions as compared to other algorithms. Climate-smart land-use practices to preserve soil quality can be adopted for sustainable food production systems.

1. Introduction

In sub-Saharan Africa, 62% of the rural population depends on agriculture as the main source of their livelihoods, and hence there is a close link between agriculture and soil health [1]. Agricultural sustainability necessitates a good understanding of soil characteristics which can inform farmers in making farming decisions and improve the practices that enhance soil quality [1,2]. Both the physical and chemical properties of soil have been used extensively to monitor soil health characteristics [3,4]; while these properties are important for farm productivity, they vary within fields and with land-use types [2,5]. If these soil properties are well-characterized, they should serve as indicators of soil health and be easy to measure using standardized methods [2]. The measurement of these soil health indicators faces significant technological difficulties due to the large number of properties involved [6]. Convectional analytical techniques such as wet chemical analysis have always been used for this purpose; however, these wet methods are time-consuming and expensive, prompting a need for a robust alternative method. Several authors have suggested near-infrared reflectance spectroscopy (12,500–4000 cm−1; 800–2500 nm) as an alternative technique to wet chemical analysis [6,7,8,9]. Near-infrared absorption bands are overtones and combinations of fundamental vibrations of XH bonding, where X can be carbon, nitrogen, oxygen, or sulfur [10]. Near-infrared spectroscopy has the advantage of being rapid, non-destructive, inexpensive, precise, and can be used to estimate water-bearing minerals, such as clay minerals and organic matter, carbon and nitrogen, and cation exchange capacity [3], as well as micro-nutrients and exchangeable cations in soil samples [1,7,11]. Additionally, the technique has been applied in precision soil management as well as regular soil analysis [12]. Soriano-Disla et al. [8] reviewed soil spectroscopic models and published and listed several soil properties that could be determined by near-infrared spectroscopy; these properties include soil water content, clay, sand, soil organic carbon (SOC), CEC, exchangeable Ca and Mg, total N and pH. These spectroscopic models used different spectral preprocessing techniques such as wavelength range selection, the scatter correction method, mean normalization, baseline offset, and derivatives [9,13,14] to increase the robustness and predictability of the models. Additionally, modeling the relationship between near-infrared spectra with soil properties requires several multivariate procedures such as principal components regression, partial least squares regression (PLSR), stepwise multiple linear regression (SMLR), Fourier regression, locally weighted regression (LWR), and artificial neural networks. None of these multivariate procedures have gained widespread adoption since a model that works well for one application may be unsuitable for another. The search for an optimum algorithm for a specific NIR-based application is difficult since no single algorithm always performs well in any domain [15]. Fortunately, the growing ‘ensemble’ concept has prompted a fundamental shift in people’s thinking [16]; rather than attempting to construct a single superb model, numerous simple models are used in tandem [15].
Ensemble modeling uses many PLS models with and without spectral preprocessing for prediction and combines several prediction models in order to improve the accuracy of weak models. Further, creating an ensemble entails two steps: (1) creating various models and (2) combining their estimates [17]. Ensembles are constructed using techniques such as Bayesian model averaging, boosting [18] and bagging [19]. The Bayesian model averages estimates from different models, weighted by their posterior evidence, while bagging bootstraps the training dataset and averages the estimates. Boosting builds models iteratively by varying case weights and employs the weighted sum of the sequence of model estimates [17]. The total ensemble algorithm technique has been little studied in situations that are common for calibration and prediction in chemistry [20]. This technique has gained increasing attention for the multivariate calibration of NIR spectra, by combining the results of multiple individual models to produce a single prediction [21]. The output of the total ensemble algorithm is computed by averaging the predicted values computed by its constituent learners [22]. The method’s key assumption is that multiple models will detect and encode more features of the relationship between independent and dependent variables than a single model [23]. To obtain a good ensemble, it is generally believed that the member models should be as accurate and diverse as possible [24].
In this study, a new approach involving total ensemble modeling of NIR spectra was used to predict numerous soil health indicators from diverse climate-smart land uses at different soil depths. We compared the prediction accuracy and performance of the total ensemble method with other five machine learning algorithms: random forest optimization (RFO), gradient boosted machines (GBM), partial least squares (PLS), Bayesian additive regression trees (BART) and Cubist. Further, we used the predicted dataset to assess how land-use practices impact on selected soil health indicators.

2. Materials and Methods

2.1. Soil Samples and Types

A total of 315 samples were collected using soil coring auger at three depths, 0–15 cm, 15–45 cm, and 45–100 cm, from six climate-smart land-use types: agroforestry, community forest, cropland with soil and water conservation (SWC), crop land without SWC, grassland and control. These land-use types are in East African Climate-Smart Villages (CSV) in Lushoto (Tanzania), Hoima (Uganda), Nyando and Wote (Kenya). Collected soil samples were oven-dried at 105 °C and then finely ground to powder and passed through a 2 mm sieve. The soil types in Lushoto are Regosols, Lithic Leptosols, Cutanic Acrisols and Ferralic Cambisols [25], while in Hoima the soil types are Vertisols, which have 30% or more clay [26]. Nyando and Wote soil types are both fluviatile and lacustrine in origin and vary from colluviums to alluvium and lacustrine clays [27].

2.2. NIR Spectroscopy and Reference Laboratory Analysis

The finely ground soil samples were loaded into a glass Petri dish and then scanned in reflectance mode using a high-intensity contact probe attached to Fourier Transform Infrared Multi-Purpose Analyzer (FTIR MPA) from 8000 to 4000 cm−1, 1250–2500 nm, at World Agroforestry Centre (ICRAF) in Nairobi, Kenya. For method validation, all the samples scanned using the NIR spectroscopy technique were taken for conventional methods for wet chemistry analysis. Selected physical and chemical properties such as pH, extractable P, K, Ca, Mg, Na, Mn, Fe, Cu, Zn, B, Mo, S, and Al, exchangeable bases (ExBas) (sum of Mehlich exch Ca, Mg, K, Na), exchangeable acidity (ExAc) and electrical conductivity (Ecd), as well as total N and C, were analyzed.

2.3. NIR Ensemble Modeling Using Spectroscopic Data

Ensemble modeling, unlike traditional single modeling methods, establishes many “weak” models then aggregates the predict results of each “weak” model through weighted average methods. In this study, we used regression modeling of the total ensemble algorithm plus other five machine learning algorithms, random forest optimization (RFO), gradient boosted machines (GBM), partial least squares (PLS), Bayesian additive regression trees (BART) and Cubist, to model the physical and chemical characteristics of soil. Using these modeling techniques, the processed spectral data were linked to laboratory-measured soil property data. RFO trains each tree separately using random sampling of the data, while GBM is a hybrid method that incorporates both boosting and bagging approaches [28]. BART is a nonparametric Bayesian regression method that employs dimensionally adaptable random basis elements to make inferences and estimate an unknown regression function [29], while the Cubist model includes boosting with training committees (typically more than one), which is comparable to the approach of “boosting” by generating a sequence of trees with successively adjusted weights [30].

2.4. Model Validation

To evaluate the accuracy of models, the coefficient of determination (R2), root mean square error (RMSE), the ratio of performance to deviation RPD and ratio of performance to interquartile distance (RPIQ) were used. Calibration models having an R2 > 0.91 are considered to be excellent, those with an R2 between 0.82 and 0.90 are good, while an R2 between 0.66 and 0.81 indicates satisfactory predictions [8]. RPD was calculated as the fraction of the standard deviation (SD) and the RMSEP (RPD = SD/RMSEP) [31], while RPIQ was calculated as a fraction of the interquartile range of the data (Q3–Q1) and the RMSEP (RPIQ = IQR/RMSEP) [32,33]. RPIQ > 1.03 indicates good predictions, 0.77–1.03 indicates reasonable prediction and <0.77 indicates non-reliable predictions. RPIQ is inversely related to R2, and so was used in isolation to rank prediction performance. Larger values of RPIQ and smaller RMSE indicate better model performance [32].

2.5. Statistical Analysis

The reference ensemble method available in the R system for statistical computing version 3.1.0 via the “caretEnsemble” add-on package [34] was used as a modeling tool. Principal component analysis (PCA) was used to visualize the variability of soil spectral signatures in the whole dataset and to identify properties explaining the greatest variability in order to select the best indicators affecting soil health. Significant differences between land-use practices were tested by analysis of variance (ANOVA). Tukey’s honest significance difference (HSD) tested the mean separation when analysis showed statistically significant differences (p < 0.05).

3. Results

3.1. Soil Properties across the Study Sites and Spectral Datasets

The laboratory-measured soil properties across the six land-use types as used for calibration are presented in Table 1. pH values range from slightly acidic to basic (4.43 to 9.08), while the average pH was 6.44 ± 1.34; this average suggests that the fields do not need lime to bring the pH to favorable levels. The total nitrogen and total carbon contents varied from 0.02% to 0.83% and 0.30–11.86%, respectively. The coefficients of variation (CV) were greater than 10% for all the elements analyzed (Table 1).
The PCA similarity maps of PC1 and PC2 made using the raw spectra show a clustering of soil samples among the spectra from various sampling counties in Kenya, Uganda and Tanzania (Figure 1). Further, there is a clear overlap of various spectra from sampling sites in Uganda. The clustering of soil samples could explain the effect of geographical origin of soil samples on NIR spectroscopy analysis when assessing both physical and chemical soil properties. Soil samples from the Mbuzii area in Tanzania were grouped away from Yamba soils samples using PC1, indicating that the largest contribution to variance in soil characteristics from Tanzania could be from these two regions from which the soil samples were taken. Although the Kenyan soil samples could not be discriminated using PC1, the samples were grouped separately according to the sampling site in PC2. There was no distinguishable variation profile amongst soil samples from Uganda to discriminate different sampling sites (Figure 1). These findings show that the primary sources of spectral variance are differences in soil types and that PCA may be used to distinguish sites but not land-use type.

3.2. Exploratory Analysis of Soils Near-Infrared Spectra

The spectra of all soils were similar in shape, with strong absorbance’s around 5200 cm−1 and three distinct absorption peaks (at roughly 7060 cm−1, 5200 cm−1, and 4520 cm−1) (Figure 2). The NIR raw spectral signature of scanned soil samples ranged from 8000 cm−1 to 4000 cm−1 regions. The NIR spectra reveal details of diverse soil constituents, including organic and inorganic components [2]. The band near 7060 cm−1 is caused by the overlapping of the first overtone of the O-H stretching vibration associated with water and hydroxyl groups in clay [2]. The band near 5200 cm−1 is related to the H-O-H bend and O-H stretch combination bands in water, while the band at around 4520 cm−1, usually associated with clay minerals, results from metal–OH bending plus O-H stretching. Based on the recorded NIR spectra, the identified spectral regions and the vibrational groups are relevant in the generation of calibration models, the overtones of O-H and H-O-H stretch vibrations of free water, and the overtones and combinations of O-H stretching and metal-OH bends in the clay lattice, and are related to the well-defined absorption patterns in soil reflectance across the NIR [35]. NIR spectra with this shape have been recorded for a variety of soils [36,37].

3.3. Comparison of Machine Learning Algorithms for Prediction of Soil Properties

The descriptive regression statistics of the predicted vs. measured soil properties using total ensemble, random forest optimization, gradient boosted machines, partial least squares, Bayesian additive regression trees and Cubist methods are provided in Table 2 and the scatterplots of the predicted vs. measured soil properties from different models are illustrated in Figure 3 and Figure 4. A good prediction was made for total carbon using the ensemble prediction (R2 of 0.87; RMSE = 0.39; RPIQ = 1.36 and RPD = 1.51) and Cubist prediction (R2 of 0.86; RMSE = 0.42; RPIQ = 1.26 and RPD = 1.40). Total nitrogen was also well-predicted using the total ensemble algorithm (R2 of 0.82; RMSE = 0.03; RPIQ = 2.00 and RPD = 1.60) and Bayesian additive regression trees (R2 of 0.78; RMSE = 0.04; RPIQ = 1.50 and RPD = 1.20). Carbon and nitrogen are two important constituents in soil organic matter, and their levels are highly connected [38]. Considering the classification by Zhao et al. [39] in terms of strong (R2 > 0.70), moderate (0.5 < R2 < 0.7), weak (0.3 < R2 < 0.5) and very weak (R2 < 0.3) prediction models, the total ensemble prediction method made: strong predictions (R2 > 0.70) for exchangeable bases (sum of Mehlich exch Ca, Mg, K, Na), Mehlich 3 extractable Cu, Mehlich 3 extractable Fe, Mehlich 3 extractable B, Mehlich 3 extractable Mn, exchangeable Na and exchangeable Ca; moderate predictions (0.5 < R2 < 0.7) for Mehlich 3 extractable Al, pH, the phosphorus sorption index (PSI) and exchangeable acidity (ExAc); weak predictions (0.3 < R2 < 0.5) for Mehlich 3 extractable Zn, Mehlich 3 extractable P, and electrical conductivity(Ecd); and a very weak (R2 < 0.3) prediction for Mehlich 3 extractable S (Table 2). The results show that total ensemble method consistently performs better than random forest optimization, gradient boosted machines, partial least squares, Bayesian additive regression trees and Cubist (lowest RMSE and highest RPIQ).
BART and Cubist made strong to moderate predictions, although cubist outperformed BART for total carbon (R2 of 0.86; RMSE = 0.42; RPIQ = 1.26 and RPD = 1.40), pH (R2 of 0.65; RMSE = 0.52; RPIQ = 4.02 and RPD = 2.58), Mehlich 3 extractable Al (R2 of 0.56; RMSE = 185.55; RPIQ = 2.58 and RPD = 1.79), Mehlich 3 extractable B (R2 of 0.71; RMSE = 0.34; RPIQ = 2.59 and RPD = 1.69), exchangeable Mg (R2 of 0.66; RMSE = 1.01; RPIQ = 1.80 and RPD = 1.76), exchangeable bases (R2 of 0.84; RMSE = 4.65; RPIQ = 6.54 and RPD = 3.94), Mehlich 3 extractable Fe (R2 of 0.69; RMSE = 32.01; RPIQ = 0.97 and RPD = 1.248) and Mehlich 3 extractable P (R2 of 0.41; RMSE = 16.79; RPIQ = 0.39 and RPD = 0.78) (Table 2). Exchangeable Ca and Mg have also been poorly predicted by other researchers [40,41], while EC and exchangeable Ca, Mg, K, and Na were poorly predicted with an RPD value < 1.3 by Pirie et al. [42]. Mehlich 3 extractable P, Mehlich 3 extractable Zn and electrical conductivity could not be predicted accurately with either model. Shepherd and Walsh [40] showed reasonable predictions for exchangeable Ca, Mg, and organic C in soils using NIR spectroscopy (R2 = 0.78–0.88), while Udelhoven et al. [43] were unable to obtain suitable predictions for organic C, plant accessible N, K, and P on land-scape scaled experiment using NIR spectroscopy.

3.4. Soil Health Indicators of Different Land Uses as Predicted by Total Ensemble Algorithm

Soil properties in Kenya, Tanzania and Uganda differed significantly among land-use types and between soil depths (Table 3). In Kenya, the effect of land use at 0–15 cm soil depth was significant for TN, TC, Mn, P, Zn, PSI, exchangeable Na, K and electrical conductivity; soils collected from community forest (CF) at 0–15 cm and 15–45 cm soil depth had significantly higher TN and TC than the control. Soil pH, which is a measure of the soil acidity and alkalinity of soil with different land uses, was only significant at 45–100 cm depth (p < 0.05). The soil pH values of grassland (8.08 ± 0.59) and the control (8.48 ± 0.30) were significantly higher than for cropland with soil and water conservation (CLSWC) (6.19 ± 0.65), community forest (6.416.94 ± 1.69) and agroforestry (6.41 ± 0.76; p < 0.05) (Table 3). The soils used for different land uses in Kenya are therefore suitable for agricultural activities (pH range of 6.62–7.46). Extremes in acidity or alkalinity of soil pH alter the available nutrients, resulting in imbalanced element uptake in plants [44]. Certain elements are more or less available at low soil acidic pH values (<5.5 units), thus resulting in nutrient constrains and unavailability in soils [45]. Maximum P availability occurs at soil pH range of 6.5–7.0; however, severe limitations with respect to P at 0–10 cm depth were observed in grassland (2.96 + 0.75 mg kg−1), cropland without soil and water conservation (5.58 + 2.35 mg kg−1) and community forest (5.50 + 3.02 mg kg−1) as compared to the critical value of 8.5 mg kg−1 suggested for most crops [46]. The effect of land use and soil depth was non-significant for exchangeable Ca and exchangeable Mg at 0–15 cm, 15–45 cm and 45–100 cm, respectively. At 0–15 cm soil depth, cropland with soil and water conservation (CLSWC) had significantly high ExK (2.55 ± 0.57 cmolc kg−1), compared to agroforestry (1.24 ± 0.47 cmolc kg−1), and cropland without soil and water conservation (0.92 ± 0.53 cmolc kg−1).
In Tanzania, the effect of land use on topsoil (0–15 cm) was significant for TN, TC, Mn, Al, Zn, PSI, exchangeable Ca, Mg, exbases and electrical conductivity (Table 3). TN was significantly high (p < 0.05) for the control (0.83 ± 0.03 %) than for cropland with soil and water conservation (0.18 ± 0.03 %) and other land uses at 0–15 cm soil depth. Soils from the community forest in Tanzania at 0–15 cm, 15–45 cm and 45–100 cm had significantly higher Fe contents than any other land use type, while community forest soils, collected at 0–15 cm and 15–45 cm soil depths had significantly higher Mn content than any other land-use type (p < 0.05; Table 3). Although Mn toxicity is ideal in acidic soils, soil pH across different land uses was not significant and ranged from 5.99 ± 0.46 to 6.47 ± 0.48, indicating that land uses have no significant effect on pH in Tanzania. Land use and soil depth were non-significant for exchangeable K and exchangeable Na, but remained significant for exchangeable Ca, exchangeable Mg at 0–15 cm, and 15–45 cm and 45–100 cm soil depths. Exchangeable K and exchangeable Ca in the soil also depends on the composition of parent rock materials [45].
In Uganda, land use and soil depth were significant for TN, TC, Mn, P, Zn, PSI, exchangeable Na and exchangeable K (Table 3). TN and TC were significantly high at the surface soil (0–15 cm) for five land use types compared to the control. Soil pH at the surface soil across different land uses were significant, with more acidic soils observed for the control (5.31 ± 0.59) and grasslands (5.79 ± 0.52). The highest soil pH (6.89 ± 0. 72) and the lowest soil pH (5.31 ± 0.59) were recorded for agroforestry and control land uses, respectively. Lower pH values for soils collected from grassland and the control could be as result of basic cation depletion in these soils. However, at 0–15 cm, 15–45 cm and 45–100 cm soil depths, soils from grassland recorded significantly higher Mn contents than the control (Table 3). The decline in soil pH across grassland could be due to encroachment into the grassland to create land for grazing and agricultural activities. Land use and soil depth did not affect exchangeable K and exchangeable Na, which could be a result of the removal of vegetation cover due to human and livestock interference [47]. A significantly high concentration of exchangeable Ca (18.21 ± 7.51 cmol/kg soil) for surface soils (0–15 cm) under agroforestry land use as compared to the control was probably due to decomposing litter that change the relative quantities of exchangeable base (Ca, Mg) and acid (Al, Fe) cations in soil [45]. Agroforestry through practices such as litter incorporation and manure application has been suggested to improve soil carbon sequestration [48].

4. Discussion

Total ensemble modeling was used as an alternative method for predicting soil quality indicators of different climate-smart land-use practices using soil NIR spectra. The ensemble modeling approaches have received little attention in contexts where calibration and prediction are common in chemistry [20]. The adoption of this method showed improved accuracy as compared to random forest optimization, gradient boosted machines, partial least squares, Bayesian additive regression trees and Cubist (Figure 3 and Figure 4). A good prediction was obtained for total carbon (R2 = 0.87), TN (R2 = 0.82), exchangeable bases (sum of Mehlich exch Ca, Mg, K, Na), Mehlich 3 extractable Cu, Mehlich 3 extractable Fe, Mehlich 3 extractable B, Mehlich 3 extractable Mn, exchangeable Na and exchangeable Ca. Exchangeable Ca and Mg have also been satisfactorily predicted (R2 > 0.70) by other authors [37,40,41]. Nevertheless, Chang et al. [38] found that some soil chemical or physical variables that do not have their primary response in the NIR region are poorly predicted, in comparison with biological indicators. However, the predictions of such properties can be achieved by correlations with biological properties [49]. Bian et al. [50] found that the linearity relationship of the total ensemble method was better than other modeling methods such as PLS.
Total ensemble modeling first generates diverse training subsets from the sample direction to build multiple sub-models [51], while gradient boosted machines (GBM) are hybrid methods that combines boosting and bagging techniques [28]. Cubist and random forests models have also been used to generate a more comprehensive predictive accuracy on NIR spectra [30,52]. The failure in calibration for Mehlich 3 extractable P and Mehlich 3 extractable Zn could probably be as a result of very narrow ranges, with extreme minimum values of 0.00 mg kg−1 for both Zn and P (Table 1). However, Minasny et al. [53] attributed the poor prediction of these properties to the fact that they are not related to the soil matrix or solid constituents. Additionally, electrical conductivity (R2 = 0.40) and pH (R2 = 0.65) do not exhibit a primary response in the NIR region and their prediction depends on their relationships with organic matter and clay contents [54]. Unsatisfactory predictions of Mehlich 3 extractable P and electrical conductivity (R2 < 0.66 and RPD < 2) were also observed by other authors [42,54]. In general, electrolyte concentration-related parameters such as electrical conductivity (ECd) are not accurately predicted [55].
The concepts of soil health and soil quality are based on the ability of a given soil to sustain production and measure both positive and negative changes in soil health [46]. Soil characteristics and overall soil health are influenced by land-use systems and intrinsic features, with soil nutrient availability and soil metal concentrations accounting for the majority of the variation seen. Additionally, land-use patterns have an impact on the soil nutrient distribution and availability by affecting soil characteristics [56]. However, the different soil pH ranges for different climate-smart land-use types observed in this study suggest that soil pH is unlikely to have a negative impact on plant development and agricultural productivity. Some of the physical-chemical qualities of soil are static across time, while others are dynamic across a range of time scales. Some are resistant to change brought about by land-use types, whereas others are quickly influenced in both positive and negative ways by land use [57]. In most of the quality indicators assessed at different soil depths, agroforestry, community forest, cropland with soil and water conservation, cropland without soil and water conservation and grassland land uses resulted in better soil quality as compared to controlled land use. As a result, the findings of this study showed the potential of different climate-smart land-use types to improve or preserve soil quality. Soil quality changes are typically the result of increased erosion processes caused by plowing, burning, overgrazing, and other land-use practices that remove the protective vegetative cover [58]. Using a threshold value of 0.2% TN [25] to assess the agricultural productivity of soils at different depths for different land uses, the majority of topsoils (0–15 cm), mid depth (15–45 cm), and at a depth of 45–100 cm had >0.2% N. The results agree with other authors who concluded that climate-smart land use can reduce the degradation of soil quality. Adopting climate-smart land-use practices is a prerequisite to ensure soil health and crop production for sustainable food production systems. These practices, such as agroforestry, community forestry, cropland with soil and water conservation, and grasslands, can increase soil microbial activities and the cation exchange capacity [59]. The soil health indicators used for the assessment of any soil are based on the intended use of the soil. Hence, the set of soil parameters used to assess soil for sustainable rangeland may differ from those used to assess the soil for sustainable crop production, grasslands, and community forests. As a result, there are no universal soil health indicators that will work in all systems. Some soil health indicators, such as water-holding capacity and soil texture, do not change rapidly over time, while other soil health indicators fluctuate more throughout a field and over time according to soil management, crop variety, precipitation, and temperature.

5. Conclusions

This study highlighted the robustness and the predictive power of using total ensemble modeling on NIR spectra for predicting the soil chemical and physical properties as soil health indicators. The results are better when using the total ensemble algorithm as opposed to the random forest optimization, gradient boosted machines, partial least squares and Bayesian additive regression trees algorithms. Further, the developed models were able to predict the physical and chemical properties of soil collected from soils used for different land-use practices at different depths. The soils of the research locations are suitable for sustainable food production systems according to the current assessment of soil health indicators; however, land-use systems and sampling depth, as well as intrinsic traits, influence the overall soil health, with soil nutrient availability accounting for the majority of the variation observed. Thus, the findings of this study highlighted the land-use categories’ potential for improving or sustaining soil quality. As a result, there is reason for adopting climate-smart agriculture (CSA) practices as a viable alternative.

Author Contributions

Conceptualization, J.W.R.; K.O.O.; A.S.; M.R.; and D.S.; methodology, J.W.R.; K.O.O.; and A.S.; software, J.W.R.; K.O.O.; and A.S.; validation, J.W.R.; K.O.O.; and A.S.; formal analysis, J.W.R.; K.O.O.; and A.S.; investigation, J.W.R.; K.O.O.; A.S.; and G.A.; resources, J.W.R.; K.O.O.; A.S.; M.R.; and D.S.; data curation, J.W.R.; K.O.O.; A.S.; and G.A.; writing—original draft preparation, J.W.R.; K.O.O.; and A.S.; writing—review and editing, J.W.R.; K.O.O.; A.S.; G.A.; M.R.; and D.S.; visualization, J.W.R.; K.O.O.; A.S.; G.A.; M.R.; and D.S.; supervision, J.W.R.; M.R; and D.S.; project administration, J.W.R.; M.R.; and D.S.; funding acquisition, J.W.R.; M.R.; and D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported through CCAFS-EC grant reference: 2000002575 for the project on Building Livelihoods and Resilience to Climate Change in East and West Africa: Agricultural Research for Development (AR4D) for large-scale implementation of Climate-Smart Agriculture. The funds are administered by the International Fund for Agricultural Development (IFAD), Rome, Italy, while the project was implemented by Alliance Bioversity-CIAT. The Article Processing Charges were paid through a grant from the World Bank “Accelerating Impacts of CGIAR Climate Research for Africa (AICCRA)” ESA regional project (Grant No. D7540) signed between CIAT and the International Development Association (IDA).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to confidentiality requirements of the funding organisations for the project.

Acknowledgments

The Kenya Agricultural and Livestock Research Organization (KALRO) and Vi Agroforestry are acknowledged for providing support during fieldwork in Nyando CSVs. Maseno University is acknowledged for the bulk density determination of Nyando CSVs soil samples. The Tanzania Agricultural Research Institute (TARI) is acknowledged for providing support during fieldwork as well as for bulk density determination of the Lushoto CSVs soil samples. The Uganda National Agricultural Research Organisation (NARO) is acknowledged for providing support during fieldwork as well as for bulk density determination of the Hoima CSVs soil samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Segnini, A.; Posadas, A.; da Silva, W.T.L.; Milori, D.M.B.P.; Gavilan, C.; Claessens, L.; Quiroz, R. Quantifying soil carbon stocks and humification through spectroscopic methods: A scoping assessment in EMBU-Kenya. J. Environ. Manag. 2019, 234, 476–483. [Google Scholar] [CrossRef] [PubMed]
  2. Comino, F.; Aranda, V.; García-Ruiz, R.; Ayora-Cañada, M.J.; Domínguez-Vidal, A. Infrared spectroscopy as a tool for the assessment of soil biological quality in agricultural soils under contrasting management practices. Ecol. Indic. 2018, 87, 117–126. [Google Scholar] [CrossRef]
  3. Chodak, M.; Khanna, P.; Horvath, B.; Beese, F. Near infrared spectroscopy for determination of total and exchangeable cations in geologically heterogeneous forest soils. J. Near Infrared Spectrosc. 2004, 12, 315–324. [Google Scholar] [CrossRef]
  4. Nosrati, K. Assessing soil quality indicator under different land use and soil erosion using multivariate statistical techniques. Environ. Monit. Assess. 2013, 185, 2895–2907. [Google Scholar] [CrossRef] [PubMed]
  5. Butkute, B.; Šlepetiene, A. Application of near infrared reflectance spectroscopy for the assessment of soil quality in a long-term pasture. Commun. Soil Sci. Plant Anal. 2006, 37, 2389–2409. [Google Scholar] [CrossRef]
  6. Cécillon, L.; Barthès, B.G.; Gomez, C.; Ertlen, D.; Genot, V.; Hedde, M.; Stevens, A.; Brun, J.J. Assessment and monitoring of soil quality using near-infrared reflectance spectroscopy (NIRS). Eur. J. Soil Sci. 2009, 60, 770–784. [Google Scholar] [CrossRef] [Green Version]
  7. Ahmadi, A.; Emami, M.; Daccache, A.; He, L. Soil Properties Prediction for Precision Agriculture Using Visible and Near-Infrared Spectroscopy: A Systematic Review and Meta-Analysis. Agronomy 2021, 11, 433. [Google Scholar] [CrossRef]
  8. Soriano-Disla, J.M.; Janik, L.J.; Rossel, R.A.V.; Macdonald, L.M.; McLaughlin, M.J. The Performance of Visible, Near-, and Mid-Infrared Reflectance Spectroscopy for Prediction of Soil Physical, Chemical, and Biological Properties. Appl. Spectrosc. Rev. 2014, 49, 139–186. [Google Scholar] [CrossRef]
  9. Pinheiro, É.F.M.; Ceddia, M.B.; Clingensmith, C.M.; Grunwald, S.; Vasques, G.M. Prediction of soil physical and chemical properties by visible and near-infrared diffuse reflectance spectroscopy in the Central Amazon. Remote Sens. 2017, 9, 293. [Google Scholar] [CrossRef] [Green Version]
  10. Ozaki, Y.; Huck, C.W.; Beć, K.B. Near-IR Spectroscopy and Its Applications. Mol. Laser Spectrosc. Adv. Appl. 2018, 11–38. [Google Scholar] [CrossRef]
  11. De Mastro, F.; Cocozza, C.; Brunetti, G.; Traversa, A. Chemical and spectroscopic investigation of different soil fractions as affected by soil management. Appl. Sci. 2020, 10, 2571. [Google Scholar] [CrossRef] [Green Version]
  12. Nduwamungu, C.; Ziadi, N.; Tremblay, G.F.; Parent, L.-É. Near-Infrared Reflectance Spectroscopy Prediction of Soil Properties: Effects of Sample Cups and Preparation. Soil Sci. Soc. Am. J. 2009, 73, 1896–1903. [Google Scholar] [CrossRef]
  13. Askari, M.S.; O’Rourke, S.M.; Holden, N.M. Evaluation of soil quality for agricultural production using visible-near-infrared spectroscopy. Geoderma 2015, 243–244, 80–91. [Google Scholar] [CrossRef]
  14. Kuang, B.; Mouazen, A.M. Non-biased prediction of soil organic carbon and total nitrogen with vis-NIR spectroscopy, as affected by soil moisture content and texture. Biosyst. Eng. 2013, 114, 249–258. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, H.; Tan, C.; Lin, Z.; Wu, T. Classification of different liquid milk by near-infrared spectroscopy and ensemble modeling. Spectrochim. Acta-Part A Mol. Biomol. Spectrosc. 2021, 251, 119460. [Google Scholar] [CrossRef]
  16. Dietterich, T.G. An Experimental Comparison of Three Methods for Constructing Ensembles of Decision Trees: Bagging, Boosting, and Randomization. Mach. Learn. 1999, 22, 1–22. [Google Scholar]
  17. Elder, J. The Apparent Paradox of Complexity in Ensemble Modeling. In Handbook of Statistical Analysis and Data Mining Applications; Academic Press: Cambridge, MA, USA, 2018; pp. 705–718. ISBN 978-0-12-416632-5. [Google Scholar]
  18. Freund, Y.; Schapire, R.E. A short introduction to Boosting. J. Jpn. Soc. Artif. Intell. 1999, 14, 771–780. [Google Scholar]
  19. Genuer, R.; Poggi, J.M.; Tuleau-Malot, C. Variable selection using random forests. Pattern Recognit. Lett. 2010, 31, 2225–2236. [Google Scholar] [CrossRef] [Green Version]
  20. Mevik, B.H.; Segtnan, V.H.; Næs, T. Ensemble methods and partial least squares regression. J. Chemom. 2004, 18, 498–507. [Google Scholar] [CrossRef]
  21. Shao, X.; Bian, X.; Cai, W. An improved boosting partial least squares method for near-infrared spectroscopic quantitative analysis. Anal. Chim. Acta 2010, 666, 32–37. [Google Scholar] [CrossRef]
  22. Reda, R.; Saffaj, T.; Ilham, B.; Saidi, O.; Issam, K.; Brahim, L.; El Hadrami, E.M. A comparative study between a new method and other machine learning algorithms for soil organic carbon and total nitrogen prediction using near infrared spectroscopy. Chemom. Intell. Lab. Syst. 2019, 195. [Google Scholar] [CrossRef]
  23. Shao, X.; Bian, X.; Liu, J.; Zhang, M.; Cai, W. Multivariate calibration methods in near infrared spectroscopic analysis. Anal. Methods 2010, 2, 1662–1666. [Google Scholar] [CrossRef]
  24. Chen, H.; Tan, C.; Lin, Z. Ensemble of extreme learning machines for multivariate calibration of near-infrared spectroscopy. Spectrochim. Acta-Part A Mol. Biomol. Spectrosc. 2020, 229, 117982. [Google Scholar] [CrossRef] [PubMed]
  25. Winowiecki, L.; Vågen, T.G.; Massawe, B.; Jelinski, N.A.; Lyamchai, C.; Sayula, G.; Msoka, E. Landscape-scale variability of soil health indicators: Effects of cultivation on soil organic carbon in the Usambara Mountains of Tanzania. Nutr. Cycl. Agroecosystems 2016, 105, 263–274. [Google Scholar] [CrossRef] [Green Version]
  26. Recha, J.; Radeny, M.; Kinyangi, J.; Kimeli, P.; Atakos, V.; Lyamchai, C.; Ngatoluwa, R.; Sayula, G. Climate-Smart Villages and Progress in Achieving Household Food Security in Lushoto, Tanzania; CCAFS Info Note; CGIAR Research Program on Climate Change, Agriculture and Food Security (CCAFS): Copenhagen, Denmark, 2015. [Google Scholar]
  27. Kinyangi, J.; Recha, J.; Kimeli, P.; Atakos, V. Climate-Smart Villages and the Hope of Food Security in Kenya; CCAFS Info Note; CGIAR Research Program on Climate Change, Agriculture and Food Security (CCAFS): Copenhagen, Denmark, 2015. [Google Scholar]
  28. Friedman, J.H. Stochastic gradient boosting. Comput. Stat. Data Anal. 2002, 38, 367–378. [Google Scholar] [CrossRef]
  29. Chipman, H.A.; George, E.I.; McCulloch, R.E. BART: Bayesian additive regression trees. Ann. Appl. Stat. 2012, 6, 266–298. [Google Scholar] [CrossRef]
  30. Zhou, J.; Li, E.; Wei, H.; Li, C.; Qiao, Q.; Armaghani, D.J. Random forests and cubist algorithms for predicting shear strengths of rockfill materials. Appl. Sci. 2019, 9, 1621. [Google Scholar] [CrossRef] [Green Version]
  31. Metzger, K.; Zhang, C.; Ward, M.; Daly, K. Mid-infrared spectroscopy as an alternative to laboratory extraction for the determination of lime requirement in tillage soils. Geoderma 2020, 364, 114171. [Google Scholar] [CrossRef]
  32. Bellon-Maurel, V.; Fernandez-Ahumada, E.; Palagos, B.; Roger, J.M.; McBratney, A. Critical review of chemometric indicators commonly used for assessing the quality of the prediction of soil attributes by NIR spectroscopy. TrAC-Trends Anal. Chem. 2010, 29, 1073–1081. [Google Scholar] [CrossRef]
  33. O’Rourke, S.M.; Stockmann, U.; Holden, N.M.; McBratney, A.B.; Minasny, B. An assessment of model averaging to improve predictive power of portable vis-NIR and XRF for the determination of agronomic soil properties. Geoderma 2016, 279, 31–44. [Google Scholar] [CrossRef]
  34. R Development Core Team. R: A Language and Environment for Statistical Computing. 2014, p. 2673. Available online: https://www.r-project.org/ (accessed on 6 October 2021).
  35. Nocita, M.; Stevens, A.; van Wesemael, B.; Aitkenhead, M.; Bachmann, M.; Barthès, B.; Dor, E.B.; Brown, D.J.; Clairotte, M.; Csorba, A.; et al. Soil Spectroscopy: An Alternative to Wet Chemistry for Soil Monitoring. Adv. Agron. 2015, 132, 139–159. [Google Scholar] [CrossRef]
  36. Ludwig, B.; Murugan, R.; Parama, V.R.R.; Vohland, M. Use of different chemometric approaches for an estimation of soil properties at field scale with near infrared spectroscopy. J. Plant Nutr. Soil Sci. 2018, 181, 704–713. [Google Scholar] [CrossRef]
  37. Pudełko, A.; Chodak, M. Estimation of total nitrogen and organic carbon contents in mine soils with NIR reflectance spectroscopy and various chemometric methods. Geoderma 2020, 368, 114306. [Google Scholar] [CrossRef]
  38. Chang, C.; Laird, D.; Mausbach, M.J. 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]
  39. Zhao, D.; Arshad, M.; Li, N.; Triantafilis, J. Predicting soil physical and chemical properties using vis-NIR in Australian cotton areas. Catena 2021, 196, 104938. [Google Scholar] [CrossRef]
  40. 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]
  41. Moron, A.; Cozzolino, D. Exploring the use of near infrared reflectance spectroscopy to study physical properties and microelements in soils. J. Near Infrared Spectrosc. 2003, 11, 145–154. [Google Scholar] [CrossRef]
  42. Pirie, A.; Singh, B.; Islam, K. Spectroscopic Techniques to Predict Several Soil Properties. Aust. J. Soil Res. 2005, 43, 713–721. [Google Scholar] [CrossRef]
  43. Udelhoven, T.; Emmerling, C.; Jarmer, T. Quantitative analysis of soil chemical properties with diffuse reflectance spectrometry and partial least-square regression: A feasibility study Thomas. Plant Soil 2003, 251, 319–329. [Google Scholar] [CrossRef]
  44. Zhao, J.; Dong, Y.; Xie, X.; Li, X.; Zhang, X.; Shen, X. Effect of annual variation in soil pH on available soil nutrients in pear orchards. Acta Ecol. Sin. 2011, 31, 212–216. [Google Scholar] [CrossRef]
  45. Takoutsing, B.; Weber, J.C.; Tchoundjeu, Z.; Shepherd, K. Soil chemical properties dynamics as affected by land use change in the humid forest zone of Cameroon. Agrofor. Syst. 2016, 90, 1089–1102. [Google Scholar] [CrossRef]
  46. Takoutsing, B.; Weber, J.; Aynekulu, E.; Rodríguez Martín, J.A.; Shepherd, K.; Sila, A.; Tchoundjeu, Z.; Diby, L. Assessment of soil health indicators for sustainable production of maize in smallholder farming systems in the highlands of Cameroon. Geoderma 2016, 276, 64–73. [Google Scholar] [CrossRef]
  47. Moges, A.; Dagnachew, M.; Yimer, F. Land use effects on soil quality indicators: A case study of Abo-Wonsho Southern Ethiopia. Appl. Environ. Soil Sci. 2013, 2013, 1–9. [Google Scholar] [CrossRef]
  48. Ambaw, G.; Recha, J.W.; Nigussie, A.; Solomon, D.; Radeny, M. Soil carbon sequestration potential of climate-smart villages in East African countries. Climate 2020, 8, 124. [Google Scholar] [CrossRef]
  49. Sanderman, J.; Savage, K.; Dangal, S.R.S. Mid-infrared spectroscopy for prediction of soil health indicators in the United States. Soil Sci. Soc. Am. J. 2020, 84, 251–261. [Google Scholar] [CrossRef] [Green Version]
  50. Bian, X.; Wang, K.; Tan, E.; Diwu, P.; Zhang, F.; Guo, Y. A selective ensemble preprocessing strategy for near-infrared spectral quantitative analysis of complex samples. Chemom. Intell. Lab. Syst. 2020, 197, 103916. [Google Scholar] [CrossRef]
  51. Wang, K.; Bian, X.; Tan, X.; Wang, H.; Li, Y. A new ensemble modeling method for multivariate calibration of near infrared spectra. Anal. Methods 2021, 13, 1374–1380. [Google Scholar] [CrossRef]
  52. Henaka Arachchi, M.P.N.K.; Field, D.J.; McBratney, A.B. Quantification of soil carbon from bulk soil samples to predict the aggregate-carbon fractions within using near- and mid-infrared spectroscopic techniques. Geoderma 2016, 267, 207–214. [Google Scholar] [CrossRef]
  53. Minasny, B.; Tranter, G.; McBratney, A.B.; Brough, D.M.; Murphy, B.W. Regional transferability of mid-infrared diffuse reflectance spectroscopic prediction for soil chemical properties. Geoderma 2009, 153, 155–162. [Google Scholar] [CrossRef]
  54. 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]
  55. Viscarrarossel, R.; Walvoort, D.; Mcbratney, A.; Janik, L.; Skjemstad, J. 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]
  56. Mulat, Y.; Kibret, K.; Bedadi, B.; Mohammed, M. Soil quality evaluation under different land use types in Kersa sub-watershed, eastern Ethiopia. Environ. Syst. Res. 2021, 10. [Google Scholar] [CrossRef]
  57. Xia, Y.; Ugarte, C.M.; Guan, K.; Pentrak, M.; Wander, M.M.; Comino, F.; Aranda, V.; García-Ruiz, R.; Ayora-Cañada, M.J.; Domínguez-Vidal, A.; et al. Assessment of soil quality of arable soils in Hungary using DRIFT spectroscopy and chemometrics. Soil Sci. Soc. Am. J. 2018, 11, 80–91. [Google Scholar]
  58. Arias, M.E.; González-Pérez, J.A.; González-Vila, F.J.; Ball, A.S. Soil health—A new challenge for microbiologists and chemists. Int. Microbiol. 2005, 8, 13–21. [Google Scholar] [CrossRef] [PubMed]
  59. Tahat, M.M.; Alananbeh, K.M.; Othman, Y.A.; Leskovar, D.I. Soil health and sustainable agriculture. Sustainability 2020, 12, 4859. [Google Scholar] [CrossRef]
Figure 1. Biplot for the principal components explaining the largest variation.
Figure 1. Biplot for the principal components explaining the largest variation.
Soilsystems 05 00069 g001
Figure 2. NIR raw spectral signatures for all soil samples.
Figure 2. NIR raw spectral signatures for all soil samples.
Soilsystems 05 00069 g002
Figure 3. Scatter plot of observed and predicted total carbon content using RFO, GBM, PLS, BART, CUBIST and Total Ensemble.
Figure 3. Scatter plot of observed and predicted total carbon content using RFO, GBM, PLS, BART, CUBIST and Total Ensemble.
Soilsystems 05 00069 g003
Figure 4. Scatter plot of observed and predicted total nitrogen content values using RFO, GBM, PLS, BART, CUBIST and total ensemble.
Figure 4. Scatter plot of observed and predicted total nitrogen content values using RFO, GBM, PLS, BART, CUBIST and total ensemble.
Soilsystems 05 00069 g004
Table 1. Summary statistics for laboratory soil properties measured as potential indicators of soil health.
Table 1. Summary statistics for laboratory soil properties measured as potential indicators of soil health.
Soil PropertyNMin.MedianMax.Mean ± SDRangeIQRSkewnessCV%Kurtosis
Total Nitrogen3150.020.110.830.13 ± 0.100.270.060.3371.300.64
Total Carbon3150.301.4611.861.79 ± 1.323.440.532.5973.419.34
Sand3150.543.2917.804.56 ± 3.711721.9181.363.49
Silt3151.686.2256.399.67 ± 8.68551031.3289.6914.81
Clay31531.7188.9496.2285.85 ± 10.726510−2.8112.4811.34
pH3154.436.359.086.44 ± 1.344.652.090.2714.97−0.83
m3.Al315456.00951.002700.001006.97 ± 323.301274.00479.000.4532.11−0.55
m3.B3150.000.654.180.77 ± 0.642.050.880.8982.38−0.26
m3.Cu3150.003.0716.003.47 ± 2.637.342.081.2575.761.44
m3.Fe31523.9092.10436.00107.55 ± 66.52238.1031.102.3061.858.20
m3.Mn3150.00214.00660.00215.29 ± 155.74390186.400.0772.34−1.15
m3.P3150.001.91166.006.98 ± 18.9685.406.524.70271.5626.77
m3.S3150.003.29226.009.24 ± 22.76151.0017.582.49246.355.61
m3.Zn3150.001.0132.302.14 ± 3.4214.001.104.80159.9724.96
PSI3150.98116.00655.00137.08 ± 87.23332.00132.050.8763.64−0.27
ExNa3150.000.0511.700.66 ± 1.6010.823.022.03240.484.03
ExCa3150.318.6044.0512.46 ± 10.6743.4923.511.0585.65−0.49
ExMg3150.073.179.833.26 ± 1.776.501.821.1054.360.09
ExK3150.000.285.170.72 ± 0.873.251.491.10119.620.16
ExBas3150.4912.2558.2617.11 ± 13.5656.7730.421.0379.24−0.58
ECd3150.010.050.770.08 ± 0.090.760.171.93108.763.84
ExAc3150.000.008.750.27 ± 0.944.870.2493.00344.639.02
m3. = Mehlich 3 extractable; PSI = Phosphorus sorption index; ExNa = Exchangeable Na; ExCa = Exchangeable Ca; ExMg = Exchangeable Mg; ExK = Exchangeable K; ExBas = Exchangeable bases (sum of Mehlich exch Ca, Mg, K, Na); Ecd = Electrical conductivity; ExAc = Exchangeable Acidity; IQR = Interquartile Range; SD = standard deviation.
Table 2. Comparison of prediction accuracy from six machine learning algorithms using soil NIR spectra.
Table 2. Comparison of prediction accuracy from six machine learning algorithms using soil NIR spectra.
Soil PropertyMethodR2RMSERPIQRPDSoil PropertyMethodR2RMSERPIQRPD
Total CarbonRFO0.830.461.151.28m3.CuRFO0.681.301.601.38
GBM0.750.590.901.00 GBM0.591.511.381.19
PLS0.780.550.961.07 PLS0.342.091.000.86
BART0.820.521.021.13 BART0.721.221.701.47
CUBIST0.860.421.261.40 CUBIST0.691.271.641.41
ENS0.870.391.361.51 ENS0.731.21.731.50
Total NitrogenRFO0.750.041.501.20m3.FeRFO0.6334.960.891.14
GBM0.760.041.501.20 GBM0.4542.550.730.93
PLS0.720.041.501.20 PLS0.5340.980.760.97
BART0.780.041.501.20 BART0.5439.800.781.00
CUBIST0.670.051.200.96 CUBIST0.6932.010.971.24
ENS0.820.032.001.60 ENS0.7329.671.051.34
pHRFO0.560.583.602.31m3.MnRFO0.65103.361.801.12
GBM0.560.603.482.23 GBM0.49125.621.480.92
PLS0.460.663.172.03 PLS0.21212.830.880.54
BART0.570.583.602.31 BART0.7292.262.021.25
CUBIST0.650.524.022.58 CUBIST0.7099.431.871.16
ENS0.660.514.102.63 ENS0.7585.302.191.35
m3.AlRFO0.49201.582.381.65m3.PRFO0.2618.650.350.70
GBM0.41212.012.261.56 GBM0.1719.660.330.67
PLS0.62169.482.831.96 PLS0.0526.780.240.49
BART0.53190.962.511.74 BART0.1619.840.330.66
CUBIST0.56185.552.581.79 CUBIST0.4116.790.390.78
ENS0.68157.123.052.11 ENS0.4116.580.390.79
m3.BRFO0.610.392.261.47m3.SRFO0.0315.51.132.23
GBM0.650.392.261.47 GBM0.0113.121.342.63
PLS0.520.481.831.19 PLS0.1112.681.392.73
BART0.620.382.321.51 BART0.0014.081.252.46
CUBIST0.710.342.591.69 CUBIST0.0212.941.362.67
ENS0.730.322.751.79 ENS0.1411.731.502.95
m3.ZnRFO0.332.590.420.86ExNaRFO0.740.585.214.34
GBM0.402.500.440.89 GBM0.570.744.083.40
PLS0.272.840.390.79 PLS0.231.092.772.31
BART0.442.490.440.90 BART0.720.65.034.19
CUBIST0.402.450.450.91 CUBIST0.750.575.304.41
ENS0.492.240.491.00 ENS0.810.506.045.03
PSIRFO0.3263.642.071.46ExCaRFO0.813.916.013.57
GBM0.3860.692.181.53 GBM0.794.295.483.25
PLS0.4457.962.281.60 PLS0.597.033.341.99
BART0.3763.322.091.46 BART0.803.975.923.52
CUBIST0.3761.382.151.51 CUBIST0.853.516.703.98
ENS0.5252.612.511.76 ENS0.853.476.784.02
ExMgRFO0.551.141.601.56ExKRFO0.400.662.261.38
GBM0.501.221.491.46 GBM0.330.712.101.28
PLS0.201.830.990.97 PLS0.220.811.841.12
BART0.541.141.601.56 BART0.470.622.401.47
CUBIST0.661.011.801.76 CUBIST0.480.622.401.47
ENS0.670.961.901.85 ENS0.510.602.481.52
ExBasRFO0.805.145.923.56ECdRFO0.360.062.832.71
GBM0.775.745.303.19 GBM0.370.053.403.25
PLS0.598.863.432.07 PLS0.300.053.403.25
BART0.795.255.793.49 BART0.230.072.432.32
CUBIST0.844.666.533.93 CUBIST0.350.062.832.71
ENS0.844.656.543.94 ENS0.400.053.401.38
ExAcRFO0.280.670.371.52
GBM0.280.680.371.50
PLS0.430.770.321.32
BART0.380.660.381.54
CUBIST0.340.640.391.59
ENS0.520.550.451.85
RFO = Random Forest Optimization; GBM = Gradient boosted machines; PLS = Partial least squares; BART = Bayesian additive regression trees; RPIQ = IQ/RMSE; RPD = SD/SEP.
Table 3. Soil quality indicators in different land uses in three countries at different depths as predicted by the ensemble algorithm on NIR spectra.
Table 3. Soil quality indicators in different land uses in three countries at different depths as predicted by the ensemble algorithm on NIR spectra.
CountryDepth (cm)Land Use (n)TN
%
T C
%
pH
Units
Al
mg kg−1
Cu
mg kg−1
Fe
mg kg−1
Mn
mg kg−1
P
mg kg−1
Zn
mg kg−1
PSI
Units
ExNa
cmolc kg−1
ExCa
cmolc kg−1
ExMg
cmolc kg−1
ExK
cmolc kg−1
ExBas
cmolc kg−1
ECd
cmolc kg−1
KENYA0–15AF(6)0.15 ± 0.06 ab2.14 ± 0.89 ab6.72 ± 0.89 a676.33 ± 205.70 a1.66 ± 0.99 a168.4 ± 80.65 a338.67 ± 75.46 a61.74 ± 61.97 ab8.71 ± 11.57 ab61.66 ± 56.70 b0.06 ± 0.05 b18.34 ± 12.80 a3.67 ± 1.62 a1.24 ± 0.47 bc23.31 ± 14.69 a0.84 ± 0.05 b
CF(6)0.24 ± 0.16 a3.88 ± 2.82 a7.09 ± 0.73 a818.83 ± 153.70 a2.86 ± 2.02 a153.00 ± 52.00 a269.50 ± 56.87 ab5.50 ± 3.02 b4.73 ± 3.91 ab136.00 ± 17.75 a1.08 ± 1.36 b28.78 ± 79.44 a5.97 ± 0.80 a2.03 ± 0.71 ab37.85 ± 9.72 a0.22 ± 0.14 ab
CLSWC(6)0.23 ± 0.05 ab3.22 ± 0.62 ab7.04 ± 0.73 a843.83 ± 103.34 a4.33 ± 0.29 a145.00 ± 12.31 a269.50 ± 56.87 ab68.95 ± 53.30 a12.97 ± 7.01 ab89.55 ± 30.15 ab0.04 ± 0.03 b18.69 ± 6.06 a4.33 ± 0.72 a2.55 ± 0.57 a25.62 ± 6.72 a0.12 ± 0.05 ab
CLNSWC(6)0.13 ± 0.06 ab1.78 ± 0.78 ab6.62 ± 0.89 a778.67 ± 48.32 a2.13 ± 1.38 a221.67 ± 78.57 a301.00 ± 81.19 b5.58 ± 2.35 b2.33 ± 0.58 b64.02 ± 32.80 ab0.67 ± 0.26 b17.32 ± 9.62 a3.09 ± 1.92 a0.92 ± 0.53 c21.98 ± 12.11 a0.06 ± 0.02 ab
GL(6)0.14 ± 0.06 ab2.11 ± 0.47 ab6.93 ± 1.08 a765.00 ± 73.52 a2.60 ± 1.86 a217.88 ± 130.59 a169.97 ± 75.45 ab2.96 ± 0.75 b1.71 ± 0.44 b99.05 ± 20.06 ab1.04 ± 0.38 b22.92 ± 15.52 a3.33 ± 2.17 a1.59 ± 0.26 c28.22 ± 18.34 a0.76 ± 0.03 b
C(6)0.09 ± 0.02 b1.27 ± 0.41 b7.46 ± 0.89 a804.17 ± 132.09 a2.25 ± 1.68 a145.40 ± 62.75 a254.33 ± 91.44 ab56.27 ± 6.32 b1.16 ± 0.28 b99.88 ± 24.82 ab3.52 ± 2.22 a21.91 ± 10.71 a3.98 ± 2.06 a1.54 ± 0.78 bc34.52 ± 10.45 a0.26 ± 0.17 a
15–45AF(6)0.10 ± 0.02 ab 1.49 ± 0.42 ab6.70 ± 0.99 a789.67 ± 172.60 ab1.79 ± 1.50 a189.17 ± 101.92 a255.17 ± 107.79 a13.91 ± 11.28 a2.99 ± 2.78 a95.08 ± 46.17 ab0.19 ± 0.17 b15.37 ± 5.48 a4.09 ± 1.24 a1.43 ± 0.81 abc 21.07 ± 7.12 a0.08 ± 0.06 b
CF(6)0.10 ± 0.04 ab 1.67 ± 0.60 ab 6.82 ± 1.53 a941.17 ± 342.97 ab2.83 ± 2.21 a170.03 ± 92.27 a301.17 ± 61.11 a3.00 ± 3.01 a2.35 ± 1.22 a149.45 ± 50.47 a1.73 ± 2.02 b23.23 ± 15.02 a4.79 ± 0.83 a1.32 ± 0.58 abc31.07 ± 17.57 a0.12 ± 1.10 ab
CLSWC(6)0.15 ± 0.04 a2.04 ± 0.65 a6.55 ± 0.59 a999.67 ± 198.04 a4.07 ± 4.17 a124.35 ± 35.66 a191.60 ± 90.40 b 17.72 ± 30.79 a4.33 ± 4.78 a128.52 ± 42.29 ab0.08 ± 0.05 b13.70 ± 3.65 a5.01 ± 1.85 a1.86 ± 0.58 ab 20.66 ± 5.07 a0.06 ± 0.02 b
CLNSWC(6)0.10 ± 0.05 ab 1.54 ± 0.61 ab 6.94 ± 0.85 a694.67 ± 89.00 ab2.00 ± 1.64 a166.50 ± 34.64 a105.85 ± 61.54 ab 2.37 ± 1.59 a2.94 ± 1.95 a68.67 ± 46.99 b1.49 ± 1.13 b18.15 ± 12.27 a3.20 ± 2.44 a0.88 ± 0.46 bc 23.71 ± 16.11 a0.12 ± 0.10 ab
GL(6)0.07 ± 0.03 b 1.23 ± 0.17 ab7.33 ± 1.13 a605.50 ± 38.40 b2.82 ± 2.52 a153.03 ± 83.50 a175.27 ± 98.26 ab2.48 ± 1.58 a4.78 ± 1.75 a78.78 ± 27.82 ab 1.89 ± 1.43 b23.48 ± 15.32 a3.32 ± 2.64 ab1.77 ± 0.35 c29.46 ± 23.72 a0.12 ± 0.07 ab
C(6)0.07 ± 0.01 b1.07 ± 0.59 b8.03 ± 0.60 a774.17 ± 153.54 ab2.73 ± 1.40 a92.47 ± 925.57 a265.50 ± 77.56 ab12.24 ± 8.91 a2.24 ± 1.16 a83.65 ± 19.92 ab 4.64 ± 2.41 a30.50 ± 6.07 a4.37 ± 2.00 a1.97 ± 0.58 a41.47 ± 8.38 a0.25 ± 0.09 a
45–100AF(6)0.05 ± 0.01 a0.74 ± 0.23 a6.41 ± 0.76 cd863.17 ± 316.21 ab1.57 ± 1.64 a163.83 ± 47.49 a221.00 ± 41.69 a15.52 ± 15.63 a1.30 ± 0.75 a111.23 ± 84.40 ab0.29 ± 0.30 b11.21 ± 6.80 b3.67 ± 1.38 a1.77 ± 0.74 ab16.94 ± 7.98 b0.06 ± 0.05 b
CF(6)0.07 ± 0.04 a1.27 ± 0.62 a6.94 ± 1.69 bcd928.50 ± 362.77 ab2.83 ± 2.07 a153.90 ± 89.84 ab304.17 ± 114.36 a3.92 ± 4.61 a1.63 ± 0.78 a132.88 ± 60.88 ab2.28 ± 2.28 b23.40 ± 15.27 ab4.48 ± 0.91 b1.16 ± 0.26 b31.31 ± 18.35 ab0.12 ± 0.11 ab
CLSWC(6)0.09 ± 0.03 a1.37 ± 0.67 a6.19 ± 0.65 d1124.33 ± 151.24 cd3.08 ± 3.63 a112.43 ± 20.22 abc132.82 ± 40.16 a2.28 ± 2.79 a0.66 ± 0.33 a165.23 ± 69.74 a0.14 ± 0.13 b10.24 ± 4.38 b5.54 ± 2.49 a1.28 ± 0.51 ab17.20 ± 6.50 b0.07 ± 0.08 b
CLNSWC(6)0.06 ± 0.01 a1.02 ± 0.28 a7.91 ± 0.48 abc689.33 ± 157.26 b2.38 ± 1.84 a83.00 ± 23.44 bc155.10 ± 146.93 a5.71 ± 4.44 a0.78 ± 0.27 a63.40 ± 17.81 b2.78 ± 1.59 ab28.71 ± 9.43 a4.82 ± 1.53 a1.90 ± 0.61 ab38.21 ± 11.82 a0.21 ± 0.13 ab
GL(6)0.05 ± 0.01 a0.88 ± 0.10 a8.08 ± 0.59 ab731.671 ± 203.55 ab2.89 ± 2.44 a80.00 ± 29.26 bc190.98 ± 168.07 a3.33 ± 1.97 a0.85 ± 0.28 a94.23 ± 34.45 ab3.21 ± 1.52 ab30.58 ± 9.84 a4.92 ± 0.96 a1.91 ± 0.75 ab40.62 ± 11.14 a0.17 ± 0.07 ab
C(6)0.06 ± 0.06 a0.78 ± 0.28 a8.48 ± 0.30 a684.67 ± 186.27 b2.89 ± 1.87 a61.20 ± 3.44 c250.00 ± 18.68 a26.22 ± 29.61 a3.50 ± 5.15 a79.27 ± 7.05 ab5.36 ± 2.90 a38.14 ± 6.29 a4.48 ± 1.71 a2.32 ± 0.74 a50.30 ± 7.06 a0.31 ± 0.17 a
Tanzania0–15AF(6)0.21 ± 0.04 b 2.64 ± 0.48 a 6.40 ± 0.25 a908.38 ± 25.28 ab 7.44 ± 5.56 a 62.35 ± 10.23 b306.33 ± 98.08 ab 2.39 ± 1.75 a7.32 ± 2.79 a121.67 ± 28.20 a0.02 ± 0.01 a10.56 ± 3.17 a3.27 ± 0.12 ab 0.27 ± 0.18 a14.11 ± 3.15 a0.07 ± 0.01 a
CF(3)0.24 ± 0.05 b2.76 ± 0.50 a 6.42 ± 0.22 a832.67 ± 44.74 b7.40 ± 1.23 a157.33 ± 41.31 a348.67 ± 87.27 a 2.53 ± 0.97 a7.24 ± 1.65 a94.47 ± 17.59 a0.02 ± 0.00 a11.92 ± 3.17 a3.88 ± 065 a 0.13 ± 0.05 a 15.95 ± 3.86 a0.09 ± 0.02 a
CLSWC(6)0.18 ± 0.03 b2.23 ± 0.33 a 5.99 ± 0.46 a995.17 ± 89.56 a5.97 ± 3.07 a78.00 ± 18.57 b129.55 ± 97.35 b5.80 ± 3.11 a2.98 ± 2.40 ab97.05 ± 6.71 a0.02 ± 0.03 a7.89 ± 2.49 ab 2.52 ± 0.92 ab 0.17 ± 0.14 a10.61 ± 3.42 ab 0.04 ± 0.01 a
CLNSWC(6)0.22 ± 0.05 b2.61 ± 0.61 a6.47 ± 0.48 a937.33 ± 86.31 ab 8.32 ± 4.70 a81.40 ± 37.31 b279.17 ± 124.67 ab 15.74 ± 22.08 a5.85 ± 1.96 a96.78 ± 30.02 a0.02 ± 0.01 a10.53 ± 3.48 a 3.30 ± 1.18 ab 0.31 ± 0.39 a14.15 ± 5.00 a0.08 ± 0.03 a
GL(6)0.21 ± 0.05 b2.59 ± 0.43 a 6.10 ± 0.42 a963.33 ± 764.66 a6.46 ± 3.72 a73.32 ± 18.43 b159.37 ± 129.13 ab 4.49 ± 2.62 a3.82 ± 3.27 ab101.33 ± 14.60 a0.01 ± 0.01 a8.68 ± 3.21 ab3.00 ± 0.50 ab 0.13 ± 0.09 a11.82 ± 3.42 ab 0.06 ± 0.01 ab
C(6)0.83 ± 0.03 a0.93 ± 0.23 b 6.28 ± 0.52 a930.50 ± 64.20 ab4.05 ± 2.01 a60.28 ± 13.19 b160.05 ± 113.86 ab 0.44 ± 0.49 a0.47 ± 3.27 b135.02 ± 39.39 a0.02 ± 0.01 a4.77 ± 1.90 b2.05 ± 1.00 b 0.08 ± 0.10 a9.92 ± 2.71 b0.03 ± 0.00 c
15–45AF(6)0.14 ± 0.05 ab1.62 ± 0.64 ab6.59 ± 0.34 a879.50 ± 43.22 a6.46 ± 5.34 a53.97 ± 11.77 b194.90 ± 127.46 ab0.93 ± 1.08 ab3.35 ± 3.31 a133.92 ± 27.43 a0.03 ± 0.03 a8.23 ± 3.88 a2.79 ± 0.66 ab0.07 ± 0.05 a11.11 ± 4.52 a0.05 ± 0.01 a
CF(3)0.13 ± 0.01 ab1.41 ± 0.11 ab6.37 ± 0.12 a844.00 ± 10.82 a4.64 ± 0.07 a116.67 ± 9.50 a392.00 ± 124.90 a0.00 ± 0.00 b3.21 ± 1.35 a107.67 ± 5.69 a0.03 ± 0.00 a7.97 ± 0.80 a3.96 ± 0.22 a0.08 ± 0.01 a12.04 ± 1.00 a0.05 ± 0.01 a
CLSWC(6)0.12 ± 0.05 ab1.42 ± 0.62 ab6.18 ± 0.52 a949.83 ± 96.80 a4.46 ± 1.83 a54.52 ± 17.50 b48.97 ± 30.32 b1.42 ± 1.81 ab0.56 ± 0.45 aa115.00 ± 12.55 a0.04 ± 0.04 a6.63 ± 2.36 a2.39 ± 0.95 ab0.05 ± 0.02 a9.11 ± 3.27 a0.04 ± 0.02 a
CLNSWC(6)0.14 ± 0.05 ab1.60 ± 0.55 ab6.38 ± 0.40 a916.67 ± 90.96 a5.48 ± 2.34 a58.82 ± 11.23 b185.50 ± 127.83 ab2.75 ± 3.11 ab1.93 ± 1.43 a107.27 ± 35.80 a0.02 ± 0.01 a7.52 ± 2.92 a2.10 ± 0.59 b0.08 ± 0.08 a9.72 ± 3.51 a0.04 ± 0.01 a
GL(6)0.20 ± 0.04 a2.43 ± 0.60 a5.98 ± 0.48 a981.83 ± 91.02 a5.94 ± 3.57 a77.03 ± 16.88 b134.30 ± 114.82 b3.63 ± 2.22 a2.79 ± 3.54 a104.40 ± 25.07 a0.02 ± 0.03 a8.40 ± 2.56 a2.28 ± 0.86 b0.07 ± 0.06 a10.77 ± 3.37 a0.05 ± 0.01 a
C(6)0.06 ± 0.03 b0.75 ± 0.28 b6.27 ± 0.71 a962.50 ± 100.41 a3.45 ± 2.21 a57.98 ± 32.52 b110.69 ± 120.97 b0.21 ± 0.34 ab0.16 ± 0.28 a141.88 ± 53.08 a0.04 ± 0.02 a4.29 ± 2.11 a1.99 ± 1.16 b0.02 ± 0.02 a6.33 ± 3.00 a0.03 ± 0.01 a
45–100AF(6)0.05 ± 0.01 b0.71 ± 0.24 ab6.49 ± 0.59 a877.17 ± 106.93 a2.65 ± 1.14 a42.85 ± 20.11 b35.22 ± 43.35 b0.05 ± 0.12 a0.22 ± 0.27 a143.93 ± 41.59 a0.04 ± 0.01 a4.99 ± 2.17 a2.06 ± 0.75 ab0.03 ± 0.03 a7.11 ± 2.76 a0.03 ± 0.01 a
CF(3)0.07 ± 0.01 ab0.69 ± 0.08 ab6.48 ± 0.17 a888.00 ± 21.66 a2.52 ± 0.25 a91.30 ± 23.13 a233.77 ± 141.97 a0.00 ± 0.00 a0.54 ± 0.63 a142.67 ± 29.67 a0.08 ± 0.01 a5.46 ± 0.99 a3.70 ± 0.58 a0.06 ± 0.01 a9.30 ± 1.56 a0.04 ± 0.01 a
CLSWC(6)0.07 ± 0.04 ab0.91 ± 0.52 ab6.18 ± 0.47 a1000.00 ± 105.50 a3.69 ± 1.22 a43.15 ± 16.81 b21.55 ± 36.55 b0.66 ± 0.84 a0.04 ± 0.09 a122.27 ± 20.81 a0.05 ± 0.07 a5.59 ± 2.62 a2.29 ± 1.26 ab0.02 ± 0.02 a7,95 ± 3.87 a0.04 ± 0.02 a
CLNSWC(6)0.06 ± 0.02 ab0.81 ± 0.13 ab6.49 ± 0.36 a928.00 ± 13.44 a4.11 ± 1.17 a39.38 ± 9.40 b59.62 ± 78.58 b0.97 ± 1.09 a0.19 ± 0.16 a118.80 ± 25.23 a0.04 ± 0.03 a5.16 ± 0.99 a1.57 ± 0.51 a0.03 ± 0.03 a6.79 ± 1.41 a0.03 ± 0.01 a
GL(6)0.12 ± 0.05 a1.45 ± 0.63 a6.21 ± 0.40 a962.67 ± 44.13 a4.38 ± 2.66 a56.42 ± 12.34 b75.33 ± 54.92 b0.78 ± 0.95 a0.98 ± 1.21 a123.17 ± 19.46 a0.04 ± 0.03 a6.65 ± 1.87 a2.04 ± 0.81 ab0.04 ± 0.02 a8.77 ± 2.59 a0.04 ± 0.02 a
C(6)0.05 ± 0.03 b0.71 ± 0.34 b6.29 ± 0.86 a920.33 ± 178.02 a3.33 ± 2.15 a46.13 ± 18.50 b73.62 ± 96.54 b0.17 ± 0.19 a0.24 ± 0.42 a148.02 ± 56.09 a0.04 ± 0.02 a3.93 ± 2.66 a1.84 ± 1.12 b0.02 ± 0.02 a5.84 ± 3.46 a0.03 ± 0.01 a
Uganda0–15AF(6)0.22 ± 0.03 a3.07 ± 0.32 ab6.89 ± 0. 72 a937.17 ± 115.00 b3.68 ± 0. 52 b75.08 ± 9.67 b 468.83 ± 74.49 a 22.52 ± 45.35 a3.25 ± 1.25 a 89.83 ± 17.08 ab0.03 ± 0.01 a 18.21 ± 7.51 a 5.50 ± 0.79 a 0.91 ± 0.74 a 24.65 ± 8.67 a 0.10 ± 0.02 a
CF(6)0.36 ± 0.13 ab4.50 ± 1.66 a6.03 ± 0.87 abc1031.50 ± 351.99 ab 3.40 ± 1.66 ab 193.50 ± 104.62 a200.38 ± 154.89 c 14.45 ± 8.43 a3.27 ± 0.59 a 94.95 ± 53.41 ab0.13 ± 0.18 a 14.23 ± 7.23 ab 5.17 ± 2.13 a 0.61 ± 0.25 a 20.37 ± 9.25 ab 0.1140.02 a
CLSWC(6)0.20 ± 0.05 ab 2.68 ± 0.76 ab6.46 ± 0.29 ab 922.50 ± 90.68 b 4.52 ± 0.32 a100.02 ± 22.71 b 557.00 ± 62.53 a 6.08 ± 4.26 a3.24 ± 0.59 a 64.78 ± 14.50 b0.03 ± 0.02 a 12.15 ± 2.49 abc3.24 ± 0.70 ab 1.17 ± 0.93 a 16.59 ± 3.66 abc 0.09 ± 0.03 a
CLNSWC(6)0.23 ± 0.04 ab3.40 ± 0.32 ab6.32 ± 0.25 abc 969.67 ± 127.33 ab 3.36 ± 0.60 ab96.62 ± 20.24 b 398.83 ± 51.47 ab2.57 ± 0.46 a 1.88 ± 1.41 ab108.05 ± 42.91 ab 0.03 ± 0.02 a 14.41 ± 1.89 ab4.13 ± 0.42 ab 0.56 ± 0.50 a 19.13 ± 2.07 ab 0.06 ± 0.02 a
GL(6)0.32 ± 0.25 ab4.82 ± 3.47 a5.79 ± 0.52 bc 1483.33 ± 483.47 a 2.19 ± 0.83 ab 141.38 ± 55.77 ab104.23 ± 122.66 c3.63 ± 0. 88 a0.76 ± 0.35 b220.35 ± 164.81 a0.06 ± 0.03 a 7.96 ± 4.72 bc 3.11 ± 2.13 ab 0.40 ± 0.47 a 11.53 ± 6.74 bc 0.04 ± 0.02 a
C(6)0.13 ± 0.03 b1.35 ± 0.73 b 5.31 ± 0.59 c 1391.67 ± 200.00 ab 2.16 ± 0.90 b91.08 ± 29.77 b 252.66 ± 140.10 bc4.91 ± 8. 09 a0.76 ± 1.17 b223.00 ± 104.40 a0.26 ± 0.45 a 4.79 ± 3.12 c2.13 ± 0.91 b 0.58 ± 0.57 a7.73 ± 3.72 c 0.22 ± 0.27 a
15–45AF(6)013. ± 0.04 a1.76 ± 0.54 a6.59 ± 0.84 a1145.00 ± 165.62 b3.54 ± 0.54 a92.27 ± 9.40 a390.83 ± 129.82 a1.28 ± 1.03 ab0.89 ± 0.52 ab134.55 ± 50.11 bc0.05 ± 0.03 a10.17 ± 3.73 a4.19 ± 1.16 a0.82 ± 1.35 a15.23 ± 4.77 a0.06 ± 0.03 b
CF(6)0.16 ± 0.05 a2.06 ± 0.73 a5.71 ± 0.54 ab1251.67 ± 252.62 b3.37 ± 1.12 a196.43 ± 136.02 ab128.95 ± 119.65 b3.34 ± 2.34 a0.68 ± 0.32 b240.35 ± 80.11 bc0.78 ± 1.71 a3.13 ± 4.42 ab3.67 ± 1.70 a047. ± 0.39 a12.05 ± 5.82 a0.05 ± 0.02 ab
CLSWC(6)0.13 ± 0.03 a1.58 ± 0.64 a6.46 ± 0.10 a1024.17 ± 40.05 b3.84 ± 0.90 a120.40 ± 28.20 ab519.17 ± 123.35 a1.61 ± 2.50 ab1.59 ± 0.83 a101.97 ± 27.29 c0.02 ± 0.01 a9.68 ± 3.28 a3.10 ± 0.82 ab0.45 ± 0.35 a13.25 ± 3.48 a0.05 ± 0.02 ab
CLNSWC(6)0.15 ± 0.03 a2.12 ± 0.20 a6.17 ± 0.53 a1144.50 ± 186.72 b3.45 ± 0.50 a109.78 ± 13.02 ab358.00 ± 90.40 a0.54 ± 0.34 b0.69 ± 0.26 b158.50 ± 53.09 bc0.05 ± 0.02 a10.33 ± 2.78 a3.07 ± 0.31 ab0.26 ± 0.23 a13.71 ± 3.16 a0.04 ± 0.01 ab
GL(6)0.23 ± 0.24 a3.28 ± 20.95 a5.13 ± 0.33 b1863.33 ± 496.41 a1.73 ± 0.80 b103.38 ± 50.45 a33.06 ± 41.91 b1.29 ± 0.84 ab0.39 ± 0.13 b346.83 ± 172.25 a0.07 ± 0.05 a2.03 ± 2.16 b0.83 ± 1.15 c0.13 ± 0.06 a3.06 ± 3.27 b0.02 ± 0.01 ab
C(6)0.11 ± 0.09 a1.26 ± 1.24 a5.05 ± 0.62 b1471.67 ± 190.62 ab1.57 ± 0.88 b70.35 ± 22.24 b152.43 ± 106.79 b0.23 ± 0.28 b0.19 ± 0.15 b275.33 ± 72.81 ab0.05 ± 0.01 a2.76 ± 1.68 b1.47 ± 0.68 bc0.15 ± 0.09 a4.42 ± 2.22 b0.07 ± 0.04 a
45–100AF(6)0.09 ± 0.05 a 1.09 ± 0.87 ab 6.28 ± 0.74 a1215.00 ± 70.92 b 2.60 ± 0.93 ab103.15 ± 23.37 ab331.83 ± 167.95 a0.61 ± 0.53 ab0.39 ± 0.23 a170.63 ± 50.87 b0.03 ± 0.01 a5.72 ± 3.38 a2.95 ± 1.03 ab0.98 ± 2.05 a9.68 ± 4.23 a0.04 ± 0.02 a
CF(6)0.11 ± 0.03 a1.18 ± 0.36 ab 5.84 ± 1.18 a1330.83 ± 294.86 b4.58 ± 2.86 ab130.12 ± 69.57 a99.99 ± 4.02 bc1.16 ± 0.94 a0.34 ± 0.17 a172.52 ± 101.91 b2,12 ± 4.71 a5.16 ± 2.71 a4.28 ± 3.11 a0.32 ± 0.23 a11.88 ± 9.48 a0.06 ± 0.06 a
CLSWC(6)0.08 ± 0.01 a0.75 ± 0.26 ab 6.11 ± 0.50 a1118.33 ± 59.81 b2.14 ± 0.66 b118.45 ± 29.32 ab 392.00 ± 143.57 a0.11 ± 0.28 b0.49 ± 0.35 a176.67 ± 48.27 b0.02 ± 0.01 a6.01 ± 2.13 a2.83 ± 0.56 ab0.23 ± 0.10 a9.09 ± 2.16 ab0.04 ± 0.01 a
CLNSWC(6)0.09 ± 0.01 a0.96 ± 0.09 ab 5.88 ± 0.68 a1278.33 ± 234.15 b2.13 ± 0.76 b100.60 ± 21.48 ab 275.00 ± 78.22 ab0.00 ± 0.00 b0.28 ± 0.04 a253.67 ± 69.94 b0.03 ± 0.01 a5.43 ± 2.87 a2.10 ± 1.11 abc0.13 ± 0.08 a7.69 ± 3.93 b0.03 ± 0.01 a
GL(6)0.17 ± 0.17 a2.15 ± 1.68 ab5.06 ± 0.26 a1928.33 ± 335.76 a 1.50 ± 0.76 b76.25 ± 33.43 ab 18.06 ± 40.90 c0.63 ± 0.47 ab0.18 ± 0.20 a391.33 ± 136.02 b0.04 ± 0.02 a0.61 ± 0.30 a0.26 ± 0.23 c0.09 ± 0.06 a1.01 ± 0.45 ab0.01 ± 0.00 a
C(6)0.07 ± 0.02 a0.70 ± 0.22 b 5.11 ± 0.66 a 1455.00 ± 169.56 b1.14 ± 0.63 b57.25 ± 8.67 b85.67 ± 140.90 c0.15 ± 0.23 b0.24 ± 0.25 a290.33 ± 66.70 ab0.06 ± 0.03 a2.58 ± 1.83 ab1.01 ± 0.60 bc0.19 ± 0.23 a3.84 ± 2.50 ab0.06 ± 0.05 a
AF = Agroforestry; CF = Community forest; CLSWC = Cropland with soil and water conservation; CLNSWC = Cropland without soil and water conservation; GL = Grassland; C = Control. Means down the groups within a country followed by different letters are significantly different (p < 0.05) with respect to land use and soil depths; n = number of samples per land-use category.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Recha, J.W.; Olale, K.O.; Sila, A.; Ambaw, G.; Radeny, M.; Solomon, D. Ensemble Modeling on Near-Infrared Spectra as Rapid Tool for Assessment of Soil Health Indicators for Sustainable Food Production Systems. Soil Syst. 2021, 5, 69. https://doi.org/10.3390/soilsystems5040069

AMA Style

Recha JW, Olale KO, Sila A, Ambaw G, Radeny M, Solomon D. Ensemble Modeling on Near-Infrared Spectra as Rapid Tool for Assessment of Soil Health Indicators for Sustainable Food Production Systems. Soil Systems. 2021; 5(4):69. https://doi.org/10.3390/soilsystems5040069

Chicago/Turabian Style

Recha, John Walker, Kennedy O. Olale, Andrew Sila, Gebermedihin Ambaw, Maren Radeny, and Dawit Solomon. 2021. "Ensemble Modeling on Near-Infrared Spectra as Rapid Tool for Assessment of Soil Health Indicators for Sustainable Food Production Systems" Soil Systems 5, no. 4: 69. https://doi.org/10.3390/soilsystems5040069

APA Style

Recha, J. W., Olale, K. O., Sila, A., Ambaw, G., Radeny, M., & Solomon, D. (2021). Ensemble Modeling on Near-Infrared Spectra as Rapid Tool for Assessment of Soil Health Indicators for Sustainable Food Production Systems. Soil Systems, 5(4), 69. https://doi.org/10.3390/soilsystems5040069

Article Metrics

Back to TopTop