Maximum entropy production and plant optimization theories

2. Chemical entropy production by plants and ecosystems

The instantaneous rate of chemical entropy production (σ_chem, J K⁻¹ s⁻¹) of a system (e.g. a plant or ecosystem) within a prescribed volume V is given by (e.g. Dewar 2003) 2.1 where F_i (a vector) is the molar flux density (mol m⁻² s⁻¹) of chemical species i, μ_i (J mol⁻¹) is the chemical potential of species i, T (K) is the temperature, ν_ir is the stoichiometric coefficient of species i in reaction r (positive for products and negative for reactants) and J_r (mol m⁻³ s⁻¹) is the rate of reaction r per unit volume. In general, F_i, μ_i, T and J_r are functions of space and time. The first term in curly brackets is the local rate of entropy production due to mass flow across chemical gradients (i.e. gradients in −μ_i/T); the second term is the local rate of entropy production due to chemical reactions.

The local mass balance of chemical species i is described by the continuity equation 2.2 where ρ_i (mol m⁻³) is the molar density. When the system is in a steady state (∂ρ_i/∂t = 0), equation (2.2) gives and equation (2.1) then simplifies to 2.3 where Ω is the system boundary, and dS is the local surface element (a vector of magnitude |dS| pointing in the direction outwardly normal to the surface). Equation (2.3) only involves contributions from species i whose mass flux across the boundary (F_i) is non-zero, and σ_chem may be interpreted as the rate of entropy export by those boundary fluxes.

In §3, I focus on the chemical entropy export from plants and ecosystems associated with three carbon species—photosynthates (sugars), carbon dioxide (CO₂) and structural biomass (proteins, cellulose, etc.). I ignore other potential contributions to σ_chem in plants (e.g. due to radiative exchange, plant water transport, evaporation of liquid water, etc.). I consider the maximization of σ_chem on different time scales. At each time scale, I identify a system of ‘fast’ carbon pools that can be considered to be in an approximate steady state on that time scale, i.e. the net carbon exchange between the ‘fast’ system and its environment is approximately zero. This approximation is useful because then σ_chem can be approximated by equation (2.3) in which Ω is the boundary of the ‘fast’ system. I will assume for simplicity that T is a constant in equation (2.3) (isothermal boundary conditions).

3. Maximum chemical entropy production at different scales

(a) MEP applied to plant photosynthates and CO₂

Figure 1 schematically depicts the carbon balance of plant photosynthates (CH₂O) and CO₂. This system may be considered to be in an approximate steady state on a time scale of the order of 1 year (i.e. approximately zero net annual accumulation of CH₂O and CO₂).

• Download figure
• Open in new tab
• Download powerpoint

Figure 1.

Carbon balance of plant photosynthates (CH₂O) and CO₂. The dashed box indicates the system boundary Ω. The internal distributions of CH₂O and CO₂ within Ω are not represented. Plant structural biomass (protein, cellulose, etc.) lies outside Ω. The direction of the arrows indicates the sense in which fluxes are taken to be positive: P, photosynthesis; R, respiration for plant maintenance and growth; S, conversion of photosynthate carbon to new plant structure. Values of chemical potentials refer to the boundary Ω: μ_P, chemical potential of source photosynthate; μ_S, chemical potential of photosynthate at the sites of growth (sinks); μ_R, chemical potential of respired CO₂. In the steady state, R = P − S.

The input flux is identified with the end products of photosynthesis (P). Some of the photosynthate is converted to CO₂ during plant respiration (R) and exported to the environment; the remainder is incorporated into various carbon products (proteins, cellulose, etc.) during structural growth (S). For each carbon flux across the system boundary, the chemical potential of the associated carbon species is indicated: μ_P, chemical potential of source photosynthate; μ_S, chemical potential of photosynthate at the sites of growth (sinks); μ_R, chemical potential of respired CO₂.

Applying equation (2.3) (noting the minus sign) to the system in figure 1 (noting the flux sign convention) gives 3.1 Substituting the steady-state flux relation R = P − S into equation (3.1) then gives 3.2 If one assumes that photosynthate is an ideal solute, then μ_P = μ_ref + R_GT ln(ρ_P/ρ_ref) and μ_S = μ_ref + R_GT ln(ρ_S/ρ_ref), where R_G is the universal gas constant, ρ_P and ρ_S are the source and sink photosynthate concentrations, and μ_ref and ρ_ref are the chemical potential and concentration of photosynthate in some reference state. Similarly, for CO₂ one may assume μ_R = μ_R,ref + R_GT ln(ρ_R/ρ_R,ref). The chemical potentials μ_P, μ_S and μ_R will vary in time to some degree due to variations in ρ_P, ρ_S and ρ_R (which depend on the fluxes P, S and R). As a first approximation, I will ignore these variations and treat μ_P, μ_S and μ_R as fixed parameters. Then, from equation (3.2) and recalling that T is also assumed to be constant, we have 3.3 where 3.4 is a constant. According to the Münch hypothesis of phloem transport (e.g. Christy & Ferrier 1973), the internal movement of photosynthates between sources and sinks—represented, respectively, by the upper and lower system boundaries in figure 1—occurs from high to low concentrations (i.e. ρ_P ≥ ρ_S), implying μ_P ≥ μ_S. Also, we have μ_S ≥ μ_R since respiration involves the dissipation of high-quality substrates (CH₂O) to low-quality products (CO₂). Therefore, μ_P ≥ μ_S ≥ μ_R and so, from equation (3.4), 0 ≤ λ_C ≤ 1. Note that the second law of thermodynamics () is satisfied so long as R ≥ 0 (since R = P − S ≥ 0 and λ_C ≤ 1 imply P ≥ S ≥ λ_CS).

When λ_C = 0, equation (3.3) implies that MEP is equivalent to maximizing P, which is a realistic goal from the perspective of natural selection; however, this lower limit is physiologically and thermodynamically unrealistic (μ_S = μ_R). The upper limit λ_C = 1 (i.e. μ_P = μ_S) is possibly a reasonable approximation for small plants (small internal gradients in CH₂O concentration). In this case, MEP is equivalent to maximizing plant respiration, R = P − S. At first sight, maximizing R might seem a less obvious goal for plant survival than maximizing P because R is often viewed negatively (especially by modellers) as a carbon ‘cost’ for plant growth (since S = P − R). This view neglects the fact that R drives all the metabolic processes that are crucial to plant function and survival (including growth, S) so that maximizing R is reasonable on fitness grounds. Moreover, in the steady state, maximizing R (= P − S) cannot be sustained without also maximizing P, because photosynthesis provides the substrate for respiration. Consideration of the hypothetical limiting cases λ_C = 0 and λ_C = 1 therefore suggests that, across the entire range 0 ≤ λ_C ≤ 1, maximization of σ_{CH2O + CO2} (equation (3.3)) is reasonable from the perspective of natural selection.

To explore this suggestion in more detail, I now consider the biological implications of maximizing σ_CH2O+CO2 ∝ P − λ_CS for the more realistic intermediate values 0 < λ_C < 1. In general, P is a saturating function of plant light absorption (I) and plant nitrogen content (N) (e.g. Dewar 1996), while S is more nearly proportional to N (e.g. Ågren & Franklin 2003). It follows that there is an optimal value of N which maximizes P − λ_CS; as N increases, the ‘benefit’ of increased P is eventually offset by the ‘cost’ of increased λ_CS. It should be remembered, however, that the benefit and cost here are being interpreted in terms of chemical entropy production (=entropy export) rather than carbon gain.

As a simple example, let us assume the rectangular hyperbolic relationship P = hαIkN/(αI + kN) (Dewar 1996) and the linear relationship S = gN (I, plant light absorption; N, plant nitrogen content; h, daylength; α, quantum yield; k, carboxylation coefficient; g, nitrogen growth efficiency). Maximization of P − λ_CS then predicts that the optimal plant nitrogen content is 3.5 where . The optimal rate of structural growth is then predicted to be directly proportional to plant light interception, S_MEP = ɛ_SI, where ɛ_S = αgθ/k can be interpreted as the plant growth ‘light-use efficiency’. A linear relationship between plant growth rate and light absorption has been observed empirically over a wide range of different plant types (e.g. Dewar 1996 and references therein). The above equations imply that a similar result also applies to photosynthesis itself: P_MEP = ɛ_PI, where ɛ_P = αhθ/(1 + θ) is the photosynthetic light-use efficiency—I will use this result in §3b.

The nitrogen-based trade-off here between P and S is mathematically equivalent to the nitrogen-based trade-off between P and maintenance respiration proposed previously under the assumption that plants maximize their net primary productivity (Dewar 1996; Haxeltine & Prentice 1996). Only a verbal justification for that assumption was given—it seems reasonable for plant survival. Here, this trade-off is given a novel thermodynamic interpretation—it is the result of MEP applied to plant photosynthates and CO₂.

(b) MEP applied to whole plants

Figure 2 depicts the carbon balance of whole plants. Here, the system boundary (Ω) in figure 1 has been extended to include plant structure. This extended system may be considered to be in an approximate steady state on a time scale of 1–10 years (i.e. of the order of the lifetime of plant structural biomass, which depends on plant type). Structural growth S is now an internal flux; litter production (L) takes the place of S as an external flux, and the associated chemical potential is μ_L.

• Download figure
• Open in new tab
• Download powerpoint

Figure 2.

Whole-plant carbon balance. Notation as in figure 1, except that the system boundary (Ω) has been extended to include plant structure (proteins, cellulose, etc.). L, plant litter production; μ_L, chemical potential of plant litter. In the steady state, R = P − L.

Analogous to equation (3.1), the plant entropy export rate is given by 3.6 Substituting the steady-state flux relation R = P − L into equation (3.6) then gives a result analogous to equation (3.3) 3.7 where 3.8 Here, again, I have ignored variations in the chemical potentials μ_P, μ_L and μ_R, and the temperature T. As before, I assume that μ_P ≥ μ_R (CH₂O → CO₂ representing dissipation of chemical free energy), and also μ_L ≥ μ_P (plant structure being more reduced than sugars), so that λ_P ≥ 1. Then, the second law of thermodynamics (σ_plant ≥ 0) is satisfied so long as R ≥ (λ_P − 1)L (since R = P − L ≥ (λ_P − 1)L implies P ≥ λ_PL); this condition reflects the fact that structural growth is an active process (i.e. μ_L ≥ μ_P) that is driven by the free energy generated by respiration (in the form of non-equilibrium ATP/ADP and NADPH/NADP ratios).

The chemical potentials μ_P and μ_R may be modelled as before in terms of the respective concentrations of CH₂O and CO₂ on the system boundary Ω. The chemical potential of litter (μ_L) is well defined theoretically as a function of its chemical composition by μ_L = ∑_kν_kμ_k (ν_k and μ_k being, respectively, the fraction and chemical potential of component k). However, the practical determination of μ_L remains challenging due to the compositional complexity of biomass (e.g. Meysman & Bruers 2007 and references therein); for the purposes of this study, I simply assume that μ_L is a given constant.

I now consider a simple example of the application of MEP to equation (3.7), which leads to the prediction of an optimal leaf biomass B_l. We have already seen from the maximization of σ_CH2O+CO2 on shorter time scales (§3a) that plant photosynthesis at the optimal nitrogen content is given by P = ɛ_PI, where ɛ_P is the photosynthetic light-use efficiency and I is plant light absorption. In general, I is a saturating function of B_l due to the effect of leaf mutual shading; a simple model of this is given by the Beer–Lambert law I = I_in(1 − e^−ksB_l), where I_in is the incident radiation at the top of the canopy, s is the leaf area per unit leaf biomass (inversely related to leaf thickness) and k (a function of leaf orientation and clumping) describes the exponential extinction of light within the canopy. In contrast, plant litter production L is more appropriately modelled as a linear function of each biomass compartment, L = ∑_j m_jB_j, where m_j is the specific mortality rate of biomass compartment j.

With reference to equation (3.7), as B_l increases (all other B_j being held fixed), the entropic ‘benefit’ of increased P (via increased light absorption I) is eventually offset by the entropic ‘cost’ of increased λ_PL. The optimal leaf biomass that maximizes σ_plant is easily calculated as 3.9 The optimal rate of photosynthesis is then found to be P_MEP = ɛ_PI_in − λ_Pm_l/ks.

The interpretation of observed leaf biomass in terms of an optimal trade-off between canopy photosynthesis and leaf litter production has been proposed previously (e.g. McMurtrie 1985 and references therein). The objective function given by equation (3.7) is also mathematically similar to that adopted by Franklin (2007). In either case, only a verbal justification for the choice of objective function was given. Here, the use of equation (3.7) as an objective function for plant optimization models is given a new thermodynamic interpretation—it is the result of MEP applied at the whole-plant scale.

(c) MEP applied to ecosystems

Figure 3 depicts the carbon balance of an ecosystem. At this scale, the system boundary now includes litter and soil organic carbon. This system may be considered to be in an approximate steady state on a time scale of the order of 10–100 years (i.e. the residence time of carbon in litter and soil organic matter).

• Download figure
• Open in new tab
• Download powerpoint

Figure 3.

Ecosystem carbon balance. Notation as in figure 2, except that the system boundary (Ω) has been extended to include litter and soil organic matter. R_A, autotrophic (plant) respiration; R_H, heterotrophic (soil) respiration. In the steady state, R_A + R_H = P.

The ecosystem entropy export rate is (notation as in figure 3) 3.10 Substituting the steady-state flux relation R_A + R_H = P into equation (3.10) then yields (cf. equations (3.3) and (3.7)) 3.11 where again I have assumed that μ_P, μ_R and T are fixed. At the ecosystem scale, therefore, the steady-state chemical entropy export is associated with the overall dissipative reaction that converts photosynthates (CH₂O) to CO₂ (μ_P ≥ μ_R), the second law (σ_ecosystem ≥ 0) is automatically satisfied since P ≥ 0 and MEP is equivalent to maximizing canopy photosynthesis. Maximization of canopy photosynthesis has been proposed previously as a plant optimization goal (e.g. Field 1983; Anten 2005; McMurtrie et al. 2008). MEP provides a thermodynamic justification for maximizing P—it is the result of MEP applied at the ecosystem scale.

The application of optimization theories is traditionally confined to the plant or canopy scales, reflecting the popular view that natural selection acts uniquely at the level of individual organisms. But as the analysis here suggests, MEP provides an alternative thermodynamic interpretation of optimization theories that can be extended beyond individual plants to whole ecosystems.

As a simple example of applying MEP at this scale, recall the result P_MEP = ɛ_PI_in − λ_Pm_l/ks obtained previously by maximizing whole-plant entropy export (§3b). The maximization of P may thus be accomplished in part by maximizing ɛ_PI_in. Kleidon (2004) has demonstrated that a maximum in ɛ_PI_in exists with respect to variations in stomatal conductance, when large-scale vegetation–atmosphere feedbacks are taken into account. Specifically, increasing stomatal conductance leads, on the one hand, to increased plant CO₂ uptake (hence increased ɛ_P) and, on the other hand, to increased transpiration (hence increased cloud cover and reduced I_in at the land surface). Optimization of stomatal conductance was found to predict realistic vegetation–climate states (Kleidon 2004).

We may envisage MEP (maximum P_MEP = ɛ_PI_in − λ_Pm_l/ks) also operating through minimization of the term λ_Pm_l/ks, involving co-adaptation of leaf lifespan (affecting m_l), leaf orientation (affecting k) and leaf thickness (affecting s). Observed correlations between leaf traits (e.g. Reich et al. 1992; Wright et al. 2004) offer a fertile testing ground for MEP and other candidate optimization theories of plant function (e.g. McMurtrie & Dewar submitted). Finally, extending MEP to ecosystems also raises the possibility of predicting optimal soil characteristics (e.g. soil depth, moisture content and nutrient cycling), since the long-term maximization of photosynthesis may involve trade-offs that depend on plant–soil feedbacks.