Tree water uptake patterns across the globe

Plant water uptake from the soil is a crucial element of the global hydrological cycle and essential for vegetation drought resilience. Yet, knowledge of how the distribution of water uptake depth (WUD) varies across species, climates, and seasons is scarce relative to our knowledge of aboveground plant functions. With a global literature review, we found that average WUD varied more among biomes than plant functional types (i.e. deciduous/evergreen broadleaves and conifers), illustrating the importance of the hydroclimate, especially precipitation seasonality, on WUD. By combining records of rooting depth with WUD, we observed a consistently deeper maximum rooting depth than WUD with the largest differences in arid regions - indicating that deep taproots act as lifelines while not contributing to the majority of water uptake. The most ubiquitous observation across the literature was that woody plants switch water sources to soil layers with the highest water availability within short timescales. Hence, seasonal shifts to deep soil layers occur across the globe when shallow soils are drying out, allowing continued transpiration and hydraulic safety. While there are still significant gaps in our understanding of WUD, the consistency across global ecosystems allows integration of existing knowledge into the next generation of vegetation process models.

Combining PSII photochemistry and hydraulics improves predictions of photosynthesis and water use from mild to lethal drought

Rising temperatures and increases in drought negatively impact the efficiency and sustainability of both agricultural and forest ecosystems. Although hydraulic limitations on photosynthesis have been extensively studied, a solid understanding of the links between whole plant hydraulics and photosynthetic processes at the cellular level under changing environmental conditions is still missing, hampering our predictive power for plant mortality. Here, we examined plant hydraulic traits and CO2 assimilation rate under progressive water limitation by implementing Photosystem II (PSII) dynamics with a whole plant process model (TREES). The photosynthetic responses to plant water status were parameterized based on measurements of chlorophyll a fluorescence, gas exchange and water potential for Brassica rapa (R500) grown in a greenhouse under fully watered to lethal drought conditions. The updated model significantly improved predictions of photosynthesis, stomatal conductance and leaf water potential. TREES with PSII knowledge predicted a larger hydraulic safety margin and a decrease in percent loss of conductivity. TREES predicted a slower decrease in leaf water potential, which agreed with measurements. Our results highlight the pressing need for incorporating PSII drought photochemistry into current process models to capture cross-scale plant water dynamics from cell to whole plant level.

The influence of increasing atmospheric CO2 , temperature, and vapor pressure deficit on seawater-induced tree mortality

Increasing seawater exposure is killing coastal trees globally, with expectations of accelerating mortality with rising sea levels. However, the impact of concomitant changes in atmospheric CO₂ concentration, temperature, and vapor pressure deficit (VPD) on seawater-induced tree mortality is uncertain. We examined the mechanisms of seawater-induced mortality under varying climate scenarios using a photosynthetic gain and hydraulic cost optimization model validated against observations in a mature stand of Sitka spruce (Picea sitchensis) trees in the Pacific Northwest, USA, that were dying from recent seawater exposure. The simulations matched well with observations of photosynthesis, transpiration, nonstructural carbohydrates concentrations, leaf water potential, the percentage loss of xylem conductivity, and stand-level mortality rates. The simulations suggest that seawater-induced mortality could decrease by c. 16.7% with increasing atmospheric CO₂ levels due to reduced risk of carbon starvation. Conversely, rising VPD could increase mortality by c. 5.6% because of increasing risk of hydraulic failure. Across all scenarios, seawater-induced mortality was driven by hydraulic failure in the first 2 yr after seawater exposure began, with carbon starvation becoming more important in subsequent years. Changing CO₂ and climate appear unlikely to have a significant impact on coastal tree mortality under rising sea levels.

Stability of tropical forest tree carbon-water relations in a rainfall exclusion treatment through shifts in effective water uptake depth

Increasing severity and frequency of drought is predicted for large portions of the terrestrial biosphere, with major impacts already documented in wet tropical forests. Using a 4-year rainfall exclusion experiment in the Daintree Rainforest in northeast Australia, we examined canopy tree responses to reduced precipitation and soil water availability by quantifying seasonal changes in plant hydraulic and carbon traits for 11 tree species between control and drought treatments. Even with reduced soil volumetric water content in the upper 1 m of soil in the drought treatment, we found no significant difference between treatments for predawn and midday leaf water potential, photosynthesis, stomatal conductance, foliar stable carbon isotope composition, leaf mass per area, turgor loss point, xylem vessel anatomy, or leaf and stem nonstructural carbohydrates. While empirical measurements of aboveground traits revealed homeostatic maintenance of plant water status and traits in response to reduced soil moisture, modeled belowground dynamics revealed that trees in the drought treatment shifted the depth from which water was acquired to deeper soil layers. These findings reveal that belowground acclimation of tree water uptake depth may buffer tropical rainforests from more severe droughts that may arise in future with climate change.

Use of hydraulic traits for modeling genotype-specific acclimation in cotton under drought

Understanding the genetic and physiological basis of abiotic stress tolerance under field conditions is key to varietal crop improvement in the face of climate variability. Here, we investigate dynamic physiological responses to water stress in silico and their relationships to genotypic variation in hydraulic traits of cotton (Gossypium hirsutum), an economically important species for renewable textile fiber production. In conjunction with an ecophysiological process-based model, heterogeneous data (plant hydraulic traits, spatially-distributed soil texture, soil water content and canopy temperature) were used to examine hydraulic characteristics of cotton, evaluate their consequences on whole plant performance under drought, and explore potential genotype × environment effects. Cotton was found to have R-shaped hydraulic vulnerability curves (VCs), which were consistent under drought stress initiated at flowering. Stem VCs, expressed as percent loss of conductivity, differed across genotypes, whereas root VCs did not. Simulation results demonstrated how plant physiological stress can depend on the interaction between soil properties and irrigation management, which in turn affect genotypic rankings of transpiration in a time-dependent manner. Our study shows how a process-based modeling framework can be used to link genotypic variation in hydraulic traits to differential acclimating behaviors under drought.

Forecasting semi-arid biome shifts in the Anthropocene

Shrub encroachment, forest decline and wildfires have caused large-scale changes in semi-arid vegetation over the past 50 years. Climate is a primary determinant of plant growth in semi-arid ecosystems, yet it remains difficult to forecast large-scale vegetation shifts (i.e. biome shifts) in response to climate change. We highlight recent advances from four conceptual perspectives that are improving forecasts of semi-arid biome shifts. Moving from small to large scales, first, tree-level models that simulate the carbon costs of drought-induced plant hydraulic failure are improving predictions of delayed-mortality responses to drought. Second, tracer-informed water flow models are improving predictions of species coexistence as a function of climate. Third, new applications of ecohydrological models are beginning to simulate small-scale water movement processes at large scales. Fourth, remotely-sensed measurements of plant traits such as relative canopy moisture are providing early-warning signals that predict forest mortality more than a year in advance. We suggest that a community of researchers using modeling approaches (e.g. machine learning) that can integrate these perspectives will rapidly improve forecasts of semi-arid biome shifts. Better forecasts can be expected to help prevent catastrophic changes in vegetation states by identifying improved monitoring approaches and by prioritizing high-risk areas for management.

Conifers depend on established roots during drought: results from a coupled model of carbon allocation and hydraulics

Trees may survive prolonged droughts by shifting water uptake to reliable water sources, but it is unknown if the dominant mechanism involves activating existing roots or growing new roots during drought, or some combination of the two. To gain mechanistic insights on this unknown, a dynamic root-hydraulic modeling framework was developed that set up a feedback between hydraulic controls over carbon allocation and the role of root growth on soil-plant hydraulics. The new model was tested using a 5 yr drought/heat field experiment on an established piñon-juniper stand with root access to bedrock groundwater. Owing to the high carbon cost per unit root area, modeled trees initialized without adequate bedrock groundwater access experienced potentially lethal declines in water potential, while all of the experimental trees maintained nonlethal water potentials. Simulated trees were unable to grow roots rapidly enough to mediate the hydraulic stress, particularly during warm droughts. Alternatively, modeled trees initiated with root access to bedrock groundwater matched the hydraulics of the experimental trees by increasing their water uptake from bedrock groundwater when soil layers dried out. Therefore, the modeling framework identified a critical mechanism for drought response that required trees to shift water uptake among existing roots rather than growing new roots.

A framework for genomics-informed ecophysiological modeling in plants

Dynamic process-based plant models capture complex physiological response across time, carrying the potential to extend simulations out to novel environments and lend mechanistic insight to observed phenotypes. Despite the translational opportunities for varietal crop improvement that could be unlocked by linking natural genetic variation to first principles-based modeling, these models are challenging to apply to large populations of related individuals. Here we use a combination of model development, experimental evaluation, and genomic prediction in Brassica rapa L. to set the stage for future large-scale process-based modeling of intraspecific variation. We develop a new canopy growth submodel for B. rapa within the process-based model Terrestrial Regional Ecosystem Exchange Simulator (TREES), test input parameters for feasibility of direct estimation with observed phenotypes across cultivated morphotypes and indirect estimation using genomic prediction on a recombinant inbred line population, and explore model performance on an in silico population under non-stressed and mild water-stressed conditions. We find evidence that the updated whole-plant model has the capacity to distill genotype by environment interaction (G×E) into tractable components. The framework presented offers a means to link genetic variation with environment-modulated plant response and serves as a stepping stone towards large-scale prediction of unphenotyped, genetically related individuals under untested environmental scenarios.

Co-occurring woody species have diverse hydraulic strategies and mortality rates during an extreme drought

From 2011 to 2013, Texas experienced its worst drought in recorded history. This event provided a unique natural experiment to assess species-specific responses to extreme drought and mortality of four co-occurring woody species: Quercus fusiformis, Diospyros texana, Prosopis glandulosa, and Juniperus ashei. We examined hypothesized mechanisms that could promote these species' diverse mortality patterns using postdrought measurements on surviving trees coupled to retrospective process modelling. The species exhibited a wide range of gas exchange responses, hydraulic strategies, and mortality rates. Multiple proposed indices of mortality mechanisms were inconsistent with the observed mortality patterns across species, including measures of the degree of iso/anisohydry, photosynthesis, carbohydrate depletion, and hydraulic safety margins. Large losses of spring and summer whole-tree conductance (driven by belowground losses of conductance) and shallower rooting depths were associated with species that exhibited greater mortality. Based on this retrospective analysis, we suggest that species more vulnerable to drought were more likely to have succumbed to hydraulic failure belowground.

Phenotypic Trait Identification Using a Multimodel Bayesian Method: A Case Study Using Photosynthesis in Brassica rapa Genotypes

Agronomists have used statistical crop models to predict yield on a genotype-by-genotype basis. Mechanistic models, based on fundamental physiological processes common across plant taxa, will ultimately enable yield prediction applicable to diverse genotypes and crops. Here, genotypic information is combined with multiple mechanistically based models to characterize photosynthetic trait differentiation among genotypes of Brassica rapa. Infrared leaf gas exchange and chlorophyll fluorescence observations are analyzed using Bayesian methods. Three advantages of Bayesian approaches are employed: a hierarchical model structure, the testing of parameter estimates with posterior predictive checks and a multimodel complexity analysis. In all, eight models of photosynthesis are compared for fit to data and penalized for complexity using deviance information criteria (DIC) at the genotype scale. The multimodel evaluation improves the credibility of trait estimates using posterior distributions. Traits with important implications for yield in crops, including maximum rate of carboxylation (V_cmax ) and maximum rate of electron transport (J_max ) show genotypic differentiation. B. rapa shows phenotypic diversity in causal traits with the potential for genetic enhancement of photosynthesis. This multimodel screening represents a statistically rigorous method for characterizing genotypic differences in traits with clear biophysical consequences to growth and productivity within large crop breeding populations with application across plant processes.

Predicting stomatal responses to the environment from the optimization of photosynthetic gain and hydraulic cost

Stomatal regulation presumably evolved to optimize CO₂ for H₂ O exchange in response to changing conditions. If the optimization criterion can be readily measured or calculated, then stomatal responses can be efficiently modelled without recourse to empirical models or underlying mechanism. Previous efforts have been challenged by the lack of a transparent index for the cost of losing water. Yet it is accepted that stomata control water loss to avoid excessive loss of hydraulic conductance from cavitation and soil drying. Proximity to hydraulic failure and desiccation can represent the cost of water loss. If at any given instant, the stomatal aperture adjusts to maximize the instantaneous difference between photosynthetic gain and hydraulic cost, then a model can predict the trajectory of stomatal responses to changes in environment across time. Results of this optimization model are consistent with the widely used Ball-Berry-Leuning empirical model (r²  > 0.99) across a wide range of vapour pressure deficits and ambient CO₂ concentrations for wet soil. The advantage of the optimization approach is the absence of empirical coefficients, applicability to dry as well as wet soil and prediction of plant hydraulic status along with gas exchange.

Plant hydraulics improves and topography mediates prediction of aspen mortality in southwestern USA

Elevated forest mortality has been attributed to climate change-induced droughts, but prediction of spatial mortality patterns remains challenging. We evaluated whether introducing plant hydraulics and topographic convergence-induced soil moisture variation to land surface models (LSM) can help explain spatial patterns of mortality. A scheme predicting plant hydraulic safety loss from soil moisture was developed using field measurements and a plant physiology-hydraulics model, TREES. The scheme was upscaled to Populus tremuloides forests across Colorado, USA, using LSM-modeled and topography-mediated soil moisture, respectively. The spatial patterns of hydraulic safety loss were compared against aerial surveyed mortality. Incorporating hydraulic safety loss raised the explanatory power of mortality by 40% compared to LSM-modeled soil moisture. Topographic convergence was mostly influential in suppressing mortality in low and concave areas, explaining an additional 10% of the variations in mortality for those regions. Plant hydraulics integrated water stress along the soil-plant continuum and was more closely tied to plant physiological response to drought. In addition to the well-recognized topo-climate influence due to elevation and aspect, we found evidence that topographic convergence mediates tree mortality in certain parts of the landscape that are low and convergent, likely through influences on plant-available water.

Pragmatic hydraulic theory predicts stomatal responses to climatic water deficits

Ecosystem models have difficulty predicting plant drought responses, partially from uncertainty in the stomatal response to water deficits in soil and atmosphere. We evaluate a 'supply-demand' theory for water-limited stomatal behavior that avoids the typical scaffold of empirical response functions. The premise is that canopy water demand is regulated in proportion to threat to supply posed by xylem cavitation and soil drying. The theory was implemented in a trait-based soil-plant-atmosphere model. The model predicted canopy transpiration (E), canopy diffusive conductance (G), and canopy xylem pressure (P_canopy ) from soil water potential (P_soil ) and vapor pressure deficit (D). Modeled responses to D and P_soil were consistent with empirical response functions, but controlling parameters were hydraulic traits rather than coefficients. Maximum hydraulic and diffusive conductances and vulnerability to loss in hydraulic conductance dictated stomatal sensitivity and hence the iso- to anisohydric spectrum of regulation. The model matched wide fluctuations in G and P_canopy across nine data sets from seasonally dry tropical forest and piñon-juniper woodland with < 26% mean error. Promising initial performance suggests the theory could be useful in improving ecosystem models. Better understanding of the variation in hydraulic properties along the root-stem-leaf continuum will simplify parameterization.

Evaluating theories of drought-induced vegetation mortality using a multimodel-experiment framework

Model-data comparisons of plant physiological processes provide an understanding of mechanisms underlying vegetation responses to climate. We simulated the physiology of a piñon pine-juniper woodland (Pinus edulis-Juniperus monosperma) that experienced mortality during a 5 yr precipitation-reduction experiment, allowing a framework with which to examine our knowledge of drought-induced tree mortality. We used six models designed for scales ranging from individual plants to a global level, all containing state-of-the-art representations of the internal hydraulic and carbohydrate dynamics of woody plants. Despite the large range of model structures, tuning, and parameterization employed, all simulations predicted hydraulic failure and carbon starvation processes co-occurring in dying trees of both species, with the time spent with severe hydraulic failure and carbon starvation, rather than absolute thresholds per se, being a better predictor of impending mortality. Model and empirical data suggest that limited carbon and water exchanges at stomatal, phloem, and below-ground interfaces were associated with mortality of both species. The model-data comparison suggests that the introduction of a mechanistic process into physiology-based models provides equal or improved predictive power over traditional process-model or empirical thresholds. Both biophysical and empirical modeling approaches are useful in understanding processes, particularly when the models fail, because they reveal mechanisms that are likely to underlie mortality. We suggest that for some ecosystems, integration of mechanistic pathogen models into current vegetation models, and evaluation against observations, could result in a breakthrough capability to simulate vegetation dynamics.

Identification of environmental covariates of West Nile virus vector mosquito population abundance

The rapid spread of West Nile virus (WNv) in North America is a major public health concern. Culex pipiens-restuans is the principle mosquito vector of WNv in the northeastern United States while Aedes vexans is an important bridge vector of the virus in this region. Vector mosquito abundance is directly dependent on physical environmental factors that provide mosquito habitats. The objective of this research is to determine landscape elements that explain the population abundance and distribution of WNv vector mosquitoes using stepwise linear regression. We developed a novel approach for examining a large set of landscape variables based on a land use and land cover classification by selecting variables in stages to minimize multicollinearity. We also investigated the distance at which landscape elements influence abundance of vector populations using buffer distances of 200, 400, and 1000 m. Results show landscape effects have a significant impact on Cx. pipiens-estuans population distribution while the effects of landscape features are less important for prediction of Ae. vexans population distributions. Cx. pipiens-restuans population abundance is positively correlated with human population density, housing unit density, and urban land use and land cover classes and negatively correlated with age of dwellings and amount of forested land.

Meteorologically conditioned time-series predictions of West Nile virus vector mosquitoes

An empirical model to forecast West Nile virus mosquito vector populations is developed using time series analysis techniques. Specifically, multivariate seasonal autoregressive integrated moving average (SARIMA) models were developed for Aedes vexans and the combined group of Culex pipiens and Culex restuans in Erie County, New York. Weekly mosquito collections data were obtained for the four mosquito seasons from 2002 to 2005 from the Erie County Department of Health, Vector and Pest Control Program. Climate variables were tested for significance with cross-correlation analysis. Minimum temperature (T(min)), maximum temperature (T(max)), average temperature (T(ave)), precipitation (P), relative humidity (R(H)), and evapotranspiration (E(T)) were acquired from the Northeast Regional Climate Center (NRCC) at Cornell University. Weekly averages or sums of climate variables were calculated from the daily data. Other climate indexes were calculated and were tested for significance with the mosquito population data, including cooling degree days base 60 degrees (C(DD_60)), cooling degree days base 63 (C(DD_63)), cooling degree days base 65 (C(DD_65)), a ponding index (I(P)), and an interactive C(DD_65)-precipitation variable (C(DD_65) x P(week_4)). Ae. vexans were adequately modeled with a (2,1,1)(1,1,0)(52) SARIMA model. The combined group of Culex pipiens-restuans were modeled with a (0,1,1)(1,1,0)(52) SARIMA model. The most significant meteorological variables for forecasting Aedes vexans abundance was the interactive C(DD_65) x P(week_4) variable at a lag of two weeks, E(T) x E(T) at a lag of five weeks, and C(DD_65) x C(DD_65) at a lag of seven weeks. The most significant predictive variables for the grouped Culex pipiens-restuans were C(DD_63) x C(DD_63) at a lag of zero weeks, C(DD_63) at a lag of eight weeks, and the cumulative maximum ponding index (I(Pcum)) at a lag of zero weeks.

Use of temporal patterns in vapor pressure deficit to explain spatial autocorrelation dynamics in tree transpiration

To quantify the relationship between temporal and spatial variation in tree transpiration, we measured sap flow in 129 trees with constant-heat sap flow sensors in a subalpine forest in southern Wyoming, USA. The forest stand was located along a soil water gradient from a stream side to near the top of a ridge. The stand was dominated by Pinus contorta Dougl. ex Loud. with Picea engelmannii Parry ex Engelm and Abies lasiocarpa (Hook.) Nutt. present near the stream and scattered individuals of Populus tremuloides Michx. throughout the stand. We used a cyclic sampling design that maximized spatial information with a minimum number of samples for semivariogram analyses. All species exhibited previously established responses to environmental variables in which the dominant driver was a saturating response to vapor pressure deficit (D). This response to D is predictable from tree hydraulic theory in which stomatal conductance declines as D increases to prevent excessive cavitation. The degree to which stomatal conductance declines with D is dependent on both species and individual tree physiology and increases the variability in transpiration as D increases. We quantified this variability spatially by calculating the spatial autocorrelation within 0.2-kPa D bins. Across 11 bins of D, spatial autocorrelation in individual tree transpiration was inversely correlated to D and dropped from 45 to 20 m. Spatial autocorrelation was much less for transpiration per unit leaf area and not significant for transpiration per unit sapwood area suggesting that spatial autocorrelation within a particular D bin could be explained by tree size. Future research should focus on the mechanisms underlying tree size spatial variability, and the potentially broad applicability of the inverse relationship between D and spatial autocorrelation in tree transpiration.

Interannual consistency in canopy stomatal conductance control of leaf water potential across seven tree species

We investigated interannual variability of canopy transpiration per unit ground area (E (C)) and per unit leaf area (E (L)) across seven tree species in northern Wisconsin over two years. These species have previously been shown to be sufficient to upscale stand-level transpiration to the landscape level during one growing season. Our objective was to test whether a simple plant hydraulic model could capture interannual variation in transpiration. Three species, wetland balsam fir (Abies balsamea (L.) Mill), basswood (Tilia Americana L.) and speckled alder (Alnus rugosa (DuRoi) Spreng), had no change in E (C) or E (L) between 2000 and 2001. Red pine (Pinus resinosa Ait) had a 57 and 19% increase in E (C) and E (L), respectively, and sugar maple (Acer saccharum Marsh) had an 83 and 41% increase in E (C) and E (L), respectively, from 2000 to 2001. Quaking aspen (Populus tremuloides Michx) had a 50 and 21% decrease in E (C) and E (L), respectively, from 2000 to 2001 in response to complete defoliation by forest tent caterpillar (Malascoma distria Hüber) and subsequent lower total leaf area index of the reflushed foliage. White cedar (Thuja occidentalis L.) had a 20% decrease in both E (C) and E (L) caused by lowered surface water in wetlands in 2001 because of lower precipitation and wetland flow management. Upland A. balsamea increased E (L) and E (C) by 55 and 53%, respectively, as a result of release from light competition of the defoliated, overstory P. tremuloides. We hypothesized that regardless of different drivers of interannual variability in E (C) and E (L), minimum leaf water potential would be regulated at the same value. Minimum midday water potentials were consistent over the two years within each of the seven species despite large changes in transpiration between years. This regulation was independently verified by the exponential saturation between daily E (C) and vapor pressure deficit (D) and the tradeoff between a reference canopy stomatal conductance (G (S)) and the sensitivity of G (S) to D, indicating that trees with high G (S) must decrease G (S) in response to atmospheric drought faster than trees with low G (S). Our results show that models of forest canopy transpiration can be simplified by incorporating G (S) regulation of minimum leaf water potential for isohydric species.

A locus for autosomal dominant posterior polar cataract on chromosome 1p

Autosomal dominant congenital cataract is a clinically and genetically heterogeneous lens disease. Here we report the linkage of a locus for autosomal dominant posterior polar cataract (CPP) to the distal short arm of chromosome 1. To map the CPP locus we performed molecular genetic linkage analysis using microsatellite markers in a three-generation pedigree. After exclusion of 13 known loci and candidate lens genes for autosomal dominant cataract, we obtained significantly positive LOD scores for markers D1S508 (Z = 3.14, theta = 0) and D1S468 (Z = 2.71, theta = 0). Multipoint analysis gave a maximum LOD score of 3.48 (theta = 0.07) between markers D1S508 and D1S468. From haplotype data, however, CPP probably lies in the telomeric interval D1S2845-1pter, which includes the locus for the clinically distinct Volkman congenital cataract (CCV). This study provides the first evidence for genetic heterogeneity of autosomal dominant posterior polar cataract for which a locus had been linked previously to chromosome 16q.