Many wetland ecosystems such as peatlands and wet tundra hold large amounts of organic carbon (C) in their soils, and are thus important in the terrestrial C cycle. We have synthesized data on the carbon dioxide (CO 2 ) exchange obtained from eddy covariance measurements from 12 wetland sites, covering 1-7 years at each site, across Europe and North America, ranging from ombrotrophic and minerotrophic peatlands to wet tundra ecosystems, spanning temperate to arctic climate zones. The average summertime net ecosystem exchange of CO 2 (NEE) was highly variable between sites. However, all sites with complete annual datasets, seven in total, acted as annual net sinks for atmospheric CO 2 . To evaluate the influence of gross primary production (GPP) and ecosystem respiration (R eco ) on NEE, we first removed the artificial correlation emanating from the method of partitioning NEE into GPP and R eco . After this correction neither R eco (P 5 0.162) nor GPP (P 5 0.110) correlated significantly with NEE on an annual basis. Spatial variation in annual and summertime R eco was associated with growing season period, air temperature, growing degree days, normalized difference vegetation index and vapour pressure deficit. GPP showed weaker correlations with environmental variables as compared with R eco , the exception being leaf area index (LAI), which correlated with both GPP and NEE, but not with R eco . Length of growing season period was found to be the most important variable describing the spatial variation in summertime GPP and R eco ; global warming will thus cause these components to increase. Annual GPP and NEE correlated significantly with LAI and pH, thus, in order to predict wetland C exchange, differences in ecosystem structure such as leaf area and biomass as well as nutritional status must be taken into account.

Light use efficiency (LUE) models are widely used to simulate gross primary production (GPP).However, the treatment of the plant canopy as a big leaf by these models can introduce large uncertainties in simulated GPP. Recently, a two-leaf light use efficiency (TL-LUE) model was developed to simulate GPP separately for sunlit and shaded leaves and has been shown to outperform the big-leaf MOD17 model at six FLUX sites in China. In this study we investigated the performance of the TL-LUE model for a wider range of biomes. For this we optimized the parameters and tested the TL-LUE model using data from 98 FLUXNET sites which are distributed across the globe. The results showed that the TL-LUE model performed in general better than the MOD17 model in simulating 8 day GPP. Optimized maximum light use efficiency of shaded leaves (ε msh ) was 2.63 to 4.59 times that of sunlit leaves (ε msu ). Generally, the relationships of ε msh and ε msu with ε max ZHOU ET AL. were well described by linear equations, indicating the existence of general patterns across biomes. GPP simulated by the TL-LUE model was much less sensitive to biases in the photosynthetically active radiation (PAR) input than the MOD17 model. The results of this study suggest that the proposed TL-LUE model has the potential for simulating regional and global GPP of terrestrial ecosystems, and it is more robust with regard to usual biases in input data than existing approaches which neglect the bimodal within-canopy distribution of PAR.

Summary• It is well established that individual organisms can acclimate and adapt to temperature to optimize their functioning. However, thermal optimization of ecosystems, as an assemblage of organisms, has not been examined at broad spatial and temporal scales.• Here, we compiled data from 169 globally distributed sites of eddy covariance and quantified the temperature response functions of net ecosystem exchange (NEE), an ecosystem-level property, to determine whether NEE shows thermal optimality and to explore the underlying mechanisms.• We found that the temperature response of NEE followed a peak curve, with the optimum temperature (corresponding to the maximum magnitude of NEE) being positively correlated with annual mean temperature over years and across sites. Shifts of the optimum temperature of NEE were mostly a result of temperature acclimation of gross primary productivity (upward shift of optimum temperature) rather than changes in the temperature sensitivity of ecosystem respiration.• Ecosystem-level thermal optimality is a newly revealed ecosystem property, presumably reflecting associated evolutionary adaptation of organisms within ecosystems, and has the potential to significantly regulate ecosystem-climate change feedbacks. The thermal optimality of NEE has implications for understanding fundamental properties of ecosystems in changing environments and benchmarking global models.

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.