- Department of Sustainable Land Management, SAPD, The School of Agriculture, Policy, and Development, University of Reading, Reading, United Kingdom
Conventional farming poses threats to sustainable agriculture in growing food demands and increasing flooding risks. This research introduces a Bayesian Belief Network (BBN) to address these concerns. The model explores tillage adaptation for flood management in soils with varying organic carbon (OC) contents for winter wheat production. Three real soils, emphasizing texture and soil water properties, were sourced from the NETMAP soilscape of the Pang catchment area in Berkshire, United Kingdom. Modified with OC content at four levels (1, 3, 5, 7%), they were modeled alongside relevant variables in a BBN. The Decision Support System for Agrotechnology Transfer (DSSAT) simulated datasets across 48 cropping seasons to parameterize the BBN. The study compared tillage effects on wheat yield, surface runoff, and GHG-CO2 emissions, categorizing model parameters (from lower to higher bands) based on statistical data distribution. Results revealed that NT outperformed CT in the highest parametric category, comparing probabilistic estimates with reduced GHG-CO2 emissions from “7.34 to 7.31%” and cumulative runoff from “8.52 to 8.50%,” while yield increased from “7.46 to 7.56%.” Conversely, CT exhibited increased emissions from “7.34 to 7.36%” and cumulative runoff from “8.52 to 8.55%,” along with reduced yield from “7.46 to 7.35%.” The BBN model effectively captured uncertainties, offering posterior probability distributions reflecting conditional relationships across variables and offered decision choice for NT favoring soil carbon stocks in winter wheat (highest among soils “NT.OC-7%PDPG8,” e.g., 286,634 kg/ha) over CT (lowest in “CT.OC-3.9%PDPG8,” e.g., 5,894 kg/ha). On average, NT released minimum GHG- CO2 emissions to “3,985 kgCO2eqv/ha,” while CT emitted “7,415 kgCO2eqv/ha.” Conversely, NT emitted “8,747 kgCO2eqv/ha” for maximum emissions, while CT emitted “15,356 kgCO2eqv/ha.” NT resulted in lower surface runoff against CT in all soils and limits runoff generations naturally for flood alleviation with the potential for customized improvement. The study recommends the model for extensive assessments of various spatiotemporal conditions. The research findings align with sustainable development goals, e.g., SDG12 and SDG13 for responsible production and climate actions, respectively, as defined by the Agriculture and Food Organization of the United Nations.
1 Introduction
Climate change and sustainability challenges have appeared as substantial threats affecting numerous aspects of life, including agriculture (Whitmee et al., 2015; Cramer et al., 2018). Climatic factors are essential as unpredictable weather events directly influence farming practices. Excessive rainfall or drought conditions can intensely affect farming actions and outcomes. For instance, prolonged drought leads to water shortage, crop yield limitations, and soil erosion. Contrarily, heavy rainfalls lead to flash flooding, delay crop planting and harvesting, and impact crop yields and losses. Excessive moisture levels affect crop quality and attract insect pest attacks and disease infestations. Similarly, extreme temperatures influence key plant development stages, affecting crop growth and yield (Nnadi et al., 2013; Rial-Lovera et al., 2017). Human-induced activities are recognized as significant contributors to climate change, often resulting from unsustainable environmental actions. These actions include deforestation, land-use practices, greenhouse gas emissions, industrial activities, burning fossil fuels, and waste management (Retnowati et al., 2014). Consequently, making the right decisions and implementing appropriate methods have become pertinent, as the ramifications will have far-reaching implications for future generations (Adger et al., 2005). The solution relies on sustainability-based actions (Gibson, 2001). Identifying and adopting the right choices among alternatives can lead to embracing sustainable solutions. The potential exists to consider one decision among others to serve the wider ambit (Walker, 2014).
The primary factors in farming are local weather conditions, soil types, and soil organic carbon levels (Mkonda and He, 2023). There is little room and limited choices for these factors to play around beyond certain points. For instance, geographical & seasonal weather conditions and extreme events can limit farming choices. However, local soilscapes exhibiting diverse health conditions should prioritize enhancing and sustaining their vigor to withstand climatic pressures and ensure sustainable yield. Hence, this study investigated soilscapes to extract real soil attributes as fundamental components (Field, 2012; Vermeulen et al., 2012). Soil organic carbon (SOC) content levels are the main indicator of the soil health in the local soils for performing farming functions, which can help build good soil structure, delivering effective crop production in return (Ryan et al., 2009). SOC serves as a reservoir of essential nutrients and is pivotal in supporting nutrient availability through ascertaining soil nutrient cycling for nitrogen, phosphorus, and sulfur. Soils with higher SOC levels improve soil fertility and reduce reliance on external inputs such as synthetic fertilizers. Hence, the SOC levels have a compelling influence on crop productivity. Higher SOC levels enhance soil water-holding capacity, which helps retain moisture and resist drought stress by helping crops withstand water scarcity periods and maintain moisture levels essential for crop growth and development (Carter, 2002). Understanding the appropriate carbon-soil organic matter (C-SOM) level is challenging, which must focus on achieving environment friendly soil carbon stocks and sustainable productivity. Testing different SOC levels in real-life situations may face spatiotemporal limitations. It can make the exploration costly and ineffective, requiring resource-intensive, prolonged field experimentations (Sierra et al., 2015; Searchinger et al., 2018; Murindangabo et al., 2023).
Mathematical modeling techniques are applied to perform this task most efficiently because the scientific approach of using modeling tools can offer scenario-based choices among alternatives to evaluate the impact of one decision over the others (Collins et al., 2013). This approach is efficient and flexible in making necessary adaptations according to the situation. Among several modeling techniques, the Bayesian modeling approach provides a comprehensive and flexible framework which can represent complex relations and diverse systems, such as the influence of Climate change and sustainable production challenges in current farming practices. Heavy rainfall events can turn into flash flooding, and extreme temperature or drought situations can cause water stress in crop growth and development. Adapting appropriate farming practices, such as tillage preference, maintaining SOC levels, and limiting GHG emissions, becomes extremely important. Bayesian networks can capture uncertainties among complicated interactions between climatology, agronomy, ecology, ecosystem dynamics and socio-economic factors to enhance resilience and sustainability. Moreover, BBN can integrate multi-domain knowledge to support decision-making choosing among alternatives (Collins et al., 2013; Stockmann et al., 2013; Parton et al., 2020; Heidari et al., 2021).
The Bayesian Belief Network (BBN) modeling technique’s rationale is its ability to handle uncertainty, integrate prior knowledge, and provide probabilistic inference. Bayesian methods offer a flexible framework for analyzing complex systems, allowing for integrating various sources of information and quantifying uncertainty in model predictions. Additionally, the BBN facilitates updating beliefs as new data becomes available, enabling iterative refinement of models and decision-making processes. Bayesian network modeling is one approach that offers the flexibility to choose several relevant factors as variables of interest to include in building a BBN model structure (Natcvetova, 2021). The BBN model can represent common and scientific interest phenomena with underlying conditional associations. This feature enables the BBN model to capture the hidden uncertainty among the complex interactions (Dal Ferro et al., 2018). These interfaces offer conditional probabilities to interrelate based on the intensity of dependencies backed by causality or conditional relationships among interacting variables. This phenomenon applies under the parent and child relationship, where a parent is independent, and a child is dependent on a parent. A variable is called a node connected to another node where the parent node is the causal or independent variable, and the child node is the dependent or effect node. A parent node connects its child node through a link arrow directed from the parent to the child node to avoid having a complete cycle. This structure develops into a network model structure, which means a directed network graph without any cycle and is also known as a directed acyclic graph (DAG) (Donald, 2011; Gopal et al., 2012).
Determining conditional relationships is also very important in this modeling process. Several approaches are used in the BBN modeling frameworks to include such variables as part of a network structure. Manual, data-driven, and hybrid approaches are known methods for constructing a BBN model using expert elicitation, data-based learning, and a combination of the earlier two approaches. There are different pros and cons for each approach. Most have limitations in reaching the sources due to Spatio-temporal factors (Neil et al., 2000; Heckerman, 2008; Koski and Noble, 2012). However, modelers or model designers can opt for one or the other approaches to identify the vital variables and their pertinent interactions for building BBN models. For instance, various modeling tools use efficient algorithms representing certain phenomena or scientific systems following ascertained relations among various variables. This study pursued established interactions for developing a new BBN model and found an opportunity to evaluate the variables’ performance with new data inputs. Parts of these could be helpful to utilize where access to limited datasets becomes a hindrance (De Oliveira and Gaudio, 2000; Pearl, 2011). The technique is also helpful where datasets are highly demanding to explore long-term impacts for predictive inferences and access to real-time empirical data is limited for other reasons. These tools can help produce synthetic datasets by executing simulations based on similar experimental setups induced with available minimum dataset applications (Drury et al., 2017; Groth et al., 2019). The hierarchical model with a process-based framework may facilitate this goal.
Process Based Models (PBMs) have been used for predicting yield and environmental regulations of plant physiological processes, which increased the use of data acquisition through sensors and technologies (Sands et al., 2000). This trend amplified the use of machine learning (ML) tools, which became popular with time. Several techniques are introduced combining different tools for attaining high prediction accuracy, such as PBM and ML were combined to develop knowledge- and data-driven. This approach presented a combination of simulation models with data science in agriculture, focusing on genetic and physiological processes involved in food production, climate change mitigation and sustainability (Hailegnaw et al., 2024). A team of researchers also introduced a comparable approach for integrating machine learning and empirical evapotranspiration modeling to explore implications for agricultural water management using the Decision Support System for Agrotechnological Transfer (DSSAT). They compared three empirical models (Hargreaves Samani – HS, Priestly Taylor – PT, and Turc - TU) with three machine learning models (Multiple linear regression – LR, Random Forest – RF, and Artificial Neural Network – NN). They found that machine-learning models outperformed the empirical models. They also informed that site specific model calibration is extremely significant for higher predictive accuracy (Korb, 2009). The potential to handle target variables and crucial processes and reflect their responses in a mechanistic model (DSSAT) is an extremely important point.
The rationale for using DSSAT to generate simulated datasets is its reputation as a robust and widely used agricultural simulation model. DSSAT integrates various modules such as crop growth, soil processes, and climate data to simulate agricultural systems precisely. Using DSSAT, researchers can simulate diverse agricultural settings under different conditions, providing valuable insights into crop performance, soil dynamics, and environmental impacts (Jones et al., 2003). This approach allows for confidence in simulated data acquisition when access to real-world data is limited. A team of scientists used the soil management data module of the DSSAT model for analyzing various types of soils with contents for their crop suitability and prediction yields by observing various input parameters such as soil, weather, rainfall, etc. (Zhang et al., 2023). Hence, the integration of synthetic datasets is targeted through executing a process-based model for data simulations using the DSSAT tool. And prototype model in the DSSAT requires a minimum level of empirical data requirement to calibrate and execute the model, which delivers the output responses of the pertinent variables aimed at the BBN model for this study. A research scientist described that the BBN models are the new generation of probabilistic systems which capture uncertainty by modeling the physical, biological and social systems (Yu et al., 2004; Cheema and Singh, 2021). Combining DSSAT with BBNs allows researchers to leverage the detailed process-based simulations provided by DSSAT while incorporating uncertainties and prior knowledge into the modeling framework using BBNs. This integration enables more robust decision-making in agricultural systems by providing probabilistic assessments of crop performance, soil dynamics, and environmental impacts while accounting for uncertainties inherent in the modeling process (Finley et al., 2011).
In rainfed agricultural systems, farming has several other associated challenges, not limited to flooding risk in the context of climate change and sustainability (Shiferaw et al., 2009; Lashford et al., 2022). Member states of the Food and Agriculture Organization (FAO) of the United Nations Organization have the understanding to achieve their targets for sustainable development goals (SDGs). In farming, SDG12 for responsible production and SDG13 for climate action measures are the most relevant targets to focus on achieving sustainability in targeted areas (Louman et al., 2019; Hales and Birdthistle, 2023; Reijneveld et al., 2023). These aspects have become extremely important in farming to consider the right choices and activities which could not endanger the environment and must not result in detrimental outcomes for others. Responsible production approaches, in combination with climate action measures, could serve this purpose (Beerling et al., 2018). For instance, whether land tilling methods (conventional or no-tillage) reduce surface runoff to facilitate flood alleviation during wet seasons without compromising the main purpose of attaining sustainable crop yield?
Tilling systems in winter wheat production are the most significant area of attention in the context of climate change and sustainability in this research study. However, tilling systems can substantially influence soil structural changes, health, carbon sequestration, moisture conservation, climate-related risks and sustainable land management in dry summers and wet winters (Hermle et al., 2008). In Berkshire, United Kingdom, farming has local challenges of flash flooding and groundwater table issues. In the winter of 2000/2001, the Pang and Lambourn catchments in Berkshire suffered from groundwater flooding, and the water table reached the land surface and produced long-duration surface flows (Hughes et al., 2011). Catchment-based, natural flood management (NFM) approaches are vital for formulating solutions around other farming objectives. Adaptation of appropriate tillage systems between no-till (NT) or conventional tillage (CT) for winter wheat can deliver in the face of Climate change and sustainability challenges. NT represents limiting the interface between soil and farm machinery or tilling implements. This practice is also known as zero tillage, conservational or non-inversion tillage. At the same time, CT presents conventional methods of tilling soils with frequent, intense tilling with regular interfaces between soils and farm machinery/ tilling implements for agronomic practices for crop establishment for production. This practice is also known as traditional inversion tillage. This study hypothesized that tillage preference could be an NFM strategy and reduce GHG-CO2 emissions to deliver sustainable crop yield. However, local soilscapes, which have diverse soil compositions and proportions of sand, silt, gravel, bulk density, and SOC content levels, remained the key focus of the study (Broadmeadow et al., 2023).
Based on the above background, this study focuses on achieving the following.
• Comparing tilling preference as an NFM strategy in winter wheat production for climate resilience.
• Quantifying the uncertainty among variables with complex interactions using long-term simulations.
• Evaluating tilling impacts for wheat production, surface runoff and greenhouse gas carbon dioxide (GHG-CO2) emissions against soil carbon variability.
2 Methods
2.1 Case study and data acquisition
The study site is associated with the River Pang adjacent to the Kennet Catchment. Site selection criteria are mainly linked with a project partner, e.g., Landwise NFM Project. They held the license for data acquisition for the catchment. This research was a part of modeling work packages for exploring NFM strategies for local farm production. Moreover, the site is near the University of Reading, Berkshire, United Kingdom. The University of Reading Weather Observatory provided the weather data used in this study. Access to resources such as data acquisition remained the key consideration for catchment selection for this research. The Pang River flows through the Berkshire Downs to join the River Thames at Pangbourne. The ground water table is the main fed water, and a little water from tributaries drains the cover of tertiary clay and sand in the lower catchment. The Pang has some challenges due to physical modifications, water quality, low flow, flooding and invasive species. Physical modifications are not limited to building structures such as culverts, bridges, tracks, infrastructure, etc., next to watercourses and water bodies. They lead to the loss of natural processes hampering the water environment (Pang River, 2020).
The Environment Agency monitors water bodies and catchment areas for their ecological, physio-chemical quality elements, hydro morphological supporting elements, chemical, priority substances and other pollutants. They classified the ecological status of the Pang as moderate and the water body type as a river, and the Pang water body administered that. The length of the river is 36.57 km2, and the catchment area is 170.531 km2 (Pang Water Body, 2024). There is no declared protected status. However, small sites of specific scientific interest are in the wider catchment, including Sulham, Tidmarsh Woods, Meadows, Coombe Wood, and Frilsham. The Pang Valley has a history of fluvial and groundwater flooding. The risk of flooding from surface water extent is about 0.1 percent annually (About the Pang, 2019). Local farming revolves around farming objectives and their interactions with businesses and lifestyles. The main land cover plus crops include grasses, wheat, oilseed rape, barley, maize, potato, beans, potatoes, etc. Winter wheat is a major cereal crop widely cultivated in arable farming systems in the catchment (Pearce, 1952; Perkin and Rehman, 1994).
This study utilized the Decision Support System for Agrotechnology Transfer (DSSAT), a software application that delivers dynamic crop growth simulation models for different crops. This tool evolved with multiple functions, incorporating weather, soil, crop management, and sample datasets for several crop models. The crop simulation models within DSSAT simulate crop development, growth, and yield based on the interaction of soil–plant-atmosphere dynamics (Jones et al., 2003). Employing three core modules, namely the Sbuild-module (for soil data file input), Weatherman-module (for climate data file input), and the Xbuild-module (for crop management data input), this tool facilitated wide-ranging crop growth modeling research (Decision Support System for Agrotechnology Transfer, n.d.). Hence, this research explored the DSSAT tool, found the most appropriate prototype experimental model for winter wheat and incorporated the available datasets (such as soils, winter wheat crop management, and climatic weather information) from the catchment area.
This research applied sample data features from the local soilscape, which describe the spatial variation and diversity of soils in a landscape. This data includes a range of soil types, textural composition, physical properties, and key features in a particular area. Various soil mapping and classification techniques delineate soilscapes and characterize spatial distribution and variability across landscapes (Farewell et al., 2011). Supplementary Table S1 shows the soilscapes explored for the Pang catchment. Soilscapes range from medium soils, deep clays/ silty, chalk and limestone soils with varying organic carbon contents. Important catchment features are highlighted in the map, shown in Supplementary Figure S1. This study used key features of data samples, which include soil features from soilscape, weather indices from climate, and winter wheat crop phenology from a prototype experiment. Data acquisition for three soils from local soilscape and weather indices of the Pang catchment area were accomplished through the relevant sources (NETMAP soilscape data) by the NFM Landwise Project and the University of Reading, United Kingdom Atmospheric Observatory, respectively. Three soilscape PDPG6, PDPG7 and PDPG8, containing key soil features, characteristics and spatial distribution, are summarized in Supplementary Table S1. These are not the entire soilscapes of the catchment but represent a considerable land in the Pang catchment area. This modeling research used detailed variability of real soils for given zones depth/ layers, soil texture, composition, SOC, bulk density and saturated hydraulic conductivity, lower limit, drained upper limit and saturated water contents are highlighted in the brown colored section of the Supplementary Table S2 for PDPG6, Supplementary Table S3 for PDPG7, and Supplementary Table S4 PDPG8, respectively. The “Weather indices” data file is available in the Supplementary material; the soilscape information is described in Supplementary Table S1, with experimental details in Supplementary Table S5 (Jamagne and King, 2002; Rameshwaram et al., 2023; Crop Cover Data Source, 2024; Geology Data Source, 2024; Land Cover Data Source, 2024; Soilscapes Data Source, 2024; Weather Data Source, 2024).
The study incorporated soils as independent variables, utilizing real soil samples with original organic carbon content from various soilscapes. Classification and characterization details are presented in Supplementary Table S1; Supplementary Figure S1. The study applied the “RB209” soil classification approach for this project. Other independent variables include weather data of seasonal rainfall (mm) and temperature (C°) from the Reading Atmospheric Observatory (Weather Data Source, 2024), as this research compared no-tillage versus conventional tillage practices in winter wheat crop production with tillage preference as treatment under a frequentist approach (Supplementary material S2) which represented as the main study input variable (independent variable) under the BBN modeling approach.
Four model soils were synthesized based on the three real soils mentioned above, with varying organic carbon contents in topsoil (1, 3, 5, 7%) (Table 1). The variations in organic carbon content within the modeled soils were adjusted in proportion to their existing levels in the real soil relative to the mineral particles. These variations were implemented using the Sbuild-module of the Decision Support System for Agrotechnology Transfer (DSSAT) version 4.8, with each layer adjusted accordingly (Decision Support System for Agrotechnology Transfer, n.d.). Fifteen soil profiles, representing each real soil from a soilscape along with its four modeled counterparts, were generated using the Sbuild-module of DSSAT. Refer to Supplementary Table S1 and Supplementary Tables S2–S4 for more information. Detailed mechanisms for synthesizing modeled soils are described as “Keynotes” under Supplementary Table S4. The Sbuild-module of the DSSAT tool applied the below regression equation on input data of soil texture to estimate the carbon in the stable organic matter based on the relationship developed by po Adiku (Porter et al., 2010; Prout et al., 2021). The Sbuild shell of the DSSAT tool calculated the bulk density of each layer of modeled soil.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-t001.jpg)
Table 1. A schematic plan comparing no-till and conventional tillage methods was designed for experimenting three real soils from the soilscape, each synthesized into four modeled soils with top layer organic carbon content levels of 1, 3, 5, and 7% using the S-module of the DSSAT tool.
The Cropping System Model (CSM) was employed within DSSAT v4.8 and contains various crop models to have simulated data for longer periods (Jones et al., 2003; Li et al., 2015). The crop models use the CERES models for wheat, and various controls & management scenarios are employed within the shell to simulate crop growth. Hence, 15 soils each for No-till and Conventional till (CT) executed in the CSM of DSSAT for having 1,440 simulations over 48 years from 1974 to 2022 and a summary of experiments is described in Table 2. The rationale for conducting simulations over an extended period is not only linked with the evaluation of historical data for improved forecasting of the variable of interest but also encompasses the provision for accessing the model with training and testing datasets. This method helped to find the responses of soils demonstrating local soilscapes containing varied carbon contents to deliver influence on wheat yield, surface runoff and GHG-CO2 emissions through related variables over 48 years in No-till (NT) and conventional till (CT) management applications using the Cropping System Model (CSM) of the DSSAT (Zhao et al., 2024).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-t002.jpg)
Table 2. The DSSAT model was used to generate synthetic datasets for 30 individual soils, e.g., one real (in orange) with four modeled (pink) soils from each of 3 real soils from local soilscapes, executing distinctive sets of simulations each over 48 years from 1974 to 2022.
2.2 Synthesis of Bayesian network model structure
The BBN model structure was synthesized by retrieving evidence reported in the published literature and practical applications available and employed in comparable software. Summary of modules for variables interactions applied in the DSSAT and the list of variables with their relationships and evidence from the literature for BBN model structure are highlighted in Supplementary Tables S6, S7. The model structure for variable relationships was also confirmed through the source manuals available with the software, e.g., the DSSAT and AquaCrop. Key considerations revolved around ascertaining and validating the model assumptions as defined and identified in the relevant sections of this manuscript (Waffa and Benoit, 2015).
This study used the Netica software interface to build a graphical network structure and represented variables as nodes and their relationships as a link or arcs representing their connection. The identified conditional (ideally causal) relationships evident through the published literature determined the direction of the link or arc (Supplementary Table S7). A single relationship between two variables helped develop a simple BBN model structure foundation. This rationale is followed for every single relationship to build a Directed Acyclic Graph (DAG), which delivers a model structure having all interrelated variables connected in a network, as shown in Figure 1 (Ni et al., 2011). Additionally, the study proved the structural morphology of the BBN model, analyzed by the Netica software, which demonstrated key interactions and is shown in Supplementary Table S9. These interactions depicted their probabilistic estimates and performance accuracy, delivering the desired model outcomes. The model was evaluated through key performance indicators outlined for model outputs in the following sections.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-g001.jpg)
Figure 1. The BBN model shows 11 nodes (variables), 17 links (relationships) and their belief based on parametric learning through simulated datasets (75%) generated in the CERES-Wheat model of the DSSAT.
This study selected the output (dependent) variables based on the study objectives for assessing tillage preferences toward winter wheat production, as tillage practices are practiced to prepare land for crop establishment through farm management to achieve crop yield. Hence, crop yield suits a major output (dependent) variable (Jug et al., 2011). Surface runoff is another vital variable that can measure flood risk reduction dependent upon tilling choices to facilitate or resist the surface water flow during cropping seasons. For instance, decreased runoff facilitates flood alleviation while tilling wheat if no-till or reduced tilling preference is opted. This approach will cover the land through vegetation and limit the interface between soil cover and farm machinery. Hence, surface runoff serves as another vital output (dependent) variable. Wheat tilling on bare land or cover (such as crop residue or mulches on land) can create a difference in resisting surface runoff during wet winters (DeLaune and Sij, 2012). Last but not least, greenhouse gas emissions for carbon dioxide (GHG-CO2) are another major output (dependent) variable because tillage preferences are linked with the disturbance of soil aggregates, impacting the decomposition rate of soil organic matter and carbon sequestration, etc. Hence, GHG-CO2 emissions were considered a vital output (dependent) variable based on tillage preferences (Galic et al., 2019).
To quantify output (dependent) variables in response to input (independent) variables of rainfall, temperature, soils from soilscapes, and tillage systems, this research explored the output variables response accessible through the DSSAT cropping system model (CSM). This study included the most applicable variables with conditional or causal influence among conditional interactions. For instance, the abovementioned parameters are focused as input variables in the DSSAT shell. The study investigated output variables of biomass, yield, precipitation interception, available water space, runoff, carbon-soil organic matter, and GHG-CO2 emissions. The investigations also explored the literature to confirm such key relationships. For instance, crop biomass accumulation is sensitive to weather indices during reproductive, flowering and grain-filling stages. A research study validated this relationship using the AquaCrop Model, which underwent Bayesian calibration involving meteorological, soil, crop, and management parameters (Zhang et al., 2018, 2019). Hence, these parameters (biomass and weather indices) were shortlisted and categorized in the BBN structure (Table 3).
The BBN model was developed based on a few assumptions, such as weather indices with SOC influencing wheat crop yield (Droste et al., 2020). Similarly, different tillage practices influence available water space and change in SOC contents over time (Sharma et al., 2016). Temperature intensity not only affects crop development & growth rate for biomass and influences GHG-CO2 emissions (Tu and Li, 2017). Rainfall affects biomass production through water stress in drought conditions and can generate surface runoff where extreme rainfall occurs. Higher organic carbon content levels in soils increase soil fertility to generate higher crop yields, but lower OC levels affect the soil’s water-holding capacity. Soil textural compositions with higher rocks limit available water space in soils to impact surface runoff. Higher biomass intercepts precipitations to reduce surface runoff and deliver higher crop yield (Refer to details in Supplementary Tables S6, S7). Hence, the most pertinent variables were included in the BBN.
The BBN model included 11 nodes (variables), 17 links (relationships), and 2,800 conditional probabilities in the final BBN model structure, as shown in Figure 1. The structure of the BBN model is informed by evidence from published literature regarding variable relationships and their practical applications, as reported and utilized in other software. For example, a few software utilize crop growth simulation methods in hierarchical models, expressing these relationships through mathematical equations in the respective tools (Ahmed et al., 2020). These software functions on their fixed models for executing crop simulations. The DSSAT software uses the Crop Environment Resource Synthesis (CERES-Wheat) model to simulate wheat crop growth to climate, soil, and management, etc. and uses the default module, namely CERES-Godwin for soil organic matter and CENTURY model, to simulate SOC dynamics (Ahmed et al., 2023). A few modules and sub-modules for Actual CO2, Mauna Loa, Hawaii (Keeling curve) are used for daily records of global atmospheric carbon dioxide concentration and others for soil conservation service for infiltration, modified soil profile for soil layer distribution, etc. (Keeling, 1986). Tillage practices impact soil carbon dynamics by altering soil structure, aeration, and microbial actions, which induce CO2 emissions. The variable interface for the primary and sub-module in DSSAT is summarized in Supplementary Table S6. Then, the model executes to generate simulations after making an experimental set-up with all modules fed with the required input datasets and files.
In DSSAT, crop models can simulate single cropping and subsequent cropping systems. The soil water balance is based on a formulation by Ritchie methods, applying the concept of a drained upper limit (DUL) and drained lower limit (LL) for available soil water. The water accounting procedure for each layer in the soil profile applies to the DSSAT by Porter and Ritchie methods. Water movement from an upper layer cascades to lower layers, resembling a series of linear reservoirs. The difference between rainfall and runoff calculates the infiltration. The drainage process between layers occurs if soil water in the layer exceeds its water-holding capacity. Root water uptake drives the upward flow due to transpiration and soil evaporation. Available soil water determines the potential root water extraction and the root length density of each layer in the soil (Ritchie et al., 1984; Parton et al., 1994; Ritchie et al., 1998; Suleiman and Ritchie, 2004; Boote et al., 2008; Porter et al., 2010; Asseng et al., 2013; Howe, 2015; Khosravi et al., 2022). The study included a parameter, “available water space,” and calculated that using figures of total available water in soil profile minus field capacity in entire simulations for all seasons (Cassel and Nielsen, 1986). This study also explored the DSSAT and Aqua Crop modeling tools manuals, confirming the underlying mechanisms and the variables’ conditional relations in delivering simulations and model outputs. Some models exhibit identical relationships between variables to give respective model outputs. For instance, the simulation of crop yield relies on biomass. The conditional relationships were identified and confirmed through literature and recorded in Supplementary Table S7. The next step was to execute the DSSAT model.
Data synthesis was achieved through the DSSAT tool that generated numerous simulations executed for 30 individual soils, each for 48 years, and a summary is described in Table 2. The DSSAT tool generated multiple variable responses based on its predefined algorithms for modular and sub-modular functions. However, this research explored only specified variables identified for their potential responses in comparing tillage preference choice in relation to wheat yield, surface runoff, and GHG-CO2 emissions. The main user interface of the DSSAT tool was used to run the adapted model for simulations through window command selection. An information alert, “Simulations are completed,” popped up shortly with an option to click the next tab to access the analysis if the process is completed successfully. Alternatively, error messages and alerts may appear to navigate the inaccuracy for model calibrations at the appropriate modules, such as Xbuild, Sbuild, Weatherman, etc., and run the program again to complete the process successfully. A flowchart illustrating the process of generating simulated datasets in the DSSAT main interface, resulting in the creation of DSSAT OUTPUT files, is included as an annexure to Supplementary Table S6. This DSSAT data served as input data for the parametrization of the BBN (Jones et al., 2011).
Hence, the DSSAT served as a data simulator only, which provided a range of variable responses to soil attributes, weather, and crop management practices depicted in the Pang catchment area. Only pertinent variables and their responses were accessed following the DSSAT simulation compilations. This phase retrieved the system-generated output files (GHG.OUT, PlantGr2.OUT, PlantGro.OUT, SoilOrg.OUT, SoilWat.OUT, Weather.OUT) using the analysis tab and GBuild module (a plotting tool for data visualization). These output files contained the simulated datasets for all variables in the BBN model. This approach resulted in retrieving a total of 1,440 simulations. The acquired datasets were then formatted into a specified case file to parametrize the BBN.
Moreover, this study utilized the Netica software and its user interface to construct and parametrize the BBN model. The structural construction started from one variable (independent) as a parent node to the second (dependent) as a child and connected them based on their conditional relationship. All parent nodes are connected to their child nodes through an arrow link directing a conditional relationship from the parent to the child nodes by avoiding having a cycle in the network. This structure developed into a directed acyclic graph (DAG), a BBN model. All 11 nodes with 17 links were connected in a BBN and were now ready to parametrize using simulated datasets. The specified case file containing (the DSSAT outputs for specified variables) is used to parametrize the BBN model. The study applied the “Holdout Method” for model training and testing for validation. Therefore, seventy-five (75%) data were used for training the model by parametrization, as exemplified in the Supplementary Table S11. The assessment used this to populate the BBN model against unseen holdout datasets for the model, and the remaining twenty-five (25%) were then used to validate the BBN model by testing. The case files for the Holdout method’s application through 25% of datasets are available as specimens in Supplementary Table S12. The study also compared model validation results against error rates among “K-Fold Cross Validation” and “Train Test” Methods (Yadav and Shukla, 2016). The summary is highlighted in Table 4.
The parametric details of the BBN model depend upon the relationships between parent and child nodes and their states defined for each variable in the network and their corresponding datasets (Supplementary Table S8). This aspect depends upon the nature and quality of datasets used to parametrize the BBN model. This study used synthetic datasets generated by the DSSAT tool, which delivered variable output responses following specified functions of certain phenomena, such as the pedotransfer function for soil water dynamics and CERES-Godwin and CENTURY-Parton SOC simulation. The probabilistic estimates for quantified variables in the BBN model were well-aligned with the known and established responses for the assessed parameters. The results highlight the adaptation of conservational tillage can lead to sustainable crop yield, reduced runoff, and GHG-CO2 emissions by offering precipitation interception and available water space in soils (Rahman et al., 2021).
2.3 BBN model outputs
The output variables of the BBN model are wheat crop yield, surface runoff and total GHG-CO2 emissions. The response of output variables is the reflection of model input variables and their vital interactions with associated variables based on their identified conditional phenomena in a hierarchical order, as shown in Figure 1. Results found model output variables have exhibited seasonal variations in their responses through long-term simulations. Weather indices have an imminent influence over crop yield, surface runoff, and total GHG- CO2 emissions. Tillage preference for no-till or conventional methods has a comparable influence over SOC levels in the soils cultivated over extended periods. Tillage practices have profound effects on carbon sequestration and preserving soil structure. This investigation evaluated the model performance and validated that using test cases. The study also analyzed posterior probability distributions of interacting variables to assess their comparative impacts on each other.
2.4 The BBN model performance and evaluation procedure
Various metrics are applied to evaluate the performance of the BBN model to measure uncertainty (Table 5). These include model complexity (e.g., the number of nodes, links, etc.) (e.g., scenario evaluations) and sensitivity analysis (e.g., sensitivity to findings, case file simulation, etc.) using Netica software (Marcot, 2012; Namdari and Li, 2019).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-g002.jpg)
Figure 2. This comparison shows an influence analysis of the BBN model evaluation based on comparing two scenarios, e.g., (A,B), one with favorable and the other with unfavorable effects outcomes.
2.4.1 Model complexity analysis
The BBN model network is a directed acyclic graph, and complexity depends upon the number of variables and their interactions included in the BBN model; e.g., many nodes with numerous links/ arrows can make a BBN model more complex and inefficient. The BBN model consists of 11 nodes (variables), 17 links (conditional relationships), and 2,800 conditional probabilities learnt through parametric learning from simulated datasets by the DSSAT (Thogmartin, 2010). The number of discretized states of continuous variables also determines model complexity, as highlighted in Supplementary Table S8. The BBN model network contains the following junction trees with member nodes, as shown in Supplementary Table S9. A Clique is a set of nodes connected to another set of nodes, and clique size represents a complete subgraph of a given graph which contains the maximum number of nodes, e.g., vertices and edges in a graph (Almond and Almond, n.d.). Hence, this BBN model is neither very complex nor very simple shown in Figure 1 and Table 3.
2.4.2 Influence analysis
Influence analysis evaluates posterior probability distributions from chosen input variables and results in good-to-best or bad-to-worst-case scenario values. This information helps in decisions by measuring individual or set input variables’ impact on outcomes, as shown in Figure 2. Probabilistic inferences are illustrated on uncertainty estimates measured in 0–100% percentages. These probabilities represent the likelihood or confidence level associated with different states or outcomes of the variables in the model (Hájek, 2002). It’s important to note that even small changes in certain influential factors or inputs can lead to variability in probabilistic estimates. Model parametric changes can also reveal fluctuations in these values. Sensitivity analysis and other evaluation metrics were performed to highlight this and described in the relevant sections.
2.4.3 Sensitivity analysis
Using Netica, the BBN model performed a sensitivity analysis called sensitivity to findings for GHG-CO2 emissions, runoff, and yield. This analysis ranked and ordered the input variables in the Supplementary material, providing a results list in Supplementary Tables S10a–c. Tabulated variables are listed in descending order regarding the above variables (output variables) as per the generated report by Netica. Variance reduction, mutual information, and variance in belief values are reduced in descending order, depicting the comparable sensitivity of each variable with that of subject (output) variables (Supplementary Tables S10a–c). This method quantitatively compares the variables’ ability to reduce uncertainty or variance (entropy reduction) and allows for mutual comparisons (Thogmartin, 2010; Marcot, 2012; Namdari and Li, 2019).
2.5 Model validation procedure
The BBN model is validated using simulated datasets due to the scarcity of real data availability. Access to real-world data sources, resource constraints, and complications linked to data protection are common challenges in most modeling works. Establishing a long-term tillage-based experiment in competing time frames with challenges that arose due to COVID-19 pandemic implications has prompted this shift toward assessing data through simulation methods. This approach brings some limitations while using the simulated dataset for model validation. These may not capture the adaptation in real-world situations, face result validation concerns, have limited capacity in capturing factual variability, and carry bias embedded in algorithms defined in tools used for data simulations. However, simulated datasets offer a great opportunity to explore complex scientific systems where access to real-world data is a valid concern. This approach presents probabilistic solutions and diverse scenarios for comparing consequences, which could support decision-making for the right choices.
The assessment employed 25% of the synthetic dataset obtained from the DSSAT-generated simulations to test the BBN model for performing model validation in this research (Morgan and Granger and Max Henrion, 1990; Netica, 2024). The test cases representing datasets are used for response variables for 2012–22, the specimen listed in Supplementary Figure S7. This approach is also known as the holdout method for model validation (Jones et al., 2011). Confusion matrices GHG-CO2 Emissions, Runoff, and Crop yield are tabulated in 6, 7, and 8, respectively. Error rate illustrated two types of error by the BBN model, e.g., type 1 error for false positive, which means rejecting a true hypothesis, and type 2 error means failing to reject a false hypothesis. The sum of these is reflected in the confusion matrix.
The study utilized Netica software to evaluate the model’s performance, which applied scoring rules based on its system-generated computations. The scoring rules are the model evaluation metrics and assess the model’s performance of probabilistic estimates or predictive distributions. These rules measure the predicted distribution, the observed outcome, and the model’s accuracy and calibration. Netica software retrieved the report with a confusion matrix carrying the scoring rule outcomes of the individual output variable. Scoring rule outcomes highlighted the actual beliefs of the states that followed the values in the case file. Results found three types of scoring rules (Logarithmic loss, Quadratic loss, and Spherical loss) performed by Netica software (Pearl, 1978; Netica, 2023).
Logarithmic loss value was calculated based on a natural log and ranged from zero to infinity, where zero represents the best performance at an acceptable level in all the above cases. For instance, the Logarithmic loss value is depicted as 0.9555 from the confusion matrix (for GHG-CO2 emissions) under Table 6. This value is closer to zero than infinity and declares an acceptable performance level. Logarithmic loss penalizes incorrect predictions more seriously when they are confidently wrong. Hence, the Logarithmic loss value was calculated based on a natural log and ranges from zero to infinity, where zero represents the best performance and having the value of 0.9555 is an acceptable level of model performance in the given example.
Quadratic loss provides another metric to evaluate the performance of a model and is also known as the Brier score and ranges from zero to 2, where zero is the best. Quadratic loss penalizes larger errors more severely than smaller errors. Therefore, this value becomes sensitive to outliers, as larger errors contribute disproportionately to the overall loss. The quadratic loss value was depicted as 0.4859 from the confusion matrix (for GHG-CO2 emissions) under Table 6. This value is closer to zero and sets an acceptable level for the model performance. Spherical payoff is another matrix to measure the performance of the BBN model by considering the accuracy and calibration of predictions. This scoring value is linked to the balanced risk–reward profile of a sphere. The spherical payoff scoring rule aims to reward the model for achieving high accuracy through making correct predictions following a well-calibrated model for accurate estimates of the likelihood of those predictions. This matrix measures the model performance using an index ranging [0,1], where 1 is the highest value for illustrating better model performance. In the above example, the Spherical loss value was depicted as 0.7122, which is closer to one and considered an acceptable level for the model performance. Similarly, this BBN model performed better than average for all output variables in dealing with the nuances of probability values under scoring rule results as depicted through the Confusion matrix (for Runoff) and Confusion matrix (for Crop yield) and shown under Tables 7, 8, respectively (Marcot, 2012; Netica, 2023).
However, the model was also validated using the K-Fold cross-validation method and tested through the train split approach. The comparative results are summarized in Table 4. The details are provided in the Supplementary material S1. Moreover, the multiple regression model was employed in this analysis and findings are detailed as Supplementary material S2. Such comparisons evaluate the modeling performance among different models, such as multiple regression and network models and promote the use of the business value approach (Bansal et al., 1993).
3 Results and discussion
3.1 Posterior probability distribution analysis
The BBN model was developed to integrate prior knowledge using simulated datasets. It highlights analyzing model inferences through probabilistic estimates represented by the modeled parameters. Salient features are described below. The BBN model captured uncertainties linked to seasonal variations by providing posterior probability data distributions reflecting probabilistic relationships across all interacting variables in the model and offering decision choices for competing tillage preferences. This study analyzed the posterior probability distribution of intermediary and output variables in the BBN model. This evaluation technique revealed the performance of the best and worst-performing variables individually and on average (Solheim, 2021).
3.1.1 Wheat crop yield
Wheat crop yield varied in quantitative results among three real soils and their four modeled soils, each under No-till versus Conventional till cultivations using simulated data generated in the DSSAT v4.8 for over 48 years. Results found the maximum product weight (kg/ha) in the “NT.OC-7%PDPG8,” e.g., 8,575 kg/ha in the year 2016 and the minimum product weight (kg/ha) in “CT.OC-1%PDPG6,” e.g., 174 kg/ha in 1996 (Refer to Figure 3; Table 9). This year is ranked among the top three wettest seasons in the entire time series (Refer to Figure 4) and correlates with lower wheat production in soils with lower levels of SOC. Similarly, “CT.OC-7%PDPG8” averaged the highest product weight at 6,008 kg/ha, while “NT.OC-1%.PDPG6” averaged the lowest at 2,087 kg/ha. Therefore, the No-till method produced the highest wheat crop yield in soils with 7% SOC from PDPG8, while the conventional till method yielded the minimum with 1% SOC from PDPG6. Results are aligned with several research findings (Soane et al., 2012; Nunes et al., 2018). Seasonal variation in Wheat yield (kg/ha) is available in Figure 3 with basic statistics highlighted in Table 9.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-g003.jpg)
Figure 3. Seasonal variation in Wheat yield (kg/ha) in 30 soils under no-till versus conventional till using simulated data from DSSAT v4.8 since 1975–2022.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-t009.jpg)
Table 9. Statistical analysis of wheat yield (kg/ha) distribution in 30 soils under No-till versus Conventional tillage - Refer to Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-g004.jpg)
Figure 4. Comparison of cumulative surface runoff (mm) in 11 wet seasons among 30 soils under no-till versus Conventional till using simulated data from DSSAT v4.8 since 1975–2022.
3.1.2 Comparison between observed data and simulated data of wheat yield (kg/ha) among 30 soils under no-till versus conventional till from 2016 to 2021
Wheat yield was compared between observed and simulated data among 30 soils under No-till versus Conventional till from 2016 to 2021. Observed data of average wheat yield from the same catchment was tracked from published data (Boote et al., 2008). The comparison showed that soils with higher carbon content levels showed higher yields overall. PDPG8 soils demonstrated strong performance, closely aligning with the observed average yield. However, the observed yield showed the highest produce level in 2016, 2017 and 2019, followed by PDPG soils with higher soil carbon content levels from 3 to 7% among NT and CT practices in all simulations. The PDPG8 soil with organic carbon contents of 7% outperformed in 2018, 2020 and 2021 by yielding 6,922 kg/ha, 6,351 kg/ha, and 8,240 kg/ha, respectively. However, the soils “CT.OC-7%PDPG8 & NT.OC-7%PDPG8” consistently produced sustainable yields over the competing period. The observed and simulated data for the competing soils are depicted by comparative wheat yields in Figure 5. Another study presented a probabilistic model for improved forecasting of crop yield under environmental uncertainty. (Supit, 1997; Mahmood et al., 2023)
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-g005.jpg)
Figure 5. Comparison between observed data and simulated data of wheat yield (Kg/ha) among 30 soils under No-till versus Conventional till since 2016–2021.
3.1.3 Surface runoff
Cumulative runoff (mm) varied in quantitative results among all real and modeled soils under No-till and Conventional tillage, simulated over 48 years with DSSAT v4.8. The highest cumulative runoff of 44 mm was recorded in “CT.OC-5%PDPG7” in 2014, a year characterized by flooding. In the same year, “NT.OC-5%.PDPG7” and “NT.OC-7%.PDPG7” both experienced a maximum runoff of 26 mm. Among 30 soils, 18 had a seasonal minimum runoff of 0 mm. On average, “CT.OC-5%PDPG8” had the highest cumulative runoff (18.15 mm), while “NT.OC-1%.PDPG6” had the lowest (2.38 mm). Seasonal variation is evident in Figure 6, with basic statistics in Table 10. These findings are comparable with other studies (DeLaune and Sij, 2012). This indicates that NT results in decreased runoff in soils with both higher and lower organic carbon levels. This response assumes that higher SOC levels could no longer help water retention after reaching saturation in a high rainfall season. This aspect highlights the level of resilience or vulnerability to challenging surface runoff characteristics among soils with high organic carbon (OC) levels under conventional tillage (CT).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-g006.jpg)
Figure 6. Seasonal variation in surface runoff (mm) in 30 soils under no-till versus conventional till using simulated data from DSSAT v4.8 since 1975–2022.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-t010.jpg)
Table 10. Statistical analysis of surface runoff (mm) distribution in 30 soils under No-till versus conventional tillage - See Figure 6.
3.1.4 Comparison of cumulative surface runoff (mm) in 11 Wet seasons among three soilscapes under no-till versus conventional till methods using simulated data generated in the DSSAT v4.8 for over 48 years
The results show that in 2014, during a flooding year, “CT.OC-5%PDPG7” had the highest cumulative runoff of 44 mm over 11 wet seasons. Conversely, “NT.OC-5%PDPG6” had the lowest cumulative runoff at 3 mm in 2001 and 2010. On average, “CT.OC-5%PDPG7” had the highest cumulative runoff (28.82 mm), while “NT.OC-1%.PDPG6” had the lowest (5.55 mm) across all wet seasons. These findings align with previous research (Gautam et al., 2015). The comparison is depictable in Figure 4. There are strong implications of higher levels of cumulative surface runoff that can lead to soil erosion and deterioration of water quality. Runoff over the soil surface removes soil particles and sediments away from their location and becomes a cause of soil loss, which results in soil fertility loss and impacts soil health and productivity. Runoff carries substances and residuals of agrochemical pollutants to the hydrological system, which endangers the ecosystem’s balance. Hence, suitable tillage approaches such as reduced to no-tillage can help conserve soil aggregates’ moisture level and facilitate land cover by reducing evapotranspiration through retaining crop residuals.
3.1.5 Carbon – soil organic matter
Soil organic matter represents all organic components of soil, while organic carbon (OC) refers to the carbon components within the SOM. Carbon in soil organic matter (C-SOM) varied in quantitative results across real and modeled soils under both No-till and Conventional tillage, simulated over 48 years with DSSAT v4.8. Seasonal changes in C-SOM were observed across all soils. Positive changes in C-SOM were found in soil “CT.OC-7%PDPG8” (e.g., +862 kg/ha in 2012), while negative changes were seen in “NT.OC-1%.PDPG7” (e.g., −200 kg/ha in 1996, a wet season). On average, “NT.OC-7%PDPG8” had the highest C-SOM (e.g., 637 kg/ha), while “NT.OC-1%.PDPG7” had the lowest (e.g., −151 kg/ha). These findings are comparable with other research findings (Pinheiro et al., 2015). Seasonal change in C-SOM is confirmed in Supplementary Figure S2. This response assumes that NT practices support higher organic carbon levels in soils. Carbon–soil organic carbon (C-SOM) is vital in maintaining soil health and carbon sequestration. C-SOM is not limited to its functional role in nutrient cycling, maintaining soil structure, enhancing water retention, facilitating microbial activity and soil biodiversity support. Tillage modifications can contribute to carbon sequestration by reducing soil carbon loss by minimizing the direct contact of soil particles and tilling equipment. Hence, with higher C-SOM levels in soils, soil resilience and resistance to soil erosion will result in better soil health conditions. Suitable C-SOM levels also support a balance between soil flora and fauna for increased soil biodiversity and facilities for ecosystem services.
3.1.6 Comparison between minimum and maximum C-SOM (kg/ha) in 30 soils under no-till versus conventional till methods using simulated data generated in the DSSAT v4.8 for over 48 years
The maximum soil carbon content (C-SOM) of 286,634 kg/ha was found in “NT.OC-7%PDPG8” in 1984, while the minimum of 5,894 kg/ha was observed in “CT.OC-3.9%PDPG8” in 2009. These results indicate higher C-SOM levels in NT compared to CT cultivation (Ferreira et al., 2020). Analysis from 1975 to 2022 using DSSAT v4.8 simulated data showed that NT had a higher average minimum C-SOM (148,752 kg/ha) than CT (130,315 kg/ha). Conversely, NT contained 148,920 kg/ha more on average for maximum C-SOM, whereas CT contained 131,995 kg/ha. This information is evident in Supplementary Figure S3. The observed differences in C-SOM levels between different soil types and tillage methods have ecological implications for carbon sequestration, erosion control, soil health, soil biodiversity, nutrient cycling, and water quality. Hence, tillage methods which can facilitate these aspects objectively should be preferred and promoted. For example, reduced tillage, non-inversion, zero or conservational tillage could endorse achieving the above phenomena for advancing ecological stability.
3.1.7 Comparison between minimum and maximum GHG-CO2 emissions (eqv. Kg/ha) among 30 soils under no-till versus conventional till methods using simulated data generated in the DSSAT v4.8 for over 48 years
Results found that “CT.OC-7%PDPG6” soil emitted the highest GHG-CO2 at 34,806 kgCO2eqv./ha in 1996, while “NT.OC-7%PDPG6” emitted 10,483 kgCO2eqv./ha in the same year. This information suggests higher emissions in soils with higher organic carbon content under CT cultivation and significantly lower under NT. Analysis from 1975 to 2022 using DSSAT v4.8 simulated data showed that NT had lower average minimum GHG-CO2 emissions (3,985 kgCO2eqv./ha) than CT (7,415 kgCO2eqv./ha). Conversely, NT emitted 8,747 kgCO2eqv./ha less average for maximum emissions, whereas CT emitted 15,356 kgCO2eqv./ha. These findings revealed that NT cultivation is better for having minimum and maximum GHG-CO2 emissions (total average scale) at reduced levels (Six et al., 2004; Mangalassery et al., 2014). The comparison is obvious in Figure 7. Soil organic matter contents, crop residue management, crop fertilization, soil moisture, soil temperature, and tillage intensity are the potential drivers of GHG-CO2 emission between different soils and tillage methods. However, suitable strategies can be adopted to mitigate them. These are not limited to adapting tillage modification (reduced to no-till), cover cropping (cropping between cash crops), crop rotation (rotating crops), agroforestry (tree planting in agri-landscapes), manuring and organic amendments (compost manuring), and precision agriculture (precision agriculture technologies).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1402098/frai-07-1402098-HTML/image_m/frai-07-1402098-g007.jpg)
Figure 7. Comparison between minimum and maximum GHG-CO2 emissions (KgCO2eq./ha) among 30 soils under no-till versus Conventional till cultivation using simulated data from DSSAT v4.8 since 1975–2022.
3.1.8 Available water space
Significant variation in available water space (cm3 cm3) was observed among three real soils and four modeled soils under No-till and Conventional till methods using simulated data from DSSAT v4.8 over 48 years. The maximum available water space was recorded in “CT.OC-7%PDPG8” at 528 (cm3 cm3) in 1996, followed by “NT.OC-7%PDPG8” at 520 (cm3 cm3) in the same year. In contrast, the minimum was observed in “NT.OC-1%.PDPG7” at 127 (cm3 cm3) in 2022. On average, “CT.OC-7%PDPG8” exhibited the highest available water space at 508 (cm3 cm3), while “NT.OC-1%.PDPG7” had the lowest average available water space at 138 (cm3 cm3). These findings show that soils with higher OC levels under CT and NT in both cultivations provide more available water space. However, the choice and management of tillage practices can particularly influence soil water dynamics, impacting water space availability in soils by altering soil structure, compaction, organic matter content, and surface cover to mitigate evaporation. That’s why Available Water Space (AWS) has the highest number of relationships (arc links) with other variables (nodes) in the BBN (see Figure 1). Additionally, AWS is sensitive to soil features and carbon-soil organic matter, as indicated in Supplementary Table S10d in the Supplementary material. Results are comparable with a similar study in the US (Hill, 1990). Seasonal variation in available water space (cm3 cm3) is illustrated in Supplementary Figure S4.
3.1.9 Precipitation interception
Among the three real soils and their four modeled soils, precipitation interception (%) revealed significant variation under both No-till and Conventional till methods using simulated data generated in DSSAT v4.8 over 48 years. The maximum “precipitation interception (%)” was noted in soils labeled “NT.OC-7%PDPG8” and “CT.OC-7%PDPG8,” up to 0.59 percent in 2016. In contrast, the minimum “precipitation interception (%)” was noticed in soils labeled “NT.OC-1%.PDPG6” and “CT.OC-1%PDPG6,” with only 0.11 percent recorded in 2010. However, considering average values, soils labeled “NT.OC-7%PDPG8 & CT.OC-7%PDPG8” exhibited the highest precipitation interception (%) at 0.44 percent. Conversely, soils labeled “NT.OC-1%.PDPG6 & CT.OC-1%PDPG6” pointed to the lowest average precipitation interception (%) at 0.27. The findings are comparable with other studies (Leuning et al., 1994). Supplementary Figure S5 shows seasonal variation in precipitation interception.
3.2 Strengths and limitations of the BBN model of tilling wheat conflicting climatic and sustainability challenges
This BBN model assessed tilling systems under climate change and sustainability challenges for winter wheat crop production. The assessment compared 48 cropping seasons of simulations for each included variable response and assessed their long-term impacts on wheat yield, surface runoff, and GHG-CO2 emissions as output variables (Nielsen and Vigil, 2010; Krauss et al., 2017; Nath et al., 2017). This model gained insight into those parameter responses where empirical data is not readily available except for soil, weather, wheat management and yield data. However, the synthetic datasets were generated using established machine-learning models in DSSAT. Then, the BBN model utilized a crop growth hierarchal model DSSAT output data from its prototype model experiment. The BBN model offers integration of identified variables using standardized responses to capture uncertainties among model parameters. These responses reflect conditional dependencies and are consistent with the established knowledge in the domain.
It is important to understand BBN, which provides a graphical probabilistic model representing a set of variables and their probabilistic dependencies using a directed acyclic graph. It modeled uncertain relationships between variables and made probabilistic inferences based on observed evidence. Contrarily, an emulator is a simplified mathematical model that approximates the behavior of a more complex, computationally expensive model. Emulators are typically used to speed up the simulation process by quickly estimating model outputs based on a subset of input values and mimicking the source (Aakula et al., 2024). However, the DSSAT is a mechanistic crop growth hierarchical model and serves as a simulator in this study. This tool operated based on its embedded machine learning capabilities and provided the DSSAT output data used as an input in this BBN. This BBN model presents probabilistic inferences and reasoning, serves different purposes, operates on different assumptions. However, this approach may have the potential to emulate functions not evaluated in this study (Liu and West, 2009). BBN models could have other challenges in building and maintaining large complex networks, and computational expansive, especially data-driven networks are highly dependent on the amount and quality of data. The BBN models can not represent feedback loops or cyclic dependencies between variables which could challenge applicability to certain domains. In such cases, alternative probabilistic models such as hidden Markov models or Markov random fields might be a better substitute. Nonetheless, this BBN model did not come across the challenges above (Wang, 2004).
The BBN model is appraised per established metrics and summarized in Table 5. No contradictory responses were observed from model performance except for a higher error rate in predictive accuracy highlighted in the confusion matrix. The error rate is higher for GHG, runoff, and yield and is reported as 29.39, 33.64, and 43.94%, respectively. This outcome could be due to a greater variation in the organic carbon contents defined with varying fluctuations among 30 soils (Tibshirani and Tibshirani, 2009). Three real soils with built in soil organic level contents have been used to simulate each into four varying OC levels. Hence, the variation in GHG, runoff, and yield must appear in the results and are also subject to a higher error rate. However, the model’s predictive accuracy is still over 56% among all output variables of interest using the holdout method of model validation. The model was trained using data from 1975 to 2011 and tested against 2012–2022. However, the model was also validated using the K-Fold Cross Validation (CV) method using K-fold splits into training and test datasets. Then, the model was trained using data from 1975 to 2014 (K1 + K2 + K3 + K4) and tested using data from 2015 to 2022 (K5) to yield consistent and reliable results, proved (Table 4; Supplementary material S1). This challenge commonly appears among models in several agroecological studies, and an example is presented in Figure 5 by comparing published crop yield datasets (2016–2021) from research covering the same catchment area (Hazarika et al., 2009; Karimi et al., 2021; Mahmood et al., 2023). Model performance could be improved by employing parameter tuning, varying (scaling) parametric categories, changing (normalization) datasets, data homogeneousness, etc (Singh and Singh, 2022). Further exploration of the DSSAT tool for its interacting functionalities could be another area of interest for comparing model parameters’ performance against those of other competing tools.
The BBN model illustrated that the No-till (NT) choice of cultivation outclassed the Conventional till (CT) by reducing GHG-CO2 emissions probability estimation in the highest band from “7.34 to 7.31%.” Probabilistic estimates are expressed in percentages (e.g., 0 to 100%), representing the likelihood or probability of different outcomes based on the available evidence and model assumptions. They can be presented as ranges or uncertainty intervals, with one bound (lower or upper) specified (Hohle and Teigen, 2018). NT also exhibited a probability of a rising trend in higher wheat yield bands, increasing from “7.46 to 7.56%.” Cumulative runoff probability decreased in higher bands, with the “high band” diminishing from “8.52 to 8.50%.”
In contrast, the CT option had a snowballing trend in GHG-CO2 emissions in the highest band from “7.34 to 7.36%.” CT resulted in probabilistic estimates of lessening in wheat yield probability in the highest band from “7.46 to 7.35%,” while cumulative runoff amplified in higher bands, with a probability of the “high band” rising from “8.52 to 8.55%” (Radford, 2007; Lu et al., 2015). The impact of tillage methods on GHG-CO2 emission is linked with the potential phenomenon of creating soil disturbance in adapting tilling frequency and intensity, which can trigger disturbing soil aggregates, impacting water retention and moisture availability for plant uptake and influencing root development. These methods also influence soil microbial communities and their activities in decomposing soil organic matter. Nutrient cycling is another aspect affected by tilling choices, which can favor or resist underlying mechanisms. Tilling preference influences SOC contents. Higher levels of OC contents favor maintaining soil health conditions for building better soil structure and more resilient soils to help achieve better crop production.
The BBN model effectively apprehended uncertainties linked to seasonal variations by providing posterior probability data distributions reflecting probabilistic relationships across all variables and offered decision choice for NT favoring soil carbon stocks (highest among soil “NT.OC-7%PDPG8,” e.g., 286,634 kg/ha) in winter wheat. Furthermore, the model favors reduced GHG-CO2 emissions through NT practices and naturally restraining surface runoff generations for flood mitigation (Maraseni and Cockfield, 2011). For additional improvement of the BBN model, examination of early and late seasonal planting, supply rates of nutrients, and the use of various crop varieties under diverse farming systems are to be explored. A few constraints are not limited to access to empirical and synthetic datasets used in this model. Simulated datasets are used to calibrate and validate the BBN model (Aliferis and Cooper, 1994; Sebastiani and Ramoni, 2001; Newlands et al., 2014; Ali, 2023), so results variation could be expected where entire real/ empirical datasets will be applied.
This study focuses only on limited output variables, e.g., GHG-CO2 emissions and does not cater to nitrous (N2O) and methane (CH4) emissions. Surface runoff measurement presents only standard soil-specific features in the S-module of the DSSAT and could consider slope and topographic features for future exploration. The study utilized simulated wheat yield and compared that with real field data from the catchment area for 2016–2021. Utilizing catchment-wise yield data could enhance yield estimations on a larger scale. While simulated datasets were generated solely using the DSSAT tool, exploring other competing tools could yield comparable results. These areas offer futuristic research on this model. Moreover, variation in model response could also be expected due to spatio-temporal factors and the extent of noisy datasets once employed. Hence, perfect datasets can deliver more reliable model outcomes (Merino et al., 2016).
This research employed simulated datasets, offering valuable research, experimentation, and decision-making prospects. However, using simulated datasets entails implications regarding predictive accuracy, cost-effectiveness, experimental flexibility, risk mitigation, and ethical considerations. Researchers must carefully consider these factors when utilizing simulated data in their studies. Seeking real-world datasets whenever possible allows for comparing predicted results with observed data, mitigating some of these concerns. This study compared wheat yield datasets for comparing against simulated datasets and found few soils exhibited comparable results.
Overall, valuable improvements were made to the BBN, which would not have been possible without data access. This was particularly relevant in this study, where dataset access for a longer period was nonexistent. This approach helps integrate DSSAT simulated data into a BBN, provides an opportunity to measure against actual wheat yield (Figure 5) and provides a good insight to compare the model and identify the higher-performing soils under competing tilling methods. Hence, the model can potentially improve the accuracy, realism, and robustness of the model’s predictions and enable more informed decision-making in agricultural systems.
The current iteration of the BBN is in its alpha stage and recommends testing in diverse real-world scenarios to progress into beta versions. Successful scalability may lead to the development of a commercial gamma-level BBN model, customizable for diverse applications. Furthermore, enrichments in decision support functionalities, informed by real-life data testing, could enable the introduction of individual predictive models as independent software tools for mobile apps and portal-based interfaces.
The following are evident advantages achieved by this BBN model.
• The BBN model uniquely compared two different tilling systems, such as No-till (NT) and Conventional till (CT), magnifying the superiority of NT practices over CT.
• This BBN model effectively captured the uncertainties of key parameters responsible for winter wheat tilling systems toward climatic and sustainability challenges (Refer to Supplementary Tables S7, S8).
• The BBN model highlighted the dynamic nature of model variables and accurately reflected the responses of output variables based on the conditionality between interacting variables.
• This BBN model highlighted the role of C-SOM levels influenced by the competing tilling systems. This helped analyze the carbon sequestration impact between competitive tilling methods, e.g., NT methodology favors C-SOM over CT.
• The BBN model provided probabilistic precipitation retention/resilience estimates, comparing tilling techniques. This analysis can aid in assessing flood alleviation impacts, with the No-till approach showing a higher capacity for precipitation interception than Conventional tillage.
• The BBN model showed that adopting the No-till method reduced GHG-CO2 emissions compared to Conventional tillage practices.
The model can be used to inform decision-making and improve agricultural sustainability. For instance, adapting suitable tillage preferences is a crucial management choice with a commercial focus on obtaining crop yield. Conventional tillage choice could deliver instant yield gains at the cost of hampering long-term soil health and structural resilience. This approach will cater only to the interest of acquiring optimum crop yields at the cost of attracting several complicated issues of soil compaction erosion, consequently attracting higher surface runoff & greater flooding risk in the long run. For example, CT promotes frequent and intense tillage application, exerting mechanical pressure through interfacing soils. Heavy and frequent tillage can expose the soil organic matter to increased oxygen levels and microbial decomposition, affecting soil functioning for crop growth and triggering GHG-CO2 emissions. Tillage methods can also influence soil nutrient cycling processes by disturbing soil physical, chemical, and biological properties. Conventional tillage practices could accelerate the mineralization of organic matter, releasing nutrients and raising nutrient availability for plant uptake with an increased risk of loss of nutrients. These changes in nutrient availability affect crop growth, yield and nutrient use efficiency, which could influence soil carbon dynamics and GHG-CO2 emissions.
However, the reduced-to-no-till choice of cultivation can bring multifaceted benefits along with attaining sustainable crop yields over time, limiting the loss of SOC, preserving soil health conditions, limiting surface runoff due to land cover vegetation preference favoring flood alleviation, and reducing the risk of soil compaction and erosion (Farahani et al., 2022). At the helm of climate change and sustainability challenges, the importance of carbon sequestration has become manifold. This model supports reduced to no-till (NT) as the right choice of tillage practices through which crop residues and re can be retained on the land to get incorporated into soils, supporting land cover and increasing SOC contributions. This aligns with sustainable crop farming practices.
4 Conclusion
This paper presents a novel contribution by introducing a BBN model for no-tilling as an NFM strategy for wheat crop cultivation. This model presents valuable insight into the effects of tillage preferences on soil carbon dynamics, greenhouse gas emissions, and crop yield. The study promotes understanding how tilling choices influence the relationships between soil health, crop productivity, and ecological sustainability. This research identifies best tilling practices for management of SOC for sustainable crop yield while minimizing environmental impacts. Reduced to no-till practices promote offset of greenhouse gas emissions and encourage agroecosystem resilience. These research findings have implications for extending policy efforts for promoting sustainable soil management practices and provide scientific evidence of reduced to no-till practices’ benefits. These conclusions inform policy decisions, support extension efforts for the farming community, and incentivize the adoption of sustainable tilling practices.
This BBN model offers decision support among choosing alternative tilling systems to achieve sustainable wheat production by reducing surface runoff, and GHG-CO2 emissions by delivering probabilistic estimates of pertinent parameters. This BBN offers unique relationships between variables of multiple interests, such as seasonal wheat crop production, alleviation of floods and GHG-CO2 emissions. This BBN model provides an opportunity to cater to precipitation interception and soil water absorption to alleviate flooding risk in winter wheat cultivation regions while aiming for sustainable wheat production and soil carbon stocks. This BBN endorses the sustainable development goals defined by the FAO of the United Nations through promoting sustainable production and climate action while practicing sustainable farming.
This comparison illustrates the impacts of No-till (NT) and conventional till (CT) cultivation methods on three real soils from the Pang catchment area, each modeled with varying soil organic carbon (SOC) content in topsoil at four levels (1, 3, 5, 7%). Synthetic datasets were produced using the DSSAT tool over 48 years to evaluate long-term impacts on wheat crop yield, surface runoff, and GHG-CO2 emissions. DSSAT output files containing variable responses were accessed, reflecting embedded crop growth models. DSSAT, with multiple interfaces, served as a data simulator for this study and provided variable responses for soil attributes, weather, and crop management practices in the Pang catchment area. Contrarily, Netica software established a single interface for a BBN model, comparing tillage preference to wheat yield, surface runoff, and GHG-CO2 emissions. The pertinent variables’ data were extracted from DSSAT simulation compilations and formatted into a specified case file for BBN parameterization. The rationale for combining DSSAT and BBN strengthens from their complementary strengths in agricultural modeling and uncertainty quantification. DSSAT presented a well-established agricultural simulation model that represented crop growth, soil processes, and management practices for comparable model outputs. By contrast, the BBNs are adept at handling uncertainty and integrating various sources of information to make probabilistic predictions.
A BBN model was developed to capture the uncertainty in data distributions and model the variables’ interactions as a single interface tool. This approach helped analyze the impacts on model outputs, e.g., crop yield, runoff, and emissions. Tillage methods indirectly affect model-output variables of crop yield, surface runoff and GHG-CO2 emissions through carbon sequestration. These effects are insightful by various mechanisms such as soil structural changes caused by conventional tilling methods where soil aggregates are broken down, increase soil compaction, and reduce soil porosity, and water infiltration rates. These changes affect water retention, root penetration, and nutrient availability, which influence soil carbon dynamics, microbial activity and crop growth. Microbial activity is another crucial mechanism in organic matter decomposition, nutrient cycling, and soil health. Tilling methods affect soil physical and chemical properties by disturbing microbial activities and community composition, such as moisture content level, pH, oxygen availability, microbial habitat disturbance, and substrate availability.
The study revealed that NT and CT do not directly impact runoff and product weight gains but rather influence through intermediate variables, as shown in Figure 1. However, these methods can indirectly influence output variables through the implications of carbon sequestration organic matter (C-SOM) by altering the rate of C-SOM reduction in the soil profile over time. NT soils exhibited lower rates of C-SOM reduction than CT soils, indicating that the NT method enhances resilience toward C-SOM contents through multiple pathways. Higher C-SOM levels in soils provided greater available water space, reducing runoff impact and contributing to higher biomass and product weight. The individual model parameter is categorized in ascending order, ranging from lower to higher bands. This study discovered from variable responses from their simulated dataset response over 48 years that probabilistic estimates for the NT option are better than CT because this reduced GHG-CO2 emissions in the highest band from “7.34 to 7.31%.” It exhibited a rising trend in higher bands of wheat yield, with the highest band rising from “7.46 to 7.56%.” The cumulative runoff also tends to reduce in higher bands, with the high band declining from “8.52 to 8.50%.” Contrarily, the CT option increased GHG-CO2 emissions in the highest band from “7.34 to 7.36%.” This option also reduced the wheat yield in the highest band from “7.46 to 7.35%,” with cumulative runoff trending to rise in higher bands with the “high band” declining from “8.52 to 8.55%”.
The research concludes that soils with higher OC contents have relatively higher GHG-CO2 emissions, but this trend decreases in soils under NT applications. Hence, promoting NT cultivation for soils with higher OC contents can help reduce GHG-CO2 emissions, while the same applies to soils with lower OC contents to sustain emissions at reduced levels. Compared with CT, NT cultivation supports precipitation interception through biomass land cover, reducing the direct impact of rainfall on bare soils and minimizing flood risks. Additionally, higher crop biomass supports sustainable wheat crop yield.
The BBN model effectively captured uncertainties, offering posterior probability distributions reflecting conditional relationships across variables and offered decision choice for NT favoring soil carbon stocks in winter wheat (highest among soils “NT.OC-7%PDPG8,” e.g., 286,634 kg/ha) over CT (lowest in “CT.OC-3.9%PDPG8,” e.g., 5,894 kg/ha). On average, NT reduced minimum GHG- CO2 emissions to “3,985 kgCO2eqv/ha,” while CT emitted “7,415 kgCO2eqv/ha.” Conversely, NT emitted “8,747 kg CO2eqv/ha” for maximum emissions, while CT emitted “15,356 kg CO2eqv/ha.” This model represents probabilistic inferences based on specified datasets from the Pang catchment area. The model is recommended for testing by farmers, growers, landowners, or other stakeholders if they have datasets for all the variables used in this BBN. They may need to tune up the parametric adjustment for respective variables (nodes) in the model to incorporate the full range of data specified into appropriate categories or bin sizes to avoid loss of information. Moreover, the model can be customized based on farm-specific data and variables of interest.
However, improvements to the BBN model could involve considering early or late seasonal planting, diverse plant seedlings and nutrient supply rates, different farming systems, and various cultivars. There is also potential for future exploration to address multiple cropping, land fallowing, and crop rotation aspects. The BBN model could face potential challenges linked to long-term real-world datasets and computational resource availability. The research points out that the implications depicted from the model findings are confined solely to this study and were derived through specified data-driven exploration. Factors such as spatiotemporal variability, the availability of high-quality datasets, computational resources, and the diversity of soil conditions and farming systems may limit the accomplishment of comparable results across different situations. Moreover, this model could be investigated further for the next step to develop into a user-experience software application package and the model could also be tested for diverse spatiotemporal conditions for comparable variable responses.
Conclusively, reduced to no-tillage proved better than conventional tillage because the former encourages retention, maintenance, and storing of carbon stocks in soil and lessens GHG emissions. The model is a decision support tool that can be customized. For further improvement of the BBN model, exploration of early and late seasonal planting, nutrient supply rates and use of different cultivars under diverse farming systems are recommended. The research points out that the implications depicted from the model findings are confined solely to this study and were derived through specified data-driven exploration. Factors such as spatiotemporal variability, the availability of high-quality datasets, computational resources, and the diversity of soil conditions and farming systems may limit the accomplishment of comparable results across different situations.
However, research findings align with relevant sustainable development goals (SDGs), e.g., SDG12 and SDG13 for responsible production and climate actions defined by the Agriculture and Food Organization of the United Nations. SDG12 defines responsible consumption & production, and this research study promotes reduced or no-tillage practices as responsible production techniques that promote minimized soil disturbance, soil health, and limiting farm mechanization involving fuel and resource depletions. These research findings also encourage responsible use of resources, maintaining soil carbon stocks through tillage promoting carbon sequestration. SDG13 defines climate action, and this research study emphasizes adopting reduced or no-till practices as recommendations for contributing to climate action.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
QA: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author declares that financial support was received for the research, authorship, and/or publication of this article. The author appreciates the funding partners of this project, namely, The Douglas Bomford Trust and The School of Agriculture, Policy, and Development, University of Reading, United Kingdom.
Acknowledgments
The author expresses gratitude to all who offered support throughout the improvement of this manuscript, including casting reviews, feedback, and comments, from fellow researchers, scientists, editorial staff and the quoted work. Special thanks to esteemed mentors Lindsay Todman (for securing funding) and Martin Lukac (for overseeing the project). The author also thanks Dr. Maleki Badjana and the NFM Landwise Project partners for data synthesis and acquisition of NATMAP soilscape for the Pang catchment area. The author acknowledges the generous support of the project funders, including the “Douglas Bomford Trust” and the “School of Agriculture, Policy, and Development, University of Reading, United Kingdom.” Furthermore, the considerable response and comments toward improving this manuscript are also greatly acknowledged at all levels. Without such valuable resources, this project would not have met quality completion.
Conflict of interest
The author declares that the research was conducted without any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frai.2024.1402098/full#supplementary-material
References
Aakula, V., Vira, J., and Petre, I. (2024). Emulation of dynamic process-based agroecosystem models using Long Short-term memory networks.
About the Pang. (2019). Action for the River Kennet. Available at: https://storymaps.arcgis.com/stories/f2019a2649274e20b8663a59c411e5bf.
Adger, W. N., Arnell, N. W., and Tompkins, E. L. (2005). Successful adaptation to climate change across scales. Glob. Environ. Chang. 15, 77–86. doi: 10.1016/j.gloenvcha.2004.12.005
Ahmed, M., Ali, S., Zahid, A., Ahmad, S., Yasin, N. A., and Hayat, R. (2023). “Climate change and process-based soil modeling” in Global agricultural production: Resilience to climate change. ed. M. Ahmed (Cham: Springer International Publishing), 73–106.
Ahmed, M., Raza, M. A., and Hussain, T. (2020). Dynamic modeling. Syst. Model. 2020, 111–150. doi: 10.1007/978-981-15-4728-7_4
Ali, Q. (2023). A trade-off between farm production and flood alleviation using land use tillage preferences as a natural flood management (NFM) strategy. Smart Agric. Technol. 6:100361. doi: 10.1016/j.atech.2023.100361
Aliferis, C. F., and Cooper, G. F. An evaluation of an algorithm for inductive learning of Bayesian belief networks using simulated data sets. Uncertain. Proc., (1994), 1: 8–14.
Asseng, S., Ewert, F., Rosenzweig, C., Jones, J. W., Hatfield, J. L., Ruane, A. C., et al. (2013). Uncertainty in simulating wheat yields under climate change. Nat. Clim. Chang. 3, 827–832. doi: 10.1038/nclimate1916
Bansal, A., Kauffman, R. J., and Weitz, R. R. (1993). 1997. Construction and assessment of classification rules. John Wiley & Sons, New York, 232 pp. Bansal a, Kauffman RJ, Weitz RR. Comparing the modelling performance of regression and neural networks as data quality varies: a business value approach. J. Manag. Inf. Syst. 10, 11–32. doi: 10.1080/07421222.1993.11517988
Beerling, D. J., Leake, J. R., Long, S. P., Scholes, J. D., Ton, J., Nelson, P. N., et al. (2018). Farming with crops and rocks to address global climate, food and soil security. Nat. Plants 4, 138–147. doi: 10.1038/s41477-018-0108-y
Boote, K. J., Sau, F., Hoogenboom, G., and Jones, J. W. (2008). “Experience with water balance, evapotranspiration, and predictions of water stress effects in the CROPGRO model” in Response of crops to limited water: Understanding and modelling water stress effects on plant growth processes. ed. L. Ahuja, vol. 1, 59–103.
Broadmeadow, S., Nisbet, T., Palmer, R., Webb, L., Short, C., Chivers, C. A., et al. (2023). Incorporating technical and farmer knowledge to improve land use and management for natural flood management in lowland catchments. Land Use Policy 128:106596. doi: 10.1016/j.landusepol.2023.106596
Carter, M. R. (2002). Soil quality for sustainable land management: organic matter and aggregation interactions that maintain soil functions. Agron. J. 94, 38–47. doi: 10.2134/agronj2002.3800
Cassel, D. K., and Nielsen, D. R. (1986). “Field capacity and available water capacity” in Methods of soil analysis: Part 1 Physical and mineralogical methods. ed. D. K. Cassel, vol. 5, 901–926.
Cheema, S. S., and Singh, A. (2021). Analysis of soil management data module of DSSAT model. Int. J. Comput. Appl. Inf. Technol. 13, 390–399.
Collins, A. J., Vegesana, K. B., Seiler, M. J., O’Shea, P., Hettiarachchi, P., and McKenzie, F. (2013). Simulation and mathematical programming decision-making support for smallholder farming. Environ. Syst. Decis. 33, 427–439. doi: 10.1007/s10669-013-9460-7
Cramer, W., Guiot, J., Fader, M., Garrabou, J., Gattuso, J. P., Iglesias, A., et al. (2018). Climate change and interconnected risks to sustainable development in the Mediterranean. Nat. Clim. Chang. 8, 972–980. doi: 10.1038/s41558-018-0299-2
Crop Cover Data Source. (2024). Available at: https://environment.data.gov.uk/dataset/f0f54bc1-b77a-42c8-b601-2f4aaf4dd851.
Dal Ferro, N., Quinn, C., and Morari, F. (2018). A Bayesian belief network framework to predict SOC dynamics of alternative management scenarios. Soil Tillage Res. 179, 114–124. doi: 10.1016/j.still.2018.01.002
De Oliveira, D. B., and Gaudio, A. C. (2000). BuildQSAR: a new computer program for QSAR analysis. Quant. Struct. Act. Relat. 19, 599–601. doi: 10.1002/1521-3838(200012)19:6<599::AID-QSAR599>3.0.CO;2-B
Decision Support System for Agrotechnology Transfer. (n.d.). Decision Support System for Agrotechnology Transfer (DSSAT) tool. Available at: https://dssat.net/.
DeLaune, P. B., and Sij, J. W. (2012). Impact of tillage on runoff in long term no-till wheat systems. Soil Tillage Res. 124, 32–35. doi: 10.1016/j.still.2012.04.009
Donald, M. (2011). Using Bayesian methods for the estimation of uncertainty in complex statistical models (doctoral dissertation, Queensland University of Technology).
Droste, N., May, W., Clough, Y., Börjesson, G., Brady, M., and Hedlund, K. (2020). Soil carbon insures arable crop production against increasing adverse weather due to climate change. Environ. Res. Lett. 15:124034. doi: 10.1088/1748-9326/abc5e3
Drury, B., Valverde-Rebaza, J., Moura, M. F., and de Andrade, L. A. (2017). A survey of the applications of Bayesian networks in agriculture. Eng. Appl. Artif. Intell. 65, 29–42. doi: 10.1016/j.engappai.2017.07.003
Farahani, E., Emami, H., and Keshavarz, P. (2022). Impacts of soil organic carbon and tillage systems on structural stability as quantified by the high energy moisture characteristic (HEMC) method. Int. Agrophys. 1, 13–26. doi: 10.31545/intagr/145805
Farewell, T. S., Truckell, I. G., Keay, C. A., and Hallett, S. H. (2011). Use and applications of the Soilscapes datasets. England: Cranfield University.
Ferreira, C. d. R., da Silva Neto, E. C., Pereira, M. G., Guedes, J. N., Rosset, J. S., and Anjos, L. H. C. (2020). Dynamics of soil aggregation and organic carbon fractions over 23 years of no-till management. Soil Tillage Res. 198:104533. doi: 10.1016/j.still.2019.104533
Field, C. B. (2012). Managing the risks of extreme events and disasters to advance climate change adaptation: Special report of the intergovernmental panel on climate change. Cambridge, UK: Cambridge University Press.
Finley, A. O., Banerjee, S., and Basso, B. (2011). Improving crop model inference through Bayesian melding with spatially varying parameters. J. Agric. Biol. Environ. Stat. 16, 453–474. doi: 10.1007/s13253-011-0070-x
Galic, M., Zgorelec, Z., and Bilandzija, D. (2019). Soil carbon dioxide emissions in winter wheat vegetation influenced by agro-ecological factors and fertilization.
Gautam, S., Mbonimpa, E. G., Kumar, S., Bonta, J. V., and Lal, R. (2015). Agricultural policy environmental eXtender model simulation of climate change impacts on runoff from a small no-till watershed. J. Soil Water Conserv. 70, 101–109. doi: 10.2489/jswc.70.2.101
Geology Data Source. (2024). Available at: https://berksgeoconservation.org.uk/geology.php.
Gibson, R. B. (2001). Specification of sustainability-based environmental assessment decision criteria and implications for determining" significance" in environmental assessment. Ottawa: Canadian Environmental Assessment Agency.
Gopal, S., Yang, Y., Bai, B., and Niculescu-Mizil, A. (2012). Bayesian models for large-scale hierarchical classification. Adv. Neural Inf. Proces. Syst. 25, 1–9.
Groth, K. M., Smith, R., and Moradi, R. (2019). A hybrid algorithm for developing third-generation HRA methods using simulator data, causal models, and cognitive science. Reliab. Eng. Syst. Saf. 191:106507. doi: 10.1016/j.ress.2019.106507
Hailegnaw, N. S., Bayabil, H. K., Berihun, M. L., Teshome, F. T., Shelia, V., and Getachew, F. (2024). Integrating machine learning and empirical evapotranspiration modeling with DSSAT: implications for agricultural water management. Sci. Total Environ. 912:169403. doi: 10.1016/j.scitotenv.2023.169403
Hájek, A. (2002). Interpretations of probability. Available at: https://plato.stanford.edu/entries/probability-interpret/?ref=jcreinhold.com.
Hales, R., and Birdthistle, N. (2023). “The sustainable development goals–SDG# 12 responsible consumption and production” in Attaining the 2030 sustainable development goal of responsible consumption and production. ed. R. Hales (Bingley, UK: Emerald Publishing Limited), 1–10.
Hazarika, S., Parkinson, R., Bol, R., Dixon, L., Russell, P., Donovan, S., et al. (2009). Effect of tillage system and straw management on organic matter dynamics. Agron. Sustain. Dev. 29, 525–533. doi: 10.1051/agro/2009024
Heckerman, D. (2008). “A tutorial on learning with Bayesian networks” in Innovations in Bayesian networks: Theory and applications. ed. D. E. Holmes (Berlin: Springer Science and Business Media), 33–82.
Heidari, M., Salaudeen, S., Arku, P., Acharya, B., Tasnim, S., and Dutta, A. (2021). Development of a mathematical model for hydrothermal carbonisation of biomass: comparison of experimental measurements with model predictions. Energy 214:119020. doi: 10.1016/j.energy.2020.119020
Hermle, S., Anken, T., Leifeld, J., and Weisskopf, P. (2008). The effect of the tillage system on soil organic carbon content under moist, cold-temperate conditions. Soil Tillage Res. 98, 94–105. doi: 10.1016/j.still.2007.10.010
Hill, R. L. (1990). Long-term conventional and no-tillage effects on selected soil physical properties. Soil Sci. Soc. Am. J. 54, 161–166. doi: 10.2136/sssaj1990.03615995005400010025x
Hohle, S. M., and Teigen, K. H. (2018). More than 50% or less than 70% chance: pragmatic implications of single-bound probability estimates. J. Behav. Decis. Mak. 31, 138–150. doi: 10.1002/bdm.2052
Howe, J. P. (2015). This is nature; this is un-nature: Reading the Keeling curve. Environ. Hist. 20, 286–293. doi: 10.1093/envhis/emv005
Hughes, A. G., Vounaki, T., Peach, D. W., Ireson, A. M., Jackson, C. R., Butler, A. P., et al. (2011). Flood risk from groundwater: examples from a chalk catchment in southern England. J. Flood Risk Manage. 4, 143–155. doi: 10.1111/j.1753-318X.2011.01095.x
Jamagne, M., and King, D. (2002). “The current French approach to a soilscapes typology” in Soil Classification. ed. M. Jamagne (Boca Raton, FL: CRC Press), 157–178.
Jones, J. W., He, J., Boote, K. J., Wilkens, P., Porter, C. H., and Hu, Z. (2011). “Estimating DSSAT cropping system cultivar-specific parameters using Bayesian techniques” in Methods Introducing System models into agricultural research. eds. L. R. Ahuja and L. Ma, vol. 2 (New York: John Wiley and Sons), 365–393.
Jones, J. W., Hoogenboom, G., Porter, C. H., Boote, K. J., Batchelor, W. D., Hunt, L. A., et al. (2003). The DSSAT cropping system model. Eur. J. Agron. 18, 235–265. doi: 10.1016/S1161-0301(02)00107-7
Jug, I., Jug, D., Sabo, M., Stipesevic, B., and Stosic, M. (2011). Winter wheat yield and yield components as affected by soil tillage systems. Turk. J. Agric. For. 35, 1–7. doi: 10.3906/tar-0909-376
Karimi, J. D., Harris, J. A., and Corstanje, R. (2021). Using Bayesian belief networks to assess the influence of landscape connectivity on ecosystem service trade-offs and synergies in urban landscapes in the UK. Landsc. Ecol. 36, 3345–3363. doi: 10.1007/s10980-021-01307-6
Keeling, C. D. (1986). Atmospheric CO2 concentrations-Mauna Loa observatory, Hawaii 1958–1986. Tennessee: Oak Ridge National Laboratory.
Khosravi, K., Phuong, T. T., Barzegar, R., Quilty, J., and Aalami, M. T. (2022). Comparing the soil conservation service model with new machine learning algorithms for predicting cumulative infiltration in semi-arid regions. Pedosphere 32, 718–732. doi: 10.1016/j.pedsph.2022.06.009
Korb, K. B. (2009). Bayesian AI introduction IEEE Computer Society and IEEE Computational Intelligence Society. Available at: https://www.abnms.org/resources.php.
Koski, T. J., and Noble, J. (2012). A review of Bayesian networks and structure learning. Math. Appl. 40:278. doi: 10.14708/ma.v40i1.278
Krauss, M., Ruser, R., Müller, T., Hansen, S., Mäder, P., and Gattinger, A. (2017). Impact of reduced tillage on greenhouse gas emissions and soil carbon stocks in an organic grass-clover ley-winter wheat cropping sequence. Agric. Ecosyst. Environ. 239, 324–333. doi: 10.1016/j.agee.2017.01.029
Land Cover Data Source. (2024). Available at: https://www.ceh.ac.uk/data/ceh-land-cover-plus-crops-2015.
Lashford, C., Lavers, T., Reaney, S., Charlesworth, S., Burgess-Gamble, L., and Dale, J. (2022). Sustainable catchment-wide flood management: a review of the terminology and application of sustainable catchment flood management techniques in the UK. Water 14:1204. doi: 10.3390/w14081204
Leuning, R., Condon, A. G., Dunin, F. X., Zegelin, S., and Denmead, O. T. (1994). Rainfall interception and evaporation from the soil below a wheat canopy. Agric. For. Meteorol. 67, 221–238. doi: 10.1016/0168-1923(94)90004-3
Li, Z. T., Yang, J. Y., Drury, C. F., and Hoogenboom, G. (2015). Evaluation of the DSSAT-CSM for simulating yield and soil organic C and N of a long-term maise and wheat rotation experiment in the loess plateau of northwestern China. Agric. Syst. 135, 90–104. doi: 10.1016/j.agsy.2014.12.006
Liu, F., and West, M. (2009). A dynamic modelling strategy for Bayesian computer model emulation. doi: 10.1214/09-BA415
Louman, B., Keenan, R. J., Kleinschmit, D., Atmadja, S., Sitoe, A. A., Nhantumbo, I., et al. (2019). “SDG 13: climate action–impacts on forests and people” in Sustainable development goals: their impacts on forests and people. ed. B. Louman (Cambridge, UK: Cambridge University Press), 419–444.
Lu, X., Lu, X., Tanveer, S. K., Wen, X., and Liao, Y. (2015). Effects of tillage management on soil CO2 emission and wheat yield under rain-fed conditions. Soil Res. 54, 38–48. doi: 10.1071/SR14300
Mahmood, S. A., Karampoiki, M., Hammond, J. P., Paraforos, D. S., Murdoch, A. J., and Todman, L. (2023). Embedding expert opinion in a Bayesian network model to predict wheat yield from spring-summer weather. Smart Agric. Technol. 4:100224. doi: 10.1016/j.atech.2023.100224
Mangalassery, S., Sjögersten, S., Sparkes, D. L., Sturrock, C. J., Craigon, J., and Mooney, S. J. (2014). To what extent can zero tillage lead to a reduction in greenhouse gas emissions from temperate soils? Sci. Rep. 4:4586. doi: 10.1038/srep04586
Maraseni, T. N., and Cockfield, G. (2011). Does the adoption of zero tillage reduce greenhouse gas emissions? An assessment for the grains industry in Australia. Agric. Syst. 104, 451–458. doi: 10.1016/j.agsy.2011.03.002
Marcot, B. G. (2012). Metrics for evaluating performance and uncertainty of Bayesian network models. Ecol. Model. 230, 50–62. doi: 10.1016/j.ecolmodel.2012.01.013
Merino, J., Caballero, I., Rivas, B., Serrano, M., and Piattini, M. (2016). A data quality in use model for big data. Futur. Gener. Comput. Syst. 63, 123–130. doi: 10.1016/j.future.2015.11.024
Mkonda, M. Y., and He, X. (2023). The influence of soil organic carbon and climate variability on crop yields in Kongwa district, Tanzania. Environ. Manag. 71, 170–178. doi: 10.1007/s00267-022-01592-0
Morgan, M.Granger and Max Henrion (1990). Uncertainty: a guide to dealing with uncertainty in quantitative risk and policy analysis. New York: Cambridge University Press.
Murindangabo, Y. T., Kopecký, M., Konvalina, P., Ghorbani, M., Perná, K., Nguyen, T. G., et al. (2023). Quantitative approaches in assessing soil organic matter dynamics for sustainable management. Agronomy 13:1776. doi: 10.3390/agronomy13071776
Namdari, A., and Li, Z. (2019). A review of entropy measures for uncertainty quantification of stochastic processes. Adv. Mech. Eng. 11:168781401985735. doi: 10.1177/1687814019857350
Natcvetova, A. (2021). The role of artificial intelligence in measuring and modelling soil organic carbon in agricultural lands. Available at: https://urn.fi/URN:NBN:fi-fe2021061838806.
Nath, C. P., Das, T. K., Rana, K. S., Bhattacharyya, R., Pathak, H., Paul, S., et al. (2017). Greenhouse gases emission, soil organic carbon and wheat yield as affected by tillage systems and nitrogen management practices. Arch. Agron. Soil Sci. 63, 1644–1660. doi: 10.1080/03650340.2017.1300657
Neil, M., Fenton, N., and Nielson, L. (2000). Building large-scale Bayesian networks. Knowl. Eng. Rev. 15, 257–284. doi: 10.1017/S0269888900003039
Netica. (2024). Testing with cases. Available at: https://www.norsys.com/tutorials/netica/secD/tut_D2.htm.
Netica. (2023). Scoring rule. Available at: https://www.norsys.com/WebHelp/NETICA/X_Scoring_Rule_Results.htm.
Newlands, N. K., Zamar, D. S., Kouadio, L. A., Zhang, Y., Chipanshi, A., Potgieter, A., et al. (2014). An integrated, probabilistic model for improved seasonal forecasting of agricultural crop yield under environmental uncertainty. Front. Environ. Sci. 2:17. doi: 10.3389/fenvs.2014.00017
Ni, Z., Phillips, L. D., and Hanna, G. B. (2011). “Exploring bayesian belief networks using netica®” in Evidence synthesis in healthcare: A practical handbook for clinicians. eds. T. Athanasiou and A. Darzi (London: Springer London), 293–318.
Nielsen, D. C., and Vigil, M. F. (2010). Precipitation storage efficiency during fallow in wheat-fallow systems. Agron. J. 102, 537–543. doi: 10.2134/agronj2009.0348
Nnadi, F. N., Chikaire, J., Echetama, J. A., Ihenacho, R. A., Umunnakwe, P. C., and Utazi, C. O. (2013). Agricultural insurance: a strategic tool for climate change adaptation in the agricultural sector. Net J. Agric. Sci. 1, 1–9.
Nunes, M. R., van Es, H. M., Schindelbeck, R., Ristow, A. J., and Ryan, M. (2018). No-till and cropping system diversification improve soil health and crop yield. Geoderma 328, 30–43. doi: 10.1016/j.geoderma.2018.04.031
Pang River. (2020). Action for the River Kennet. Available at: https://www.riverkennet.org/about-the-river/pang-catchment.
Pang Water Body. (2024). Available at: https://environment.data.gov.uk/catchment-planning/WaterBody/GB106039023300.
Parton, W. J., Ojima, D. S., Cole, C. V., and Schimel, D. S. (1994). A general model for soil organic matter dynamics: sensitivity to litter chemistry, texture and management. Quant. Model. Soil Form. Process. 39, 147–167. doi: 10.2136/sssaspecpub39.c9
Parton, W. J., Ojima, D. S., and Schimel, D. S. (2020). Models to evaluate soil organic matter storage and dynamics. InStructure and organic matter storage in agricultural soils. Boca Raton, FL: CRC Press, 421–448.
Pearce, J. (1952). Farming in Berkshire. Grass Forage Sci. 7, 37–42. doi: 10.1111/j.1365-2494.1952.tb01301
Pearl, J. (1978). An economic basis for certain methods of evaluating probabilistic forecasts. Int. J. Man Machine Stud. 10, 175–183.
Pearl, J. (2011). Bayesian networks. Available at: https://escholarship.org/uc/item/53n4f34m.
Perkin, P., and Rehman, T. (1994). Farmers' objectives and their interactions with business and lifestyles: Evidence from Berkshire, England. Rural and farming systems analysis: European perspectives, pp. 193–212.
Pinheiro, É. F., de Campos, D. V., de Carvalho, B. F., dos Anjos, L. H., and Pereira, M. G. (2015). Tillage systems effects soil carbon stock and physical fractions of soil organic matter. Agric. Syst. 132, 35–39. doi: 10.1016/j.agsy.2014.08.008
Porter, C. H., Jones, J. W., Adiku, S., Gijsman, A. J., Gargiulo, O., and Naab, J. B. (2010). Modelling organic carbon and carbon-mediated soil processes in DSSAT v4. 5. Oper. Res. 10, 247–278. doi: 10.1007/s12351-009-0059-1
Prout, J. M., Shepherd, K. D., McGrath, S. P., Kirk, G. J., and Haefele, S. M. (2021). What is a good level of soil organic matter? An index based on organic carbon to clay ratio. Eur. J. Soil Sci. 72, 2493–2503. doi: 10.1111/ejss.13012
Radford, B. J. (2007). No-tillage and conservation farming practices in grain growing areas of Queensland–a review of 40 years of development. Aust. J. Exp. Agric. 47, 887–898. doi: 10.1071/EA06204
Rahman, M. M., Aravindakshan, S., Hoque, M. A., Rahman, M. A., Gulandaz, M. A., Rahman, J., et al. (2021). Conservation tillage (CT) for climate-smart sustainable intensification: assessing the impact of CT on soil organic carbon accumulation, greenhouse gas emission and water footprint of wheat cultivation in Bangladesh. Environ. Sustain. Indic. 10:100106. doi: 10.1016/j.indic.2021.100106
Rameshwaram, P., Trill, E., O’Brien, A., Robotham, J., Scarlett, P., and Blake, J. (2023). LANDWISE NFM: How soil management and land use affects soil properties and flood risk: Results from the broad-scale and detailed field surveys. Available at: https://nora.nerc.ac.uk/id/eprint/533842.
Reijneveld, J. A., van Oostrum, M. J., Brolsma, K. M., and Oenema, O. (2023). Soil carbon check: a tool for monitoring and guiding soil carbon sequestration in farmer fields. Front. Agric. Sci. Eng. 2023:499. doi: 10.15302/J-FASE-2023499
Retnowati, A., Anantasari, E., Marfai, M. A., and Dittmann, A. (2014). Environmental ethics in local knowledge responding to climate change: an understanding of seasonal traditional calendar pranotomongso and its phenology in karst area of Gunungkidul, Yogyakarta, Indonesia. Procedia Environ. Sci. 20, 785–794. doi: 10.1016/j.proenv.2014.03.095
Rial-Lovera, K., Davies, W. P., and Cannon, N. D. (2017). Implications of climate change predictions for UK cropping and prospects for possible mitigation: a review of challenges and potential responses. J. Sci. Food Agric. 97, 17–32. doi: 10.1002/jsfa.7767
Ritchie, J. T., Godwin, D. C., and Otter, S. (1984). CERES-wheat: A user-oriented wheat yield model. Preliminary documentation. AGRISTARS publication no. YM-U3-04442-JSC-18892, East Lansing, Michigan, pp. 252.
Ritchie, J. T., Singh, U., Godwin, D. C., and Bowen, W. T. (1998). Cereal growth, development and yield. Understanding options for agricultural production. Dordrecht, The Netherlands: Kluwer Academic Publication, 79–98.
Ryan, J., Masri, S., and Singh, M. (2009). Seasonal changes in soil organic matter and biomass and labile forms of carbon as influenced by crop rotations. Commun. Soil Sci. Plant Anal. 40, 188–199. doi: 10.1080/00103620802625617
Sands, P. J., Battaglia, M., and Mummery, D. (2000). Application of process-based models to forest management: experience with PROMOD, a simple plantation productivity model. Tree Physiol. 20, 383–392. doi: 10.1093/treephys/20.5-6.383
Searchinger, T. D., Wirsenius, S., Beringer, T., and Dumas, P. (2018). Assessing the efficiency of changes in land use for mitigating climate change. Nature 564, 249–253. doi: 10.1038/s41586-018-0757-z
Sebastiani, P., and Ramoni, M. (2001). On the use of Bayesian networks to analyse survey data. Res. Official Stat. 4, 53–64.
Sharma, P., Abrol, V., Sharma, K. R., Sharma, N., Phogat, V. K., and Vikas, V. (2016). Impact of conservation tillage on soil organic carbon and physical properties–a review. Int. J. Bioresour. Stress Manage. 7, 151–161. doi: 10.23910/IJBSM/2016.7.1.1387
Shiferaw, B., Okello, J., and Reddy, V. R. (2009). Challenges of adoption and adaptation of land and water management options in smallholder agriculture: Synthesis of lessons and experiences. InRainfed agriculture: Unlocking the potential. Wallingford UK: CABI, 258–275.
Sierra, C. A., Trumbore, S. E., Davidson, E. A., Vicca, S., and Janssens, I. (2015). Sensitivity of decomposition rates of soil organic matter with respect to simultaneous changes in temperature and moisture. J. Adv. Mod. Earth Syst. 7, 335–356. doi: 10.1002/2014MS000358
Singh, D., and Singh, B. (2022). Feature wise normalization: an effective way of normalizing data. Pattern Recogn. 122:108307. doi: 10.1016/j.patcog.2021.108307
Six, J., Ogle, S. M., Jay Breidt, F., Conant, R. T., Mosier, A. R., and Paustian, K. (2004). The potential to mitigate global warming with no-tillage management is only realised when practised in the long term. Glob. Chang. Biol. 10, 155–160. doi: 10.1111/j.1529-8817.2003.00730.x
Soane, B. D., Ball, B. C., Arvidsson, J., Basch, G., Moreno, F., and Roger-Estrade, J. (2012). No-till in northern, western and South-Western Europe: a review of problems and opportunities for crop production and the environment. Soil Tillage Res. 118, 66–87. doi: 10.1016/j.still.2011.10.015
Soilscapes Data Source. (2024). Available at: https://www.landis.org.uk/soilscapes/.
Solheim, S. (2021). Timing the emergence and development of arable farming in southeastern Norway by using the summed probability distribution of radiocarbon dates and a Bayesian age model. Radiocarbon 63, 1503–1524. doi: 10.1017/RDC.2021.80
Stockmann, U., Adams, M. A., Crawford, J. W., Field, D. J., Henakaarchchi, N., Jenkins, M., et al. (2013). The knowns, known unknowns and unknowns of sequestration of soil organic carbon. Agric. Ecosyst. Environ. 164, 80–99. doi: 10.1016/j.agee.2012.10.001
Suleiman, A. A., and Ritchie, J. T. (2004). Modifications to the DSSAT vertical drainage model for more accurate soil water dynamics estimation. Soil Sci. 169, 745–757. doi: 10.1097/01.ss.0000148740.90616.fd
Supit, I. (1997). Predicting national wheat yields using a crop simulation and trend models. Agric. For. Meteorol. 88, 199–214. doi: 10.1016/S0168-1923(97)00037-3
Thogmartin, W. E. (2010). Sensitivity analysis of north American bird population estimates. Ecol. Model. 221, 173–177. doi: 10.1016/j.ecolmodel.2009.09.013
Tibshirani, R. J., and Tibshirani, R. (2009). A bias correction for the minimum error rate in cross-validation. Ann. Appl. Stat. 3, 822–829. doi: 10.1214/08-AOAS224
Tu, C., and Li, F. (2017). Responses of greenhouse gas fluxes to experimental warming in wheat season under conventional tillage and no-tillage fields. J. Environ. Sci. 54, 314–327. doi: 10.1016/j.jes.2016.09.016
Vermeulen, S. J., Aggarwal, P. K., Ainslie, A., Angelone, C., Campbell, B. M., Challinor, A. J., et al. (2012). Options for support to agriculture and food security under climate change. Environ. Sci. Pol. 15, 136–144. doi: 10.1016/j.envsci.2011.09.003
Waffa, R., and Benoit, G. (2015). Simulation of climate change impact on wheat production in the Tiaret region of Algeria using the DSSAT model. Eur. Sci. J. 9, 249–259.
Walker, S. (2014). Designing sustainability: Making radical changes in a material world. Abingdon: Routledge.
Wang, P. (2004). The limitation of Bayesianism. Artif. Intell. 158, 97–106. doi: 10.1016/j.artint.2003.09.003
Weather Data Source. (2024). Available at: https://research.reading.ac.uk/meteorology/atmospheric-observatory/.
Whitmee, S., Haines, A., Beyrer, C., Boltz, F., Capon, A. G., de Souza Dias, B. F., et al. (2015). Safeguarding human health in the Anthropocene epoch: report of the Rockefeller Foundation–lancet commission on planetary health. Lancet 386, 1973–2028. doi: 10.1016/S0140-6736(15)60901-1
Yadav, S., and Shukla, S. (2016). Analysis of k-fold cross-validation over hold-out validation on colossal datasets for quality classification. In: 2016 IEEE 6th international conference on advanced computing (IACC). IEEE. pp. 78–83.
Yu, J., Smith, V. A., Wang, P. P., Hartemink, A. J., and Jarvis, E. D. (2004). Advances to Bayesian network inference for generating causal networks from observational biological data. Bioinformatics 20, 3594–3603. doi: 10.1093/bioinformatics/bth448
Zhang, T., Su, J., Liu, C., and Chen, W. H. (2018). Bayesian calibration of the AquaCrop model. In: 2018 37th Chinese control conference (CCC). IEEE, pp. 10334–10339.
Zhang, T., Su, J., Liu, C., and Chen, W. H. (2019). Bayesian calibration of AquaCrop model for winter wheat by assimilating UAV multi-spectral images. Comput. Electron. Agric. 167:105052. doi: 10.1016/j.compag.2019.105052
Zhang, N., Zhou, X., Kang, M., Hu, B. G., Heuvelink, E., and Marcelis, L. F. (2023). Machine learning versus crop growth models: an ally, not a rival. AoB Plants 15:61. doi: 10.1093/aobpla/plac061
Keywords: Bayesian model, climate change, sustainable challenges, tillage preferences, NFM strategies, synthetic datasets, DSSAT simulations, GHG-CO2 emissions
Citation: Ali Q (2024) Bayesian model of tilling wheat confronting climatic and sustainability challenges. Front. Artif. Intell. 7:1402098. doi: 10.3389/frai.2024.1402098
Edited by:
Tarun Belwal, Texas A and M University, United StatesReviewed by:
Murat Kirisci, Istanbul University-Cerrahpasa, TürkiyeZiyi Li, University of Illinois at Urbana-Champaign, United States
Copyright © 2024 Ali. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qaisar Ali, q.ali@pgr.reading.ac.uk