Global maps of soil temperature.

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-km2 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-km2 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.

Canopy photosynthesis of six major arable crops is enhanced under diffuse light due to canopy architecture.

Diffuse radiation generally increases photosynthetic rates if total radiation is kept constant. Different hypotheses have been proposed to explain this enhancement of photosynthesis, but conclusive results over a wide range of diffuse conditions or about the effect of canopy architecture are lacking. Here, we show the response of canopy photosynthesis to different fractions of diffuse light conditions for five major arable crops (pea, potato, wheat, barley, rapeseed) and cover crops characterized by different canopy architecture. We used 13 years of flux and microclimate measurements over a field with a typical 4 year crop rotation scheme in Switzerland. We investigated the effect of diffuse light on photosynthesis over a gradient of diffuse light fractions ranging from 100% diffuse (overcast sky) to 11% diffuse light (clear-sky conditions). Gross primary productivity (GPP) increased with diffuse fraction and thus was greater under diffuse than direct light conditions if the absolute photon flux density per unit surface area was kept constant. Mean leaf tilt angle (MTA) and canopy height were found to be the best predictors of the diffuse versus direct radiation effect on photosynthesis. Climatic factors, such as the drought index and growing degree days (GDD), had a significant influence on initial quantum yield under direct but not diffuse light conditions, which depended primarily on MTA. The maximum photosynthetic rate at 2,000 µmol m-2  s-1 photosynthetically active radiation under direct conditions strongly depended on GDD, MTA, leaf area index (LAI) and the interaction between MTA and LAI, while under diffuse conditions, this parameter depended mostly on MTA and only to a minor extent on canopy height and their interaction. The strongest photosynthesis enhancement under diffuse light was found for wheat, barley and rapeseed, whereas the lowest was for pea. Thus, we suggest that measuring canopy architecture and diffuse radiation will greatly improve GPP estimates of global cropping systems.

Joint control of terrestrial gross primary productivity by plant phenology and physiology.

Terrestrial gross primary productivity (GPP) varies greatly over time and space. A better understanding of this variability is necessary for more accurate predictions of the future climate-carbon cycle feedback. Recent studies have suggested that variability in GPP is driven by a broad range of biotic and abiotic factors operating mainly through changes in vegetation phenology and physiological processes. However, it is still unclear how plant phenology and physiology can be integrated to explain the spatiotemporal variability of terrestrial GPP. Based on analyses of eddy-covariance and satellite-derived data, we decomposed annual terrestrial GPP into the length of the CO2 uptake period (CUP) and the seasonal maximal capacity of CO2 uptake (GPPmax). The product of CUP and GPPmax explained >90% of the temporal GPP variability in most areas of North America during 2000-2010 and the spatial GPP variation among globally distributed eddy flux tower sites. It also explained GPP response to the European heatwave in 2003 (r(2) = 0.90) and GPP recovery after a fire disturbance in South Dakota (r(2) = 0.88). Additional analysis of the eddy-covariance flux data shows that the interbiome variation in annual GPP is better explained by that in GPPmax than CUP. These findings indicate that terrestrial GPP is jointly controlled by ecosystem-level plant phenology and photosynthetic capacity, and greater understanding of GPPmax and CUP responses to environmental and biological variations will, thus, improve predictions of GPP over time and space.

N2O emission factors for cattle urine: effect of patch characteristics and environmental drivers.

Urine patches from grazing cattle are hotspots of nitrous oxide (N2O) emissions. The default IPCC emission factor for urine patches (EFurine) is 0.77% for wet climates and 0.32% for dry climates. However, literature reports a considerable range of cattle urine EF values and urine characteristics used in experimental studies, revealing contrary results on the effects of urine patch characteristics and seasonal pattern. Therefore, we examined N2O emissions and corresponding EFurine values in relation to urine patch characteristics (urine N concentration, urine volume, patch area, urine composition) and environmental drivers (precipitation, water filled pore space, soil temperature). Ten artificial urine application experiments were performed from July 2020 to June 2022 on a pasture located in Eastern Switzerland. Urine N concentration, patch area, volume and urine N composition showed no significant effects on the EFurine value (p > 0.05). EFurine varied, however, strongly over time (0.17-2.05%). A large part of the variation could be predicted either by cumulative precipitation 20 days after urine application using a second order polynomial model (Adj. R2 = 0.60) or average WFPS 30 days after urine application using a linear model (Adj. R2 = 0.45). The derived precipitation model was used to simulate EFurine weekly over the last 20 years showing no significant differences between the seasons of a year. The resulting overall average EFurine was 0.67%. More field studies are needed across sites/regions differing in climate and soil properties to implement a country-specific EF3 for Switzerland and to improve the quantification of N2O emissions at the national scales.

Thermal optimality of net ecosystem exchange of carbon dioxide and underlying mechanisms.

• 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.

Large regional differences of soil water limitation effect on ozone induced yield loss for wheat and potato in Switzerland.

The accumulated stomatal ozone (O3) uptake over a threshold (Phytotoxic Ozone Dose POD6), calculated by an ozone deposition model, has been shown to be the most appropriate metric to quantify negative effects of O3 on food crops. In this study we used data of 13 sites in different regions of Switzerland with multiple years of O3 measurements to quantify the stomatal O3 uptake and the related yield loss of wheat and potato. Flux patterns for different years were calculated with the DO3SE model to disentangle the influence of contrasting seasonal environmental conditions. Regional and inter-annual differences in meteorological conditions led to considerable variations in soil water conditions and the POD6 values for wheat. Potato stomatal uptake was much less influenced by soil water and showed a more even distribution of POD6 values across sites and years. The estimated nationally and temporally average yield loss was 3.2 ± 1.2% for wheat and 2.4 ± 0.8% for potato, calculated based on an area weighting. It was found that soil water deficit, observed frequently in the western part of Switzerland, had a large attenuation effect on stomatal O3 uptake by wheat and on corresponding yield losses. This highlights the importance of including soil moisture limitation in O3 uptake modelling even in moist climatic regions. The comparison of modelled evapotranspiration with water flux measurements over a wheat field showed a reasonable agreement concerning the temporal pattern and the magnitude. But it is also concluded that the DO3SE soil moisture module will need further testing and adaptation to improve accuracy of the model in dryer conditions.

The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data.

The FLUXNET2015 dataset provides ecosystem-scale data on CO2, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Application and test of a simple tool for operational footprint evaluations.

We present a user-friendly tool for footprint calculations of flux measurements in the surface layer. The calculations are based on the analytical footprint model by Kormann, R. and Meixner, F.X. [2001. An analytical footprint model for Non-neutral Stratification. Boundary-Layer Meteorology 99, 207-224]. The footprint density function of a flux sensor is determined using readily available data from standard eddy covariance measurements. This footprint density function is integrated over defined surface areas given as quadrangular polygons representing e.g. agricultural fields. We illustrate the use and performance of the tool by applying it to CO2 flux measurements with three eddy covariance system at the Swiss CarboEurope grassland site. Two flux towers were positioned in the centre of two neighbouring fields, respectively, that showed a very different CO2 flux during the study period. The third tower was located near the border of the two fields and was frequently influenced by both fields to a similar degree. The calculated footprint fractions were used to simulate the latter flux from the other two systems. The measured and simulated fluxes showed a good agreement and thus support the reliability of the footprint calculation. The presented simple footprint tool can be used as a routine quality check for flux monitoring stations influenced by surface areas with varying vegetation covers and/or land-use.

High-resolution modelling of AOT40 and stomatal ozone uptake in wheat and grassland: a comparison between 2000 and the hot summer of 2003 in Switzerland.

The aim was to compare the ozone risk for agricultural crops in Switzerland during the hot and dry year 2003 with the more 'normal' situation in 2000. An improved version of the Ozone DEposition Model ODEM was used at a 2 x 2 km resolution. The distribution of the index AOT40 was compared with the accumulated stomatal ozone flux, AF(st). Averaged AOT40 at 2 m and at canopy height was much higher in 2003 than in 2000, but inter-annual differences in AF(st) for wheat and grasslands were small due to the limiting effect of low soil water contents in 2003. AOT40 suggested larger potential yield losses in wheat in 2003, while using AF(st) with a threshold of 6 nmol m(-2) s(-1) (AF(st)6) yielded similar estimates for both years. The data show that modelling of AF(st) can be used to differentiate ozone risks between regions and years at a national scale.

A modeling study on mitigation of N2O emissions and NO3 leaching at different agricultural sites across Europe using LandscapeDNDC.

The identification of site-specific agricultural management practices in order to maximize yield while minimizing environmental nitrogen losses remains in the center of environmental pollution research. Here, we used the biogeochemical model LandscapeDNDC to explore different agricultural practices with regard to their potential to reduce soil N2O emissions and NO3 leaching while maintaining yields. In a first step, the model was tested against observations of N2O emissions, NO3 leaching, soil micrometeorology as well as crop growth for eight European cropland and grassland sites. Across sites, LandscapeDNDC predicts very well mean N2O emissions (r(2)=0.99) and simulates the magnitude and general temporal dynamics of soil inorganic nitrogen pools. For the assessment of site-specific mitigation potentials of environmental nitrogen losses a Monte Carlo optimization technique considering different agricultural management options (i.e., timing of planting, harvest and fertilization, amount of applied fertilizer as well as residue management) was used. The identified optimized field management practices reduce N2O emissions and NO3 leaching from croplands on average by 21% and 31%, respectively. Likewise, average reductions of 55% for N2O emissions and 16% for NO3 leaching are estimated for grasslands. For mitigating environmental loss - while maintaining yield levels - it was most important to reduce fertilizer application rates by in average 10%. Our analyses indicate that yield scaled N2O emissions and NO3 leaching indicate possible improvements of nitrogen use efficiencies in European farming systems. Moreover, the applied optimization approach can be used also in a prognostic way to predict optimal timings and fertilization options (rates and splitting) upon accurate weather forecasts combined with the knowledge of modeled soil nutrient availability and plant nitrogen demand.