Potential pre-industrial-like new particle formation induced by pure biogenic organic vapors in Finnish peatland.

Pure biogenic new particle formation (NPF) induced by highly oxygenated organic molecules (HOMs) could be an important mechanism for pre-industrial aerosol formation. However, it has not been unambiguously confirmed in the ambient due to the scarcity of truly pristine continental locations in the present-day atmosphere or the lack of chemical characterization of NPF precursors. Here, we report ambient observations of pure biogenic HOM-driven NPF over a peatland in southern Finland. Meteorological decoupling processes formed an "air pocket" (i.e., a very shallow surface layer) at night and favored NPF initiated entirely by biogenic HOM from this peatland, whose atmospheric environment closely resembles that of the pre-industrial era. Our study sheds light on pre-industrial aerosol formation, which represents the baseline for estimating the impact of present and future aerosol on climate, as well as on future NPF, the features of which may revert toward pre-industrial-like conditions due to air pollution mitigation.

Radiation and temperature drive diurnal variation of aerobic methane emissions from Scots pine canopy.

Methane emissions from plant foliage may play an important role in the global methane cycle, but their size and the underlying source processes remain poorly understood. Here, we quantify methane fluxes from the shoots of Scots pine trees, a dominant tree species in boreal forests, to identify source processes and environmental drivers, and we evaluate whether these fluxes can be constrained at the ecosystem-level by eddy covariance flux measurements. We show that shoot-level measurements conducted in forest, garden, or greenhouse settings; on mature trees and saplings; manually and with an automated CO2-, temperature-, and water-controlled chamber system; and with multiple methane analyzers all resulted in comparable daytime fluxes (0.144 ± 0.019 to 0.375 ± 0.074 nmol CH4 g-1 foliar d.w. h-1). We further find that these emissions exhibit a pronounced diurnal cycle that closely follows photosynthetically active radiation and is further modulated by temperature. These diurnal patterns indicate that methane production is associated with diurnal cycle of sunlight, indicating that this production is either a byproduct of photosynthesis-associated biochemical reactions (e.g., the methionine cycle) or produced through nonenzymatic photochemical reactions in plant biomass. Moreover, we identified a light-dependent component in stand-level methane fluxes, which showed order-of-magnitude agreement with shoot-level measurements (0.968 ± 0.031 nmol CH4 g-1 h-1) and which provides an upper limit for shoot methane emissions.

A single-point modeling approach for the intercomparison and evaluation of ozone dry deposition across chemical transport models (Activity 2 of AQMEII4).

A primary sink of air pollutants and their precursors is dry deposition. Dry deposition estimates differ across chemical transport models, yet an understanding of the model spread is incomplete. Here, we introduce Activity 2 of the Air Quality Model Evaluation International Initiative Phase 4 (AQMEII4). We examine 18 dry deposition schemes from regional and global chemical transport models as well as standalone models used for impact assessments or process understanding. We configure the schemes as single-point models at eight Northern Hemisphere locations with observed ozone fluxes. Single-point models are driven by a common set of site-specific meteorological and environmental conditions. Five of eight sites have at least 3 years and up to 12 years of ozone fluxes. The interquartile range across models in multiyear mean ozone deposition velocities ranges from a factor of 1.2 to 1.9 annually across sites and tends to be highest during winter compared with summer. No model is within 50 % of observed multiyear averages across all sites and seasons, but some models perform well for some sites and seasons. For the first time, we demonstrate how contributions from depositional pathways vary across models. Models can disagree with respect to relative contributions from the pathways, even when they predict similar deposition velocities, or agree with respect to the relative contributions but predict different deposition velocities. Both stomatal and nonstomatal uptake contribute to the large model spread across sites. Our findings are the beginning of results from AQMEII4 Activity 2, which brings scientists who model air quality and dry deposition together with scientists who measure ozone fluxes to evaluate and improve dry deposition schemes in the chemical transport models used for research, planning, and regulatory purposes.

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.

Modeled production, oxidation, and transport processes of wetland methane emissions in temperate, boreal, and Arctic regions.

Wetlands are the largest natural source of methane (CH4 ) to the atmosphere. The eddy covariance method provides robust measurements of net ecosystem exchange of CH4 , but interpreting its spatiotemporal variations is challenging due to the co-occurrence of CH4 production, oxidation, and transport dynamics. Here, we estimate these three processes using a data-model fusion approach across 25 wetlands in temperate, boreal, and Arctic regions. Our data-constrained model-iPEACE-reasonably reproduced CH4 emissions at 19 of the 25 sites with normalized root mean square error of 0.59, correlation coefficient of 0.82, and normalized standard deviation of 0.87. Among the three processes, CH4 production appeared to be the most important process, followed by oxidation in explaining inter-site variations in CH4 emissions. Based on a sensitivity analysis, CH4 emissions were generally more sensitive to decreased water table than to increased gross primary productivity or soil temperature. For periods with leaf area index (LAI) of ≥20% of its annual peak, plant-mediated transport appeared to be the major pathway for CH4 transport. Contributions from ebullition and diffusion were relatively high during low LAI (<20%) periods. The lag time between CH4 production and CH4 emissions tended to be short in fen sites (3 ± 2 days) and long in bog sites (13 ± 10 days). Based on a principal component analysis, we found that parameters for CH4 production, plant-mediated transport, and diffusion through water explained 77% of the variance in the parameters across the 19 sites, highlighting the importance of these parameters for predicting wetland CH4 emissions across biomes. These processes and associated parameters for CH4 emissions among and within the wetlands provide useful insights for interpreting observed net CH4 fluxes, estimating sensitivities to biophysical variables, and modeling global CH4 fluxes.

Atmospheric observations suggest methane emissions in north-eastern China growing with natural gas use.

The dramatic increase of natural gas use in China, as a substitute for coal, helps to reduce CO2 emissions and air pollution, but the climate mitigation benefit can be offset by methane leakage into the atmosphere. We estimate methane emissions from 2010 to 2018 in four regions of China using the GOSAT satellite data and in-situ observations with a high-resolution (0.1° × 0.1°) inverse model and analyze interannual changes of emissions by source sectors. We find that estimated methane emission over the north-eastern China region contributes the largest part (0.77 Tg CH4 yr-1) of the methane emission growth rate of China (0.87 Tg CH4 yr-1) and is largely attributable to the growth in natural gas use. The results provide evidence of a detectable impact on atmospheric methane observations by the increasing natural gas use in China and call for methane emission reductions throughout the gas supply chain and promotion of low emission end-use facilities.

Intercomparison of methods to estimate gross primary production based on CO2 and COS flux measurements.

Separating the components of ecosystem-scale carbon exchange is crucial in order to develop better models and future predictions of the terrestrial carbon cycle. However, there are several uncertainties and unknowns related to current photosynthesis estimates. In this study, we evaluate four different methods for estimating photosynthesis at a boreal forest at the ecosystem scale, of which two are based on carbon dioxide (CO2) flux measurements and two on carbonyl sulfide (COS) flux measurements. The CO2-based methods use traditional flux partitioning and artificial neural networks to separate the net CO2 flux into respiration and photosynthesis. The COS-based methods make use of a unique 5-year COS flux data set and involve two different approaches to determine the leaf-scale relative uptake ratio of COS and CO2 (LRU), of which one (LRUCAP) was developed in this study. LRUCAP was based on a previously tested stomatal optimization theory (CAP), while LRUPAR was based on an empirical relation to measured radiation. For the measurement period 2013-2017, the artificial neural network method gave a GPP estimate very close to that of traditional flux partitioning at all timescales. On average, the COS-based methods gave higher GPP estimates than the CO2-based estimates on daily (23% and 7% higher, using LRUPAR and LRUCAP, respectively) and monthly scales (20% and 3% higher), as well as a higher cumulative sum over 3 months in all years (on average 25% and 3% higher). LRUCAP was higher than LRU estimated from chamber measurements at high radiation, leading to underestimation of midday GPP relative to other GPP methods. In general, however, use of LRUCAP gave closer agreement with CO2-based estimates of GPP than use of LRUPAR. When extended to other sites, LRUCAP may be more robust than LRUPAR because it is based on a physiological model whose parameters can be estimated from simple measurements or obtained from the literature. In contrast, the empirical radiation relation in LRUPAR may be more site-specific. However, this requires further testing at other measurement sites.

Does growing atmospheric CO2 explain increasing carbon sink in a boreal coniferous forest?

The terrestrial net ecosystem productivity (NEP) has increased during the past three decades, but the mechanisms responsible are still unclear. We analyzed 17 years (2001-2017) of eddy-covariance measurements of NEP, evapotranspiration (ET) and light and water use efficiency from a boreal coniferous forest in Southern Finland for trends and inter-annual variability (IAV). The forest was a mean annual carbon sink (252 [ ± 42] gC m-2a-1 ), and NEP increased at rate +6.4-7.0 gC m-2a-1 (or ca. +2.5% a-1 ) during the period. This was attributed to the increasing gross-primary productivity GPP and occurred without detectable change in ET. The start of annual carbon uptake period was advanced by 0.7 d a-1 , and increase in GPP and NEP outside the main growing season contributed ca. one-third and one-fourth of the annual trend, respectively. Meteorological factors were responsible for the IAV of fluxes but did not explain the long-term trends. The growing season GPP trend was strongest in ample light during the peak growing season. Using a multi-layer ecosystem model, we showed that direct CO2 fertilization effect diminishes when moving from leaf to ecosystem, and only 30-40% of the observed ecosystem GPP increase could be attributed to CO2 . The increasing trend in leaf-area index (LAI), stimulated by forest thinning in 2002, was the main driver of the enhanced GPP and NEP of the mid-rotation managed forest. It also compensated for the decrease of mean leaf stomatal conductance with increasing CO2 and LAI, explaining the apparent proportionality between observed GPP and CO2 trends. The results emphasize that attributing trends to their physical and physiological drivers is challenged by strong IAV, and uncertainty of LAI and species composition changes due to the dynamic flux footprint. The results enlighten the underlying mechanisms responsible for the increasing terrestrial carbon uptake in the boreal zone.

Uncovering the critical soil moisture thresholds of plant water stress for European ecosystems.

Understanding the critical soil moisture (SM) threshold (θcrit ) of plant water stress and land surface energy partitioning is a basis to evaluate drought impacts and improve models for predicting future ecosystem condition and climate. Quantifying the θcrit across biomes and climates is challenging because observations of surface energy fluxes and SM remain sparse. Here, we used the latest database of eddy covariance measurements to estimate θcrit across Europe by evaluating evaporative fraction (EF)-SM relationships and investigating the covariance between vapor pressure deficit (VPD) and gross primary production (GPP) during SM dry-down periods. We found that the θcrit and soil matric potential threshold in Europe are 16.5% and -0.7 MPa, respectively. Surface energy partitioning characteristics varied among different vegetation types; EF in savannas had the highest sensitivities to SM in water-limited stage, and the lowest in forests. The sign of the covariance between daily VPD and GPP consistently changed from positive to negative during dry-down across all sites when EF shifted from relatively high to low values. This sign of the covariance changed after longer period of SM decline in forests than in grasslands and savannas. Estimated θcrit from the VPD-GPP covariance method match well with the EF-SM method, showing this covariance method can be used to detect the θcrit . We further found that soil texture dominates the spatial variability of θcrit while shortwave radiation and VPD are the major drivers in determining the spatial pattern of EF sensitivities. Our results highlight for the first time that the sign change of the covariance between daily VPD and GPP can be used as an indicator of how ecosystems transition from energy to SM limitation. We also characterized the corresponding θcrit and its drivers across diverse ecosystems in Europe, an essential variable to improve the representation of water stress in land surface models.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.

Statistical upscaling of ecosystem CO2 fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties.

The regional variability in tundra and boreal carbon dioxide (CO2 ) fluxes can be high, complicating efforts to quantify sink-source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO2 fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO2 fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990-2015 from 148 terrestrial high-latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO2 fluxes and test the accuracy and uncertainty of different statistical models. CO2 fluxes were upscaled at relatively high spatial resolution (1 km2 ) across the high-latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE-focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO2 sink strength was larger in the boreal biome (observed and predicted average annual NEE -46 and -29 g C m-2  yr-1 , respectively) compared to tundra (average annual NEE +10 and -2 g C m-2  yr-1 ). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO2 budget, estimated using the annual NEE ensemble prediction, suggests the high-latitude region was on average an annual CO2 sink during 1990-2015, although uncertainty remains high.

Identifying dominant environmental predictors of freshwater wetland methane fluxes across diurnal to seasonal time scales.

While wetlands are the largest natural source of methane (CH4 ) to the atmosphere, they represent a large source of uncertainty in the global CH4 budget due to the complex biogeochemical controls on CH4 dynamics. Here we present, to our knowledge, the first multi-site synthesis of how predictors of CH4 fluxes (FCH4) in freshwater wetlands vary across wetland types at diel, multiday (synoptic), and seasonal time scales. We used several statistical approaches (correlation analysis, generalized additive modeling, mutual information, and random forests) in a wavelet-based multi-resolution framework to assess the importance of environmental predictors, nonlinearities and lags on FCH4 across 23 eddy covariance sites. Seasonally, soil and air temperature were dominant predictors of FCH4 at sites with smaller seasonal variation in water table depth (WTD). In contrast, WTD was the dominant predictor for wetlands with smaller variations in temperature (e.g., seasonal tropical/subtropical wetlands). Changes in seasonal FCH4 lagged fluctuations in WTD by ~17 ± 11 days, and lagged air and soil temperature by median values of 8 ± 16 and 5 ± 15 days, respectively. Temperature and WTD were also dominant predictors at the multiday scale. Atmospheric pressure (PA) was another important multiday scale predictor for peat-dominated sites, with drops in PA coinciding with synchronous releases of CH4 . At the diel scale, synchronous relationships with latent heat flux and vapor pressure deficit suggest that physical processes controlling evaporation and boundary layer mixing exert similar controls on CH4 volatilization, and suggest the influence of pressurized ventilation in aerenchymatous vegetation. In addition, 1- to 4-h lagged relationships with ecosystem photosynthesis indicate recent carbon substrates, such as root exudates, may also control FCH4. By addressing issues of scale, asynchrony, and nonlinearity, this work improves understanding of the predictors and timing of wetland FCH4 that can inform future studies and models, and help constrain wetland CH4 emissions.

Warming homogenizes apparent temperature sensitivity of ecosystem respiration.

Warming-induced carbon loss through terrestrial ecosystem respiration (Re) is likely getting stronger in high latitudes and cold regions because of the more rapid warming and higher temperature sensitivity of Re (Q 10). However, it is not known whether the spatial relationship between Q 10 and temperature also holds temporally under a future warmer climate. Here, we analyzed apparent Q 10 values derived from multiyear observations at 74 FLUXNET sites spanning diverse climates and biomes. We found warming-induced decline in Q 10 is stronger at colder regions than other locations, which is consistent with a meta-analysis of 54 field warming experiments across the globe. We predict future warming will shrink the global variability of Q 10 values to an average of 1.44 across the globe under a high emission trajectory (RCP 8.5) by the end of the century. Therefore, warming-induced carbon loss may be less than previously assumed because of Q 10 homogenization in a warming world.

Substantial hysteresis in emergent temperature sensitivity of global wetland CH4 emissions.

Wetland methane (CH4) emissions ([Formula: see text]) are important in global carbon budgets and climate change assessments. Currently, [Formula: see text] projections rely on prescribed static temperature sensitivity that varies among biogeochemical models. Meta-analyses have proposed a consistent [Formula: see text] temperature dependence across spatial scales for use in models; however, site-level studies demonstrate that [Formula: see text] are often controlled by factors beyond temperature. Here, we evaluate the relationship between [Formula: see text] and temperature using observations from the FLUXNET-CH4 database. Measurements collected across the globe show substantial seasonal hysteresis between [Formula: see text] and temperature, suggesting larger [Formula: see text] sensitivity to temperature later in the frost-free season (about 77% of site-years). Results derived from a machine-learning model and several regression models highlight the importance of representing the large spatial and temporal variability within site-years and ecosystem types. Mechanistic advancements in biogeochemical model parameterization and detailed measurements in factors modulating CH4 production are thus needed to improve global CH4 budget assessments.

Assessing model performance via the most limiting environmental driver in two differently stressed pine stands.

Climate change will impact forest productivity worldwide. Forecasting the magnitude of such impact, with multiple environmental stressors changing simultaneously, is only possible with the help of process-based models. In order to assess their performance, such models require careful evaluation against measurements. However, direct comparison of model outputs against observational data is often not reliable, as models may provide the right answers due to the wrong reasons. This would severely hinder forecasting abilities under unprecedented climate conditions. Here, we present a methodology for model assessment, which supplements the traditional output-to-observation model validation. It evaluates model performance through its ability to reproduce observed seasonal changes of the most limiting environmental driver (MLED) for a given process, here daily gross primary productivity (GPP). We analyzed seasonal changes of the MLED for GPP in two contrasting pine forests, the Mediterranean Pinus halepensis Mill. Yatir (Israel) and the boreal Pinus sylvestris L. Hyytiälä (Finland) from three years of eddy-covariance flux data. Then, we simulated the same period with a state-of-the-art process-based simulation model (LandscapeDNDC). Finally, we assessed if the model was able to reproduce both GPP observations and MLED seasonality. We found that the model reproduced the seasonality of GPP in both stands, but it was slightly overestimated without site-specific fine-tuning. Interestingly, although LandscapeDNDC properly captured the main MLED in Hyytiälä (temperature) and in Yatir (soil water availability), it failed to reproduce high-temperature and high-vapor pressure limitations of GPP in Yatir during spring and summer. We deduced that the most likely reason for this divergence is an incomplete description of stomatal behavior. In summary, this study validates the MLED approach as a model evaluation tool, and opens up new possibilities for model improvement.

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.

Effects of drought and meteorological forcing on carbon and water fluxes in Nordic forests during the dry summer of 2018.

The Nordic region was subjected to severe drought in 2018 with a particularly long-lasting and large soil water deficit in Denmark, Southern Sweden and Estonia. Here, we analyse the impact of the drought on carbon and water fluxes in 11 forest ecosystems of different composition: spruce, pine, mixed and deciduous. We assess the impact of drought on fluxes by estimating the difference (anomaly) between year 2018 and a reference year without drought. Unexpectedly, the evaporation was only slightly reduced during 2018 compared to the reference year at two sites while it increased or was nearly unchanged at all other sites. This occurred under a 40 to 60% reduction in mean surface conductance and the concurrent increase in evaporative demand due to the warm and dry weather. The anomaly in the net ecosystem productivity (NEP) was 93% explained by a multilinear regression with the anomaly in heterotrophic respiration and the relative precipitation deficit as independent variables. Most of the variation (77%) was explained by the heterotrophic component. Six out of 11 forests reduced their annual NEP with more than 50 g C m-2 yr-1 during 2018 as compared to the reference year. The NEP anomaly ranged between -389 and +74 g C m-2 yr-1 with a median value of -59 g C m-2 yr-1. 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 CO2 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 CO2 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'.