Vapor pressure deficit controls on fire ignition and fire spread in boreal forest ecosystems

Introduction Conclusions References


Introduction
Fire plays an important role in regulating carbon and energy fluxes in boreal forest ecosystems on multiple time scales.Carbon emissions during combustion vary considerably as a function of fire weather, topography, permafrost status, and past fire history (Harden et al., 2000;Amiro et al., 2001;Balshi et al., 2007;de Groot et al., 2009;Kasischke and Hoy, 2012).Over a period of decades, fire is one of several types of disturbance that control the distribution of forest stand ages and thus is an important factor regulating the magnitude of terrestrial carbon sinks and sources across regions (Harden et al., 2000;Bond-Lamberty et al., 2007;Balshi et al., 2009;McGuire et al., 2009).On longer time scales, fires influence the carbon cycle through their impacts on permafrost levels (Jorgerson et al., 2010;O'Donell et al., 2011;Nossov et al., 2013), thermokarst dynamics (Myers-Smith et al., 2008), and the stability of deep soil organic carbon layers (Harden et al., 2006).In parallel, fires influence surface energy fluxes by means of several interrelated pathways.By killing overstory trees and initiating succession, fires modify vegetation composition over a period of decades (Johnstone and Kasischke 2005;Goetz et al., 2007;Verbyla, 2008;Johnstone et al., 2010) that, in turn, influences surface albedo (Randerson et al., 2006;Lyons et al., 2006) and sensible and latent heat fluxes (Liu et al., 2005;Amiro et al., 2006;Lee et al., 2011).These fire-ecosystem interactions have the potential to modify regional climate (Rogers et al., 2013) and are likely to influence boreal forest species composition in the future (Goetz et al., 2007;Johnstone et al., 2011;Barret et al., 2011;Mann et al., 2012;Hollingsworth et al., 2013).
Major uncertainties remain in understanding how boreal fire regimes are likely to respond to climate warming.Burned area in boreal North America has increased over the last several decades (Gillet et al., 2004;Kasischke and Turetsky, 2006).Changes in fire disturbance regimes have already occurred (Turetsky et al., 2011;Beck and Goetz, 2011) and future climate change scenarios point towards additional changes (Flannigan et al., 2005;Amiro et al., 2009;Balshi et al., 2009;de Groot et al., 2013) Variations in fire frequency, intensity, size, and spatial and seasonal distributions may increase fire severity and carbon emissions (Kasischke et al., 2010;Turetsky et al., 2011) and affect post-fire vegetation regeneration and succession (Chapin et al., 2004;Goetz et al., 2007;Johnstone et al., 2011;Hollingsworth et al., 2013).In this context, more information about climate controls on various aspects of fire ignition, spread, and extinction is needed to develop more realistic models of future change.
Several studies have explored the relationship between large-scale atmospheric processes and wildfire activity in boreal ecosystems.Skinner et al. (1999) studied circulation anomalies in the mid-troposphere and identified statistically significant correlations between regional annual burned area in Canada and 5000 hPa geopotential heights.Hess et al. (2001) reported linkages between El Niño events, weather patterns in Alaska, and annual burned area in Alaska.Duffy et al. (2005) quantified the impact of variations in atmospheric circulation and weather variables on fire season severity in Alaska, concluding that the average June temperature could explain one-third of the variability in the logarithm of annual burned area.Subsequent studies have combined meteorological variables and fire weather indices with step-wise regression (Flannigan et al., 2005) and multivariate adaptive regression splines (Balshi et al., 2009b) to estimate burned area at large spatial scales under future climate warming scenarios.
While these studies provide insights about the controls behind large fire years in boreal ecosystems at regional level, there is still a lack of information about climate controls on different aspects of fire dynamics at finer scales.This information is needed to develop more mechanistic models describing the ignition and spread rates of fires under scenarios of future climate warming derived from Earth system model simulations.Recently Abatzoglou and Kolden (2011) worked to close this gap by increasing the spatial and temporal resolution of their climate-fire analysis.The authors found that while antecedent weather conditions have negligible importance on fire growth, postignition conditions and the timing of precipitation events are critical determinants of fire size in boreal forests.Introduction

Conclusions References
Tables Figures

Back Close
Full Building on this work, here we developed daily maps of fire progression with a 1 km resolution for approximately 700 wildfires in interior Alaska during 2002-2011 using MODIS active fire observations and fire perimeters from the Alaska Fire Service.To our knowledge this represents the largest dataset of daily wildland fire spread rates in boreal ecosystems.We then used these fire-spread maps with North American Regional Reanalysis data (Mesinger et al., 2006) to investigate how climate influences lightning-driven ignition, and separately, how climate influences fire spread rates and burned area on synoptic, seasonal, and interannual time scales.Our analyses identified consistent relationships between vapor pressure deficit and daily burned area for individual fires, between vapor pressure deficit and daily regional burned area, and between mean summer vapor pressure deficit and annual burned area.The analyses of fire spread provided an improved understanding about temporal and spatial dynamics of wildfires in boreal forest ecosystems and have implications for modeling the terrestrial carbon cycle and the design mitigation strategies.

Study area
The area of study covered an extent of circa 725 000 km 2 in the Alaskan Intermontane Boreal Interior (Nowacki et al., 2001) (Fig. 1).Spanning regions from longitudes 141-159 • West and latitudes 62-69 • North this area included the land between the Canadian border to the east, the Innoko area to the west, the south slope of the Brooks Range to the north and the north slope of the Alaska Range to the south.The vast majority of historical wildland fires in Alaska occurred within this study domain.

Data products
We used active fires information from the MODIS Active Fire Global Monthly Fire Location Product (MCD14ML) (Giglio, 2010) to develop our maps of daily fire spread rates.Figures

Back Close
Full The MCD14ML product consists of an ASCII list of all of the fire thermal anomalies detected in 1 km pixels at the time of overpass of the Terra or Aqua satellites.The fire detection algorithm uses combination of contextual information, a spectral threshold in the middle-infrared, and thermal infrared brightness temperatures to identify fires and separate them from bright clouds and other variations in the background imagery (Giglio et al., 2003).The MCD14ML product includes geographic coordinates for each active fire detection, together with information about the time and date, satellite, scan angle, brightness temperatures, detection confidence and fire radiative power.For this study, we used active fires with low, medium, and high detection confidence levels from both the Aqua and Terra satellites.
To refine fire perimeter identification and identify unburned islands within the fire perimeters we created difference normalized burn ratio (dNBR) spectral index images using the MODIS Vegetation Indices product (MOD13A1).MOD13A1 provides near infrared and medium infrared 16 day composites at 500 m spatial resolution with a sinusoidal projection.For each year the images from the second June composite in the preceding and following years were acquired for the dNBR calculation.Three tiles (h10v02, h11v02 and h12v02) were required to cover to study area.The images were reprojected into a geographic projection (WGS84) using the MODIS Reprojection Tool (https://lpdaac.usgs.gov/tools/modis_reprojection_tool).The Alaska Large Fire Database (LFDB) for the State of Alaska (Kasischke et al., 2002) was used to identify the final perimeters of large wildland fires for our fire spread analysis.LFDB covers the period from 1940-2012 and includes information on a variety of attributes, including fire name, start date, out date, estimated size, and cause.
We obtained meteorological data from the North American Regional Reanalysis (NARR) archived at the National Climatic Data Center (NCDC) (Kalnay et al., 1996;Mesinger et al., 2006).The NARR database combines observational data and models to provide sub-daily information (3 hourly) for 45 meteorological variables at a 32 km spatial resolution from 1979 to 2012.From this dataset we extracted information on relative humidity and air temperature at 2 m heights above surface for the midday interval Introduction

Conclusions References
Tables Figures

Back Close
Full (15.00-18.00p.m. Alaska time).The NARR daily raster images also were projected to the WGS84 geographic projection.
Lightning data was obtained from the Alaska Lightning Detection Network (ALDN) maintained by the Bureau of Land Management Alaska Fire Service (AFS) (Reap, 1991).The network includes 9 sensors distributed in the Alaska interior.Previous studies have documented spatial accuracies for lightning detections processed after 1995 on the order of 2-4 km.Detection efficiencies are assumed to be between 60-80 % in interior Alaska (Dissing and Vervyla, 2003).Upgraded during the 2000 fire season, new IMPACT sensors provide a detection efficiency of 99 % and a median location accuracy of 500 m or better (Vaisala Oyj, 2011).The ALDN provides a vector spatial representation of lightning strikes for each year during 1986-2012.Additionally it also contains attribute fields on location (latitude and longitude) and date (time, day and year), polarity and number of recorded strokes.For this study we extracted data for the period 2002-2011.

Fire spread rate estimation method
We interpolated MCD14ML active fires detection times within LFDB perimeters to retrieve a spatially continuous (raster) representation of the time at which each pixel burned.The output is a raster layer resampled at 500 m spatial resolution of the date of burning for each year during 2002-2011 (Fig. 2).Satellite sensors with high revisit frequencies have proven useful for obtaining information on wildfire occurrence and progression (Boschetti et al., 2010).Previous studies have used active fire information from MODIS to obtain information about fire spread.Loboda and Csiszar (2007), for example, estimated fire spread rates for wildfires in Northern Eurasia solely based on active fire information.More recently Thorsteinsson et al. (2011)  (2012) combined active fires and fire perimeters in an attempt to estimate daily carbon emissions for a set of 41 fires in Alaska.
To retrieve a consistent and reliable fire spread map the number of active fires for a single wildfire must be sufficiently large and their distribution within the perimeter must be relatively homogeneous.While this can be a limitation at lower latitudes, the shorter MODIS revisit period at higher latitudes provides dense active fire distributions within the perimeter of many large fires.Still, the spatiotemporal pattern of active fires can be discontinuous because of factors such as cloud and smoke coverage, and variable fire spread rates.The spatial interpolation approach must be robust to minimize the influence of discontinuous active fires spatial patterns and produce reliable results under a wide variety of conditions.We chose an inverse distance weight interpolation (IDW) for the spatial interpolation of active fire detection dates.IDW is an exact deterministic interpolation approach that calculates values at unknown locations as a weighted average of the values at known locations within a fixed spatial context (Watson and Philip, 1985).The weights are assigned following an inverse relationship to the spatial distances between unknown and known locations.Despite of its simplicity, IDW is a robust approach that requires a limited number of initial parameters and enables a straightforward physical interpretation.This interpolation method has been extensively used to retrieve continuous representations of a variety of environmental variables such as precipitation, temperature, or elevation from point data (Chen and Liu 2012;Didari et al., 2012;Schwendel et al., 2012).While more sophisticated statistical interpolation methods, such as kriging, exist, their more severe assumptions cannot be always met and their increasing complexity commonly requires a larger number of empirical parameters.
The IDW spatial interpolation consists of two steps.The first calculates the weights for each unknown location based on its distance to known locations (Eq.1).The second combines the weights and the values at the known locations to solve for the value at Figures

Back Close
Full the unknown location (Eq.2): where x is the unknown location; x i are known locations; w i are the pair weights; d are the distances between two locations; p is the power parameter; N is the number of known locations used for the calculation; and u is the calculated value at an unknown location.
The interpolation output is sensitive to the initial parameter selection.While small search radii can result in gaps in the final product, very large search areas result in excessively smooth surfaces.Larger p values assign higher weights to neighboring locations.After a preliminary analysis of active fires spatial patterns and the potential parameter combinations for the study area, we selected a circular neighborhood 1 km radius, a power of 2 and a maximum of 10 active fires for the number of known locations.This combination of parameters retained the spatial structure of active fires burning dates while minimizing the number of gaps in the final map.
The day of the year was selected as interpolation variable.Sub-daily estimates at unknown locations were truncated to the day.When only a single active fire existed the fire progression could not be calculated.The annual fire spread maps also provided estimates of fire ignition and extinction dates and locations.Ignition and extinction dates were defined for each fire as the date of the first and last active fire detected within a fire perimeter respectively.However the precise detection of the ignition and extinction dates can be affected by time lags with the time of overpass of satellites and the absence of cloud free imagery.
The IDW interpolation of active fires can potentially retrieve values for unburned islands within the fire perimeters, resulting in an overestimation of the burned area for Figures

Back Close
Full individual fires and aggregated burned area statistics at regional level.To remove these unburned pixels a mask based on the difference normalized burned ratio (dNBR) index was applied.The dNBR is a spectral index that compares pre and post fire values of the normalized burned ratio index (NBR) taking advantage of the differences in spectral signatures before and after fire to identify burned areas and burned area severity.The dNBR has been extensively tested as a proxy for fire severity (Epting et al., 2005;Loboda et al., 2007;De Santis and Chuvieco 2007;Murphy et al., 2008;Verbyla et al., 2008).In this study the dNBR index was calculated from 500 m spatial resolution MODIS imagery (MOD13A1), using June images from the preceding and posterior year of each fire.The corresponding MODIS band 2 (841-876 nm) and band 7 (2105-2155 nm) were used as NIR and SWIR bands for the calculation.A dNBR-based mask was created by applying a threshold between burned and unburned surface, with dNBR values smaller than 0.15 considered as unburned pixels.Subsequently the dNBR mask was used to remove unburned pixels from the fire-spread maps.
Finally fire perimeters were overlaid on the annual fire spread maps in order to assign the fire patches generated in the interpolation to fires recorded in the LFDB dataset.The LFDB included a total of 966 fires for the study area in the period 2002-2011.Of the 966 fires in the database, 289 fire perimeters, with an average size of 74 ha, did not have any active fire within their perimeter.These fires represented 0.35 % of the burned area during 2002-2011.Because of the lack of active fires the information on fire progression for these fires could not be generated.One hundred and seventy fire patches generated from the spatial interpolation of the active fires did not correspond to any fire perimeter in the LFDB.The average size of these patches was 200 ha and they accounted for 0.57 % of the burned area during 2002-2011.Since these fires could not be associated to any event in the LFDB they also were not considered in the analysis.

Back Close
Full Given the limited size of fires without active fires and the small amount of burned area corresponding to fires not registered in the LFDB, their exclusion is unlikely to affect the outcome of the analyses described below.Daily burned area was computed as the number of pixels with the same pixel value (date).The total burned area for individual fires was computed as the sum of the daily burned areas.The annual burned area for the study area was computed as the sum of burned area of individual fires during a given year.Burned areas were calculated after subtracting the unburned island areas as defined by the dNBR mask.For the set of largest fires accounting for 95 % of the Alaska LFDB burned area during the study period, unburned islands accounted for 15 % of the total perimeter area.

Vapor pressure deficit
Vapor pressure deficit (VPD) represents the difference between the saturation vapor pressure and the actual vapor pressure.The VPD compresses meteorological information about the drivers of evaporative demand by the atmosphere into a single variable and can be easily calculated from meteorological variables in existing reanalysis datasets (Bonan, 2008).VPD has been previously used to characterize tree physiology in boreal ecosystems (Zimmermann et al., 2000;Angstmann et al., 2012).It has also been used to analyze fire regimes under various climatic conditions in forest ecosystems (Balch et al., 2008) and as predictor variable in regional and global fire models (Pechony and Shindell, 2009;Silvestrini et al., 2011).We calculated VPD from the air temperature and relative humidity (%) using the two meters variables in the NARR dataset.VPD was calculated with the following expressions (Bonan, 2008): SVP = 610.7 × 10 7.5•T/(237.3+T ) (5) Where SPV is the saturation vapor pressure, T is surface air temperature ( • C), and RH is the relative humidity.To reduce the data volume a single 3 h NARR record corresponding to the 00:00-03:00 UTC time interval (15.00-18.00p.m. Alaska time) was 1319 Introduction

Conclusions References
Tables Figures

Back Close
Full used for each day.Since the NARR records provide variable averages for 3 h cycles and the longitudinal range of the area of study accounts for variations in local time of less than two hours, local time was not dynamically adjusted.
For each VPD pixel, we computed a daily average for the 33 yr period to serve as the baseline value for the calculation of VPD anomalies.The daily VPD maps were further compressed to produce a time series of area-weighted average daily VPD within the study domain shown in Fig. 1.Subsequently the regional time series also were combined to produce a 33 yr average to enable the calculation of regional VPD anomalies within our study region.
For each fire a time series of daily VPD values was constructed by extracting the pixels overlaying the fire perimeter.If more than one VPD pixel overlaid the fire perimeter, the average value of all VPD pixels was calculated.Since the spatial resolution of the daily VPD maps was significantly coarser than that of the fire perimeters, these VPD time series described only the general meteorological conditions near the fire front, and did not capture not finer scale variations caused by topography, fire-driven weather, or other processes.

Time and location of lightning-induced ignitions
To identify lightning-related fire ignitions the ALDN, we used the following approach.First, for each fire, we defined a 1.

Meteorological controls on ignition and initial fire spread rates
The number of lightning strikes varied considerably from year to year (Fig. 3a).The years with the highest number of recorded lightning strikes corresponded to years with the largest burned area, with an R 2 of 0.52 (P < 0.01).Most lightning strikes occurred during June and July, with these months combined accounting for 89 % of the lightning strikes during 2002-2012.Most of the lightning outside of this period occurred during May (3.8 %) or August and September (6.9 %) (Fig. 3a).Fire ignitions in boreal forest ecosystems of Alaska required the coincidence in time and space of lightning strikes and high VPD values (Fig. 3b).For each year between 2002 and 2011 daily VPD values during lightning strikes leading to fire ignitions were significantly higher than VPD values for all other lightning strikes (z test, 0.05 significance level) (Fig. 3b).The seasonal pattern of fire start dates follows a similar pattern to that of lightning strikes, with fires less clustered in the central summer months (Fig. 3d).The majority (72 %) of the lightning-induced fires occurred in June and July, 14 % in May and 13 % in August and September (Fig. 3c).The total number of annual lightning strikes was significantly related to the number of lightning-induced fires (R 2 = 0.63, n = 10, p = 0.01) and the total number of fires (R 2 = 0.44, n = 10, p = 0.01) (Fig. 3d).
To identify differences in meteorological controls on fires of various sizes we grouped fires in four categories: small fires (0-1000 ha); medium fires (1000-10 000 ha); large fires (10 000-50 000 ha); and very large fires (> 50 000 ha).The analysis revealed that fire lengths varied from a single day to nearly 90 days.A positive relationship existed between fire duration and fire size (Fig. 4).Many small fires persisted for a few days and were often characterized by an initial period of limited growth (low daily burned area) followed by an extinction during the first period of negative VPD anomalies.In a few cases small fires were dormant for longer periods, and then began to grow after several weeks.For all of the different fire size classes, mean VPD anomalies averaged 1-5 days prior to ignition were all positive, indicating that lower fuel moisture levels appeared to be a necessary pre-requisite for lightning triggered ignition.These positive VPD anomalies, however, were not significantly different as a function of fire size, providing evidence that pre-ignition conditions did not play an important role in determining the size of the final fire perimeter (Fig. 5a).
In contrast, VPD anomalies in the first five days after ignition were a significant driver of fire size.The VPD anomalies in this period were negative for small fires, and were positive and considerably higher both for large and very large fires (Fig. 5b).These results imply that fire weather immediately following ignition is crucial for the growth of large wildland fires.Small fires tended to occur during oscillating VPD conditions under which the fires were unable to grow and eventually were extinguished during periods of wet weather.In contrast, sustained periods of positive VPD anomalies immediately after ignition may have contributed to the resilience of fires that eventually grew to very large sizes.One possible mechanism for this diverging fire behavior might be the establishment of larger fire perimeters for fires with a higher initial VPD, thus creating the potential for a larger number of surface or subterranean hotspots.These hotspots may allow these fires to survive during relatively long periods of wet weather, and then expand later when the climate is more favorable for fire spread.Qualitative evidence for this hypothesis is provided in Fig. 6.For many of the larger fires during 2004, 2005, and 2009, long quiescent period preceded periods of rapid burned area growth in mid or late summer.

Meteorological controls on the burned area of individual fires
Vapor pressure deficit regulated both daily fire spread rates and the final size of individual fires.Time series of daily burned areas for large fires were significantly correlated with daily VPD anomalies.For the 160 fires greater than 5000 ha in our 2002-2011 dataset, correlation coefficients between VPD anomalies and daily burned area varied between 0.03 and 0.77, with a mean of 0.14.This set of fires had a mean duration 1322 Figures

Back Close
Full .Daily fire spread rates were often highly variable within the same fire, alternating between periods of large daily fire spread rates and long periods of quiescence (Figs. 6 and 7).Positive VPD anomalies averaged over fire lifetimes was positive and without significant differences for the different group sizes (Fig. 5c), indicating that drier than normal atmospheric conditions were needed to sustain fire spread rates.Cumulative VPD anomalies during the lifetime of each fire in this subset explained 60 % of the variance in final fire size (Fig. 8).These results provide evidence that after ignition, the persistence of extended periods with positive VPD anomalies allowed fires to grow to larger sizes.Fire extinctions also occurred for many fires during periods with negative VPD anomalies.Mean VPD anomalies during the fire extinction, defined as the five days after the last recorded active fire, were negative for all fire sizes, with lower (more negative) values for smaller fires (Fig. 5d).Since larger fires often persisted into August when monthly mean VPD levels were lower, their extinction may require a less intense period of wet weather than smaller fires that are extinguished in mid-summer.

Meteorological controls on regional burned area
For many fires in interior Alaska, the timing of dry and wet periods during the fire season was similar.For example, for the four fires shown in Fig. 7, many of the high and low VPD periods occurred at similar time intervals, even though these fires were located in different parts of the state.As a consequence of the regional coherence in meteorology, ignition, daily fire spread rates, and fire extinction periods were often synchronized across different fires.In Fig. 6a-c distinct temporal clusters in different years (Fig. 6).The correlation between daily VPD and daily burned area, averaged at a regional level, varied from year to year, with stronger correlations observed in general during high fire years and weaker correlations observed in low fire years (Table 2).Large positive VPD anomalies that persisted for 10 days or more during 2002, 2004, 2005, and 2009 fire seasons were responsible for most of the burned area during the decade (Fig. 9).Periods with positive VPD anomalies accounted for 79 % of the total regional burned area during 2002-2011, yet accounted for only 40 % of the time during the fire season (Fig. 9).

Meteorological controls on interannual variability of regional burned area
The meteorological controls identified above for individual fires and for regional burned area during individual fire seasons allowed us to develop hypotheses about climate drivers of interannual variability in wildfire activity.We hypothesized specifically, that years with more positive summer VPD anomalies would have higher levels of cumulative burned area because these drier conditions would enable more ignitions and allow for longer intervals when fire spread rates were high.During our study period of 2002-2011, we found a significant positive correlation between the sum of positive VPD anomalies during the fire season (DOY 150-240) and annual burned area (R 2 = 0.71, n = 10, p < 0.01) for the 2002-2011 period.The correlation was weaker but still highly significant for the longer 1979-2011 period of the NARR dataset (R 2 = 0.45, n = 33, p < 0.01) (Fig. 10).These relationships supported our hypothesis and provided evidence that VPD was a key driver of regional fire dynamics on longer timescales.We also investigated how year-to-year changes in the seasonal timing of fire ignitions and fire extinctions contributed to interannual variability in regional burned area.The mean annual ignition dates varied from early-June (DOY 175) to mid-July (DOY 191) (Table 1).Early annual ignition dates did not correspond to large annual burned area years (R 2 = 0.39, n = 10, p = 0.01).In contrast, the mean extinction dates ranged from early-June (DOY 172) to early-September (DOY 221), and were significantly correlated defined as the difference between annual mean extinction and ignition dates varied from 6 to 30 days (Table 1), and there was a positive relationship between mean fire duration and annual burned area (R 2 = 0.60, n = 10, p = 0.01) (Fig. 11).

Implications of meteorological controls for fire management
Each year, some of the fires ignited in late spring or early summer remained dormant for extended periods until more extreme periods of high fire weather caused these perimeters to rapidly grow (Fig. 6).Other fires were ignited during periods of sustained dry conditions, and rapidly expanded within a few days of the initial lightning strike.The former set of fires may represent a target of opportunity for fire suppression efforts designed to maintain the current fire regime in future decades, given large unidirectional changes in burned area expected from climate warming.To investigate the importance of these different fire types in contributing to annual burned area, we conducted the following analysis.First, we quantified the number of fires that stayed below a certain size threshold (5000 ha) for the first 10 days after ignition.From the 160 fires that consumed more than 5000 ha during 2002-2011, 61 fires required more than 10 days after ignition to reach this area threshold.These fires accounted for 30 % of the total burned area during the period.We also analyzed the 50 largest fires within the study domain during this interval.These fires accounted for 69 % of the total burned area during 2002-2011.For this set, 23 fires remained below a threshold of 5000 ha during the first 10 days after ignition, and 32 fires burned less than 20 % of their total size during the first 10 days after ignition (Fig. 12).Fires that initially remain in a quiescent/smoldering phase could be targeted for suppression, although as described below more information is needed to identify landscape characteristics that enable fires within this class to grow to very large sizes.Introduction

Conclusions References
Tables Figures

Back Close
Full

VPD, satellite-based information and fire danger estimation
We found a strong positive relationship between VPD and burned area across different temporal and spatial scales.This relationship stands at fine temporal resolution (daily VPD against daily burned area) for individual fires, at a regional level within a single fire season, and across different years for the study domain as a whole.VPD also was implicated as an important climate regulator during multiple fire stages including, specifically, the probability that a lightning strike triggered ignition, during periods of initial fire spread, for daily burned area variations in larger fires, and the timing of fire extinction.Our findings suggest that by combining satellite estimates of fire spread with high-resolution climate information, there is potential to considerably improve fire weather systems and fire behavior models.Fire danger and fire behavior in boreal forests have normally been estimated with weather-based empirical indices developed before satellite information was widely accessible.For instance the Canadian Forest Fire Weather Index (FFWI) System (Van Wagner, 1987)

Implications for Earth system models
Knowledge of fire-spread dynamics boreal forests is fundamental for developing realistic projections of burned area for different climate change scenarios.Models of terrestrial carbon dynamics in boreal ecosystems rely on future burned area estimates at a regional level (Mann et al., 2012;Balshi et al., 2007Balshi et al., , 2009)).These estimates are commonly based on linear relationships between meteorological variables and burned area.Duffy et al. (2005), for example, applied the average June temperature as a predictor of burned area at a regional level.Bachelet et al. (2005) modeled burned area as a function of a drought index within a global vegetation model to study carbon dynamics in Alaska.Rupp et al. (2007) estimated burned area using an empirical relationship driven by growing-season average temperature and cumulative precipitation, and integrated this information within a landscape-level model of vegetation dynamics.Several studies have explored the relationships between fire weather indices and burned area to predict future burned area using climate model simulations of the 20th and 21st centuries (Flannigan et al., 2005;Balshi et al., 2009;De Groot et al., 2013).
Our results show that multiple interacting mechanisms likely govern the emergent relationships observed between climate and burned area at coarse spatial and temporal scales.Large fires commonly last for several weeks, alternating periods of high activity and dormancy associated with VPD dynamics.For extreme high fire years, our analysis provides evidence that several environmental factors must align.First, extensive lightning storms and drier than normal conditions need to occur during May and June to enable the creation of hundreds of fire ignitions during the early part of the fire season.Many of these storms need to be followed by several days of dry weather to enable small fires to reach some minimum size that allows them to survive during subsequent wet periods.In the middle of the growing season, high pressure systems need to establish persistent periods of dry and hot weather.Large, multi-day positive VPD anomalies occurred during both the 2004 and 2005 fire seasons (Fig. 9).Finally, near the end of the fire season in, climate conditions need to be drier than average to Introduction

Conclusions References
Tables Figures

Back Close
Full reduce the probability of an early, region-wide extinction event.This complex behavior is unlikely to be captured by monthly mean climate or weather indices, and thus an important direction for future work is to develop more sophisticated predictive models that take into account daily weather patterns and atmospheric processes that regulate convection and lightning storms.The inclusion of these processes in fire models, along with better representation of ecosystem and topographic controls, will lead to improved estimates of regional burned area changes and carbon dynamics in Earth system models.

Implications for fire management
Understanding the meteorological controls behind fires in boreal forests opens the possibility for the design and implementation of fire suppression operations to maintain historical or current fire regimes.Today, the extent of the area and limited accessibility to many remote fires in Alaska increases costs for suppression.Fire suppression and management resources in Alaska are allocated based on risks to human settlements and ecosystems; the cost of suppression activities must be commensurate with the value of the resources needing protection (Alaska Interagency Wildland Fire Management Plan, 2010).Fire suppression activities are categorized in different management options in which areas with settlements and infrastructure receive the highest priority.However, given the magnitude of the estimated carbon pool in Alaskan boreal forests, the potential future implementation of international carbon tax mechanisms or other climate policies could justify more extensive fire suppression activities in the near future.Fire suppression may be even more valuable from the perspective of adaptive resource management: maintaining current burned area levels may slow changes in species composition and other changes in ecosystem function ignited in late spring or early summer that remained dormant for extended periods until more extreme periods of high fire weather caused these perimeters to rapidly grow (Fig. 6).The discontinuous progression of some wildfires offers the opportunity to target selected fires that have the potential to grow to very large sizes and identify time windows in which suppression efforts could be implemented.
Large fires in Alaska account for a very significant share of the annual burned area.For instance the twenty largest fires of 2004 accounted for 80 % of the annual burned area.Our analyses indicate that for more than half of the largest fires (32 out of 50) during the period 2002-2011, less than 20 % of their burned area occurred in the first 10 days after the ignition (Fig. 12).Concentrating fire suppression operations on this limited number of fires with slow initial growth might have considerably reduced the total burned area during this period.An important challenge for the future is to more effectively identify fires that represent ideal suppression opportunities based on their location and timing of ignition.

Conclusions
We combined MODIS active fires (MCD14ML) information, MODIS imagery (MOD13A1) and ancillary historic fire perimeter database (Alaska Large Fire Dataset) to produce a dataset of fire-spread maps (daily burned area) in Alaska for the period of 2002-2011.This product allowed us to identify fire ignition location and dates and fire progression for individual wildfires.We used the fire spread maps with spatially and temporally explicit climate records to examine relationships between the driving gradient for vegetation water loss (VPD) and wildfire activity in boreal forest ecosystems.We found that VPD exerted controls on burned area during multiple stages of fire growth, including ignition probabilities, initial rates of fire spread, daily burned area, and likelihood of extinction.Lightning strikes were more likely to trigger fire ignitions when they occurred during periods of high VPD.The probability of a fire growing into large or very large sizes was enhanced when VPD levels were above average in the first five days Figures following ignition.Our analyses underline the importance of fine temporal patterns of VPD and lightning strikes and fire season length in regulating year-to-year variations in regional burned area.The findings of this study improve the spatial and temporal characterization of the climate controls regulating the contemporary fire disturbance regime and highlight in the importance of developing improved models to predict future changes in fire regimes.Information on fire progression of large wildfires may create potential opportunities for adaptive fire management in the future.The link between daily burned area from the fire spread maps and meteorological and landscape variables should allow for a more precise estimation of fire severity and fuel consumption and refine estimates of daily carbon emissions from boreal forest wildfires.Introduction

Conclusions References
Tables Figures

Back Close
Full       Relationship fire duration and fire size.Small fires were defined as having sizes between (0-1000 ha, medium fires between 1000-10 000 ha, large fires between 10 000-50 000 ha), and very large fires as greater than 50 000 ha.The percentage of number of fires and burned area for each group was respectively 63.5 and 2 % for small fires; 24.3 % and 18 % for medium fires; 10 and 39.9 % for large fires; and 2.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | used MODIS active fires to monitor the evolution of a single fire event in Iceland and Kasischke and Hoy Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 5 km radius neighborhood as around the 1 km location of the first active fire detected within each fire perimeter.Any lightning strikes occurring during the same day or the day before of the active fire observation was classified as lightning that triggered a fire ignition.All other lightning observations were classified as those that did not trigger an ignition.The annual average VPD values at the location and day of the lightning-related fire ignitions were computed and compared against the annual average VPD values of all others lightning strikes through a test of significance (z test).
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | of 39 days and accounted for 90 % of the burned area during 2002-2011.Examples of these relationships for 4 large fires during 2004, including the Boundary, Dall City, Winter Trail, and Billy Creek fires are shown in Fig. 7 , for example, the end of the 2005 fire season occurred abruptly as a negative VPD anomaly extinguished almost all of the fires across the state between days 230 and 240.In 2009 the end of the fire season also was abrupt in the northern part of the state, and occurred about 15 days earlier than in 2005.Ignitions appeared less synchronized than extinctions, but also showed 1323 Discussion Paper | Discussion Paper | Discussion Paper | with annual burned areas (R 2 = 0.68, n = 10, p = 0.01).Mean fire season duration, Discussion Paper | Discussion Paper | Discussion Paper | uses weather variables interpolated from a number of weather stations to estimate six empirical indices related to soil moisture and fire behavior.The connection of spatially explicit VPD information and remote sensing-based fire spread rates shows potential to optimize parameters and develop more mechanistic models of fire spread.Vapor pressure deficit estimated from NARR variables can provide spatially explicit information related to fuel moisture, which improves the interpolation of basic weather variables in regions with a limited network of weather stations.Furthermore, the combination of active fire points from moderate spatial resolution satellite remote sensing and spatially explicit weather variables provides an approach to improve the understanding of the relationships between weather and fire spread as a first step toward developing more advanced fire spread indices.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | (i.e.permafrost melt) expected with climate warming.The development of decision systems including information on fire progression drivers, weather forecasts, environmental variables (including land cover, topography, and landscape fragmentation) have potential to improve the efficiency of fire suppression operations.As previously presented, some fires were 1328 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 2 .
Fig. 2. Active fires and daily fire spread rates for the Billy Creek fire that burned during the summer of 2004.(A) Active fires within the LFDB perimeter; (B) the fire spread map after spatial interpolation; (C) the difference normalized burn ratio index (dNBR) estimated from 500 m MODIS data from June 2003 (pre-fire) and June 2005 (post-fire; D) fire spread map after applying dNBR-based mask to remove unburned islands.

Fig. 3 .
Fig. 3. (A) Annual number of lightning strikes from ALDN disaggregated by month in the period 2002-2011.(B): Average VPD values for lightning strikes leading to fire ignitions (light grey) and average VPD values for all other lightning strikes (dark grey).(C) Annual number of lightninginduced fires disaggregated by month in the period 2002-2011.(D) Annual number of fires disaggregated by month during 2002-2011.

Fig. 5 .
Fig.4.Relationship fire duration and fire size.Small fires were defined as having sizes between (0-1000 ha, medium fires between 1000-10 000 ha, large fires between 10 000-50 000 ha), and very large fires as greater than 50 000 ha.The percentage of number of fires and burned area for each group was respectively 63.5 and 2 % for small fires; 24.3 % and 18 % for medium fires; 10 and 39.9 % for large fires; and 2.7 and 40 % very large fires.

Fig. 6a .)Fig. 8 .Fig. 10 .Fig. 11 .
Fig. 6a.Top plate: length and daily burned area of individual fires in 2004, 2005 and 2009.Each line represents a recorded fire event in each year.Fires are organized from north (top) to south (bottom) within the study domain.The color bar represents daily burned area, with values normalized by annual regional daily burned area values between 0 and 100.Light blue (color code −20) represents days between ignition and extinction in which burned area was not registered.Bottom plate: daily VPD values averaged for interior Alaska (blue), 1979-2011 daily VPD average (red) and 75th percentile values in the period 1979-2011 (green).Interior Alaska was defined by the longitudes 142-154 • W and latitudes 63.5-67.5 • N. The horizontal axes in both plates represent the day of the year.
Fig. 12. (A)Percentage of total burned area within ten and twenty days after fire ignition for fires larger than 50 000 ha during the period 2002-2011.(B) Time series of daily (red line) and cumulative burned area (blue line) for three large fires in 2004.The fire in the lower middle panel had burned less than 35 % of its final burned area within the first forty days after its ignition, which made it a suitable candidate for fire suppression.

Table 1 .
Annual burned area, number of fires, mean ignition and extinction dates, and fire season duration.

Table 2 .
Annual VPD summer average (June-August), slope of the burned area-VPD relationship and r squared.Figures