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.

The increase in aridity, mainly by decreases in precipitation but also by higher temperatures, is likely the main threat to the diversity and survival of Mediterranean forests. Changes in land use, including the abandonment of extensive crop activities, mainly in mountains and remote areas, and the increases in human settlements and demand for more resources with the resulting fragmentation of the landscape, hinder the establishment of appropriate management tools to protect Mediterranean forests and their provision of services and biodiversity. Experiments and observations indicate that if changes in climate, land use and other components of global change, such as pollution and overexploitation of resources, continue, the resilience of many forests will likely be exceeded, altering their structure and function and changing, mostly decreasing, their capacity to continue to provide their current services. A consistent assessment of the impacts of the changes, however, remains elusive due to the difficulty of obtaining simultaneous and complete data for all scales of the impacts in the same forests, areas and regions. We review the impacts of climate change and other components of global change and their interactions on the terrestrial forests of Mediterranean regions, with special attention to their impacts on ecosystem services. Management tools for counteracting the negative effects of global change on Mediterranean ecosystem- services are finally discussed.

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.

Widespread population declines have been reported for diverse Mediterranean butterflies over the last three decades, and have been significantly associated with increased global change impacts. The specific landscape and climatic drivers of these declines remain uncertain for most declining species. Here, we analyse whether plastic phenotypic traits of a model butterfly species (Pieris napi) perform as reliable biomarkers of vulnerability to extreme temperature impacts in natural populations, showing contrasting trends in thermally exposed and thermally buffered populations. We also examine whether improved descriptions of thermal exposure of insect populations can be achieved by combining multiple information sources (i.e., integrating measurements of habitat thermal buffering, habitat thermal amplification, host plant transpiration, and experimental assessments of thermal death time (TDT), thermal avoidance behaviour (TAB) and thermally induced trait plasticity). These integrative analyses are conducted in two demographically declining and two non‐declining populations of P. napi. The results show that plastic phenotypic traits (butterfly body mass and wing size) are reliable biomarkers of population vulnerability to extreme thermal conditions. Butterfly wing size is strongly reduced only in thermally exposed populations during summer drought periods. Laboratory rearing of these populations documented reduced wing size due to significant negative effects of increased temperatures affecting larval growth. We conclude that these thermal biomarkers are indicative of the population vulnerability to increasing global warming impacts, showing contrasting trends in thermally exposed and buffered populations. Thermal effects in host plant microsites significantly differ between populations, with stressful thermal conditions only effectively ameliorated in mid‐elevation populations. In lowland populations, we observe a sixfold reduction in vegetation thermal buffering effects, and larval growth occurs in these populations at significantly higher temperatures. Lowland populations show reduced host plant quality (C/N ratio), reduced leaf transpiration rates and complete above‐ground plant senescence during the peak of summer drought. Amplified host plant temperatures are observed in open microsites, reaching thermal thresholds that can affect larval survival. Overall, our results suggest that butterfly population vulnerability to long‐term drought periods is associated with multiple co‐occurring and interrelated ecological factors, including limited vegetation thermal buffering effects at lowland sites, significant drought impacts on host plant transpiration and amplified leaf surface temperature, as well as reduced leaf quality linked to the seasonal advance of plant phenology. Our results also identify multiannual summer droughts affecting larval growing periods as a key driver of the recently reported butterfly population declines in the Mediterranean biome.