Search for: All records

We analyze a spatially extended version of a well-known model of forest-savanna dynamics, which presents as a system of nonlinear partial integro-differential equations, and study necessary conditions for pattern-forming bifurcations. Homogeneous solutions dominate the dynamics of the standard forest-savanna model, regardless of the length scales of the various spatial processes considered. However, several different pattern-forming scenarios are possible upon including spatial resource limitation, such as competition for water, soil nutrients, or herbivory effects. Using numerical simulations and continuation, we study the nature of the resulting patterns as a function of system parameters and length scales, uncovering subcritical pattern-forming bifurcations and observing significant regions of multistability for realistic parameter regimes. Finally, we discuss our results in the context of extant savanna-forest modeling efforts and highlight ongoing challenges in building a unifying mathematical model for savannas across different rainfall levels.

The frequency distributions can characterize the population-potential landscape related to the stability of ecological states. We illustrate the practical utility of this approach by analyzing a forest–savanna model. Savanna and forest states coexist under certain conditions, consistent with past theoretical work and empirical observations. However, a grassland state, unseen in the corresponding deterministic model, emerges as an alternative quasi-stable state under fluctuations, providing a theoretical basis for the appearance of widespread grasslands in some empirical analyses. The ecological dynamics are determined by both the population-potential landscape gradient and the steady-state probability flux. The flux quantifies the net input/output to the ecological system and therefore the degree of nonequilibriumness. Landscape and flux together determine the transitions between stable states characterized by dominant paths and switching rates. The intrinsic potential landscape admits a Lyapunov function, which provides a quantitative measure of global stability. We find that the average flux, entropy production rate, and free energy have significant changes near bifurcations under both finite and zero fluctuation. These may provide both dynamical and thermodynamic origins of the bifurcations. We identified the variances in observed frequency time traces, fluctuations, and time irreversibility as kinematic measures for bifurcations. This framework opens the way to characterize ecological systems globally, to uncover how they change among states, and to quantify the emergence of quasi-stable states under stochastic fluctuations.

Tree clusters in savannas are commonly found in sizes that follow power laws with well-established exponents. We show that their size distributions could result from the space–time probabilistic structure of soil moisture, estimated over the range of rainfall observed in semiarid savannas; patterns of soil moisture display islands whose size, for moisture thresholds above the mean, follows power laws. These islands are the regions where trees are expected to exist and they have a fractal structure whose perimeter–area relationship is the same as observed in field data for the clustering of trees. When the impact of fire and herbivores is accounted for, as acting through the perimeter of the tree clusters, the power law of the soil moisture islands is transformed into a power law with the same exponents observed in the tree cluster data.

Understory fires represent an accelerating threat to Amazonian tropical forests and can, during drought, affect larger areas than deforestation itself. These fires kill trees at rates varying from < 10 to c. 90% depending on fire intensity, forest disturbance history and tree functional traits. Here, we examine variation in bark thickness across the Amazon. Bark can protect trees from fires, but it is often assumed to be consistently thin across tropical forests. Here, we show that investment in bark varies, with thicker bark in dry forests and thinner in wetter forests. We also show that thinner bark translated into higher fire‐driven tree mortality in wetter forests, with between 0.67 and 5.86 gigatonnes CO₂lost in Amazon understory fires between 2001 and 2010. Trait‐enabled global vegetation models that explicitly include variation in bark thickness are likely to improve the predictions of fire effects on carbon cycling in tropical forests.