next up previous
Next: About the productivity function Up: Climatevegetation, and global Previous: Introduction

The model: Formulation and simplification

The proposed model is called the ``simplest'' model because it contains the minimal set of state variables to realize a global vegetation model incorporating the global carbon cycle. The climate of our hypothetical planet is described by one variable, namely, the annual average temperature of its surface. The planetary atmosphere is an isotropic one, containing one ``greenhouse'' gas - tex2html_wrap_inline1275 - ,with carbon concentration C(t). The equation for temperature will be [3, 8]:
equation36
where k is the surface heat capacity, S is the solar radiation, tex2html_wrap_inline1293 is the surface albedo, and tex2html_wrap_inline1295 is the Stephan-Boltzmann constant. The ``greenhouse'' effect is described by the function tex2html_wrap_inline1297, which is a monotonous decreasing function with saturation for tex2html_wrap_inline1299: tex2html_wrap_inline1301 , tex2html_wrap_inline1303. A good approximation is a hyperbola.

We consider the ``point'' planet to be without ocean, covered by vegetation with the density tex2html_wrap_inline1305 (in carbon units per unit surface). If tex2html_wrap_inline1307 is the total area of the planet, then tex2html_wrap_inline1309 is the total amount of carbon contained in the vegetation. (Without loss of generality we can put tex2html_wrap_inline1311).

As in the previous model we assume that the albedo tex2html_wrap_inline1293 depends on the density of vegetation N, so that tex2html_wrap_inline1293 is a monotonous decreasing function of N (see Fig. 1).

  figure39
Figure 1: Albedo tex2html_wrap_inline1293 as a function of the carbon mass in the vegetation N.

We include in the model the simplest two-compartmental submodel of the global carbon cycle, when the total carbon is allocated to the compartments: atmosphere and biota (vegetation), with corresponding concentrations C and N.

The equation for biota is:
 equation46
where P is the annual net-production of vegetation, m is the value which is inverse to the residence time of carbon in the biota, tex2html_wrap_inline1333.

The equation for atmospheric carbon is
 equation51
where e(t) is the annual anthropogenic emission.

The total amount of carbon A(t)=C(t)+N(t) has the following evolution in time:
equation56
and integrating the equation one obtains
equation64
where tex2html_wrap_inline1339. Using the equality to exclude the variable C(t) from the equations (2) and (3), we get
 eqnarray70
where
 equation79
The system is the simplest model of the biosphere of our hypothetical planet. Due to the explicit time dependence of A(t) it is a non-autonomous dynamical system. If we suppose that A(t) changes quasi-stationarily with t, i.e. is a slow process in relation to T and N then
equation82
Therefore we analyse the stationary solutions of the system with tex2html_wrap_inline1353.



next up previous
Next: About the productivity function Up: Climatevegetation, and global Previous: Introduction

Werner von Bloh (Data & Computation)
Thu Jul 13 11:24:58 MEST 2000