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 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.

Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO₂ fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (V_cmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r² = 0.76; Nash–Sutcliffe modeling efficiency, MEF = 0.76) and ecosystem respiration (ER, r² = 0.78, MEF = 0.75), with lesser accuracy for latent heat fluxes (LE, r² = 0.42, MEF = 0.14) and and net ecosystem CO₂ exchange (NEE, r² = 0.38, MEF = 0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r² values (0.57–0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r² values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r² < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized V_cmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average V_cmax value.

One of the negative consequences of climate change is the also increase in the severity, frequency and length of droughts appearing in Europe. The effects of meteorological drought are often substantial, not only for the natural environment but also for humans. Hence, the main purpose of this research was to determine the trends in the severity and occurrence of droughts in Europe during the period 1951–2015 using the standardized precipitation index and the standardized precipitation and evapotranspiration index (SPEI). For six European sites located in mid latitudes, the number of dry months was determined and the trend of their occurrence was examined. Moreover, for the summer months in which the indicators fell below 0, the trend related to the severity of the drought was determined for each site. Despite the absence of a statistically significant trend of an increase in the occurrence of dry months in general, an increase in the severity of droughts occurring in summer was observed for the SPEI over a 6-month scale for all the investigated stations. The highest number of dry months since the 1970s appeared during the last 5 years of analysis.

In the global methane budget, the largest natural source is attributed to wetlands, which encompass all ecosystems composed of waterlogged or inundated ground, capable of methane production. Among them, northern peatlands that store large amounts of soil organic carbon have been functioning, since the end of the last glaciation period, as long-term sources of methane (CH₄) and are one of the most significant methane sources among wetlands. To reduce uncertainty of quantifying methane flux in the global methane budget, it is of significance to understand the underlying processes for methane production and fluxes in northern peatlands. A methane model that features methane production and transport by plants, ebullition process and diffusion in soil, oxidation to CO₂, and CH₄ fluxes to the atmosphere has been embedded in the ORCHIDEE-PEAT land surface model that includes an explicit representation of northern peatlands. ORCHIDEE-PCH₄ was calibrated and evaluated on 14 peatland sites distributed on both the Eurasian and American continents in the northern boreal and temperate regions. Data assimilation approaches were employed to optimized parameters at each site and at all sites simultaneously. Results show that methanogenesis is sensitive to temperature and substrate availability over the top 75 cm of soil depth. Methane emissions estimated using single site optimization (SSO) of model parameters are underestimated by 9 g CH₄ m⁻² yr⁻¹ on average (i.e., 50 % higher than the site average of yearly methane emissions). While using the multi-site optimization (MSO), methane emissions are overestimated by 5 g CH₄ m⁻² yr⁻¹ on average across all investigated sites (i.e., 37 % lower than the site average of yearly methane emissions).

The effects of atmospheric nitrogen deposition (Ndep) on carbon (C) sequestration in forests have often been assessed by relating differences in productivity to spatial variations of Ndep across a large geographic domain. These correlations generally suffer from covariation of other confounding variables related to climate and other growth-limiting factors, as well as large uncertainties in total (dry + wet) reactive nitrogen (Nr) deposition. We propose a methodology for untangling the effects of Ndep from those of meteorological variables, soil water retention capacity and stand age, using a mechanistic forest growth model in combination with eddy covariance CO2 exchange fluxes from a Europe-wide network of 22 forest flux towers. Total Nr deposition rates were estimated from local measurements as far as possible. The forest data were compared with data from natural or semi-natural, non-woody vegetation sites. The response of forest net ecosystem productivity to nitrogen deposition (dNEP ∕ dNdep) was estimated after accounting for the effects on gross primary productivity (GPP) of the co-correlates by means of a meta-modelling standardization procedure, which resulted in a reduction by a factor of about 2 of the uncorrected, apparent dGPP ∕ dNdep value. This model-enhanced analysis of the C and Ndep flux observations at the scale of the European network suggests a mean overall dNEP ∕ dNdep response of forest lifetime C sequestration to Ndep of the order of 40–50 g C per g N, which is slightly larger but not significantly different from the range of estimates published in the most recent reviews. Importantly, patterns of gross primary and net ecosystem productivity versus Ndep were non-linear, with no further growth responses at high Ndep levels (Ndep > 2.5–3 g N m−2 yr−1) but accompanied by increasingly large ecosystem N losses by leaching and gaseous emissions. The reduced increase in productivity per unit N deposited at high Ndep levels implies that the forecast increased Nr emissions and increased Ndep levels in large areas of Asia may not positively impact the continent's forest CO2 sink. The large level of unexplained variability in observed carbon sequestration efficiency (CSE) across sites further adds to the uncertainty in the dC∕dN response.