Variability of midday net ecosystem CO2 exchange (NEE) and respiration was measured using a transect of closed system chambers spanning transitions from channel fen, permafrost plateau, and ombrotrophic flat bog land cover types during the spring melt season (26 April—6 June 2008). The primary objective was to compare fluxes from different land cover types and topographic variability within zones adjacent to and including rapid permafrost thaw. During this period, the bog was the greatest net source of CO2 to the atmosphere, followed by plateau, and fen. NEE was slightly positive (indicating CO2 loss to the atmosphere) during the snowmelt period (average = 0.009 ± 0.004 mg CO2 m-2 s-1), and increased to 0.025 ± 0.012 mg CO2 m-2 s—1, on average, possibly due to soil thaw and increased microbial activity within two days of completely snow-free conditions. Near surface soil temperature and depth to the water table were the most significant controls of soil and ground cover CO2 fluxes within chambers at all sites (p < 0.05). Analysis of historical aerial photographs and satellite imagery of the area from 1947 to 2008 indicates that plateaus are converting more rapidly into bogs than fen, where 73% of plateau areas (since 1970) that thawed had become bogs (as opposed to 27% conversion into fen). Future research requires establishment of a full ecosystem or land cover greenhouse gas and soil nutrient exchange/transfer program, including CO2 and water fluxes as well as dissolved organic and inorganic C, and CH4 losses from the soil. These results contribute to a better understanding of northern soil and ground-cover carbon exchanges as greater areas of permafrost plateaus collapse and form bogs.
Introduction
The discontinuous permafrost zone of North America is expanding as a result of a northward shift of the southern boundary of continuous permafrost, while the overall area coverage and density of permafrost has been shrinking (e.g. Quinton et al., 2010; Chasmer et al., 2011a). Ansimov and Reneva (1996) predicted that permafrost area will reduce by 15%–30% by 2050 as a result of warming air temperatures. Because northern peatlands contain approximately one-third of the global carbon (C) pool, and almost half of the atmospheric C pool (Gorham, 1991), small changes to land surface morphology, hydrology, and vegetation communities could have significant cumulative impacts to northern C exchanges and global climate change (Camill et al., 2001).
Several studies predict that northern (boreal/tundra) ecosystems will become a stronger source of C with increased air temperatures (e.g. Oechel et al., 1994; Goulden et al., 1998), while others suggest the opposite (e.g. Oechel et al., 2000; Turetsky et al., 2002; Myers-Smith et al., 2007, 2008). These tend to follow two opposing thoughts: (a) higher air temperatures will increase rates of permafrost thaw. Improved soil drainage and aridity due to creation of drainage channels may increase rates of decomposition, which will provide feedback to the global climate system by increasing global atmospheric CO2 concentrations and air temperatures (e.g. Gorham, 1991). (b) Changes in vegetation community composition could result in increased C sequestration over many years, which will exceed C released via soil decomposition (e.g. Turetsky et al., 2002). Observations have been made to support both hypotheses; however, their relevance to climate change will depend on the spatial extent of these observations within highly heterogeneous ecosystems and feedbacks between local energy, water and nutrient flows, and vegetation succession (Petrone et al., 2011).
The discontinuous permafrost region of Canada is exceedingly complex. Interactions between soil moisture, vegetation functional type and succession, and microclimate are difficult to separate and often confound relationships with C fluxes. Slight changes to local hydrology and energy receipt vary biotic and abiotic controls on the net directionality of C fluxes from sink to net source and vice versa (Moore et al., 1998; Trumbore et al., 1999; Bubier et al., 2003; Sulman et al., 2010; Wilson and Humphreys, 2010; Humphreys and Lafleur, 2011). Generally, microtopographic highs (hummocks, etc.) are dryer than microtopographic lows (hollows, depressions, etc.) (Waddington and Roulet, 1996; Petrone et al., 2005), and these often contain different vegetation community distributions and variable CO2 fluxes (e.g. Oechel et al., 1994; Goulden et al., 1998; Strack et al., 2006; Johansson et al., 2006; Petrone et al., 2011; Wilson and Humphreys, 2010). Within the fragmented discontinuous permafrost zone, the spatial variability of water storage and drainage is even more crucial because of preferential thaw and the obstruction and redirection of water from plateaus through channel fens and bogs (Quinton et al., 2010). On plateaus, excess water is shed to adjacent fens and bogs as subsurface flow over the relatively impermeable frost table (Quinton and Marsh, 1999), which descends through the active layer during thaw (Hayashi et al., 2007). This often creates wetter (not drier) peatlands as permafrost thaws, resulting in changes to plant communities (Camill, 1999). Local increases in radiation loading in areas of reduced canopy cover (as a result of tree mortality due to soil saturation) create preferential thaw of the active layer, thaw depressions, and changes to local hydrology and plant functional types (Wright et al., 2009; Hayashi et al., 2007; Quinton et al., 2010; Chasmer et al., 2011b), thereby affecting the C balance (e.g. Myers-Smith et al., 2007).
Rates of C exchange within northern peatlands, including those underlain with permafrost support (a) soil drying, increased decomposition, increased CO2 efflux; and (b) soil drying, changes to plant functional types, and increased CO2 sequestration hypotheses, depending on local environmental/topographical conditions and land cover type. Goulden et al. (1998) found that respiration increased during warm summer months in a 120-year-old black spruce forest in northern Manitoba, while photosynthesis remained fairly constant. They suggested that the depth to the frost table and increasing volume of unfrozen soil throughout the summer had the greatest influence on CO2 efflux, while saturated wetlands were not prone to seasonal increases in respiration due to anoxic soil conditions after thaw. During the winter to spring shoulder period, Bubier et al. (2002) found that the largest release of CO2 occurred when ground temperatures reached 0 °C in a temperate peatland. Degassing of soils likely occurred as a result of the buildup and release of stored CO2 as soon as the ground thawed, and also likely increased microbial activity. They note that emissions of between 3% and 50% of annual CO2 efflux in northern ecosystems vary due to soil moisture status, the length of the nongrowing season, snow cover, and temperature. Further, increased CO2 effluxes were also associated with deeper snow packs, likely as a result of soil insulation from freezing air temperatures and minimal microbial activity.
Vegetation succession and changes in plant functional type that occur as trends in local soil moisture change can also have significant influences on C fluxes. For example, Turetsky et al. (2000, 2002), found that wet internal lawns (depressions) had greater C uptake likely as a result of increased bryophyte primary production, while Petrone et al. (2011) found that Sphagnum-dominant lawns had greatest C uptake compared with feather moss-dominant lawns and Sphagnum- and feather moss-dominant depressions. Reduced overstory canopy cover and increased photosynthetically active radiation (PAR) incident on Sphagnum-dominant lawns increased photosynthesis. This is supported by greater maximum CO2 assimilation (Amax) found for Sphagnum-dominant lawns compared with those of feather moss (Petrone et al., 2011). Within peatlands in northern England, Ward et al. (2009) found that changes to plant functional groups and composition affected short-term C fluxes, especially if rapid CO2 uptake and turnover via graminoids is altered by the migration and growth of dwarf shrubs. These reduce CO2 uptake because of shading by dwarf shrubs and litter accumulation. Artificial application of N to vegetation within plots also increased biomass and leaf area, but reduced net ecosystem exchanges (NEE) due to reduced moss photosynthesis (Bubier et al., 2007). Over longer time periods, Camill et al. (2001) found that changes in biomass and decomposition were more strongly related to vegetation succession than climatic changes within the discontinuous permafrost region of northern Manitoba. However, warmer air temperatures and predicted drought (due to increased evaporation and evapotranspiration in these areas) are hypothesized to reduce the success of black spruce and Sphagnum vegetation communities, which could lead to decreases in the efficiency of ecosystem C storage over the long term (Myers-Smith et al., 2008). Further, because C accumulation is greater in areas of permafrost collapse (Myers-Smith et al., 2008), climate warming and forest fires (inducing further permafrost collapse) could result in greater C storage, while increased drying and decomposition will have the opposite effect (Camill and Clark, 1998; Robinson and Moore, 2000). Myers-Smith et al. (2007) suggested that expansion of peatland vegetation with warming air temperatures will likely result in increased C uptake, but also increased CH4 emissions as permafrost degrades.
The sensitivity of permafrost degradation to increasing air temperatures and topographic morphology (as a mechanism for water movement and transfer of heat through the soil column) has been documented in several studies (e.g. Camill and Clark, 1998; Quinton et al., 2009), but has not been related to changes in C fluxes at the plateau/peatland margin where the greatest changes are taking place. In this study, we examine for the first time relationships between the 3D topographical characteristics of the land surface derived from very high resolution airborne Light Detection And Ranging (LiDAR), soil and ground cover CO2 fluxes, and soil chemistry. The objectives of this study are to: (1) compare CO2 fluxes within and between plateau, bog, and fen sites, which will provide evidence to support soil drying/efflux and/or vegetation succession/uptake hypotheses per land cover type; (2) quantify environmental influences on fluxes, including meteorological, microtopographic, and soil property; and (3) discuss possible implications of climate change in the discontinuous permafrost zone by quantifying plateau conversion to bogs and fens from historic aerial photography and satellite imagery. This study focuses on the critical final transition from the early snowmelt to snow-free period. This is the most hydrologically dynamic period in the discontinuous permafrost zone, initiating microbial activity, photosynthesis, and gas exchanges between the biosphere and the atmosphere (Griffis et al., 2000; Humphreys and Lafleur, 2011). Feedbacks and feedback strength related to land cover type, topography, and vegetation structure will be quantified.
Study Site
CO2 exchange was measured within the Scotty Creek watershed (61°18′N, 121°18′W), 50 km south of Fort Simpson, Northwest Territories, Canada. An ∼160-m-long transect was established through channel fen, permafrost plateau, and flat ombrotrophic bog land cover types (Fig. 1). The channel fen is located on the northwest side of the permafrost plateau and contains a ground cover of primarily Carex spp., Helodium blandowii, Brachytheuum rivulane, Aulacommium palustre, and patches of shrubs with some Picea mariana (Table 1; Hayashi et al., 2004). Soils are comprised of a dense organic layer with mineral soils extending to a depth of 3 m below the water surface (Hayashi et al., 2004). The permafrost plateau rises ∼0.9 m above the surrounding peatlands and is rapidly degrading at an average rate of 0.3 m per year on the northwest side (from 1977 to 2008 determined from aerial photography; Chasmer et al., 2011a). Signs of ground surface collapse are evident along the outer perimeter of the plateau, with additional subsidence in areas of isolated surficial permafrost loss towards the center. The plateau supports an open tree canopy of Picea mariana with average canopy heights of 3.1 m (±2.2 m) and a density of 1 stem m-2 (Wright et al., 2008). Ground cover is mainly Betula glandulosa, Rubus chamaemorus, and Sphagnum capilofolium (Table 1). An organic peat layer extends to a depth of up to 8 m, under which a silt-sand layer and thick clay to siltclay layer exists (Aylesworth and Kettles, 2000). The upper plateau soil layer (0.2–0.5 m) comprises living vegetation, and a newly decomposed fibric peat layer, under which lies a layer of dense sylvic peat containing wood material, lichen, moss, rootlets, and needles at a more advanced state of decomposition (Quinton et al., 2003).
TABLE 1
Vegetation canopy characteristics (±range) within 2 m radius of center of chamber measurement location, and vegetation species types within chambers (only).
Peat deposit within bogs and fens in this region have an average thickness of 2.4 m and are underlain by a mineral silt-sand layer ranging in thickness up to 1 m. The mineral silt-sand layer overlays a thick layer (6 m on average) of glaciolacustrine clay (Rutter et al., 1973; Aylesworth et al., 1993). Both bogs and fens are seasonally frozen wetlands with the water table remaining near the ground surface throughout the year. The bog surface is covered by Sphagnum fuscum and Pleurozium schreberi. Ericaceous shrubs, such as leatherleaf (Chamaedaphne calyculata), bog birch (Betula glandulosa), and northern bog-laurel (Kalmia polifolia), grow on slightly raised hummocks. The bog has a few stunted and scattered black spruce and tamarack (Larix lancina), and is surrounded on three sides by permafrost plateaus. The fen is a large linear feature which conveys water received from surrounding permafrost plateaus to the basin outlet (Quinton et al., 2003). The fen surface is covered by a buoyant Sphagnum riparium—dominated peat mat, which is approximately 0.5 to 1.0 m thick. The peat mat responds to changes in water table height, resulting in a variable surface elevation. The mat, which sits just below the water surface, supports the development of sedges, feather mosses, grasses, and various herbs and shrubs. The establishment of trees is minimal due to ground surface instability. Beneath the peat mat (∼3 m) lies a dense organic layer with mineral soils (Hayashi et al., 2004).
Mean annual air temperature and precipitation measured at Fort Simpson is -3.2 °C and 369 mm, respectively. Approximately ∼46% of total precipitation falls as snow (from 1971 to 2000; MSC, 2002). The year studied (2008) was cooler with more snow than normal. Average annual air temperature was -4.1 °C, and ∼72% of the cumulative precipitation for that year (361 mm) fell as snow. In 2008, snowmelt began in late March and ended during the first week of May.
Materials and Methods
CO2 FLUXES
CO2 fluxes were measured each day between 10:00 and 17: 00 (CST) from 26 April to 6 June 2008 using a closed chamber system, which measured all fluxes within the chamber (c.f. Solondz et al., 2008; Petrone et al., 2011). This resulted in 44 measurements per chamber site throughout the spring snowmelt and the early plant growth (“green-up”) period (Fig. 1). Total respiration (Rtot) and net ecosystem exchange (NEE) from the soil and ground cover plants within the chamber were measured using the Vaisala CARBOCAP® GMP343 single-beam, dual-wavelength, nondispersive infrared Carbon Dioxide Probe (Vaisala Oyj, Vantaa, Finland). Terminology for NEE and Rtot refers to within-chamber fluxes and does not extend to the ecosystem (as is the case with tower-based eddy covariance methods). The probe was inserted into the side of the chamber at a height of approximately half the height of the chamber. Clear Polyvinylchloride (PVC) collars (interior radius = 12.3 cm, 40 cm tall) were inserted into the moss at each location on 12 April, two weeks prior to sampling (to minimize disturbance) and remained permanently in place for the duration of the study. As soils thawed, collars were cut deeper into the moss to a depth of 20 cm, while making note of headspace. This minimized damage to within chamber vegetation. NEE was measured within clear Lexan chambers, whereas Rtot was measured during darkened conditions within a few minutes of NEE measurements after a two-layer opaque shroud was placed over the chamber system. The enclosed surface area of the tubular chamber was 0.56 m2, with a volume of 0.03 m3. Grooves were cut along the top of each collar and were filled with water before the chamber was placed on top of the collar. This ensured an airtight seal, limiting interaction between the chamber environment and the surrounding (outside) atmosphere. The interior air temperature was cooled by a flexible tube that was used to circulate water, thereby maintaining an interior chamber atmosphere within ± 2% of ambient air temperature (Petrone et al., 2011). A small fan mounted to the inside of the chamber at a height of approximately two-thirds the height of the chamber was used to circulate air for equilibration of gas concentrations and to avoid an artificial decrease in atmospheric gradient. During each chamber measurement CO2 concentrations were measured every 30 s for a 5 min period (Table 2), and CO2 flux was calculated as:
where F is the gas flux (mg CO2 m2 s-1), Δ is the linear change in CO2 concentration with time (mg S-1), M is the molar mass of CO2 (44,010 mg mol-1), N is the molar volume of a gas (0.224 m3 mol—1) at standard temperature and pressure (STP), V is the temperature corrected volume within the chamber (m3) (using temperature measured on the GMP343 head), and A is the chamber area (m2). To ensure the quality control of exchanges per chamber, each regression slope (change in CO2 over the duration of each sampling interval, 30 s) was checked for significant deviation from a linear slope, and any measurements that were significantly different were removed, thereby removing any bias from the data set.
Timing of measurements per site were deliberately and randomly varied from day to day to minimize preferential sampling at specific hours within the day. A Kolmogorov-Smirnov (K-S) test used to determine if the timing of flux measurements throughout the day (e.g. biased grouping of measurements in the morning, late afternoon, etc.) varied significantly from one site to the next. The results indicate that the frequency and timing of data collection per site did not vary significantly (p > 0.2, in all cases), and therefore chamber data are comparable.
ENVIRONMENTAL MEASUREMENTS
Meteorology at Scotty Creek was described from measurements made at a meteorological tower located in the center of the bog, and detailed environmental measurements were made at each chamber location. These are described in Table 2.
Six representative soil cores (two from each site) were sampled near chamber locations at the end of the growing season (24 August 2008). Cores were inserted into the ground to a depth of 50 cm, and a depth of 20 cm was used (due to saturated conditions at bog/fen and compaction at plateau). Cores were cut horizontally into 5 cm subsamples, each with a volume of approximately 210 cm3. Degree of decomposition (VonPost), total N (TN), total C (TC), and carbon-to-nitrogen (C:N) ratio of each subsample was determined. For C:N ratios, the samples were oven dried at 100 °C for 24 hours and ground in a tumbling ball mill for 2–5 min until homogenized into a powder for sampling. The percentage of TC and TN was determined through combustion utilizing an Isochrom—elemental analysis, Carlo—Erba Isotope Ratio Mass Spectrometry, autocombustion carbon—nitrogen analyzer (Micromass UK, Ltd., Environmental Isotope Laboratory, Department of Earth Sciences, University of Waterloo, Waterloo, Ontario, Canada).
TABLE 2
Daily and instantaneous measurements of flux and hydrometeorology at ecosystem level and at each individual chamber.
AIRBORNE LIDAR ANALYSIS OF LOCAL TOPOGRAPHIC VARIABILITY AND SITE REPRESENTATION
Local topographic variability surrounding CO2 measurement chambers was determined from airborne LiDAR data collected on 6 August 2008 (Chasmer et al., 2011b). Laser pulse reflections from the ground surface were classified into ‘ground’ returns in TerraScan software (Terrasolid, Finland). Returns within a 5 m radius of the chamber location were used to create a digital elevation model (DEM) at a resolution of 0.5 m based on an inverse distance weighting algorithm (O'Sullivan and Unwin, 2003) in Surfer (Golden Software, Inc., Colorado). An x,y search radius of 1 m and a power of one was used so that the return integrity (and ground surface morphology from a number of points) was maintained. The representation of chamber sites within the wider study area was also determined by applying a Boolean classification to ± 1 standard deviation of topographic attributes (slope, aspect, highs, lows), and vegetation attributes (vegetation height and vegetation fractional cover) within a 2 m radius of each chamber site.
STATISTICAL ANALYSIS
Quantification of differences in measured CO2 fluxes within and between land cover types, assuming that replication exists between all chambers per land cover type, was tested using Friedman repeated measures Analysis of Variance (RM ANOVA) on ranks where fluxes are not normally distributed at P < 0.05. Where they are normally distributed, a one-way RM ANOVA (with Chi-square test) is used instead. A time series analysis of land cover differences in NEE and Rtot was also performed based on the differences between plateau vs. bog, plateau vs. fen, and bog vs. fen land cover types. The purpose was to quantify days when fluxes between land cover types were significantly different. Those days where land cover differences in average flux exceeded twice the standard deviation (approximately 95% confidence level) of the maximum land cover variance in NEE or Rtot were considered significantly different. Multiple linear regression with an ordinary least squares estimation was used to determine correlation between multiple environmental controls on CO2 fluxes. Relationships between fluxes and flux controlling mechanisms were determined using coefficient of determination (R2) adjusted for covariance between controls. In all tests, a significance level of 0.05 (α) is used. Here we assume that chambers are representative of processes occurring at plateau, fen, and bog sites, but in reality, we cannot truly represent flux variability at the same land covers throughout the basin because of limited numbers of chambers. Therefore, caution is required when interpreting results.
Results
VARIABILITY OF AVERAGE CO2 EXCHANGE BETWEEN LAND COVER TYPES
Average midday NEE at the bog, plateau, and fen chamber sites were 0.034 ± 0.025, 0.019 ± 0.013, and 0.014 ± 0.016 mg CO2 m-2 s-1, respectively, from 8 May to 6 June. A summary of fluxes and environmental conditions are presented in Table 3. Comparisons between chamber locations and land cover types based on one-way RM ANOVA indicate that there were no significant differences in NEE measured between plateau chambers (F = 2.5, p = 0.09, df = 2), but there were significant differences in NEE measured between fen chambers (RM ANOVA on ranks: Chi-square = 8.72, p = 0.013, df = 2) and bog chambers (RM ANOVA on ranks: Chi-square = 50.72, p < 0.001, df = 2). A time series analysis comparing daily average plateau, bog, and fen land cover types indicates that significant differences (p = 0.05) between plateau and bog average NEE occur only once during the entire time period examined, while significant differences between plateau and fen, and bog and fen occur on the first and last four to five days of the measurement period.
TABLE 3
Average meteorological, hydrological, and fluxes of bog, fen, and plateau chamber sites. Numbers in parentheses represent ± 1 standard deviation. Incident PAR is at chamber height.
Average midday Rtot at bog, plateau, and fen chamber sites was 0.036 ± 0.026, 0.035 ± 0.026, and 0.042 ± 0.022 mg CO2 m-2 s-1, respectively. The bog was the greatest net source of CO2 to the atmosphere (indicated by positive NEE and Rtot), ranging from 0.002 to 0.05 ± 0.02 mg CO2 m-2 s—1 despite greater Rtot from the fen, which was offset by uptake as the fen progressed from a source of CO2 to a sink throughout the period studied (Fig. 2). As soils warmed in early June, Rtot at the plateau site increased rapidly from 0.03 to 0.09 ± 0.02 mg CO2 m-2 s-1 over an 8-day period. Rtot was significantly different between chambers within land cover types (Chi-square = 24.261 [plateau]; 13.04 [bog]; and 28.58 [fen], p < 0.001, df = 2), and increases in site average PAR corresponded with increases in average NEE and Rtot at plateau and bog sites (although average Tsoil does not correspond with average incident PAR). Daily differences between average Rtot, depending on land cover type, occurred most frequently between plateau and bog, where 6 of 34 days were significantly different and between bog and fen sites where 8 of 27 days were significantly different. Greatest differences occurred at the end of the measurement period. Differences between plateau and fen land cover types were not significant.
ENVIRONMENTAL CONDITIONS AND INFLUENCE ON CO2 FLUXES
CO2 fluxes were strongly influenced by the timing of snowmelt (Fig. 3). Average snow cover at the bog and fen sites persisted until DOY 121 (30 April) and 123 (2 May), respectively, and until DOY 126 (5 May) at the plateau sites. Deeper snowpack found at the plateau sites (45% greater than fen and bog) was caused by wind-driven drift and entrapment of snow within trees, reduced sublimation, and lower melt rates. Average Tsoil at the permafrost plateau remained near 0 °C due to shallow frost table depth, but by DOY 153 (2 June) frost table depth had expanded, corresponding with an increase in Tsoil (r2 = 0.39, p < 0.001). Average Tair and Tsoil increased linearly from 0 °C to ∼14 °C (varying slightly per day). During the snow-cover period, fluxes remained below 0.02 (average = 0.015 ± 0.018) mg CO2 m-2 s-1 at all sites, and were dominated by Rtot.
Increased rates of CO2 flux into the atmosphere began on DOY 124 (3 May), increasing from 0.008 to 0.04 mg CO2 m-2 s—1 one day after the snow pack disappeared at the bog site (on average for sites, including five days prior to and following snow disappearance). Rtot at the plateau sites increased starting on DOY 128 (7 May) (from 0.01 to 0.03 mg CO2 m-2 s-1 on average), approximately two days following snow-free conditions. By DOY 128, CO2 uptake increased at both the fen and permafrost plateau sites, while they both remained small net sources (Figure 2) (0.02 and 0.03 mg CO2 m-2 S-1, respectively).
Of the environmental drivers for CO2 flux examined (including Tsoil, depth to water table, PAR, and the number of days where the depth to the frost table was less than 1 m), daily average Tsoil most strongly contributed to variability in Rtot and NEE at bog and fen sites, and depth to water table most strongly contributed to Rtot and NEE variability at plateau sites (Tables 4 and 5). PAR was a secondary driving mechanism of Rtot and NEE at all sites except at Fen chamber site 1 (FC1) (Tables 4 and 5) due to a mismatch in temporal scales between rapid variations in PAR and slower response to soil heating, ground thaw, microbial activity, and CO2 exchanges. Fen and bog sites were saturated with water such that the water table rose above the ground surface. Therefore depth to water table was not included in this comparison at these sites.
MICROTOPOGRAPHIC INFLUENCES ON FLUXES
Topographic position of chamber sites did not necessarily influence ground vegetation species cover. Two plateau sites (PC1 and PC2) were dominated by Sphagnum, while PC3 contained mostly Cladina rangiferina (reindeer lichen). Significant differences in Rtot between chambers were found at plateau sites (discussed above). PC1 and PC3 were located within depressions of up to 1.5 m deep within a 5 m radius of the site location, while PC2 was located approximately near the middle of a 2.3° westfacing slope (Fig. 4). Aboveground foliage cover was greater at PC1 and PC2 than at PC3 (Table 1), resulting in greater PAR incident at the ground surface at PC3. Soil and ground cover efflux were greatest at PC1 and PC3, corresponding with topographic depression and surprisingly lower average soil VMC at these sites. Earlier active layer thaw and exposure of belowground biomass to oxygen at PC1 and PC3 likely dominated differences in Rtot compared with PC2, where active layer remained near the soil surface for a longer period of time.
TABLE 4
Explanation of variance in Rtot determined from multiple linear regression based on least squares method at plateau (PC), bog (BC), and fen (FC) sites. Contributions to variance are listed in order of importance (1 = more important, 2 = less important).
All three bog sites were dominated by equal proportions of Sphagnum capillifolium. and S. fuscum. Similarity in vegetation ground cover and minimal range of topographic variability at these sites (maximum 40 cm) makes it difficult to link CO2 fluxes to topographic or vegetation communities. Yet significant differences in NEE and Rtot were found between sites (discussed above). BC1 had the lowest average NEE and Rtot and was dominated by an overstory of Carex spp., Andromeda polifolia, and Chamaedaphne calyculata. Increased overstory (sedge and shrub) species may have formed on this slightly sloped (0.06°, east-facing) surface (Fig. 4); however, these have not resulted in lower average Tsoil compared with other sites (Table 3). In this case, higher average site Tsoil was not related to greater CO2 flux.
Significant differences in both NEE and Rtot were found at fen chamber sites using RM ANOVA (discussed above); however, differences in topographic variability were minimal (maximum range ∼25 cm, and 5 cm within 1 m of chamber locations). Each site was located near the apex of the slope of small upraised areas within less than a meter from the chamber site location (Fig. 4). Average NEE was the same at FC1 and FC3, whilst average NEE and Rtot were less at FC2. Average Tsoil was greatest at FC2 (Table 3), which was dominated by Aulacomnium palustre (bog moss), while FC1 and FC3 were dominated by feather moss and bryophytes. Understory species coverage at FC2 was greater than the other two sites, but also had greatest average Tsoil and incident PAR.
TABLE 5
Explanation of variance in NEE (from multiple linear regression) at plateau (PC), bog (BC), and fen (FC) sites. Contributions to variance are listed in order of importance (1 = most important, 2 = less important, 3 = least important). NA = not applicable.
SOIL CHARACTERISTICS
C:N followed trophic differences between sites where highest C:N was found at the plateau and lowest at the channel fen (e.g. Trumbore et al., 1999). The ratio of C:N within the soil may also be used as an indicator of substrate quality, which may affect soil respiration through soil microbial activity. High N is often related to greater rates of decomposition of organic matter and higher rates of respiration (Raich and Schlesinger, 1992; Moore et al., 1998). Rtot was greatest at the plateau sites (average = 0.042 ± 0.02 mg CO2 m-2 s-1), despite a relatively high soil C:N of 55.9, low N (0.80%) (compared with fen and bog sites), and low VonPost humification scale (H1 to H3). Average Rtot at the fen site was slightly less than that measured at the plateau (0.041 ± 0.02 mg CO2 m-2 s-1), despite high total N (1.2%) and low C:N ratio (40.6) measured within the top 0.2 m of the soil column (Bekku et al., 2003). Further, porous and drier soils found at the plateau site (porosity average = 80.6% ± 4.1% within the upper 0.5 m) likely aided in movement and diffusion of CO2 to the atmosphere.
Discussion
To date, variability of CO2 fluxes between land cover types within northern discontinuous permafrost ecosystems is not well understood. Highly variable microtopographical controls on soil moisture, soil temperature, and plant functional types can have significant influences on soil and ground cover CO2 fluxes. At Scotty Creek, soil and ground cover fluxes were strongly related to soil moisture, soil temperature and the position of the frost table within the soil profile at daily time scales. Wright et al. (2009) found that depth to frost table on permafrost plateaus is highly variable over short distances of up to 1 m, and is strongly correlated to soil moisture distribution and lateral flow into subsurface depressions. When the soils are frozen, permafrost limits groundwater flow and storage capacity and decreases hydraulic conductivity (Wright et al., 2009). As soils thaw, the frost table and the watersaturated zone above it descends down through the soil profile, reducing saturated hydraulic conductivity and lateral flow (Quinton et al., 2000).
Topographical variability over very small distances also can have significant influences on plant functional types, which can alter the C sequestration potential of land cover types over broad areas, although these may be confounded by soil moisture (Petrone et al., 2011). Within plateaus, Wright et al. (2009) found significant correlation between depth to frost table, soil moisture, and ground cover vegetation species at the same permafrost plateau examined in this study. They found that spatial ground coverage of mosses vs. lichens affected soil temperature and the rate of soil thaw in spring. Sphagnum mosses are able to retain capillary water far above the water table due to storage capacity of hyaline cells (Vitt, 2000; Lafleur et al., 2005), while lichens are not able to retain water, and act as thermal insulation by separating wet soil from the atmosphere (Rouse and Kershaw, 1971). Sulman et al. (2010) and Bubier et al. (1998) found that soil/ground cover efflux increased in topographically upraised areas within fens, whilst fluxes remained near zero in saturated hollows. The the opposite was true for bogs. Sulman et al. (2010) hypothesized that this may be related to the relative contributions of vascular and moss vegetation; where bogs contained mostly Sphagnum ground cover, fen sites contained more shrubs and sedges. They suggested that vascular species are more productive during dry periods, because they have greater maximum productivity than mosses and typically have higher GEP. This becomes more pronounced during drying periods. Petrone et al. (2011) found that Sphagnum-dominant lawns had the greatest rate of uptake, while Sphagnum-dominant depressions, and feather moss-dominant depressions and lawns were net sources of CO2. Differences in dominant plant functional types found among plateau, fen, and bog sites in the current study may also contribute to differences in fluxes, although linkages cannot be clearly defined (Table 3). Chamber sites contain up to 100% moss ground cover, while coverage of aboveground vascular species often varies among land cover types: up to 27% and 50% of chamber planar area contain vascular species (herbs and shrubs) within plateaus and fens, respectively, while vascular species are contained within up to 12% of the chamber planar area within bogs. Coverage of vascular vs. moss vegetation at fen and bog sites do not appear to influence NEE or Rtot; however, plateau and fen chamber sites containing increased coverage of vascular vegetation (PC1 and FC1) also have higher average Rtot. This may be due to litter production and decomposition from deciduous plants within the chamber, although the complexity of soil and vegetation interactions, microtopography, and biogeochemistry (e.g. losses due to dissolved organic C, methane, variability of soil nutrients, etc., not examined here) confound these results.
Microtopography may result in increased Rtot in areas between runoff channels, via aerobic conditions (Moore et al., 1998; Aerts et al., 2001; Bubier et al., 2003; Glenn et al., 2006). In this study, PC2 was located on a sloped surface and had the highest VMC compared with PC1 and PC2, which were located in slight depressions. Although the other two plateau chamber sites contributed more CO2 to the atmosphere, on average, the proportion of Rtot to NEE was the lowest at this site. Chambers located at topographically high positions within the bog (e.g. BC1) released the lowest amount of CO2 to the atmosphere, but in terms of the proportion of Rtot to NEE, respiration at this site was 20% greater than sequestration and was a much larger proportion of the flux than that measured at the other two (topographically lower) chamber sites. All chambers within the fen were located on a very slight slope, and therefore topographical influences were difficult to determine (although fluxes were significantly different at these sites). Petrone et al. (2011) found that differences in GEP were observed between Sphagnum-dominated lawns and depressions, while differences in GEP were not observed at feather moss sites of varying topographical elevation.
During the snowmelt period, fluxes remained below 0.02 (average = 0.015 ± 0.018) mg CO2 m-2 s—1 at all sites (e.g. Zimov et al., 1996; Corradi et al., 2005; Elberling, 2007), and were dominated by Rtot due to limited photosynthesis by vegetation within chambers. A prolonged period of about 10 days of high Rtot is also found at the bog sites lagging the date of completely snow-free conditions and the start of frost table thaw by 2–3 days at plateau sites. This corresponds with an increase in air temperature from 8 °C to ∼15 °C and warming soils (Fig. 3). CO2 efflux continued at all flux sites throughout the period studied, except for the final four measurement days at the fen sites, where sites began to sequester CO2. Photosynthesis resulted in a decrease in NEE at fen sites (Fig. 2), while NEE was relatively constant at bog and plateau sites (Griffis et al., 2000; Rouse et al., 2002). Low C:N ratio and greater productivity in the form of vigorous shrub development within the fen were evidence of increased nutrient availability due to plateau runoff and the movement of water through the channel fen (NWWG, 1988). At fen chamber sites, increased vegetation production and growth correspond with greater rates of CO2 sequestration (compared with bog and plateau). This resulted in higher maximum levels of gross ecosystem production (GEP), on average at fen sites (land cover average GEP peaked at 0.05 mg CO2 m-2 s-1 [0.03 mg CO2 m-2 s-1 on average throughout the measurement period]), than that for mosses, which are dominant in bogs (average landcover GEP = 0.0008 mg CO2 m-2 s-1). Further, mosses are more sensitive to water table height than vascular vegetation found within fens (e.g. Talbot et al., 2010). Average fluxes compare well with those of Myers-Smith et al. (2007), who used a similar methodology and data set. They found that growing season average midday NEE was 0.06 (±0.017) mg CO2 m-2 s—1 along a transect containing a collapsed (previously burned) plateau, moat, and bog in Alaskan discontinuous permafrost.
AERIALLY WEIGHTING MIDDAY SPRING CO2 FLUXES BASED ON LAND COVER CHANGE
Increases in annual average air temperature of 1.6 °C and a lengthening of the snow-free period by approximately 14 days over the past 40 years (linear regression, p = 0.07, measured at Fort Simpson, Meteorological Service of Canada) has resulted in dramatic permafrost degradation and an increase in peatland areas at Scotty Creek. For example, Quinton et al. (2011) have shown that plateau area has decreased by approximately 11% between 1970 and 2008, and bog and fen areas have increased by 8% and 3%, respectively. Land cover conversion alters C fluxes, hydrology, and vegetation type/characteristics. Further, localized snow distribution and freeze-thaw cycles vary the timing of microbial activity, respiration, and photosynthesis. Historic aerial photography and high resolution topographical data provides a context from which these results may be further examined.
At Scotty Creek, the ratio of upraised areas to depressions indicates that 77% of the total watershed plateau area was above the planar surface, whilst 23% fell below (from LiDAR data). Ratios within fen and bogs were 37% to 63% and 22% to 78% (upraised:depression), respectively. Saturated land surface depressions within plateaus collapse into small, isolated bogs (Trumbore et al., 1999) and within fen and bog land cover types, depressions are more prevalent. At plateau edges, perimeter to area ratios have increased linearly from 0.03 (1947) to 0.08 (2008) (Chasmer et al., 2011a), and the proportion of isolated bogs within permafrost plateaus is ∼15% of the cumulative plateau area (in 2008). Bogs are the greatest source of CO2 to the atmosphere during the spring snowmelt and early green-up period, and 73% of the degraded plateau area has been converted into either isolated or connected bogs.
If chamber measurements are scaled according to their classified areas in 1970, a decrease in plateau area has the potential to reduce CO2 fluxes by 26% during spring snowmelt. Increased fen and bog areas may result in increased CO2 emissions by ∼9% and 23%, respectively (all else being equal to 2008, including timing of snowmelt and onset of photosynthesis). Large-scale conversion of land covers may increase CO2 effluxes by 7% based on land cover conversion alone. However, these results must be viewed with caution because only one plateau, one bog, and one fen are measured. Replication of this experiment at other sites within the basin is required to add greater certainty to these results. Further, because both the bog and fen areas were saturated with water at the time of study, the extrapolation procedure implies that all bogs and fens within the basin must be completely saturated (adding even more uncertainty). Finally, this does not include total ecosystem flux at plateau sites (due to exclusion of trees and shrubs), and within some areas containing shrubs at fen sites. The characteristics (vegetation and topographic) found within chambers represent characteristics found within other plateaus, bogs, and fens 53%, 72%, and 59% of the time within a broader area, although important plant functional types are not included in this classification.
Conclusions
During the snowmelt and early spring green-up period, midday CO2 exchanges observed at plateau, fen, and bog sites suggest that the ombrotrophic flat bog was the greatest source of CO2 to the atmosphere, followed by permafrost plateau and channel fen. Relationships between CO2 fluxes and soil/ground cover warming and drying in topographically upraised areas vs. establishment of localized plant communities as a result of shifting soil moisture regimes within plateau, fen, and bog sites are not clear during the spring thaw period at Scotty Creek. Interacting environmental driving mechanisms of C flux, including soil temperature, soil moisture, and PAR are confounded by cryospheric processes, the longevity of the frost table at the surface and movement through the soil profile, snow depth, cover, and timing of complete melt. The spatial distribution of plant functional types and biological mechanisms used by individual species to retain water and process C add additional complexity. Finally, biogeochemical processes, including nutrient flows and DOC and CH4 losses were not considered, but also likely contributed to spatial and temporal variability in C uptake and respiration. In this study, the complexity of interacting processes influencing CO2 exchanges over varying temporal and spatial scales makes quantification of a single descriptive proxy for NEE or Rtot difficult. For example, the timing of complete snow melt was strongly related to CO2 exchanges at all sites, lagging between one day (bog and fen sites) and two days (plateau sites), on average. At plateau and bog sites, springtime average soil warming as a result of increased air temperatures and PAR incident on the ground surface corresponded with increased NEE and Rtot. However, at daily time periods Tsoil and Rtot were not related to PAR due in part to a mismatch in temporal scales between rapidly varying PAR and slower rates of soil thaw and soil heating. Unraveling end member sites that have contrasting environmental, cryospheric, biogeochemical, and ecological traits will improve understanding of within land cover drivers of CO2 fluxes. This is especially important as permafrost plateaus become increasingly fragmented, forming connected bog depressions, and greater volumes of water are conveyed through fens as a result of permafrost thaw.
Acknowledgments
Support for this research was provided by the Improved Processes and Parameterisation for Prediction in Cold Regions network (Dr. John Pomeroy and Dr. Masaki Hayashi), Liidlii Kue First Nation. Funding for this project was provided by the Canadian Foundation for Climate and Atmospheric Sciences, and the Natural Sciences and Engineering Research Support. Logistical support was provided by the Aurora Research Institute and Water Survey of Canada. Field support was provided by Dr. Nicole Wright and Tyler Veness. Subsidized LiDAR data were organized, collected, and processed by Dr. Chris Hopkinson. We would like to thank Dr. Oliver Sonnentag for discussion and two anonymous reviewers for helpful comments and insights.