Recent studies from mountainous areas of small spatial extent (<2500 km(2) ) suggest that fine-grained thermal variability over tens or hundreds of metres exceeds much of the climate warming expected for the coming decades. Such variability in temperature provides buffering to mitigate climate-change impacts. Is this local spatial buffering restricted to topographically complex terrains? To answer this, we here study fine-grained thermal variability across a 2500-km wide latitudinal gradient in Northern Europe encompassing a large array of topographic complexities. We first combined plant community data, Ellenberg temperature indicator values, locally measured temperatures (LmT) and globally interpolated temperatures (GiT) in a modelling framework to infer biologically relevant temperature conditions from plant assemblages within <1000-m(2) units (community-inferred temperatures: CiT). We then assessed: (1) CiT range (thermal variability) within 1-km(2) units; (2) the relationship between CiT range and topographically and geographically derived predictors at 1-km resolution; and (3) whether spatial turnover in CiT is greater than spatial turnover in GiT within 100-km(2) units. Ellenberg temperature indicator values in combination with plant assemblages explained 46-72% of variation in LmT and 92-96% of variation in GiT during the growing season (June, July, August). Growing-season CiT range within 1-km(2) units peaked at 60-65°N and increased with terrain roughness, averaging 1.97 °C (SD = 0.84 °C) and 2.68 °C (SD = 1.26 °C) within the flattest and roughest units respectively. Complex interactions between topography-related variables and latitude explained 35% of variation in growing-season CiT range when accounting for sampling effort and residual spatial autocorrelation. Spatial turnover in growing-season CiT within 100-km(2) units was, on average, 1.8 times greater (0.32 °C km(-1) ) than spatial turnover in growing-season GiT (0.18 °C km(-1) ). We conclude that thermal variability within 1-km(2) units strongly increases local spatial buffering of future climate warming across Northern Europe, even in the flattest terrains.

Current analyses and predictions of spatially explicit patterns and processes in ecology most often rely on climate data interpolated from standardized weather stations. This interpolated climate data represents long-term average thermal conditions at coarse spatial resolutions only. Hence, many climate-forcing factors that operate at fine spatiotemporal resolutions are overlooked. This is particularly important in relation to effects of observation height (e.g. vegetation, snow and soil characteristics) and in habitats varying in their exposure to radiation, moisture and wind (e.g. topography, radiative forcing or cold-air pooling). Since organisms living close to the ground relate more strongly to these microclimatic conditions than to free-air temperatures, microclimatic ground and near-surface data are needed to provide realistic forecasts of the fate of such organisms under anthropogenic climate change, as well as of the functioning of the ecosystems they live in. To fill this critical gap, we highlight a call for temperature time series submissions to SoilTemp, a geospatial database initiative compiling soil and near-surface temperature data from all over the world. Currently, this database contains time series from 7,538 temperature sensors from 51 countries

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.