Abstract. Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20∘ S to 20∘ N) the spring onset of elevated CH4 emissions starts 3 d earlier, and the CH4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

Global and regional projections of climate change by Earth system models are limited by their uncertain estimates of terrestrial ecosystem productivity. At the middle to low latitudes, the East Asian monsoon region has higher productivity than forests in Europe‐Africa and North America, but its estimate by current generation of terrestrial biosphere models (TBMs) has seldom been systematically evaluated. Here, we developed a traceability framework to evaluate the simulated gross primary productivity (GPP) by 15 TBMs in the East Asian monsoon region. The framework links GPP to net primary productivity, biomass, leaf area and back to GPP via incorporating multiple vegetation functional properties of carbon‐use efficiency (CUE), vegetation C turnover time (τveg), leaf C fraction (Fleaf), specific leaf area (SLA), and leaf area index (LAI)‐level photosynthesis (PLAI), respectively. We then applied a relative importance algorithm to attribute intermodel variation at each node. The results showed that large intermodel variation in GPP over 1901–2010 were mainly propagated from their different representation of vegetation functional properties. For example, SLA explained 77% of the intermodel difference in leaf area, which contributed 90% to the simulated GPP differences. In addition, the models simulated higher CUE (18.1 ± 21.3%), τveg (18.2 ± 26.9%), and SLA (27.4±36.5%) than observations, leading to the overestimation of simulated GPP across the East Asian monsoon region. These results suggest the large uncertainty of current TBMs in simulating GPP is largely propagated from their poor representation of the vegetation functional properties and call for a better understanding of the covariations between plant functional properties in terrestrial ecosystems.