Quick Start

Ecological-economic model for optimal control of fire-driven, semi-arid rangelands

Model Type
SES Library XPP model

This is an ecological-economic model that endogenizes discontinuous change between states of fire-driven, semi-arid rangelands which may exist in varying degrees as grassland and woodlands depending partly on the given soil conditions. More sandy soils will result in woodlands, whereas clay soils will more likely result in grasslands. One of the critical features of this model is the existence and impact of fires on these different types of rangelands. Grasslands may build up biomass and create a fireload that once ignited may destroy grass and woodlands. The natural recovery rates of these different rangelands are important in determining the system dynamics.  Without human control through the use of livestock, the equilibrium of the grasslands is determined by their intrinsic growth rates, competition coefficients, as well as their susceptibility to the effects of fire. Accordingly, the model also considers the behavior of fire and herbivores and the interaction between these elements in order to determine optimally controlled management parameters. The time horizon of management planning can heavily influence the pathways and conditions of the system.  Specifically, by controlling the offtake of livestock that grazes on the grasslands, planners may influence the interaction between the grasslands and woodlands to optimal paths which depend not only on natural parameters but also on the time horizons considered and discount rates. 

 

 

No scenarios available.

No animations available.

Name Initial Value Description
\(x\) - livestock biomass at time t
Name Default Value Description
\(alpha\) - growth rate of herd at time t
\(beta\) - rate of regeneration of grasses at time t
\(gamma\) - rate of depletion of grasses due to grazing
\(u\) - herd offtake at time t
$$x_{t+1} - x_t = \alpha_t x_t (1 - \gamma x_t/g_t) - u_t \\\r \r $$
$$g_{t+1}-g_{t} = \beta_t g_t (1 - c_{gg} g_t/g_{max} - c_{wg}w_t/w_{max}) - \gamma_t x_t - b(g_t) g_t\\$$
$$w_{t+1} - w{t} = \phi_t w_t (1 - c_{ww} w_t/w_{max} - c_{gw} g_t/g_{max}) - f(g_t)w_t \\$$
$$b(g_t)= \frac{\sigma g_t[|(\frac{(g_t-\gamma x_t)}{g_f} - 1)| + (\frac{(g_t - \gamma_t x_t)}{g_f}-1)]}{2(\frac{g_t-\gamma_tx_t}{g_f}-1)}\\$$
$$f(g_t)= \frac{\overline{\omega} w_t[|(\frac{(g_t-\gamma x_t)}{g_f} - 1)| + (\frac{(g_t - \gamma_t x_t)}{g_f}-1)]}{2(\frac{g_t-\gamma_tx_t}{g_f}-1)}$$
fireland.ode
Loading source code...

Download Source File (fireland.ode)

The source code file can be used with XPP/XPPAUT or other simulation tools.

Loading model data...