loading page

Hybrid Modeling of Evapotranspiration: Inferring Stomatal and Aerodynamic Resistances Using Combined Physics-Based and Machine Learning
  • +4
  • Reda ElGhawi,
  • Basil Kraft,
  • Christian Reimers,
  • Markus Reichstein,
  • Marco KΓΆrner,
  • Pierre Gentine,
  • Alexander J Winkler
Reda ElGhawi
Max Planck Institute for Biogeochemistry, Max Planck Institute for Biogeochemistry

Corresponding Author:[email protected]

Author Profile
Basil Kraft
Max Planck Institute for Biogeochemistry, Max Planck Institute for Biogeochemistry
Author Profile
Christian Reimers
Max Planck Institute for Biogeochemistry, Max Planck Institute for Biogeochemistry
Author Profile
Markus Reichstein
Max Planck Institute for Biogeochemistry, Max Planck Institute for Biogeochemistry
Author Profile
Marco KΓΆrner
Technical University of Munich, Technical University of Munich
Author Profile
Pierre Gentine
Columbia University, Columbia University
Author Profile
Alexander J Winkler
Max Planck Institute for Biogeochemistry, Max Planck Institute for Biogeochemistry
Author Profile

Abstract

The process of evapotranspiration transfers water vapour from vegetation and soil surfaces to the atmosphere, the so-called latent heat flux (𝑄 LE), and thus crucially modulates Earth’s energy, water, and carbon cycles. Vegetation controls 𝑄 LE through regulating the leaf stomata (i.e., surface resistance π‘Ÿ s) and through altering surface roughness (aerodynamic resistance π‘Ÿ a). Estimating π‘Ÿ s and π‘Ÿ a across different vegetation types proves to be a key challenge in predicting 𝑄 LE. Here, we propose a hybrid modeling approach (i.e., combining mechanistic modeling and machine learning) for 𝑄 LE where neural networks independently learn the resistances from observations as intermediate variables. In our hybrid modeling setup, we make use of the Penman-Monteith equation based on the Big Leaf theory in conjunction with multi-year flux measurements across different forest and grassland sites from the FLUXNET database. We follow two conceptually different strategies to constrain the hybrid model to control for equifinality arising when estimating the two resistances simultaneously. One strategy is to impose an a priori constraint on π‘Ÿ a based on our mechanistic understanding (theory-driven strategy), while the other strategy makes use of more observational data and adds a constraint in predicting π‘Ÿ a through multi-task learning of the latent as well as the sensible heat flux (𝑄 H ; data-driven strategy). Our results show that all hybrid models exhibit a fairly high predictive skill for the target variables with 𝑅 2 = 0.82-0.89 for grasslands and 𝑅 2 = 0.70-0.80 for forests sites at the mean diurnal scale. The predictions of π‘Ÿ s and π‘Ÿ a show physical consistency across the two regularized hybrid models, but are physically implausible in the under-constrained hybrid model. The hybrid models are robust in reproducing consistent results for energy fluxes and resistances across different scales (diurnal, seasonal, interannual), reflecting their ability to learn the physical dependence of the target variables on the meteorological inputs. As a next step, we propose to test these heavily observation-informed parameterizations derived through hybrid modeling as a substitute for overly simple ad hoc formulations in Earth system models.