Savannas are dynamical systems where grasses and trees can either dominate or coexist. Fires are known to be central in the functioning of the savanna biome though their characteristics are expected to vary along the rainfall gradients as observed in
Sub-Saharan Africa. In this paper, we model the tree-grass dynamics using impulsive differential equations that consider fires as discrete events. This framework allows us to carry out a comprehensive qualitative mathematical analysis that revealed more diverse possible outcomes than the analogous continuous model. We investigated local and global properties of the equilibria and show that various states exist for the physiognomy of vegetation. Though several abrupt shifts between vegetation states appeared determined by fire periodicity, we showed that direct shading of grasses by trees is also an influential process embodied in the model by a competition parameter leading to bifurcations. Relying on a suitable nonstandard finite difference scheme, we carried out numerical simulations in reference to three main climatic zones as observable in Central Africa.
Fires and rainfall are major mechanisms that regulate woody and grassy biomasses in savanna ecosystems. Conditions of long-lasting coexistence of trees and grasses have been mainly studied using continuous-time modelling of tree-grass competition. In
these frameworks, fire is a time-continuous forcing while the relationship between woody plant size and fire-sensitivity is not systematically considered. In this paper, we propose a new mathematical framework to model tree-grass interaction that takes into account both the discrete nature of fire occurrence and size-dependent fire sensitivity (via two classes of woody plants). We carry out a qualitative analysis that highlights ecological thresholds and bifurcations parameters that shape the dynamics of the savanna-like systems within the main ecological zones. Moreover, through a qualitative analysis, we show that the impulsive modelling of fire occurrences leads to more diverse behaviors and a more realistic array of solutions than the analogous time-continuous fire models. Numerical simulations are provided to illustrate the theoretical results and to support a discussion about the bifurcation parameters and future developments.
