Abstract: Terrestrial evapotranspiration is the second‐largest component of the land water cycle, linking the water, energy, and carbon cycles and influencing the productivity and health of ecosystems. The dynamics of ET across a spectrum of spatiotemporal scales and their controls remain an active focus of research across different science disciplines. Here, we provide an overview of the current state of ET science across in situ measurements, partitioning of ET, and remote sensing, and discuss how different approaches complement one another based on their advantages and shortcomings. We aim to facilitate collaboration among a cross‐disciplinary group of ET scientists to overcome the challenges identified in this paper and ultimately advance our integrated understanding of ET.

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.

Fundamental axes of variation in plant traits result from trade-offs between costs and benefits of resource-use strategies at the leaf scale. However, it is unclear whether similar trade-offs propagate to the ecosystem level. Here, we test whether trait correlation patterns predicted by three well-known leaf- and plant-level coordination theories - the leaf economics spectrum, the global spectrum of plant form and function, and the least-cost hypothesis - are also observed between community mean traits and ecosystem processes. We combined ecosystem functional properties from FLUXNET sites, vegetation properties, and community mean plant traits into three corresponding principal component analyses. We find that the leaf economics spectrum (90 sites), the global spectrum of plant form and function (89 sites), and the least-cost hypothesis (82 sites) all propagate at the ecosystem level. However, we also find evidence of additional scale-emergent properties. Evaluating the coordination of ecosystem functional properties may aid the development of more realistic global dynamic vegetation models with critical empirical data, reducing the uncertainty of climate change projections.

The leaf economics spectrum^1,2 and the global spectrum of plant forms and functions³ revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species². Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities⁴. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability^4,5. Here we derive a set of ecosystem functions⁶ from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems^7,8.