Analytical Solutions
Last updated
Last updated
The balanced equation for in a single layer (Eq. [3:5.1.2a]) can be combined with Eq. [3:5.1.3a] and Eq. [3:5.1.7] and with a few assumptions allow for an analytical solution to the differential equation. The soil clay content and the C input rate were assumed to be constant and the powers and set to 1, so that the balance equation can be solved explicitly for . The time step was considered to be a year and we assumed that residues and manure fully decompose in that time frame so that . Since the influence of residues and manure on the balance is similar we assumed that C inputs are only through residues. With these substitutions Eq. [3:5.1.2a] becomes:
3:5.3.1
The constant k substitutes for . This differential equation can be solved analytically, with integration rendering the following solution:
3:5.3.2
The integration constant depends on the initial . The steady state solution for Eq. 3:5.3.2 is:
3:5.3.3
The ratio is the equilibrium that would be obtained if neither nor had a dependence on . As increases, the value of the fraction will get smaller. Therefore, the squared root term in Eq. 3:5.3.3 can be approximated as by preserving just the first two terms of the binomial expansion, from which Eq. 3:5.3.4 can be re-written as:
3:5.3.4
Therefore, as increases approaches asymptotically (Figure 1). Taking as a reference a soil layer with = 25 Mg C ha, = 0.2, and , it can be seen in Figure 1 that doubling and have a similar effect but of opposite sign such that the equilibrium increases with increasing and decreases with increasing . In both cases the increase and the decrease in are less than proportional to the increase in these two parameters. The equilibrium , however, is very sensitive to changes in , which makes this variable critical for a correct representation of dynamics. This formulation is a mathematical representation of the concept of saturation (Hassink and Whitmore, 1997; Six et al., 2002), enhanced with a control of the decomposition rate by .
The transient trajectory of is controlled by the quotient of the two exponential terms in Eq. [3:5.3.5].
[3.5.3.5]
For a given and initial , increasing the inputs () changes the steady state with decreasing marginal increments as approaches , yet the steady state condition is approached faster with higher inputs (Figure 3:5-1). For a given , changing has a substantial impact on the rate of change of when the inputs are medium to high (Figure 3:5-1) but a minor effect if inputs are too low. This formulation strongly suggests that soils with higher carbon storage capacity (higher ) that are currently depleted of should be the primary targets for storing , or that soils with low may store carbon quickly for a few years but the rate of gains will decrease earlier than in soils with higher .
The conditions for which the can be modeled analytically as shown here are very restrictive. The numerical solution implemented in the model is more flexible as the constants and are allowed to differ from 1. The model can be expanded to accommodate saturation of different SOM pools, instead of just one uniform pool, as strongly suggested by the results and analysis of Stewart et al. (2008). Yet, this will require a level of parameterization for which we consider there is simply not sufficient information for a realistic implementation in numerical models.
Figure 3:5-1. Equilibrium soil organic carbon (SC, Mg C ha) for the steady state condition (Eq. [3:5.3.3]) with different values for humification (, kg kg), SOM apparent turnover rate (, yr), and saturation soil organic carbon (, Mg C ha). The line without a symbol in both panels was arbitrarily chosen as a reference. The linear, no asymptotic line in Panel A shows the equilibrium for the case in which and do not depend on so that at equilibrium = , where is the residue carbon input rate (Mg C ha yr). Panel A shows the equilibrium when is doubled and Panel B shows the equilibrium when or is doubled.