Report

Changes in Climatic Water Balance Drive Downhill Shifts in Plant Species’ Optimum Elevations

See allHide authors and affiliations

Science  21 Jan 2011:
Vol. 331, Issue 6015, pp. 324-327
DOI: 10.1126/science.1199040

Abstract

Uphill shifts of species’ distributions in response to historical warming are well documented, which leads to widespread expectations of continued uphill shifts under future warming. Conversely, downhill shifts are often considered anomalous and unrelated to climate change. By comparing the altitudinal distributions of 64 plant species between the 1930s and the present day within California, we show that climate changes have resulted in a significant downward shift in species’ optimum elevations. This downhill shift is counter to what would be expected given 20th-century warming but is readily explained by species’ niche tracking of regional changes in climatic water balance rather than temperature. Similar downhill shifts can be expected to occur where future climate change scenarios project increases in water availability that outpace evaporative demand.

Climatic warming during the 20th century has led to a variety of responses from biota (1), including changes in phenology for plants (2, 3) and animals (4, 5), and shifts in the geographic distributions of species poleward in latitude (6, 7) or uphill in elevation (8, 9). Because climate-change scenarios project warming during the 21st century (10), changes in the distributions of species are predicted to continue, with shifts toward higher elevations and latitudes projected for species that are able to track changes in temperature (11).

The assumption that temperature is the principal factor defining species’ distributions ignores the fact that many species, including plants, are constrained by energy and water availability (12, 13). Consequently, considering changes in temperature alone may be inadequate for understanding distributional shifts of plant species. Consistent with this, there is ample evidence of stasis and even downhill shifts in species’ distributions despite climatic warming (1416). These seemingly inconsistent responses are attributed to changes in competitive interactions at range margins (17), human-induced landscape modification (14), lags in the effect of climatic warming (15), or idiosyncratic responses of individual species (18). As such, there has been limited progress in understanding the drivers of shifts in species’ distributions and their associated mechanisms.

We sought to determine how the altitudinal distributions of vascular plant species in California, USA, have changed during the 20th century and whether changes, if any, could be attributed to changing climatic parameters. We defined our study area based on the availability of long-term climate station data and spatial coverage of vegetation survey data in both historical (1930–1935) and modern (2000–2005) time periods (19). This region includes approximately one-half of the state (177,000 km2) and encompasses most of the major mountain ranges north of 35° latitude (Fig. 1). Mean annual temperatures in California increased by ~0.6°C during this time period, with warming occurring across the entire state (19). Changes in precipitation exhibited regional-scale variability, with increases across the northern portion of the state (19) that have resulted in a net decrease in the climatic water deficit (difference between potential evapotranspiration and precipitation) relative to the early 20th century (Fig. 2).

Fig. 1

Portion of California used in study (~177,000 km2).

Fig. 2

Change in annual precipitation (A), potential evapotranspiration (B), and climatic water deficit (C) between historical (1920–1949) and modern (1976–2005) time periods at long-term meteorological stations in California, USA (study area outlined in black). Values expressed in mm and precipitation changes inverted for consistency with water deficit changes. Stations that exhibited statistically significant differences are shown by circles. Anomalies (red, positive; blue, negative) for mean annual temperature (D), potential evapotranspiration (E), annual precipitation (F), and climatic water deficit (G) were averaged across 33 climate stations from our study area. Solid curve represents a running mean of 11 years. Linear trend is indicated in each plot and is significant at the 90th percentile.

We quantified species’ responses by analyzing shifts in species’ optimum elevations, the elevation at which there is a maximum probability of occurrence within the species’ altitudinal distribution (16). Using modern climate records and records from the early 20th century, we tested whether changes in optimum elevation were consistent with patterns of observed changes in climate (19). We quantified climatic water balance using climatic water deficit, a measure of the evaporative demand of a reference crop that is not met by precipitation and thus estimates drought at a site independent of that site’s vegetation (12). Climate data from weather stations in our study area indicated a decrease in climatic water deficit of 100.0 mm [63.3 to 136.7 mm, 95% confidence interval (CI)] between the early and latter years of the 20th century (Fig. 2) and are consistent with a decreasing linear trend of 19.4 mm per decade from 1920 to 2005 (19) (Fig. 2 and fig. S1).

We used two temporally independent data sets of survey plots collected from the major mountain ranges of the state. The first data set comprised 13,746 survey plots sampled in the 1930s (1921); the second was compiled from a variety of sources and comprised ~33,000 plots sampled from 2000 to 2005 (19). We conducted a subsampling procedure to remove potential sampling bias between the data sets (19), which left a balanced sample of 8747 plots in each period. We used interpolated climate surfaces representing 30-year monthly normals to extract climatic data from historical (1905–1935) and modern (1975–2005) plot locations (19, 22). We used logistic regression to estimate species’ optimum elevations in each period, focusing on species with at least 50 occurrences that exhibited unimodal bell-shaped response functions in each period (16). Altitudinal shifts were measured as the difference in optimum elevation between periods.

We found significant downhill shifts in optimum elevations (mean difference = –88.2 m, t = –2.49, df = 63, P = 0.016), with a higher proportion of species shifting their distributions downhill [proportion (p) = 0.72, 95% CI = 0.59 to 0.82] than uphill (p = 0.28, 95% CI = 0.18 to 0.41) (Fig. 3). We found no changes in optimum climatic water deficit (mean difference = –11.1 mm, t = –0.20, df = 45, P = 0.843), with similar proportions shifting their distributions toward drier (p = 0.41, 95% CI = 0.27 to 0.57) and wetter conditions (p = 0.59, 95% CI = 0.43 to 0.73) (Fig. 3). There was moderate evidence of increases in optimum temperatures experienced by species (mean difference = 0.37°C, t = 1.87, df = 65, P = 0.067), with a higher proportion of species shifting toward warmer conditions (p = 0.67, 95% CI = 0.54 to 0.78) than cooler conditions (p = 0.33, 95% CI = 0.22 to 0.46) (Fig. 3). Species traits were unrelated to shifts in optimums (Kruskal-Wallis tests, P > 0.05 for all comparisons) (19).

Fig. 3

Scatterplots of optimum elevation (A), climatic water deficit (B), and mean annual temperature (C) between historical (circa 1935) and modern (circa 2005) time periods. Each point represents a single species. Species below the diagonal represent species with decreasing optimums. Species with nonoverlapping 95% CIs around optimums between time periods are displayed as filled triangles (▲); species with overlapping 95% CIs are displayed as open circles (○). Insets represent number of species with increasing or decreasing optimums.

Plant species in our study area appear to be tracking their climatic niche by shifting their altitudinal distributions downhill in response to decreased climatic water deficit. These shifts represent niche tracking of changes in water availability rather than changes in temperature. Assuming constant temperature and precipitation lapse rates for our study area, we estimated that the observed historical changes in temperature and precipitation would result in a mean downward shift in optimum elevation of 85 m (95% CI = 64 to 106 m) (19), an interval that covers the observed change in mean optimum elevation (88.2 m downhill). The consistency of our results across a wide array of taxa (46 species from 17 families) (table S1) and broad geographic scale suggests that these are the effects of large-scale climatic drivers of species’ distributions. One line of evidence to support this conclusion is that species with higher optimum elevations exhibited larger downhill shifts (linear model, β1 = –0.162, F = 11.8, df = 62, P = 0.001) (Fig. 4). Regional changes in precipitation and temperature differentially alter the local climatic water balance across complex topography owing to the nonlinear relationship between temperature and vapor pressure deficit. Cooler and wetter sites at higher elevations would be expected to exhibit greater sensitivity to changes in precipitation, whereas warmer and drier sites at lower elevations would be expected to exhibit greater sensitivity to changes in temperature. Given the observed study area–wide increase in precipitation, we would thus expect downhill shifts to be more likely at wetter sites (higher elevations) (Fig. 4). Our results are also consistent with observational studies within the region showing an increase in the density of younger cohorts of tree species at lower elevations during the past century (18, 23). Further, we found no evidence of changes in the width of species’ response curves between periods, which suggests that the observed shifts in elevation are occurring across their entire altitudinal distribution rather than at range margins (fig. S4). This, coupled with evidence of downhill shifts across the entire elevation gradient (Fig. 4), supports a widespread driver such as climate.

Fig. 4

Scatterplot of shift in optimum elevation (m) versus historical (circa 1935) altitudinal position (m) for plant species in California. Solid line represents linear regression model with 95% CI bands (as dashed lines).

The assumption that plant species will move uphill to track increases in temperature does not account for coupled environmental constraints on their growth and persistence. We show that plant species in our study area are experiencing an increase in their optimum temperature (0.36°C) due to both climatic warming and downhill shifts. For these organisms, the effects of warming on moisture balance are being offset by enhanced water availability. This combination of factors has likely increased net primary productivity in the region (24). Our findings highlight the importance of characterizing climatic water balance for understanding plant responses to climate change. Regional increases in climatic water balance have been identified elsewhere in the northern hemisphere (25), which suggests that downhill shifts by plant species may not be isolated to our study region. Most locations north of 45° latitude have experienced increased precipitation over the past century (10), and at continental scales, global climate models generally predict those locations to become wetter over the next century (10). If increases in precipitation outpace increases in evaporative demand in these regions, we can expect downhill shifts by species to occur.

Supporting Online Material

www.sciencemag.org/cgi/content/full/331/6015/324/DC1

Materials and Methods

Figs. S1 to S4

Tables S1 to S3

References

References and Notes

  1. Materials and methods are available as supporting material on Science Online.
  2. Supported by NSF (BCS-0819430 to S.Z.D. and BCS-0819493 to J.A.G.). J.T.A. was supported by NSF’s Experimental Program to Stimulate Competitive Research (EPSCoR) (EPS-0814387). A.R.M. was supported by U.S. Department of Agriculture Forest Service’s (USFS’s) Rocky Mountain Research Station (JV11221635-201). We thank C. Daly for providing interpolated climate data; M. Kelly for Vegetation Type Mapping project data; C. Ramirez (Region 5 USFS Remote Sensing Lab) for the USFS Forest Inventory and Analysis data (plot locations were not provided to the authors); J. Thorne for compiling modern data sources; and S. Running for helpful comments on the manuscript.
View Abstract

Navigate This Article