Atmospheric warming is expected to cause shifts in arctic tundra vegetation composition, especially in the abundance and distribution of shrub species. Greater shrub abundance will impact the carbon exchanges between tundra ecosystems and the atmosphere, including ecosystem respiration. Here, total respiration under the shrub canopy (RT) and its components soil respiration (RS) and respiration from the ground cover vegetation (RG) were investigated at three tundra sites in the Canadian Low Arctic with varying shrub coverage. Seasonal RT and RS mean values were significantly greater (P < 0.05) at the site with greatest shrub abundance; mean values were 3.70 and 3.22 μmol m-2 s-1, respectively. Mean RG did not differ among sites; mean values ranged from 0.45 to 0.52 μmol m-2 s-1. Soil temperature exerted a stronger control on RT and RS compared to soil moisture. Differences in RT and RS among sites were attributed to differences in soil properties, such as soil total N content and bulk density. These findings suggest that belowground sources of respired carbon dioxide in Low Arctic tundra may vary with long-term shrub expansion as soil microclimate conditions and physiochemical properties adjust to changes in shrub coverage.
Warming in the Arctic is projected to alter tundra vegetation communities, including an increase in the relative abundance and cover of deciduous shrub species (Walker et al., 2003; Tape et al., 2006; Pearson et al., 2013). Such changes will ultimately result in feedbacks on carbon cycling in the Arctic. In recent decades, satellite observations of greening trends in the Arctic, which suggest increased productivity and photosynthetic uptake of carbon dioxide (CO2), have been linked to increasing temperatures (Bunn et al., 2005; Beck and Goetz, 2011). However, broad-scale changes in ecosystem respiration (ER) are more difficult to detect, and ER response to changing arctic climate and vegetation cover remains uncertain. To understand these impacts, it is necessary to separate ER into its component sources. Different sources of ER may respond differently to environmental controls (Yan et al., 2009; Subke et al., 2011) and, therefore, may be influenced differently by increasing shrub cover and density. Although tundra ER can be partitioned in various ways in terms of changing shrub cover effects, a useful approach is to divide ER into two main streams: (1) autotrophic respiration from the shrub canopy, and (2) respiration from below the canopy. The latter includes the ground cover vegetation respiration (RG) and soil respiration (RS) representing the combined respiration from roots and microbial communities. While it is expected that autotrophic respiration from the shrub canopy and root-derived RS should increase with greater plant biomass and root productivity (Hobbie and Chapin, 1998), the expected variations in RS and RG with increasing shrub cover are relatively unexplored.
There are large and potentially labile carbon (C) pools stored in arctic soils (Ping et al., 2008; Tarnocai et al., 2009; Hugelius et al., 2014) and microbe-derived RS is the major pathway for this stored C to enter the atmosphere. In most models, respiration from soil sources is determined by soil microclimate and other properties, including soil temperature, soil moisture, and aeration (Raich and Potter, 1995; Shaver et al., 2007, 2013), substrate supply (Ryan and Law, 2005), and substrate lability (Buckeridge et al., 2010). Thus, how these variables change with increased shrub cover and density will determine future trends in RS of arctic tundra.
Climate warming generally favors growth of deciduous shrubs at the expense of forbs, graminoids, bryophytes, and lichens in arctic tundra ecosystems (Walker et al., 2006) because of the interspecific competition for light and nutrients (Totland and Esaete, 2002; Fargione et al., 2003). Differences in ground cover and overstory vegetation could affect RS through changes in rhizosphere conditions, soil temperature, moisture, and substrate quantity and quality. Changes in aboveground vegetation could also influence root respiration through changes in root biomass, productivity, and turnover (Parker et al., 2015). Ground cover, such as lichens, increase surface albedo, and mosses tend to insulate the soil from summer heat, resulting in cooler soil temperatures (Gornall et al., 2007; Blok et al., 2011). In winter, taller shrubs can trap blowing snow, resulting in greater snow depth, which insulates the soil, leading to warmer winter soil temperatures (Sturm et al., 2001; Bewley et al., 2010). Conversely, by experimentally removing Betula nana canopies, Blok et al. (2010) found that soils were warmer and thaw depths were deeper in summer, which suggests that shading from shrubs in summer resulted in cooler soil temperatures and less thawing within the active layer. However, model studies of the impacts of increasing shrub abundance suggest that when climate feedbacks such as surface albedo changes are considered, large-scale expansion of shrubs could lead to warmer soil temperature and deeper permafrost thaw (Lawrence and Swenson, 2011; Bonfils et al., 2012).
While there is uncertainty about the impact of shrub expansion in regulating soil temperatures, there is also continued uncertainty regarding the temperature sensitivity of soil C decomposition (Davidson and Janssens, 2006; Karhu et al., 2014). Short-term experiments have demonstrated that the rate of soil respiration increases with temperature (Friedlingstein et al., 2006; Kirschbaum, 2006); however, long-term responses are not well known (Bradford, 2013). For instance, acclimation of plant respiration to temperature could greatly reduce the sensitivity of respiration to changes in environmental temperature resulting from climate warming (Atkin and Tjoelker, 2003). Warming-induced changes in the soil microbial community have also been observed to amplify the instantaneous increase in the rates of CO2 production in arctic soils (Hartley et al., 2008) and in other ecosystems (Hartley et al., 2007; Nie et al., 2013).
Apart from temperature, soil moisture may also change with long-term shrub expansion. The accumulation of wind-blown snow in tall shrubs can lead to greater soil moisture content in spring (Sturm et al., 2005). Conversely, some authors suggest that greater deciduous shrub cover may increase evapotranspiration rates and lead to soil drying in summer (Pearson et al., 2013). It is well established that soil respiration can vary significantly among tundra ecosystems of different hydrological regimes (Oberbauer et al., 1996; Grogan and Jonasson, 2005). Soil incubation studies have demonstrated that respiration responds positively to moisture at low levels, but negatively at saturating levels because of oxygen limitation (Johnson et al., 1996).
Besides its influences on soil temperature and moisture, expansion of shrubs will lead to more deciduous litter production, which decomposes 60% faster than evergreen species litter (Cornwell et al., 2008). Buckeridge et al. (2010) also found soil nitrogen cycling rates were ∼3 times faster in tall deciduous birch tundra compared to surrounding dwarf birch hummock vegetation. In contrast, other studies have suggested deciduous shrub litter is more recalcitrant relative to herbaceous litter resulting in slower decomposition rates (Hobbie, 1996; Baptist et al., 2010) because woody plants store more C as lignin (De Deyn et al., 2008). Thus, the net effect of shrub expansion on RS will likely depend on the vegetation that was replaced.
The objectives of this study were to quantify RS and RG across three tundra sites with varying levels of shrub cover and to examine the importance of soil temperature, moisture, and soil physiochemical properties in governing differences in respiration rates among the three sites. We hypothesized that (1) soil respiration would increase with increases in shrub cover, (2) respiration from ground cover vegetation would decrease with greater shrub cover, and (3) soil temperature and moisture would drive temporal variations in RS and RG. In testing these hypotheses, our aim was to develop a better understanding of the potential changes in respiration rates that may result through climate-induced changes in arctic tundra ecosystem vegetation.
The study area was located in the Low Arctic near the Tundra Ecosystem Research Station (TERS) at Daring Lake in Canada's Northwest Territories (64°52′N, 111°34′W). It is approximately 70 km north of treeline. The area is underlain by continuous permafrost with the maximum active layer depth typically ranging from 0.3 m to 1.0 m. Local elevation ranges from 414 to 470 m above sea level. Soils in the region typically have a thin surface organic horizon in drier areas (less than 0.1 m) and a thicker organic horizon in wetter areas (up to 0.7 m), overlying coarse mineral soil (Humphreys and Lafleur, 2011). Climate in this region is characterized by long cold winters with a short growing season (mid-June to early September). The mean annual temperature is -9 °C and mean annual precipitation ranges between 200 and 300 mm (Lafleur et al., 2012).
Three research sites, located within 3 km of each other and where the vegetation varied in percent cover and height of the dominant birch shrub, Betula glandulosa (Michx.), were selected. Sites were described and named according to the Circumpolar Arctic Vegetation Map classification (CAVM; Walker et al., 2005). The mixed tundra site (CAVM designations S1 and G4) included10- to 30-cm-tall B. glandulosa shrubs with coverage of ∼17%. The dwarf shrub site (S1) was located approximately 1 km northeast of the mixed tundra site, had 40- to 60-cm-tall B. glandulosa shrubs covering ∼45% percent of the site. The low shrub site (S2), located approximately 2.5 km southwest of the mixed tundra site, had ∼60% coverage of taller (50 to 100 cm) B. glandulosa shrubs. The other common species at all three sites included Vaccinium uliginosum (L.), Rhododendron subsect. Ledum, Empetrum nigrum (L.), Vaccinium vitis-idaea (L.), Andromeda polifolia (L.), Rubus chamaemorus (L.), some graminoids, and various mosses and lichens. All three sites were located on slightly sloping terrain (<1°) with generally coarse textured mineral soils varying from sand to silt loam with a trend toward more silt at the dwarf and low shrub sites. All sites were characterized by a large range of near-surface soil moisture conditions with no significant difference among sites (see Results section below).
At each site, six plots (1 m × 1 m) were randomly selected that contained at least one B. glandulosa shrub. In each plot, three replicate locations (∼15 cm in diameter) were established for the measurement of RT, including respiration from ground cover species such as mosses, lichen, and small stature vascular plants, and RS only. To achieve this, two collars constructed of opaque PVC tubing with an inner diameter of 11 cm and a height of 8 cm, were installed in each replicate location (Fig. 1). The collars were inserted into the soil to a depth of 5 cm, leaving approximately 3 cm above the surface. In the RT collars (labeled A, B, and C in Fig. 1, part a) surface vegetation and litter were left undisturbed (Fig. 1, part b). The second set of collars was a plant removal treatment (labeled a, b, and c in Fig. 1, part a) designed to measure RS. After installation, all green moss and aboveground vascular vegetation with associated surface roots (top ∼1–2 cm) were removed within these collars (Fig. 1, part b). Any senesced leaf litter was replaced on the soil surface. This treatment reduced the total root biomass in the collars to some extent such that RS represented respiration from roots below ∼2 cm and from soil microbial communities. The difference between RT and RS was designated as ground cover vegetation respiration (RG), which also included contributions from surface roots. We realize the potential for spatial variability and disturbance confounding results in this type of manipulation. To reduce these impacts, we set the paired collars as close together as possible and were careful to minimize the disturbance. This approach is similar to that used in past experiments (Blok et al., 2011; Nobrega and Grogan, 2008). Using this design of six plots and three replications, each respiration variable had an n of 18 for each site. Respiration measurements began 1 week after establishment of the treatments to eliminate any disturbance effects from inserting the collars (Nobrega and Grogan, 2008).
Respiration measurements were conducted weekly during the growing season (27 June–25 August 2014) at all sites. Respiration rates for each collar location were measured using a LI-COR 6400-09 soil respiration chamber linked to the 6400XT Portable Photosynthesis System (LI-COR, Lincoln, Nebraska, U.S.A.). All measurements took place between 10:00 a.m. and 4:00 p.m. local time, and selection of the site, plots, and collars was randomized for each week of measurements. Measurements were performed according to the protocol suggested by LI-COR (2003). Every measurement cycle automatically repeated three times at each collar and final soil CO2 flux for each soil collar was calculated as the average of the three measurement cycles. During the measurement of each set of paired collars, RT was measured first, then the chamber was lifted for approximately 20 s to restore ambient conditions, and then it was placed on the vegetation removal collar to measure RS.
Soil Microclimate Measurements
Soil temperature was measured at 5 cm depth (Thermo Recorder TR-51, T & D Co. Ltd., Japan) adjacent to the three paired collars and at the center of each plot (n = 4 for each plot). Soil moisture was assessed as volumetric water content (VWC) measured using a Hydrosense probe (Campbell Scientific Inc., Logan, Utah, U.S.A.) outfitted with a 12 cm probe inserted vertically into the soil beside the paired collars and at the center of each plot (n = 4). Because of the vegetation manipulation in the RS collars, near surface soil moisture and temperature may have differed from the surroundings and the RT collars. Whereas this disturbance was not assessed, we assumed the effects would be minimal.
Thaw depth beside each collar pair and at the center of each plot (n = 4) was measured by inserting a steel rod into the soil. Thaw depth was defined as the distance between the top of the moss layer and the permafrost table. These data were collected every week immediately after the respiration measurements in each plot.
Soil Properties Measurements
Surface soils were sampled for carbon and nitrogen content analysis from the edge of each plot at a weekly interval in 2014. Soil samples were taken with an open-sided soil auger (5 cm in diameter) at weekly intervals. The auger was inserted to a depth of ∼15 cm at least 50 cm from the central shrub to minimize disturbance to the shrub. Loose leaf litter, living vegetation including mosses (and up to ∼5 cm of moss stems where Sphagnum moss was present) and large roots were removed from the samples in an attempt to characterize the near-surface rooting medium of the shrubs. After collection, samples were dried at 70 °C for 48 h, and samples from each site were mixed together in sealed plastic bags and stored in a dry environment for later analysis. At the end of the field season, surface litter and soil samples were transported to Trent University, prepared, and analyzed for total carbon and nitrogen content (Elementar CNS analyzer, Hanau, Germany).
In 2014, soil bulk density for surface soil (below the litter) was determined by collecting a soil core with a metal ring (2.5 cm in radius and 7.3 cm in length = volume 143.3 cm3) at the edge of each semi-permanent plot. The soil sample was dried at 60 °C for 48 h and weighed to determine the dry mass.
Ten vegetation survey plots had already been randomly established at each site, where the point quadrat method (PQM) (Groeneveld, 1997) was used to provide estimates of species leaf area index (LAI). A vegetation survey was conducted three times during late July through early August of 2014. The PQM was performed using a frame (0.5 m × 0.5 m) with evenly strung fishing line creating 25 intersects within the grid. At each survey plot, the frame was positioned above the vegetation and a long, thin metal pin was inserted vertically through the grid at each intersect (n = 25). The number of pin contacts with a green leaf for each species was recorded. The LAI for a given species was computed as the total number of contacts for all intersections divided by the number of grid intersects. As well, at each intersect the ground cover (moss, lichen, or bare soil) was recorded and used to estimate percent cover of moss and lichen in each plot. LAI at the ten vegetation survey plots was also measured using a plant canopy analyzer (LAI-2000, LI-COR Inc.) (LAIoptical) four times from mid-July to mid-August in 2014. This method is unable to represent the leaf area associated with prostrate vegetation and other ground cover plant species and, as a result, is a measure of the leaf (and stem) area of the overstory vegetation only.
A mixed effects model was used to test for differences in respiration rates and their relationships with soil microclimate variables among the three sites. Plots and collars were treated as random effects, site as a fixed effect, and soil temperature and soil moisture measured directly beside each collar (n = 3/plot) as covariates. All soil microclimate measurements at each plot (n = 4/plot) were also tested separately for differences among sites with plots and collars treated as random effects, weeks as a repeated effect, and site as a fixed effect. Differences in soil properties and vegetation among the three sites were determined using one-way ANOVA. When a site was found to be a significant model effect, a least significant difference (LSD) test was then used to identify differences between the sites. Relationships between respiration rates and various soil properties for a given site (n = 8 for soil nitrogen, n = 48 for soil temperature and soil moisture) were also tested using simple linear regression and non-linear regression. However, given the narrow range of soil temperature variation (4 to 10 °C) and the poorer fit for the non-linear regressions, such as exponential regression of Rs on soil temperature, only linear regression analyses are presented in this paper. Prior to statistical analyses all data were checked for normality. Respiration rates at the low shrub site were non-normal, and these data were log-transformed for all sites to achieve normality for the statistical testing. We used P = 0.05 for significance level for all tests. All statistical analyses were conducted using SPSS version 16.0 software (SPSS, Chicago, Illinois, U.S.A.).
Air Temperature and Precipitation
The 10-year average (2005–2014) air temperature and rainfall at Daring Lake for July and August was 12.0 ± 1.1 °C and 89 ± 27 mm (±1 SD). In comparison, the 2014 study period (26 June–28 August) was relatively warm (13.7 °C).The distribution of precipitation during the period was extremely uneven. Two heavy rain events occurred on 16 July (57.5 mm) and 16 August (16.0 mm) (Fig. 2), which together accounted for 74% of the total precipitation (99.3 mm) during the two-month period.
Soil Temperature, Moisture, and Thaw Depth
There were significant effects of site, weeks, and their interaction on soil temperature at 5 cm depth among the three sites (Table1). Soil temperature at the low shrub site was significantly lower than at the other two sites. Soil temperature varied over the growing season, and despite a significant effect of site ´ weeks (Table 1), the seasonal patterns were somewhat similar for all sites (Fig. 3). However, there was no significant difference in soil moisture among the three sites where seasonal mean values ranged from 0.42 to 0.59 m3 m-3. Thaw depth increased over the growing season at all sites and was significantly shallower at the low shrub site compared to the mixed tundra and dwarf shrub sites (Table 1, Fig. 3).
Mean 5 cm soil temperature (T), volumetric water content (VWC) in the upper 12 cm of the soil and thaw depth in 2014. Values in parentheses are standard errors. Different letters indicate a significant difference based on one-way ANOVA followed by the LSD test.
The 0–7.3 cm soil dry bulk density ranged from 0.21 g cm-3 to 0.50 g cm-3 with significantly lower bulk density at low shrub site compared to the other sites. Bulk density at all three sites indicated the presence of substantial organic matter in the near-surface soils with values well below the root growth limiting bulk density of 1.6 g cm-3 in sandy soils. Average soil C and N were similar at the mixed tundra and dwarf shrub sites, but greater soil C and significantly greater soil total N content were found at the low shrub site (Table 2). Soil C/N also differed significantly among three sites, with greatest C/N at the mixed tundra and lowest C/N at the low shrub site.
The leaf area index (LAIoptical) of overstory shrubs was significantly greater at the low shrub site; however, ground cover vascular LAI was not significantly different among sites (Table 3). The low shrub site also had the greatest moss cover (55.5%) and the lowest lichen cover (9.5%) among the three sites.
Respiration Fluxes and Controls
RT, RS, and RG varied over the growing season (Fig. 4) with the highest values at all sites measured in the first week of the study (29 June—5 July). After the heavy rain event on 16 July and accompanying cool temperatures (Fig. 2), respiration rates decreased and did not recover entirely by 24 July when soil and air temperatures reached their highest values of the season and were coincident with peak leaf area. Following this period, respiration rates declined as soil and air temperatures decreased.
Mean surface soil dry bulk density (BD) (n = 6), soil total N and C content (n = 8) and soil C/N (n = 8) at the three sites. Parentheses indicate standard errors of the mean. Values within a column with different superscript letters were significantly different (One-way ANOVA followed by LSD test).
Mean total shrub leaf area index as measured using a plant canopy analyzer (LAIoptical) (n = 10), other ground cover vegetation LAI, and percent cover of moss and lichen as measured using the PQM (n = 10) at the three sites. Parentheses indicate standard errors of the mean. Values within a column with different superscript letters were significantly different (one-way ANOVA followed by LSD test).
RT and RS differed among the three sites with significantly greater respiration rates at the low shrub site and no differences between the other two sites (Table 4). There was no significant difference in mean RG among the three sites. All respiration rates were significantly related to soil temperature whereas only RT and RS were significantly related to soil moisture (Table 4). The significant interactions between site and these variables indicated the relationships between RT and RS and soil temperature and/or soil moisture differed among sites (Table 4). Again, RG did not relate to soil moisture and there were no significant site interactions with this environmental variable.
In general, our results show that respiration rates increased with increasing soil temperature and decreased with increasing soil moisture. Respiration rates showed stronger relationships with soil temperature than with soil moisture at all sites (Fig. 5). However, the proportion of variation in respiration explained by temperature was always less than 30% (Table 5). Using both soil moisture and soil temperature to explain variations in respiration produced modest improvements in the coefficient of determination, especially at the dwarf shrub and low shrub sites. RT and RS variations were significantly related to soil temperature at all three sites, and the response (slope) tended to increase with increasing shrub cover (Fig. 5). RT and RS were significantly negatively related to soil moisture at dwarf shrub and low shrub sites but not at the mixed tundra site. There was no significant relationship between soil temperature or moisture and RG at the low shrub site.
Mixed effects model results for respiration rates (μmol m-2 s-1) of total understory (RT), bulk soil (RS), and ground cover respiration (RG) with the covariates of soil temperature (T) and soil volumetric water content (VWC) at mixed tundra, dwarf shrub, and low shrub sites. Different superscript letters indicate a significant difference among the three sites based on LSD test.
To analyze the controls on spatial variations in respiration rates across the three sites, we regressed average respiration rates from each plot during the peak growing week, 24 July–26 July, against soil temperature, moisture, and bulk density for that same period. Spatial variations in soil temperature (Fig. 6, part a) and soil volumetric water content (Fig. 6, part b) did not relate well to spatial variations in respiration rates among the three sites. However, 61% and 66% of the variation in RT and RS, respectively, could be explained by variations in surface soil bulk density with respiration rates increasing with decreasing bulk density (Fig. 6, part c). Among the three sites, the mixed tundra site had the highest bulk density and lowest RT and RS rates, whereas the low shrub site has the lowest soil bulk density and highest RT and RS rates.
Seasonal mean RT and RS rates had similar trends as soil total N content across the three sites (Fig. 6, part d). Greater RT and RS rates at the low shrub site were associated with relatively high soil total N content, and lower RT and RS rates at the mixed tundra site were associated with relatively low soil total N content. However, with only three values it is uncertain how well these relationships may explain spatial variations in respiration at Daring Lake.
Respiration from Soil
We hypothesized that greater RS would be associated with greater shrub coverage. This is supported by the evidence that the low shrub site had significantly greater RS compared to the other sites. There was little difference in RS at the dwarf shrub and mixed tundra sites where there was also no significant difference in most vegetation, soil, and soil microclimate characteristics. Only at the low shrub site, with its taller and denser shrub cover, did soil conditions differ substantially. Although not assessed directly, the near-surface soil contained greater organic matter content (inferred from the significantly lower soil bulk density) at the low shrub site. The low shrub site also had significantly lower litter C/N compared with the other two sites (29.8 at low shrub site, 35.8 at the dwarf shrub site, and 39.7 at low shrub site). The low litter C/N generally represents relatively easily decomposable (high quality) substrate matter (Quested et al., 2003; Cornwell et al., 2008). The lowest soil C/N at the low shrub site (Table 2) may reflect highly decomposed soil organic material, high turnover rates, and/or less recalcitrant substrates (Gundersen et al., 1998; Ollinger et al., 2002; Björk et al., 2007). In contrast to our results, previ ous studies have shown that soil decomposition rates decreased with increasing shrubs due to a more recalcitrant substrate supply (Hobbie, 1996; Cornelissen et al., 2007; Baptist et al., 2010). However, this is likely to depend on the nature of the vegetation replaced, with the previous studies finding deciduous shrub litter to be more recalcitrant compared to herbaceous litter (De Deyn et al., 2008). At the Daring Lake sites, ericaceous and evergreen dwarf shrubs such as E. nigrum and R. decumbens are common and decrease in abundance in areas where B. glandulosa dominates. As a result, a shift to deciduous over evergreen litter may be driving faster decomposition and greater soil respiration rates. Parker et al. (2015) also found greater soil respiration rates for higher productivity deciduous shrub vegetation than for lower productivity tundra heath, and this might be caused by the increased growth rates of fungal hyphae across the transition from heath to shrub. Although we cannot conclude that soil total N or C/N are the dominating controls on soil respiration rates with our limited sample size, these observations at least highlight the importance of substrate characteristics in controlling soil respired CO2 with increased shrub cover and density.
Coefficient of determination (R2) for soil temperature (T), soil volumetric water content (VWC), and T × VWC in least squares linear regression analysis to predict variation in total understory (RT), soil (RS), and ground cover respiration (RG). Levels of significance were represented as “*” P < 0.05, “**” P < 0.001, and “***” P < 0.0001.
Although not assessed directly, contributions from root respiration might also differ among sites. The variation in root respiration may vary in response to aboveground plant phenology and the differences in community composition among sites (Rey et al., 2002; Yuste et al., 2004). The low shrub site likely had greater root productivity because of the greater shrub biomass, thereby resulting in greater root respiration. Significantly higher shrub photosynthetic capacity at the low shrub site (data not shown) also suggests stronger root respiration at this site because root respiration has been observed to be tightly coupled with allocation of photosynthate to the rhizosphere and aboveground productivity (Kuzyakov and Cheng, 2001; Vargas et al., 2011).
Respiration from Ground Cover Vegetation
Increased shrub cover in arctic tundra ecosystems might reduce the abundance of lichens and mosses and other ground cover plants (forbs and graminoids) because of the dominating role of taller deciduous shrubs in obtaining light and nutrients (Walker et al., 2006; Dahlgren et al., 2006). In this study, we found relatively low respiration rates from ground cover vegetation at the three sites (Fig. 4). The contribution of RG to RT varied from 23% at the mixed tundra site, and 21% at the dwarf shrub site, to 15% at the low shrub site. This might be caused by the low rates of respiration from mosses and lichens, as observed in a previous study (Sveinbjörnsson and Sonesson, 1997). We also found no significant differences in RG among sites (Table 4) despite an inverse trend of moss cover and lichen cover across the sites, as the site with more moss cover had less lichen cover (Table 3). Thus, an increase in the contribution of respired CO2 from more mosses might be offset by a decrease in respired CO2 from fewer lichens, resulting in little net effect on total RG. In addition, the lack of a major difference in the understory vascular plant LAI between the low shrub site and the other two sites may have also been a reason for why RG was similar among sites.
Spatial and Temporal Variation of Respiration
Instantaneous RS rates measured at our sites ranged from 1.67 to 4.45 μmol m-2 s-1 over the growing seasons and were comparable with soil respiration values reported for other arctic ecosystems (Elberling, 2007; Nobrega and Grogan, 2008; Cahoon et al., 2016). Soil temperature is the dominant factor regulating temporal variations in CO2 effluxes in the Arctic (Kim et al., 2013; Jensen et al., 2014) and other terrestrial ecosystems (Davidson et al., 1998; Davidson and Janssens, 2006). In general, the role of soil moisture in controlling soil respiration has not been as clearly distinguished as soil temperature. Incubation studies in tundra ecosystems (Johnson et al., 1996) and temperate ecosystems (Davidson et al., 1998) indicate that soil respiration positively correlates with soil moisture at low levels (<0.2 m3 m-3), but negatively correlates with soil moisture at high levels (>0.2 m3 m-3) due to oxygen limitation. We found weak negative relationships between RS and soil moisture at the dwarf shrub and low shrub sites where soil moisture content exceeded 0.2 m3 m-3 during the growing season (Fig. 5, Table 5). Jensen et al. (2014) reported variations in annual CO2 effluxes in a tundra ecosystem as a result of changes in soil moisture because of severe rainfall events. A substantial input of rainwater can inhibit soil production of CO2 and limit diffusion rates to the atmosphere. We observed a decline in RS and RT at all three sites immediately following a heavy precipitation event (and cool air and soil temperatures) on 16 July 2014. Respiration rates did not fully recover despite a return to warm conditions in the following week. This may have been because the soils remained relatively wet; limiting both CO2 production and transport rates (Oechel et al., 2000; Oberbauer et al., 2007). These findings highlight the importance of unusually large rain events in controlling respired CO2 rates, which could be important considering expected increases in precipitation in the Arctic associated with climate change (Macdonald et al., 2005; Koenigk et al., 2013).
Although our hypothesis that soil temperature and moisture would be the main drivers of respiration variability was supported by the temporal variations in these variables, spatial variations in RT and RS were better correlated with soil physiochemical properties as noted above. The low shrub site with the coldest soil had the highest RS and RT, therefore soil temperature was clearly not the primary control on respiration differences among sites. Nobrega and Grogan (2008) found soil moisture to be the dominating factor controlling patterns of RS among tundra ecosystems across a moisture gradient, where summer soil respiration was greatest at moderate soil moisture conditions (∼0.3 m3 m-3). There was no relationship between peak season respiration and soil moisture across the three sites (Fig. 6b). Instead, we also found RS and RT increased with increasing soil moisture content up to ∼0.45 m3 m-3 and decreased with wetter conditions. It is notable, however, that there was large variability in respiration rates across the range of soil moisture values.
The three sites examined in this study experience the same climate conditions, yet have developed differences in shrub cover due to other differences in habitat characteristics. Although we recognize this fact tempers possible interpretations from our study, if we consider the three sites to reflect the impact of long-term shrub expansion in similar tundra of the Low Arctic, there are important implications arising from this study. No significant differences in mean RT, RS, and RG rates between the mixed tundra and dwarf shrub sites suggests a weak response of respiration rates in the early stages of shrub expansion (i.e., as invading shrubs start to become dominant on herbaceous tundra). As shrub expansion approaches a relative mature phase (i.e., similar to our low shrub site) continued inputs of larger amounts and relatively higher quality litter may substantially improve the decomposability of soil organic matter, to support increased Rs and RT.. Our results suggest that total understory respiration in Low Arctic tundra may not be fixed with long-term shrub expansion. Presently, these notions are untested, but ecosystem models could be used to disentangle the linkages between shrub expansion, soil nutrients, and ecosystem respiration.
This study documented the temporal and spatial variations in soil and ground cover vegetation respiration at three Low Arctic sites with varying birch shrub cover. Despite colder soils, increased soil-respired CO2 at the taller and denser deciduous shrub site was likely caused by greater plant productivity and higher soil nitrogen content. This highlights the importance of shifts in litter input and soil nitrogen content in driving the future soil respiration trend during long-term shrub expansion. Future shrub expansion in arctic tundra will alter the exchanges of energy and carbon that can impact the trajectory of future atmospheric warming. Hence, accurately predicting these changes is important, and our study shows that the amount and nature of shrub cover change combined with direct and indirect effects on soil climate, vegetation structure, and biogeochemistry are all important factors that should be considered when predicting the response of tundra ecosystem respiration to continued climate warming.
This work was supported by the NSERC Discovery grants to P. M. Lafleur and E. R. Humphreys. We thank Sean Arruda, Kristyn Foster, Mary-Claire Buell, and Claire Elliott for their help in the field and laboratory. We also thank Karin Clark and the Daring Lake TERS staff for assistance in the field. Finally, we are thankful for the comments of two anonymous reviewers, which greatly improved this manuscript.
- Atkin, O. K., and Tjoelker, M. G., 2003: Thermal acclimation and the dynamic response of plant respiration to temperature. Trends in Plant Science, 8(7): 343–351. Google Scholar
- Baptist, F., Yoccoz, N. G., and Choler, P., 2010: Direct and indirect control by snow cover over decomposition in alpine tundra along a snowmelt gradient. Plant and Soil, 328(1-2): 397–410. Google Scholar
- Beck, P. S. A., and Goetz, S. J., 2011: Satellite observations of high northern latitude vegetation productivity changes between 1982 and 2008: ecological variability and regional differences. Environmental Research Letters , 6: 045501. Google Scholar
- Bewley, D., Essery, R., Pomeroy, J., and Menard, C., 2010: Measurements and modelling of snowmelt and turbulent heat fluxes over shrub tundra. Hydrology and Earth System Sciences, 14(7): 1331–1340. Google Scholar
- Björk, R. G., Klemedtsson, L., Molau, U., Harndorf, J., Odman, A., and Giesler, R., 2007: Linkages between N turnover and plant community structure in a tundra landscape. Plant and Soil, 294(1-2): 247–261. Google Scholar
- Blok, D., Heijmans, M., Schaepman-Strub, G., Konono, A. V., Maximov, T. C., and Berendse, F., 2010: Shrub expansion may reduce summer permafrost thaw in Siberian tundra. Global Change Biology , 16(4):1296–1305. Google Scholar
- Blok, D., Heijmans, M., Schaepman-Strub, G., van Ruijven, J., Parmentier, F. J. W., Maximov, T. C., and Berendse, F., 2011: The cooling capacity of mosses: controls on water and energy fluxes in a Siberian tundra site. Ecosystems , 14: 1055–1065. Google Scholar
- Bonfils, C. J. W., Phillips, T. J., Lawrence, D. M., Cameron-Smith, P., Riley, W. J., and Subin, Z. M., 2012: On the influence of shrub height and expansion on northern high latitude climate. Environmental Research Letters, 7: 015503. Google Scholar
- Bradford, M. A., 2013: Thermal adaptation of decomposer communities in warming soils. Frontiers in Microbiology, 4: 00333. Google Scholar
- Buckeridge, K. M., Zufelt, E., Chu, H., and Grogan, P., 2010: Soil nitrogen cycling rates in low arctic shrub tundra are enhanced by litter feedbacks. Plant and Soil, 330(1-2): 407–421. Google Scholar
- Bunn, A. G., Goetz, S. J., and Fiske, G. J., 2005: Observed and predicted responses of plant growth to climate across Canada. Geophysical Research Letters , 32: L16710. Google Scholar
- Cahoon, S. M. P., Sullivan, P. F., Gamm, C., Welker, J. M., Eissenstat, D., and Post, E., 2016: Limited variation in proportional contributions of auto- and heterotrophic soil respiration, despite large differences in vegetation structure and function in the Low Arctic. Biogeochemistry, 127(2-3): 339–351. Google Scholar
- Cornelissen, J. H. C., van Bodegom, P. M., Aerts, R., Callaghan, T. V., van Logtestijn, R. S. P., Alatalo, J., Chapin, F. S., Gerdol, R., Gudmundsson, J., Gwynn-Jones, D., Hartley, A. E., Hik, D. S., Hofgaard, A., Jonsdottir, I. S., Karlsson, S., Klein, J. A., Laundre, J., Magnusson, B., Michelsen, A., Molau, U., Onipchenko, V. G., Quested, H. M., Sandvik, S. M., Schmidt, I. K., Shaver, G. R., Solheim, B., Soudzilovskaia, N. A., Stenstrom, A., Tolvanen, A., Totland, O., Wada, N., Welker, J. M., Zhao, X. Q., and Team, M. O. L., 2007: Global negative vegetation feedback to climate warming responses of leaf litter decomposition rates in cold biomes. Ecology Letters, 10(7): 619–627. Google Scholar
- Cornwell, W. K., Cornelissen, J. H. C., Amatangelo, K., Dorrepaal, E., Eviner, V. T., Godoy, O., Hobbie, S. E., Hoorens, B., Kurokawa, H., Perez-Harguindeguy, N., Quested, H. M., Santiago, L. S., Wardle, D. A., Wright, I. J., Aerts, R., Allison, S. D., van Bodegom, P., Brovkin, V., Chatain, A., Callaghan, T. V., Diaz, S., Garnier, E., Gurvich, D. E., Kazakou, E., Klein, J. A., Read, J., Reich, P. B., Soudzilovskaia, N. A., Victoria Vaieretti, M., and Westoby, M., 2008: Plant species traits are the predominant control on litter decomposition rates within biomes worldwide. Ecology Letters, 11(10): 1065–1071. Google Scholar
- Dahlgren, J. P., Eriksson, O., Bolmgren, K., Strindell, M., and Ehrlen, J., 2006: Specific leaf area as a superior predictor of changes in field layer abundance during forest succession. Journal of Vegetation Science, 17(5): 577–582. Google Scholar
- Davidson, E. A., and Janssens, I. A., 2006: Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature, 440(7081): 165–173. Google Scholar
- Davidson, E. A., Belk, E., and Boone, R. D., 1998: Soil water content and temperature as independent or confounded factors controlling soil respiration in a temperate mixed hardwood forest. Global Change Biology, 4(2): 217–227. Google Scholar
- De Deyn, G. B., Cornelissen, J. H. C., and Bardgett, R. D., 2008: Plant functional traits and soil carbon sequestration in contrasting biomes. Ecology Letters, 11(5): 516–531. Google Scholar
- Elberling, B., 2007: Annual soil CO2 effluxes in the High Arctic: the role of snow thickness and vegetation type. Soil Biology & Biochemistry, 39(2): 646–654. Google Scholar
- Fargione, J., Brown, C. S., and Tilman, D., 2003: Community assembly and invasion: an experimental test of neutral versus niche processes. Proceedings of the National Academy of Sciences of the United States of America, 100(15): 8916–8920. Google Scholar
- Friedlingstein, P., Cox, P., Betts, R., Bopp, L., Von Bloh, W., Brovkin, V., Cadule, P., Doney, S., Eby, M., Fung, I., Bala, G., John, J., Jones, C., Joos, F., Kato, T., Kawamiya, M., Knorr, W., Lindsay, K., Matthews, H. D., Raddatz, T., Rayner, P., Reick, C., Roeckner, E., Schnitzler, K. G., Schnur, R., Strassmann, K., Weaver, A. J., Yoshikawa, C., and Zeng, N., 2006: Climate-carbon cycle feedback analysis: results from the (CMIP)-M-4 model intercomparison. Journal of Climate, 19(14): 3337–3353. Google Scholar
- Gornall, J. L., Jonsdottir, I. S., Woodin, S. J., and van der Wal, R., 2007: Arctic mosses govern below-ground environment and ecosystem processes. Oecologia, 153(4): 931–941. Google Scholar
- Groeneveld, D. P., 1997: Vertical point quadrat sampling and an extinction factor to calculate leaf area index. Journal of Arid Environments , 36(3): 475–485. Google Scholar
- Grogan, P., and Jonasson, S., 2005: Temperature and substrate controls on intra-annual variation in ecosystem respiration in two subarctic vegetation types. Global Change Biology, 11(3): 465–475. Google Scholar
- Gundersen, P., Emmett, B. A., Kjonaas, O. J., Koopmans, C. J., and Tietema, A., 1998: Impact of nitrogen deposition on nitrogen cycling in forests: a synthesis of NITREX data. Forest Ecology and Management, 101(1-3): 37–55. Google Scholar
- Hartley, I. P., Heinemeyer, A., Evans, S. P., and Ineson, P., 2007: The effect of soil warming on bulk soil vs. rhizosphere respiration. Global Change Biology, 13(12): 2654–2667. Google Scholar
- Hartley, I. P., Hopkins, D. W., Garnett, M. H., Sommerkorn, M., and Wookey, P. A., 2008: Soil microbial respiration in arctic soil does not acclimate to temperature. Ecology Letters, 11(10): 1092–1100. Google Scholar
- Hobbie, S. E., 1996: Temperature and plant species control over litter decomposition in Alaskan tundra. Ecological Monographs, 66(4): 503–522. Google Scholar
- Hobbie, S. E. and Chapin, F. S., 1998: Response of tundra plant biomass, aboveground production, nitrogen, and CO2 flux to experimental warming. Ecology, 79(5): 1526–1544. Google Scholar
- Hugelius, G., Strauss, J., Zubrzycki, S., Harden, J. W., Schuur, E. A. G., Ping, C. L., Schirrmeister, L., Grosse, G., Michaelson, G. J., Koven, C. D., O'Donnell, J. A., Elberling, B., Mishra, U., Camill, P., Yu, Z., Palmtag, J., and Kuhry, P., 2014: Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps. Biogeosciences, 11(23): 6573–6593. Google Scholar
- Humphreys, E. R., and Lafleur, P. M., 2011: Does earlier snowmelt lead to greater CO2 sequestration in two low Arctic tundra ecosystems? Geophysical Research Letters , 38: L09703. Google Scholar
- Jensen, A. E., Lohse, K. A., Crosby, B. T., and Mora, C. I., 2014: Variations in soil carbon dioxide efflux across a thaw slump chronosequence in northwestern Alaska. Environmental Research Letters, 9(2): 025001. Google Scholar
- Johnson, L. C., Shaver, G. R., Giblin, A. E., Nadelhoffer, K. J., Rastetter, E. R., Laundre, J. A., and Murray, G. L., 1996: Effects of drainage and temperature on carbon balance of tussock tundra microcosms. Oecologia, 108(4): 737–748. Google Scholar
- Karhu, K., Auffret, M. D., Dungait, J. A. J., Hopkins, D. W., Prosser, J. I., Singh, B. K., Subke, J. A., Wookey, P. A., Agren, G. I., Sebastia, M. T., Gouriveau, F., Bergkvist, G., Meir, P., Nottingham, A. T., Salinas, N., and Hartley, I. P., 2014: Temperature sensitivity of soil respiration rates enhanced by microbial community response. Nature, 513(7516): 81–84. Google Scholar
- Kim, Y., Kim, S. D., Enomoto, H., Kushida, K., Kondoh, M., and Uchida, M., 2013: Latitudinal distribution of soil CO2 efflux and temperature along the Dalton Highway, Alaska. Polar Science, 7(2): 162–173. Google Scholar
- Kirschbaum, M. U. F., 2006: The temperature dependence of organic-matter decomposition—Still a topic of debate. Soil Biology & Biochemistry, 38(9): 2510–2518. Google Scholar
- Koenigk, T., Brodeau, L., Graversen, R. G., Karlsson, J., Svensson, G., Tjernstrom, M., Willen, U., and Wyser, K., 2013: Arctic climate change in 21st century CMIP5 simulations with EC-Earth. Climate Dynamics, 40(11-12): 2719–2743. Google Scholar
- Kuzyakov, Y., and Cheng, W., 2001: Photosynthesis controls of rhizosphere respiration and organic matter decomposition. Soil Biology & Biochemistry, 33(14): 1915–1925. Google Scholar
- Lafleur, P. M., Humphreys, E. R., St Louis, V. L., Myklebust, M. C., Papakyriakou, T., Poissant, L., Barker, J. D., Pilote, M., and Swystun, K. A., 2012: Variation in peak growing season net ecosystem production across the Canadian Arctic. Environmental Science & Technology, 46(15): 7971–7977. Google Scholar
- Lawrence, D. M., and Swenson, S. C., 2011: Permafrost response to increasing arctic shrub abundance depends on the relative influence of shrubs on local soil cooling versus large-scale climate warming. Environmental Research Letters, 6: 024014. Google Scholar
- LI-COR, 2003: 6400-09 Soil CO2 Flux Chamber Instruction Manual. LI-COR Inc., Lincoln, Nebraska, U.S.A. Google Scholar
- Macdonald, R. W., Harner, T., and Fyfe, J., 2005: Recent climate change in the Arctic and its impact on contaminant pathways and interpretation of temporal trend data. Science of the Total Environment, 342(1-3): 5–86. Google Scholar
- Nie, M., Pendall, E., Bell, C., Gasch, C. K., Raut, S., Tamang, S., and Wallenstein, M. D., 2013: Positive climate feedbacks of soil microbial communities in a semi-arid grassland. Ecology Letters, 16(2): 234–241. Google Scholar
- Nobrega, S., and Grogan, P., 2008: Landscape and ecosystem-level controls on net carbon dioxide exchange along a natural moisture gradient in Canadian Low Arctic tundra. Ecosystems, 11(3): 377–396. Google Scholar
- Oberbauer, S. F., Gillespie, C. T., Cheng, W. X., Sala, A., Gebauer, R., and Tenhunen, J. D., 1996: Diurnal and seasonal patterns of ecosystem CO2 efflux from upland tundra in the foothills of the Brooks Range, Alaska, USA. Arctic and Alpine Research, 28(3): 328–338. Google Scholar
- Oberbauer, S. F., Tweedie, C. E., Welker, J. M., Fahnestock, J. T., Henry, G. H. R., Webber, P. J., Hollister, R. D., Walker, M. D., Kuchy, A., Elmore, E., and Starr, G., 2007: Tundra CO2 fluxes in response to experimental warming across latitudinal and moisture gradients. Ecological Monographs, 77(2): 221–238. Google Scholar
- Oechel, W. C., Vourlitis, G. L., Verfaillie, J., Crawford, T., Brooks, S., Dumas, E., Hope, A., Stow, D., Boynton, B., Nosov, V., and Zulueta, R., 2000: A scaling approach for quantifying the net CO2 flux of the Kuparuk River Basin, Alaska. Global Change Biology, 6: 160–173. Google Scholar
- Ollinger, S. V., Smith, M. L., Martin, M. E., Hallett, R. A., Goodale, C. L., and Aber, J. D., 2002: Regional variation in foliar chemistry and N cycling among forests of diverse history and composition. Ecology, 83(2): 339–355. Google Scholar
- Parker, T. C., Subke, J., and Wookey, P. A., 2015: Rapid carbon turnover beneath shrub and tree vegetation is associated with low soil carbon stocks at a subarctic treeline. Global Change Biology, 21(5): 2070–2081. Google Scholar
- Pearson, R. G., Phillips, S. J., Loranty, M. M., Beck, P. S. A., Damoulas, T., Knight, S. J., and Goetz, S. J., 2013: Shifts in arctic vegetation and associated feedbacks under climate change. Nature Climate Change, 3(7): 673–677. Google Scholar
- Ping, C.L., Michaelson, G. J., Jorgenson, M. T., Kimble, J. M., Epstein, H., Romanovsky, V. E., and Walker, D. A., 2008: High stocks of soil organic carbon in the North American arctic region. Nature Geoscience, 1(9): 615–619. Google Scholar
- Quested, H. M., Press, M. C., and Callaghan, T. V., 2003: Litter of the hemiparasite Bartsia alpina enhances plant growth: evidence for a functional role in nutrient cycling. Oecologia, 135(4): 606–614. Google Scholar
- Raich, J. W., and Potter, C. S., 1995: Global patterns of carbon dioxide emissions from soils. Global Biogeochemical Cycles, 9(1): 23–36. Google Scholar
- Rey, A., Pegoraro, E., Tedeschi, V., De Parri, I., Jarvis, P. G., and Valentini, R., 2002: Annual variation in soil respiration and its components in a coppice oak forest in central Italy. Global Change Biology, 8(9): 851–866. Google Scholar
- Ryan, M. G., and Law, B. E., 2005: Interpreting, measuring, and modeling soil respiration. Biogeochemistry, 73(1): 3–27. Google Scholar
- Shaver, G. R., Street, L. E., Rastetter, E. B., Van Wijk, M. T., and Williams, M., 2007: Functional convergence in regulation of net CO2 flux in heterogeneous tundra landscapes in Alaska and Sweden. Journal of Ecology, 95(4): 802–817. Google Scholar
- Shaver, G. R., Rastetter, E. B., Salmon, V., Street, L. E., van de Weg, M. J., Rocha, A., van Wijk, M. T., and Williams, M., 2013: Pan-Arctic modelling of net ecosystem exchange of CO2. Philosophical Transactions of the Royal Society B—Biological Sciences, 368: 20120485. Google Scholar
- Sturm, M., McFadden, J. P., Liston, G. E., Chapin, F. S., Racine, C. H., and Holmgren, J., 2001: Snow-shrub interactions in arctic tundra: a hypothesis with climatic implications. Journal of Climate, 14(3): 336–344. Google Scholar
- Sturm, M., Schimel, J., Michaelson, G., Welker, J. M., Oberbauer, S. F., Liston, G. E., Fahnestock, J., and Romanovsky, V. E., 2005: Winter biological processes could help convert arctic tundra to shrubland. Bioscience, 55(1): 17–26. Google Scholar
- Subke, J. A., Voke, N. R., Leronni, V., Garnett, M. H., and Ineson, P., 2011: Dynamics and pathways of autotrophic and heterotrophic soil CO2 efflux revealed by forest girdling. Journal of Ecology, 99(1): 186–193. Google Scholar
- Sveinbjörnsson, B., and Sonesson, M., 1997: Photosynthesis and respiration in mosses and lichens. In Oechel, W. C., Callaghan, T., Gilmanov, T., Holten, J. I., Maxwell, B., Molau, U., and Sveinbjörnsson, B. (eds.), Global Change and Arctic Terrestrial Ecosystems. New York: Springer-Verlag, 113–128. Google Scholar
- Tape, K., Sturm, M., and Racine, C., 2006: The evidence for shrub expansion in northern Alaska and the Pan-Arctic. Global Change Biology, 12(4): 686–702. Google Scholar
- Tarnocai, C., Canadell, J. G., Schuur, E. A. G., Kuhry, P., Mazhitova, G., and Zimov, S., 2009: Soil organic carbon pools in the northern circumpolar permafrost region. Global Biogeochemical Cycles, 23: GB2023. Google Scholar
- Totland, O., and Esaete, J., 2002: Effects of willow canopies on plant species performance in a low-alpine community. Plant Ecology, 161(2): 157–166. Google Scholar
- Vargas, R., Carbone, M. S., Reichstein, M., and Baldocchi, D. D, 2011: Frontiers and challenges in soil respiration research: from measurements to model-data integration. Biogeochemistry , 102(1): 1–13. Google Scholar
- Walker, D. A., Epstein, H. E., Jia, G. J., Balser, A., Copass, C., Edwards, E. J., Gould, W. A., Hollingsworth, J., Knudson, J., Maier, H. A., Moody, A., and Raynolds, M. K., 2003: Phytomass, LAI, and NDVI in northern Alaska: relationships to summer warmth, soil pH, plant functional types, and extrapolation to the circumpolar Arctic. Journal of Geophysical Research-Atmospheres, 108(D2): 000986. Google Scholar
- Walker, D. A., Raynolds, M. K., Daniels, F. J. A., Einarsson, E., Elvebakk, A., Gould, W. A., Katenin, A. E., Kholod, S. S., Markon, C. J., Melnikov, E. S., Moskalenko, N. G., Talbot, S. S., Yurtsev, B. A., and Team, C., 2005: The circumpolar arctic vegetation map. Journal of Vegetation Science, 16(3): 267–282. Google Scholar
- Walker, M. D., Wahren, C. H., Hollister, R. D., Henry, G. H. R., Ahlquist, L. E., Alatalo, J. M., Bret-Harte, M. S., Calef, M. P., Callaghan, T. V., Carroll, A. B., Epstein, H. E., Jonsdottir, I. S., Klein, J. A., Magnusson, B., Molau, U., Oberbauer, S. F., Rewa, S. P., Robinson, C. H., Shaver, G. R., Suding, K. N., Thompson, C. C., Tolvanen, A., Totland, O., Turner, P. L., Tweedie, C. E., Webber, P. J., and Wookey, P. A., 2006: Plant community responses to experimental warming across the tundra biome. Proceedings of the National Academy of Sciences of the United States of America, 103(5): 1342–1346. Google Scholar
- Yan, J. X., Li, H. J., Tang, Y., and Zhang, Y. H., 2009: Effect of vegetation types on soil respiration characteristics on a smaller scale. Huan Jing KeXue, 30(11): 3121–3129. Google Scholar
- Yuste, J. C., Janssens, I. A., Carrara, A., and Ceulemans, R., 2004: Annual Q10 of soil respiration reflects plant phenological patterns as well as temperature sensitivity. Global Change Biology, 10(2): 161–169. Google Scholar