The purpose of the agricultural system development was the identification of the some of the most relevant causal relationships in the system plant-soil-micro-biota within the frame of agricultural production for a subsequent simulation. The agricultural system was considered based on previous agricultural soil models [21]. The fundamental elements that would permit the definition of a “generic agricultural system” were extracted. With these elements and their interactions, the microbiota element was incorporated and disaggregated from a “black box” (Figure 1). Given that the concept of agricultural management deals with a great quantity of interactions between the farmer and the culture,

the management is defined based on the input of organic matter to soil, because it allows the development of the rhizospheric soil microbiota and its expression in a simulation. Microbiota is developed inside the system around two fundamental processes, cellulose transformation into glucose and organic Nitrogen processing into mineral Nitrogen bioavailable for plants (Figure 1). Carbon metabolism or the cellulose transformation is based in the work of Choi et al. (2018). Nitrogen metabolism is based on a general conceptual frame in terms of mineralization and Kegg. These processes are defined form one aspect of agricultural management, the input of organic matter that could be defined as complex organic Carbon molecules in the form of cellulose and organic Nitrogen in the form of proteins [22].

Individual based modelling (IBM) applies the concept of characterize the simulation different individuals, characterized by variables, in this case, microorganisms characterized by attributes. Specifically, rhizospheric soil microorganisms defined or characterized by biological variables accordingly to the biology of soil microorganisms. The computational simulation of the IBM was developed and implemented in the NetLogo Platform (Wilensky, 1999); the performed experiments in the Behavior Space tool of NetLogo considered interactions between variables.

Community process meaning synergistic associations between microorganisms defined Carbon and Nitrogen transformations. Those transformations are community process, nor a unique population process. Carbon and Nitrogen functional microorganisms’ community assemblages were defined from a functional and dynamical perspective in the context of one agricultural management aspect that is the input or organic matter (Figure 1). For this reason, the Carbon consortium represents a group of microorganisms that synergistically converts cellulose into glucose considering that cellulose is one of the dominant molecules in green amendments or green soil covers. On the other hand, the input of Nitrogen, in terms of amendments is defined as organic Nitrogen, proteins and aminoacids that are mineralized by the second consortium that looks for represent a group of microorganisms that synergistically process organic Nitrogen and converts it to mineral Nitrogen; the system also considers the input of atmospheric Nitrogen.

The Carbon microbiota assemblage characteristics and attributes permits reflect the principal biological attributes (Table 1). The differentiation between the microorganisms that are part of the Carbon consortium is essentially functional, but the context is also physiological and biological. This specialized microbiota is

defined by three different kinds of microorganisms considered as the simulation individuals or agents and these microorganisms populations has its own functional dynamics that allows that the community transform recalcitrant Carbon to produce labile available Carbon for all microorganisms and is expressed as the system glucose level. This Carbon metabolism is an emergent property of the system, resembling the interaction between the Carbon consortia, their exoenzymes and the sequential transformation of the cellulose to glucose (Table 1) in the form cellulose Cellodextrinàcellobioseàglucose (Figure 1) [23].

Nitrogen consortium is proposed considering the different stages in the transformation of organic Nitrogen to mineral Nitrogen forms. This biotransformation has the potential to generate bioavailable Nitrogen for the plants in the form of nitrate. Additionally, N can enter the soil by the fixation process, that is developed by a specialized group of microorganisms called FixingBacts and resemble Nitrogen fixating bacteria (Table 2).

Ammonia and ammonium ions suffer a posterior transformation that is mediated by another group of microorganisms performing nitrification process. The microorganisms that perform the nitrification process are defined in Table 3.

Each group of agents or microorganism’s populations are characterized by attributes (Table 4) that has been expressed as NetLogo variables associated to a specific biological function. Emergent community patterns arise from the possibility for each microorganism of take labile Carbon that in turn defines an internal energy level. The internal energy level defines the actions considered for microorganisms, movement, enzyme production, glucose assimilation, reproduction and mortality. Each action is characterized by one or more variables; the definition of this biological variables (Table 4) has been considered from the work of Moore [19].

The principal characteristic that differentiates each microorganism or agent group is functional. The exoenzyme produced by each microbial group implies specificity in terms of enzyme-substrate-product (Tables 1-3). This scheme allows representing the diversity in a way that rules out the functional redundancy as new works has expressed the “taxonomic distribution” of some enzymatic groups. A group of attributes that aim to reflect the biological reality of each represented microbial community from a functional perspective defines the regulation or control of the exoenzymes by the microbiota functional groups. Exoenzymes properties are defined by the duration of the enzyme in the soil, the diffusion rate of the enzyme, the number of product particles generated per unit of enzyme, the maximum enzyme speed, the number of particles of substrate that saturates the enzyme and the cost of producing an exoenzyme unit (Table 5).

Considering the development of an agricultural system context, we defined the plant as an agent characterized by the uptake of energy based on the availability of nitrate, NO₃, in the soil. The values refer to an experiment developed in BehaviorSpace. Nitrite oxidoreductase cost is the energy that the microorganism that produces this specific enzyme loses per unit of enzyme generated. Maximum nitrate taken plants is the level of nitrate that the plant converts in new plant cells or plant growth, the nitrate energy value is the level of energy gained by the plant cells per unit of nitrate taken and the initial plants refers to the initial number of plant cells. These variables look for a representation of the variables that could affect plant growth in an agricultural system and the plant adaptation to its environment.

In Netlogo, the environment consists of a grid composed of patches; these patches, represents for the system, soil fragments or plant root cells; the importance of the patches is that emulates different conditions in the proposed environment. In the case of this system, each grid is defined as a soil fragment, and characterized by the level of the enzymes and substrates including the glucose level. The microorganisms walk around the patches looking for glucose and the specific substrate in the metabolic scheme defined. The principal factors affecting soil microbiota that were considered in the subsequent simulations were the root exudates and conservation agriculture in terms of organic matter inputs (Table 6). The patches in this simulation could also take values of the plant root cells, which means that each root cell delivers labile Carbon for microorganisms’ populations.

In addition to Table 5 and Table 6, the experiments designed in BehaviorSpace of Netlogo, and the correspondent results presented in Figures 2-6, the community composition was not modified, presenting an initial number of 250 individuals for all of the runs. Considering the impotance of interspecific competition in an ecological context represented by the rhizosphere, the composition were maintained equal but other values observed (Table 5 and Table 6) were modified. So the composition patterns should be attributed to interspecific competition.

The microbiota interacts in the agricultural system in terms of the input of organic matter, assuming that this practice allows microbiota development. This is represented as cellulose that is an abundant molecule and proteins, substrates processed by the microbiota (Figure 1).

Microbiota community dynamics

Eleven different populations (equivalent to agents in terms of IBM) have been proposed in the model, defining a community with two specialized consortia, one consortium for Carbon and the other for Nitrogen transformations. All the

populations that conform each consortium present its own dynamics, that altogether generates emergent patterns in terms of the abundance and hence the community structure. The population growth patterns suggests that the rhizospheric

soil conditions allow the community sustainable growth and that the input of only cellulose (Figures 2(a)-(c)) is not enough to maintain microorganisms’ populations.

Some of the different populations that conforms the Carbon and Nitrogen consortia get a point of exponential growth but sharply declines, nonetheless under rhizospheric conditions (Figures 3(a)-(c)), meaning presence of roots in the simulation and glucose input from the roots, a slightly different tendency is observed.

As the populations present differences, the model’s population dynamics suggest that these results are in concordance with the rhizospheric environment theory, in which plant root supplies microbiota with nutrients, allowing defined groups to improve the fitness and hence the population patterns.

To verify the behavior of the system organic Carbon substrates, cellulose and glucose, variation in the input rate of glucose and cellulose were developed along with variation in the total root length of the rhizospheric soil. Varying the input of Carbon substrates cellulose and glucose, the number of cells defined as roots in the soil and the time between substrate inputs gives the substrate levels patterns (Table 6). The results suggest that the Carbon consortia processed the cellulose considering that this substrate doesn’t accumulate in the system (Figure 4) and the microorganisms constantly use that glucose.

In the system simulation, enzymes could be considered as the principal response variables even over microorganism’s population dynamics. The enzyme behavior in the system presents a different emergent pattern as the microorganism’s dynamics (Figure 5), suggesting that the soil functional capacity depends on the microorganism population level but, depends too on other variables associated to the enzymes own identity. Variables considered in this sense are half time of the enzyme, the patch concentration of the enzymes, the biological cost of producing the enzyme for microorganisms, an energy level under which microorganisms doesn’t produce the enzymes and the diffusion rate of the enzyme [19].

The plant depends on the system nitrate generation for growth. At the same time, Nitrate is related to the process of cellulose metabolism to glucose for the microbiota energy and is related to the organic Nitrogen transformation proposed as input of proteins. Cellulose and organic Nitrogen are structural components of organic amendments, so the system considers conservative agricultural practices and its effect on the rhizospheric microbiota articulating ultimately to plant biomass generation. The variables used for plant definition (Table 5) defines different biomass production levels, nonetheless the nitrate generation is fundamental for the appearance of the observed patterns (Figure 6).

The system was validated from the plant perspective in terms of vegetal growth patterns. These patterns indicate that in general, agricultural plants presents a lineal growth followed by a stationary phase, aspects that has been observed in Chrysanthemun [24] and other plants [25]. The growth patterns have been defining as expo linear growth, is recognized as a theoretical pattern in the context of agricultural plant development, and could be affirmed that the simulation (Figure 6) also shows this biomass increment.