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 similar to 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.Peer reviewe

Gap-filling eddy covariance methane fluxes : Comparison of machine learning model predictions and uncertainties at FLUXNET-CH4 wetlands

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 halfhourly 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).Peer reviewe

AmeriFlux FLUXNET-1F US-Tw5 East Pond Wetland

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Tw5 East Pond Wetland. This is the FLUXNET version of the carbon flux data for the site US-Tw5 East Pond Wetland produced by applying the standard ONEFlux (1F) software. Site Description - The Twitchell Wetland site is a 6.5 acre restored wetland on Twitchell Island, that is managed by the California Department of Water Resources (DWR) and the U.S. Geological Survey (USGS). In the fall of 1997, the site was permanently flooded to a depth of approximately 55 cm. The wetland remained fairly unvegetated in patches increasing in size towards the east. The site underwent a major disturbance in 2013 when the vegetation was removed to seed a nearby restored wetland. A flux tower equipped to analyze energy, H2O, CO2, and CH4 fluxes was installed on April 17, 2018

Remotely sensed phenological heterogeneity of restored wetlands: linking vegetation structure and function

Seasonal phenological dynamics of vegetation hold important clues on ecosystem performance towards management goals, such as carbon uptake, and thus should be considered in projections of their targeted services. However, in wetlands spatio-temporal heterogeneity due to mixing of open water, soil, green and dead vegetation makes it difficult to generalize ecosystem functioning across different regions. Remote sensing observations can provide spatially-explicit, cost-effective phenology indicators; however, little is known about their capacity to indicate the links between wetland ecosystem structure and function. Here we assessed this potential by comparing one-year Enhanced Vegetation Index (EVI) from satellite products at high (5m; RapidEye) and low (30m; Landsat) spatial resolutions with eddy covariance time series of net carbon exchange, field digital camera (phenocam) greenness and water temperature among three floristically similar restored wetlands in California, USA. Phenological timing differed by wetland site: depending on satellite, the range in site-median start of greening was up to 28 days, end of greening – up to 73 days, start of senescence – up to 79 days, and end of senescence – up to 10 days. Key transition dates from satellite inputs agreed with seasonal changes in net carbon exchange, phenocam greenness and water temperatures, suggesting that phenological contrasts could result in part from site differences in vegetation configuration and litter affecting the exposure of canopy, soil and water to sunlight and thus sub-canopy microclimate and ecosystem functioning. Yet, the agreement between satellite inputs was non-systematic, with the greatest disparities at the more heterogeneous, less vegetated site. Phenological model fitting uncertainty increased with greater spatial resolution, highlighting the tradeoff between the accuracy of representing vegetation and the complexity of local seasonal variation. These findings highlight the sensitivity of satellite-derived phenology to structural and functional heterogeneity of ecosystems and call for more rigorous spatially-explicit analyses to inform assessments of restoration and management outcomes

Productive wetlands restored for carbon sequestration quickly become net CO2 sinks with site-level factors driving uptake variability.

Inundated wetlands can potentially sequester substantial amounts of soil carbon (C) over the long-term because of slow decomposition and high primary productivity, particularly in climates with long growing seasons. Restoring such wetlands may provide one of several effective negative emission technologies to remove atmospheric CO2 and mitigate climate change. However, there remains considerable uncertainty whether these heterogeneous ecotones are consistent net C sinks and to what degree restoration and management methods affect C sequestration. Since wetland C dynamics are largely driven by climate, it is difficult to draw comparisons across regions. With many restored wetlands having different functional outcomes, we need to better understand the importance of site-specific conditions and how they change over time. We report on 21 site-years of C fluxes using eddy covariance measurements from five restored fresh to brackish wetlands in a Mediterranean climate. The wetlands ranged from 3 to 23 years after restoration and showed that several factors related to restoration methods and site conditions altered the magnitude of C sequestration by affecting vegetation cover and structure. Vegetation established within two years of re-flooding but followed different trajectories depending on design aspects, such as bathymetry-determined water levels, planting methods, and soil nutrients. A minimum of 55% vegetation cover was needed to become a net C sink, which most wetlands achieved once vegetation was established. Established wetlands had a high C sequestration efficiency (i.e. the ratio of net to gross ecosystem productivity) comparable to upland ecosystems but varied between years undergoing boom-bust growth cycles and C uptake strength was susceptible to disturbance events. We highlight the large C sequestration potential of productive inundated marshes, aided by restoration design and management targeted to maximise vegetation extent and minimise disturbance. These findings have important implications for wetland restoration, policy, and management practitioners

Productive wetlands restored for carbon sequestration quickly become net CO2 sinks with site-level factors driving uptake variability.

Inundated wetlands can potentially sequester substantial amounts of soil carbon (C) over the long-term because of slow decomposition and high primary productivity, particularly in climates with long growing seasons. Restoring such wetlands may provide one of several effective negative emission technologies to remove atmospheric CO2 and mitigate climate change. However, there remains considerable uncertainty whether these heterogeneous ecotones are consistent net C sinks and to what degree restoration and management methods affect C sequestration. Since wetland C dynamics are largely driven by climate, it is difficult to draw comparisons across regions. With many restored wetlands having different functional outcomes, we need to better understand the importance of site-specific conditions and how they change over time. We report on 21 site-years of C fluxes using eddy covariance measurements from five restored fresh to brackish wetlands in a Mediterranean climate. The wetlands ranged from 3 to 23 years after restoration and showed that several factors related to restoration methods and site conditions altered the magnitude of C sequestration by affecting vegetation cover and structure. Vegetation established within two years of re-flooding but followed different trajectories depending on design aspects, such as bathymetry-determined water levels, planting methods, and soil nutrients. A minimum of 55% vegetation cover was needed to become a net C sink, which most wetlands achieved once vegetation was established. Established wetlands had a high C sequestration efficiency (i.e. the ratio of net to gross ecosystem productivity) comparable to upland ecosystems but varied between years undergoing boom-bust growth cycles and C uptake strength was susceptible to disturbance events. We highlight the large C sequestration potential of productive inundated marshes, aided by restoration design and management targeted to maximise vegetation extent and minimise disturbance. These findings have important implications for wetland restoration, policy, and management practitioners

Restoring wetlands on intensive agricultural lands modifies nitrogen cycling microbial communities and reduces N2O production potential.

The concentration of nitrous oxide (N2O), an ozone-depleting greenhouse gas, is rapidly increasing in the atmosphere. Most atmospheric N2O originates in terrestrial ecosystems, of which the majority can be attributed to microbial cycling of nitrogen in agricultural soils. Here, we demonstrate how the abundance of nitrogen cycling genes vary across intensively managed agricultural fields and adjacent restored wetlands in the Sacramento-San Joaquin Delta in California, USA. We found that the abundances of nirS and nirK genes were highest at the intensively managed organic-rich cornfield and significantly outnumber any other gene abundances, suggesting very high N2O production potential. The quantity of nitrogen transforming genes, particularly those responsible for denitrification, nitrification and DNRA, were highest in the agricultural sites, whereas nitrogen fixation and ANAMMOX was strongly associated with the wetland sites. Although the abundance of nosZ genes was also high at the agricultural sites, the ratio of nosZ genes to nir genes was significantly higher in wetland sites indicating that these sites could act as a sink of N2O. These findings suggest that wetland restoration could be a promising natural climate solution not only for carbon sequestration but also for reduced N2O emissions