Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla.

BACKGROUND: Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla. RESULTS: We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions. CONCLUSIONS: Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Pushing the envelope: do narrowly and widely distributed Eucalyptus species differ in response to climate warming?

Contemporary climate change will push many tree species into conditions that are outside their current climate envelopes. Using the Eucalyptus genus as a model, we addressed whether species with narrower geographical distributions show constrained ability to cope with warming relative to species with wider distributions, and whether this ability differs among species from tropical and temperate climates. We grew seedlings of widely and narrowly distributed Eucalyptus species from temperate and tropical Australia in a glasshouse under two temperature regimes: the summer temperature at seed origin and +3.5°C. We measured physical traits and leaf-level gas exchange to assess warming influences on growth rates, allocation patterns, and physiological acclimation capacity. Warming generally stimulated growth, such that higher relative growth rates early in development placed seedlings on a trajectory of greater mass accumulation. The growth enhancement under warming was larger among widely than narrowly distributed species and among temperate rather than tropical provenances. The differential growth enhancement was primarily attributable to leaf area production and adjustments of specific leaf area. Our results suggest that tree species, including those with climate envelopes that will be exceeded by contemporary climate warming, possess capacity to physiologically acclimate but may have varying ability to adjust morphology.

The importance of species-specific and temperature-sensitive parameterisation of A/Ci models: A case study using cotton (Gossypium hirsutum L.) and the automated 'OptiFitACi' R-package.

Leaf gas exchange measurements are an important tool for inferring a plant's photosynthetic biochemistry. In most cases, the responses of photosynthetic CO2 assimilation to variable intercellular CO2 concentrations (A/Ci response curves) are used to model the maximum (potential) rate of carboxylation by ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco, Vcmax) and the rate of photosynthetic electron transport at a given incident photosynthetically active radiation flux density (PAR; JPAR). The standard Farquhar-von Caemmerer-Berry model is often used with default parameters of Rubisco kinetic values and mesophyll conductance to CO2 (gm) derived from tobacco that may be inapplicable across species. To study the significance of using such parameters for other species, here we measured the temperature responses of key in vitro Rubisco catalytic properties and gm in cotton (Gossypium hirsutum cv. Sicot 71) and derived Vcmax and J2000 (JPAR at 2000 µmol m-2 s-1 PAR) from cotton A/Ci curves incrementally measured at 15°C-40°C using cotton and other species-specific sets of input parameters with our new automated fitting R package 'OptiFitACi'. Notably, parameterisation by a set of tobacco parameters produced unrealistic J2000:Vcmax ratio of <1 at 25°C, two- to three-fold higher estimates of Vcmax above 15°C, up to 2.3-fold higher estimates of J2000 and more variable estimates of Vcmax and J2000, for our cotton data compared to model parameterisation with cotton-derived values. We determined that errors arise when using a gm,25 of 2.3 mol m-2 s-1 MPa-1 or less and Rubisco CO2-affinities in 21% O2 (KC 21%O2) at 25°C outside the range of 46-63 Pa to model A/Ci responses in cotton. We show how the A/Ci modelling capabilities of 'OptiFitACi' serves as a robust, user-friendly, and flexible extension of 'plantecophys' by providing simplified temperature-sensitivity and species-specificity parameterisation capabilities to reduce variability when modelling Vcmax and J2000.

The double-edged sword of potassium and sodium fertilization in xylem embolism resistance of two Eucalyptus species under drought stress.

Sodium (Na+) is a beneficial element for most plants and may replace potassium (K+) in osmoregulatory process to a certain extent, increasing plant water use efficiency. Thus, understanding coordinated mechanisms underlying the combined use of K+ and Na+ in tree drought tolerance is a key challenge for forestry in dealing with productivity and water limitations. A pot experiment with three ratios of K/Na (K-supplied, partial K replacement by Na, and K-deficient plants) and two water regimes, well-watered (W+) and water-stressed (W-), was conducted on saplings of two Eucalyptus species with contrasting drought sensitivities. We evaluated the point of stomatal closure (Pgs90), xylem water potential at 12, 50, and 88% embolized xylem area (P12, P50, P88), hydraulic safety margin, leaf gas exchange (A, E, gs, and dark respiration), pre-dawn and midday leaf water potential (ΨPD and ΨMD), long-term water use efficiency (WUEL) and total dry mass. Partial K replacement by Na increased leaf gas exchange, WUEL, and total dry mass, while Pgs90, P12, P50, P88, and ΨMD decreased (were more negative), compared with plants exclusively supplied with K and K-deficient plants of both species. Fertilized plants had narrower hydraulic safety margins than K-deficient plants, indicating that these Eucalyptus species adopt the functional adaptive strategy of operating close to their hydraulic limits to maximize carbon uptake while increasing the risk of hydraulic failure under drought stress.

Does the response of Rubisco and photosynthesis to elevated [CO2] change with unfavourable environmental conditions?

Climate change due to anthropogenic CO2 emissions affects plant performance globally. To improve crop resilience, we need to understand the effects of elevated CO2 concentration (e[CO2]) on CO2 assimilation and Rubisco biochemistry. However, the interactive effects of e[CO2] and abiotic stress are especially unclear. This study analyses the CO2 effect on photosynthetic capacity under different water availability and temperature conditions in 42 different crop species, varying in functional group, photosynthetic pathway and phenological stage. We analysed close to 3000 data points extracted from 120 published manuscripts. For C3 species, e[CO2] increases net photosynthesis and intercellular [CO2], while reducing stomatal conductance and transpiration. Vmaxc, Rubisco in vitro extractable maximal activity and content also decrease with e[CO2] in C3 species, while C4 crops are less responsive to e[CO2]. The interaction with drought and/or heat stress does not significantly alter these photosynthetic responses, indicating that the photosynthetic capacity of stressed plants responds to e[CO2]. Moreover, e[CO2] has strong effect on the photosynthetic capacity of grasses mainly in the final stages of development. This study provides insight into the intricate interactions within the plant photosynthetic apparatus under the influence of climate change, enhancing the understanding of mechanisms governing plant responses to environmental parameters.

Long-term warming increased carbon sequestration capacity in a humid subtropical forest.

Tropical and subtropical forests play a crucial role in global carbon (C) pools, and their responses to warming can significantly impact C-climate feedback and predictions of future global warming. Despite earth system models projecting reductions in land C storage with warming, the magnitude of this response varies greatly between models, particularly in tropical and subtropical regions. Here, we conducted a field ecosystem-level warming experiment in a subtropical forest in southern China, by translocating mesocosms (ecosystem composed of soils and plants) across 600 m elevation gradients with temperature gradients of 2.1°C (moderate warming), to explore the response of ecosystem C dynamics of the subtropical forest to continuous 6-year warming. Compared with the control, the ecosystem C stock decreased by 3.8% under the first year of 2.1°C warming; but increased by 13.4% by the sixth year of 2.1°C warming. The increased ecosystem C stock by the sixth year of warming was mainly attributed to a combination of sustained increased plant C stock due to the maintenance of a high plant growth rate and unchanged soil C stock. The unchanged soil C stock was driven by compensating and offsetting thermal adaptation of soil microorganisms (unresponsive soil respiration and enzyme activity, and more stable microbial community), increased plant C input, and inhibitory C loss (decreased C leaching and inhibited temperature sensitivity of soil respiration) from soil drying. These results suggest that the humid subtropical forest C pool would not necessarily diminish consistently under future long-term warming. We highlight that differential and asynchronous responses of plant and soil C processes over relatively long-term periods should be considered when predicting the effects of climate warming on ecosystem C dynamics of subtropical forests.

Root relative water content is a potential signal for impending mortality of a subtropical conifer during extreme drought stress.

Adaptation to future climates characterized by more frequent severe droughts requires enhanced mechanistic understanding of tree mortality. However, our knowledge of the physiological limits to withstand extreme drought, and how the coordination between water and carbon traits enhances survival, is still limited. Potted seedlings of Pinus massoniana were dehydrated to three target droughts (percentage loss of stem hydraulic conductivity of ca. 50%, 85%, and 100%; PLC50 , PLC85 and PLC100 ) and then relieved from these target droughts by fully rewatering. Predawn and midday water potentials (Ψ), relative water content (RWC), PLC and nonstructural carbohydrates (NSC) were monitored. During drought, Ψ and RWC declined as PLC increased. Root RWC declined more rapidly than other organ RWCs, particularly after PLC50 stress. All organ NSC concentrations were above predrought values. During rewatering, water trait recovery declined as drought increased, with no mortality at PLC50 but 75% mortality at PLC85 . The observed stem hydraulic recovery at PLC50 following rewatering was not correlated to NSC dynamics. Collectively, our results highlighted the primary role of hydraulic failure in Pinus massoniana seedling mortality by assessing mortality threshold and links among water status and water supply. Root RWC can be considered as a potential warning signal of P. massoniana mortality.

Response of the plant core microbiome to Fusarium oxysporum infection and identification of the pathobiome.

Plant core microbiomes consist of persistent key members that provide critical host functions, but their assemblages can be interrupted by biotic and abiotic stresses. The pathobiome is comprised of dynamic microbial interactions in response to disease status of the host. Hence, identifying variation in the core microbiome and pathobiome can significantly advance our understanding of microbial-microbial interactions and consequences for disease progression and host functions. In this study, we combined glasshouse and field studies to analyse the soil and plant rhizosphere microbiome of cotton plants (Gossypium hirsutum) in the presence of a cotton-specific fungal pathogen, Fusarium oxysporum f. sp. vasinfectum (FOV). We found that FOV directly and consistently altered the rhizosphere microbiome, but the biocontrol agents enabled microbial assemblages to resist pathogenic stress. Using co-occurrence network analysis of the core microbiome, we identified the pathobiome comprised of the pathogen and key associate phylotypes in the cotton microbiome. Isolation and application of some negatively correlated pathobiome members provided protection against plant infection. Importantly, our field survey from multiple cotton fields validated the pattern and responses of core microbiomes under FOV infection. This study advances key understanding of core microbiome responses and existence of plant pathobiomes, which provides a novel framework to better manage plant diseases in agriculture and natural settings.

Mechanical stress acclimation in plants: Linking hormones and somatic memory to thigmomorphogenesis.

A single event of mechanical stimulation is perceived by mechanoreceptors that transduce rapid transient signalling to regulate gene expression. Prolonged mechanical stress for days to weeks culminates in cellular changes that strengthen the plant architecture leading to thigmomorphogenesis. The convergence of multiple signalling pathways regulates mechanically induced tolerance to numerous biotic and abiotic stresses. Emerging evidence showed prolonged mechanical stimulation can modify the baseline level of gene expression in naive tissues, heighten gene expression, and prime disease resistance upon a subsequent pathogen encounter. The phenotypes of thigmomorphogenesis can persist throughout growth without continued stimulation, revealing somatic-stress memory. Epigenetic processes regulate TOUCH gene expression and could program transcriptional memory in differentiating cells to program thigmomorphogenesis. We discuss the early perception, gene regulatory and phytohormone pathways that facilitate thigmomorphogenesis and mechanical stress acclimation in Arabidopsis and other plant species. We provide insights regarding: (1) the regulatory mechanisms induced by single or prolonged events of mechanical stress, (2) how mechanical stress confers transcriptional memory to induce cross-acclimation to future stress, and (3) why thigmomorphogenesis might resemble an epigenetic phenomenon. Deeper knowledge of how prolonged mechanical stimulation programs somatic memory and primes defence acclimation could transform solutions to improve agricultural sustainability in stressful environments.

The roles of divergent and parallel molecular evolution contributing to thermal adaptive strategies in trees.

Local adaptation is a driver of biological diversity, and species may develop analogous (parallel evolution) or alternative (divergent evolution) solutions to similar ecological challenges. We expect these adaptive solutions would culminate in both phenotypic and genotypic signals. Using two Eucalyptus species (Eucalyptus grandis and Eucalyptus tereticornis) with overlapping distributions grown under contrasting 'local' temperature conditions to investigate the independent contribution of adaptation and plasticity at molecular, physiological and morphological levels. The link between gene expression and traits markedly differed between species. Divergent evolution was the dominant pattern driving adaptation (91% of all significant genes); but overlapping gene (homologous) responses were dependent on the determining factor (plastic, adaptive or genotype by environment interaction). Ninety-eight percent of the plastic homologs were similarly regulated, while 50% of the adaptive homologs and 100% of the interaction homologs were antagonistical. Parallel evolution for the adaptive effect in homologous genes was greater than expected but not in favour of divergent evolution. Heat shock proteins for E. grandis were almost entirely driven by adaptation, and plasticity in E. tereticornis. These results suggest divergent molecular evolutionary solutions dominated the adaptive mechanisms among species, even in similar ecological circumstances. Suggesting that tree species with overlapping distributions are unlikely to equally persist in the future.

Testing the limits of plant drought stress and subsequent recovery in four provenances of a widely distributed subtropical tree species.

Drought-induced tree mortality may increase with ongoing climate change. Unraveling the links between stem hydraulics and mortality thresholds, and the effects of intraspecific variation, remain important unresolved issues. We conducted a water manipulation experiment in a rain-out shelter, using four provenances of Schima superba originating from a gradient of annual precipitation (1124-1796 mm) and temperature (16.4-22.4°C). Seedlings were droughted to three levels of percentage loss of hydraulic conductivity (i.e., P50 , P88  and P99) and subsequently rewatered to field capacity for 30 days; traits related to water and carbon relations were measured. The lethal water potential associated with incipient mortality was between P50 and P88 . Seedlings exhibited similar drought responses in xylem water potential, hydraulic conductivity and gas exchange. Upon rehydration, patterns of gas exchange differed among provenances but were not related to the climate at the origin. The four provenances exhibited a similar degree of stem hydraulic recovery, which was correlated with the magnitude of antecedent drought and stem soluble sugar at the end of the drought. Results suggest that there were intraspecific differences in the capacity of S. superba seedlings for carbon assimilation during recovery, indicating a decoupling between gas exchange recovery and stem hydraulics across provenances.

Mechanisms of xylem hydraulic recovery after drought in Eucalyptus saligna.

The mechanisms by which woody plants recover xylem hydraulic capacity after drought stress are not well understood, particularly with regard to the role of embolism refilling. We evaluated the recovery of xylem hydraulic capacity in young Eucalyptus saligna plants exposed to cycles of drought stress and rewatering. Plants were exposed to moderate and severe drought stress treatments, with recovery monitored at time intervals from 24 h to 6 months after rewatering. The percentage loss of xylem vessels due to embolism (PLV) was quantified at each time point using microcomputed tomography with stem water potential (Ψx ) and canopy transpiration (Ec ) measured before scans. Plants exposed to severe drought stress suffered high levels of embolism (47.38% ± 10.97% PLV) and almost complete canopy loss. No evidence of embolism refilling was observed at 24 h, 1 week, or 3 weeks after rewatering despite rapid recovery in Ψx . Recovery of hydraulic capacity was achieved over a 6-month period by growth of new xylem tissue, with canopy leaf area and Ec recovering over the same period. These findings indicate that E. saligna recovers slowly from severe drought stress, with potential for embolism to persist in the xylem for many months after rainfall events.

High safety margins to drought-induced hydraulic failure found in five pasture grasses.

Determining the relationship between reductions in stomatal conductance (gs ) and leaf water transport during dehydration is key to understanding plant drought responses. While numerous studies have analysed the hydraulic function of woody species, minimal research has been conducted on grasses. Here, we sought to characterize hydraulic vulnerability in five widely-occurring pasture grasses (including both C3 and C4 grasses) and determine whether reductions in gs and leaf hydraulic conductance (Kleaf ) during dehydration could be attributed to xylem embolism. Using the optical vulnerability (OV) technique, we found that all species were highly resistant to xylem embolism when compared to other herbaceous angiosperms, with 50% xylem embolism (PX50 ) occurring at xylem pressures ranging from -4.4 to -6.1 MPa. We observed similar reductions in gs and Kleaf under mild water stress for all species, occurring well before PX50 . The onset of xylem embolism (PX12 ) occurred consistently after stomatal closure and 90% reduction of Kleaf . Our results suggest that factors other than xylem embolism are responsible for the majority of reductions in gs and Kleaf during drought and reductions in the productivity of pasture species under moderate drought may not be driven by embolism.

Warming drives sustained plant phosphorus demand in a humid tropical forest.

Phosphorus (P) is often one of the most limiting nutrients in highly weathered soils of humid tropical forests and may regulate the responses of carbon (C) feedback to climate warming. However, the response of P to warming at the ecosystem level in tropical forests is not well understood because previous studies have not comprehensively assessed changes in multiple P processes associated with warming. Here, we detected changes in the ecosystem P cycle in response to a 7-year continuous warming experiment by translocating model plant-soil ecosystems across a 600-m elevation gradient, equivalent to a temperature change of 2.1°C. We found that warming increased plant P content (55.4%) and decreased foliar N:P. Increased plant P content was supplied by multiple processes, including enhanced plant P resorption (9.7%), soil P mineralization (15.5% decrease in moderately available organic P), and dissolution (6.8% decrease in iron-bound inorganic P), without changing litter P mineralization and leachate P. These findings suggest that warming sustained plant P demand by increasing the biological and geochemical controls of the plant-soil P-cycle, which has important implications for C fixation in P-deficient and highly productive tropical forests in future warmer climates.

To what extent can rising [CO2 ] ameliorate plant drought stress?

Plant responses to elevated atmospheric carbon dioxide (eCO2 ) have been hypothesized as a key mechanism that may ameliorate the impact of future drought. Yet, despite decades of experiments, the question of whether eCO2 reduces plant water use, yielding 'water savings' that can be used to maintain plant function during periods of water stress, remains unresolved. In this Viewpoint, we identify the experimental challenges and limitations to our understanding of plant responses to drought under eCO2 . In particular, we argue that future studies need to move beyond exploring whether eCO2 played 'a role' or 'no role' in responses to drought, but instead more carefully consider the timescales and conditions that would induce an influence. We also argue that considering emergent differences in soil water content may be an insufficient means of assessing the impact of eCO2 . We identify eCO2 impact during severe drought (e.g. to the point of mortality), interactions with future changes in vapour pressure deficit and uncertainty about changes in leaf area as key gaps in our current understanding. New insights into CO2 × drought interactions are essential to better constrain model theory that governs future climate model projections of land-atmosphere interactions during periods of water stress.

Drought by CO2 interactions in trees: a test of the water savings mechanism.

Elevated atmospheric CO2 (eCa ) may benefit plants during drought by reducing stomatal conductance (gs ) but any 'water savings effect' could be neutralized by concurrent stimulation of leaf area. We investigated whether eCa enhanced water savings, thereby ameliorating the impact of drought on carbon and water relations in trees. We report leaf-level gas exchange and whole-plant and soil water relations during a short-term dry-down in two Eucalyptus species with contrasting drought tolerance. Plants had previously been established for 9 to 11 months in steady-state conditions of ambient atmospheric CO2 (aCa ) and eCa , with half of each treatment group exposed to sustained drought for 5 to 7 months. The lower stomatal conductance under eCa did not lead to soil moisture savings during the dry-down due to the counteractive effect of increased whole-plant leaf area. Nonetheless, eCa -grown plants maintained higher photosynthetic rates and leaf water potentials, making them less stressed during the dry-down, despite being larger. These effects were more pronounced in the xeric species than the mesic species, and in previously water-stressed plants. Our findings indicate that eCa may enhance plant performance during drought despite a lack of soil water savings, especially in species with more conservative growth and water-use strategies.

Climate and stomatal traits drive covariation in nighttime stomatal conductance and daytime gas exchange rates in a widespread C4 grass.

Nighttime stomatal conductance (gsn ) varies among plant functional types and species, but factors shaping the evolution of gsn remain unclear. Examinations of intraspecific variation in gsn as a function of climate and co-varying leaf traits may provide new insight into the evolution of gsn and its adaptive significance. We grew 11 genotypes of Panicum virgatum (switchgrass) representing differing home-climates in a common garden experiment and measured nighttime and daytime leaf gas exchange, as well as stomatal density (SD) and size during early-, mid-, and late-summer. We used piecewise structural equation modelling to determine direct and indirect relationships between home-climate, gas exchange, and stomatal traits. We found no direct relationship between home-climate and gsn . However, genotypes from hotter climates possessed higher SD, which resulted in higher gsn . Across genotypes, higher gsn was associated with higher daytime stomatal conductance and net photosynthesis. Our results indicate that higher gsn may arise in genotypes from hotter climates via increased SD. High SD may provide benefits to genotypes from hotter climates through enhanced daytime transpirational cooling or by permitting maximal gas exchange when conditions are suitable. These results highlight the role of climate and trait coordination in shaping genetic differentiation in gsn .

Mesophyll conductance in two cultivars of wheat grown in glacial to super-elevated CO2 concentrations.

Mesophyll conductance (gm) is an important factor limiting photosynthesis. However, gm response to long-term growth in variable [CO2] is not well understood, particularly in crop plants. Here, we grew two cultivars of wheat (Halberd and Cranbrook), known to differ in gm under current environmental conditions, in four [CO2] treatments: glacial (206 µmol mol-1), pre-industrial (344 µmol mol-1), current ambient (489 µmol mol-1), and super-elevated (1085 µmol mol-1), and two water treatments (well-watered and moderate water limitation), to develop an evolutionary and future climate perspective on gm control of photosynthesis and water-use efficiency (WUE). In the two wheat genotypes, gm increased with rising [CO2] from glacial to ambient [CO2], but declined at super-elevated [CO2]. The responses of gm to different growth [CO2] also depend on water stress; however, the specific mechanism of gm response to [CO2] remains unclear. Although gm and gm/gsc (mesophyll conductance/stomatal conductance) were strongly associated with the variability of photosynthetic rates (A) and WUE, we found that plants with higher gm may increase A without increasing gsc, which increased WUE. These results may be useful to inform plant breeding programmes and cultivar selection for Australian wheat under future environmental conditions.

Smart glass impacts stomatal sensitivity of greenhouse Capsicum through altered light.

Optical films that alter light transmittance may reduce energy consumption in high-tech greenhouses, but their impact on crop physiology remains unclear. We compared the stomatal responses of Capsicum plants grown hydroponically under control glass (70% diffuse light) or the smart glass (SG) film ULR-80, which blocked >50% of short-wave radiation and ~9% of photosynthetically active radiation (PAR). SG had no significant effects on steady-state (gs) or maximal (gmax) stomatal conductance. In contrast, SG reduced stomatal pore size and sensitivity to exogenous abscisic acid (ABA), thereby increasing rates of leaf water loss, guard cell K+ and Cl- efflux, and Ca2+ influx. SG induced faster stomatal closing and opening rates on transition between low (100 µmol m-2 s-1) and high PAR (1500 µmol m-2 s-1), which compromised water use efficiency relative to control plants. The fraction of blue light (0% or 10%) did not affect gs in either treatment. Increased expression of stomatal closure and photoreceptor genes in epidermal peels of SG plants is consistent with fast stomatal responses to light changes. In conclusion, stomatal responses of Capsicum to SG were more affected by changes in light intensity than spectral quality, and re-engineering of the SG should maximize PAR transmission, and hence CO2 assimilation.

Increasing aridity will not offset CO2 fertilization in fast-growing eucalypts with access to deep soil water.

Rising atmospheric [CO2 ] (Ca ) generally enhances tree growth if nutrients are not limiting. However, reduced water availability and elevated evaporative demand may offset such fertilization. Trees with access to deep soil water may be able to mitigate such stresses and respond more positively to Ca . Here, we sought to evaluate how increased vapor pressure deficit and reduced precipitation are likely to modify the impact of elevated Ca (eCa ) on tree productivity in an Australian Eucalyptus saligna Sm. plantation with access to deep soil water. We parameterized a forest growth simulation model (GOTILWA+) using data from two field experiments on E. saligna: a 2-year whole-tree chamber experiment with factorial Ca (ambient =380, elevated =620 µmol mol-1 ) and watering treatments, and a 10-year stand-scale irrigation experiment. Model evaluation showed that GOTILWA+ can capture the responses of canopy C uptake to (1) rising vapor pressure deficit (D) under both Ca treatments; (2) alterations in tree water uptake from shallow and deep soil layers during soil dry-down; and (3) the impact of irrigation on tree growth. Simulations suggest that increasing Ca up to 700 µmol mol-1 alone would result in a 33% increase in annual gross primary production (GPP) and a 62% increase in biomass over 10 years. However, a combined 48% increase in D and a 20% reduction in precipitation would halve these values. Our simulations identify high D conditions as a key limiting factor for GPP. They also suggest that rising Ca will compensate for increasing aridity limitations in E. saligna trees with access to deep soil water under non-nutrient limiting conditions, thereby reducing the negative impacts of global warming upon this eucalypt species. Simulation models not accounting for water sources available to deep-rooting trees are likely to overestimate aridity impacts on forest productivity and C stocks.