AUTHOR=Vidale P. L. , Egea G. , McGuire P. C. , Todt M. , Peters W. , Müller O. , Balan-Sarojini B. , Verhoef A. 

TITLE=On the Treatment of Soil Water Stress in GCM Simulations of Vegetation Physiology

JOURNAL=Frontiers in Environmental Science

VOLUME=9

YEAR=2021

URL=https://www.frontiersin.org/journals/environmental-science/articles/10.3389/fenvs.2021.689301

DOI=10.3389/fenvs.2021.689301

ISSN=2296-665X

ABSTRACT=<p>Current land surface schemes in weather and climate models make use of the so-called coupled photosynthesis–stomatal conductance (<italic>A–g</italic><sub><italic>s</italic></sub>) models of plant function to determine the surface fluxes that govern the terrestrial energy, water and carbon budgets. Plant physiology is controlled by many environmental factors, and a number of complex feedbacks are involved, but soil moisture control on root water uptake is primary, particularly in sub-tropical to temperate ecosystems. Land surface models represent plant water stress in different ways, but most implement a water stress factor, <italic>β</italic>, which ranges linearly (more recently also curvilinearly) between <italic>β</italic> = 1 for unstressed vegetation and <italic>β</italic> = 0 at the wilting point, expressed in terms of volumetric water content (<inline-formula id="inf1"><mml:math id="m1" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>θ</mml:mi></mml:math></inline-formula>). <inline-formula id="inf2"><mml:math id="m2" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mo> </mml:mo><mml:mi>β</mml:mi></mml:mrow></mml:math></inline-formula> is most commonly used to either limit <italic>A</italic> or <italic>g</italic><sub>s</sub>, and hence carbon and water fluxes, and a pertinent research question is whether these treatments are in fact interchangeable. Following Egea et al. (Agricultural and Forest Meteorology, 2011, 151 (10), 1,370–1,384) and Verhoef et al. (Agricultural and Forest Meteorology, 2014, 191, 22–32), we have implemented new <italic>β</italic> treatments, reflecting higher levels of biophysical complexity in a state-of-the-art LSM, Joint UK Land Environment Simulator, by allowing root zone soil moisture to limit plant function non-linearly and via individual routes (carbon assimilation, stomatal conductance, or mesophyll conductance) as well as any (non-linear) combinations thereof. The treatment of <italic>β</italic> does matter to the prediction of water and carbon fluxes: this study demonstrates that it represents a key structural uncertainty in contemporary LSMs, in terms of predictions of gross primary productivity, energy fluxes and soil moisture evolution, both in terms of climate means and response to a number of European droughts, including the 2003 heat wave. Treatments allowing <italic>ß</italic> to act on vegetation fluxes via stomatal and mesophyll routes are able to simulate the spatiotemporal variability in water use efficiency with higher fidelity during the growing season; they also support a broader range of ecosystem responses, e.g., those observed in regions that are radiation limited or water limited. We conclude that current practice in weather and climate modelling is inconsistent, as well as too simplistic, failing to credibly simulate vegetation response to soil water stress across the typical range of variability that is encountered for current European weather and climate conditions, including extremes of land surface temperature and soil moisture drought. A generalized approach performs better in current climate conditions and promises to be, based on responses to recently observed extremes, more trustworthy for predicting the impacts of climate change.</p>