Search

International audience; 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 (

We present 49 metagenome assemblies of the microbiome associated with Sphagnum (peat moss) collected from ambient, artificially warmed, and geothermally warmed conditions across Europe. These data will enable further research regarding the impact of climate change on plant-microbe symbiosis, ecology, and ecosystem functioning of northern peatland ecosystems.

Water table constitutes a master control of the general biogeochemistry in northern peatlands. The performance of peatland simulations in global ecosystem models is strongly hampered by the accuracy of the water table predictions. We examined the applicability of the Optical TRApezoid Model (OPTRAM) to monitor the temporal fluctuations in water table over 53 intact, restored, and drained northern peatlands in Finland, Estonia, Sweden, Canada, and the USA from 2018 through 2021. Various OPTRAM were computed based on Sentinel-2 data with the Google Earth Engine cloud platform. We found that (i) the choice of vegetation index utilised in OPTRAM does not significantly affect OPTRAM performance; (ii) the tree cover density is a significant factor controlling the sensitivity of OPTRAM to water table dynamics; (iii) the relationship between water table and OPTRAM often disappears for deep water tables. Based on an anomaly analysis, we further found that OPTRAM seems to be in particular suitable to monitor long-term (i.e., interannual) water table variability while the performance for short-term changes (e.g., response to individual rain events) was lower. Overall, our results support the application of OPTRAM to monitor water table dynamics in intact and restored northern peatlands with low tree cover density when the water table is shallow to moderately deep.

Ecosystem water use efficiency (WUE) is a key characteristic that describes the coupling of carbon and water exchange and can be used as an indicator of a forest's adaptability to varying climatic conditions. Mixed forests, characterized by the coexistence of two or more dominant tree species, may potentially exhibit higher productivity and greater resistance to extreme weather events due to possible niche differentiation among dominant species, leading to more efficient nutrient utilization. However, the increased productivity may also result in higher evapotranspiration demand, resulting in lower WUE compared to monospecific forests. In this study, we aim to assess the variation in WUE of mixed and monospecific boreal forests in response to different environmental factors using eddy-covariance measurements. The two study sites are represented by forest stands of similar age, growing under the same climatic conditions and located in close proximity (~10km distance) in Northern Sweden. The Rosinedalsheden site is a ~100-year-old monospecific pine (Pinus sylvestris) forest stand with sandy soils. The Svartberget site is a mixed ~110-year-old forest featuring pine (Pinus sylvestris, 61%), spruce (Picea abies, 34%), and birch (Betula sp., 5%) species, with soils dominated by till and sorted sediments. Our study spans a period of seven years (2014-2020) and covers a wide range of weather conditions, including the 2018 heatwave.

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

Peatland vegetation constantly takes up mercury (Hg) from the atmosphere, typically contributing to net production and export of neurotoxic Methyl-Hg to downstream ecosystems. Chemical reduction processes can slow down Methyl-Hg production by releasing Hg from peat back to the atmosphere. The extent of these processes remain, however, unclear. Here we present results from a comprehensive study covering concentrations and isotopic signatures of Hg in an open boreal peatland system to identify post-depositional Hg transformation processes. Isotope mass balances suggest photoreduction of HgII is the predominant process by which 30% of annually deposited Hg is emitted back to the atmosphere. Isotopic analyses indicate that dark abiotic oxidation decreases peat soil gaseous Hg0 concentrations above the water table, and that supersaturated gaseous Hg below the water table was likely created more by photoreduction of rainfall rather than release from peat soil. These findings constrain the possibilities for mobilizing previously deposited Hg. Identification and quantification of these processes further advance our understanding of the fate of Hg in peatlands and their export of Methyl-Hg to downstream ecosystems.

Little is known about the long-term expansion of mire ecosystems, despite their importance in the global carbon and hydrogeochemical cycles. It has been firmly established that mires do not expand linearly over time. Despite this, mires are often assumed to have expanded at a constant rate after initiation simply for lack of a better understanding. There has not yet been a serious attempt to determine the rate and drivers of mire expansion at the regional, or larger spatial scales. Here we make use of a natural chronosequence, spanning the Holocene, which is provided by the retreating coastline of Northern Sweden. By studying an isostatic rebound area we can infer mire expansion dynamics by looking at the portion of the landscape where mires become progressively scarce as the land becomes younger. Our results confirms that mires expanded non-linearly across the landscape and that their expansion is related to the availability of suitably wet areas, which, in our case, depends primarily on the hydro-edaphic properties of the landscape. Importantly, we found that mires occupied the wettest locations in the landscape within only one to two thousand years, while it took mires three to four thousand years to expand into slightly drier areas. Our results imply that the lateral expansion of mires, and thus peat accumulation is a non-linear process, occurring at different rates depending, above all else, on the wetness of the landscape.

Northern peatland stores a large amount of organic soil carbon and is considered to be one of the most significant CH4 sources among wetlands. The default wetland CH4 emission scheme in JULES (land surface model of the UK Earth System model) only takes into account the CH4 emissions from inundated areas in a simple way. However, it is known that the processes for peatland CH4 emission are complex. In this work, we coupled the process-based peatland CH4 emission model HIMMELI (HelsinkI Model of MEthane buiLd-up and emIssion for peatlands) with JULES (JULES-HIMMELI) by taking the HIMMELI input data from JULES simulations. Firstly, the soil temperature, water table depth (WTD) and soil carbon simulated by JULES, as well as the prescribed maximum leaf area index (LAI) in JULES were evaluated against available datasets at the studied northern wetland sites. Then, the simulated CH4 emissions from JULES and JULES-HIMMELI simulations were compared against the observed CH4 emissions at these sites. Moreover, sensitivities of CH4 emissions to the rate of anoxic soil respiration (anoxic Rs), surface soil temperature and WTD were investigated. Results show that JULES can well represent the magnitude and seasonality of surface (5–10 cm) and relatively deep (34–50 cm) soil temperatures, whereas the simulated WTD and soil carbon density profiles show large deviations from the site observations. The prescribed maximum LAI in JULES was within one standard deviation of the maximum LAIs derived from the Sentinel-2 satellite data for Siikaneva, Kopytkowo and Degerö sites, but lower for the other three sites. The simulated CH4 emissions by JULES have much smaller inter-annual variability than the observations. However, no specific simulation setup of the coupled model can lead to consistent improvements in the simulated CH4 emissions for all the sites. When using observed WTD or modified soil decomposition rate, there were only improvements in simulated CH4 fluxes at certain sites or years. Both simulated and observed CH4 emissions at sites strongly depend on the rate of anoxic Rs, which is the basis of CH4 emission estimates in HIMMELI. By excluding the effect from the rate of anoxic Rs on CH4 emissions, it is found that the Rs-log-normalized CH4 emissions (log normalization of the ratio of CH4 emission to anoxic Rs rate) show similar increasing trends with increased surface soil temperature from both observations and simulations, but different trends with raised WTD which may due to the uncertainty in simulated O2 concentration in HIMMELI. In general, we consider the JULES-HIMMELI model is more appropriate in simulating the wetland CH4 emissions than the default wetland CH4 emission scheme in JULES. Nevertheless, in order to improve the accuracy of simulated wetland CH4 emissions with the JULES-HIMMELI model, it is still necessary to better represent the peat soil carbon and hydrologic processes in JULES and the CH4 production and transportation processes in HIMMELI, such as plant transportation of gases, seasonality of parameters controlling oxidation and production, and adding microbial activities.

Peatlands store carbon in the form of dead organic residues. Climate change and human impact impose risks on the sustainability of the peatlands carbon balance due to increased peat decomposition. Here, we investigated molecular changes in the upper peat layers (0-40 cm), inferred from high-resolution vertical depth profiles, from a boreal peatland using two-dimensional H-1-C-13 nuclear magnetic resonance (NMR) spectroscopy, and comparison to delta C-13, delta N-15, and carbon and nitrogen content. Effects of hydrological conditions were investigated at respective sites: natural moist, drainage ditch, and natural dry. The molecular characterization revealed preferential degradation of specific side-chain linkages of xylan-type hemicelluloses within 0-14 cm at all sites, indicating organic matter losses up to 25%. In contrast, the xylan backbone, galactomannan-type hemicelluloses, and cellulose were more resistant to degradation and accumulated at the natural moist and drainage site. delta C-13, delta N-15, and carbon and nitrogen content did not correlate with specific hemicellulose structures but reflected changes in total carbohydrates. Our analysis provides novel insights into peat carbohydrate decomposition and indicates substantial organic matter losses in the acrotelm due to the degradation of specific hemicellulose structures. This suggests that variations in hemicellulose content and structure influence peat stability, which may have important implications with respect to climate change.

Sediments represent the main reservoir of mercury (Hg) in aquatic environments and may act as a source of Hg to aquatic food webs. Yet, accumulation routes of Hg from the sediment to benthic organisms are poorly constrained. We studied the bioaccumulation of inorganic and methylmercury (HgII and MeHg, respectively) from different geochemical pools of Hg into four groups of benthic invertebrates (amphipods, polychaetes, chironomids, and bivalves). The study was conducted using mesocosm experiments entailing the use of multiple isotopically enriched Hg tracers and simulation of estuarine systems with brackish water and sediment. We applied different loading regimes of nutrients and terrestrial organic matter and showed that the vertical localization and the chemical speciation of HgII and MeHg in the sediment, in combination with the diet composition of the invertebrates, consistently controlled the bioaccumulation of HgII and MeHg into the benthic organisms. Our results suggest a direct link between the concentration of MeHg in the pelagic planktonic food web and the concentration of MeHg in benthic amphipods and, to some extent, in bivalves. In contrast, the quantity of MeHg in benthic chironomids and polychaetes seems to be driven by MeHg accumulation via the benthic food web. Accounting for these geochemical and dietary drivers of Hg bioaccumulation in benthic invertebrates will be important to understand and predict Hg transfer between the benthic and the pelagic food web, under current and future environmental scenarios.

Mercury (Hg) in peatlands remains a problem of global interest. To mitigate the risks of this neurotoxin, accurate assessments of Hg in peat are needed. Treatment of peat that will be analysed for Hg is, however, not straightforward due to the volatile nature of Hg. The drying process is of particular concern since Hg evasion increases with the temperature. Samples are, therefore, often freeze-dried to limit Hg loss during the drying processes. A problem with freeze-drying is that cost and equipment resources can limit the number of samples analysed in large projects. To avoid this bottleneck, we tested if drying in a 60 °C-degree oven could be an acceptable alternative to freeze-drying. We both freeze-dried and oven-dried (60 °C) 203 replicate pairs of peat samples, and then examined the differences in total Hg concentration. The Hg concentration differed significantly between the two drying methods with a median Hg deficit in oven-dried samples of 4.2%. Whether a 4.2% deficit of Hg depends on one’s purpose. The lower median Hg concentration in oven-dried samples has to be weighed against the upside efficiently drying large sets of peat samples. By freeze-drying a subset of the samples, we fitted a function to correct for Hg loss during oven-drying ($$y=0.96x+0.08)$$ y = 0.96 x + 0.08 ) . By applying this correction, the freeze-drying bottleneck could oven-dry large-scale inventories of total Hg in peatlands with results equivalent to freeze-drying, but only have to freeze-dry a subset.

Ditch cleaning (DC) is increasingly applied to facilitate forest regeneration following clear-cutting in Fennoscandinavia. However, its impact on the ecosystem carbon and greenhouse gas (GHG) balances is poorly understood. We conducted chamber measurements to assess the initial DC effects on carbon dioxide (CO2) and methane (CH4) fluxes in a recent forest clear-cut on wet mineral soil in boreal Sweden. Measurements were conducted in two adjacent areas over two pre-treatments (2018/19) and two years (2020/21) after conducting DC in one area. We further assessed the spatial variation of fluxes at three distances (4, 20, 40 m) from ditches. We found that DC lowered the water table level by 12 ± 2 cm (mean ± standard error) and topsoil moisture by 0.12 ± 0.01 m3 m−3. DC had a limited initial effect on the net CO2 exchange and its component fluxes. CH4 emissions were low during the dry pre-treatment years but increased particularly in the control area during the wet years of 2020/21. Distance to ditch had no consistent effects on CO2 and CH4 fluxes. Model extrapolations suggest that annual carbon emissions decreased over the four years from 6.7 ± 1.4 to 1.6 ± 1.6 t-C ha−1 year−1, without treatment differences. Annual CH4 emissions contributed

To meet the sustainable development goals and enable protection of peatlands, there is a strong need to plan and align land-use management with the needs of the environment. The most critical tool to succeed in sustainable spatial planning is accurate and detailed maps. Here we present a novel approach to map mineral and peat soils based on a high-resolution digital soil moisture map. This soil moisture map was produced by combining LIDAR-derived terrain indices and machine learning to model soil moisture at 2 m spatial resolution across the Swedish landscape with high accuracy (Kappa = 0.69, MCC = 0.68). We used field data from about 20,000 sites across Sweden to train an extreme gradient boosting model to predict soil moisture. The predictor features included a suite of terrain indices generated from national LIDAR digital elevation model and other ancillary environmental features, including surficial geology, climate, land use information, allowing for adjustment of soil moisture maps to regional/local conditions. As soil moisture is an important control on peat formation, we investigated if this map can be used to improve the mapping of peatlands. In this study, we included a total of 5 479 soil pit data for organic layer thickness from the Swedish Forest Soil Inventory. Peat was defined as organic layer thickness > 50 cm. The data was split into a calibration dataset and a validation dataset using a randomized 50% split. An empirical relationship between the thickness of the organic layer and the continuous SLU soil moisture map (R2 = 0.66, p < 0.001) was used to generate both a categorical map (of mineral soil and peat) and continuous map (of organic layer thickness) to demonstrate how these two mapping approaches can be useful for different management objectives. The peat coverage on the new categorical map, the quaternary deposits map and topographical map was 17.3%, 14.1% and 13.5%, respectively. Map quality measures from the evaluation dataset showed that the newly developed peat map had higher recall and MCC (80.4, 0.73) than quaternary deposits map (68.5, 0.65) and topographical map (49.8, 0.61). The continuous map of the organic layer ranged 6-95 cm with an RMSE of 4 cm. Using Sweden as a test case, this study provides a guide to improved mapping of mineral and peat soils from Lidar data in other boreal forest regions for effective ecosystem management. The map of organic soils was developed to support the need for land use management optimization by incorporating landscape sensitivity and hydrological connectivity into a framework that promotes the protection of soil and water quality. The organic soil map can be used to address fundamental considerations, such as;guiding the restoration of drained wetlands; designing riparian protection zones to optimize the protection of water quality and biodiversity as the new map also include riparian peats.

Ditch cleaning (DC) is a well-established forestry practice across Fennoscandia to lower water table levels (WTL) and thereby facilitate the establishment of tree seedlings following clear-cutting. However, the implications from these activities for ecosystem-atmosphere greenhouse gas (GHG) exchanges are poorly understood at present. We assessed the initial DC effects on the GHG fluxes in a forest clear-cut on a drained fertile peatland in hemiboreal Sweden, by comparing chamber measurements of carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O) fluxes from soil and ditches in DC and uncleaned (UC) areas over the first two post-harvest years. We also evaluated spatial effects by comparing fluxes at 4 m and 40 m from ditches. We found that 2 years after DC, mean (+/- standard error) WTL of-65 +/- 2 cm was significantly lower in the DC area compared to-56 +/- 2 cm in the UC area. We further observed lower gross primary production and ecosystem respiration in the first year after DC which coincided with delayed development of herbaceous ground vegetation. We also found higher CH4 uptake but no difference in N2O fluxes after DC. Greater CH4 uptake occurred at 4 m compared to 40 m away from both cleaned and uncleaned ditches. Model extrapolation suggests that total annual GHG emissions in the second year were reduced from 49.4 +/- 17.0 t-CO2-eq-ha(-1) -year(-1) in the UC area to 27.8 +/- 10.3 t-CO2-eq-ha(-1) -year(-1) in the DC area. A flux partitioning approach suggested that this was likely caused by decreased heterotrophic respiration, possibly because of enhanced soil dryness following DC during the dry meteorological conditions. CH(4 )and N2O fluxes from clear-cut areas contributed < 2 % to the total (soil, ditches) GHG budget. Similarly the area -weighted contributions by CO2 and CH4 emissions from both cleaned and uncleaned ditches were < 2 %. Thus, our study highlights that DC may considerably alter the post-harvest GHG fluxes of drained peatland forests. However, long-term observations under various site conditions and forest rotation stages are warranted to better understand DC effects on the forest GHG balance.

Wetland CH4 emissions are among the most uncertain components of the global CH4 budget. The complex nature of wetland CH4 processes makes it challenging to identify causal relationships for improving our understanding and predictability of CH4 emissions. In this study, we used the flux measurements of CH4 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 CH4 emissions at sub-seasonal scale. We found that soil temperature is the dominant factor for CH4 emissions in all studied wetland types. Ecosystem respiration (CO2) 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 CH4 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 CH4 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 CH4 emissions within earth system land models.

To meet the sustainable development goals and enable sustainable management and protection of peatlands, there is a strong need for improving the mapping of peatlands. Here we present a novel approach to identify peat soils based on a high-resolution digital soil moisture map that was produced by combining airborne laser scanning-derived terrain indices and machine learning to model soil moisture at 2 m spatial resolution across the Swedish landscape. As soil moisture is a key factor in peat formation, we fitted an empirical relationship between the thickness of the organic layer (measured at 5479 soil plots across the country) and the continuous SLU (Swedish University of Agricultural Science) soil moisture map (R2= 0.66, p < 0.001). We generated categorical maps of peat occurrence using three different definitions of peat (30, 40, and 50 cm thickness of the organic layer) and a continuous map of organic layer thickness. The predicted peat maps had a higher overall quality (MCC = 0.69–0.73) compared to traditional Quaternary deposits maps (MCC = 0.65) and topographical maps (MCC = 0.61) and captured the peatlands with a recall of ca. 80 % compared to 50 %–70 % on the traditional maps. The predicted peat maps identified more peatland area than previous maps, and the areal coverage estimates fell within the same order as upscaling estimates from national field surveys. Our method was able to identify smaller peatlands resulting in more accurate maps of peat soils, which was not restricted to only large peatlands that can be visually detected from aerial imagery – the historical approach of mapping. We also provided a continuous map of the organic layer, which ranged 6–88 cm organic layer thickness, with an R2 of 0.67 and RMSE (root mean square error) of 19 cm. The continuous map exhibits a smooth transition of organic layers from mineral soil to peat soils and likely provides a more natural representation of the distribution of soils. The continuous map also provides an intuitive uncertainty estimate in the delineation of peat soils, critically useful for sustainable spatial planning, e.g., greenhouse gas or biodiversity inventories and landscape ecological research.

Natural peatlands contribute significantly to global carbon sequestration and storage of biomass, most of which derives from Sphagnum peat mosses. Atmospheric CO2 levels have increased dramatically during the twentieth century, from 280 to > 400 ppm, which has affected plant carbon dynamics. Net carbon assimilation is strongly reduced by photorespiration, a process that depends on the CO2 to O2 ratio. Here we investigate the response of the photorespiration to photosynthesis ratio in Sphagnum mosses to recent CO2 increases by comparing deuterium isotopomers of historical and contemporary Sphagnum tissues collected from 36 peat cores from five continents. Rising CO2 levels generally suppressed photorespiration relative to photosynthesis but the magnitude of suppression depended on the current water table depth. By estimating the changes in water table depth, temperature, and precipitation during the twentieth century, we excluded potential effects of these climate parameters on the observed isotopomer responses. Further, we showed that the photorespiration to photosynthesis ratio varied between Sphagnum subgenera, indicating differences in their photosynthetic capacity. The global suppression of photorespiration in Sphagnum suggests an increased net primary production potential in response to the ongoing rise in atmospheric CO2, in particular for mire structures with intermediate water table depths.

In-vitro incubation of environmental samples is a common approach to estimate CH

The boreal forest is an important global carbon (C) sink. Since low soil nitrogen (N) availability is commonly a key constraint on forest productivity, the prevalent view is that increased N input enhances its C sink-strength. This understanding however relies primarily on observations of increased aboveground tree biomass and soil C stock following N fertilization, whereas empirical data evaluating the effects on the whole ecosystem-scale C balance are lacking. Here we use a unique long-term experiment consisting of paired forest stands with eddy covariance measurements to explore the effect of ecosystem-scale N fertilization on the C balance of a managed boreal pine forest. We find that the annual C uptake (i.e. net ecosystem production, NEP) at the fertilized stand was 16 +/- 2% greater relative to the control stand by the end of the first decade of N addition. Subsequently, the ratio of NEP between the fertilized and control stand remained at a stable level during the following five years with an average NEP to N response of 7 & PLUSMN; 1 g C per g N. Our study reveals that this non-linear response of NEP to long-term N fertilization was the result of a cross-seasonal feedback between the N-induced increases in both growing-season C uptake and subsequent winter C emission. We further find that one decade of N addition altered the sensitivity of ecosystem C fluxes to key environmental drivers resulting in divergent responses to weather patterns. Thus, our study highlights the need to account for ecosystem-scale responses to perturbations to improve our understanding of nitrogen-carbon-climate feedbacks in boreal forests.

We measured CO2 and CH4 fluxes using chambers and eddy covariance (only CO2) from a moist moss tundra in Svalbard. The average net ecosystem exchange (NEE) during the summer (June–August) was −0.40 g C m−2 day−1 or −37 g C m−2 for the whole summer. Including spring and autumn periods the NEE was reduced to −6.8 g C m−2 and the annual NEE became positive, 24.7 gC m−2 due to the losses during the winter. The CH4 flux during the summer period showed a large spatial and temporal variability. The mean value of all 214 samples was 0.000511 ± 0.000315 µmol m−2s−1 which corresponds to a growing season estimate of 0.04 to 0.16 g CH4 m−2. We find that this moss tundra emits about 94–100 g CO2-equivalents m−2 yr−1 of which CH4 is responsible for 3.5–9.3 % using GWP100 of 27.9 respectively GWP20. Air temperature, soil moisture and greenness index contributed significantly to explain the variation in ecosystem respiration (Reco) while active layer depth, soil moisture and greenness index were the variables that best explained CH4 emissions. Estimate of temperature sensitivity of Reco and gross primary productivity showed that a modest increase in air temperature of 1 degree did not significantly change the NEE during the growing season but that the annual NEE would be even more positive adding another 8.5 g C m−2 to the atmosphere. We tentatively suggest that the warming of the Arctic that has already taken place is partly responsible for the fact that the moist moss tundra now is a source of CO2 to the atmosphere.

Sphagnum mosses account for most accumulated dead organic matter in peatlands. Therefore, understanding their responses to increasing atmospheric CO

The regional variability in tundra and boreal carbon dioxide (CO

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, including their seasonality, due to quasi-continuous and high temporal resolution of flux measurements, coincident measurements of carbon, water, and energy fluxes, 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 global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4- community-product/). 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 and because freshwater wetlands are a substantial source of total atmospheric CH4 emissions; and 3) 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 three days earlier, and the CH4 emission season lasts 4 days longer, for each degree C increase in mean annual air temperature. On average, the onset of increasing CH4 emissions lags soil warming by one 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). The FLUXNET-CH4 dataset provides an open-access resource for CH4 flux synthesis, has a range of applications, and is unique in that it includes coupled measurements of important CH4 drivers such as GPP and temperature. Although FLUXNET-CH4 could certainly be improved by adding more sites in tropical ecosystems and by increasing the number of site-years at existing sites, it is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4408468. Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/, and a complete list of the 79 individual site data DOIs is provided in Table 2 in the Data Availability section of this document., Earth System Science Data Discussions, ISSN:1866-3591

The forest-floor represents an important interface for various carbon dioxide (CO2) fluxes, however, our knowledge of their variability and drivers across a managed boreal forest landscape is limited. Here, we used a three-year (2016-2018) data set of biometric- and chamber-based flux measurements to investigate the net forest-floor CO2 exchange (NEff) and its component fluxes across 50 forest stands spanning different soil types, tree species, and age classes within a 68 km2 boreal catchment in Sweden. We found that the forest-floor acted as a net CO2 source with the 10th-90th percentile (used hereafter for describing reported variations) ranging from 149 to 399 g C m- 2 yr-1. Among the key landscape attributes, stand age strongly affected most NEff component fluxes, whereas tree species and soil type effects were weak and absent, respectively. Specifically, forest-floor net CO2 emissions increased with stand age due to declining understory gross and net primary production, ranging between 77-275 and 49-163 g C m- 2 yr- 1, respectively. Furthermore, we observed higher understory production rates in pine than in spruce stands. Across the 50 stands, the total forest-floor respiration ranged between 340 and 549 g C m- 2 yr-1 and its spatial variation was primarily regulated by its autotrophic components, i.e., understory and tree root respiration, which displayed divergent increasing and decreasing age-related trends, respectively. Furthermore, heterotrophic soil respiration remained within a relatively narrow range (154-290 g C m- 2 yr- 1), possibly owing to compensating gradients in forest-floor properties. We further identified tree biomass as the major driver of the landscape-scale variations of CO2 fluxes, likely attributable to modulating effects on forest-floor resource availability and growing conditions. This implies that tree growth responses to forest management and global change will be particularly important for regulating magnitudes and spatial variations of forest-floor CO2 fluxes in boreal forests.

We investigated the influence of sulfate (SO

Peatlands are globally important ecosystems where inorganic mercury is converted to bioaccumulating and highly toxic methylmercury, resulting in high risks of methylmercury exposure in adjacent aquatic ecosystems. Although biological mercury methylation has been known for decades, there is still a lack of knowledge about the organisms involved in mercury methylation and the drivers controlling their methylating capacity. In order to investigate the metabolisms responsible for mercury methylation and methylmercury degradation as well as the controls of both processes, we studied a chronosequence of boreal peatlands covering fundamentally different biogeochemical conditions. Potential mercury methylation rates decreased with peatland age, being up to 53 times higher in the youngest peatland compared to the oldest. Methylation in young mires was driven by sulfate reduction, while methanogenic and syntrophic metabolisms became more important in older systems. Demethylation rates were also highest in young wetlands, with a gradual shift from biotic to abiotic methylmercury degradation along the chronosequence. Our findings reveal how metabolic shifts drive mercury methylation and its ratio to demethylation as peatlands age.

The net ecosystem CO

Peatlands store substantial amount of carbon, are vulnerable to climate change. To predict the fate of carbon stored in peatlands, the complex interactions between water, peat and vegetations need more attention. This study describes a modified version of the ORCHIDEE land surface model for simulating the hydrology, surface energy and CO2 fluxes of peatlands on daily to annual time scales. The model, referred to as ORCHIDEE-PEAT, includes a separate soil tile in each 0.5° grid-cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation with a grid-cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model is evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (Vcmax) being optimized at each site to match the peak of growing season gross primary productivity (GPP), derived from direct EC measurements. Regarding short-term variations from day to day, the model performance was good for the variations in GPP (r2 = 0.76, Nash-Sutcliff modeling efficiency, MEF = 0.76), with lesser accuracy for latent heat fluxes (LE, r2 = 0.42, MEF = 0.14) and Net ecosystem CO2 exchange (NEE, r2 = 0.38, MEF = 0.26). Seasonal variations in GPP, NEE and energy fluxes on monthly scales showed moderate to high r2 values ranging from 0.57 to 0.86. For spatial across-sites gradients of annual mean GPP, NEE and LE, r2 of 0.93, 0.27, and 0.71, respectively, were achieved. The water table variations are not well predicted (r2 < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, when using the observed water table in the carbon module to define the fraction of oxic and anoxic decomposition instead of the modeled water table, ORCHIDEE-PEAT shows a small improvement in reproducing NEE. Moreover, we found a significant relationship between optimized Vcmax and the latitude (temperature), which can better reflect the spatial gradients of annual NEE than using an average Vcmax value. In a future version of ORCHIDEE-PEAT, the influences of water table on photosynthesis and depth-dependent influences of soil temperature on respiration may be included.

At high latitudes, winter climate change alters snow cover and, consequently, may cause a sustained change in soil frost dynamics. Altered winter soil conditions could influence the ecosystem exchange of carbon dioxide (CO

In high-latitude regions, carbon dioxide (CO

Terrestrial vegetation currently absorbs approximately a third of anthropogenic CO2 emissions, mitigating the rise of atmospheric CO2. However, terrestrial net primary production is highly sensitive to atmospheric CO2 levels and associated climatic changes. In C-3 plants, which dominate terrestrial vegetation, net photosynthesis depends on the ratio between photorespiration and gross photosynthesis. This metabolic flux ratio depends strongly on CO2 levels, but changes in this ratio over the past CO2 rise have not been analyzed experimentally. Combining CO2 manipulation experiments and deuterium NMR, we first establish that the intramolecular deuterium distribution (deuterium isotopomers) of photosynthetic C-3 glucose contains a signal of the photorespiration/photosynthesis ratio. By tracing this isotopomer signal in herbarium samples of natural C-3 vascular plant species, crops, and a Sphagnum moss species, we detect a consistent reduction in the photorespiration/photosynthesis ratio in response to the similar to 100-ppm CO2 increase between similar to 1900 and 2013. No difference was detected in the isotopomer trends between beet sugar samples covering the 20th century and CO2 manipulation experiments, suggesting that photosynthetic metabolism in sugar beet has not acclimated to increasing CO2 over >100 y. This provides observational evidence that the reduction of the photorespiration/photosynthesis ratio was ca. 25%. The Sphagnum results are consistent with the observed positive correlations between peat accumulation rates and photosynthetic rates over the Northern Hemisphere. Our results establish that isotopomers of plant archives contain metabolic information covering centuries. Our data provide direct quantitative information on the "CO2 fertilization" effect over decades, thus addressing a major uncertainty in Earth system models.

The effects of wood ash application on soil microbial processes were investigated in three drained forested peatlands, which differed in nutrient status and time since application. Measured variables included the concentrations of soil elements and phospholipid fatty acids (PLFAs), net nitrogen (N) mineralization, nitrification and denitrification enzyme activity, potential methane (CH(4)) oxidation, CH(4) production and microbial respiration kinetics. Wood ash application had a considerable influence on soil element concentrations. This mirrored a decrease in the majority of the microbial biomarkers by more than one-third in the two oligotrophic peatlands, although the microbial community composition was not altered. The decreases in PLFAs coincided with reduced net ammonification and net N mineralization. Other measured variables did not change systematically as a result of wood ash application. No significant changes in microbial biomass or processes were found in the mesotrophic peatland, possibly because too little time (1 year) had elapsed since the wood ash application. This study suggests that oligotrophic peatlands can be substantially affected by wood ash for a period of at least 4 years after application. However, within 25 years of the wood ash application, the microbial biomass seemed to have recovered or adapted to enhanced element concentrations in the soil.