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.

• We evaluate methane flux gap-filling methods across 17 boreal-to-tropical wetlands • New methods for generating realistic artificial gaps and uncertainties are proposed • Decision tree algorithms perform slightly better than neural networks on average • Soil temperature and generic seasonality are the most important predictors • Open-source code is released for gap-filling steps and uncertainty evaluation Time series of wetland methane fluxes measured by eddy covariance require gap-filling to estimate daily, seasonal, and annual emissions. Gap-filling methane fluxes is challenging because of high variability and complex responses to multiple drivers. To date, there is no widely established gap-filling standard for wetland methane fluxes, with regards both to the best model algorithms and predictors. This study synthesizes results of different gap-filling methods systematically applied at 17 wetland sites spanning boreal to tropical regions and including all major wetland classes and two rice paddies. Procedures are proposed for: 1) creating realistic artificial gap scenarios, 2) training and evaluating gap-filling models without overstating performance, and 3) predicting half-hourly methane fluxes and annual emissions with realistic uncertainty estimates. Performance is compared between a conventional method (marginal distribution sampling) and four machine learning algorithms. The conventional method achieved similar median performance as the machine learning models but was worse than the best machine learning models and relatively insensitive to predictor choices. Of the machine learning models, decision tree algorithms performed the best in cross-validation experiments, even with a baseline predictor set, and artificial neural networks showed comparable performance when using all predictors. Soil temperature was frequently the most important predictor whilst water table depth was important at sites with substantial water table fluctuations, highlighting the value of data on wetland soil conditions. Raw gap-filling uncertainties from the machine learning models were underestimated and we propose a method to calibrate uncertainties to observations. The python code for model development, evaluation, and uncertainty estimation is publicly available. This study outlines a modular and robust machine learning workflow and makes recommendations for, and evaluates an improved baseline of, methane gap-filling models that can be implemented in multi-site syntheses or standardized products from regional and global flux networks (e.g., FLUXNET).

• 514 sites-years of flux data were used to analyze the potential of physiological and phenological metrics in explaining the variability of forest NEP; • Summer physiological metrics performed better than phenological metrics in explaining IAV of NEP; • Ecosystem respiration played an important role in controlling the variability of NEP in forest ecosystem; • MCUI exhibited a great potential in explaining both IAV and SV of NEP. Understanding the feedback of ecosystem carbon uptake on climate change at temporal and spatial scales is crucial for developing ecosystem models. Previous studies have focused on the role of spring and autumn phenology in regulating carbon sequestration in forest stands, but few on the impact of physiological status in summer. However, plant accumulated the most carbon in summer compared with spring and autumn, therefore, it is of great significance to explore the role of summer phenological metrics on the variability of carbon sequestration. Using 514 site-years of flux data obtained at 40 FLUXNET sites including three forest ecosystems (i.e. evergreen needleleaf forest (ENF), deciduous broadleaf forest (DBF) and mixed forest (MF)) in Europe and North America, we compared the potential of physiological and phenological metrics of Gross Primary Production (GPP) and Ecosystem Respiration (RECO) in explaining the interannual and spatial variability (IAV and SV) of forest net ecosystem production (NEP). In view of the better performance of physiological metrics, we developed the maximum carbon uptake index (MCUI), which integrated the physiology metrics of photosynthesis and respiration in summer, and further explored its ability in explaining the IAV and SV of NEP. The results suggest that the MCUI had a better ability than respiration-growth length ratio (RGR) in predicting NEP for all three forest types. The interpretation of MCUI based on meteorological variables illustrated that the controlling meteorological factors of MCUI differed substantially among ecosystems. The summer shortwave radiation had the greatest influence on MCUI at DBF sites, while the soil water content played an important but opposite role at ENF and DBF sites, and no significant meteorological driver was found at MF sites. The higher potential of MCUI in explaining IAV and SV of NEP highlights the importance of summer physiology in controlling the forest carbon sequestration, and further confirms the significant role of peak plant growth in regulating carbon cycle of forest ecosystems. Understanding the drivers of peak plant growth is therefore of a great significance for further improving the precious of ecosystem model in the future.

A Correction to this paper has been published: https://doi.org/10.1038/s41597-021-00851-9.

Abstract The FLUXNET2015 dataset provides ecosystem-scale data on CO 2 , water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Abstract Terrestrial evapotranspiration (ET) is thermodynamically expected to increase with increasing atmospheric temperature; however, the actual constraints on the intensification of ET remain uncertain due to a lack of direct observations. Based on the FLUXNET2015 Dataset, we found that relative humidity (RH) is a more important driver of ET than temperature. While actual ET decrease at reduced RH, potential ET increases, consistently with the complementary relationship (CR) framework stating that the fraction of energy not used for actual ET is dissipated as increased sensible heat flux that in turn increases potential ET. In this study, we proposed an improved CR formulation requiring no parameter calibration and assessed its reliability in estimating ET both at site-level with the FLUXNET2015 Dataset and at basin-level. Using the ERA-Interim meteorological dataset for 1979–2017 to calculate ET, we found that the global terrestrial ET showed an increasing trend until 1998, while the trend started to decline afterwards. Such decline was largely associated with a reduced RH, inducing water stress conditions that triggered stomatal closure to conserve water. For the first time, this study quantified the global-scale implications of changes in RH on terrestrial ET, indicating that the temperature-driven acceleration of the terrestrial water cycle will be likely constrained by terrestrial vegetation feedbacks.

Forests dominate carbon (C) exchanges between the terrestrial biosphere and the atmosphere on land. In the long term, the net carbon flux between forests and the atmosphere has been significantly impacted by changes in forest cover area and structure due to ecological disturbances and management activities. Current empirical approaches for estimating net ecosystem productivity (NEP) rarely consider forest age as a predictor, which represents variation in physiological processes that can respond differently to environmental drivers, and regrowth following disturbance. Here, we conduct an observational synthesis to empirically determine to what extent climate, soil properties, nitrogen deposition, forest age and management influence the spatial and interannual variability of forest NEP across 126 forest eddy-covariance flux sites worldwide. The empirical models explained up to 62% and 71% of spatio-temporal and across-site variability of annual NEP, respectively. An investigation of model structures revealed that forest age was a dominant factor of NEP spatio-temporal variability in both space and time at the global scale as compared to abiotic factors, such as nutrient availability, soil characteristics and climate. These findings emphasize the importance of forest age in quantifying spatio-temporal variation in NEP using empirical approaches.

Abstract Research infrastructures play a key role in launching a new generation of integrated long-term, geographically distributed observation programmes designed to monitor climate change, better understand its impacts on global ecosystems, and evaluate possible mitigation and adaptation strategies. The pan-European Integrated Carbon Observation System combines carbon and greenhouse gas (GHG; CO 2 , CH 4 , N 2 O, H 2 O) observations within the atmosphere, terrestrial ecosystems and oceans. High-precision measurements are obtained using standardised methodologies, are centrally processed and openly available in a traceable and verifiable fashion in combination with detailed metadata. The Integrated Carbon Observation System ecosystem station network aims to sample climate and land-cover variability across Europe. In addition to GHG flux measurements, a large set of complementary data (including management practices, vegetation and soil characteristics) is collected to support the interpretation, spatial upscaling and modelling of observed ecosystem carbon and GHG dynamics. The applied sampling design was developed and formulated in protocols by the scientific community, representing a trade-off between an ideal dataset and practical feasibility. The use of open-access, high-quality and multi-level data products by different user communities is crucial for the Integrated Carbon Observation System in order to achieve its scientific potential and societal value.

Abstract The Integrated Carbon Observation System is a Pan-European distributed research infrastructure that has as its main goal to monitor the greenhouse gas balance of Europe. The ecosystem component of Integrated Carbon Observation System consists of a multitude of stations where the net greenhouse gas exchange is monitored continuously by eddy covariance measurements while, in addition many other measurements are carried out that are a key to an understanding of the greenhouse gas balance. Amongst them are the continuous meteorological measurements and a set of non-continuous measurements related to vegetation. The latter include Green Area Index, aboveground biomass and litter biomass. The standardized methodology that is used at the Integrated Carbon Observation System ecosystem stations to monitor these vegetation related variables differs between the ecosystem types that are represented within the network, whereby in this paper we focus on forests, grasslands, croplands and mires. For each of the variables and ecosystems a spatial and temporal sampling design was developed so that the variables can be monitored in a consistent way within the ICOS network. The standardisation of the methodology to collect Green Area Index, above ground biomass and litter biomass and the methods to evaluate the quality of the collected data ensures that all stations within the ICOS ecosystem network produce data sets with small and similar errors, which allows for inter-comparison comparisons across the Integrated Carbon Observation System ecosystem network.