Arctic tundra has the potential to generate significant climate feedbacks, but spatial complexity makes it difficult to quantify the impacts of climate on ecosystem-atmosphere fluxes, particularly in polygonal tundra comprising wetter and drier polygon types on the scale of tens of meters. We measured CO2, CH4, and energy fluxes using eddy covariance for 7 yr (April to November, 2013–2019) in polygonal tundra near Utqiagvik, Alaska. This period saw the earliest snowmelt, latest snow accumulation, and hottest summer on record. To estimate fluxes by polygon type, we combined a polygon classification with a flux-footprint model. Methane fluxes were highest in the summer months but were also large during freeze-up and increased with the warming trend in August–November temperatures. While CO2 respiration had a consistent, exponential relationship with temperature, net ecosystem exchange was more variable among years. CO2 and CH4 exchange (June–September) ranged between −0.83 (Standard error [SE] = 0.03) and −1.32 (SE = 0.04) μmol m−2 s−1 and 13.92 (SE = 0.26)—23.42 (SE = 0.45) nmol m−2 s−1, respectively, and varied interannually (p ≤ 0.05). The maximum-influence method effectively attributed fluxes to polygon types. Areas dominated by low-centered polygons had higher CO2 fluxes except in 2016–2017. Methane fluxes were highest in low-centered polygons 2013–2015 and in flat-centered polygons in subsequent years, possibly due to increasing temperature and precipitation. Sensible and latent heat fluxes also varied significantly among polygon types. Accurate characterization of Arctic fluxes and their climate dependencies requires spatial disaggregation and long term observations.

Net photosynthesis is the largest single flux in the global carbon cycle, but controls over its variability are poorly understood because there is no direct way of measuring it at the ecosystem scale. We report observations of ecosystem carbonyl sulfide (COS) and CO2 fluxes that resolve key gaps in an emerging framework for using concurrent COS and CO2 measurements to quantify terrestrial gross primary productivity. At a wheat field in Oklahoma we found that in the peak growing season the flux-weighted leaf relative uptake of COS and CO2 during photosynthesis was 1.3, at the lower end of values from laboratory studies, and varied systematically with light. Due to nocturnal stomatal conductance, COS uptake by vegetation continued at night, contributing a large fraction (29%) of daily net ecosystem COS fluxes. In comparison, the contribution of soil fluxes was small (1-6%) during the peak growing season. Upland soils are usually considered sinks of COS. In contrast, the well-aerated soil at the site switched from COS uptake to emissions at a soil temperature of around 15 °C. We observed COS production from the roots of wheat and other species and COS uptake by root-free soil up to a soil temperature of around 25 °C. Our dataset demonstrates that vegetation uptake is the dominant ecosystem COS flux in the peak growing season, providing support of COS as an independent tracer of terrestrial photosynthesis. However, the observation that ecosystems may become a COS source at high temperature needs to be considered in global modeling studies.

The Southern Great Plains are characterized by a fine-scale mixture of different land-cover types, predominantly winter-wheat and grazed pasture, with relatively small areas of other crops, native prairie, and switchgrass. Recent droughts and predictions of increased drought in the Southern Great Plains, especially during the summer months, raise concern for these ecosystems. We measured ecosystem carbon and water fluxes with eddy-covariance systems over cultivated cropland for 10 years, and over lightly grazed prairie and new switchgrass fields for 2 years each. Growing-season precipitation showed the strongest control over net carbon uptake for all ecosystems, but with a variable effect: grasses (prairie and switchgrass) needed at least 350mm of precipitation during the growing season to become net carbon sinks, while crops needed only 100mm. In summer, high temperatures enhanced evaporation and led to higher likelihood of dry soil conditions. Therefore, summer-growing native prairie species and switchgrass experienced more seasonal droughts than spring-growing crops. For wheat, the net reduction in carbon uptake resulted mostly from a decrease in gross primary production rather than an increase in respiration. Flux measurements suggested that management practices for crops were effective in suppressing evapotranspiration and decomposition (by harvesting and removing secondary growth), and in increasing carbon uptake (by fertilizing and conserving summer soil water). In light of future projections for wetter springs and drier and warmer summers in the Southern Great Plains, our study indicates an increased vulnerability in native ecosystems and summer crops over time.

Land-atmosphere interactions are important to climate prediction, but the underlying effects of surface forcing of the atmosphere are not well understood. In the U.S. Southern Great Plains, grassland/pasture and winter wheat are the dominant land covers but have distinct growing periods that may differently influence land-atmosphere coupling during spring and summer. Variables that influence surface flux partitioning can change seasonally, depending on the state of local vegetation. Here we use surface observations from multiple sites in the U.S. Department of Energy Atmospheric Radiation Measurement Southern Great Plains Climate Research Facility and statistical modeling at a paired grassland/agricultural site within this facility to quantify land cover influence on surface energy balance and variables controlling evaporative fraction (latent heat flux normalized by the sum of sensible and latent heat fluxes). We demonstrate that the radiative balance and evaporative fraction are closely related to green leaf area at both winter wheat and grassland/pasture sites and that the early summer harvest of winter wheat abruptly shifts the relationship between evaporative fraction and surface state variables. Prior to harvest, evaporative fraction of winter wheat is strongly influenced by leaf area and soil-atmosphere temperature differences. After harvest, variations in soil moisture have a stronger effect on evaporative fraction. This is in contrast with grassland/pasture sites, where variation in green leaf area has a large influence on evaporative fraction throughout spring and summer, and changes in soil-atmosphere temperature difference and soil moisture are of relatively minor importance.

The few prethaw observations of tundra carbon fluxes suggest that there may be large spring releases, but little is known about the scale and underlying mechanisms of this phenomenon. To address these questions, we combined ecosystem eddy flux measurements from two towers near Barrow, Alaska, with mechanistic soil-core thawing experiment. During a 2 week period prior to snowmelt in 2014, large fluxes were measured, reducing net summer uptake of CO2 by 46% and adding 6% to cumulative CH4 emissions. Emission pulses were linked to unique rain-on-snow events enhancing soil cracking. Controlled laboratory experiment revealed that as surface ice thaws, an immediate, large pulse of trapped gases is emitted. These results suggest that the Arctic CO2 and CH4 spring pulse is a delayed release of biogenic gas production from the previous fall and that the pulse can be large enough to offset a significant fraction of the moderate Arctic tundra carbon sink.

Wetlands are the largest global natural methane (CH4/ source, and emissions between 50 and 70° N latitude contribute 10-30% to this source. Predictive capability of land models for northern wetland CH4 emissions is still low due to limited site measurements, strong spatial and temporal variability in emissions, and complex hydrological and biogeochemical dynamics. To explore this issue, we compare wetland CH4 emission predictions from the Community Land Model 4.5 (CLM4.5-BGC) with siteto regional-scale observations. A comparison of the CH4 fluxes with eddy flux data highlighted needed changes to the model's estimate of aerenchyma area, which we implemented and tested. The model modification substantially reduced biases in CH4 emissions when compared with CarbonTracker CH4 predictions. CLM4.5 CH4 emission predictions agree well with growing season (May-September) CarbonTracker Alaskan regional-level CH4 predictions and sitelevel observations. However, CLM4.5 underestimated CH4 emissions in the cold season (October-April). The monthly atmospheric CH4 mole fraction enhancements due to wetland emissions are also assessed using the Weather Research and Forecasting-Stochastic Time-Inverted Lagrangian Transport (WRF-STILT) model coupled with daily emissions from CLM4.5 and compared with aircraft CH4 mole fraction measurements from the Carbon in Arctic Reservoirs Vulnerability Experiment (CARVE) campaign. Both the tower and aircraft analyses confirm the underestimate of cold-season CH4 emissions by CLM4.5. The greatest uncertainties in predicting the seasonal CH4 cycle are from the wetland extent, coldseason CH4 production and CH4 transport processes. We recommend more cold-season experimental studies in highlatitude systems, which could improve the understanding and parameterization of ecosystem structure and function during this period. Predicted CH4 emissions remain uncertain, but we show here that benchmarking against observations across spatial scales can inform model structural and parameter improvements.

While the eddy covariance (EC) technique is a well-established method for measuring water fluxes (i.e., evaporation or 'evapotranspiration’, ET), the measurement is susceptible to many uncertainties. One such issue is the potential underestimation of ET when relative humidity (RH) is high (>70%), due to low-pass filtering with some EC systems. Yet, this underestimation for different types of EC systems (e.g. open-path or closed-path sensors) has not been characterized for synthesis datasets such as the widely used FLUXNET2015 dataset. Here, we assess the RH-associated underestimation of latent heat fluxes (LE, or ET) from different EC systems for 163 sites in the FLUXNET2015 dataset. We found that the LE underestimation is most apparent during hours when RH is higher than 70%, predominantly observed at sites using closed-path EC systems, but the extent of the LE underestimation is highly site-specific. We then propose a machine learning based method to correct for this underestimation, and compare it to two energy balance closure based LE correction approaches (Bowen ratio correction, BRC, and attributing all errors to LE). Our correction increases LE by 189% for closed-path sites at high RH (>90%), while BRC increases LE by around 30% for all RH conditions. Additionally, we assess the influence of these corrections on ET-based transpiration (T) estimates using two different ET partitioning methods. Results show opposite responses (increasing vs. slightly decreasing T-to-ET ratios, T/ET) between the two methods when comparing T based on corrected and uncorrected LE. Overall, our results demonstrate the existence of a high RH bias in water fluxes in the FLUXNET2015 dataset and suggest that this bias is a pronounced source of uncertainty in ET measurements to be considered when estimating ecosystem T/ET and WUE.

Large datasets of greenhouse gas and energy surface-atmosphere fluxes measured with the eddy-covariance technique (e.g., FLUXNET2015, AmeriFlux BASE) are widely used to benchmark models and remote-sensing products. This study addresses one of the major challenges facing model-data integration: To what spatial extent do flux measurements taken at individual eddy-covariance sites reflect model- or satellite-based grid cells? We evaluate flux footprints—the temporally dynamic source areas that contribute to measured fluxes—and the representativeness of these footprints for target areas (e.g., within 250–3000 m radii around flux towers) that are often used in flux-data synthesis and modeling studies. We examine the land-cover composition and vegetation characteristics, represented here by the Enhanced Vegetation Index (EVI), in the flux footprints and target areas across 214 AmeriFlux sites, and evaluate potential biases as a consequence of the footprint-to-target-area mismatch. Monthly 80% footprint climatologies vary across sites and through time ranging four orders of magnitude from 103 to 107 m2 due to the measurement heights, underlying vegetation- and ground-surface characteristics, wind directions, and turbulent state of the atmosphere. Few eddy-covariance sites are located in a truly homogeneous landscape. Thus, the common model-data integration approaches that use a fixed-extent target area across sites introduce biases on the order of 4%–20% for EVI and 6%–20% for the dominant land cover percentage. These biases are site-specific functions of measurement heights, target area extents, and land-surface characteristics. We advocate that flux datasets need to be used with footprint awareness, especially in research and applications that benchmark against models and data products with explicit spatial information. We propose a simple representativeness index based on our evaluations that can be used as a guide to identify site-periods suitable for specific applications and to provide general guidance for data use.

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.

The following authors were omitted from the original version of this Data Descriptor: Markus Reichstein and Nicolas Vuichard. Both contributed to the code development and N. Vuichard contributed to the processing of the ERA-Interim data downscaling. Furthermore, the contribution of the co-author Frank Tiedemann was re-evaluated relative to the colleague Corinna Rebmann, both working at the same sites, and based on this re-evaluation a substitution in the co-author list is implemented (with Rebmann replacing Tiedemann). Finally, two affiliations were listed incorrectly and are corrected here (entries 190 and 193). The author list and affiliations have been amended to address these omissions in both the HTML and PDF versions.