Search

Soil moisture plays a key role in soil nutrient and carbon cycling; plant productivity; and energy, water, and greenhouse gas exchanges between the land and the atmosphere. The knowledge on drivers of spatiotemporal soil moisture dynamics in subarctic landscapes is limited. In this study, we used the Spatial Forest Hydrology (SpaFHy) model, in situ soil moisture data, and Sentinel-1 synthetic aperture radar (SAR)-based soil moisture estimates to explore spatiotemporal controls of soil moisture in a subarctic headwater catchment in northwestern Finland. The role of groundwater dynamics and lateral flow in soil moisture was studied through three groundwater model conceptualizations: (i) omission of groundwater storage and lateral flow, (ii) conceptual TOPMODEL approach based on topographic wetness index, and (iii) explicit 2D lateral groundwater flow. The model simulations were compared against continuous point soil moisture measurements, distributed manual measurements, and novel SAR-based soil moisture estimates available at high spatial and temporal resolutions. Based on model scenarios and model–data comparisons, we assessed when and where the lateral groundwater flow shapes shallow soil moisture and under which conditions soil moisture variability is driven more by local ecohydrology, i.e., the balance of infiltration, drainage, and evapotranspiration. The choice of groundwater flow model was shown to have a strong impact on modeled soil moisture dynamics within the catchment. All model conceptualizations captured the observed soil moisture dynamics in the upland forests, but accounting for the lateral groundwater flow was necessary to reproduce the saturated conditions common in the peatlands and occasionally in lowland forest grid cells. We further highlight the potential of integrating multi-scale observations with land surface and hydrological models. The results have implications for ecohydrological and biogeochemical processes, as well as for modeling hydrology and Earth system feedbacks in subarctic and boreal environments.

Wetlands are responsible for 20%–31% of global methane (CH4) emissions and account for a large source of uncertainty in the global CH4 budget. Data-driven upscaling of CH4 fluxes from eddy covariance measurements can provide new and independent bottom-up estimates of wetland CH4 emissions. Here, we develop a six-predictor random forest upscaling model (UpCH4), trained on 119 site-years of eddy covariance CH4 flux data from 43 freshwater wetland sites in the FLUXNET-CH4 Community Product. Network patterns in site-level annual means and mean seasonal cycles of CH4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash-Sutcliffe Efficiency ∼0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH4 emissions of 146 ± 43 TgCH4 y−1 for 2001–2018 which agrees closely with current bottom-up land surface models (102–181 TgCH4 y−1) and overlaps with top-down atmospheric inversion models (155–200 TgCH4 y−1). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH4 fluxes has the potential to produce realistic extra-tropical wetland CH4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid-to-arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions. The monthly wetland methane products gridded at 0.25° from UpCH4 are available via ORNL DAAC (https://doi.org/10.3334/ORNLDAAC/2253).

Abstract Expanding and intensifying agriculture has led to a loss of soil carbon. As agroecosystems cover over 40% of Earth’s land surface, they must be part of the solution put in action to mitigate climate change. Development of efficient management practices to maximize soil carbon retention is currently limited, in part, by a poor understanding of how plants, which input carbon to soil, and microbes, which determine its fate there, interact. Here we implement a diversity gradient by intercropping undersown species with barley in a large field trial, ranging from one to eight undersown species. We find that increasing plant diversity strengthens positive associations within the rhizosphere soil microbial community in relation to negative associations. These associations, in turn, enhance community carbon use efficiency. Jointly, our results highlight how increasing plant diversity in agriculture can be used as a management strategy to enhance carbon retention potential in agricultural soils.

The urgent need to mitigate climate change has evoked a broad interest in better understanding and estimating nitrous oxide (N2O) emissions from different ecosystems. Part of the uncertainty in N2O emission estimates still comes from an inadequate understanding of the temporal and small-scale spatial variability of N2O fluxes. Using 4.5 years of N2O flux data collected in a drained peatland forest with six automated chambers, we explored temporal and small-scale spatial variability of N2O fluxes. A random forest with conditional inference trees was used to find immediate and delayed relationships between N2O flux and environmental conditions across seasons and years. The spatiotemporal variation of the N2O flux was large, with daily mean N2O flux varying between −10 and +1760 µgN2Om-2h-1 and annual N2O budgets of different chambers between +60 and +2110 mgN2Om-2yr-1. Spatial differences in fluxes persisted through years of different environmental conditions. Soil moisture, water table level, and air temperature were the most important variables explaining the temporal variation of N2O fluxes. N2O fluxes responded to precipitation events with peak fluxes measured on average 4 d after peaks in soil moisture and water table level. The length of the time lags varied in space and between seasons indicating possible interactions with temperature and other soil conditions. The high temporal variation in N2O flux was related to (a) temporal variation in environmental conditions, with the highest N2O fluxes measured after summer precipitation events and winter soil freezing, and (b) to annually varying seasonal weather conditions, with the highest N2O emissions measured during wet summers and winters with discontinuous snow cover. Climate change may thus increase winter N2O emissions, which may be offset by lower summer N2O emissions in dry years. The high sensitivity of N2O fluxes to seasonal weather conditions suggests increasing variability in annual peatland forest N2O budgets as the frequency of extreme weather events, such as droughts, is predicted to increase.

The snowpack has a major influence on the land surface energy budget. Accurate simulation of the snowpack energy and radiation budget is challenging due to, e.g., effects of vegetation and topography, as well as limitations in the theoretical understanding of turbulent transfer in the stable boundary layer. Studies that evaluate snow, hydrology and land surface models against detailed observations of all surface energy balance components at high latitudes are scarce. In this study, we compared different configurations of the SURFEX land surface model against surface energy flux, snow depth and soil temperature observations from four eddy-covariance stations in Finland. The sites cover two different climate and snow conditions, representing the southern and northern subarctic zones, as well as the contrasting forest and peatland ecosystems typical for the boreal landscape. We tested different turbulent flux parameterizations implemented in the Crocus snowpack model. In addition, we examined common alternative approaches to conceptualize soil and vegetation, and we assessed their performance in simulating surface energy fluxes, snow conditions and soil thermal regime. Our results show that a stability correction function that increases the turbulent exchange under stable atmospheric conditions is imperative to simulate sensible heat fluxes over the peatland snowpacks and that realistic peat soil texture (soil organic content) parameterization greatly improves the soil temperature simulations. For accurate simulations of surface energy fluxes, snow and soil conditions in forests, an explicit vegetation representation is necessary. Moreover, we demonstrate the high sensitivity of surface fluxes to a poorly documented parameter involved in snow cover fraction computation. Although we focused on models within the SURFEX platform, the results have broader implications for choosing suitable turbulent flux parameterization and model structures depending on the potential use cases for high-latitude land surface modeling.

Wetland CH₄ emissions are among the most uncertain components of the global CH₄ budget. The complex nature of wetland CH₄ processes makes it challenging to identify causal relationships for improving our understanding and predictability of CH₄ emissions. In this study, we used the flux measurements of CH₄ from eddy covariance towers (30 sites from 4 wetlands types: bog, fen, marsh, and wet tundra) to construct a causality-constrained machine learning (ML) framework to explain the regulative factors and to capture CH₄ emissions at sub-seasonal scale. We found that soil temperature is the dominant factor for CH₄ emissions in all studied wetland types. Ecosystem respiration (CO₂) and gross primary productivity exert controls at bog, fen, and marsh sites with lagged responses of days to weeks. Integrating these asynchronous environmental and biological causal relationships in predictive models significantly improved model performance. More importantly, modeled CH₄ emissions differed by up to a factor of 4 under a +1°C warming scenario when causality constraints were considered. These results highlight the significant role of causality in modeling wetland CH₄ emissions especially under future warming conditions, while traditional data-driven ML models may reproduce observations for the wrong reasons. Our proposed causality-guided model could benefit predictive modeling, large-scale upscaling, data gap-filling, and surrogate modeling of wetland CH₄ emissions within earth system land models.

We compiled published peer-reviewed CO2, CH4, and N2O data on managed drained organic forest soils in boreal and temperate zones to revisit the current Tier 1 default emission factors (EFs) provided in the IPCC (2014) Wetlands Supplement: to see whether their uncertainty may be reduced; to evaluate possibilities for breaking the broad categories used for the IPCC EFs into more site-type-specific ones; and to inspect the potential relevance of a number of environmental variables for predicting the annual soil greenhouse gas (GHG) balances, on which the EFs are based. Despite a considerable number of publications applicable for compiling EFs being added, only modest changes were found compared to the Tier 1 default EFs. However, the more specific site type categories generated in this study showed narrower confidence intervals compared to the default categories. Overall, the highest CO2 EFs were found for temperate afforested agricultural lands and boreal forestry-drained sites with very low tree stand productivity. The highest CH4 EFs in turn prevailed in boreal nutrient-poor forests with very low tree stand productivity and temperate forests irrespective of nutrient status, while the EFs for afforested sites were low or showed a sink function. The highest N2O EFs were found for afforested agricultural lands and forestry-drained nutrient-rich sites. The occasional wide confidence intervals could be mainly explained by single or a few highly deviating estimates rather than the broadness of the categories applied. Our EFs for the novel categories were further supported by the statistical models connecting the annual soil GHG balances to site-specific soil nutrient status indicators, tree stand characteristics, and temperature-associated weather and climate variables. The results of this synthesis have important implications for EF revisions and national emission reporting, e.g. by the use of different categories for afforested sites and forestry-drained sites, and more specific site productivity categories based on timber production potential.

Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO2) fluxes in terrestrial ecosystems across the rapidly warming Arctic-boreal zone (ABZ) have provided valuable information but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic-boreal CO2 fluxes (ABCflux) that aggregates in situ measurements of terrestrial net ecosystem CO2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures), and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19% of the monthly observations), snow diffusion (3%) and eddy covariance (78%) techniques. The largest number of observations were collected during the climatological summer (June-August; 32%), and fewer observations were available for autumn (September-October; 25%), winter (December-February; 18%), and spring (March-May; 25%). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO2 fluxes and to better estimate the terrestrial ABZ CO2 budget. ABCflux is openly and freely available online (Virkkala et al., 2021b, 10.3334/ORNLDAAC/1934). Copyright:

Time series of wetland methane fluxes measured by eddy covariance require gap-filling to estimate daily, seasonal, and annual emissions. Gap-filling methane fluxes is challenging because of high variability and complex responses to multiple drivers. To date, there is no widely established gap-filling standard for wetland methane fluxes, with regards both to the best model algorithms and predictors. This study synthesizes results of different gap-filling methods systematically applied at 17 wetland sites spanning boreal to tropical regions and including all major wetland classes and two rice paddies. Procedures are proposed for: 1) creating realistic artificial gap scenarios, 2) training and evaluating gap-filling models without overstating performance, and 3) predicting half-hourly methane fluxes and annual emissions with realistic uncertainty estimates. Performance is compared between a conventional method (marginal distribution sampling) and four machine learning algorithms. The conventional method achieved similar median performance as the machine learning models but was worse than the best machine learning models and relatively insensitive to predictor choices. Of the machine learning models, decision tree algorithms performed the best in cross-validation experiments, even with a baseline predictor set, and artificial neural networks showed comparable performance when using all predictors. Soil temperature was frequently the most important predictor whilst water table depth was important at sites with substantial water table fluctuations, highlighting the value of data on wetland soil conditions. Raw gap-filling uncertainties from the machine learning models were underestimated and we propose a method to calibrate uncertainties to observations. The python code for model development, evaluation, and uncertainty estimation is publicly available. This study outlines a modular and robust machine learning workflow and makes recommendations for, and evaluates an improved baseline of, methane gap-filling models that can be implemented in multi-site syntheses or standardized products from regional and global flux networks (e.g., FLUXNET).

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

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.

Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20g g€¯S to 20g g€¯N) the spring onset of elevated CH4 emissions starts 3g€¯d earlier, and the CH4 emission season lasts 4g€¯d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at 10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

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.

Abstract The water table depth (WTD) in peatlands determines the soil carbon decomposition rate and influences vegetation growth, hence the above‐ground carbon assimilation. Here, we used satellite‐observed Solar‐Induced chlorophyll Fluorescence (SIF) as a proxy of Gross Primary Production (GPP) to investigate water‐related vegetation stress over northern peatlands. A linear model with interaction effects was used to relate short‐ and long‐term anomalies in SIF with WTD anomalies and the absolute WTD. Most locations showed the occurrence of drought and waterlogging stress though regions with exclusively waterlogging or drought stress were also detected. As a spatial median, minimal water‐related vegetation stress was found for a WTD of −0.22 m (short‐term) and −0.20 m (long‐term) (±0.01 m, 95% confidence interval of statistical uncertainty). The stress response observed with SIF is supported by an analysis of in situ GPP data. Our findings provide insight into how changes in WTD of northern peatlands could affect GPP under climate change.

Abstract Wetlands are responsible for 20%–31% of global methane (CH4) emissions and account for a large source of uncertainty in the global CH4 budget. Data‐driven upscaling of CH4 fluxes from eddy covariance measurements can provide new and independent bottom‐up estimates of wetland CH4 emissions. Here, we develop a six‐predictor random forest upscaling model (UpCH4), trained on 119 site‐years of eddy covariance CH4 flux data from 43 freshwater wetland sites in the FLUXNET‐CH4 Community Product. Network patterns in site‐level annual means and mean seasonal cycles of CH4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash‐Sutcliffe Efficiency ∼0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH4 emissions of 146 ± 43 TgCH4 y−1 for 2001–2018 which agrees closely with current bottom‐up land surface models (102–181 TgCH4 y−1) and overlaps with top‐down atmospheric inversion models (155–200 TgCH4 y−1). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH4 fluxes has the potential to produce realistic extra‐tropical wetland CH4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid‐to‐arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions. The monthly wetland methane products gridded at 0.25° from UpCH4 are available via ORNL DAAC (https://doi.org/10.3334/ORNLDAAC/2253).

Accurate estimates of net ecosystem CO2 exchange (NEE) would improve the understanding of natural carbon sources and sinks and their role in the regulation of global atmospheric carbon. In this work, we use and compare the random forest (RF) and the gradient boosting (GB) machine learning (ML) methods for predicting year-round 6 h NEE over 1996–2018 in a pine-dominated boreal forest in southern Finland and analyze the predictability of NEE. Additionally, aggregation to weekly NEE values was applied to get information about longer term behavior of the method. The meteorological ERA5 reanalysis variables were used as predictors. Spatial and temporal neighborhood (predictor lagging) was used to provide the models more data to learn from, which was found to improve considerably the accuracy of both ML approaches compared to using only the nearest grid cell and time step. Both ML methods can explain temporal variability of NEE in the observational site of this study with meteorological predictors, but the GB method was more accurate. Only minor signs of overfitting could be detected for the GB algorithm when redundant variables were included. The accuracy of the approaches, measured mainly using cross-validated R2 score between the model result and the observed NEE, was high, reaching a best estimate value of 0.92 for GB and 0.88 for RF. In addition to the standard RF approach, we recommend using GB for modeling the CO2 fluxes of the ecosystems due to its potential for better performance.

The subarctic landscape consists of a mosaic of forest, peatland, and aquatic ecosystems and their ecotones. The carbon (C) exchange between ecosystems and the atmosphere through carbon dioxide (CO2) and methane (CH4) fluxes varies spatially and temporally among these ecosystems. Our study area in Kaamanen in northern Finland covered 7 km2 of boreal subarctic landscape with upland forest, open peatland, pine bogs, and lakes. We measured the CO2 and CH4 fluxes with eddy covariance and chambers between June 2017 and June 2019 and studied the C flux responses to varying meteorological conditions. The landscape area was an annual CO2 sink of −45 ± 22 and −33 ± 23 g C m−2 and a CH4 source of 3.0 ± 0.2 and 2.7 ± 0.2 g C m−2 during the first and second study years, respectively. The pine forest had the largest contribution to the landscape-level CO2 sink, −126 ± 21 and −101 ± 19 g C m−2, and the fen to the CH4 emissions, 7.8 ± 0.2 and 6.3 ± 0.3 g C m−2, during the first and second study years, respectively. The lakes within the area acted as CO2 and CH4 sources to the atmosphere throughout the measurement period, and a lake located downstream from the fen with organic sediment showed 4-fold fluxes compared to a mineral sediment lake. The annual C balances were affected most by the rainy peak growing season in 2017, the warm summer in 2018, and a heatwave and drought event in July 2018. The rainy period increased ecosystem respiration (ER) in the pine forest due to continuously high soil moisture content, and ER was on a level similar to the following, notably warmer, summer. A corresponding ER response to abundant precipitation was not observed for the fen ecosystem, which is adapted to high water table levels, and thus a higher ER sum was observed during the warm summer 2018. During the heatwave and drought period, similar responses were observed for all terrestrial ecosystems, with decreased gross primary productivity and net CO2 uptake, caused by the unfavourable growing conditions and plant stress due to the soil moisture and vapour pressure deficits. Additionally, the CH4 emissions from the fen decreased during and after the drought. However, the timing and duration of drought effects varied between the fen and forest ecosystems, as C fluxes were affected sooner and had a shorter post-drought recovery time in the fen than forest. The differing CO2 flux response to weather variations showed that terrestrial ecosystems can have a contrasting impact on the landscape-level C balance in a changing climate, even if they function similarly most of the time.

Abstract Climate change mitigation requires, besides reductions in greenhouse gas emissions, actions to increase carbon sinks in terrestrial ecosystems. A key measurement method for quantifying such sinks and calibrating models is the eddy covariance technique, but it requires imputation, or gap-filling, of missing data for determination of annual carbon balances of ecosystems. Previous comparisons of gap-filling methods have concluded that commonly used methods, such as marginal distribution sampling (MDS), do not have a significant impact on the carbon balance estimate. By analyzing an extensive, global data set, we show that MDS causes significant carbon balance errors for northern (latitude $$>60^\circ$$ > 60 ∘ ) sites. MDS systematically overestimates the carbon dioxide (CO $$_2$$ 2 ) emissions of carbon sources and underestimates the CO $$_2$$ 2 sequestration of carbon sinks. We also reveal reasons for these biases and show how a machine learning method called extreme gradient boosting or a modified implementation of MDS can be used to substantially reduce the northern site bias.

Vegetation phenology, which refers to the seasonal changes in plant physiology, biomass and plant cover, is affected by many abiotic factors, such as precipitation, temperature and water availability. Phenology is also associated with the carbon dioxide (CO2) exchange between ecosystems and the atmosphere. We employed digital cameras to monitor the vegetation phenology of three northern boreal peatlands during five growing seasons. We derived a greenness index (green chromatic coordinate, GCC) from the images and combined the results with measurements of CO2 flux, air temperature and high-resolution satellite data (Sentinel-2). From the digital camera images it was possible to extract greenness dynamics on the vegetation community and even species level. The highest GCC and daily maximum gross photosynthetic production (GPPmax) were observed at the site with the highest nutrient availability and richest vegetation. The short-term temperature response of GCC depended on temperature and varied among the sites and months. Although the seasonal development and year-to-year variation in GCC and GPPmax showed consistent patterns, the short-term variation in GPPmax was explained by GCC only during limited periods. GCC clearly indicated the main phases of the growing season, and peatland vegetation showed capability to fully compensate for the impaired growth resulting from a late growing season start. The GCC data derived from Sentinel-2 and digital cameras showed similar seasonal courses, but a reliable timing of different phenological phases depended upon the temporal coverage of satellite data.

We describe a new coordination activity and initial results for a global synthesis of eddy covariance CH4 flux measurements.

The occurrence of new particle formation (NPF) events detected in a coastal agricultural site, at Qvidja, in Southwestern Finland, was investigated using the data measured with a nitrate ion-based chemical-ionization atmospheric-pressure-interface time-of-flight (CI-APi-TOF) mass spectrometer. The binned positive matrix factorization method (binPMF) was applied to the measured spectra. It resulted in eight factors describing the time series of ambient gas and cluster composition at Qvidja during spring 2019. The most interesting factors related to the observed NPF events were the two factors with the highest mass-to-charge ratios, numbered 7 and 8, both having profiles with patterns of highly oxygenated organic molecules with one nitrogen atom. It was observed that factor 7 had elevated intensities during the NPF events. A variable with an even better connection to the observed NPF events is fF7, which denotes the fraction of the total spectra within the studied mass-to-charge ratio range between 169 and 450 Th being in a form of factor 7. Values of fF7 higher than 0.50±0.05 were observed during the NPF events, of which durations also correlated with the duration of fF7 exceeding this critical value. It was also observed that factor 8 acts like a precursor for factor 7 with solar radiation and that the formation of factor 8 is associated with ozone levels.

Boreal upland forests are generally considered methane (CH4) sinks due to the predominance of CH4 oxidizing bacteria over the methanogenic archaea. However, boreal upland forests can temporarily act as CH4 sources during wet seasons or years. From a landscape perspective and in annual terms, this source can be significant as weather conditions may cause flooding, which can last a considerable proportion of the active season and because often, the forest coverage within a typical boreal catchment is much higher than that of wetlands. Processes and conditions which change mineral soils from acting as a weak sink to a strong source are not well understood. We measured soil CH4 fluxes from 20 different points from regularly irrigated and control plots during two growing seasons. We also estimated potential CH4 production and oxidation rates in different soil layers and performed a laboratory experiment, where soil microcosms were subjected to different moisture levels and glucose addition simulating the fresh labile carbon (C) source from root exudates. The aim was to find the key controlling factors and conditions for boreal upland soil CH4 production. Probably due to long dry periods in both summers, we did not find occasions of CH4 production following the excess irrigation, with one exception in July 2019 with emission of 18 200 µg CH4 m−2 h−1. Otherwise, the soil was always a CH4 sink (median CH4 uptake rate of 260–290 and 150–170 µg CH4 m−2 h−1, in control and irrigated plots, respectively). The median soil CH4 uptake rates at the irrigated plot were 88 % and 50 % lower than at the control plot in 2018 and 2019, respectively. Potential CH4 production rates were highest in the organic layer (0.2–0.6 nmol CH4 g−1 d−1), but some production was also observed in the leaching layer, whereas in other soil layers, the rates were negligible. Potential CH4 oxidation rates varied mainly within 10–40 nmol CH4 g−1 d−1, except in deep soil and the organic layer in 2019, where potential oxidation rates were almost zero. The laboratory experiment revealed that high soil moisture alone does not turn upland forest soil into a CH4 source. However, a simple C source, e.g., substrates coming from root exudates with high moisture, switched the soil into a CH4 source. Our unique study provides new insights into the processes and controlling factors on CH4 production and oxidation, and the resulting net efflux that should be incorporated in process models describing global CH4 cycling.

In the global methane budget, the largest natural source is attributed to wetlands, which encompass all ecosystems composed of waterlogged or inundated ground, capable of methane production. Among them, northern peatlands that store large amounts of soil organic carbon have been functioning, since the end of the last glaciation period, as long-term sources of methane (CH4) and are one of the most significant methane sources among wetlands. To reduce uncertainty of quantifying methane flux in the global methane budget, it is of significance to understand the underlying processes for methane production and fluxes in northern peatlands. A methane model that features methane production and transport by plants, ebullition process and diffusion in soil, oxidation to CO2, and CH4 fluxes to the atmosphere has been embedded in the ORCHIDEE-PEAT land surface model that includes an explicit representation of northern peatlands. ORCHIDEE-PCH4 was calibrated and evaluated on 14 peatland sites distributed on both the Eurasian and American continents in the northern boreal and temperate regions. Data assimilation approaches were employed to optimized parameters at each site and at all sites simultaneously. Results show that methanogenesis is sensitive to temperature and substrate availability over the top 75 cm of soil depth. Methane emissions estimated using single site optimization (SSO) of model parameters are underestimated by 9 g CH4 m−2 yr−1 on average (i.e., 50 % higher than the site average of yearly methane emissions). While using the multi-site optimization (MSO), methane emissions are overestimated by 5 g CH4 m−2 yr−1 on average across all investigated sites (i.e., 37 % lower than the site average of yearly methane emissions).

Soil moisture plays a key role in soil nutrient and carbon cycling; plant productivity; and energy, water, and greenhouse gas exchanges between the land and the atmosphere. The knowledge on drivers of spatiotemporal soil moisture dynamics in subarctic landscapes is limited. In this study, we used the Spatial Forest Hydrology (SpaFHy) model, in situ soil moisture data, and Sentinel-1 synthetic aperture radar (SAR)-based soil moisture estimates to explore spatiotemporal controls of soil moisture in a subarctic headwater catchment in northwestern Finland. The role of groundwater dynamics and lateral flow in soil moisture was studied through three groundwater model conceptualizations: (i) omission of groundwater storage and lateral flow, (ii) conceptual TOPMODEL approach based on topographic wetness index, and (iii) explicit 2D lateral groundwater flow. The model simulations were compared against continuous point soil moisture measurements, distributed manual measurements, and novel SAR-based soil moisture estimates available at high spatial and temporal resolutions. Based on model scenarios and model–data comparisons, we assessed when and where the lateral groundwater flow shapes shallow soil moisture and under which conditions soil moisture variability is driven more by local ecohydrology, i.e., the balance of infiltration, drainage, and evapotranspiration. The choice of groundwater flow model was shown to have a strong impact on modeled soil moisture dynamics within the catchment. All model conceptualizations captured the observed soil moisture dynamics in the upland forests, but accounting for the lateral groundwater flow was necessary to reproduce the saturated conditions common in the peatlands and occasionally in lowland forest grid cells. We further highlight the potential of integrating multi-scale observations with land surface and hydrological models. The results have implications for ecohydrological and biogeochemical processes, as well as for modeling hydrology and Earth system feedbacks in subarctic and boreal environments. [ABSTRACT FROM AUTHOR]

Better monitoring, reporting, and verification (MRV) of the amount, additionality, and persistence of the sequestered soil carbon is needed to understand the best carbon farming practices for different soils and climate conditions, as well as their actual climate benefits or cost efficiency in mitigating greenhouse gas emissions. This paper presents our Field Observatory Network (FiON) of researchers, farmers, companies, and other stakeholders developing carbon farming practices. FiON has established a unified methodology towards monitoring and forecasting agricultural carbon sequestration by combining offline and near-real-time field measurements, weather data, satellite imagery, modeling, and computing networks. FiON's first phase consists of two intensive research sites and 20 voluntary pilot farms testing carbon farming practices in Finland. To disseminate the data, FiON built a web-based dashboard called the Field Observatory (v1.0, https://www.fieldobservatory.org/, last access: 3 February 2022). The Field Observatory is designed as an online service for near-real-time model–data synthesis, forecasting, and decision support for the farmers who are able to monitor the effects of carbon farming practices. The most advanced features of the Field Observatory are visible on the Qvidja site, which acts as a prototype for the most recent implementations. Overall, FiON aims to create new knowledge on agricultural soil carbon sequestration and effects of carbon farming practices as well as provide an MRV tool for decision support.

Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO2) fluxes in terrestrial ecosystems across the rapidly warming Arctic–boreal zone (ABZ) have provided valuable information but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic–boreal CO2 fluxes (ABCflux) that aggregates in situ measurements of terrestrial net ecosystem CO2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures), and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19 % of the monthly observations), snow diffusion (3 %) and eddy covariance (78 %) techniques. The largest number of observations were collected during the climatological summer (June–August; 32 %), and fewer observations were available for autumn (September–October; 25 %), winter (December–February; 18 %), and spring (March–May; 25 %). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO2 fluxes and to better estimate the terrestrial ABZ CO2 budget. ABCflux is openly and freely available online (Virkkala et al., 2021b, https://doi.org/10.3334/ORNLDAAC/1934).

Research infrastructures play a key role in launching a new generation of integrated long-Term, geographically distributed observation programmes designed to monitor climate change, better understand its impacts on global ecosystems, and evaluate possible mitigation and adaptation strategies. The pan-European Integrated Carbon Observation System combines carbon and greenhouse gas (GHG; CO 2 , CH 4 , N 2 O, H 2 O) observations within the atmosphere, terrestrial ecosystems and oceans. High-precision measurements are obtained using standardised methodologies, are centrally processed and openly available in a traceable and verifiable fashion in combination with detailed metadata. The Integrated Carbon Observation System ecosystem station network aims to sample climate and land-cover variability across Europe. In addition to GHG flux measurements, a large set of complementary data (including management practices, vegetation and soil characteristics) is collected to support the interpretation, spatial upscaling and modelling of observed ecosystem carbon and GHG dynamics. The applied sampling design was developed and formulated in protocols by the scientific community, representing a trade-off between an ideal dataset and practical feasibility. The use of open-Access, high-quality and multi-level data products by different user communities is crucial for the Integrated Carbon Observation System in order to achieve its scientific potential and societal value.

Organic soil amendments such as manure, biochar and compost are among the most efficient and widely used methods to increase soil carbon sequestration in agricultural soils. Even though their benefits are well known, many wood-derived materials are not yet utilized in Nordic agriculture due to a lack of incentives and knowledge of their effects in the local climate. We studied greenhouse gas exchange, plant growth and soil properties of a clay soil cultivated with oat in southern Finland in an extremely dry year. Two years earlier, the field was treated with three ligneous soil amendments—lime-stabilized fiber from the pulp industry, willow biochar and spruce biochar—which we compared against fertilized and non-fertilized controls. We found that the soil amendments increased porosity and the mean soil water holding capacity, which was most noticeable in plots amended with spruce biochar. There was a trend indicating that the mean yield and overall biomass production were larger in plots with soil amendments; however, the difference to unamended control was seldom significant due to the high variance among replicates. Manual chamber measurements revealed that carbon dioxide and methane exchange rates were reduced most probably by the exceptionally hot and dry weather conditions, but no differences could be found between the amended and unamended treatments. The nitrous oxide emissions were significantly smaller from the vegetated soil amended with willow biochar compared with the unamended control. Emissions from non-vegetated soil, representing heterotrophic respiration, were similar but without significant differences between treatments. Overall, the studied soil amendments indicated positive climatic impact two years after their application, but further research is needed to conclusively characterize the specific effects of organic soil amendments on processes affecting greenhouse gas exchange and plant growth.

Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20∘ S to 20∘ N) the spring onset of elevated CH4 emissions starts 3 d earlier, and the CH4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

A significant proportion of the global carbon emissions to the atmosphere originate from agriculture. Therefore, continuous long-term monitoring of CO2 fluxes is essential to understand the carbon dynamics and balances of different agricultural sites. Here we present results from a new eddy covariance flux measurement site located in southern Finland. We measured CO2 and H2O fluxes at this agricultural grassland site for 2 years, from May 2018 to May 2020. In particular the first summer experienced prolonged dry periods, which affected the CO2 fluxes, and substantially larger fluxes were observed in the second summer. During the dry summer, leaf area index (LAI) was notably lower than in the second summer. Water use efficiency increased with LAI in a similar manner in both years, but photosynthetic capacity per leaf area was lower during the dry summer. The annual carbon balance was calculated based on the CO2 fluxes and management measures, which included input of carbon as organic fertilizers and output as yield. The carbon balance of the field was −57 ± 10 and −86 ± 12 g C m−2 yr−1 in the first and second study years, respectively.

Water flows in peatland margins is an under-researched topic. This study examines recharge from a peatland to an esker aquifer in an aapa mire complex of northern Finland. Our objective was to study how the aapa mire margin is hydrogeologically connected to the riverside aquifer and spatial and temporal variations in the recharge of peatland water to groundwater (GW). Following geophysical studies and monitoring of the saturated zone, a GW model (MODFLOW) was used in combination with stable isotopes to quantify GW flow volumes and directions. Peatland water recharge to the sandy aquifer indicated a strong connection at the peatland–aquifer boundary. Recharge volumes from peatland to esker were high and rather constant (873 m3 d−1) and dominated esker recharge at the study site. The peat water recharging the esker boundary was rich in dissolved organic carbon (DOC). Stable isotope studies on water (δ18O, δ2H, and d-excess) from GW wells verified the recharge of DOC-rich water from peatlands to mineral soil esker. Biogeochemical analysis revealed changes from DOC to dissolved inorganic carbon in the flow pathway from peatland margin to the river Kitinen. This study highlights the importance of careful investigation of aapa mire margin areas and their potential role in regional GW recharge patterns. HIGHLIGHTS Peatland water recharge to aquifer showed connection at the peatland–aquifer boundary.; Analysis revealed changes from dissolved organic carbon to dissolved inorganic carbon in groundwater flow pathway from peatland.; Connection between an aapa mire margin and riverside esker was documented.;

Wetland methane emissions contribute to global warming, and are oversimplified in climate models. Here the authors use eddy covariance measurements from 48 global sites to demonstrate seasonal hysteresis in methane-temperature relationships and suggest the importance of microbial processes.

The patterned microtopography of subarctic mires generates a variety of environmental conditions, and carbon dioxide (CO2) and methane (CH4) dynamics vary spatially among different plant community types (PCTs). We studied the CO2 and CH4 exchange between a subarctic fen and the atmosphere at Kaamanen in northern Finland based on flux chamber and eddy covariance measurements in 2017–2018. We observed strong spatial variation in carbon dynamics between the four main PCTs studied, which were largely controlled by water table level and differences in vegetation composition. The ecosystem respiration (ER) and gross primary productivity (GPP) increased gradually from the wettest PCT to the drier ones, and both ER and GPP were larger for all PCTs during the warmer and drier growing season 2018. We estimated that in 2017 the growing season CO2 balances of the PCTs ranged from −20 g C m−2 (Trichophorum tussock PCT) to 64 g C m−2 (string margin PCT), while in 2018 all PCTs were small CO2 sources (10–22 g C m−2). We observed small growing season CH4 emissions (< 1 g C m−2) from the driest PCT, while the other three PCTs had significantly larger emissions (mean 7.9, range 5.6–10.1 g C m−2) during the two growing seasons. Compared to the annual CO2 balance (−8.5 ± 4.0 g C m−2) of the fen in 2017, in 2018 the annual balance (−5.6 ± 3.7 g C m−2) was affected by an earlier onset of photosynthesis in spring, which increased the CO2 sink, and a drought event during summer, which decreased the sink. The CH4 emissions were also affected by the drought. The annual CH4 balance of the fen was 7.3 ± 0.2 g C m−2 in 2017 and 6.2 ± 0.1 g C m−2 in 2018. Thus, the carbon balance of the fen was close to zero in both years. The PCTs that were adapted to drier conditions provided ecosystem-level resilience to carbon loss due to water level drawdown.

Northern peatlands are projected to be crucial in future atmospheric methane (CH4) budgets and have a positive feedback on global warming. Fens receive nutrients from catchments via inflowing water and are more sensitive than bogs to variations in their ecohydrology. Yet, due to a lack of data detailing the impacts of moving water on microhabitats and CH4 fluxes in fens, large uncertainties remain with respect to predicting CH4 emissions from these sites under climate changes. We measured CH4 fluxes with manual chambers over three growing seasons (2017–2019) at a northern boreal fen. To address the spatial variation at the site where a stream flows through the long and narrow valley fen, we established sample plots at varying distances from the stream. To link the variations in CH4 emissions to environmental controls, we quantified water levels, peat temperature, dissolved oxygen concentration, vegetation composition, and leaf area index in combination with flux measurements during the growing season in 2019. We found that due to the flowing water, there was a higher water level, cooler peat temperatures, and more oxygen in the peat close to the stream, which also had the highest total leaf area and gross primary production (GPP) values but the lowest CH4 emissions. CH4 emissions were highest at an intermediate distance from the stream where the oxygen concentration in the surface peat was low but GPP was still high. Further from the stream, the conditions were drier and produced low CH4 emissions. Our results emphasize the key role of ecohydrology in CH4 dynamics in fens and, for the first time, show how a stream controls CH4 emissions in a flow-through fen. As valley fens are common peatland ecosystems from the Arctic to the temperate zones, future projections of global CH4 budgets need to take flowing water features into account.

We have studied biogenic volatile organic compound (BVOC) emissions and their ambient concentrations at a sub-Arctic wetland (Lompolojänkkä, Finland), which is an open, nutrient-rich sedge fen and a part of the Pallas-Sodankylä Global Atmosphere Watch (GAW) station. Measurements were conducted during the growing season in 2018 using an in situ thermal-desorption–gas-chromatograph–mass-spectrometer (TD-GC-MS). Earlier studies have shown that isoprene is emitted from boreal wetlands, and it also turned out to be the most abundant compound in the current study. Monoterpene (MT) emissions were generally less than 10 % of the isoprene emissions (mean isoprene emission over the growing season, 44 µg m−2 h−1), but sesquiterpene (SQT) emissions were higher than MT emissions all the time. The main MTs emitted were α-pinene, 1,8-cineol, myrcene, limonene and 3Δ-carene. Of SQTs cadinene, β-cadinene and α-farnesene had the major contribution. During early growing season the SQT∕MT emission rate ratio was ∼10, but it became smaller as summer proceeded, being only ∼3 in July. Isoprene, MT and SQT emissions were exponentially dependent on temperature (correlation coefficients (R2) 0.75, 0.66 and 0.52, respectively). Isoprene emission rates were also found to be exponentially correlated with the gross primary production of CO2 (R2=0.85 in July). Even with the higher emissions from the wetland, ambient air concentrations of isoprene were on average > 100, > 10 and > 6 times lower than MT concentrations in May, June and July, respectively. This indicates that wetland was not the only source affecting atmospheric concentrations at the site, but surrounding coniferous forests, which are high MT emitters, contribute as well. Daily mean MT concentrations had high negative exponential correlation (R2=0.96) with daily mean ozone concentrations indicating that vegetation emissions can be a significant chemical sink of ozone in this sub-Arctic area.

The impact of atmospheric reactive nitrogen (Nr) deposition on carbon (C) sequestration in soils and biomass of unfertilized, natural, semi-natural and forest ecosystems has been much debated. Many previous results of this dC∕dN response were based on changes in carbon stocks from periodical soil and ecosystem inventories, associated with estimates of Nr deposition obtained from large-scale chemical transport models. This study and a companion paper (Flechard et al., 2020) strive to reduce uncertainties of N effects on C sequestration by linking multi-annual gross and net ecosystem productivity estimates from 40 eddy covariance flux towers across Europe to local measurement-based estimates of dry and wet Nr deposition from a dedicated collocated monitoring network. To identify possible ecological drivers and processes affecting the interplay between C and Nr inputs and losses, these data were also combined with in situ flux measurements of NO, N2O and CH4 fluxes; soil NO3- leaching sampling; and results of soil incubation experiments for N and greenhouse gas (GHG) emissions, as well as surveys of available data from online databases and from the literature, together with forest ecosystem (BASFOR) modelling. Multi-year averages of net ecosystem productivity (NEP) in forests ranged from −70 to 826 g C m−2 yr−1 at total wet + dry inorganic Nr deposition rates (Ndep) of 0.3 to 4.3 g N m−2 yr−1 and from −4 to 361 g C m−2 yr−1 at Ndep rates of 0.1 to 3.1 g N m−2 yr−1 in short semi-natural vegetation (moorlands, wetlands and unfertilized extensively managed grasslands). The GHG budgets of the forests were strongly dominated by CO2 exchange, while CH4 and N2O exchange comprised a larger proportion of the GHG balance in short semi-natural vegetation. Uncertainties in elemental budgets were much larger for nitrogen than carbon, especially at sites with elevated Ndep where Nr leaching losses were also very large, and compounded by the lack of reliable data on organic nitrogen and N2 losses by denitrification. Nitrogen losses in the form of NO, N2O and especially NO3- were on average 27 % (range 6 %–54 %) of Ndep at sites with Ndep −2 yr−1 versus 65 % (range 35 %–85 %) for Ndep > 3 g N m−2 yr−1. Such large levels of Nr loss likely indicate that different stages of N saturation occurred at a number of sites. The joint analysis of the C and N budgets provided further hints that N saturation could be detected in altered patterns of forest growth. Net ecosystem productivity increased with Nr deposition up to 2–2.5 g N m−2 yr−1, with large scatter associated with a wide range in carbon sequestration efficiency (CSE, defined as the NEP ∕ GPP ratio). At elevated Ndep levels (> 2.5 g N m−2 yr−1), where inorganic Nr losses were also increasingly large, NEP levelled off and then decreased. The apparent increase in NEP at low to intermediate Ndep levels was partly the result of geographical cross-correlations between Ndep and climate, indicating that the actual mean dC∕dN response at individual sites was significantly lower than would be suggested by a simple, straightforward regression of NEP vs. Ndep.

The effects of atmospheric nitrogen deposition (Ndep) on carbon (C) sequestration in forests have often been assessed by relating differences in productivity to spatial variations of Ndep across a large geographic domain. These correlations generally suffer from covariation of other confounding variables related to climate and other growth-limiting factors, as well as large uncertainties in total (dry + wet) reactive nitrogen (Nr) deposition. We propose a methodology for untangling the effects of Ndep from those of meteorological variables, soil water retention capacity and stand age, using a mechanistic forest growth model in combination with eddy covariance CO2 exchange fluxes from a Europe-wide network of 22 forest flux towers. Total Nr deposition rates were estimated from local measurements as far as possible. The forest data were compared with data from natural or semi-natural, non-woody vegetation sites. The response of forest net ecosystem productivity to nitrogen deposition (dNEP ∕ dNdep) was estimated after accounting for the effects on gross primary productivity (GPP) of the co-correlates by means of a meta-modelling standardization procedure, which resulted in a reduction by a factor of about 2 of the uncorrected, apparent dGPP ∕ dNdep value. This model-enhanced analysis of the C and Ndep flux observations at the scale of the European network suggests a mean overall dNEP ∕ dNdep response of forest lifetime C sequestration to Ndep of the order of 40–50 g C per g N, which is slightly larger but not significantly different from the range of estimates published in the most recent reviews. Importantly, patterns of gross primary and net ecosystem productivity versus Ndep were non-linear, with no further growth responses at high Ndep levels (Ndep > 2.5–3 g N m−2 yr−1) but accompanied by increasingly large ecosystem N losses by leaching and gaseous emissions. The reduced increase in productivity per unit N deposited at high Ndep levels implies that the forecast increased Nr emissions and increased Ndep levels in large areas of Asia may not positively impact the continent's forest CO2 sink. The large level of unexplained variability in observed carbon sequestration efficiency (CSE) across sites further adds to the uncertainty in the dC∕dN response.

Accounting for temporal changes in carbon dioxide (CO _2 ) effluxes from freshwaters remains a challenge for global and regional carbon budgets. Here, we synthesize 171 site-months of flux measurements of CO _2 based on the eddy covariance method from 13 lakes and reservoirs in the Northern Hemisphere, and quantify dynamics at multiple temporal scales. We found pronounced sub-annual variability in CO _2 flux at all sites. By accounting for diel variation, only 11% of site-months were net daily sinks of CO _2 . Annual CO _2 emissions had an average of 25% (range 3%–58%) interannual variation. Similar to studies on streams, nighttime emissions regularly exceeded daytime emissions. Biophysical regulations of CO _2 flux variability were delineated through mutual information analysis. Sample analysis of CO _2 fluxes indicate the importance of continuous measurements. Better characterization of short- and long-term variability is necessary to understand and improve detection of temporal changes of CO _2 fluxes in response to natural and anthropogenic drivers. Our results indicate that existing global lake carbon budgets relying primarily on daytime measurements yield underestimates of net emissions.