Many boreal forests grow in regions where climate is now warming rapidly. Changes in these vast, cold forests have the potential to affect global climate because they store huge amounts of carbon and because the relative abundances of their different tree species influence how much solar radiation reflects back to space. Both the carbon cycling and albedo of boreal forests are strongly affected by wildland fires, which in turn are closely controlled by summer climate. Here we use a forest disturbance model in both a retrospective and predictive manner to explore how the forests of Interior Alaska respond to changing climate. Results suggest that a widespread shift from coniferous to deciduous vegetation began around A.D. 1990 and will continue over the next several decades. This ecological regime shift is being driven by old, highly flammable spruce stands encountering a warmer climate conducive to larger and more frequent fires. Increased burning promotes the spread of early successional, deciduous species at the expense of spruce. These striking changes in the vegetation composition and fire regime are predicted to alter the biophysics of Alaska's forests. The ground will warm, and a surge of carbon emission is likely. Our modeling results support previous inferences that Alaska's boreal forest is now shifting to a new ecological state and that positive feedbacks to global warming will accompany this change.
Introduction
The responses of ecosystems to changing climate can be sudden, nonlinear, and sometimes lead to very different self-sustaining states (Higgins et al., 2002; Scheffer, 2009). Such ecological regime shifts involve sweeping changes in ecosystem structure, function, and species composition that entail changes in chemical and biophysical exchanges with the rest of the biosphere (Scheffer and Carpenter, 2003). Ecological regime shifts are of particular concern when they involve entire biomes and threaten to trigger positive feedbacks to human-caused warming (Lenton et al., 2008).
The circumboreal forest is of special concern in regard to climate change and ecological regime shifts because of the enormous amounts of carbon (C) it stores. Boreal forests cover some 16 million km2 globally, and every year about 8% of the atmospheric pool of carbon (C) cycles through them (McGuire et al., 2009). Altogether, boreal forests currently contain 10–20% of the world's vegetation C and are underlain by about 30% of the world's soil C (McGuire et al., 2009). This abundant soil C has accumulated in boreal forests because decomposition rates there are low due to the wet, cold soils.
Temperatures are rising rapidly in many boreal forest regions today (Hinzman et al., 2005). As soil temperatures warm and decomposition rates increase, more CO2 and methane (CH4) will be released to the atmosphere (Schuur et al., 2009). Wildland fires are the other important recycler of vegetation and soil C in boreal forests (Soja et al., 2006; Boby et al., 2010). By removing shading vegetation, insulating moss, and soil organic horizons, fires can cause soils to warm dramatically (Viereck et al., 2008), which speeds the decomposition of soil organic matter. If the boreal forest “freezers” defrost through a combination of rising air and soil temperatures caused by a combination of warmer climate and increased burning, the resulting release of trace gases could trigger significant positive feedback to global warming.
Climate is now changing rapidly in Alaska (Hinzman et al., 2005). Over the last century, mean annual temperature in Interior Alaska rose 1.4 °C, and precipitation decreased by 11% (Wendler and Shulski, 2009; Wendler et al., 2010). Growing season length has increased by about three days per decade since 1970 (Euskirchen et al., 2010). Temperatures in the top meter of the ground have risen in Alaska by 0.7 °C per decade over the last 40 years, and widespread thawing of permafrost seems inevitable if this trend continues (Jorgenson et al., 2010). How stable is Interior Alaska's forest in the face of these rapid changes in climate?
PORTENTS OF IMPENDING CHANGE?
Previous research suggests that boreal forests in general (Soja et al., 2006) and Interior Alaska's forests in particular are approaching a major ecological threshold. White spruce in Interior Alaska is sensitive to temperature-induced drought (Barber et al., 2000), and recent warming has caused reductions in its growth throughout the region (McGuire et al., 2010). Chapin et al. (2004) speculated that white spruce forests in Interior Alaska might be replaced by aspen parkland or grassland with a future temperature rise of just 2 °C. Similarly, Calef et al. (2005) inferred that a combination of warmer drier summers and more frequent fires could cause a significant expansion of deciduous stands at the expense of white spruce and tundra vegetation. Based on a review of the key processes that structure the boreal forest and that are believed to be sensitive to warming temperature (e.g., fire regime, tree regeneration, thermokarst), Chapin et al., (2010) inferred that “…the Alaskan boreal forest is on the cusp of potentially large nonlinear changes in structure and functioning.”
Focusing on the fire regime, Kasischke et al. (2010) suggested Interior Alaska's forest is vulnerable to sweeping ecological changes triggered by two fire effects. The first is a shortened firefree period that could reduce the opportunity for the recruitment of black spruce seedlings (Johnstone and Chapin, 2006) and so allow an increase in deciduous vegetation. The second is the increased combustion of organic soils, which could significantly alter the soil-legacy effects (Harden et al., 2006) that today favor the self-replacement of black spruce. This might further encourage the spread of short-lived deciduous vegetation at the expense of black spruce (Johnstone et al., 2010). Similarly, Kurkowski et al. (2008) suggested that a warming climate could trigger changes in the relative importance of different post-fire successional pathways in such a way as to cause the widespread replacement of spruce by deciduous trees.
Beck et al. (2011) presented compelling evidence that forest productivity has declined over the last several decades in Interior Alaska as summers have warmed and trees become more drought stressed. Based on this evidence, Beck et al. (2011) warned of an impending ecological regime shift in Interior Alaska. These threads of inference and evidence converge on the conclusion that a sudden shift in the structure and functioning of Interior Alaska's boreal forest is possible as climate continues to warm. What has been missing is a way to combine these different threads into a unified, testable hypothesis of forest function and forest response to climatic drivers.
RESEARCH QUESTIONS AND MODELING APPROACH
Our primary research questions are these: How sensitive is the boreal forest in Interior Alaska to ongoing climate forcing? What key processes mediate this sensitivity? Does an integrative, modeling analysis of this forest substantiate the suggestions that an ecological regime shift is imminent?
To answer these questions, we use the ALFRESCO computer model, a forest-disturbance model that integrates fire, forest succession, and climate change across spatially explicit landscapes of regional extent. By “spatially explicit” we mean that the ALFRESCO model tracks the changing vegetation and fire history of a myriad of separate, contiguous sites on real topography through time. Its historical, climate-influenced, and spatially explicit structure endows ALFRESCO with the ability to estimate transient change in albedo, heat fluxes, and carbon budgets on a forested landscape experiencing rapid climate changes.
The main virtue of a computer model is its ability to meld inferences about diverse ecological processes into a single, unified hypothesis about how an entire forest functions and consequently how it might respond to changing climate. The relative simplicity of boreal forests makes modeling their responses to climate change more straightforward than for many other biomes. Many boreal forests still exist in primeval states because of their remoteness from human activities, and a typical boreal forest contains just a handful of tree species. Disturbance-regeneration cycles in these forests are triggered mainly by lightning-ignited fires (Kasischke and Turetsky, 2006; Kasischke et al., 2010) whose frequency and size are controlled by relatively few weather parameters (Duffy et al., 2005). Here we use ALFRESCO in both a retrospective and predictive manner to explore how boreal forests of Interior Alaska respond to changing climate.
Materials and Methods
STUDY AREA
Interior Alaska includes approximately 47 million ha between the Alaska and Brooks Range. The regional topography consists of a series of large tectonic basins separated by low mountain ranges and interconnected by meandering rivers. The combination of a diverse topography with low sun angles in summer sets up a complex mosaic of radiative microclimates (Kurkowski et al., 2008), which give rise to a corresponding mosaic of soil microclimates. These soil microclimates influence both the distribution of permafrost and the thickness of the active layer, which together strongly influence vegetation distribution, decomposition rate, and soilorganic horizon thickness (Harden et al., 2006). Interior Alaska lies within the zone of discontinuous permafrost (Osterkamp and Romanovsky, 1999), and stand-replacing fires are the most frequent forest disturbance (Viereck, 1983). The most abundant tree species in the region is black spruce (Picea mariana), which usually grows on the coldest, wettest soils with the thickest organic horizons. Aspen (Populus tremuloides) occupies the warmest, driest sites where the organic mat is thin or absent. White spruce (Picea glauca) and birch (Betula neoalaskana) occupy sites that have microclimates and soil characteristics intermediate between black spruce and aspen. Detailed descriptions of this forest and its physical environment can be found in Chapin et al. (2006). Kasischke et al. (2010) described fire regimes in Alaskan forests, Johnstone et al. (2010) described fire effects on tree regeneration there, and Euskirchen et al. (2010) reviewed forest-climate interactions.
THE ALFRESCO MODEL
Overview
ALFRESCO is a vegetation disturbance model to which biogeochemical fluxes can be appended ( http://www.snap.uaf.edu/resource_page.php?resourceid=11). Because vegetation cover in ALFRESCO is a function of topographic position and time-since-last-fire (TSLF), the overall species composition and stand-age distribution of the study region emerges from the continually changing responses of the overall forest mosaic to local fire histories. In this way the model incorporates both the emergent properties of forests and the crucial role played by climate-fire linkages in structuring them. ALFRESCO's strength lies in its ability to simulate rapid changes in the forest's overall species composition and tree-age distribution, which can then be used to quantify biophysical properties using estimates of the rates and magnitude of heat fluxes, and carbon emissions within specific vegetation types (e.g., 40-year-old aspen versus 100-year-old black spruce).
ALFRESCO quantifies interactions among climate, topography, vegetation distribution, post-fire succession, and age-dependent flammability specifically for Interior Alaska's boreal forest (Rupp et al., 2000, 2006, 2007). The model simulates the transient states of vegetation at a spatial scale of 1 × 1 km and at annual time steps in response to climate-driven changes in the fire regime (fire frequency, fire severity, and fire size). The fire regime is simulated stochastically and constrained by climate and vegetation type. Ignition is determined randomly and is a function of pixel flammability. The model uses a cellular automaton approach in which an ignited pixel can spread fire to any neighboring pixel (Brubaker et al., 2009). Fire spread depends on the flammability of the adjacent pixels as modified by the presence of firebreaks such as non-vegetated mountain slopes, rivers, and lakes. Fireclimate rules come from a statistical model (Duffy et al., 2005) that explains 79% of the interannual variability in area burned in Interior Alaska as a function of monthly climate parameters over the period 1950–2005. Post-fire successional pathways are determined by the burn severity in each pixel, which is modeled as a function of fire size and topography (Duffy et al., 2007).
We used the age structure, species composition, and fire regime of the present-day forest to establish model rules and to calibrate ALFRESCO within reasonable values of observed parameters taken from the extensive literature on post-fire succession and its relation to topography and pre-fire stand type in Interior Alaska (Chapin et al., 2006; Kurkowski et al., 2008; Johnstone et al, 2010). For more details on model calibration, the reader is referred to the model's user manual ( http://www.snap.uaf.edu/resourcepage.php?resourceid=11). We estimated annual area burned back to A.D. 1860 using the climate-area burned relationship of Duffy et al. (2005) in conjunction with reconstructed and downscaled data sets of monthly temperature and precipitation. Climate reconstructions for 1860–1900 are from the Potsdam Institute for Climate Impact Research (Leemans and Cramer, 1991); spatially interpolated, observational climate data for 1901–2002 are from the Climate Research Unit (Mitchell and Jones, 2005). Further details on the climate inputs to the model are given below.
We tuned ALFRESCO's pixel-scale rules governing flammability, fire spread, and hence annual area burned using estimates of annual area burned from two sources. For 1950–2007, we used fire management records (Kasischke et al., 2002, 2010; Bureau of Land Management—Alaska Fire Service, 2008). For 1860–1949, we used back-casts of annual area burned made using the fire-climate relationship (Duffy et al., 2005). We assessed model accuracy by iteratively running the model forward from 1860 and comparing the simulated fire-regime metrics and vegetation composition with the back-cast fire history (1860–1949), with he observed fire history (1950–2007), with regional forest composition in 2001 (Homer et al., 2004), and with comparisons between observed and modeled power law exponents. Power-law exponents, β, are a sensitive metric useful in characterizing fire regimes across multiple spatial and temporal scales (Malamud et al., 2005; Song et al., 2006; Pueyo, 2007). We applied ALFRESCO to the future by inputting the downscaled climate predictions of the five global circulation models found to work best over Alaska (Walsh et al., 2008).
ALFRESCO Vegetation
The model uses the minimum number of vegetation types needed to realistically depict the stand mosaic that comprises the boreal forest of Interior Alaska today. The vegetation data used in model spin-up were reclassified from the 1990 AVHRR vegetation classification of Interior Alaska ( http://agdcftpl.wr.usgs.gov/pub/projects/fhm/vegcls.tar.gz) and the 2001 National Land Cover Database vegetation classification ( http://www.mrlc.gov) (Homer et al., 2004) into tundra, black spruce, white spruce, or deciduous vegetation. Deciduous vegetation includes aspen, cottonwood (Populus balsamifera), and birch. In the first few decades after a fire, deciduous vegetation includes seedlings and saplings of the abovementioned tree species plus herbaceous species like grass, sedge, horsetail (Equisetum spp.), and fireweed (Epilobium angustifolium) (Mann and Plug, 2002; Johnstone and Chapin, 2006). Remote sensing is unable to distinguish black spruce from white spruce (Kurkowski et al., 2008), so ALFRESCO uses aspect and topographic position as a predictor of spruce species, with black spruce located on north-facing slopes, poorly drained flats, and toe slopes. These are all sites where organic soils are common, so in this way the model incorporates the role that soil organic matter plays in maintaining ecological inertia on this landscape (Harden et al., 2006).
ALFRESCO Climate
Two spatially explicit (0.5° × 0.5°) data sets provide ALFRESCO with monthly averages of temperature and monthly totals of precipitation between A.D. 1860 and 2002. The Potsdam Institute for Climate Impact Research (PIK) data set is a modified version of that presented in Leemans and Cramer (1991) (McGuire et al., 2001). We used PIK data for the years 1860–1900 and CRU data ( http://www.cru.uea.ac.uk/∼timm/grid/CRU_TS_2_0.html) for the years 1901–2002. Predictions of future climate come from the general circulation models (GCMs) that best match historical weather records in Alaska (Walsh et al., 2008). These five are the MPI_ECHAM5, the GFDL_CM2_1, MIROCC3_2_MEDRES, UKMO_HADCM3, and CCCMA_CGCM3_1. We used the midrange emission scenario (A1B) with each of these models. To span the several years missing between the CRU and GCM-predicted weather data, we used output of whichever of the five models best predicted the observed area burned in Interior Alaska in that year. Details about how we down-scaled GCM-predicted climate to the 1×1 km scale used in ALFRESCO can be found at http://www.snap.uaf.edu/.
Inputting Climate and Calibrating Flammability
Climate data are used in ALFRESCO to define the flammability coefficient of each pixel. Pixel flammability also varies according to its vegetation cover, which depends on TSLF and topographic position. Monthly fields of average temperature and total precipitation at a 1 × 1 km scale are superimposed over vegetation mapped at the same scale. The flammability coefficients of different vegetation types under varying climatic conditions are iteratively tuned within ecologically realistic bounds so that the model's output matches the observational records of several key parameters ( http://www.snap.uaf.edu/resource_page.php?resourceid=11). The first of these is annual area burned as recorded by the Alaska Fire Service of the Bureau of Land Management ( http://afsmaps.blm.gov/) (Kasischke et al., 2002; Kasischke and Turetsky, 2006). These historical records provide a target for tuning ALFRESCO's flammability coefficients between 1950 and 2007. We extended the annual area burned record into the past by using the climate-area burned relationship inferred by Duffy et al. (2005) in conjunction with the CRU and PIK climate data sets to back-cast annual area burned between 1860 and 1949. While comparing ALFRESCO's simulated annual area burned to these historical and back-cast estimates of area burned, we iteratively adjusted the rules governing pixel flammability within ecologically reasonable limits until the model yielded predictions of annual area burned that were consistent with observed and back-cast values (Fig. 1).
Biophysical Effects in ALFRESCO
Landscape-scale changes in vegetation and fire regime strongly influence surface radiation budgets and C cycling in Alaskan forests (McGuire et al., 2006; Balshi et al., 2007; Chapin et al., 2009; Euskirchen et al., 2009, 2010). By combining ALFRESCO's predictions with published estimates of the effects of TSLF and vegetation type on various biophysical parameters, we can predict the effects of changes in the Alaskan boreal forest at local, regional, and global scales. We used the measurements of Randerson et al. (2006) to define a TSLF-albedo function. Similar trends in albedo change over time result from using the TSLF-albedo estimates of Amiro et al. (2006), Lyons et al. (2008), and Euskirchen et al., (2007, 2010). We based estimates of net solar radiation, latent heat flux, ground heat flux, and sensible heat flux on the chronosequence data of Liu et al. (2005) and Liu and Randerson (2008). To estimate C emissions from burning boreal forests, we mainly followed Turquety et al. (2007). Estimates of the amount of C in dry matter come from Harden et al. (2000) and Soja et al. (2004).
Approximately 40% of Interior Alaska is classified as poorly drained and is covered by thick (>25 cm), often water-saturated, organic soil horizons (Harden et al., 2001). Because water content controls the maximum amount of biomass combusted during a fire (Soja et al., 2004; French et al., 2004; Turquety et al., 2007), carbon release from organic soils is strongly affected by seasonal variations in water-table height. To depict this in the model, we added rules that scaled C output from burning to seasonal changes in moisture. The seasonal pattern of soil moisture is fairly predictable in Interior Alaska. Fires occurring early in the summer typically encounter relatively high water tables in areas of thick organic soils, while the rare late summer fires usually encounter low water tables.
Thus for the 40% of the annual area burned that is underlain by thick organic deposits (peat), if a fire occurs before July 1:
where 0.7 is a scaling factor applying to burning at times of high water tables. If a fire occurs between July 1 and August 1:
where 1.0 is the scaling factor depicting the “average” height of the water table. If a fire occurs after August 1, conditions are deemed unusually dry, so:
where 1.5 is the scaling factor accounting for unusually low water tables in organic soils. AU the above estimates are divided by 2 to convert dry matter to C. For the other 60% of the landscape that burns and is not peatland, we assumed that 3.7 kg of Dry Matter m-2 (1.85 kg of C m-2) is combusted regardless of the forest successional stage or the preceding fire history.
Model Spin-Up
To establish realistic initial conditions for the model, we performed spin-ups that were > 10 × the length of a typical fire cycle in Interior Alaska, which is estimated at 105 years since A.D. 1920 (Kasischke et al., 2010). For the first 1000 years of each spin-up, pairs of temperature and precipitation values were chosen at random from the CRU-PIK-derived maps of monthly mean temperature and precipitation during summers between A.D. 1860 and 1990. This was followed by a 392-year additional spin-up that used 1860–1990 repeated in order three times. To check for possible “sequencing” artifacts, we compared these results against the result of a single 392-year, randomized climate scenario drawn from the same climate data and found no difference. The resultant standage, vegetation, and burn-severity maps were then used as the model's depiction of the region's forested landscape in the year A.D. 1860. To assess the model's sensitivity to initial conditions, we made multiple model runs using a variety of flammability coefficients. Spin-up under conditions of heightened flammability resulted in an 1860 forest with approximately equal spruce and deciduous trees. Reducing flammability resulted in more spruce trees. Regardless of initial flammability, the modeled vegetation converged towards the spruce/deciduous ratio actually observed in the region in 2001 (See Results section), demonstrating that ALFRESCO–s behavior after A.D. 1860 is quite robust with respect to initial conditions.
Model Verification
After calibrating the model to agree with observed and backcast estimates of annual area burned, we assessed its performance in two ways. First, we used ALFRESCO to predict the age structure of the present-day forest and then compared this prediction to the age structure inferred from 4800 tree ages obtained along an eastwest transect across the region (Mann, Rupp, and Duffy, unpublished data). The close agreement of observed and predicted age structures demonstrates that the model accurately simulates the interactions among climate, fire, and post-fire succession that underlie the age structure of the present-day forest (Fig. 2).
The second way we checked the model's realism was by comparing the exponents of the power-law relationships between frequency and magnitude of annual area burned. As mentioned above, power-law exponents, β, are a sensitive metric for characterizing fire regimes across multiple spatial and temporal scales. ALFRESCO produces estimates of β that are statistically similar to those estimated from the observed and back-cast data between 1988 and 2007 (Fig. 3). We did not use the observed area-burned data collected prior to 1988 in these calculations because of irregularities in the way fire records were kept in those years (Lyons et al., 2008).
Results
FIRE REGIME
ALFRESCO predicts that the annual area burned in Interior Alaska will increase markedly after 2010 compared to the annual areas burned between 1860 and 2010 that are estimated from historical records and from our back-casts of annual area burned (Fig. 4). This steepening occurs under all five GCMs we used as predictors of future climate. In the model, a region-wide increase in burning started around 1990 (Fig. 4), and is predicted to continue until at least 2040. Much of the increase in annual area burned is caused by more frequent mega-fire seasons like those occurring in 2004 and 2005 when > 15,000 km2 burned each summer.
As annual area burned increases, other aspects of the fire regime are predicted to change as well. Mean fire size will increase, continuing a trend that started in the 1990s (Fig. 5, part a). An increased frequency of sporadic, mega-fire summers like those in 2004 and 2005 will probably typify Interior Alaska's fire regime over the next 40 years. ALFRESCO also predicts that the variance in the annual area burned will increase until ca. 2020 when the forest begins to stabilize under its new fire regime (Fig. 5, part b).
FOREST STRUCTURE AND COMPOSITION
According to the ALFRESCO model, climate-driven changes in the fire regime are causing striking changes in the structure of Alaska's forests already today, and these changes are predicted to continue over the next 30 years. Increased burning is removing old-aged stands, which is causing a decline in both the mean stand age and the spatial variability of stand ages across the landscape (Fig. 6). In 2001, white and black spruce together comprised more than half of all forest stands in Interior Alaska (Homer et al., 2004). Our retrospective modeling suggests spruce were even more abundant between 1920 and 1990 when they comprised roughly two-thirds of all stands (Fig. 7). Looking ahead, ALFRESCO predicts that changes in the fire regime will cause deciduous plants (aspen, birch, willow, and herbaceous taxa) to replace spruce at many sites. In the model, the magnitude of fire-conducive climatic changes predicted by the GCMs overwhelms the potential, negative feedback between the increased spatial coverage of less flammable, deciduous vegetation and annual area burned. The replacement of spruce by deciduous species is predicted to slow by ca. 2020, by which time deciduous vegetation cover is predicted to be nearly twice as abundant as spruce.
ENERGY FLUXES
Because ALFRESCO is spatially explicit and tracks fire history and successional stages, we can quantify the biophysical effects of the predicted changes in fire regime and forest composition across the entire shifting mosaic of fires and post-fire vegetation that covers Interior Alaska (Fig. 8). The predicted net increase in summertime albedo caused by the proliferation of deciduous vegetation at the expense of spruce across Interior Alaska could lower net solar radiation in summer by 3–5 W m-2 by 2020 (Fig. 8, part a), a decline of 6–8% from present values. This albedo modification is comparable in magnitude (though opposite in sign) to the estimated effects of earlier snowmelt in springtime (Euskirchen et al., 2009). In response to changes in summertime albedo,ALFRESCO predicts that sensible heat fluxes will decline (Fig. 8, part b). Note that these results do not address how changes in vegetation cover could affect snow-season albedos and energy fluxes (Euskirchen et al., 2010). In contrast to net radiation, summertime fluxes of latent heat are predicted to increase because deciduous leaves have higher evapotranspiration rates than conifer needles (Beringer et al., 2005; Amiro et al., 2006) (Fig. 8, part c). An increase in burned-over, blackened ground will likely cause a transient increase of 3–6% in heat flux into the ground in summer (Fig. 8, part d).
CARBON EMISSIONS
The altered fire regime and vegetation cover that ALFRESCO predicts could have dramatic effects on the carbon balance of Interior Alaska's forest. Burning releases an immediate pulse of C to the atmosphere, followed by a decade of declining C emission from decaying vegetation (Welp et al., 2006) (Fig. 9). C sequestration occurs most rapidly in the first 10–40 years after a fire during the initial stages of revegetation (O'Neill et al., 2003; Welp et al., 2006). By integrating C emission and sequestration across the disturbance/recovery mosaic of the entire forest, ALFRESCO is able to estimate the net impacts of changes in fire regime and vegetation composition on C balance (Fig. 9). Results suggest that starting ca. 1990, Interior Alaska became a net source of C to the atmosphere. This was accentuated in the mega-fire years of 2004 and 2005, and Interior Alaska is predicted to continue to be a C source until sometime after 2020 (Fig. 9). This predicted 25- to 50-year pulse of heightened C emissions results from the release of C that accumulated in Interior Alaska's vegetation and soils over preceding centuries when the climate was less favorable for burning.
Discussion
These modeling results support previous inferences based on assessments of ecological sensitivities (Chapin et al., 2004) and on tree-growth responses (Beck et al., 2011) that a critical ecological threshold is now being crossed in Interior Alaska's boreal forest. The key processes that have created this threshold event can be inferred from the interactions between fire, forest succession, and climate change depicted in the ALFRESCO model. As previously suggested by Duffy et al. (2005) and Kasischke et al. (2010) and others, the first and largest impact of changing climate in Interior Alaska is occurring through the climate-fire linkage. This is happening because slight changes in summer weather are capable of triggering large changes in the fire regime. Specifically, annual area burned increases nonlinearly as summer weather becomes warmer and drier.
Another important threshold that emerges from the ALFRESCO model is the outcome of interactions between fire frequency and post-fire succession. Increasing annual area burned shortens the temporal window available for post-fire succession, which favors plant species with shorter juvenile periods. Once the time between successive fires drops below a threshold of 60–80 years, spruce trees are increasingly excluded from this forest. Fires intervene before the slower growing, longer-lived spruce can replace the faster growing deciduous species in the forest canopy (Kurkowski et al., 2008).
The impacts of both the climate-area burned relationship and its domino effect on forest succession are accentuated by the ecological legacy of fuels accumulated under preceding fire regimes that were characterized by less frequent and smaller fires. The more fuel that is available, the greater is the surge in burning that occurs as climate warms, and the less likely it is that a spruce canopy will have time to develop. According to the ALFRESCO model, the ecological threshold now being crossed in Interior Alaska is the result of a fuel-load legacy encountering rapid summertime warming, which together are triggering the nonlinear climate-fire and fire-succession processes just described.
Why doesn't a negative feedback kick in to reduce flammability as deciduous vegetation becomes more widespread on the landscape? As every firefighter knows, deciduous plants tend to be less flammable than spruce trees. In the ALFRESCO model, the scenarios of future climate are so warm and dry that they overwhelm the potential negative feedback caused by increases in deciduous vegetation cover. This result is not surprising given that much of the increase in burning is predicted to occur during extreme fire seasons like those of 2004 and 2005 when fires were able to burn spruce and deciduous stands at similar frequencies (Kasischke et al., 2010).
WARNING SIGNS OF IMPENDING THRESHOLDS?
Our results identify several warning signs that might prove useful for detecting ecological regime shifts elsewhere. The predicted shift to a deciduous state in Interior Alaska was preceded by decades of accumulation of old (> 100 years) conifer stands on the landscape. This created high levels of spatial variation in stand ages (Fig. 6, part b). Because of their heightened flammabilities, the presence of these old stands predisposed the forest to extensive fires and lowered its overall resistance (the capacity of a system to absorb stresses and continue functioning) to climate-induced perturbations to its fire regime. When combined with the ongoing climate trend, the 1990 age structure of this forest was unstable because annual area burned inevitably increased as summers became warmer and drier. Old spruce stands were first to burn because of their higher flammability, and their loss increased the coverage of younger stands, which caused the spatial variability of stand ages to decline. Based on these results, we speculate that in ecosystems where the probability of stand-replacing disturbances increases with stand age, heightened variability in stand age is a warning sign of lowered resistance and an impending threshold.
Increasing variance sometimes precedes ecological transitions (Carpenter and Brock, 2006; Contamin and Ellison, 2009). In Interior Alaska, interannual variance in area burned appears to be peaking now as the forest shifts towards a new, deciduous-dominated state (Fig. 5, part b). This is happening because fires characteristic of both the old and the new fire regime occur. In effect, the forest is now straddling two different attractor basins as it transits into a new disturbance regime.
EXISTING ANALOGS TO THE FUTURE FOREST?
The future forest predicted by the ALFRESCO model has both prehistoric and present-day analogs. A poplar-dominated, parkland vegetation was widespread in Interior Alaska during the early Holocene (Higuera et al., 2009) when a precession-driven summer insolation anomaly caused summertime temperatures 1–2 °C higher than today (Kaufman et al., 2004).
A present-day analogue for Alaska's future forests may be the boreal mixedwoods of south-central Canada (Hogg and Hurdle, 1995; Cumming, 2001; Schneider et al., 2009) that grow under a climate slightly warmer and drier than Interior Alaska is today. Aspen dominates the boreal mixedwoods, white spruce is restricted to mesic sites (Alberta Natural Regions Committee, 2006), and lightning-caused fires are common (Weir et al., 2000).
How close is the analogy between Alberta's mixedwoods and Interior Alaska's future forests? To address this question we compared the power-law exponents describing the frequency-magnitude relationships of predicted fires in Interior Alaska and recent fires in the boreal mixedwoods (Fig. 10). ALFRESCO predicts that the β exponent declines rapidly in Interior Alaska after 1990, falling to values between 1.4 and 1.5 by 2030. In the mixedwoods of northern Alberta between 1961 and 2007, β was 1.48 (r2 = 0.99) for all fires >0.5 km2 in size (Government of Alberta Sustainable Resources Development Wildfire Information, 2009). Alaska's declining β is the result of an increase in the relative frequency of larger fires, which is consistent with both a predicted increase in fire size (Fig. 5) and with the trajectory of the historical fire data (Kasischke et al., 2010). The close match between the frequencymagnitude relationship of fires in Alberta's boreal mixedwoods between 1961 and 2007 and that predicted by ALFRESCO for Interior Alaska between 2010 and 2040 supports the idea that Interior Alaska's forest and its fire regime are changing to resemble Alberta's present mixedwoods.
LOCAL IMPLICATIONS
Besides the prospect of more mega-fire seasons like those in 2004 and 2005, one of the most troubling implications of our modeling results concerns the predicted increase of heat flux into the ground (Fig. 8, part d). The transient spike that is predicted in ground heat flux during summer implies an increase in thermokarsting throughout the region. Thermokarst is the alteration of existing topography by the thawing of underlying permafrost. In Interior Alaska, thermokarsting often occurs after fires because burning removes organic matter that previously shaded or insulated the ground (Jorgenson and Osterkamp, 2005; Jorgenson et al., 2010). Thermokarst may continue for decades to a century after a fire (Viereck et al., 2008) and it can increase trace gas emissions by thawing organic-rich sediment and disrupting stream networks (Schuur et al., 2008, 2009). Thermokarsting can also enhance Hgmethylation rates by creating waterlogged conditions (Grigal, 2003; Turetsky et al., 2006).
GLOBAL IMPLICATIONS
Our results suggest that Interior Alaska switched from being a C sink to a C source starting ca. 1990 accompanied by a 20- to 30-year pulse of heightened C emissions (Fig. 9). These findings are consistent with the results of other studies using different models that suggest boreal North America shifted from being a carbon sink to a source at the end of the 20th century (Balshi et al., 2007, 2009; McGuire et al., 2009). If boreal forests in other parts of the subarctic undergo similarly radical shifts in their fire and successional processes, a significant positive effect on global warming seems likely, especially if these changes in fire regimes and vegetation composition are accompanied by decreases in annually averaged albedo caused by lengthening of the snow-free season (Euskirchen et al., 2010). We have resisted applying our regionally summed estimates of future C emissions, albedo, and heat fluxes to estimating Alaska's contributions to global radiative forcing. While interesting, these back-of-the-envelope calculations have enormous error terms. For example, Randerson et al. (2006) estimated the total forcing derived from a fire in Interior Alaska to be 18 ± 42 W m-2 in the year following the fire and -2.4 ± 2.3 W m-2 during years 0–80 post-fire.
HOW GOOD IS ALFRESCO?
Models are deliberate abstractions of reality. A good model is one that quantifies what we think we know about how a natural system works in the simplest possible manner leading to meaningful results. Meaningful results either expose major gaps in our current understanding and/or make plausible predictions about what could happen to the system in the future. One can add realistic embellishments to a model ad infinitum without improving on these two main purposes.
The ALFRESCO model has several strengths and several weaknesses. Its greatest strength is its ability to synthesize the ecological states of myriads of separate, post-fire successional disturbance patches across a vegetation mosaic of regional extent. ALFRESCO's other strength is its ability to depict climate-firevegetation interactions, including feedbacks caused by the vegetation mosaic's changing flammability due to successional stage.
ALFRESCO has two possible weaknesses at present. The first one involves the stationarity of the relationship between summer climate and annual area burned as deduced by Duffy et al. (2005). This relationship may change in the future as the region's vegetation and climate both change. Twenty years hence it is possible that snowmelt may have advanced a month, June mean temperature may be significantly warmer, and the onset of summer rains may routinely be delayed into September. Sizable changes in the climate of Interior Alaska are indicated by tree-ring studies that cover the last 200 years (Barber et al., 2004). Such changes could alter how ALFRESCO now predicts annual area burned from climate parameters.
ALFRESCO's second potential weakness involves the way it depicts post-fire succession. The model now assumes that all vegetated sites, regardless of their pre-fire vegetation cover, undergo “species-dominance relay” in which all plants establish within a few decades after a fire. The faster growing, deciduous species then dominate the canopy until either the next fire destroys them or spruce trees slowly emerge from the understory. Recent work by Kurkowski et al. (2008) and Johnstone et al. (2010) shows that another successional pathway, self-replacement, is important in certain settings. One such setting is hilly terrain where low summer sun angles cause steep gradients in incoming solar radiation and so create persistent safe sites for species possessing different physiological tolerances (Kurkowski et al., 2008). Another setting where self-replacement may be important is peatland where fires are unable to alter the pre-existing soil and drainage conditions, and thus black spruce communities immediately reestablish (Johnstone et al., 2010). In these two settings, it could be that more frequent and severe fires will have little impact on the landscape's overall conifer/deciduous ratio. It is possible that incorporating a self-replacement pathway of post-fire succession into the ALFRESCO model will reveal a larger inertia in the response of vegetation to climate change. This might reduce the speed and magnitude of the ecological changes predicted here.
Conclusions
Results of the ALFRESCO forest-disturbance model support previous suggestions that a major ecological regime shift is now underway in Interior Alaska. The formerly spruce-dominated vegetation seems to be responding to rapid climate warming by shifting back towards its early Holocene state, which was similar to the boreal mixedwoods now growing in southern Canada. Interior Alaska's forest started across this ecological threshold ca. 1990, and the transition is predicted to be completed by 2040. Our modeling results suggest that larger and more frequent wildland fires are an important driver of this ecological regime shift. The changing fire regime is the result of the intersection of warming summer climate with the legacy of fuel accumulated over preceding decades of less frequent fires. Crossing the ecological threshold between conifer forest and boreal mixedwoods will have large impacts on the biophysical properties of Interior Alaska, including its heat fluxes and carbon budget, with possible positive feedbacks to global warming.
Acknowledgments
Partial funding came from National Science Foundation grant ARC-0902169. We thank Randy Jandt, Jane Wolken, Chris Fastie, Pieter Beck, Glenn Juday, and Janet Wilmshurst for useful discussions in the course of this project.