The math behind the daisyworld modelby Ralph H. Abraham, abraham@vismath.orgThis is simple account of the mathematical analysis behind the daisyworld model, as originally published in Andrew J. Watson and James E. Lovelock, "Biological homeostasis of the global environment: the parable of Daisyworld", Tellus (1983), 35B, 284-289. We will refer to this article as "WL". The science behind the model is discussed in WL and elsewhere, see the bibliography. As indicated in the title of WL, the heart of the model is a point attractor of a dynamical scheme. In this case, the main control parameter is
The state variables are:
A = alphaGAG + alphaBAB + alphaWAWThus we have a two-dimensional dynamical system, given in equation (1) of WL, for the rates of change of alphaB and alphaW, alphaW' = alphaW(x beta - gamma)where x = alphaG, gamma is the death rate of all daisies, taken as 0.3 in the simulations, and beta is a quadratic function of the local temperature, equation (3) of WL, beta(T) = max {0, 1 - 0.003265 (22.5 - T)^2}Now we look for the critical points. Assuming that both daisy areas are positive (a zero value means the game is over) we find the conditions for a critical point, as given in equations (14) of WL, Twhich are constants independent of L, a surprising and hopeful result. From these equilibrium conditions, we find beta But to find them individually, it is necessary to proceed with numerical integration. The results of these simulations occupy the bulk of tje WL paper. Revised 06 May 2000 by Ralph Abraham. |