Persistence of fire-killed conifer snags in California, USA

Background: Coarse woody debris has numerous functions in forest ecosystems, including wildlife habitat, fuel loading, and nutrient cycling. Standing dead trees, or snags, are particularly important resources for wildlife, and guidelines for snag recruitment and retention are integral parts of forest management. However, snag persistence is not well understood and existing models are poorly validated. We monitored fire-killed conifer snag persistence after four California wildfires and one prescribed fire. Species sampled included white fir ( Abies concolor [Gord. & Glend.] Lindl. ex Hildebr.), incense-cedar ( Calocedrus decurrens [Torr.] Florin), Jeffrey pine ( Pinus jeffreyi [Balf.]), sugar pine ( Pinus lambertiana Douglas), and ponderosa pine ( Pinus ponderosa Lawson & C. Lawson var. ponderosa C. Lawson). Results: Fire-injured trees that survived at least one year post fire were monitored for 10 years after they died. We developed two sets of species-specific snag persistence rates: a simplified model based only on time since death, and a full model that included tree characteristics. The primary factors affecting persistence rates in all models were snag age (years since death) and species. Secondary factors included diameter at breast height, time between fire and tree death, and crown injury. Snag persistence rates decreased with snag age. Yellow pine snags (a combined group of ponderosa pine and Jeffrey pine) fell significantly faster than any other species. Larger-diameter snags persisted longer than smaller-diameter snags for all species. Incense-cedar trees that survived longer post fire persisted longer as snags. Conclusions: This study provides land managers with guidance on snag longevity, which can aid in planning for long-term wildlife habitat, fuel loading estimates, and forest structural diversity.


Background
Dead trees are an important component of forest ecosystems. Either as standing snags or downed logs, they play vital roles in nutrient cycling and in carbon sequestration (Harmon et al. 1986;Krankina and Harmon 1995;Brais and Paré 2006). They also provide habitat and foraging sites for vertebrates and invertebrates, as well as growth sites for fungi, lichen, and bryophytes (Thomas 1979;Harmon et al. 1986;Bunnell and Kremsater 1990;Bull and Parks 1997). In the northwestern United States, about 9% of all wildlife species and over 25% of bird species utilize dead trees (R. Harger, USDA Forest Service, Northern Region, Missoula, Montana, USA, unpublished report; McClelland et al. 1979, Hicks 1983. At least 169 wildlife species in California utilize snags (California Department of Fish and Wildlife 2014). Nearly all snags larger than 10 cm DBH (diameter at breast height; 1.37 m above ground) can be used, although some wildlife species prefer specific tree species or sizes (Thomas 1979;Parish et al. 2010). In general, larger snags can support more wildlife species, though a diversity of snags with varying DBH, height, decay class, and species is important to manage for a wide variety of wildlife (Raphael and White 1984;Bunnell 2013).
Predicting snag longevity is important in many aspects of forest management, from quantifying habitat requirements, carbon stores, and fuel loading, to identifying safety hazards. Despite this importance, snag dynamics are still poorly understood. Studies have identified many different and often discrepant factors affecting snag persistence. Most studies have found snag age (i.e., time since death), and species significantly affect snag persistence, with some overlap between similar species. While most studies report that snag diameter affects longevity, with larger snags persisting longer than smaller snags (Keen 1955;Bull 1983;Raphael and Morrison 1987;Chambers and Mast 2005), some studies have found no impact of DBH (Schmid et al. 1985;Harrington 1996). Other factors including cause of death, tree height, soil type, site aspect, slope, microclimate, species composition, landform, salvage harvesting, and stand density have been identified as significant in some studies, but not significant in others (Keen 1955;Lyon 1977;Mitchell and Preisler 1998;Russell et al. 2006;Parish et al. 2010). Large-scale factors such as climate may also play an important role, but lack empirical evidence to describe their effects (Keen 1955;Harmon et al. 1986;Garber et al. 2005). In addition to these studies, the Fire and Fuels Extension of the Forest Vegetation Simulator (FFE-FVS), a commonly used forest growth model, includes submodels of snag dynamics (Reinhardt and Crookston 2003). These snag submodels are largely not validated, although a few have incorporated work from the Snag Recruitment Simulator and the Coarse Woody Dynamics Model (Marcot 1992;Mellen and Ager 2002).
Wildland fire burns millions of forested hectares annually in the western United States (National Interagency Fire Center: https://www.nifc.gov/fireInfo/nfn.htm) and is a major cause of snag recruitment (Raphael and White 1984). Several studies indicate that the snag longevity of fire-killed trees is different than other causes of mortality (Raphael and Morrison 1987;Everett et al. 1999;Garber et al. 2005), although few exist that directly compare snags created by fire versus other means (Raphael and Morrison 1987;Morrison and Raphael 1993;Chambers and Mast 2005). Most suggest that fire-killed snags persist for less time than snags created by other mortality agents. These fire-killed snags persist for only a short time, with very few remaining after 10 years (Harrington 1996;Everett et al. 1999;Ritchie et al. 2013;Ritchie and Knapp 2014). Current post-fire management practices often include harvest activities to recoup economic losses from dead trees, to reduce long-term coarse woody fuel loading, and for hazard-tree abatement (Dunn and Bailey 2015;Peterson et al. 2015). Snag retention guidelines frequently focus on tree sizes and species anticipated to persist the longest. Improved understanding of post-fire snag persistence is needed to better predict long-term fire effects and to ensure that forest management objectives for snag density are met.
We monitored trees from five species for 10 years after death in four wildfires and one prescribed fire in California, USA, which burned between 2000 and 2004. There are limited studies on the snag dynamics of fire-killed trees, and even fewer that include the species monitored here or that account for tree-level fire injuries (Table 1). Our objectives were: (1) to quantify post-fire snag persistence, (2) to determine which snag characteristics most affect persistence, and (3) to evaluate the predictive accuracy of the snag submodel in FFE-FVS.

Study area and sampling
The data and sampling methods used here are largely a subset of those described in detail in Hood et al. (2010), which were used to develop post-fire tree mortality models. Briefly, we sampled in five burned areas (Cone, McNally, Power, and Storrie wildfires, and Prattville prescribed burn) across a wide geographical area in California extending from the southern end of the Cascade Range to the southern end of the Sierra Nevada mountain range. All fires except Prattville, a mixed-severity prescribed fire, were reported in Hood et al. (2010).
These areas burned between 2000 and 2004. All fires were in the Sierra Nevada mixed-conifer forest type (SAF Type 243;Tappeiner 1980) with the exception of the Cone Fire, which occurred in the interior ponderosa pine type (SAF Type 237; Barrett et al. 1980). Personnel from local forest districts identified areas of mixedseverity fire within each fire boundary for possible sampling sites. These areas burned primarily as surface fires, with the majority of trees retaining some green needles. We selected trees during the summer of the year following each fire based on categorical bins in an effort to assess a full range of tree and injury characteristics. Specific categories included DBH classes from 25 to 50 cm, 51 to 76 cm, and ≥77 cm; crown kill ranges ≤10%, 11 to 50%, 51 to 60%, 61 to 70%, 71 to 80%, 81 to 90%, and 91 to 100%; and cambium kill categorical ratings of 0, 1, 2, 3, and 4. Any tree that survived at least one year post fire was monitored annually for five years or until death, then for 10 years post mortality. Only trees with 10 years of post-mortality observations were used in this analysis. One year post fire, we assessed tree and injury characteristics including species, DBH, tree height, pre-fire crown length, number of red turpentine beetle (RTB; Dendroctonus valens LeConte) pitch tubes for pines, and fire injury. We quantified fire injury by recording cambium kill rating (CKR; ordinal variable between 0 and 4 describing the number of dead cambium samples at each quadrant of the tree base), crown length killed (all species; observed portion of the crown with dead tissue from fire injury), and crown length scorched (Jeffrey pine and ponderosa pine only; observed portion of the crown with dead needles, but living buds and branches). Lengths were measured using a clinometer or laser hypsometer to the nearest 0.3 m. We calculated pre-fire crown length and crown length killed from the estimated prefire crown base height, post-fire crown base height, and tree height measurements. The pre-fire crown was estimated using the methods described in Ryan (1982), by visually reconstructing the crown to its pre-fire state based on crown scorch and fine branch structure. We converted crown length killed into percent crown length killed, by dividing the crown length killed by the pre-fire crown length to create a crown injury variable (CI). When crown length killed was not available, we used an estimated live crown measurement to determine percent crown length killed. Snag persistence may be influenced by whether it is in isolation or in a group of living trees (Mitchell and Preisler 1998;Russell et al. 2006). We (1) = factor had a significant relationship with snag longevity, which varied by factor level; (0) = factor had no significant relationship with snag longevity therefore also assessed whether each snag was in a stand that had been salvage logged (Sal) and if the snag was surrounded by live conifers (protected by green conifers; PGC) as binary categories (0 = no and 1 = yes). Each year after this initial assessment, trees were monitored for live or dead status and break height, if applicable. Snags were monitored until they fell or broke below 3.65 m. This threshold was chosen as a combination of the taller snag requirements in the Lassen National Forest Resource Management Plan (USDA Forest Service 1992), which recommends that suitable wildlife snags be at least 6.1 m tall, and studies of snag nesting heights and foraging use that suggest shorter snags are poorly utilized. For example, Ganey (2016) found that snags of 3 m or more held 95% of nesting cavities. In a monograph of 18 cavity-nesting bird species in the Sierra Nevada, California, mean nest snag height was 12.3 m for all species combined and the lowest individual species mean nest snag height was 3.8 m (Raphael and White 1984). In a review of snag use by mammals and birds, only one bird species had a mean nest height shorter than 5 m (Bunnell et al. 2002). Ponderosa pine and Jeffrey pine were combined into one group during data collection and are referred to as yellow pine hereafter.

Analysis
Descriptive analyses by standing and downed snag categories were performed for each species, including range, mean, and standard error of the mean for DBH, CI, and YTD (time between the year of the fire and the year of tree death). These were compared within and among species using a Wilcoxon rank-sum test to determine if there was a difference between the characteristics of downed and standing snags. With this and all other analyses, a P-value of less than 0.05 was used to define statistical significance unless otherwise noted. Preliminary results showed marked species differences; we therefore evaluated species separately. Analyses were performed with the R Statistical Program v. 3.3.1 (R Development Core Team. 2016).

Data censoring
Snags were evaluated for 10 years after death. Over 35% of the snags were still standing by the end of the observation period. These snags were considered rightcensored (i.e., snags were followed for a set amount of time rather than until all snags in our sample population fell). In addition to the right-censoring, several snags had years in which they were not measured. If these snags were standing before the missed observation, but down at the next observation, the censoring was assumed to be non-informative and those snags were considered to be standing for the missing years. This is the standard protocol for such interval censoring (Kleinbaum and Klein 2010).

Snag longevity trends
We first fit unique models to each species for both overall persistence rates and individual fall risk using time since fire as the only independent variable (time-only model). The time-only models can be used to estimate stand-level snag persistence. We fit five common distributions-Weibull, exponential, logistic, log-logistic, and Gaussian-to each species' persistence time to determine the general shape before developing final time-only models. When comparing maximum log-likelihood, Weibull and log-logistic were the top models for all species. However, the likelihood ratio was not significantly different for the two distributions in the case of incensecedar. We chose to use the Weibull function because the heavier tail predicted with the log-logistic distribution cannot be evaluated within the timeframe of our observations and because the hazard function, or the probability that a snag will fall within the next time unit (in this case, one year) given it has survived until a certain time, is monotonically increasing for the Weibull distribution (where the shape parameter is greater than 1), but is nonmonotonic (first increasing, then decreasing) for log-logistic distributions. Ecologically, the older a snag is, the more likely it is to fall before the next observation, making the Weibull function a more logical fit. To ensure that differences in average DBH between species in our sample were not confounding the timeonly models, we also fitted species-specific models using a subset of snags centered on a mean of 52 cm DBH and compared to the time-only models with the full datasets. The plotted curves exhibited similar trends (data not shown).
Weibull curves can be parameterized in several different ways to achieve equivalent results. For the survival analysis to predict the standing time of a snag after death, we used the accelerated failure time (AFT) model (Eqs. 1-3; Qi 2009) rather than the proportional hazards model. Although both forms use the same standard notation, the parameter interpretations are not equivalent. The resulting AFT parameters can be used to determine the hazard of falling from one year to the next, the probability of persistence to a certain age, and the age of snag-fall given a certain probability threshold. In the AFT model, the covariate accelerates or decelerates the life course (i.e., how much time passes before a snag transitions from standing to fallen) by a constant. In addition, parameter coefficient estimates in AFT regressions are robust to unknown, unaccounted for covariates and are not strongly affected by probability distribution choice.
Equation 2 is the hazard rate (h) at a given time (t) using the same variables as Eq. 1: Equation 3 is survival time (t) given a percentile (p) of the survival distribution: Factors affecting snag longevity We next developed individual snag-persistence equations with time since fire and additional covariates as independent variables (covariate model) through a combination of purposeful variable selection, model averaging, and information theoretic selection. Variables tested as covariates included CKR, DBH, CI, RTB, PGC, YTD, snag height, and species. The covariate models can be used to identify specific individual snags within a stand that are likely to persist the longest. Preliminary main effects variables were chosen by univariate significance (P < 0.25) or known ecological importance, then narrowed to those that were statistically and logically conditionally independent from other variables. We additionally assessed the impact of random effects of individual variability and fire (e.g., location) on snag persistence. A full model was then developed using model averaging to determine the suite of models that had the best Akaike's information criterion (AIC). Model-averaged importance of terms and model weights were calculated. All variables with importance greater than 0.8 were included, and any terms included in the models as a covariate but not as a variate were included as the latter for interpretability. Although these averaged models provided the likely best predictive model, they are less interpretable and tend to have many more terms than a more practical model. Terms were thus selected from these full equations using both a backward and a stepwise AIC procedure until the difference in AIC of the simpler model and the full model fell to less than 7 (Anderson and Burnham 2002;Burnham et al. 2011). We chose a model from this suite of best models for each species that maximized accuracy, goodness of fit, and simplicity, while utilizing the most reliable and easy to collect variables. Following term selection, the models were bootstrapped with a set seed 2500 times with 20% replacement to estimate coefficients. Estimates are the median bootstrap value minus the bias estimate. Estimated stand-level persistence was calculated using the median values of each covariate. The estimated probability of persistence at year 10 was calculated for each individual snag to assess the accuracy of the models. Snags with a probability of persistence greater than or equal to 0.5 (i.e., ≥50% chance of snag persisting through a given year) at year 10 were considered predicted to be standing and those with a probability of persistence less than 0.5 were considered predicted to be down.

FFE-FVS snag submodel evaluation
The Forest Vegetation Simulator (FVS) is widely used to predict stand dynamics in response to time, management, disturbance, and natural succession. The FFE-FVE integrates FVS with models of fuel dynamics, fire behavior, and fire effects. In addition, it includes a snag submodel to predict the breakage, decay, and persistence rates of snags based on snag height, decay status, diameter class, and a species-and variant-specific fall rate multiplier. The last 5% of snags in a stand are predicted to fall at a rate of 0.01 per year, with a maximum persistence of 100 years (Reinhardt and Crookston 2003). Other than this change, FFE-FVS does not account for snag age when calculating persistence rates. FFE-FVS does not take cause of mortality into account, although there is an option to simulate a fire in the snags. The snags in this study were from two variants of FVS: Inland California and Southern Cascades (ICASCA) and Western Sierra Nevada. The fall rate multipliers for these variants are the same when rounded to two decimal places, at 0.69 for incense-cedar, 0.88 for white fir, and 1.24 for yellow pine and sugar pine. We compared observed yearly cumulative snag persistence to those predicted by FFE-FVS.

Results
A total of 1140 snags were followed for 10 years after death: 96 incense-cedar, 189 sugar pine, 308 white fir, and 547 yellow pine. At 10 years post mortality, 13% of incense-cedar, 42% of sugar pine, 46% of white fir, and 95% of yellow pine had fallen ( Figs. 1 and 2). Average initial break heights were similar across species examined: 0.6 m for incense-cedar, 3.5 m for sugar pine, 3.2 m for white fir, and 1.0 m for yellow pine. Of the 869 snags that broke during the 10 years that they were followed, 613 snags (71%) broke at groundline. All species datasets covered the full range of fire-caused injuries and a wide range of diameters (Table 2). We only report summary statistics for variables that were significant in the best models; these included DBH, CI, and YTD (Table 2), and PGC. Standing sugar pine (DBH: P < 0.0001), white fir (DBH: P < 0.0001), and yellow pine (DBH: P = 0.001) were larger than downed snags, but DBH did not differ between standing and downed incense-cedar (P = 0.298). PGC did not differ between standing and downed snags of any species (P > 0.05): 95% of standing and 100% of downed incense-cedar, 77% of standing and 79% of downed sugar pine, 89% of standing and 93% of downed white fir, and 53% of standing and 45% of downed yellow pine snags were protected by green conifers. Crown injury did not differ for any species between standing and downed snags (CI: P > 0.05). Incense-cedar (YTD: P = 0.026) and yellow pine (YTD: P = 0.028) trees that survived longer post fire persisted longer as snags, while YTD was not different

Persistence trends
There were clear differences in species in snag persistence ( Fig. 1; symbols show observed values). Yellow pine fell fastest and incense-cedar persisted the longest. Sugar pine and white fir persistence rates were similar to each other and fell in between yellow pine and incense-cedar. The models fit for each species all differed from each other in scale, shape, or both, indicating that snag dynamics required species-specific equations ( Fig. 1; Table 3). The first set of models predicted persistence rates by only accounting for time since fire ( Fig. 1; dashed lines). The models estimated that 50% of snags would fall by 7 years after dying for yellow pine, 10 years for white fir, 11 years for sugar pine, and 14 years for incense-cedar. They predicted that 99% of snags would fall by 21 years for incense-cedar, 25 years for sugar pine, 20 years for white fir, and 12 years for yellow pine after dying. Although incense-cedar displayed the greatest lag time in persistence drop-off, the snag fall after this initial lag increased rapidly. In contrast, sugar pine snags initially fell faster than incense-cedar and at a similar rate to white fir, but the model estimated the rate of fall to slow over time, leading to a small portion of the snags persisting for >25 years after dying from fire. White fir displayed a moderate lag at both the beginning and end of the life cycle. Yellow pine showed almost no lag at the beginning or end of the snag life cycle; snag persistence probability decreased rapidly and consistently.  Table 3 Coefficient estimates of regressions (from Eqs. 1-3: α 1 , α 2 , etc., are the covariate coefficients corresponding to x 1 , x 2 , etc.; μ is the intercept; and σ is the scale); coefficient standard error estimate in parentheses. Variable codes are: diameter at breast height (DBH; cm); crown injury percent (CI); years to mortality post fire (YTD), and protected by green conifers (PGC; 0 = no and 1 = yes) by species and status. Species codes are: IC = incense cedar; SP = sugar pine; WF = white fir; and YP = yellow pine. Superscript numbers indicate model-averaged variable importance rank. ΔAIC is the difference between this model and the best AIC model. Those terms that are less than 0.8 in model-averaged importance are noted with an asterisk (*)

Factors affecting persistence
Several factors affected individual snag persistence besides time since death ( Fig. 1; solid lines). While time since death was the most important factor for all species, accounting for DBH, CI, and PGC allowed us to assess snag persistence on an individual basis. DBH was the most important covariate for white fir, sugar pine, and yellow pine; larger snags persisted longer than smaller snags (Table 3). White fir and sugar pine snags over 50 cm DBH and yellow pine snags over 90 cm DBH were likely to persist for at least 10 years post mortality (Fig. 3). CI and its interaction with DBH were also significant for white fir and yellow pine. Snags with lower crown injury persisted longer than those with higher crown injury in white fir greater than 50 cm DBH. Larger white fir and yellow pine snags with lower CI stood longest compared to similarly sized snags with high CI. Smaller white fir and unprotected yellow pine snags with higher CI persisted longer than similarly sized snags with low CI (Figs. 3c and 4). These interactions in pattern occurred around 50 cm DBH and 50% CI for white fir and 30 cm DBH and 50% CI for unprotected yellow pine. Yellow pine smaller than 30 cm DBH that were protected by green conifers lasted longer than unprotected snags (Fig. 4a, b); however, protection by green conifers dramatically decreased the probability of persistence of large yellow pine snags as CI increased (Fig.  4a, b). Years since death after fire was the only significant covariate for incense-cedar. Incense-cedar that survived longer post fire persisted longer as snags. The random effects of fire or individual variability were not significant in any species. The AIC of the recommended models for all species were within 7 ΔAIC of the best AIC model. All covariate models predicted that snags would stand longer than what the time-only models predicted. Sugar pine was particularly affected, with the final 1% of snags remaining to year 44 in the covariate model versus year 25 in the time-only model (Fig. 1). The covariate model overestimated yellow pine persistence by 31%, but was highly accurate (i.e., within 3%) for incense-cedar, sugar pine, and white fir in this dataset (Fig. 2). Model accuracy was assessed on an individual snag level at year 10 and threshold for probability of persistence of 0.5 (P < 0.5 downed, P ≥ 0.5 standing). The covariate models for yellow pine, sugar pine, white fir, and incense-cedar predicted 36.1%, 58.6%, 56.9%, and 87.4% of snags standing Fig. 3 Sample predicted snag persistence once a tree dies from wildfire (1 = standing; 0 = downed or <3.65 m tall) from covariate models for (a) incense-cedar, (b) sugar pine, and (c) white fir. Year to death = number of years trees survives a wildfire. CI = crown injury from fire. Fire-killed tree data were collected from five fires in California that were monitored for 10 years after death at year 10, respectively. The covariate models accurately predicted individual snag status at year 10 and probability of persistence of 0.5 for 79.2% of incense-cedar, 75.7% of sugar pine, 64% of white fir, and 93.1% of yellow pine (Table 4).

Discussion
Many factors interact to affect wood decay and snag longevity (DeNitto et al. 2000;Lowell et al. 2010;Dunn and Bailey 2015). Snag persistence was species specific in our study, with yellow pine snags falling at a much faster rate than the other species studied. Our results agree with previous findings (Bull 1983;Chambers and Mast 2005;Russell et al. 2006;Dunn and Bailey 2012;Ritchie and Knapp 2014). The differences between species are likely due to a combination of species-specific wood characteristics, xylogphagous insects, and wood decay fungi (Scheffer 1966;Lowell et al. 2010;Dunn and Bailey 2015), although interspecific differences in tree characteristics and injury could play a role (Table 2). Incensecedar was very stable for the first 8 years, as also reported in Ritchie and Knapp (2014) and Ritchie et al. (2013). This is likely attributable to the inherent rot resistance of incense-cedar (Scheffer 1966). Yellow pine is host to many bark and wood-boring beetle species, which can vector fungi and facilitate decay by perforating bark and tunneling into sapwood (Furniss and Carolin 1977). Foraging by woodpeckers can also lead to sapwood cracking and further decay. Hadfield and Magelssen (2006) reported more than three quarters of ponderosa pine fire-killed snag volume was affected by decay within 5 years after death.
Snag age was the strongest predictor of persistence for all species, followed by DBH for yellow pine, sugar pine, and white fir. This concurs with other studies (Keen 1955, Russell et al. 2006, Parish et al. 2010, Dunn and Bailey 2015 Table 1). In those species for which DBH was a significant factor, larger-diameter snags persisted much longer than smaller-diameter snags. This is likely a result of the higher volume of rot-resistant heartwood, the greater percentage of heartwood to sapwood, and the sheer increase in volume in larger snags (Harmon et al. 1986;Lowell et al. 2010). White fir sapwood and heartwood deteriorate quickly after dying, followed by Year to death = number of years trees survives a wildfire. CI = crown injury from fire. Fire-killed tree data were collected from five fires in California that were monitored for 10 years after death. In (b), CI = 75% lines overlap for the plotted values Table 4 Accuracy rates at year 10 and probability of snag persistence cutoff equal to 0.5 for the covariate equations by species. True positive rate is the percentage of snags with a probability of falling greater than 0.5 by year 10 that were observed to fall by year 10. True negative rate is the percentage of snags with a probability of falling less than 0.5 by year 10 that were observed to be standing at year 10. Overall accuracy is the percentage of snags that had their status at year 10 correctly predicted by the covariate model at P = 0.5. ponderosa and Jeffrey pine (Kimmey 1955). Therefore, the tendency of larger snags to persist longer than smaller snags, all else equal, may just be due to more volume and the time it takes for that extra volume to decay. Factors such as aspect, elevation, and how long bark is retained on the snag can also influence snag persistence (Kimmey 1955;Harmon et al. 1986). Crown injury was a significant factor in snag persistence in white fir and yellow pine. Trees with lower crown injury persisted longer as snags than those with higher crown injury in white fir greater than 50 cm DBH and in yellow pine greater than 30 cm DBH. Higher crown injury in smaller white fir and unprotected yellow pine led to increased persistence. Additional research is needed to provide mechanistic explanations of snag decay and persistence patterns. Several studies have shown that snags of some species, such as white fir, may persist longer than other species due to a tendency to break above the downed height threshold and gradually lose height over time (Raphael and White 1984, Everett et al. 1999, Landram et al. 2002, Chambers and Mast 2005 Table 1). This pattern is inconsistent with our results, with average initial break heights below 3.5 m and the majority of snags breaking at groundline. However, Kimmey (1955) also reported that it was common for ponderosa pine and white fir fire-killed snags to break at ground level. We used 3.65 m as the height criterion for a snag to move from standing to downed; therefore, most snags converted to downed upon the initial break and we were not able to detect further trends in snag height loss over time. While most other studies of fire-killed snags use shorter height criteria of between 1.5-1.8 m (Table 1), our findings of species differences and the importance of DBH in determining persistence are still in agreement with these studies.
Choosing a snag persistence model will depend on the intended use and tolerance of accuracy and cost. The FFE-FVS models over-predicted snag longevity for sugar pine, white fir, and yellow pine, and under-predicted snag persistence for incense-cedar at year 10. The large over-prediction of yellow pine snag persistence by FFE-FVS suggests that caution is needed when using this software to estimate snag longevity after fire. The covariate models were more accurate than FFE-FVS, but not as accurate in predicting overall snag persistence as the time-only models. This is due to the increased difficulty in modeling behavior at the individual snag level rather than at the stand level. Only 12 incense-cedars had fallen by the end of the observation period. Although this is enough to develop a reasonable model, in future studies, incense-cedar should be followed for a longer period of time to assess the accuracy of these models past 10 years.
When choosing a basis for predicting snag persistence, it is important to note the limitations of a particular study. We summarized published studies evaluating firekilled snags of the species we assessed for comparison with our results (Table 1). Here, we had right-censored data, only assessed trees that survived at least one year post fire, did not assess unburned trees, and sampled a small number of downed incense-cedar. All studies in Table 1 also contain censored data. Many did their initial snag survey several years post fire, assuming all standing (and in some cases, downed) snags at the time of the measurement were killed in the fire (Dahms 1949;Fig. 5 Observed (solid lines) and FFE-FVS (dashed lines) predicted snag persistence once a tree dies from wildfire (1 = standing; 0 = downed or < 3.65 m tall) by species. Fire-killed tree data were collected from five fires in California that were monitored for 10 years after death Raphael and White 1984;Raphael and Morrison 1987;Morrison and Raphael 1993;Everett et al. 1999;Smith 1999;Chambers and Mast 2005;Russell et al. 2006;Bagne et al. 2008;Dunn and Bailey 2012). Only five studies compared burned and unburned snags (Raphael and White 1984;Raphael and Morrison 1987;Morrison and Raphael 1993;Chambers and Mast 2005;Bagne et al. 2008). Several studies had low sample sizes or inadequate range distributions in one or more subsets of the data (Raphael and White 1984;Everett et al. 1999;Smith 1999;Chambers and Mast 2005;Bagne et al. 2008). Some studies did not follow individual snags, but rather inventoried or censused snags in the study area (Dahms 1949;Raphael and White 1984). Lastly, longitudinal data of natural phenomena are almost certain to have some confounding factors as well as variables that are strongly correlated; these can cause statistical use and interpretation challenges. As with virtually any model application, it is important to choose a model that was developed from data that most closely matches the tree and site conditions in question.

Conclusions and management implications
Snags and downed logs are essential to forest ecosystems, providing numerous ecosystem benefits, but also contributing to fuel loading (Dunn and Bailey 2015). This study provides post-fire snag persistence information for five conifer species to aid land managers in meeting snag-related forest objectives. Snag age and DBH were the strongest predictors of fire-killed snag persistence. For all species, smaller-diameter snags fell more quickly than larger snags, with white fir and sugar pine snags over 50 cm DBH and yellow pine over 90 cm DBH predicted to persist much longer than smaller snags. Incense-cedar snags were fairly stable, but they are not as widely used by wildlife as the other species studied (Thomas 1979;Morrison et al. 1987;Bagne et al. 2008) and therefore may not be the preferred tree species for snag retention for wildlife objectives in a postfire management scenario. However, incense-cedar could be useful for other objectives such as carbon storage. White fir and sugar pine snags both persisted for an extended time on the landscape and are moderately used by wildlife (Thomas 1979;Morrison et al. 1987;Ganey and Vojta 2004;Bagne et al. 2008). Although yellow pine snags are widely used by wildlife, virtually all yellow pine fell within 10 years of dying. In addition, we found FFE-FVS greatly over-predicted yellow pine fire-killed snag persistence and recommend revising the equations in the software to reflect the consistently reported finding of fast fall rates in yellow pine. Snags and logs provide numerous ecosystem benefits (e.g., wildlife habitat, nutrient cycling, germination microsites), but also can potentially increase future fire severity. Our results can be used to inform post-fire snag management considerations.