Search

Wetlands are responsible for 20%–31% of global methane (CH4) emissions and account for a large source of uncertainty in the global CH4 budget. Data-driven upscaling of CH4 fluxes from eddy covariance measurements can provide new and independent bottom-up estimates of wetland CH4 emissions. Here, we develop a six-predictor random forest upscaling model (UpCH4), trained on 119 site-years of eddy covariance CH4 flux data from 43 freshwater wetland sites in the FLUXNET-CH4 Community Product. Network patterns in site-level annual means and mean seasonal cycles of CH4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash-Sutcliffe Efficiency ∼0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH4 emissions of 146 ± 43 TgCH4 y−1 for 2001–2018 which agrees closely with current bottom-up land surface models (102–181 TgCH4 y−1) and overlaps with top-down atmospheric inversion models (155–200 TgCH4 y−1). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH4 fluxes has the potential to produce realistic extra-tropical wetland CH4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid-to-arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions. The monthly wetland methane products gridded at 0.25° from UpCH4 are available via ORNL DAAC (https://doi.org/10.3334/ORNLDAAC/2253).

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

The continued warming of the Arctic could release vast stores of carbon into the atmosphere from high-latitude ecosystems, especially from thawing permafrost. Increasing uptake of carbon dioxide (CO2) by vegetation during longer growing seasons may partially offset such release of carbon. However, evidence of significant net annual release of carbon from site-level observations and model simulations across tundra ecosystems has been inconclusive. To address this knowledge gap, we combined top-down observations of atmospheric CO2 concentration enhancements from aircraft and a tall tower, which integrate ecosystem exchange over large regions, with bottom-up observed CO2 fluxes from tundra environments and found that the Alaska North Slope is not a consistent net source nor net sink of CO2 to the atmosphere (ranging from -6 to +6TgCyr-1 for 2012-2017). Our analysis suggests that significant biogenic CO2 fluxes from unfrozen terrestrial soils, and likely inland waters, during the early cold season (September-December) are major factors in determining the net annual carbon balance of the North Slope, implying strong sensitivity to the rapidly warming freeze-up period. At the regional level, we find no evidence of the previously reported large late-cold-season (January-April) CO2 emissions to the atmosphere during the study period. Despite the importance of the cold-season CO2 emissions to the annual total, the interannual variability in the net CO2 flux is driven by the variability in growing season fluxes. During the growing season, the regional net CO2 flux is also highly sensitive to the distribution of tundra vegetation types throughout the North Slope. This study shows that quantification and characterization of year-round CO2 fluxes from the heterogeneous terrestrial and aquatic ecosystems in the Arctic using both site-level and atmospheric observations are important to accurately project the Earth system response to future warming.

Arctic warming is affecting snow cover and soil hydrology, with consequences for carbon sequestration in tundra ecosystems. The scarcity of observations in the Arctic has limited our understanding of the impact of covarying environmental drivers on the carbon balance of tundra ecosystems. In this study, we address some of these uncertainties through a novel record of 119 site-years of summer data from eddy covariance towers representing dominant tundra vegetation types located on continuous permafrost in the Arctic. Here we found that earlier snowmelt was associated with more tundra net CO2 sequestration and higher gross primary productivity (GPP) only in June and July, but with lower net carbon sequestration and lower GPP in August. Although higher evapotranspiration (ET) can result in soil drying with the progression of the summer, we did not find significantly lower soil moisture with earlier snowmelt, nor evidence that water stress affected GPP in the late growing season. Our results suggest that the expected increased CO2 sequestration arising from Arctic warming and the associated increase in growing season length may not materialize if tundra ecosystems are not able to continue sequestering CO2 later in the season.

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

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.

Cold air pooling (CAP) occurs in low-lying areas where cold, dense air collects during nighttime hours, producing colder temperatures than surrounding higher elevations. Devils Postpile National Monument (DEPO) is confined by steep mountain ridges, which promote cold air drainage into lower-elevation meadows and river valleys. These low lying areas where CAP occurs could help facilitate a potential refugium from some of the greatest impacts of regional climate warming. A strong focus on CAP occurrence is part of a seven-step climate change refugia conservation cycle, outlined in Morelli et al. (2016), that DEPO has instituted, wherein resource management seeks to identify and focus on parts of the landscape that may be sheltered from the intensity and pace of rapid climate change. Locales that harbor persistent CAP may provide vulnerable species and ecosystems a sort of refuge, allowing more time to adapt to new conditions. Central to DEPO’s strategy is better monitoring and understanding of CAP occurrence. Based on observations from 10 years of temperature loggers in and around DEPO, CAP operates reliably throughout the year, and especially during summer. Importantly, CAP has occurred strongly during recent unusually warm years. These findings reinforce the value of monitoring and ongoing analysis as a way to guide conservation and adaptation using potential climate change refugia. As with this co-developed park–university investigation, other land managers could consider climate refugia-oriented management as a viable conservation and adaptation strategy.

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 due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, 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 open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). 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 which are a substantial source of total atmospheric CH4 emissions; and (3) we 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 20g g€¯S to 20g g€¯N) the spring onset of elevated CH4 emissions starts 3g€¯d earlier, and the CH4 emission season lasts 4g€¯d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 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). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at 10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

California’s Central Valley (CV) is one of the most productive agricultural regions in the world, enabled by the conjunctive use of surface water and groundwater. We investigated variations in the CV’s managed surface water diversions relative to climate variability. Using a historical record (1979−2010) of diversions from 531 sites, we found diversions are largest in the wetter Sacramento basin to the north, but most variable in the drier Tulare basin to the south. A rotated empirical orthogonal function (REOF) analysis finds 72% of the variance of diversions is captured by the first three REOFs. The leading REOF (35% of variance) exhibited strong positive loadings in the Tulare basin, and the corresponding principal component time-series (RPC1) was strongly correlated (ρ > 0.9) with contemporaneous hydrologic variability. This pattern indicates larger than average diversions in the south, with neutral or slightly less than average diversions to the north during wet years, with the opposite true for dry years. The second and third REOFs (20% and 17% of variance, respectively), were strongest in the Sacramento basin and San Francisco Bay−Delta. RPC2 and RPC3 were associated with variations in agricultural- and municipal-bound diversions, respectively. RPC2 and RPC3 were also moderately correlated with 7-year cumulative precipitation based on lagged correlation analysis, indicating that diversions in the north and central portions of the CV respond to longer-term hydrologic variations. The results illustrate a dichotomy of regimes wherein diversions in the more arid Tulare are governed by year-to-year hydrologic variability, while those in wetter northern basins reflect land-use patterns and low-frequency hydrologic variations.

Recent warming in the Arctic, which has been amplified during the winter1-3, greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO2)4. However, the amount of CO2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates5,6. Here we synthesize regional in situ observations of CO2 flux from arctic and boreal soils to assess current and future winter carbon losses from the northern permafrost domain. We estimate a contemporary loss of 1662 Tg C yr-1 from the permafrost region during the winter season (October through April). This loss is greater than the average growing season carbon uptake for this region estimated from process models (-1032 Tg C yr-1). Extending model predictions to warmer conditions in 2100 indicates that winter CO2 emissions will increase 17% under a moderate mitigation scenario-Representative Concentration Pathway (RCP) 4.5-and 41% under business-as-usual emissions scenario-RCP 8.5. Our results provide a new baseline for winter CO2 emissions from northern terrestrial regions and indicate that enhanced soil CO2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.

Arctic terrestrial ecosystems are major global sources of methane (CH4); hence, it is important to understand the seasonal and climatic controls on CH4 emissions from these systems. Here, we report year-round CH4 emissions from Alaskan Arctic tundra eddy flux sites and regional fluxes derived from aircraft data. We find that emissions during the cold season (September to May) account for ≥ 50% of the annual CH4 flux, with the highest emissions from noninundated upland tundra. A major fraction of cold season emissions occur during the zero curtain period, when subsurface soil temperatures are poised near 0 °C. The zero curtain may persist longer than the growing season, and CH4 emissions are enhanced when the duration is extended by a deep thawed layer as can occur with thick snow cover. Regional scale fluxes of CH4 derived from aircraft data demonstrate the large spatial extent of late season CH4 emissions. Scaled to the circumpolar Arctic, cold season fluxes from tundra total 12 ± 5 (95% confidence interval) Tg CH4 y(-1), ∼ 25% of global emissions from extratropical wetlands, or ∼ 6% of total global wetland methane emissions. The dominance of late-season emissions, sensitivity to soil environmental conditions, and importance of dry tundra are not currently simulated in most global climate models. Because Arctic warming disproportionally impacts the cold season, our results suggest that higher cold-season CH4 emissions will result from observed and predicted increases in snow thickness, active layer depth, and soil temperature, representing important positive feedbacks on climate warming.

Understanding the sources and sinks of CH4 is critical to both predicting and mitigating future climate change. There are large uncertainties in the global budget of atmospheric CH4, but natural emissions are estimated to be of a similar magnitude to total anthropogenic emissions. The largest sources of uncertainty in scaling bottom-up CH4 estimates stem from limited ground-based measurements and the misalignment between drivers of CH4 fluxes and current land use classifications. To understand the CH4 flux potential of natural ecosystems and agricultural lands in the United States (US) of America, a multi-scale CH4 observation network focused on CH4 flux rates, processes, and scaling methods is required. This can be achieved with a network of ground-based observations that are distributed based on climatic regions and landcover. To determine the gaps in physical infrastructure for developing this network, we need to understand the representativeness of current measurements. We focus here on eddy covariance (EC) flux towers because they are essential for a bottom-up framework that bridges the gap between point-based chamber measurements and airborne or satellite platforms, informing the remote sensing and modelling communities and policy decisions, all the way to IPCC reports. Using multidimensional scaling and a cluster analysis, the US was divided into 10 clusters that were distributed across temperature and wetness gradients. We evaluated the distance to the medoid condition within each cluster for research sites with EC tower infrastructure to identify the gaps in existing infrastructure that limit our ability to constrain the contribution of US biogenic CH4 emissions to the global budget. These gaps occurred across all EC flux tower networks and independently managed sites as well as in some environmental clusters. Through our analysis using climate, land cover, and location variables, we have identified priority areas to target for research infrastructure to provide a more complete understanding of the CH4 flux potential of ecosystem types across the US. [ABSTRACT FROM AUTHOR]

Highlights: [•] We measured net CO2 exchange (NEE) in a remnant raised peat bog for two years. [•] Mean annual NEE was −234gCm−2 yr−1 indicating substantial CO2 sink strength. [•] Year-round growing conditions accounted for large component CO2 fluxes and NEE. [•] Summertime light response of NEE was similar to northern hemisphere peatlands. [•] Small differences in cool season NEE accounted for differences in annual NEE. [Copyright &y& Elsevier]

Drained peatlands are major sources of CO 2 to the atmosphere, yet the effects of land management and hydrological extremes have been little-studied at spatial scales relevant to agricultural enterprises. We measured fluxes of CO 2 using the eddy covariance (EC) technique at two adjacent dairy farms on a drained peatland in Aotearoa New Zealand with remaining peat depths 5.5-8 m. One site (SD) had shallow surface drains and mean water table depth (WTD) -657 mm, while the other site (BD) had deep field border drains and mean WTD -838 mm. Net ecosystem CO 2 production (NEP) was similar at the two sites when the soils were moist but diverged during late-summer drying, with site BD having 4.56 t C ha -1 greater CO 2 emission than site SD over the four-month dry period. Soil drying reduced gross primary production (GPP) at both sites, while ecosystem respiration (ER) was reduced at site SD but not at site BD. The low dry season respiration rates at site SD contributed to near-zero annual NEP, while higher respiration rates at site BD led to annual CO 2 loss of -4.95 ± 0.59 t C ha -1  yr -1 . Accounting for other imports and exports of carbon, annual net ecosystem carbon balances were -2.23 and -8.47 t C ha -1  yr -1 at sites SD and BD, respectively. It is likely that the contrasting dry season respiration rates resulted from differences in soil physical properties affecting soil moisture vertical redistribution and availability to plants and microbes rather than from the relatively small differences in WTD. These differences could be caused by soil physical disturbances during pasture renewal or paddock recontouring, or time since initial drainage. Therefore, improved soil management might provide practical mitigation against excessive CO 2 emissions during dry conditions, including droughts., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Recent warming in the Arctic, which has been amplified during the winter 1-3 , greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO 2 ) 4 . However, the amount of CO 2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates 5,6 . Here we synthesize regional in situ observations of CO 2 flux from arctic and boreal soils to assess current and future winter carbon losses from the northern permafrost domain. We estimate a contemporary loss of 1662 Tg C yr -1 from the permafrost region during the winter season (October through April). This loss is greater than the average growing season carbon uptake for this region estimated from process models (-1032 Tg C yr -1 ). Extending model predictions to warmer conditions in 2100 indicates that winter CO 2 emissions will increase 17% under a moderate mitigation scenario-Representative Concentration Pathway (RCP) 4.5-and 41% under business-as-usual emissions scenario-RCP 8.5. Our results provide a new baseline for winter CO 2 emissions from northern terrestrial regions and indicate that enhanced soil CO 2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.