Tree-based Modeling for Large-scale Management in Agriculture: Explaining Organic Matter Content in Soil

1. Introduction

Machine learning (ML) has emerged as a crucial domain in science and technology, exerting a substantial socioeconomic–environmental influence on various aspects of human and natural systems [1,2]. ML allows us to learn from vast amounts of data and improve the predictive performance of models. However, ML often employs complex algorithms, which results in black-box models due to their intricate internal processes that are not readily interpretable [3,4,5]. Such opaqueness may lead stakeholders to overlook meaningful patterns or issues arising from hidden biases in the data. This hinders the handling of effective predictive resource management, mainly when based on large-scale ML [6]. For example, the accuracy problem of inaccurate yield mapping is well-known due to errors inherent to a high data volume and algorithms' opacity [7].

Considerable concern has been expressed about relying on opaque models that may result in decisions that are not fully comprehended or, even worse, violate ethical principles in terms of business and the environment or legal norms [1,8]. These risks are particularly relevant for critical decision making in real-life scenarios and for access to public benefits [9], for example, digitalization in agriculture [10] and terrestrial conservation [11]. This partly explains the relatively low adoption rate of current ML-based decision support systems in many areas. Land managers, government agencies, and companies that incorporate black-box ML models into their practices, products, and applications potentially face efficiency, safety, and trust issues [12]. Therefore, the lack of interpretability must be addressed, and increasing the explainability of predictive modeling for data analysis is becoming increasingly important for mitigating these unintended risks and promoting the correct application of ML models in critical domains.

In 2018, the European Parliament implemented the General Data Protection Regulation, which establish provisions regarding automated decision making [1]. These regulations aim to ensure that individuals have the right to receive "comprehensible explanations of the underlying reasoning" when automated decision-making processes are used. Additionally, in 2019, the European Union’s High-Level Expert Group on Artificial Intelligence (AI) introduced ethical guidelines for trustworthy AI, and one of the critical requirements is explainability [1]. This requirement has been incorporated into the proposed EU regulation known as the AI Act [13], which establishes standardized rules for AI, thus affecting ML as a subfield of AI. Similar but nonregulatory proposals exist for AI risk management in the U.S., such as the "Identifying Outputs of Generative Adversarial Networks Act" and the "National Artificial Intelligence Initiative Act of 2020" [14]. Overall, the consensus on the importance of developing practical explanation tools is growing. Meaningful explanations are critical for describing data, testing models, identifying potential biases, addressing risks, and fostering trust and collaboration between humans and their AI assistants. However, this remains an ongoing scientific challenge [5].

An optimum model should be highly accurate and easy to interpret [2]; however, despite the rising interest in these models, achieving both interpretable and highly accurate model outputs has presented a considerable challenge [15], particularly in response to the abovementioned concerns at the management and policy scales. Consequently, the development of various explanation methods has increased for black-box models in both academia and industry [16,17,18,19,20]. Explainable ML emerged in the late 2010s for prediction in different systems to better explain black-box models and comprehensively address diverse aspects of the food and agriculture sector [21,22]. Explainable ML seeks to enhance the interpretability of complex algorithms while still maintaining their accuracy. Prioritizing interpretable predictions is more important than solely focusing on accurate predictions when using ML models, especially for decision making [9,23].

The objective of the study was to showcase the potential of explainable ML algorithms in analyzing large-scale data. Agricultural system data and models allow us to explore the biophysical, practical, and social aspects of food production [24,25]. However, black-box machine learning often determines false relationships between components in the system, making it unsuitable for predicting and explaining [26]. It is important to address these issues in this field because agricultural production is influenced by the management decisions of growers in response to changes in our climate and environment. We specifically focused on how different tree-based ML models can uncover novel patterns of organic matter in soil from data from Korea's field cropland, compiled on a national scale. Here, we targeted soil organic matter, as soils are integral to food production and have the potential to mitigate greenhouse gas emissions [27], while the soil pools are the most vulnerable to land degradation and climate change [17,18] and are constrained by various social, economic, and political factors [28]. This approach highlights the dominant controls of soil organic matter across fields in which its distribution interacts with the environmental state and the sociocultural matrix [29,30].

2. Materials and Methods

2.1. Data compilation and processing

We compiled environmental variables representing soil, terrain, climate, and vegetation to describe diverse field conditions in South Korea. We categorized explanatory variables into five groups: 1) soil chemical properties; 2) soil map; 3) terrain; 4) vegetation; and 5) climate (Table 1).

The National Institute of Agricultural Sciences (NAS) provides data on soil chemical properties (0–0.15 m depth) for agricultural fields [31]. Specifically, the data include organic matter (g kg^-1), pH (1:5 H₂O), available P₂O₅ (mg kg^-1), available SiO₂ (mg kg^-1), exchangeable magnesium (cmol⁺ kg^-1), exchangeable potassium (cmol⁺ kg^-1), exchangeable calcium (cmol⁺ kg^-1), and electric conductivity (dS m^-1) [31]. The soil data are updated annually across fields to recommend a crop management plan and fertilizer application rates. The dataset is open-source and contains field-specific data for the last three years, accessible from the administrative division. Soil chemical data within each division are organized by location and sampling date. The Rural Development Administration (RDA) provides thematic soil maps of 30 physical properties at a 125 m resolution (1:25,000 scale) [32]; we selected topsoil (0–20 cm) texture, drainage class, erosion grade, soil order, soil structure, and parent material in this study.

We obtained the digital elevation model (DEM) data (90 m resolution), which represented a continuous topographic elevation surface, from the National Spatial Data Infrastructure (NSDI) portal [33]. We used the DEM data to derive slope, aspect, roughness, topographic position index, terrain ruggedness index, and flow direction. Hill shade was computed from slope and aspect [36], assuming sun elevation and direction (azimuth) angles of 45° and 0°, respectively. Lastly, we estimated the upslope contributing area and topographic wetness index according to Quinn et al. [37]. We used a global estimate of annual net primary production in TIFF format from 2019 to 2021 to represent vegetation from the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard NASA's Terra and Aqua satellites [35]. The spatial resolution of the data is 0.1° (approximately 11 km). To represent the current climate, we used mean annual air temperature, maximum/minimum air temperature, mean annual precipitation, solar irradiation, relative humidity, and wind speed data based on the Modified Korean-Parameter-elevation Regressions on Independent Slopes Model (MK-PRISM, version 1.2). All data were available from 2000 to 2019, except for solar irradiation (2014–2019) [34]. The MK-PRISM data were in netCDF format. All climatic variables were averaged annually and then over the entire period, whereas the sum of precipitation values was calculated and averaged over the same period.

According to the land cover data from EGIS [38], agricultural land subclasses include: 1) fields (including rearranged fields); 2) rice paddies; 3) facility cultivation; 4) orchards; and 5) pastures and nurseries. Based on the survey of arable land in 2021 [39], the total field area was 766,000 ha. We first compiled valid soil chemical data and the coordinates of the fields (n = 263,790). We then extracted the values of the climate, terrain, and soil map variables with the coordinates of the samples. We removed outliers from the complete dataset based on Cook’s distance (n = 9,584) and scaled the data for modeling. All data were compiled and processed in R [40].

2.2. Explainable Tree-based Models

We chose decision tree (DT), random forest (RF), and gradient boosting (GB) due to their ease of interpretation and proven success in handling structured datasets [41,42]. The DT algorithm is a recursive partitioning method that iteratively generates child nodes, which are further divided into pairs of nodes. DT serves as an interpretable ML model, as the decision path from the root to the leaf or terminal nodes keeps track of the features used for predictions. However, DT learning is difficult to comprehend and suffers from overfitting when the models become more complex (i.e., with a larger maximum depth). Ensemble methods, such as RF and GB, are used instead to address this issue. By incorporating multiple decision trees, these approaches produce robust and accurate models. Leveraging the diversity of individual trees, they mitigate overfitting while enhancing overall performance. RF uses feature randomization to create each tree. At each split, a random subset of features is considered to determine the best splitting point for the node. This randomness increases the diversity among trees and enhances the generalization capability of the ensemble model.

GB is an alternative ensemble learning technique in which an additive model is constructed by sequentially combining predictions from multiple decision trees, thus creating a robust predictive model. RF primarily aims to reduce variance through subset (i.e., bagging) and feature randomization. In contrast, GB focuses on diminishing bias and enhancing predictive accuracy through iterative optimization processes. Among the models, the DT algorithm is relatively simple and its interpretability is high. The last two models are composed of many trees that have been combined, which, therefore, must be supplemented with interpretable methods for understanding model behavior.

2.2.1. Model evaluation

We integrated the Shapley additive explanation (SHAP) method with the optimized ML to understand the prediction criteria of the models and the level of contribution of the individual features that strongly correlate in organic matter prediction. Feature importance provides a crucial reference for feature selection and model interpretability [43]. A higher feature importance value, compared with another, implies the greater importance of the feature for generating a prediction. The SHAP method, developed by Lundberg and Lee [44], is a method for interpreting the predictions of ML models. Shapley values quantify the contribution of each explanatory variable in each instance. The term "additive" signifies that the corresponding Shapley value for each explanatory variable of an instance can be additively combined [45]. This approach provides quantitative information to explain how individual explanatory variables either positively or negatively impact the target variable of interest in the model. The following formula represents the Shapley value:

$$\mathit{g}(\mathit{z}\mathit{\text{'}})={\mathit{\varphi }}_0+\sum _{\mathit{i}=1}^{\mathit{M}}{\mathit{\varphi }}_{\mathit{i}}\mathit{z}{\mathit{\text{'}}}_{\mathit{i}\mathit{\text{'}}}$$

where ${\varphi }_i$ is the feature attribution value (${\varphi }_i\in \mathbb{R}$), and $z{\text{'}}_i$ represents the binary variables ($z{\text{'}}_i\in {\left\{\mathrm{0,1}\right\}}^M$). When an explanatory variable is included and has a value of $z{\text{'}}_i=1$, it becomes an excluded input feature if its value is $z{\text{'}}_i=0$. M represents the total number of input features, while g explains the model. Furthermore, the importance of feature i is defined by the Shapley value:

$${\mathit{\varphi }}_{\mathit{i}}=\frac{1}{\left|\mathit{N}\right|!}{\sum }_{\mathit{S}\subseteq \left|\mathit{N}\right|\backslash \left\{\mathit{i}\right\}}\left|\mathit{S}\right|!(\left|\mathit{N}\right|-\left|\mathit{S}\right|-1)!\left[\mathit{f}\left(\mathit{S}\cup \left\{\mathit{i}\right\}\right)-\mathit{f}(\mathit{S})\right],$$

where f(S) represents the output of the ML model being explained with a set S of features, and N refers to the entire set of available features. The Shapley value or contribution of a specific feature i(${\varphi }_i$) is calculated as the average contribution across all possible permutations of the feature set [46].

Although the Shapley-based feature importance method effectively identifies important features, it does not provide insights into two crucial aspects: 1) how changes in these important features affect predictions and 2) whether a specific threshold for these key features can enhance the accuracy of predicting the target variable. To address these concerns, we employed partial dependence (PD) plotting, which enabled us to analyze how specific input features contribute to variations in the expected response of the target variable. One-dimensional PD plots illustrate how changes in the chosen independent variables impact the expected value of their dependent variable while other independent variables remain constant. Furthermore, two-dimensional PD plots are used to assess predicted values for their dependent variable as both independent variables simultaneously vary.

PD plots are a tool for comprehending the impacts of features in any machine learning model [47], such as$\widehat{f}$(x). The explanatory variables x = ($x_s$, $x_c$) can be divided into two subsets to emphasize the concept and effectively convey this relationship and visualization: $x_s$ and $x_c$. $x_s$ represents the "chosen" independent variable, and $x_c$ comprises the set of other independent variables. The function that represents PD is defined as:

$${ PD}_S\left(\mathrm{X}\right)={PD}_S\left(X_S\right)=\int \widehat{f}(X_S, X_C)dP(X_C).$$

The PD curve is estimated by calculating the average prediction when a specific variable subset, denoted as $X_S$, takes on specific values, represented as $x_S$ [47].

$$\widehat{{PD}_S}\left(x\right)=\widehat{{PD}_S}\left(x_s\right)=\frac{1}{n}\sum _{i=1}^n\widehat{f}\left(x_s,x_{ic}\right).$$

Note that, when $X_S$ is equal to 𝑋, the value of the partial dependence function ${PD}_S\left(\mathrm{x}\right)$ corresponds to the estimated output function $\widehat{f}$(x). In a hypothetical scenario where s is an empty set, such as when $X_C$ = 𝑋, this would result in:

$${PD}_{\mathrm{\varnothing }}\left(X\right)={PD}_{\mathrm{\varnothing }}=\int \widehat{f}\left(X\right)dP\left(X\right),$$

which does not depend on x and represents the constant average prediction of the model calculated as:

$$\widehat{{PD}_{\mathrm{\varnothing }}}\left(x\right) =\widehat{{ PD}_{\mathrm{\varnothing }}}=\frac{1}{n}{\sum }_{i=1}^n\widehat{f}\left(x_i\right)$$

2.2.2. Assessment statistics

To train the model, we performed five-times repeated five-fold cross validation. Grid search was employed for exhaustive searching to systematically tune the hyperparameters in ML from a predefined set of values to enhance performance and mitigate the risk of overfitting (Table 2).

To assess the performance of the models, we computed the root mean squared error (RMSE) to quantify the inaccuracy of the estimates, the mean absolute error (MAE) to quantify the errors between modeled and observed values, and the coefficient of determination (R²). The RMSE accounts for both the bias and the imprecision of the analysis. We report the mean and standard deviation of the assessment statistics from the cross validations.

3. Results

3.1. Comparison of Selected Tree-based Models

Among the models, GB exhibited the superior performance, as indicated by the highest R² of 0.59 (Figure 1). Following closely, RF demonstrated reasonable predictive ability with an R² of 0.57. Lastly, the DT model displayed the weakest performance, with an R² of 0.45. Furthermore, the RF and GB models exhibited higher prediction accuracy with fewer errors. Moreover, these two models showed a more concentrated distribution of estimates around the observations, suggesting increased consistency and reliability in their predictions compared with those of the DT model (Figure 1).

The optimal values for the hyperparameters of each model used in the study are summarized in Table 2. For the DT model, a maximum depth of six was selected to prevent overfitting caused by excessive expansion of the tree structure. The minimum number of samples per node was set to 12. In the RF model, a maximum depth of 12 and a minimum number of eight samples per leaf were chosen. Additionally, an internal node requires at least eight samples for splitting to manage complexity and minimize overfitting concerns. Our RF and GB models comprised 100 trees to ensure accurate results. However, adding more trees beyond this point yielded no further improvements.

3.2. Feature Overall Importance Analysis

The global feature explanation obtained from the SHAP method for the three tree-based models is shown in Figure 2 and Figure 3. We represent the average impact of each variable on soil organic matter based on the global feature importance, which is the mean absolute representative SHAP value for that feature over all the given samples (Figure 2). Figure 3 illustrates the impact of the features on the model's predictions. Higher SHAP values for precipitation, indicated in red, signify a large and positive average impact on predictions, suggesting a higher likelihood of increasing the prediction accuracy due to the high mean absolute SHAP values.

The features are ordered from top to bottom by their predictive importance. In this case, we found that all models consistently identified the mean annual precipitation (RN) as the primary predictor. This was closely trailed by exchangeable K content in the DT and RF models (Figure 3). Similarly, the GB model acknowledged exchangeable K as one of the top five predictors; however, it attributed relatively less importance to it (Figure 3).

3.2. Feature Partial Dependence Analysis

Subsequently, we undertook a more in-depth investigation into the impacts of precipitation and exchangeable K on soil organic matter content. This choice was driven by the results of the DT and RF analyses, which highlighted these variables as being within the top two regarding Shapley-based feature importance (Figure 2 and Figure 3). In Figure 4, we present PD plots illustrating how the three models portrayed the connections between precipitation and organic matter content; all models consistently showed a positive correlation between them. Nonetheless, we identified disparities in the strength and form of the relationships across the models. Specifically, the DT model exhibited a distinct, step-like association centered around a precipitation value of 1400 mm (Figure 4). However, both RF and GB models suggested a positive yet nonlinear relationship.

When scrutinizing the association between soil organic matter content and other variables, all three models concurred on the positive nature of these relationships for precipitation, K, and available P. Notably, only the DT model revealed a conspicuous stepwise linkage, discernible at the values of precipitation (1400 mm), K (1 cmol⁺ kg^-1), and available P (600 mg kg^-1).

Two-dimensional PD plots were employed to visually determine the interaction effects between precipitation and exchangeable K (Figure 5). All models indicated that organic matter is contingent upon both precipitation and K. The observed patterns across all models displayed distinct divisions in the patterns along the K contents around the value of 1 cmol⁺ kg^-1 and along the amount of precipitation of approximately 1400 mm. This observation strongly suggests the presence of an interaction effect between these variables. The presentation of PD plots focusing on precipitation while conditioning on varying K thresholds effectively demonstrates this interaction effect in Figure 5. A discernible pattern emerges from these plots, showing that the association with precipitation exhibits heightened strength when the level of exchangeable K surpasses 1 cmol⁺ kg^-1.

4. Discussion

Based on the results of the model performance analysis, we found that GB yielded the highest accuracy among the models evaluated, followed by RF and then DT. The precipitation variable (RN) consistently emerged as the most influential predictor in all these models under Korean conditions. Subsequently, our analysis showed that soil organic matter levels under field crops also depend on soil fertility properties, such as available K, available P, and Ca levels. In addition, analyzing the bee swarm summary plots provided a more comprehensive understanding of how other variables influenced the prediction (Figure 2). All models showed positive relationships between nutrients and soil organic matter, but the intensity and shape of these relationships varied among the models. For example, the two-dimensional PD plots illustrated the divergent interaction effects of precipitation and exchangeable K (Figure 5). The findings suggested that the relationship with organic matter is more robust when exchangeable K exceeds 1 cmol⁺ kg^-1 and annual precipitation surpasses 1400 mm. Thus, organic matter content is contingent upon both precipitation and K on a large scale. This assessment provides a comprehensive understanding of the varying importance levels attributed to key predictors using distinct modeling approaches. However, the impact of these features on the model is still not fully understood [48].

The SHAP method was also employed to assess the importance of variables at a specific local site (data not shown). The results indicated that pH was more important than RN, identified as the most influential factor in controlling global model behavior. Only the DT model revealed a conspicuous stepwise linkage, discernible at certain values of precipitation (1400 mm), available P (600 mg kg^-1), and K (1 cmol⁺ kg^-1) (Figure 4). This PD analysis underscores the nuances in how different models represent the intricate interplay between soil variables and organic matter. This insight provides valuable clarity regarding these two variables' complex interplay and collective impact on soil organic matter. Up to this point, our focus was on elucidating the overarching behavior of global models, aiming to comprehend the insights derived from the data. However, this approach falls short in explaining the nuances of local model behavior, a critical aspect in discerning the factors deemed important by the models when predicting values for specific instances.

Tree-based modeling can lead to the reliable prediction of organic matter content in soil and the identification of vital environmental factors that affect organic matter content. This knowledge is essential for managing the soil’s health and the sustainability of cropping systems, and for fine-tuning fertilizer and water use management. Moreover, when addressing agricultural challenges, especially with adopting alternative production systems [49,50], the results can support the selection of a farm location to improve not only soil organic matter and fertility conditions but also the marketability and profitability of crop harvests. For example, in Kenya, agricultural practices include the use of fertilizers, pesticides, and irrigation to enhance soil organic matter, potentially producing premium-market-priced organic products [51]. When premium prices are available for organic produce, the organic system yields significantly higher net returns than conventionally managed systems, achieving a gross margin that is 1.3 to 4.1 times higher. Additionally, intercropping various crops enhanced overall productivity and profitability [51]. Thus, these results provide insights into the economic implications of choosing an alternative farming system based on soil organic matter levels and related conditions. In enhancing soil organic matter content, several viable management approaches may offer additional means to boosting agricultural productivity, including soil amendment; yet, the economic evaluation of each of these methods has been limited, and a gap exists in terms of the comprehensive evaluation of the socioeconomic impacts, necessitating further research in this area. Petersen and Hoyle [52] modeled the benefits of soil organic carbon, mainly focusing on the increased availability of nitrogen and increased plant-available water-holding capacity (PAWC). The value of soil organic carbon is estimated to be between AUD 7.1 and 8.7 Mg^-1 ha^-1 annually [134]. This valuation includes approximately 75% for carbon sequestration and smaller proportions for productivity improvements. The enhancements in PAWC (~5%) and nitrogen replacement value (~20%) contribute to higher agricultural productivity. An increased PAWC allows for increased water retention in the soil, and increased nitrogen availability supports healthier plant growth. Both soil quality and profit result in higher-quality land having a higher market value, providing an additional incentive to adopt a soil-conserving crop production system.

Mikhailova et al. [53] evaluated the monetary value of soil organic carbon stocks in the U.S., considering various factors such as soil order, depth, and geographic region. They estimated the total value of soil organic carbon storage to range from USD 4.64 trillion to USD 23.1 trillion. This valuation highlights the critical role of its management in delivering environmental and economic benefits. Similarly, Dube et al. [54] offered a detailed financial analysis of ecosystem services from healthy soils in Vermont, highlighting benefits such as increased carbon storage (USD 19 acre^-1 year^-1 in climate mitigation), reduced phosphorus losses (USD 8 acre^-1 year^-1 in water quality), erosion control (USD 2 acre^-1 year^-1 in waterway damage reduction), and enhanced water retention (USD 2 acre^-1 year^-1 in flood damage reduction), cumulatively valued at over USD 25 million annually. This emphasizes the economic importance of soil health investment and preservation. Hacisalihoglu et al. [55] used the “market value of soil” method to calculate the cost of soil erosion, considering nutrient loss and fertilizer market prices. They estimated an average economic loss of USD 59.54 per hectare per year in pasture lands and USD 102.36 in agricultural lands in Turkey due to soil erosion. This erosion, which tends to remove the nutrient- and organic-matter-rich topsoil, diminishes soil fertility, provides nutrients, sustains structure and moisture, and affects economic values by depleting a crucial soil component. Additionally, Kane et al. [56] found that counties with higher soil organic matter levels had increased yields, lower yield losses, and lower crop insurance payout rates. A 1% increase in soil organic matter corresponded to a yield boost of 2.2 ± 0.33 Mg ha⁻¹ and a notable reduction of 36 ± 4.76% in the average proportion of liabilities paid. Sparling et al. [57] quantified the monetary value of soil organic matter in enhancing crop production in New Zealand soils. This value was determined by estimating the worth of dairy milk solids, derived from a computer simulation modeling the yield of dry pasture matter and the accumulation of organic matter. The findings revealed that soils with lower organic matter levels yielded between 8.5 to 47.7 kg fewer milk solids per hectare annually, translating to a financial impact of NZD 27 to NZD 150 per hectare. Over recovery periods of 36, 90, and 125 years, the cumulative loss per hectare at Pukekohe, due to reduced productivity, was estimated at NZD 1239 with a 3.5% discount factor and NZD 772 with a 10% discount.

Finally, Fan et al. [58] showed that field practices with varying organic matter inputs could affect the total ecosystem service valuation in organic cereal crop production systems. This impact likely stems from altered soil properties due to long-term diverse field management. The authors estimated the economic value of ES in these systems under different management strategies, finding that the economic values ranged from USD 1492 to USD 1969 per hectare per year. Reyes and Elias [59] demonstrated that drought and excess precipitation were the primary causes of crop losses in the U.S. from 2001 to 2016, leading to over USD 440 billion in economic damage. These studies emphasize the importance of understanding and predicting environmental factors in agriculture. Therefore, tree-based modeling for predicting organic matter in soil and identifying key factors can help improve soil health and contribute to more sustainable and profitable agricultural practices.

5. Conclusions

We addressed the lack of interpretability of ML modeling to gain more insight into agricultural management. Specifically, we selected three tree-based models (decision tree, random forest, and gradient boosting) and focused on soil organic matter, which can provide multiple benefits as a valuable resource but represents challenges commonly faced in agriculture. We found that tree-based explainable ML enables the reliable prediction of the organic matter content in soil and identifies vital environmental factors that are relevant to organic matter management. This knowledge is essential for managing soil health and the sustainability of cropping systems and for fine-tuning fertilizer and water use management. We further examined the potential interactions between selected variables within the top four as determined by the Shapley-based feature importance. Based on our results, incorporating deep learning into our analytical framework would provide considerable benefit. Deep learning has demonstrated remarkable capabilities in handling complex data patterns and may provide enhanced insights when applied to agricultural science. Our prior emphasis on overall feature importance analysis should complement the value of local study. Whereas overall feature importance analysis provides valuable insights at the global level, local analysis allows for a more granular understanding of how specific features impact the system under investigation. This perspective can lead to more precise and actionable recommendations. Lastly, employing feature elimination represents a meaningful avenue for analysis. This approach can aid with dimensionality reduction and identifying critical features, potentially simplifying the model while preserving its predictive power. Such streamlining can enhance model interpretability and efficiency, which are paramount in agricultural science.