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.

Although variation in effect sizes and predicted values among studies of similar phenomena is inevitable, such variation far exceeds what might be produced by sampling error alone. One possible explanation for variation among results is differences among researchers in the decisions they make regarding statistical analyses. A growing array of studies has explored this analytical variability in different fields and has found substantial variability among results despite analysts having the same data and research question. Many of these studies have been in the social sciences, but one small “many analyst” study found similar variability in ecology. We expanded the scope of this prior work by implementing a large-scale empirical exploration of the variation in effect sizes and model predictions generated by the analytical decisions of different researchers in ecology and evolutionary biology. We used two unpublished datasets, one from evolutionary ecology (blue tit, Cyanistes caeruleus, to compare sibling number and nestling growth) and one from conservation ecology (Eucalyptus, to compare grass cover and tree seedling recruitment). The project leaders recruited 174 analyst teams, comprising 246 analysts, to investigate the answers to prespecified research questions. Analyses conducted by these teams yielded 141 usable effects (compatible with our meta-analyses and with all necessary information provided) for the blue tit dataset, and 85 usable effects for the Eucalyptus dataset. We found substantial heterogeneity among results for both datasets, although the patterns of variation differed between them. For the blue tit analyses, the average effect was convincingly negative, with less growth for nestlings living with more siblings, but there was near continuous variation in effect size from large negative effects to effects near zero, and even effects crossing the traditional threshold of statistical significance in the opposite direction. In contrast, the average relationship between grass cover and Eucalyptus seedling number was only slightly negative and not convincingly different from zero, and most effects ranged from weakly negative to weakly positive, with about a third of effects crossing the traditional threshold of significance in one direction or the other. However, there were also several striking outliers in the Eucalyptus dataset, with effects far from zero. For both datasets, we found substantial variation in the variable selection and random effects structures among analyses, as well as in the ratings of the analytical methods by peer reviewers, but we found no strong relationship between any of these and deviation from the meta-analytic mean. In other words, analyses with results that were far from the mean were no more or less likely to have dissimilar variable sets, use random effects in their models, or receive poor peer reviews than those analyses that found results that were close to the mean. The existence of substantial variability among analysis outcomes raises important questions about how ecologists and evolutionary biologists should interpret published results, and how they should conduct analyses in the future.

Climate change and other global change drivers threaten plant diversity in mountains worldwide. A widely documented response to such environmental modifications is for plant species to change their elevational ranges. Range shifts are often idiosyncratic and difficult to generalize, partly due to variation in sampling methods. There is thus a need for a standardized monitoring strategy that can be applied across mountain regions to assess distribution changes and community turnover of native and non‐native plant species over space and time. Here, we present a conceptually intuitive and standardized protocol developed by the Mountain Invasion Research Network (MIREN) to systematically quantify global patterns of native and non‐native species distributions along elevation gradients and shifts arising from interactive effects of climate change and human disturbance. Usually repeated every five years, surveys consist of 20 sample sites located at equal elevation increments along three replicate roads per sampling region. At each site, three plots extend from the side of a mountain road into surrounding natural vegetation. The protocol has been successfully used in 18 regions worldwide from 2007 to present. Analyses of one point in time already generated some salient results, and revealed region‐specific elevational patterns of native plant species richness, but a globally consistent elevational decline in non‐native species richness. Non‐native plants were also more abundant directly adjacent to road edges, suggesting that disturbed roadsides serve as a vector for invasions into mountains. From the upcoming analyses of time series, even more exciting results can be expected, especially about range shifts. Implementing the protocol in more mountain regions globally would help to generate a more complete picture of how global change alters species distributions. This would inform conservation policy in mountain ecosystems, where some conservation policies remain poorly implemented.

Pinus contorta is considered one of the most invasive tree species worldwide, generating significant impacts on biodiversity and ecosystems. In several Patagonian ecosystems in southern Chile, it has escaped from plantations established mainly in the 1970s, and is now invading both forests and treeless environments. In this study, we evaluated the impact of the invasion of P. contorta on microenvironmental conditions in Araucaria araucana forest and Patagonian steppe ecosystems, and assessed how these changes related to the richness and abundance of native and non-native plant species. In each ecosystem, 24 plots of 100 m² were established along a gradient of P. contorta biomass, where 18 environmental variables and the composition of native and non-native vegetation were measured at a local scale. Our results indicated that increased pine biomass was associated with differences in microclimatic conditions (soil and air temperature, photosynthetically active radiation (PAR), and soil moisture) and soil properties (potassium, nitrate, pH, and litter accumulation). These changes were ecosystem dependent, however, as well as associated with the level of invasion. Finally, the reduction in the richness and abundance of native plants was associated with the changes in soil properties (accumulation of leaf litter, pH, and organic matter) as well as in the microclimate (minimum air temperature, PAR) generated by the invasion of P. contorta. Overall, our results confirm that the invasion of P. contorta impacts microenvironmental conditions (i.e., canopy cover, litter accumulation, minimum air temperature, and maximum soil temperature) and reduces native plant diversity. For future restoration plans, more emphasis should be given to how environmental changes can influence the recovery of invaded ecosystems even after the removal of the living pine biomass (i.e., legacy of the invasion).