Recent increases in annual, seasonal, and extreme methane fluxes driven by changes in climate and vegetation in boreal and temperate wetland ecosystems.

Climate warming is expected to increase global methane (CH4 ) emissions from wetland ecosystems. Although in situ eddy covariance (EC) measurements at ecosystem scales can potentially detect CH4 flux changes, most EC systems have only a few years of data collected, so temporal trends in CH4 remain uncertain. Here, we use established drivers to hindcast changes in CH4 fluxes (FCH4 ) since the early 1980s. We trained a machine learning (ML) model on CH4 flux measurements from 22 [methane-producing sites] in wetland, upland, and lake sites of the FLUXNET-CH4 database with at least two full years of measurements across temperate and boreal biomes. The gradient boosting decision tree ML model then hindcasted daily FCH4 over 1981-2018 using meteorological reanalysis data. We found that, mainly driven by rising temperature, half of the sites (n = 11) showed significant increases in annual, seasonal, and extreme FCH4 , with increases in FCH4 of ca. 10% or higher found in the fall from 1981-1989 to 2010-2018. The annual trends were driven by increases during summer and fall, particularly at high-CH4 -emitting fen sites dominated by aerenchymatous plants. We also found that the distribution of days of extremely high FCH4 (defined according to the 95th percentile of the daily FCH4 values over a reference period) have become more frequent during the last four decades and currently account for 10-40% of the total seasonal fluxes. The share of extreme FCH4 days in the total seasonal fluxes was greatest in winter for boreal/taiga sites and in spring for temperate sites, which highlights the increasing importance of the non-growing seasons in annual budgets. Our results shed light on the effects of climate warming on wetlands, which appears to be extending the CH4 emission seasons and boosting extreme emissions.

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.

Observational constraints reduce model spread but not uncertainty in global wetland methane emission estimates.

The recent rise in atmospheric methane (CH4 ) concentrations accelerates climate change and offsets mitigation efforts. Although wetlands are the largest natural CH4 source, estimates of global wetland CH4 emissions vary widely among approaches taken by bottom-up (BU) process-based biogeochemical models and top-down (TD) atmospheric inversion methods. Here, we integrate in situ measurements, multi-model ensembles, and a machine learning upscaling product into the International Land Model Benchmarking system to examine the relationship between wetland CH4 emission estimates and model performance. We find that using better-performing models identified by observational constraints reduces the spread of wetland CH4 emission estimates by 62% and 39% for BU- and TD-based approaches, respectively. However, global BU and TD CH4 emission estimate discrepancies increased by about 15% (from 31 to 36 TgCH4 year-1 ) when the top 20% models were used, although we consider this result moderately uncertain given the unevenly distributed global observations. Our analyses demonstrate that model performance ranking is subject to benchmark selection due to large inter-site variability, highlighting the importance of expanding coverage of benchmark sites to diverse environmental conditions. We encourage future development of wetland CH4 models to move beyond static benchmarking and focus on evaluating site-specific and ecosystem-specific variabilities inferred from observations.

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.

Where old meets new: An ecosystem study of methanogenesis in a reflooded agricultural peatland.

Reflooding formerly drained peatlands has been proposed as a means to reduce losses of organic matter and sequester soil carbon for climate change mitigation, but a renewal of high methane emissions has been reported for these ecosystems, offsetting mitigation potential. Our ability to interpret observed methane fluxes in reflooded peatlands and make predictions about future flux trends is limited due to a lack of detailed studies of methanogenic processes. In this study we investigate methanogenesis in a reflooded agricultural peatland in the Sacramento Delta, California. We use the stable-and radio-carbon isotopic signatures of wetland sediment methane, ecosystem-scale eddy covariance flux observations, and laboratory incubation experiments, to identify which carbon sources and methanogenic production pathways fuel methanogenesis and how these processes are affected by vegetation and seasonality. We found that the old peat contribution to annual methane emissions was large (~30%) compared to intact wetlands, indicating a biogeochemical legacy of drainage. However, fresh carbon and the acetoclastic pathway still accounted for the majority of methanogenesis throughout the year. Although temperature sensitivities for bulk peat methanogenesis were similar between open-water (Q10  = 2.1) and vegetated (Q10  = 2.3) soils, methane production from both fresh and old carbon sources showed pronounced seasonality in vegetated zones. We conclude that high methane emissions in restored wetlands constitute a biogeochemical trade-off with contemporary carbon uptake, given that methane efflux is fueled primarily by fresh carbon inputs.

Gap-filling approaches for eddy covariance methane fluxes: A comparison of three machine learning algorithms and a traditional method with principal component analysis.

Methane flux (FCH4 ) measurements using the eddy covariance technique have increased over the past decade. FCH4 measurements commonly include data gaps, as is the case with CO2 and energy fluxes. However, gap-filling FCH4 data are more challenging than other fluxes due to its unique characteristics including multidriver dependency, variabilities across multiple timescales, nonstationarity, spatial heterogeneity of flux footprints, and lagged influence of biophysical drivers. Some researchers have applied a marginal distribution sampling (MDS) algorithm, a standard gap-filling method for other fluxes, to FCH4 datasets, and others have applied artificial neural networks (ANN) to resolve the challenging characteristics of FCH4 . However, there is still no consensus regarding FCH4 gap-filling methods due to limited comparative research. We are not aware of the applications of machine learning (ML) algorithms beyond ANN to FCH4 datasets. Here, we compare the performance of MDS and three ML algorithms (ANN, random forest [RF], and support vector machine [SVM]) using multiple combinations of ancillary variables. In addition, we applied principal component analysis (PCA) as an input to the algorithms to address multidriver dependency of FCH4 and reduce the internal complexity of the algorithmic structures. We applied this approach to five benchmark FCH4 datasets from both natural and managed systems located in temperate and tropical wetlands and rice paddies. Results indicate that PCA improved the performance of MDS compared to traditional inputs. ML algorithms performed better when using all available biophysical variables compared to using PCA-derived inputs. Overall, RF was found to outperform other techniques for all sites. We found gap-filling uncertainty is much larger than measurement uncertainty in accumulated CH4 budget. Therefore, the approach used for FCH4 gap filling can have important implications for characterizing annual ecosystem-scale methane budgets, the accuracy of which is important for evaluating natural and managed systems and their interactions with global change processes.

Effects of seasonality, transport pathway, and spatial structure on greenhouse gas fluxes in a restored wetland.

Wetlands can influence global climate via greenhouse gas (GHG) exchange of carbon dioxide (CO2 ), methane (CH4 ), and nitrous oxide (N2 O). Few studies have quantified the full GHG budget of wetlands due to the high spatial and temporal variability of fluxes. We report annual open-water diffusion and ebullition fluxes of CO2 , CH4 , and N2 O from a restored emergent marsh ecosystem. We combined these data with concurrent eddy-covariance measurements of whole-ecosystem CO2 and CH4 exchange to estimate GHG fluxes and associated radiative forcing effects for the whole wetland, and separately for open-water and vegetated cover types. Annual open-water CO2 , CH4 , and N2 O emissions were 915 ± 95 g C-CO2  m-2  yr-1 , 2.9 ± 0.5 g C-CH4  m-2  yr-1 , and 62 ± 17 mg N-N2 O m-2  yr-1 , respectively. Diffusion dominated open-water GHG transport, accounting for >99% of CO2 and N2 O emissions, and ~71% of CH4 emissions. Seasonality was minor for CO2 emissions, whereas CH4 and N2 O fluxes displayed strong and asynchronous seasonal dynamics. Notably, the overall radiative forcing of open-water fluxes (3.5 ± 0.3 kg CO2 -eq m-2  yr-1 ) exceeded that of vegetated zones (1.4 ± 0.4 kg CO2 -eq m-2  yr-1 ) due to high ecosystem respiration. After scaling results to the entire wetland using object-based cover classification of remote sensing imagery, net uptake of CO2 (-1.4 ± 0.6 kt CO2 -eq yr-1 ) did not offset CH4 emission (3.7 ± 0.03 kt CO2 -eq yr-1 ), producing an overall positive radiative forcing effect of 2.4 ± 0.3 kt CO2 -eq yr-1 . These results demonstrate clear effects of seasonality, spatial structure, and transport pathway on the magnitude and composition of wetland GHG emissions, and the efficacy of multiscale flux measurement to overcome challenges of wetland heterogeneity.

Salinity causes widespread restriction of methane emissions from small inland waters.

Inland waters are one of the largest natural sources of methane (CH4), a potent greenhouse gas, but emissions models and estimates were developed for solute-poor ecosystems and may not apply to salt-rich inland waters. Here we combine field surveys and eddy covariance measurements to show that salinity constrains microbial CH4 cycling through complex mechanisms, restricting aquatic emissions from one of the largest global hardwater regions (the Canadian Prairies). Existing models overestimated CH4 emissions from ponds and wetlands by up to several orders of magnitude, with discrepancies linked to salinity. While not significant for rivers and larger lakes, salinity interacted with organic matter availability to shape CH4 patterns in small lentic habitats. We estimate that excluding salinity leads to overestimation of emissions from small Canadian Prairie waterbodies by at least 81% ( ~ 1 Tg yr-1 CO2 equivalent), a quantity comparable to other major national emissions sources. Our findings are consistent with patterns in other hardwater landscapes, likely leading to an overestimation of global lentic CH4 emissions. Widespread salinization of inland waters may impact CH4 cycling and should be considered in future projections of aquatic emissions.

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.