Multiple microbial guilds mediate soil methane cycling along a wetland salinity gradient.

Estuarine wetlands harbor considerable carbon stocks, but rising sea levels could affect their ability to sequester soil carbon as well as their potential to emit methane (CH4). While sulfate loading from seawater intrusion may reduce CH4 production due to the higher energy yield of microbial sulfate reduction, existing studies suggest other factors are likely at play. Our study of 11 wetland complexes spanning a natural salinity and productivity gradient across the San Francisco Bay and Delta found that while CH4 fluxes generally declined with salinity, they were highest in oligohaline wetlands (ca. 3-ppt salinity). Methanogens and methanogenesis genes were weakly correlated with CH4 fluxes but alone did not explain the highest rates observed. Taxonomic and functional gene data suggested that other microbial guilds that influence carbon and nitrogen cycling need to be accounted for to better predict CH4 fluxes at landscape scales. Higher methane production occurring near the freshwater boundary with slight salinization (and sulfate incursion) might result from increased sulfate-reducing fermenter and syntrophic populations, which can produce substrates used by methanogens. Moreover, higher salinities can solubilize ionically bound ammonium abundant in the lower salinity wetland soils examined here, which could inhibit methanotrophs and potentially contribute to greater CH4 fluxes observed in oligohaline sediments.IMPORTANCELow-level salinity intrusion could increase CH4 flux in tidal freshwater wetlands, while higher levels of salinization might instead decrease CH4 fluxes. High CH4 emissions in oligohaline sites are concerning because seawater intrusion will cause tidal freshwater wetlands to become oligohaline. Methanogenesis genes alone did not account for landscape patterns of CH4 fluxes, suggesting mechanisms altering methanogenesis, methanotrophy, nitrogen cycling, and ammonium release, and increasing decomposition and syntrophic bacterial populations could contribute to increases in net CH4 flux at oligohaline salinities. Improved understanding of these influences on net CH4 emissions could improve restoration efforts and accounting of carbon sequestration in estuarine wetlands. More pristine reference sites may have older and more abundant organic matter with higher carbon:nitrogen compared to wetlands impacted by agricultural activity and may present different interactions between salinity and CH4. This distinction might be critical for modeling efforts to scale up biogeochemical process interactions in estuarine wetlands.

Practical Guide to Measuring Wetland Carbon Pools and Fluxes.

Wetlands cover a small portion of the world, but have disproportionate influence on global carbon (C) sequestration, carbon dioxide and methane emissions, and aquatic C fluxes. However, the underlying biogeochemical processes that affect wetland C pools and fluxes are complex and dynamic, making measurements of wetland C challenging. Over decades of research, many observational, experimental, and analytical approaches have been developed to understand and quantify pools and fluxes of wetland C. Sampling approaches range in their representation of wetland C from short to long timeframes and local to landscape spatial scales. This review summarizes common and cutting-edge methodological approaches for quantifying wetland C pools and fluxes. We first define each of the major C pools and fluxes and provide rationale for their importance to wetland C dynamics. For each approach, we clarify what component of wetland C is measured and its spatial and temporal representativeness and constraints. We describe practical considerations for each approach, such as where and when an approach is typically used, who can conduct the measurements (expertise, training requirements), and how approaches are conducted, including considerations on equipment complexity and costs. Finally, we review key covariates and ancillary measurements that enhance the interpretation of findings and facilitate model development. The protocols that we describe to measure soil, water, vegetation, and gases are also relevant for related disciplines such as ecology. Improved quality and consistency of data collection and reporting across studies will help reduce global uncertainties and develop management strategies to use wetlands as nature-based climate solutions. Supplementary Information: The online version contains supplementary material available at 10.1007/s13157-023-01722-2.

AmeriFlux BASE data pipeline to support network growth and data sharing.

AmeriFlux is a network of research sites that measure carbon, water, and energy fluxes between ecosystems and the atmosphere using the eddy covariance technique to study a variety of Earth science questions. AmeriFlux's diversity of ecosystems, instruments, and data-processing routines create challenges for data standardization, quality assurance, and sharing across the network. To address these challenges, the AmeriFlux Management Project (AMP) designed and implemented the BASE data-processing pipeline. The pipeline begins with data uploaded by the site teams, followed by the AMP team's quality assurance and quality control (QA/QC), ingestion of site metadata, and publication of the BASE data product. The semi-automated pipeline enables us to keep pace with the rapid growth of the network. As of 2022, the AmeriFlux BASE data product contains 3,130 site years of data from 444 sites, with standardized units and variable names of more than 60 common variables, representing the largest long-term data repository for flux-met data in the world. The standardized, quality-ensured data product facilitates multisite comparisons, model evaluations, and data syntheses.

Carbon-sink potential of continuous alfalfa agriculture lowered by short-term nitrous oxide emission events.

Alfalfa is the most widely grown forage crop worldwide and is thought to be a significant carbon sink due to high productivity, extensive root systems, and nitrogen-fixation. However, these conditions may increase nitrous oxide (N2O) emissions thus lowering the climate change mitigation potential. We used a suite of long-term automated instrumentation and satellite imagery to quantify patterns and drivers of greenhouse gas fluxes in a continuous alfalfa agroecosystem in California. We show that this continuous alfalfa system was a large N2O source (624 ± 28 mg N2O m2 y-1), offsetting the ecosystem carbon (carbon dioxide (CO2) and methane (CH4)) sink by up to 14% annually. Short-term N2O emissions events (i.e., hot moments) accounted for ≤1% of measurements but up to 57% of annual emissions. Seasonal and daily trends in rainfall and irrigation were the primary drivers of hot moments of N2O emissions. Significant coherence between satellite-derived photosynthetic activity and N2O fluxes suggested plant activity was an important driver of background emissions. Combined data show annual N2O emissions can significantly lower the carbon-sink potential of continuous alfalfa agriculture.

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.

The global spectrum of plant form and function: enhanced species-level trait dataset.

Here we provide the 'Global Spectrum of Plant Form and Function Dataset', containing species mean values for six vascular plant traits. Together, these traits -plant height, stem specific density, leaf area, leaf mass per area, leaf nitrogen content per dry mass, and diaspore (seed or spore) mass - define the primary axes of variation in plant form and function. The dataset is based on ca. 1 million trait records received via the TRY database (representing ca. 2,500 original publications) and additional unpublished data. It provides 92,159 species mean values for the six traits, covering 46,047 species. The data are complemented by higher-level taxonomic classification and six categorical traits (woodiness, growth form, succulence, adaptation to terrestrial or aquatic habitats, nutrition type and leaf type). Data quality management is based on a probabilistic approach combined with comprehensive validation against expert knowledge and external information. Intense data acquisition and thorough quality control produced the largest and, to our knowledge, most accurate compilation of empirically observed vascular plant species mean traits to date.

Detecting Hot Spots of Methane Flux Using Footprint-Weighted Flux Maps.

In this study, we propose a new technique for mapping the spatial heterogeneity in gas exchange around flux towers using flux footprint modeling and focusing on detecting hot spots of methane (CH4) flux. In the first part of the study, we used a CH4 release experiment to evaluate three common flux footprint models: the Hsieh model (Hsieh et al., 2000), the Kljun model (Kljun et al., 2015), and the K & M model (Kormann and Meixner, 2001), finding that the K & M model was the most accurate under these conditions. In the second part of the study, we introduce the Footprint-Weighted Flux Map, a new technique to map spatial heterogeneity in fluxes. Using artificial CH4 release experiments, natural tracer approaches and flux chambers we mapped the spatial flux heterogeneity, and detected and validated a hot spot of CH4 flux in a oligohaline restored marsh. Through chamber measurements during the months of April and May, we found that fluxes at the hot spot were on average as high as 6589 ± 7889 nmol m-2 s-1 whereas background flux from the open water were on average 15.2 ± 7.5 nmol m-2 s-1. This study provides a novel tool to evaluate the spatial heterogeneity of fluxes around eddy-covariance towers and creates important insights for the interpretation of hot spots of CH4 flux, paving the way for future studies aiming to understand subsurface biogeochemical processes and the microbiological conditions that lead to the occurrence of hot spots and hot moments of CH4 flux.

Cross-biome synthesis of source versus sink limits to tree growth.

Uncertainties surrounding tree carbon allocation to growth are a major limitation to projections of forest carbon sequestration and response to climate change. The prevalence and extent to which carbon assimilation (source) or cambial activity (sink) mediate wood production are fundamentally important and remain elusive. We quantified source-sink relations across biomes by combining eddy-covariance gross primary production with extensive on-site and regional tree ring observations. We found widespread temporal decoupling between carbon assimilation and tree growth, underpinned by contrasting climatic sensitivities of these two processes. Substantial differences in assimilation-growth decoupling between angiosperms and gymnosperms were determined, as well as stronger decoupling with canopy closure, aridity, and decreasing temperatures. Our results reveal pervasive sink control over tree growth that is likely to be increasingly prominent under global climate change.

Confronting the water potential information gap.

Water potential directly controls the function of leaves, roots, and microbes, and gradients in water potential drive water flows throughout the soil-plant-atmosphere continuum. Notwithstanding its clear relevance for many ecosystem processes, soil water potential is rarely measured in-situ, and plant water potential observations are generally discrete, sparse, and not yet aggregated into accessible databases. These gaps limit our conceptual understanding of biophysical responses to moisture stress and inject large uncertainty into hydrologic and land surface models. Here, we outline the conceptual and predictive gains that could be made with more continuous and discoverable observations of water potential in soils and plants. We discuss improvements to sensor technologies that facilitate in situ characterization of water potential, as well as strategies for building new networks that aggregate water potential data across sites. We end by highlighting novel opportunities for linking more representative site-level observations of water potential to remotely-sensed proxies. Together, these considerations offer a roadmap for clearer links between ecohydrological processes and the water potential gradients that have the 'potential' to substantially reduce conceptual and modeling uncertainties.

A novel approach to partitioning evapotranspiration into evaporation and transpiration in flooded ecosystems.

Reliable partitioning of micrometeorologically measured evapotranspiration (ET) into evaporation (E) and transpiration (T) would greatly enhance our understanding of the water cycle and its response to climate change related shifts in local-to-regional climate conditions and rising global levels of vapor pressure deficit (VPD). While some methods on ET partitioning have been developed, their underlying assumptions make them difficult to apply more generally, especially in sites with large contributions of E. Here, we report a novel ET partitioning method using artificial neural networks (ANNs) in combination with a range of environmental input variables to predict daytime E from nighttime ET measurements. The study uses eddy covariance data from four restored wetlands in the Sacramento-San Joaquin Delta, California, USA, as well as leaf-level T data for validation. The four wetlands vary in their vegetation make-up and structure, representing a range of ET conditions. The ANNs were built with increasing complexity by adding the input variable that resulted in the next highest average value of model testing R2 across all sites. The order of variable inclusion (and importance) was: VPD > gap-filled sensible heat flux (H_gf) > air temperature (Tair ) > friction velocity (u∗ ) > other variables. The model using VPD, H_gf, Tair , and u∗ showed the best performance during validation with independent data and had a mean testing R2  value of 0.853 (averaged across all sites, range from 0.728 to 0.910). In comparison to other methods, our ANN method generated T/ET partitioning results which were more consistent with CO2 exchange data especially for more heterogeneous sites with large E contributions. Our method improves the understanding of T/ET partitioning. While it may be particularly suited to flooded ecosystems, it can also improve T/ET partitioning in other systems, increasing our knowledge of the global water cycle and ecosystem functioning.

A remote sensing-based three-source energy balance model to improve global estimations of evapotranspiration in semi-arid tree-grass ecosystems.

It is well documented that energy balance and other remote sensing-based evapotranspiration (ET) models face greater uncertainty over water-limited tree-grass ecosystems (TGEs), representing nearly 1/6th of the global land surface. Their dual vegetation strata, the grass-dominated understory and tree-dominated overstory, make for distinct structural, physiological and phenological characteristics, which challenge models compared to more homogeneous and energy-limited ecosystems. Along with this, the contribution of grasses and trees to total transpiration (T), along with their different climatic drivers, is still largely unknown nor quantified in TGEs. This study proposes a thermal-based three-source energy balance (3SEB) model, accommodating an additional vegetation source within the well-known two-source energy balance (TSEB) model. The model was implemented at both tower and continental scales using eddy-covariance (EC) TGE sites, with variable tree canopy cover and rainfall (P) regimes and Meteosat Second Generation (MSG) images. 3SEB robustly simulated latent heat (LE) and related energy fluxes in all sites (Tower: LE RMSD ~60 W/m2 ; MSG: LE RMSD ~90 W/m2 ), improving over both TSEB and seasonally changing TSEB (TSEB-2S) models. In addition, 3SEB inherently partitions water fluxes between the tree, grass and soil sources. The modelled T correlated well with EC T estimates (r > .76), derived from a machine learning ET partitioning method. The T/ET was found positively related to both P and leaf area index, especially compared to the decomposed grass understory T/ET. However, trees and grasses had contrasting relations with respect to monthly P. These results demonstrate the importance in decomposing total ET into the different vegetation sources, as they have distinct climatic drivers, and hence, different relations to seasonal water availability. These promising results improved ET and energy flux estimations over complex TGEs, which may contribute to enhance global drought monitoring and understanding, and their responses to climate change feedbacks.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.

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.

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.

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.

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.

An Ecosystem-Scale Flux Measurement Strategy to Assess Natural Climate Solutions.

Eddy covariance measurement systems provide direct observation of the exchange of greenhouse gases between ecosystems and the atmosphere, but have only occasionally been intentionally applied to quantify the carbon dynamics associated with specific climate mitigation strategies. Natural climate solutions (NCS) harness the photosynthetic power of ecosystems to avoid emissions and remove atmospheric carbon dioxide (CO2), sequestering it in biological carbon pools. In this perspective, we aim to determine which kinds of NCS strategies are most suitable for ecosystem-scale flux measurements and how these measurements should be deployed for diverse NCS scales and goals. We find that ecosystem-scale flux measurements bring unique value when assessing NCS strategies characterized by inaccessible and hard-to-observe carbon pool changes, important non-CO2 greenhouse gas fluxes, the potential for biophysical impacts, or dynamic successional changes. We propose three deployment types for ecosystem-scale flux measurements at various NCS scales to constrain wide uncertainties and chart a workable path forward: "pilot", "upscale", and "monitor". Together, the integration of ecosystem-scale flux measurements by the NCS community and the prioritization of NCS measurements by the flux community, have the potential to improve accounting in ways that capture the net impacts, unintended feedbacks, and on-the-ground specifics of a wide range of emerging NCS strategies.