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.

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.

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.

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.

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.

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.

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.

How is water-use efficiency of terrestrial ecosystems distributed and changing on Earth?

A better understanding of ecosystem water-use efficiency (WUE) will help us improve ecosystem management for mitigation as well as adaption to global hydrological change. Here, long-term flux tower observations of productivity and evapotranspiration allow us to detect a consistent latitudinal trend in WUE, rising from the subtropics to the northern high-latitudes. The trend peaks at approximately 51°N, and then declines toward higher latitudes. These ground-based observations are consistent with global-scale estimates of WUE. Global analysis of WUE reveals existence of strong regional variations that correspond to global climate patterns. The latitudinal trends of global WUE for Earth's major plant functional types reveal two peaks in the Northern Hemisphere not detected by ground-based measurements. One peak is located at 20° ~ 30°N and the other extends a little farther north than 51°N. Finally, long-term spatiotemporal trend analysis using satellite-based remote sensing data reveals that land-cover and land-use change in recent years has led to a decline in global WUE. Our study provides a new framework for global research on the interactions between carbon and water cycles as well as responses to natural and human impacts.

High tolerance of subalpine grassland to long-term ozone exposure is independent of N input and climatic drivers.

In a seven-year study, we tested effects of increased N and O3 deposition and climatic conditions on biomass of subalpine grassland. Ozone risk was assessed as exposure (AOT40) and as stomatal flux (POD0,1). We hypothesized that productivity is higher under N- and lower under O3 deposition, with interactions with climatic conditions. Aboveground biomass was best correlated with growing-degree days for May (GDDMay). Nitrogen deposition increased biomass up to 60% in the highest treatment, and 30% in the lowest addition. Also belowground biomass showed a positive N-response. Ozone enrichment had no effect on biomass, and no interaction between O3 and N was observed. Growth response to N deposition was not correlated to GDDMay or precipitation, but indicated a cumulative effect over time. Productivity of subalpine grassland is tolerant to increasing ozone exposure, independent of N input and climatic drivers. N deposition rates at current critical loads, strongly increase the grassland yield.

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.

Global convergence in the temperature sensitivity of respiration at ecosystem level.

The respiratory release of carbon dioxide (CO(2)) from the land surface is a major flux in the global carbon cycle, antipodal to photosynthetic CO(2) uptake. Understanding the sensitivity of respiratory processes to temperature is central for quantifying the climate-carbon cycle feedback. We approximated the sensitivity of terrestrial ecosystem respiration to air temperature (Q(10)) across 60 FLUXNET sites with the use of a methodology that circumvents confounding effects. Contrary to previous findings, our results suggest that Q(10) is independent of mean annual temperature, does not differ among biomes, and is confined to values around 1.4 +/- 0.1. The strong relation between photosynthesis and respiration, by contrast, is highly variable among sites. The results may partly explain a less pronounced climate-carbon cycle feedback than suggested by current carbon cycle climate models.

Soil Respiration in European Grasslands in Relation to Climate and Assimilate Supply.

Soil respiration constitutes the second largest flux of carbon (C) between terrestrial ecosystems and the atmosphere. This study provides a synthesis of soil respiration (R(s)) in 20 European grasslands across a climatic transect, including ten meadows, eight pastures and two unmanaged grasslands. Maximum rates of R(s) (R(s(max) )), R(s) at a reference soil temperature (10°C; R(s(10) )) and annual R(s) (estimated for 13 sites) ranged from 1.9 to 15.9 µmol CO(2) m(-2) s(-1), 0.3 to 5.5 µmol CO(2) m(-2) s(-1) and 58 to 1988 g C m(-2) y(-1), respectively. Values obtained for Central European mountain meadows are amongst the highest so far reported for any type of ecosystem. Across all sites R(s(max) ) was closely related to R(s(10) ).Assimilate supply affected R(s) at timescales from daily (but not necessarily diurnal) to annual. Reductions of assimilate supply by removal of aboveground biomass through grazing and cutting resulted in a rapid and a significant decrease of R(s). Temperature-independent seasonal fluctuations of R(s) of an intensively managed pasture were closely related to changes in leaf area index (LAI). Across sites R(s(10) ) increased with mean annual soil temperature (MAT), LAI and gross primary productivity (GPP), indicating that assimilate supply overrides potential acclimation to prevailing temperatures. Also annual R(s) was closely related to LAI and GPP. Because the latter two parameters were coupled to MAT, temperature was a suitable surrogate for deriving estimates of annual R(s) across the grasslands studied. These findings contribute to our understanding of regional patterns of soil C fluxes and highlight the importance of assimilate supply for soil CO(2) emissions at various timescales.

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.

Bi-directional soil/atmosphere N2 O exchange over two mown grassland systems with contrasting management practices.

Nitrous oxide (N2 O) fluxes from soil under mown grassland were monitored using static chambers over three growing seasons in intensively and extensively managed systems in Central Switzerland. Emissions were largest following the application of mineral (NH4 NO3 ) fertilizer, but there were also substantial emissions following cattle slurry application, after grass cuts and during the thawing of frozen soil. Continuous flux sampling, using automatic chambers, showed marked diurnal patterns in N2 O fluxes during emission peaks, with highest values in the afternoon. Net uptake fluxes of N2 O and subambient N2 O concentrations in soil open pore space were frequently measured on both fields. Flux integration over 2.5 years yields a cumulated emission of +4.7 kgN2 O-N ha-1 for the intensively managed field, equivalent to an average emission factor of 1.1%, and a small net sink activity of -0.4 kg N2 O-N ha-1 for the unfertilized system. The data suggest the existence of a consumption mechanism for N2 O in dry, areated soil conditions, which cannot be explained by conventional anaerobic denitrification. The effect of fertilization on greenhouse gas budgets of grassland at the ecosystem level is discussed.