The National Ecological Observatory Network (NEON) is a multidecadal and continental-scale observatory with sites across the United States. Having entered its operational phase in 2018, NEON data products, software, and services become available to facilitate research on the impacts of climate change, land-use change, and invasive species. An essential component of NEON are its 47 tower sites, where eddy-covariance (EC) sensors are operated to determine the surface–atmosphere exchange of momentum, heat, water, and CO₂. EC tower networks such as AmeriFlux, the Integrated Carbon Observation System (ICOS), and NEON are vital for providing the distributed observations to address interactions at the soil–vegetation–atmosphere interface. NEON represents the largest single-provider EC network globally, with standardized observations and data processing explicitly designed for intersite comparability and analysis of feedbacks across multiple spatial and temporal scales. Furthermore, EC is tightly integrated with soil, meteorology, atmospheric chemistry, isotope, phenology, and rich contextual observations such as airborne remote sensing and in situ sampling bouts. Here, we present an overview of NEON’s observational design, field operation, and data processing that yield community resources for the study of surface–atmosphere interactions. Near-real-time data products become available from the NEON Data Portal, and EC and meteorological data are ingested into AmeriFlux and FLUXNET globally harmonized data releases. Open-source software for reproducible, extensible, and portable data analysis includes the eddy4R family of R packages underlying the EC data product generation. These resources strive to integrate with existing infrastructures and networks, to suggest novel systemic solutions, and to synergize ongoing research efforts across science communities.

Large differences in instrumentation, site setup, data format, and operating system stymie the adoption of a universal computational environment for processing and analyzing eddy-covariance (EC) data. This results in limited software applicability and extensibility in addition to often substantial inconsistencies in flux estimates. Addressing these concerns, this paper presents the systematic development of portable, reproducible, and extensible EC software achieved by adopting a development and systems operation (DevOps) approach. This software development model is used for the creation of the eddy4R family of EC code packages in the open-source R language for statistical computing. These packages are community developed, iterated via the Git distributed version control system, and wrapped into a portable and reproducible Docker filesystem that is independent of the underlying host operating system. The HDF5 hierarchical data format then provides a streamlined mechanism for highly compressed and fully self-documented data ingest and output. The usefulness of the DevOps approach was evaluated for three test applications. First, the resultant EC processing software was used to analyze standard flux tower data from the first EC instruments installed at a National Ecological Observatory (NEON) field site. Second, through an aircraft test application, we demonstrate the modular extensibility of eddy4R to analyze EC data from other platforms. Third, an intercomparison with commercial-grade software showed excellent agreement (R² = 1.0 for CO₂ flux). In conjunction with this study, a Docker image containing the first two eddy4R packages and an executable example workflow, as well as first NEON EC data products are released publicly. We conclude by describing the work remaining to arrive at the automated generation of science-grade EC fluxes and benefits to the science community at large. This software development model is applicable beyond EC and more generally builds the capacity to deploy complex algorithms developed by scientists in an efficient and scalable manner. In addition, modularity permits meeting project milestones while retaining extensibility with time.

Continuous measurements of net ecosystem CO₂ exchange (NEE) using the eddy-covariance method were made over an agricultural ecosystem in the southeastern US. During optimum environmental conditions, photosynthetically active radiation (PAR) was the primary driver controlling daytime NEE, accounting for as much as 67 to 89% of the variation in NEE. However, soil water content became the dominant factor limiting the NEE-PAR response during the peak growth stage. NEE was significantly depressed when high PAR values coincided with very low soil water content. The presence of a counter-clockwise hysteresis of daytime NEE with PAR was observed during periods of water stress. This is a result of the stomatal closure control of photosynthesis at high vapor pressure deficit and enhanced respiration at high temperature. This result is significant since this hysteresis effect limits the range of applicability of the Michaelis-Menten equation and other related expressions in the determination of daytime NEE as a function of PAR. The systematic presence of hysteresis in the response of NEE to PAR suggests that the gap-filling technique based on a non-linear regression approach should take into account the presence of water-limited field conditions. Including this step is therefore likely to improve current evaluation of ecosystem response to increased precipitation variability arising from climatic changes.

Flux measurements of nitrogen oxides (NO_x) were made over London using airborne eddy covariance from a low-flying aircraft. Seven low-altitude flights were conducted over Greater London, performing multiple overpasses across the city during eight days in July 2014. NO_x fluxes across the Greater London region (GLR) exhibited high heterogeneity and strong diurnal variability, with central areas responsible for the highest emission rates (20–30 mg m⁻² h⁻¹). Other high-emission areas included the M25 orbital motorway. The complexity of London's emission characteristics makes it challenging to pinpoint single emissions sources definitively using airborne measurements. Multiple sources, including road transport and residential, commercial and industrial combustion sources, are all likely to contribute to measured fluxes. Measured flux estimates were compared to scaled National Atmospheric Emissions Inventory (NAEI) estimates, accounting for monthly, daily and hourly variability. Significant differences were found between the flux-driven emissions and the NAEI estimates across Greater London, with measured values up to 2 times higher in Central London than those predicted by the inventory. To overcome the limitations of using the national inventory to contextualise measured fluxes, we used physics-guided flux data fusion to train environmental response functions (ERFs) between measured flux and environmental drivers (meteorological and surface). The aim was to generate time-of-day emission surfaces using calculated ERF relationships for the entire GLR; 98 % spatial coverage was achieved across the GLR at 400 m² spatial resolution. All flight leg projections showed substantial heterogeneity across the domain, with high emissions emanating from Central London and major road infrastructure. The diurnal emission structure of the GLR was also investigated, through ERF, with the morning rush hour distinguished from lower emissions during the early afternoon. Overall, the integration of airborne fluxes with an ERF-driven strategy enabled the first independent generation of surface NO_x emissions, at high resolution using an eddy-covariance approach, for an entire city region.