# The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data

Gilberto Pastorello, Carlo Trotta, Eleonora Canfora, Housen Chu, Danielle Christianson, You-Wei Cheah, Cristina Poindexter, Jiquan Chen, Abdelrahman Elbashandy, Marty Humphrey, et al.

► To cite this version:

Gilberto Pastorello, Carlo Trotta, Eleonora Canfora, Housen Chu, Danielle Christianson, et al. The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data. Scientific Data, 2020, 7 (1), pp.1-27. 10.1038/s41597-020-0534-3. hal-03778635

HAL Id: hal-03778635

https://hal.science/hal-03778635

Submitted on 16 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC DATA 110110
0111101
11011110
011101101

<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

# The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data

Gilberto Pastorello et al. #

The FLUXNET2015 dataset provides ecosystem-scale data on CO₂ 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.

## Background & Summary

For over three decades, the eddy covariance technique¹ has been used to measure land-atmosphere exchanges of greenhouse gases and energy at sites around the world to study and determine the function and trajectories of both ecosystems and the climate system. The technique allows nondestructive measurement of these fluxes at a high temporal resolution and ecosystem level, making it a unique tool. Based on high frequency (10–20 Hz) measurements of vertical wind velocity and a scalar (CO₂, H₂O, temperature, etc.), it provides estimates of the net exchange of the scalar over a source footprint area that extends up to hundreds of meters around the point of measurement. Soon after the first consistent measurement sites were operational, regional networks of sites were formed in Europe² and the US³,⁴, followed by similar initiatives in other continents⁵⁻⁷. Networks enabled the use of eddy covariance data beyond a single site or ecosystem for cross-site comparisons and regional-to-global studies⁸⁻¹⁵. These regional networks have evolved into long-term research infrastructures or monitoring activities, such as ICOS, AmeriFlux, NEON, AsiaFlux, ChinaFLUX, and TERN-OzFlux.

FLUXNET was created as a global network of networks¹⁶⁻¹⁸, a joint effort among regional networks to harmonize and increase standardization of the data being collected. It made possible the creation of global eddy covariance datasets. The first gap-filled, global FLUXNET dataset, which included derived, partitioned fluxes like photosynthesis and respiration, was the Marconi dataset¹⁹ in 2000, with 97 site-years of data, followed by the 2007 LaThuile dataset²⁰ with 965 site-years of data, and finally in 2015 the FLUXNET2015 dataset¹⁸,²¹ (hereafter FLUXNET2015) with 1532 site-years of data. Two main factors limited the numbers of sites and years included in each dataset: data policy and data quality. Willingness to share data under the selected data policy is a major reason why FLUXNET2015 likely only includes between 10–20% of existing sites globally—the total number of existing sites is still unknown. Then there is the evolution of processing pipelines and quality controls, leading to new issues being identified in the data that, if not solved in time, led to leaving out that data. The LaThuile dataset had a more restrictive policy and, in a few cases, previously undiscovered data issues, leading to fewer sites being included in FLUXNET2015.

<footnote># A full list of authors and their affiliations appears at the end of the paper.</footnote>
<figure><img src="image_3.png" /></figure>

**Fig. 1** Map of 206 tower sites included in this paper from the 212 sites in the February 2020 release of the FLUXNET2015 dataset. The size of the circle indicates the length of the data record. The color of the circles represents the ecosystem type based on the International Geosphere-Biosphere Programme (IGBP) definition. When overlapping, locations are offset slightly to improve readability. Numbers in parentheses indicate the number of sites in each IGBP group. The inset shows the distribution of data record lengths. See also Supplementary Fig. SM4 for continental scale maps of Australia, Europe, and North America.

FLUXNET2015 includes for the first time sites with records over two decades long (Fig. 1). The dataset was created through collaborations among many regional networks, with data preparation efforts happening at site, regional network, and global network levels. The global coordination of data preparation activities and data processing was done by a team from the AmeriFlux Management Project (AMP), the European Ecosystem Fluxes Database, and the ICOS Ecosystem Thematic Centre (ICOS-ETC). This team was responsible for the coding efforts, quality checks, and execution of the data processing pipeline. These combined efforts led to a dataset that is standardized with respect to the (1) data products themselves, (2) data distribution formatting, and (3) data quality across sites. The wide application of these datasets in global synthesis and modeling activities highlights their value. At the same time, however, heterogeneity in the data-caused mainly by differences in data collection, flux calculations, and data curation before submission-highlights the need for estimates of data uncertainty and uniform evaluation of data quality.

The data processing pipeline uses well-established and published methods, with new code implemented for this release as well as code adapted from implementations by the community. The main products in this pipeline are: (1) thorough data quality control checks; (2) calculation of a range of friction velocity thresholds to filter low turbulence periods, allowing an estimate of the uncertainty from this filtering along with the random uncertainty; (3) gap-filling of meteorological and flux measurements, including the use of a downscaled reanalysis data product to fill long gaps in meteorological variables; (4) partitioning of CO₂ fluxes into respiration and photosynthesis (gross primary productivity) components using three distinct methods; and (5) calculation of a correction factor for energy fluxes estimating the deviation from energy balance closure for the site. Two features of this pipeline are the ranges of friction velocity thresholds, and the multiple methods for partitioning CO₂ fluxes. Both features support a more thorough evaluation of the uncertainty introduced by the processing steps themselves. Our implementation of this pipeline is available as an open-source code package called ONEFlux (Open Network-Enabled Flux processing pipeline)²². The goal of this paper is to describe FLUXNET2015 and additional products, present the details about this processing pipeline, and document the methods used to generate the dataset. Doing so will provide the community of FLUXNET end-users with the technical and practical knowledge necessary to harness the full potential of the FLUXNET data, including data from the FLUXNET2015 release, and data submitted to the network since.

## Data Processing Methods

The data contributed by site teams for inclusion in FLUXNET2015 encompassed fluxes, meteorological, environmental, and soil time series at half-hourly or hourly resolutions. Contributed data underwent a uniform data quality control process, with issues addressed in consultation with site teams. Data were then processed using the pipeline (Fig. 2) described in this section. The resulting data products were distributed through the FLUXNET-Fluxdata web portal²³, where the usage of the dataset is tracked through a registration of all the
<figure><img src="image_4.png" /><figcaption>**Fig. 2** The logic of the data processing steps for FLUXNET2015 (details about the different steps and meaning of abbreviations in the text).</figcaption></figure>

requests and details about the user and the data use plan. This information is crucial to better understand the user needs and the impact of the dataset.

**Data sources.** The first level of processing was completed by the site teams, including the calculation of half-hourly or hourly turbulent fluxes from high-frequency wind and concentration measurements, the averaging of meteorological variables sampled at shorter temporal resolutions, and the site team's own quality control procedures. The contributed fluxes were required to be submitted separately as turbulent and storage fluxes (components to be added for total flux), and not gap-filled or filtered for low turbulence conditions-see Aubinet et al.²⁴ for details. Control checks were implemented to ensure that this was the case, and to detect additional inconsistencies that site teams were asked to address, for example coherence of the timestamps, consistency across correlated variables, etc²⁵. Starting from these data, we applied the same set of methods for all the remaining processing steps (from filtering to gap-filling and partitioning), increasing uniformity, and allowing quantification of the uncertainty introduced by the processing. The harmonization of data-especially of data quality-was a high priority while creating the dataset, and therefore extensive interaction with the site teams was necessary. Data were contributed through a regional network, and the formats provided by the networks were converted into a standard input for processing. This FLUXNET dataset led to the creation of a new cross-network specification for standard site data and metadata submission formats, now adopted by different regional networks.

**Data processing pipeline overview.** The data processing pipeline (Fig. 2) is divided into four main processing blocks. The first one is the data quality assurance and quality control (QA/QC) activities our team applied to data from all sites. This part was done with a combination of automated procedures and manual checks for all the variables in the contributed data. The next three blocks were all part of an automated pipeline that was executed separately for each site: the Energy & Water Fluxes processing block encompasses sensible and latent heat variables; the Carbon Fluxes processing block handles the variables for CO₂ fluxes like net ecosystem exchange (NEE); and, the Meteorological Variables processing block deals with all the meteorological measurements that are also used in the processing of fluxes and other products. At each processing step, a set of automated pre- and post-conditions are enforced, making sure the inputs and outputs of each step are within the expected behavior. The final step involved merging all the products generated at previous steps, and adding daily through yearly temporal aggregations of most of the variables in the dataset and related quality flags. At this step, automated checks were performed on all the variables to ensure consistency, and the final files with all the contributed and derived data products to be distributed are created. Supplementary Fig. SM2 shows the general steps involved in the processing in the sequence organized in the code, as available in the ONEFlux package²². All the steps have been implemented in the shared code except the Sundown partitioning (see Implementation Approach section for details).

**Ensemble of results.** The adoption of multiple methods for the same step (e.g., two methods for USTAR threshold calculation, or two or three methods for CO₂ flux partitioning, see below in this section for details) is motivated by the existence of different methods in the literature, using different assumptions and potentially
having diverging results²⁶,²⁷. On the one hand, this lack of uniformity can represent a problem for synthesis studies. On the other hand, adopting a single method could lead to biases and underestimation of the uncertainty in the methodology. The approach taken here, which simultaneously adopts multiple methods, allows the creation of an ensemble, helping assess uncertainty, and also the suitability of individual methods to a site's conditions.

**Data quality assurance and quality control.** Prior to the processing that generated the derived data products (hereafter called *post-processing*), the data for each site went through QA/QC checks following Pastorello *et al.*²⁵. All variables included in the dataset underwent checks, and critical variables underwent further scrutiny. These additional checks targeted variables critical to the processing, e.g., flux variables, meteorological variables used in gap-filling, and variables used by uncertainty estimation procedures. The processing did not proceed for sites with pending issues in critical variables.

**Critical metadata variables for the post-processing:**

* The site FluxID in the form CC-SSS (two character country code, three character site identifier within country) – e.g., US-Ha1

* The latitude and longitude for the site in the WGS 84 decimal format with at least four decimal points resolution – e.g., 42.5378/-72.1715

* Time zone of the site (time series, if time zone changed; timestamps are all local standard time, no daylight savings) – e.g., UTC-5

* Height of the gas analyzer – e.g., 30.0 m

**Critical data variables for the post-processing, averaged or integrated over 30 or 60 minutes (*required*):**

* *CO₂ (µmolCO₂ mol⁻¹): Carbon Dioxide (CO₂) mole fraction in moist air

* *FC (µmolCO₂ m⁻² s⁻¹): Carbon Dioxide (CO₂) turbulent flux (without storage component)

* *SC (µmolCO₂ m⁻² s⁻¹): Carbon Dioxide (CO₂) storage flux measured with a vertical profile system, optional if tower shorter than 3 m

* *H (W m⁻²): sensible heat turbulent flux, without storage correction

* *LE (W m⁻²): latent heat turbulent flux, without storage correction

* *WS (m s⁻¹): horizontal wind speed

* *USTAR (m s⁻¹): friction velocity

* *TA (deg C): air temperature

* *RH (%): relative humidity (range 0–100%)

* *PA (kPa): atmospheric pressure

* *G (W m⁻²): ground heat flux, not mandatory, but needed for the energy balance closure calculations

* *NETRAD (W m⁻²): net radiation, not mandatory, but needed for the energy balance closure calculations

* *SW_IN (W m⁻²): incoming shortwave radiation

* *SW_IN_POT (W m⁻²): potential incoming shortwave radiation (top of atmosphere theoretical maximum radiation), calculated based on the site coordinates²²

* *PPFD_IN (µmolPhotons m⁻² s⁻¹): incoming photosynthetic photon flux density

* *P (mm): precipitation total of each 30 or 60 minute period

* *LW_IN (W m⁻²): incoming (down-welling) longwave radiation

* *SWC (%): soil water content (volumetric), range 0–100%

* *TS (deg C): soil temperature

**File format standardization.** To process hundreds of sites, we needed consistent file formats that supported the input data and metadata. This led to multi-network agreements and creation of formats for data and metadata contribution to the regional networks²⁸,²⁹. These formats have now been adopted by networks in Europe and the Americas and by some instrument manufacturers, and are under consideration by other regional networks. In addition, automated extraction and conversion tools for direct format translation were implemented to work with data in older formats.

**Data QA/QC steps.** Data quality was checked before the processing started. If issues were identified that could not be resolved by the network-level data team, the site team was asked to suggest a course of action or send a new version of the data addressing the quality issue identified. The main data QA/QC steps were: single-variable checks, multi-variable checks, specialized checks, and automatic checks. Single-variable checks look at patterns in the time series of one variable at a time, for long- and short-term trends and other issues. Multi-variable checks look at the relationships among correlated variables (e.g., different radiation variables) to identify periods with disagreements. Specialized checks test for common issues in EC and meteorological data, like timestamp shifts or sensor deterioration. During this phase, a time series of top-of-the-atmosphere potential radiation (SW_IN_POT) is also computed, using latitude/longitude coordinates and time²². These three types of checks are detailed in Pastorello *et al.*²⁵. The automated checks apply variable-specific despiking routines adapted from Papale *et al.*³⁰ and apply a set of range controls per variable. This last step creates a series of flags that were discussed with the site managers for corrections and resubmissions and then used to filter the data in subsequent steps.

**Meteorological products.** The main processing applied to meteorological data was gap-filling by two independent methods: Marginal Distribution Sampling³¹ (MDS) and ERA-Interim³². Data gap-filled by MDS (applied to all variables that are gap-filled) are identified by the _F_MDS suffix. Data gap-filled using ERA-Interim
downscaling (six variables that are available in the reanalysis dataset) have a _F_ERA suffix. The final gap-filled time series for variables combines both of these methods (indicated by an unqualified _F_ suffix), following a data-quality-based selection approach (see below). For SW_IN, in case of gaps or in case the variable was not measured, we performed a calculation from PPFD_IN when available, calculating the conversion factor from the periods of overlap of the two measurements (and assuming a factor 0.48 J (µmol photon)⁻¹ when the sensors did not run in parallel at the site).

**MDS.** The MDS method, introduced in Reichstein et al.³¹, is applied to all variables that may be gap-filled. It works by seeking meteorological conditions physically and temporally similar to the ones for the missing data point(s). The restrictions on the size of the time window and which variables must be available are incrementally relaxed until a suitable set of records is found to fill the gap in the target variable. All values of the target variable satisfying the current set of conditions are averaged to generate the fill value. The method was applied as described in the original implementation³¹ using SW_IN, TA and VPD as drivers. The basic three scenarios for the time when the target variable is missing are: (i) all three drivers are available; (ii) only SW_IN is available; and (iii) all three drivers are missing. Based on the available co-located variables, a search for similar conditions is started, keeping the searching window as small as possible to avoid changes in other slow-changing drivers (phenology, water availability, etc.). The more variables missing and the larger the time window, the lower the confidence indicated by the _F_MDS_QC_ flag. The values for this flag are (0-3): _F_MDS_QC=0_ (measured); _F_MDS_QC=1_ (filled with high confidence); _F_MDS_QC=2_ (filled with medium confidence); _F_MDS_QC=3_ (filled with low confidence). For details on the implementation, see the original paper³¹ and the ONEFlux source code²².

**ERA-Interim.** This method is based on the ERA-Interim (ERA-I) Reanalysis global atmospheric product³³,³⁴ created by the European Centre for Medium-Range Weather Forecasts (ECMWF)-ERA stands for ECMWF Re-Analysis. Applied to the subset of variables that are also available in the ERA-Interim product, the method involves a spatial and temporal downscaling process using the measured variable at the site. The ERA-I variables that were used are: air temperature at 2 m (t2m, K), incoming shortwave solar radiation at the surface (Sw, W m⁻²), dew point temperature at 2 m (dt2m, K), wind speed horizontal components at 10 m (u10 and v10, m s⁻¹), total precipitation (Pr, m of water per time step), and incoming longwave solar radiation at the surface (Lw, W m⁻²). The gap-filling procedure harmonizes units, identifies periods that are long enough to allow a linear relationship to be built, a simple debiasing of the linear relationship, evaluation of the diurnal cycle in the subset of variables, and other evaluations of the results to identify potential missing or incorrect information (e.g., coordinate or temporal mismatches). The linear relationships are built taking into account instantaneous and averaged variables, and then applied to the whole ERA-I record, generating the spatially (coordinate-based) and temporally (diurnal cycle-based) downscaled version of each variable. The method was applied as in the original implementation; additional details can be found in Vuichard and Papale³².

**Final gap-filled product.** Measured or high quality gap-filled records using MDS (_F_MDS_QC<2_) are used in the final gap-filled products (_F_ suffixed variables, without _MDS_ or _ERA_). If the variable has a low quality gap-fill flag (2 or 3), the ERA-I product is used instead. The final quality flag (_F_QC_) is 0 for measured, 1 for high quality fill using MDS, and 2 for data gap-filled with the ERA-I downscaled product. A gap-filled version of CO₂ concentration is also generated (CO₂_F_MDS) using the MDS method as described above, including the corresponding quality flag.

**Energy and water products.** The main data products associated with energy and water fluxes are the gap-filled versions of the data and the estimation of a version ensuring the energy balance closure and estimating its uncertainty—for a description of the issue see Stoy et al.³⁵ Turbulent energy fluxes (sensible and latent heat, H and LE, respectively) are gap-filled using the MDS method³¹ described above. From LE, it is possible to calculate the water flux (evapotranspiration) using the latent heat of vaporization. An energy balance corrected version of LE and H is also created, a data product often needed when data are used in model parameterization and validation for which the closure of the energy balance is prescribed. There is no general agreement on the reasons and approaches to correct the imbalance in the energy budget within EC measurements. In this product, the methodology used to calculate the energy balance corrected fluxes is based on the assumption that the Bowen ratio is correct³⁶. Fluxes are corrected by multiplying the original, gap-filled LE and H data by an energy balance closure correction factor (EBC_CF, in the dataset). The correction factor is calculated starting from the half-hours where all the variables needed to estimate energy balance closure are available (measured NETRAD and G, and measured or good-quality gap-filled H and LE). The correction factor for each single half-hour is calculated as in Eq. (1), but is not applied directly.

$$EBC\_CF = (NETRAD - G)/(H + LE) \qquad (1)$$

First, to avoid transient conditions, the calculated EBC_CF time series is filtered by removing values outside of 1.5 times its own interquartile range. Then, the correction factor used in the calculations is obtained using one of three methods, applied hierarchically (see also diagram in Supplementary Fig. SM3):

* **EBC_CF Method 1:** For each half-hour, a sliding window of ±15 days (31 days total) is used to select half-hours between time periods 22:00-02:30 and 10:00-14:30 (local standard time). These time-of-day restrictions aim at removing sunrise and sunset time periods, when changes in ecosystem heat storage (not measured) are more significant, preventing energy balance closures. For all half-hours meeting these criteria, the corresponding EBC_CFs are selected and used to calculate the corrected values of H and LE for the half-hour processed (center of the sliding window), generating a pool of values for each of these two variables.
From each of these two pools, the 25th, 50th (median), and 75th percentiles are extracted for their corresponding variables, generating the values for H_CORR25, H_CORR, H_CORR75 and LE_CORR25, LE CORR, LE_CORR75. If fewer than five EBC_CF values are present in the sliding window, Method 2 is used for the half-hour. (Note on temporal aggregations: for DD the sliding window size is ±7 days and the EBC_CF are calculated from the daily average values of G, NETRAD, H and LE. For WW, MM, and YY the EBC_CFs are calculated from corresponding average fluxes of the period analysed, but no percentiles are computed. For WW, MM, and YY, Method 1 fails if less than 50% of half-hours within the window have measured values for all four component variables.)

* **EBC_CF Method 2:** For the current half-hour, EBC_CF is calculated as the average of the EBC_CF values used to calculate the H_CORR and LE_CORR with Method 1 within a sliding window of ±5 days and ±1 hour of the time-of-day of the current timestamp. H_CORR and LE_CORR are calculated and the corresponding CORR25 and CORR75 percentiles are not generated. If no EBC_CF is available, Method 3 is used for the current half-hour. (Note on temporal aggregations: differing sliding windows are: DD: ±2 weeks, WW: ±2 weeks, MM: ±1 month, and YY: ±1 year.)

* **EBC_CF Method 3:** An approach like Method 2 is applied but using a sliding window of ±5 days for the same half-hour in the previous and next years, with the current EBC_CF being calculated from the average of the EBC_CF values used to calculate the H_CORR and LE_CORR. H_CORR and LE_CORR are calculated and the corresponding CORR25 and CORR75 percentiles are not generated. In case this method also cannot be applied due to missing values, the energy balance closure corrected fluxes are not calculated. (Note on temporal aggregations: differing sliding windows are: DD: ±2 weeks, WW: ±2 weeks, MM: ±1 month, and YY: ±2 years.)

**H and LE Random Uncertainty.** The random uncertainty for H and LE is also estimated at half-hourly resolution, based on the method introduced by Hollinger & Richardson³⁷ and then aggregated at the other temporal resolutions. The random uncertainty (indicated by the suffix _RANDUNC) in the measurements is estimated using one of two methods, applied hierarchically:

* **H-LE-RANDUNC Method 1 (direct standard deviation method):** For a sliding window of ±5 days and ±1 hour of the time-of-day of the current timestamp, the random uncertainty is calculated as the standard deviation of the measured fluxes. The similarity in the meteorological conditions evaluated as in the MDS gap-filling method³¹ and a minimum of five measured values must be present; otherwise, method 2 is used.

* **H-LE-RANDUNC Method 2 (median standard deviation method):** For the same sliding window of ±5 days and ±1 hour of the time-of-day of the current timestamp, random uncertainty is calculated as the median of the random uncertainty (calculated with H-LE-RANDUNC Method 1) of similar fluxes, i.e., within the range of ±20% and not less than 10 W m⁻².

The joint uncertainty for H and LE is computed from the combination of the uncertainty from the energy balance closure correction factor and random uncertainty.

$$H\_CORR\_JOINTUNC = \sqrt{H\_RANDUNC^2 + \left( \frac{H\_CORR75 - H\_CORR25}{1.349} \right)^2} \quad (2)$$

These variables are identified by the _JOINTUNC suffix and are computed for H as in Eq. (2), and similarly for LE. (Note on temporal aggregations: joint uncertainties for H and LE are recomputed at HH and DD resolutions separately, and not generated for WW, MM, and YY resolutions.)

**CO₂ products.** The processing steps applied to CO₂ fluxes were: calculation of net ecosystem exchange (NEE) from CO₂ turbulent and storage fluxes, applying a spike detection algorithm, filtering for low turbulence conditions using multiple friction velocity (USTAR) thresholds, gap-filling of all NEE time series generated by an ensemble of USTAR thresholds, estimation of random uncertainty, and partitioning of NEE into its ecosystem respiration (RECO) and gross primary production (GPP) components.

**Calculation of NEE.** CO₂ storage fluxes (SC) express the change of CO₂ concentration below the measurement level of the eddy covariance system within the half-hour. NEE was calculated as the sum of the CO₂ turbulent fluxes (FC) and SC. Both FC and SC are part of the required data contributed by site teams. SC is usually estimated using a profile system³⁸. If SC was not provided or missing, two cases were implemented: for measurement heights lower than 3 m and short canopies, the SC term was considered to be 0; for taller towers/canopies, a discrete estimation based on the top measurement of CO₂ concentration was used to compute SC³⁹.

**Despiking of NEE.** Although the processing of high frequency data into half-hourly fluxes usually includes steps to remove spikes from instantaneous measurements, spikes can also occur in the half-hourly data. The method described in Papale et al.³⁰, based on the median absolute deviation (MAD) with z = 5.5, was applied to filter NEE for residual spikes that were removed.

**USTAR threshold estimation and filtering.** Filtering for low turbulence conditions is necessary when there is not enough turbulence, causing the ecosystem flux to be transported by advective flows and missed by both the eddy covariance system and the storage profile, resulting in underestimated fluxes. Despite different approaches having been tested to measure and quantify horizontal and vertical advection⁴⁰, the most often used method to avoid the
underestimation of fluxes is removing the data points potentially affected by strong advection¹. These points are identified using the friction velocity (USTAR) as an indicator of turbulence strength, defining a threshold value under which NEE measurements are discarded and replaced by gap-filled estimates.

This USTAR threshold is linked to the canopy structure, measurement height, wind regimes, and other factors specific to an individual site. It is estimated using nighttime NEE measurements (only ecosystem respiration), based on the dependency between USTAR and NEE at similar temperatures and periods of the year (main drivers of ecosystem respiration). Under these conditions, NEE is assumed (and expected) to be independent from USTAR, which is not a driver of respiration. However, in most sites below a certain USTAR threshold value, NEE is found to increase with USTAR; this USTAR value is selected as the threshold to define conditions with reduced risk of flux underestimation. Different methods have been proposed to estimate the USTAR thresholds and the related uncertainty as to how the approach works at a specific site¹.

**CP and MP USTAR threshold methods.** Two methods to calculate USTAR thresholds were used: change-point-detection (CP) proposed by Barr *et al.*⁴¹ and a modified version of the moving-point-transition (MP) described originally by Reichstein *et al.*³¹ and Papale *et al.*³⁰. Both methods are similar in terms of data selection, preparation and grouping and aim to estimate the USTAR threshold value. Measurements collected when USTAR is below the threshold are removed. The difference between these methods is in how this threshold value is estimated. For both methods, the nighttime data of a full year are divided in four three-month periods (seasons) and 7 temperature classes (of equal size in terms of number of observations). For each season/temperature group the data are divided into 20 USTAR classes (also with equal number of observations) and the average NEE for each USTAR class is computed. The calculation of the threshold uses each of the methods (see below for details on their differences). For each season, the median value of the 7 temperature classes is calculated and a final threshold is defined by selecting the maximum of the 4 seasonal values.

The CP method uses two linear regressions between NEE and USTAR, the second with an imposed zero slope. The change point is defined as where the two lines cross, i.e., constraining the shape of the NEE-USTAR dependency. The method is extensively used to detect temporal discontinuities in climatic data. Details can be found in Barr *et al.*⁴¹.

For the MP method³⁰,³¹, the mean NEE value in each of the 20 USTAR classes is compared to the mean NEE measured in the 10 higher USTAR classes. The threshold selected is the USTAR class in which the average nighttime NEE reaches more than 99% of the average NEE at the higher USTAR classes. An improvement of the MP method was implemented here for robustness over noisy data, by adding a second step to the original MP implementation: when a threshold is selected, it was tested to ensure it was also valid for the following USTAR class. In other words, assuming that Eq. (3) holds, where *x* is one of the 20 USTAR classes and *NEE_USTAR*(*x*) is the average NEE for that USTAR class.

$$NEE_{USTAR}(x) > 0.99 \times \text{MEAN}(NEE_{USTAR}(x + 1), NEE_{USTAR}(x + 2), \dots, NEE_{USTAR}(x + 10)) \quad (3)$$

The USTAR value associated to the *x*th-class was selected as threshold only if Eq. (4) also holds, to confirm that the plateau where NEE is USTAR-independent was reached. If not, the search for the plateau and threshold continued toward higher USTAR values.

$$NEE_{USTAR}(x + 1) > 0.99 \times \text{MEAN}(NEE_{USTAR}(x + 2), NEE_{USTAR}(x + 3), \dots, NEE_{USTAR}(x + 11)) \quad (4)$$

**Bootstrapping USTAR threshold estimation.** For each of the two methods, a bootstrapping technique was used. The full dataset (year of measurement) was re-sampled 100 times with the possibility to select the same data point multiple times (i.e., with replacement), creating 100 versions of the dataset. The threshold values were calculated for each of them, obtaining 100 threshold values per method (CP and MP) and year, for a total of 200 USTAR threshold estimates for each year. This process and next steps are illustrated in Fig. 3. These 200 threshold values represent the uncertainty in the threshold estimation that could also impact the uncertainty of NEE. It is worth noting that there is not always a direct relationship between the threshold and NEE uncertainties. It is possible, for instance, that a small variability in the thresholds has a strong effect on NEE or, conversely, with NEE almost insensitive to the threshold value. This is related to the site characteristics (USTAR variability) and to the level of difficulty in filling the gaps created by the filtering.

There are cases where not enough data are present to calculate a USTAR threshold (for both the CP and MP methods) or where it is not possible to identify a clear change point (CP method only). This leads to the uncertainty being underestimated (fewer or no USTAR threshold values available). This should be considered as a general indication of difficulties in the application of the USTAR filtering for the specific sites or years. Sites and years where these conditions occurred are reported in the SUCCESS_RUN variable in the AUXNEE product (values 1: threshold found, 0: failed/no threshold found).

**Variable and constant variants of the USTAR threshold methods.** To calculate the uncertainty in NEE due to the uncertainty in the selected USTAR threshold, all the threshold values obtained with the two methods and the bootstrapping were pooled together, from which 40 representative values were extracted: from the percentile 1.25 of the series to the percentile 98.75, with a step of 2.5, i.e., [1.25:2.5:98.75]. When long time series (multi-years) are processed, it is possible to extract the 40 representative thresholds for each of the years. The threshold is
<figure><img src="image_5.png" /><figcaption>**Fig. 3** To identify and remove data collected under low turbulence conditions, under which advective fluxes could lead to an underestimation of fluxes, filtering based on the USTAR threshold was used. In order to estimate the uncertainty in the USTAR threshold calculation, a bootstrapping approach was implemented, with a selection of values representative of the distribution included in the final data products. From the (up to) 200 thresholds from the combined bootstrapping of the two methods, 40 percentiles are extracted. All the subsequent steps of the pipeline are applied to all 40 versions. For each of the final output products (e.g., NEE, as illustrated here), seven percentiles representative of the distribution are included.</figcaption></figure>

a function of slow-changing dynamics (height of canopy, height of measurement, roughness), but a threshold changing every year could introduce false interannual variability. On the other hand, a constant threshold across all the years would not represent changes in the ecosystem structure and EC system setup. For this reason, two approaches were implemented:

*   **Variable USTAR Threshold (VUT):** The thresholds found for each year and the years immediately before and after (if available) have been pooled together, and from their joint population, the final 40 thresholds extracted. With that, the USTAR thresholds vary from year to year; however, they are still influenced by neighboring years. This is identified in FLUXNET2015 variables by the “_VUT” suffix;

*   **Constant USTAR Threshold (CUT):** Across years, all the thresholds found have been pooled together and the final 40 thresholds extracted from this dataset. With that, all years were filtered with the same USTAR threshold. This is identified in FLUXNET2015 variables by the “_CUT” suffix.

If the dataset includes up to two years of data, the two methods give the same result, and only the _VUT is generated.

For both the VUT and CUT approaches, 40 NEE datasets have been created, filtering the original NEE time series using 40 different USTAR values estimated as explained above. The values of the thresholds are reported in the AUXNEE product file. These 40 NEE versions have been used as the basis for all the derived variables provided. An example of the variability of the two methods (CP and MP) is shown in Fig. 4, contrasting the distribution of the bootstrapped results for each method, showing comparable values for some years and divergent values for other years (of the same site). This highlights the importance of applying both methods in this ensemble-like way.

**Filtering NEE based on USTAR thresholds.** The USTAR thresholds are applied to daytime and nighttime data, removing NEE values collected when USTAR is below the threshold and removing also the first half-hour with high turbulence after a period of low turbulence to avoid false emission pulses due to CO₂ accumulated under the canopy and not detected by the storage system (in particular, when a profile is not available at the site). The USTAR filtering is not applied to H and LE, because it has not been proved that when there are CO₂ advective fluxes, these also impact energy fluxes, specifically due to the fact that when advection is in general large (night-time), energy fluxes are small. Figure 5 shows the range of thresholds found (interquartile ranges) across sites in FLUXNET2015. While some sites had low thresholds and low variability in the USTAR thresholds, others show large ranges of values in some more extreme cases (indicating difficulties in estimating the “real” threshold).
<figure><img src="image_6.png" /></figure>

**Fig. 4** Example of the distribution of USTAR thresholds calculated for each year using the MP³⁰ method in blue and CP⁴¹ method in green for the US-UMB site (dark green where they overlap). All these thresholds were pulled together to extract the CUT final 40 thresholds, while for the VUT thresholds, each year was pulled with the two immediately before and after (e.g., 2005 + 2006 + 2007 to extract the 40 thresholds to be used to filter 2006). Note that the level of agreement between methods and between subsequent years is variable, justifying the approach that propagates this variability into uncertainty in NEE.

<figure><img src="image_7.png" /></figure>

**Fig. 5** Ranked USTAR thresholds based on median threshold and error bars showing 25th to 75th percentiles of the 40 thresholds calculated with the Constant USTAR Threshold (CUT) method – only computed for sites with 3 or more years, so only 177 sites out of the 206 are shown. Colors show different ecosystem classes based on the site's IGBP.

*Gap-filling of NEE.* Existing gaps from instrument or power failures are further increased after QC and USTAR filtering. The time series with gaps need to be filled, especially before aggregated values can be calculated (from daily to annual). Moffat *et al.*²⁶ compared different gap-filling methods for CO₂ concluding that most of the methods currently available perform sufficiently well with respect to the general uncertainty associated with the measurements. The method implemented here is the Marginal Distribution Sampling (MDS) method already described in the meteorological products.

*Selection of reference NEE variables.* After filtering NEE using the 40 USTAR thresholds and gap-filling, 40 complete (gap-free) NEE time series were available for each site. For each half hour, it is possible to use the 40 values to estimate the NEE uncertainty resulting from the USTAR threshold estimation (reported as percentiles
of the NEE distribution, identified by the “_XX” numeric suffix) and the average value (identified in the dataset by the “_MEAN” suffix). Since the average value has a smoothing effect on the time series, an additional reference value of NEE was selected and identified in FLUXNET2015 variables by the “_REF” suffix, in an attempt to identify which of the 40 NEE realizations was the most representative of the ensemble. The “_REF” NEE was selected among the 40 different NEE instances in this way: (1) the Nash–Sutcliffe model efficiency coefficient⁴² was calculated between each NEE instance and the remaining 39; (2) the reference NEE has been selected as the one with the highest model efficiency coefficients sum, i.e., the most similar to the other 39. Note that determining the reference NEE is done independently for variables using VUT and CUT USTAR thresholds, as well as for each temporal aggregation. Therefore, the version selected as REF could be different for different temporal resolutions. For instance, NEE_VUT_REF at half-hourly resolution might have been generated using a different USTAR threshold than NEE_VUT_REF at daily resolution. Information on which threshold values were used for each version and temporal aggregation can be found in the auxiliary products for NEE processing (AUXNEE). In addition to the reference NEE, the NEE instance obtained by filtering the data with the median value of the USTAR thresholds distribution is also included. This NEE is identified in FLUXNET2015 variables by the “_USTAR50” suffix (for both CP and MP methods, and both VUT and CUT approaches) and is stable across temporal aggregation resolutions. Individual percentiles of the USTAR thresholds distribution are reported in the AUXNEE file (40 instances for CUT and 40 per year for VUT).

Random uncertainty for NEE. In addition to the uncertainty estimates based on multiple thresholds for USTAR filtering, the random uncertainty for NEE is also estimated based on the method used by Hollinger & Richardson³⁷. Variables expressing random uncertainty are identified by the suffix _RANDUNC. One of two methods are used to estimate random uncertainty, applied hierarchically:

* NEE-RANDUNC Method 1 (direct standard deviation method): For a sliding window of ±7 days and ±1 hour of the time-of-day of the current timestamp, the random uncertainty is calculated as the standard deviation of the measured fluxes. The similarity in the meteorological conditions evaluated as in the MDS gap-filling method³¹ and a minimum of five measured values must be present; otherwise, method 2 is used.

* NEE-RANDUNC Method 2 (median standard deviation method): For a sliding window of ±5 days and ±1 hour of the time-of-day of the current timestamp, random uncertainty is calculated as the median of the random uncertainty (calculated with NEE-RANDUNC Method 1) of similar fluxes, i.e., within the range of ±20% and not less than 2 µmolCO₂ m⁻² s⁻¹. (Note on temporal aggregations: differing sliding windows are: WW: ±2 weeks, MM: ±1 month, and YY: ±2 years.)

The joint uncertainty for NEE is computed from the combination of the uncertainty from multiple USTAR thresholds and random uncertainty. These variables identified by the _JOINTUNC suffix and are computed for NEE filtered using the VUT method as in Eq. (5), and similarly for NEE filtered with the CUT method.

$$NEE\_VUT\_REF\_JOINTUNC = \sqrt{NEE\_VUT\_REF\_RANDUNC^2 + \left( \frac{NEE\_VUT\_84 - NEE\_VUT\_16}{2} \right)^2} \quad (5)$$

The 16th and 84th percentiles are used because they are equivalent to ±1 Standard Deviation in case of a normal distribution. (Note on temporal aggregations: joint uncertainties for NEE are recomputed at all temporal resolutions.)

CO₂ flux partitioning in GPP and RECO. Partitioning CO₂ fluxes from NEE into estimates of its two main components, Gross Primary Production (GPP) and Ecosystem Respiration (RECO), was done by parameterizations of models using measured data. All sites were partitioned with the nighttime fluxes method³¹ (_NT suffixes) and the daytime fluxes method⁴³ (_DT suffixes), while a third method, sundown reference respiration⁴⁴ (_SR suffixes), was applied to all sites meeting the method's requirements (e.g., high quality storage measurement).

The nighttime method uses nighttime data to parameterize a respiration-temperature model that is then applied to the whole dataset to estimate RECO. GPP is then calculated as the difference between RECO and NEE. The parameterization uses short windows of time (14 days) to account for the dynamic of other important respiration drivers such as water, substrate availability, and phenology (see Reichstein et al.³¹ for details on the implementation and ONEFlux²² for the code).

The daytime method uses daytime and nighttime data to parameterize a model with one component based on a light-response curve and vapor pressure deficit for GPP, and a second component using a respiration-temperature relationship similar to the nighttime method. In this case, NEE becomes a function of both GPP and RECO, both of which are estimated by the model. Similarly to the nighttime method, the parameterization is done for short windows (8 days) to take into consideration other slower-changing factors (see Lasslop et al.⁴³ for details on the implementation and ONEFlux²² for the code).

For forest sites where a CO₂ concentration profile for storage fluxes was available, an additional RECO estimate was calculated using the method from van Gorsel et al.⁴⁴, with variables identified by the _SR suffix. In this method, the parameterization of a respiration-temperature model is based solely on data acquired just after sundown, aiming at excluding the measurements potentially affected by advection and also assuming that in the first hour of the evening the advective transport is not yet established.

The sundown partitioning method requires that the NEE is not filtered for low turbulence conditions (USTAR), and for this reason it was applied only to the original time series. The nighttime and daytime methods instead require NEE filtered for low turbulence conditions. For this reason they were applied to all the 40 NEE
versions resulting from the 40 USTAR thresholds, obtaining 40 versions of GPP and RECO for each of the two partitioning methods, propagating the uncertainty from NEE to GPP and RECO. This has been done for both the CUT and VUT filtering methods.

Similarly to NEE, the 40 GPP and RECO estimates (for each method and for CUT and VUT) have been used to calculate the percentiles of their distribution for each timestep (describing their uncertainty due to the NEE uncertainty). The average value (_MEAN) and the reference value use the same model efficiency approach used for NEE for each temporal aggregation. Similarly, NEE filtered with the median USTAR value (_USTAR50) has been partitioned into GPP and RECO. Information on the threshold values used for all versions of GPP and RECO (_NT and _DT, _VUT and _CUT, HH to YY resolutions) are in the auxiliary files for NEE processing (AUXNEE). Variables for reference GPP and RECO are also identified by a _REF suffix. The two methods for the partitioning (three for the cases in which the sundown method is applied) are not merged in any way, because their difference is informative with the respect to the uncertainty of the methods, as in the case of a model comparison exercise.

**Implementation approach.** To increase the traceability of changes between versions of datasets and reduce uncertainty stemming from choices made at implementation time, we favored using original code implementations or thoroughly validated re-implementations of original codes. Thus, our code organization strings together loosely coupled components which implement each step, with clear-cut interfaces between steps. This modular approach eases the maintenance and change efforts for any individual step, but adds complexity to evaluating changes for the entire pipeline. Different programming languages (Python, C, MATLAB and IDL, plus PV-WAVE for FLUXNET2015) were used to implement the different steps, all connected using a controller code that makes appropriate calls in the correct order. The ONEFlux<sup>22</sup> code collection replaced the PV-WAVE code with a re-implementation in Python, and also collates most of these steps into a cohesive pipeline (see also the Code Availability section). The IDL code, which applies the sundown partitioning method<sup>44</sup>, is not yet currently implemented in ONEFlux, because some additional testing and development are needed to make it robust and more suitable for general application. Implementation details of individual steps are discussed next, with references to the outputs each step identified by an execution sequential number and the step name-e.g., 01_qc visual contains the results of the first processing step, the visual check step. Each of these steps correspond to a code module. Supplementary Fig. SM2 shows the steps and their inter-dependencies.

*Steps implemented in python.* The main controller code for ONEFlux is implemented in Python. Besides being the glue code that executes each step, pre- and post-checks are also executed before and after each step. These checks guarantee that the input data meet the minimum requirements to run the step, that the minimum expected outputs were generated by the execution of the step, and that any errors or exception conditions were handled correctly. Information about execution is recorded in a log for the entire pipeline, along with logs for individual steps. Besides the controller code, two of the three flux partitioning steps were re-implemented in Python (the nighttime and daytime methods, 10_nee_partition_nt and 11_nee_partitionDt), together with other specific steps such as data preparation for the uncertainty estimates (12_ure_input), and the creation and checking of final products (99_fluxnet2015). The original flux partitioning implementation in PV-WAVE was used for the LaThuile2007 and FLUXNET2015 datasets. Also, the tool for the downscaling of the ERA-I meteorological data is implemented in Python and runs on a server connected to the ERA data.

*Steps implemented in C.* Several steps are implemented in the C programming language, allowing better control over execution performance of these steps. These steps include:

automated QA/QC flagging (02_qc_auto), USTAR threshold estimation using the MP method (04_ustarPop), the filtering and gap-filling of meteorological data, including the merging with the ERA-I downscaled data (07_meteo_proc), the filtering and gap-filling of CO₂ fluxes (08_nee_proc), the filtering, gap-filling, and energy corrections of energy fluxes (09_energy_proc), and the computation of uncertainty products (12_ure). The source codes and the compiled executables are provided for steps implemented in C, as well as build procedures in make/Makefile format.

*Steps implemented in MATLAB.* The estimation of USTAR thresholds using the CP method (05_ustar cp) is the only step implemented in MATLAB. It is distributed both as source code and compiled code to be used with the MATLAB Runtime Environment, such that it does not require a license purchase.

## Data Records

The FLUXNET2015 portion presented in this paper contains 1496 site-years of data from 206 sites<sup>45-250</sup>, characterizing ecosystem-level carbon and energy fluxes in diverse ecosystems across the globe (Fig. 1, Supplementary Fig. SM1<sup>251,252</sup>), spanning from the early 1990s to 2014, with 69 sites having decade-long records. The dataset covers the distribution of ecosystem fluxes as reported in the recent meta-analyses<sup>253,254</sup> (Fig. 6).

The dataset is distributed in files separated by sites, by temporal aggregation resolutions (e.g., hourly, weekly), and by data products (e.g., FULLSET with all the variables and SUBSET designed for less experienced users). All data files for a site are available for download as a single ZIP file archive with site-specific DOI. The file-naming conventions details these options for each file (Table 1). Site metadata are also available as a single file containing metadata for all sites, detailed later in this section and Supplementary Table SM8. Note that DOIs are assigned at the site level, one DOI per site for all of that site's products. A DOI was not assigned to the whole FLUXNET2015 dataset, since this would make citation and assigning credit imprecise and hard to track.

The FLUXNET2015 dataset provides data at five temporal resolutions. Site teams contribute either half-hourly (HH) or hourly (HR) datasets, depending on the integration/aggregation time decided by the site managers and
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Distribution of the yearly (a) net ecosystem exchange (NEE), (b) gross primary production (GPP), and (c) ecosystem respiration (RECO) in FLUXNET2015. Only data with QC flag (NEE_VUT_REF_QC) higher than 0.5 are shown here. The values are reference NEE, GPP, and RECO based on the Variable USTAR Threshold (VUT) and selected reference for model efficiency (REF). GPP and RECO are based on the nighttime partitioning (NT) method. The grey histogram (bin width 100 gC m⁻² y⁻¹) shows the flux distribution in 1224 of the available site-years; negative GPP and RECO values are kept to preserve distributions, see Data processing methods section for details. Black lines show the distribution curves based on published data²⁵³,²⁵⁴. The boxplots show the flux distribution (i.e., 25th, 50th, and 75th percentiles) for vegetation types defined and color-coded according to IGBP (International Geosphere-Biosphere Programme) definitions. Circles represent data points beyond the 1.5-times interquartile range (25th to 75th percentile) plus the 75th percentile or minus 25th percentile (whisker). Numbers in parentheses indicate the number of site-years used in each IGBP group. The NO-Blv site from the snow/ice IGBP group is not shown in the boxplots.</figcaption></figure>

function of the characteristics of the turbulence. References to half-hourly in this paper also apply to hourly data, unless explicitly stated otherwise. Half-hourly data are the basis of all the processing done for this dataset and are the finest grained temporal resolution provided. Coarser aggregations are generated uniformly from half-hourly data within the data processing pipeline. The other standard temporal aggregations are: daily (DD), weekly (WW), monthly (MM), and yearly (YY).

The complete output from the data processing pipeline includes over 200 variables—among which are measured and derived data, quality flags, uncertainty quantification variables, and results from intermediate data processing steps. The variable names follow the naming conventions of `<BASENAME>_<QUALIFIER>`, where `BASENAME` describes the physical quantities (e.g., TA, NEE, Table 2) and `QUALIFIER` describes the information of processing methods (e.g., VUT, CUT), uncertainties (e.g., RANDUNC), and quality flags (e.g., QC) (see Supplementary Table SM1).

To serve the users with an easier-to-use data product, we created two variants with different selections of variables for data distribution: the FULLSET with all the results and variables; and the SUBSET, designed to help non-expert users, with a reduced set of variables that should fit most needs.

* **FULLSET**: variables generated by the processing such as uncertainty quantification variables, all variants of the data products, all quality information flags, and many variables generated by intermediate processing steps to allow in-depth understanding of individual processing steps and their effect in the final data products. A summary of the main variable basenames is in Table 2, while a full list of variables is provided in Supplementary Table SM1. Key features of the FULLSET version are:
File Name Conventions

FLX_[SITE_ID]_FLUXNET2015_[DATA_PRODUCT]_[RESOLUTION]_[FIRST YEAR]_[LAST YEAR]_[SITE_VERSION]_[CODE_VERSION]_[EXT]

<table><thead><tr><th>Field</th><th>Definition</th><th>Possible options</th></tr></thead><tbody><tr><td>SITE_ID</td><td>FLUXNET site ID in the format CC-SSS (CC is two-letter country code, SSS is three-character site-level identifier)</td><td></td></tr><tr><td>DATA_PRODUCT</td><td>Grouping of variables from release included in file.</td><td><ul><li>SUBSET: Core set of variables with quality and uncertainty information needed for general uses of the data</li><li>FULLSET: All variables, including all quality and uncertainty information, and key variables from intermediate processing steps</li><li>AUXMETEO: Auxiliary variables related to meteorological downscaling</li><li>AUXNEE: Auxiliary variables related to NEE, RECO, and GPP processing</li><li>ERAI: Full record (1989–2014) of ERA-Interim downscaled meteorological variables for the site</li></ul></td></tr><tr><td>RESOLUTION</td><td>Temporal resolution of data product</td><td><ul><li>HH: Half-Hourly time steps</li><li>HR: Hourly time steps</li><li>DD: Daily time steps</li><li>WW: Weekly time steps</li><li>MM: Monthly time steps</li><li>YY: Yearly time steps</li></ul></td></tr><tr><td>FIRST YEAR LAST YEAR</td><td>First and last years of eddy covariance flux data</td><td></td></tr><tr><td>SITE_VERSION CODE_VERSION</td><td>Version string in integer. SITE_VERSION indicates the version of the original dataset for the site used; CODE_VERSION indicates the version of the code of the data processing pipeline used to process the dataset for the site</td><td></td></tr><tr><td>EXT</td><td>File extension</td><td><ul><li>csv: Comma-separated values in a text file (ASCII)</li><li>zip: Archive file with all temporal resolutions for the same site and data product</li></ul></td></tr></tbody></table>

Examples of file names and structures:

FLX_US-Ha1 FLUXNET2015_FULLSET_HH_1992-2012_1-3.zip

- FLX_US-Ha1 FLUXNET2015_FULLSET_HH_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015_FULLSET_DD_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015_FULLSET_WW_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015_FULLSET_MM_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015_FULLSET_YY_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015 ERAI_HH_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015 ERAI_DD_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015 ERAI_WW_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015 ERAI_MM_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015 ERAI_YY_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015 AUXMETEO_1992-2012_1-3.csv

- FLX_US-Ha1 FLUXNET2015 AUXNEE_1992-2012_1-3.csv

<caption>**Table 1.** The template of file naming conventions, including the field, field definition, and the possible options. Examples of file names from a zipped file of a single site are provided.</caption>

* Meteorological variables filled with multiple gap-filling methods (e.g., MDS, ERA) are provided separately.

* NEE versions filtered with two different methods of extracting the USTAR thresholds (i.e., CUT, VUT) are provided. Multiple percentiles and reference NEE are also provided.

* GPP and RECO partitioned from NEE filtered with VUT and CUT methods, using both daytime and nighttime partitioning methods (i.e., NT, DT). Multiple percentiles and reference GPP and RECO are provided.

* LE and H gap-filled, adjusted and non-adjusted for energy balance closure, are both provided.

* Random, methodological, and joint uncertainties for NEE, GPP, RECO, LE, and H are provided.

* **SUBSET**: Includes a subset of the data product. The selection of the variables for this data product was done based on the expected usage for most users and to help less experienced users. Although the number of variables used is reduced, they are still accompanied by a set of quality flags and uncertainty quantification variables essential to correctly interpret the data. Key features of the SUBSET version are:

* Only the consolidated gap-filled meteorological variables are provided.

* Only the REF version of NEE filtered with the VUT method is provided. Selected percentiles and reference NEE are also provided.

* GPP and RECO (only REF versions) partitioned from NEE filtered with only the VUT method, using both daytime and nighttime partitioning methods, are provided. Selected percentiles and reference GPP and RECO are also provided.
<table><caption>**Table 2.** List of the variable basenames, descriptions, available resolutions and units. Separate units are listed if different units are used in different temporal aggregation resolutions. n/a indicates a variable is not provided at the specified resolution.</caption><thead><tr><th rowspan="2">Basename</th><th rowspan="2">Description</th><th colspan="4">Units by Resolution</th></tr><tr><th>HH/HR</th><th>DD</th><th>WW</th><th>MM</th><th>YY</th></tr></thead><tbody><tr><td>TA</td><td>Air temperature</td><td colspan="4">deg C</td></tr><tr><td>SW_IN_POT</td><td>Shortwave radiation, incoming, potential (top of atmosphere)</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>SW_IN</td><td>Shortwave radiation</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>LW_IN</td><td>Longwave radiation, incoming</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>VPD</td><td>Vapor Pressure saturation Deficit</td><td colspan="4">hPa</td></tr><tr><td>PA</td><td>Atmospheric pressure</td><td colspan="4">kPa</td></tr><tr><td>P</td><td>Precipitation</td><td>mm</td><td>mm d<sup>-1</sup></td><td>mm y<sup>-1</sup></td><td></td></tr><tr><td>WS</td><td>Wind speed</td><td colspan="4">m s<sup>-1</sup></td></tr><tr><td>WD</td><td>Wind direction</td><td>Decimal degrees</td><td colspan="3">n/a</td></tr><tr><td>RH</td><td>Relative humidity</td><td>%</td><td colspan="3">n/a</td></tr><tr><td>USTAR</td><td>Friction velocity</td><td colspan="4">m s<sup>-1</sup></td></tr><tr><td>NETRAD</td><td>Net radiation</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>PPFD_IN</td><td>Photosynthetic photon flux density, incoming</td><td colspan="4">µmolPhoton m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>PPFD_DIF</td><td>Photosynthetic photon flux density, diffuse incoming</td><td colspan="4">µmolPhoton m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>PPFD_OUT</td><td>Photosynthetic photon flux density, outgoing</td><td colspan="4">µmolPhoton m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>SW_DIF</td><td>Shortwave radiation, diffuse incoming</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>SW_OUT</td><td>Shortwave radiation, outgoing</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>LW_OUT</td><td>Longwave radiation, outgoing</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>CO<sub>2</sub></td><td>CO<sub>2</sub> mole fraction</td><td colspan="4">µmolCO<sub>2</sub> mol<sup>-1</sup></td></tr><tr><td>TS</td><td>Soil temperature</td><td colspan="4">deg C</td></tr><tr><td>SWC</td><td>Soil water content</td><td colspan="4">%</td></tr><tr><td>G</td><td>Soil heat flux</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>LE</td><td>Latent heat flux</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>H</td><td>Sensible heat flux</td><td colspan="4">W m<sup>-2</sup></td></tr><tr><td>NEE</td><td>Net Ecosystem Exchange</td><td>µmolCO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup></td><td>gC m<sup>-2</sup> d<sup>-1</sup></td><td>gC m<sup>-2</sup> y<sup>-1</sup></td><td></td></tr><tr><td>RECO</td><td>Ecosystem Respiration</td><td>µmolCO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup></td><td>gC m<sup>-2</sup> d<sup>-1</sup></td><td>gC m<sup>-2</sup> y<sup>-1</sup></td><td></td></tr><tr><td>GPP</td><td>Gross Primary Production</td><td>µmolCO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup></td><td>gC m<sup>-2</sup> d<sup>-1</sup></td><td>gC m<sup>-2</sup> y<sup>-1</sup></td><td></td></tr></tbody></table>

* LE and H gap-filled, adjusted and non-adjusted for energy balance closure, are both provided.

* Random and methodological uncertainties for NEE, GPP, RECO, LE, and H are provided.

The variable proposed in the SUBSET product is NEE_VUT_REF since it maintains the temporal variability (as opposed to the MEAN NEE), it is representative of the ensemble, and the VUT method is sensitive to possible changes of the canopy (density and height) and site setup, which can have an impact on the turbulence and consequently on the USTAR threshold. The RECO and GPP products in SUBSET are calculated from the corresponding NEE variables filtered with the VUT method, generating RECO_NT_VUT_REF and RECO_DT_VUT_REF for RECO, and GPP_NT_VUT_REF and GPP_DT_VUT_REF for GPP. It is important to use both daytime (DT) and nighttime (NT) variables, and consider their difference as uncertainty.

Auxiliary data products provide extra information on specific parameters of the data processing pipeline. The groups of products are:

* **AUXMETEO:** Auxiliary data product containing information about the downscaling of meteorological variables using the ERA-Interim reanalysis data product (TA, PA, VPD, WS, P, SW_IN, and LW_IN). Variables in these files relate to the linear regression and error/correlation estimates for each data variable used in the downscaling.
* Parameters:

* ERA_SLOPE: the slope of linear regression

* ERA_INTERCEPT: intercept point of linear regression

* ERA_RMSE: root mean square error between site data and downscaled data

* ERA_CORRELATION: correlation coefficient of linear fit

* AUXNEE: Auxiliary data product with variables resulting from the processing of NEE (mainly related to USTAR filtering) and generation of RECO and GPP. Variables in this product include success/failure of execution of USTAR filtering methods, USTAR thresholds applied to different versions of variables, and percentile/threshold pairs with best model efficiency results.

* Variables:

* USTAR MP METHOD: Moving Point Test USTAR threshold method run

* USTAR CP METHOD: Change Point Detection USTAR threshold method run

* NEE_USTAR50_[UT]: 50th percentile of USTAR thresholds obtained from bootstrapping and sed to generate NEE_USTAR50_[UT] (with UT either CUT or VUT)

* NEE_[UT]_REF: USTAR threshold used to calculate the reference NEE, using model efficiency approach (with UT either CUT or VUT)

* [PROD]_[ALG]_[UT]_REF: USTAR threshold used to filter the NEE that was used to produce the reference product PROD (RECO or GPP), selected using model efficiency approach, using algorithm ALG (NT, DT) (with UT either CUT or VUT)

* Parameters:

* SUCCESS_RUN: 1 if a run of a method (USTAR MP METHOD or USTAR CP METHOD) was successful, 0 otherwise

* USTAR_PERCENTILE: percentile of USTAR thresholds from bootstrapping at USTAR filtering step

* USTAR_THRESHOLD: USTAR threshold value corresponding to USTAR_PERCENTILE

* [RR]_USTAR_PERCENTILE: percentile of USTAR thresholds from bootstrapping at USTAR filtering step at resolution RR (HH, DD, WW, MM, YY)

* [RR]_USTAR_THRESHOLD: USTAR threshold value corresponding to USTAR_PERCENTILE at resolution RR (HH, DD, WW, MM, YY)

* ERAI: Auxiliary data product containing full record (1989–2014) of downscaled meteorological variables using the ERA-Interim reanalysis data product, including TA, PA, VPD, WS, P, SW_IN, and LW_IN.

The FLUXNET2015 metadata are included in a single file (FLX_AA-Flx_BIF_[RESOLUTION]_[YYYYMMDD].xlxs) for all sites for each data product resolution (see Table 1 for resolution options). The metadata follow the Biological, Ancillary, Disturbance, and Metadata (BADM²⁵⁵,²⁵⁶) standards and are provided in the BADM Interchange Format²⁹ (BIF). Table 3 illustrates the type of metadata included with selected metadata variables (See full lists and descriptions of the metadata in Supplementary Tables SM2–SM7). Height and instrument models for the flux variables, as well as soil temperature and moisture depths, are reported in the Variable Information metadata.

## Technical Validation

Eddy covariance measurements offer a direct method to estimate trace gas or energy exchanges between surface and atmosphere at an ecosystem scale (approximately up to 1 km around the measurement point). This makes eddy covariance difficult to compare with other methods. Nonetheless, eddy covariance data have been extensively used in numerous scientific papers and studies that indirectly validate their reliability and usefulness. Hundreds of articles have been published based on eddy covariance measurements; examples of multi-site studies using FLUXNET2015 data include Jung *et al.*¹⁵, Tramontana *et al.*²⁵⁷, and Keenan *et al.*²⁵⁸.

Eddy covariance data were evaluated with respect to other methods such as inventory and chambers by Campioli *et al.*²⁵⁹, who showed that “EC [eddy covariance] biases are not apparent across sites, suggesting the effectiveness of standard post-processing procedures. Our results increase confidence in EC …”. The approach of Campioli *et al.*²⁵⁹ requires sites that have several additional (and rare) pieces of information; therefore, it is not generally applicable, particularly not across the sites used in this study. However, the eddy covariance site teams co-authoring this paper have compared and technically validated their measurements with respect to knowledge of their site. Unavoidably, measurement and processing uncertainties exist, and can be large for certain sites and ecosystem conditions. However, in general, flux values provided in this dataset are consistent with expectations, and eddy covariance remains one of the more reliable techniques for assessing land-air exchanges at ecosystem scales.

## Usage Notes

Detailed documentation on how to use and interpret FLUXNET2015 is available online at https://fluxnet.flux-data.org/data/fluxnet2015-dataset/. Here, we present some of the main points to guide the usage of the data.

**Risks in the application of standard procedures.** When standardized procedures are applied across different sites, the possible differences owing to data treatment are avoided or minimized; this is one of the main
<table><caption>**Table 3.** Metadata types and selected variables. See Supplementary Tables SM2–SM7 for a full list of metadata with descriptions. Variables collected from or generated for all sites are in bold.</caption><thead><tr><th>Metadata Type</th><th>Selected Metadata Variables</th></tr></thead><tbody><tr><td>Site General Info<br>(25 variables)</td><td><table><tr><td>SITE_ID:</td><td>Unique site identifier (CC-sss, where CC is the country code)</td></tr><tr><td>SITE_NAME:</td><td>Site name</td></tr><tr><td>SITE_DESC:</td><td>Site description</td></tr><tr><td>LOCATION_LAT:</td><td>Latitude of site</td></tr><tr><td>LOCATION_LONG:</td><td>Longitude of site</td></tr><tr><td>FLUX_MEASUREMENTS_VARIABLE:</td><td>Flux variables measured at the site</td></tr><tr><td>IGBP:</td><td>Vegetation type based on International Geosphere-Biosphere Programme classification</td></tr><tr><td>UTC_OFFSET:</td><td>Offset from UTC of site data</td></tr><tr><td>TEAM_membership_NAME:</td><td>Team member name</td></tr><tr><td>MAT, MAP:</td><td>Mean annual temperature and precipitation</td></tr><tr><td>TOWER_TYPE, TOWER_POWER:</td><td>Type of tower and power type</td></tr></table></td></tr><tr><td>DOI<br>(12 variables)</td><td><table><tr><td>DOI:</td><td>Digital Object Identifier (DOI) for the flux-met data product</td></tr><tr><td>DOI_CONTRIBUTOR_NAME:</td><td>Name of contributor to the development of data (and associated info)</td></tr><tr><td>DOI_ORGANIZATION:</td><td>Organization contributing to the data</td></tr></table></td></tr><tr><td>Reference publications<br>(4 variables)</td><td><table><tr><td>REFERENCE_PAPER:</td><td>Reference for understanding the site</td></tr><tr><td>REFERENCE_DOI:</td><td>DOI of the reference</td></tr><tr><td>REFERENCE_USAGE:</td><td>Suggested usage of the reference</td></tr></table></td></tr><tr><td>Canopy Height<br>(2 variables)</td><td><table><tr><td>HEIGHTC:</td><td>Canopy height. In a forest ecosystem, canopy height is the distribution of overstory trees that see light at the top of the canopy.</td></tr><tr><td>Note:</td><td>The reported value is representative of the mean of such a distribution.</td></tr><tr><td>HEIGHTC_DATE:</td><td>Date of canopy height observation</td></tr></table></td></tr><tr><td>Variable Information<br>(5 variables)<br>Note: Variable Information groups are only reported for variables with data.</td><td><table><tr><td>VAR_INFO_VARNAME:</td><td>Variable name</td></tr><tr><td>VAR_INFO_UNIT:</td><td>Variable unit</td></tr><tr><td>VAR_INFO_DATE:</td><td>Start date for reported variable information</td></tr><tr><td>VAR_INFOHeight:</td><td>Height/depth of observation (meters)</td></tr><tr><td>VAR_INFO MODEL:</td><td>Model(s) used to collect observation.</td></tr></table></td></tr></tbody></table>

goals of FLUXNET2015 and ONEFlux. However, there is also the risk and possibility that the standard methods don't work properly or as expected at specific sites and under certain conditions. This is particularly true for the CO₂ flux partitioning, which as with all models is based on assumptions that could not always be not valid. For this reason, it could be necessary to contact the site PIs that are listed in Supplementary Table SM9.

**Using the QC flags.** There are quality-control flag variables in the dataset to help users filter and interpret variables, especially for gap-filled and process knowledge-based variables. These flags are described in the variable documentation (Supplementary Table SM1). It is highly recommended that one carefully considers the QC flags when using the data.

**Percentile variants for fluxes and reference values.** For most flux variables, there are reference values and percentile versions of the variables to help understand some of the uncertainty in the record. For NEE, RECO, and GPP, the percentiles are generated from the bootstrapping of the USTAR threshold estimation step, i.e., they characterize the variability from a range of values obtained as USTAR thresholds. In addition, three different reference values are provided (“_MEAN”, “_USTAR50” and “_REF”) in order to cover different user needs. In general the “_REF” version should be the most representative, particularly if related to the percentiles. It is, however, important to clearly refer to which NEE version is used in order to ensure reproducibility. For the energy balance corrected H and LE variables, the percentiles indicate the variability due to the uncertainty in the correction factor applied. Similarly to NEE, there are gap-filled and energy balance corrected versions of H and LE variables; therefore, it is also important to clearly refer to which version is used. The SUBSET version of the dataset includes a reduced number of variables, selected for non-expert users. We encourage users to carefully evaluate their requirements and options in the dataset, and if needed to contact regional networks, site teams, or even co-authors of this article for help and recommendations. For more detail, see the Methods section above.

**Temporally aggregated resolutions.** All data products are provided at multiple temporal resolutions where feasible. The finest resolution is either hourly or half-hourly (indicated by the filename tags HR and HH, respectively). These data are then aggregated into daily (DD), 7-day weekly (WW), monthly (MM), and yearly (YY) resolutions, with appropriate aggregations for each variable, such as averaging for TA and summation for P.

**Timestamps.** Timestamps in the data and metadata files use the format YYYYMMDDHHMM, truncated at the appropriate resolution (e.g., YYYYMMDD for a date or YYYYMM for a month). Two formats of time associated with a record are used: (1) single timestamp, where the meaning of the timestamp is unambiguous, and, (2) a pair of timestamps, indicating the start and end of the period the timestamps represent.

**Time zones.** To allow more direct site comparability, all time variables are reported in local standard time (i.e., without daylight saving time). The time zone information with respect to UTC time is reported in the site metadata.

**Numeric resolution.** The floating point numbers are maintained at their original resolution throughout processing steps, using double precision for the majority of cases, and are truncated at up to nine decimal places in the distributed files for numbers between 0.0 and 1.0, and at up to five decimal places for larger numbers.
**Column ordering.** The order of columns is not always the same in different files (e.g., different sites). User data-processing routines should use the variable label (which is always consistent) and not the order of occurrence of that variable in the file. Timestamps are the only exception and will always be the first variable(s)/ column(s) of the data file. This applies to text file data representations (i.e., CSV formatted).

**Missing data.** Missing data values are indicated with -9999, without decimal points, independent of the cause of the missing value.

**Known issues.** A list of known issues and limitations relevant to the dataset is maintained online: http://fluxnet.fluxdata.org/data/fluxnet2015-dataset/known-issues/.

**Releases of the FLUXNET2015 dataset.** The original FLUXNET2015 release was in December 2015, followed by incremental releases in July 2016 and November 2016, and, finally, a release in February 2020 with fixes and additional metadata as described in this paper. More information on the releases can be found in the online change log: http://fluxnet.fluxdata.org/change-log/. A newer release replaces all previous ones, and only the newest release is available for direct download. Access to previous versions can be obtained upon request.

**Support to FLUXNET2015 data users.** Scientists and staff responsible for the creation of the dataset offer support to data users and can be reached at fluxdata-support@fluxdata.org.

**Updates and future versions.** There is strong interest and engagement in order to ensure the availability of new data (new sites and new years), keeping the open policy and the high quality data that we tried to reach with this work. We expect that the processing pipeline and QA/QC procedures will continue evolving, in support of new products. However, the amount of both technical and coordination work, along with difficulty securing long-term international funding, hamper creation of new versions of the dataset. There are ongoing discussions among regional networks and FLUXNET on this coordination, but currently there is no plan for a follow-up version of FLUXNET2015.

## Code availability

The ONEFlux collection of codes used to create data intercomparable with FLUXNET2015 has been packaged to be executed as a complete pipeline and is available in both source-code and executable forms under a 3-clause BSD license on GitHub: https://github.com/AmeriFlux/ONEFlux. The complete environment to run this pipeline requires a GCC compatible C compiler (or capability to run pre-compiled Windows, Linux, and/or Mac executables), a MATLAB Runtime Environment, and a Python interpreter with a few numeric and scientific packages installed. All of these can be obtained at no cost.

Received: 5 March 2020; Accepted: 20 May 2020;

Published: 9 July 2020

## References

1. Aubinet, M., Vesala, T. & Papale, D. (Eds.). *Eddy Covariance: A Practical Guide to Measurement and Data Analysis* (Springer Netherlands, Dordrecht, 2012).

2. Aubinet, M. et al. Estimates of the Annual Net Carbon and Water Exchange of Forests: The EUROFLUX Methodology, in *Advances in Ecological Research* (ed. Fitter, A. H., Raffaelli, D. G.), pp. 113–175 (Elsevier, 1999).

3. Law, B. AmeriFlux Network aids global synthesis. *Eos Trans. AGU* **88**, 286–286 (2007).

4. Novick, K. A. et al. The AmeriFlux network: A coalition of the willing. *Agr. Forest Meteorol.* **249**, 444–456 (2018).

5. Yamamoto, S. et al. Findings through the AsiaFlux network and a view toward the future. *J. Geogr. Sci.* **15**, 142–148 (2005).

6. Yu, G.-R. et al. Overview of ChinaFLUX and evaluation of its eddy covariance measurement. *Agr. Forest Meteorol.* **137**, 125–137 (2006).

7. Beringer, J. et al. An introduction to the Australian and New Zealand flux tower network - OzFlux. *Biogeosciences* **13**, 5895–5916 (2016).

8. Williams, C.A. et al. Climate and vegetation controls on the surface water balance: Synthesis of evapotranspiration measured across a global network of flux towers. *Water Resour. Res.* **48**, W06523 (2012).

9. Vargas, R. et al. Drought Influences the Accuracy of Simulated Ecosystem Fluxes: A Model-Data Meta-analysis for Mediterranean Oak Woodlands. *Ecosystems* **16**, 749–764 (2013).

10. Shao, J. et al. Biotic and climatic controls on interannual variability in carbon fluxes across terrestrial ecosystems. *Agr. Forest Meteorol.* **205**, 11–22 (2015).

11. Balzarolo, M. et al. On the relationship between ecosystem-scale hyperspectral reflectance and CO₂ exchange in European mountain grasslands. *Biogeosciences* **12**, 3089–3108 (2015).

12. Wolf, S. et al. Warm spring reduced carbon cycle impact of the 2012 US summer drought. *Proc Natl Acad Sci USA* **113**, 5880–5885 (2016).

13. Jung, M. et al. Compensatory water effects link yearly global land CO₂ sink changes to temperature. *Nature* **541**, 516–520 (2017).

14. Besnard, S. et al. Quantifying the effect of forest age in annual net forest carbon balance. *Environ. Res. Lett.* **13**, 124018 (2018).

15. Jung, M. et al. The FLUXCOM ensemble of global land-atmosphere energy fluxes. *Sci. Data* **6**, 74 (2019).

16. Baldocchi, D. et al. FLUXNET: A New Tool to Study the Temporal and Spatial Variability of Ecosystem-Scale Carbon Dioxide, Water Vapor, and Energy Flux Densities. *Bull. Amer. Meteor. Soc.* **82**, 2415–2434 (2001).

17. Chu, H., Baldocchi, D. D., John, R., Wolf, S. & Reichstein, M. Fluxes all of the time? A primer on the temporal representativeness of FLUXNET. *J. Geophys. Res. Biogeosci.* **122**, 289–307 (2017).

18. Pastorello, G. et al. A New Data Set to Keep a Sharper Eye on Land-Air Exchanges. *Eos* **98**, 2017 Apr 17 (2017).

19. Falge, E. et al. FLUXNET Marconi Conference Gap-Filled Flux and Meteorology Data, 1992–2000. ORNL DAAC, https://doi.org/10.3334/ornlaac/811 (2000).

20. FLUXNET-Fluxdata Portal. LaThuile Dataset, https://fluxnet.fluxdata.org/data/la-thuile-dataset/ (2007).

21. FLUXNET-Fluxdata Portal. FLUXNET2015 Dataset, https://fluxnet.fluxdata.org/data/fluxnet2015-dataset/ (2020).
22. Pastorello, G. et al. ONEFlux: Open Network-Enabled Flux processing pipeline. *AmeriFlux Management Project, European Ecosystem Fluxes Database, ICOS Ecosystem Thematic Center*, https://github.com/AmeriFlux/ONEFlux/ (2019).

23. FLUXNET-Fluxdata Portal. Fluxdata - The Data Portal serving the FLUXNET community, https://fluxnet.fluxdata.org/ (2020).

24. Aubinet, M. et al. Nighttime Flux Correction, in Eddy Covariance (eds. Aubinet, M., Vesala, T. & Papale, D.), pp. 133–157, (Springer Netherlands, Dordrecht, 2012).

25. Pastorello, G. et al. Observational Data Patterns for Time Series Data Quality Assessment, in Proc. IEEE 10th International Conference on e-Science, Sao Paulo, Brazil, pp. 271–278 (IEEE, 2014).

26. Moffat, A. M. et al. Comprehensive comparison of gap-filling techniques for eddy covariance net carbon fluxes. *Agr. Forest Meteorol.* **147**, 209–232 (2007).

27. Desai, A. R. et al. Cross-site evaluation of eddy covariance GPP and RE decomposition techniques. *Agr. Forest Meteorol.* **148**, 821–838 (2008).

28. AmeriFlux. Data Variables and Format, https://ameriflux.lbl.gov/data/aboutdata/data-variables/ (2020).

29. AmeriFlux. BADM Data Templates, https://ameriflux.lbl.gov/data/badm-data-templates/ (2020).

30. Papale, D. et al. Towards a standardized processing of Net Ecosystem Exchange measured with eddy covariance technique: algorithms and uncertainty estimation. *Biogeosciences* **3**, 571–583 (2006).

31. Reichstein, M. et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm. *Global Change Biol.* **11**, 1424–1439 (2005).

32. Vuichard, N. & Papale, D. Filling the gaps in meteorological continuous data measured at FLUXNET sites with ERA-Interim reanalysis. *Earth Syst. Sci. Data* **7**, 157–171 (2015).

33. Dee, D. P. et al. The ERA-Interim reanalysis: configuration and performance of the data assimilation system. *Q.J.R. Meteorol. Soc.* **137**, 553–597 (2011).

34. European Centre for Medium-Range Weather Forecasts. ERA-Interim is a global atmospheric reanalysis dataset, https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era-interim (2020).

35. Stoy, P. C. et al. A data-driven analysis of energy balance closure across FLUXNET research sites: The role of landscape scale heterogeneity. *Agr. Forest Meteorol.* **171–172**, 137–152 (2013).

36. Foken, T. The Energy Balance Closure Problem: An Overview. *Ecol. Appl.* **18**, 1351–1367 (2008).

37. Hollinger, D. Y. & Richardson, A. D. Uncertainty in eddy covariance measurements and its application to physiological models. *Tree Physiol.* **25**, 873–885 (2005).

38. Finnigan, J. The storage term in eddy flux calculations. *Agr. Forest Meteorol.* **136**, 108–113 (2006).

39. Nicolini, G. et al. Impact of CO₂ storage flux sampling uncertainty on net ecosystem exchange measured by eddy covariance. *Agr. Forest Meteorol.* **248**, 228–239 (2018).

40. Aubinet, M. et al. Direct advection measurements do not help to solve the night-time CO₂ closure problem: Evidence from three different forests. *Agr. Forest Meteorol.* **150**, 655–664 (2010).

41. Barr, A. G. et al. Use of change-point detection for friction–velocity threshold evaluation in eddy-covariance studies. *Agr. Forest Meteorol.* **171–172**, 31–45 (2013).

42. Nash, J. E. & Sutcliffe, J. V. River flow forecasting through conceptual models part I — A discussion of principles. *J. Hydrol.* **10**, 282–290 (1970).

43. Lasslop, G. et al. Separation of net ecosystem exchange into assimilation and respiration using a light response curve approach: critical issues and global evaluation. *Glob. Change Biol.* **16**, 187–208 (2010).

44. van Gorsel, E. et al. Estimating nocturnal ecosystem respiration from the vertical turbulent flux and change in storage of CO₂. *Agr. Forest Meteorol.* **149**, 1919–1930 (2009).

45. Garcia, A. et al. FLUXNET2015 AR-SLu San Luis. FLUXNET; Instituto Nacional de Tecnología Agropecuaria (INTA) https://doi.org/10.18140/flx/1440191 (2016).

46. Posse, G., Lewczuk, N., Richter, K. & Cristiano, P. FLUXNET2015 AR-Vir Virasoro. FLUXNET; Instituto Nacional de Tecnología Agropecuaria https://doi.org/10.18140/flx/1440192 (2016).

47. Wohlfahrt, G., Hammerle, A. & Hörtnagl, L. FLUXNET2015 AT-Neu Neustift. FLUXNET; University of Innsbruck https://doi.org/10.18140/flx/1440121 (2016).

48. Beringer, J. & Hutley, L. FLUXNET2015 AU-Ade Adelaide River. FLUXNET; Monash University; Charles Darwin University https://doi.org/10.18140/flx/1440193 (2016).

49. Cleverly, J. & Eamus, D. FLUXNET2015 AU-ASM Alice Springs. FLUXNET; University of Technology Sydney https://doi.org/10.18140/flx/1440194 (2016).

50. Meyer, W. et al. FLUXNET2015 AU-Cpr Calperum. FLUXNET; University of Adelaide https://doi.org/10.18140/flx/1440195 (2016).

51. Pendall, E. & Griebel, A. FLUXNET2015 AU-Cum Cumberland Plains. FLUXNET; Western Sydney University https://doi.org/10.18140/flx/1440196 (2016).

52. Beringer, J. & Hutley, L. FLUXNET2015 AU-DaP Daly River Savanna. FLUXNET; Monash University; Charles Darwin University https://doi.org/10.18140/flx/1440123 (2016).

53. Beringer, J. & Hutley, P. FLUXNET2015 AU-DaS Daly River Cleared. FLUXNET; University of Western Australia; Charles Darwin University; Monash University https://doi.org/10.18140/flx/1440122 (2016).

54. Beringer, J. & Hutley, L. FLUXNET2015 AU-Dry Dry River. FLUXNET; Monash University; University of Western Australia; Charles Darwin University https://doi.org/10.18140/flx/1440197 (2016).

55. Schroder, I., Zegelin, S., Palu, T. & Feitz, A. FLUXNET2015 AU-Emr Emerald. FLUXNET; CSIRO; Geoscience Australia https://doi.org/10.18140/flx/1440198 (2016).

56. Beringer, J. & Hutley, L. FLUXNET2015 AU-Fog Fogg Dam. FLUXNET; Monash University; Charles Darwin University https://doi.org/10.18140/flx/1440124 (2016).

57. Macfarlane, C. et al. FLUXNET2015 AU-Gin Gingin. FLUXNET; Edith Cowan University (Centre for Ecosystem Management) https://doi.org/10.18140/flx/1440199 (2016).

58. Macfarlane, C., Prober, S. & Wiehl, G. FLUXNET2015 AU-GWW Great Western Woodlands, Western Australia, Australia. FLUXNET; CSIRO https://doi.org/10.18140/flx/1440200 (2016).

59. Beringer, J. & Hutley, L. FLUXNET2015 AU-How Howard Springs. FLUXNET; Charles Darwin University; University of Western Australia; Monash University https://doi.org/10.18140/flx/1440125 (2016).

60. Ewenz, C., Stevens, R. & Grigson, G. FLUXNET2015 AU-Lox Loxton. FLUXNET; South Australian Research and Development Institute (SARDI) https://doi.org/10.18140/flx/1440247 (2016).

61. Beringer, J. & Hutley, L. FLUXNET2015 AU-RDF Red Dirt Melon Farm, Northern Territory. FLUXNET; Monash University; Charles Darwin University https://doi.org/10.18140/flx/1440201 (2016).

62. Beringer, J. et al. FLUXNET2015 AU-Rig Riggs Creek. FLUXNET; Monash University https://doi.org/10.18140/flx/1440202 (2016).

63. Liddell, M. FLUXNET2015 AU-Rob Robson Creek, Queensland, Australia. FLUXNET; James Cook University https://doi.org/10.18140/flx/1440203 (2016).

64. Beringer, J. & Hutley, L. FLUXNET2015 AU-Stp Sturt Plains. FLUXNET; University of Western Australia; Charles Darwin University; Monash University https://doi.org/10.18140/flx/1440204 (2016).

65. Cleverly, J. & Eamus, D. FLUXNET2015 AU-TTE Ti Tree East. FLUXNET; University of Technology Sydney https://doi.org/10.18140/flx/1440205 (2016).
66. Woodgate, W., Van Gorsel, E. & Leuning, R. FLUXNET2015 AU-Tum Tumbarumba. FLUXNET; CSIRO https://doi.org/10.18140/flx/1440126 (2016).

67. Beringer, J., Hutley, L., McGuire, D. & U, P. FLUXNET2015 AU-Wac Wallaby Creek, FLUXNET; Monash University; University of California Davis; Charles Darwin University; University of Alaska Fairbanks; University of Melbourne, https://doi.org/10.18140/flx/1440127 (2016).

68. Beringer, J. et al. FLUXNET2015 AU-Whr Whroo. FLUXNET; Monash University https://doi.org/10.18140/flx/1440206 (2016).

69. Arndt, S., Hinko-Najera, N. & Griebel, A. FLUXNET2015 AU-Wom Wombat. FLUXNET; University of Melbourne, School of Ecosystem and Forest Sciences https://doi.org/10.18140/flx/1440207 (2016).

70. Beringer, J. & Walker, J. FLUXNET2015 AU-Ync Jaxa. FLUXNET; University of Western Australia; Monash University https://doi.org/10.18140/flx/1440208 (2016).

71. Neirynck, J. et al. FLUXNET2015 BE-Bra Brasschaat. FLUXNET; University of Antwerp https://doi.org/10.18140/flx/1440128 (2016).

72. De Ligne, A. et al. FLUXNET2015 BE-Lon Lonzee. FLUXNET; University of Liege - Gembloux Agro-Bio Tech https://doi.org/10.18140/flx/1440129 (2016).

73. De Ligne, A. et al. FLUXNET2015 BE-Vie Vielsalm. FLUXNET; University of Liege - Gembloux Agro-Bio Tech; University catholic of Louvain-la-Neuve, https://doi.org/10.18140/flx/1440130 (2016).

74. Saleska, S. FLUXNET2015 BR-Sa1 Santarem-Km67-Primary Forest. FLUXNET; University of Arizona https://doi.org/10.18140/flx/1440032 (2016).

75. Goulden, M. FLUXNET2015 BR-Sa3 Santarem-Km83-Logged Forest. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440033 (2016).

76. McCaughey, H. FLUXNET2015 CA-Gro Ontario - Groundhog River, Boreal Mixedwood Forest. FLUXNET; Queen's University https://doi.org/10.18140/flx/1440034 (2016).

77. Amiro, B. FLUXNET2015 CA-Man Manitoba - Northern Old Black Spruce (former BOREAS Northern Study Area). FLUXNET; University of Manitoba https://doi.org/10.18140/flx/1440035 (2016).

78. Goulden, M. FLUXNET2015 CA-NS1 UCI-1850 burn site. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440036 (2016).

79. Goulden, M. FLUXNET2015 CA-NS2 UCI-1930 burn site. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440037 (2016).

80. Goulden, M. FLUXNET2015 CA-NS3 UCI-1964 burn site. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440038 (2016).

81. Goulden, M. FLUXNET2015 CA-NS4 UCI-1964 burn site wet. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440039 (2016).

82. Goulden, M. FLUXNET2015 CA-NS5 UCI-1981 burn site. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440040 (2016).

83. Goulden, M. FLUXNET2015 CA-NS6 UCI-1989 burn site. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440041 (2016).

84. Goulden, M. FLUXNET2015 CA-NS7 UCI-1998 burn site. FLUXNET; University of California - Irvine https://doi.org/10.18140/flx/1440042 (2016).

85. Black, T. FLUXNET2015 CA-Oas Saskatchewan - Western Boreal, Mature Aspen. FLUXNET; The University of British Columbia https://doi.org/10.18140/flx/1440043 (2016).

86. Black, T. FLUXNET2015 CA-Obs Saskatchewan - Western Boreal, Mature Black Spruce. FLUXNET; The University of British Columbia https://doi.org/10.18140/flx/1440044 (2016).

87. Margolis, H. FLUXNET2015 CA-Qfo Quebec - Eastern Boreal, Mature Black Spruce. FLUXNET; Université Laval https://doi.org/10.18140/flx/1440045 (2016).

88. Amiro, B. FLUXNET2015 CA-SF1 Saskatchewan - Western Boreal, forest burned in 1977. FLUXNET; University of Manitoba https://doi.org/10.18140/flx/1440046 (2016).

89. Amiro, B. FLUXNET2015 CA-SF2 Saskatchewan - Western Boreal, forest burned in 1989. FLUXNET; University of Manitoba https://doi.org/10.18140/flx/1440047 (2016).

90. Amiro, B. FLUXNET2015 CA-SF3 Saskatchewan - Western Boreal, forest burned in 1998. FLUXNET; University of Manitoba; Canadian Forest Service https://doi.org/10.18140/flx/1440048 (2016).

91. Arain, M. FLUXNET2015 CA-TP1 Ontario - Turkey Point 2002 Plantation White Pine. FLUXNET; McMaster University https://doi.org/10.18140/flx/1440050 (2016).

92. Arain, M. FLUXNET2015 CA-TP2 Ontario - Turkey Point 1989 Plantation White Pine. FLUXNET; McMaster University https://doi.org/10.18140/flx/1440051 (2016).

93. Arain, M. FLUXNET2015 CA-TP3 Ontario - Turkey Point 1974 Plantation White Pine. FLUXNET; McMaster University https://doi.org/10.18140/flx/1440052 (2016).

94. Arain, M. FLUXNET2015 CA-TP4 Ontario - Turkey Point 1939 Plantation White Pine. FLUXNET; McMaster University https://doi.org/10.18140/flx/1440053 (2016).

95. Arain, M. FLUXNET2015 CA-TPD Ontario - Turkey Point Mature Deciduous. FLUXNET; McMaster University https://doi.org/10.18140/flx/1440112 (2016).

96. Nouvellon, Y. FLUXNET2015 CG-Tch Tchizalamou. FLUXNET; Centre de coopération internationale en recherche agronomique pour le développement https://doi.org/10.18140/flx/1440142 (2016).

97. Hörtnagl, L. et al. FLUXNET2015 CH-Cha Chamau. FLUXNET; ETH Zurich https://doi.org/10.18140/flx/1440131 (2016).

98. Hörtnagl, L. et al. FLUXNET2015 CH-Dav Davos. FLUXNET; ETH Zurich https://doi.org/10.18140/flx/1440132 (2016).

99. Hörtnagl, L. et al. FLUXNET2015 CH-Fru Früebüel. FLUXNET; ETH Zurich https://doi.org/10.18140/flx/1440133 (2016).

100. Hörtnagl, L. et al. FLUXNET2015 CH-Lae Laegern. FLUXNET; ETH Zurich https://doi.org/10.18140/flx/1440134 (2016).

101. Ammann, C. FLUXNET2015 CH-Oe1 Oensingen grassland. FLUXNET; Agroscope Zuerich https://doi.org/10.18140/flx/1440135 (2016).

102. Hörtnagl, L. et al. FLUXNET2015 CH-Oe2 Oensingen crop. FLUXNET; ETH Zurich https://doi.org/10.18140/flx/1440136 (2016).

103. Zhang, J. & Han, S. FLUXNET2015 CN-Cha Changbaishan. FLUXNET; IAE Chinese Academy of Sciences https://doi.org/10.18140/flx/1440137 (2016).

104. Dong, G. FLUXNET2015 CN-Cng Changling. FLUXNET; Shanxi University https://doi.org/10.18140/flx/1440209 (2016).

105. Shi, P., Zhang, X. & He, Y. FLUXNET2015 CN-Dan Dangxiong. FLUXNET; IGSNRR Chinese Academy of Sciences https://doi.org/10.18140/flx/1440138 (2016).

106. Zhou, G. & Yan, J. FLUXNET2015 CN-Din Dinghushan. FLUXNET; SCIB Chinese Academy of Sciences https://doi.org/10.18140/flx/1440139 (2016).

107. Chen, S. FLUXNET2015 CN-Du2 Duolun_grassland (D01). FLUXNET; Institute of Botany, Chinese Academy of Sciences https://doi.org/10.18140/flx/1440140 (2016).

108. Shao, C. FLUXNET2015 CN-Du3 Duolun Degraded Meadow. FLUXNET https://doi.org/10.18140/flx/1440210 (2016).

109. Li, Y. FLUXNET2015 CN-Ha2 Haibei Shrubland. FLUXNET; NWIPB Chinese Academy of Sciences https://doi.org/10.18140/flx/1440211 (2016).
110. Tang, Y., Kato, T. & Du, M. FLUXNET2015 CN-HaM Haibei Alpine Tibet site. FLUXNET; National Institute for Environmental Studies https://doi.org/10.18140/flx/1440190 (2016).

111. Wang, H. & Fu, X. FLUXNET2015 CN-Qia Qianyanzhou. FLUXNET; IGSNRR Chinese Academy of Sciences https://doi.org/10.18140/flx/1440141 (2016).

112. Shao, C. FLUXNET2015 CN-Sw2 Siziwang Grazed (SZWG). FLUXNET https://doi.org/10.18140/flx/1440212 (2016).

113. Sigut, L. et al. FLUXNET2015 CZ-BK1 Bily Kriz forest. FLUXNET; Global Change Research Institute CAS https://doi.org/10.18140/flx/1440143 (2016).

114. Sigut, L. et al. FLUXNET2015 CZ-BK2 Bily Kriz grassland. FLUXNET; Global Change Research Institute CAS https://doi.org/10.18140/flx/1440144 (2016).

115. Dušek, J., Janouš, D. & Pavelka, M. FLUXNET2015 CZ-wet Trebon (CZECHWET). FLUXNET; Global Change Research Institute CAS https://doi.org/10.18140/flx/1440145 (2016).

116. Bernhofer, C. et al. FLUXNET2015 DE-Akm Ankłam. FLUXNET; TU Dresden https://doi.org/10.18140/flx/1440213 (2016).

117. Brümmer, C., Lucas-Moffat, A., Herbst, M. & Kolle, O. FLUXNET2015 DE-Geb Gebesee. FLUXNET; Thünen Institute of Climate-Smart Agriculture, Braunschweig https://doi.org/10.18140/flx/1440146 (2016).

118. Bernhofer, C. et al. FLUXNET2015 DE-Gri Grillenburg. FLUXNET; TU Dresden https://doi.org/10.18140/flx/1440147 (2016).

119. Knohl, A. et al. FLUXNET2015 DE-Hai Hainich. FLUXNET; University of Goettingen, Bioclimatology https://doi.org/10.18140/flx/1440148 (2016).

120. Bernhofer, C. et al. FLUXNET2015 DE-Kli Klingenberg. FLUXNET; TU Dresden https://doi.org/10.18140/flx/1440149 (2016).

121. Lindauer, M. et al. FLUXNET2015 DE-Lkb Lackenberg. FLUXNET; Karlsruhe Institute of Technology, IMK-IFU https://doi.org/10.18140/flx/1440214 (2016).

122. Knohl, A. et al. FLUXNET2015 DE-Lnf Leinefelde. FLUXNET; University of Goettingen, Bioclimatology https://doi.org/10.18140/flx/1440150 (2016).

123. Bernhofer, C. et al. FLUXNET2015 DE-Obe Oberbärenburg. FLUXNET; TU Dresden https://doi.org/10.18140/flx/1440151 (2016).

124. Schmidt, M. & Graf, A. FLUXNET2015 DE-RuR Rollesbroich. FLUXNET; IBG-3 Agrosphäre, Research Centre Jülich GmbH https://doi.org/10.18140/flx/1440215 (2016).

125. Schmidt, M. & Graf, A. FLUXNET2015 DE-RuS Selhausen Juelich. FLUXNET; IBG-3 Agrosphäre, Research Centre Jülich GmbH https://doi.org/10.18140/flx/1440216 (2016).

126. Schneider, K. & Schmidt, M. FLUXNET2015 DE-Seh Selhausen. FLUXNET; University of Cologne https://doi.org/10.18140/flx/1440217 (2016).

127. Klatt, J., Schmid, H., Mauder, M. & Steinbrecher, R. FLUXNET2015 DE-SfN Schechenfilz Nord. FLUXNET; Karlsruhe Institute of Technology, IMK-IFU https://doi.org/10.18140/flx/1440219 (2016).

128. Bernhofer, C. et al. FLUXNET2015 DE-Spw Spreewald. FLUXNET; TU Dresden https://doi.org/10.18140/flx/1440220 (2016).

129. Bernhofer, C. et al. FLUXNET2015 DE-Tha Tharandt. FLUXNET; TU Dresden https://doi.org/10.18140/flx/1440152 (2016).

130. Sachs, T., Wille, C., Larmanou, E. & Franz, D. FLUXNET2015 DE-Zrk Zarnekow. FLUXNET; GFZ German Research Centre for Geosciences https://doi.org/10.18140/flx/1440221 (2016).

131. Pilegaard, K. & Ibrom, A. FLUXNET2015 DK-Eng Enghave. FLUXNET; Technical University of Denmark (DTU) https://doi.org/10.18140/flx/1440153 (2016).

132. Olesen, J. FLUXNET2015 DK-Fou Foulum. FLUXNET; Danish Institute of Agricultural Sciences https://doi.org/10.18140/flx/1440154 (2016).

133. Ibrom, A. & Pilegaard, K. FLUXNET2015 DK-Sor Soroe. FLUXNET; Technical University of Denmark (DTU) https://doi.org/10.18140/flx/1440155 (2016).

134. Poveda, F. et al. FLUXNET2015 ES-Amo Amoladeras. FLUXNET; Estación Experimental de Zona Áridas (EEZA, CSIC) https://doi.org/10.18140/flx/1440156 (2016).

135. Reverter, B., Perez-Cañete, E. & Kowalski, A. FLUXNET2015 ES-LgS Laguna Seca. FLUXNET; Universidad de Granada https://doi.org/10.18140/flx/1440225 (2016).

136. Cañete, E. et al. FLUXNET2015 ES-LJ Llano de los Juanes. FLUXNET; University of Granada https://doi.org/10.18140/flx/1440157 (2016).

137. Reverter, B., Perez-Cañete, E. & Kowalski, A. FLUXNET2015 ES-Ln2 Lanjaron-Salvage logging. FLUXNET; Universidad de Granada https://doi.org/10.18140/flx/1440226 (2016).

138. Mammarella, I. et al. FLUXNET2015 FI-Hyy Hytiälä. FLUXNET; University of Helsinki https://doi.org/10.18140/flx/1440158 (2016).

139. Lohila, A. et al. FLUXNET2015 FI-Jok Jokioinen. FLUXNET; Finnish Meteorological Institute https://doi.org/10.18140/flx/1440159 (2016).

140. Lohila, A. et al. FLUXNET2015 FI-Let Lettosuo. FLUXNET; Finnish Meteorological Institute https://doi.org/10.18140/flx/1440227 (2016).

141. Aurela, M. et al. FLUXNET2015 FI-Lom Lompolojankka. FLUXNET; Finnish Meteorological Institute https://doi.org/10.18140/flx/1440228 (2016).

142. Aurela, M. et al. FLUXNET2015 FI-Sod Sodankyla. FLUXNET; Finnish Meteorological Institute https://doi.org/10.18140/flx/1440160 (2016).

143. Berveiller, D. et al. FLUXNET2015 FR-Fon Fontainebleau-Barbeau. FLUXNET; CNRS https://doi.org/10.18140/flx/1440161 (2016).

144. Buysse, P. et al. FLUXNET2015 FR-Gri Grignon. FLUXNET; French National Institute for Agricultural Research https://doi.org/10.18140/flx/1440162 (2016).

145. Berbigier, P. & Loustau, D. FLUXNET2015 FR-LBr Le Bray. FLUXNET; INRA - UMR ISPA https://doi.org/10.18140/flx/1440163 (2016).

146. Ourcival, J. FLUXNET2015 FR-Pue Puechabon. FLUXNET; CNRS https://doi.org/10.18140/flx/1440164 (2016).

147. Bonal, D. & Burban, B. FLUXNET2015 GF-Guy Guyaflux (French Guiana). FLUXNET; INRA https://doi.org/10.18140/flx/1440165 (2016).

148. Valentini, R. et al. FLUXNET2015 GH-Ank Ankasa. FLUXNET; Euro Mediterranean Center for Climate Change - Viterbo; University of Tuscia - Vietrbo, https://doi.org/10.18140/flx/1440229 (2016).

149. Hansen, B. FLUXNET2015 GL-NuF Nuuk Fen. FLUXNET; University of Copenhagen; University of Aarhus; Asiaq - Greenland Survey https://doi.org/10.18140/flx/1440222 (2016).

150. Lund, M., Jackowicz-Korczyński, M. & Abermann, J. FLUXNET2015 GL-ZaF Zackenberg Fen. FLUXNET; Aarhus University https://doi.org/10.18140/flx/1440223 (2016).

151. Lund, M., Jackowicz-Korczyński, M. & Abermann, J. FLUXNET2015 GL-ZaH Zackenberg Heath. FLUXNET; Aarhus University https://doi.org/10.18140/flx/1440224 (2016).

152. Magliulo, V. et al. FLUXNET2015 IT-BCi Borgo Cioffi. FLUXNET; CNR ISAFOM https://doi.org/10.18140/flx/1440166 (2016).

153. Sabbatini, S., Arriga, N. & Papale, D. FLUXNET2015 IT-CA1 Castel d'Asso1. FLUXNET; University of Tuscia - Vietrbo https://doi.org/10.18140/flx/1440230 (2016).

154. Sabbatini, S., Arriga, N., Gioli, B. & Papale, D. FLUXNET2015 IT-CA2 Castel d'Asso2. FLUXNET; CNR IBIMET; University of Tuscia - Vietrbo https://doi.org/10.18140/flx/1440231 (2016).
155. Sabbatini, S., Arriga, N., Matteucci, G. & Papale, D. FLUXNET2015 IT-CA3 Castel d'Asso 3. FLUXNET; University of Tuscia - Vietrbo; CNR IBAF https://doi.org/10.18140/flx/1440232 (2016).

156. Matteucci, G. FLUXNET2015 IT-Col Collelongo. FLUXNET; Istituto di Ecologia e Idrologia Forestale CNR https://doi.org/10.18140/flx/1440167 (2016).

157. Fares, S., Savi, F. & Conte, A. FLUXNET2015 IT-Cp2 Castelporziano 2. FLUXNET; Council for Agricultural Research and Economics https://doi.org/10.18140/flx/1440233 (2016).

158. Valentini, R. et al. FLUXNET2015 IT-Cp2 Castelporziano. FLUXNET; University of Tuscia - Vietrbo https://doi.org/10.18140/flx/1440168 (2016).

159. Gruening, C., Goded, I., Cescatti, A. & Pokorska, O. FLUXNET2015 IT-Isp Ispra ABC-IS. FLUXNET; European Commission - Joint Research Centre https://doi.org/10.18140/flx/1440234 (2016).

160. Cescatti, A., Marcolla, B., Zorer, R. & Gianelle, D. FLUXNET2015 IT-La2 Lavarone2. FLUXNET; Centro di Ecologia Alpina https://doi.org/10.18140/flx/1440235 (2016).

161. Gianelle, D., Zampedri, R., Cavagna, M. & Sottocornola, M. FLUXNET2015 IT-Lav Lavarone. FLUXNET; Edmund Mach Foundation https://doi.org/10.18140/flx/1440169 (2016).

162. Gianelle, D., Cavagna, M., Zampedri, R. & Marcolla, B. FLUXNET2015 IT-MBo Monte Bondone. FLUXNET; Edmund Mach Foundation https://doi.org/10.18140/flx/1440170 (2016).

163. Spano, D. et al. FLUXNET2015 IT-Noe Arca di Noe - Le Prigionette. FLUXNET; University of Sassari; CNR-Ibimet Sassari https://doi.org/10.18140/flx/1440171 (2016).

164. Manca, G. & Goded, I. FLUXNET2015 IT-PT1 Parco Ticino forest. FLUXNET; European Commission - DG Joint Research Centre https://doi.org/10.18140/flx/1440172 (2016).

165. Montagnani, L. & Minerbi, S. FLUXNET2015 IT-Ren Renon. FLUXNET; Autonomous Province of Bolzano, Forest Services https://doi.org/10.18140/flx/1440173 (2016).

166. Valentini, R. et al. FLUXNET2015 IT-Ro1 Roccarespampani 1. FLUXNET; University of Tuscia - Vietrbo https://doi.org/10.18140/flx/1440174 (2016).

167. Papale, D. et al. FLUXNET2015 IT-Ro2 Roccarespampani 2. FLUXNET; University of Tuscia - Vietrbo https://doi.org/10.18140/flx/1440175 (2016).

168. Gruening, C., Goded, I., Cescatti, A. & Pokorska, O. FLUXNET2015 IT-SR2 San Rossore 2. FLUXNET; European Commission - Joint Research Centre https://doi.org/10.18140/flx/1440236 (2016).

169. Gruening, C. et al. FLUXNET2015 IT-SRo San Rossore. FLUXNET; European Commission - Joint Research Centre https://doi.org/10.18140/flx/1440176 (2016).

170. Cremonese, E., Galvagno, M., Di Cella, U. & Migliavacca, M. FLUXNET2015 IT-Tor Torgnon. FLUXNET; Environmental Protection Agency of Aosta Valley https://doi.org/10.18140/flx/1440237 (2016).

171. Kotani, A. FLUXNET2015 JP-MBF Moshiri Birch Forest Site. FLUXNET; Nagoya University https://doi.org/10.18140/flx/1440238 (2016).

172. Kotani, A. FLUXNET2015 JP-SMF Seto Mixed Forest Site. FLUXNET; Nagoya University https://doi.org/10.18140/flx/1440239 (2016).

173. Kosugi, Y. & Takanashi, S. FLUXNET2015 MY-PSO Pasoh Forest Reserve (PSO). FLUXNET; FRIM (Forest Research Institute of Malaysia); Kyoto University https://doi.org/10.18140/flx/1440240 (2016).

174. Dolman, H. et al. FLUXNET2015 NL-Hor Horstermeer. FLUXNET; Vrije Universiteit Amsterdam https://doi.org/10.18140/flx/1440177 (2016).

175. Moors, E. & Elbers, J. FLUXNET2015 NL-Loo Loobos. FLUXNET; ALTERRA/Wageningen Environmental Research https://doi.org/10.18140/flx/1440178 (2016).

176. Wolf, S., Eugster, W. & Buchmann, N. FLUXNET2015 PA-SPn Sardinilla Plantation. FLUXNET; ETH Zurich https://doi.org/10.18140/flx/1440180 (2016).

177. Wolf, S., Eugster, W. & Buchmann, N. FLUXNET2015 PA-SPs Sardinilla-Pasture. FLUXNET; ETH Zurich https://doi.org/10.18140/flx/1440179 (2016).

178. Merbold, L., Rebmann, C. & Corradi, C. FLUXNET2015 RU-Che Cherski. FLUXNET; Max-Planck Institute for Biogeochemistry https://doi.org/10.18140/flx/1440181 (2016).

179. Dolman, H. et al. FLUXNET2015 RU-Cok Chokurdakh. FLUXNET; Vrije Universiteit Amsterdam https://doi.org/10.18140/flx/1440182 (2016).

180. Varlagin, A., Kurbatova, J. & Vygodskaya, N. FLUXNET2015 RU-Fyo Fyodorovskoye. FLUXNET; A.N. Severtsov Institute of Ecology and Evolution https://doi.org/10.18140/flx/1440183 (2016).

181. Belelli, L., Papale, D. & Valentini, R. FLUXNET2015 RU-Ha1 Hakasia steppe. FLUXNET; University of Tuscia - Vietrbo https://doi.org/10.18140/flx/1440184 (2016).

182. Ardö, J., El Tahir, B. & ElKhidir, H. FLUXNET2015 SD-Dem Demokeya. FLUXNET; LUND UNIVERSITY https://doi.org/10.18140/flx/1440186 (2016).

183. Christensen, T. FLUXNET2015 SJ-Adv Adventdalen. FLUXNET; NATEKO; Lund University https://doi.org/10.18140/flx/1440241 (2016).

184. Boike, J. et al. FLUXNET2015 SJ-Blv Bayelva, Spitsbergen. FLUXNET; University of Oslo, Department of Geosciences, 0316 OSLO, Norway; Universität Bayreuth, Department of Earth Sciences, 95440 Bayreuth, Germany; Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Periglacial Research Unit, 14473 Potsdam, Germany, https://doi.org/10.18140/flx/1440242 (2016).

185. Tagesson, T., Ardö, J. & Fensholt, R. FLUXNET2015 SN-Dhr Dahra. FLUXNET; Lund University https://doi.org/10.18140/flx/1440246 (2016).

186. Billesbach, D., Bradford, J. & Torn, M. FLUXNET2015 US-AR1 ARM USDA UNL OSU Woodward Switchgrass 1. FLUXNET; Lawrence Berkeley National Lab; U.S. Department of Agriculture; University of Nebraska https://doi.org/10.18140/flx/1440103 (2016).

187. Billesbach, D., Bradford, J. & Torn, M. FLUXNET2015 US-AR2 ARM USDA UNL OSU Woodward Switchgrass 2. FLUXNET; Lawrence Berkeley National Lab; U.S. Department of Agriculture; University of Nebraska https://doi.org/10.18140/flx/1440104 (2016).

188. Torn, M. FLUXNET2015 US-ARb ARM Southern Great Plains burn site- Lamont. FLUXNET; Lawrence Berkeley National Laboratory https://doi.org/10.18140/flx/1440064 (2016).

189. Torn, M. FLUXNET2015 US-ARc ARM Southern Great Plains control site- Lamont. FLUXNET; Lawrence Berkeley National Laboratory https://doi.org/10.18140/flx/1440065 (2016).

190. Biraud, S., Fischer, M., Chan, S. & Torn, M. FLUXNET2015 US-ARM ARM Southern Great Plains site- Lamont. FLUXNET; Lawrence Berkeley National Laboratory https://doi.org/10.18140/flx/1440066 (2016).

191. Zona, D. & Oechel, W. FLUXNET2015 US-Atq Atqasuk. FLUXNET; San Diego State University https://doi.org/10.18140/flx/1440067 (2016).

192. Goldstein, A. FLUXNET2015 US-Blo Blodgett Forest. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440068 (2016).

193. Bowling, D. FLUXNET2015 US-Cop Corral Pocket. FLUXNET; University of Utah https://doi.org/10.18140/flx/1440100 (2016).

194. Chen, J. & Chu, H. FLUXNET2015 US-CRT Curtice Walter-Berger cropland. FLUXNET; University of Toledo/Michigan State University https://doi.org/10.18140/flx/1440117 (2016).
195. Massman, B. FLUXNET2015 US-GBT GLEES Brooklyn Tower. FLUXNET; USDA Forest Service https://doi.org/10.18140/flx/1440118 (2016).

196. Massman, B. FLUXNET2015 US-GLE GLEES. FLUXNET; USDA Forest Service https://doi.org/10.18140/flx/1440069 (2016).

197. Meyers, T. FLUXNET2015 US-Goo Goodwin Creek. FLUXNET; NOAA/ARL https://doi.org/10.18140/flx/1440070 (2016).

198. Munger, J. FLUXNET2015 US-Ha1 Harvard Forest EMS Tower (HFR1). FLUXNET; Harvard University https://doi.org/10.18140/flx/1440071 (2016).

199. Matamala, R. FLUXNET2015 US-IB2 Fermi National Accelerator Laboratory- Batavia (Prairie site). FLUXNET; Argonne National Laboratory https://doi.org/10.18140/flx/1440072 (2016).

200. Zona, D. & Oechel, W. FLUXNET2015 US-Ivo Ivtotuk. FLUXNET; San Diego State University https://doi.org/10.18140/flx/1440073 (2016).

201. Drake, B. & Hinkle, R. FLUXNET2015 US-KS1 Kennedy Space Center (slash pine). FLUXNET; Smithsonian Environmental Research Center; University of Central Florida https://doi.org/10.18140/flx/1440074 (2016).

202. Drake, B. & Hinkle, R. FLUXNET2015 US-KS2 Kennedy Space Center (scrub oak). FLUXNET; Smithsonian Environmental Research Center; University of Central Florida https://doi.org/10.18140/flx/1440075 (2016).

203. Fares, S. FLUXNET2015 US-Lin Lindcove Orange Orchard. FLUXNET; Entecra https://doi.org/10.18140/flx/1440107 (2016).

204. Desai, A. FLUXNET2015 US-Los Lost Creek. FLUXNET; University of Wisconsin https://doi.org/10.18140/flx/1440076 (2016).

205. Meyers, T. FLUXNET2015 US-LWW Little Washita Watershed. FLUXNET; NOAA/ARL https://doi.org/10.18140/flx/1440077 (2016).

206. Law, B. FLUXNET2015 US-Me1 Metolius - Eyerly burn. FLUXNET; Oregon State University https://doi.org/10.18140/flx/1440078 (2016).

207. Law, B. FLUXNET2015 US-Me2 Metolius mature ponderosa pine. FLUXNET; Oregon State University https://doi.org/10.18140/flx/1440079 (2016).

208. Law, B. FLUXNET2015 US-Me3 Metolius-second young aged pine. FLUXNET; Oregon State University https://doi.org/10.18140/flx/1440080 (2016).

209. Law, B. FLUXNET2015 US-Me4 Metolius-old aged ponderosa pine. FLUXNET; Oregon State University https://doi.org/10.18140/flx/1440081 (2016).

210. Law, B. FLUXNET2015 US-Me5 Metolius-first young aged pine. FLUXNET; Oregon State University https://doi.org/10.18140/flx/1440082 (2016).

211. Law, B. FLUXNET2015 US-Me6 Metolius Young Pine Burn. FLUXNET; Oregon State University https://doi.org/10.18140/flx/1440099 (2016).

212. Novick, K. & Phillips, R. FLUXNET2015 US-MMS Morgan Monroe State Forest. FLUXNET; Indiana University https://doi.org/10.18140/flx/1440083 (2016).

213. Sturtevant, C. et al. FLUXNET2015 US-Myb Mayberry Wetland. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440105 (2016).

214. Suyker, A. FLUXNET2015 US-Ne1 Mead - irrigated continuous maize site. FLUXNET; University of Nebraska - Lincoln https://doi.org/10.18140/flx/1440084 (2016).

215. Suyker, A. FLUXNET2015 US-Ne2 Mead - irrigated maize-soybean rotation site. FLUXNET; University of Nebraska - Lincoln https://doi.org/10.18140/flx/1440085 (2016).

216. Suyker, A. FLUXNET2015 US-Ne3 Mead - rainfed maize-soybean rotation site. FLUXNET; University of Nebraska - Lincoln https://doi.org/10.18140/flx/1440086 (2016).

217. Blanken, P. et al. FLUXNET2015 US-NR1 Niwot Ridge Forest (LTER NWT1). FLUXNET; University of Colorado https://doi.org/10.18140/flx/1440087 (2016).

218. Chen, J., Chu, H. & Noormets, A. FLUXNET2015 US-Oho Oak Openings. FLUXNET; University of Toledo/Michigan State University https://doi.org/10.18140/flx/1440088 (2016).

219. Bohrer, G. FLUXNET2015 US-ORv Olentangy River Wetland Research Park. FLUXNET; The Ohio State University https://doi.org/10.18140/flx/1440102 (2016).

220. Desai, A. FLUXNET2015 US-PFa Park Falls/WLEF. FLUXNET; University of Wisconsin https://doi.org/10.18140/flx/1440089 (2016).

221. Kobayashi, H. & Suzuki, R. FLUXNET2015 US-Prr Poker Flat Research Range Black Spruce Forest. FLUXNET; Japan Agency for Marine-Earth Science and Technology https://doi.org/10.18140/flx/1440113 (2016).

222. Kurc, S. FLUXNET2015 US-SRC Santa Rita Creosote. FLUXNET; University of Arizona https://doi.org/10.18140/flx/1440098 (2016).

223. Scott, R. FLUXNET2015 US-SRG Santa Rita Grassland. FLUXNET; United States Department of Agriculture https://doi.org/10.18140/flx/1440114 (2016).

224. Scott, R. FLUXNET2015 US-SRM Santa Rita Mesquite. FLUXNET; United States Department of Agriculture https://doi.org/10.18140/flx/1440090 (2016).

225. Ewers, B. & Pendall, E. FLUXNET2015 US-Sta Saratoga. FLUXNET; University of Wyoming https://doi.org/10.18140/flx/1440115 (2016).

226. Desai, A. FLUXNET2015 US-Syv Sylvania Wilderness Area. FLUXNET; University of Wisconsin https://doi.org/10.18140/flx/1440091 (2016).

227. Baldocchi, D. & Ma, S. FLUXNET2015 US-Ton Tonzi Ranch. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440092 (2016).

228. Szutu, D., Baldocchi, D., Eichelmann, E. & Knox, S. FLUXNET2015 US-Tw1 Twitchell Wetland West Pond. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440108 (2016).

229. Baldocchi, D. FLUXNET2015 US-Tw2 Twitchell Corn. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440109 (2016).

230. Szutu, D. & Baldocchi, D. FLUXNET2015 US-Tw3 Twitchell Alfalfa. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440110 (2016).

231. Sanchez, C. et al. FLUXNET2015 US-Tw4 Twitchell East End Wetland. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440111 (2016).

232. Baldocchi, D. FLUXNET2015 US-Twt Twitchell Island. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440106 (2016).

233. Gough, C., Bohrer, G. & Curtis, P. FLUXNET2015 US-UMB Univ. of Mich. Biological Station. FLUXNET; Ohio State University; Virginia Commonwealth University https://doi.org/10.18140/flx/1440093 (2016).

234. Gough, C., Bohrer, G. & Curtis, P. FLUXNET2015 US-UMd UMBS Disturbance. FLUXNET; Ohio State University; Virginia Commonwealth University https://doi.org/10.18140/flx/1440101 (2016).

235. Baldocchi, D., Ma, S. & Xu, L. FLUXNET2015 US-Var Vaira Ranch- Ione. FLUXNET; University of California, Berkeley https://doi.org/10.18140/flx/1440094 (2016).

236. Desai, A. FLUXNET2015 US-WCr Willow Creek. FLUXNET; University of Wisconsin https://doi.org/10.18140/flx/1440095 (2016).

237. Scott, R. FLUXNET2015 US-Whs Walnut Gulch Lucky Hills Shrub. FLUXNET; United States Department of Agriculture https://doi.org/10.18140/flx/1440097 (2016).

238. Chen, J. FLUXNET2015 US-Wi0 Young red pine (YRP). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440055 (2016).
239. Chen, J. FLUXNET2015 US-Wi1 Intermediate hardwood (IHW). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440054 (2016).

240. Chen, J. FLUXNET2015 US-Wi2 Intermediate red pine (IRP). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440056 (2016).

241. Chen, J. FLUXNET2015 US-Wi3 Mature hardwood (MHW). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440057 (2016).

242. Chen, J. FLUXNET2015 US-Wi4 Mature red pine (MRP). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440058 (2016).

243. Chen, J. FLUXNET2015 US-Wi5 Mixed young jack pine (MYJP). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440059 (2016).

244. Chen, J. FLUXNET2015 US-Wi6 Pine barrens #1 (PB1). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440060 (2016).

245. Chen, J. FLUXNET2015 US-Wi7 Red pine clearcut (RPCC). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440061 (2016).

246. Chen, J. FLUXNET2015 US-Wi8 Young hardwood clearcut (YHW). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440062 (2016).

247. Chen, J. FLUXNET2015 US-Wi9 Young Jack pine (YJP). FLUXNET; Michigan State University https://doi.org/10.18140/flx/1440063 (2016).

248. Scott, R. FLUXNET2015 US-Wkg Walnut Gulch Kendall Grasslands. FLUXNET; United States Department of Agriculture https://doi.org/10.18140/flx/1440096 (2016).

249. Chen, J. & Chu, H. FLUXNET2015 US-WPT Winous Point North Marsh. FLUXNET; University of Toledo/Michigan State University https://doi.org/10.18140/flx/1440116 (2016).

250. Kutsch, W., Merbold, L. & Kolle, O. FLUXNET2015 ZM-Mon Mongu. FLUXNET; Max-Planck Institute for Biogeochemistry https://doi.org/10.18140/flx/1440189 (2016).

251. Harris, I., Jones, P. D., Osborn, T. J. & Lister, D. H. Updated high-resolution grids of monthly climatic observations - the CRU TS3.10 Dataset. *Int. J. Climatol.* **34**, 623–642 (2014).

252. Channan, S., Collins, K. & Emanuel, W. Global mosaics of the standard MODIS land cover type data. University of Maryland and the Pacific Northwest National Laboratory, College Park, Maryland, USA (2014).

253. Baldocchi, D., Chu, H. & Reichstein, M. Inter-annual variability of net and gross ecosystem carbon fluxes: A review. *Agr. Forest Meteorol.* **249**, 520–533 (2018).

254. Baldocchi, D. D. How eddy covariance flux measurements have contributed to our understanding of Global Change Biology. *Glob. Change Biol.* **26**, 242–260 (2019).

255. Law, B. E. et al. Terrestrial Carbon Observations: Protocols for Vegetation Sampling and Data Submission. Report No.55 (Global Terrestrial Observing System, FAO, Rome, 2008).

256. Papale, D., Canfora, E. & Polidori, D. ICOS Ecosystem Instructions for Use the ICOS BADM (Version 20171013), ICOS Ecosystem Thematic Centre, https://doi.org/10.18160/6m8s-fy7m (2017).

257. Tramontana, G. et al. Predicting carbon dioxide and energy fluxes across global FLUXNET sites with regression algorithms. *Biogeosciences* **13**, 4291–4313 (2016).

258. Keenan, T. F. et al. Widespread inhibition of daytime ecosystem respiration. *Nat. Ecol. Evol.* **3**, 407–415 (2019).

259. Campioli, M. et al. Evaluating the convergence between eddy-covariance and biometric methods for assessing carbon budgets of forests. *Nat. Commun.* **7**, 13717 (2016).

## Acknowledgements

We thank the many people who helped with generating high quality data through the years at all sites, and the funding sources to the sites and networks, for the data collection, curation, and sharing that make this dataset possible. The eddy covariance data were acquired and shared by the following networks: AmeriFlux, AfriFlux, AsiaFlux, CarboAfrica, CarboEuropeIP, CarboItaly, CarboMont, ChinaFlux, Fluxnet-Canada, GreenGrass, ICOS, KoFlux, LBA, NECC, OzFlux-TERN, TCOS-Siberia, Swiss FluxNet and USCCC. The ERA-Interim reanalysis data were provided by ECMWF and processed by Laboratoire des sciences du climat et de l'environnement (LSCE). The FLUXNET eddy covariance data processing and harmonization was carried out by the European Fluxes Database Cluster and the AmeriFlux Management Project (with support by European Union H2020 projects and U.S. Department of Energy Office of Science, respectively), with contributions from the Carbon Dioxide Information Analysis Center, ICOS Ecosystem Thematic Centre, and OzFlux, ChinaFlux, and AsiaFlux offices.

## Author contributions

Dario Papale, Gilberto Pastorello, Margaret Torn, and Deb Agarwal conceived of and organized the FLUXNET2015 dataset. Eleonora Canfora, You-Wei Cheah, Danielle Christianson, Housen Chu, Dario Papale, Gilberto Pastorello, Diego Polidori, Cristina Poindexter, and Carlo Trotta were responsible for quality checking, post-processing, and creating data products. You-Wei Cheah, Danielle Christianson, Abdelrahman Elbashandy, Marty Humphrey, Gilberto Pastorello, Diego Polidori, Markus Reichstein, Alessio Ribeca, Carlo Trotta, Catharine van Ingen and Nicolas Vuichard were responsible for code and software preparation and implementation (processing pipeline and data distribution platform). Eleonora Canfora, You-Wei Cheah, Jiquan Chen, Danielle Christianson, Housen Chu, Peter Isaac, Dario Papale and Leiming Zhang managed the flux data and metadata collections. Gilberto Pastorello, Dario Papale, Deb Agarwal, Sebastien Biraud, You-Wei Cheah, Danielle Christianson, Housen Chu, and Margaret Torn conceived of the paper and prepared the first draft that was reviewed and commented on by all the coauthors. All the coauthors collected, processed and contributed the data, also participated in the quality assessment and correction of errors. The link between sites and coauthors is provided in Supplementary Table SM9.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41597-020-0534-3.

Correspondence and requests for materials should be addressed to G.P. or D.P.
Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files associated with this article.

This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply 2020, corrected publication 2021

Gilberto Pastorello<sup>1,∗</sup>, Carlo Trotta<sup>2</sup>, Eleonora Canfora<sup>2,3</sup>, Housen Chu<sup>4</sup>,
Danielle Christianson<sup>1</sup>, You-Wei Cheah<sup>1</sup>, Cristina Poindexter<sup>5</sup>, Jiquan Chen<sup>6</sup>,
Abdelrahman Elbashandy<sup>1</sup>, Marty Humphrey<sup>7</sup>, Peter Isaac<sup>8</sup>, Diego Polidori<sup>2,3</sup>, Markus Reichstein<sup>9</sup>,
Alessio Ribeca<sup>2,3</sup>, Catharine van Ingen<sup>1</sup>, Nicolas Vuichard<sup>10</sup>, Leiming Zhang<sup>11</sup>, Brian Amiro<sup>12</sup>,
Christof Ammann<sup>13</sup>, M. Altaf Arain<sup>14</sup>, Jonas Ardö<sup>15</sup>, Timothy Arkebauer<sup>16</sup>, Stefan K. Arndt<sup>17</sup>,
Nicola Arriga<sup>18,19</sup>, Marc Aubinet<sup>20</sup>, Mika Aurela<sup>21</sup>, Dennis Baldocchi<sup>22</sup>, Alan Barr<sup>23,24</sup>,
Eric Beamesderfer<sup>14</sup>, Luca Belelli Marchesini<sup>25,26</sup>, Onil Bergeron<sup>27</sup>, Jason Beringer<sup>28</sup>,
Christian Bernhofer<sup>29</sup>, Daniel Berveiller<sup>30</sup>, Dave Billesbach<sup>31</sup>, Thomas Andrew Black<sup>32</sup>,
Peter D. Blanken<sup>33</sup>, Gil Bohrer<sup>34</sup>, Julia Boike<sup>35,36</sup>, Paul V. Bolstad<sup>37</sup>, Damien Bonal<sup>38</sup>,
Jean-Marc Bonnefond<sup>39</sup>, David R. Bowling<sup>40</sup>, Rosvel Bracho<sup>41</sup>, Jason Brodeur<sup>42</sup>,
Christian Brümmer<sup>43</sup>, Nina Buchmann<sup>44</sup>, Benoit Burban<sup>45</sup>, Sean P. Burns<sup>33,46</sup>,
Pauline Buysse<sup>47</sup>, Peter Cale<sup>48</sup>, Mauro Cavagna<sup>25</sup>, Pierre Cellier<sup>47</sup>, Shiping Chen<sup>49</sup>, Isaac Chini<sup>25</sup>,
Torben R. Christensen<sup>50</sup>, James Cleverly<sup>51,52</sup>, Alessio Collalti<sup>2,53</sup>, Claudia Consalvo<sup>2,54</sup>,
Bruce D. Cook<sup>55</sup>, David Cook<sup>56</sup>, Carole Coursolle<sup>57,58</sup>, Edoardo Cremonese<sup>59</sup>,
Peter S. Curtis<sup>60</sup>, Ettore D'Andrea<sup>53</sup>, Humberto da Rocha<sup>61</sup>, Xiaoqin Dai<sup>11</sup>, Kenneth J. Davis<sup>62</sup>,
Bruno De Cinti<sup>63</sup>, Agnes de Grandcourt<sup>64</sup>, Anne De Ligne<sup>20</sup>, Raimundo C. De Oliveira<sup>65</sup>,
Nicolas Delpierre<sup>30</sup>, Ankur R. Desai<sup>66</sup>, Carlos Marcelo Di Bella<sup>67</sup>, Paul di Tommasi<sup>53</sup>,
Han Dolman<sup>68</sup>, Francisco Domingo<sup>69</sup>, Gang Dong<sup>70</sup>, Sabina Dore<sup>71</sup>, Pierpaolo Duce<sup>72</sup>,
Eric Dufrêne<sup>30</sup>, Allison Dunn<sup>73</sup>, Jiří Dušek<sup>74</sup>, Derek Eamus<sup>51</sup>, Uwe Eichelmann<sup>29</sup>,
Hatim Abdalla M. ElKhidir<sup>75</sup>, Werner Eugster<sup>44</sup>, Cacilia M. Ewenz<sup>76</sup>, Brent Ewers<sup>77</sup>,
Daniela Famulari<sup>53</sup>, Silvano Fares<sup>78,79</sup>, Iris Feigenwinter<sup>44</sup>, Andrew Feitz<sup>80</sup>, Rasmus Fensholt<sup>81</sup>,
Gianluca Filippa<sup>59</sup>, Marc Fischer<sup>82</sup>, John Frank<sup>83</sup>, Marta Galvagno<sup>59</sup>, Mana Gharun<sup>44</sup>,
Damiano Gianelle<sup>25</sup>, Bert Gielen<sup>18</sup>, Beniamino Gioli<sup>84</sup>, Anatoly Gitelson<sup>85</sup>, Ignacio Goded<sup>19</sup>,
Mathias Goeckede<sup>86</sup>, Allen H. Goldstein<sup>22</sup>, Christopher M. Gough<sup>87</sup>, Michael L. Goulden<sup>88</sup>,
Alexander Graf<sup>89</sup>, Anne Griebel<sup>17</sup>, Carsten Gruening<sup>19</sup>, Thomas Grünwald<sup>29</sup>,
Albin Hammerle<sup>90</sup>, Shijie Han<sup>91,92</sup>, Xingguo Han<sup>49</sup>, Birger Ulf Hansen<sup>81</sup>, Chad Hanson<sup>93</sup>,
Juha Hatakka<sup>21</sup>, Yongtao He<sup>11,94</sup>, Markus Hehn<sup>29</sup>, Bernard Heinesch<sup>20</sup>, Nina Hinko-Najera<sup>95</sup>,
Lukas Hörtnagl<sup>44</sup>, Lindsay Hutley<sup>96</sup>, Andreas Ibrom<sup>97</sup>, Hiroki Ikawa<sup>98</sup>, Marcin Jackowicz-
Korczynski<sup>15,50</sup>, Dalibor Janouš<sup>74</sup>, Wilma Jans<sup>99</sup>, Rachhpal Jassal<sup>32</sup>, Shicheng Jiang<sup>100</sup>,
Tomomichi Kato<sup>101,102</sup>, Myroslava Khomik<sup>14,103</sup>, Janina Klatt<sup>104</sup>, Alexander Knohl<sup>105,106</sup>,
Sara Knox<sup>107</sup>, Hideki Kobayashi<sup>108</sup>, Georgia Koerber<sup>109</sup>, Olaf Kolle<sup>86</sup>, Yoshiko Kosugi<sup>110</sup>,
Ayumi Kotani<sup>111</sup>, Andrew Kowalski<sup>112</sup>, Bart Kruijt<sup>113</sup>, Julia Kurbatova<sup>114</sup>, Werner L. Kutsch<sup>115</sup>,
Hyojung Kwon<sup>93</sup>, Samuli Launiainen<sup>116</sup>, Tuomas Laurila<sup>21</sup>, Bev Law<sup>93</sup>, Ray Leuning<sup>197</sup>,
Yingnian Li<sup>117</sup>, Michael Liddell<sup>118</sup>, Jean-Marc Limousin<sup>119</sup>, Marryanna Lion<sup>120</sup>, Adam J. Liska<sup>31</sup>,
Annalea Lohila<sup>21,121</sup>, Ana López-Ballesteros<sup>122</sup>, Efrén López-Blanco<sup>50</sup>, Benjamin Loubet<sup>47</sup>,
Denis Loustau<sup>39</sup>, Antje Lucas-Moffat<sup>43,123</sup>, Johannes Lüers<sup>124,125</sup>, Siyan Ma<sup>22</sup>,
Craig Macfarlane<sup>126</sup>, Vincenzo Magliulo<sup>53</sup>, Regine Maier<sup>44</sup>, Ivan Mammarella<sup>121</sup>,
Giovanni Manca<sup>19</sup>, Barbara Marcolla<sup>25</sup>, Hank A. Margolis<sup>58</sup>, Serena Marras<sup>3,127</sup>, William Massman<sup>83</sup>,
Mikhail Mastepanov<sup>50,128</sup>, Roser Matamala<sup>56</sup>, Jaclyn Hatala Matthes<sup>129</sup>,
Francesco Mazzenga¹³⁰, Harry McCaughey¹³¹, Ian McHugh¹⁷, Andrew M. S. McMillan¹³²,
Lutz Merbold¹³³, Wayne Meyer¹⁰⁹, Tilden Meyers¹³⁴, Scott D. Miller¹³⁵, Stefano Minerbi¹³⁶,
Uta Moderow²⁹, Russell K. Monson¹³⁷, Leonardo Montagnani¹³⁶,¹³⁸, Caitlin E. Moore¹³⁹,
Eddy Moors¹⁴⁰,¹⁴¹, Virginie Moreaux³⁹,¹⁴², Christine Moureaux²⁰, J. William Munger¹⁴³,¹⁴⁴,
Taro Nakai¹⁴⁵,¹⁴⁶, Johan Neirynck¹⁴⁷, Zoran Nesic³², Giacomo Nicolini²,³, Asko Noormets¹⁴⁸,
Matthew Northwood¹⁴⁹, Marcelo Nosetto¹⁵⁰,¹⁵¹, Yann Nouvellon⁶⁴,¹⁵², Kimberly Novick¹⁵³,
Walter Oechel¹⁵⁴,¹⁵⁵, Jørgen Eivind Olesen¹⁵⁶,¹⁵⁷, Jean-Marc Ourcival¹¹⁹, Shirley A. Papuga¹⁵⁸,
Frans-Jan Parmentier¹⁵,¹⁵⁹, Eugenie Paul-Limoges¹⁶⁰, Marian Pavelka⁷⁴, Matthias Peichl¹⁶¹,
Elise Pendall¹⁶², Richard P. Phillips¹⁶³, Kim Pilegaard⁹⁷, Norbert Pirk¹⁵,¹⁶⁴, Gabriela Posse¹⁶⁵,
Thomas Powell⁴, Heiko Prasse²⁹, Suzanne M. Prober¹⁶⁴, Serge Rambal¹¹⁹, Üllar Rannik¹²¹,
Naama Raz-Yaseef⁴, Corinna Rebmann¹⁶⁶, David Reed¹⁶⁷, Victor Resco de Dios¹⁶²,¹⁶⁸,
Natalia Restrepo-Coupe¹³⁷, Borja R. Reverter¹⁶⁹, Marilyn Roland¹⁸, Simone Sabbatini²,
Torsten Sachs¹⁷⁰, Scott R. Saleska¹³⁷, Enrique P. Sánchez-Cañete¹¹²,¹¹⁷, Zulia M. Sanchez-
Mejia¹⁷², Hans Peter Schmid¹⁰⁴, Marius Schmidt⁸⁹, Karl Schneider¹⁷³, Frederik Schrader¹⁴³,
Ivan Schroder¹⁷⁴, Russell L. Scott¹⁷⁵, Pavel Sedlák⁷⁴,¹⁷⁶, Penélope Serrano-Ortíz¹⁷¹,¹⁷⁷,
Changliang Shao¹⁷⁸, Peili Shi¹¹, Ivan Shironya¹¹⁴, Lukas Siebicke¹⁰⁵, Ladislav Šigut⁷⁴,
Richard Silberstein²⁸,¹⁷⁹, Costantino Sirca³,¹²⁷, Donatella Spano³,¹²⁷, Rainer Steinbrecher¹⁰⁴,
Robert M. Stevens¹⁸⁰, Cove Sturtevant¹⁸¹, Andy Suyker⁸⁵, Torbern Tagesson¹⁵,⁸¹,
Satoru Takanashi¹⁸², Yanhong Tang¹⁸³, Nigel Tapper¹⁸⁴, Jonathan Thom¹⁸⁵,
Michele Tomassucci²,¹⁸⁶, Juha-Pekka Tuovinen²¹, Shawn Urbanski¹⁸⁷, Riccardo Valentini²,³,
Michiel van der Molen¹⁸⁸, Eva van Gorsel¹⁸⁹, Ko van Huissteden⁶⁸, Andrej Varlagin¹¹⁴,
Joseph Verfaillie²², Timo Vesala¹²¹, Caroline Vincke¹⁹⁰, Domenico Vitale²,³,
Natalia Vygodskaya¹¹⁴, Jeffrey P. Walker¹⁹¹, Elizabeth Walter-Shea⁸⁵, Huimin Wang¹¹,
Robin Weber²², Sebastian Westermann¹⁶⁵, Christian Wille¹⁷⁰, Steven Wofsy¹⁴³,¹⁴⁴,
Georg Wohlfahrt⁹⁰, Sebastian Wolf⁴⁴, William Woodgate¹⁹²,¹⁹³, Yuelin Li¹⁹⁴,
Roberto Zampedri²⁵, Junhui Zhang⁹², Guoyi Zhou¹⁹⁵, Donatella Zona¹⁵⁴,¹⁹⁶, Deb Agarwal¹,
Sebastien Biraud⁴, Margaret Torn⁴ & Dario Papale²,³✉

¹Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA. ²DIBAF, University of Tuscia, Viterbo, 01100, Italy. ³Euro-Mediterranean Centre on Climate Change Foundation (CMCC), Lecce, 73100, Italy. ⁴Climate & Ecosystem Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA. ⁵Department of Civil Engineering, California State University, Sacramento, CA, 95819, USA. ⁶Department of Geography, Environment, and Spatial Sciences, Michigan State University, East Lansing, MI, 48823, USA. ⁷Department of Computer Science, University of Virginia, Charlottesville, VA, 22904, USA. ⁸TERN Ecosystem Processes, Menzies Creek, VIC3159, Australia. ⁹Max Planck Institute for Biogeochemistry, Jena, 07701, Germany. ¹⁰Laboratoire des Sciences du Climat et de l'Environnement (LSCE), CEA CNRS, UVSQ UPSACLAY, Gif sur Yvette, 91190, France. ¹¹Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, 100101, China. ¹²Department of Soil Science, University of Manitoba, Winnipeg, MB, R3T2N2, Canada. ¹³Department of Agroecology and Environment, Agroscope Research Institute, Zürich, 8046, Switzerland. ¹⁴School of Geography and Earth Sciences, McMaster University, L8S4K1, Hamilton, ON, Canada. ¹⁵Department of Physical Geography and Ecosystem Science, Lund University, Lund, 22362, Sweden. ¹⁶Department of Agronomy and Horticulture, University of Nebraska-Lincoln, Lincoln, NE, 68583, USA. ¹⁷School of Ecosystem and Forest Sciences, The University of Melbourne, Richmond, VIC3121, Australia. ¹⁸Department of Biology, Research Group PLECO, University of Antwerp, Antwerp, 2610, Belgium. ¹⁹Joint Research Centre, European Commission, Ispra, 21027, Italy. ²⁰TERRA Teaching and Research Center, University of Liege, Gembloux, B-5030, Belgium. ²¹Finnish Meteorological Institute, Helsinki, 00560, Finland. ²²ESPM, University of California Berkeley, Berkeley, CA, 94720, USA. ²³Global Institute for Water Security, University of Saskatchewan, Saskatoon, SK, S7N3H5, Canada. ²⁴Climate Research Division, Environment and Climate Change Canada, Saskatoon, SK, S7N3H5, Canada. ²⁵Department of Sustainable Agro-ecosystems and Bioresources, Research and Innovation Centre, Fondazione Edmund Mach, San Michele All'adige, 38010, Italy. ²⁶Department of Landscape Design and Sustainable Ecosystems, Agrarian-Technological Institute, RUDN University, Moscow, 117198, Russia. ²⁷Direction du marché du carbone, Ministère du Développement durable de l'Environnement et de la Lutte contre les changements climatiques, Québec, QC, G1R5V7, Canada. ²⁸School of Agriculture and Environment, University of Western Australia, Crawley, 6009, Australia. ²⁹Institute of Hydrology and Meteorology, Technische Universität Dresden, Tharandt, 01737, Germany. ³⁰Université Paris-Saclay, CNRS, AgroParisTech, Ecologie Systématique et Evolution, Orsay, 91405, France. ³¹Biological Systems Engineering, University of Nebraska-Lincoln, Lincoln, NE, 68583, USA. ³²Faculty of Land and Food Systems, University of British Columbia, Vancouver, BC, V6T1Z4, Canada. ³³Department of Geography, University of Colorado, Boulder, CO, 80309, USA. ³⁴Department of Civil, Environmental & Geodetic Engineering, Ohio State University, Columbus, OH, 43210, USA. ³⁵Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Potsdam, 14482, Germany. ³⁶Geography Department, Humboldt-Universität zu Berlin, Berlin, Germany. ³⁷Forest Resources, University of Minnesota, St Paul, MN, 55108, USA. ³⁸Université de Lorraine, AgroParisTech, INRAE, UMR Silva, Nancy, 54000, France. ³⁹ISPA, Bordeaux Sciences
Agro, INRAE, Villenave d'Ornon, 33140, France. ⁴⁰School of Biological Sciences, University of Utah, Salt Lake City, UT, 84112, USA. ⁴¹School of Forest Resources and Conservation, University of Florida, Gainesville, FL, 32611, USA. ⁴²McMaster University Library, McMaster University, Hamilton, ON, L8S4L6, Canada. ⁴³Thünen Institute of Climate-Smart Agriculture, Federal Research Institute of Rural Areas, Forestry and Fisheries, Braunschweig, 38116, Germany. ⁴⁴Department of Environmental Systems Science, ETH Zurich, Zurich, 8092, Switzerland. ⁴⁵INRAE UMR ECOFOG, Kourou, 97387, French Guiana. ⁴⁶Mesoscale and Microscale Meteorology Laboratory, National Center for Atmospheric Research, Boulder, CO, 80301, USA. ⁴⁷Université Paris-Saclay, INRAE, AgroParisTech, UMR ECOSYS, Thiverval-Grignon, 78850, France. ⁴⁸Australian Landscape Trust, Renmark, SA5341, Australia. ⁴⁹State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, China. ⁵⁰Department of Bioscience, Arctic Research Center, Aarhus University, Roskilde, 4000, Denmark. ⁵¹School of Life Sciences, University of Technology Sydney, Sydney, 2007, Australia. ⁵²Terrestrial Ecosystem Research Network TERN, University of Technology, Sydney, 2007, Australia. ⁵³Institute for Agricultural and Forestry Systems in the Mediterranean, National Research Council of Italy, Ercolano, 80056, Italy. ⁵⁴Research Institute on Terrestrial Ecosystems, National Research Council of Italy, Porano, 05010, Italy. ⁵⁵Biospheric Sciences Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD, 20771, USA. ⁵⁶Environmental Science Division, Argonne National Laboratory, Lemont, IL, 60439, USA. ⁵⁷Canadian Forest Service, Natural Resources Canada, Québec, QC, G1V4C7, Canada. ⁵⁸Centre d'étude de la forêt, Faculté de foresterie, de géographie et de géomatique, Université Laval, Québec, QC, G1V0A6, Canada. ⁵⁹Climate Change Unit, Environmental Protection Agency of Aosta Valley, Saint Christophe, 11020, Italy. ⁶⁰Department of Evolution, Ecology, and Organismal Biology, Ohio State University, Columbus, OH, 43210, USA. ⁶¹Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, São Paulo, SP, 01000-000, Brazil. ⁶²Department of Meteorology and Atmospheric Science, The Pennsylvania State University, University Park, PA, 16802, USA. ⁶³Institute of Research on Terrestrial Ecosystems, National Research Council of Italy, Montelibretti, 00010, Italy. ⁶⁴UMR Eco&Sols, CIRAD, Montpellier, 34060, France. ⁶⁵Pedology, Embrapa Amazonia Oriental, Belém, PA, 68020640, Brazil. ⁶⁶Atmospheric and Oceanic Sciences, University of Wisconsin-Madison, Madison, WI, 53706, USA. ⁶⁷Departamento de Métodos Cuantitativos y Sistemas de Información, Facultad de Agronomía, UBA, Buenos Aires, 1417, Argentina. ⁶⁸Department of Earth Sciences, Vrije Universiteit Amsterdam, Amsterdam, 1081 HV, The Netherlands. ⁶⁹Desertification and Geoecology Department, Experimental Station of Arid Zones, CSIC, Almería, 04120, Spain. ⁷⁰School of Life Science, Shanxi University, Taiyuan, 030006, China. ⁷¹HydroFocus, Davis, CA, 95618, USA. ⁷²Institute of BioEconomy, National Research Council of Italy, Sassari, 07100, Italy. ⁷³Department of Earth, Environment, and Physics, Worcester State University, Worcester, MA, 01602, USA. ⁷⁴Department of Matter and Energy Fluxes, Global Change Research Institute of the Czech Academy of Sciences, Brno, 60300, Czech Republic. ⁷⁵EIObeid Research Station, Agricultural Research Corporation, EIObeid, 51111, Sudan. ⁷⁶Airborne Research Australia, TERN Ecosystem Processes Central Node, Parafield, 5106, Australia. ⁷⁷Department of Botany, Program in Ecology, University of Wyoming, 1000 E. Univ. Ave, Laramie, WY, 82071, USA. ⁷⁸Institute of BioEconomy, National Research Council of Italy, Rome, 00100, Italy. ⁷⁹Research Centre for Forestry and Wood, Council for Agricultural Research and Economics, Rome, 00166, Italy. ⁸⁰Geoscience Australia, Canberra, 2601, Australia. ⁸¹Department of Geosciences and Natural Resource Management, University of Copenhagen, Copenhagen, 1350, Denmark. ⁸²Energy Analysis & Environmental Impacts Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA. ⁸³USDA Forest Service, Rocky Mountain Research Station, Fort Collins, CO, 80526, USA. ⁸⁴Institute of BioEconomy, National Research Council of Italy, Firenze, 50145, Italy. ⁸⁵School of Natural Resources, University of Nebraska-Lincoln, Lincoln, NE, 68583, USA. ⁸⁶Max Planck Institute for Biogeochemistry, Jena, 03641, Germany. ⁸⁷Department of Biology, Virginia Commonwealth University, Richmond, VA, 23284, USA. ⁸⁸Department of Earth System Science, University of California, Irvine, CA, 92697, USA. ⁸⁹Agrosphere (IBG3), Forschungszentrum Jülich, Jülich, 52428, Germany. ⁹⁰Department of Ecology, University of Innsbruck, Innsbruck, 6020, Austria. ⁹¹International Joint Research Laboratory for Global Change Ecology, School of Life Sciences, Henan University, Kaifeng, 450000, China. ⁹²Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang, 110016, China. ⁹³Department of Forest Ecosystems and Society, Oregon State University, Corvallis, OR, 97333, USA. ⁹⁴College of Resources and Environment, University of Chinese Academy of Sciences, Beijing, 100190, China. ⁹⁵School of Ecosystem and Forest Sciences, The University of Melbourne, Creswick, VIC3363, Australia. ⁹⁶Research Institute for the Environment and Livelihoods, Charles Darwin University, Darwin, 0909, Australia. ⁹⁷Department of Environmental Engineering, Technical University of Denmark (DTU), Kongens Lyngby, 2800, Denmark. ⁹⁸Institute for Agro-Environmental Sciences, National Agriculture and Food Research Organization, Tsukuba, 305-8604, Japan. ⁹⁹Wageningen Environmental Research, Wageningen University and Research, Wageningen, 6708PB, The Netherlands. ¹⁰⁰Key Laboratory of Vegetation Ecology, Ministry of Education, Northeast Normal University, Changchun, 130024, China. ¹⁰¹Research Faculty of Agriculture, Hokkaido University, Sapporo, 060-8589, Japan. ¹⁰²GI-Core, Hokkaido University, Sapporo, 060-0808, Japan. ¹⁰³Geography and Environmental Management, Waterloo, ON, N2L3G1, Canada. ¹⁰⁴Institute of Meteorology and Climate Research, Karlsruhe Institute of Technology, Garmisch-Partenkirchen, 82467, Germany. ¹⁰⁵Bioclimatology, University of Goettingen, Goettingen, 37077, Germany. ¹⁰⁶Centre of Biodiversity and Sustainable Land Use (CBL), University of Goettingen, Goettingen, 37077, Germany. ¹⁰⁷Department of Geography, The University of British Columbia, Vancouver, BC, V6T1Z2, Canada. ¹⁰⁸Research Institute for Global Change, Institute of Arctic Climate and Environment Research, Japan Agency for Marine-Earth Science and Technology, Yokoama, 236-0001, Japan. ¹⁰⁹Biological Sciences, University of Adelaide, Adelaide, SA5064, Australia. ¹¹⁰Graduate School of Agriculture, Kyoto University, Kyoto, 606-8502, Japan. ¹¹¹Graduate School of Bioagricultural Sciences, Nagoya University, Nagoya, 4648601, Japan. ¹¹²Department of Applied Physics, University of Granada, Granada, 18071, Spain. ¹¹³Water systems and Global Change group, Wageningen University, Wageningen, 6500, The Netherlands. ¹¹⁴A.N. Severtsov Institute of Ecology and Evolution, Russian Academy of Sciences, Moscow, 119071, Russia. ¹¹⁵Head Office, Integrated Carbon Observation System (ICOS ERIC), Helsinki, 00560, Finland. ¹¹⁶Natural Resources Institute Finland, Helsinki, 00790, Finland. ¹¹⁷Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy
of Sciences, Xining, 810008, China. ¹¹⁸Centre for Tropical Environmental Sustainability Studies, James Cook University, Cairns, 4878, Australia. ¹¹⁹CEFE, CNRS, Univ Montpellier, Montpellier, 34293, France. ¹²⁰Forestry and Environment Division, Forest Research Institute Malaysia (FRIM), Kepong, 52109, Malaysia. ¹²¹Institute for Atmosphere and Earth System Research/Physics, University of Helsinki, Helsinki, 00560, Finland. ¹²²Department of Botany, School of Natural Sciences, Trinity College Dublin, Dublin, D02PN40, Ireland. ¹²³German Meteorological Service (DWD), Centre for Agrometeorological Research, Braunschweig, 38116, Germany. ¹²⁴Micrometeorology, University of Bayreuth, Bayreuth, 95440, Germany. ¹²⁵Bayreuth Center of Ecology and Environmental Research, 95448, Bayreuth, Germany. ¹²⁶CSIRO Land and Water, Floreat, 6014, Australia. ¹²⁷Department of Agriculture, University of Sassari, Sassari, 07100, Italy. ¹²⁸Oulanka research station, University of Oulu, Kuusamo, 93900, Finland. ¹²⁹Dept. Biological Sciences, Wellesley College, Wellesley, MA, 02481, USA. ¹³⁰Research Institute on Terrestrial Ecosystems, National Research Council of Italy, Monterotondo Scalo, 00015, Italy. ¹³¹Department of Geography and Planning, Queen's University, Kingston, ON, K7L3N6, Canada. ¹³²Environmental Analytics NZ, Ltd. Raumati South, Paraparaumu, 5032, New Zealand. ¹³³Mazingira Centre, International Livestock Research Institute (ILRI), Nairobi, 00100, Kenya. ¹³⁴NOAA/OAR/Air Resources Laboratory, 325 Broadway, Boulder, CO, 80303, USA. ¹³⁵Atmospheric Sciences Research Center, State University of New York at Albany, Albany, NY, 12203, USA. ¹³⁶Forest Department of South Tyrol, Bolzano, 39100, Italy. ¹³⁷Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, 85721, USA. ¹³⁸Faculty of Science and Technology, Free University of Bolzano, Bolzano, 39100, Italy. ¹³⁹Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA. ¹⁴⁰IHE Delft, Delft, 2611, The Netherlands. ¹⁴¹Faculty of Science, VU Amsterdam, Amsterdam, 1081, The Netherlands. ¹⁴²University Grenoble Alpes, IRD, CNRS, IGE, Grenoble, 38000, France. ¹⁴³School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA. ¹⁴⁴Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA, 02138, USA. ¹⁴⁵School of Forestry and Resource Conservation, National Taiwan University, Taipei, 0617, Taiwan. ¹⁴⁶International Arctic Research Center, University of Alaska Fairbanks, Fairbanks, AK, 99775, USA. ¹⁴⁷Environment and Climate, Research Institute for Nature and Forest, Geraardsbergen, 9500, Belgium. ¹⁴⁸Department of Ecosystem Science and Management, Texas A&M University, College Station, TX, 77843, USA. ¹⁴⁹Research Institute for the Environment and Livelihoods, Charles Darwin University, Casuarina, 0810, Australia. ¹⁵⁰Grupo de Estudios Ambientales, Instituto de Matemática Aplicada San Luis (UNSL & CONICET), San Luis, D5700HHW, Argentina. ¹⁵¹Facultad de Ciencias Agropecuarias (UNER), Oro Verde, 3100, Argentina. ¹⁵²Eco&Sols, Univ Montpellier-CIRAD-INRA-IRD-Montpellier SupAgro, Montpellier, 34060, France. ¹⁵³O'Neill School of Public and Environmental Affairs, Indiana University Bloomington, Bloomington, IN, 47405, USA. ¹⁵⁴Global Change Research Group, Dept. Biology, San Diego State University, San Diego, CA, 92182, USA. ¹⁵⁵Department of Geography, College of Life and Environmental Sciences, University of Exeter, Exeter, EX44RJ, United Kingdom. ¹⁵⁶Department of Agroecology, Aarhus University, Tjele, 8830, Denmark. ¹⁵⁷iCLIMATE, Aarhus University, Tjele, 8830, Denmark. ¹⁵⁸Department of Geology, Wayne State University, Detroit, MI, 48202, USA. ¹⁵⁹Department of Geosciences, University of Oslo, Oslo, 0315, Norway. ¹⁶⁰Department of Geography, University of Zurich, Zurich, 8057, Switzerland. ¹⁶¹Department of Forest Ecology and Management, Swedish University of Agricultural Sciences, Umeå, 90183, Sweden. ¹⁶²Hawkesbury Institute for the Environment, Western Sydney University, Penrith, 2751, Australia. ¹⁶³Department of Biology, Indiana University Bloomington, Bloomington, IN, 47401, USA. ¹⁶⁴CSIRO Land and Water, Wembley, 6913, Australia. ¹⁶⁵Instituto de Clima y Agua, Instituto Nacional de Tecnología Agropecuaria (INTA), Buenos Aires, 1686, Argentina. ¹⁶⁶Department Computational Hydrosystems, Helmholtz Centre for Environmental Research UFZ, Leipzig, 04318, Germany. ¹⁶⁷Center for Global Change & Earth Observations, Michigan State University, East Lansing, MI, 48823, USA. ¹⁶⁸School of Life Science and Engineering, Southwest University of Science and Technology, Mianyang, 621010, China. ¹⁶⁹Departamento de Química e Física, Universidade Federal da Paraíba, Areia, PB, 58397-000, Brazil. ¹⁷⁰Remote Sensing and Geoinformatics, GFZ German Research Centre for Geosciences, Potsdam, 14473, Germany. ¹⁷¹Andalusian Institute for Earth System Research (CEAMA-IISTA), Granada, 18006, Spain. ¹⁷²Ciencias del Agua y Medioambiente, Instituto Tecnológico de Sonora, Ciudad Obregón, 85000, Mexico. ¹⁷³Geographical Institute, University of Cologne, Cologne, 50923, Germany. ¹⁷⁴Department of Industry, Innovation and Science, Geoscience Australia, Canberra, 2609, Australia. ¹⁷⁵Southwest Watershed Research Center, USDA-ARS, Tucson, AZ, 85719, USA. ¹⁷⁶Institute of Atmospheric Physics of the Czech Academy of Sciences, Prague, 14100, Czech Republic. ¹⁷⁷Department of Ecology, University of Granada, Granada, 18071, Spain. ¹⁷⁸National Hulunber Grassland Ecosystem Observation and Research Station & Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, 100081, China. ¹⁷⁹School of Science, Edith Cowan University, Joondalup, 6027, Australia. ¹⁸⁰Sentek Pty Ltd, Stepney, SA5069, Australia. ¹⁸¹National Ecological Observatory Network Program, Boulder, CO, 80301, USA. ¹⁸²Kansai Research Center, Forestry and Forest Products Research Institute, Kyoto, 612-0855, Japan. ¹⁸³College of Urban and Environmental Sciences, Peking University, Beijing, 100871, China. ¹⁸⁴School of Earth, Atmosphere and Environment, Monash University, Clayton, 3800, Australia. ¹⁸⁵Space Science and Engineering Center, University of Wisconsin-Madison, Madison, WI, 53706, USA. ¹⁸⁶Terrasystem srl, Viterbo, 01100, Italy. ¹⁸⁷USDA Forest Service, Rocky Mountain Research Station, Missoula, MT, 59808, USA. ¹⁸⁸Meteorology and Air Quality group, Wageningen University, 6500, Wageningen, The Netherlands. ¹⁸⁹Fenner School of Environment and Society, Australian National University Canberra, Canberra, ACT, 2600, Australia. ¹⁹⁰Earth and Life Institute, Université Catholique de Louvain, Louvain, 1348, Belgium. ¹⁹¹Department of Civil Engineering, Monash University, Clayton, 3800, Australia. ¹⁹²CSIRO Land and Water, Canberra, 2601, Australia. ¹⁹³School of Earth and Environmental Sciences, The University of Queensland, St Lucia, 4072, Australia. ¹⁹⁴South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, 510650, China. ¹⁹⁵College of Applied Meteorology, Nanjing University of Information Science & Technology, Nanjing, 210044, China. ¹⁹⁶Department of Animal and Plant Sciences, University of Sheffield, Sheffield, S102TN, United Kingdom. ¹⁹⁷Deceased: Ray Leuning.

e-mail: gzpastorello@lbl.gov; darpap@unitus.it