HUNGARIAN JOURNAL OF INDUSTRY AND CHEMISTRY Vol. 50 pp. 57–65 (2022) hjic.mk.uni-pannon.hu DOI: 10.33927/hjic-2022-09 DEVELOPING PLANT MODELS OF REDUCED COMPLEXITY BY CHEMI- CAL PROCESS ENGINEERING WAY OF THINKING MÓNIKA VARGA1 1Hungarian University of Agriculture and Life Sciences, Kaposvár Campus, Guba Sándor u. 40, Kaposvár, 7400, HUNGARY Given the increasing complexity of agricultural systems within the broader context of the bio-based circular economy, simplified and unified plant models are needed that represent the primary biomass production by solar-driven carbon- dioxide sequestration. Utilizing experiences from process systems engineering, which was originally inspired by chemical engineering, a suitable plant model is proposed. The structure of the model is generated from the process net of the underlying state and transition elements. Two special-state elements are introduced for the short-term storage of the supplied biomass to be distributed and of the uptake of nutrient-containing water, necessary for evapotranspiration and photosynthesis. The transition-oriented description of functionalities follows the essential causalities and balances of natural self-control. Implementation of the model is illustrated by a simple example. Keywords: plant model, reduced complexity, stoichiometric processes, process network, supply/demand lo- gistics, natural self-control 1. Introduction The bio-based circular economy is crucial to secure the supply of food and materials for mankind given the burden of depleting non-renewable resources and finite reservoirs. The replacement of open process systems with circular ones needs conscious engineering planning and operations while systemically overviewing the underly- ing processes. Multisectoral process networks require the coupling of sub-models from various disciplines on dif- ferent scales. In the bio-based circular economy, the pho- tosynthetic biosystems (plants) represent biomass pro- duction from solar energy, i.e. from the only external, unlimited energy resource for our planet. Accordingly, the model-based analysis, planning and operation of cul- tivated and natural plants plays an important role. Recently, motivated also by the increasingly inte- grated engineering of natural and man-made systems, in- tensive bidirectional learning has commenced between the computational modeling of natural and man-made processes whereby: • the principles of more sustainable and resilient natu- ral ecosystems can be applied in the design and con- trol of man-made systems on the one hand, while Recieved: 26 May 2022; Revised: 7 June 2022; Accepted: 7 June 2022 Correspondence: varga.monika@uni-mate.hu • the experiences of industrial systems designed by engineers can be taken into account with regard to the model-based analysis, planning and operation of agricultural cultivation on the other. Although systems of chemical process engineering played a unique role in this knowledge transfer because the underlying multidisciplinary processes were complex enough to represent various features, they were not too complex for the application of formerly applied com- putational tools. The lessons learnt from chemical pro- cess systems are still clearly important in terms of the rapidly developing model-supported problem-solving of complex agri-food and agro-environmental process sys- tems. This paper shows how chemical process engineer- ing can be used to develop plant (crop) models of limited complexity. First, some available plant models will be overviewed in brief. Two different approaches are available, namely (i) empirical (statistical) models to calculate various spe- cific sources of biomass production and (ii) the mecha- nistic (biophysical) models that describe the underlying physical, chemical and biological processes, representing causalities and balances. This overview focuses on the mechanistic crop models. Physiological models are important to further our understanding of metabolism, growth and how plants respond to environmental conditions, e.g. climate change [1]. The detailed dynamic modeling of physio- logical processes is not a novel endeavor, in fact it was https://doi.org/10.33927/hjic-2022-09 mailto:varga.monika@uni-mate.hu 58 VARGA already applied at the end of the previous century [2–5]. However, because of the increasing complexity of Water- Energy-Food-Ecosystems Nexus [6–8], besides detailed biophysical models, the systemic coupling of these bio- physical models from various disciplines is also neces- sary, which also requires flexible models to couple eco- nomic considerations [9]. Regarding the level of detail, the improved under- standing of physiological characteristics and the expected response to environmental changes require mechanistic models at both the cellular and organ level [10]. How- ever, their practical applicability, considering the avail- able data and knowledge, also requires the development of advanced coupling models. Various complex crop models are available. In a re- cent paper, eight kinds of crop models are classified and compared [11], the most important of which are as fol- lows: Agricultural Production Systems sIMulator (APSIM) is an actively evolving tool for the modeling and simu- lation of a wide range of agricultural systems, includ- ing plants, animals and soil, which also takes into con- sideration management actions and climatic effects [12]. The crop-related parts contain the detailed biophysical description concerning the phenology, biomass accumu- lation and distribution of newly synthesized biomass as well as the uptake of water and components by taking into account the related limitations and dependencies on environmental conditions. STICS is a detailed biophysical modeling tool that considers water, carbon as well as thermal and radia- tion energy balances for many (ca. 20) different crops ([13, 14]). It also clearly represents the phenological stages and the most important biophysical processes, e.g. light interception, transpiration, uptake of water and nu- trients, etc. Cropping Systems Simulation Model [15] is also a frequently used simulation model that takes into con- sideration the soil water budget, soil-plant nitrogen bud- get, crop-canopy and root growth, dry matter produc- tion, yield, residue production and decomposition as well as several management options, e.g. cultivar selection, crop rotation, irrigation, fertilization, tillage operations, residue management, etc. Recently, the rediscovery of advantages regarding the coupling of tree and crop systems by combining their models has also come to the fore. Both APSIM and STICS follow this direction. A few validated tree mod- els are available that are integrated in APSIM to be used in combination with plants [12]. On the other hand, a STICS crop model is embedded in the Hi-sAFe agro- forestry tool [16]. 2. Challenges and objective Plant models for well-defined important crops and trees are available. These detailed, specific models require a considerable set of parameters to be identified, more- over, the increasing design space of the bio-based circular economy needs simplified, approximate, unified, flexible, extensible and connectible plant models. The objective of this paper is to introduce the con- ceptual framework and experimental implementation of a unified plant model of reduced complexity. To develop the model, the following chemical process engineering- based principles were applied: • generalized unit operations; • specific stoichiometric composition of pseudo- components and other entities; • stoichiometric conservation processes based on the model-specific conservation laws, e.g. conservation of atoms in chemistry; and • demand-supply chain-like representation of material flows driven by the underlying push or pull logistics. 3. Materials and methods 3.1 Data and calculation formulae for a typical example of a plant to be modelled As an illustrative example of a man-made and operational plant ecosystem, a cultivated field of maize was modeled, where 9600 individual plants were cultivated over an area of 1600 m2. The maize-related specific data were derived from the literature [17–19]. Within the contours of the outlined system, this cul- tivated field was associated with the connected layers of soil and the compartment of air. The environmental con- ditions were taken into consideration in accordance with the data from the respective meteorological database. The initial data of the plants refer to the stage fol- lowing the sowing of the seeds when the initial biomass of the plants is contained within the seeds and sprouting has not yet occurred. The initial conditions of the seed biomass (based on estimations by experts) and its com- ponents [20, 21] are the following: Biomass = 0.0003 kg/plant C = 0.03747 kmol/kg H = 0.06999 kmol/kg O = 0.02846 kmol/kg N = 0.00154 kmol/kg P = 0.00011 kmol/kg X = 0 kmol/kg H2O = 0.118 kmol/kg O2 = 0.00066 kmol/kg Germination is an event-driven process that occurs af- ter the time-driven sowing. Sowing, which is a manage- ment process, is modelled by putting the seeds into down- flow material storage. Afterwards, in the event of the ap- propriate environmental (meteorological and hydrologi- cal) conditions, the seeds germinate resulting in the re- lease of these stored materials in accordance with the fol- lowing parameters (based on estimations by experts): Seed biomass = 0.0001 kg/pc Hungarian Journal of Industry and Chemistry DEVELOPING PLANT MODELS OF REDUCED COMPLEXITY 59 Seed rate = 6.579 × 10−7 kg/h Proportion of leaves = 0.853 Proportion of roots = 0.147 Surface area ratio of leaves = 6.667 m2/kg Surface area ratio of roots = 3.003 m2/kg In our example model, the germination period is from April 20th until May 9th. After the leaves and roots appear, resulting from the event-driven process of germination, the life processes of plants, that is, photosynthesis, growth, respiration, evap- otranspiration and uptake, start. The rate of photosynthesis is calculated by the follow- ing simplified equations [19]: ∆BiomassDry = Num Rad Ft Et ρ DT (1) where ∆BiomassDry = biomass produced, kg, Num = number of plants, pc. Rad = radiation, W/mJ2 Ft = proportion of radiation absorbed by the plants Et = radiation-use efficiency, kgdry matter/Wabsorbed radiation ρ = density, number of plants/m2, where Ft = 1 − e−ktLAIact (2) LAIact = LeavesSurf Num Aact (3) Aact = LAIratio LandSurf (4) ρ = Num Aact (5) LandSurf = surface area of land, m2 LAIratio = 1 m2 area of leaf / m2 area of land kt = 0.8 - light extinction coefficient [19] Et = 0.01409 kg/MJ, radiation-use efficiency [19] The process of evapotranspiration is calculated from the reference evapotranspiration (ET0, mm/day), deter- mined from meteorological data according to the well- known Penman-Monteith combination equation [22]. Based on this equation, the evaporation from the land and plants during a time step are calculated separately based on the following equations: ETland = Ke ETM LandSurf DT (6) ETplants = Kcb ETM Surf Num DT (7) where ETM = reference evapotranspiration, recalculated from ET0 / kmol/h ETland = land-related evapotranspiration, kmol ETplant = plant-related evapotranspiration, kmol LandSurf = surface area of the land, m2 Surf = surface area of the leaves, m2 Num = number of plants, pc Ke = 1 m−2, soil-related part of the dual crop coefficient [22] Kcb = 1 m−2, basal plant-related part of the dual crop coefficient [22] DT = the time step in hours, distinguished for the changing actual daylight or night period Growth is calculated after germination and is inter- preted as the distribution of the photosynthetic biomass between the parts of the plants. Before the time- or event- driven appearance of the product, the following ratios are applied in line with estimations by experts: Proportion of leaves = 0.853, Proportion of roots = 0.147. Afterwards: Proportion of leaves = 0.637, Proportion of products = 0.253, Proportion of roots = 0.110. Respiration is calculated for the individual parts of the plants. For all the parts, two kinds of respiration is sim- ulated: one as a given proportion of the biomass synthe- sized and the other as a given proportion of the already existing biomass. The applied equations are the follow- ing: Rleaves = K DMleaves + C Mleaves Num DT (8) Rprod = K DMprod + C Mprod Num DT (9) Rroot = K DMroots + C Mroots Num DT (10) where Rleaves = respired biomass of leaves, kg Rprod = respired biomass of product, kg Rroot = respired biomass of roots, kg DMleaves = synthesized biomass of leaves, kg DMprod = synthesized biomass of product, kg DMroots = synthesized biomass of roots, kg Mleaves = existing biomass of leaves, kg Mprod = existing biomass of product, kg Mroots = existing biomass of roots, kg K = the constant of 0.1 h−1 [22] C = the constant of 0.0001 h−1 [22] Num = number of plants, pc. DT = the time step in hours, distinguished for the changing actual daylight or night period The uptake of water, nitrogen and phosphorus (or of other optional elements) is calculated as the minimum amount: • required for evapotranspiration and photosynthesis together and • available in the soil. 50 pp. 57–65 (2022) 60 VARGA 3.2 Non-conventional methodology of Pro- grammable Process Structures Programmable Process Structures (PPS, [23–26]) have developed from its antecedent, that is, Direct Computer Mapping [27]. PPS offers automatic generation of easily extensible, connectible and combined dynamic balance- and rule-based models for the analysis, planning and op- eration of complex process systems, even beyond indus- tries that apply CIM (Chemical Integrated Manufactur- ing). These models consist of unified state and transition elements, transition-oriented representations of structures as well as locally programmable functional prototypes. The main sources of inspiration behind PPS are: 1) the general functional definition of process systems in Kalman’s State Space Model [28]; 2) the structural rep- resentation of General Net Theory [29]; 3) the concept of communicating autonomous programs in terms of Agent- Based Modeling [30]. Accordingly, PPS models are derived from two gen- eral (state and transition) "meta-agents," namely the structure of the generated state and transition elements form a net structure, moreover, the locally programmed state and transition prototypes represent the distributed functionalities in terms of Kalman’s model. In addition, PPS can consciously make a distinction between model- specific conservation laws based on additive measures and signals. In fact, PPS models can be generated from two gen- eral meta-prototypes and from the corresponding descrip- tion of the process net. The local program containing prototype elements (that are responsible for the case- specific calculations), are also derived from the same meta-prototypes. The simulation can be executed accord- ing to the connections between the actual state and tran- sition elements, accompanied by the data transfer be- tween the actual elements and their calculated prototypes. This architecture and its AI programming language-based (SWI-Prolog to be exact) implementation strongly sup- port the integration of various field- and task-specific models. 4. Results and discussion 4.1 Chemical process engineering-inspired principles and hypotheses of plant mod- els Natural and cultivated plants, including crops, vegetables, herbs, grasses and bushes, trees, etc. in a broader context, cover a wide variety of biological species embedded in a naturally occurring and partly human-controlled pro- cess system. Since many different species exist, at first glance this resembles the early stages of chemical en- gineering when individual technologies were interpreted as a system of various case-specific reactors, separators, etc. The essential invariant elements were later general- ized according to the concept of unit operations. Simi- larly, the unified, essential features of agricultural models can be formulated as “biological engineering unit opera- tions” within the complex system of the connected agro- technological, ecological and environmental systems. Considering the need for unified and simplified bi- ological, ecological and environmental engineering pro- cess units, these systems must be represented by the nec- essary and sufficient types as well as numbers of state and transition elements calculated by a limited set of gener- ally usable program prototypes. The coherent and connectible set of the underlying ‘first principles’-based mechanistic (physical, chemical and biological) models must be based on causally correct, model-specific conservation laws-based material and en- ergy balances. Considering the numerous biochemical compounds and biological objects, that is, organs, etc., synthesized from these compounds, the various typical biological units can be characterized by their specific sto- ichiometric composition that facilitates the representation of balances in accordance with the conservation of atoms and in line with chemical principles. Besides the conservation-based balances, the causally determined (driving force-based) transformations and transportations are the second pillar of process engineer- ing models. In this regard, plant-like biological process units represent a special case because the major driving force is solar radiation originating from beyond the con- tours of the system. This feature determines the unique position of plants in the bio-based circular carbon econ- omy. In fact, solar radiation-driven photosynthesis pro- duces a stoichiometric composition of biomass that sup- plies biomass in the various state elements of plants through downflow transporting short-term storage. More- over, the forces of solar radiation-driven evapotranspi- ration result in the uptake of water and dissolved nutri- ents through upflow transporting short-term storage that supplies the additional resources required for photosyn- thesis as well as removes the by-products of the energy- producing respiration. Accordingly, the essential self-control of plant life is organized by the solar radiation-driven push logistics of downflow as well as by the solar radiation-driven pull logistics of upflow. The daily and seasonal changes in plant behavior are determined by the temporally changing environmental functionalities, while human intervention can be taken into consideration by the respective manage- rial events. The hypotheses for the simplified and unified plant model can be summarized as follows: • The state elements are described by the specific biomass (or mass); the stoichiometric amounts of C, H, O, N, P and optional X atoms; as well as those originating from H2O, O2 and CO2. • The transition elements, e.g. photosynthesis, growth, respiration, evapotranspiration, uptake, etc., determine the functionalities resulting in stoi- chiometric changes in the aforementioned sources Hungarian Journal of Industry and Chemistry https://www.swi-prolog.org/ DEVELOPING PLANT MODELS OF REDUCED COMPLEXITY 61 of biomass, mass, atoms and components in the respective state elements. • The life processes of plants as self-controlled liv- ing systems can be characterized by (i) the supply logistics of the photosynthesis-driven utilization of CO2 from air and H2O, N, P, etc. (from top soil) to produce O2 which is emitted into the atmosphere in addition to stoichiometric pools of C, H, O, N and P that is incorporated into downflow material storage, as well as by (ii) the demand logistics of the solar energy-driven evapotranspiration-controlled uptake of H2O, N and P from the soil and the emission of CO2 and H2O into the air. 4.2 Structure of the investigated process sys- tem The process net structure of the simplified plant model, embedded in its natural environment and extended with human managerial interventions, is illustrated in Fig. 1. In the net model, the dots and bars represent the state and transition elements, respectively. The state and transition elements of the simplified plant model are the following: Plant-related model elements: • State elements: – roots (responsible for the water uptake, trans- portation of dissolved nutrients and long-term biomass storage that is also capable of gener- ating useful products); – leaves (including stems which are responsible for solar radiation-driven photosynthesis and evapotranspiration); – products (which facilitate the storage of biomass for reproduction that also generates useful products); – downflow of material (short-term storage of photosynthesized biomass to be distributed amongst the roots, leaves and products); – upflow of material (short-term storage of up- taken water and nutrients as well as of respired components to be distributed between evapo- transpiration and photosynthesis). • Transition elements: – photosynthesis: utilizes solar radiation to syn- thesize biomass from atmospheric carbon dioxide, uptaken water and dissolved compo- nents, e.g. nitrogen, phosphorus, etc.; – growth: distributes the photosynthesized biomass between the parts of the plant ac- cording to the phenological phase-specific stoichiometry; – respiration: creates energy to synthesize tis- sues from already synthesized biomass and in part maintain existing plant biomass; – evapotranspiration: which is determined by the atmospheric conditions, i.e. level of solar en- ergy, generates the driving force for the uptake of water and dissolved nutrients from the soil as well as releases the CO2 and H2O produced by respiration; – uptake: supplies the necessary water and dis- solved nutrients from the soil. Soil-related model elements: • State elements: – residue (only in the topsoil): contains organic residues, e.g. from leaf littering or the plough- ing of roots; – humus (only in the topsoil): transformed or- ganic biomass in the soil; – solution containing water and dissolved com- ponents; – inorganic solid phase. • Transition elements: – transform (only in the topsoil): describes the production of humus and dissolved nutrients from the residues; – air_land (only in the topsoil): calculates the levels of precipitation and nitrogen fixation from the atmosphere into the soil as well as those of evaporation and CO2 emission from the soil into the atmosphere; – miner_deminer: determines the degrees of mineralization and demineralization of dis- solved components in the soil; – seepage: calculates the vertical downflow of water and dissolved components between the layers of soil. Human interventions in terms of cultivation: • Typical state elements: manure, seeds, harvested products, etc. • Typical time- and/or event-driven transition ele- ments: manuring, sowing, harvesting, ploughing, etc. Other environmental state elements: the atmosphere (air), solar radiation-related meteorology, ground layer below the soil that absorbs water and nutrients. 4.3 Solar radiation-driven “natural supply lo- gistics” of plant biomass generation The essential functionalities of the investigated process system can be represented as solar energy-driven as well as predominantly self-regulated, natural supply-and- demand logistics. Moreover, the respective supply-and- demand processes are causally connected that establish a natural (basically cooperative) feedback between each other. Solar radiation facilitates the synthesis of biomass, the latter is calculated according to Eqs. 1-5. The related subprocesses of the plant model are denoted in green lines in Fig. 1. The synthesis of biomass is driven by solar radiation but limited by the available amounts of water, nitrogen and phosphorus with regard to the upflow material to be 50 pp. 57–65 (2022) 62 VARGA Figure 1: Process network of the simplified plant model stored. The rate of photosynthesis also depends on the surface area of the leaves and stem. Controlled by these conditions and limitations, the synthesis sequestrates the calculated amount of the practically unlimited supply of atmospheric CO2 in the biomass. Simultaneously, pho- tosynthesis emits O2 into the atmosphere as a result of the oxidation of uptaken water, while the associated hy- drogen is incorporated into the synthesized biomass. The photosynthesized biomass supplies the downflow of ma- terial to be stored according to the plant-specific atomic stoichiometry. In the case of the studied maize, the sto- ichiometry of the synthesized dry matter was as fol- lows [20, 21]: [C, H, O, N, P, X] = [0.037, 0.062, 0.029, 0.00087, 4.8 × 10−5, 0] The dry matter is supplied by an additional amount of H2O, according to the average water content of the plant (in our case 0.682). The synthesized biomass is divided between the leaves, roots and product state elements of the plant, ac- cording to the plant-specific ratios that also depend on the phenological condition of the model. Furthermore, differ- ent stoichiometries can be used by the various parts of the plant. In addition to the increase in the amount of synthe- sized biomass, a given proportion of that which is stored in the short term is utilized to meet the energy demand supplied by respiration for the synthesis of plant biomass according to Eqs. 8-10. 4.4 Solar energy-driven “natural demand lo- gistics” of water and nutrient uptake Considering the “natural demand logistics,” the solar en- ergy facilitates the uptake of water and dissolved nutri- ents from the soil as well as of H2O and CO2, the by- products of respiration. The respective evapotranspiration is calculated according to the Penman-Monteith combi- nation equation [22] as well by Eqs. 6-7. The related sub- processes of the plant model are denoted in red lines in Fig. 1. Depending on the meteorological conditions, the H2O and CO2 content of the uptaken material to be stored is emitted into the atmosphere by the process of evapotran- spiration. Simultaneously, the by-products of respiration, namely H2O and CO2, accumulate here, moreover, the necessary amounts of H2O as well as of dissolved nitro- gen and phosphorus are utilized for photosynthesis from this store. In the knowledge of the resultant actual condi- tions with regard to the uptaken material to be stored, the uptake is controlled by demand-determined pull actions according to the concentration bounds and uptake rules as follows: Lower and upper bounds for H2O as well as N and P atoms are known. The rule that is applied is the following: if Actual ≤ Lower, then Demand = min((Upper-Actual), Available), otherwise Demand = 0, where Actual = actual amount of uptaken material to be stored, Lower = lower bound, Hungarian Journal of Industry and Chemistry DEVELOPING PLANT MODELS OF REDUCED COMPLEXITY 63 Upper = upper bound, Demand = quantity to be uptaken, Available = available amount in the topsoil. 4.5 Generation and simulation of the PPS model The files describing the respective model are found in the “Plant” directory of a Mendeley database [31]. The process network of the example model is defined in the text file named Plant_N.pl. The initial conditions (mass, biomass, components, etc.) and parameters (coef- ficients of equations, bounds, etc.) are described in the text file entitled Plant_D.pl, which was derived using an appropriately configured MS Excel spreadsheet. Starting from these case-specific files and the general definition of state and transition meta-prototypes, the general-purpose kernel program generates the editable graphical model Plant_G.graphml. In parallel, by utilizing the meta-prototypes-based templates, modeling experts have to prepare the local programs for the respective elements of the prototypes. The locally executable programs of the prototypes are de- scribed in the file entitled Plant_G_prot.graphml. In the knowledge of the prototype elements, the sec- ond generating algorithm of the PPS kernel prepares the dynamic databases of the simulation, namely • Plant_Exp.pl containing the declaration of the Pro- log clauses describing the local programs; and • Plant_Use.pl containing the declaration of the Pro- log facts describing the case-specific elements of the model along with their initial values and parameters. The actually selected simulation results are saved in the file named Plant_Out.csv, while the data can be visual- ized using a case-specific MS Excel spreadsheet. Some examples of the simulated results are illustrated in Figs. 2-5. Fig. 2 shows the change in the total amount of biomass produced over one hectare during consecutive half days in the simulation. Time = 0 indicates when the land was sowed. Biomass begins to be produced follow- ing germination and stops at the end of the vegetation period. Fluctuations in biomass production indicate that its rate is higher during the daytime compared to at night- time. In Fig. 3, water that evaporated from the plant (blue) and from the land (red) are illustrated. Evaporation from the land follows the weather conditions much more closely and dominates during the early stages of the growing season, particularly when the leaves develop. Similarly to the production of biomass, evaporation rates are also higher in this simulation during the daytime. In Fig. 4, the biomass of the parts of one plant (leaves, stem, roots, products) can be seen. The leaves, stem and roots begin to develop at the point of germination, while products appear later on. Although the leaves, stem and products are removed at harvest time, the roots are trans- formed into residues as a result of ploughing. Of course, Figure 2: Total amount of photosynthesized biomass per hectare Figure 3: Dynamic changes in the evaporation rate of H2O from plants and the surface of the land Figure 4: Biomass of plant parts throughout the growing season the stepwise increase in biomass is insignificant in this integrated illustration. Fig. 5 shows the sequestrated CO2 over one hectare during consecutive half days of the simulation. The neg- ative values refer to the reduction in its concentration in the surrounding atmosphere. 5. Conclusion Although detailed plant models for well-defined impor- tant crops are available that require hundreds of parame- 50 pp. 57–65 (2022) 64 VARGA Figure 5: Amount of sequestrated CO2 ters to be identified, the increasing complexity of agricul- tural systems within the broader context of the bio-based circular carbon economy requires simplified and unified plant models, which can describe the primary production of biomass. The suggested conceptual framework and its experimental implementation show how the chemical process engineering principles of process units, stoichio- metric conservational processes, process networks, driv- ing force-controlled functionalities and supply/demand processes can result in a reduced, unified plant model. The structure of the model can be generated from the process net of the underlying state and transition ele- ments. The model contains two special state elements for the short-term storage of the synthesized biomass to be distributed amongst the state elements of the plant as well as for the short-term storage of the uptaken aqueous nutri- ents required for evapotranspiration and photosynthesis. These two forms of logistical storage used represent the roles of phloem and xylem in detailed biophysical mod- els. The transition-related dynamic models follow the es- sential causalities and balances of natural self-control. Solar radiation-driven photosynthesis produces a stoi- chiometric composition of biomass that, as a result of the downflow of products to be stored in the short term, in- creases the biomass in the various plant elements. More- over, solar radiation-driven evapotranspiration forces wa- ter and dissolved nutrients to be uptaken through the up- flow of products to be stored in the short term, which increases the amount of water and the additional re- sources for photosynthesis required, as well as removes the by-products of energy-producing respiration. The nat- ural self-control of plant life is organized by the solar radiation-driven push logistics of downflow as well as by the solar radiation-driven pull logistics of upflow. The suggested stoichiometric approach underlines the increasing importance of elemental analysis with regard to agriculture-related raw materials and products. The experimental PPS implementation of a simple ex- ample model illustrates the possible application of the re- duced plant model prepared according to the suggested principles inspired by process systems engineering. Acknowledgements This research was partly supported by the 2019-2.1.11- TÉT-2020-00252 program. The author is especially grateful to Béla Csukás for his valuable advice. REFERENCES [1] Marin, F. R.; Ribeiro, R. V.; Marchiori, P. E. R.: How can crop modeling and plant physiology help to un- derstand the plant responses to climate change? A case study with sugarcane, Theor. Exp. Plant Phys- iol., 2014, 26, 49–63 DOI: 10.1007/s40626-014-0006-2 [2] Whisler, F.; Acock, B,; Baker, D.; et al: Crop sim- ulation models in agronomic systems, Adv. Agron., 1986, 40, 141–208 DOI: 10.1016/S0065-2113(08)60282-5 [3] Penning de Vries, F.; Jansen, D.; Ten Berge, H.; Bakema, A.: Simulation of ecophysiological pro- cesses of growth in several annual crops. Wagenin- gen: Centre for Agricultural Publishing and Docu- mentation (Pudoc); 1989. [4] McMaster, G. S.; Morgan, J. A.; Wilhelm, W. W.: Simulating winter wheat spike development and growth. Agric. For. Meteorol., 1992, 60, 193–220 DOI: 10.1016/0168-1923(92)90038-6 [5] Salminen, H.; Saarenmaa, H.; Perttunen, J.; Sievä- nen, R.; Väkevä, J.; Nikinmaa, E.: Modelling trees using an object-oriented scheme, Math. Com- put. Model., 1994, 20, 49–64 DOI: 10.1016/0895- 7177(94)90230-5 [6] Garcia, D. J., You, F.: The water-energy-food nexus and process systems engineering: A new fo- cus, Comput. Chem. Eng., 2016, 91, 49–67 DOI: 10.1016/j.compchemeng.2016.03.003 [7] Fouladi, J.; AlNouss, A.; Al-Ansari, T.: Sus- tainable energy-water-food nexus integration and optimisation in eco-industrial parks, Com- put. Chem. Eng., 2021, 146, 107229 DOI: 10.1016/j.compchemeng.2021.107229 [8] Yoon, P. R.; Lee, S-H.; Choi, J-Y.; Yoo, S-H.; Hur, S-O.: Analysis of climate change impact on resource intensity and carbon emission in pro- tected farming systems using water-energy-food- carbon nexus, SSRN Electron. J., 2022, 184, 106394 DOI: 10.2139/ssrn.4054485 [9] Henderson, J. D.; Parajuli, R.; Abt, R. C.: Biologi- cal and market responses of pine forests in the US Southeast to carbon fertilization, Ecol. Econ., 2020, 169, 106491 DOI: 10.1016/j.ecolecon.2019.106491 [10] Poorter, H.; Anten, N. P. R.; Marcelis, L. F. M.: Physiological mechanisms in plant growth models: Do we need a supra-cellular systems biology ap- proach? Plant, Cell Environ., 2013, 36, 1673–1690 DOI: 10.1111/pce.12123 [11] Tao, F.; Palosuo, T.; Rötter, R. P.; et al. Why do crop models diverge substantially in climate impact pro- jections? A comprehensive analysis based on eight barley crop models, Agric. For. Meteorol., 2020, 281, 107851 DOI: 10.1016/j.agrformet.2019.107851 Hungarian Journal of Industry and Chemistry https://doi.org/10.1007/s40626-014-0006-2 https://doi.org/10.1016/S0065-2113(08)60282-5 https://doi.org/10.1016/0168-1923(92)90038-6 https://doi.org/10.1016/0895-7177(94)90230-5 https://doi.org/10.1016/0895-7177(94)90230-5 https://doi.org/10.1016/j.compchemeng.2016.03.003 https://doi.org/10.1016/j.compchemeng.2016.03.003 https://doi.org/10.1016/j.compchemeng.2021.107229 https://doi.org/10.1016/j.compchemeng.2021.107229 https://doi.org/10.2139/ssrn.4054485 https://doi.org/10.1016/j.ecolecon.2019.106491 https://doi.org/10.1111/pce.12123 https://doi.org/10.1016/j.agrformet.2019.107851 DEVELOPING PLANT MODELS OF REDUCED COMPLEXITY 65 [12] Holzworth, D. P.; Huth, N. I.; deVoil, P. G.; et al. APSIM – Evolution towards a new gen- eration of agricultural systems simulation, En- viron. Model. Softw., 2014, 62, 327–350 DOI: 10.1016/J.ENVSOFT.2014.07.009 [13] Brisson, N.; Mary, B.; Ripoche, D.; et al. STICS: A generic model for the simulation of crops and their water and nitrogen balances. I. Theory and param- eterization applied to wheat and corn, Agronomie, 1998, 18, 311–346 DOI: 10.1051/agro:19980501 [14] Brisson, N.; Launay, M.; Mary, B.; Beaudoin, N. (Eds.): Conceptual basis, formalisations and param- eterization of the STICS crop model (Éditions Quæ, Versailles, France) 2009 ISBN: 978-2-7592-0290-4 [15] Stöckle, C. O.; Kemanian, A. R.; Nelson, R. L.; Adam, J. C.; Sommer, R.; Carlson, B.: CropSyst model evolution: From field to regional to global scales and from research to decision support sys- tems, Environ. Model. Softw., 2014, 62, 361–369 DOI: 10.1016/j.envsoft.2014.09.006 [16] Dupraz, C.; Wolz, K. J.; Lecomte, I.; et al.: Hi- sAFe: A 3D agroforestry model for integrating dy- namic tree-crop interactions, Sustain., 2019, 11(8), 2293, DOI: 10.3390/su11082293 [17] Bergez, J. E.; Raynal, H.; Launay, M.; et al.: Evo- lution of the STICS crop model to tackle new en- vironmental issues: New formalisms and integration in the modelling and simulation platform RECORD, Environ. Model. Softw., 2014, 62, 370–384 DOI: 10.1016/J.ENVSOFT.2014.07.010 [18] Varga, M.; Gyalog, G.; Raso, J.; Kucska, B.; Csukas, B.: Programmable process structures of uni- fied elements for model-based planning and op- eration of complex agri-environmental processes. In: Bochtis D. D.; Sørensen, C. G.; Fountas, S.; Moysiadis, V.; Pardalos, P. Inf. Commun. Technol. Agric. III Decis., (Cham.: Springer) 2022, pp. 223– 249 DOI: 10.1007/978-3-030-84152-2_11 [19] van der Werf, W.; Keesman, K.; Burgess, P.; et al. Yield-SAFE: A parameter-sparse, process-based dy- namic model for predicting resource capture, growth, and production in agroforestry systems, Ecol. Eng. 2007, 29 419–433 DOI: 10.1016/j.ecoleng.2006.09.017 [20] TNO Biobased and Circular Technologies, The Netherlands. Database for the physico-chemical composition of (treated) lignocellulosic biomass, micro- and macroalgae, various feedstocks for biogas production and biochar n.d. https://phyllis.nl/ (accessed January 4, 2022) [21] Antal, M.; Allen, S. G.; Dai, X.; Shimizu, B.; Tam, M. S.; Grønli, M.: Attainment of the theoretical yield of carbon from biomass, Ind. Eng. Chem. Res., 2000, 39, 4024–4031 DOI: 10.1021/ie000511u [22] Allen, R.; Pereira, L.; Raes, D.; Smith, M.: Crop evapotranspiration - Guidelines for computing crop water requirements - FAO Irrigation and drainage pa- per 56, Chapter 2 - FAO Penman-Monteith equation, Rome: FAO; 1998. [23] Varga, M.; Prokop, A.; Csukas, B.: Biosystem mod- els, generated from a complex rule/reaction/influence network and from two functionality proto- types, BioSystems, 2017, 152, 24–43 DOI: 10.1016/j.biosystems.2016.12.005 [24] Varga, M.; Csukas, B.: Generation of extensible ecosystem models from a network structure and from locally executable programs, Ecol. Modell., 2017, 364, 25–41 DOI: 10.1016/J.ECOLMODEL.2017.09.014 [25] Varga, M.; Csukas, B.; Kucska, B.: Implementation of an easily reconfigurable dynamic simulator for re- circulating aquaculture systems, Aquac. Eng., 2020, 90, 102073 DOI: 10.1016/J.AQUAENG.2020.102073 [26] Varga, M.; Berzi-Nagy, L.; Csukas, B.; Gyalog, G.: Long-term dynamic simulation of environmen- tal impacts on ecosystem-based pond aquaculture, Environ. Model. Softw., 2020, 134, 104755 DOI: 10.1016/j.envsoft.2020.104755 [27] Csukas, B.: Simulation by direct mapping of the structural models onto executable programs, AIChE Annu. Meet., Miami, FL: AIChE, 1998, Paper 239/9. [28] Kalman, R.; Falb P.; Arbib, M.: Topics in mathemat- ical system theory (McGraw-Hill, New York, USA) 1969 ISBN: 978-0-0703-3255-3 [29] Petri, C.: Introduction to general net theory, In: Brauer, W. (Ed.): Net Theory Appl. Lect. Notes Com- put. Sci. (Heidelberg: Springer, Berlin, Germany), 1980, DOI: 10.1007/3-540-10001-6_21 [30] Abar, S.; Theodoropoulos, G. K.; Lemarinier, P.; O’Hare, G. M. P.: Agent based modelling and sim- ulation tools: A review of the state-of-art soft- ware, Comput. Sci. Rev., 2017, 24, 13–33 DOI: 10.1016/J.COSREV.2017.03.001 [31] Varga, M.: Mendeley Data, V1, 2022, linked to this paper DOI: 10.17632/nw3cjv75j5.1 50 pp. 57–65 (2022) https://doi.org/10.1016/J.ENVSOFT.2014.07.009 https://doi.org/10.1016/J.ENVSOFT.2014.07.009 https://doi.org/10.1051/agro:19980501 https://doi.org/10.1016/j.envsoft.2014.09.006 https://doi.org/10.1016/j.envsoft.2014.09.006 https://doi.org/10.3390/su11082293 https://doi.org/10.1016/J.ENVSOFT.2014.07.010 https://doi.org/10.1016/J.ENVSOFT.2014.07.010 https://doi.org/10.1007/978-3-030-84152-2_11 https://doi.org/10.1016/j.ecoleng.2006.09.017 https://phyllis.nl/ https://doi.org/10.1021/ie000511u https://doi.org/10.1016/j.biosystems.2016.12.005 https://doi.org/10.1016/j.biosystems.2016.12.005 https://doi.org/10.1016/J.ECOLMODEL.2017.09.014 https://doi.org/10.1016/J.AQUAENG.2020.102073 https://doi.org/10.1016/j.envsoft.2020.104755 https://doi.org/10.1016/j.envsoft.2020.104755 https://doi.org/10.1007/3-540-10001-6_21 https://doi.org/10.1016/J.COSREV.2017.03.001 https://doi.org/10.1016/J.COSREV.2017.03.001 https://doi.org/10.17632/nw3cjv75j5.1 Introduction Challenges and objective Materials and methods Data and calculation formulae for a typical example of a plant to be modelled Non-conventional methodology of Programmable Process Structures Results and discussion Chemical process engineering-inspired principles and hypotheses of plant models Structure of the investigated process system Solar radiation-driven ``natural supply logistics'' of plant biomass generation Solar energy-driven ``natural demand logistics'' of water and nutrient uptake Generation and simulation of the PPS model Conclusion