Fire frequency, area burned, and severity: a ... - Fire Ecology Journal

34 downloads 4695 Views 3MB Size Report
burned from a combination of field- and satellite-derived data. We used the continuous ...... Carbon recovery rates following different wildfire risk miti- .... H.D. Safford, W.B. Cohen, R.S.H. Kennedy, E.E. Knapp, J.K. Agee, and M. Moeur. 2010.
Lutz et al.: Defining Fire Normals Page 51

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051 Research Article

FIRE FREQUENCY, AREA BURNED, AND SEVERITY: A QUANTITATIVE APPROACH TO DEFINING A NORMAL FIRE YEAR James A. Lutz1*, Carl H. Key2, Crystal A. Kolden3, Jonathan T. Kane1, and Jan W. van Wagtendonk4 College of the Environment, Box 352100, University of Washington, Seattle, Washington 98195-2100, USA 1

US Geological Survey, Northern Rocky Mountain Science Center, Glacier Field Station, P.O. Box 128, West Glacier, Montana 59936, USA

2

3

Department of Geography, University of Idaho, P.O. Box 443021, Moscow, Idaho 83844, USA

US Geological Survey (retired), Western Ecological Research Center, Yosemite Field Station, P.O. Box 700, El Portal, California 95318-0700, USA

4

Corresponding author: Tel.: 001-206-616-3827; e-mail: [email protected]

*

ABSTRACT Fire frequency, area burned, and fire severity are important attributes of a fire regime, but few studies have quantified the interrelationships among them in evaluating a fire year. Although area burned is often used to summarize a fire season, burned area may not be well correlated with either the number or ecological effect of fires. Using the Landsat data archive, we examined all 148 wildland fires (prescribed fires and wildfires) >40 ha from 1984 through 2009 for the portion of the Sierra Nevada centered on Yosemite National Park, California, USA. We calculated mean fire frequency and mean annual area burned from a combination of field- and satellite-derived data. We used the continuous probability distribution of the differenced Normalized Burn Ratio (dNBR) values to describe fire severity. For fires >40 ha, fire frequency, annual area burned, and cumulative severity were consistent in only 13 of 26 years (50 %), but all pair-wise comparisons among these fire regime attributes were significant. Borrowing from long-established practice in climate science, we defined “fire normals” to be the 26 year means of fire frequency, annual area burned, and the area under the cumulative probability distribution of dNBR. Fire severity normals were significantly lower when they were aggregated by year compared to aggregation by area. Cumulative severity distributions for each year were best modeled with Weibull functions (all 26 years, r2 ≥ 0.99; P < 0.001). Explicit modeling of the cumulative severity distributions may allow more comprehensive modeling of climate-severity and area-severity relationships. Together, the three metrics of number of fires, size of fires, and severity of fires provide land managers with a more comprehensive summary of a given fire year than any single metric. Keywords: differenced Normalized Burn Ratio, fire severity, fire severity normals, Sierra Nevada, Weibull distribution, Yosemite National Park

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Lutz et al.: Defining Fire Normals Page 52

Citation: Lutz, J.A., C.H. Key, C.A. Kolden, J.T. Kane, and J.W. van Wagtendonk. 2011. Fire frequency, area burned, and severity: a quantitative approach to defining a normal fire year. Fire Ecology 7(2): 51-65. doi: 10.4996/fireecology0702051 INTRODUCTION Fire frequency, fire extent, and fire severity are three of the seven fire regime attributes of particular importance to ecologists and land managers (Sugihara et al. 2006) who closely monitor changes and trends in natural processes. Changes in fire regimes can alter vegetation type (van Wagtendonk and Fites-Kauffman 2006), species composition (Lutz et al. 2009b), forest structure (Peterson et al. 2005), and regeneration patterns (Swanson et al. 2011), which in turn all affect carbon cycling (Hurteau and North 2010), smoke production (Tarnay and Lutz 2011), vertebrate habitat (Roberts et al. 2008), and recreation amenity values (Boxall and Englin 2008). Changing fire regimes also impact efforts to manage (Kolden and Brown 2010) and suppress fire. Projections of increased fire activity associated with climate change (e.g., Littell et al. 2010, Wotton et al. 2010) have increased the need to develop historical fire regime baselines and methods to monitor change. However, few frameworks have been proposed to help land managers quantitatively examine these changes over time or between areas. We propose a set of fire regime metrics—annual fire frequency, annual area burned, and cumulative fire severity distribution—that can be monitored on a purely statistical basis using data available from the Monitoring Trends in Burn Severity (MTBS) database (Eidenshink et al. 2007). The MTBS database provides Landsatbased, 30 m resolution fire severity data and fire perimeters for wildfires in the US from 1984 to present. From the MTBS data, it is possible to quantify the number of large fires (i.e., those above a specified size threshold) for a region of interest, and extract its area burned and severity. The quantity of fires above a certain size threshold (likely to vary based on the

specifics of the ecosystem) is important to land managers dealing with logistics and costs of fire and air quality monitoring and, if necessary, suppression. The quantity of fires may also give some indication of the distribution of burned areas within the study region. Area burned is the most commonly monitored wildfire metric and represents that portion of the landscape that has been burned by fire. Area burned is most often quantified by the mapped fire perimeter area rather than the actual area burned within the fire perimeter (e.g., Westerling et al. 2006, Morgan et al. 2008), which may overestimate the actual burned area (Kolden and Weisberg 2007). Using satellite measurements of burned area and severity is significantly more accurate than using mapped perimeters (Kolden and Weisberg 2007), especially if fire perimeters were approximated by aerial observers. Fire severity approximates the immediate ecological effect of a fire on vegetation and soil. It is distinct from fire intensity, which is a measure of the energy released by the fire (in kW m-1) (van Wagtendonk 2006). Correlations between fire intensity and fire severity vary depending on specifics of vegetation, fuel loading, and fire behavior (Sugihara et al. 2006). Fire severity gives an indication of post-fire conditions within the area burned, regardless of whether that area is large or small (Lentile et al. 2006). Fire severity, although important ecologically, has been more difficult to quantify and standardize due to varying objectives and perceptions, and most analyses have either ignored it or considered broad severity classifications (i.e., low versus high severity) (see Lentile et al. 2006 for examples in different ecosystems). Area burned can be considered for an individual fire or for an entire fire year (the superimposition of all fires for that year). Within a burned area, some portions burn at compara-

Lutz et al.: Defining Fire Normals Page 53

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Reference severity Lower severity Higher severity

1.00

Cumulative distribution

tively higher and some at comparatively lower severities. Fire severity can be measured through ground observations at localized sites (Key 2006, Key and Benson 2006), but is more commonly measured at broader spatial scales through satellite-derived proxy indices of fire severity, most commonly, the differenced Normalized Burn Ratio, dNBR (Key and Benson 2006). Using data acquired by Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper-plus (ETM+), dNBR is calculated from the near-infrared (Landsat band 4) and the mid-infrared (Landsat band 7), providing 30 m spatial resolution (one pixel = 0.09 ha). The dNBR compares the difference between pre-fire and post-fire Landsat scenes (dNBR = NBRPRE – NBRPOST); where NBR = ([B4 – B7] ÷ [B4 + B7]), and where B4 and B7 represent reflectance values). The dNBR can range between –2.0 and 2.0. Over natural landscapes, however, non-anomalous dNBR values typically have a more limited range of about –0.6 to 1.2., which is usually scaled by a factor of 1000. High values of dNBR represent a combination of a decrease in the reflectance of Landsat B4, indicating primarily a decrease in photosynthetic materials, and an increase in the reflectance of Landsat B7, indicating an increase in ash, carbon, and soil, as well as a decrease in surface materials holding water. In most burned areas, the dNBR value of all pixels follows a smooth cumulative probability distribution (Figure 1). Fire severity can also be derived from other satellite indexes such as the relative differenced Normalized Burn Ratio (RdNBR) (Miller and Thode 2007). The RdNBR provides a means of nomalizing satellite-derived fire severity values across vegetation types or between stands of different ages or productivities. The RdNBR value of all pixels also follows a smooth cumulative probability distribution (Thode et al. 2011). Satellite-derived dNBR is most commonly stratified into four burn severity categories (high, moderate, and low severities; and no detected change), with subsequent analysis per-

0.75

0.50

0.25

0.00 -200

0

200

400

600

800 1000 1200

dNBR

Figure 1. Cumulative distributions of fire severity for a reference area, and for areas that burned with lower and higher severities. Compared to the curve of reference severity (black), the lower severity curve (green) has few high severity pixels. Conversely, the higher severity curve (red) shows an area that burned with more higher severity pixels. These model curves are of the ͳ form: ‫ ݕ‬ൌ . ͳ ൅ ݁ ିሺௗே஻ோି௔ሻ

formed on the four classification levels or upon the simple area burned (Miller and Thode 2007, Miller et al. 2008, Lutz et al. 2009a). Demarcations between severity levels are ideally determined after extensive ground verification of the effect of fire on each vegetation type (e.g., Thode et al. 2011). Classifying severity data without detailed ground information has the potential to distort the interpretation of fire severity, particularly in the highest and lowest classifications, or where vegetation types or stand ages are closely intermixed. Furthermore, classification of any continuous variable into discrete bins sacrifices information and, thereby, increases noise in subsequent analyses. Classification is also sensitive to potentially differing interpretations made by remote sensing analysts and field ecologists (Spies et al. 2010), suggesting that a statistical

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Lutz et al.: Defining Fire Normals Page 54

approach would be more robust, at least for sufficiently large areas (Hudak et al. 2007). Evaluation of the fire severity distribution itself (Figure 1) can serve as a quantification of the severity of an entire burned area, whether that area represents a single fire or all the fires in a given fire year. The area under the cumulative severity distribution curve decreases as the number of pixels burned at relatively higher severity (i.e., have higher dNBR values) increases (Figure 1). A quantitative metric that both represents the overall severity of the area burned and that increases with severity is one minus the area under the curve, with higher values of this metric indicating a higher cumulative severity distribution. The discrete nature of dNBR values and the variability inherent in real world data suggest that data be aggregated between reasonable possible values of dNBR, limits that will vary between ecoregions. METHODS Study Area

Yosemite National Park is a contiguous management unit of 3027 km2 located in the central Sierra Nevada, California, USA. Yosemite experiences multiple wildland fires (prescribed fires and wildfires) each year and, since 1972, many naturally ignited fires have been allowed to burn under prescribed conditions (van Wagtendonk 2007). In adjacent lands managed by the Forest Service, most fires are suppressed. Our study area (hereafter referred to as Yosemite) included the area of the park and a buffer radiating 6.5 km from the park boundaries (4771 km2) (Figure 2). The buffer was selected to maximize the size of the study area while minimizing non-forest cover and developed areas. The study area is in federal ownership (Yosemite National Park, Stanislaus National Forest, Toiyabe National Forest, Inyo National Forest, and Sierra National Forest) with the exception of small inholdings. Yosemite’s climate is Mediterranean. July

Figure 2. Land management and forest cover within the Yosemite study area. The 4771 km2 study area (green outline) comprises Yosemite National Park and portions of the Stanislaus, Toiyabe, Inyo, and Sierra national forests—areas with low levels of anthropogenic disturbance. Perimeters for fires >40 ha between 1984 and 2009 are shown in red. Landsat 5 Thematic Mapper scene from 16 August 2010 (bands 3, 2, 1).

mean minimum and maximum temperatures are 2 °C to 13 °C at higher elevations and 16 °C to 35 °C at lower elevations. Annual precipitation ranges from 800 mm to 1720 mm, with most precipitation falling in the winter as snow (Lutz et al. 2010). The forest vegetation of Yosemite comprises a mosaic of forest types, species, and structural stages (van Wagtendonk and FitesKaufman 2006, Fites-Kaufman et al. 2007). The lower elevation portions of the park include shrub patches, shrub fields, and woodland as well as forests. The Sierra Nevada fire regime is mixed; fires burn with patches of high, moderate, and low severities at intervals ranging from years to centuries (Agee 1993,

Lutz et al.: Defining Fire Normals Page 55

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

van Wagtendonk et al. 2002, Sugihara et al. 2006a, van Wagtendonk and Fites-Kaufman 2006, van Wagtendonk and Lutz 2007). Fire intensity and severity vary by forest type (van Wagtendonk et al. 2002, Thode et al. 2011). The natural fire return interval for the forested ecosystems of Yosemite National Park ranges from 4 yr to 187 yr (Caprio and Swetnam 1995, Caprio and Lineback 1997, van Wagtendonk et al. 2002, Collins and Stephens 2007). Fire Severity Data

Satellite fire severity data (dNBR) were obtained from MTBS. We first obtained fire perimeters from Yosemite National Park (National Park Service 2010) and from the Forest Service (USDA 2010). Using the area burned from the fire perimeter data, we determined that 84 % of area burned is accounted for by fires ≥400 ha (the standard minimum fire size for MTBS in the western US), while 97 % of the area burned is accounted for by fires ≥40 ha. We therefore set 40 ha as the minimum fire size for analysis. We examined all Landsat scenes with 150, we extended fire perimeters to include those pixels.

We tabulated area burned in multiples of the 0.09 ha Landsat pixels. Data Reduction and Analysis

We constructed cumulative distributions of burn severity using all the Landsat pixels within the fire perimeters for each year between 1984 and 2009. We then compared each year’s cumulative severity distribution with an average severity distribution for all years. We calculated the average severity distributions in two different ways: one that gave more weight to years with greater area burned, and one that gave equal weight to each year in the study period. To determine the overall cumulative severity distribution of all fires, we aggregated all burned pixels for the entire study period (areal averaging). To determine an annual normal distribution, we aggregated all fires >40 ha in each year, and then combined them with equal weight (yearly averaging). Preliminary analysis showed that >99 % of dNBR pixel values fell in a dNBR range of –200 to 1200. Therefore, we constructed severity distributions using only pixels with dNBR values between –200 and 1200, using the actual dNBR values (no binning). We approximated an integration of these cumulative severity distributions by calculating the area under each distribution for each year to generate a single number representative of the shape of the severity distribution. We represented our severity metric (SM) as one minus the area under the curve so that higher numbers reflect higher severity. We calculated: ଵଶ଴଴

ܵ‫ ܯ‬ൌ ͳ െ  ෍

௜ୀିଶ଴଴

ܲ‫ ݏ݈݁ݔ݅݌ ݂݋ ݊݋݅ݐݎ݋݌݋ݎ‬൑ ݅ ͳͶͲͳ

(1)

In addition to the yearly and areal averages, we also used the dNBR range of –200 to 1200 to fit sigmoid (Equation 2), logistic (Equation 3), and Weibull (Equation 4) functions (Weibull 1951) to each year’s burned area. Because the Weibull distribution is sin-

Lutz et al.: Defining Fire Normals Page 56

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

gle-sided (having a zero value at some lower limit [Equation 4]), we hypothesized that it might better model the transition between unburned and burned. Data were fit to the equations with SigmaPlot Version 11.0 (Systat Software, Chicago, Illinois, USA) with x0, b, and c being the shape parameters of the equations. We bounded the data range to less than the maximum range of dNBR so as to avoid overstating the goodness of fit.

‫ݕ‬ൌ

ͳ൅

‫ ݕ‬ൌ

ͳ

௫ି௫బ  ିቀ ቁ ௕ ݁

ͳ

‫ ݔ‬௕ ͳ൅ቀ ቁ ‫ݔ‬଴



(2)

(3)

Table 1. Annual fire statistics for fires >40 ha in the Yosemite study area, 1984 to 2009, and comparisons among fire regime metrics. Area burned represents the area of Landsat pixels (0.09 ha) either within the reported fire perimeter or outside and adjacent to the reported perimeter with dNBR ≥ 150. The frequency, area burned, and severity of fires are classified by quartile (L = lowest quartile, M = middle 50 %, H = upper quartile). The years with the most extensive annual area burned are associated with high cumulative severity, and the years with low annual area burned are associated with low cumulative fire severity. However, the years with an annual area burned in the middle 50 % are associated with all levels of cumulative burn severity. The three metrics of frequency, area, and severity were consistent in 13 of 26 years (bold). Comparison with Area of fires Severity§ middle 50 % † years >40 ha No. Area Severity (ha) fires (ha) (pixel) 795 0.20 M L L 3 395 0.24 H M M 1 789 0.25 M M M 12 925 0.35 H H H 6 550 0.25 H H M 698 0.17 L L L 10 909 0.36 H H H 2 718 0.27 M M M 538 0.21 M L L 488 0.20 L L L 2 412 0.21 M M L 407 0.25 L L M 24 812 0.34 M H H 1 656 0.26 L M M 1 008 0.21 L L L 7 257 0.25 H H M 144 0.16 L L L 3 018 0.28 L M H 2 348 0.22 M M M 5 253 0.31 H H H 3 846 0.29 M M H 1 935 0.24 H M M 3 124 0.25 M M M 1 228 0.21 M M M 2 933 0.28 M M M 5 546 0.29 M H H † L = annual value lower than the middle 50 % of years, M = within the middle 50 % of years, H = annual value higher than the middle 50 % of years.

No. fires Year >40 భ ೎ ha ೣషೣబ శ್೗೙మ೎ ିቌ ቍ భ ್ 1984 6 భ ೎ ‫ ݕ‬ൌͳെ݁ ;‫ ݔ‬൐  ‫ݔ‬଴ (4) െ ܾ ൈ Žሺʹሻ  1985 ೎11 ೣషೣబ శ್೗೙మ೎ ିቌ ቍ 1986 4 భ ್ 1987 10 ೎ ൌͳെ݁ ;‫ݔ‬  ൐  ‫ݔ‬଴ െ ܾ ൈ Žሺʹሻ  1988 15 1989 1 1990 8 For each value of fire frequency, area, and 1991 4 severity, we classified each year into quartiles, 1992 4 considering the first quartile as low, the middle 1993 1 1994 4 two quartiles as normal, and the upper quartile 1995 3 as high. We compared our purely statistical 1996 6 division of dNBR fire severities with the clas1997 3 1998 3 sified fire atlas based on RdNBR from van 1999 11 Wagtendonk and Lutz (2007). We limited the 2000 1 comparison to fires covered by both studies. 2001 1 2002 5 The comparison therefore was limited to the 2003 10 years 1984 through 2005, to fires within or 2004 4 crossing the Yosemite National Park boundar2005 10 2006 7 ies, and to the 6.5 km buffer around the park. 2007 6 2008 6 RESULTS 2009 4

From 1984 through 2009, there were 148 fires >40 ha in the Yosemite study area. The average area burned per year in fires >40 ha was 4144 ha (Table 1). Cumulative distributions of dNBR severity values between –200 and 1200 were best modeled by Weibull func-

§

Severity metric (SM) calculated from the cumulative distribution of severity for each year as: ଵଶ଴଴

ܵ‫ ܯ‬ൌ ͳ െ  ෍

௜ୀିଶ଴଴

ܲ‫ ݏ݈݁ݔ݅݌ ݂݋ ݊݋݅ݐݎ݋݌݋ݎ‬൑ ݅ ͳͶͲͳ

ሺͳሻ

Lutz et al.: Defining Fire Normals Page 57

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

tion (Equation 4). All years, as well as areal and yearly aggregations, could be described by this function, with every year and the areal and yearly averages having an r2 ≥ 0.993 and P < 0.001. The logistic function (Equation 3) was not a good fit when x0 approached zero (in 1992 and 1994). The sigmoid functions had lower adjusted r2 than the Weibull function in every year. Cumulative distribution of fire severity varied among years (Figure 3, Table 2). The average cumulative distribution of severity for all years differed significantly (KS-test, P < 0.001) whether the average was calculated on a yearly or areal basis, confirming the potential

1.0

differences between annual and areal averaging. On an areal basis, the mean dNBR severity in Yosemite was 158, with lower and upper quartiles at 65 and 313 (Figure 4). On a yearly basis, the mean dNBR severity in Yosemite was 102, with lower and upper quartiles at 33 and 214. Calculating SM between dNBR values of –200 and 1200 yielded values of SM from a minimum of 0.16 in 2000 to a maximum of 0.36 in 1990. The areal SM was 0.30 and the yearly SM was 0.25. Years with low severity distributions were characterized by lower area burned, and years with high severity distributions were characterized by greater area burned (Figure 3, Table 1). However,

1989 2000

Cumulative distribution

0.8 1987 1990 1996

0.6

0.4

0.2 By years By pixels

0.0 -200

0

200

400

600

800

1000

1200

dNBR Figure 3. Cumulative fire severity for all fires >40 ha that burned in each year between 1984 and 2009 (gray lines). Each gray line represents the cumulative severity distribution for one fire year. The annual average cumulative severity (each fire year weighted equally) is shown in red. The cumulative severity for the entire area burned by all fires >40 ha in the Yosemite study area between 1984 and 2009 is shown in orange. The two years with the lowest cumulative severity (2000, SM = 0.16; and 1989, SM = 0.17) are at left; the three fire years with the highest cumulative severity (1987, SM = 0.35; 1990, SM = 0.36; and 1996, SM = 0.34) are at right. The year with the least burned area (2000) also had the lowest cumulative severity, but the year with the most area burned (1996) was third in cumulative severity. The years 1987 and 1990 had higher cumulative severities reflecting the lower elevation and more severe fire regimes of the vegetation (mixed forest and shrubs).

Lutz et al.: Defining Fire Normals Page 58

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Table 2. Weibull parameters for area burned by fires >40 ha in Yosemite between 1984 and 2009. The data illustrated in Figure 3 were fit to a curve భ ೎ of the form: ೣషೣబ శ್೗೙మ೎ All P values < 0.001. 





By pixels By years



;‫ ݔ‬൐  ‫ݔ‬଴ െ ܾ ൈ Žሺʹሻ೎ 

Year b c x0 1984 118.7810 1.0810 54.2297 1985 146.1845 0.8969 75.6685 1986 153.8405 0.9766 92.2235 1987 387.9411 1.4239 243.9727 1988 200.7272 1.0391 97.7872 1989 67.3086 1.0719 9.5990 1990 416.7735 1.4848 253.9185 1991 227.3193 1.1967 135.0312 1992 123.2930 0.9445 52.3423 1993 123.8919 1.6866 74.1267 1994 168.1034 1.2922 59.4338 1995 193.2867 1.0225 94.6221 1996 291.1345 1.1743 212.7715 1997 211.6579 1.3631 131.3029 1998 138.9025 1.6552 75.9263 1999 220.8271 1.3292 121.1346 2000 83.7842 1.5352 1.5352 2001 186.6766 1.0194 140.9227 2002 137.7264 1.4858 90.4649 2003 282.6585 1.3200 187.2860 2004 293.9449 1.2722 153.6713 2005 161.8939 1.3377 111.8003 2006 205.0641 1.6230 123.5403 2007 123.0021 1.2231 75.7596 2008 238.3653 1.5045 160.7112 2009 305.4238 1.4179 165.9796 All pixels† 262.1238 1.1511 160.3395 All years‡ 185.5725 1.0499 101.7855 † All burned pixels from all years aggregated. ‡ All burned pixels within one year aggregated, and then years averaged.

Proportion of pixles

‫ ݕ‬ൌͳെ݁

ିቌ

-200

0

200 400 600 800 1000 1200

dNBR

Figure 4. Distribution of fire severity for the Yosemite study area. All dNBR pixels from fires >40 ha from 1984 through 2009 are included. Vertical lines delineate the severity quartiles for all area burned between 1984 and 2009.

those fire regime attributes relating frequency, area, and severity were only consistent in 13 of 26 years (Table 1). Although the three metrics were not consistently related, the individual pair-wise comparisons among variables were significant (Figure 5). Greater numbers of fires were associated with larger annual area burned (r2 = 0.34, P = 0.001) and with higher

Figure 5. Relationships among number of fires >40 ha, area burned by fires >40 ha, and severity distributions of fires >40 ha in the Yosemite study area for each year between 1984 and 2009. Points shown in red indicate years with low snowpack (sensu Lutz et al. 2009a) and points shown in blue represent years with high snowpack.

Lutz et al.: Defining Fire Normals Page 59

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

annual burn severity as measured by SM (r2 = 0.12, P = 0.04). The strongest relationship was between annual area burned and SM (r2 = 0.69, P < 0.001). Between 1984 and 2009, there was no trend in the frequency, area, or severity of fires >40 ha (all r2 < 0.01, all P > 0.5) in the Yosemite study area (Figure 6). The purely statistical approach to severity classification (quartiles) was consistent with the classified RdNBR values for the fires in common between the MTBS and previous fire atlas. In comparing the two fire atlases, the MTBS processing of fire severity between 1984 and 2005 yielded 92 374 ha of burned area. The previous Yosemite fire atlas yielded 94 405 ha. The simple quartile numbers understated severity (Table 3). However, the correspondence between classes was high, and the

spatial structure of severity patches was consistent between the two approaches (Figure 7). DISCUSSION Fire Severity and Annual Area Burned

The years with the most extensive annual area burned are associated with high cumulative severity, and the years with low annual area burned are associated with low cumulative fire severity (Table 1). However, the association often depends on the size of large fires and the vegetation type in which the fire burns. The lower elevation portions of the study area with more shrub cover burn at higher aggregate severity (Thode et al. 2011). Fires originating in these lower elevation forests

Figure 6. There is no temporal trend in number of fires >40 ha, area burned by fires >40 ha, and severity distributions of fires >40 ha in the Yosemite study area between the years of 1984 and 2009. Points shown in red indicate years with low snowpack (sensu Lutz et al. 2009a) and points shown in blue represent years with high snowpack. The scale of the study area (4771 km2) is likely too small to discern statistically significant temporal trends. Table 3. Comparison between severity quartiles and classified severity between 1984 and 2005. The classified severity (from van Wagtendonk and Lutz 2007) uses a fire atlas originally developed by Thode (Miller and Thode 2007, Thode et al. 2011). The present study used the MTBS fire atlas. The comparison is for the burned areas of fires covered by both studies. Outside perimeter Outside perimeter Unchanged Low severity Moderate severity High severity

0.0 307.5 224.9 81.5 12.7

Lowest Low-mid High-mid quartile quartile quartile ····································ha···································· 225.3 111.2 62.7 728.9 287.6 135.7 887.6 962.9 617.5 388.4 781.0 1195.3 68.9 109.6 279.4

Highest quartile 24.4 38.9 25.8 916.0 1228.4

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Lutz et al.: Defining Fire Normals Page 60

severity, presumably because conditions not conducive to fire spread are also not conducive to more complete vegetation consumption. The two years with the lowest SM (1989 and 2000) featured only one fire >40 ha (Table 1, Figure 3). Furthermore, these fires were management-ignited prescribed burns, which usually have lower severities (van Wagtendonk and Lutz 2007). Although we did not detect temporal trends in fire frequency, area, or severity, we cannot discount that such trends may exist. Rather, the high interannual variability of fires in the immediate Yosemite area (Table 1, Figures 5 and 6), partially driven by spring snowpack, masks any gradual trend (Lutz et al. 2009a). Larger spatial scales (e.g., Miller et al. 2008) or a longer data record would be necessary to confirm fire regime trends. Figure 7. Fires in the Yosemite study area >40 ha between 1984 and 2009 stratified by dNBR quartiles. Quartiles were determined from all burned area over the 26 yr study period. Where areas were reburned, the severity of the most recent fire is shown.

have often been large (1987, 1990, and 1996), and the effect of vegetation type on fire severity distribution seems more important than the effect of fire area (details in Thode et al. 2011). The year with the greatest area burned (1996) had SM = 0.35. However, fire years 1987 and 1990, which had approximately half the area burned of 1996, featured equivalent or higher severity (SM = 0.38 and 0.35). These two years had very low levels of spring snowpack. All three of these fires burned in lower mixed conifer forest types, but the 1996 fire year featured higher area burned in closed canopy forest, and the 1987 and 1990 fire years burned in areas vegetated with mixed forest, woodland, and chaparral. Lower quantities of area burned were more closely associated with lower cumulative burn

Modeling Continuous Fire Severity Distributions

We found the cumulative severity distributions to be very amenable to approximation by continuous functions (Equation 4, Table 2). Although these equations had r2 > 0.99, the residuals were not normally distributed (as would be expected from a cumulative distribution). Therefore, Equation 4 may not model what may be important secondary fire regime attributes, nor may it model ecosystems with multimodal cumulative severity distributions. Our objective in seeking a continuous equation was to represent severity distributions in a way that could be easily manipulated with parameters that could be correlated with climate and vegetation. Although Equation 4 represents severity well, all the Weibull shape parameters (Table 2) are needed to characterize the severity distributions (Figure 3). The single-number approximation from Equation 1 may be as useful, and because it is easier to calculate (see Appendix), may be a convenient and useful tool for land managers.

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Local Land Management through Time

We borrowed the “fire normals”—the longterm average values of multiple fire parameters from a geographic area—approach from climatology, where characterizations of climate conditions are most frequently framed as departures from a multidecadal mean, or normal. Our approach in delineating fire severity distributions into quartiles responds to a statistic commonly used in a number of disciplines. Rather than quartiles, it could just as effectively be based on some other percentile if there is ecological justification to do so, but whatever the case, the approach remains well defined and repeatable. There are several important advantages and limitations. The principal advantage is that it can be done automatically, without need to analyze the specific vegetation related correlates of fire severity; the approach is entirely statistical and does not rely on expert judgment. A purely quantitative approach avoids any possibility of bias in the selection of burn severity thresholds. As fire effects information becomes increasingly integrated by third parties into ecological research, resource management, and policy, it will be ever more important to ensure that severity data are not misused by advocates for particular policies. However, this statistical approach is only relevant within a particular study area. Accurate interpretation of satellite-derived fire severity requires calibration with local field data. The results permit a comparison through time, but there can be little, if any, extrapolation to other geographic areas (Table 3, Figure 7). On the other hand, such aspects of fire regimes in one region can be compared to other regions with some quantifiable and meaningful interpretation of the differences or similarities. A fire normals approach also requires a sufficient period of time so that the fire regime can be adequately represented. The climatological community commonly uses 30 yr periods to describe climate normals for a region. The Landsat-derived dNBR record of 26 yr, used

Lutz et al.: Defining Fire Normals Page 61

here, approximates that period. A fire normals approach also requires judgment in selecting an area and a set of fires that is sufficiently characteristic of the area for the summary statistics to be of immediate use to land managers. Provided that these conditions are met, land managers tasked with understanding local changes can use these fire normals to compare individual fires and fire years to the historical record. In this study area, the purely statistical approach gives equivalent results to fire-byfire classification (Table 3, Figure 6). The lack of this fire-by-fire analysis of fire effects also constitutes the greatest weakness to this approach. Without fire-by-fire analysis of pre-fire vegetation and fire effects, inferences from the raw dNBR values could be incorrect. There is no inherently ecological reason for considering a certain dNBR population percentile to represent a particular level of severity. Inferring ecological effects of satellitederived severities depends on vegetation type (Thode et al. 2011), time since previous fire (Miller and Thode 2007, Thompson et al. 2007, Larson et al. 2008), and the nature of post-fire regeneration (Key 2006; Kane et al. 2008, 2010). In addition, while dNBR is well accepted as a measure of fire severity, in the Sierra Nevada, the relativized version of dNBR (RdNBR) has been shown to be robust (Miller and Thode 2007, but see Soverel et al. 2010). This statistical approach could be applied to RdNBR values as well. Choosing a Severity Metric

The severity metric (Equation 1) and the continuous approximation (Equation 4) provide additional methods for quantifying severity. Severity is most often quantified by classifying into four, or sometimes six, levels of severity. The severity of entire burned areas has previously been quantified by averages (i. e., Roberts et al. 2008) or proportions of high severity area (Lutz et al. 2009b), but SM represents a continuum of severities and distribu-

Lutz et al.: Defining Fire Normals Page 62

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

tions. No matter what characteristic severity distribution exists for a particular management unit, questions related to a fire year (i.e., “Was this year an active fire year?”) can be examined by comparing the SM for the year in question to the historical average (Figure 3), along with the comparisons of number of fires and annual area burned. The Weibull shape parameters (Table 2), as well as the SM, maintain the continuous distribution of dNBR, and they may also reveal relationships between fire severity and abiotic predictors such as climate conditions, elevation, or vegetation type that are not evident in analyses of classified fire severity. When examining questions related to an area burned, either a single fire or a fire year, (i.e., “Did this fire burn more severely than normal for the area?”), a better comparison is with the severity mean calculated from all pixels over the period of record. An areal metric decreases the influence of those years where there was little burned area. Programs such as MTBS help managers track the dynamics of the landscapes they are responsible for stewarding by using consistent

processing methods. The MTBS methodology (usually limited to fires >400 ha in the western US and >200 ha in the eastern US) can be extended to smaller fires so that essentially all of the burned area (97 % in this work) can be analyzed. Using the Landsat Thematic Mapper period of record, the historical fire severity can be considered in purely statistical terms, allowing comparisons of severity between years. The current no-cost availability of the entire Landsat archive allows many scenes to be examined for the highest quality pre-fire and post-fire scene pair, and the uniform processing procedures of MTBS provide a standardized set of information. The normal fire regime attributes of frequency, area, and severity can all be easily calculated from the data provided by the MTBS program so that land managers can easily calculate the fire normals specific to their location (see Appendix for a calculation example from MTBS data). These fire normals provide a framework for examining how fires or fire years differ from recent means, especially as fire regimes may be affected by future climate change.

ACKNOWLEDGMENTS We thank Stephen Howard for assistance with MTBS processing, Jay Miller for discussions about Sierra Nevada fire history and data, and Kent van Wagtendonk and Yosemite National Park for data and data management. Funding was provided by the US Geological Survey Global Change Research Program (Climate change impacts on burn severity in three forest ecoregions of the US). Previous versions of this manuscript were improved by comments from Nate Benson, Kim Keating, Janine Powell, and two anonymous reviewers. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the US government. LITERATURE CITED Boxall, P.C., and J.E. Englin. 2008. Fire and recreation values in fire-prone forests: exploring an intertemporal amenity function using pooled RP-SP data. Journal of Agricultural and Resource Economics 33(1): 19-33. . Accessed 22 February 2011. Caprio, A., and P. Lineback. 1997. Pre-twentieth century fire history of Sequoia and Kings Canyon national parks: a review and evaluation of our knowledge. Pages 180-199 in: N.G. Sugihara, M.A. Morales, and T.J Morales, editors. Proceedings of the conference on fire in California ecosystems: integrating ecology, prevention, and management. Association for Fire Ecology Miscellaneous Publication 1. Sacramento, California, USA

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Lutz et al.: Defining Fire Normals Page 63

Caprio, A.C., and T.W. Swetnam. 1995. Historic fire regimes along an elevational gradient on the west slope of the Sierra Nevada, California. Pages 173-179 in: J.K. Brown, R.W. Mutch, C.W. Spoon, and R.H. Wakimoto, technical coordinators. Proceedings of the Symposium on Fire in Wilderness and Park Management. USDA Forest Service General Technical Report INT-GTR-320. Intermountain Research Station, Ogden, Utah, USA. Collins, J.B., and S.L. Stephens. 2007. Fire scarring patterns in Sierra Nevada wilderness areas burned by multiple wildland fire use fires. Fire Ecology 3(2): 53-67. doi: 10.4996/ fireecology.0302053 Eidenshink, J., B. Schwind, K. Brewer, Z. Zhu, B. Quayle, and S. Howard. 2007. A project for monitoring trends in burn severity. Fire Ecology 3(1): 3-21. doi: 10.4996/fireecology.0301003 Fites-Kaufman, J., P. Rundel, N. Stephenson, and D.A. Weixelman. 2007. Montane and subalpine vegetation of the Sierra Nevada and Cascade Ranges. Pages 456-501 in: M. Barbour, T. Keeler-Wolf, and A.A. Schoenherr, editors. Terrestrial vegetation of California. University of California Press, Berkeley, USA. Hudak, A.T., P. Morgan, M. Bobbitt, and L. Lentile. 2007. Characterizing stand-replacing harvest and fire disturbance patches in a forested landscape: a case study from Cooney Ridge, Montana. Pages 209-231 in: M.A. Wulder and S.E. Franklin, editors. Understanding forest disturbance and spatial patterns: remote sensing and GIS approaches. Taylor and Francis, London, United Kingdom. Hurteau, M.D., and M. North. 2010. Carbon recovery rates following different wildfire risk mitigation treatments. Forest Ecology and Management 260: 930-937. doi: 10.1016/j. foreco.2010.06.015 Kane, V.R., A.R. Gillespie, R. McGaughey, J.A. Lutz, K. Ceder, and J.F. Franklin. 2008. Interpretation and topographic correction of conifer forest canopy self-shadowing using spectral mixture analysis. Remote Sensing of Environment 112: 3820-3832. doi: 10.1016/j. rse.2008.06.001 Kane, V.R., R. McGaughey, J.D. Bakker, R. Gersonde, J.A. Lutz, J.F. Franklin. 2010. Comparisons between field- and LiDAR-based measures of stand structural complexity. Canadian Journal of Forest Research 40: 761-773. doi: 10.1139/X10-024 Kane, V.R., J.D. Bakker, R.J. McGaughey, J.A. Lutz, R. Gersonde, and J.F. Franklin. 2010. Examining conifer canopy structural complexity across forest ages and zones with LiDAR data. Canadian Journal of Forest Research 40: 774-787. doi: 10.1139/X10-064 Key, C.H. 2006. Ecological and sampling constraints on defining landscape fire severity. Fire Ecology 2(2): 178-203. doi: 10.4996/fireecology.0202178 Key, C.H., and N.C. Benson. 2006. Landscape assessment: ground measure of severity, the Composite Burn Index, and remote sensing of severity, the Normalized Burn Ratio. Pages LA1-LA55 in FIREMON: fire effects monitoring and inventory system. D.C. Lutes, R.E. Keane, J.F. Caratti, C.H. Key, N.C. Benson, S. Sutherland, and L.J. Gangi. USDA Forest Service General Technical Report RMRS-GTR-164CD. Rocky Mountain Research Station, Fort Collins, Colorado, USA. Kolden, C.A., and T.J. Brown. 2010. Beyond wildfire: perspectives of climate, managed fire and policy in the USA. International Journal of Wildland Fire 19: 364-373. doi: 10.1071/ WF08111 Kolden, C.A., and P.W. Weisberg. 2007. Assessing accuracy of manually mapped wildfire perimeters in topographically dissected areas. Fire Ecology 3(1): 22-31. doi: 10.4996/ fireecology.0301022

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Lutz et al.: Defining Fire Normals Page 64

Larson, A.J., J.A. Lutz, R.F. Gersonde, J.F. Franklin, and F.F. Hietpas. 2008. Potential site productivity influences the rate of forest structural development. Ecological Applications 18: 899-910. doi: 10.1890/07-1191.1 Lentile, L.B., Z.A. Holden, A.M.S. Smith, M.J. Falkowski, A.T. Hudak, P. Morgan, S.A. Lewis, P.E. Gessler, and N.C. Benson. 2006. Remote sensing techniques to assess active fire characteristics and post-fire effects. International Journal of Wildland Fire 15: 319-345. doi: 10.1071/WF05097 Littell, J.S., E.E. Oneil, D. McKenzie, J.A. Hicke, J.A. Lutz, R.A. Norheim, and M.M. Elsner. 2010. Forest ecosystems, disturbance, and climatic change in Washington state, USA. Climatic Change 102: 129-158. doi: 10.1007/s10584-010-9858-x Lutz, J.A., J.W. van Wagtendonk, and J.F. Franklin. 2010. Climatic water deficit, tree species ranges, and climate change in Yosemite National Park. Journal of Biogeography 37: 936-950. doi: 10.1111/j.1365-2699.2009.02268.x Lutz, J.A., J.W. van Wagtendonk, A.E. Thode, J.D. Miller, and J.F. Franklin. 2009a. Climate, lightning ignitions, and fire severity in Yosemite National Park, California, USA. International Journal of Wildland Fire 18: 765-774. doi: 10.1071/WF08117 Lutz, J.A., J.W. van Wagtendonk, and J.F. Franklin. 2009b. Twentieth-century decline of largediameter trees in Yosemite National Park, California, USA. Forest Ecology and Management 257: 2296-2307. doi: 10.1016/j.foreco.2009.03.009 Lutz, J.A., and C.B. Halpern. 2006. Tree mortality during early forest development: a long-term study of rates, causes, and consequences. Ecological Monographs 76: 257-275. doi: 10.1890/0012-9615(2006)076[0257:TMDEFD]2.0.CO;2 Miller, J.D., and A.E. Thode. 2007. Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR). Remote Sensing of Environment 109: 66-80. doi: 10.1016/j.rse.2006.12.006 Miller J.D., H.D. Safford, M. Crimmins, and A.E. Thode. 2008. Quantitative evidence for increasing forest fire severity in the Sierra Nevada and southern Cascade Mountains, California and Nevada, USA. Ecosystems 12: 16-32. doi: 10.1007/s10021-008-9201-9 Morgan, P., E.K. Heyerdahl, and C.E. Gibson. 2008. Multi-season climate synchronized forest fires throughout the 20th century, northern Rockies, USA. Ecology 89: 717-728. doi: 10.1890/06-2049.1 National Park Service. 2010. Fire history for Yosemite National Park 1930-2009. . Accessed 1 May 2010. Peterson, D.L., M.C. Johnson, J.K. Agee, T.B. Jain, D. McKenzie, and E.D. Reinhardt. 2005. Forest structure and fire hazard in dry forests of the western United States. USDA Forest Service General Technical Report PNW-GTR-628. Pacific Northwest Research Station, Portland, Oregon, USA. Roberts, S.L., J.W. van Wagtendonk, D.A. Kelt, A.K. Miles, and J.A. Lutz. 2008. Modeling the effects of fire severity and spatial complexity on small mammals in Yosemite National Park, California. Fire Ecology 4(2): 83-104. doi: 10.4996/fireecology.0402083 Soverel, N.O., D.D.B Perrakis, and N.C. Coops. 2010. Estimating burn severity from Landsat dNBR and RdNBR indices across western Canada. Remote Sensing of Environment 114: 1896-1909. doi: 10.1016/j.rse.2010.03.013 Spies, T.A., J.D. Miller, J.B. Buchanan, J.F. Lehmkuhl, J.F. Franklin, S.P. Healey, P.F. Hessburg, H.D. Safford, W.B. Cohen, R.S.H. Kennedy, E.E. Knapp, J.K. Agee, and M. Moeur. 2010. Underestimating risks to the northern spotted owl in fire-prone forests: response to Hanson et al. Conservation Biology 24: 330-333. doi: 10.1111/j.1523-1739.2009.01414.x

Fire Ecology Volume 7, Issue 2, 2011 doi: 10.4996/fireecology.0702051

Lutz et al.: Defining Fire Normals Page 65

Sugihara, N.G., J.W. van Wagtendonk, and J. Fites-Kaufman. 2006. Fire as an ecological process. Pages 58-74 in: N.G. Sugihara, J.W. van Wagtendonk, K.E. Shaffer, J. Fites-Kaufman, and A.E. Thode, editors. Fire in California’s ecosystems. University of California Press, Berkeley, USA. Swanson, M.E., J.F. Franklin, R.L. Beschta, C.M. Crisafulli, D.A. DellaSala, R.L. Hutto, D.B. Lindenmayer, and F.J. Swanson. 2011. The forgotten stage of forest succession: early-successional ecosystems on forest sites. Frontiers in Ecology and the Environment 9: 117-125. doi: 10.1890/090157 Tarnay, L.W., and J.A. Lutz. 2011. Sustainable fire: preserving carbon stocks and protecting air quality as the Sierra Nevada warm. Park Science 28(1): 48-55. Thode, A.E., J.W. van Wagtendonk, J.D. Miller, and J.F. Quinn. 2011. Quantifying the fire regime distributions for severity in Yosemite National Park, California, USA. International Journal of Wildland Fire 20: 223-239. doi: 10.1071/WF09060 Thompson, J.R., T.A. Spies, and L.M. Ganio. 2007. Reburn severities in managed and unmanaged vegetation in a large wildfire. Proceedings of the National Academy of Sciences 104: 10743-10748. doi: 10.1073/pnas.0700229104 US Department of Agriculture [USDA]. 2010. CA_R5_FireHistory09_2. . Accessed 1 August 2010. van Wagtendonk, J.W. 2006. Fire as a physical process. Pages 38-57 in: N.G. Sugihara, J.W. van Wagtendonk, K.E. Shaffer, J. Fites-Kaufman, and A.E. Thode, editors. Fire in California’s ecosystems. University of California Press, Berkeley, USA. van Wagtendonk, J.W. 2007. The history and evolution of wildland fire use. Fire Ecology 3(2): 3-17. doi: 10.4996/fireecology.0302003 van Wagtendonk, J.W., and J. Fites-Kaufman. 2006. Sierra Nevada bioregion. Pages 264-294 in: N.G. Sugihara, J.W. van Wagtendonk, K.E. Shaffer, J. Fites-Kaufman, and A.E. Thode, editors. Fire in California’s ecosystems. University of California Press, Berkeley, USA. van Wagtendonk, J.W., and J.A. Lutz. 2007. Fire regime attributes of wildland fires in Yosemite National Park, USA. Fire Ecology 3(2): 34-52. doi: 10.4996/fireecology.0302034 van Wagtendonk, J.W., K.A. van Wagtendonk, J.B. Meyer, and K. Painter. 2002. The use of geographic information for fire management in Yosemite National Park. George Wright Forum 19: 19-39. Weibull, W. 1951. A statistical distribution function of wide applicability. Journal of Applied Mechanics 18: 293-297. Westerling, A.L., H.G. Hidalgo, D.R. Cayan, and T.W. Swetnam. 2006. Warming and earlier spring increases western US Forest wildfire activity. Science 313: 940-943. doi: 10.1126/ science.1128834 Wotton, B.J., C.A. Nock, and M.D. Flannigan. 2010. Forest fire occurrence and climate change in Canada. International Journal of Wildland Fire 19(3): 253-271. doi: 10.1071/WF09002