Experimental warming differentially affects vegetative and reproductive phenology of tundra plants

Rapid climate warming is altering Arctic and alpine tundra ecosystem structure and function, including shifts in plant phenology. While the advancement of green up and flowering are well-documented, it remains unclear whether all phenophases, particularly those later in the season, will shift in unison or respond divergently to warming. Here, we present the largest synthesis to our knowledge of experimental warming effects on tundra plant phenology from the International Tundra Experiment. We examine the effect of warming on a suite of season-wide plant phenophases. Results challenge the expectation that all phenophases will advance in unison to warming. Instead, we find that experimental warming caused: (1) larger phenological shifts in reproductive versus vegetative phenophases and (2) advanced reproductive phenophases and green up but delayed leaf senescence which translated to a lengthening of the growing season by approximately 3%. Patterns were consistent across sites, plant species and over time. The advancement of reproductive seasons and lengthening of growing seasons may have significant consequences for trophic interactions and ecosystem function across the tundra.publishedVersio

The tundra phenology database: more than two decades of tundra phenology responses to climate change

Observations of changes in phenology have provided some of the strongest signals of the effects of climate change on terrestrial ecosystems. The International Tundra Experiment (ITEX), initiated in the early 1990s, established a common protocol to measure plant phenology in tundra study areas across the globe. Today, this valuable collection of phenology measurements depicts the responses of plants at the colder extremes of our planet to experimental and ambient changes in temperature over the past decades. The database contains 150 434 phenology observations of 278 plant species taken at 28 study areas for periods of 1\u201326 years. Here we describe the full data set to increase the visibility and use of these data in global analyses and to invite phenology data contributions from underrepresented tundra locations. Portions of this tundra phenology database have been used in three recent syntheses, some data sets are expanded, others are from entirely new study areas, and the entirety of these data are now available at the Polar Data Catalogue (https://doi.org/10.21963/13215)

Global maps of soil temperature

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2 m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km2 resolution for 0–5 and 5–15 cm soil depth. These maps were created by calculating the difference (i.e. offset) between in situ soil temperature measurements, based on time series from over 1200 1-km2 pixels (summarized from 8519 unique temperature sensors) across all the world\u27s major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10°C (mean = 3.0 ± 2.1°C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6 ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler (−0.7 ± 2.3°C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications

Global maps of soil temperature

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2 m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km² resolution for 0–5 and 5–15 cm soil depth. These maps were created by calculating the difference (i.e., offset) between in-situ soil temperature measurements, based on time series from over 1200 1-km² pixels (summarized from 8500 unique temperature sensors) across all the world’s major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10°C (mean = 3.0 ± 2.1°C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6 ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler (-0.7 ± 2.3°C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in-situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications

Global maps of soil temperature.

Research in global change ecology relies heavily on global climatic grids derived from estimates of air temperature in open areas at around 2 m above the ground. These climatic grids do not reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions occur and most terrestrial species reside. Here, we provide global maps of soil temperature and bioclimatic variables at a 1-km2 resolution for 0-5 and 5-15 cm soil depth. These maps were created by calculating the difference (i.e. offset) between in situ soil temperature measurements, based on time series from over 1200 1-km2 pixels (summarized from 8519 unique temperature sensors) across all the world's major terrestrial biomes, and coarse-grained air temperature estimates from ERA5-Land (an atmospheric reanalysis by the European Centre for Medium-Range Weather Forecasts). We show that mean annual soil temperature differs markedly from the corresponding gridded air temperature, by up to 10°C (mean = 3.0 ± 2.1°C), with substantial variation across biomes and seasons. Over the year, soils in cold and/or dry biomes are substantially warmer (+3.6 ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are on average slightly cooler (-0.7 ± 2.3°C). The observed substantial and biome-specific offsets emphasize that the projected impacts of climate and climate change on near-surface biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments. The global soil-related bioclimatic variables provided here are an important step forward for any application in ecology and related disciplines. Nevertheless, we highlight the need to fill remaining geographic gaps by collecting more in situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications

Overview of the nine calculated cumulative temperature variables (temperature sums) based on daily average temperatures.

<p>Melt-out dates were recorded for each subplot.</p

Estimates of treatment effects on plant sizes throughout the growing season, as well as peak size and number of thawing degree days needed to reach peak size.

<p>Effect values other than the intercept (here: Normal treatment) are deviations from the latter. Given are standard deviation (sd), t- and p-value and degrees of freedom (df) of the model. Modeled maximal plant sizes and TDD (thawing degree days, i.e. positive air temperature sums in °C) to reach the maximum size are based on the model shown in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0086281#pone-0086281-t003" target="_blank">Table 3</a>, and are shown as values for each treatment, rather than deviation from the Normal treatment. Normal = unmanipulated snow cover; Removed = snow removal on 1. May; Shallow = naturally early snowmelt; Medium = intermediately increased snow; Deep = maximally increased snow. Measures for plant species: average sum of leaf lengths (in mm) for <i>Alopecurus magellanicus</i>, <i>Dryas octopetala</i>, <i>Luzula arcuata</i>, <i>Salix polaris</i>; average sum of leaf areas (in mm²) for <i>Bistorta vivipara</i>; average annual shoot increment (in mm) for <i>Cassiope tetragona</i>; average plant sizes (in mm) for <i>Pedicularis hirsuta</i>, <i>Stellaria crassipes</i>.</p

Overview of the species-specific parameters per individual, growth form and habitat association.

<p>Overview of the species-specific parameters per individual, growth form and habitat association.</p

Estimates of treatment effects on melt-out dates.

<p>Melt-out dates in days of year (doy) and on TDD (thawing degree days, i.e. positive air temperature sums in °C). Effect values other than the intercept (here: Normal treatment) are deviations from the latter. Normal = unmanipulated snow cover; Removed = snow removal on 1. May; Shallow = naturally early snowmelt; Medium = intermediately increased snow; Deep = maximally increased snow. Given are standard deviation (sd), t- and p-values and degrees of freedom (df) of the model.</p