With climatic warming, wildfire occurrence is increasing in the boreal forest of interior Alaska. Loss of catchment vegetation during fire can impact streams directly through altered solute and debris inputs and changed light and temperature regimes. Over longer time scales, fire can accelerate permafrost degradation, altering catchment hydrology and stream nutrient dynamics. In 2004, the 217,000 ha Boundary Fire burned 65% of an established study site in the Caribou-Poker Creeks Research Watershed. We used this opportunity to investigate the impact of wildfire on stream chemistry and metabolism in boreal forest catchments. Wildfire impacts on chemistry were evaluated by examining solute chemistry in four catchments from 2002 to 2007. Ecosystem metabolism was measured over the summer of 2005 in one burned and two unburned catchments. Wildfire led to stream nitrate concentration increasing up to threefold, whereas dissolved organic carbon (DOC) and dissolved organic nitrogen concentrations decreased post-fire. Average stream gross primary production in the burned catchment was double that of the unburned sites (2.4 and 1.2 g O2 m−2 day−1, respectively). Respiration rate was also elevated in the burned stream (6.6 g O2 m−2 day−1) compared with the control streams (1.2 and 4.5 g O2 m−2 day−1). Climatic warming has the potential to impact boreal forest streams through permafrost thaw and increased fire frequency, leading to altered solute inputs and production and respiration rates.
In the boreal forest, wildfire is a widespread disturbance and fundamental in shaping forest structure, function, and successional processes (Chapin et al., 2006). With climatic warming, the northern boreal forest is predicted to experience earlier onset of the wildfire season combined with a potential increase in the number and size of fires (Kasischke et al., 1995; Stocks et al., 1998). In regions with permafrost, the effects of fire on surface organic layer depth, albedo, and vegetation are likely to initiate a positive feedback on soil temperature, contributing to decreased permafrost extent and increased depth of seasonal thaw (Burn, 1998; Yoshikawa et al., 2003).
Wildfire has the potential to impact boreal forest streams over short time scales via change in catchment vegetation and soil characteristics, and over longer time scales through reduced permafrost extent, altered vegetation successional trajectories, and changes in hydrological pathways. Immediate impacts of fire on stream ecosystems include changes in temperature and light regimes through loss of shading riparian canopy (Gresswell, 1999; Hitt, 2003). Streams are closely coupled to their catchments with allochthonous inputs serving as an important source of organic matter and nutrients. Combustion of terrestrial organic matter may alter delivery of solutes to streams, in particular carbon, nitrogen, phosphorus, and base cations (Bisson et al., 2003; Gresswell, 1999). In addition, loss of riparian vegetation can alter inputs of coarse and fine particulate organic matter and woody debris to streams (Benfield, 1997). A reduced rate of catchment evapotranspiration following wildfire may also lead to changes in stream discharge and hydrologic flowpaths through watersheds, thereby altering both the magnitude and timing of flows and solute delivery to streams (Earl and Blinn, 2003; Petrone et al., 2007; Williams and Melack, 1997).
Vegetation stand type and permafrost extent are intimately linked and generally maintained through negative feedbacks in boreal forest watersheds. Soils underlying black spruce stands are generally wetter and more acidic than soil underlying deciduous stands (Baldocchi et al., 2000; Ping et al., 2005), which keeps organic soil mineralization low, and promotes maintenance of a thick organic layer and underlying permafrost. When fire in black spruce forest is accompanied by removal of organic soils and disruption of the soil thermal regime, it can lead to establishment of deciduous stands, which is most likely to occur in areas where permafrost is unstable (Chapin et al., 2004). Over successional time scales, change in vegetation type and the loss of underlying permafrost have potential to alter soil nutrient chemistry, soil active layer depth, and catchment hydrologic flowpaths, with associated change in allochthonous inputs to streams.
Alteration of stream physical and chemical properties by fire will likely impact stream ecosystem functioning. Whole stream metabolism is a measure of instream primary production and ecosystem respiration and, as such, is an integrative measure of stream function. The rate of metabolism can have significant ramifications for nutrient uptake in streams, downstream transport of nutrients, and stream food webs (Amon, 2002; Bunn et al., 1999; Hall and Meyer, 1998; Hall and Tank, 2003; Thorp, 2002). Temperature, light, dissolved organic matter, nutrient availability, and discharge are important determinants of production and respiration rates in streams (Elser et al., 1990; Elser et al., 2000; Gausch et al., 1995; Uehlinger, 2000). Post-fire alteration of discharge regime, canopy loss, and nutrient enrichment have been shown to alter invertebrate assemblages and stream food webs (Minshall et al., 1997; Spencer et al., 2003; Vieira et al., 2004).
The objectives of this research were twofold. First, we examined the impact of wildfire on stream nutrient and dissolved organic matter (DOM) chemistry over a six-year period. Second, we examined the impacts of wildfire on stream functioning by measuring whole stream metabolism in one stream draining a fire impacted catchment and in two control streams. We hypothesized that wildfire in the surrounding catchment would lead to increased inorganic nutrient input to streams, thereby alleviating any nutrient limitation and resulting in increased rate of GPP. In contrast, we hypothesized that combustion of soil carbon stocks in the surrounding catchment and reduced allochthonous DOM input to streams would result in a lower rate of stream ecosystem respiration.
Our study was conducted in the Caribou-Poker Creeks Research Watershed (CPCRW; 65.15°N, 147.5°W), located approximately 50 km NE of Fairbanks, Alaska. The CPCRW is a pristine, 104-km2 catchment reserved for hydrological and ecological research. The climate of CPCRW is continental, with warm summers (mean = 16.4 °C in July), cold winters (mean = −29 °C in January), and low precipitation (411 mm, of which 31% falls as snow). The research site is located in the zone of discontinuous permafrost, with the extent of permafrost underlying sub-catchments of the CPCRW ranging from 3 to 53%. The distribution of permafrost is largely determined by aspect and winter temperature, with permafrost found predominately on north-facing slopes and in valley bottoms (Haugen et al., 1982). Uplands on south-facing slopes are dominated by hardwood forests of Alaskan paper birch (Betula neoalaskana) and quaking aspen (Populus tremuloides), whereas north-facing slopes generally have black spruce (Picea mariana), and feathermosses (Pleurozium schreberi and Hylocomium schreberi). Several alder species have a patchy distribution in both uplands (Alnus viridis) and valley bottoms (A. incana). Valley bottoms typically have saturated soil and support growth of mosses (e.g., Sphagnum spp., Hylocomium) and dwarf shrubs (e.g., Betula nana, Salix spp., Vaccinium uliginosum).
We studied a total of four streams that drained catchments that varied in size from 5.2 to 10.0 km2 and permafrost extents from 4 to 26% (C1, C2, C4, and P6 catchments; Table 1, Fig. 1). Within the CPCRW, the C4 catchment was partially burned in 1999 by a prescribed fire to assess fire and climate feedbacks in the boreal forest (FROSTFIRE project; Hinzman et al., 2003). The prescribed fire was of moderate intensity, burned 28% of the catchment area, and was largely restricted to the black spruce–dominated north-facing slopes with the riparian zone left mostly unburned. Stream nitrogen was elevated in C4 during a single post-fire storm flow, after which stream nutrient concentrations were not significantly different from pre-fire levels (Petrone et al., 2007).
Selected physical characteristics of study catchments at the Caribou-Poker Creeks Research Watershed. Area burned from D. Verbyla (personal communication). Temperature, conductivity, pH, and discharge are summer (June–September) averages (±SE). Riparian canopy is classified as open, closed, or a combination of open and closed (mixed), based on visual assessment.
The Boundary Fire burned approximately 217,000 ha of boreal forest in interior Alaska between mid-June and late August 2004. Within the CPCRW, the Boundary Fire burned to the south of Poker Creek and much of the P6 catchment (Fig. 1). The fire was of moderate to severe intensity and burned approximately 65% of the P6 catchment. Burning was primarily in upslope, black spruce–dominated stands. Following the fire, organic soil was largely intact and there was considerable charred litter and dead vegetation remaining. The riparian corridor, which had a width of approximately 40 m, was largely untouched by the fire. Precipitation was unusually low in 2004 and gradually increased from 2005 through 2007 (Fig. 2).
The impacts of fire on stream chemistry were assessed by comparing pre- and post-fire stream chemistry in two paired streams (C1 and P6). P6 was the only catchment in the CPCRW with comprehensive pre-fire data that was burned during the wildfire. The C1 stream was selected as the control based on the similarity in catchment size, permafrost extent, dominant vegetation, and similarity in response of stream chemistry with changes in discharge (Table 1). Although the P6 and C4 watersheds had essentially the same extent of permafrost (18 and 19%, respectively), the streams have substantially different responses in streamflow and chemistry during storms. In addition to monitoring stream chemistry, we measured whole-stream metabolism in three streams (C2, C4, and P6) during summer 2005 to evaluate nutrient controls over instream respiration and gross primary production (GPP). All studies were conducted during the ice-free period when streams are free-flowing (typically June to September).
STREAM CHEMISTRY AND PHYSICAL MEASUREMENTS
Stream stage height was measured continuously every 15 minutes using pressure transducers at Parshall flumes at the C2, C4, and P6 streams. Rating curves were developed for all sites from velocity and stream cross-section measurements. Water chemistry was measured throughout the summers of 2002–2007 in the C1, C2, and C4 streams. The P6 stream was sampled from June to September in all years, except 2004 when water samples were collected from July to September (post–Boundary Fire). Stream water samples were collected daily using ISCO autosamplers and an additional grab sample was collected every 1 to 2 weeks. Previous analyses at CPCRW have shown that for most solutes, including nitrate and DOC, the concentration measured in autosampler-collected samples (without use of preservatives) stored for two weeks is not significantly different from grab sample concentrations (r2 = 0.95; slope = 0.99; Jones et al., 2005). Samples were collected in acid washed high density polyethylene bottles, stored in a cooler for transport to the lab and filtered within 12 h using glass fiber filters (Gelman A/E). Samples were stored at 4 °C until analysis (usually within 48 h). When samples could not be analyzed immediately, samples were filtered using glass fiber filters (Gelman A/E) and frozen. Several sample dates were missed at each stream due to mechanical failure of the autosamplers.
In-stream GPP, respiration, and net ecosystem productivity (NEP) were calculated from the diel change in dissolved oxygen concentration using the single-station, whole-stream technique (Bott, 1996; Odum, 1956). Metabolism was measured over 24 h approximately biweekly in 2005 from June through August in C2 and C4 streams, and approximately biweekly from July through August in the P6 stream. Metabolism was not measured in the P6 stream in June due to equipment failures. Because interior Alaska experiences up to 20 h of daylight during the summer, we calculated daily respiration from the minimum hourly change in dissolved oxygen over the day. On each date, dissolved oxygen was measured every 15 minutes using CS511–L Sensorex dissolved oxygen probes connected to Campbell Scientific CR10X dataloggers (Campbell Scientific). Dissolved oxygen probes were calibrated on each metabolism measurement date using a spectrophotometric dissolved oxygen method (Roland et al., 1999). Triplicate 250-mL water samples were collected in BOD bottles from the thalweg of each stream, and immediately fixed in the field with 2.5 mL MnSO4 and 2.5 mL KOH–KI. Stream temperature probes were calibrated at the time of metabolism measurements using a high-precision thermometer. To express metabolism on an areal basis, stream widths were measured once at all sites in 2005. All of the study streams have incised stream channels, and, as a consequence, width does not vary appreciably with discharge. Width was measured at 5 m intervals over a 200 m reach and averaged.
Oxygen exchange with the atmosphere was determined using the steady-state injection technique (Hibbs et al., 1998; Mulholland et al., 2001). We conducted eight injections across a range of flows to develop a model relating evasion rate to streamflow. Evasion was quantified in the C2 and C4 streams, which have the lowest and highest streamflows of the streams studied. In the lab prior to injections, SF6 was injected into collapsible containers of conservative tracer (dissolved Cl− or Br−) and allowed to equilibrate for 24 h. The resulting solution was injected into the study stream 300 m upstream from the dissolved oxygen probe using a metering pump (flow rate = 50 mL min−1) until steady-state of the conservative tracer was achieved at the downstream station (approximately 3 h). Prior to each injection, triplicate water and triplicate gas samples (30 mL stream water collected in 60 mL polyethylene syringes) were collected from the thalweg at sampling stations every 50 m downstream from the injection site (300 m total reach length).
Stream discharge (Q) at the time of injection was determined using the dilution gauging technique and calculated asStream Solute Workshop, 1990). Stream velocity (v) was calculated as t0 is the injection start time, and t1 is the time when solute concentration at the downstream station reached half plateau concentration. Stream depth (z) was calculated by dividing discharge by velocity and average stream width. τ is the hydrologic travel time between the up- and downstream sampling points, G is the steady-state tracer gas concentration, C is the conservative tracer concentration, and the subscripts 1 and 2 are up- and downstream sampling locations, respectively (Generaux and Hemond, 1992). SF6 evasion was converted to O2 evasion using the equation is evasion rate at ambient stream temperature (TS; Elmore and West, 1961), and are the Schmidt numbers for SF6 are O2, respectively (Wanninkhof, 1992), and n = 0.7 (Generaux and Hemond, 1992). Oxygen evasion was standardized to evasion at 20 °C () using the formula
To develop an evasion model for CPCRW streams, we compared evasion against velocity and depth via regression analysis using the general form of the equation by Bennett and Rathburn (1972)
Anions (NO3−, SO42−) and cations (NH4+, Ca2+, Mg2+, Na+, K+) were quantified using a Dionex DX-320 Ion Chromatograph. Soluble reactive phosphorus (SRP) was measured using the molybdate blue method using a Beckman DU 640B spectrophotometer (American Public Health Association, 1998). DOC concentration was measured using a Shimadzu 5000 total organic carbon analyzer, which was plumbed to an Antek 7050 nitric oxide chemoluminescent detector to quantify total dissolved nitrogen (TDN). Dissolved organic nitrogen (DON) was calculated as TDN − (NH4+ + NO3−). Dissolved oxygen samples were acidified with 2.5 mL 12N H3PO4 and the light absorbance at 430 nm read using a Beckman DU 640B spectrophotometer according to the method described by Roland et al. (1999). SF6 was measured by allowing syringes to warm to room temperature, adding a 30 mL helium headspace, and equilibrating for 2 h. The headspace was then injected into a Varian CP 3800 gas chromatograph (1 mL) and measured with an electron capture detector.
Differences in stream metabolism among streams were tested using one-way ANOVA (α = 0.05). The effects of fire on stream chemistry were tested in two manners. First, we used a two-way ANOVA to test for significant differences in chemistry among streams, and prior and post-fire. The interaction term, in particular, was evaluated to determine if the degree of change in stream chemistry before and after the fire in the burned watershed differed from the other streams. Second, we further evaluated interannual patterns in stream chemistry using a one-way ANOVA with the differences in stream chemistry between the stream draining the burned watershed and a control stream. This second approach was used to maximize our ability to detect interannual changes following the fire. The control was selected based on the greatest similarity in the daily change in stream chemistry between the control and treatment stream, which was the C1 stream. Significant differences among streams (p < 0.05) were further evaluated using Tukey's multiple-comparison test.
STREAM SOLUTE CHEMISTRY
Stream nitrate concentration varied significantly among CPCRW streams (Table 2). Prior to the fire, the mean nitrate concentration in the streams in the lower permafrost watersheds was 0.51 and 0.59 mg N L−1 (C2 and C4 streams, respectively) compared with 0.24 and 0.28 mg N L−1 in the streams draining the higher permafrost watersheds (C1 and P6, respectively; Table 2). Stream water DOC concentration, in contrast, exhibited the opposite pattern with greatest concentration in the streams draining the higher permafrost watersheds. Prior to the fire, the mean DOC concentration in streams in the lower permafrost watersheds was 3.3 and 3.4 mg C L−1 (C2 and C4 streams, respectively) compared with 4.7 and 5.7 mg C L−1 in streams draining the higher permafrost watersheds (C1 and P6, respectively). While the mean DOC concentration among streams varied nearly twofold, DON concentration prior to the fire only varied from 0.28 to 0.32 mg N L−1 (Table 2). Coupled to the patterns in DOC, the C:N of dissolved organic matter tended to be greatest in higher permafrost streams with values of 20.8 in P6, 19.0 in C1, 15.1 in C4 and 13.2 in C2 streams. Similar to the pattern in nitrate, sulfate and base cation concentrations tended to be greatest in the streams draining the lower permafrost watersheds for streams in the Caribou Creek sub-watershed (Table 2). In the P6 stream, sulfate and most base cation concentrations were always equal to or greater than in the tributaries of Caribou Creek, suggesting a difference in parent geology (Table 2). SRP concentration was almost always below detection limits (<1 µg L−1).
Mean summer (June–September) stream solute concentrations (±SE) prior and post fire for the four study streams. Superscript letters indicate significantly different concentrations of solute.
WILDFIRE IMPACTS ON STREAM CHEMISTRY
Following the wildfire, nitrate, sulfate, and potassium concentrations increased in the impacted stream, whereas DOC and calcium concentrations declined relative to the three undisturbed streams (Tables 2 and 3). Wildfire had a dramatic effect on stream nitrate concentration (Fig. 3, Tables 2 and 3). During 2002–2003, nitrate concentration was on average only 0.08 mg N L−1 lower in the P6 (burned) stream compared with the C1 (control) stream. Nitrate concentration post-fire in 2004 did not significantly change in the P6 stream compared with the control, although 2004 was a dry summer with little rainfall (Fig. 2) and low streamflow. In subsequent summers following the fire (2005–2007), nitrate concentration in the P6 stream increased substantially compared with the pre-fire difference between burned and control streams (mean difference 0.30 mg N L−1; Table 3, Fig. 3). Post-fire, nitrate concentration in the P6 stream was on average over twofold greater with 0.65 mg N L−1 compared with 0.28 mg N L−1 in years 2002–2004 (Table 3, Fig. 3).
Mean summer (June–September) stream solute concentrations and the mean difference in concentrations between streams (±SE) for the C1 (control) and P6 (burned) sub-catchments of the Caribou-Poker Creeks Research Watershed from 2002 to 2007. For differences in solute concentrations, superscript letters indicate years with significantly different differences between the C1 and P6 streams.
DOC and DON exhibited the opposite pattern to nitrate with the concentration declining post-fire. Prior to the fire the DON concentration in the P6 stream was on average 0.05 mg N L−1 greater than in the control stream (2002 and 2003; Table 3, Fig. 4). Post-fire, DON concentration in the burned stream declined compared with the control to the point that DON concentration in the control stream exceeded the burned stream in 2005 (Table 3, Fig. 4). Similar to the pattern in DON, DOC concentration declined in the burned stream post-fire, although the response was delayed until the wetter years of 2006 and 2007. From 2002 to 2005, the DOC concentration in P6 stream was on average 1.2 mg C L−1 greater than in the control stream (Table 3, Fig. 5). In 2006 and 2007, however, the DOC concentration in the P6 stream significantly declined relative to the control stream to the point at which the burned stream had on average 0.5 mg C L−1 less than the unburned stream (Table 3). Coupled to the change in DOC and DON concentrations, the C:N of DOM significantly declined in the burned stream compared with the reference. Prior to the fire, the C:N of the DOM in the P6 stream averaged 22.9 (2002 and 2003), but tended to decline post-fire (2004, 2006, 2007).
Gas evasion rate was quantified across streamflows ranging from 19 to 60 L s−1 and velocities from 0.13 to 0.24 m s−1. Using our field-based measurements of gas evasion and streamflow with a generalized form of the Bennett and Rathburn (1972) model, the equation that best described oxygen evasion was (R2 = 0.72)
GPP and respiration were greatest in the burn-impacted P6 stream compared with the control C2 and C4 streams (p < 0.05; Fig. 6). GPP in the burned stream was double that of the unburned streams with a summer mean of 2.4 g O2 m−2 d−1 compared with an average rate of 1.2 g O2 m−2 d−1 (Fig. 6). Similarly, respiration rate was higher in the P6 stream with a summer average of 6.6 g O2 m−2 d−1, whereas average respiration rate in the unburned C2 and C4 streams was 1.2 and 4.5 g O2 m−2 d−1, respectively (Fig. 6). GPP and respiration were not statistically correlated with stream temperature. However, GPP may be related to canopy cover. With a closed riparian canopy, the C4 stream had the lowest GPP (average 0.9 g O2 m−2 d−1) whereas C2, with an open canopy, had a higher average GPP of 1.5 g O2 m−2 d−1. P6, however, had a mixed canopy cover that was undamaged by the wildfire and substantially higher production than either C2 or C4 (Figure 6). Similarly, GPP and respiration were unrelated to stream nitrate or DOC concentrations.
CPCRW STREAM CHEMISTRY
In this study, stream solute concentrations varied widely among sub-catchments and tended to vary in response to extent of permafrost (MacLean et al., 1999), although other factors such as permafrost distribution, catchment size, mineral geology, and vegetation composition also undoubtedly affect solute chemistry. Permafrost serves as a barrier to the infiltration of precipitation into deep soil horizons and is therefore an important regulator of stream solute concentrations in the region of discontinuous permafrost (Carey, 2003; Jones et al., 2005; MacLean et al., 1999; Petrone et al., 2006). In catchments underlain by extensive permafrost, hydrological flowpaths are restricted to the upper organic soil horizon, with streams receiving greater DOM inputs. In contrast, in low permafrost catchments precipitation can infiltrate through organic horizons to the deeper mineral soils, producing streams characterized by reduced DOM inputs and greater nutrient concentrations (Jones et al., 2005; MacLean et al., 1999; Petrone et al., 2006). High DOC concentration in surface soil has been attributed to DOC leaching from leaf litter and high biological activity, whereas reduction of DOC as flow passes through mineral soil is thought due to abiotic sorption of DOC to soil particles (McDowell and Likens, 1988). In CPCRW, sulfate, sodium, and calcium concentrations are considerably higher in groundwater than in shallow soil, suggesting recently thawed permafrost soil, mineral soil, and weathering as the primary sources of these solutes to streamflow (Petrone et al., 2006).
WILDFIRE AND SOLUTE DYNAMICS
Comparison of pre- and post-fire control (C1) and burned (P6) stream chemistry revealed that wildfire resulted in significant change in stream solute concentrations in the burned catchment. A significant decline in stream DOC concentration in P6 was observed two years following the wildfire. While soil organic matter was likely altered immediately post fire, the effect on stream chemistry was not realized until the wetter year of 2006 (2004 and 2005 had particularly low summer precipitation; Fig. 2). Stream DOC concentration is linked to rainfall and flushing of organic matter from soil organic horizons in CPCRW streams (Petrone et al., 2006). The flush of DOC in the C1 catchment in 2007 highlights this relationship, when DOC that had likely accumulated in soil during the dry period was delivered to the stream. The lack of a corresponding flush of DOC in P6 reinforces the interpretation that soil carbon stocks had declined in this catchment post-fire. The delayed decline in stream DOC concentration we observed may be a natural response to fire, but is more likely a product of low summer precipitation during and in the year following the fire.
The reduction in DOC concentration observed in 2006 and 2007 may be caused by several processes, including combustion of soil carbon stocks, transformation of soil carbon to more recalcitrant and less soluble carbon compounds, or decreased microbial mineralization of organic carbon due to a post-fire reduction in microbial biomass (Certini, 2005). A reduction in allochthonous inputs of coarse and fine particulate organic matter were likely not contributing to the observed reduction in DOC as riparian vegetation remained intact along the length of the P6 stream. Following the prescribed burning of the C4 catchment in 1999 (FROSTFIRE; Hinzman et al., 2003), both soil water DOC and water extractable organic carbon from the organic horizon of burned slopes declined (Shibata et al., 2003). Ping et al. (2005) proposed that in boreal forest soils, fire leads to oxidative combustion of immobile and labile soil organic matter, which is lost to the atmosphere through and conversion to CO species. The results from research following the FROSTFIRE, coupled with the importance of the organic horizon in DOC delivery to streams in the CPCRW, suggests the observed reduction in stream DOC concentration in the P6 catchment was due to combustion and conversion of soil organic carbon stocks in surface soils.
Elevated stream nitrogen concentration, as observed in this study, is a common response to fire (Bayley and Schindler, 1991; Hauer and Spencer, 1998; Minshall et al., 1997; Williams and Melack, 1997). Increased soil ammonium concentration, commonly observed following fire, is potentially due to increased microbial or thermal decomposition of vegetation and soil organic matter (Chorover et al., 1994; Wan et al, 2001). Any increase in available ammonium has the potential to be rapidly nitrified and, in turn, increase soil nitrate availability (Certini, 2005; Wan et al., 2001). Despite evidence that soil nitrification is not stimulated in black spruce stands of CPCRW following fire (Lyle, 2006; Smith et al., 2000), stream nitrate concentration increased in this study, suggesting there may have been decreased microbial competition with vegetation for ammonium and increased nitrification. Alternatively, the pattern of increased nitrate concentration may be due to increased soil active depth, with two possible mechanisms leading to nitrate concentration increases. First, increased thaw depth initiated by wildfire (latent heat transfer and lower surface albedo) may have improved soil drainage and increased aeration, thus reducing nitrate removal via denitrification. Second, increased active layer depth and improved drainage may have increased percolation of precipitation to deep mineral soils, similar to the high nitrate levels in streams draining low permafrost watersheds.
Comparing the results of our study with a previous study conducted in CPCRW on fire impacts on stream chemistry (Hinzman et al., 2003), stream nitrate in the former study increased for a short period during a single post-fire storm, but was not elevated at any other time (Petrone et al., 2007). The difference in the response of stream chemistry between the fires is likely due to extent of catchment burned and the fire intensity. The 2004 wildfire in P6 burned 65% of the catchment, compared with 28% of C4 burned during FROSTFIRE. In addition, although the riparian zone was not burned in the P6 catchment, the fire burned closer to the stream channel than during the FROSTFIRE burn. The type of vegetation and soil burned may also play an important role in nutrient export to streams. In a study of soil chemistry following wildfire in interior Alaska, soil solution nitrogen was slightly elevated in black spruce forests, in contrast to soil in aspen, birch, and white spruce stands where soil nitrogen concentrations declined following fire (Dyrness et al., 1989). The 1999 and 2004 fires both burned black spruce predominately, indicating that the difference in catchment area burned and proximity to the stream are the likely causes of the sustained elevation of stream nitrate at P6.
Stream phosphorus concentration has been shown to either dramatically increase (Earl and Blinn, 2003; Prepas et al., 2003), albeit for a brief period, or be unchanged by fire (Minshall et al., 1997; Stephens et al., 2004). Prior to the 2004 wildfire, SRP concentration in the C1 and P6 streams was below detection limit (<1 μg P L−1), and we observed no detectable post-fire increase in stream SRP concentration in P6. Other studies in the boreal forest, however, have reported a change in stream phosphorus concentration after fire. In a study in the boreal forest of Canada, total dissolved phosphorus was unchanged after fire, but particulate phosphorus increased significantly during periods of high discharge (Prepas et al., 2003). In northwest Montana, stream phosphorus concentration was over fivefold higher in burned than control streams, and this difference was sustained for a period of several years (Hauer and Spencer, 1998). The observed increase in nitrate, sulfate, and cations in P6 may have been accompanied by increased total phosphorus concentration, but given our monitoring program we were unable to detect a change.
METABOLISM, WILDFIRE, AND PERMAFROST
Stream respiration and productivity have been shown to be influenced by various environmental factors, including light (Acuña et al., 2004, Mulholland et al., 2001; Young and Huryn, 1996), temperature (Lamberti and Steinman, 1997; Sinsabaugh, 1997), nutrient concentrations, organic matter supply (Acuña et al., 2004, Kreutzweiser and Capell, 2003), discharge (Lamberti and Steinman, 1997; Uehlinger, 2000; Young and Huryn, 1996), and catchment disturbance (Houser et al., 2005). In this study, both GPP and respiration were higher in the stream draining the burned than unburned catchments, suggesting that fire did indeed have some indirect impact on ecosystem metabolism.
GPP and respiration were not correlated with temperature, although average daily stream temperature varied little over CPCRW sub-catchments in the summer of 2005. Average daily temperature ranged from 4.8 to 7.9 °C across metabolism measurement dates in all streams, and is unlikely the cause of differences in metabolism at these times. Given the high nitrate concentration in streams of CPCRW, GPP is not likely nitrogen limited. Moreover, nitrate concentrations were higher in the C2 and C4 streams than the P6 stream, but GPP rates were lower. The relationship between canopy cover and production was inconclusive, with the highest production occurring in the stream with a mixed riparian cover. Additionally, although algal standing stocks were not measured, we did not observe any obvious difference in substrate biomass among streams. Rather, production may be controlled by phosphorus availability. Given the very low concentration of SRP in CPCRW streams, limitation of production by this nutrient is highly likely, as was found in another study in interior Alaska (Slavik et al., 2004). Considering the low ambient concentration of SRP in CPCRW streams, even a minor SRP increase may have been sufficient to stimulate GPP in streams draining burned catchments.
Stream respiration rates in CPCRW were greatest in the stream impacted by fire, perhaps due to increased input of labile DOC (Findlay et al., 2003, Kreutzweizer and Capell, 2003) or extent of the hyporheic and transient storage zones (Fellows et al., 2001; Mulholland et al., 2001). Carbon bioavailability has been found to be both positively (Striegl et al., 2005) and negatively (Kawahigashi et al., 2004) correlated with permafrost extent. In CPCRW streams, significant differences in stream DOC concentration with catchment stand type suggest differences in carbon lability among streams. Moreover, White et al. (2008) found that biologically labile DOC increased with discharge in CPCRW streams, when flow bypasses mineral soils. We observed a significant decrease in C:N ratio of DOM following wildfire in P6, which may be contributing to the higher respiration rate observed in this stream. Similarly, the variability in catchment permafrost extent and soil saturation among CPCRW sub-catchments may result in differences in size of the hyporheic zone and transient storage among streams.
WARMING, WILDFIRE, PERMAFROST, AND ECOSYSTEM METABOLISM
Our ability to draw conclusions from this study regarding drivers on ecosystem metabolism in boreal forest streams is limited by the lack of treatment replication. Due to the unexpected nature of the 2004 Boundary Fire and confounding effect of permafrost coverage on stream chemistry, we were constrained by availability of pre-existing study sites. Despite this, we have established that wildfire leads to sustained change in stream solute chemistry. The role of stream nutrient concentrations and environmental parameters in driving ecosystem metabolism remains unclear and warrants further investigation, particularly during this period of rapid change in the region of discontinuous permafrost.
Understanding wildfire impacts on stream metabolism and nutrient chemistry in boreal forest streams needs to be viewed in light of discontinuous permafrost extent and vegetation feedbacks and their effect on catchment hydrology. Wildfire impacts on stream chemistry may be exacerbated by climatic warming and permafrost degradation. Over the last century, climatic warming has resulted in a 2.1% reduction in permafrost extent at CPCRW, and over one-third of the catchment now has either unstable or actively thawing permafrost (Hinzman et al., 2006). In the boreal forest, wildfire most often burns through the highly flammable black spruce stands (Kasischke et al., 2006). Permafrost typically underlies black spruce stands, and with fire dramatically altering the thermal balance of soils, permafrost is likely to become unstable and/or begin to degrade (Hinzman et al., 2006; Ping et al., 2005). After fire, black spruce may reestablish, leading to stabilization of permafrost, or, alternatively, black spruce forests may be replaced with deciduous stands (Chapin et al., 2004) under which further degradation of the underlying permafrost is likely. Fire effects on permafrost extent could be substantial under a changing climate, and where permafrost is already near the point of thawing, wildfire may accelerate the rate of degradation.
From this and previous research, it is clear that climatic change and increased wildfire occurrence may impact stream chemistry and function in the region of discontinuous permafrost via two pathways. Short- to medium-term changes to stream solute concentrations after fire may alter nutrient limitation and increase rates of whole-stream metabolism. Over longer time scales, reduction in the extent of underlying permafrost within watersheds from wildfire and/or warming will lead to change in watershed hydrologic flowpaths and nutrient inputs into streams, and again may alter ecosystem functioning.
We would like to thank Kelly Balcarczyk, Hannah Clilverd, Jessica Eichmiller, Jonathon O'Donnell, Amanda Rinehart, and Emily Schwing for their assistance in the field and laboratory. Thanks to Bob Busey, Prathap Kodial, and the Water and Environmental Resources Center at the University of Alaska Fairbanks for assistance with field data collection. We would also like to thank Richard Boone and Stuart Chapin for providing valuable comments on this research and manuscript, and to Dave Verbyla and Jamie Hollingsworth for data on Boundary Fire burn extent and mapping.
This research was part of the Bonanza Creek Long Term Ecological Research program (funded jointly by NSF grant DEB–0423442 and USDA Forest Service, Pacific Northwest Research Station grant PNW01–JV11261952–231). Additional assistance was provided by the Center for Global Change and Arctic Systems Research, University of Alaska Fairbanks.