|Home | About | Journals | Submit | Contact Us | Français|
A snow addition experiment in moist acidic tussock tundra at Toolik Lake, Alaska, increased winter snow depths 2–3 m, and resulted in a doubling of the summer active layer depth. We used radiocarbon (14C) to (1) determine the age of C respired in the deep soils under control and deepened active layer conditions (deep snow drifts), and (2) to determine the impact of increased snow and permafrost thawing on surface CO2 efflux by partitioning respiration into autotrophic and heterotrophic components. 14C signatures of surface respiration were higher in the deep snow areas, reflecting a decrease in the proportion of autotrophic respiration. The radiocarbon age of soil pore CO2 sampled near the maximum mid-July thaw depth was approximately 1,000 years in deep snow treatment plots (45–55 cm thaw depth), while CO2 from the ambient snow areas was ~100 years old (30-cm thaw depth). Heterotrophic respiration 14C signatures from incubations were similar between the two snow depths for the organic horizon and were extremely variable in the mineral horizon, resulting in no significant differences between treatments in either month. Radiocarbon ages of heterotrophically respired C ranged from <50 to 235 years BP in July mineral soil samples and from 1,525 to 8,300 years BP in August samples, suggesting that old soil C in permafrost soils may be metabolized upon thawing. In the surface fluxes, this old C signal is obscured by the organic horizon fluxes, which are significantly higher. Our results indicate that, as permafrost in tussock tundra ecosystems of arctic Alaska thaws, carbon buried up to several thousands of years ago will become an active component of the carbon cycle, potentially accelerating the rise of CO2 in the atmosphere.
The online version of this article (doi:10.1007/s00442-009-1556-x) contains supplementary material, which is available to authorized users.
Thawing permafrost, as evidenced by thermokarst formation (Jorgenson et al. 2001, 2006) and observations of warmer permafrost temperatures across the arctic (Osterkamp and Romanovsky 1999), is one of the most important climate feedback concerns today (Romanovsky and Osterkamp 2000; IPCC 2007; Ping et al. 2008; Schuur et al. 2009) This carbon (C) pool contains 250–455 Pg of C (Miller et al. 1983; Post et al. 1985; Gorham 1991) that, upon thawing, may be released as CO2 and CH4 (McKane et al. 1997; Ping et al. 2008) or as dissolved organic C (Neff et al. 2006), creating a major positive feedback scenario. Today, with the arctic warming at least 0.5°C per decade (Serreze et al. 2000), there is a strong likelihood that permafrost will become an increasingly important source of C to the atmosphere, as active layers deepen, soils warm and microbial degradation of ancient organic matter releases CO2 and nutrients (Schuur et al. 2008). For instance, recent estimates indicate that up to 40 Pg of C may be released to the atmosphere due to thawing in Siberia over the next 40 years (Dutta et al. 2006). However, the magnitude of this release may be ameliorated by the stimulation of vegetation growth by nutrients released during thaw, particularly in the early stages (Schuur et al. 2009).
Tundra field experiments have focused primarily on quantifying the extent to which changes in summer conditions (warmer air and soil, added fertilizer, reduced light, added water) alter C biogeochemistry and feedback dynamics (Chapin and Shaver 1985; Welker et al. 1993, 1999, 2000, 2004, 2006; Wookey et al. 1993; Chapin et al. 1995, 2005; Robinson et al. 1995; Shaver et al. 1998; Callaghan et al. 2004; Mack et al. 2004; Aerts et al. 2007). However, it is now becoming increasingly clear that changes in winter snow depth may have direct and indirect effects on winter soil CO2 efflux (e.g., Welker et al. 2000), soil C and N mineralization (Schimel et al. 2004; Borner et al. 2008), and subsequent effects on annual C budgets and summer leaf mineral nutrition (Welker et al. 2005), culminating in shifts in tundra vegetation composition and canopy structure (Sturm et al. 2005; Wahren et al. 2005; Tape et al. 2006).
Deep snowpacks (~3 m) insulate soil from cold wintertime atmospheric temperatures, resulting in warmer soil temperatures and increased summer thaw, as well as shortening the growing season by up to 3 weeks (Wahren et al. 2005). Recently, we have observed that, after 13 years of experimental snow pack manipulations (Welker et al. 2005), the soil active layer under the deep snow treatment has doubled as a result of permafrost thawing (Fig. 1). In order to better understand the effects of winter snow manipulation and increased active layer depth on soil C storage and turnover, we used radiocarbon measurements of CO2 to determine the age of C respired by microbes and to partition surface CO2 efflux into autotrophic and heterotrophic sources of respiration. Radiocarbon allows us to understand C turnover on annual to millennial timescales. Conventional radiocarbon dating is used to determine ages of preserved organic materials on centennial to millennial timescales, while tracking the fate of radiocarbon from nuclear weapons testing in the mid-twentieth century can be used to ascertain turnover times on annual to decadal timescales (Broecker and Peng 1982).
This study was designed to address two questions: (1) is the Δ14C of microbial respiration in areas with permafrost thawing (deep winter snow) significantly different than areas without experimental thawing (ambient snow depth); and (2) are changes in vegetation cover or microbial respiration in response to snow depth manipulation evident in the radiocarbon signatures of the surface CO2 flux?
A 3 × 60 m snow fence was erected in 1994, in moist acidic tussock tundra at Toolik Lake, Alaska (68°38′ N, 149°34′ W), in the foothills of the Brooks Range (760 m) as part of the International Tundra Experiment (Henry and Molau 1997; Welker et al. 1997). This fence has allowed an additional 2–3 m of snow to accumulate each winter within 20 m of the leeward side of the fence (Walker et al. 1999), shortening the growing season by 2–3 weeks (Jones et al. 1998; Welker et al. 2000, 2005; Wahren et al. 2005) and increasing permafrost thaw an extra 30 cm by the end of the growing season (Fig. 1). One important consequence of this constrained growing season is that vegetation cover has been reduced by 50% and total live aboveground biomass of Eriophorum vaginatum and Betula nana has been reduced from 550 to 100 gm−2 (Welker et al., unpublished data). Future climate scenarios do not predict a snow depth increase of this magnitude, consequently the aim of this study is to examine the effects of permafrost thaw.
Our samples were derived from two areas, the snow fence area (SF) (deep winter snow depth ~3 m) area behind the snow fence and the ambient snow control (C) areas (winter snow depth ~50 cm). In early July 2007, eight 25-cm-diameter PVC collars were installed per treatment into the soil to a depth of approximately 10 cm (including litter/moss). In mid-July and mid-August, four collars with a representative mix of vegetation were capped using opaque lids fitted with a rubber closed cell foam gasket in each treatment, resulting in chamber height of 10–15 cm. CO2 was scrubbed from the chambers using soda lime, preceded by a Drierite column, and then CO2 was allowed to accumulate for approximately 30–60 min until there was at least 2 mg of C in the chamber as determined by previously measured respiration rates. This CO2 was then collected and saved using molecular sieve traps (preceded by Drierite columns to prevent the sorbtion of water), which release CO2 when heated to ~600°C. At the same time, a sample of atmospheric CO2 was also obtained. One collar in the July C treatment was excluded from subsequent analysis due to air contamination.
Gas sampling tubes were installed in mid-July near the maximum thaw depth at that time, four in the control plots, where the thaw depth was ~30 cm, and four in the treatment plots, where the maximum thaw depth was 45–55 cm. The maximum thaw depths in July 2007 were greater than the average for that time of year. One gas sample from each probe was obtained in July by slowly filling pre-evacuated canisters and brought back to the laboratory for 14C analysis (Gaudinski et al. 2000). One sample from the control and two from the treatment plots contained too much water to be analyzed.
Soils were collected in mid-July and mid-August within 2–3 days of the field respiration measurements. Four cores were collected from the surface of the organic horizon to the top of the permafrost in each treatment using a 5.5-cm-diameter corer (live surface vegetation was excluded except moss, as it was difficult to separate living and dead portions). They were frozen after collection and thawed just prior to heterotrophic respiration and C content analysis. After transport to the laboratory, soil core samples were thawed and divided into two layers representing the litter/organic and mineral horizons.
Organic and mineral soils were put into separate glass incubation jars for heterotrophic respiration analysis (four per horizon/treatment). Organic soils (including litter) were only measured from the July sampling, as the respiration rates and signatures should not change (Czimczik et al., 2006). CO2-free air was pumped into the jars, after which they were sealed, stored in darkness at 7°C, and CO2 from the soils was allowed to accumulate. For a 1-week period, CO2 concentrations were measured daily using a Licor 6252 Infrared Gas Analyzer and average respiration rates were calculated (Davidson and Trumbore 1995). Air in the jars was replaced with CO2-free air whenever CO2 concentrations exceeded 3%. At the end of a 2-week period, the CO2 in the jars was collected into evacuated canisters and stored until extraction (Dioumaeva et al. 2002; Czimczik et al. 2006).
To establish a relationship between CO2 evolution rate and temperature, the temperature of incubation for half the July jars was increased to 15°C and rates were measured as above (Dioumaeva et al. 2002; Czimczik et al. 2006). The rates were then used to calculate separate Q10s for the organic and mineral horizons using the equation:
Following incubation, soil cores were dried at 60°C for several days, and the difference in weight was used to calculate the gravimetric moisture content.
Low fluxes and leaks in chambers meant that trapped samples included air that entered the chambers following the scrubbing procedure. In order to correct for air in the field respiration samples, the typical δ13C signatures of C3 vegetation (δ13Csamp.), the atmosphere (δ13Cair), and the measured value (δ13Cmeas.) were used to calculate the air fraction (Fair) using a mass balance equation.
where δ13Cair = −8‰, δ13Csamp. = −24‰.
Subsequently, the air fraction, the atmospheric 14C signature (14Cair) and the measured 14C signature (14Cmeas.), were used to determine the 14C signature of the CO2 derived from surface respiration (14Csamp.). Samples were excluded when Fair exceeded 0.5.
where 14Cair was 50‰ (±1‰) based on three air samples collected at the time of field respiration measurement.
Finally, the proportion of C derived from autotrophic respiration was calculated using the following equation:
where FAR is the fraction of autotrophic respiration and Δ14CHR is the average flux-weighted radiocarbon signature of heterotrophic respiration (July) derived from incubations, using fluxes that were temperature corrected to 3°C for the organic and 1°C for the mineral soil layers using Q10s. Autotrophic respiration is generally made up of recent photosynthates; therefore, it was assumed to have atmospheric 14C values but C3 plant δ13C values (Ekblad and Högberg 2001). The heterotrophic respiration values measured from soil cores were extremely variable, consequently an average of the July incubations was used as the end member for each respective treatment. Field measurements likely show less spatial variability because the area measured by one sample is equivalent to 20 soil cores, consequently the average is the most realistic representation of the heterotrophic respiration component. The July heterotrophic respiration value was used for partitioning in both cases, as the flux of old C from the deep soils was likely over-estimated in the August incubations to a greater degree. Additionally, it is not clear whether CO2 from the very deep soils would leave via the atmosphere or dissolved in water under natural conditions.
C contents were analyzed on soil cores taken in August at the same time as the cores used for incubation were collected. After thawing, sub-samples were taken every 10 cm down the depth of each of the three cores per treatment and analyzed for %C in a Costech elemental analyzer.
In the laboratory, CO2 samples were purified cryogenically from air samples stored in canisters and from the molecular sieve traps by heating to 600°C to release the CO2. All CO2 was purified on a vacuum line and reduced to graphite using sealed tube zinc reduction with an iron catalyst (Xu et al. 2007). The graphite was then analyzed for 14C and 13C at the W.M. Keck AMS facility at UC-Irvine (Southon et al. 2004). Δ14C units are used to report these data, and are defined as:
In other words, the deviation (in parts per thousand) of the 14C/12C ratio of the sample from that of Oxalic Acid I, the international standard, was corrected for radioactive decay of Δ14C between the year of sampling and 1950 using the real radiocarbon half-life (Stuiver and Polach 1977). Δ14C measurements do not reflect mass dependent fractionation, because samples were normalized to a δ13C of −25‰, assuming Δ14C fractionation is twice that of δ13C (Southon et al. 2004). An aliquot of purified CO2 was retained and used for higher precision δ13C measurements made on a continuous flow stable isotope ratio mass spectrometer (Delta-Plus CFIRMS), operated in Gasbench mode.
Conventional radiocarbon ages (years BP) (Stuiver and Polach 1977) are reported on samples with Δ14C less than 0‰ and were calculated with the following equation, which uses the Libby half-life (5,568 years) and 1950 as year 0.
Prior to nuclear weapons testing and the advent of fossil fuel use, the atmospheric Δ14C signature was approximately 0‰, although slight variations existed. Δ14C signatures increased dramatically due to nuclear weapons testing in the mid-twentieth century and have subsequently been declining due to biological uptake and fossil fuel inputs (−1000‰). Samples with negative Δ14C values were fixed long enough prior to 1950 for significant radioactive decay of radiocarbon to have occurred; for these samples, we used the radiocarbon age. For C fixed since 1960, the highest Δ14C signatures are associated with decadally cycling carbon pools (i.e., most of the C was fixed closer in time to the 1960s maximum in Δ14C), while the most recently fixed C will have values close to the 2007 atmosphere (48‰).
Statistical analyses were performed using SPSS statistical software. Group means were compared using t tests. Standard deviations reported reflect spatial variability.
The Δ14C signature of CO2 in soil pores in the deepened active layer area was significantly lower at the maximum July depth of thaw by up to 70‰ when compared to the control areas (P = 0.007, t test). In the deep snow (SF) area the Δ14C values averaged −106‰ ± 39 SD (~50 cm depth) versus 41‰ ± 18 SD (30 cm depth) in the control (C) area. This corresponds to mean ages of approximately 600–1,100 years BP in the SF area compared to <100 years BP (i.e., containing ‘bomb’ radiocarbon) in the C areas.
Although permafrost thawing resulted in the decomposition of old C and the subsequent production of older CO2 in the soil profile, older C did not dominate the surface CO2 efflux (Fig. 2a). The radiocarbon values (Δ14C) of the surface CO2 efflux were on average 25‰ higher (P = 0.034, July; P = 0.038, August) in the SF area than the C area (Fig. 2a). Because the heterotrophic fluxes are dominated by decadal-aged C (higher Δ14C) and the autotrophic respiration is derived from recently fixed photosynthetic products (lower Δ14C), the higher Δ14C values in the SF area indicate that ~60% of respiration (sample range 40–80%) was derived from the vegetation in the C area and ~27% (sample range 0–55%) in the SF area (Fig. 2b).
Heterotrophic respiration rates were much higher in the organic soils than in the mineral soils and were particularly variable in the mineral soils (data not shown), but showed few other trends. There were no differences between treatments or months in the flux-weighted total heterotrophic respiration rates (org. + min., corrected for temperature) (data not shown). Q10s in mineral soil tended to be higher than the organic soil, although they were also quite variable (Q10 = 1.5 ± 0.4, organic; Q10 = 2.7 ± 1.3, mineral).
There were no differences between the permafrost thaw conditions in the Δ14C values of heterotrophic respiration from the organic layer (Fig. 3). The observed 14C values correspond to organic matter turnover times of 9–17 years in this horizon. These measurements were only made on the soils collected in July, as the results should not change between July and August (Czimczik et al. 2006).
Ages of heterotrophically respired CO2 from the mineral soil were modern (14C > 0‰) in the control (C) and modern–235 years BP in the SF soils in July (Fig. 3). Although some samples from the SF treatment respired CO2 with 14C < 0‰ (i.e., pre-1950 values), the means between treatments were not significantly different. In August, C respired from the mineral soil—which integrated a deeper thaw depth in both treatments—was older than in the July samples (P < 0.001, t test 14C), with the 14C ages ranging from 1,720 to 5,960 years BP in the control areas and from 1,525 to 8,300 years in the SF areas. There was no relationship between 14C and thaw depth in August, despite substantial differences in thaw depth. C contents were variable and showed few trends (Fig. 4), and there were no significant differences between treatments. There appeared to be an increase in C contents with depth in the deeper soils, although this was not linear.
At the surface, total (autotrophic + heterotrophic) C fluxes were dominated by C fixed 1950–present (>0‰ 14C) (Fig. 2a). Autotrophic respiration was responsible for a smaller fraction of total respiration in the SF areas (Fig. 2b), which could reflect either increased heterotrophic or decreased autotrophic respiration fluxes. Weighted heterotrophic respiration fluxes were dominated by the organic horizon and did not appear to change with the depth of permafrost thaw (data not shown). Therefore, the difference between treatments likely reflects a reduction in the autotrophic component in the increased snow treatment, which had a shortened growing season and less vegetation cover where snow depths were greatest (Walker et al. 1999; Wahren et al. 2005). Debris also accumulates in the deep drift region and is associated with reduced photosynthetically active radiation, possibly interfering with autotrophic plant activity (Fahnestock et al. 2000).
July soil pore 14CO2 values indicate C aged 600–1,100 years BP was being decomposed in the SF area soil profile, at least 8 times older than the C in soil pores from the C areas. July heterotrophic respiration measurements did not reflect the influence of this older CO2 at depth, likely because soil pore gas measurements were taken at the maximum thaw depth, which was 15–25 cm deeper than the average thaw depth represented in our soil core samples. Older C was not detectable in the soil surface respiration measurements in either treatment, due primarily to the overwhelming fluxes associated with plants and heterotrophic respiration fluxes from the organic horizon. Increased efflux from mineral soils could be detected at the surface if respiration rates from these depths were similar to the highest individual rates observed in our (isothermal 7°C) incubations; however, we did not observe this in the field. This may be because the C is not escaping via the surface flux route or because these rates are unrealistic in the field. Laboratory incubations result in changes in moisture regimes, which may increase respiration and the lack of temperature gradients in the laboratory result in the over-representation of respiration from the deep soils. It is also possible that thawing exposes small pockets of C that have very high decomposition rates, which may dominate a small sample (i.e., a soil core), but not a large one (i.e., a surface flux measurement). Despite all of this, it is clear from the soil pore gas results that some old C is respired at depth, although perhaps less than would be estimated by extrapolating incubation measurements directly to field conditions.
Our observation that the oldest ages of respired CO2 occurs at depth in mineral soils is not surprising. Measurements of heterotrophic respiration of the deeper mineral soils would be expected to be older than the shallower soils. Incubations of mineral soil samples taken in July showed that the C respired was mostly >0‰ 14C, indicating the presence of post-bomb carbon. However, in August, when mineral soils were thawed to greater depths, the samples were much older, suggesting there was indeed a relationship between thaw depth and age of C respired, but not a linear one (Fig. 3). The lack of a linear relationship may be due to cryoturbation, which results in soils that are very heterogeneous both in terms of C contents and decomposability (Fig. 4) (Bockheim and Tarnocai 1998). Other studies suggest that this mechanism is extremely common and is likely to result in the high spatial heterogeneity observed in C stocks observed in tussock tundra (Michaelson et al. 1996; Bockheim 2007).
Although this was a small study with limited statistical power, these findings indicate that exposing long frozen C to thawing has the potential to allow very old C to be decomposed, potentially accelerating the rate of CO2 increase in the atmosphere. Our findings complement those recently reported by Schuur et al. (2009). In their study of tussock tundra near Denali National Park, they find C exceeding 1,000 years BP is being lost from sites with degrading permafrost, the magnitude of this loss being dependent upon the time since thawing begun. Early on, total C sequestration increases due to increases in soil nutrient pools and the stimulation of plant growth. However, as the period of degradation increases, respiration losses increase and the landscapes become net C sources to the atmosphere. As with this study, old, deep C mobilized by thaw is a relatively minor contribution to surface fluxes.
These findings suggest that while large C stocks exist in permafrost soil (Miller et al. 1983; Post et al. 1985; Gorham 1991), upon thaw, they will initially contribute only a minor amount of CO2 to atmosphere relative to overall fluxes. However, as the area of thawed soil increases, along with nutrient availability, permafrost C will become an increasingly important CO2 source, although the future fluxes will depend on the interaction of biology with physical changes in the soil due to warming. Differences in C losses are not discernable at the surface due to the overwhelming influence of the organic horizon and plant respiration; however, it is clear that very old C at depth is labile and, when exposed to thawing, can be metabolized. How much C will be metabolized depends mainly on the thaw depth, the amount and composition of organic C stored in currently frozen soil, and physical changes in the soil environment. Although C losses from deep permafrost soils appear to be insignificant in the short term, the cumulative C loss due to decades of thawing is likely to have a significant impact on the global C cycle.
Below is the link to the electronic supplementary material.
We thank Toolik Long Term Ecological Research Station for logistical field support, Jeremy Chingell for field assistance, and Claudia Czimczik, Monica Castenada, Matt Khosh, and Xiaomei Xu for laboratory assistance. This research was supported by NSF OPP grants 0632184 and 0612534 awarded to J.M.W. and the WM Keck Carbon Cycle AMS Facility. This experiment complied with all laws of the United States.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.