Abstract

Little is known about the short-term impacts of warming on native plant community dynamics in the northern Canadian prairies. This study examined the immediate effects of elevated temperature and defoliation on plant community diversity, composition, and biomass within a native rough fescue (Festuca hallii (Vasey) Piper) grassland over two growing seasons. We used open-top chambers to simulate climate change and defoliated vegetation in midsummer of the first year to simulate biomass loss associated with periodic ungulate grazing. Warming marginally increased plant species diversity and changed community composition shortly after treatment, but effects were not apparent the second year, and there were no apparent impacts on plant biomass. Nonetheless, warming may have impacted community diversity indirectly through reduced soil moisture content, a pattern that persisted into the second year. Overall, this northern temperate grassland demonstrated limited community-level changes to warming even in the presence of defoliation.

1. Introduction

Over long time scales, the plant community composition within savanna habitats is sensitive to climate change [1, 2]. Over the past 50 years, western Canadian savannas, such as the Aspen Parkland ecoregion, have experienced an increase in mean annual temperature of approximately 1.3°C [3]. Mean temperatures for this region are predicted to further increase over the next century, leading to higher evaporation rates and reduced soil moisture availabilities [4, 5], as well as a longer frost-free period [6]. The consequences for native vegetation within these regions are unclear.

Under anticipated climate change scenarios, native savanna ecosystems, existing as a mosaic of community types, may be particularly prone to ecological change [7]. Moreover, land-use practices such as livestock grazing may exacerbate any climate-driven changes in vegetation [8]. Plant community changes arise when species vary in their fundamental response to environmental and management factors, including temperature and defoliation in the form of grazing. Warming can impact plants both directly by altering the potential for photosynthesis throughout the growing season [9] or indirectly by altering complex secondary biological processes in the ecosystem [10]. Collectively, warming and defoliation may lead to reductions in the diversity of native species or the invasion of undesirable species [11].

The Aspen Parkland of western Canada is historically a distinct savanna habitat comprised of a mosaic of forest and grassland communities. Currently, few native grasslands persist in the region [12] due to human land-use practices, making ecosystems in this region one of the most threatened in western Canada. Plant communities in the Aspen Parkland are known to have varied greatly over time [1], presumably due to the marked impact of ongoing fluctuations in climate coupled with secondary disturbance [13]. Plant growth in this region is typically limited by low temperatures during winter and spring and by water stress during late summer [5], but recent projections suggest that the Aspen Parkland may expand by as much as 21% due to climate change [14].

Grazing by large ungulates is both historically and currently a common disturbance in grasslands of the Aspen Parkland [15, 16]. Removal of plant biomass may exacerbate soil warming through reduced shading of the soil surface and altered soil moisture availability [17]. Moreover, warming and defoliation induced changes to plant growth and associated litter biomass may further impact plant community composition, including plant species diversity [18]. Ultimately, there is reason to suspect that nonadditive effects may exist between the disturbances of warming and defoliation, with the magnitude of final changes dependent on whether plant species adapted to one stressor are adapted to the second.

Climate change and the associated changes in growing conditions (e.g., moisture availability) have been widely implicated for their potential to negatively affect agricultural ecosystems [19]. Thus, investigation of the impact of elevated temperatures in conjunction with grazing as the prevailing land use will improve our understanding of the response of these ecosystems to climate change. To understand and mitigate the potential impacts of climate change on the biodiversity and sustainable production of Canada’s remaining native Parkland grasslands, a mechanistic understanding is needed of the relationship between temperature and defoliation within these areas.

We examined the immediate impacts of warming, with and without defoliation, on plant community composition, diversity, and biomass, within a native rough fescue grassland in the Aspen Parkland of Alberta, Canada. Moreover, we predicted that defoliation would augment the impacts of warming with the removal of litter, as the latter is known to be an important driver of plant growth in these grasslands [20]. We also examined how plant community characteristics were associated with key abiotic environmental variables, including light availability, soil moisture content, soil nitrogen mineralization rates, and soil temperature.

2. Materials and Methods

2.1. Field Site

We conducted a field experiment within a native rough fescue grassland in the Aspen Parkland natural subregion of Alberta, Canada, in 2006 and 2007. The study was conducted in a 20 ha area at the University of Alberta Kinsella Research Station (53°05′N; 111°33′W), situated 150 km southeast of Edmonton. The area climate is dry subhumid, with mean annual precipitation of 431 mm and mean annual temperature of 1.5°C [20]. Mean annual growing season (May–September) precipitation over the past 30 years was 312 mm, with mean annual growing season temperature of 13.9°C. In 2006 and 2007, annual growing season precipitation was 269 and 287 mm, respectively; although below the 30-year norm, these conditions are not unlike those of the last decade during which growing season precipitation averaged 284 mm [21]. Mean annual growing season temperature was 16.1°C in 2006, 2.2°C warmer than in 2007 (13.9°C).

Dominant grasses at the study site included Hesperostipa curtiseta (A.S. Hitchc) Barkworth, Festuca hallii (Vasey) Piper, Koeleria macrantha (Ledeb.) Schult., and dryland Carex species. Dominant forbs included Achillea millefolium L., Commandra umbellata L., Aster falcatus (Lindl.), Artemisia ludoviciana (Nutt.), and Solidago missouriensis (Nutt.). Rosa arkensana (Porter) was also prevalent [22]. A more detailed description is available in Lamb [18].

Soils at the study site are Dark Brown Chernozems and are generally well drained. On average, soils are loamy in texture (22% clay, 33% silt, and 45% clay) have a pH of 5.7, Ah depth of 13.1 cm, and approximately 13% organic matter. Although the dominant native vegetation found at this site has been characterized as deep rooted [23], most root biomass in the study area remains situated within 20 cm of the ground surface as rocks and hard packed clay limit the ability of roots to reach lower depths [24]. While the study area had a long history of fall grazing by cattle, it had not been grazed for 3 years prior to the start of this study.

2.2. Experimental Design

The experiment was established in a randomized complete block design with five replicate blocks each 5×5 m in size. Each block was more homogeneous with respect to soil type, slope, aspect, and drainage than between blocks. Within each block, four 2 m diameter circular plots were established and four treatment combinations were randomly assigned. Treatments included a control, warming, defoliation, and warming + defoliation.

We used open-top chambers (OTCs) situated directly over the plots to increase air temperatures. The OTCs are a low-cost, proven method to increase air temperature in field studies [25] with previous studies demonstrating an increase of 1 to 2°C in average daily temperature [26, 27]. Our OTCs were 40 cm high and cone shaped, with 2 m diameter at the bottom and 1.6 m diameter at the top (Figure 1). The OTCs were made from Sunlite HP solar glazing material (Solar Components Corporation, Manchester, NH, USA), which allows visible light transmission but reradiates light in the infrared range, creating a greenhouse effect within the OTC. OTCs were placed on the plots in April right after snowmelt and removed in October after the current years growth had senesced. Initial data collected in June 2006 indicated that this design provided an increase of up to 4°C in air temperature near midday, although average daytime temperature increases were closer to 1.6°C.

Defoliation was achieved through manual clipping of all standing vegetation in and around each plot (including a 50 cm nonsampled buffer). Clipping was done in mid-June 2006 to a stubble height of 2.5 cm above ground and is consistent in timing and intensity with the common land-use practice of midsummer cattle grazing across the region. Although actual ungulate grazing would have been more representative of the disturbance these communities experience, ungulates could not be used due to their tendency to feed in a spatially heterogeneous (i.e., patchy) manner as well as pose risk of damaging OTCs. Defoliation was not done in 2007 in order to examine the residual impacts of defoliation during the previous year and assess community recovery from grazing under both warmed and unwarmed conditions.

2.3. Vegetation Sampling

Plant community composition in each plot was assessed in early August 2006 and 2007 within a centrally located 0.5×0.5 m quadrat. Percent canopy cover for each species was visually assessed in each quadrat, from which we calculated plant species richness (number of species) and Shannon’s diversity index as 𝐻=Σ(𝑃𝑖ln𝑃𝑖),(1) where 𝑃𝑖 is the proportion of total species cover comprised of each individual species “𝑖.” Shannon’s diversity was used as it accounts for both the relative abundance and density of plant species in a plot.

Destructive biomass sampling occurred within a separate 0.25×0.5 m area within each plot at peak biomass (early August) of each year, within which total live shoot biomass was sorted from plant litter. Root cores (5 cm wide × 20 cm deep) were also removed from each quadrat at the time of harvest to assess root biomass at the community level [28]. Only one core was removed from each plot due to limited space for destructive sampling. Cores were washed free of soil through a 0.5 mm sieve. Both above and below ground vegetation samples were dried at 65°C to constant mass and weighed.

2.4. Environmental Sampling

To determine the predominant environmental factors that were influenced by our experimental treatments and associated with community composition, we measured several parameters during the growing season. Air temperatures were recorded every 10 minutes in the summer of 2006 within warmed and unwarmed plots using Boxcar HOBO data loggers (Hobo H8 Pro Series, Onset Computer Co. Procasset, Mass, USA) installed in PVC radiation shields 2 cm above ground level. Temperatures were averaged daily, as well as from 0600 until 1800 hr to assess daytime air temperatures within each plot. Soil moisture and temperature readings were taken 4–7 days after rainfall in both June and July, near midday to assess treatment-induced differences in the microenvironment under field conditions. Moisture content was measured in the top 10 cm of mineral soil using a TDR-MLX2 moisture probe (Delta-T Devices, Cambridge, UK), with a minimum of 4 replicate measurements per plot in each of 2006 and 2007: measurements were then averaged to obtain a per plot estimate of moisture content.

Light interception by vegetation was measured coincident with soil moisture content using a Li-Cor ceptometer (Li-Cor Inc, Lincoln, Neb, USA) and reported as photosynthetically active radiation (PAR). Values of PAR were measured above and below the vegetation canopy using an average of 10 segmented sensors to assess light penetration to the soil surface. Readings were taken in the center of the plot and outside the shading influence of the OTC. The proportion of light intercepted was then calculated using the following equation: (abovecanopyPARbelowcanopyPAR)abovecanopyPAR100%.(2)

In order to assess in situ net nitrogen mineralization rates, the buried-bag method [29] was used. Duplicate soil cores (5×15 cm) were taken monthly at 0–5 cm depth and combined for subsequent analysis. Additionally, 𝑁 supply rates were measured monthly during 2007 using Plant Root Simulator (PRS) probes (Western Ag Innovations Inc, Saskatoon, SK, Canada). Detailed methods are provided by Attaeian [30].

2.5. Data Analysis

To examine the impact of warming on the microenvironment of plots prior to plant responses, we conducted separate general linear mixed models using the June 2006 (the start of the experiment) measures of soil moisture and light interception as response variables. In each analysis, the presence or absence of an OTC served as the fixed effect, with block as a random factor. Soil temperature was tested in July 2006 after soil temperatures were assumed to have equilibrated with the presence of the OTCs. Analyses were conducted using Proc Mixed in SAS 9.1.3 (SAS Institute Inc., Cary, NC, USA).

To test the impact of warming and defoliation on plant community structure, a permutation-based MANOVA (perMANOVA) [31] and a multi-response permutation procedure (MRPP) [32] were run using a Sorenson distance measure on the community composition data for August 2006 and 2007 in PC-ORD version 5 (MJM Software Design, Gleneden Beach, Ore, USA). For the perMANOVA, 4999 randomizations were used and significance was based on the proportion of randomized trials with an indicator cover greater than or equal to the observed cover value, such that 𝑃=(1+numberofrunsobserved)(1+numberofrandomizedruns).(3) For the MRPP, significance was based on the probability of the observed delta statistic being smaller or equal to the expected delta.

To further examine plant species responses to the warming treatment, a nonmetric multidimensional scaling (NMS) [33, 34] ordination was run on vascular plant cover data from early August in both 2006 and 2007 using a Sorensen distance measure in PC-ORD. Real data were run 250 times, as were the randomized data for the Monte Carlo test. A total of 500 iterations were used to obtain the final stable solution with instability of 0.00001. Axes scores were interpreted based on Pearson and Kendall correlations with all species found across the site. An indicator species analysis [35] using 4999 permutations was also performed on the August 2006 community data to examine specific plant species responses to warming.

To test the impact of warming and defoliation on overall plant species diversity and biomass, general linear mixed models were run on August 2006 and 2007 measures of plant species richness, diversity, and shoot, litter, and root biomass. We used an alpha of 5% to indicate significance and minimize the potential of a type II error, although relationships with an alpha of 10% were also examined.

Finally, to identify abiotic variables associated with variation in community composition, Shannon’s diversity (plant species), and biomass, we performed a series of multiple regression analyses that included soil temperature and moisture content, 𝑁 mineralization rates, and light interception from the preceding two months (June and July) within each year as independent factors. A stepwise multiple regression was used with the minimum significance level of entry into the model set at 10%.

3. Results

3.1. Microenvironment and Relation to Vegetation

The OTCs had the anticipated effects on abiotic conditions, increasing daily surface air temperatures by 1.0°C and average daytime air temperatures by 1.6°C during 2006. Warming did not impact soil temperature in July 2006 (𝐹1,12=0.12, 𝑃=0.73), but showed a potential impact on June 2006 soil moisture content (𝐹1,14=3.78, 𝑃=0.07), with warmed plots drier (24.7±1.4%) than unwarmed plots (26.7±1.4%). There was no evidence that the OTCs altered light interception by vegetation in June 2006 (𝐹1,14=0.97, 𝑃=0.34), with warmed plots (31.8±2.7%) similar to unwarmed plots (28.9±2.7%).

Environmental variables showed a minor relationship with vegetation during both years of the study, with the strongest association between soil moisture content and the plant community, particularly Shannon’s diversity and richness (Table 1). Additionally, shoot biomass was significantly related to light interception (Table 1).

3.2. Influence of Warming

Our perMANOVA analyses of plant community composition revealed a weak effect of warming in 2006 (perMANOVA pseudo-𝐹=2.10, 𝑃=0.07; MRPP 𝐴=0.021, 𝑃=0.12). The 2006 NMS also indicated an effect of warming on community composition, however, with Axis 1 showing evidence of groupings of plots based on the warming treatment (Figure 2(a)). NMS Axis 1 was closely associated with Hesperostipa curtiseta (𝑟=0.82), Koeleria macrantha (𝑟=0.60), and Artemisia frigida (𝑟=0.54) (𝑃<0.05), with results further confirmed in the indicator species analysis through the response of K. macrantha (𝑃=0.05) and A. frigida (𝑃=0.10) to warming. The 2006 NMS analysis resulted in a two-dimensional final ordination solution (Axis 1: 𝑃=0.004, Axis 2: 𝑃=0.004), with final stress of 10.34. The proportion of species variance represented by Axis 1 was 65.6% and by Axis 2 was 24.1%, based on the correlations between ordination distances and distances in the original 𝑛-dimensional space.

By 2007, there was no evidence of any lasting effect of warming on community composition, regardless of the methods of analysis used (perMANOVA pseudo-𝐹=0.91, 𝑃=0.52; MRPP 𝐴=0.006, 𝑃=0.64). The 2007 ordination with NMS also showed no distinct patterns associated with the warming treatment (Figure 2(b)). The 2007 NMS resulted in a three-dimensional final ordination solution (all axes 𝑃=0.004), with final stress of 10.26. Axis 1 represented 27.2% of the variance, while axis 2 represented 5.6% and Axis 3, 56.2%.

Direct impacts of warming were found through a general linear mixed model on Shannon’s species diversity, but only in 2006 (𝐹1,12=4.68, 𝑃=0.05), during which warmed plots (𝐻=1.91±0.11) were more diverse than unwarmed plots (𝐻=1.69±0.11). In contrast, warming showed no impact on species richness, shoot biomass, root biomass, or litter biomass (Table 2).

3.3. Influence of Defoliation

Defoliation strongly affected plant community composition in 2006 (perMANOVA pseudo 𝐹=6.62, 𝑃=0.0006), though these effects also did not carry over into 2007 (pseudo 𝐹=0.78, 𝑃=0.64). The lack of a warming by defoliation interaction in 2006 (pseudo 𝐹=0.55, 𝑃=0.77) and 2007 (pseudo 𝐹=0.82, 𝑃=0.62) indicated that the effects of defoliation on community composition were not dependent on warming.

Defoliation impacted species diversity and peak shoot biomass in 2006, as well as peak shoot and litter biomass in 2007 (Table 2). Species diversity was lower in the defoliated plots (𝐻=1.61±0.12) compared to the nondefoliated plots (𝐻=2.00±0.11). Total (accumulated) season-long shoot biomass was 25% greater in defoliated plots (259.2±20.3 g m−2) compared to non-defoliated plots (206.7±20.3 g m−2) in 2006. However, defoliated plots remained 19% lower in shoot biomass relative to non-defoliated plots in 2007 (199.6±17 versus 246.1±17 g m−2, resp.), the year of no defoliation. Despite this, no interaction of warming with defoliation was found on any of the response variables, including species richness, peak shoot biomass, peak root biomass, and litter biomass in either year of the study (Table 2), again indicating that warming did not alter the response of plants to the defoliation treatment.

4. Discussion

Warmer conditions produced by the OTCs (~1.6°C in daytime) were similar to the increases in mean annual air temperature seen in western Canada over the past 50 years [4, 36] and reinforce the utility of this method of imposing a modest increase in air temperature for examining the biological consequences of climate change. The observed reduction in soil moisture availability associated with the warming treatment suggests either that the OTCs indirectly reduced ambient rainfall into plots or, more likely, that the modest increase in temperature increased evapo-transpiration. While limited in magnitude, the reduction in measured soil moisture availability (by 5.4%) provides support to the notion that a warmer microenvironment may increase the potential for water stress in this northern temperate grassland and could also explain why some climate change models [14] predict a shift northward in the boundaries of the Aspen Parkland. However, this observation must be tempered by the fact that rainfall was below normal during the study period and must therefore be interpreted in this context. Although it is important to recognize that temperature fluctuations, both seasonally and interannually, are typical of most grasslands, to which this is no exception, making it possible that the increase here was within the normal range of variation for this ecosystem, growing conditions during the two study years examined here remained consistent with recent climatic trends for the Aspen Parkland, which has experienced an increased frequency of summer drought during the 1990s [5].

Warming had a limited impact on plant community composition immediately following its initiation, but impacts were not apparent in the second year of study. The marginal increase in diversity rather than richness under warming in 2006 suggests that the combination of increased temperature and lower soil moisture favored the development of a more even plant community rather than a greater number of species. Richness is often slower to respond than Shannon’s diversity, as the former entails a change in the absolute presence of individual plant species. Several plant species, including H. curtiseta, K. macrantha, and A. frigida, showed evidence of an initial positive response to warming in our study. Notably, these species are common dominants of grasslands across the Mixed-grass Prairie region to the south of the Aspen Parkland [37] and a shift in dominance to such species could suggest a potential change in prairie structure under warmer and potentially drier conditions. Niu and Wan [38] documented changes in plant competitive hierarchy arising from differential responses to warming in northern China and may account for some of the short term changes observed here.

Plant shoot and root biomass representing the entire plant community appeared resistant to the influence of warming in the short-term. Given that changes in aggregate community biomass are unlikely to be seen until reductions in vigor occur within plants, changes in biomass may be expected to take much longer than two years to develop. Nonetheless, this does not discount potential changes in biomass of individual species. For example, working in wet meadows of northern Alaska, Hollister and Flaherty [39] found an increase in sedge biomass under 1-2°C warming, but similar to our findings total community biomass did not change. In the present study, warming also failed to interact with defoliation to alter plant composition or biomass production, suggesting that warming induced by the OTCs did not influence tolerance of this community to defoliation, including recovery the following year. In general, defoliation had a more acute and consistent impact on the plant community, leading to immediate reductions in diversity as well as extended declines in biomass production.

In general, few studies exist to directly compare our results to, particularly given that no studies have examined the direct effects of warming on northern temperate grasslands dominated by cool season plant species. Nevertheless, our results indicate that defoliation (i.e., the predominant land use in the area) has a greater immediate impact on Parkland grasslands than warming and is consistent with projections of global biodiversity responses for terrestrial ecosystems by others [11]. These same authors, however, suggested that northern temperate ecosystems may experience the least biodiversity change in the future, largely due to the acute changes already produced by previous changes in land use activities [11]. Our results provide at least some evidence that contrasts this perspective, as warming of Parkland grasslands did induce changes in species composition, albeit minor and potentially short lived. Interpretation of defoliation impacts must also be placed in the context of the clipping treatments imposed rather than the use of grazing animals. As animals are selective during defoliation and have significant impact on soils (e.g., through compaction and nutrient cycling), the impacts of large herbivores under warming may differ from the current results in these landscapes. More research is needed on these questions in this region of North America.

Long-term studies are recommended to more definitively establish the mechanisms behind the observed changes in species composition, particularly those xeric adapted plant species. For example, community dynamics in this study were associated primarily with light and soil moisture availabilities rather than soil temperature. This observation suggests that the predominant driver of climate-induced vegetation change within these ecotonal grasslands may be soil moisture availability, which, in turn, may be indirectly altered by thermodynamics within the plant canopy and rhizosphere. Plant litter accumulation has been found elsewhere to have an overriding impact on soil moisture content [40] as well as associated diversity responses [18] within Parkland grasslands. Moreover, although species richness had no direct response to our experimental treatments, it is worth noting that richness was positively related to soil moisture content in both years, suggesting that long-term declines in effective moisture availability, brought about either through reduced precipitation or increased evapo-transpiration, may lead to a parallel loss of plant species. This finding is supported by other studies documenting strong effects of interannual variation in precipitation on plant communities [41].

Not surprisingly, light interception was most closely associated with Axis 1 in the NMS from 2006, and directly reflects the prominent impact of leaf area removal during defoliation on the plant community during that year. The negative relationship of diversity with 𝑁 mineralization in 2006 suggests that greater nutrient cycling may lead to a less uniform plant community, potentially due to suppression of those species less adapted to increased 𝑁 availability and an increase in those with favorable responses to 𝑁. Notably, these results support the work by Dijkstra et al. [42] who documented a more open 𝑁 cycle (i.e., increased soil inorganic 𝑁 and total plant 𝑁) in a semiarid grassland under warming. Root biomass was also positively associated with 𝑁 mineralization in 2006 and may reflect opportunistic root foraging by plant species adapted to high 𝑁 availability. A common invasive grass of the area, Poa pratensis L., is well adapted to 𝑁 addition [43]. Finally, this response during 2006 may also be tied to defoliation in that year, as this disturbance is known to lead to root death of cool temperate grasses in the region [44], which, in turn, would lead to mineralization and the aforementioned increase in nitrophilic species.

5. Conclusions

Despite the short-term nature of this ecological investigation, our results provide useful preliminary information on the initial response of native northern temperate Parkland grasslands to warming and defoliation. Although these grasslands demonstrated prompt but minor changes in diversity and plant species composition in response to warming, these changes disappeared in the second year and were independent of defoliation. Instead, most plant community and biomass responses in this study appeared to depend directly on defoliation and changes to soil moisture availability, light interception, and 𝑁 mineralization, rather than soil temperature. This grassland ecosystem showed resilience to short-term fluctuations in mean annual air temperatures under normal grazing intensities, though long-term impacts require further study.

Acknowledgments

This study was funded by the Natural Sciences and Engineering Research Council of Canada (Discovery Grant to J. F. Cahill Jr. and CGS-M scholarship to E. S. Deutsch), the Alberta Conservation Association (Biodiversity Challenge Grants Program), and the University of Alberta (Range Management Post-Graduate Endowment). The authors thank Behnaz Attaeian for assistance in setting up the study and for provision of the nitrogen mineralization data. They also thank Jennifer Caudron for field assistance.