Tropical tree mortality has increased with rising atmospheric water stress.

Evidence exists that tree mortality is accelerating in some regions of the tropics1,2, with profound consequences for the future of the tropical carbon sink and the global anthropogenic carbon budget left to limit peak global warming below 2 °C. However, the mechanisms that may be driving such mortality changes and whether particular species are especially vulnerable remain unclear3-8. Here we analyse a 49-year record of tree dynamics from 24 old-growth forest plots encompassing a broad climatic gradient across the Australian moist tropics and find that annual tree mortality risk has, on average, doubled across all plots and species over the last 35 years, indicating a potential halving in life expectancy and carbon residence time. Associated losses in biomass were not offset by gains from growth and recruitment. Plots in less moist local climates presented higher average mortality risk, but local mean climate did not predict the pace of temporal increase in mortality risk. Species varied in the trajectories of their mortality risk, with the highest average risk found nearer to the upper end of the atmospheric vapour pressure deficit niches of species. A long-term increase in vapour pressure deficit was evident across the region, suggesting that thresholds involving atmospheric water stress, driven by global warming, may be a primary cause of increasing tree mortality in moist tropical forests.

Mesophyll airspace unsaturation drives C4 plant success under vapor pressure deficit stress.

A fundamental assumption in plant science posits that leaf air spaces remain vapor saturated, leading to the predominant view that stomata alone control leaf water loss. This concept has been pivotal in photosynthesis and water-use efficiency research. However, recent evidence has refuted this longstanding assumption by providing evidence of unsaturation in the leaf air space of C3 plants under relatively mild vapor pressure deficit (VPD) stress. This phenomenon represents a nonstomatal mechanism restricting water loss from the mesophyll. The potential ubiquity and physiological implications of this phenomenon, its driving mechanisms in different plant species and habitats, and its interaction with other ecological adaptations have. In this context, C4 plants spark particular interest for their importance as crops, bundle sheath cells' unique anatomical characteristics and specialized functions, and notably higher water-use efficiency relative to C3 plants. Here, we confirm reduced relative humidities in the substomatal cavity of the C4 plants maize, sorghum, and proso millet down to 80% under mild VPD stress. We demonstrate the critical role of nonstomatal control in these plants, indicating that the role of the CO2 concentration mechanism in CO2 management at a high VPD may have been overestimated. Our findings offer a mechanistic reconciliation between discrepancies in CO2 and VPD responses reported in C4 species. They also reveal that nonstomatal control is integral to maintaining an advantageous microclimate of relatively higher CO2 concentrations in the mesophyll air space of C4 plants for carbon fixation, proving vital when these plants face VPD stress.

Impacts of elevated temperature and vapour pressure deficit on leaf gas exchange and plant growth across six tropical rainforest tree species.

Elevated air temperature (Tair) and vapour pressure deficit (VPDair) significantly influence plant functioning, yet their relative impacts are difficult to disentangle. We examined the effects of elevated Tair (+6°C) and VPDair (+0.7 kPa) on the growth and physiology of six tropical tree species. Saplings were grown under well-watered conditions in climate-controlled glasshouses for 6 months under three treatments: (1) low Tair and low VPDair, (2) high Tair and low VPDair, and (3) high Tair and high VPDair. To assess acclimation, physiological parameters were measured at a set temperature. Warm-grown plants grown under elevated VPDair had significantly reduced stomatal conductance and increased instantaneous water use efficiency compared to plants grown under low VPDair. Photosynthetic biochemistry and thermal tolerance (Tcrit) were unaffected by VPDair, but elevated Tair caused Jmax25 to decrease and Tcrit to increase. Sapling biomass accumulation for all species responded positively to an increase in Tair, but elevated VPDair limited growth. This study shows that stomatal limitation caused by even moderate increases in VPDair can decrease productivity and growth rates in tropical species independently from Tair and has important implications for modelling the impacts of climate change on tropical forests.

Harnessing photosynthetic C18O16O discrimination dynamics under leaf water nonsteady state to estimate mesophyll conductance: a new, regression-based method.

Mesophyll conductance (gm) is a crucial plant trait that can significantly limit photosynthesis. Measurement of photosynthetic C18O16O discrimination (Δ18O) has proved to be the only viable means of resolving gm in both C3 and C4 plants. However, the currently available methods to exploit Δ18O for gm estimation are error prone due to their inadequacy in constraining the degree of oxygen isotope exchange (θ) during mesophyll CO2 hydration. Here, we capitalized on experimental manipulation of leaf water isotopic dynamics to establish a novel, nonsteady state, regression-based approach for simultaneous determination of gm and θ from online Δ18O measurements. We demonstrated the methodological and theoretical robustness of this new Δ18O-gm estimation approach and showed through measurements on several C3 and C4 species that this approach can serve as a benchmark method against which to identify previously-unrecognized biases of the existing Δ18O-gm methods. Our results highlight the unique value of this nonsteady state-based approach for contributing to ongoing efforts toward quantitative understanding of mesophyll conductance for crop yield improvement and carbon cycle modeling.

Dry inside: progressive unsaturation within leaves with increasing vapour pressure deficit affects estimation of key leaf gas exchange parameters.

Climate change not only leads to higher air temperatures but also increases the vapour pressure deficit (VPD) of the air. Understanding the direct effect of VPD on leaf gas exchange is crucial for precise modelling of stomatal functioning. We conducted combined leaf gas exchange and online isotope discrimination measurements on four common European tree species across a VPD range of 0.8-3.6 kPa, while maintaining constant temperatures without soil water limitation. In addition to applying the standard assumption of saturated vapour pressure inside leaves (ei), we inferred ei from oxygen isotope discrimination of CO2 and water vapour. ei desaturated progressively with increasing VPD, consistently across species, resulting in an intercellular relative humidity as low as 0.73 ± 0.11 at the highest tested VPD. Assuming saturation of ei overestimated the extent of reductions in stomatal conductance and CO2 mole fraction inside leaves in response to increasing VPD compared with calculations that accounted for unsaturation. In addition, a significant decrease in mesophyll conductance with increasing VPD only occurred when the unsaturation of ei was considered. We suggest that the possibility of unsaturated ei should not be overlooked in measurements related to leaf gas exchange and in stomatal models, especially at high VPD.

Uncoupling of stomatal conductance and photosynthesis at high temperatures: mechanistic insights from online stable isotope techniques.

The strong covariation of temperature and vapour pressure deficit (VPD) in nature limits our understanding of the direct effects of temperature on leaf gas exchange. Stable isotopes in CO2 and H2 O vapour provide mechanistic insight into physiological and biochemical processes during leaf gas exchange. We conducted combined leaf gas exchange and online isotope discrimination measurements on four common European tree species across a leaf temperature range of 5-40°C, while maintaining a constant leaf-to-air VPD (0.8 kPa) without soil water limitation. Above the optimum temperature for photosynthesis (30°C) under the controlled environmental conditions, stomatal conductance (gs ) and net photosynthesis rate (An ) decoupled across all tested species, with gs increasing but An decreasing. During this decoupling, mesophyll conductance (cell wall, plasma membrane and chloroplast membrane conductance) consistently and significantly decreased among species; however, this reduction did not lead to reductions in CO2 concentration at the chloroplast surface and stroma. We question the conventional understanding that diffusional limitations of CO2 contribute to the reduction in photosynthesis at high temperatures. We suggest that stomata and mesophyll membranes could work strategically to facilitate transpiration cooling and CO2 supply, thus alleviating heat stress on leaf photosynthetic function, albeit at the cost of reduced water-use efficiency.

The determining factors of hydrogen isotope offsets between plants and their source waters.

A fundamental assumption when using hydrogen and oxygen stable isotopes to understand ecohydrological processes is that no isotope fractionation occurs during plant water uptake/transport/redistribution. A growing body of evidence has indicated that hydrogen isotope fractionation occurs in certain environments or for certain plant species. However, whether the plant water source hydrogen isotope offset (δ2 H offset) is a common phenomenon and how it varies among different climates and plant functional types remains unclear. Here, we demonstrated the presence of positive, negative, and zero offsets based on extensive observations of 12 plant species of 635 paired stable isotopic compositions along a strong climate gradient within an inland river basin. Both temperature and relative humidity affected δ2 H offsets. In cool and moist environments, temperature mainly affected δ2 H offsets negatively due to its role in physiological activity. In warm and dry environments, relative humidity mainly affected δ2 H offsets, likely by impacting plant leaf stomatal conductance. These δ2 H offsets also showed substantial linkages with leaf water 18 O enrichment, an indicator of transpiration and evaporative demand. Further studies focusing on the ecophysiological and biochemical understanding of plant δ2 H dynamics under specific environments are essential for understanding regional ecohydrological processes and for conducting paleoclimate reconstructions.

Hydrogen isotope fractionation in plants with C3, C4, and CAM CO2 fixation.

Measurements of stable isotope ratios in organic compounds are widely used tools for plant ecophysiological studies. However, the complexity of the processes involved in shaping hydrogen isotope values (δ2H) in plant carbohydrates has limited its broader application. To investigate the underlying biochemical processes responsible for 2H fractionation among water, sugars, and cellulose in leaves, we studied the three main CO2 fixation pathways (C3, C4, and CAM) and their response to changes in temperature and vapor pressure deficit (VPD). We show significant differences in autotrophic 2H fractionation (εA) from water to sugar among the pathways and their response to changes in air temperature and VPD. The strong 2H depleting εA in C3 plants is likely driven by the photosynthetic H+ production within the thylakoids, a reaction that is spatially separated in C4 and strongly reduced in CAM plants, leading to the absence of 2H depletion in the latter two types. By contrast, we found that the heterotrophic 2H-fractionation (εH) from sugar to cellulose was very similar among the three pathways and is likely driven by the plant's metabolism, rather than by isotopic exchange with leaf water. Our study offers new insights into the biochemical drivers of 2H fractionation in plant carbohydrates.

Shifts in internal stem damage along a tropical precipitation gradient and implications for forest biomass estimation.

Woody biomass is a large carbon store in terrestrial ecosystems. In calculating biomass, tree stems are assumed to be solid structures. However, decomposer agents such as microbes and insects target stem heartwood, causing internal wood decay which is poorly quantified. We investigated internal stem damage across five sites in tropical Australia along a precipitation gradient. We estimated the amount of internal aboveground biomass damaged in living trees and measured four potential stem damage predictors: wood density, stem diameter, annual precipitation, and termite pressure (measured as termite damage in downed deadwood). Stem damage increased with increasing diameter, wood density, and termite pressure and decreased with increasing precipitation. High wood density stems sustained less damage in wet sites and more damage in dry sites, likely a result of shifting decomposer communities and their differing responses to changes in tree species and wood traits across sites. Incorporating stem damage reduced aboveground biomass estimates by > 30% in Australian savannas, compared to only 3% in rainforests. Accurate estimates of carbon storage across woody plant communities are critical for understanding the global carbon budget. Future biomass estimates should consider stem damage in concert with the effects of changes in decomposer communities and abiotic conditions.

Unsaturation in the air spaces of leaves and its implications.

Modern plant physiological theory stipulates that the resistance to water movement from plants to the atmosphere is overwhelmingly dominated by stomata. This conception necessitates a corollary assumption-that the air spaces in leaves must be nearly saturated with water vapour; that is, with a relative humidity that does not decline materially below unity. As this idea became progressively engrained in scientific discourse and textbooks over the last century, observations inconsistent with this corollary assumption were occasionally reported. Yet, evidence of unsaturation gained little traction, with acceptance of the prevailing framework motivated by three considerations: (1) leaf water potentials measured by either thermocouple psychrometry or the Scholander pressure chamber are largely consistent with the framework; (2) being able to assume near saturation of intercellular air spaces was transformational to leaf gas exchange analysis; and (3) there has been no obvious mechanism to explain a variable, liquid-phase resistance in the leaf mesophyll. Here, we review the evidence that refutes the assumption of universal, near saturation of air spaces in leaves. Refining the prevailing paradigm with respect to this assumption provides opportunities for identifying and developing mechanisms for increased plant productivity in the face of increasing evaporative demand imposed by global climate change.

Ecotypic Variation in Leaf Thermoregulation and Heat Tolerance but Not Thermal Safety Margins in Tropical Trees.

To avoid reaching lethal temperatures during periods of heat stress, plants may acclimate either their biochemical thermal tolerance or leaf morphological and physiological characteristics to reduce leaf temperature (Tleaf). While plants from warmer environments may have a greater capacity to regulate Tleaf, the extent of intraspecific variation and contribution of provenance is relatively unexplored. We tested whether upland and lowland provenances of four tropical tree species grown in a common garden differed in their thermal safety margins by measuring leaf thermal traits, midday leaf-to-air temperature differences (∆Tleaf) and critical leaf temperatures defined by chlorophyll fluorescence (Tcrit). Provenance variation was species- and trait-specific. Higher ∆Tleaf and Tcrit were observed in the lowland provenance for Terminalia microcarpa, and in the upland provenance for Castanospermum australe, with no provenance effects in the other two species. Within-species covariation of Tcrit and ∆Tleaf led to a convergence of thermal safety margins across provenances. While future studies should expand the number of provenances and species investigated, our findings suggest that lowland and upland provenances may not differ substantially in their vulnerability to heat stress, as determined by thermal safety margins, despite differences in operating temperatures and Tcrit.

Assessing the CO2 concentration at the surface of photosynthetic mesophyll cells.

We present a robust estimation of the CO2 concentration at the surface of photosynthetic mesophyll cells (cw ), applicable under reasonable assumptions of assimilation distribution within the leaf. We used Capsicum annuum, Helianthus annuus and Gossypium hirsutumas model plants for our experiments. We introduce calculations to estimate cw using independent adaxial and abaxial gas exchange measurements, and accounting for the mesophyll airspace resistances. The cw was lower than adaxial and abaxial estimated intercellular CO2 concentrations (ci ). Differences between cw and the ci of each surface were usually larger than 10 µmol mol-1 . Differences between adaxial and abaxial ci ranged from a few µmol mol-1 to almost 50 µmol mol-1 , where the largest differences were found at high air saturation deficits (ASD). Differences between adaxial and abaxial ci and the ci estimated by mixing both fluxes ranged from -30 to +20 µmol mol-1 , where the largest differences were found under high ASD or high ambient CO2 concentrations. Accounting for cw improves the information that can be extracted from gas exchange experiments, allowing a more detailed description of the CO2 and water vapor gradients within the leaf.

Covariation between oxygen and hydrogen stable isotopes declines along the path from xylem water to wood cellulose across an aridity gradient.

Oxygen and hydrogen isotopes of cellulose in plant biology are commonly used to infer environmental conditions, often from time series measurements of tree rings. However, the covariation (or the lack thereof) between δ18 O and δ2 H in plant cellulose is still poorly understood. We compared plant water, and leaf and branch cellulose from dominant tree species across an aridity gradient in Northern Australia, to examine how δ18 O and δ2 H relate to each other and to mean annual precipitation (MAP). We identified a decline in covariation from xylem to leaf water, and onwards from leaf to branch wood cellulose. Covariation in leaf water isotopic enrichment (Δ) was partially preserved in leaf cellulose but not branch wood cellulose. Furthermore, whilst δ2 H was well-correlated between leaf and branch, there was an offset in δ18 O between organs that increased with decreasing MAP. Our findings strongly suggest that postphotosynthetic isotope exchange with water is more apparent for oxygen isotopes, whereas variable kinetic and nonequilibrium isotope effects add complexity to interpreting metabolic-induced δ2 H patterns. Varying oxygen isotope exchange in wood and leaf cellulose must be accounted for when δ18 O is used to reconstruct climatic scenarios. Conversely, comparing δ2 H and δ18 O patterns may reveal environmentally induced shifts in metabolism.

Unenriched xylem water contribution during cellulose synthesis influenced by atmospheric demand governs the intra-annual tree-ring δ18 O signature.

The oxygen isotope composition (δ18 O) of tree-ring cellulose is used to evaluate tree physiological responses to climate, but their interpretation is still limited due to the complexity of the isotope fractionation pathways. We assessed the relative contribution of seasonal needle and xylem water δ18 O variations to the intra-annual tree-ring cellulose δ18 O signature of larch trees at two sites with contrasting soil water availability in the Swiss Alps. We combined biweekly δ18 O measurements of soil water, needle water, and twig xylem water with intra-annual δ18 O measurements of tree-ring cellulose, xylogenesis analysis, and mechanistic and structural equation modeling. Intra-annual cellulose δ18 O values resembled source water δ18 O mean levels better than needle water δ18 O. Large parts of the rings were formed under high proportional exchange with unenriched xylem water (pex ). Maximum pex values were achieved in August and imprinted on sections at 50-75% of the ring. High pex values were associated with periods of high atmospheric evaporative demand (VPD). While VPD governed needle water δ18 O variability, we estimated a limited Péclet effect at both sites. Due to a variable pex , source water has a strong influence over large parts of the intra-annual tree-ring cellulose δ18 O variations, potentially masking signals coming from needle-level processes.

Vapour pressure deficit modulates hydraulic function and structure of tropical rainforests under nonlimiting soil water supply.

Atmospheric conditions are expected to become warmer and drier in the future, but little is known about how evaporative demand influences forest structure and function independently from soil moisture availability, and how fast-response variables (such as canopy water potential and stomatal conductance) may mediate longer-term changes in forest structure and function in response to climate change. We used two tropical rainforest sites with different temperatures and vapour pressure deficits (VPD), but nonlimiting soil water supply, to assess the impact of evaporative demand on ecophysiological function and forest structure. Common species between sites allowed us to test the extent to which species composition, relative abundance and intraspecific variability contributed to site-level differences. The highest VPD site had lower midday canopy water potentials, canopy conductance (gc ), annual transpiration, forest stature, and biomass, while the transpiration rate was less sensitive to changes in VPD; it also had different height-diameter allometry (accounting for 51% of the difference in biomass between sites) and higher plot-level wood density. Our findings suggest that increases in VPD, even in the absence of soil water limitation, influence fast-response variables, such as canopy water potentials and gc , potentially leading to longer-term changes in forest stature resulting in reductions in biomass.

Coordination of photosynthetic traits across soil and climate gradients.

"Least-cost theory" posits that C3 plants should balance rates of photosynthetic water loss and carboxylation in relation to the relative acquisition and maintenance costs of resources required for these activities. Here we investigated the dependency of photosynthetic traits on climate and soil properties using a new Australia-wide trait dataset spanning 528 species from 67 sites. We tested the hypotheses that plants on relatively cold or dry sites, or on relatively more fertile sites, would typically operate at greater CO2 drawdown (lower ratio of leaf internal to ambient CO2 , Ci :Ca ) during light-saturated photosynthesis, and at higher leaf N per area (Narea ) and higher carboxylation capacity (Vcmax 25 ) for a given rate of stomatal conductance to water vapour, gsw . These results would be indicative of plants having relatively higher water costs than nutrient costs. In general, our hypotheses were supported. Soil total phosphorus (P) concentration and (more weakly) soil pH exerted positive effects on the Narea -gsw and Vcmax 25 -gsw slopes, and negative effects on Ci :Ca . The P effect strengthened when the effect of climate was removed via partial regression. We observed similar trends with increasing soil cation exchange capacity and clay content, which affect soil nutrient availability, and found that soil properties explained similar amounts of variation in the focal traits as climate did. Although climate typically explained more trait variation than soil did, together they explained up to 52% of variation in the slope relationships and soil properties explained up to 30% of the variation in individual traits. Soils influenced photosynthetic traits as well as their coordination. In particular, the influence of soil P likely reflects the Australia's geologically ancient low-relief landscapes with highly leached soils. Least-cost theory provides a valuable framework for understanding trade-offs between resource costs and use in plants, including limiting soil nutrients.

Forest system hydraulic conductance: partitioning tree and soil components.

Soil-leaf hydraulic conductance determines canopy-atmosphere coupling in vegetation models, but it is typically derived from ex-situ measurements of stem segments and soil samples. Using a novel approach, we derive robust in-situ estimates for whole-tree conductance (ktree ), 'functional' soil conductance (ksoil ), and 'system' conductance (ksystem , water table to canopy), at two climatically different tropical rainforest sites. Hydraulic 'functional rooting depth', determined for each tree using profiles of soil water potential (Ψsoil ) and sap flux data, enabled a robust determination of ktree and ksoil . ktree was compared across species, size classes, seasons, height above nearest drainage (HAND), two field sites, and to alternative representations of ktree ; ksoil was analysed with respect to variations in site, season and HAND. ktree was lower and changed seasonally at the site with higher vapour pressure deficit (VPD) and rainfall; ktree differed little across species but scaled with tree circumference; rsoil (1/ksoil ) ranged from 0 in the wet season to 10× less than rtree (1/ktree ) in the dry season. VPD and not rainfall may influence plot-level k; leaf water potentials and sap flux can be used to determine ktree , ksoil and ksystem ; Ψsoil profiles can provide mechanistic insights into ecosystem-level water fluxes.

A meta-analysis of responses of C3 plants to atmospheric CO2 : dose-response curves for 85 traits ranging from the molecular to the whole-plant level.

Generalised dose-response curves are essential to understand how plants acclimate to atmospheric CO2 . We carried out a meta-analysis of 630 experiments in which C3 plants were experimentally grown at different [CO2 ] under relatively benign conditions, and derived dose-response curves for 85 phenotypic traits. These curves were characterised by form, plasticity, consistency and reliability. Considered over a range of 200-1200 µmol mol-1 CO2 , some traits more than doubled (e.g. area-based photosynthesis; intrinsic water-use efficiency), whereas others more than halved (area-based transpiration). At current atmospheric [CO2 ], 64% of the total stimulation in biomass over the 200-1200 µmol mol-1 range has already been realised. We also mapped the trait responses of plants to [CO2 ] against those we have quantified before for light intensity. For most traits, CO2 and light responses were of similar direction. However, some traits (such as reproductive effort) only responded to light, others (such as plant height) only to [CO2 ], and some traits (such as area-based transpiration) responded in opposite directions. This synthesis provides a comprehensive picture of plant responses to [CO2 ] at different integration levels and offers the quantitative dose-response curves that can be used to improve global change simulation models.

Towards species-level forecasts of drought-induced tree mortality risk.

Predicting species-level responses to drought at the landscape scale is critical to reducing uncertainty in future terrestrial carbon and water cycle projections. We embedded a stomatal optimisation model in the Community Atmosphere Biosphere Land Exchange (CABLE) land surface model and parameterised the model for 15 canopy dominant eucalypt tree species across South-Eastern Australia (mean annual precipitation range: 344-1424 mm yr-1 ). We conducted three experiments: applying CABLE to the 2017-2019 drought; a 20% drier drought; and a 20% drier drought with a doubling of atmospheric carbon dioxide (CO2 ). The severity of the drought was highlighted as for at least 25% of their distribution ranges, 60% of species experienced leaf water potentials beyond the water potential at which 50% of hydraulic conductivity is lost due to embolism. We identified areas of severe hydraulic stress within-species' ranges, but we also pinpointed resilience in species found in predominantly semiarid areas. The importance of the role of CO2 in ameliorating drought stress was consistent across species. Our results represent an important advance in our capacity to forecast the resilience of individual tree species, providing an evidence base for decision-making around the resilience of restoration plantings or net-zero emission strategies.

Do 2 H and 18 O in leaf water reflect environmental drivers differently?

We compiled hydrogen and oxygen stable isotope compositions (δ2 H and δ18 O) of leaf water from multiple biomes to examine variations with environmental drivers. Leaf water δ2 H was more closely correlated with δ2 H of xylem water or atmospheric vapour, whereas leaf water δ18 O was more closely correlated with air relative humidity. This resulted from the larger proportional range for δ2 H of meteoric waters relative to the extent of leaf water evaporative enrichment compared with δ18 O. We next expressed leaf water as isotopic enrichment above xylem water (Δ2 H and Δ18 O) to remove the impact of xylem water isotopic variation. For Δ2 H, leaf water still correlated with atmospheric vapour, whereas Δ18 O showed no such correlation. This was explained by covariance between air relative humidity and the Δ18 O of atmospheric vapour. This is consistent with a previously observed diurnal correlation between air relative humidity and the deuterium excess of atmospheric vapour across a range of ecosystems. We conclude that 2 H and 18 O in leaf water do indeed reflect the balance of environmental drivers differently; our results have implications for understanding isotopic effects associated with water cycling in terrestrial ecosystems and for inferring environmental change from isotopic biomarkers that act as proxies for leaf water.