A data-driven analysis of energy balance closure across FLUXNET research sites: The role of landscape scale heterogeneity

The energy balance at most surface-atmosphere flux research sites remains unclosed. The mechanisms underlying the discrepancy between measured energy inputs and outputs across the global FLUXNET tower network are still under debate. Recent reviews have identified exchange processes and turbulent motions at large spatial and temporal scales in heterogeneous landscapes as the primary cause of the lack of energy balance closure at some intensively-researched sites, while unmeasured storage terms cannot be ruled out as a dominant contributor to the lack of energy balance closure at many other sites. We analyzed energy balance closure across 173 ecosystems in the FLUXNET database and explored the relationship between energy balance closure and landscape heterogeneity using MODIS products and GLOBEstat elevation data. Energy balance closure per research site (C_EB,s) averaged 0.84±0.20, with best average closures in evergreen broadleaf forests and savannas (0.91-0.94) and worst average closures in crops, deciduous broadleaf forests, mixed forests and wetlands (0.70-0.78). Half-hourly or hourly energy balance closure on a percent basis increased with friction velocity (u_*) and was highest on average under near-neutral atmospheric conditions. C_EB,s was significantly related to mean precipitation, gross primary productivity and landscape-level enhanced vegetation index (EVI) from MODIS, and the variability in elevation, MODIS plant functional type, and MODIS EVI. A linear model including landscape-level variability in both EVI and elevation, mean precipitation, and an interaction term between EVI variability and precipitation had the lowest Akaike's information criterion value. C_EB,s in landscapes with uniform plant functional type approached 0.9 and C_EB,s in landscapes with uniform EVI approached 1. These results suggest that landscape-level heterogeneity in vegetation and topography cannot be ignored as a contributor to incomplete energy balance closure at the flux network level, although net radiation measurements, biological energy assimilation, unmeasured storage terms, and the importance of good practice including site selection when making flux measurements should not be discounted. Our results suggest that future research should focus on the quantitative mechanistic relationships between energy balance closure and landscape-scale heterogeneity, and the consequences of mesoscale circulations for surface-atmosphere exchange measurements.