Some Applications of Thermodynamics for Ecological Systems

Now ecologists feel necessary to construct the theoretical building of system ecology, to break strong reductionistic tradition of ecology and to include the use of thermodynamics in a new holistic approach to study ecosystems, their structure, functioning and natural history. We tried to present here the current state of thermodynamic view on ecosystems. The first law of thermodynamics proclaims constancy of the total energy of isolated system for all changes, taking place in this system: energy cannot be created or destroyed. According to the second law of thermodynamics in isolated system entropy is always increasing or remaining constant. All processes in the Universe are oriented to the equilibrium state. Nevertheless, biological systems, and, consequently, ecological systems create order from disorder, they create and support chemical and physical non-equilibrium state – the basis they live on. In this chapter the general overview of ecosystem as thermodynamic system is given and the concept of Eco-Exergy is introduced. The use of this concept in ecology is demonstrated to be very fruitful. To make it easy for other researchers to use the Eco-Exergy the procedure of exergy evaluation for ecosystems is followed with special attention to dimensions used. The main applications of exergy in modern ecology are reviewed with special focus on practical use of Eco-Exergy, exergy index, and structural exergy for real ecosystem assessment and estimation of their health and disturbance. Another application of irreversible thermodynamics (Prigogine’s inventions) is discussed. The theory of hypercycles, developed for cycles of autocatalytic reactions and widely accepted in biochemistry and molecular biology can also be applied for ecological systems. The model of conjugated hypercycles, applied to ecological systems explains many aspects of their non-linear dynamics and can be used for analysis of oscillating processes in ecological systems.

surrounding and to thermodynamic equilibrium (The First Ecological Law of Thermodynamics).10.If the ecosystem is offered more pathways or combinations of pathways to move away from thermodynamic equilibrium, then the combinations of pathways that move the system most away from thermodynamic equilibrium (=yield the highest Eco-Exergy of the ecosystem) will win (The Second Ecological Law of Thermodynamics).

Eco-exergy
The exergy of a system is a measure of its deviation from thermodynamic equilibrium with the environment, and represents the maximum capacity of energy to perform useful work as the system proceeds to equilibrium, with irreversibility increasing its entropy at the expense of exergy (Ludovisi, 2009).Taken by itself, the total exergy of an ecosystem is a measure of the change in entropy content from the equilibrium and the actual state (Svirezhev, 2000).We may distinguish between technological exergy and Eco-Exergy: technological exergy uses the environment as reference state and is useful to find the first class energy (work) that a power plant can produce, Eco-Exergy uses as reference state the same ecosystem with the same temperature and pressure but at thermodynamic -chemical equilibrium (Fig. 1).Below we use the terms exergy and Eco-Exergy as synonymic.
The development and maintenance of the far-from-equilibrium condition of ecosystems is due to the steady storage of free energy into complex organic structures, biosynthesized from simple inorganic compounds.Accordingly, the total exergy of an ecosystem actually reflects the accumulation of biomass into the system, irrespective of the distribution of biogenic matter among ecosystem components.Exergy is a measure of the free energy of a system with contributions from all components including the energy of organisms.The measure for exergy in ecology also includes a factor to weigh the "complexity" of the ecological species.Moving from macroscopic to microscopic information storage, the exergetic contribution due to information grows and becomes even three orders of magnitude higher than the physical one in the more complex living systems.The capacity of packaging information at the molecular level (DNA) that differs from one organism to another can be taken into account using Eco-Exergy function.
Fig. 1.Exergy is calculated for the system relatively to reference environment, Eco-Exergy relatively to the same system at the same temperature and pressure, but as inorganic solution without life and even organic molecules.
Here, we accept the following definition of exergy (according to Jørgensen, 1992;Svirezhev, 2000;Jørgensen & Svirezhev, 2004): Exergy is the distance between the present state of the system and the state of it in thermodynamic equilibrium with the environment, measured in the units of energy.It demonstrates the amount of work performed to create a given system from its primary components (in the case of ecological systems -from primary chemical compounds).Exergy related to the total biomass (structural, specific or normalized exergy) measures the possibility of the ecosystem to accept and utilize external fluxes of energy.It reflects the degree of ecosystem development or complexity, and has advantages in comparison with the total exergy such as independence from the total biomass of the ecosystem and possibility to serve as an indicator, demonstrating the level of evolutionary development of organisms in the ecosystem.

Eco-exergy calculation
According to Prigogine's theoreme, an entropy production in every linear system without external influences is decreasing until it reaches minimum at steady state of dynamic equilibrium.Every living system is thermodynamic open system (EL 3) continuously converting potential chemical energy of organic matter into useful energy of creative processes (EL 9) and, in the end of ends, dispose to environment in the form of heat (EL 2).
As a result of it, there is no thermodynamic equilibrium in living system.At temperatures, normal for life (see EL 7), living structures are labile and are destructed constantly.To compensate this destruction the permanent internal work in the form of synthesis is fulfilled.These working synthetic processes are processes producing negative entropy (negentropy), they create order with the use of chemical energy of low-entropy energy-rich compounds (consumers and reducers) or low-entropy energy-rich solar radiation (producers).Termination of these processes causes the loss of order, death.Dead body is in thermodynamic equilibrium with maximal entropy.Exergy is the useful part of energy involved in some process, i.e. the maximal work fulfilled by the system during transit to thermodynamic equilibrium with environment state.This equilibrium means all components to be: 1) inorganic, 2) oxidized to maximum degree, 3) distributed homogenously (there is no gradients in the system).So, if we shall transfer the system into thermodynamic equilibrium with its environment, temperature and pressure will be equal for system and environment, so the only component exergy consists of will be chemical energy.Differences in temperature and pressure between system and environment are small, so we can ignore them in our calculations.The maximal input to exergetic constituent of ecological system will be done due to chemical energy, stored in organic matter and biological structures (Jørgensen et al., 2000;Jørgensen, Fath, 2011;Jørgensen, 2011).Taking into account this concept Eco-Exergy index can be calculated basing on chemical energy: state.Eco-Exergy index for the system is calculated by reference of this system to the same system in the form of inorganic soup (i.e.-without life, structure, information, organic matter).The equation for exergy calculation was proposed by S.E.Jørgensen (Mejer & Jørgensen, 1977): where Ех -exergy, J; R -gas constant, J Mol -1 •К -1 ; T -temperature, К; с i -concentration of component i, Mol; с i,eq -cncentration of the same component in the state of thermodynamic equilibrium with environment.Mol; N -number of components.The problem is to find value of с i,eq .From one hand, exergy of compounds can be calculated on the basis of elementary composition.The disadvantage of this approach consists in that, firstly, input of inorganic components into total exergy of ecosystem is too low, secondly, higher and lower organisms have approximately the same stochiometry (EL 6), so their exergy will be equal, that is in contradiction with information constituent of exergy.From another hand, the value с i,eq can be assessed from the probability P i,eq of discovering of component i in thermodynamic equilibrium state: , , , 0 If we could find this probability, we can find the ratio of с i,eq to current concentration..As inorganic component с 0 prevails in thermodynamic equilibrium state, we can rewrite (2) as: / 0, ,, Pcc eq ie q ie q  Chemical potential of dead organic matter (i =1) can be found from classic thermodynamics as: where  -chemical potential.The difference μ i -μi ,eq is known for detritus.From (5) we see: For biological components (i =2, 3, 4, …, N) probability P i,eq is composed from probability of detritus production P 1,eq , and probability P i,а of genetic information collection to determine protein structure.Supposing this events independent: 1, ,, PPP eq ie q iа   (i ≥2).
[g m -3 ] (18) Now we are to find P i,а -the probability of creation of specific genetic information, characteristic for the organism given.Originally (Jørgensen, 1992(Jørgensen, , 2002;;Jørgensen et al., www.intechopen.comSome Applications of Thermodynamics for Ecological Systems 325 2000) P i,а was determined basing on number of informative (structural) genes (each gene codes the sequence of 700 aminoacids in average) for various taxonomic groups: where g -числ е .Then exergy of typical single cell alga (850 genes approximately) can be calculated as: If we relate values of different components of ecosystem exergy to one of detritus (7,34•10 5 ), we shall get relative recalculation coefficient β i .Corresponding coefficients were calculated for many systematic groups and published (Jørgensen, 1992;Bendoricchio & Jørgensen, 1997;Jørgensen et al., 2000).These coefficients reflect relative complexity of organisms (simpler organisms have lower β values).Later, with the use of new genetic data and some indirect methods of β values assessment, ratio of non-informative genes to total genes number and others, new list of β was composed and published (Jørgensen et al., 2005;Jørgensen, 2007).New β values are added every year (table 1).Therefore, total exergy of ecosystem, based on chemical energy of organic matter (biomass) and information, stored in living organisms (recalculating coefficient β), can be calculated as: This exergy now is often called Eco-Exergy (sometimes -exergy index) to distinguish it from physical or technological exergy (Marques et al, 2003;Jørgensen, 2006Jørgensen, , 2007)).
Another indicator of ecosystem state, based on Eco-Exergy, was proposed -structural or specific exergy (structural or specific Eco-Exergy).Structural exergy (Ex str ) is the exergy related to total biomass (Silow, 1998(Silow, , 1999(Silow, , 2006;;Xu et al., 1999Xu et al., , 2001Xu et al., , 2004Xu et al., , 2005;;Marques et al., 2003;Jørgensen, 2006a).Unlike total exergy it does not depends on biomass and it reflects the ability of ecosystem to accept and utilize the flow of energy from external sources, serving simultaneously as indicator of ecosystem development, its complexity and level of evolutionaty development of biological species composed in it.
We can measure the following aspects of an ecosystem state with the Eco-Exergy: 1) the distance from thermodynamic equilibrium, i.e. general measure of total complexity of ecosystem; 2) structure (biomass and network size) and functions (available information) of ecosystem; 3) ability of ecosystem to survive (expressed via biomass and information of system).Structural exergy reflects: 1) efficiency of energy use by organisms; 2) relative information content of ecosystem and, 3) consequently, the ability of ecosystem to regulate interactions between organisms or groups of organisms.We have seen above exergy approach was demonstrated to be very fruitful during the analysis of the application of thermodynamic principles and laws to the main fundamental concepts of ecology at the end of the XX century.The analysis of three thermodynamic laws expressions in ecological rules together with exergy analysis led to formulation of the 10 Ecological Laws, in particular the Fourth (Ecological) Law of Thermodynamics, EL9 (Patten et al., 1997;Jørgensen et al., 1999;Straškraba et al., 1999;Jørgensen, 2006b).Non-equilibrium thermodynamics models based on the concept of exergy provided a common basis for representing many aspects of ecosystem development and response to environmental impacts as a single measure (Pykh et al., 2000).The use of exergy made possible the investigation of the flows of an ecosystem in terms of exergy and to arrange the system as a hierarchically ordered sequence of systems, thermodynamically embedded in each other (Nielsen, 2000).Experiments with mathematical models supported the hypothesis that an ecosystem can coordinate the most complex behaviour in the case of high level of exergy of the systems at the edge of oscillation before entering into the chaotic situation (Mandal et al., 2007).The thermodynamic notion of exergy was shown to give better insight both to the patterns of nonlinear ecosystem behaviour and to comparison of the patterns in ecological modelling (Svirezhev & Steinbom, 2001).

Eco-exergy and aquatic ecology
There are very few researches devoted to analysis of plankton communities with the aid of exergy.The implications of body sizes of phytoplankton and zooplankton for total system dynamics by optimizing exergy as a goal function for system performance indicator with mathematical models have been analyzed (Ray et al., 2001).A structurally dynamic model based on phosphorus nutrient limitation has been developed for Lake Mogan located nearby Ankara, Turkey.Exergy was applied as a goal function to consider the dynamic adaptation and the seasonality of plankton species (e.g., size shifts) (Zhang et al, 2003a, b).
The ecosystem of the North Sea integrity was approved to be reflected in exergy capture, storage capacity, cycling, matter losses, and heterogeneity (the diatom/non-diatom ratio of planktonic algae was used) with ecosystem model.Its feasibility was assessed as an ecosystem model of the North Sea, for the Elbe plume, after prior satisfactory calibration.
The modeling effort suggested that drastic nutrient load reduction from the Elbe alone would have a limited effect on the larger German Bight: even a 60% reduction scenario would only lead to moderate changes in all five indicators (Windhorst et al., 2005).More representative and multiple are applications of exergy to benthos communities.Exergy was used in optimization models of phytobenthos (Nielsen, 1997).Exergy concept allowed the finding of the best adapted water plants species in a given environmental condition and to explain in a satisfactory way the observed distributions of them in the Lagoon of Venice, Italy (Coffaro et al., 1997).
Exergy storage was estimated for benthic communities of sandy and muddy bottoms of the North Adriatic Sea subjected to experimental disturbance, induced by means of a controlled trawl fishing haul.The results showed a decrease of local exergy content in the disturbed area, with the minimum, both in sandy and muddy bottom, one month after the experimental disturbance.The exergy of the benthic community increased to the reference level, i.e., the surrounding control area, in accordance with the proposed hypothesis on the dynamics of exergy storage during a systems' development (Libralato et al., 2006).The changes of exergy and specific exergy were studied with data of benthic macrofauna in the Mondego estuary (Western Portugal).Estimates for the exergy indices provided useful indications for the evaluation of environmental impact due to the eutrophication process (Fonseka et al., 2002).Export of exergy was estimated for benthic communities on the South-Western Atlantic Coast of France.This export was mainly composed of the migration of grazing fish during the warm season, and of cultivated bivalves during the cold season (Leguerrier et al., 2007).
In the following study a self-organizing map for patterning exergy of benthic macroinvertebrate communities of 650 sampling sites in the Netherlands, including 855 species was implemented.Using these datasets, authors have calculated exergy of five trophic functional groups for each sampling site on the basis of the biomass data.Exergy of different trophic groups responded differently to different water types reflecting characteristics of target ecosystems (Park et al., 2006).Eco-Exergy and Specific Eco-Exergy were used to characterize the state of the community during the recovery process after damage to the benthic communities caused by ecological engineering Yangtze River, China (Zhang et al., 2009).Changes of the macro-benthic community structure (Venice lagoon, Italy) over almost 70 years were pictured, showing a sharp decrease in its diversity and system efficiency, estimated with the use of exergy (Pranovi et al., 2008).

Eco-exergy as an ecosystem health indicator
The idea to use exergy as an indicator of ecosystem health was proposed by S.E.Jørgensen (1992Jørgensen ( , 1999Jørgensen ( , 2002Jørgensen ( , 2006a,b),b), but the first applications of exergy as an ecosystem health indicator were fulfilled with mathematical models.Mathematical modelling with the use of exergy was shown to be applicable to explain ecosystem reactions (Jørgensen & Padisak, 1996), and to facilitate the estimation of parameters of models.The first pioneer papers describing the application of exergy indicators for natural aquatic ecosystems were published in 1997 (Xu et al., 1997;Marques et al., 1997).In 1998 the first application of exergy analysis to the results of field and laboratory experiments with real aquatic ecosystem was published (Silow, 1998).This work was continued by few more publications (Xu et al., 1999a(Xu et al., , 1999b;;Silow & Oh, 2004).The possibility to use such parameters as structural exergy and exergy for estimation of ecosystem state and its changes under various external influences was demonstrated for real natural and experimental ecosystems.These parameters were shown to reflect the state of ecosystem and they can indicate the degree of ecosystem adaptation, decreasing when important for ecosystem functioning components were depressed or eliminated.S.E.Jørgensen (2006a) proposed to use Eco-Exergy, specific Eco-Exergy and ecological buffer capacities as Ecological indicators for ecosystem development and health assessment.

Case study; eco-exergy analysis of lake Baikal state
The first works were fulfilled with the use of mathematic models.Different sensitivity of under-ice and open water plankton communities to contaminant additions was demonstrated.This can be related to different structural exergy content in plankton community.Exergy buffer capacity seems to be a more realistic measure for pliability of ecosystem reaction to external factors than biomass buffer capacity (Silow, 1998(Silow, , 1999)).In field researches the structural exergy of benthic communities at control (pristine) site, and in the region of Baikalsk Pulp and Paper Combine wastewaters discharge region at the same depths and kind of sediments was shown to differ strongly (structural exergy in polluted area was much lower than in pristine one), while biomass and total exergy behaved in not such an expressive way (Silow & Oh, 2004;Silow. 2006).The next step was the analysis of exergy and structural exergy of plankton community response to different chemical stressors analyzed in mesocosms experiments.Results obtained with mesocosms and microcosms demonstrate structural exergy decrease in experiments proportionally to a value of the added toxicant concentration, while other parameters (biomasses of components, total biomass of community, total exergy) fluctuated (Silow & Oh, 2004;Silow. 2006).Here we present the results of exergy calculations for natural plankton community of the lake Baikal.
Yearly average values of structural exergy during 1951-1999 fluctuated around their longterm average within the limits "long-term average ± mean square deviation" (154,9±26,0) without any trends.More interesting is the picture for total eco-exergy for the same period.
It demonstrates well expressed linear trend of increase with r 2 = 0.31 (Fig. 2).
We have also analysed the long-term dynamics of exergetic parameters for four limnological seasons at Baikal: inverted stratification (limnological Winter, under-ice season, February -April), spring overturn (limnological Spring, ice melting, May -June), direct stratification (limnological Summer, July -October), fall overturn (limnological Autumn, November -January)1 .Analysis of eco-exergy and structural exergy behaviour during different seasons cleared that the positive trend of eco-exergy is observed during limnological Summer (Fig. 3, 4).Dynamics of pelagic plankton biomass in Baikal for 1951-1999 is given in Fig. 2.There is neither expressed directional change of total biomass, nor changes of biomasses of different components (only slight positive trend of zooplankton biomass).Long-term oscillations of individual components are easily observed.Taking into account all discussed above and remembering the relative constancy of the total biomass of pelagic plankton, we can try to explain the tendency of its exergy to increase according to three listed above strategies (EL8 -increase of biomass, increase of network, increase of information).According to the first strategy it is the primary production increase, based on the mass development of small sized alga in summer period.Actually it is observed in the lake (Izmesyeva & Silow, 2010).
According to the second strategy it might be some recently observed structural changes in the plankton community (Hampton et al., 2008;Moore et al., 2009;Silow, 2010), and the third startegy is realized through the growth of share of larger zooplankton, like Cladocerans (Pislegina & Silow, 2010).Total biomass of plankton community and its individual components remains constant, while the total exergy of the community tends to increase.This increase can be explained with the principles of S.E.Jørgensen (section 2.1)the principles of exergy maximization (EL9 and EL10) via the growth of solar exergy consuming capacity, sophistication of ecosystem networking and increase of ecosystem information storage (EL8).
The calculated values of structural exergy for different seasons in the lake Baikal plankton for the second half of XX century, on the basis of long-term monitoring data, fluctuate within their natural limits (long-term average ± mean square deviation) and do not demonstrate any positive or negative trends (Fig. 2, Fig. 4).It points to the lack of expressed unfavourable changes in the lake Baikal pelagic.Fig. 2. Long-term dynamics of year-average biomasses of components (mg•m -3 ), exergy (g detritus eq•m -3 ) and structural exergy of lake Baikal plankton.Dotted line -long-term average.Fig. 3. Long-term dynamics of total exergy (g detritus eq•m -3 ) of lake Baikal plankton for different seasons.As we know from above (sections 3.1, 3.2.4)structural exergy reflects ecosystem health and ability of it to withstand to external influences, e.g.human impact.It is seen from Fig. 4 average long-term under-ice structural exergy (157,2) is practically equal to long-term yearaverage (154,9), summer structural exergy (175,9) is sufficiently higher than year-average , www.intechopen.comwhile overturn seasons exergy is lower -142,5 and 139,4 for spring and fall overturns.It means summer community is much more resistant to external disturbances than under-ice one, and during overturns Baikal is especially sensitive to pollution and other kinds of human impact.It is in good concordance with previously obtained results of our experiments with mathematical models (Silow et al., 1995(Silow et al., , 2001;;Silow, 1999) and mesocosms (Silow et al., 1991;Silow & Oh, 2004).

Hypercycles
One of the characteristic features of ecosystems behaviour are cyclic changes of component biomasses and numbers.They are observed both in natural objects and in artificial ecosystems.Application of Lotka -Volterra and energy flow models to systems containing more than two components demonstrates reducing of oscillation amplitude and stabilization of components parameters at fixed levels (Limburg, 1985;Ruan, 2001;Mougi Nishimura, 2007).Some researchers make their systems to oscillate via chaotic (Stone & He, 2007) or stochastic (Abta et al., 2008) behaviour of its components.
Another fact hardly simulated with the use of existing approaches is "paradox of plankton"coexisting of two or more species in the same ecological niche (Hutchinson, 1961).To explain it researchers are forced to find any, though very small differences in ecological characteristics of these species, such as optimal temperatures or oxygen contents, growth rates, nutrient thresholds for growth, mortality (Ebenhöh, 1988), time and duration of mass development (Nikolaev, 1986), albeit the last can be not the cause but effect of coexistence.Nicolis & Prigozhin (1977) basing on the following suggections:

Description of models
where A -food, X -component biomass, N -general organic matter content in closed system, k, d -growth and death rate coefficients, obtained This equation is identical to those used for description of autocatalytic processes, where component X serves as catalyst for self-creation from substance A. Such autocatalytic and self-reproduction units can be regarded in cycles called hypercycles (Eigen & Schuster, 1979).Ecosystem also can be described as hypercycle, where every next trophic level obtains material from previous level to reproduce itself as autocatalyst.Example of such structure is given in Fig. 5. Phytoplankton obtains nutrients from bacteria to create organic matter with the use of external energy (solar irradiation).Zooplankton feeding on phytoplankton obtains organic matter for growth.Bacteriae in turn get food supply from zooplankton corpses and faeces.Of course, this scheme is idealized as bacteriae obtain food from phytoplankton extracellular products and dead phytoplankton cells, zooplankton can consume not only phytoplankton but also bacteriae etc. Nevertheless this scheme represents the most important ways of energy and matter transfer in closed ecosystem including producers, consumers and reducers.Dynamics of components is determined by following system of equations where x i -biomasses, f -growth functions, m -death functions, g -grazing function, effectiveness of energy and matter conversion (for phytoplankton -relation between bacteriae concentration and nutrients availability) between components, φ -effectiveness of grazing,  -maximum growth rate.Indices i mean: 1 -phytoplankton, 2 -zooplankton, 3bacteriae.Function parameters were calculated at the ecosystem stability condition dx i /dt=0.Also we have investigated a system including two species of phytoplankton (x 11 , x 12 ) competing for nutrient supply and two species of zooplankton (x 21 , x 22 ), and bacteriae (Fig. 6).Starting biomasses for these newly introduced species at stability state were x 12 =0,33 where  -competition for nutrients function.

Behaviour of models
Dynamics of model ( 25) after external influence shows its returning to the stability point.Such behaviour is characteristic for stable non-linear systems (Gnauck, 1979).
There are no oscillations in this system, it always returns to stable state after initial biomasses changes.To make the system oscillate it is necessary to imitate input of nutrients or toxicants into it (Fig. 7).We can remind similar situation obtained by group of R. Pal et al. (2009).Their phytoplankton -zooplankton -nutrients model (with much more sophisticated mathematics, than ours) demonstrated oscillations under toxification.In other works oscillations of ecosystem components were caused by externally driven forces (oscillating environment) (Eladydi & Sacker, 2006;Koszalka et al., 2007).Hypercycles with not more than three components are shown to remain stable with equilibrium concentrations of components regardless with initial concentrations (Köppers, 1985).Model (26) demonstrates oscillation behaviour around stability point but never reaches it (Fig. 8).It is in good accordance, e.g. with the biomass fluctuations and species population fluctuations, cased by increase of community size and complexity (Fowler, 2009).It may be connected with the competition for resources (in our case -for nutrients, released by bacteria), as in many works similar results were obtained.For example, in system of two plant populations, competing for one nutrient (Damgaard, 2004), two predators, competing for one prey (Saleem et al., 2003;Yu et al., 2011), three microbial populations, competing for resources (Li, 2001).

Conclusion
It is now becoming clear that the movement away from thermodynamic equilibrium, and the subsequent increase in organization during ecosystem growth and development, is a result of system components and configurations that maximize the flux of useful energy and the amount of stored exergy.Both empirical data, as well as theoretical models, support these conclusions.Exergy is widely used in ecology to analyze theoretical problems and to solve applied tasks.The most perspective use of exergy parameters in recent ecology is the use of them as ecosystem health indicators.

Bacterio plankton
Exergy, and, especially structural exergy, is shown to be a good health indicator for ecosystems in many model, experimental, field and complex case studies.The application of exergy calculations to long-term dynamics of the lake Baikal plankton demonstrates the steady state of plankton community structural exergy and well observed increase of its total exergy.Presented model of two competing hypercycles shows: 1) simultaneous coexistence of two ecologically equivalent phytoplankton species obtaining nutrients from the one source; 2) auto-oscillations of all the components included in model in constant environmental conditions, similar to those observed in real ecosystems.

Fig. 4 .
Fig. 4. Long-term dynamics of structural exergy of lake Baikal plankton for different seasons.

Fig. 5 .
Fig.5.Idealized scheme of ecosystem as hypercycle.Blue arrows -flows of energy and matter between ecosystem components, green -self-reproduction of components, yellow arrow -income of low entropy (high exergy) energy from external source (Sun), dotted arrows -dissipation of energy in the form of high entropy (low exergy) heat energy.

Eco-exergy and structural exergy applications in ecology and environmental science 3.2.1 Eco-exergy in theoretical ecology and in aquatic ecology
•x 11 , x 22 =0,1•x 21 .System was described by the following equations: