Got shrubs? Precipitation mediates long-term shrub and introduced grass dynamics in chaparral communities after fire

Background: Short-term post-fire field studies have shown that native shrub cover in chaparral ecosystems negatively affects introduced cover, which is influenced by burn severity, elevation, aspect, and climate. Using the southern California 2003 Old and Simi fires and the 2008 Sesnon Fire, we investigated the role of native shrubs in post-fire ecosystem responses across gradients of elevation, aspect, climate, burn severities, fire histories, and time. We collected field estimates of species cover in 2004 and 2015 at nested sampling sites. We used structural equation models with introduced and shrub cover as dependent variables. Results: Shrub cover in 2004 was most influenced by the number of reburns, while shrub cover in 2015 was most influenced by the time between the two most recent fires. In 2004, introduced cover was most influenced by burn severity in 2003; similarly, in 2015, introduced cover was most influenced by burn severity in 2008. In both one and twelve years post fire, average precipitation increased the length of time between fires and decreased the number of times a site burned. This direct reduction in the number of times a site had burned due to average precipitation resulted in lower shrub cover one and twelve years post fire. Additionally, mean annual precipitation increased burn severity one year post fire, which resulted in lower introduced cover. However, this indirect relationship between precipitation and introduced cover through burn severity was no longer present twelve years post fire. Shrub cover increased with a longer average time between fires twelve years after fire. Shrub cover did not mediate any indirect relationships between burn severity or fire history metrics and introduced cover in either year, suggesting competitive exclusion of introduced species by shrubs. Conclusion: Our research found that significant fire effects on shrub and introduced species are often mediated by precipitation. Precipitation trends are likely to change fire regimes and thus alter plant community dynamics.


Background
Plant community dynamics after fire are influenced by a suite of factors, including climate, fire characteristics, and non-living groundcover. Post-fire precipitation can affect plant community responses to fire by altering plant growth and survival, resulting in changes in plant community composition and successional trajectories (Keeley et al. 2005b). Additionally, plant establishment after fire can be influenced by non-living groundcover, such as litter, rock, and bare ground. For example, high litter and rock cover may limit post-fire soil erosion (Shakesby et al. 1994) and provide microsites for seed germination (Schlesinger and Gill 1978). Local conditions, such as burn severity, most recent fire return interval, as well as plant canopy, can also affect soil heating, reducing the seed bank and fuel load for subsequent fires (Odion and Davis 2000). Fire characteristics, such as burn severity, can impact future plant communities by affecting fire behavior and effects within subsequent fires (Parks et al. 2014a). Furthermore, repeat burns may affect a site's invasibility and thus ecosystem trajectory (Keeley and Brennan 2012).
Historically, wildfire in chaparral ecosystems burned with high severity through closed shrub canopies (Borchert andOdion 1995, Keeley andFotheringham 2001) every 30 to 150 years (Van de Water and Safford 2011), generally from July to August, with lightning as the primary ignition source (Keeley and Fotheringham 2001). At 15 years post fire, leaf area index is maximized, and shrubs can burn with high intensity, resulting in stand-replacing crown fires (McMichael et al. 2004). The chaparral biotic community evolved with these strong selective pressures under its historical fire regime (Keeley et al. 2011). However, the historical fire regime has changed due to fire suppression, the introduction of non-native plant species, urban growth into previously shrub-dominated landscapes, both humancaused ignitions and climate change extending the fire season to year-round, and increasing fire size and frequency (D'Antonio and Vitousek 1992;Moritz 2003;Keeley et al. 2005c;Van de Water and Safford 2011;Molinari et al. 2018;Syphard et al. 2018aSyphard et al. , 2018b. These changes to the historical fire regime have led to changes in the structure, function, and range of the chaparral plant community Keeley 2005, Molinari et al. 2018).
Pre-fire plant communities may also influence how a post-fire plant community responds to wildfire. In chaparral, shrubs are highly flammable due to high resin content, woody stems, standing dead branches, and the accumulation of litter beneath their bases. Some shrub species are adapted to survive high-severity fires either by resprouting (obligate resprouter), having high seed survival in the seed bank (obligate reseeder), or both (facultative seeder). Areas that burn more frequently are at a higher risk of extirpating the native plant community in favor of introduced species, as many native shrubs cannot survive frequent fires. For example, Adenostoma fasciculatum Hook. & Arn. will not resprout without an adequate growth period between fires (Zedler et al. 1983). Increases in fire frequency may lead to lower burn severity during subsequent fires and thus higher introduced species invasion (Keeley et al. 2005d, Keeley et al. 2008. Following a fire, succession in the plant community is dominated by residual species that survived the disturbance either by resprouting or germinating from the seed bank. Colonizers also play an essential role in the biodiversity of early post-fire communities and typically arrive onto a site in the first or second growing season. Short-term field studies after fire in chaparral ecosystems have shown that native shrub cover is negatively associated with both introduced forb and graminoid cover, and that introduced species are influenced by burn severity, elevation, aspect, and precipitation (Keeley et al. 2005a, Grace and Keeley 2006). Specifically, Keeley et al. (2005a) found that alien cover is negatively associated with woody plant cover five years after fire. In early post-fire years, shrubs are the primary functional group to both resprout and germinate in chaparral communities (Keeley et al. 2005a). However, introduced annual grasses are also rapid post-fire colonizers (D' Antonio and Vitousek 1992). Therefore, the dynamic between shrub and annual grass establishment is likely to determine long-term successional trajectories.
Few to no studies have examined the direct and indirect relationships between fire history, burn severity, years between fires, non-living groundcover, time, and precipitation that mediate post-fire chaparral community dynamics. Using hypothesized relationships, as well as other studies, we created a structural equation model (SEM) to understand the relationships between these variables (Fig. 1, Table 1). SEM's are a powerful statistical tool designed to determine causal multivariate relationships and uniquely allow for examining simultaneous multivariate relationships (Grace andKeeley 2006, Dell et al. 2017). We examined these relationships in a relatively intact montane chaparral ecosystem that burned typical of past chaparral fire regimes, as well as in an interior chaparral ecosystem that has had more frequent fires and a greater invasion of non-native species. We hypothesized that: 1) the amount of precipitation would directly affect burn severity and fire history, as well as shrub and introduced species cover; 2) fire history metrics (time between the most recent two fires and number of reburns) would indirectly affect shrub and introduced species cover; and 3) shrub species would adversely affect introduced species, as mediated by precipitation, burn severity, fire history metrics, and non-living groundcover variables. Using SEM to evaluate species responses across interacting burn severity, fire return interval, precipitation, and invasion gradients may provide insight into post-fire successional trajectories.

Field measurements
Immediately after the 2003 fires, Burn Area Reflectance Classification maps (Clark et al. 2003) were used to select six sites across burn severity levels (low, moderate,  (Parks et al. 2014b).
2 (+) Precipitation increases the number of fires at a site by increasing site productivity and thus fuel loads (Davis and Michaelsen 1995). 6 (+) Rocky topography provides favorable microsites for young shrub stands (Schlesinger and Gill 1978).

7
( −) Persistent negative relationships between introduced species richness and rock cover were found over five years (Keeley et al. 2005d).

8
(+) Fire severity is high where plant cover is high (Grace and Keeley 2006).

9
( −) The relationship between dNBR and non-native cover is negative ).

10
(+) Burn severity increases with the length of time since the last burn (Parks et al. 2014a).

11
(−) Fire severity is significantly lower when fires burn within a previously recorded fire perimeter (i.e., it reburned; Parks et al. 2015).

13
(−) Sites that burn twice in four years have higher introduced cover then sites that burn once in four years (Keeley and Brennan 2012).

15
(−) Introduced cover increases with younger pre-fire stand age (Keeley et al. 2005d or high) at each fire, with designations verified by field observations of vegetation canopy condition; predominantly green, brown, or black vegetation crowns were considered to be burned at low, moderate, or high severity, respectively ). Prior to 2015 field sampling, the burn areas were stratified by burn severity level, elevation, and transformed aspect (Roberts and Cooper 1989); two additional sites were added at the Simi Fire to fill two strata not represented by the initial field sites. Five 1 m 2 plots were sampled within each 60 × 60 m site; four plots were placed around a central plot in orthogonal directions oriented with the prevailing slope and plots were spaced 30 m apart to include three to four adjacent Landsat satellite image pixels. The center plot was monumented, and the sites were revisited in 2004 and 2015 (Bright et al. 2019). Understory vascular plant cover was ocularly estimated at each plot; cover included tree seedlings and saplings less than 1.37 m in height. All plants were identified to species, genus, or family. Each plant identified was assigned to a functional group (forbs, graminoids, shrubs, or trees) and designated as either native or introduced based on the USDA PLANTS Database (USDA NRCS 2016).

Statistics
A conceptual model to integrate the relationships among burn severity, fire history metrics, climate, plant community, and abiotic variables was developed from preliminary analyses and previous studies (Fig. 1). Based on this conceptual model, we tested direct and mediated (indirect) relationships between our collected variables (Grace 2006). Variables included in the hypothesized model were first examined for normality using the Shapiro-Wilks test of normality and then examined for skewness and kurtosis. All variables failed the Shapiro-Wilks test for normality (P < 0.05). To account for non-normality, the variable standard error and structural equation model test statistic were calculated by bootstrapping with repetitions set to 1000. Estimates of direct relationships were calculated using the lavaan package version 0.5-23.1097 (Rosseel 2012) in R version 3.5.0 (R Core Team 2018). Estimates and significance of the mediated relationships were also obtained using the lavaan package by explicitly defining interactions of interest using the test for mediation as outlined by Grace (2006). To account for the two sites that were added to the Simi Fire in 2015, we used full information maximum likelihood estimation. Autocorrelation among plots within sites was tested using Moran's I (Plant 2012); we found no autocorrelation, so the plot was used as the replicate throughout this study.
Correlations in the model residuals were computed from the SEM. Correlations in model residuals (< 0.10) were then evaluated for exclusion by comparing the Akaike information criterion (AIC) of models with and without the variable. If no significant reduction in AIC was found by excluding the variable, the more parsimonious model was kept. Next, correlations between variable errors were evaluated using modification indices; errors with indices greater than 3.84 were included in the model. Competing models were evaluated for a significant reduction in AIC. The final model fit was evaluated for acceptability based on no significant difference between the expected and observed covariance, a comparative fit index of 0.90 to 1.00, and root mean square error of 0.00 to 0.10. Direct relationships were significant when P < 0.05. Significant bivariate relationships found in SEM were further examined using Spearman rank correlation coefficients.
For fire history information, we obtained CalFire's fire history GIS layer. This was used to determine how many times a site had burned since 1878, which is when the database records began (CalFire 2017). From this we calculated the time between the two most recent fires and the number of times a site had burned. These two metrics were used in the subsequent analysis to represent fire history and are referred to as fire history metrics throughout text. We tallied the number of fires and time between the two most recent fires one year after the Old and Simi fires (2004), and again in 2015, which captured the Sesnon Fire reburning some of the Simi Fire sites. To indicate burn severity, we used differenced normalized burn ratio (dNBR), which was generated from Landsat 5 Thematic Mapper imagery ). The dNBR was calculated in 2003 for the Old and Simi fires and again in 2008 for the Sesnon Fire.
Climate variables were interpolated at 1 km resolution using a thin-plate spline model (Hutchinson 2000) of the western United States (Rehfeldt 2006, Rehfeldt et al. 2015. To capture finer scale variation in climate (30 m), we applied the model, which was based on a resolution of~1 km, to annual climatic data for the years 1981 to 2010 using the digital elevation model from the Shuttle Radar Topographic Mission (Farr et al. 2007). At each plot, climate normals were extracted and calculated from 2003 to 2010, including mean annual temperature, mean annual precipitation, mean maximum temperature in the warmest month, mean minimum temperature in the coldest month, growing season precipitation (Apr to Sep), winter precipitation (Nov to Feb), summer precipitation (Jul to Aug), and spring precipitation (Apr to May).
We calculated residuals for the climate normals at each fire separately by subtracting the mean at each fire to account for site differences. Spearman correlation coefficients and corresponding P-values were obtained for the relationships between all climate variables. All climate variables were highly correlated. Unlike temperature, chaparral fire occurrence is directly dependent on precipitation variability (Dennison and Moritz 2010), and this variability can mitigate modeled temperature changes under future climate scenarios (Bachelet et al. 2016). Thus, mean annual precipitation residuals were selected to broadly represent climate (Fig. 2).
Three field-collected abiotic variables were evaluated using Spearman's correlation tests for inclusion in the model: percentage cover of litter, rock, and bare ground. Litter and bare ground (2015: r = −0.43, P = 0.0002), rock and bare ground (2015: r = 0.43, P < 0.0001), and rock and litter (2015: r = −0.45, P < 0.0001) were significantly correlated. Bare ground was removed from consideration due to its correlation with litter and rock. Litter cover was removed from consideration as it was significantly correlated with other plant metrics (introduced cover 2015: r = 0.41, P = 0.0004). Rock cover was thus selected as the non-living ground cover variable to include in the model.
Throughout the text and figures, variables are labeled with the year the data was collected. For example, time between the two most recent fires in 2004 used a data set that included all fires up to 2004, while time between the two most recent fires in 2015 used a dataset that included all fires up to 2015.

Fire history and abiotic factors
Mean dNBR in 2003 was 301 (range = −215 to 758), while mean dNBR in 2008 was 158 (range = −84 to 670). In 2003, the Old Fire had a mean dNBR of 469 (range = 130 to 758), while the Simi Fire had a mean dNBR of 135 (range = −215 to 401). After the Sesnon Fire in 2008, the mean dNBR at the Simi sites was 302 (range = −41 to 670), while at the Old sites, which did not burn that year, the mean dNBR was −33 (range = −84 to 62). The range of time between the two most recent fires was the same in 2004 and 2015 (4 to 125 years), although the median dropped from 22 years in 2004 to just 5 years in 2015, primarily due to the Sesnon Fire reburning some of the Simi Fire sites only five years after the Simi Fire. The number of reburns in 2004 ranged from 1 to 5, with a mean of 2.9 burns, and in 2015 ranged from 1 to 6 burns with a mean of 3.4 burns. Rock cover was similar at the Old and Simi fire sites in 2015 with a mean of 7.3% and 7.2%, respectively, and ranged from 0 to 60% at the Old Fire sites and 0 to 75% at the Simi Fire sites (Additional file 1, Fig. 1).

SEM results
Our initial SEM (Fig. 1) yielded a highly discrepant, singular model (comparative fit index = 0.972, root mean square error approximation = 0.186, chi-square = 30.8, degrees of freedom = 9, and P < 0.0001). Using the residual covariance matrix from the SEM, we were able to retain relationships within this model with high residual covariance. Our final inferential model yielded no significant discrepancies between the model and the data (comparative fit index = 0.997, root mean square error approximation = 0.036, chi-square = 30.52, degrees of freedom = 28, and P = 0.34; Fig. 3, Additional file 3, Fig.  1). Many of the direct relationships originally hypothesized for exploration using SEM were not supported by the collected data. The unsupported relationships were excluded from the final inferential model; these excluded relationships were mean annual precipitation not significantly predicting shrub cover in either years nor introduced cover in 2004, rock cover not predicting shrub cover, shrub cover and introduced cover not predicting dNBR or fire history metrics, and, with the exception of 2003, dNBR predicting 2008 dNBR, all relationships between dNBR and fire history metrics were accounted for in the model's residual covariance structure due to high correlations. Our model explained 8% of the variation in shrub cover in 2004, which increased to 21% by 2015 (Fig. 3). The model explained 35% and 63% of the variability in introduced cover in 2004 and 2015, respectively. Based on standardized estimates (Table 2, Fig. 3), the number of reburns was the most influential variable affecting shrub cover in 2004, while the time between fires was most influential variable affecting shrub cover in 2015. The most influential variable affecting introduced cover in 2004 was dNBR in 2003; similarly, introduced cover in 2015 was most influenced by dNBR in 2008 (Table 2, Fig. 3).
Mean annual precipitation, an exogenous variable not predicted by other variables in the model, directly affected dNBR and fire history metrics. Overall, sites that were wetter tended to burn with higher burn severity ( Fig. 4A and B), had longer times between the two most recent fires ( Fig. 4C and D), burned less frequently (  Table 2). The relationship between mean annual precipitation and introduced cover in 2004 was mediated by dNBR in 2003; sites that were wetter had higher dNBR and lower introduced cover (Fig. 3, Table 3). The relationship between mean annual precipitation and shrub cover in 2004 and 2015 was mediated by the number of reburns in 2004 and 2015; wetter sites tended to burn less often, resulting in lower shrub cover (Fig. 3, Table 3). Rock cover in 2015 had a negative, direct effect on introduced cover in 2015 (Fig. 5), but there were no other significant direct or indirect relationships with rock cover in the model (Fig. 3, Table 2).
We found a direct, negative effect of 2003 dNBR on introduced cover in 2004 (Fig. 6A), and a direct positive effect of 2008 dNBR on introduced cover in 2015 (Fig.  6B). There was a direct negative relationship between 2003 dNBR and 2008 dNBR (Fig. 6C). We also found that this relationship between 2003 dNBR and 2008 dNBR indirectly affected introduced cover in 2015 (Fig.  3, Table 3). Overall, sites in 2003 had higher dNBR than sites in 2008, which resulted in higher introduced cover by 2015.
Time between the two most recent fires had a positive effect on shrub cover in 2015 (Fig. 7A) but not in 2004 (Fig. 3, Table 2). In both 2004 and 2015, there was not a significant direct or indirect relationship between the introduced cover and time between the two most recent fires (Fig. 3, Table 3). The number of reburns had positive effects on shrub cover in both 2015 (Figs. 3 and 7B, C; Table 2). We did not find a significant direct or indirect relationship between the number of reburns in 2004 and 2015 and introduced cover in 2004 and 2015 (Fig. 3, Table 3). For plant-based variables, introduced cover in either year was not predictive of any of the variables in the model (Table 2), although shrub cover in 2004 and 2015 had a significant negative relationship with introduced cover in 2004 (Fig. 8A) and 2015 (Fig. 8B)(Fig. 3, Table 2). Shrub cover did not significantly mediate any indirect relationships between dNBR, fire history, and introduced cover (Table 3).

Discussion
We predicted that climate would have direct effects on burn severity (dNBR) and fire history metrics, resulting in changes in plant communities. Indeed, mean annual precipitation had direct effects on burn severity and fire history metrics in this study. Wetter sites tended to have higher burn severity, burn less often, and have a longer time between fires, which directly and indirectly affected shrub and introduced cover.
Abiotic and biotic factors influence burn severity, including climate, fuel loads, pre-fire vegetation composition, and time since the last fire (Parks et al. 2015). In many ecosystems, including chaparral, if fire re-occurs before the ecosystem has had adequate time to recover, the severity of the subsequent fire is limited (Parks et al. 2014a, Parks et al. 2015. Indeed, higher burn severity in the first fire resulted in lower burn severity in the subsequent fire, and the maximum burn severity at the 2008 sites was reduced, which we presume was due to lower fuel loads. Also, high mean annual precipitation increased burn severity in both 2003 and 2008, likely due to higher fuel loads in those areas. Although time between the two most recent fires and the number of reburns are strongly correlated, they have slightly different implications for ecosystems. Time between the two most recent fires impacts stand age, species composition, and invasibility, while the number of reburns affects a site's resiliency to fire (Parks et al. 2014a, Parks et al. 2015. We found that mean annual precipitation increased the time between the two most recent fires. This may be due to increases in live fuel moisture, which has a limiting effect on a site's ability to burn (Dennison and Moritz 2010). Time between the two most recent fires determines how much fuel is available to burn, which in turn significantly influences fire intensity and burn severity. Many shrub species, such as the dominant and perennial Ceanothus oliganthus, Ceanothus megacarpus Nutt., and Adenostoma fasciculatum, can be entirely or nearly eliminated from sites with a shorter than historic fire return interval (Zedler et al. 1983, Jacobsen et al. 2004). This potential extirpation can be exacerbated by the fact that native obligate seeders such as Ceanothus megacarpus and Ceanothus crassifolius thrive in arid locations where fires have the potential to burn more frequently (Kolb andDavis 1994, Pratt et al. 2014). While short time intervals between fires can be a detriment to the native plant community, introduced grasses with an annual life history strategy can typically survive short intervals between fires (Whisenant 1990).
Lower mean annual precipitation also decreased the number of reburns, acting as a top-down control by increasing live fuel moisture, which limits its ability to burn (Dennison and Moritz 2010), and decreases its invasibility by the more arid-dependent, introduced Bromus species (Germino et al. 2016). In chaparral ecosystems, fine fuels created by high introduced-species cover can increase fire frequency (Keeley 2000). Frequent ground fires can impact native species, as many native shrubs are unable to survive during short fire return intervals caused by frequent burns (Keeley et al. 2003,  Keeley 2006, Keeley and Brennan 2012, Meng et al. 2014). Thus, introduced species drive increases in fire frequency that may result in community state changes (Zedler et al. 1983).
Few studies have considered long-term plant community dynamics in relation to fire. Although shrub cover remained similar over time, the number of shrub species doubled between 2004 and 2015. The most common shrub species one and twelve years post fire were facultative seeders, with Adenostoma fasciculatum being the most frequently encountered. Shrub diversity, in the first growing season after fire, is primarily attributed to  species that survived the fire and those that benefit from fire, such as facultative and obligate seeders and obligate resprouters (Moreno and Oechel 1991). Shrub diversity in chaparral after the second growing season is a dynamic process reliant on spatial patterns of resources, especially precipitation (Keeley et al. 2005b). The number of introduced species doubled over twelve years, with Bromus diandrus consistently common over that period. Keeley et al. (2005a) found that introduced species diversity peaked two to five years post fire, although they only studied the first five years post fire. Our results suggest that regardless of when introduced species diversity peaks, numerous introduced species can maintain a population over time.
As predicted, we found that the direct relationships between mean annual precipitation and fire had indirect consequences on plant communities, but these were not consistent across time, burn severity, fire history, or plant functional groups. First, the number of reburns was lower in areas of high precipitation, which resulted in lower shrub cover both one and twelve years post fire. Furthermore, we found that areas with high precipitation also had high burn severity, which resulted in lower introduced cover, but this indirect relationship was only significant one year after fire. We also found that shrub species and not introduced species benefited when sites had a longer than average time between fires, although this was only evident twelve years post fire, not one year post fire.
Rock cover affects percolation and thus site hydrology (Urbanek and Shakesby 2009); creates a more heterogeneous surface that may limit post-fire erosion (Shakesby et al. 1994); provides microsites for seed germination (Schlesinger and Gill 1978); and alters soil nutrients (Houlton et al. 2018). Additionally, species richness has been found to be positively associated with rock cover in chaparral ecosystems (Grace and Keeley 2006). There were no significant relationships between rock cover and burn severity or fire history metrics in this study. However, similar to Keeley et al. (2005d), who found a negative relationship between rock and introduced cover one, two, and five years post fire, we found that rock cover decreased introduced cover but only twelve years post fire. We also found that this relationship between rock and introduced cover was most likely masked by a significant relationship between burn severity and introduced cover, as did Keeley et al. (2005d). In our study, during the first growing season, burn severity had such a significant impact on introduced cover that there was little room for other influences. These results suggest that rock cover can have an impact on introduced cover, and the timing and magnitude of these effects depend on other environmental variables.
Introduced cover was also affected by mean annual precipitation, burn severity, fire history metrics, and shrub cover, resulting in complex interactions. Higher mean annual precipitation decreased introduced cover in year twelve, which is not surprising given that our sites were predominantly invaded by Bromus species, which flourish in more arid ecosystems (Germino et al. 2016). However, we did not find that mean annual precipitation influenced introduced cover during the first growing season after fire, which we attribute to low introduced cover during the first year, and colonization, which usually peaks two to five years post fire (Keeley et al. 2005a).
The positive effects of mean annual precipitation on burn severity immediately after fire resulted in lower introduced cover one year after fire. Others have found that higher severity favors shrub species, which can then outcompete introduced species (Keeley et al. 2005d). However, the positive effects of mean annual precipitation on 2003 burn severity reduced 2008 burn severity, and thus increased introduced cover in year twelve. Higher burn severity in year one likely reduced fuel. Combined with a shorter fire return interval of five years, which was likely detrimental to native shrubs, this allowed introduced species to increase after two successive burns. Surprisingly, time between the two most recent fires and number of reburns had no direct effects on introduced cover. Climate and fire effects on introduced cover were mediated through shrub cover. Like others, we found that shrub cover reduced introduced cover (Keeley and Brennan 2012, Condon and Pyke 2018), often through complex interactions with precipitation, burn severity, and fire history metrics. Overall, sites that had longer time between fires and more reburns had higher shrub cover, simliar to what others have found. Zedler et al. (1983) showed that longer time between fires benefit native shrubs, and McMichael et al. (2004) used remote sensing to show that leaf area index is maximized at 15 years post fire, a timeline similar to our study. These results suggest that the sites in our study were within the normal range of variability for time between fires and number of reburns. This study's results do not refute the cause for concern that Zedler et al. (1983) addressed, which was that reburning frequency higher than observed in this study could lead to a state change, due to shrubs being less tolerant than introduced species to more frequent fires.

Conclusions
Post-fire ecology studies typically focus on the direct relationship between fire metrics and plant community within five years after fire. Our study was unique in that it compared these direct relationships, as well as indirect relationships, between burn severity, time between the two most recent fires, number of reburns, climate, abiotic variables, and the plant community during the first and twelfth growing season after fire, and across variable fire histories. We found that mean annual precipitation was a driving factor in chaparral fire regimes, generally resulting in positive effects on native shrubs and adverse effects on introduced species. Thus, changes to climate will result in changes in chaparral fire regimes, which may have negative consequences for native species. We found that these fire effects on native and introduced species changed over time, implying that long-term monitoring is essential for understanding the successional trajectory of a chaparral site. This may be especially true at sites that receive precipitation outside (more likely less than) the normal range of variability. Finally, we found that shrub cover was negatively associated with introduced cover, which are both influenced by burn severity, either directly or indirectly, suggesting dynamic competitive exclusion between shrubs and introduced species.