A drained nutrient-poor peatland forest in boreal Sweden constitutes a net carbon sink after integrating terrestrial and aquatic fluxes.

Northern peatlands provide a globally important carbon (C) store. Since the beginning of the 20th century, however, large areas of natural peatlands have been drained for biomass production across Fennoscandia. Today, drained peatland forests constitute a common feature of the managed boreal landscape, yet their ecosystem C balance and associated climate impact are not well understood, particularly within the nutrient-poor boreal region. In this study, we estimated the net ecosystem carbon balance (NECB) from a nutrient-poor drained peatland forest and an adjacent natural mire in northern Sweden by integrating terrestrial carbon dioxide (CO2 ) and methane (CH4 ) fluxes with aquatic losses of dissolved organic C (DOC) and inorganic C based on eddy covariance and stream discharge measurements, respectively, over two hydrological years. Since the forest included a dense spruce-birch area and a sparse pine area, we were able to further evaluate the effect of contrasting forest structure on the NECB and component fluxes. We found that the drained peatland forest was a net C sink with a 2-year mean NECB of -115 ± 5 g C m-2 year-1 while the adjacent mire was close to C neutral with 14.6 ± 1.7 g C m-2 year-1 . The NECB of the drained peatland forest was dominated by the net CO2 exchange (net ecosystem exchange [NEE]), whereas NEE and DOC export fluxes contributed equally to the mire NECB. We further found that the C sink strength in the sparse pine forest area (-153 ± 8 g C m-2 year-1 ) was about 1.5 times as high as in the dense spruce-birch forest area (-95 ± 8 g C m-2 year-1 ) due to enhanced C uptake by ground vegetation and lower DOC export. Our study suggests that historically drained peatland forests in nutrient-poor boreal regions may provide a significant net ecosystem C sink and associated climate benefits.

Mercury deposition and redox transformation processes in peatland constrained by mercury stable isotopes.

Peatland vegetation 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 remains, 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 redox 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 above the water table, dark abiotic oxidation decreases peat soil gaseous Hg0 concentrations. Below the water table, supersaturation of gaseous Hg is likely created more by direct photoreduction of rainfall rather than by reduction and release of Hg from the peat soil. Identification and quantification of these light-driven and dark redox processes advance our understanding of the fate of Hg in peatlands, including the potential for mobilization and methylation of HgII.

Catchment characteristics control boreal mire nutrient regime and vegetation patterns over ~5000 years of landscape development.

Vegetation holds the key to many properties that make natural mires unique, such as surface microtopography, high biodiversity values, effective carbon sequestration and regulation of water and nutrient fluxes across the landscape. Despite this, landscape controls behind mire vegetation patterns have previously been poorly described at large spatial scales, which limits the understanding of basic drivers underpinning mire ecosystem services. We studied catchment controls on mire nutrient regimes and vegetation patterns using a geographically constrained natural mire chronosequence along the isostatically rising coastline in Northern Sweden. By comparing mires of different ages, we can partition vegetation patterns caused by long-term mire succession (<5000 years) and present-day vegetation responses to catchment eco-hydrological settings. We used the remote sensing based normalized difference vegetation index (NDVI) to describe mire vegetation and combined peat physicochemical measures with catchment properties to identify the most important factors that determine mire NDVI. We found strong evidence that mire NDVI depends on nutrient inputs from the catchment area or underlying mineral soil, especially concerning phosphorus and potassium concentrations. Steep mire and catchment slopes, dry conditions and large catchment areas relative to mire areas were associated with higher NDVI. We also found long-term successional patterns, with lower NDVI in older mires. Importantly, the NDVI should be used to describe mire vegetation patterns in open mires if the focus is on surface vegetation, since the canopy cover in tree-covered mires completely dominated the NDVI signal. With our study approach, we can quantitatively describe the connection between landscape properties and mire nutrient regime. Our results confirm that mire vegetation responds to the upslope catchment area, but importantly, also suggest that mire and catchment aging can override the role of catchment influence. This effect was clear across mires of all ages, but was strongest in younger mires.

On the uncertainty in estimates of the carbon balance recovery time after forest clear-cutting.

An essential metric for describing carbon dynamics in managed forest landscapes is the recovery time of the carbon balance after clear-cutting. Here, we demonstrate how the age-dependent mathematical trajectory is affected by both the selected model and data availability, leading to considerable uncertainty in the modelling of the net ecosystem production (NEP) over stand age. We further show that the initial carbon loss estimates associated with the timing of the source-sink transition (SST) are significant, but may have a limited effect on the total carbon sequestration at the end of the standard (RP, 120 years) or optimal (OCS) rotation periods.

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

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

Microbial communities mediating net methylmercury formation along a trophic gradient in a peatland chronosequence.

Peatlands are generally important sources of methylmercury (MeHg) to adjacent aquatic ecosystems, increasing the risk of human and wildlife exposure to this highly toxic compound. While microorganisms play important roles in mercury (Hg) geochemical cycles where they directly and indirectly affect MeHg formation in peatlands, potential linkages between net MeHg formation and microbial communities involving these microorganisms remain unclear. To address this gap, microbial community composition and specific marker gene transcripts were investigated along a trophic gradient in a geographically constrained peatland chronosequence. Our results showed a clear spatial pattern in microbial community composition along the gradient that was highly driven by peat soil properties and significantly associated with net MeHg formation as approximated by MeHg concentration and %MeHg of total Hg concentration. Known fermentative, syntrophic, methanogenic and iron-reducing metabolic guilds had the strong positive correlations to net MeHg formation, while methanotrophic and methylotrophic microorganisms were negatively correlated. Our results indicated that sulfate reducers did not have a key role in net MeHg formation. Microbial activity as interpreted from 16S rRNA sequences was significantly correlated with MeHg and %MeHg. Our findings shed new light on the role of microbial community in net MeHg formation of peatlands that undergo ontogenetic change.

Landscape-variability of the carbon balance across managed boreal forests.

Boreal forests are important global carbon (C) sinks and, therefore, considered as a key element in climate change mitigation policies. However, their actual C sink strength is uncertain and under debate, particularly for the actively managed forests in the boreal regions of Fennoscandia. In this study, we use an extensive set of biometric- and chamber-based C flux data collected in 50 forest stands (ranging from 5 to 211 years) over 3 years (2016-2018) with the aim to explore the variations of the annual net ecosystem production (NEP; i.e., the ecosystem C balance) across a 68 km2 managed boreal forest landscape in northern Sweden. Our results demonstrate that net primary production rather than heterotrophic respiration regulated the spatio-temporal variations of NEP across the heterogeneous mosaic of the managed boreal forest landscape. We further find divergent successional patterns of NEP in our managed forests relative to naturally regenerating boreal forests, including (i) a fast recovery of the C sink function within the first decade after harvest due to the rapid establishment of a productive understory layer and (ii) a sustained C sink in old stands (131-211 years). We estimate that the rotation period for optimum C sequestration extends to 138 years, which over multiple rotations results in a long-term C sequestration rate of 86.5 t C ha-1 per rotation. Our study highlights the potential of forest management to maximize C sequestration of boreal forest landscapes and associate climate change mitigation effects by developing strategies that optimize tree biomass production rather than heterotrophic soil C emissions.

Organochemical Characterization of Peat Reveals Decomposition of Specific Hemicellulose Structures as the Main Cause of Organic Matter Loss in the Acrotelm.

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 1H-13C nuclear magnetic resonance (NMR) spectroscopy, and comparison to δ13C, δ15N, 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. δ13C, δ15N, 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.

Draft Metagenome Sequences of the Sphagnum (Peat Moss) Microbiome from Ambient and Warmed Environments across Europe.

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.

Geochemical and Dietary Drivers of Mercury Bioaccumulation in Estuarine Benthic Invertebrates.

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.

A novel belowground in-situ gas labeling approach: CH4 oxidation in deep peat using passive diffusion chambers and 13C excess.

In-vitro incubation of environmental samples is a common approach to estimate CH4 oxidation potential. Here we developed and verified an in-situ method utilizing passive diffusion chambers (PDC, silicone tubes) to deliver 13C-labeled CH4 into peat for the determination of the CH4 oxidation potential based on 13C excess of CO2. To target CH4 oxidation under semi-aerobic and anaerobic conditions, we installed 20 PDCs (30 ml volume) below the water table in profiles from 35-cm to 2-m depths of a peatland in north-eastern Sweden in July 2017 using a peat auger. 13C-labeled CH4 was injected into PDCs through tubing twice during 12 days (day 0 and 6) and samples were collected at days 1, 3, 6, 8 and 11. Background (non-labeled) δ13C of CO2 ranged from -7.3 (35 cm) to +5.7‰ (200 cm) with depth. These δ13C values rose to +110 and + 204‰ after the CH4 injection. The estimated CH4-derived C in CO2 was the lowest at the bottom of the profile (0.3 µmol L-1), whereas the maximum was at 100 cm (6.1 µmol L-1) at five days after the second labeling. This corresponded to 1.5-7.2% of the total CH4 pool to be oxidized, depending on depth. This novel approach with belowground in-situ 13C labeling of gases demonstrated the suitability of tracing the transformations of these gases in soil depth by PDCs and for the first time verified the in-situ occurrence of a deep-peat CH4 oxidation.

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

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

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

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

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

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

CO2 fertilization of Sphagnum peat mosses is modulated by water table level and other environmental factors.

Sphagnum mosses account for most accumulated dead organic matter in peatlands. Therefore, understanding their responses to increasing atmospheric CO2 is needed for estimating peatland C balances under climate change. A key process is photorespiration: a major determinant of net photosynthetic C assimilation that depends on the CO2 to O2 ratio. We used climate chambers to investigate photorespiratory responses of Sphagnum fuscum hummocks to recent increases in atmospheric CO2 (from 280 to 400 ppm) under different water table, temperature, and light intensity levels. We tested the photorespiratory variability using a novel method based on deuterium isotopomers (D6S /D6R ratio) of photosynthetic glucose. The effect of elevated CO2 on photorespiration was highly dependent on water table. At low water table (-20 cm), elevated CO2 suppressed photorespiration relative to C assimilation, thus substantially increasing the net primary production potential. In contrast, a high water table (~0 cm) favored photorespiration and abolished this CO2 effect. The response was further tested for Sphagnum majus lawns at typical water table levels (~0 and -7 cm), revealing no effect of CO2 under those conditions. Our results indicate that hummocks, which typically experience low water table levels, benefit from the 20th century's increase in atmospheric CO2 .

Chronic Atmospheric Reactive Nitrogen Deposition Suppresses Biological Nitrogen Fixation in Peatlands.

Biological nitrogen fixation (BNF) represents the natural pathway by which mosses meet their demands for bioavailable/reactive nitrogen (Nr) in peatlands. However, following intensification of nitrogen fertilizer and fossil fuel use, atmospheric Nr deposition has increased exposing peatlands to Nr loading often above the ecological threshold. As BNF is energy intensive, therefore, it is unclear whether BNF shuts down when Nr availability is no longer a rarity. We studied the response of BNF under a gradient of Nr deposition extending over decades in three peatlands in the U.K., and at a background deposition peatland in Sweden. Experimental nitrogen fertilization plots in the Swedish site were also evaluated for BNF activity. In situ BNF activity of peatlands receiving Nr deposition of 6, 17, and 27 kg N ha-1 yr-1 was not shut down but rather suppressed by 54, 69, and 74%, respectively, compared to the rates under background Nr deposition of ∼2 kg N ha-1 yr-1. These findings were corroborated by similar BNF suppression at the fertilization plots in Sweden. Therefore, contribution of BNF in peatlands exposed to chronic Nr deposition needs accounting when modeling peatland's nitrogen pools, given that nitrogen availability exerts a key control on the carbon capture of peatlands, globally.

Global CO2 fertilization of Sphagnum peat mosses via suppression of photorespiration during the twentieth century.

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.

Bimodal diel pattern in peatland ecosystem respiration rebuts uniform temperature response.

Accurate projections of climate change impacts on the vast carbon stores of northern peatlands require detailed knowledge of ecosystem respiration (ER) and its heterotrophic (Rh) and autotrophic (Ra) components. Currently, however, standard flux measurement techniques, i.e. eddy covariance and manual chambers, generate empirical ER data during only night- or daytime, respectively, which are extrapolated to the daily scale based on the paradigm that assumes a uniform diel temperature response. Here, using continuous autochamber measurements, we demonstrate a distinct bimodal pattern in diel peatland ER which contrasts the unimodal pattern inherent to the classical assumption. This feature results from divergent temperature dependencies of day- and nighttime ER due to varying contributions from Rh and Ra. We further find that disregarding these bimodal dynamics causes significant bias in ER estimates across multiple temporal scales. This calls for improved process-based understanding of ER to advance our ability to simulate peatland carbon cycle-climate feedbacks.

Formation and mobilization of methylmercury across natural and experimental sulfur deposition gradients.

We investigated the influence of sulfate (SO42-) deposition and concentrations on the net formation and solubility of methylmercury (MeHg) in peat soils. We used data from a natural sulfate deposition gradient running 300 km across southern Sweden to test the hypothesis posed by results from an experimental field study in northern Sweden: that increased loading of SO42- both increases net MeHg formation and redistributes methylmercury (MeHg) from the peat soil to its porewater. Sulfur concentrations in peat soils correlated positively with MeHg concentrations in peat porewater, along the deposition gradient similar to the response to added SO42- in the experimental field study. The combined results from the experimental field study and deposition gradient accentuate the multiple, distinct and interacting roles of SO42- deposition in the formation and redistribution of MeHg in the environment.

Opposing spatial trends in methylmercury and total mercury along a peatland chronosequence trophic gradient.

Peatlands are abundant elements of boreal landscapes where inorganic mercury (IHg) can be transformed into bioaccumulating and highly toxic methylmercury (MeHg). We studied fifteen peatlands divided into three age classes (young, intermediate and old) along a geographically constrained chronosequence to determine the role of biogeochemical factors and nutrient availability in controlling the formation of MeHg. In the 10 cm soil layer just below the average annual growing season water table, concentrations of MeHg and %MeHg (of total Hg) were higher in younger, more mesotrophic peatlands than in older, more oligotrophic peatlands. In contrast, total mercury (THg) concentrations were higher in the older peatlands. Partial least squares (PLS) analysis indicates that the net MeHg production was positively correlated to trophic demands of vegetation and an increased availability of potential electron acceptors and donors for Hg methylating microorganisms. An important question for further studies will be to elucidate why there is less THg in the younger peatlands compared to the older peatlands, even though the age of the superficial peat itself is similar for all sites. We hypothesize that ecosystem features which enhance microbial processes involved in Hg methylation also promote Hg reduction that makes previously deposited Hg more available for evasion back to the atmosphere.