Articles | Volume 19, issue 9
Research article
 | Highlight paper
13 May 2022
Research article | Highlight paper |  | 13 May 2022

Gaps in network infrastructure limit our understanding of biogenic methane emissions for the United States

Sparkle L. Malone, Youmi Oh, Kyle A. Arndt, George Burba, Roisin Commane, Alexandra R. Contosta, Jordan P. Goodrich, Henry W. Loescher, Gregory Starr, and Ruth K. Varner

Understanding the sources and sinks of methane (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 anthropogenic emissions. To understand CH4 flux from biogenic sources 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 land cover. To determine the gaps in physical infrastructure for developing this network, we need to understand the landscape representativeness of the current infrastructure. 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 that inform policy decisions and global climate agreements. Using dissimilarity, multidimensional scaling, and cluster analysis, the US was divided into 10 clusters distributed across temperature and precipitation gradients. We evaluated dissimilarity within each cluster for research sites with active CH4 EC towers to identify gaps in existing infrastructure that limit our ability to constrain the contribution of US biogenic CH4 emissions to the global budget. Through our analysis using climate, land cover, and location variables, we identified priority areas for research infrastructure to provide a more complete understanding of the CH4 flux potential of ecosystem types across the US. Clusters corresponding to Alaska and the Rocky Mountains, which are inherently difficult to capture, are the most poorly represented, and all clusters require a greater representation of vegetation types.

1 Introduction

The 21st century is characterized by ongoing changes in Earth's climate system that result from increasing concentrations of radiatively important trace gases in the atmosphere. Unlike the relatively steady increases of atmospheric carbon dioxide (CO2) and nitrous oxide (N2O), atmospheric methane (CH4) concentrations show dynamic trends with a rapid increase of  10 ppb yr−1 since 2014 (Nisbet et al., 2019). The annual increase of atmospheric CH4 in 2020 was the largest on record at  15 ppb yr−1 (Dlugokencky, 2021), despite the global pandemic reducing energy demand (Le Quéré et al., 2021). Increasing atmospheric CH4 concentrations (Nisbet et al., 2019) is of concern because CH4 is 34 times more effective at trapping heat in the atmosphere compared to an equivalent mass of CO2 over a 100-year timeframe and accounts for  42 % of warming since the pre-industrial period (IPCC, 2021). These rapid increases in atmospheric CH4 challenge us to reach the goals of the Paris Agreement (Nisbet et al., 2019) but also provide an opportunity given the relatively short atmospheric residence time ( 9 years) of CH4. Understanding the sources and sinks of CH4 is therefore critical in predicting and mitigating future climate change.

Quantifying the national CH4 budget is important for assessing realistic pathways to mitigate climate change, yet uncertainties in the magnitude, size, and location of sources and sinks limit budget development (Saunois et al., 2020; Bruhwiler et al., 2021). Methane is emitted from a variety of often co-located biogenic, thermogenic, and pyrogenic sources (IPCC, 2013; Nisbet et al., 2019). Biogenic emissions are thought to be of a similar magnitude to total anthropogenic emissions, yet biogenic CH4 emissions remain the most uncertain source of the global CH4 budget (Saunois et al., 2020). Surface–atmosphere exchange from biogenic sources and sinks, the biological and environmental processes driving these fluxes (e.g., ebullition, aerenchyma pumping), and how CH4 sources and sinks change over space and time, including interannual variability (Michalak et al., 2009; Kirschke et al., 2013; Knox et al., 2019; Nisbet et al., 2019), are not well constrained. Finally, both the vast areas with relatively small uptake and emission rates (e.g., deserts, grasslands, forests) and the lake–ocean water continuum that transports CH4 (e.g., fens, streams, and rivers) have been largely understudied but could contribute significantly to regional and global budgets (Hutchins et al., 2019; Rosentreter et al., 2021; Zhou et al., 2021). These unknowns hinder our ability to predict future climate change due to the complex feedbacks between biological processes (e.g., microbial production and consumption) (Sherwood et al., 2017; Zhang et al., 2017; Oh et al., 2020), climate change (Zhang et al., 2017), and land cover change (Kirschke et al., 2013; Knox et al., 2019; Saunois et al., 2020).

To understand the biogenic CH4 flux potential of the United States of America (US), a multi-scale CH4 observation network focused on CH4 flux rates, processes, and scaling methods is required. When scaling bottom-up measurements to the landscape and regional scale, measurements of CH4 from existing infrastructure tend not to be sufficiently geographically distributed to capture the true spatial variation that is innate to the production and consumption of CH4, and is compounded by large source/sink strengths in small areas (e.g., periodic wetting/drying of seasonal wetlands, saturated soils) (IPCC, 2013; Knox et al., 2019; Thornton et al., 2016) and by very small source/sink strengths in very large areas. In addition, bottom-up biogenic CH4 process-level estimates have historically been limited to short periods (< 1–2 years), are discontinuous (grab sampling), and/or occur only during the growing season at middle and high latitudes (though see Groffman et al., 2006, and Arndt et al., 2019, for notable exceptions).

There is a pressing need to assess the capacity of existing infrastructure for current and future applications (Lovett et al., 2007; Kumar et al., 2016; Jongman et al., 2017; Novick et al., 2018; Villarreal et al., 2018; Chu et al., 2021). The representativeness of research infrastructure is often described in terms of the extent to which the measurements collected at any given location and time represent the conditions at any other location and time, and this is often driven by ecological and climatic conditions (Sulkava et al., 2011; Chu et al., 2021). Representativeness is also measured across landscapes, and studies have evaluated how well tower infrastructure captures the variability observed at specific sites (Chu et al., 2021). These approaches seek to understand the representativeness of the measurements for a broader landscape, which is critical for upscaling point measurements to regional and global scales. These types of assessments inform the scientific community on how to increase their utility and are often designed to support network design, upscaling, and bias estimation (Chen et al., 2011; Ciais et al., 2014; Jongman et al., 2017; Schimel and Keller, 2015; Villarreal et al., 2018; Kumar et al., 2016). There have been many attempts to assess the representativeness of existing eddy covariance (EC) tower networks for various purposes. To date, no study has focused on CH4 infrastructure across the US, though many studies have used clustering and ecoregions (Sulkava et al., 2011; Hargrove and Hoffman, 2003), dissimilarity (Yang et al., 2008), and distance measures (Hargrove and Hoffman, 2003; Yang et al., 2008; He et al., 2015; Hoffman et al., 2013) on climatic (Novick et al., 2018) and vegetation type structure and function (Chu et al., 2021) to measure the representativeness of existing research infrastructure. The primary goal of this work is to fill this key knowledge gap by determining the regions where biogenic CH4 infrastructure is needed within the US in order to constrain both the national and global CH4 budget.

To determine key regions where biogenic CH4 infrastructure is needed within the US, we statistically identify gaps in active research infrastructure and evaluate areas where infrastructure can be augmented to include new CH4 measurements. We use a combination of climate data and dominant land cover to guide the scientific community on how we can develop a distributed observational network for the US by leveraging existing infrastructure. While this analysis does not capture the heterogeneity of the conditions that drive CH4 fluxes at the ecosystem scale, it is designed to evaluate the sampling intensity of research sites at the landscape scale. This coarse resolution influences the capacity to scale ecosystem-level results to the landscape, regional, and national level, which is required for the development of CH4 budgets and emission reduction strategies.

2 Methods

2.1 Overview

To determine the gaps in physical research infrastructure for ecosystem-scale CH4 fluxes, we need to understand how the current infrastructure is distributed across the US. We focus here on EC flux towers given their capabilities for continuous measurements and use in upscaling flux estimates and are therefore a useful basis for identifying gaps in the current network of CH4 observations. The AmeriFlux network of EC towers was launched in 1996 and grew from about 15 sites in 1997 to more than 110 active sites registered today. It was originally a network of PI-managed sites measuring ecosystem CO2, H2O, and energy fluxes. The network was established to connect research on field sites representing major climatic and ecological biomes, including tundra, grasslands, savanna, crops, and coniferous, deciduous, and tropical forests. The AmeriFlux community tailored instrumentation to suit each unique ecosystem but now also includes towers that are a part of the standardized network, the National Ecological Observatory Network (NEON). In 2012, the US Department of Energy established the AmeriFlux Management Project (AMP) at Lawrence Berkeley National Laboratory (LBNL) to support the broad AmeriFlux community and the AmeriFlux sites. The AMP standardizes, post-processes, and makes flux data available to the research community. More recently, flux towers began measuring CH4 in freshwater, coastal, upland, natural, and managed ecosystems. Although we have information on the location of existing EC tower infrastructure that is a part of AmeriFlux (n=223), NEON (n=47), and known, independent PI-managed sites (n=141), we focus this analysis on the towers measuring CH4 (n=100) and we distinguish between towers providing data to AmeriFlux (yes = 49, no = 51) and tower activity (active = 70; inactive = 30). We understand that additional towers exist within the US, but because these towers are not reporting or providing data to the flux community, we cannot include them in this analysis.

To understand the landscape representativeness across geographic clusters, we measured dissimilarity based on climate and land cover type, as these two factors together are characteristic of regional resource availability and disturbance regimes. First, we developed a dissimilarity matrix that was condensed down to a two-dimensional ordination to determine regional clusters and calculate cluster dissimilarity for each location within a cluster (Fig. 1). It is important to note that a tower should be representative of the ecosystem type and the region where it is stationed (Desai, 2010; Jung et al., 2011; Xiao et al., 2012; Chu et al., 2021); however, the landscape representativeness analysis done here uses a coarser classification of land cover types that are more emblematic of resource availability and factors that influence how ecosystems function, not the specific ecosystem type where the tower is situated. Chu et al. (2021) examined the land-cover composition and vegetation characteristics of 214 AmeriFlux tower site footprints. They found that most sites do not represent the dominant land-cover type of the ecosystems they exist within, and when paired with common model–data integration approaches this mismatch introduces biases on the order of 4 %–20 % for the enhanced vegetation index (EVI) and 6 %–20 % for the dominant land cover percentage (Chu et al., 2021), making it essential to consider landscape characteristics in the design and evaluation of network infrastructure. Infrastructure representativeness at the landscape scale is indicative of the capacity to upscale information by climate and the dominant ecosystems of locations within a landscape.

Figure 1To determine the gaps in physical research infrastructure for CH4 fluxes we measured landscape cover and climate dissimilarity across the US and evaluated the current distribution of CH4 tower infrastructure.


2.2 Climate and dominant land cover types

We used the National Land Cover Database (NLCD;, last access: 1 October 2021) to create a land cover layer for the contiguous US (Jin et al., 2019). The NLCD has a 30 m resolution with a 16-class legend based on a modified Anderson Level II classification system. We reclassified the NLCD into eight major land cover types (water, developed, barren, forest, scrub, herbaceous, crop, and wetland). Where the NLCD was not available (Alaska, Hawaii, and Puerto Rico), we used the Moderate Resolution Imaging Spectroradiometer (MODIS; 1 km) land cover (type 5 – vegetation functional types) for vegetation functional type (MCD12Q1.006) (Sulla-Menashe and Friedl, 2018), which was also reclassified to the eight major land cover types (Table 1). The crop land cover type was expanded to non-irrigated and irrigated classes using agricultural information from the US Department of Agriculture's CropScape and Cropland Data layer (Boryan et al., 2011), and the wetland class was expanded using information from the US Fish and Wildlife Service's National Wetland Inventory. Expanded wetland classes were emergent coastal, emergent freshwater, and forest freshwater (Wilen and Bates, 1995). Climate data were obtained from DAYMET (Thornton et al., 2017). We used five climate variables to characterize the climatic conditions across the US: annual mean daily minimum, daily average, and daily maximum temperatures, annual total precipitation, and mean annual daily vapor pressure deficit from 2010–2020. Understanding that these patterns are changing with climate change, we chose a shorter time period than the commonly used 30-year climate normal to better represent current conditions (Bessembinder et al., 2021). Land cover was resampled to match the DAYMET climate data (1 km), and all pre-processing was done in R version 4.0.4 (R Core Team, 2021) with the raster package (Hijmans, 2021). This approach allowed us to create a land cover layer of the dominant land cover types at 1 km resolution that was expanded in categories of interest for CH4. The land cover and climate layers were chosen to represent the primary environmental conditions that are often indicative of a combination of resource availability and disturbance regimes. These coarse layers are essential for considering the landscape and large-scale climate effects that can influence how ecosystems within landscapes function. While the available land cover information is appropriate for the coarse, landscape-scale analysis done here, it is important to note that the products used here are not designed to estimate the potential CH4 source/sink status, particularly from the aquatic, wetland, and agricultural land cover types.

Table 1Land cover and data sources. The blended land cover product comprises the National Land Cover Database (NLCD) and Moderate Resolution Imaging Spectroradiometer (MODIS). The crop category is enhanced with CropScape and the wetland category with the National Wetland Inventory (NWI) to identify areas dominated by land cover types with additional classes added for types with expected CH4 source potential.

n/a: not applicable.

Download Print Version | Download XLSX

2.3 Measuring landscape dissimilarity across clusters within the US

Climate, land cover, and location (latitude/longitude) were used in a multivariate distance analysis (Venables and Ripley, 2002; Ripley, 2007; Cox and Cox, 2008) to measure the dissimilarity across the US (all 50 states and Puerto Rico) at the landscape scale and divide it into ecological clusters. The purpose of this is to identify the interrelatedness of ecological components within a landscape (Ippoliti et al., 2019). We included location (latitude/longitude) to incorporate the interaction between climate, land cover, and most importantly, seasonality. The US was subsampled because of limitations in the maximum number of points that can be evaluated in the cluster analysis. To measure dissimilarity, we first randomly sampled (n= 20 000 1 km pixels) the US, maintaining the distribution of land cover and climate to define dissimilarity between observations. Although there were more than 8 million 1 km pixels available for the US, there are limits to the number of samples that can be analyzed by the functions used for the multidimensional scaling (MDS) analysis. We first developed a dissimilarity matrix by calculating Gower dissimilarity (Gower, 1971; Huang, 1997; Podani, 1999; Ahmad and Dey, 2007; Harikumar and Pv, 2015) using the function distmix from the package kmed in R. We used Gower dissimilarity because it can handle mixed data types. For each variable type in the data set, the dissimilarity metric that works well for that type is used and scaled to fall between 0 and 1. Then, a linear combination featuring user-specified weights (most simply an average) is calculated to create the final dissimilarity matrix. This approach measures the dissimilarity for each location within the US using land cover, climate, and location information (land cover, five climate variables, and location) and creates a dissimilarity matrix (20 000 × 20 000) that indicates dissimilarity for a location to every other location in the US.

Once we created the dissimilarity matrix, we used MDS to generate a two-dimensional ordination showing landscape dissimilarity with the MASS package in R (Venables and Ripley, 2002). The MDS makes it possible to evaluate dissimilarity in two dimensions, which is essential to our goal to evaluate representativeness. We used the Kruskal method of non-metric scaling with the IsoMDS function in the MASS package (Venables and Ripley, 2002). IsoMDS works best when applied to metric variables (Torgerson, 1958). Torgerson (1958) initially developed this method, which assumes that the data obey distance axioms. It uses eigendecomposition of the dissimilarity to identify major components and axes and represents any point as a linear combination of dimensions. This is very similar to principal component analysis (PCA) or factor analysis, but it uses the dissimilarity matrix rather than a correlation matrix as input. Furthermore, the included dimensions are the most important dimensions produced, like PCA which is able to identify all of the dimensions that exist in the original data up to N−1, but will retain only the most important ones.

Knowing that regional patterns in climate and land cover will be important for scaling CH4 to the regional and national scale, we divided the US into clusters to evaluate representativeness using the first and second dimension from the MDS. Cluster analysis has been used to assess the spatial representativeness of network infrastructure and to suggest arrangements of study sites (Sulkava et al., 2011; Kumar et al., 2016). It is an objective method of producing meaningful, mutually exclusive groups based on similarities among entities (Balijepally et al., 2011). This approach is descriptive, a-theoretical, and non-inferential with sound mathematical support (Balijepally et al., 2011). Clustering outcomes are driven by large effect sizes or the accumulation of many smaller effects across features, and they are mostly unaffected by differences in covariance structure (Dalmaijer et al., 2020). Sufficient statistical power is achieved with relatively small samples (Dalmaijer et al., 2020), provided cluster separation is sufficient. Traditional notions about statistical power only partially apply to cluster analysis (Dalmaijer et al., 2020). Increasing the number of sample points above a sufficient sample size does not improve power, but effect size is important (Dalmaijer et al., 2020). Clustering is useful when large subgroup separation is expected and when MDS improves cluster separation (Dalmaijer et al., 2020).

We determined the optimal number of clusters using the library cluster and the function pam in R (Reynolds et al., 2006; Schubert and Rousseeuw, 2019, 2021). This approach uses the k-medoids algorithm, which partitions a data set into k groups or clusters and is a robust alternative to k-means clustering (Kaufman and Rousseeuw, 2009). The k-medoid algorithm is less sensitive to noise and outliers, compared to k-means, because it uses medoids as cluster centers. The k-medoids algorithm requires the user to specify k, the number of clusters to be generated. A useful approach to determine the optimal number of clusters is the silhouette method. We fit an increasing number of clusters from 2 to 50 to construct a silhouette plot and choose the number of clusters that maximized the average silhouette width (Fig. S2).

While useful, there are limitations to cluster analysis that can affect cluster patterns and the stability of clusters. The final cluster solution is dependent upon the clustering variables, the similarity/dissimilarity measure used, the clustering algorithm, and the data used to estimate clusters. Therefore, varying elements of clustering methods can lead to many alternative cluster solutions (Balijepally et al., 2011). Cluster solutions can also be produced in the absence of natural structure in the data, and there is no statistical basis to reject the null hypothesis that there are no natural groupings in the data (Balijepally et al., 2011). Cluster algorithms also cannot differentiate between relevant versus irrelevant variables. Therefore, only the variables expected to be influential should be used (Balijepally et al., 2011) and should emanate from past research or explicit theory and be consistent with the objectives of the study.

Due to the limitations of this approach, it is important to validate the cluster solution to ensure its meaningfulness and utility (Punj and Stewart, 1983; Balijepally et al., 2011). Consistency is established by checking the stability of cluster solutions obtained by using multiple algorithms (Punj and Stewart, 1983) or through splitting a sample, analyzing the cluster solutions for the two halves separately, and checking their consistency. After checking for reliability, the validity of a cluster solution is established through external validity and criterion-related validity. External validity ensures that clusters are representative of the actual population (Cook and Campbell, 1979) and can be verified by clustering on a hold-out sample using the same variables and assessing the similarity of the two solutions. This analysis was repeated five times to ensure that the 20 000 pixel subsample would produce similar results in the dimensions and clustering. For simplicity, we show the results of the first analysis, and a comparison of clustering methods and measures of stability are available in the Supplement.

To measure dissimilarity across the cluster once defined, each cluster was represented by one of the data points in the cluster named the cluster medoid. The medoid had the lowest average dissimilarity between it and all other objects in the cluster. The medoid can be considered a representative example of the members of that cluster. We calculated the dissimilarity between every location within the cluster to the medoid to create a measure of how different each location was from the medoid condition of each cluster. We utilized the pointDistance function in the raster package, which provided a unit-less relative measure of dissimilarity that was determined by measuring the difference between the first and second dimensions produced by the isoMDS of each point in a cluster to the dimensions of the medoid.

To extrapolate the cluster and dissimilarity layers across the entire US beyond the 20 000-pixel subsample and to show the predictive validity (Kerlinger, 1986), we employed the machine learning algorithm random forest (RF) with the package randomForest (Liaw and Wiener, 2002) to model the first and second dimensions using the land cover and climate layers as predictors. We then created a random forest model of the cluster layer using the first and second dimension as the explanatory variables. All models were then projected spatially to produce a spatially explicit cluster layer and a dissimilarity layer beyond the 20 000 sample points that were used in the MDS analysis. The RF algorithm was first introduced by Breiman (2001) and uses an ensemble of regression trees to predict target values. In RF, a series of bootstrapped data sets are used to generate independent regression trees; at each node, a random sample of predictor variables is selected for use. The RF prediction is the ensemble of multiple individual trees. We created 500 trees for each year and site, using 80 % of the data for model fitting and 20 % for model validation. The fit of each RF was evaluated with the out-of-bag mean square error (OOB MSE), and variable importance was computed as the amount of the prediction error increased when a particular predictor was permuted. Initially, 500 RF trees were generated. Overall model fit was evaluated with the average of the 500 OOB MSEs from the final model for each year and site, and variable importance was calculated as the average rank of each predictor variable for the 500 models. This approach allowed us to measure the importance of the original data on the first and second dimensions defined by the MDS and how the MDS leads to cluster and dissimilarity patterns. This step was essential to producing a spatially explicit cluster and dissimilarity layers for the entire US, since the MDS analysis limits the number of observations that can be analyzed. This is also important for evaluating the meaningfulness of the cluster by using the original variables used in the development of the distance matrix to predict clusters.

2.4 Measuring the landscape representativeness of research infrastructure

Representativeness studies discern when, where, and at what frequency networks are measuring ecological processes (Baldocchi et al., 2012; Jongman et al., 2017; Vaughan et al., 2001; Villarreal et al., 2018). To understand the representativeness of current CH4 infrastructure, we defined clusters (Sulkava et al., 2011) and measured the dissimilarity between each location in a cluster to the medoid. We extracted the cluster and dissimilarity for all active tower sites measuring CH4 that were distributed across the US and measured the tower cluster representativeness (TRcluster) as the percent overlap between the range of dissimilarity sampled by the infrastructure (rcluster) divided by the range of dissimilarity observed in the entire cluster (r; Eq. 1).

(1) TR cluster = r cluster r × 100

We recognize that it is essential to capture the distribution of dissimilarity across an entire cluster to upscale ecosystem measurements. We also report the sampling intensity of the major ecosystem types within the cluster and report the ecosystem representativeness (TRIGBP) by the International Geosphere–Biosphere Programme (IGBP) vegetation types of the towers (Eq. 2).

(2) TR IGBP = r IGBP r × 100

This approach allows the evaluation of representativeness that is not based on a specific research site, but on the dissimilarity of a location to other locations in the landscape, and we use the range to indicate a capacity to scale within a cluster which is based on both the effects of landscape dominant land cover, climate, and the specific ecosystems measured (IGBP).

3 Results

3.1 Measuring landscape dissimilarity across clusters within the US

Land cover, climate, and location were condensed down to two dissimilarity dimensions (Fig. 2a). Both climatic factors and location were the most important variables for determining dimensions and explained 99 % of the variance in dimensions (Fig. S1). Using the first and second dimensions, the US was divided into 10 clusters (Fig. S2) that were distributed across temperature and wetness gradients (Fig. 2; Table 2). The coldest zones were in Alaska and included clusters Na and Nb. Cool to temperate clusters in the midwestern and western US include NW, W, and NEa. Temperate clusters extend from the midwestern to the eastern US and include clusters NEb and Ea. Warm regions were distributed across clusters Eb, SW, and SE. Dry clusters (Na, SW, W, and Nb) were distributed across the western US and Alaska, and wet clusters (Ea, Eb, and SE) were in the south-eastern US and Hawaii. Individual clusters represented 7 %–16 % of the US each by area (Table 2) with cluster NW as the largest cluster in the Pacific Northwest, and the smallest cluster being cluster Nb in the northern half of Alaska.

Figure 2(a) Multidimensional scaling across the United States (US) produced 10 clusters using ecotype (Table 1), climate, and location (latitude/longitude). (b) Spatial distribution of the identified clusters.

Across all clusters, dissimilarity ranged from 0.01 to 0.33 (Fig. 3). The mean dissimilarity was 0.04, and most areas within a cluster were less than or equal to the mean. Southern Alaska (cluster Na), Hawaii (clusters SE and Eb), Florida (cluster SE), Puerto Rico (cluster SE), and the northeast (cluster NEa) had greater than average dissimilarity in their respective clusters.

Dominant landscape land cover types also varied across clusters, with forests, scrub, and herbaceous ecosystems dominating clusters (> 20 % coverage; Table 2). Although irrigated croplands did not have high coverage rates across any cluster, non-irrigated croplands had high coverage rates in NEb, Ea, and NEb. Wetlands did not have high coverage rates in any cluster.

Table 2The land cover and climate of the 10 clusters in the US. Crops were divided into irrigated (CropI) and non-irrigated (CropNI) and wetlands into emergent coastal (WetEC), emergent freshwater (WetEF), and freshwater forest (WetFF). Percent coverage (% Cov) is the percent area occupied by a cluster and R is the range in dissimilarity for each cluster.

Download Print Version | Download XLSX

Figure 3Cluster dissimilarity for the US. Inset: the distributions of dissimilarity across all clusters shown in a histogram, in which the line denotes the mean dissimilarity across all clusters.

3.2 Landscape representativeness of existing CH4 tower infrastructure

There were 70 active EC towers measuring CH4 distributed across forest (3 towers), grasslands (4 towers), shrublands (1 tower), agriculture (19 towers), wetlands (37 towers), barren (2 towers), and aquatic (4 towers) IGBP vegetation classes. Less than half of the active towers (43 %) were providing data to the community through AmeriFlux, limiting the development of CH4-derived products. For this reason, we will first focus this analysis on the active towers providing data to AmeriFlux. Although CH4 EC tower infrastructure was not a part of a single organized network designed to be representative of the climate, landscape, and dominant IGBP vegetation classes that exist within the US, EC tower infrastructure that was providing data to AmeriFlux was distributed across 8 of the 10 clusters (Table 3), with clusters NW and SE without any active towers providing data to the community. Tower representativeness (TRcluster) of clusters ranged from 0 %–88 %. The greatest TRcluster was for Eb and NEa, and the lowest TRcluster was for NW and SE, which had no towers. TRcluster was low (< 50 %) for most clusters, and high coverage was not associated with a higher frequency of towers. A high TRcluster was found in clusters where towers were dispersed across IGBP vegetation classes and where towers in wetlands, forests, or the arctic tundra (barren) were distributed across the observed range in the dissimilarity of clusters. Most clusters were substantially under-sampled (Table 3, Fig. 4) due to an insufficient number of towers measuring CH4 and poor distribution across the cluster.

The representativeness of IGBP vegetation types within clusters was poor for all vegetation types, excluding forests in the NEa. TRIGBP ranged from 0 %–79 %, and wetlands were the only IGBP class to be sampled across eight clusters. Ideally, IGBP classes should be distributed both within and across clusters where the classes exist. There was not a single cluster with towers in all of the IGBP classes (forest, scrub, aquatic ecosystems, crops, wetlands, barren tundra, and grasslands) that are found within that cluster.

Table 3The total number of eddy covariance (EC) towers measuring CH4 and providing data to AmeriFlux. The tower frequency by dominant landscape type, the total cluster representativeness (TRcluster), and cluster representativeness by major ecosystem types are shown (TRIGBP). For TRcluster and TRIGBP values of 0.01 were assigned where a single tower was present.

Download Print Version | Download XLSX

Figure 4The range in dissimilarity for clusters (black bar), active CH4 towers providing CH4 data to AmeriFlux (cyan), all active CH4 towers (magenta), and for NEON towers (blue). The black lines show the range in dissimilarity observed for a cluster and greater overlap between the cluster range and the tower range is important for landscape representativeness.


Table 4The TRcluster for CH4 towers that are active and providing data to AmeriFlux, the TRcluster for all active CH4 towers, and the TRcluster for all active towers in addition to NEON towers.

Download Print Version | Download XLSX

There were important gains in the TRcluster when considering all CH4 towers regardless of if they were providing data to AmeriFlux (Table 4 and Fig. 4). The clusters with substantial gains in representativeness (> 10 %) include Na, NEb, Ea, and the SE. The TRcluster of the NW, Ea, SW, W, and the SE would be further enhanced by more than 10 % with the addition of CH4 instrumentation at NEON tower sites.

4 Discussion

To determine key regions where biogenic CH4 infrastructure is needed within the US, we identified gaps in active research infrastructure. We found that there is an insufficient number of towers measuring CH4, and the distribution of these sites across the range in dissimilarity observed is poor for all clusters. Current EC towers measuring CH4 are in ecosystems known to be sources of CH4. This is extremely limiting when trying to upscale CH4 fluxes because it leads to a serious bias towards CH4 emissions in model results and constrains our capacity to appropriately model ecosystems that are CH4 sinks. In this analysis, we include NEON towers because they are purposefully distributed across climate zones and ecosystem types, they provide consistent and standardized measurements, existing infrastructure at these sites could be quickly adapted to measure CH4, and all data are publicly available. We understand that for PI-managed infrastructure the placement of towers is driven by the scientific question being asked and research funding priorities (Papale et al., 2015; Mahecha et al., 2017; Villarreal et al., 2018; Knox et al., 2019), but as the number of towers measuring CH4 fluxes continues to grow, consideration for key underrepresented regions where towers are needed or where more efforts are needed for existing but nonreporting towers to contribute to AmeriFlux is of utmost importance. Making all data available must become the standard of the trace gas flux and biogeochemistry communities. Notable infrastructure gaps were in clusters Na, NW, SW, W, Nb, and the SE, and all clusters require a greater representation of IGBP vegetation types. Our analysis shows that the Na, W, and Nb clusters are the most poorly represented regions, corresponding to Alaska (Na and Nb) and the Rocky Mountains (W), where large elevational changes in the landscape are inherently difficult to capture.

One reason for gaps in CH4 flux tower infrastructure may be the lag in technological capability behind that of CO2 flux measurements. Methane gas analyzers with sufficient measurement frequency for EC were not common before the late 1990s and early 2000s (Shurpali et al., 1993; Billesbach et al., 1998; Rinne et al., 2007), and the number of commercial options has expanded only more recently (Peltola et al., 2013; Nemitz et al., 2018; Burba et al., 2019; Burba, 2021). Therefore, as the flux tower infrastructure has expanded to measure CH4, decisions on the locations of measurement sites have largely been tied to CO2 and water vapor exchange research (Baldocchi, 2014) and to the availability of suitable infrastructure (McDermitt et al., 2011), and not necessarily to address CH4 hypotheses. In addition to technological limitations, the environments where we expect CH4 fluxes to be highest complicate considerations for where best to place instrumentation. Large sources of natural biogenic CH4 can sometimes originate from small, heterogeneous components within a landscape, such as patchy wetlands within an otherwise upland forested region, causing the area to be a net source of CH4 (Desai et al., 2015). In contrast, some systems covering large areas that are known to be important CH4 sources, such as arctic tundra ecosystems and shallow lakes (Wik et al., 2016; Elder et al., 2020), are simply too remote and difficult to instrument. When they are instrumented, towers are often clustered together regionally, resulting in incremental changes in landscape representativeness. A non-negligible portion of the existing CH4 measurements, including both towers and chambers, are not placed where CH4 sources or sinks are but where the grid power is available to run such measurements. The likely incomplete quantification of CH4 fluxes within heterogeneous sites and the measurement of CH4 fluxes at sites that were established to measure CO2 and energy fluxes together introduce an inherent source of site-level bias in existing data and our analyses. Hence, we interpret our results as a best-case scenario, as this bias likely would reduce even further our reported degree of representativeness.

Gaps in our US infrastructure and current capability to measure CH4 were most noted when considering only the AmeriFlux sites that provide CH4 data. When evaluating all sites with CH4 infrastructure with the addition of the measurement capability from NEON sites, there were great improvements in landscape representativeness. Still, the largest gaps in infrastructure capability to measure expected CH4 sources were from aquatic sites. These gaps in representation have been noted in other investigations of CH4 flux and budget studies, as a part of larger global CH4 analyses (Saunois et al., 2020) and FLUXNET CH4 flux syntheses (Knox et al., 2019; Delwiche et al., 2021). In fact, the call for more measurements of CH4 from natural sites is not new (Matthews and Fung, 1987; Bartlett and Harriss, 1993; Dlugokencky et al., 2011; Nisbet et al., 2014) and has been advocated as necessary to reduce the uncertainty in CH4 budget estimates from natural ecosystems (Peltola et al., 2019), which is among the largest uncertainty in the global CH4 budget (Saunois et al., 2020). Even areas that have been traditionally thought to have negligible CH4 emission or consumption rates should be monitored because their contribution to CH4 budgets may be significant when considering their large spatial extent. There is also a strong need for a continental CH4 observatory to aid in reducing these uncertainties in the natural CH4 sources and sinks.

A large source of uncertainty in scaling bottom-up CH4 estimates are in the current land use classification (LUC) products (Kirschke et al., 2013; Knox et al., 2019; Saunois et al., 2020), which are not designed to estimate the potential CH4 source/sink status, particularly from aquatic, wetland, and agricultural land cover. Aquatic ecosystems contribute significantly to global CH4 emissions, with emissions increasing from natural to impacted aquatic ecosystems and from coastal to freshwater ecosystems (Rosentreter et al., 2021). Specific ecosystems within the landscape can contribute significantly to landscape-level and regional CH4 source/sink estimates. Aquatic emissions are likely to change in the future due to an increase in urbanization, eutrophication, and positive climate feedbacks (IPCC, 2021). Yet current wetland classifications from land use data products are not suitable to capture these potential changes, or the potential feedbacks they may have on CH4 processes. Wetland classifications are often generalized too broadly in current LUC schemas to accurately scale and predict CH4 flux rates and processes. Small changes in the delineation or characterization of LUC can result in changing the source/sink status of whole regions (Kirschke et al., 2013; Barkley et al., 2017; Knox et al., 2019). For wetlands these include (i) delineation of wetland area, the largest natural CH4 source, especially in regions like Alaska and Florida, (ii) conflation of fluxes from wetlands and fresh waters leading to double counting (Thornton et al., 2016), and (iii) classification of saturated soils as non-wetland, possibly missing strong CH4 emission potential. For agricultural lands, we must also consider (iv) deforestation for agricultural use, which reduces the soil CH4 sink potential (Robertson et al., 2000), or (v) accurate representation of agricultural land CH4 potential when land use includes a complex mixture of ruminants feedlots, manure, and pastures (Lassey, 2008). These potential large sources of uncertainties in biogenic CH4 flux estimates cannot be addressed with the land cover maps currently used to scale CH4 fluxes and the existing distribution of CH4 observation sites (Rosentreter et al., 2021). Hence, if we are to build a US CH4 budget using a scaled-up land use classification scheme (as is done for CO2), we need both better representation of CH4 measurement sites and better identification and quantification of the CH4 source/sink potential of the land use classes themselves, i.e., specific development of land use classes based on CH4 potential.

Ideally, CH4 measurement infrastructure should have representation of all IGBP vegetation classes within and across clusters, where appropriate, and address the scale of spatial heterogeneity that reduces uncertainty in a national CH4 budget with confidence limits that can inform both research objectives and mitigation policy. Thus, the incorporation of representative CH4 sources and sink strength is essential to develop national CH4 budgets. Neglecting sinks would further bias models that suggest sources occur where we are confident they do not. Advancing research and our process-level understanding of biogenic CH4, we need to determine the measurement scales to assess the degree of spatial heterogeneity required to reduce uncertainty within and among sites. One means to address the within-site scale of spatial uncertainty is from automated chamber measurements within flux tower footprints, such as that found in soils, or first- and second-order streams. This would also allow the scientific community to determine the within site CH4 source/sink strength from local (chamber; < 1 m2), ecosystem (EC flux tower;  1 km2), and landscape scales (tower concentrations;  100 s km2). At even a larger scale, airborne observations of atmospheric CH4 concentrations can be used to estimate boundary-level surface–atmosphere CH4 fluxes and potentially provide greater spatial coverage than towers (Chang et al., 2014; Zona et al., 2016) and provide a mechanistic link between tower-based and satellite-derived CH4 estimates.

The rate of global climate change re-enforces the urgency to establish a continental-scale CH4 observatory network that can enable the first national CH4 budget. As it stands, we currently do not know the scale, location, or the magnitude of site-based biogenic CH4 source/sinks to estimate a national budget. For example, we lack the quantitative information about specific processes (particularly those that are stochastic, e.g., temperature sensitivity, susceptibility to drought and flooding, tipping points) from representative ecosystems that would scale and inform a national CH4 budget. In addition to the current uncertainty in basic ecosystem-level CH4 processes and the way they spatially scale, the backdrop of climate change is also changing the rates of CH4 production and consumption, as well as the CH4 transport pathways. For example, arctic regions are warming faster than most other regions of the world (Serreze and Barry, 2011), turning permafrost into wetlands and changing traditional CH4 sinks to sources on short timescales (Chadburn et al., 2017; Schaefer, 2019; Yumashev et al., 2019). In temperate areas, higher climate-change-induced variability in precipitation (e.g., higher moisture of upland forested soils, prolonged droughts) results in a reduction of soil CH4 uptake and a reduced global CH4 sink (Ni and Groffman, 2018). Sea-level rise, which leads to the inundation of coastal regions turning previously dry upland environments into saturated, anoxic areas, can in some cases increase CH4 production and emission rates (Lu et al., 2018). Hence, we do not have a baseline US CH4 budget to establish a starting point now and to compare to in the future, and as a baseline to estimate the efficacy of any mitigation decision (policy) made today. As scientists, we are often asked what the most likely future state of an ecological system is and what the most likely state of a system is, given a decision or action is made today. The current state of CH4 research and its ability to inform these questions are still nascent.

5 Conclusions

We used landscape dissimilarity to assess gaps in current CH4 infrastructure at the landscape scale in the US. Evaluating the strengths and limitations of existing measurement infrastructure is critical for strategic augmentation to provide the most valuable information toward reducing uncertainties in future large-scale budget estimations. This analysis complements previous studies based on climatic or vegetation characteristics (Hargrove and Hoffman, 2003; Yang et al., 2008; Villarreal et al., 2018) and identifies regions within the US where gaps are limiting the development of upscaling techniques. To accurately understand the impact of climate and land cover change on biogenic CH4 emissions, we need a long-term, calibrated, and strategic continental-scale CH4 observatory network. Current gaps in existing measurement infrastructure limit our ability to capture the spatial and temporal variations of biogenic CH4 fluxes and therefore limit our ability to predict future CH4 emissions. Maps of potential CH4 emissions require land cover classification targeted at land cover types like wetlands that are important sources of CH4 to the atmosphere. Aquatic ecosystems like streams and lakes as well as coastal ecosystems are significant and variable sources of CH4 not well studied on a long-term basis. Through our analysis using climate, land cover, and location variables, we have identified priority areas to enhance research infrastructure to provide a more complete understanding of the CH4 flux potential of ecosystem types in the US. For EC tower locations, dissimilarity coverage was lacking for clusters Na, W, and Nb, and currently clusters Na, W, Eb, and Nb are substantially undersampled. All aquatic sites were undersampled within each cluster. An enhanced network would allow for us to monitor both the response of CH4 fluxes to climate and land use change as well as to assess the impact of future policy and mitigation strategies.

Code and data availability

The data products produced are available on the Knowledge Network for Biocomplexity (; Malone, 2021).


The supplement related to this article is available online at:

Author contributions

All authors contributed to conceptualization; SLM, KAA, RC, ARC, JPG, and RKV designed the manuscript; SLM, KAA, JPG, GS, and RKV wrote parts of the manuscript; RKV, SLM, and JPG supervised the project/manuscript; SLM, KAA, and YO contributed to data curation and formal analysis, and all authors performed critical reviews.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors appreciate the substantial contributions of reviewers. The authors would also like to acknowledge the researchers and support staff for their contributions to discussions that led to the idea for this article: Melissa Genazzio, Amy Lafreniere, Kim Nitschke, Lori Bruhwiler, Julia Bryce, Patrick Crill, Amarnath Gupta, Ilya Zaslavsky, Ilkay Altintas, Stephen Hale, Mike Stewart, Michael Thomson, and Mark Milutinovich. Henry W. Loescher acknowledges the National Science Foundation (NSF) for ongoing support. NEON is a project sponsored by the NSF and managed under a cooperative support agreement (EF-1029808) to Battelle. Ruth K. Varner and Alexandra R. Contosta acknowledge UNH's Collaborative Research Excellence (CoRE) grant. Sparkle L. Malone acknowledges support provided by NSF grant no. 2047687. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of our sponsoring agencies. This is contribution no. 1432 from the Institute of Environment at Florida International University.

Financial support

This research has been supported by the National Science Foundation (NSF) through a cooperative support agreement (EF-1029808) to Battelle and through grant number 2047687. Support was also provided by the University of New Hampshire’s Collaborative Research Excellence (CoRE) grant.

Review statement

This paper was edited by Ben Bond-Lamberty and reviewed by Jitendra Kumar and one anonymous referee.


Ahmad, A. and Dey, L.: A k-mean clustering algorithm for mixed numeric and categorical data, Data Knowl. Eng., 63, 503–527,, 2007.  

Arndt, K. A., Oechel, W. C., Goodrich, J. P., Bailey, B. A., Kalhori, A., Hashemi, J., Sweeney, C., and Zona, D.: Sensitivity of methane emissions to later soil freezing in arctic tundra ecosystems, J. Geophys. Res.-Biogeosci., 124, 2595–2609,, 2019. 

Baldocchi, D.: Measuring fluxes of trace gases and energy between ecosystems and the atmosphere – the state and future of the eddy covariance method, Glob. Change Biol., 20, 3600–3609,, 2014. 

Baldocchi, D., Reichstein, M., and Papale, D.: The role of trace gas flux networks in the biogeosciences, Eos Trans. Am. Geophys. Union,, 2012. 

Balijepally, V., Mangalaraj, G., and Iyengar, K.: Are We Wielding this Hammer Correctly? A Reflective Review of the Application of Cluster Analysis in Information Systems Research, J. Assoc. Inf. Syst., 12, 375–413,, 2011. 

Barkley, Z. R., Lauvaux, T., Davis, K. J., Deng, A., Miles, N. L., Richardson, S. J., Cao, Y., Sweeney, C., Karion, A., Smith, M., Kort, E. A., Schwietzke, S., Murphy, T., Cervone, G., Martins, D., and Maasakkers, J. D.: Quantifying methane emissions from natural gas production in north-eastern Pennsylvania, Atmos. Chem. Phys., 17, 13941–13966,, 2017. 

Bartlett, K. B. and Harriss, R. C.: Review and assessment of methane emissions from wetlands, Chemosphere, 26, 261–320,, 1993. 

Bessembinder, J., Overbeek, B., and Siegmund, P.: Climate normals and climate change: how to communicate these together?, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-4032,, 2021. 

Billesbach, D. P., Kim, J., Clement, R. J., Verma, S. B., and Ullman, F. G.: An Intercomparison of Two Tunable Diode Laser Spectrometers Used for Eddy Correlation Measurements of Methane Flux in a Prairie Wetland, J. Atmos. Ocean. Technol., 15, 197–206,<0197:aiottd>;2, 1998. 

Boryan, C., Yang, Z., Mueller, R., and Craig, M.: Monitoring US agriculture: the US Department of Agriculture, National Agricultural Statistics Service, Cropland Data Layer Program, Geocarto Int., 26, 341–358,, 2011. 

Breiman, L.: Random Forests, Mach. Learn., 45, 5–32,, 2001. 

Bruhwiler, L., Parmentier, F.-J. W., Crill, P., Leonard, M., and Palmer, P. I.: The Arctic Carbon Cycle and Its Response to Changing Climate, 7, 14–34,, 2021. 

Burba, G.: 9 – Atmospheric flux measurements, in: Advances in Spectroscopic Monitoring of the Atmosphere, edited by: Chen, W., Venables, D. S., and Sigrist, M. W., Elsevier, 443–520,, 2021. 

Burba, G., Anderson, T., and Komissarov, A.: Accounting for spectroscopic effects in laser-based open-path eddy covariance flux measurements, Glob. Change Biol., 25, 2189–2202,, 2019. 

Chadburn, S. E., Burke, E. J., Cox, P. M., Friedlingstein, P., Hugelius, G., and Westermann, S.: An observation-based constraint on permafrost loss as a function of global warming, Nat. Clim. Change, 7, 340–344,, 2017. 

Chang, R. Y.-W., Miller, C. E., Dinardo, S. J., Karion, A., Sweeney, C., Daube, B. C., Henderson, J. M., Mountain, M. E., Eluszkiewicz, J., Miller, J. B., Bruhwiler, L. M. P., and Wofsy, S. C.: Methane emissions from Alaska in 2012 from CARVE airborne observations, P. Natl. Acad. Sci. USA, 111, 16694–16699,, 2014. 

Chen, B., Coops, N. C., Fu, D., Margolis, H. A., Amiro, B. D., Barr, A. G., Black, T. A., Arain, M. A., Bourque, C. P.-A., Flanagan, L. B., Lafleur, P. M., McCaughey, J. H., and Wofsy, S. C.: Assessing eddy-covariance flux tower location bias across the Fluxnet-Canada Research Network based on remote sensing and footprint modelling, Agr. Forest Meteorol., 151, 87–100,, 2011. 

Chu, H., Luo, X., Ouyang, Z., Chan, W. S., Dengel, S., Biraud, S. C., Torn, M. S., Metzger, S., Kumar, J., Arain, M. A., Arkebauer, T. J., Baldocchi, D., Bernacchi, C., Billesbach, D., Black, T. A., Blanken, P. D., Bohrer, G., Bracho, R., Brown, S., Brunsell, N. A., Chen, J., Chen, X., Clark, K., Desai, A. R., Duman, T., Durden, D., Fares, S., Forbrich, I., Gamon, J. A., Gough, C. M., Griffis, T., Helbig, M., Hollinger, D., Humphreys, E., Ikawa, H., Iwata, H., Ju, Y., Knowles, J. F., Knox, S. H., Kobayashi, H., Kolb, T., Law, B., Lee, X., Litvak, M., Liu, H., Munger, J. W., Noormets, A., Novick, K., Oberbauer, S. F., Oechel, W., Oikawa, P., Papuga, S. A., Pendall, E., Prajapati, P., Prueger, J., Quinton, W. L., Richardson, A. D., Russell, E. S., Scott, R. L., Starr, G., Staebler, R., Stoy, P. C., Stuart-Haëntjens, E., Sonnentag, O., Sullivan, R. C., Suyker, A., Ueyama, M., Vargas, R., Wood, J. D., and Zona, D.: Representativeness of Eddy-Covariance flux footprints for areas surrounding AmeriFlux sites, Agr. Forest Meteorol., 301/302, 108350,, 2021. 

Ciais, P., Dolman, A. J., Bombelli, A., Duren, R., Peregon, A., Rayner, P. J., Miller, C., Gobron, N., Kinderman, G., Marland, G., Gruber, N., Chevallier, F., Andres, R. J., Balsamo, G., Bopp, L., Bréon, F.-M., Broquet, G., Dargaville, R., Battin, T. J., Borges, A., Bovensmann, H., Buchwitz, M., Butler, J., Canadell, J. G., Cook, R. B., DeFries, R., Engelen, R., Gurney, K. R., Heinze, C., Heimann, M., Held, A., Henry, M., Law, B., Luyssaert, S., Miller, J., Moriyama, T., Moulin, C., Myneni, R. B., Nussli, C., Obersteiner, M., Ojima, D., Pan, Y., Paris, J.-D., Piao, S. L., Poulter, B., Plummer, S., Quegan, S., Raymond, P., Reichstein, M., Rivier, L., Sabine, C., Schimel, D., Tarasova, O., Valentini, R., Wang, R., van der Werf, G., Wickland, D., Williams, M., and Zehner, C.: Current systematic carbon-cycle observations and the need for implementing a policy-relevant carbon observing system, Biogeosciences, 11, 3547–3602,, 2014. 

Cook, T. D. and Campbell, D. T.: Quasi-experimentation: Design and Analysis Issues for Field Settings, Rand McNally College, 405 pp., ISBN 9780528686948, 1979. 

Cox, M. A. A. and Cox, T. F.: Multidimensional Scaling, in: Handbook of Data Visualization, edited by: Chen, C.-H., Härdle, W., and Unwin, A., Springer Berlin Heidelberg, Berlin, Heidelberg, 315–347,, 2008. 

Dalmaijer, E. S., Nord, C. L., and Astle, D. E.: Statistical power for cluster analysis, arXiv [stat.ML], arXiv,, 2020. 

Delwiche, K. B., Knox, S. H., Malhotra, A., Fluet-Chouinard, E., McNicol, G., Feron, S., Ouyang, Z., Papale, D., Trotta, C., Canfora, E., Cheah, Y.-W., Christianson, D., Alberto, Ma. C. R., Alekseychik, P., Aurela, M., Baldocchi, D., Bansal, S., Billesbach, D. P., Bohrer, G., Bracho, R., Buchmann, N., Campbell, D. I., Celis, G., Chen, J., Chen, W., Chu, H., Dalmagro, H. J., Dengel, S., Desai, A. R., Detto, M., Dolman, H., Eichelmann, E., Euskirchen, E., Famulari, D., Fuchs, K., Goeckede, M., Gogo, S., Gondwe, M. J., Goodrich, J. P., Gottschalk, P., Graham, S. L., Heimann, M., Helbig, M., Helfter, C., Hemes, K. S., Hirano, T., Hollinger, D., Hörtnagl, L., Iwata, H., Jacotot, A., Jurasinski, G., Kang, M., Kasak, K., King, J., Klatt, J., Koebsch, F., Krauss, K. W., Lai, D. Y. F., Lohila, A., Mammarella, I., Belelli Marchesini, L., Manca, G., Matthes, J. H., Maximov, T., Merbold, L., Mitra, B., Morin, T. H., Nemitz, E., Nilsson, M. B., Niu, S., Oechel, W. C., Oikawa, P. Y., Ono, K., Peichl, M., Peltola, O., Reba, M. L., Richardson, A. D., Riley, W., Runkle, B. R. K., Ryu, Y., Sachs, T., Sakabe, A., Sanchez, C. R., Schuur, E. A., Schäfer, K. V. R., Sonnentag, O., Sparks, J. P., Stuart-Haëntjens, E., Sturtevant, C., Sullivan, R. C., Szutu, D. J., Thom, J. E., Torn, M. S., Tuittila, E.-S., Turner, J., Ueyama, M., Valach, A. C., Vargas, R., Varlagin, A., Vazquez-Lule, A., Verfaillie, J. G., Vesala, T., Vourlitis, G. L., Ward, E. J., Wille, C., Wohlfahrt, G., Wong, G. X., Zhang, Z., Zona, D., Windham-Myers, L., Poulter, B., and Jackson, R. B.: FLUXNET-CH4: a global, multi-ecosystem dataset and analysis of methane seasonality from freshwater wetlands, Earth Syst. Sci. Data, 13, 3607–3689,, 2021. 

Desai, A. R.: Climatic and phenological controls on coherent regional interannual variability of carbon dioxide flux in a heterogeneous landscape, J. Geophys. Res., 115, G00J02,, 2010. 

Desai, A. R., Xu, K., Tian, H., Weishampel, P., Thom, J., Baumann, D., Andrews, A. E., Cook, B. D., King, J. Y., and Kolka, R.: Landscape-level terrestrial methane flux observed from a very tall tower, Agr. Forest Meteorol., 201, 61–75,, 2015. 

Dlugokencky, E.: Trends in Atmospheric Methane Global CH4 Monthly Means, NOAA, (last access: 5 January 2022), 2021. 

Dlugokencky, E. J., Nisbet, E. G., Fisher, R., and Lowry, D.: Global atmospheric methane: budget, changes and dangers, Philos. Trans. A Math. Phys. Eng. Sci., 369, 2058–2072,, 2011. 

Elder, C. D., Thompson, D. R., Thorpe, A. K., Hanke, P., Walter Anthony, K. M., and Miller, C. E.: Airborne mapping reveals emergent power law of arctic methane emissions, Geophys. Res. Lett., 47, e2019GL085707,, 2020. 

Gower, J. C.: A General Coefficient of Similarity and Some of Its Properties, Biometrics, 27, 857–871,, 1971. 

Groffman, P. M., Hardy, J. P., Driscoll, C. T., and Fahey, T. J.: Snow depth, soil freezing, and fluxes of carbon dioxide, nitrous oxide and methane in a northern hardwood forest, Glob. Change Biol., 12, 1748–1760,, 2006. 

Hargrove, W. W. and Hoffman, F. M.: New analysis reveals representativeness of the AmeriFlux network, Eos Trans. Amer. Geophys. Union, 84, 529–544, 2003. 

Harikumar, S. and Pv, S.: K-Medoid Clustering for Heterogeneous DataSets, Procedia Comput. Sci., 70, 226–237,, 2015. 

He, H., Zhang, L., Gao, Y., Ren, X., Zhang, L., Yu, G., and Wang, S.: Regional representativeness assessment and improvement of eddy flux observations in China, Sci. Total Environ., 502, 688–698,, 2015. 

Hijmans, R. J.: Geographic Data Analysis and Modeling [R package raster version 3.4-13], Comprehensive R Archive Network (CRAN), last access: 12 August 2021. 

Hoffman, F. M., Kumar, J., Mills, R. T., and Hargrove, W. W.: Representativeness-based sampling network design for the State of Alaska, 28, 1567–1586,, 2013. 

Huang, Z.: Clustering large data sets with mixed numeric and categorical values, in: Proceedings of the 1st pacific-asia conference on knowledge discovery and data mining (PAKDD), PAKDD, Singapore, 21–34, 23–24 February, (last access: 9 May 2022), 1997. 

Hutchins, D. A., Jansson, J. K., Remais, J. V., Rich, V. I., Singh, B. K., and Trivedi, P.: Climate change microbiology - problems and perspectives, Nat. Rev. Microbiol., 17, 391–396,, 2019. 

IPCC: The physical science basis, Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change, USA, Cambridge University Press, 1535 pp., (last access: 9 May 2022), 2013. 

IPCC: Climate Change 2021: The Physical Science Basis, Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press,, 2021. 

Ippoliti, C., Candeloro, L., Gilbert, M., Goffredo, M., Mancini, G., Curci, G., Falasca, S., Tora, S., Di Lorenzo, A., Quaglia, M., and Conte, A.: Defining ecological regions in Italy based on a multivariate clustering approach: A first step towards a targeted vector borne disease surveillance, PLoS One, 14, e0219072,, 2019. 

Jin, S., Homer, C., Yang, L., Danielson, P., Dewitz, J., Li, C., Zhu, Z., Xian, G., and Howard, D.: Overall Methodology Design for the United States National Land Cover Database 2016 Products, Remote Sens., 11, 2971,, 2019. 

Jongman, R. H. G., Skidmore, A. K., Mücher, C. A. S., Bunce, R. G. H., and Metzger, M. J.: Global terrestrial ecosystem observations: why, where, what and how?, in: The GEO handbook on biodiversity observation networks, Springer, Cham, 19–38, ISBN 978-3-319-27288-7, 2017. 

Jung, M., Reichstein, M., Margolis, H. A., Cescatti, A., Richardson, A. D., Altaf Arain, M., Arneth, A., Bernhofer, C., Bonal, D., Chen, J., Gianelle, D., Gobron, N., Kiely, G., Kutsch, W., Lasslop, G., Law, B. E., Lindroth, A., Merbold, L., Montagnani, L., Moors, E. J., Papale, D., Sottocornola, M., Vaccari, F., and Williams, C.: Global patterns of land-atmosphere fluxes of carbon dioxide, latent heat, and sensible heat derived from eddy covariance, satellite, and meteorological observations, J. Geophys. Res., 116, G00J07,, 2011. 

Kaufman, L. and Rousseeuw, P. J.: Finding Groups in Data: An Introduction to Cluster Analysis, John Wiley & Sons, 342 pp.,, 2009. 

Kerlinger, F. N.: Foundations of Behavioral Research, Holt, Rinehart Winston, New York, NY, ISBN 9780030417610, 1986. 

Kirschke, S., Bousquet, P., Ciais, P., Saunois, M., Canadell, J. G., Dlugokencky, E. J., Bergamaschi, P., Bergmann, D., Blake, D. R., Bruhwiler, L., Cameron-Smith, P., Castaldi, S., Chevallier, F., Feng, L., Fraser, A., Heimann, M., Hodson, E. L., Houweling, S., Josse, B., Fraser, P. J., Krummel, P. B., Lamarque, J.-F., Langenfelds, R. L., Le Quéré, C., Naik, V., O'Doherty, S., Palmer, P. I., Pison, I., Plummer, D., Poulter, B., Prinn, R. G., Rigby, M., Ringeval, B., Santini, M., Schmidt, M., Shindell, D. T., Simpson, I. J., Spahni, R., Steele, L. P., Strode, S. A., Sudo, K., Szopa, S., van der Werf, G. R., Voulgarakis, A., van Weele, M., Weiss, R. F., Williams, J. E., and Zeng, G.: Three decades of global methane sources and sinks, Nat. Geosci., 6, 813–823,, 2013. 

Knox, S. H., Jackson, R. B., Poulter, B., McNicol, G., Fluet-Chouinard, E., Zhang, Z., Hugelius, G., Bousquet, P., Canadell, J. G., Saunois, M., Papale, D., Chu, H., Keenan, T. F., Baldocchi, D., Torn, M. S., Mammarella, I., Trotta, C., Aurela, M., Bohrer, G., Campbell, D. I., Cescatti, A., Chamberlain, S., Chen, J., Chen, W., Dengel, S., Desai, A. R., Euskirchen, E., Friborg, T., Gasbarra, D., Goded, I., Goeckede, M., Heimann, M., Helbig, M., Hirano, T., Hollinger, D. Y., Iwata, H., Kang, M., Klatt, J., Krauss, K. W., Kutzbach, L., Lohila, A., Mitra, B., Morin, T. H., Nilsson, M. B., Niu, S., Noormets, A., Oechel, W. C., Peichl, M., Peltola, O., Reba, M. L., Richardson, A. D., Runkle, B. R. K., Ryu, Y., Sachs, T., Schäfer, K. V. R., Schmid, H. P., Shurpali, N., Sonnentag, O., Tang, A. C. I., Ueyama, M., Vargas, R., Vesala, T., Ward, E. J., Windham-Myers, L., Wohlfahrt, G., and Zona, D: FLUXNET-CH 4 Synthesis Activity: Objectives, Observations, and Future Directions, Bull. Am. Meteorol. Soc., 100, 2607–2632, 2019. 

Kumar, J., Hoffman, F. M., Hargrove, W. W., and Collier, N.: Understanding the representativeness of FLUXNET for upscaling carbon flux from eddy covariance measurements, Earth Syst. Sci. Data Discuss. [preprint],, 2016. 

Lassey, K. R.: Livestock methane emission and its perspective in the global methane cycle, Aust. J. Exp. Agr., 48, 114–118,, 2008. 

Le Quéré, C., Peters, G. P., Friedlingstein, P., Andrew, R. M., Canadell, J. G., Davis, S. J., Jackson, R. B., and Jones, M. W.: Fossil CO2 emissions in the post-COVID-19 era, Nat. Clim. Change, 11, 197–199,, 2021. 

Liaw, A. and Wiener, M.: Classification and regression by randomForest, R news, 2, 18–22, 2002. 

Lovett, G. M., Burns, D. A., Driscoll, C. T., Jenkins, J. C., Mitchell, M. J., Rustad, L., Shanley, J. B., Likens, G. E., and Haeuber, R.: Who needs environmental monitoring?, Front. Ecol. Environ., 5, 253–260,[253:WNEM]2.0.CO;2, 2007. 

Lu, X., Zhou, Y., Zhuang, Q., Prigent, C., Liu, Y., and Teuling, A.: Increasing methane emissions from natural land ecosystems due to sea-level rise, J. Geophys. Res.-Biogeo., 123, 1756–1768,, 2018. 

Mahecha, M. D., Gans, F., Sippel, S., Donges, J. F., Kaminski, T., Metzger, S., Migliavacca, M., Papale, D., Rammig, A., and Zscheischler, J.: Detecting impacts of extreme events with ecological in situ monitoring networks, Biogeosciences, 14, 4255–4277,, 2017. 

Malone, S.: Gaps in Network Infrastructure limit our understanding of biogenic methane emissions in the United States, knb [data set],, 2021. 

Matthews, E. and Fung, I.: Methane emission from natural wetlands: Global distribution, area, and environmental characteristics of sources, Global Biogeochem. Cy., 1, 61–86,, 1987. 

McDermitt, D., Burba, G., Xu, L., Anderson, T., Komissarov, A., Riensche, B., Schedlbauer, J., Starr, G., Zona, D., Oechel, W., Oberbauer, S., and Hastings, S.: A new low-power, open-path instrument for measuring methane flux by eddy covariance, Appl. Phys. B, 102, 391–405,, 2011. 

Michalak, A. M., Jackson, R., Marland, G., and Sabine, C.: A U.S. Carbon Cycle Science Plan:, First Meeting of the Carbon Cycle Science Working Group, Eos Transactions American Geophysical Union, Washington, D. C, 102–103,, 2009. 

Nemitz, E., Mammarella, I., Ibrom, A., Aurela, M., Burba, G. G., Dengel, S., Gielen, B., Grelle, A., Heinesch, B., Herbst, M., Hörtnagl, L., Klemedtsson, L., Lindroth, A., Lohila, A., McDermitt, D. K., Meier, P., Merbold, L., Nelson, D., Nicolini, G., Nilsson, M. B., Peltola, O., Rinne, J., and Zahniser, M.: Standardisation of eddy-covariance flux measurements of methane and nitrous oxide, Int. Agrophys., 32, 517–549,, 2018. 

Ni, X. and Groffman, P. M.: Declines in methane uptake in forest soils, P. Natl. Acad. Sci. USA, 115, 8587–8590,, 2018. 

Nisbet, E. G., Dlugokencky, E. J., and Bousquet, P.: Methane on the Rise – Again, Science, 343, 493–495,, 2014. 

Nisbet, E. G., Manning, M. R., Dlugokencky, E. J., Fisher, R. E., Lowry, D., Michel, S. E., Myhre, C. L., Platt, S. M., Allen, G., Bousquet, P., Brownlow, R., Cain, M., France, J. L., Hermansen, O., Hossaini, R., Jones, A. E., Levin, I., Manning, A. C., Myhre, G., Pyle, J. A., Vaughn, B. H., Warwick, N. J., and White, J. W. C.: Very strong atmospheric methane growth in the 4 years 2014–2017: Implications for the Paris agreement, Global Biogeochem. Cy., 33, 318–342,, 2019. 

Novick, K. A., Biederman, J. A., Desai, A. R., Litvak, M. E., Moore, D. J. P., Scott, R. L., and Torn, M. S.: The AmeriFlux network: A coalition of the willing, Agr. Forest Meteorol., 249, 444–456,, 2018. 

Oh, Y., Zhuang, Q., Liu, L., Welp, L. R., Lau, M. C. Y., Onstott, T. C., Medvigy, D., Bruhwiler, L., Dlugokencky, E. J., Hugelius, G., D'Imperio, L., and Elberling, B.: Reduced net methane emissions due to microbial methane oxidation in a warmer Arctic, Nat. Clim. Change, 10, 317–321,, 2020. 

Papale, D., Black, T. A., Carvalhais, N., Cescatti, A., Chen, J., Jung, M., Kiely, G., Lasslop, G., Mahecha, M. D., Margolis, H., Merbold, L., Montagnani, L., Moors, E., Olesen, J. E., Reichstein, M., Tramontana, G., Gorsel, E., Wohlfahrt, G., and Ráduly, B.: Effect of spatial sampling from European flux towers for estimating carbon and water fluxes with artificial neural networks, J. Geophys. Res.-Biogeo., 120, 1941–1957,, 2015. 

Peltola, O., Mammarella, I., Haapanala, S., Burba, G., and Vesala, T.: Field intercomparison of four methane gas analyzers suitable for eddy covariance flux measurements, Biogeosciences, 10, 3749–3765,, 2013. 

Peltola, O., Vesala, T., Gao, Y., Räty, O., Alekseychik, P., Aurela, M., Chojnicki, B., Desai, A. R., Dolman, A. J., Euskirchen, E. S., Friborg, T., Göckede, M., Helbig, M., Humphreys, E., Jackson, R. B., Jocher, G., Joos, F., Klatt, J., Knox, S. H., Kowalska, N., Kutzbach, L., Lienert, S., Lohila, A., Mammarella, I., Nadeau, D. F., Nilsson, M. B., Oechel, W. C., Peichl, M., Pypker, T., Quinton, W., Rinne, J., Sachs, T., Samson, M., Schmid, H. P., Sonnentag, O., Wille, C., Zona, D., and Aalto, T.: Monthly gridded data product of northern wetland methane emissions based on upscaling eddy covariance observations, Earth Syst. Sci. Data, 11, 1263–1289,, 2019. 

Podani, J.: Extending Gower's general coefficient of similarity to ordinal characters, Taxon, 48, 331–340,, 1999. 

Punj, G. and Stewart, D. W.: Cluster Analysis in Marketing Research: Review and Suggestions for Application, J. Mark. Res., 20, 134–148,, 1983. 

R Core Team: R: A language and environment for statistical computing, Version 4.0.4, R Foundation for Statistical Computing, (last access: 9 May 2022), 2021. 

Reynolds, A. P., Richards, G., de la Iglesia, B., and Rayward-Smith, V. J.: Clustering Rules: A Comparison of Partitioning and Hierarchical Clustering Algorithms, J. Math. Model. Algor., 5, 475–504,, 2006. 

Rinne, J., Riutta, T., Pihlatie, M., Aurela, M., Haapanala, S., Tuovinen, J.-P., Tuittila, E.-S., and Vesala, T.: Annual cycle of methane emission from a boreal fen measured by the eddy covariance technique, Tellus B, 59, 449–457,, 2007. 

Ripley, B. D.: Pattern Recognition and Neural Networks, Cambridge University Press, 403 pp., ISBN 9780521717700, 2007. 

Robertson, G. P., Paul, E. A., and Harwood, R. R.: Greenhouse gases in intensive agriculture: contributions of individual gases to the radiative forcing of the atmosphere, Science, 289, 1922–1925,, 2000. 

Rosentreter, J. A., Borges, A. V., Deemer, B. R., Holgerson, M. A., Liu, S., Song, C., Melack, J., Raymond, P. A., Duarte, C. M., Allen, G. H., Olefeldt, D., Poulter, B., Battin, T. I., and Eyre, B. D.: Half of global methane emissions come from highly variable aquatic ecosystem sources, Nat. Geosci., 14, 225–230,, 2021. 

Saunois, M., Stavert, A. R., Poulter, B., Bousquet, P., Canadell, J. G., Jackson, R. B., Raymond, P. A., Dlugokencky, E. J., Houweling, S., Patra, P. K., Ciais, P., Arora, V. K., Bastviken, D., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Carlson, K. M., Carrol, M., Castaldi, S., Chandra, N., Crevoisier, C., Crill, P. M., Covey, K., Curry, C. L., Etiope, G., Frankenberg, C., Gedney, N., Hegglin, M. I., Höglund-Isaksson, L., Hugelius, G., Ishizawa, M., Ito, A., Janssens-Maenhout, G., Jensen, K. M., Joos, F., Kleinen, T., Krummel, P. B., Langenfelds, R. L., Laruelle, G. G., Liu, L., Machida, T., Maksyutov, S., McDonald, K. C., McNorton, J., Miller, P. A., Melton, J. R., Morino, I., Müller, J., Murguia-Flores, F., Naik, V., Niwa, Y., Noce, S., O'Doherty, S., Parker, R. J., Peng, C., Peng, S., Peters, G. P., Prigent, C., Prinn, R., Ramonet, M., Regnier, P., Riley, W. J., Rosentreter, J. A., Segers, A., Simpson, I. J., Shi, H., Smith, S. J., Steele, L. P., Thornton, B. F., Tian, H., Tohjima, Y., Tubiello, F. N., Tsuruta, A., Viovy, N., Voulgarakis, A., Weber, T. S., van Weele, M., van der Werf, G. R., Weiss, R. F., Worthy, D., Wunch, D., Yin, Y., Yoshida, Y., Zhang, W., Zhang, Z., Zhao, Y., Zheng, B., Zhu, Q., Zhu, Q., and Zhuang, Q.: The global methane budget 2000–2017, Earth Syst. Sci. Data, 12, 1561–1623,, 2020. 

Schaefer, H.: On the Causes and Consequences of Recent Trends in Atmospheric Methane, , Current Climate Change Reports, 5, 259–274,, 2019. 

Schimel, D. and Keller, M.: Big questions, big science: meeting the challenges of global ecology, Oecologia, 177, 925–934,, 2015. 

Schubert, E. and Rousseeuw, P. J.: Faster k-Medoids Clustering: Improving the PAM, CLARA, and CLARANS Algorithms, V1, 171–187,, 2019. 

Schubert, E. and Rousseeuw, P. J.: Fast and eager k-medoids clustering: O(k) runtime improvement of the PAM, CLARA, and CLARANS algorithms, Inf. Syst., 101, 101804,, 2021. 

Serreze, M. C. and Barry, R. G.: Processes and impacts of Arctic amplification: A research synthesis, Glob. Planet. Change, 77, 85–96,, 2011. 

Sherwood, O. A., Schwietzke, S., Arling, V. A., and Etiope, G.: Global inventory of gas geochemistry data from fossil fuel, microbial and burning sources, version 2017, Earth Syst. Sci. Data, 9, 639–656,, 2017. 

Shurpali, N. J., Verma, S. B., Clement, R. J., and Billesbach, D. P.: Seasonal distribution of methane flux in a Minnesota peatland measured by eddy correlation, J. Geophys. Res., 98, 20649,, 1993. 

Sulkava, M., Luyssaert, S., Zaehle, S., and Papale, D.: Assessing and improving the representativeness of monitoring networks: The European flux tower network example, J. Geophys. Res., 116,, 2011. 

Sulla-Menashe, D. and Friedl, M. A.: User guide to collection 6 MODIS land cover (MCD12Q1 and MCD12C1) product, NASA, 2018. 

Thornton, B. F., Wik, M., and Crill, P. M.: Double-counting challenges the accuracy of high-latitude methane inventories, Geophys. Res. Lett., 43, 12569–12577,, 2016. 

Thornton, M. M., Thornton, P. E., Wei, Y., Vose, R. S., and Boyer, A. G.: Daymet: Station-level inputs and model predicted values for North America, Version 3, 2017. 

Torgerson, W. S.: Theory and methods of scaling, Wiley, Oxford, England, ISBN 195907320000, 1958. 

Vaughan, H., Brydges, T., Fenech, A., and Lumb, A.: Monitoring long-term ecological changes through the Ecological Monitoring and Assessment Network: science-based and policy relevant, Environ. Monit. Assess., 67, 3–28,, 2001. 

Venables, W. N. and Ripley, B. D.: Modern applied statistics with S, Springer, New York, NY,, 2002. 

Villarreal, S., Guevara, M., Alcaraz-Segura, D., Brunsell, N. A., Hayes, D., Loescher, H. W., and Vargas, R.: Ecosystem functional diversity and the representativeness of environmental networks across the conterminous United States, Agr. Forest Meteorol., 262, 423–433,, 2018. 

Wik, M., Thornton, B. F., Bastviken, D., Uhlbäck, J., and Crill, P. M.: Biased sampling of methane release from northern lakes: A problem for extrapolation, Geophys. Res. Lett., 43, 1256–1262,, 2016. 

Wilen, B. O. and Bates, M. K.: The US fish and wildlife service's national wetlands inventory project, in: Classification and Inventory of the World's Wetlands, Springer Netherlands, Dordrecht, 153–169,, 1995. 

Xiao, J., Chen, J., Davis, K. J., and Reichstein, M.: Advances in upscaling of eddy covariance measurements of carbon and water fluxes, J. Geophys. Res., 117, G00J01,, 2012. 

Yang, F., Zhu, A.-X., Ichii, K., White, M. A., Hashimoto, H., and Nemani, R. R.: Assessing the representativeness of the AmeriFlux network using MODIS and GOES data, J. Geophys. Res., 113, G04036,, 2008.  

Yumashev, D., Hope, C., Schaefer, K., Riemann-Campe, K., Iglesias-Suarez, F., Jafarov, E., Burke, E. J., Young, P. J., Elshorbany, Y., and Whiteman, G.: Climate policy implications of nonlinear decline of Arctic land permafrost and other cryosphere elements, Nat. Commun., 10, 1900,, 2019. 

Zhang, Z., Zimmermann, N. E., Stenke, A., Li, X., Hodson, E. L., Zhu, G., Huang, C., and Poulter, B.: Emerging role of wetland methane emissions in driving 21st century climate change, P. Natl. Acad. Sci. USA, 114, 9647–9652,, 2017. 

Zhou, X., Zhang, M., Krause, S. M. B., Bu, X., Gu, X., Guo, Z., Jia, Z., Zhou, X., Wang, X., Chen, X., and Wang, Y.: Soil aeration rather than methanotrophic community drives methane uptake under drought in a subtropical forest, Sci. Total Environ., 792, 148292,, 2021. 

Zona, D., Gioli, B., Commane, R., Lindaas, J., Wofsy, S. C., Miller, C. E., Dinardo, S. J., Dengel, S., Sweeney, C., Karion, A., Chang, R. Y.-W., Henderson, J. M., Murphy, P. C., Goodrich, J. P., Moreaux, V., Liljedahl, A., Watts, J. D., Kimball, J. S., Lipson, D. A., and Oechel, W. C.: Cold season emissions dominate the Arctic tundra methane budget, P. Natl. Acad. Sci. USA, 113, 40–45,, 2016. 

Short summary
To understand the CH4 flux potential of natural ecosystems and agricultural lands in the United States 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 land cover.
Final-revised paper