This paper describes the formation of, and initial results for, a new FLUXNET coordination network for ecosystem-scale methane (CH4) measurements at 60 sites globally, organized by the Global Carbon Project in partnership with other initiatives and regional flux tower networks. The objectives of the effort are presented along with an overview of the coverage of eddy covariance (EC) CH4 flux measurements globally, initial results comparing CH4 fluxes across the sites, and future research directions and needs. Annual estimates of net CH4 fluxes across sites ranged from −0.2 ± 0.02 g C m–2 yr–1 for an upland forest site to 114.9 ± 13.4 g C m–2 yr–1 for an estuarine freshwater marsh, with fluxes exceeding 40 g C m–2 yr–1 at multiple sites. Average annual soil and air temperatures were found to be the strongest predictor of annual CH4 flux across wetland sites globally. Water table position was positively correlated with annual CH4 emissions, although only for wetland sites that were not consistently inundated throughout the year. The ratio of annual CH4 fluxes to ecosystem respiration increased significantly with mean site temperature. Uncertainties in annual CH4 estimates due to gap-filling and random errors were on average ±1.6 g C m–2 yr–1 at 95% confidence, with the relative error decreasing exponentially with increasing flux magnitude across sites. Through the analysis and synthesis of a growing EC CH4 flux database, the controls on ecosystem CH4 fluxes can be better understood, used to inform and validate Earth system models, and reconcile differences between land surface model- and atmospheric-based estimates of CH4 emissions.
We describe a new coordination activity and initial results for a global synthesis of eddy covariance CH4 flux measurements.
Atmospheric methane (CH4) is the second-most important anthropogenic greenhouse gas following carbon dioxide (CO2) (Myhre et al. 2013). The concentration of CH4 in the atmosphere today is about 2.5 times higher than in 1750 (Saunois et al. 2016a). The increase in atmospheric CH4 has arisen from human activities in agriculture, energy production, and waste disposal, and from changes in natural CH4 sources and sinks (Saunois et al. 2016a,b, 2017; Turner et al. 2019). Based on top-down atmospheric inversions, global CH4 emissions for the decade of 2003–12 were an estimated ∼420 Tg C yr–1 (range 405–426 Tg C yr–1) (Saunois et al. 2016a). However, some analyses suggest that uncertainties in global CH4 sources and sinks are higher than those for CO2, and uncertainties from natural sources exceed those from anthropogenic emissions (Saunois et al. 2016a). In particular, the largest source of uncertainty in the global CH4 budget is related to emissions from wetlands and inland waters (Saunois et al. 2016a; Melton et al. 2013; Bastviken et al. 2011). Wetland CH4 emissions may contribute as much as 25%–40% of the global total and are a leading source of interannual variability in total atmospheric CH4 concentrations (Bousquet et al. 2006; Chen and Prinn 2006; Saunois et al. 2016a).
Direct, ground-based measurements of in situ CH4 fluxes with high measurement frequency are important for understanding the responses of CH4 fluxes to environmental factors including climate, for providing validation datasets for the land surface models used to infer global CH4 budgets, and for constraining CH4 budgets. Eddy covariance (EC) flux towers measure real-time exchange of gases such as CO2, CH4, water vapor, and energy between the land surface and the atmosphere. The EC technique has emerged as a widespread means of measuring trace gas exchange because it provides direct and near-continuous ecosystem-scale flux measurements without disturbing the soil or vegetation (Baldocchi 2003; Aubinet et al. 2012). There are more than 900 reported active and historical flux tower sites globally and approximately 7,000 site years of data collected (Chu et al. 2017). While most of these sites measure CO2, water vapor, and energy exchange, the development of new and robust CH4 sensors has resulted in a rapidly growing number of CH4 EC measurements (Baldocchi 2014; Morin 2018), primarily in natural and agricultural wetlands (Petrescu et al. 2015).
Since the late 1990s, with a growing number of long-term, near-continuous EC measurements, the EC community has been well coordinated for integrating and synthesizing CO2, water vapor and energy fluxes. This cross-site coordination resulted in the development of regional flux networks for Europe [EuroFlux, CarboEurope, and Integrated Carbon Observing System (ICOS)], Australia (OzFlux), North and South America (AmeriFlux, Large Biosphere Amazon, Fluxnet-Canada/Canadian Carbon Program, and MexFlux), Asia [AsiaFlux, ChinaFlux, KoFlux, and U.S.–China Carbon Consortium (USCCC)], and globally, FLUXNET (Papale et al. 2012; Baldocchi 2014). The resulting FLUXNET database (http://fluxnet.fluxdata.org/) has been used extensively to evaluate satellite measurements, inform Earth system models, generate data-driven CO2 flux products, and provide answers to a broad range of questions about atmospheric fluxes related to ecosystems, land use and climate (Pastorello et al. 2017). FLUXNET has grown steadily over the past 25 years, enhancing our understanding of carbon, water and energy cycles in terrestrial ecosystems (Chu et al. 2017).
Similar community efforts and syntheses for CH4 remain limited in part because EC measurements for CH4 fluxes were rarer until recently. Whereas the earliest EC measurements of CO2 fluxes date back to the late 1970s and early 1980s (Desjardins 1974; Anderson et al. 1984), the first EC CH4 flux measurements only began in the 1990s (Verma et al. 1992; Shurpali and Verma 1998; Fan et al. 1992; Kim et al. 1999), with reliable, easy-to-deploy field sensors only becoming available in the past decade or so. EC CH4 flux measurements became more feasible with advances in sensor development, such as tunable diode laser absorption spectrometers, that allowed researchers to measure previously undetectable trace gas fluxes with higher signal to noise ratios (Rinne et al. 2007; McDermitt et al. 2011). After these new sensors were commercialized, and low-power, low-maintenance open-path sensors were developed that could be operated by solar panels in remote locations, the number of CH4 flux tower measurements increased substantially (Baldocchi 2014; Morin 2018). The rapidly growing number of EC CH4 flux measurements presents new opportunities for FLUXNET-type analyses and syntheses of ecosystem-scale CH4 flux observations.
This manuscript describes initial results from a new coordination activity for flux tower CH4 measurements organized by the Global Carbon Project (GCP) in collaboration with regional flux networks and FLUXNET. The goal of the activity is to develop a global database for EC CH4 observations to answer regional and global questions related to CH4 cycling. Here, we describe the objectives of the FLUXNET-CH4 activity, provide an overview of the current geographic and temporal coverage of CH4 flux measurements globally, present initial analyses exploring time scales of variability, uncertainty, trends, and drivers of CH4 fluxes across 60 sites, and discuss future research opportunities for examining controls on CH4 emissions and reducing uncertainties in the role of wetlands in the global CH4 cycle.
FLUXNET-CH4 SYNTHESIS OBJECTIVES AND TASKS.
This activity is part of a larger GCP effort to establish and better constrain the global methane budget (www.globalcarbonproject.org/methanebudget/index.htm), and is designed to develop a CH4 database component in FLUXNET for a global synthesis of CH4 flux tower data. To this end, we are surveying, assembling, and synthesizing data from the EC community, in coordination with regional networks, including AmeriFlux’s 2019 “Year of Methane” (http://ameriflux.lbl.gov/year-of-methane/year-of-methane/), FLUXNET initiatives, and other complementary activities. In particular, this work is being carried out in parallel with the EU’s Readiness of ICOS for Necessities of Integrated Global Observations (RINGO) project, which is working to standardize protocols for flux calculations, quality control and gap-filling for CH4 fluxes (Nemitz et al. 2018). Methane-specific protocols are needed because of the added complexities and high variability of CH4 flux measurements and dynamics (Nemitz et al. 2018).
Our approach is to include all currently available and future CH4 flux tower observations in a global CH4 database, including freshwater, coastal, natural, and managed ecosystems, as well as upland ecosystems that may be measuring CH4 uptake by soils. The initiative is open to all members of the EC community. Database compilation began in 2017 and is ongoing. Data from sites in the Americas can be submitted to AmeriFlux (http://ameriflux.lbl.gov/data/how-to-uploaddownload-data/); otherwise, data can be submitted to the European Fluxes Database Cluster (www.europe-fluxdata.eu/home/sites-list).
In addition to many applications, an ultimate goal of the FLUXNET-CH4 activity is to generate a publicly available, open-access, data-driven global CH4 emissions product using similar machine-learning-based approaches used for CO2 fluxes (Jung et al. 2009; Tramontana et al. 2016). The product will be based on mechanistic factors associated with CH4 emissions and new spatiotemporal information on wetland area and dynamics for constraining CH4-producing areas. This gridded product will provide an independent bottom-up estimate of global wetland CH4 emissions to compare with estimates of global CH4 emissions from land surface models and atmospheric inversions. Recent work has shown the potential to upscale EC CH4 flux observations across northern wetlands, with predictive performance comparable to previous studies upscaling net CO2 exchange (Peltola et al. 2019); however, our focus is on a globally gridded product.
The near-continuous, high-frequency nature of EC measurements also offers significant promise for improving our understanding of ecosystem-scale CH4 flux dynamics. As such, this synthesis also aims to investigate the dominant controls on net ecosystem-scale CH4 fluxes from hourly to interannual time scales across wetlands globally, and to characterize scale-emergent, nonlinear, and lagged processes of CH4 exchange.
Methane is produced during decomposition under anaerobic or reducing conditions and is transported to the atmosphere via plant-mediated transport, ebullition, and diffusion (Bridgham et al. 2013). During transport, CH4 can pass through unsaturated soil layers and be consumed or oxidized by aerobic bacteria (Wahlen 1993). Process-based biogeochemical models developed and applied at site, regional, and global scales simulate these individual processes with varying degrees of complexity (Bridgham et al. 2013; Melton et al. 2013; Poulter et al. 2017; Castro-Morales et al. 2018; Grant and Roulet 2002). The large range in predicted wetland CH4 emissions rates suggests that there is both substantial parameter and structural uncertainty in large-scale CH4 flux models, even after accounting for uncertainties in wetland areas (Poulter et al. 2017; Saunois et al. 2016a; Melton et al. 2013; Riley et al. 2011). A global EC CH4 database and associated environmental variables can help constrain the parameterization of process-based biogeochemistry models (Saunois et al. 2016a; Bridgham et al. 2013; Oikawa et al. 2017). Furthermore, a key challenge is evaluating globally applicable process-based CH4 models at a spatial scale comparable to model grid cells (Melton et al. 2013; Riley et al. 2011). A globally gridded wetland CH4 emissions product upscaled from EC fluxes can help resolve this issue by providing a scale-appropriate model evaluation dataset. As such, the global CH4 database and gridded product will also be used to parameterize and benchmark the performance of land surface models of global CH4 emissions, providing a unique opportunity for informing and validating biogeochemical models.
Based on a survey of the EC community (announced via the firstname.lastname@example.org and AmeriFlux-Community@lbl.gov listservs), information available in regional networks and FLUXNET, and the scientific literature, we estimate that at least 200 sites worldwide are currently applying the EC method for CH4 flux measurements (Fig. 1). Here we focus on findings from across 60 of the ∼110 sites currently committed to participating in our FLUXNET-CH4 activity [Table A1 in the appendix and Table ES1 in the online supplemental material]. Data from this initial set of sites were selected because they were publicly available or were contributed directly by site principal investigators (PIs). We will continue to engage the EC community more broadly and expand the database in the future.
Data standardization, gap-filling, and partitioning.
We used similar data processing procedures as FLUXNET to standardize and gap-fill measurements, and in the case of net CO2 exchange, partition fluxes across sites (http://fluxnet.fluxdata.org/data/aboutdata/data-processing-101-pipeline-and-procedures/). Standard quality assurance and quality control of the data were first performed by site PIs. In nearly all cases, data collected by the local tower teams were first submitted to the data archives hosted by the regional flux networks, where data are prescreened and formatted based on the regional network data protocols. Data from the regional networks then entered our flux processing procedure.
Within our processing procedure, data were first checked for obvious problems including unit errors, spikes, and out-of-range values based on visualization of the data and statistical metrics. Next, the data were filtered, gap-filled, and partitioned. Friction velocity () filtering, based on relating nighttime CO2 fluxes to , was implemented using the REddyProc package (Wutzler et al. 2018) for R statistical software (R Core Team 2018, version 3.5.0), although in a few cases filtering was performed by the site PIs. Gaps in meteorological variables including air temperature (TA), incoming shortwave (SWIN) and longwave (LWIN) radiation, vapor pressure deficit (VPD), pressure (PA), precipitation (P), and wind speed (WS) were filled with ERA-Interim (ERA-I) reanalysis data (Vuichard and Papale 2015). Gaps in CO2 and latent and sensible heat fluxes were filled using the marginal distribution sampling method (Reichstein et al. 2005) using the REddyProc package (Wutzler et al. 2018). Net CO2 fluxes were partitioned into gross primary production (GPP) and ecosystem respiration (ER) using both the nighttime (Reichstein et al. 2005) and daytime (Lasslop et al. 2010) approaches also implemented in REddyProc (Wutzler et al. 2018).
There are as yet no standards for gap-filling CH4 flux measurements and this is an active and ongoing area of research (Nemitz et al. 2018). Gaps in CH4 fluxes were filled using artificial neural networks (ANNs), as they have shown good performance for gap-filling CH4 flux data (Dengel et al. 2013; Knox et al. 2015; Morin et al. 2014a; Nemitz et al. 2018; Goodrich et al. 2015). Details of the ANN routine are provided in Knox et al. (2016) and are summarized here briefly. The ANN routine was optimized for both generalizability and representativeness. To facilitate representativeness, explanatory data were divided into a maximum of 15 data clusters using the k-means algorithm. To avoid biasing toward conditions with better flux data coverage (e.g., summer and daytime), data used to train, test, and validate the ANN were proportionately sampled from these clusters. Several neural network architectures of increasing complexity were tested, ranging from one hidden layer with the number of nodes equal to the number of explanatory data variables (N) to two hidden layers with 1.5N and 0.75N nodes, respectively. The architecture of each neural network was initialized 10 times with random starting weights, and the initialization resulting in the lowest mean sampling error was selected. The simplest architecture, whereby additional increases in complexity resulted in <5% reduction in mean squared error, was chosen and the prediction saved. This procedure was repeated with 20 resamplings of the data, and missing half hours were filled using the median prediction. A standard set of variables available across all sites were used to gap-fill CH4 fluxes (Dengel et al. 2013), including TA, SWIN, WS, PA, and sine and cosine functions to represent seasonality. These meteorological variables were selected since they are relevant to CH4 exchange and were gap-filled using the ERA-I reanalysis data. Other variables related to CH4 exchange such as water table depth (WTD) or soil temperature (TS) were not included as explanatory variables as they were not available across all sites or had large gaps that could not be filled using the ERA-I reanalysis data. These missing data for variables highlight some of the key challenges in standardizing CH4 gap-filling methods across sites and emphasize the need for standardized protocols of auxiliary measurements across sites (cf. “Future research directions and needs” section) (Nemitz et al. 2018; Dengel et al. 2013). ANN gap-filling was performed using MATLAB (MathWorks 2018, version 9.4.0).
Annual CH4 budgets represent gap-filled, half-hourly fluxes integrated over an entire year or growing season. If fluxes were only measured during the growing season, we assumed that fluxes outside of this period were negligible, although we acknowledge that cold season fluxes can account for as much as ∼13%–50% of the annual CH4 emissions in some locations (Zona et al. 2016; Treat et al. 2018b; Helbig et al. 2017a; Kittler et al. 2017).
ANNs were also used to estimate annual gap-filled and random uncertainty in CH4 flux measurements (Richardson et al. 2008; Moffat et al. 2007; Anderson et al. 2016; Knox et al. 2018). Here, we focus on assessing the random error, but a full assessment of total flux measurement error also requires quantifying systematic error or bias (Baldocchi 2003). Systematic errors, due to incomplete spectral response, lack of nocturnal mixing, submesoscale circulations, and other factors are discussed elsewhere (Baldocchi 2003; Peltola et al. 2015) and are the focus of other ongoing initiatives.
Random errors in EC fluxes follow a double exponential (Laplace) distribution with a standard deviation varying with flux magnitude (Richardson et al. 2012, 2006). Model residuals of gap-filling algorithms such as ANNs provide a reliable, and conservative “upper limit,” estimate of the random flux uncertainty (Moffat et al. 2007; Richardson et al. 2008). For half-hourly CH4 flux measurements, random error was estimated using the residuals of the median ANN predictions. At each site, the probability density function (PDF) of the random flux measurement error more closely followed a double-exponential (Laplace) rather than normal (Gaussian) distribution, with the root-mean-square error (RMSE) for the Laplace distribution fitted to the PDF of random errors consistently lower than the normal distributed error. From half-hourly flux measurements, random error can also be estimated using the daily differencing approach (Richardson et al. 2012). Random error estimates [σ(δ)], as expressed as the standard deviation of the double-exponential distribution with scaling parameter β, where (Richardson et al. 2006), were found to be nearly identical using the two approaches [σ(δ)model_residual = 1.0 × σ(δ)daily_differencing + 1.21; r2 = 0.97, p < 0.001], supporting the use of the model residual approach for estimating random error. As discussed below, σ(δ) scaled linearly with the magnitude of CH4 fluxes at nearly all sites. To quantify random uncertainty of cumulative fluxes, we used a Monte Carlo simulation that randomly draws 1,000 random errors for every original measurement using σ(δ) binned by flux magnitude, and then computed the variance of the cumulative sums (Anderson et al. 2016). For gap-filled values, the combined gap-filling and random uncertainty was calculated from the variance of the cumulative sums of the 20 ANN predictions (Anderson et al. 2016; Oikawa et al. 2017; Knox et al. 2015). The annual cumulative uncertainty at 95% confidence was estimated by adding the cumulative gap-filling and random measurement uncertainties in quadrature (Richardson and Hollinger 2007; Anderson et al. 2016). Note that when reporting mean or median annual CH4 fluxes across sites, error bars represent the standard error.
Wavelet-based time-scale decomposition.
Methane fluxes are highly dynamic and vary across a range of time scales (Sturtevant et al. 2016; Koebsch et al. 2015). For example, in wetlands with permanent inundation, the seasonal variation of CH4 exchange is predominantly controlled by temperature and plant phenology (Chu et al. 2014; Sturtevant et al. 2016). Ecosystem CH4 exchange also varies considerably at both longer (e.g., interannual; Knox et al. 2016; Rinne et al. 2018) and shorter (e.g., weeks, days, or hours; Koebsch et al. 2015; Hatala et al. 2012; Schaller et al. 2018) time scales. Wavelet decomposition is a particularly useful tool for investigating scale in geophysical and ecological analysis (Cazelles et al. 2008; Torrence and Compo 1998), because it can characterize both the time scale and location of patterns and perturbations in the data. Partitioning variability across temporal scales can help to isolate and characterize important processes (Schaller et al. 2018).
The maximal overlap discrete wavelet transform (MODWT) was used to decompose the time scales of variability in gap-filled CH4 flux measurements, as described in Sturtevant et al. (2016). The MODWT allows the time series to be decomposed into the detail added from progressively coarser to finer scales and either summed or treated individually to investigate patterns across scales. We reconstructed the detail in the fluxes for dyadic scales 1 (21 measurements = 1 h) to 14 (214 measurements = 341 days). Since patterns generated by ecological processes tend to occur over a scale range rather than at one individual scale, the detail over adjacent scales were summed to analyze four general time scales of variation (Sturtevant et al. 2016). These time scales included the “hourly” scale (1–2 h) representing perturbations such the passage of clouds overhead and turbulent scales up to the spectral gap, the “diel” scale (4 h to 1.3 days) encompassing the diel cycles in sunlight and temperature, the “multiday” scale (2.7 to 21.3 days) reflecting synoptic weather variability or fluctuations in water levels, and the “seasonal” scale (42.7 to 341 days) representing the annual solar cycle and phenology. Data were wavelet decomposed into the hourly, diel, and multiday scales using the Wavelet Methods for Time Series Analysis (WMTSA) Wavelet Toolkit in MATLAB.
We tested for significant relationships between log-transformed annual CH4 emissions and a number of covariates using linear mixed-effects models as described in Treat et al. (2018b). The predictor variables of CH4 flux we evaluated included: biome or ecosystem type (categorical variables), and continuous biophysical variables including mean seasonal WTD, mean annual soil and air temperature (TMST and TMAT, respectively), net ecosystem exchange (NEE), GPP, and ER. When considering continuous variables, we focused on freshwater wetlands for comparison with previous CH4 synthesis activities. Soil temperature was measured between 2 and 25 cm below the surface in different studies. The results below are presented for GPP and ER covariates that are partitioned using the nighttime flux partitioning algorithm (Wutzler et al. 2018; Reichstein et al. 2005), although similar findings were obtained using daytime partitioned estimates. Additionally, individual sites or site years were excluded when gaps in measurements exceeded two consecutive months, which explains the differences in the number of sites and site years in the “Environmental controls on annual CH4 emissions across freshwater wetland sites” section below.
Mixed-effects modeling was used because of the potential bias of having measurements over several years, with site included as a random effect in the analysis (Treat et al. 2018b). The significance of individual predictor variables was evaluated using a χ2 test against a null model using only site as a random variable (Bates et al. 2015), with both models fit without reduced maximum likelihood. For multiple linear regression models, we used the model selection process outlined in Zuur et al. (2009). To incorporate annual cumulative uncertainty when assessing the significance of trends and differences in annual CH4 fluxes across biomes and ecosystem types, we used a Monte Carlo simulation that randomly draws 1,000 annual cumulative uncertainties for each estimate of annual CH4 flux. For each random draw the significance of the categorical variable was tested using a χ2 test against the null model with only site as a random variable. We report the marginal r2 (r2m), which describes the proportion of variance explained by the fixed factors alone (Nakagawa and Schielzeth 2013). The mixed-effects modeling was implemented using the lmer command from the lme4 package (Bates et al. 2015) for R statistical software.
RESULTS AND DISCUSSION.
Geographic and temporal coverage of eddy covariance CH4 flux measurements.
We identified 200 sites worldwide that are applying the EC method for CH4 (Fig. 1; Table ES1); wetlands (including natural, managed, and restored wetlands) comprise the majority of sites (59%), with rice agriculture (10%) as the second-most represented vegetation type. The predominance of wetland and rice paddy sites in the database is unsurprising because many studies are designed to target ecosystems expected to have relatively large CH4 emissions. However, there are also sites in ecosystems that are typically smaller sources or even sinks of CH4 such as upland forests (13%) and grasslands (8%). Additionally, six sites (∼3%) are urban, with another five sites measuring CH4 fluxes from open water bodies. Although identified sites span all continents except Antarctica, the majority are concentrated in North America and Europe, with a growing number of sites in Asia (Fig. 1; Table ES1).
Measurements of CH4 fluxes cover a broad range of climates and a large fraction of wetland habitats (Fig. 2), with the tropics and tropical wetlands notably underrepresented. As discussed below (see “Future research directions and needs” section), one important goal of FLUXNET and the regional networks is to increase site representativeness and extend measurements in undersampled regions. Increasing the number of tropical sites is particularly important for CH4 because more than half of global CH4 emissions are thought to come from this region (Saunois et al. 2016a; Dean et al. 2018). Furthermore, compared to northern wetlands, their biogeochemistry remains relatively poorly understood (Mitsch et al. 2009; Pangala et al. 2017). We expect the number of CH4 flux sites and their geographic and temporal coverage to continue to increase, as has occurred through time for CO2, water vapor, and energy flux measurements in FLUXNET (Pastorello et al. 2017; Chu et al. 2017).
Long-term CH4 flux time series are key to understanding the causes of year-to-year variability and trends in fluxes (Chu et al. 2017; Euskirchen et al. 2017; Pugh et al. 2018). The longest continuous record of CH4 flux measurements, from a fen in Finland (Rinne et al. 2018), is now ∼14 years and ongoing (Table ES1). Three other sites have measurements exceeding 10 years; however, the median length is 5 years, with most sites established from 2013 onward (Table ES1). Longer time series are also important for both exploring the short- and long-term effects of extreme events on fluxes and tracking the response of disturbed or restored ecosystems over time (Pastorello et al. 2017). Furthermore, they can help address new and emerging science questions, such as quantifying CH4 feedbacks to climate with rising temperatures and associated changes in ecosystem composition, structure and function (Helbig et al. 2017a,b; Dean et al. 2018), and the role of wetland emissions in atmospheric CH4 variability (McNorton et al. 2016; Poulter et al. 2017).
CH4 fluxes and trends across biomes and ecosystem types.
Half-hourly and annual net CH4 fluxes for the 60 sites currently included in the database exhibited strong variability across sites (Figs. 3 and 4). Across the dataset, the mean half-hourly CH4 flux was greater than the median flux, indicating a positively skewed distribution with infrequent, large emissions (Fig. 3a), similar to findings from chamber-based syntheses (Olefeldt et al. 2013; Turetsky et al. 2014). Mean and median CH4 fluxes were smaller at higher latitudes and larger at lower latitudes (Fig. 3b), comparable again to trends in CH4 fluxes observed in predominantly chamber-based syntheses (Bartlett and Harriss 1993; Turetsky et al. 2014; Treat et al. 2018b).
The continuous nature of EC flux measurements is well suited for quantifying annual ecosystem-scale CH4 budgets, along with accumulated uncertainty (cf. “Gap-filling performance and uncertainty quantification” section). Annual estimates of net CH4 flux for each of the 60 sites in the flux tower database ranged from −0.2 ± 0.02 g C m–2 yr–1 for an upland forest site to 114.9 ± 13.4 g C m–2 yr–1 for an estuarine freshwater marsh (Rey-Sanchez et al. 2018), with fluxes exceeding 40 g C m–2 yr–1 at multiple sites (Fig. 4b). These emissions are of a considerably broader range and have much higher annual values than in an earlier synthesis by Baldocchi (2014), which included published values from 13 sites (Fig. 4a); median annual CH4 fluxes (±SE) in that study were 6.4 ± 1.9 g C m–2 yr–1, compared with 10.0 ± 1.6 g C m–2 yr–1 for our expanded database. Annual CH4 sums in our database were positively skewed, with skewness increasing with additional observations due largely to the inclusion of high CH4-emitting freshwater marsh sites (Fig. 4).
As suggested from Fig. 3b, annual wetland CH4 emissions differed significantly among biomes, even when considering accumulated uncertainty [average Monte Carlo χ2 = 13.4 (12.1–14.7, 95% confidence interval), degrees of freedom (df) = 3, p < 0.05] (Table 1). Median CH4 emissions were significantly lower for tundra wetlands (2.9 ± 1.3 g C m–2 yr–1) than temperate wetlands (27.4 ± 3.4 g C m–2 yr–1). Higher CH4 emissions were observed from subtropical/tropical wetlands (43.2 ± 11.2 g C m–2 yr–1), based on only three site years of data; however, emphasizing the need for additional flux tower measurements in the tropics.
Whereas annual boreal/taiga wetland CH4 emissions were comparable to values reported in a recent synthesis of predominantly chamber-based CH4 flux measurements (Treat et al. 2018b), our tower-based measurements are ∼50% lower and over 6 times higher for tundra and temperate wetlands, respectively (Table 1). The inconsistencies highlighted in Table 1 not only reflect the differences in the number and location of sites between datasets, but also the discrepancies resulting from different measurement techniques. Several studies have noted considerable differences in CH4 emissions measured using EC and chamber techniques, with estimates from chambers often higher than those from the EC measurements (Schrier-Uijl et al. 2010; Hendriks et al. 2010; Meijide et al. 2011; Krauss et al. 2016). This distinction highlights the need for additional studies investigating the systematic differences caused by the different spatial and temporal sampling footprints of these methods (Krauss et al. 2016; Morin et al. 2017; Windham-Myers et al. 2018; Xu et al. 2017). Characterizing discrepancies between measurement techniques may also help constrain bottom-up estimates of CH4 emissions and reduce the disagreement of ∼15 Tg C yr–1 between bottom-up (139 Tg CH4 yr–1) and top-down (125 Tg CH4 yr–1) estimates of CH4 emissions from natural wetlands (Saunois et al. 2016a).
Annual CH4 emissions also differed significantly across ecosystems [average Monte Carlo χ2 = 45.5 (39.3–50.1), df = 9, p < 0.001; Fig. 5], with median fluxes highest for freshwater marshes (43.2 ± 4.2 g C m–2 yr–1) and lowest for upland ecosystems (1.3 ± 0.7 g C m–2 yr–1). Treat et al. (2018b) also observed the highest annual emissions in marshes and reported a similar median value for temperate marshes (49.6 g C m–2 yr–1). Wet tundra and bogs had significantly lower annual emissions than marshes (Fig. 5), which in part reflects their presence in colder boreal and tundra systems, as well as differences in vegetation type, nutrient status, and hydrological regime (Treat et al. 2018b). Low median CH4 emission was observed from salt marshes in our dataset (0.8 ± 2.9 g C m–2 yr–1), because high sulfate concentrations inhibit methanogenesis (Poffenbarger et al. 2011; Holm et al. 2016). Even drained wetlands converted to agricultural land can be large sources of CH4 associated with seasonal flooding (Fig. 5). Median annual CH4 flux from rice was 12.6 ± 1.6 g C m–2 yr–1, which is slightly lower than the IPCC default value of 15 g C m–2 yr–1 (Sass 2003).
Environmental controls on annual CH4 emissions across freshwater wetland sites.
Using an integrated CH4 flux database, we can begin to investigate the factors associated with varying CH4 emissions across sites. We explored the effects of WTD, TMST or TMAT, NEE, GPP, and ER on annual CH4 flux. At global scales, TMAT and TMST were the most important predictors of annual CH4 flux across wetland sites (p < 0.001 for each), with the fixed factor of TMAT or TMST explaining ∼65% of the variation in log transformed annual CH4 emission (Figs. 6a,b). Previous synthesis studies also observed a significant, but weaker, relationship between soil temperature and average CH4 emissions across sites, explaining <15% of the variation in CH4 flux in those studies (Olefeldt et al. 2013; Yvon-Durocher et al. 2014). However, our findings are consistent with numerous site-level studies that report a strong correlation between wetland CH4 emissions and temperature, with nearly 95% of all EC studies reporting a significant relationship between temperature and CH4 flux (Morin 2018). Across sites, Peltola et al. (2019) found that the most important predictor in a random forest model used to upscale EC CH4 emissions across northern latitudes was temperature, again highlighting the importance of temperature in regulating CH4 emissions within and across sites.
Water table depth has also commonly been identified as a key control on CH4 emissions (Turetsky et al. 2014; Bubier et al. 2005), because higher water levels often inhibit oxygen availability and lower the soil reduction potential, making methanogenesis more thermodynamically favorable. Although predominantly chamber-based wetland CH4 syntheses have found a positive relationship between WTD and average or annual CH4 emissions across sites (Olefeldt et al. 2013; Turetsky et al. 2014; Treat et al. 2018b), we observed no significant relationship between mean WTD and annual CH4 flux across all sites (χ2 = 0.2, df = 1, p = 0.66, Nsites = 20, Nsite_yr = 46), even when considering WTD2 or WTD3 (Olefeldt et al. 2013). However, if we consider only sites where WTD was below the soil surface for part or all of the year (Fig. 6c, solid circles), we did observe a significant relationship with WTD (p < 0.05). Conversely, CH4 emissions for permanently inundated sites showed no significant relationship with WTD (Fig. 6c, open circles) (χ2 = 0.5, df = 1, p = 0.50, Nsites = 13, Nsite_yr = 19). This result supports the finding that wetlands that are permanently inundated or exhibit little variation in WTD tend to show weak to no correlation between WTD and CH4 emissions (Chu et al. 2014; Jackowicz-Korczyński et al. 2010; Rinne et al. 2007; Christensen et al. 2003); in contrast, wetlands with lower and more variable water levels often have a significant relationship between WTD and CH4 emissions (Bubier et al. 2005; Treat et al. 2007). However, only half of the sites currently included in the database report water table position, and given the importance of WTD in regulating CH4 exchange, it is critical to ensure that WTD is measured across all sites.
Gross primary production and ER were both significant positive predictors of annual CH4 flux (χ2 = 21.3, df = 1, p < 0.001, r2m = 0.29 and χ2 = 17.1, df = 1, p < 0.001, r2m = 0.25, respectively, Nsites = 26, Nsite_yr = 64), although there was no significant relationship between NEE and annual CH4 flux (χ2 = 0.9, df = 1, p = 0.33, Nsites = 2, Nsite_yr = 64). However, when considering GPP or ER in a multiple linear regression model with TMST, including interaction terms (Chu et al. 2014), neither GPP nor ER were significant, suggesting that the observed relationship with GPP or ER was due to covariation with soil temperature and, possibly, other environmental drivers.
The strong temperature dependence of ecosystem-scale CH4 emissions we observed across wetland sites is in line with the high temperature sensitivity of CH4 emissions found across microbial to ecosystem scales (Yvon-Durocher et al. 2014). CH4 emissions also have a higher temperature dependence than ER, such that the ratio of CH4 to CO2 emissions was found to increase markedly with seasonal increases in temperature (Yvon-Durocher et al. 2014). Similarly, we observed a significant increase in the ratio of annual CH4 to ER along geographic temperature gradients, ranging from 0.4% to 7.9%, with a median value of 2.8% across the dataset (Fig. 6d). This relationship suggests that warming may result in a greater relative contribution of CH4 to total carbon emissions from wetland ecosystems. With a growing FLUXNET CH4 database, it will be possible to further explore the dominant controls on CH4 fluxes within and across ecosystem types, as well as further investigate the temperature dependence of ecosystem-scale CH4 exchange (Schipper et al. 2014; Arcus et al. 2016; Yvon-Durocher et al. 2014).
Time scales of variability.
Methane fluxes exhibited strong variability over a range of time scales, with the variation across time scales differing between wetland types (Fig. 7). As observed previously (Sturtevant et al. 2016), the seasonal time scale dominated CH4 flux variability across wetland types, but was most pronounced in rice paddies, which have a distinct growing season, and least pronounced in bogs (Fig. 7). Across ecosystem types, variation was lowest at the multiday scale, although multiday CH4 flux variation was slightly greater in rice paddies and wet tundra, potentially indicating greater water table fluctuations (Sturtevant et al. 2016), particularly at rice paddy sites, which are subject to seasonal drainage (Knox et al. 2016; Runkle et al. 2019). Whereas some studies report a strong diel pattern in CH4 emissions from wetlands and rice paddies (Knox et al. 2016; Chu et al. 2014; Morin et al. 2014b; Kim et al. 1999), others have found little or no diel variation (Rinne et al. 2018; Jackowicz-Korczyński et al. 2010; Yagi and Minami 1990; Nadeau et al. 2013). Across wetland types, diel variation was greatest in freshwater marshes (Fig. 7), which is consistent with the observations that the vegetation at sites with a strong diel cycle of CH4 emissions is typically dominated by species with convective gas flow such as Phragmites australis or Typha spp. (Brix et al. 1992; Chanton et al. 1993). Bogs, fens, and wet tundra showed the greatest variation at the hourly scale (Fig. 7). This is likely in part due to typically lower fluxes at these sites as hourly perturbations of turbulent time series are largely dominated by noise (Hollinger and Richardson 2005), as well as the fact that near-surface turbulence and short-term pressure fluctuations can strongly influence CH4 exchange in these peat dominated ecosystems (Nadeau et al. 2013; Sachs et al. 2008).
Gap-filling performance and uncertainty quantification.
The performance of the neural networks varied strongly across sites (Fig. 8). Model r2, calculated from the median ANN prediction and observed fluxes at each site, ranged from ∼0 to 0.92 across sites, with a median value of 0.41. Across sites, ANN performance was strongly linked to the percentage of total variance at diel and seasonal scales (r2 = 0.69, p < 0.001), indicating that across the wide range of observed flux magnitudes, sites with a more distinct seasonal and diel pattern tended to be more predictable (Fig. 8). There was also a significant negative relationship between model r2 and the percentage of total variance at the hourly scale across sites (r2 = 0.72, p < 0.001), because, as noted previously, hourly perturbations are largely dominated by noise (Hollinger and Richardson 2005).
Knowledge of the random errors in half-hourly flux measurements is not only important for evaluating the uncertainty in cumulative fluxes (e.g., daily, monthly, or annual) and comparing fluxes across tower sites, but it also needed to incorporate information about random flux errors in model-data synthesis activities (Richardson et al. 2006). As noted above, random flux error more closely followed a Laplace rather than Gaussian distribution. Within sites, σ(δ) was not constant, but rather nearly always scaled with the magnitude of CH4 fluxes (Fig. 9a), as predicted from theory (Richardson et al. 2006). As observed for other fluxes (Richardson et al. 2006), both the slope and intercept of this relationship varied among sites, and depending on the sign of the flux (Fig. 9a). Across sites, random flux error therefore scaled linearly with the magnitude of mean CH4 flux (r2 = 0.86, p < 0.001), even when excluding the two highest CH4-emitting sites (r2 = 0.46, p < 0.001) (Fig. 9b). Whereas closed-path CH4 analyzers have been found to have lower random errors and instrument noise compared with open-path sensors (Peltola et al. 2014), there was no clear evidence of a systematic effect of the influence of closed- versus open-path sensors on random errors across sites (Fig. 9).
The total annual cumulative uncertainty in CH4 fluxes, including both random and gap-filling errors, ranged from ±0.01 to ±13.4 g C m–2 yr–1, with a median value of ±1.0 g C m–2 yr–1 at 95% confidence (Fig. 10a). Relative error decreased exponentially with flux magnitude, ranging from 1.5% to 60% in most cases (Fig. 10b), although a few sites where annual CH4 sums were near zero had relative errors exceeding 200% (data not shown). The highest relative errors therefore tended to be associated with low CH4-emitting sites, such as upland sites and bogs, and the lowest relative errors were generally associated with high CH4-emitting sites such as freshwater marshes (Fig. 10b).
FUTURE RESEARCH DIRECTIONS AND NEEDS.
Better quantification of CH4 sources and sinks will improve estimates of regional and global CH4 budgets and reduce uncertainties in the CH4 cycle. In this general context, high-frequency observations of ecosystem-scale CH4 emissions should help constrain bottom-up CH4 budgets, improve our understanding of the environmental factors controlling CH4 fluxes, and inform and validate land surface models used to estimate global CH4 fluxes. Unlike well-established efforts synthesizing CO2, water vapor, and energy observations, no such global data synthesis or initiative previously existed for CH4. The database presented here addresses this gap with the EC community by organizing the collection and aggregation of a global EC CH4 database through FLUXNET.
EC flux data quality assessment.
Much of what has been learned within FLUXNET for CO2, water vapor, and energy measurements is informing, and should continue to inform, new efforts for CH4. Reliable EC measurements of CO2 and water vapor fluxes have been conducted at hundreds of sites across broad regional networks (Papale et al. 2012), and substantial efforts have focused on developing best practices and harmonizing approaches across sites to ensure consistent, high-quality flux measurements (Aubinet et al. 1999; Reichstein et al. 2005; Moffat et al. 2007). CH4 fluxes are often characterized by small fluxes with episodic spikes, and additional research is needed to ensure reliable measurements (Peltola et al. 2014, 2013), and refine and standardize methods and routines for data processing and quality checking (Nemitz et al. 2018; Schaller et al. 2018). Recent efforts provided guidance on instrument selection, setup and maintenance, and data processing for EC CH4 flux measurements (Nemitz et al. 2018). However, with respect to instrument setup and data processing, more research is needed in best practices for storage flux quantification, despiking, and u* filtering (Nemitz et al. 2018).
Whereas neural networks have shown strong performance for gap-filling CH4 fluxes (Dengel et al. 2013; Knox et al. 2016), our results reveal some of the challenges of gap-filling CH4 fluxes at sites with low fluxes and/or a lack of seasonal and diel variation (Fig. 8). More research is therefore needed for best practices for gap-filling to estimate annual CH4 budgets (Nemitz et al. 2018). For example, there has yet to be a comprehensive analysis comparing a wide range of gap-filling approaches for CH4 fluxes similar to the study by Moffat et al. (2007) for CO2 exchange. While ANNs are one gap-filling method (Dengel et al. 2013; Shoemaker et al. 2014; Morin et al. 2014a), numerous other gap-filling approaches exist, including nonlinear regression techniques, mean diurnal variation, lookup tables, marginal distribution sampling, and the multiple imputation method (Moffat et al. 2007; Vitale et al. 2019). Future efforts should focus on systematically investigating these approaches across a range of sites to provide best practices for gap-filling CH4 exchange.
Along with research that addresses the challenges of measuring and processing EC CH4 fluxes, key ancillary variables to help gap-fill, predict, and scale CH4 fluxes should also be measured more comprehensively across sites. For instance, although WTD is known to strongly influence CH4 emissions (Turetsky et al. 2014; Treat et al. 2018b), as noted above, only half of the sites currently included in the database report water table position. Generally, EC CH4 measurements are implemented at sites also collecting CO2 fluxes and common meteorological measurements used in the flux community; however, guidelines are only beginning to emerge for which additional supporting variables should be collected at sites measuring CH4 fluxes (Nemitz et al. 2018).
Measurements of variables beyond those relevant for CO2 are needed to better understand and predict the complex and interacting processes of CH4 production, consumption, and transport, the latter of which includes diffusion, ebullition, and plant-mediated transport. Guidance on the description of some basic variables affecting these processes is available through new protocols in the flux community detailing soil meteorological measurements, ancillary vegetation measurements, and site description, management and disturbance (Saunders et al. 2018; Op De Beeck et al. 2018; Gielen et al. 2018). These protocols provide guidance on variables such as soil temperature and soil moisture profiles, water table depth and snow depth, soil pH and soil type, bulk density, and livestock density. However, although WTD is an easily measured proxy for anaerobic conditions, direct and continuous measurement of redox potential and oxygen content in particular would be valuable additional measurements (Nemitz et al. 2018). Similarly, measuring variables such as conductivity, below-ground CH4 concentrations, dissolved organic carbon concentrations, and the presence of alternative electron acceptors such as nitrate, iron, sulfate, and humic substances in the water and soil column would provide useful information for the interpretation of CH4 emissions. Stable isotope analyses of CH4 are also valuable as they provide important information on mechanisms of CH4 production, transport, and oxidation (Chanton et al. 1997; Marushchak et al. 2016). Detailed information on soil microbial communities driving CH4 production and consumption could also be helpful (Kwon et al. 2017). Vegetation biomass, species composition, and phenology are also important variables to consider, because plants are a primary source of carbon substrates for methanogenic metabolism, and they mediate CH4 transport through aerenchymous tissue (Kwon et al. 2017; Joabsson et al. 1999; Carmichael et al. 2014). New guidance is now available for such measurements at flux tower locations (Gielen et al. 2018; Hufkens et al. 2018). Continuing to develop a consensus on best practices for ancillary measurements is important for interpreting, gap-filling, and upscaling CH4 flux measurements.
Characterizing spatial variability.
CH4 fluxes exhibit finescale spatial variability that can span orders of magnitude within a landscape (Peltola et al. 2015; Marushchak et al. 2016; Desai et al. 2015; Treat et al. 2018a; Iwata et al. 2018), attributable to heterogeneous soil properties and moisture conditions, vegetation composition, and land use (Davidson et al. 2016; Parmentier et al. 2011; Chamberlain et al. 2018). Furthermore, there is evidence that traditionally unmeasured surfaces (i.e., tree stems) are important sources of CH4 to the atmosphere and could explain spatial heterogeneity within ecosystems (Barba et al. 2019). Accurately representing spatial heterogeneity and the relative fraction of uplands and wetlands is imperative for interpreting and predicting CH4 emissions within many ecosystems, and for upscaling flux measurements regionally and globally as wetlands are hot spots for carbon cycling (Treat et al. 2018a; Tuovinen et al. 2019; Rößger et al. 2019). Flux footprint analysis characterizing the fractional coverage of the dominant surface types, particularly the fraction of open water and aerenchymatous plants, is important for interpreting EC CH4 flux measurements and quantifying annual CH4 budgets at spatially heterogeneous sites (Franz et al. 2016; Helbig et al. 2017a; Jammet et al. 2017) (Fig. 11). This integration can be achieved by combining CH4 measurements, flux footprint analysis, and near-surface (e.g., phenocams) and/or high-resolution drone or satellite remote sensing data, and should be common practice for all sites measuring CH4 fluxes.
Spatial variability in ecosystem-scale CH4 flux can further be examined by combining chamber and EC measurements, including manual and autochambers, multitower approaches, and airborne flux measurements (Peltola et al. 2015; Zona et al. 2016; Helbig et al. 2017a; Wolfe et al. 2018; Kohnert et al. 2018; Lai et al. 2014; McNicol et al. 2017). Integrating additional observations such as information on microbial communities, isotopic measurements, and laboratory incubation observations along with chamber and EC CH4 flux measurements can further help explain CH4 dynamics across scales (Angle et al. 2017; Chamberlain et al. 2018; Yang et al. 2017). However, as discussed above, additional research is needed to reconcile differences in fluxes measured across scales (Gioli et al. 2004; Holm et al. 2016; Meijide et al. 2011). Explicitly considering source area composition and spatial heterogeneity will provide enhanced processed-based understanding of CH4 fluxes and improve upscaled regional and global estimates of CH4 emissions, which can help reconcile the discrepancy between bottom-up and top-down budgets (Saunois et al. 2016a; Morin et al. 2017; Davidson et al. 2016).
More sites in key regions.
We expect the number of flux towers measuring CH4 fluxes will continue to grow (Chu et al. 2017; Pastorello et al. 2017; Morin 2018), but our compilation of EC CH4 flux sites highlights key underrepresented regions where future flux towers are needed or where more efforts are needed for existing but nonreporting towers to contribute to FLUXNET (Fig. 1). As noted previously, notable gaps include both tropical and subtropical regions, as well as eastern Canada, and the boreal forests of Russia. Figure 1 also provides guidance on where new towers could be strategically located to help reconcile differences between top-down and bottom-up estimates of wetland CH4 emissions. In particular, substantial disagreements between top-down and bottom-up estimates are found over the Congo basin, the Inner Niger delta, the Orinoco River delta, the Maranon–Ucayali palm swamps, the Pantanal, the Ganges–Brahmaputra delta, Sumatra, the western Siberian lowlands, and the Hudson Bay lowlands (Fig. 1). However, the placement of new towers is a strong function of the scientific question being asked and research funding priorities, and therefore the optimal tower network could be different for different applications (Mahecha et al. 2017; Papale et al. 2015; Villarreal et al. 2018).
Better understanding and representing processes.
One of the biggest challenges for understanding ecosystem functioning is resolving overlapping, asynchronous (i.e., lagged) and nonlinear processes (Sturtevant et al. 2016). This challenge is particularly relevant for interpreting continuous, ecosystem-scale measurements of CH4 exchange where scale-specific, nonlinear, and lagged processes may dominate (Franz et al. 2016; Sturtevant et al. 2016; Knox et al. 2018). For instance, CH4 emission responses to water table fluctuation can be nonlinear and lagged on the order of days to months (Goodrich et al. 2015; Sturtevant et al. 2016). CH4 flux has also been observed to lag GPP by hours to days (Rinne et al. 2018; Hatala et al. 2012). Adequately representing these dynamics in process models is important, and further research is needed to better characterize the complex and nonlinear processes influencing ecosystem-scale CH4 exchange across time scales.
The complex nature of CH4 flux dynamics requires moving beyond traditional linear correlation and regression, and using methods such as wavelets, information theory, and Granger causality that are more tailored to address scale, nonlinearity, and lags directly (Stoy et al. 2005; Vargas et al. 2011; Schäfer et al. 2014; Knox et al. 2016; Detto et al. 2012). Through a USGS Powell Center working group activity, we will continue to investigate controls on CH4 emissions within and across wetland types. To further explore interactions between ecosystem-scale CH4 exchange and drivers across time scales, wavelet analysis will be combined with information theory to explore biosphere–atmosphere interactions regardless of form or asynchrony (Sturtevant et al. 2016; Knox et al. 2018; Chamberlain et al. 2018). By coupling wavelet decomposition with information theory, future research will investigate key controls on CH4 fluxes across time scales, as well as the importance of nonlinearities and lags in predicting CH4 flux dynamics. Future research will also use the global CH4 database to parameterize and benchmark the performance of land surface models of global CH4 emissions, providing a unique opportunity for informing and validating biogeochemical models.
Coordinating, organizing and improving the integration of CH4 fluxes in regional networks and ultimately FLUXNET will bring us one step closer to achieving the goal of providing flux information “everywhere and all of the time” (Baldocchi 2008). In the long term, we hope to integrate the global eddy covariance CH4 database with other methods for measuring CH4 fluxes, such as chamber, aircraft, and satellite measurements. By integrating CH4 flux measurements, remote sensing, and modeling, we aim to better characterize CH4 emissions from terrestrial ecosystems and ultimately reduce uncertainties in the global CH4 cycle.
This study was supported by the Gordon and Betty Moore Foundation through Grant GBMF5439 “Advancing Understanding of the Global Methane Cycle” to Stanford University supporting the Methane Budget activity for the Global Carbon Project (globalcarbonproject.org). This work was also conducted as a part of the Wetland FLUXNET Synthesis for Methane Working Group supported by the John Wesley Powell Center for Analysis and Synthesis of the U.S. Geological Survey. We acknowledge the following AmeriFlux sites for their data records: CA-SCC, CA-SCB, U.S.-Beo, U.S.-Bes, U.S.-Atq, U.S.-Ivo, U.S.-ICs, U.S.-NGB, U.S.-NGC, U.S.-Los, U.S.-NC4, U.S.-HRA, U.S.-HRC, U.S.-Myb, U.S.-Sne, U.S.-Tw1, U.S.-Tw4, U.S.-Twt, U.S.-Bi2, U.S.-Bi1, U.S.-Snd, U.S.-OWC, U.S.-ORv, U.S.-WPT, U.S.-CRT, U.S.-MRM, U.S.-Srr, U.S.-Uaf, U.S.-Ho1, U.S.-StJ, U.S.-LA2, U.S.-LA1. In addition, funding for AmeriFlux data resources were provided by the U.S. Department of Energy’s Office of Science. BRKR was supported by NSF Award 1752083. WCO and DZ acknowledge support from NSF 1204263 and 1702797; NASA ABoVE NNX15AT74A, NNX16AF94A, and NNX17AC61A; NOAA EPP NA16SEC4810008; EU Horizon 2020 No. 629 727890; and NERC UK NE/P002552/1. TFK and HC acknowledge support from the AmeriFlux Management Project. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. government. DP, IM, and TV thank the European Union for supporting the RINGO project funded by the Horizon 2020 Research and Innovation Programme under Grant Agreement 730944. The authors also thank Claire Treat, Nigel Roulet, and another anonymous reviewer for their helpful comments to improve the manuscript.
APPENDIX: SITES CURRENTLY INCLUDED IN THE DATABASE.
Table A1 presents characteristics of sites currently included in the database.
A supplement to this article is available online (10.1175/BAMS-D-18-0268.2)