Arctic-Alpine Plants Decline over Two Decades in Glacier National Park, Montana, U.S.A.

Abstract

Climate-based models predict global warming will cause mountain plants to migrate upward with local extirpation of species currently restricted to the highest elevations. Alpine monitoring studies have generally documented changes in alpine plant distributions consistent with these predictions; however, relatively few such studies have been reported from North America. I estimated canopy cover in 71 permanent microplots at two alpine moist-turf sites in Glacier National Park, Montana, U.S.A., three times between 1988 and 2011. Mean annual and summer temperatures were approximately 0.6 and 0.7 °C higher than the previous four decades, respectively. Plants more restricted to high elevations declined more than species with a broader elevational amplitude during this time. Dicots were more likely to have declined than monocots. These data support predictions and provide information that may help refine climate-envelope models.

Introduction

There is widespread agreement that increasing levels of atmospheric carbon dioxide and other greenhouse gases will result in global warming during the next 50 years (Henderson-Sellers, 1990; IPCC, 2007). Less rapid changes of comparable amplitude in the past altered the ranges of many species (Graham and Grimm, 1990; Webb, 1992; Jackson and Overpeck, 2000), and there is evidence that biota are currently shifting their geographic ranges upward in elevation or latitude (Root et al., 2003). High-mountain ecosystems are thought to be particularly vulnerable to the effects of warming (Grabherr et al., 1995; Cannone et al., 2007; Nagy and Grabherr, 2009; Engler et al., 2011), especially at higher latitudes (Körner, 1999). In general, climate envelope models that employ the relation between current climate variables and present-day species distributions predict that species will migrate upward with local extirpation of species currently restricted to the highest elevations (Gottfried et al., 1998; Walker et al., 2006). However, there are reasons to believe that some alpine species will persist by shifting microsites rather than elevational zones (Scherrer and Körner, 2010).

Many monitoring studies have documented changes in alpine plant distributions consistent with climate-warming predictions. The majority of these studies are based on presence/absence data in which a study area is surveyed for the presence of species at one point in time and then resurveyed at a subsequent time (Grabherr et al., 1994; Klanderud and Birks, 2003; Walther et al., 2005; Parolo and Rossi, 2008; Erschbamer et al., 2012). In most cases, researchers have documented increases in species richness resulting from upward migration of lower-elevation species but have only infrequently found a loss of high-elevation species. Presence/absence measurements can be more repeatable than canopy cover estimates (Elzinga et al., 2001), but change is more likely to be detected earlier by accurate measures of canopy cover. Canopy-cover monitoring of alpine vegetation was begun on a large scale in Europe in 1994. One ten-year study in the Austrian Alps that monitored plant cover in permanent plots found an increase in the cover of “low-alpine” (= alpine) species but also documented a decline in “high-alpine” (i.e., nival and subnival) species (Pauli et al., 2007).

Several studies have indicated that a species' response to climatic warming will often be associated with its functional or life-form group (Arft et al., 1999). For example, experimental warming studies suggest that shrubs will increase at lower alpine zones in Europe (Cannone et al., 2007), Alaska (Epstein et al., 2000; Walker et al., 2006), and Asia (Klein et al., 2007) at the expense of graminoids. Unfortunately, results of long-term monitoring studies have rarely been examined for life-form correlations. Although a good deal of experimental work has been done in North America, there have been relatively few published, long-term monitoring studies of alpine vegetation (Walker et al., 1994). Here I report on a 24-year monitoring study of moist alpine turf in Glacier National Park in the Rocky Mountains of western North America. Canopy cover data are analyzed to determine if there is a signal similar to those found in Europe and which species and life-form groups have increased or declined.

Study Sites

Permanent monitoring transects were established at two locations above timberline in Glacier National Park in the Logan Pass area: Lunch Creek and Mount Reynolds (Fig. 1). Both sites are on gentle south-facing slopes, subirrigated throughout at least part of the growing season by snowmelt originating upslope. These unusual sites mimic arctic tundra by having a long growing season due to the south exposure but perennially cold and moist soil. These moist turf sites harbor several arctic plant species that are at the southern margin of their range (Lesica and Steele, 1996; Lesica and McCune, 2004). Parent materials of both sites are limestones and argillites of the Precambrian Belt Series, and both sites have coarse to fine colluvium overlain by highly organic soil of varying thickness (sometimes absent) that is saturated or nearly so for much of the growing season. The Mount Reynolds site is on the south-facing side of Mount Reynolds, 4 km SSW of the Logan Pass visitor center at an elevation of 2285 m. There is a late-season snowfield at the base of the near-vertical slope ∼50 m above the transects. The Lunch Creek site is in Lunch Creek Cirque in the valley created by Polock Mountain and Piegan Mountain, 2 km NNE of the Logan Pass visitor center at an elevation of 2165 m. There are permanent snowfields ∼150 m farther up the cirque.

Climate data are available for several sites at both high and low elevations in Glacier National Park, but only discontinuous data are available from high-elevation recording stations. Annual precipitation at Sperry Glacier, at 2440 m and 5–10 km south of the study sites exceeds 250 cm. Precipitation occurs as rainfall during June through September and as snowfall the remainder of the year (Finklin, 1986). Summer climate is cool and moist; during the 1959 growing season, average daily temperature and relative humidity were 10 °C and 74% (Choate and Habeck, 1967). Mean annual daily high and low temperatures were estimated to be 6.6 and -4.9 °C, respectively, between 1980 and 2003, with a mean annual precipitation of 159 cm (Malanson et al., 2012). Continuous long-term (i.e., >50 yr) climatic data have been collected only at West Glacier at 975 m. Climate at West Glacier over the two-decade course of the study was warmer but not wetter than the previous four decades (Table 1). Snowpack across the Rocky Mountains has been melting earlier over the past three decades due to warmer spring temperatures (Pederson et al., 2011), suggesting that climate trends at West Glacier mirror those at higher elevations.

FIGURE 1. Location of Lunch Creek and Mount Reynolds study areas in Glacier National Park, Montana, U.S.A. Aerial photograph was taken in late summer, 1990. Note the snowfields just above the two sites.

TABLE 1 Mean annual and summer (June through September) temperature and mean annual precipitation at West Glacier.

Methods

FIELD METHODS

Two (Mount Reynolds) or three (Lunch Creek) permanent transects were established at the study sites in 1988. Transects were subjectively located to pass through the middle of longitudinal strips of moist turf habitat running perpendicular to the slope. I estimated canopy cover for all vascular plant species in 0.1 m² plots (50 × 20 cm) placed 1 m apart along the transect lines. There was a total of 35 and 36 sample plots at Mount Reynolds and Lunch Creek sites, respectively. Canopy cover (Daubenmire, 1959) was estimated to the nearest 5%, but estimates less than 2% cover were recorded as 1%. Canopy cover of vascular plants was recorded at Mount Reynolds in 1988, 2000, 2011, and at Lunch Creek in 1988, 2000, and 2010. There is a subjective aspect to canopy cover estimation (Elzinga et al., 2001) that can introduce unwanted error when data collected by different observers are compared. All canopy-cover estimates were made by the author for all three recording periods. Although the estimates may not be perfectly accurate, the change between estimates at different times by the same observer will usually result in acceptable precision (Elzinga et al., 2001). Vascular plant nomenclature follows Lesica (2002).

DATA ANALYSIS

I used repeated-measures analysis of variance (ANOVA) followed by a linear polynomial contrast test to ascertain the presence of trend for each of the common species (recorded in >10 sample plots). Plots were considered independent samples. The data met requirements for normality of residuals and homogeneity of variance. Site was used as a factor in the ANOVA model if a species occurred at both Lunch Creek and Mount Reynolds, and sites were analyzed separately if there was a significant Year × Site interaction (P ≤ 0.05).

For each species, I used the slope of the best-fit linear regression model of canopy cover across the three years of the study (see previous section) divided by the mean canopy cover for the three years as a measure of trend. This adjusted regression slope (hereafter “trend index”) allowed me to compare trend for common species with those that were less common. I used a recent floristic treatment for Montana (Lesica, 2012) to assign an “alpine index number” to each of the 33 common species found at the study sites. This index number reflects the degree to which each species is restricted to high elevations in Montana. Species that occur only above treeline were assigned an alpine index of (4), both alpine and subalpine (3), montane through alpine (2), and valleys through alpine (1). I used linear regression analysis to assess the relationship between the trend index and the alpine index number, that is, whether species restricted to high elevations declined more or increased less than less restricted species. I used a two-sample t-test to determine whether adjusted regression slope (trend) differed between monocots and dicots over the course of the study.

Results

Vegetation at the Lunch Creek and Mount Reynolds sites was similar. The plant community at Lunch Creek was dominated by the dwarf shrubs Salix reticulata (mean canopy cover across all 3 years = 11%), and Dryas octopetala (9%). Carex scirpoidea (18%) and C. capillaris (7%) were the dominant graminoids. The most common forbs were Solidago multiradiata (7%), Hedysarum sulphurescens (7%), Polygonum viviparum (6%), Silene acaulis (5%), and Anemone parviflora (3%). Tundra at the Mount Reynolds site was also dominated by Dryas octopetala (26%) and Salix nivalis (8%). Carex scirpoidea (18%), Kobresia simpliciuscula (17%), and C. rupestris (5%) were the dominant graminoids. The most common forbs were Polygonum viviparum (9%), Silene acaulis (3%), and Solidago multiradiata (3%). The cover of bare rock was 5% at the Lunch Creek site and 3% at Mount Reynolds.

A total of 63 vascular plant species were recorded at the Lunch Creek and Mount Reynolds sites over the course of the study. Of these, 33 species were recorded in more than ten plots (Table 2). Canopy cover of Carex scirpoidea, Potentilla fruticosa, and Sedum rosea (at Lunch Creek) showed statistically significant increases over the course of the study, while cover of Gentianella propinqua, Kobresia simpliciuscula (at Mount Reynolds), Polygonum viviparum, Salix arctica, S. reticulata, and Sedum rosea (at Mount Reynolds) declined significantly (Table 2). All six of the significantly declining species, but only one of the three increasing species, are restricted to alpine and subalpine habitats in Montana (Table 2).

There was a marginally significant negative association between alpine index and trend index indicating a tendency for species restricted to higher elevations (greater alpine index) to have declined more in canopy cover (more negative trend index) than species with greater ecological amplitude over the two decades of the study (R ² = 0.09, P = 0.089, Table 2). In addition, canopy cover of the 21 dicots declined at a greater rate compared to the 11 monocots (t = 2.1, P = 0.045, Table 2). Only three species of woody plants were common enough for analysis. All three were in the genus Salix (Table 2). The two dwarf shrubs (S. arctica and S. reticulata) both declined, while the taller S. vestita remained stable.

Discussion

There was a declining trend in the canopy cover of species restricted to high elevations compared to species with a broader elevational amplitude over the past two decades in the moist-turf study sites in Glacier National Park. During this same period, mean annual daily and summer daily temperatures increased 0.6 and 0.7 °C, respectively, at the West Glacier recording station compared to the previous 40 years, while annual precipitation was virtually unchanged. Results from this study are descriptive in nature rather than experimental, so I cannot rigorously infer that climatic warming caused the observed changes. However, my results are consistent with forecasts for alpine areas in Europe and elsewhere (Gottfried et al., 2002; Nagy and Grabherr, 2009).

Although the observed changes in species canopy cover were associated with increased temperature rather than declines in precipitation, I believe that changes in snowpack mediated by increased temperatures are the most likely cause of the observed vegetation shifts at my study sites because the duration of snowmelt and disposition of meltwater strongly affect the pattern of alpine plant communities (Walker et al., 1993; Malanson et al., 2012). Tundra plants often have detrimentally high base metabolic rates in abnormally warm environments (Billings and Mooney, 1968; Spomer and Salisbury, 1968; McNulty et al., 1988; Körner, 1999, p. 192), and there is evidence that species from cold, moist sites have less ability to adapt to warming temperatures (Cooper, 2004). Moist tundra vegetation at the two study sites is subirrigated by cold meltwater, which presumably keeps the roots and, to some extent, the above-ground vegetation cooler than ambient temperatures. Warmer temperatures in spring and summer in Glacier National Park have resulted in a diminution of summer snowpack since the mid-1980s (Pederson et al., 2010; P. Lesica, observation), presumably leading to drier and warmer soil downslope. I suggest that the warmer and drier soil of these sites later in the growing season is causing a loss of species adapted to cold, moist habitats in favor of those adapted to warmer, drier, lower-elevation environments.

Studies of plant distributions in alpine sites have frequently documented upward migration of low-elevation species but less commonly declines in high-elevation species (Grabherr et al., 1995; Cannone et al., 2007; Nagy and Grabherr, 2009; Engler et al., 2011), while my study suggests both processes are occurring. Most alpine dominants are long-lived perennials, so detecting a decline of species will take a longer time than detecting immigration of lower-elevation plants (Grabherr, 2003). Furthermore, many monitoring studies have been conducted on exposed alpine summits with drier soils and hydrology that is less affected by persistent snowfields. Alpine plants in these environments are expected to be less sensitive to drying, so the main effect on the plants will be a longer growing season and warmer ambient temperatures that allow immigration by lower-elevation species. Most alpine monitoring has been done using presence/absence data from large plots or zones (Grabherr et al., 1994; Klanderud and Birks, 2003; Walther et al., 2005; Parolo and Rossi, 2008; Erschbamer et al., 2012). Monitoring canopy cover in numerous small sample plots will usually be more sensitive than monitoring presence/absence over large areas because extirpation in a large area will take more time than a diminution of cover. Pauli et al. (2007), who also measured canopy cover in the Austrian Alps, documented an increase in low-elevation species but also a decline in high-elevation species, similar to my results.

TABLE 2 Common species (# of plots > 10) at Lunch Creek and Mount Reynolds sites and their mean canopy cover over the course of the study, trend index (see Methods), alpine index (see Methods), and significant change as determined by repeated-measures ANOVA followed by a linear polynomial contrast test.

Some researchers have suggested that the response of arctic species to climatic warming will depend on the plant's growth form (i.e., functional group) (Chapin and Shaver, 1996; Chapin et al., 1996; Arft et al., 1999; Yang et al., 2011; but see Klanderud, 2008). Common alpine graminoids increased during the warmest decade on record in the Austrian Alps (Bahn and Körner, 2003), and Walker et al. (2006) found that arctic graminoids increased under experimental warming. My study suggests that the same may be true for alpine plant communities. Populations of graminoids and other monocots (5 of 11 species in the Cyperaceae) tended to increase or remain stable more than broad-leaved forbs and shrubs over the course of the study. Greater resilience to drier conditions by monocots compared to dicots may be the result of anatomical differences such as structural root mass or size of conductive tissues (Dickison, 2000; Carlquist, 2012) or differences in physiology (Edwards et al., 2007). Similar observations were made on Great Plains grasslands during the great drought of the 1930s. Weaver (1968) reported that drought caused only a shift in the dominance of the grass species, but resulted in the disappearance of many species of forbs (Albertson and Weaver, 1944).

Climate-change experimental research on arctic plants also suggests that shrubs will increase with a warming climate, but dwarf (i.e., prostrate) shrubs are likely to decline (Klanderud and Totland, 2005; Wipf et al., 2005; Walker et al., 2006; but see Bahn and Körner, 2003). Although the sample size is small, my results are consistent with these reports. The two prostrate willow species (Salix arctica, S. reticulata) declined significantly, but there was no detectable change in the cover of the taller shrub, Salix vestita (Table 2), which grows 20–100 cm tall in Montana (Lesica, 2012).

Conclusion

The results of my study support the general prediction for alpine vegetation that lower-elevation species will expand their range, while species restricted to the highest elevations will decline in abundance with climatic warming. Furthermore, it appears that dicots may be more prone to decline than monocots. Finally, four of the common recorded species, Pinguicula vulgaris, Tofieldia pusilla, Kobresia simpliciuscula, and Gentianella propinqua, are at or near the southern edge of their arctic-alpine distribution. Of these, the latter two (50%) declined compared to only 18% of the common species overall, suggesting that peripheral populations may also be at greater risk from climate change (Lesica and McCune, 2004). These results may help make climate-based models more biologically realistic and accurate in predicting the fate of alpine plants on a warming planet.

Acknowledgments

Kathy Ahlenslager, Teri Allendorf, Jennifer Asebrook, Levi Besaw, Tara Carolin, Ann DeBolt, Pete Del Zotto, Jennifer Hintz, Laurie Kurth, Bruce McCune, and Diane Taylor assisted in the field. Rober Pal provided helpful comments on the manuscript. This study was funded by the U.S. National Park Service.