A new method for extracting seasonality information from time-series of satellite sensor data is presented. The method is based on nonlinear least squares fits of asymmetric Gaussian model functions to the time-series. The smooth model functions are then used for defining key seasonality parameters, such as the number of growing seasons, the beginning and end of the seasons, and the rates of growth and decline. The method is implemented in a computer program TIMESAT and tested on Advanced Very High Resolution Radiometer (AVHRR) normalized difference vegetation index (NDVI) data over Africa. Ancillary cloud data [clouds from AVHRR (CLAVR)] are used as estimates of the uncertainty levels of the data values. Being general in nature, the proposed method can be applied also to new types of satellite-derived time-series data.

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.