A theoretical and empirical assessment of stomatal optimization modeling

Global convergence in terrestrial gross primary production response to atmospheric vapor pressure deficit

Atmospheric vapor pressure deficit (VPD) increases with climate warming and may limit plant growth. However, gross primary production (GPP) responses to VPD remain a mystery, offering a significant source of uncertainty in the estimation of global terrestrial ecosystems carbon dynamics. In this study, in-situ measurements, satellite-derived data, and Earth System Models (ESMs) simulations were analysed to show that the GPP of most ecosystems has a similar threshold in response to VPD: first increasing and then declining. When VPD exceeds these thresholds, atmospheric drought stress reduces soil moisture and stomatal conductance, thereby decreasing the productivity of terrestrial ecosystems. Current ESMs underscore CO2 fertilization effects but predict significant GPP decline in low-latitude ecosystems when VPD exceeds the thresholds. These results emphasize the impacts of climate warming on VPD and propose limitations to future ecosystems productivity caused by increased atmospheric water demand. Incorporating VPD, soil moisture, and canopy conductance interactions into ESMs enhances the prediction of terrestrial ecosystem responses to climate change.

Incorporating photosynthetic acclimation improves stomatal optimisation models

Stomatal opening in plant leaves is regulated through a balance of carbon and water exchange under different environmental conditions. Accurate estimation of stomatal regulation is crucial for understanding how plants respond to changing environmental conditions, particularly under climate change. A new generation of optimality-based modelling schemes determines instantaneous stomatal responses from a balance of trade-offs between carbon gains and hydraulic costs, but most such schemes do not account for biochemical acclimation in response to drought. Here, we compare the performance of six instantaneous stomatal optimisation models with and without accounting for photosynthetic acclimation. Using experimental data from 37 plant species, we found that accounting for photosynthetic acclimation improves the prediction of carbon assimilation in a majority of the tested models. Photosynthetic acclimation contributed significantly to the reduction of photosynthesis under drought conditions in all tested models. Drought effects on photosynthesis could not accurately be explained by the hydraulic impairment functions embedded in the stomatal models alone, indicating that photosynthetic acclimation must be considered to improve estimates of carbon assimilation during drought.

The poorly-explored stomatal response to temperature at constant evaporative demand

Changes in leaf temperature are known to drive stomatal responses, because the leaf-to-air water vapour gradient (Δw) increases with temperature if ambient vapour pressure is held constant, and stomata respond to changes in Δw. However, the direct response of stomata to temperature (DRST; the response when Δw is held constant by adjusting ambient humidity) has been examined far less extensively. Though the meagre available data suggest the response is usually positive, results differ widely and defy broad generalisation. As a result, little is known about the DRST. This review discusses the current state of knowledge about the DRST, including numerous hypothesised biophysical mechanisms, potential implications of the response for plant adaptation, and possible impacts of the DRST on plant-atmosphere carbon and water exchange in a changing climate.

Different roads, same destination: The shared future of plant ecophysiology and ecohydrology

Terrestrial water fluxes are substantially mediated by vegetation, while the distribution, growth, health, and mortality of plants are strongly influenced by the availability of water. These interactions, playing out across multiple spatial and temporal scales, link the disciplines of plant ecophysiology and ecohydrology. Despite this connection, the disciplines have provided complementary, but largely independent, perspectives on the soil-plant-atmosphere continuum since their crystallization as modern scientific disciplines in the late 20th century. This review traces the development of the two disciplines, from their respective origins in engineering and ecology, their largely independent growth and maturation, and the eventual development of common conceptual and quantitative frameworks. This common ground has allowed explicit coupling of the disciplines to better understand plant function. Case studies both illuminate the limitations of the disciplines working in isolation, and reveal the exciting possibilities created by consilience between the disciplines. The histories of the two disciplines suggest opportunities for new advances will arise from sharing methodologies, working across multiple levels of complexity, and leveraging new observational technologies. Practically, these exchanges can be supported by creating shared scientific spaces. This review argues that consilience and collaboration are essential for robust and evidence-based predictions and policy responses under global change.

Sap flow through partially embolized xylem vessel networks

Sap is transported through numerous conduits in the xylem of woody plants along the path from the soil to the leaves. When all conduits are functional, vessel lumen diameter is a strong predictor of hydraulic conductivity. As vessels become embolized, sap movement becomes increasingly affected by factors operating at scales beyond individual conduits, creating resistances that result in hydraulic conductivity diverging from diameter-based estimates. These effects include pit resistances, connectivity, path length, network topology, and vessel or sector isolation. The impact of these factors varies with the level and distribution of emboli within the network, and manifest as alterations in the relationship between the number and diameter of embolized vessels with measured declines in hydraulic conductivity across vulnerability to embolism curves. Divergences between measured conductivity and diameter-based estimates reveal functional differences that arise because of species- and tissue-specific vessel network structures. Such divergences are not uniform, and xylem tissues may diverge in different ways and to differing degrees. Plants regularly operate under nonoptimal conditions and contain numerous embolized conduits. Understanding the hydraulic implications of emboli within a network and the function of partially embolized networks are critical gaps in our understanding of plants occurring within natural environments.

The multidimensionality of plant drought stress: The relative importance of edaphic and atmospheric drought

Drought threatens plant growth and related ecosystem services. The emergence of plant drought stress under edaphic drought is well studied, whilst the importance of atmospheric drought only recently gained momentum. Yet, little is known about the interaction and relative contribution of edaphic and atmospheric drought on the emergence of plant drought stress. We conducted a gradient experiment, fully crossing gravimetric water content (GWC: maximum water holding capacity-permanent wilting point) and vapour pressure deficit (VPD: 1-2.25 kPa) using five wheat varieties from three species (Triticum monococcum, T. durum & T. aestivum). We quantified the occurrence of plant drought stress on molecular (abscisic acid), cellular (stomatal conductance), organ (leaf water potential) and stand level (evapotranspiration). Plant drought stress increased with decreasing GWC across all organizational levels. This effect was magnified nonlinearly by VPD after passing a critical threshold of soil water availability. At around 20%GWC (soil matric potential 0.012 MPa), plants lost their ability to regulate leaf water potential via stomata regulation, followed by the emergence of hydraulic dysfunction. The emergence of plant drought stress is characterized by changing relative contributions of soil versus atmosphere and their non-linear interaction. This highly non-linear response is likely to abruptly alter plant-related ecosystem services in a drying world.

Preserving isohydricity: vertical environmental variability explains Amazon forest water-use strategies

Increases in hydrological extremes, including drought, are expected for Amazon forests. A fundamental challenge for predicting forest responses lies in identifying ecological strategies which underlie such responses. Characterization of species-specific hydraulic strategies for regulating water-use, thought to be arrayed along an 'isohydric-anisohydric' spectrum, is a widely used approach. However, recent studies have questioned the usefulness of this classification scheme, because its metrics are strongly influenced by environments, and hence can lead to divergent classifications even within the same species. Here, we propose an alternative approach positing that individual hydraulic regulation strategies emerge from the interaction of environments with traits. Specifically, we hypothesize that the vertical forest profile represents a key gradient in drought-related environments (atmospheric vapor pressure deficit, soil water availability) that drives divergent tree water-use strategies for coordinated regulation of stomatal conductance (gs) and leaf water potentials (ΨL) with tree rooting depth, a proxy for water availability. Testing this hypothesis in a seasonal eastern Amazon forest in Brazil, we found that hydraulic strategies indeed depend on height-associated environments. Upper canopy trees, experiencing high vapor pressure deficit (VPD), but stable soil water access through deep rooting, exhibited isohydric strategies, defined by little seasonal change in the diurnal pattern of gs and steady seasonal minimum ΨL. In contrast, understory trees, exposed to less variable VPD but highly variable soil water availability, exhibited anisohydric strategies, with fluctuations in diurnal gs that increased in the dry season along with increasing variation in ΨL. Our finding that canopy height structures the coordination between drought-related environmental stressors and hydraulic traits provides a basis for preserving the applicability of the isohydric-to-anisohydric spectrum, which we show here may consistently emerge from environmental context. Our work highlights the importance of understanding how environmental heterogeneity structures forest responses to climate change, providing a mechanistic basis for improving models of tropical ecosystems.

Theory and tests for coordination among hydraulic and photosynthetic traits in co-occurring woody species

Co-occurring plants show wide variation in their hydraulic and photosynthetic traits. Here, we extended 'least-cost' optimality theory to derive predictions for how variation in key hydraulic traits potentially affects the cost of acquiring and using water in photosynthesis and how this, in turn, should drive variation in photosynthetic traits. We tested these ideas across 18 woody species at a temperate woodland in eastern Australia, focusing on hydraulic traits representing different aspects of plant water balance, that is storage (sapwood capacitance, CS), demand vs supply (branch leaf : sapwood area ratio, AL : AS and leaf : sapwood mass ratio and ML : MS), access to soil water (proxied by predawn leaf water potential, ΨPD) and physical strength (sapwood density, WD). Species with higher AL : AS had higher ratio of leaf-internal to ambient CO2 concentration during photosynthesis (ci : ca), a trait central to the least-cost theory framework. CS and the daily operating range of tissue water potential (∆Ψ) had an interactive effect on ci : ca. CS, WD and ΨPD were significantly correlated with each other. These results, along with those from multivariate analyses, underscored the pivotal role leaf : sapwood allocation (AL : AS), and water storage (CS) play in coordination between plant hydraulic and photosynthetic systems. This study uniquely explored the role of hydraulic traits in predicting species-specific photosynthetic variation based on optimality theory and highlights important mechanistic links within the plant carbon-water balance.

Dynamic regulation of water potential in Juniperus osteosperma mediates ecosystem carbon fluxes

Some plants exhibit dynamic hydraulic regulation, in which the strictness of hydraulic regulation (i.e. iso/anisohydry) changes in response to environmental conditions. However, the environmental controls over iso/anisohydry and the implications of flexible hydraulic regulation for plant productivity remain unknown. In Juniperus osteosperma, a drought-resistant dryland conifer, we collected a 5-month growing season time series of in situ, high temporal-resolution plant water potential ( Ψ ) and stand gross primary productivity (GPP). We quantified the stringency of hydraulic regulation associated with environmental covariates and evaluated how predawn water potential contributes to empirically predicting carbon uptake. Juniperus osteosperma showed less stringent hydraulic regulation (more anisohydric) after monsoon precipitation pulses, when soil moisture and atmospheric demand were high, and corresponded with GPP pulses. Predawn water potential matched the timing of GPP fluxes and improved estimates of GPP more strongly than soil and/or atmospheric moisture, notably resolving GPP underestimation before vegetation green-up. Flexible hydraulic regulation appears to allow J. osteosperma to prolong soil water extraction and, therefore, the period of high carbon uptake following monsoon precipitation pulses. Water potential and its dynamic regulation may account for why process-based and empirical models commonly underestimate the magnitude and temporal variability of dryland GPP.

Beyond the visible: Accounting for ultraviolet and far-red radiation in vegetation productivity and surface energy budgets

Photosynthetically active radiation (PAR) is typically defined as light with a wavelength within 400-700 nm. However, ultra-violet (UV) radiation within 280-400 nm and far-red (FR) radiation within 700-750 nm can also excite photosystems, though not as efficiently as PAR. Vegetation and land surface models (LSMs) typically do not explicitly account for UV's contribution to energy budgets or photosynthesis, nor FR's contribution to photosynthesis. However, whether neglecting UV and FR has significant impacts remains unknown. We explored how canopy radiative transfer (RT) and photosynthesis are impacted when explicitly implementing UV in the canopy RT model and accounting for UV and FR in the photosynthesis models within a next-generation LSM that can simulate hyperspectral canopy RT. We validated our improvements using photosynthesis measurements from plants under different light sources and intensities and surface reflection from an eddy-covariance tower. Our model simulations suggested that at the whole plant level, after accounting for UV and FR explicitly, chlorophyll content, leaf area index (LAI), clumping index, and solar radiation all impact the modeling of gross primary productivity (GPP). At the global scale, mean annual GPP within a grid would increase by up to 7.3% and the increase is proportional to LAI; globally integrated GPP increases by 4.6 PgC year-1 (3.8% of the GPP without accounting for UV + FR). Further, using PAR to proxy UV could overestimate surface albedo by more than 0.1, particularly in the boreal forests. Our results highlight the importance of improving UV and FR in canopy RT and photosynthesis modeling and the necessity to implement hyperspectral or multispectral canopy RT schemes in future vegetation and LSMs.

Amphistomy increases leaf photosynthesis more in coastal than montane plants of Hawaiian 'ilima (Sida fallax)

PREMISE

The adaptive significance of amphistomy (stomata on both upper and lower leaf surfaces) is unresolved. A widespread association between amphistomy and open, sunny habitats suggests the adaptive benefit of amphistomy may be greatest in these contexts, but this hypothesis has not been tested experimentally. Understanding amphistomy informs its potential as a target for crop improvement and paleoenvironment reconstruction.

METHODS

We developed a method to quantify "amphistomy advantage" (AA $\text{AA}$ ) as the log-ratio of photosynthesis in an amphistomatous leaf to that of the same leaf but with gas exchange blocked through the upper surface (pseudohypostomy). Humidity modulated stomatal conductance and thus enabled comparing photosynthesis at the same total stomatal conductance. We estimatedAA $\text{AA}$ and leaf traits in six coastal (open, sunny) and six montane (closed, shaded) populations of the indigenous Hawaiian species 'ilima (Sida fallax).

RESULTS

Coastal 'ilima leaves benefit 4.04 times more from amphistomy than montane leaves. Evidence was equivocal with respect to two hypotheses: (1) that coastal leaves benefit more because they are thicker and have lower CO2 conductance through the internal airspace and (2) that they benefit more because they have similar conductance on each surface, as opposed to most conductance being through the lower surface.

CONCLUSIONS

This is the first direct experimental evidence that amphistomy increases photosynthesis, consistent with the hypothesis that parallel pathways through upper and lower mesophyll increase CO2 supply to chloroplasts. The prevalence of amphistomatous leaves in open, sunny habitats can partially be explained by the increased benefit of amphistomy in "sun" leaves, but the mechanistic basis remains uncertain.

Dynamically optimizing stomatal conductance for maximum turgor-driven growth over diel and seasonal cycles

Stomata have recently been theorized to have evolved strategies that maximize turgor-driven growth over plants' lifetimes, finding support through steady-state solutions in which gas exchange, carbohydrate storage and growth have all reached equilibrium. However, plants do not operate near steady state as plant responses and environmental forcings vary diurnally and seasonally. It remains unclear how gas exchange, carbohydrate storage and growth should be dynamically coordinated for stomata to maximize growth. We simulated the gas exchange, carbohydrate storage and growth that dynamically maximize growth diurnally and annually. Additionally, we test whether the growth-optimization hypothesis explains nocturnal stomatal opening, particularly through diel changes in temperature, carbohydrate storage and demand. Year-long dynamic simulations captured realistic diurnal and seasonal patterns in gas exchange as well as realistic seasonal patterns in carbohydrate storage and growth, improving upon unrealistic carbohydrate responses in steady-state simulations. Diurnal patterns of carbohydrate storage and growth in day-long simulations were hindered by faulty modelling assumptions of cyclic carbohydrate storage over an individual day and synchronization of the expansive and hardening phases of growth, respectively. The growth-optimization hypothesis cannot currently explain nocturnal stomatal opening unless employing corrective 'fitness factors' or reframing the theory in a probabilistic manner, in which stomata adopt an inaccurate statistical 'memory' of night-time temperature. The growth-optimization hypothesis suggests that diurnal and seasonal patterns of stomatal conductance are driven by a dynamic carbon-use strategy that seeks to maintain homeostasis of carbohydrate reserves.

Optimal balancing of xylem efficiency and safety explains plant vulnerability to drought

In vast areas of the world, forests and vegetation are water limited and plant survival depends on the ability to avoid catastrophic hydraulic failure. Therefore, it is remarkable that plants take hydraulic risks by operating at water potentials (ψ) that induce partial failure of the water conduits (xylem). Here we present an eco-evolutionary optimality principle for xylem conduit design that explains this phenomenon based on the hypothesis that conductive efficiency and safety are optimally co-adapted to the environment. The model explains the relationship between the tolerance to negative water potential (ψ50 ) and the environmentally dependent minimum ψ (ψmin ) across a large number of species, and along the xylem pathway within individuals of two species studied. The wider hydraulic safety margin in gymnosperms compared to angiosperms can be explained as an adaptation to a higher susceptibility to accumulation of embolism. The model provides a novel optimality-based perspective on the relationship between xylem safety and efficiency.

Observed forest trait velocities have not kept pace with hydraulic stress from climate change

The extent to which future climate change will increase forest stress and the amount to which species and forest ecosystems can acclimate or adapt to increased stress is a major unknown. We used high-resolution maps of hydraulic traits representing the diversity in tree drought tolerance across the United States, a hydraulically enabled tree model, and forest inventory observations of demographic shifts to quantify the ability for within-species acclimation and between-species range shifts to mediate climate stress. We found that forests are likely to experience increases in both acute and chronic hydraulic stress with climate change. Based on current species distributions, regional hydraulic trait diversity was sufficient to buffer against increased stress in 88% of forested areas. However, observed trait velocities in 81% of forested areas are not keeping up with the rate required to ameliorate projected future stress without leaf area acclimation.

Plant water use theory should incorporate hypotheses about extreme environments, population ecology, and community ecology

Plant water use theory has largely been developed within a plant-performance paradigm that conceptualizes water use in terms of value for carbon gain and that sits within a neoclassical economic framework. This theory works very well in many contexts but does not consider other values of water to plants that could impact their fitness. Here, we survey a range of alternative hypotheses for drivers of water use and stomatal regulation. These hypotheses are organized around relevance to extreme environments, population ecology, and community ecology. Most of these hypotheses are not yet empirically tested and some are controversial (e.g. requiring more agency and behavior than is commonly believed possible for plants). Some hypotheses, especially those focused around using water to avoid thermal stress, using water to promote reproduction instead of growth, and using water to hoard it, may be useful to incorporate into theory or to implement in Earth System Models.

Daytime stomatal regulation in mature temperate trees prioritizes stem rehydration at night

Trees remain sufficiently hydrated during drought by closing stomata and reducing canopy conductance (G_c ) in response to variations in atmospheric water demand and soil water availability. Thresholds that control the reduction of G_c are proposed to optimize hydraulic safety against carbon assimilation efficiency. However, the link between G_c and the ability of stem tissues to rehydrate at night remains unclear. We investigated whether species-specific G_c responses aim to prevent branch embolisms, or enable night-time stem rehydration, which is critical for turgor-dependent growth. For this, we used a unique combination of concurrent dendrometer, sap flow and leaf water potential measurements and collected branch-vulnerability curves of six common European tree species. Species-specific G_c reduction was weakly related to the water potentials at which 50% of branch xylem conductivity is lost (P₅₀ ). Instead, we found a stronger relationship with stem rehydration. Species with a stronger G_c control were less effective at refilling stem-water storage as the soil dries, which appeared related to their xylem architecture. Our findings highlight the importance of stem rehydration for water-use regulation in mature trees, which likely relates to the maintenance of adequate stem turgor. We thus conclude that stem rehydration must complement the widely accepted safety-efficiency stomatal control paradigm.

Do stomata optimize turgor-driven growth? A new framework for integrating stomata response with whole-plant hydraulics and carbon balance

Every existing optimal stomatal model uses photosynthetic carbon assimilation as a proxy for plant evolutionary fitness. However, assimilation and growth are often decoupled, making assimilation less ideal for representing fitness when optimizing stomatal conductance to water vapor and carbon dioxide. Instead, growth should be considered a closer proxy for fitness. We hypothesize stomata have evolved to maximize turgor-driven growth, instead of assimilation, over entire plants' lifetimes, improving their abilities to compete and reproduce. We develop a stomata model that dynamically maximizes whole-stem growth following principles from turgor-driven growth models. Stomata open to assimilate carbohydrates that supply growth and osmotically generate turgor, while stomata close to prevent losses of turgor and growth due to negative water potentials. In steady state, the growth optimization model captures realistic stomatal, growth, and carbohydrate responses to environmental cues, reconciles conflicting interpretations within existing stomatal optimization theories, and explains patterns of carbohydrate storage and xylem conductance observed during and after drought. Our growth optimization hypothesis introduces a new paradigm for stomatal optimization models, elevates the role of whole-plant carbon use and carbon storage in stomatal functioning, and has the potential to simultaneously predict gross productivity, net productivity, and plant mortality through a single, consistent modeling framework.

Kinetic factors of physiology and the dynamic light environment influence the economic landscape of short-term hydraulic risk

Optimality-based models of stomatal conductance unify biophysical and evolutionary constraints and can improve predictions of land-atmosphere carbon and water exchange. Recent models incorporate hydraulic constraints by penalizing excessive stomatal opening in relation to hydraulic damage caused by low water potentials. We used simulation models to test whether penalties based solely on vulnerability curves adequately represent the optimality hypothesis, given that they exclude the effects of kinetic factors on stomatal behavior and integrated carbon balance. To quantify the effects of nonsteady-state phenomena on the landscape of short-term hydraulic risk, we simulated diurnal dynamics of leaf physiology for 10 000 patches of leaf in a canopy and used a ray-tracing model, Helios, to simulate realistic variation in sunfleck dynamics. Our simulations demonstrated that kinetic parameters of leaf physiology and sunfleck properties influence the economic landscape of short-term hydraulic risk, as characterized by the effect of stomatal strategy (gauged by the water potential causing a 50% hydraulic penalty) on both aggregated carbon gain and the aggregated carbon cost of short-term hydraulic risk. Hydraulic penalties in optimization models should be generalized to allow their parameters to account for kinetic factors, in addition to parameters of hydraulic vulnerability.

Constant hydraulic supply and ABA dynamics facilitate the trade-offs in water and carbon

Carbon-water trade-offs in plants are adjusted through stomatal regulation. Stomatal opening enables carbon uptake and plant growth, whereas plants circumvent drought by closing stomata. The specific effects of leaf position and age on stomatal behavior remain largely unknown, especially under edaphic and atmospheric drought. Here, we compared stomatal conductance (g_s ) across the canopy of tomato during soil drying. We measured gas exchange, foliage ABA level and soil-plant hydraulics under increasing vapor pressure deficit (VPD). Our results indicate a strong effect of canopy position on stomatal behavior, especially under hydrated soil conditions and relatively low VPD. In wet soil (soil water potential > -50 kPa), upper canopy leaves had the highest g_s (0.727 ± 0.154 mol m^-2 s^-1) and assimilation rate (A; 23.4 ± 3.9 µmol m^-2 s^-1) compared to the leaves at a medium height of the canopy (g_s : 0.159 ± 0.060 mol m² s^-1; A: 15.9 ± 3.8 µmol m^-2 s^-1). Under increasing VPD (from 1.8 to 2.6 kPa), g_s , A and transpiration were initially impacted by leaf position rather than leaf age. However, under high VPD (2.6 kPa), age effect outweighed position effect. The soil-leaf hydraulic conductance was similar in all leaves. Foliage ABA levels increased with rising VPD in mature leaves at medium height (217.56 ± 85 ng g^-1 FW) compared to upper canopy leaves (85.36 ± 34 ng g^-1 FW). Under soil drought (< -50 kPa), stomata closed in all leaves resulting in no differences in g_s across the canopy. We conclude that constant hydraulic supply and ABA dynamics facilitate preferential stomatal behavior and carbon-water trade-offs across the canopy. These findings are fundamental in understanding variations within the canopy, which helps in engineering future crops, especially in the face of climate change.

Stomata conductance as a goalkeeper for increased photosynthetic efficiency

100-120 words. References should not be included. Abbreviations should be avoided as far as possible. Low stomatal conductance (g_s) poses a major constraint for improving photosynthetic efficiency for greater yield. Options at the molecular, leaf, canopy, and even the whole plant scales can be developed to enhance g_s for greater light and water use efficiencies. Among these, many genes regulating stomatal development and stomatal movement have been discovered and manipulated to increase light and water use efficiencies under well-watered, drought, or facility agriculture conditions with the manual-controlled growth environmental. Optimization of canopy conductance to increase whole plant photosynthesis with full consideration of the heterogeneities in g_s, microclimates and leaf ontology inside the canopy represents a largely uncharted area to improve crop efficiency.

Towards a unified theory of plant photosynthesis and hydraulics

The global carbon and water cycles are governed by the coupling of CO₂ and water vapour exchanges through the leaves of terrestrial plants, controlled by plant adaptations to balance carbon gains and hydraulic risks. We introduce a trait-based optimality theory that unifies the treatment of stomatal responses and biochemical acclimation of plants to environments changing on multiple timescales. Tested with experimental data from 18 species, our model successfully predicts the simultaneous decline in carbon assimilation rate, stomatal conductance and photosynthetic capacity during progressive soil drought. It also correctly predicts the dependencies of gas exchange on atmospheric vapour pressure deficit, temperature and CO₂. Model predictions are also consistent with widely observed empirical patterns, such as the distribution of hydraulic strategies. Our unified theory opens new avenues for reliably modelling the interactive effects of drying soil and rising atmospheric CO₂ on global photosynthesis and transpiration.

Simple Models Outperform More Complex Big-Leaf Models of Daily Transpiration in Forested Biomes

Transpiration makes up the bulk of total evaporation in forested environments yet remains challenging to predict at landscape-to-global scales. We harnessed independent estimates of daily transpiration derived from co-located sap flow and eddy-covariance measurement systems and applied the triple collocation technique to evaluate predictions from big leaf models requiring no calibration. In total, four models in 608 unique configurations were evaluated at 21 forested sites spanning a wide diversity of biophysical attributes and environmental backgrounds. We found that simpler models that neither explicitly represented aerodynamic forcing nor canopy conductance achieved higher accuracy and signal-to-noise levels when optimally configured (rRMSE = 20%; R 2 = 0.89). Irrespective of model type, optimal configurations were those making use of key plant functional type dependent parameters, daily LAI, and constraints based on atmospheric moisture demand over soil moisture supply. Our findings have implications for more informed water resource management based on hydrological modeling and remote sensing.

Dynamic Energy Budget models: fertile ground for understanding resource allocation in plants in a changing world

Climate change is having dramatic effects on the diversity and distribution of species. Many of these effects are mediated by how an organism's physiological patterns of resource allocation translate into fitness through effects on growth, survival and reproduction. Empirically, resource allocation is challenging to measure directly and so has often been approached using mathematical models, such as Dynamic Energy Budget (DEB) models. The fact that all plants require a very similar set of exogenous resources, namely light, water and nutrients, integrates well with the DEB framework in which a small number of variables and processes linked through pathways represent an organism's state as it changes through time. Most DEB theory has been developed in reference to animals and microorganisms. However, terrestrial vascular plants differ from these organisms in fundamental ways that make resource allocation, and the trade-offs and feedbacks arising from it, particularly fundamental to their life histories, but also challenging to represent using existing DEB theory. Here, we describe key features of the anatomy, morphology, physiology, biochemistry, and ecology of terrestrial vascular plants that should be considered in the development of a generic DEB model for plants. We then describe possible approaches to doing so using existing DEB theory and point out features that may require significant development for DEB theory to accommodate them. We end by presenting a generic DEB model for plants that accounts for many of these key features and describing gaps that would need to be addressed for DEB theory to predict the responses of plants to climate change. DEB models offer a powerful and generalizable framework for modelling resource allocation in terrestrial vascular plants, and our review contributes a framework for expansion and development of DEB theory to address how plants respond to anthropogenic change.

GriddingMachine, a database and software for Earth system modeling at global and regional scales

Land and Earth system modeling is moving towards more explicit biophysical representations, requiring increasing variety of datasets for initialization and benchmarking. However, researchers often have difficulties in identifying and integrating non-standardized datasets from various sources. We aim towards a standardized database and one-stop distribution method of global datasets. Here, we present the GriddingMachine as (1) a database of global-scale datasets commonly used to parameterize or benchmark the models, from plant traits to vegetation indices and geophysical information and (2) a cross-platform open source software to download and request a subset of datasets with only a few lines of code. The GriddingMachine datasets can be accessed either manually through traditional HTTP, or automatically using modern programming languages including Julia, Matlab, Octave, Python, and R. The GriddingMachine collections can be used for any land and Earth modeling framework and ecological research at the regional and global scales, and the number of datasets will continue to grow to meet the increasing needs of research communities.

Stomatal responses of Eucalyptus spp. under drought can be predicted with a gain-risk optimization model

The frequency and severity of drought events are expected to increase due to climate change, with optimal environmental conditions for forestry likely to shift. Modeling plant responses to a changing climate is therefore vital. We tested the process-based gain-risk model to predict stomatal responses to drought of two Eucalyptus hybrids. The process-based gain-risk model has the advantage that all the parameters used within the model are based on measurable plant traits. The gain-risk model proposes that plants optimize photosynthetic gain while minimizing a hydraulic cost. Previous versions of the model used hydraulic risk as a cost function; however, they did not account for delayed or reduced hydraulic recovery rates from embolism post-drought. Hydraulic recovery has been seen in many species, however it is still unclear how this inclusion of a partial or delayed hydraulic recovery would affect the predictive power of the gain-risk model. Many hydraulic parameters required by the model are also difficult to measure and are not freely available. We therefore tested a simplified gain-risk model that includes a delayed or reduced hydraulic recovery component post-drought. The simplified gain-risk model performed well at predicting stomatal responses in both Eucalyptus grandis × camaldulensis (GC) and Eucalyptus urophylla × grandis (UG). In this study two distinct strategies were seen between GC and UG, with GC being more resistant to embolism formation, however it could not recover hydraulic conductance compared with UG. The inclusion of a delayed or reduced hydraulic recovery component slightly improved model predictions for GC, however not for UG, which can be related to UG being able to recover lost hydraulic conductance and therefore can maintain stomatal conductance regardless of hydraulic risk. Even though the gain-risk model shows promise in predicting plant responses, more information is needed regarding hydraulic recovery after drought.

Integrating plant physiology and community ecology across scales through trait-based models to predict drought mortality

Forests are a critical carbon sink and widespread tree mortality resulting from climate-induced drought stress has the potential to alter forests from a carbon sink to a source, causing a positive feedback on climate change. Process-based vegetation models aim to represent the current understanding of the underlying mechanisms governing plant physiological and ecological responses to climate. Yet model accuracy varies across scales, and regional-scale model predictive skill is frequently poor when compared with observations of drought-driven mortality. I propose a framework that leverages differences in model predictive skill across spatial scales, mismatches between model predictions and observations, and differences in the mechanisms included and absent across models to advance the understanding of the physiological and ecological processes driving observed patterns drought-driven mortality.

CO2 fertilization of terrestrial photosynthesis inferred from site to global scales

The magnitude of the CO₂ fertilization effect on terrestrial photosynthesis is uncertain because it is not directly observed and is subject to confounding effects of climatic variability. We apply three well-established eco-evolutionary optimality theories of gas exchange and photosynthesis, constraining the main processes of CO₂ fertilization using measurable variables. Using this framework, we provide robust observationally inferred evidence that a strong CO₂ fertilization effect is detectable in globally distributed eddy covariance networks. Applying our method to upscale photosynthesis globally, we find that the magnitude of the CO₂ fertilization effect is comparable to its in situ counterpart but highlight the potential for substantial underestimation of this effect in tropical forests for many reflectance-based satellite photosynthesis products.

Global photosynthesis is increasing with elevated atmospheric CO₂ concentrations, a response known as the CO₂ fertilization effect (CFE), but the key processes of CFE are not constrained and therefore remain uncertain. Here, we quantify CFE by combining observations from a globally distributed network of eddy covariance measurements with an analytical framework based on three well-established photosynthetic optimization theories. We report a strong enhancement of photosynthesis across the observational network (9.1 gC m⁻² year⁻²) and show that the CFE is responsible for 44% of the gross primary production (GPP) enhancement since the 2000s, with additional contributions primarily from warming (28%). Soil moisture and specific humidity are the two largest contributors to GPP interannual variation through their influences on plant hydraulics. Applying our framework to satellite observations and meteorological reanalysis data, we diagnose a global CO₂-induced GPP trend of 4.4 gC m⁻² year⁻², which is at least one-third stronger than the median trends of 13 dynamic global vegetation models and eight satellite-derived GPP products, mainly because of their differences in the magnitude of CFE in evergreen broadleaf forests. These results highlight the critical role that CFE has played in the global carbon cycle in recent decades.

Stomatal closure during water deficit is controlled by below-ground hydraulics

BACKGROUND AND AIMS

Stomatal closure allows plants to promptly respond to water shortage. Although the coordination between stomatal regulation, leaf and xylem hydraulics has been extensively investigated, the impact of below-ground hydraulics on stomatal regulation remains unknown.

METHODS

We used a novel root pressure chamber to measure, during soil drying, the relation between transpiration rate (E) and leaf xylem water pressure (ψleaf-x) in tomato shoots grafted onto two contrasting rootstocks, a long and a short one. In parallel, we also measured the E(ψleaf-x) relation without pressurization. A soil-plant hydraulic model was used to reproduce the measurements. We hypothesize that (1) stomata close when the E(ψleaf-x) relation becomes non-linear and (2) non-linearity occurs at higher soil water contents and lower transpiration rates in short-rooted plants.

KEY RESULTS

The E(ψleaf-x) relation was linear in wet conditions and became non-linear as the soil dried. Changing below-ground traits (i.e. root system) significantly affected the E(ψleaf-x) relation during soil drying. Plants with shorter root systems required larger gradients in soil water pressure to sustain the same transpiration rate and exhibited an earlier non-linearity and stomatal closure.

CONCLUSIONS

We conclude that, during soil drying, stomatal regulation is controlled by below-ground hydraulics in a predictable way. The model suggests that the loss of hydraulic conductivity occurred in soil. These results prove that stomatal regulation is intimately tied to root and soil hydraulic conductances.

Exploring optimal stomatal control under alternative hypotheses for the regulation of plant sources and sinks

Experimental evidence that nonstomatal limitations to photosynthesis (NSLs) correlate with leaf sugar and/or leaf water status suggests the possibility that stomata adjust to maximise photosynthesis through a trade-off between leaf CO₂ supply and NSLs, potentially involving source-sink interactions. However, the mechanisms regulating NSLs and sink strength, as well as their implications for stomatal control, remain uncertain. We used an analytically solvable model to explore optimal stomatal control under alternative hypotheses for source and sink regulation. We assumed that either leaf sugar concentration or leaf water potential regulates NSLs, and that either phloem turgor pressure or phloem sugar concentration regulates sink phloem unloading. All hypotheses led to realistic stomatal responses to light, CO₂ and air humidity, including conservative behaviour for the intercellular-to-atmospheric CO₂ concentration ratio. Sugar-regulated and water-regulated NSLs are distinguished by the presence/absence of a stomatal closure response to changing sink strength. Turgor-regulated and sugar-regulated phloem unloading are distinguished by the presence/absence of stomatal closure under drought and avoidance/occurrence of negative phloem turgor. Results from girdling and drought experiments on Pinus sylvestris, Betula pendula, Populus tremula and Picea abies saplings are consistent with optimal stomatal control under sugar-regulated NSLs and turgor-regulated unloading. Our analytical results provide a simple representation of stomatal responses to above-ground and below-ground environmental factors and sink activity.

Hydraulic-stomatal coordination in tree seedlings: tight correlation across environments and ontogeny in Acer pseudoplatanus

Hydraulic conductance is recognized as a major determinant of gas exchange and productivity. However, whether this also applies to seedlings, a critically important stage for vegetation regeneration, has been largely unknown. We analyzed the hydraulic and stomatal conductance of leaves and shoots for 6-wk-old Acer pseudoplatanus seedlings emerging in different lowland and treeline habitats and under glasshouse conditions, respectively, as well as on 9-, 15- and 18-wk-old plants, and related findings to leaf and xylem anatomical traits. Treeline seedlings had higher leaf area-specific shoot hydraulic conductance (Kshoot-L ), and stomatal conductance (gs ), associated with wider xylem conduits, lower leaf area and higher stomatal density than lowland and glasshouse-grown plants. Across the first 18 wk of development, seedlings increased four-fold in absolute shoot hydraulic conductance (Kshoot ) and declined by half in Kshoot-L , with correlated shifts in xylem and leaf anatomy. Distal leaves had higher leaf hydraulic conductance (Kleaf ) and gs compared to basal leaves. Seedlings show strong variation across growth environments and ontogenetic shifts in hydraulic and anatomical parameters. Across growth sites, ontogenetic stages and leaf orders, gs was tightly correlated with Kshoot-L and Kleaf , balancing hydraulic supply with demand for the earliest stages of seedling establishment.

Spatial variations in stomatal traits and their coordination with leaf traits in Quercus variabilis across Eastern Asia

The stomatal traits influence ecosystem carbon-water fluxes and play essential roles that enable plants to adapt to changing environmental conditions. However, how stomatal traits vary along a large climate gradient and whether stomatal traits coordinated with other leaf functional traits in response to environmental changes remain unclear. We investigated the stomatal density (SD), stomatal size (SS), and leaf traits (leaf area (LA), leaf mass per area (LMA), and vein density (VD)) of 44 in situ Quercus variabilis populations across Eastern Asia (24 to 51.8°N, 99 to 137°E) and 15 populations grown in a common garden, and evaluated their relationships with environmental factors. Stepwise multiple regression showed that the SD was significantly associated with mean annual precipitation (MAP), LMA, and VD, and the SS with latitude, mean annual temperature (MAT), mean monthly solar radiation (MMSR), and VD. The SD was positively correlated with the LMA, while the SS was negatively correlated with the VD. The SD and LMA increased with decreasing precipitation, which indicated that they may coordinate to commonly enhance plant resistance against drought. The SS decreased; however, the VD increased with temperature. This implied that plants might further reduce their SS by increasing VD limitations under global warming. In the common garden, plants exhibited a higher SD and VD and lower SS and LA compared to those in the field; however, no relation between the stomatal and leaf traits was observed. Our results suggested that stomatal traits have high environmental plasticity and are highly coordinated with other leaf functional traits in response to environmental changes. Nevertheless, this coordination may have been formed through long-term adaptations, rather than over short time spans.

Eco-evolutionary optimality as a means to improve vegetation and land-surface models

Global vegetation and land-surface models embody interdisciplinary scientific understanding of the behaviour of plants and ecosystems, and are indispensable to project the impacts of environmental change on vegetation and the interactions between vegetation and climate. However, systematic errors and persistently large differences among carbon and water cycle projections by different models highlight the limitations of current process formulations. In this review, focusing on core plant functions in the terrestrial carbon and water cycles, we show how unifying hypotheses derived from eco-evolutionary optimality (EEO) principles can provide novel, parameter-sparse representations of plant and vegetation processes. We present case studies that demonstrate how EEO generates parsimonious representations of core, leaf-level processes that are individually testable and supported by evidence. EEO approaches to photosynthesis and primary production, dark respiration and stomatal behaviour are ripe for implementation in global models. EEO approaches to other important traits, including the leaf economics spectrum and applications of EEO at the community level are active research areas. Independently tested modules emerging from EEO studies could profitably be integrated into modelling frameworks that account for the multiple time scales on which plants and plant communities adjust to environmental change.

Soil hydraulics affect the degree of isohydricity

A soil–plant hydraulic model shows that the degree of isohydricity is constrained by below-ground hydraulic limitations and that the shape of plant water potential (WP) curve depends on soil hydraulics.

Climate and functional traits jointly mediate tree water-use strategies

Tree water use is central to plant function and ecosystem fluxes. However, it is still unknown how organ-level water-relations traits are coordinated to determine whole-tree water-use strategies in response to drought, and whether this coordination depends on climate. Here we used a global sap flow database (SAPFLUXNET) to study the response of water use, in terms of whole-tree canopy conductance (G), to vapour pressure deficit (VPD) and to soil water content (SWC) for 142 tree species. We investigated the individual and coordinated effect of six water-relations traits (vulnerability to embolism, Huber value, hydraulic conductivity, turgor-loss point, rooting depth and leaf size) on water-use parameters, also accounting for the effect of tree height and climate (mean annual precipitation, MAP). Reference G and its sensitivity to VPD were tightly coordinated with water-relations traits rather than with MAP. Species with efficient xylem transport had higher canopy conductance but also higher sensitivity to VPD. Moreover, we found that angiosperms had higher reference G and higher sensitivity to VPD than did gymnosperms. Our results highlight the need to consider trait integration and reveal the complications and challenges of defining a single, whole-plant resource use spectrum ranging from 'acquisitive' to 'conservative'.

Parsimony vs predictive and functional performance of three stomatal optimization principles in a big-leaf framework

Stomatal optimization models can improve estimates of water and carbon fluxes with relatively low complexity, yet there is no consensus on which formulations are most appropriate for ecosystem-scale applications. We implemented three existing analytical equations for stomatal conductance, based on different water penalty functions, in a big-leaf comparison framework, and determined which optimization principles were most consistent with flux tower observations from different biomes. We used information theory to dissect controls of soil water supply and atmospheric demand on evapotranspiration in wet to dry conditions and to quantify missing or inadequate information in model variants. We ranked stomatal optimization principles based on parameter uncertainty, parsimony, predictive accuracy, and functional accuracy of the interactions between soil moisture, vapor pressure deficit, and evapotranspiration. Performance was high for all model variants. Water penalty functions with explicit representation of plant hydraulics did not substantially improve predictive or functional accuracy of ecosystem-scale evapotranspiration estimates, and parameterizations were more uncertain, despite having physiological underpinnings at the plant level. Stomatal optimization based on water use efficiency thus provided more information about ecosystem-scale evapotranspiration compared to those based on xylem vulnerability and proved more useful in improving ecosystem-scale models with less complexity.

Coupled whole-tree optimality and xylem hydraulics explain dynamic biomass partitioning

Trees partition biomass in response to resource limitation and physiological activity. It is presumed that these strategies evolved to optimize some measure of fitness. If the optimization criterion can be specified, then allometry can be modeled from first principles without prescribed parameterization. We present the Tree Hydraulics and Optimal Resource Partitioning (THORP) model, which optimizes allometry by estimating allocation fractions to organs as proportional to their ratio of marginal gain to marginal cost, where gain is net canopy photosynthesis rate, and costs are senescence rates. Root total biomass and profile shape are predicted simultaneously by a unified optimization. Optimal partitioning is solved by a numerically efficient analytical solution. THORP's predictions agree with reported tree biomass partitioning in response to size, water limitations, elevated CO₂ and pruning. Roots were sensitive to soil moisture profiles and grew down to the groundwater table when present. Groundwater buffered against water stress regardless of meteorology, stabilizing allometry and root profiles as deep as c. 30 m. Much of plant allometry can be explained by hydraulic considerations. However, nutrient limitations cannot be fully ignored. Rooting mass and profiles were synchronized with hydrological conditions and groundwater even at considerable depths, illustrating that the below ground shapes whole-tree allometry.

Understanding and predicting forest mortality in the western United States using long-term forest inventory data and modeled hydraulic damage

Global warming is expected to exacerbate the duration and intensity of droughts in the western United States, which may lead to increased tree mortality. A prevailing proximal mechanism of drought-induced tree mortality is hydraulic damage, but predicting tree mortality from hydraulic theory and climate data still remains a major scientific challenge. We used forest inventory data and a plant hydraulic model (HM) to address three questions: can we capture regional patterns of drought-induced tree mortality with HM-predicted damage thresholds; do HM metrics improve predictions of mortality across broad spatial areas; and what are the dominant controls of forest mortality when considering stand characteristics, climate metrics, and simulated hydraulic stress? We found that the amount of variance explained by models predicting mortality was limited (R² median = 0.10, R² range: 0.00-0.52). HM outputs, including hydraulic damage and carbon assimilation diagnostics, moderately improve mortality prediction across the western US compared with models using stand and climate predictors alone. Among factors considered, metrics of stand density and tree size tended to be some of the most critical factors explaining mortality, probably highlighting the important roles of structural overshoot, stand development, and biotic agent host selection and outbreaks in mortality patterns.

Importance of hydraulic strategy trade-offs in structuring response of canopy trees to extreme drought in central Amazon

Plant ecophysiological trade-offs between different strategies for tolerating stresses are widely theorized to shape forest functional diversity and vulnerability to climate change. However, trade-offs between hydraulic and stomatal regulation during natural droughts remain under-studied, especially in tropical forests. We investigated eleven mature forest canopy trees in central Amazonia during the strong 2015 El Niño. We found greater xylem embolism resistance ([Formula: see text] = - 3.3 ± 0.8 MPa) and hydraulic safety margin (HSM = 2.12 ± 0.57 MPa) than previously observed in more precipitation-seasonal rainforests of eastern Amazonia and central America. We also discovered that taller trees exhibited lower embolism resistance and greater stomatal sensitivity, a height-structured trade-off between hydraulic resistance and active stomatal regulation. Such active regulation of tree water status, triggered by the onset of stem embolism, acted as a feedback to avoid further increases in embolism, and also explained declines in photosynthesis and transpiration. These results suggest that canopy trees exhibit a conservative hydraulic strategy to endure drought, with trade-offs between investment in xylem to reduce vulnerability to hydraulic failure, and active stomatal regulation to protect against low water potentials. These findings improve our understanding of strategies in tropical forest canopies and contribute to more accurate prediction of drought responses.

Optimization theory explains nighttime stomatal responses

Nocturnal transpiration is widely observed across species and biomes, and may significantly impact global water, carbon, and energy budgets. However, it remains elusive why plants lose water at night and how to model it at large scales. We hypothesized that plants optimize nighttime leaf diffusive conductance (g_wn ) to balance potential daytime photosynthetic benefits and nocturnal transpiration benefits. We quantified nighttime benefits from respiratory reductions due to evaporative leaf cooling. We described nighttime costs in terms of a reduced carbon gain during the day because of water use at night. We measured nighttime stomatal responses and tested our model with water birch (Betula occidentalis) saplings grown in a glasshouse. The g_wn of water birch decreased with drier soil, higher atmospheric CO₂ , wetter air, lower leaf temperature, and lower leaf respiration rate. Our model predicted all these responses correctly, except for the response of g_wn to air humidity. Our results also suggested that the slow decrease in g_wn after sunset could be associated with decreasing leaf respiration. The optimality-based nocturnal transpiration model smoothly integrates with daytime stomatal optimization approaches, and thus has the potential to quantitatively predict nocturnal transpiration across space and time.

On the inter- and intra-annual variability of ecosystem evapotranspiration and water use efficiency of an oak savanna and annual grassland subjected to booms and busts in rainfall

Whether annual evapotranspiration of native ecosystems is increasing or decreasing with time as CO₂ concentrations are rising, the climate is warming and rainfall experiences booms and busts, remains an unanswered question in the field of global change biology. To answer this question, we measured evapotranspiration and carbon dioxide exchange over and under an oak savanna and over an annual grassland in the Mediterranean climate of California, USA, from 2001 through 2019 with the eddy covariance method; during this 19-year period, CO₂ rose 40 ppm, air temperature increased by 1°C and annual rainfall ranged between 133 and 890 mm/year. No temporal trend in evapotranspiration or water use efficiency was observed over this time duration. Many competing positive and negative feedbacks among stomatal sensitivity to carbon dioxide concentrations, soil moisture, and vapor pressure deficit, the impact of temperature on saturation vapor pressure and access to groundwater muted the response of evapotranspiration to its changing world when integrated to the ecosystem scale and annual time steps. At the intra-annual time scale, we found that plants transmit information on soil moisture status through their influence on the vapor pressure deficit of the atmospheric boundary layer. The inter-annual variations in evaporative water use by the savanna and annual grassland were relatively decoupled from the booms and busts in rainfall. Instead, variations in length of growing season and access to groundwater explained much of this year-to-year variation in annual evapotranspiration. The access of groundwater by the oak savanna may make these ecosystems more robust in a warmer world, than was previously thought. This is a scale emergent property that needs better consideration in coupled climate-ecosystem models.

Optimization can provide the fundamental link between leaf photosynthesis, gas exchange and water relations

Tight coordination in the photosynthetic, gas exchange and water supply capacities of leaves is a globally conserved trend across land plants. Strong selective constraints on leaf carbon gain create the opportunity to use quantitative optimization theory to understand the connected evolution of leaf photosynthesis and water relations. We developed an analytical optimization model that maximizes the long-term rate of leaf carbon gain, given the carbon costs in building and maintaining stomata, leaf hydraulics and osmotic pressure. Our model demonstrates that selection for optimal gain should drive coordination between key photosynthetic, gas exchange and water relations traits. It also provides predictions of adaptation to drought and the relative costs of key leaf functional traits. Our results show that optimization in terms of carbon gain, given the carbon costs of physiological traits, successfully unites leaf photosynthesis and water relations and provides a quantitative framework to consider leaf functional evolution and adaptation.