Premium
Comparing optimal and empirical stomatal conductance models for application in Earth system models
Author(s) -
Franks Peter J.,
Bonan Gordon B.,
Berry Joseph A.,
Lombardozzi Danica L.,
Holbrook N. Michele,
Herold Nicholas,
Oleson Keith W.
Publication year - 2018
Publication title -
global change biology
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 4.146
H-Index - 255
eISSN - 1365-2486
pISSN - 1354-1013
DOI - 10.1111/gcb.14445
Subject(s) - stomatal conductance , earth system science , environmental science , earth (classical element) , empirical modelling , computer science , geology , mathematics , oceanography , biology , simulation , photosynthesis , botany , mathematical physics
Earth system models ( ESM s) rely on the calculation of canopy conductance in land surface models ( LSM s) to quantify the partitioning of land surface energy, water, and CO 2 fluxes. This is achieved by scaling stomatal conductance, g w , determined from physiological models developed for leaves. Traditionally, models for g w have been semi‐empirical, combining physiological functions with empirically determined calibration constants. More recently, optimization theory has been applied to model g w in LSM s under the premise that it has a stronger grounding in physiological theory and might ultimately lead to improved predictive accuracy. However, this premise has not been thoroughly tested. Using original field data from contrasting forest systems, we compare a widely used empirical type and a more recently developed optimization‐type g w model, termed BB and MED , respectively. Overall, we find no difference between the two models when used to simulate g w from photosynthesis data, or leaf gas exchange from a coupled photosynthesis‐conductance model, or gross primary productivity and evapotranspiration for a FLUXNET tower site with the CLM 5 community LSM . Field measurements reveal that the key fitted parameters for BB and MED , g 1B and g 1M, exhibit strong species specificity in magnitude and sensitivity to CO 2 , and CLM 5 simulations reveal that failure to include this sensitivity can result in significant overestimates of evapotranspiration for high‐ CO 2 scenarios. Further, we show that g 1B and g 1M can be determined from mean c i / c a (ratio of leaf intercellular to ambient CO 2 concentration). Applying this relationship with c i / c a values derived from a leaf δ 13 C database, we obtain a global distribution of g 1B and g 1M , and these values correlate significantly with mean annual precipitation. This provides a new methodology for global parameterization of the BB and MED models in LSM s, tied directly to leaf physiology but unconstrained by spatial boundaries separating designated biomes or plant functional types.