Next Article in Journal
A Method of Detecting Candidate Regions and Flames Based on Deep Learning Using Color-Based Pre-Processing
Next Article in Special Issue
Tree Species Composition and Diversity in Fire-Affected Areas of Miombo Woodlands, Central Mozambique
Previous Article in Journal
Leadership at the Extreme: A Longitudinal Study of Transformational Leadership Style and Well-Being in Firefighters
Previous Article in Special Issue
Topographic Factors Drive Short-Term Understory Revegetation in Burned Areas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Demographic Response of Grass Species to Fire Treatments in a Guinean Savanna

by
Kouamé Fulgence Koffi
1,2,*,
Aya Brigitte N’Dri
1,
Tharaniya Srikanthasamy
2,
Jean-Christophe Lata
2,3,
Souleymane Konaté
1,
Marcel Konan
1 and
Sébastien Barot
2
1
Research Training Unit in Natural Sciences, Environment and Sustainable Development Research Pole, Laboratory of Ecology and Sustainable Development, Ecology Station of Lamto/CRE, Nangui Abrogoua University, 02 BP 801 Abidjan 02, Côte d’Ivoire
2
Institute of Ecology and Environmental Sciences-Paris, Sorbonne Université, CNRS, IRD, INRAE, Université de Paris, UPEC, 4 place Jussieu, 75005 Paris, France
3
Department of Geoecology and Geochemistry, Institute of Natural Resources, Tomsk Polytechnic University, 30 Lenin Street, Tomsk 634050, Russia
*
Author to whom correspondence should be addressed.
Fire 2022, 5(6), 193; https://doi.org/10.3390/fire5060193
Submission received: 15 September 2022 / Revised: 28 October 2022 / Accepted: 8 November 2022 / Published: 16 November 2022
(This article belongs to the Special Issue Effects of Wildfire on Biodiversity)

Abstract

:
Fighting tree encroachment using fire promotes C4 perennial grasses but likely affects their demography according to the fire date during the dry season. We analyzed the impact of four fire treatments (early, middle, late and no fire) on the demography of the four dominant perennial grasses (Andropogon canaliculatus, Andropogon schirensis, Hyparrhenia diplandra and Loudetia simplex) in a Guinean savanna of West Africa (Lamto, Côte d’Ivoire). We carried out a yearly demographic monitoring of each grass individual during five years (2015–2019) on three plots by treatment and parametrized a size-classified matrix model with five circumference classes. The results showed that A. schirensis, H. diplandra and L. simplex declined (λ < 1.0) under late fire and will disappear after 10 years, as did L. simplex under the middle fire. Stasis influenced the most λ values and stable class distribution was nearly achieved in all species under all treatments. The size of L. simplex should increase under early and late fires. Our results suggest that late fire is the most detrimental fire regime for grasses in this Guinean area, contrary to early and middle fires, which could be recommended to savanna managers.

1. Introduction

Savannas are defined by the coexistence of trees and grasses maintained by disturbances such as fires, when the annual rainfall is above 650 mm [1]. In fact, fire permits the fighting of tree encroachment which threaten savannas and their biodiversity, by preventing an increase in tree cover and colonization by forest plants. To assess the effectiveness of fire in fighting tree encroachment, many studies have addressed its impact on tree demography [2,3,4], while studies on grass demography are rare (but see Garnier and Dajoz [5]). However, in humid savannas, perennial grasses dominate the grass stratum because they represent 70–98% of the total above-ground biomass [6], which make them the main source of fuel for fire [7,8]. The biomass of grass and its dryness change throughout the dry season; fire characteristics also change depending on its timing. Therefore, there are currently debates on the impact of this timing on savanna functioning and the best fire period to manage vegetation and biodiversity [9,10,11].
It is well known that savanna trees and grasses are adapted to fire [12,13] with fairly well-known underlying mechanisms. For example, trees have developed several mechanisms allowing them to persist as avoiders (species that grow quickly in the absence of fire), resisters (species that develop thick bark) or resprouters (species that regenerate from their underground parts protected from fire) [14,15,16]. Grasses persist by allocating reserves to their belowground system and producing aerial buds that insure their regrowth [17]. Thus, the mechanism of grass adaptation to fire is similar to the strategy employed by resprouting trees, as they both lose their aerial part during fire while keeping their underground parts, which allows growth to resume after fire [18].
Because fire intensity depends on the period of burning [19,20], this period should have contrasting effects on the demography of savanna trees [21,22] and perennial grasses [5,23,24]. Indeed, despite their adaptation to fire, fire regularly kills them, especially if it occurs at the beginning of their development cycle at the beginning of the wet season [24,25]. Moreover, fires occurring at the beginning of the wet season burn a very dry biomass and are more intense [26,27]. Thus, fighting tree encroachment may involve fire management, which implies choosing the best fire period to maintain the abundances and diversity of grasses and the desired tree density. Because tree encroachment is a worldwide issue in savannas [28], people could advise fire regimes such as late fires that are detrimental to trees [22,25]. However, are these regimes favorable to grasses? For these reasons and to determine the best timing for fires to conserve savanna biodiversity, a team of savanna managers started test four fire treatments (the absence of fire and three fire timings) in the Lamto savanna in 2013. These treatments are the early fire set at the beginning of the dry season in November, the middle fire in January (in the middle of the dry season), the late fire in March (at the end of the dry season) and the no-fire treatment. Using matrix models, we studied here the impact of these four fire treatments on the demography of the four dominant perennial grass species.
Koffi et al. [29] have already compared the demography of the four species for the middle fire. This study showed a difference in demography between species and years, depending on variations in stasis in different size-classes and on the species. One species declined with λ < 1.0 over three one-year transitions while the others were apparently stable. Our goal here was to determine the impact of the four fire treatments on the demography of the savanna perennial grasses. We parameterized the matrix demographic models for the four dominant grass species for which enough data was available to estimate transition rates [30] in the four fire treatments to compare the response of the grass species to these fire treatments. These models allowed us to compare various integrative demographic features (such as the asymptotic growth rate or the mean age within a stage) between the fire treatments. We tested the following hypotheses: (1) The late fire has a negative effect on the demography of all species because it occurs when grasses have restarted their annual development so that their leaves and buds are vulnerable to fire. (2) The no-fire treatment negatively affects the demography of the species by killing the seedlings and some adults through shading by the dry biomass (that has not been removed by any fire). (3) The early fire has a low negative effect on grass species because it has a low intensity and occurs when biomass is still wet and burns poorly. (4) As each species is effectively the same life-form but with different characteristics (e.g., tussocks of the different species have different sizes [24]), they have contrasting responses to each fire treatment.

2. Material and Methods

2.1. Study Site

Field data were collected at the Lamto Reserve in Ivory Coast (West Africa: 6°9′–6°18′ N; 5°15′–4°57′ W) between semi-deciduous forests and humid savannas [31] in the wettest end of the Guinean savanna domain. The weather is composed of a long rainy season from March to July, a short dry season in August, a short rainy season from September to November, and a long dry season from December to February. The average annual precipitation reaches 1200 mm and the mean annual temperature is about 27 °C. During the study (from 2015 to 2019), the mean annual rainfall was 1167 mm and the main temperature was 29 °C. This shows that the study years were overall slightly hotter and drier than the average values, which could influence fire characteristics and grass demography.
The vegetation of Lamto Reserve is composed of forests and savannas, with variable shrub and tree densities and a herbaceous stratum dominated by perennial grasses [6]. About 10 perennial grass species of the Andropogoneae tribe coexist in this savanna [32]. All these species are C4 and the dominant ones are: Andropogon ascinodis C.B.Cl., Andropogon canaliculatus Schumach., Andropogon schirensis Hochst. ex A. Rich., Hyparrhenia diplandra (Hack.) Stapf, Hyparrhenia smithiana (Hook.f.) Stapf and Loudetia simplex (Nees) C.E. Hubbard [6].

2.2. Study Plots

The four studied fire treatments (early, middle, late and no fire) were implemented on three replicated 3.72 ha blocks, each divided into four 50 × 100 m plots, one per fire treatment since September 2013 [20]. We monitored grass demography on a 5 × 5 m subplot delimited on each early, late and no-fire plots, and a 5 × 10 m subplot delimited on each middle fire plot (the usually used fire regime). We used 5 years of census (from 2015 to 2019) encompassing four one-year transitions (2015–2016; 2016–2017; 2017–2018 and 2018–2019).

2.3. Data Collection

The four dominant perennial grass species (Andropogon canaliculatus, Andropogon schirensis, Hyparrhenia diplandra and Loudetia simplex) were considered in this study because they had enough individuals on all plots to carry out demographic studies except L. simplex, which does not maintain in the absence of fire. However, this species was overall the third most abundant species before A. schirensis [30]. We conducted a systematic sampling of all the individuals of each species with more than five tillers (to facilitate species identification) on all plots. We permanently marked all individuals from 2015 onwards between April and May using metal labels tied to metal pegs. We measured the circumference of each tussock using a measuring tape at the ground level. All labeled tussocks were censused every year until 2019, to determine their status (dead or alive, fragmented or not, retrogressed or recruited) as in Koffi et al. [24], and the circumference (cm) at the ground level of all living tussocks was measured. The new seedlings with at least 5 tillers were also labeled from 2016 onwards during each field campaign and integrated into the data. These seedlings were distinguished from small fragments of tussocks by the absence of remaining burnt stems at the base.

2.4. Matrix Models

For these perennial grasses, size-classified matrix model is appropriate for several reasons: (1) their survival, growth and reproduction depend more on their size than their age [33,34]; (2) they exhibit clonal reproduction; and (3) their age is difficult to determine [35,36].
Thus, a size-classified matrix model was designed for each combination of the four grass species and the four fire treatments. This model was parameterised separately for the 60 combinations of four grass species, four fire treatments and four one-year transitions, minus the 4 matrices of L. simplex missing in the no-fire treatment. Average matrices of the four one-year transitions were used to derive an average matrix for each species in each fire treatment and to compare species and fire treatments (overall 15 matrices, 4 species × 4 fire treatments—1 matrix for L. simplex under the no-fire treatment). Individuals were classified into five circumference classes: 3–10 cm (the smallest tussocks measuring 3 cm), 10–20 cm, 20–35 cm, 35–50 cm and 50 cm and more. These 5 size-classes were chosen to have a good enough description of the life history and robust parameter estimation, which requires having enough individuals in each size-class. Thus, the two first size-classes, where more individuals were found, are smaller than the following three size-classes [30].
Koffi et al. [29] have already established the life-cycle diagram of the perennial grasses (Figure S1). It shows that all transitions are possible between all size-classes in these species. The transitions lead to matrices of stasis, growth, retrogression, fragmentation and birth. These matrices contain, respectively, the probabilities of staying in the same class, growing from a small class to a larger one, retrogressing from a large class to a smaller class, producing clones in a small class by individuals in a larger class and the number of seedlings produced by each class. These matrices are summed to three types of matrices: the survival matrix (stasis, growth and retrogression), the fragmentation matrix and the birth matrix, because the calculation of the age-based parameters required these three types of matrix. All these matrix types were then summed to the transition matrix, describing the contribution of each size-class to the others. Then, the survival, fragmentation and birth matrices were used for the calculation of age parameters and the transition matrix was used for all other demographic analyses. The number of individuals in each class during each field campaign was described by a five-entry vector.

2.5. Demographic Analyses

Demography is generally summarized by the equation N(t + 1) = AN(t) [37], where N(t) and N(t + 1) are vectors of abundances in each size-class at time t and t + 1, respectively. This enables us to predict the asymptotic behavior of the population, i.e., asymptotic growth rate of the population (λ, the dominant eigenvalue of the transition matrix) and the stable or predicted size structure (the left eigenvector of the matrix; [37]). The standard error and the 95% confidence intervals of λ were calculated according to Caswell [34] to test whether λ was significantly different from 1.0 or not.
The elasticity of λ in relation to the matrix parameters [37,38] was also computed as an index of the relative influence of each matrix entry on λ. In other words, this measures the contribution of each matrix entry to changes in λ [37,39].
We conducted a loglinear analysis [40] to determine whether the transition matrices depended on the species, the transition year, the fire treatment or the previous state of the individuals. The response variable was the fate (death or living in a new size-class), while the explanatory variables were the species, the fire treatments, the transition years (2015–2016; 2016–2017; 2017–2018 and 2018–2019), and the state (the initial size-classes).
We used life-table response experiment (LTRE) analyses to quantify the contribution of each transition matrix entry to differences of λ between fire treatments. Contributions were calculated using Caswell’s method [37]. To better highlight the contribution of each demographic parameter, the entries of the contribution matrix were summed in: stasis, i.e., all entries on the diagonal; fecundity, i.e., all the matrix entries on the first line minus the first; retrogression, i.e., the entries under the fecundities and above the stasis; and growth, i.e., under the first diagonal.
We used the mean matrices of the four one-year transitions for all demographic analysis and to calculate the age-based parameters for each fire treatment and species. We used the method of Cochran and Ellner [41] to calculate the age-based parameters and their standard deviation because perennial grasses exhibited clonal reproduction. We calculated (1) the mean age of individuals in each size-class, (2) the mean age of residence in each size-class, (3) the mean time to first reach a given class from the class 1, (4) the conditional remaining life-span of individuals in the given class and (5) the total conditional life-span if a class was reached [41,42].
In addition, we used Keyfitz’s metric [43] to assess the difference between the observed size-class distribution and the predicted size-class distribution, i.e., the stable size-class distribution. This metric is a measure of the distance between any two probability vectors, i.e., it provides a measure of how far a population is from the predicted distribution at the time the distribution was observed. Keyfitz’s formula allowing such analysis is:
Δ = 1 2 i n 0 , i w i
where n 0 , i is the observed proportion of individuals in size-class i and w i is the proportion expected at stable size-class distribution. This formula provides values ranging from 0.0 (no difference) to 1.0 (high difference). The ‘popbio’ package [44] was used in R software version 4.0.5 [45] for all these parameter calculations and analyses.

3. Results

3.1. Average Matrix of Species in the Fire Treatments

The total number of grass tussocks recorded for the different fire treatments were 952, 1784, 873 and 531 in 2015 (after 2 years of fire experimentation), and 1020, 1842, 500 and 456 in 2019 (after 6 years of fire experimentation), respectively, for the early, middle, late and no-fire treatments. L. simplex had a single individual in the no-fire plots during the whole period, so the demography of this species was not studied in this treatment.

3.2. Asymptotic Growth Rate (Λ) of Grass Species under the Fire Treatments

The growth rates of all species were not significantly different from 1.0 under the early fire, as well as for A. canaliculatus regardless of the three other fire treatments (Table 1). The λs of A. schirensis and H. diplandra were not significantly different from 1.0 under the early and no-fire treatments. The following three species, A. schirensis, H. diplandra and L. simplex, exhibited low growth rates (λ < 1.0) under the late fire, as well as L. simplex under the middle fire (Table 1). Overall, the late fire caused the decline of three species (A. schirensis, H. diplandra and L. simplex). Projecting the population size over time (Figure 1) showed a reduction in the total population of 13%, 32% and 67%, respectively, in 10 years, and a trend towards extinction after 10 years. L. simplex also declined under the middle fire by 67% of the total population in 10 years. The population of this species grew exponentially over time under the early fire (Figure 1).

3.3. Loglinear Analysis

There was a significant effect of the treatment × species interaction and main effects of both treatment and species on the fate of the grass individuals (log-linear test p < 0.001, Table 2). This effect of the interaction was a clear increase in the death rate and decrease in the probability of transition in the two largest size-classes for A. schirensis, H. diplandra and L. simplex under the late fire treatment (Table 3). The effect of the middle fire was a reduction in the mortality of A. canaliculatus, A. schirensis and H. diplandra. The effect of the middle fire on L. simplex was an increase in mortality and a decrease of the probability of transition in class 5. The early fire exhibited a noticeable low mortality for H. diplandra and L. simplex and low growth of L. simplex toward the largest size-classes. The main effect of the no-fire treatment was a reduction of the transition probability in the smallest size-class for A. canaliculatus and H. diplandra (Table 3).

3.4. Elasticity Analysis

The proportional sensitivity of λ to the variations in the different categories of demographic parameters followed the same pattern for all species and fire treatments (fecundity < retrogression < growth < stasis) except for A. schirensis and H. diplandra that exhibited higher elasticity to growth than to stasis under the late fire (Table 4): the λs of A. schirensis and H. diplandra showed the highest elasticity to variations in growth under the late fire, whereas in other combinations of treatment and species, this highest elasticity of λs was seen in variations in stasis. Considering the size-classes, the proportional sensitivity of λs to variations in demographic parameters was highest for stasis whatever the species and the fire treatments (Figure 2). The highest elasticity of λ of H. diplandra to variations in stasis was observed in class 3 under all treatments, unlike the other species which exhibited highest elasticity of λ in different classes according to the fire treatments. This highest elasticity of λ to variations of stasis were found in the class 3 for A. canaliculatus (0.166) and A. schirensis (0.131) under the no-fire treatment and for A. canaliculatus in the classes 5, 2 and 3, respectively, in the early fire (0.134), the middle fire (0.260) and the late fire (0.127). A. schirensis showed highest elasticity of λ to variations of stasis in the classes 2 and 3, respectively, in the early fire (0.155) and the middle fire (0.259), while the highest elasticity of λ to variations of stasis was obtained in class 2 for L. simplex in the early and late fires (0.21 and 0.193, respectively) and in class 1 under the middle fire (0.231).

3.5. Stable Size Distribution Analysis

The distance between the observed and predicted size-class distributions given by the Keyfitz metric ranged from 0.025 (for A. canaliculatus under the late fire) to 0.23 (for L. simplex under the middle fire) showed that the two size-class distributions were not much different for all species under each fire treatment (Table S1). The predicted size-class distribution changed with the fire treatment (Figure 3). The predicted distribution exhibited the same shape for all species under the late fire treatment, with the highest proportions of individuals in class 2. The highest proportions of individuals in the predicted distribution were found for H. diplandra in class 3 under the other treatments. In the other cases, there were no clear patterns describing the predicted size-class distribution: the highest proportions were obtained in class 2 for A. canaliculatus, A. schirensis and L. simplex for the early fire. These proportions were found in classes 2, 3 and 1 for the middle fire for the same species, respectively. For the no-fire treatment, A. canaliculatus, and A. schirensis exhibited the highest proportions of individuals in classes 3 and 2, respectively. The lowest proportions (sometimes zero) were found in class 5 in all cases (Figure 3).

3.6. LTRE Analysis

Table 5 shows the difference of λ between the fire treatments, the overall contribution of each treatment to this difference and the influence of each demographic parameter on this difference. The population of A. canaliculatus grew faster under the early fire compared to the other treatments. In this species, growth was the parameter that contributed the most to the difference between this treatment and the middle and late fires, with the highest growth for the early fire. The difference between the early fire and no fire was marked by a higher retrogression for the early fire than the no-fire treatment. The λ of A. canaliculatus was higher for the middle fire than for the late fire due to higher stasis for the middle fire. Higher growth and lower fecundity for the no-fire treatment than for the middle and late fires, respectively, explained the difference of demography between no fire and these treatments.
In A. schirensis, H. diplandra and L. simplex, λ values were ordered as follows for the fire treatments: early fire > middle fire > no fire > late fire. The difference of demography between the early fire and the middle fire was mostly due to a higher retrogression for the early fire than the middle fire for A. schirensis and a higher growth for the early fire than the middle fire for H. diplandra and L. simplex. The demographic differences between the middle fire and the late fire, the middle fire and no fire, and no fire and the late fire were mostly due to stasis. This stasis was higher under the middle fire than no fire, then the middle fire than for the late fire, and no fire than the late fire.

3.7. Age-Based Life-History Parameters

The average age of individuals was between 12.2 y and 72.2 y; they resided from 3 y to 34 y in each size-class and their total conditional life-span was between 6 y and 46 y. Overall, these species can live between 24 y and 120 y. The mean age of A. canaliculatus in the different size-classes was high under all fire treatments (Table S2) with values ordered as follows: early fire (57.3 y) > late fire and no-fire treatments (45.7 y for both) > middle fire (38.9 y). In the class 1, this species was older under the late fire (36.1 y) than the middle fire (30.7 y) and the no-fire treatment (28.3 y). A. canaliculatus resided longer in the size-classes for the early fire (25.4 y) than for the middle fire (18.6 y), the late fire (13.7 y) and the no-fire treatments (9.6 y). All other age-based life-history parameters of this species followed the same trend except for mean time to first reach the other classes from class 1, which did not change much between the early fire (10.2 y) and the middle fire (10.4 y). This parameter showed lower values in classes 1, 2 and 3 and higher values in classes 4 and 5 under the middle fire than the early fire.
Andropogon schirensi had older individuals for the middle fire in all size-classes (72.2 y) than for the other treatments except the late fire, under which the age in the first class was the highest (Table S2). This species can live more than 100 y under the middle fire condition, since its mean remaining life-span was about 34.5 y and it resided longer in each size-class (33.9 y). The other age-based life-history parameters exhibited the following trend: middle fire > early fire > no fire > late fire. The mean age in each size-class of A. schirensis under the late fire (59.9 y) did not reflect the values of the other parameters, as the species spent less time in each size-class (3.7 y), quickly reached the other size-classes from class 1 (≈3.2 y) and had a low remaining life-span (≈2.9 y and 6.1 y).
All age-based life-history parameters of H. diplandra exhibited higher values for the early fire than for the middle fire, followed by no fire and the late fire (Table S2). This species was about two or three times older in the early fire (61.7 y) than in the other treatments, which did not show clear differences in the mean age (29.3 y, 21.9 y and 17.9 y), respectively, for the middle fire, the no fire and the late fire treatment. However, individuals of this species in the class 1 were older for the late fire (13.4 y) than for the no-fire treatment (12.8 y).
The mean age of L. simplex was higher in all size-classes for the middle fire (37.6 y) than for the early fire (19.2 y) and the late fire (16.6 y). The other parameters did not match the mean age values in the size-classes since they exhibited too-low values compared to the mean age (Table S2). The calculations gave unrealistic values for the mean age of residence in the size-classes, the mean time to first reach the classes from the class 1, the conditional remaining life-span in the classes, and the total conditional life-span under the early fire treatment for this species.
Moreover, there was no common pattern between the species regarding the variation of age-based life-history parameters according to the fire treatments. Similarly, the values of age-based life-history parameters of these species did not always follow the same order according to the fire treatments. For example, A. schirensis exhibited a higher value of mean age for the late fire treatment compared to the other species, whereas for the other parameters and this same fire treatment, it presented the lowest values. Finally, no general effect of a specific fire treatment was noticed on the age-based life-history parameters of these grass species.

4. Discussion

4.1. General Effect of the Fire Treatments on the Demography of the Grass Species

The late fire: The late fire caused a strong decline in the populations of three of the four species (A. schirensis, H. diplandra and L. simplex). This also leads to a decrease in the proportion of individuals in the two largest size-classes (stable size distribution). These results confirm the studies by Koffi et al. [24], who suggested that the late fire had high negative impacts on the demography of the dominant grasses. These negative effects can be explained by the fact that the late fire, with its high intensity [27], increases mortality, fragmentation and retrogression [24]. In turn, fragmentation and retrogression decrease the size of individuals, which likely make them more vulnerable to fire [24]. This general pattern is likely due to the late fire occurring at the beginning of the rainy season, when grasses have started to grow, so that their terminal buds become vulnerable to fire. This may lead to the death of parts of the tussocks (therefore increasing fragmentation and retrogression) or the death of whole tussocks (therefore increasing mortality). Indeed, dry-season fires have been shown to lead to lower risk of population decline than rainy-season fires [46], considering that early and middle fires are dry-season fires and the late fire is a rainy season fire. During the period of the late fire, grass phenological state could be crucial because fire could kill parts of the tillering of larger tussocks that have already been resprouted, as explained by Koffi et al. [24], which would also promote fragmentation and retrogression and increase death rate. The tussocks of L. simplex have a naturally small size and are therefore easily killed by the late fire. Moreover, on our late fire plots, Imperata cylindrica, one of the 100 most harmful invasive species in the world according to Lowe et al. [47], has become largely dominant over the years. Perhaps competition of this species for nutrients is partly responsible of the decline of A. schirensis, H. diplandra and L. simplex under this fire treatment.
The middle fire: The middle fire induced a decline of L. simplex by increasing death rate and reducing the probability of transition to class 5. These effects could be explained by the low capacity of this species to compete with other species for light [48]. Indeed, this species is characterized by frail and rather short tillers that are dominated by the other species at the end of their vegetation cycle, which could increase the death of individuals. This is confirmed by Raventos and Silva [49], who showed that competition with neighbours decreases the tillering capacity of perennial grasses. Thus, the synergy between the small size of L. simplex tussocks and the fairly high intensity of the middle fire could increase the death of L. simplex individuals [27]. The other species were not significantly declining under the middle fire, confirming their adaptation to this fire regime that has always been applied since at least 50 years in the area. The low mortality rate of these species under the middle fire revealed by the loglinear analysis confirmed this result. This agrees with Hiers et al. [50], who consider that the middle fire is the best management strategy for savanna species because it coincides with the historical fire regimes under which species have always evolved.
The no-fire treatment: This fire treatment did not induce the growth or decline of A. schirensis, H. diplandra and A. canaliculatus (the λ values were not significantly different from 1.0). However, we had no demographic data for L. simplex under the no-fire treatment because only one individual was recorded. In fact, this treatment leads to the accumulation of dead leaves and tillers [51], including from species that have higher tillers than L. simplex. This likely strongly reduces the species’ capacity to have access to light and would increase strongly its mortality. This agrees with studies by Kaye et al. [52] and Silva et al. [53], who reported a similar decline in populations of perennial grasses from other areas protected from fire. Contrary to L. simplex, the other species did not decline under the no-fire treatment, certainly because of their generally large tussock size that allows them to better resist shading by the accumulation of dead aerial biomass. Conversely, the absence of fire would have favorable effects (in comparison to all fire treatments) by suppressing the increase in mortality, retrogression and fragmentation due to fire. This would contribute, as shown by the observed size-class distribution, to a low proportion of individuals in the small size-classes.
The early fire: The early fire did not induce the growth or decline of any species. This treatment induces a low fire intensity due to the abundant but very humid fuel [27], which decreases fire capacity to increase mortality, fragmentation and retrogression [24]. Particularly in L. simplex, seed dispersal takes place early (a few weeks before the early fire), so the seeds are disseminated before this fire treatment at the base of the tussocks [54] and are more likely to settle deeper in the litter or soil before the fire is set [46]. Since early fire often has a weak impact on vegetation, these seeds can be spared. After fire, nutrient availability, microbial activity and nitrogen fixation are increased [55,56]. This can induce abundant seed germination and rapid seedling growth [57] and increase population growth rate [58]. Thus, seeds of L. simplex, spared from the early fire, can germinate, especially when soil is fairly wet (as it was the case at the time of this fire treatment), and then take advantage of resource availability after fire to grow and form new tussocks. Higher growth, resulting in higher fecundity and recruitment in annually burnt populations [53], could explain the high growth rate of 1.061 observed in this species under the early fire.

4.2. Responses of the Grass Species to Fire Treatments

4.2.1. Similarities in the Response of Species to the Fire Treatments

Elasticity, stable size distribution, LTRE and age-based life-history parameter analyzes showed that perennial grass species had similarities in their responses to fire treatments. Indeed, all species were more sensitive to variations in stasis under all fire treatments. This high proportional sensitivity of λ to variations in stasis agrees with the fact that, for most long-lived plants in arid and semi-arid environments, stasis is the single most important demographic process [59,60]. In our case, all species we studied were long-lived plants (up to 120 y) but in a humid savanna [29]. Thus, the high elasticity of λ to stasis can also be a characteristic of long-lived plants in humid environments. Though the four species are adapted to their fire-dominated environment [61], increasing the survival of individuals would be the best way to increase the growth of their populations.. Another common point between the four species is that they were all close to the stable size-class distribution predicted by the matrix models under all fire treatments. Similar results were obtained on a perennial grass species: Andropogon semiberbis in Venezuela [53], and on a range of herb, shrub, and tree species [62]. One could predict, as Clutton-Brock and Coulson [63] and Ujvari et al. [64] did, that environmental disturbances such as fire lead to a difference between the predicted and the observed size-class distributions. This suggests that the four studied species adjust quickly their size distributions to disturbances. In any case, this proximity between the two distributions suggests, according to Williams et al. [62], that relying on the equilibrium assumption to make demographic analyzes of populations was valid.

4.2.2. Contrasting Responses of the Grass Species to the Fire Treatments

The demographic parameters of each species were impacted differently by the fire treatments, showing dissimilarity between the demography of the species. The fact that population growth was higher for the early and the middle fires than for the late fire for A. canaliculatus agrees with Liu and Menges [65], who found that recruitment was stronger with dry season fires for Chamaecrista keyensis, our late fire sometimes coinciding with the rainy season. Changes in the shape of the predicted size-class distribution of each species according to the fire treatment confirmed that the fire treatments acted slightly differently on the size distribution of each species [24]. The highest proportion of A. canaliculatus individuals was predicted in class 2 for all fire regimes, but in class 3 for the no-fire treatment, indicating that the demography of this species is maintained when submitted to fire but changes in absence of fire. The competitive ability of this species for resources would certainly increase for the no-fire treatment compared to the other treatments because individual size (the highest proportion of individuals in class 3 under the no-fire treatment) and competitive ability increases with plant size [66,67]. This could explain the difference of demography of A. canaliculatus between fire and no-fire treatments and confirms the fact that A. canaliculatus is a very fire-adapted species. The response of H. diplandra to the fire treatments was a strong reduction of the size of individuals predicted for the late fire. The same result was obtained for A. schirensis for the early fire, the late fire and the no-fire treatments compared to the middle fire. Indeed, the middle fire is the fire regime that is usually applied in the Lamto savanna and to which all species are already adapted. According to Miller et al. [68], changes from the fire treatment to which species are adapted tend to generate a mismatch between different fire-adaptation traits, which negatively affect population dynamics. For example, the highest proportion of individuals of A. schirensis was predicted in class 2 for early and late fire régimes and in class 1 for the absence of fire, whereas for the middle fire, the highest proportion of individuals was predicted in class 3. For the same type of reason, L. simplex, contrary to A. schirensis, exhibited an increase of individual size when submitted to other fire regimes than the middle fire: the predicted proportion of individuals was the highest in size-class 2 instead of size-class 1 for the middle fire.
The life-history parameters of all species differed by fire treatment. Some species had older individuals in the early fire and others in the middle fire treatment, while some species lived longer in the late fire and others in the no-fire treatment. The other age-based parameters also were characteristic. This difference in species response to fire treatments was more marked between size classes as shown by Schemske et al. [69] and Kesler et al. [70]. This probably reflects both the discussed differences in their demographic responses and the complexity in plants of the links between size and age. This complexity is probably increased here by the regression and fragmentation of tussocks.

5. Conclusions

Our results demonstrated that only the A. canaliculatus population declined, but not significantly (λ not significantly lower than 1.0), under all fire treatments, while the late fire caused decline in all other species. This suggests that applying the late fires alone should be avoided to maintain grass diversity. The early and middle fire treatments could be recommended to savanna managers, as they are less detrimental to almost all grass species and have less effect on the stasis, which is the most sensitive parameter in the demography of grasses. However, this implies that these fire treatments must also be efficient in reducing or maintaining tree density. If not, it would be appropriate to test other methods such as the combination of these fire treatments.
In fact, applying a mixture of middle, early, late and no-fire treatments in space or time could favor grass diversity and differences in grass size structure of their tussocks, as suggested by differences in their responses to these fire regimes. For example, in the case of the Lamto savanna, the late fire could be applied, for instance, every 5 years, instead of the usual annual middle fire, to regulate the populations of woody species as suggested by N’Dri et al. [22], probably without being detrimental to perennial grass biodiversity. The no-fire treatment alone cannot normally be considered as a treatment that helps to preserve savanna biodiversity because it leads to a rapid declining of some grass species and the transformation of savanna into forest [71]. However, a fire regime based on the number of years between two fires could be applied as studied in other savanna types [72]. This should not be detrimental to perennial grass diversity, even to L. simplex, which tends to disappear without fire.
We have so far studied the demography of the four dominant grass species as independent phenomena, while facilitation (which could occur through the local enhancement of soil fertility [73]) and competition (e.g., for light) could likely lead to interactions between these species. It would be thus interesting to study the competition/facilitation relationships between the perennial grass species using an experimental approach, by artificially reducing the density of individuals of one or more species and monitoring the demographic effect of this reduction on their demography. This could allow parameterising a matrix population model in which demographic rates would depend on tussock densities [34].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/fire5060193/s1.

Author Contributions

Conceptualization, K.F.K., A.B.N., J.-C.L. and S.B.; methodology, A.B.N., J.-C.L. and S.B.; formal analysis, K.F.K. and S.B.; investigation, K.F.K., J.-C.L., M.K. and S.B.; resources, A.B.N., T.S., J.-C.L., S.K., M.K. and S.B.; data curation, K.F.K., M.K. and S.B.; writing—original draft preparation, K.F.K. and S.B.; writing—review and editing, K.F.K., A.B.N., T.S., J.-C.L., S.K. and S.B.; supervision, A.B.N., J.-C.L., S.K. and S.B.; project administration, A.B.N., S.K. and S.B.; funding acquisition, K.F.K., A.B.N., T.S., J.-C.L. and S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the PhD grant provided by the Direction of Research and Training in the South of the Institute of Research for Development (IRD-DPF) to the first author from 2016 to 2019. The study was conducted within the project “Integrated and Sustainable Management of Savanna Ecosystems’’ funded by the IRD within the “Young Associated Team” program. This project funded the management of the fire regimes on the experimental plots.

Data Availability Statement

Our raw data set will be made available to the scientific community on the data repository of the Institute of Research for Development (IRD) (https://dataverse.ird.fr, accessed on 15 November 2022) as soon as the paper is accepted for publication.

Acknowledgments

We thank Yéo Kolo, the Director of the Lamto station, for allowing us to work at the station of ecology of Lamto. This research work was funded by the GainGrass Project (Global Assessment of Nitrification Inhibition by tropical Grasses Project, https://anr.fr/Project-ANR-19-CE02-0009 accessed on 10 September 2022) funded by the French National Research Agency (ANR).

Conflicts of Interest

The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study.

References

  1. Sankaran, M.; Hanan, N.P.; Scholes, R.J.; Ratnam, J.; Augustine, D.J.; Cade, B.S.; Gignoux, J.; Higgins, S.I.; Le Roux, X.; Ludwig, F.; et al. Determinants of Woody Cover in African Savannas. Nature 2005, 438, 846–849. [Google Scholar] [CrossRef] [PubMed]
  2. Menaut, J.-C.; Abbadie, L. Vegetation. In Lamto; Abbadie, L., Gignoux, J., Le Roux, X., Lepage, M., Eds.; Ecological Studies; Springer: New York, NY, USA, 2006; Volume 179, pp. 63–74. [Google Scholar]
  3. Williams, R.J.; Müller, W.J.; Wahren, C.-H.; Setterfield, S.A.; Cusack, J. Vegetation. In Fire in Tropical Savannas; Andersen, A.N., Cook, G.D., Williams, R.J., Eds.; Ecological Studies; Springer: New York, NY, USA, 2003; pp. 79–106. [Google Scholar]
  4. Savadogo, P.; Zida, D.; Sawadogo, L.; Tiveau, D.; Tigabu, M.; Oden, P.C. Fuel and Fire Characteristics in Savanna-Woodland of West Africa in Relation to Grazing and Dominant Grass Type. Int. J. Wildland Fire 2007, 16, 531–539. [Google Scholar] [CrossRef]
  5. Gignoux, J.; Lahoreau, G.; Julliard, R.; Barot, S. Establishment and Early Persistence of Tree Seedlings in an Annually Burned Savanna. J. Ecol. 2009, 97, 484–495. [Google Scholar] [CrossRef]
  6. Smit, I.P.J.; Asner, G.P.; Govender, N.; Vaughn, N.R.; van Wilgen, B.W. An Examination of the Potential Efficacy of High-Intensity Fires for Reversing Woody Encroachment in Savannas. J. Appl. Ecol. 2016, 53, 1623–1633. [Google Scholar] [CrossRef]
  7. Stevens, N.; Lehmann, C.E.R.; Murphy, B.P.; Durigan, G. Savanna Woody Encroachment Is Widespread across Three Continents. Glob. Chang. Biol. 2017, 23, 235–244. [Google Scholar] [CrossRef] [Green Version]
  8. Barot, S.; Gignoux, J.; Legendre, S.; Vuattoux, R. Demography of a Savanna Palm Tree in Ivory Coast (Lamto): Population Persistence, and Life History. J. Trop. Ecol. 2000, 16, 637–655. [Google Scholar] [CrossRef]
  9. Jurisch, K.; Hahn, K.; Wittig, R.; Bernhardt-Römermann, M. Population Structure of Woody Plants in Relation to Land Use in a Semi-Arid Savanna, West Africa. Biotropica 2012, 44, 744–751. [Google Scholar] [CrossRef]
  10. O’Connor, T.G. Demography of Woody Species in a Semi-Arid African Savanna Reserve Following the Re-Introduction of Elephants. Acta Oecologica 2017, 78, 61–70. [Google Scholar] [CrossRef]
  11. Garnier, L.K.M.; Dajoz, I. The Influence of Fire on the Demography of a Dominant Grass Species of West African Savannas, Hyparrhenia diplandra. J. Ecol. 2001, 89, 200–208. [Google Scholar] [CrossRef]
  12. Keeley, J.E.; Pausas, J.G.; Rundel, P.W.; Bond, W.J.; Bradstock, R.A. Fire as an Evolutionary Pressure Shaping Plant Traits. Trends Plant Sci. 2011, 16, 406–411. [Google Scholar] [CrossRef]
  13. Pilon, N.A.L.; Hoffmann, W.A.; Abreu, R.C.R.; Durigan, G. Quantifying the Short-Term Flowering after Fire in Some Plant Communities of a Cerrado Grassland. Plant Ecol. Divers. 2018, 11, 259–266. [Google Scholar] [CrossRef]
  14. Gignoux, J.; Clobert, J.; Menaut, J.-C. Alternative Fire Resistance Strategies in Savanna Trees. Oecologia 1997, 110, 576–583. [Google Scholar] [CrossRef] [PubMed]
  15. Bond, W.J. What Limits Trees in C4 Grasslands and Savannas? Annu. Rev. Ecol. Evol. Syst. 2008, 39, 641–659. [Google Scholar] [CrossRef] [Green Version]
  16. Accatino, F.; Wiegand, K.; Ward, D.; De Michele, C. Trees, Grass, and Fire in Humid Savannas: The Importance of Life History Traits and Spatial Processes. Ecol. Model. 2016, 320, 135–144. [Google Scholar] [CrossRef]
  17. Ott, J.P.; Klimešová, J.; Hartnett, D.C. The Ecology and Significance of Below-Ground Bud Banks in Plants. Ann. Bot. 2019, 123, 1099–1118. [Google Scholar] [CrossRef]
  18. Bond, W.J.; Wilgen, B.W. Fire and Plants; Chapman & Hall: London, UK; Springer: London, UK, 1996; pp. 1–263. [Google Scholar]
  19. Simpson, K.J.; Ripley, B.S.; Christin, P.-A.; Belcher, C.M.; Lehmann, C.E.R.; Thomas, G.H.; Osborne, C.P. Determinants of Flammability in Savanna Grass Species. J. Ecol. 2016, 104, 138–148. [Google Scholar] [CrossRef] [Green Version]
  20. N’Dri, A.B.; Fongbe, M.; Soro, T.D.; Gignoux, J.; Kone, M.; Dosso, K.; N’dri, J.K.; Kone, N.A.; Barot, S. Principaux Indices de l’intensité Du Feu Dans Une Savane Guinéenne d’Afrique de l’ouest. Int. J. Biol. Chem. Sci. 2018, 12, 266–274. [Google Scholar] [CrossRef]
  21. Werner, P.A.; Prior, L.D. Demography and Growth of Subadult Savanna Trees: Interactions of Life History, Size, Fire Season, and Grassy Understory. Ecol. Monogr. 2013, 83, 67–93. [Google Scholar] [CrossRef]
  22. N’Dri, A.B.; Kpangba, K.P.; Werner, P.A.; Koffi, K.F.; Bakayoko, A. The Response of Sub-Adult Savanna Trees to Six Successive Annual Fires: An Experimental Field Study on the Role of Fire Season. J. Appl. Ecol. 2022, 59, 1347–1361. [Google Scholar] [CrossRef]
  23. Flory, S.L.; Bauer, J.; Phillips, R.P.; Clay, K. Effects of a Non-Native Grass Invasion Decline over Time. J. Ecol. 2017, 105, 1475–1484. [Google Scholar] [CrossRef]
  24. Koffi, K.F.; N’Dri, A.B.; Lata, J.-C.; Konaté, S.; Srikanthasamy, T.; Konaré, S.; Konan, M.; Barot, S. Effect of Fire Regimes on the Demographic Parameters of the Perennial Tussock Grasses of a Humid Savanna. J. Veg. Sci. 2019, 30, 950–962. [Google Scholar] [CrossRef]
  25. Amoako, E.; Misana, S.; Kranjac-Berisavljevic, G.; Zizinga, A.; Duwieja, B.A. Effect of the Seasonal Burning on Tree Species in the Guinea Savanna Woodland, Ghana: Implications for Climate Change Mitigation. Appl. Ecol. Environ. Res. 2018, 16, 1935–1949. [Google Scholar] [CrossRef]
  26. Oliveira, S.L.J.; Maier, S.W.; Pereira, J.M.C.; Russell-Smith, J. Seasonal Differences in Fire Activity and Intensity in Tropical Savannas of Northern Australia Using Satellite Measurements of Fire Radiative Power. Int. J. Wildland Fire 2015, 24, 249. [Google Scholar] [CrossRef]
  27. N’Dri, A.B.; Kpré, A.J.-N.; Kpangba, K.P.; Soro, T.D.; Kouassi, K.V.; Koffi, K.F.; Kouamé, Y.A.G.; Koffi, A.B.; Konan, L.N. Experimental Study of Fire Behavior in Annually Burned Humid Savanna of West Africa in the Context of Bush Encroachment. In Sustainable Development in Africa; Leal Filho, W., Pretorius, R., de Sousa, L.O., Eds.; World Sustainability Series; Springer International Publishing: Berlin/Heidelberg, Germany, 2021; pp. 491–505. [Google Scholar]
  28. Moleele, N.M.; Ringrose, S.; Matheson, W.; Vanderpost, C. More Woody Plants? The Status of Bush Encroachment in Botswana’s Grazing Areas. J. Environ. Manag. 2002, 64, 3–11. [Google Scholar] [CrossRef]
  29. Koffi, K.F.; N’Dri, A.B.; Konaré, S.; Srikanthasamy, T.; Lata, J.-C.; Konaté, S.; Konan, M.; Barot, S. Demography of the Dominant Perennial Grass Species of a Humid African Savanna. Acta Oecologica 2022, 114, 103816. [Google Scholar] [CrossRef]
  30. Koffi, K.F.; N’Dri, A.B.; Lata, J.-C.; Konaté, S.; Srikanthasamy, T.; Konan, M.; Barot, S. Effect of Fire Regime on the Grass Community of the Humid Savanna of Lamto, Ivory Coast. J. Trop. Ecol. 2019, 35, 1–7. [Google Scholar] [CrossRef]
  31. Lamotte, M.; Tireford, J.-L. Le Climat de la Savane de Lamto (Côte d’Ivoire) et sa Place Dans les Climats de l’Ouest Africain; Université d’Abidjan, Station d’Ecologie Tropicale de Lamto: N’douci, Côte d’Ivoire, 1988; pp. 1–146. [Google Scholar]
  32. Abbadie, L.; Gignoux, J.; Lepage, M.; Le Roux, X. Environmental Constraints on Living Organisms. In Lamto; Abbadie, L., Gignoux, J., Le Roux, X., Lepage, M., Eds.; Ecological Studies; Springer: New York, NY, USA, 2006; Volume 179, pp. 45–61. [Google Scholar]
  33. Harvell, C.D.; Caswell, H.; Simpson, P. Density Effects in a Colonial Monoculture: Experimental Studies with a Marine Bryozoan (Membranipora membranacea L.). Oecologia 1990, 82, 227–237. [Google Scholar] [CrossRef]
  34. Caswell, H. Matrix Population Models: Construction, Analysis, and Interpretation; Sinauer Associates Inc.: Sunderland, UK, 2001; pp. 1–770. [Google Scholar]
  35. Crouse, D.T.; Crowder, L.B.; Caswell, H. A Stage-Based Population Model for Loggerhead Sea Turtles and Implications for Conservation. Ecology 1987, 68, 1412–1423. [Google Scholar] [CrossRef]
  36. Huenneke, L.F.; Marks, P.L. Stem Dynamics of the Shrub Alnus incana ssp. rugosa: Transition Matrix Models. Ecology 1987, 68, 1234–1242. [Google Scholar]
  37. Caswell, H. Matrix Population Models; Sinauer Associates Inc.: Sunderland, UK, 1989; pp. 1–328. [Google Scholar]
  38. De Kroon, H.; Plaisier, A.; Van Groenendael, J.; Caswell, H. Elasticity: The Relative Contribution of Demographic Parameters to Population Growth Rate. Ecology 1986, 67, 1427–1431. [Google Scholar] [CrossRef]
  39. Van Groenendael, J.; De Kroon, H.; Caswell, H. Projection Matrices in Population Biology. Trends Ecol. Evol. 1988, 3, 264–269. [Google Scholar] [CrossRef]
  40. Fingleton, B. Models of Category Counts; Cambridge University Press: Cambridge, UK, 1984; pp. 1–204. [Google Scholar]
  41. Cochran, M.E.; Ellner, S. Simple Methods for Calculating Age-Based Life History Parameters for Stage-Structured Populations. Ecol. Monogr. 1992, 62, 345–364. [Google Scholar] [CrossRef]
  42. Barot, S.; Gignoux, J.; Legendre, S. Stage-Structured Matrix Models and Age Estimates. Oikos 2002, 96, 56–61. [Google Scholar] [CrossRef] [Green Version]
  43. Keyfitz, N. Introduction to the Mathematics of Populations; Addison-Wesley Educational Publishers Inc: Boston, MA, USA, 1968; pp. 1–450. [Google Scholar]
  44. Stubben, C.J.; Milligan, B.G. Estimating and Analyzing Demographic Models Using the Popbio Package in R. J. Stat. Softw. 2007, 22, 1–23. [Google Scholar] [CrossRef] [Green Version]
  45. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2021. [Google Scholar]
  46. Liu, H.; Menges, E.S.; Quintana-Ascencio, P.F. Population Viability Analyses of Chamaecrista keyensis: Effects of Fire Season and Frequency. Ecol. Appl. 2005, 15, 210–221. [Google Scholar] [CrossRef] [Green Version]
  47. Lowe, S.; Browne, M.; Boudjelas, M.; De Poorter, M. 100 of the World’s Worst Invasive Alien Species: A Selection from the Global Invasive Species Database; IUCN: Gland, Switzerland; Hollands Printing Ltd: Auckland, New Zealand, 2000; pp. 1–12. [Google Scholar]
  48. D’Antonio, C.M.; Vitousek, P.M. Biological Invasions by Exotic Grasses, the Grass/Fire Cycle, and Global Change. Annu. Rev. Ecol. Syst. 1992, 2, 63–87. [Google Scholar] [CrossRef]
  49. Raventós, J.; Silva, J.F. Competition Effects and Responses to Variable Numbers of Neighbours in Two Tropical Savanna Grasses in Venezuela. J. Trop. Ecol. 1995, 11, 39–52. [Google Scholar] [CrossRef]
  50. Hiers, J.K.; Wyatt, R.; Mitchell, R.J. The Effects of Fire Regime on Legume Reproduction in Longleaf Pine Savannas: Is a Season Selective? Oecologia 2000, 125, 521–530. [Google Scholar] [CrossRef]
  51. Zimmermann, J.; Higgins, S.I.; Grimm, V.; Hoffmann, J.; Linstädter, A. Grass Mortality in Semi-Arid Savanna: The Role of Fire, Competition and Self-Shading. Perspect. Plant Ecol. Evol. Syst. 2010, 12, 1–8. [Google Scholar] [CrossRef]
  52. Kaye, T.N.; Pendergrass, K.L.; Finley, K.; Kauffman, J.B. The Effect of Fire on the Population Viability of an Endangered Prairie Plant. Ecol. Appl. 2001, 11, 15. [Google Scholar] [CrossRef]
  53. Silva, J.F.; Raventos, J.; Caswell, H.; Trevisan, M.C. Population Responses to Fire in a Tropical Savanna Grass, Andropogon semiberbis: A Matrix Model Approach. J. Ecol. 1991, 79, 345–355. [Google Scholar] [CrossRef]
  54. Garnier, L.K.; Durand, J.; Dajoz, I. Limited Seed Dispersal and Microspatial Population Structure of an Agamospermous Grass of West African Savannahs, Hyparrhenia diplandra (Poaceae). Am. J. Bot. 2002, 89, 1785–1791. [Google Scholar] [CrossRef] [PubMed]
  55. Hulbert, L.C. Causes of Fire Effects in Tallgrass Prairie. Ecology 1988, 69, 46–58. [Google Scholar] [CrossRef]
  56. Srikanthasamy, T.; Barot, S.; Koffi, F.K.; Tambosco, K.; Marcangeli, Y.; Carmignac, D.; N’Dri, A.B.; Gervaix, J.; Roux, X.L.; Lata, J.-C. Short-Term Impact of Fire on the Total Soil Microbial and Nitrifier Communities in a Wet Savanna. Ecol. Evol. 2021, 11, 9958–9969. [Google Scholar] [CrossRef] [PubMed]
  57. Baruch, Z.; Bilbao, B. Effects of Fire and Defoliation on the Life History of Native and Invader C4 Grasses in a Neotropical Savanna. Oecologia 1999, 119, 510–520. [Google Scholar] [CrossRef]
  58. Lesica, P. Effects of Fire on the Demography of the Endangered, Geophytic Herb Silene spaldingii (Caryophyllaceae). Am. J. Bot. 1999, 86, 996–1002. [Google Scholar] [CrossRef]
  59. Jiménez-Sierra, C.; Mandujano, M.C.; Eguiarte, L.E. Are Populations of the Candy Barrel Cactus (Echinocactus platyacanthus) in the Desert of Tehuacán, Mexico at Risk? Population Projection Matrix and Life Table Response Analysis. Biol. Conserv. 2007, 135, 278–292. [Google Scholar] [CrossRef]
  60. Jiménez-Valdés, M.; Godínez-Alvarez, H.; Caballero, J.; Lira, R. Population Dynamics of Agave marmorata Roezl. under Two Contrasting Management Systems in Central Mexico. Econ. Bot. 2010, 64, 149–160. [Google Scholar] [CrossRef]
  61. Morris, W.F.; Tuljapurkar, S.; Haridas, C.V.; Menges, E.S.; Horvitz, C.C.; Pfister, C.A. Sensitivity of the Population Growth Rate to Demographic Variability within and between Phases of the Disturbance Cycle. Ecol. Lett. 2006, 9, 1331–1341. [Google Scholar] [CrossRef]
  62. Williams, J.L.; Ellis, M.M.; Bricker, M.C.; Brodie, J.F.; Parsons, E.W. Distance to Stable Stage Distribution in Plant Populations and Implications for Near-Term Population Projections: Distance to Stable Stage Distribution. J. Ecol. 2011, 99, 1171–1178. [Google Scholar] [CrossRef]
  63. Clutton–Brock, T.H.; Coulson, T. Comparative Ungulate Dynamics: The Devil Is in the Detail. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 2002, 357, 1285–1298. [Google Scholar] [CrossRef] [PubMed]
  64. Ujvari, B.; Andersson, S.; Brown, G.; Shine, R.; Madsen, T. Climate-Driven Impacts of Prey Abundance on the Population Structure of a Tropical Aquatic Predator. Oikos 2010, 119, 188–196. [Google Scholar] [CrossRef]
  65. Liu, H.; Menges, E.S. Winter Fires Promote Greater Vital Rates in the Florida Keys than Summer Fires. Ecology 2005, 86, 1483–1495. [Google Scholar] [CrossRef]
  66. Violle, C.; Richarte, J.; Navas, M.-L. Effects of Litter and Standing Biomass on Growth and Reproduction of Two Annual Species in a Mediterranean Old-Field. J. Ecol. 2006, 94, 196–205. [Google Scholar] [CrossRef]
  67. Wang, P.; Stieglitz, T.; Zhou, D.W.; Cahill, J.F., Jr. Are Competitive Effect and Response Two Sides of the Same Coin, or Fundamentally Different? Funct. Ecol. 2010, 24, 196–207. [Google Scholar] [CrossRef]
  68. Miller, R.G.; Tangney, R.; Enright, N.J.; Fontaine, J.B.; Merritt, D.J.; Ooi, M.K.J.; Ruthrof, K.X.; Miller, B.P. Mechanisms of Fire Seasonality Effects on Plant Populations. Trends Ecol. Evol. 2019, 34, 1104–1117. [Google Scholar] [CrossRef] [PubMed]
  69. Schemske, D.W.; Husband, B.C.; Ruckelshaus, M.H.; Goodwillie, C.; Parker, I.M.; Bishop, J.G. Evaluating Approaches to the Conservation of Rare and Endangered Plants. Ecology 1994, 75, 584–606. [Google Scholar] [CrossRef] [Green Version]
  70. Kesler, H.C.; Trusty, J.L.; Hermann, S.M.; Guyer, C. Demographic Responses of Pinguicula ionantha to Prescribed Fire: A Regression-Design LTRE Approach. Oecologia 2008, 156, 545–557. [Google Scholar] [CrossRef]
  71. Douffi, K.G.-C.; Yao, A.C.; Koffi, K.J.; Traore, A.S.; Kone, M. Afforestation in Response to Thermal Change in the Forest-Savannah Transition of the Lamto Scientific Reserve, Côte d’Ivoire. Eur. J. For. Eng. 2021, 7, 45–56. [Google Scholar]
  72. Van Wilgen, B.W.; Govender, N.; Biggs, H.C.; Ntsala, D.; Funda, X.N. Response of Savanna Fire Regimes to Changing Fire Management Policies in a Large African National Park. Conserv. Biol. 2004, 18, 1533–1540. [Google Scholar] [CrossRef]
  73. Srikanthasamy, T.; Leloup, J.; N’Dri, A.B.; Barot, S.; Gervaix, J.; Koné, A.W.; Koffi, K.F.; Le Roux, X.; Raynaud, X.; Lata, J.-C. Contrasting Effects of Grasses and Trees on Microbial N-Cycling in an African Humid Savanna. Soil Biol. Biochem. 2018, 117, 153–163. [Google Scholar] [CrossRef]
Figure 1. Total population size projected over 50 years by species and fire treatment using the average transition matrices.
Figure 1. Total population size projected over 50 years by species and fire treatment using the average transition matrices.
Fire 05 00193 g001
Figure 2. Elasticity matrices of all species in the four fire treatments. The redder the color, the higher the value.
Figure 2. Elasticity matrices of all species in the four fire treatments. The redder the color, the higher the value.
Fire 05 00193 g002
Figure 3. Comparison of observed and predicted distributions of individuals in size-classes for each combination of species (Andropogon canaliculatus, Andropogon schirensis, Hyparrhenia diplandra and Loudetia simplex) and fire treatments. The predicted size distributions were calculated using the mean matrices over the four one-year transitions.
Figure 3. Comparison of observed and predicted distributions of individuals in size-classes for each combination of species (Andropogon canaliculatus, Andropogon schirensis, Hyparrhenia diplandra and Loudetia simplex) and fire treatments. The predicted size distributions were calculated using the mean matrices over the four one-year transitions.
Fire 05 00193 g003
Table 1. Asymptotic grow rates (λs), the associated standard errors (SE) and 95% confidence intervals calculated for the mean transition matrix of the four species in each fire treatment. The (−) denote that there was no value for L. simplex in the no-fire treatment because there was almost no individual of this species in this treatment.
Table 1. Asymptotic grow rates (λs), the associated standard errors (SE) and 95% confidence intervals calculated for the mean transition matrix of the four species in each fire treatment. The (−) denote that there was no value for L. simplex in the no-fire treatment because there was almost no individual of this species in this treatment.
SpeciesFire TreatmentsλsSE95% Confidence Interval
Andropogon canaliculatusearly fire0.9780.040[0.899, 1.057]
middle fire0.9690.031[0.907, 1.030]
late fire0.9450.047[0.853, 1.037]
no fire0.9150.055[0.807, 1.024]
Andropogon schirensisearly fire0.9890.099[0.794, 1.185]
middle fire0.9850.047[0.893, 1.078]
late fire0.6700.103[0.467, 0.872]
no fire0.9200.123[0.677, 1.162]
Hyparrhenia diplandraearly fire0.9510.043[0.865, 1.036]
middle fire0.9530.036[0.881, 1.025]
late fire0.7830.065[0.655, 0.912]
no fire0.9330.054[0.826, 1.040]
Loudetia simplexearly fire1.0610.060[0.943, 1.179]
middle fire0.8720.042[0.788, 0.955]
late fire0.7580.064[0.632, 0.885]
no fire---
Table 2. Loglinear analysis of the transition frequency table showing the effects of the treatments (T), the species (E) and the interaction on the demographic fate (D) of the four perennial grass species, conditional on the initial state (S). The DS,SET is the null model. This model contains the main effects of species and treatment and the interaction. It also excludes all interaction between fate and any combination of species and treatment. Comparing the null model to the model obtained after adding the species or treatment or the interaction in the null model give the corresponding effect. The log-likelihood ratio (G2) calculated after adding 0.5 to each table cell, as suggested by Fingleton [40], measures the goodness-of-fit of the loglinear model. The variation of G2 (ΔG2) assesses the significance of the observed effects. The degree of freedom (df) is the difference between the number of cells in the table and the number of parameters in the model.
Table 2. Loglinear analysis of the transition frequency table showing the effects of the treatments (T), the species (E) and the interaction on the demographic fate (D) of the four perennial grass species, conditional on the initial state (S). The DS,SET is the null model. This model contains the main effects of species and treatment and the interaction. It also excludes all interaction between fate and any combination of species and treatment. Comparing the null model to the model obtained after adding the species or treatment or the interaction in the null model give the corresponding effect. The log-likelihood ratio (G2) calculated after adding 0.5 to each table cell, as suggested by Fingleton [40], measures the goodness-of-fit of the loglinear model. The variation of G2 (ΔG2) assesses the significance of the observed effects. The degree of freedom (df) is the difference between the number of cells in the table and the number of parameters in the model.
ModelEffectsG2ΔG2dfΔdfp-Value
DS,SET 1573.7 375 <0.001
DSE,SET 1176.1 300 <0.001
DE,DSESpecies 397.6 75<0.001
DST,SET 900.9 300 <0.001
DT,DSTTreatment 672.8 75<0.001
DSE,DST,SET 515.0 225 <0.001
DT,DSTTreatment 661.1 75<0.001
DE,DSESpecies 385.9 75<0.001
DSET 0.0 0.0 1.000
DET,DSETSpecies × Treatment 515.0 225<0.001
Table 3. Effect of the treatment × species interaction on the fate of grass individuals shown by the loglinear analysis. Table entries give the proportions (%) of transitions from individual initial states to six possible fates (dead, classes 1–5). The total number of individual in each fate by species and treatment was used in the calculations.
Table 3. Effect of the treatment × species interaction on the fate of grass individuals shown by the loglinear analysis. Table entries give the proportions (%) of transitions from individual initial states to six possible fates (dead, classes 1–5). The total number of individual in each fate by species and treatment was used in the calculations.
SpeciesTreatments Fate
DeadClass 1Class 2Class 3Class 4Class 5
Andropogon canaliculatusearly fire 11.211.424.723.614.614.6
middle fire7.015.229.525.111.012.2
late fire14.716.430.821.910.16.1
no fire16.39.120.825.816.911.2
Andropogon schirensisearly fire 11.99.933.824.814.65.0
middle fire5.37.223.331.419.313.5
late fire36.810.132.014.54.81.8
no fire17.312.819.620.712.317.3
Hyparrhenia diplandraearly fire 8.67.225.128.316.314.6
middle fire7.29.725.629.917.110.5
late fire33.712.829.018.74.71.1
no fire13.28.424.027.416.310.7
Loudetia simplexearly fire 5.114.040.933.45.41.3
middle fire10.324.330.222.58.14.6
late fire27.314.730.023.14.01.0
no fire------
Table 4. Elasticity as the contribution (%) of fecundity, stasis, retrogression and growth to the λs of all species under each fire treatment. The (-) denote that there was no value for L. simplex in the no-fire treatment.
Table 4. Elasticity as the contribution (%) of fecundity, stasis, retrogression and growth to the λs of all species under each fire treatment. The (-) denote that there was no value for L. simplex in the no-fire treatment.
SpeciesTreatmentsMatrix Entry Categories
FecundityStasisRetrogressionGrowth
Andropogon canaliculatusearly fire 7.546.419.026.8
middle fire9.455.912.022.5
late fire8.249.313.428.9
no fire4.451.317.426.7
Andropogon schirensisearly fire 6.950.215.627.1
middle fire3.159.316.421.0
late fire10.533.020.735.6
no fire12.142.512.632.5
Hyparrheniadiplandraearly fire 4.955.116.223.5
middle fire5.552.216.525.7
late fire15.034.913.036.9
no fire8.141.017.633.0
Loudetiasimplexearly fire 13.344.713.028.8
middle fire16.848.99.025.0
late fire16.239.911.332.3
no fire----
Table 5. Contribution of the grouped demographic parameters (Fec.: fecundity, Sta.: stasis, Ret.: retrogression and Gro.: growth) to the demographic difference between treatments in each species, analysed by the life-table response experiment. When the demographic parameters are higher in treatment 2 than in treatment 1, contribution values are positive and vice versa. Treatment 1: the first of the two compared fire treatments; treatment 2: the second treatment; ∆𝜆 = λ(treatment 2) − 𝜆(treatment 1): the difference between the asymptotic growth rates (λ) of treatment 1 and treatment 2; C-value: the global contribution value calculated by summing all the entries of the contribution matrix from the comparison between the corresponding treatment. The (-) denote that there was no possible comparison between the other treatments and no fire in Loudetia simplex.
Table 5. Contribution of the grouped demographic parameters (Fec.: fecundity, Sta.: stasis, Ret.: retrogression and Gro.: growth) to the demographic difference between treatments in each species, analysed by the life-table response experiment. When the demographic parameters are higher in treatment 2 than in treatment 1, contribution values are positive and vice versa. Treatment 1: the first of the two compared fire treatments; treatment 2: the second treatment; ∆𝜆 = λ(treatment 2) − 𝜆(treatment 1): the difference between the asymptotic growth rates (λ) of treatment 1 and treatment 2; C-value: the global contribution value calculated by summing all the entries of the contribution matrix from the comparison between the corresponding treatment. The (-) denote that there was no possible comparison between the other treatments and no fire in Loudetia simplex.
SpeciesTreatment 1 × Treatment 2∆𝜆C-ValueDemographic Parameters
Fec.Sta.Ret.Gro.
Andropogon canaliculatusmiddle fire × early fire
middle fire × late fire
middle fire × no fire
early fire × late fire
early fire × no fire
late fire × no fire
0.0091
−0.0238
−0.0531
−0.0329
−0.0623
−0.0293
0.0108
−0.0249
−0.0545
−0.0312
−0.0631
−0.0290
−0.016
0.006
−0.019
0.028
−0.003
−0.026
−0.095
−0.086
−0.085
0.020
0.008
−0.012
0.008
−0.013
−0.037
−0.015
−0.038
−0.014
0.114
0.069
0.088
−0.064
−0.029
0.024
Andropogonschirensismiddle fire × early fire
middle fire × late fire
middle fire × no fire
early fire × late fire
early fire × no fire
late fire × no fire
0.0041
−0.3157
−0.0656
−0.3198
−0.0697
0.2500
0.0036
−0.3150
−0.0630
−0.3172
−0.0681
0.2484
0.035
0.020
0.076
−0.013
0.035
0.066
−0.069
−0.361
−0.180
−0.281
−0.113
0.155
0.072
0.059
−0.017
−0.006
−0.084
−0.055
−0.034
−0.033
0.058
−0.015
0.094
0.082
Hyparrhenia diplandramiddle fire × early fire
middle × late fire
middle fire × no fire
early fire × late fire
early fire × no fire
late fire × no fire
−0.0026
−0.1698
−0.0199
−0.1672
−0.0173
0.1499
−0.0028
−0.1681
−0.0198
−0.1744
−0.0171
0.1519
−0.026
0.047
0.003
0.092
0.033
−0.053
0.011
−0.208
−0.108
−0.223
−0.122
0.108
−0.036
0.014
0.001
0.047
0.041
−0.017
0.048
−0.022
0.084
−0.091
0.029
0.113
Loudetia simplexmiddle fire × early fire
middle fire × late fire
middle fire × no fire
early fire × late fire
early fire × no fire
late fire × no fire
0.1890
−0.1138
-
−0.3028
-
-
0.1862
−0.1142
-
−0.3012
-
-
−0.051
−0.037
-
0.010
-
-
0.043
−0.120
-
−0.166
-
-
0.020
−0.010
-
−0.034
-
-
0.173
0.054
-
−0.110
-
-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Koffi, K.F.; N’Dri, A.B.; Srikanthasamy, T.; Lata, J.-C.; Konaté, S.; Konan, M.; Barot, S. The Demographic Response of Grass Species to Fire Treatments in a Guinean Savanna. Fire 2022, 5, 193. https://doi.org/10.3390/fire5060193

AMA Style

Koffi KF, N’Dri AB, Srikanthasamy T, Lata J-C, Konaté S, Konan M, Barot S. The Demographic Response of Grass Species to Fire Treatments in a Guinean Savanna. Fire. 2022; 5(6):193. https://doi.org/10.3390/fire5060193

Chicago/Turabian Style

Koffi, Kouamé Fulgence, Aya Brigitte N’Dri, Tharaniya Srikanthasamy, Jean-Christophe Lata, Souleymane Konaté, Marcel Konan, and Sébastien Barot. 2022. "The Demographic Response of Grass Species to Fire Treatments in a Guinean Savanna" Fire 5, no. 6: 193. https://doi.org/10.3390/fire5060193

Article Metrics

Back to TopTop