Contrasting impact of extreme soil and atmospheric dryness on the functioning of trees and forests

Recent studies indicate an increase in the frequency of extreme compound dryness days (days with both extreme soil AND air dryness) across central Europe in the future, with little information on their impact on the functioning of trees and forests. This study aims to quantify and assess the impact of extreme soil dryness, extreme air dryness, and extreme compound dryness on the functioning of trees and forests. For this, >15 years of ecosystem-level (carbon dioxide and water vapor fluxes) and 6-10 years of tree-level measurements (transpiration and growth) each from a montane mixed deciduous forest (CH-Lae) and a subalpine evergreen coniferous forest (CH-Dav) in Switzerland, is used. The results showed extreme air dryness limitation on CO2 fluxes and extreme soil dryness limitations on water vapor fluxes. Additionally, CH-Dav was mainly affected by extreme air dryness whereas CH-Lae was affected by both extreme soil dryness and extreme air dryness. The impact of extreme compound dryness on net CO2 uptake (about 75 % decrease) was more due to higher increased ecosystem respiration (40 % and 70 % increase at CH-Dav and CH-Lae, respectively) than decreased gross primary productivity (10 % and 40 % decrease at CH-Dav and CH-Lae, respectively). A significant negative impact on evapotranspiration and transpiration was only observed at CH-Lae during extreme soil and compound dryness (about 25 % decrease). Furthermore, with some differences, the tree-level impact on tree water deficit, transpiration, and growth were consistent with the ecosystem-level impact on carbon uptake and evapotranspiration. Finally, the impact of extreme dryness showed no significant relationship with tree allometry (diameter and height) but across different tree species. The projected future is likely to expose these forest areas to more extreme and frequent dryness conditions, thus compromising the functioning of trees and forests, thereby calling for management interventions to increase the adaptive capacity and resistance of these forests.

N2O and CH4 fluxes from intensively managed grassland: The importance of biological and environmental drivers vs. management

Agriculture is the main contributor to anthropogenic nitrous oxide (N2O) and methane (CH4) emissions. Therefore, mitigation options are urgently needed. In contrast to carbon dioxide, eddy covariance measurements of N2O and CH4 fluxes are still scarce, and thus little is known how environmental and biotic drivers as well as management affect the net N2O and CH4 exchange in grasslands. Thus, we investigated the most important drivers of net ecosystem N2O and CH4 fluxes in a temperate grassland, and continued a N2O mitigation experiment (increased clover proportion vs. fertilization with slurry). Random forest gap-filling models were able to capture intermittent emission peaks, performing better for half-hourly N2O than for CH4 fluxes. The unfertilized clover parcel (parcel B) continued to show lower N2O emissions (4.4 and 2.7 kg N2O-N ha-1 yr-1) compared to the fertilized parcel (parcel A; 6.9 and 5.9 kg N2O-N ha-1 yr-1) for 2019 and 2020, respectively. Tier 1 nitrogen (N) emission factors of 2.6 % and 1.9 % were observed at the fertilized parcel during the study period. Lower soil N concentrations indicated a lower N leaching risk at the clover than at the fertilized parcel. Annual CH4 emissions (including periods with sheep grazing) were similar from both parcels, and ranged from 25 to 38.5 kg CH4-C ha-1. The most important drivers of both N2O and CH4 fluxes were lagged precipitation and water filled pore space, but also management (for N2O from parcel B; CH4 from parcel A). Biotic variables such as vegetation height and leaf area index were important predictors for the N2O exchange, while grazing temporarily increased CH4 emissions. Overall, reducing N fertilization and increasing the legume proportion were effective N2O reduction measures. In particular, adjusting N fertilization to plant N demands can help to avoid high N2O emissions from grasslands.

Detection and attribution of an anomaly in terrestrial photosynthesis in Europe during the COVID-19 lockdown

Carbon dioxide (CO2) uptake by plant photosynthesis, referred to as gross primary production (GPP) at the ecosystem level, is sensitive to environmental factors, including pollutant exposure, pollutant uptake, and changes in the scattering of solar shortwave irradiance (SWin) - the energy source for photosynthesis. The 2020 spring lockdown due to COVID-19 resulted in improved air quality and atmospheric transparency, providing a unique opportunity to assess the impact of air pollutants on terrestrial ecosystem functioning. However, detecting these effects can be challenging as GPP is influenced by other meteorological drivers and management practices. Based on data collected from 44 European ecosystem-scale CO2 flux monitoring stations, we observed significant changes in spring GPP at 34 sites during 2020 compared to 2015-2019. Among these, 14 sites showed an increase in GPP associated with higher SWin, 10 sites had lower GPP linked to atmospheric and soil dryness, and seven sites were subjected to management practices. The remaining three sites exhibited varying dynamics, with one experiencing colder and rainier weather resulting in lower GPP, and two showing higher GPP associated with earlier spring melts. Analysis using the regional atmospheric chemical transport model (LOTOS-EUROS) indicated that the ozone (O3) concentration remained relatively unchanged at the research sites, making it unlikely that O3 exposure was the dominant factor driving the primary production anomaly. In contrast, SWin increased by 9.4 % at 36 sites, suggesting enhanced GPP possibly due to reduced aerosol optical depth and cloudiness. Our findings indicate that air pollution and cloudiness may weaken the terrestrial carbon sink by up to 16 %. Accurate and continuous ground-based observations are crucial for detecting and attributing subtle changes in terrestrial ecosystem functioning in response to environmental and anthropogenic drivers.

Monitoring of carbon-water fluxes at Eurasian meteorological stations using random forest and remote sensing

Simulating the carbon-water fluxes at more widely distributed meteorological stations based on the sparsely and unevenly distributed eddy covariance flux stations is needed to accurately understand the carbon-water cycle of terrestrial ecosystems. We established a new framework consisting of machine learning, determination coefficient (R2), Euclidean distance, and remote sensing (RS), to simulate the daily net ecosystem carbon dioxide exchange (NEE) and water flux (WF) of the Eurasian meteorological stations using a random forest model or/and RS. The daily NEE and WF datasets with RS-based information (NEE-RS and WF-RS) for 3774 and 4427 meteorological stations during 2002-2020 were produced, respectively. And the daily NEE and WF datasets without RS-based information (NEE-WRS and WF-WRS) for 4667 and 6763 meteorological stations during 1983-2018 were generated, respectively. For each meteorological station, the carbon-water fluxes meet accuracy requirements and have quasi-observational properties. These four carbon-water flux datasets have great potential to improve the assessments of the ecosystem carbon-water dynamics.

Long-term changes in forest response to extreme atmospheric dryness

Atmospheric dryness, as indicated by vapor pressure deficit (VPD), has a strong influence on forest greenhouse gas exchange with the atmosphere. In this study, we used long-term (10-30 years) net ecosystem productivity (NEP) measurements from 60 forest sites across the world (1003 site-years) to quantify long-term changes in forest NEP resistance and NEP recovery in response to extreme atmospheric dryness. We tested two hypotheses: first, across sites differences in NEP resistance and NEP recovery of forests will depend on both the biophysical characteristics (i.e., leaf area index [LAI] and forest type) of the forest as well as on the local meteorological conditions of the site (i.e., mean VPD of the site), and second, forests experiencing an increasing trend in frequency and intensity of extreme dryness will show an increasing trend in NEP resistance and NEP recovery over time due to emergence of long-term ecological stress memory. We used a data-driven statistical learning approach to quantify NEP resistance and NEP recovery over multiple years. Our results showed that forest types, LAI, and median local VPD conditions explained over 50% of variance in both NEP resistance and NEP recovery, with drier sites showing higher NEP resistance and NEP recovery compared to sites with less atmospheric dryness. The impact of extreme atmospheric dryness events on NEP lasted for up to 3 days following most severe extreme events in most forests, indicated by an NEP recovery of less than 100%. We rejected our second hypothesis as we found no consistent relationship between trends of extreme VPD with trends in NEP resistance and NEP recovery across different forest sites, thus an increase in atmospheric dryness as it is predicted might not increase the resistance or recovery of forests in terms of NEP.

Joint optimization of land carbon uptake and albedo can help achieve moderate instantaneous and long-term cooling effects

Both carbon dioxide uptake and albedo of the land surface affect global climate. However, climate change mitigation by increasing carbon uptake can cause a warming trade-off by decreasing albedo, with most research focusing on afforestation and its interaction with snow. Here, we present carbon uptake and albedo observations from 176 globally distributed flux stations. We demonstrate a gradual decline in maximum achievable annual albedo as carbon uptake increases, even within subgroups of non-forest and snow-free ecosystems. Based on a paired-site permutation approach, we quantify the likely impact of land use on carbon uptake and albedo. Shifting to the maximum attainable carbon uptake at each site would likely cause moderate net global warming for the first approximately 20 years, followed by a strong cooling effect. A balanced policy co-optimizing carbon uptake and albedo is possible that avoids warming on any timescale, but results in a weaker long-term cooling effect.

Greenhouse gas fluxes (CO2, N2O and CH4) of pea and maize during two cropping seasons: Drivers, budgets, and emission factors for nitrous oxide

Agriculture contributes considerably to the increase of global greenhouse gas (GHG) emissions. Hence, magnitude and drivers of temporal variations in carbon dioxide (CO2), nitrous oxide (N2O) and methane (CH4) fluxes in croplands are urgently needed to develop sustainable, climate-smart agricultural practices. However, our knowledge of GHG fluxes from croplands is still very limited. The eddy covariance technique was used to quantify GHG budgets and N2O emission factors (EF) for pea and maize in Switzerland. The random forest technique was applied for gap-filling N2O and CH4 fluxes as well as to determine the relevance of environmental, vegetation vs. management drivers of the GHG fluxes during two cropping seasons. Environmental (i.e., net radiation, soil water content, soil temperature) and vegetation drivers (i.e., vegetation height) were more important drivers for GHG fluxes at field scale than time since management for the two crop species. Both crops acted as GHG sinks between sowing and harvest, clearly dominated by net CO2 fluxes, while CH4 emissions were negligible. However, considerable N2O emissions occurred in both crop fields early in the season when crops were still establishing. N2O fluxes in both crops were small later in the season when vegetation was tall, despite high soil water contents and temperatures. Results clearly show a strong and highly dynamic microbial-plant competition for N driving N2O fluxes at the field scale. The total loss was 1.4 kg N2O-N ha-1 over 55 days for pea and 4.8 kg N2O-N ha-1 over 127 days for maize. EFs of N2O were 1.5 % (pea) and 4.4 % (maize) during the cropping seasons, clearly exceeding the IPCC Tier 1 EF for N2O. Thus, sustainable, climate-smart agriculture needs to consider crop phenology and better adapt N supply to crop N demand for growth, particularly during the early cropping season when competition for N between establishing crops and soil microorganisms modulates N2O losses.

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.

Ecosystem transpiration and evaporation: Insights from three water flux partitioning methods across FLUXNET sites

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate.

Physiological response of Swiss ecosystems to 2018 drought across plant types and elevation

Using five eddy covariance flux sites (two forests and three grasslands), we investigated ecosystem physiological responses to the 2018 drought across elevational gradients in Switzerland. Flux measurements showed that at lower elevation sites (below 1000 m.a.s.l.; grassland and mixed forest) annual ecosystem productivity (GPP) declined by approximately 20% compared to the previous 2 years (2016 and 2017), which led to a reduced annual net ecosystem productivity (NEP). At the high elevation sites, however, GPP increased by approximately 14% and as a result NEP increased in the alpine and montane grasslands, but not in the subalpine coniferous forest. There, increased ecosystem respiration led to a reduced annual NEP, despite increased GPP and lengthening of the growing period. Among all ecosystems, the coniferous forest showed the most pronounced negative stomatal response to atmospheric dryness (i.e. vapour pressure deficit, VPD) that resulted in a decline in surface conductance and an increased water-use efficiency during drought. While increased temperature enhanced the water-use efficiency of both forests, de-coupling of GPP from evapotranspiration at the low-elevation grassland site negatively affected water-use efficiency due to non-stomatal reductions in photosynthesis. Our results show that hot droughts (such as in 2018) lead to different responses across plants types, and thus ecosystems. Particularly grasslands at lower elevations are the most vulnerable ecosystems to negative impacts of future drought in Switzerland. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Altered energy partitioning across terrestrial ecosystems in the European drought year 2018

Drought and heat events, such as the 2018 European drought, interact with the exchange of energy between the land surface and the atmosphere, potentially affecting albedo, sensible and latent heat fluxes, as well as CO₂ exchange. Each of these quantities may aggravate or mitigate the drought, heat, their side effects on productivity, water scarcity and global warming. We used measurements of 56 eddy covariance sites across Europe to examine the response of fluxes to extreme drought prevailing most of the year 2018 and how the response differed across various ecosystem types (forests, grasslands, croplands and peatlands). Each component of the surface radiation and energy balance observed in 2018 was compared to available data per site during a reference period 2004-2017. Based on anomalies in precipitation and reference evapotranspiration, we classified 46 sites as drought affected. These received on average 9% more solar radiation and released 32% more sensible heat to the atmosphere compared to the mean of the reference period. In general, drought decreased net CO₂ uptake by 17.8%, but did not significantly change net evapotranspiration. The response of these fluxes differed characteristically between ecosystems; in particular, the general increase in the evaporative index was strongest in peatlands and weakest in croplands. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

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.

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

The FLUXNET2015 dataset provides ecosystem-scale data on CO₂, 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.

Measurement(s)	net ecosystem exchange • carbon dioxide • water • energy
Technology Type(s)	eddy covariance • measurement device
Sample Characteristic - Environment	terrestrial biome • atmosphere
Sample Characteristic - Location	Earth (planet)

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12295910

Assessment of spatial variability of multiple ecosystem services in grasslands of different intensities

Grasslands provide multiple Ecosystem Services (ES) such as forage provision, carbon sequestration or habitat provision. Knowledge about the trade-offs between these ES is of great importance for grassland management. Yet, the outcome of different management strategies on ES provision is highly uncertain due to spatial variability. We aim to characterize the provision (level and spatial variability) of grassland ES under various management strategies. To do so, we combine empirical data for multiple ES with spatially explicit census data on land use intensities. We analyzed the variations of five ES (forage provision, climate regulation, pollination, biodiversity conservation and outdoor recreation) using data from biodiversity fieldwork, experimental plots for carbon as well as social network data from Flickr. These data were used to calculate the distribution of modelled individual and multiple ES values from different grassland management types in a Swiss case study region using spatial explicit information for 17,383 grassland parcels. Our results show that (1) management regime and intensity levels play an important role in ES provision but their impact depends on the ES. In general, extensive management, especially in pastures, favors all ES but forage provision, whereas intensive management favors only forage provision and outdoor recreation; (2) ES potential provision varies between parcels under the same management due to the influence of environmental drivers, related to topography and landscape structure; (3) there is a trade-offs between forage provision and other ES at the cantonal level but a synergy between forage provision and biodiversity conservation within the grassland categories, due to the negative impact of elevation on both ES. Information about multiple ES provision is key to support effective agri-environmental measures and information about the spatial variability can prevent uncertain outputs of decision-making processes.

Greenhouse gas fluxes over managed grasslands in Central Europe

Central European grasslands are characterized by a wide range of different management practices in close geographical proximity. Site-specific management strategies strongly affect the biosphere-atmosphere exchange of the three greenhouse gases (GHG) carbon dioxide (CO₂ ), nitrous oxide (N₂ O), and methane (CH₄ ). The evaluation of environmental impacts at site level is challenging, because most in situ measurements focus on the quantification of CO₂ exchange, while long-term N₂ O and CH₄ flux measurements at ecosystem scale remain scarce. Here, we synthesized ecosystem CO₂ , N₂ O, and CH₄ fluxes from 14 managed grassland sites, quantified by eddy covariance or chamber techniques. We found that grasslands were on average a CO₂ sink (-1,783 to -91 g CO₂  m^-2  year^-1 ), but a N₂ O source (18-638 g CO₂ -eq. m^-2  year^-1 ), and either a CH₄ sink or source (-9 to 488 g CO₂ -eq. m^-2  year^-1 ). The net GHG balance (NGB) of nine sites where measurements of all three GHGs were available was found between -2,761 and -58 g CO₂ -eq. m^-2  year^-1 , with N₂ O and CH₄ emissions offsetting concurrent CO₂ uptake by on average 21 ± 6% across sites. The only positive NGB was found for one site during a restoration year with ploughing. The predictive power of soil parameters for N₂ O and CH₄ fluxes was generally low and varied considerably within years. However, after site-specific data normalization, we identified environmental conditions that indicated enhanced GHG source/sink activity ("sweet spots") and gave a good prediction of normalized overall fluxes across sites. The application of animal slurry to grasslands increased N₂ O and CH₄ emissions. The N₂ O-N emission factor across sites was 1.8 ± 0.5%, but varied considerably at site level among the years (0.1%-8.6%). Although grassland management led to increased N₂ O and CH₄ emissions, the CO₂ sink strength was generally the most dominant component of the annual GHG budget.

Quantifying deforestation and forest degradation with thermal response

Deforestation and forest degradation cause the deterioration of resources and ecosystem services. However, there are still no operational indicators to measure forest status, especially for forest degradation. In the present study, we analysed the thermal response number (TRN, calculated by daily total net radiation divided by daily temperature range) of 163 sites including mature forest, disturbed forest, planted forest, shrubland, grassland, savanna vegetation and cropland. TRN generally increased with latitude, however the regression of TRN against latitude differed among vegetation types. Mature forests are superior as thermal buffers, and had significantly higher TRN than disturbed and planted forests. There was a clear boundary between TRN of forest and non-forest vegetation (i.e. grassland and savanna) with the exception of shrubland, whose TRN overlapped with that of forest vegetation. We propose to use the TRN of local mature forest as the optimal TRN (TRN_opt). A forest with lower than 75% of TRN_opt was identified as subjected to significant disturbance, and forests with 66% of TRN_opt was the threshold for deforestation within the absolute latitude from 30° to 55°. Our results emphasized the irreplaceable thermal buffer capacity of mature forest. TRN can be used for early warning of deforestation and degradation risk. It is therefore a valuable tool in the effort to protect forests and prevent deforestation.

Eddy covariance flux measurements of gaseous elemental mercury using cavity ring-down spectroscopy

A newly developed pulsed cavity ring-down spectroscopy (CRDS) system for measuring atmospheric gaseous elemental mercury (GEM) concentrations at high temporal resolution (25 Hz) was used to successfully conduct the first eddy covariance (EC) flux measurements of GEM. GEM is the main gaseous atmospheric form, and quantification of bidirectional exchange between the Earth's surface and the atmosphere is important because gas exchange is important on a global scale. For example, surface GEM emissions from natural sources, legacy emissions, and re-emission of previously deposited anthropogenic pollution may exceed direct primary anthropogenic emissions. Using the EC technique for flux measurements requires subsecond measurements, which so far has not been feasible because of the slow time response of available instrumentation. The CRDS system measured GEM fluxes, which were compared to fluxes measured with the modified Bowen ratio (MBR) and a dynamic flux chamber (DFC). Measurements took place near Reno, NV, in September and October 2012 encompassing natural, low-mercury (Hg) background soils and Hg-enriched soils. During nine days of measurements with deployment of Hg-enriched soil in boxes within 60 m upwind of the EC tower, the covariance of GEM concentration and vertical wind speed was measured, showing that EC fluxes over an Hg-enriched area were detectable. During three separate days of flux measurements over background soils (without Hg-enriched soils), no covariance was detected, indicating fluxes below the detection limit. When fluxes were measurable, they strongly correlated with wind direction; the highest fluxes occurred when winds originated from the Hg-enriched area. Comparisons among the three methods showed good agreement in direction (e.g., emission or deposition) and magnitude, especially when measured fluxes originated within the Hg-enriched soil area. EC fluxes averaged 849 ng m(-2) h(-1), compared to DFC fluxes of 1105 ng m(-2) h(-1) and MBR fluxes of 1309 ng m(-2) h(-1). This study demonstrated that a CRDS system can be used to measure GEM fluxes over Hg-enriched areas, with a conservative detection limit estimate of 32 ng m(-2) h(-1).

An ecosystem-scale perspective of the net land methanol flux: synthesis of micrometeorological flux measurements

Methanol is the second most abundant volatile organic compound in the troposphere and plays a significant role in atmospheric chemistry. While there is consensus about the dominant role of living plants as the major source and the reaction with OH as the major sink of methanol, global methanol budgets diverge considerably in terms of source/sink estimates reflecting uncertainties in the approaches used to model, and the empirical data used to separately constrain these terms. Here we compiled micrometeorological methanol flux data from eight different study sites and reviewed the corresponding literature in order to provide a first cross-site synthesis of the terrestrial ecosystem-scale methanol exchange and present an independent data-driven view of the land-atmosphere methanol exchange. Our study shows that the controls of plant growth on the production, and thus the methanol emission magnitude, and stomatal conductance on the hourly methanol emission variability, established at the leaf level, hold across sites at the ecosystem-level. Unequivocal evidence for bi-directional methanol exchange at the ecosystem scale is presented. Deposition, which at some sites even exceeds methanol emissions, represents an emerging feature of ecosystem-scale measurements and is likely related to environmental factors favouring the formation of surface wetness. Methanol may adsorb to or dissolve in this surface water and eventually be chemically or biologically removed from it. Management activities in agriculture and forestry are shown to increase local methanol emission by orders of magnitude; they are however neglected at present in global budgets. While contemporary net land methanol budgets are overall consistent with the grand mean of the micrometeorological methanol flux measurements, we caution that the present approach of simulating methanol emission and deposition separately is prone to opposing systematic errors and does not allow taking full advantage of the rich information content of micrometeorological flux measurements.

An ecosystem-scale perspective of the net land methanol flux: synthesis of micrometeorological flux measurements

Methanol is the second most abundant volatile organic compound in the troposphere and plays a significant role in atmospheric chemistry. While there is consensus about the dominant role of living plants as the major source and the reaction with OH as the major sink of methanol, global methanol budgets diverge considerably in terms of source/sink estimates reflecting uncertainties in the approaches used to model, and the empirical data used to separately constrain these terms. Here we compiled micrometeorological methanol flux data from eight different study sites and reviewed the corresponding literature in order to provide a first cross-site synthesis of the terrestrial ecosystem-scale methanol exchange and present an independent data-driven view of the land-atmosphere methanol exchange. Our study shows that the controls of plant growth on the production, and thus the methanol emission magnitude, and stomatal conductance on the hourly methanol emission variability, established at the leaf level, hold across sites at the ecosystem-level. Unequivocal evidence for bi-directional methanol exchange at the ecosystem scale is presented. Deposition, which at some sites even exceeds methanol emissions, represents an emerging feature of ecosystem-scale measurements and is likely related to environmental factors favouring the formation of surface wetness. Methanol may adsorb to or dissolve in this surface water and eventually be chemically or biologically removed from it. Management activities in agriculture and forestry are shown to increase local methanol emission by orders of magnitude; they are however neglected at present in global budgets. While contemporary net land methanol budgets are overall consistent with the grand mean of the micrometeorological methanol flux measurements, we caution that the present approach of simulating methanol emission and deposition separately is prone to opposing systematic errors and does not allow taking full advantage of the rich information content of micrometeorological flux measurements.

Methane and nitrous oxide exchange over a managed hay meadow

The methane (CH4) and nitrous oxide (N2O) exchange of a temperate mountain grassland near Neustift, Austria, was measured during 2010-2012 over a time period of 22 months using the eddy covariance method. Exchange rates of both compounds at the site were low, with 97% of all half-hourly CH4 and N2O fluxes ranging between ±200 and ±50 ng m-2 s-1, respectively. The meadow acted as a sink for both compounds during certain time periods, but was a clear source of CH4 and N2O on an annual timescale. Therefore, both gases contributed to an increase of the global warming potential (GWP), effectively reducing the sink strength in terms of CO2 equivalents of the investigated grassland site. In 2011, our best guess estimate showed a net greenhouse gas (GHG) sink of -32 g CO2 equ. m-2 yr-1 for the meadow, whereby 55% of the CO2 sink strength of -71 g CO2m-2 yr-1 was offset by CH4 (N2O) emissions of 7 (32) g CO2 equ. m-2 yr-1. When all data were pooled, the ancillary parameters explained 27 (42)% of observed CH4 (N2O) flux variability, and up to 62 (76)% on shorter timescales in-between management dates. In the case of N2O fluxes, we found the highest emissions at intermediate soil water contents and at soil temperatures close to 0 or above 14 °C. In comparison to CO2, H2O and energy fluxes, the interpretation of CH4 and N2O exchange was challenging due to footprint heterogeneity regarding their sources and sinks, uncertainties regarding post-processing and quality control. Our results emphasize that CH4 and N2O fluxes over supposedly well-aerated and moderately fertilized soils cannot be neglected when evaluating the GHG impact of temperate managed grasslands.

Tradeoffs between global warming and day length on the start of the carbon uptake period in seasonally cold ecosystems

It is well established that warming leads to longer growing seasons in seasonally cold ecosystems. Whether this goes along with an increase in the net ecosystem carbon dioxide (CO₂) uptake is much more controversial. We studied the effects of warming on the start of the carbon uptake period (CUP) of three mountain grasslands situated along an elevational gradient in the Alps. To this end we used a simple empirical model of the net ecosystem CO₂ exchange, calibrated and forced with multi-year empirical data from each site. We show that reductions in the quantity and duration of daylight associated with earlier snowmelts were responsible for diminishing returns, in terms of carbon gain, from longer growing seasons caused by reductions in daytime photosynthetic uptake and increases in nighttime losses of CO₂. This effect was less pronounced at high, compared to low, elevations, where the start of the CUP occurred closer to the summer solstice when changes in day length and incident radiation are minimal.

Acetaldehyde exchange above a managed temperate mountain grassland

An overview of acetaldehyde exchange above a managed temperate mountain grassland in Austria over four growing seasons is presented. The meadow acted as a net source of acetaldehyde in all four years, emitting between 7 and 28 mg C m-2 over the whole growing period. The cutting of the meadow resulted in huge acetaldehyde emission bursts on the day of harvesting or one day later. During undisturbed conditions, both uptake and emission fluxes were recorded. The bidirectional nature of acetaldehyde fluxes was also reflected by clear diurnal cycles during certain time periods, indicating strong deposition processes before the 1st cut and emission towards the end of the growing season. The analysis of acetaldehyde compensation points revealed a complex relationship between ambient acetaldehyde mixing ratios and respective fluxes, significantly influenced by multiple environmental parameters and variable throughout the year. As a major finding of this study, we identified both a positive and negative correlation between concentration and flux on a daily scale, where soil temperature and soil water content were the most significant factors in determining the direction of the slope. In turn, this bidirectional relationship on a daily scale resulted in compensation points between 0.40 ppbv and 0.54 ppbv, which could be well explained by collected ancillary data. We conclude that in order to model acetaldehyde fluxes at the site in Neustift on a daily scale over longer time periods, it is crucial to know the type of relationship, i.e. the direction of the slope, between mixing ratios and fluxes on a given day.

Qualitative and quantitative characterization of volatile organic compound emissions from cut grass

Mechanical wounding of plants triggers the release of a blend of reactive biogenic volatile organic compounds (BVOCs). During and after mowing and harvesting of managed grasslands, significant BVOC emissions have the potential to alter the physical and chemical properties of the atmosphere and lead to ozone and aerosol formation with consequences for regional air quality. We show that the amount and composition of BVOCs emitted per unit dry weight of plant material is comparable between laboratory enclosure measurements of artificially severed grassland plant species and in situ ecosystem-scale flux measurements above a temperate mountain grassland during and after periodic mowing and harvesting. The investigated grassland ecosystem emitted annually up to 130 mg carbon m(-2) in response to cutting and drying, the largest part being consistently represented by methanol and a blend of green leaf volatiles (GLV). In addition, we report the plant species-specific emission of furfural, terpenoid-like compounds (e.g., camphor), and sesquiterpenes from cut plant material, which may be used as tracers for the presence of given plant species in the ecosystem.

Carbonyl sulfide (COS) as a tracer for canopy photosynthesis, transpiration and stomatal conductance: potential and limitations

The theoretical basis for the link between the leaf exchange of carbonyl sulfide (COS), carbon dioxide (CO(2)) and water vapour (H(2)O) and the assumptions that need to be made in order to use COS as a tracer for canopy net photosynthesis, transpiration and stomatal conductance, are reviewed. The ratios of COS to CO(2) and H(2)O deposition velocities used to this end are shown to vary with the ratio of the internal to ambient CO(2) and H(2)O mole fractions and the relative limitations by boundary layer, stomatal and internal conductance for COS. It is suggested that these deposition velocity ratios exhibit considerable variability, a finding that challenges current parameterizations, which treat these as vegetation-specific constants. COS is shown to represent a better tracer for CO(2) than H(2)O. Using COS as a tracer for stomatal conductance is hampered by our present poor understanding of the leaf internal conductance to COS. Estimating canopy level CO(2) and H(2)O fluxes requires disentangling leaf COS exchange from other ecosystem sources/sinks of COS. We conclude that future priorities for COS research should be to improve the quantitative understanding of the variability in the ratios of COS to CO(2) and H(2)O deposition velocities and the controlling factors, and to develop operational methods for disentangling ecosystem COS exchange into contributions by leaves and other sources/sinks. To this end, integrated studies, which concurrently quantify the ecosystem-scale CO(2), H(2)O and COS exchange and the corresponding component fluxes, are urgently needed.

Leaf and ecosystem response to soil water availability in mountain grasslands

Climate change is expected to affect the Alps by increasing the frequency and intensity of summer drought events with negative impacts on ecosystem water resources. The response of CO₂ and H₂O exchange of a mountain grassland to natural fluctuations of soil water content was evaluated during 2001-2009. In addition, the physiological performance of individual mountain forb and graminoid plant species under progressive soil water shortage was explored in a laboratory drought experiment. During the 9-year study period the natural occurrence of moderately to extremely dry periods did not lead to substantial reductions in net ecosystem CO₂ exchange and evapotranspiration. Laboratory drought experiments confirmed that all the surveyed grassland plant species were insensitive to progressive soil drying until very low soil water contents (<0.01 m³ m^-3) were reached after several days of drought. In field conditions, such a low threshold was never reached. Re-watering after a short-term drought event (5±1 days) resulted in a fast and complete recovery of the leaf CO₂ and H₂O gas exchange of the investigated plant species. We conclude that the present-day frequency and intensity of dry periods does not substantially affect the functioning of the investigated grassland ecosystem. During dry periods the observed "water spending" strategy employed by the investigated mountain grassland species is expected to provide a cooling feedback on climate warming, but may have negative consequences for down-stream water users.

Biotic, abiotic and management controls on methanol exchange above a temperate mountain grassland

Methanol (CH3OH) fluxes were quantified above a managed temperate mountain grassland in the Stubai Valley (Tyrol, Austria) during the growing seasons 2008 and 2009. Half-hourly methanol fluxes were calculated by means of the virtual disjunct eddy covariance (vDEC) method using 3-dimensional wind data from a sonic anemometer and methanol volume mixing ratios measured with a proton-transfer-reaction mass spectrometer (PTR-MS). During (undisturbed) mature and growing phases methanol fluxes exhibited a clear diurnal cycle with close-to-zero fluxes during nighttime and emissions, up to 10 nmol m-2 s-1, which followed the diurnal course of radiation and air temperature. Management events were found to represent the largest perturbations of methanol exchange at the studied grassland ecosystem: Peak emissions of 144.5 nmol m-2 s-1 were found during/after cutting of the meadow reflecting the wounding of the plant material and subsequent depletion of the leaf internal aqueous methanol pools. After the application of organic fertilizer, elevated methanol emissions of up to 26.7 nmol m-2 s-1 were observed, likely reflecting enhanced microbial activity associated with the applied manure. Simple and multiple linear regression analyses revealed air temperature and radiation as the dominant abiotic controls, jointly explaining 47 % and 70 % of the variability in half-hourly and daily methanol fluxes. In contrast to published leaf-level laboratory studies, the surface conductance and the daily change in the amount of green plant area, used as ecosystem-scale proxies for stomatal conductance and growth, respectively, were found to exert only minor biotic controls on methanol exchange.

Deposition Fluxes of Terpenes over Grassland

Eddy covariance flux measurements were carried out for two subsequent vegetation periods above a temperate mountain grassland in an alpine valley using a proton-transfer-reaction - mass spectrometer (PTR-MS) and a PTR-time of flight - mass spectrometer (PTR-TOF). In 2008 and during the first half of the vegetation period 2009 the volume mixing ratios (VMRs) for the sum of monoterpenes (MTs) were typically well below 1 ppbv and neither MT emission nor deposition was observed. After a hailstorm in July 2009 an order of magnitude higher amount of terpenes was transported to the site from nearby coniferous forests causing elevated VMRs. As a consequence, deposition fluxes of terpenes to the grassland, which continued over a time period of several weeks without significant re-emission, were observed. For days without precipitation the deposition occurred at velocities close to the aerodynamic limit. In addition to monoterpene uptake, deposition fluxes of the sum of sesquiterpenes (SQTs) and the sum of oxygenated terpenes (OTs) were detected. Considering an entire growing season for the grassland (i.e., 1^st of April to 1^st of November), the cumulative carbon deposition of monoterpenes reached 276 mg C m^-2. This is comparable to the net carbon emission of methanol (329 mg C m^-2), which is the dominant non methane volatile organic compound (VOC) emitted from this site, during the same time period. It is suggested that deposition of monoterpenes to terrestrial ecosystems could play a more significant role in the reactive carbon budget than previously assumed.

Ecosystem-scale biosphere-atmosphere interactions of a hemiboreal mixed forest stand at Järvselja, Estonia

During two measurement campaigns, from August to September 2008 and 2009, we quantified the major ecosystem fluxes in a hemiboreal forest ecosystem in Järvselja, Estonia. The main aim of this study was to separate the ecosystem flux components and gain insight into the performance of a multi-species multi-layered tree stand. Carbon dioxide and water vapor fluxes were measured using the eddy covariance method above and below the canopy in conjunction with the microclimate. Leaf and soil contributions were quantified separately by cuvette and chamber measurements, including fluxes of carbon dioxide, water vapor, nitrogen oxides, nitrous oxide, methane, ozone, sulfur dioxide, and biogenic volatile organic compounds (isoprene and monoterpenes). The latter have been as well characterized for monoterpenes in detail. Based on measured atmospheric trace gas concentrations, the flux tower site can be characterized as remote and rural with low anthropogenic disturbances. Our results presented here encourage future experimental efforts to be directed towards year round integrated biosphere-atmosphere measurements and development of process-oriented models of forest-atmosphere exchange taking the special case of a multi-layered and multi-species tree stand into account. As climate change likely leads to spatial extension of hemiboreal forest ecosystems a deep understanding of the processes and interactions therein is needed to foster management and mitigation strategies.

Eddy covariance VOC emission and deposition fluxes above grassland using PTR-TOF

Eddy covariance (EC) is the preferable technique for flux measurements since it is the only direct flux determination method. It requires a continuum of high time resolution measurements (e.g. 5-20 Hz). For volatile organic compounds (VOC) soft ionization via proton transfer reaction has proven to be a quantitative method for real time mass spectrometry; here we use a proton transfer reaction time of flight mass spectrometer (PTR-TOF) for 10 Hz EC measurements of full mass spectra up to m/z 315. The mass resolution of the PTR-TOF enabled the identification of chemical formulas and separation of oxygenated and hydrocarbon species exhibiting the same nominal mass. We determined 481 ion mass peaks from ambient air concentration above a managed, temperate mountain grassland in Neustift, Stubai Valley, Austria. During harvesting we found significant fluxes of 18 compounds distributed over 43 ions, including protonated parent compounds, as well as their isotopes and fragments and VOC-H⁺ - water clusters. The dominant BVOC fluxes were methanol, acetaldehyde, ethanol, hexenal and other C₆ leaf wound compounds, acetone, acetic acid, monoterpenes and sequiterpenes. The smallest reliable fluxes we determined were less than 0.1 nmol m^-2 s^-1, as in the case of sesquiterpene emissions from freshly cut grass. Terpenoids, including mono- and sesquiterpenes, were also deposited to the grassland before and after the harvesting. During cutting, total VOC emission fluxes up to 200 nmolC m^-2 s^-1 were measured. Methanol emissions accounted for half of the emissions of oxygenated VOCs and a third of the carbon of all measured VOC emissions during harvesting.

Multiple constraints on grassland evapotranspiration: implications for closing the energy balance

When using the eddy covariance (EC) method for measuring the ecosystem-atmosphere exchange of sensible and latent heat, it is not uncommon to find that these two energy fluxes fall short of available energy by 20-30 %. As the causes for the energy imbalance are still under discussion, it is currently not clear how the energy balance should be closed. The objective of the present paper is to use independent measurements of evapotranspiration (ET) for empirically devising on how to best close the energy balance. To this end ET of a temperate mountain grassland was quantified during two measurement campaigns using both an open- and a closed-path EC system, lysimeters and an approach scaling up leaf-level stomatal conductance to canopy level transpiration. Our study showed that both EC systems underestimated ET measured independently by lysimeters and the up-scaling approach. Best correspondence to independently measured ET was achieved by assigning the entire energy imbalance to ET and by adjusting ET according to the average energy balance ratio during the first and second measurement campaign, respectively. Due to a large spatial variability in ET during the first measurement campaign and given large differences in spatial scale between the EC and the independent methods, we are more confident with the comparison of approaches during the second measurement campaign and thus recommend forcing energy balance closure by adjusting for the average energy balance ratio.

Dealing with disjunct concentration measurements in eddy covariance applications: a comparison of available approaches

Using proton transfer reaction mass spectrometry equipped with a quadrupol mass analyser to quantify the biosphere-atmosphere exchange of volatile organic compounds (VOC), concentrations of different VOC are measured sequentially. Depending on how many VOC species are targeted and their respective integration times, each VOC is measured at repeat rates on the order of a few seconds. This represents an order of magnitude longer sample interval compared to the standard eddy covariance (EC) method (5-20 Hz sampling rates). Here we simulate the effect of disjunct sampling on EC flux estimates by decreasing the time resolution of CO2 and H2O concentrations measured at 20 Hz above a temperate mountain grassland in the Austrian Alps. Fluxes for one month are calculated with the standard EC method and compared to fluxes calculated based on the disjunct data (1, 3 and 5 s sampling rates) using the following approaches: i) imputation of missing concentrations based on the nearest neighbouring samples (iDECnn), ii) imputation by linear interpolation (iDECli), and iii) virtual disjunct EC (vDEC), i.e. flux calculation based solely on the disjunct concentrations. It is shown that the two imputation methods result in additional low-pass filtering, longer lag times (as determined with the maximum cross-correlation method) and a flux loss of 3-30 % as compared to the standard EC method. A novel procedure, based on a transfer function approach, which specifically corrects for the effect of data treatment, was developed, resulting in improved correspondence (to within 2 %). The vDEC method yields fluxes which approximate the true (20 Hz) fluxes to within 3-7 % and it is this approach we recommend because it involves no additional empirical corrections. The only drawback of the vDEC method is the noisy nature of the cross-correlations, which poses problems with lag determination - practical approaches to overcome this limitation are discussed.

First eddy covariance flux measurements by PTR-TOF

The recently developed PTR-TOF instrument was evaluated to measure methanol fluxes emitted from grass land using the eddy covariance method. The high time resolution of the PTR-TOF allowed storing full mass spectra up to m/z 315 with a frequency of 10 Hz. Three isobaric ions were found at a nominal mass of m/z 33 due to the high mass resolving power of the PTR-TOF. Only one of the three peaks contributed to eddy covariance fluxes. The exact mass of this peak agrees well with the exact mass of protonated methanol (m/z 33.0335). The eddy covariance methanol fluxes measured with PTR-TOF were compared to virtual disjunct eddy covariance methanol fluxes simultaneously measured with a conventional PTR-MS. The methanol fluxes from both instruments show excellent agreement.

Estimating carbon dioxide fluxes from temperate mountain grasslands using broad-band vegetation indices

The broad-band normalised difference vegetation index (NDVI) and the simple ratio (SR) were calculated from measurements of reflectance of photosynthetically active and short-wave radiation at two temperate mountain grasslands in Austria and related to the net ecosystem CO₂ exchange (NEE) measured concurrently by means of the eddy covariance method. There was no significant statistical difference between the relationships of midday mean NEE with narrow- and broad-band NDVI and SR, measured during and calculated for that same time window, respectively. The skill of broad-band NDVI and SR in predicting CO₂ fluxes was higher for metrics dominated by gross photosynthesis and lowest for ecosystem respiration, with NEE in between. A method based on a simple light response model whose parameters were parameterised based on broad-band NDVI allowed to improve predictions of daily NEE and is suggested to hold promise for filling gaps in the NEE time series. Relationships of CO₂ flux metrics with broad-band NDVI and SR however generally differed between the two studied grassland sites indicting an influence of additional factors not yet accounted for.

Measuring eddy covariance fluxes of ozone with a slow-response analyser

Ozone (O₃) fluxes above a temperate mountain grassland were measured by means of the eddy covariance (EC) method using a slow-response O₃ analyser. The resultant flux loss was corrected for by a series of transfer functions which model the various sources of high- and, in particular, low-pass filtering. The resulting correction factors varied on average between 1.7 and 3.5 during night and day time, respectively. A cospectral analysis confirmed the accuracy of this approach. O₃ fluxes were characterised by a comparatively large random uncertainty, which during daytime typically amounted to 60 %. EC O₃ fluxes were compared against O₃ flux measurements made concurrently with the flux-gradient (FG) method. The two methods generally agreed well, except for a period between sun rise and early afternoon, when the FG method was suspected of being affected by the presence of photochemical sources/sinks. O₃ flux magnitudes and deposition velocities determined with the EC method compared nicely with the available literature from grassland studies. We conclude that our understanding of the causes and consequences of various sources of flux loss (associated with any EC system) has sufficiently matured so that also less-than-ideal instrumentation may be used in EC flux applications, albeit at the cost of relatively large empirical corrections.

On the consequences of the energy imbalance for calculating surface conductance to water vapour

The Penman-Monteith combination equation, which is most frequently used to derive the surface conductance to water vapour (G_s), implicitly assumes the energy balance to be closed. Any energy imbalance (positive or negative) will thus affect the calculated G_s. Using eddy covariance energy flux data from a temperate grassland and a desert shrub ecosystem we explored five possible approaches of closing the energy imbalance and show that calculated G_s may differ considerably between these five approaches depending on the relative magnitudes of sensible and latent heat fluxes, and the magnitude and sign of the energy imbalance. Based on our limited understanding of the nature of the energy imbalance, we tend to favour an approach which preserves the Bowen-ratio and closes the energy balance on a larger time scale.