Many wetland ecosystems such as peatlands and wet tundra hold large amounts of organic carbon (C) in their soils, and are thus important in the terrestrial C cycle. We have synthesized data on the carbon dioxide (CO 2 ) exchange obtained from eddy covariance measurements from 12 wetland sites, covering 1-7 years at each site, across Europe and North America, ranging from ombrotrophic and minerotrophic peatlands to wet tundra ecosystems, spanning temperate to arctic climate zones. The average summertime net ecosystem exchange of CO 2 (NEE) was highly variable between sites. However, all sites with complete annual datasets, seven in total, acted as annual net sinks for atmospheric CO 2 . To evaluate the influence of gross primary production (GPP) and ecosystem respiration (R eco ) on NEE, we first removed the artificial correlation emanating from the method of partitioning NEE into GPP and R eco . After this correction neither R eco (P 5 0.162) nor GPP (P 5 0.110) correlated significantly with NEE on an annual basis. Spatial variation in annual and summertime R eco was associated with growing season period, air temperature, growing degree days, normalized difference vegetation index and vapour pressure deficit. GPP showed weaker correlations with environmental variables as compared with R eco , the exception being leaf area index (LAI), which correlated with both GPP and NEE, but not with R eco . Length of growing season period was found to be the most important variable describing the spatial variation in summertime GPP and R eco ; global warming will thus cause these components to increase. Annual GPP and NEE correlated significantly with LAI and pH, thus, in order to predict wetland C exchange, differences in ecosystem structure such as leaf area and biomass as well as nutritional status must be taken into account.

Ecosystems are increasingly prone to climate extremes, such as drought, with long-lasting effects on both plant and soil communities and, subsequently, on carbon (C) cycling. However, recent studies underlined the strong variability in ecosystem's response to droughts, raising the issue of nonlinear responses in plant and soil communities. The conundrum is what causes ecosystems to shift in response to drought. Here, we investigated the response of plant and soil fungi to drought of different intensities using a water table gradient in peatlands-a major C sink ecosystem. Using moving window structural equation models, we show that substantial changes in ecosystem respiration, plant and soil fungal communities occurred when the water level fell below a tipping point of -24 cm. As a corollary, ecosystem respiration was the greatest when graminoids and saprotrophic fungi became prevalent as a response to the extreme drought. Graminoids indirectly influenced fungal functional composition and soil enzyme activities through their direct effect on dissolved organic matter quality, while saprotrophic fungi directly influenced soil enzyme activities. In turn, increasing enzyme activities promoted ecosystem respiration. We show that functional transitions in ecosystem respiration critically depend on the degree of response of graminoids and saprotrophic fungi to drought. Our results represent a major advance in understanding the nonlinear nature of ecosystem properties to drought and pave the way towards a truly mechanistic understanding of the effects of drought on ecosystem processes.

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 thus fail to reflect conditions below vegetation canopies and near the ground surface, where critical ecosystem functions are controlled 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 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 of 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 2 m 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 warmer than gridded air temperature), whereas soils in warm and humid environments are on average slightly cooler (0.7 ± 2.3°C cooler). The observed substantial and biome-specific offsets underpin that the projected impacts of climate and climate change on 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 global gaps by collecting more in-situ measurements of microclimate conditions to further enhance the spatiotemporal resolution of global soil temperature products for ecological applications.

Abstract. Natural wetlands constitute the largest and most uncertain source of methane (CH4) to the atmosphere and a large fraction of them are found in the northern latitudes. These emissions are typically estimated using process (“bottom-up”) or inversion (“top-down”) models. However, estimates from these two types of models are not independent of each other since the top-down estimates usually rely on the a priori estimation of these emissions obtained with process models. Hence, independent spatially explicit validation data are needed. Here we utilize a random forest (RF) machine-learning technique to upscale CH4 eddy covariance flux measurements from 25 sites to estimate CH4 wetland emissions from the northern latitudes (north of 45∘ N). Eddy covariance data from 2005 to 2016 are used for model development. The model is then used to predict emissions during 2013 and 2014. The predictive performance of the RF model is evaluated using a leave-one-site-out cross-validation scheme. The performance (Nash–Sutcliffe model efficiency =0.47) is comparable to previous studies upscaling net ecosystem exchange of carbon dioxide and studies comparing process model output against site-level CH4 emission data. The global distribution of wetlands is one major source of uncertainty for upscaling CH4. Thus, three wetland distribution maps are utilized in the upscaling. Depending on the wetland distribution map, the annual emissions for the northern wetlands yield 32 (22.3–41.2, 95 % confidence interval calculated from a RF model ensemble), 31 (21.4–39.9) or 38 (25.9–49.5) Tg(CH4) yr−1. To further evaluate the uncertainties of the upscaled CH4 flux data products we also compared them against output from two process models (LPX-Bern and WetCHARTs), and methodological issues related to CH4 flux upscaling are discussed. The monthly upscaled CH4 flux data products are available at https://doi.org/10.5281/zenodo.2560163 (Peltola et al., 2019).

Research in environmental science relies heavily on global climatic grids derived from estimates of air temperature at around 2 meter above ground1-3. These climatic grids however fail to reflect conditions near and below the soil surface, where critical ecosystem functions such as soil carbon storage are controlled and most biodiversity resides4-8. By using soil temperature time series from over 8500 locations across all of the world’s terrestrial biomes4, we derived global maps of soil temperature-related variables at 1 km resolution for the 0–5 and 5–15 cm depth horizons. Based on these maps, we show that mean annual soil temperature differs markedly from the corresponding 2 m gridded air temperature, by up to 10°C, with substantial variation across biomes and seasons. Soils in cold and/or dry biomes are annually substantially warmer (3.6°C ± 2.3°C) than gridded air temperature, whereas soils in warm and humid environments are slightly cooler (0.7 ± 2.3°C). As a result, annual soil temperature varies less (by 17%) across the globe than air temperature. The effect of macroclimatic conditions on the difference between soil and air temperature highlights the importance of considering that macroclimate warming may not result in the same level of soil temperature warming. Similarly, changes in precipitation could alter the relationship between soil and air temperature, with implications for soil-atmosphere feedbacks9. Our results underpin that the impacts of climate and climate change on biodiversity and ecosystem functioning are inaccurately assessed when air rather than soil temperature is used, especially in cold environments.