tealeaves: an R package for modelling leaf temperature using energy budgets

Beyond a single temperature threshold: Applying a cumulative thermal stress framework to plant heat tolerance

Most plant thermal tolerance studies focus on single critical thresholds, which limit the capacity to generalise across studies and predict heat stress under natural conditions. In animals and microbes, thermal tolerance landscapes describe the more realistic, cumulative effects of temperature. We tested this in plants by measuring the decline in leaf photosynthetic efficiency (FV/FM) following a combination of temperatures and exposure times and then modelled these physiological indices alongside recorded environmental temperatures. We demonstrate that a general relationship between stressful temperatures and exposure durations can be effectively employed to quantify and compare heat tolerance within and across plant species and over time. Importantly, we show how FV/FM curves translate to plants under natural conditions, suggesting that environmental temperatures often impair photosynthetic function. Our findings provide more robust descriptors of heat tolerance in plants and suggest that heat tolerance in disparate groups of organisms can be studied with a single predictive framework.

Leaf thermal safety margins decline at hotter temperatures in a natural warming 'experiment' in the Amazon

The threat of rising global temperatures may be especially pronounced for low-latitude, lowland plant species that have evolved under stable climatic conditions. However, little is known about how these species may acclimate to elevated temperatures. Here, we leveraged a strong, steep thermal gradient along a natural geothermal river to assess the ability of woody plants in the Amazon to acclimate to elevated air temperatures. We measured leaf traits in six common tropical woody species along the thermal gradient to investigate whether individuals of these species: acclimate their thermoregulatory traits to maintain stable leaf temperatures despite higher ambient temperatures; acclimate their photosynthetic thermal tolerances to withstand hotter leaf temperatures; and whether acclimation is sufficient to maintain stable leaf thermal safety margins (TSMs) across different growth temperatures. Individuals of three species acclimated their thermoregulatory traits, and three species increased their thermal tolerances with growth temperature. However, acclimation was generally insufficient to maintain constant TSMs. Notwithstanding, leaf health was generally consistent across growth temperatures. Acclimation in woody Amazonian plants is generally too weak to maintain TSMs at high growth temperatures, supporting previous findings that Amazonian plants will be increasingly vulnerable to thermal stress as temperatures rise.

Parameterizing mechanistic niche models in biophysical ecology: a review of empirical approaches

Mechanistic niche models are computational tools developed using biophysical principles to address grand challenges in ecology and evolution, such as the mechanisms that shape the fundamental niche and the adaptive significance of traits. Here, we review the empirical basis of mechanistic niche models in biophysical ecology, which are used to answer a broad array of questions in ecology, evolution and global change biology. We describe the experiments and observations that are frequently used to parameterize these models and how these empirical data are then incorporated into mechanistic niche models to predict performance, growth, survival and reproduction. We focus on the physiological, behavioral and morphological traits that are frequently measured and then integrated into these models. We also review the empirical approaches used to incorporate evolutionary processes, phenotypic plasticity and biotic interactions. We discuss the importance of validation experiments and observations in verifying underlying assumptions and complex processes. Despite the reliance of mechanistic niche models on biophysical theory, empirical data have and will continue to play an essential role in their development and implementation.

Seasonal trends in leaf-level photosynthetic capacity and water use efficiency in a North American Eastern deciduous forest and their impact on canopy-scale gas exchange

Vegetative transpiration (E) and photosynthetic carbon assimilation (A) are known to be seasonally dynamic, with changes in their ratio determining the marginal water use efficiency (WUE). Despite an understanding that stomata play a mechanistic role in regulating WUE, it is still unclear how stomatal and nonstomatal processes influence change in WUE over the course of the growing season. As a result, limited understanding of the primary physiological drivers of seasonal dynamics of canopy WUE remains one of the largest uncertainties in earth system model projections of carbon and water exchange in temperate deciduous forest ecosystems. We investigated seasonal patterns in leaf-level physiological, hydraulic, and anatomical properties, including the seasonal progress of the stomatal slope parameter (g1 ; inversely proportional to WUE) and the maximum carboxylation rate (Vcmax ). Vcmax and g1 were seasonally variable; however, their patterns were not temporally synchronized. g1 generally showed an increasing trend until late in the season, while Vcmax peaked during the midsummer months. Seasonal progression of Vcmax was primarily driven by changes in leaf structural, and anatomical characteristics, while seasonal changes in g1 were most strongly related to changes in Vcmax and leaf hydraulics. Using a seasonally variable Vcmax and g1 to parameterize a canopy-scale gas exchange model increased seasonally aggregated A and E by 3% and 16%, respectively.

The effect of the vertical gradients of photosynthetic parameters on the CO2 assimilation and transpiration of a Panamanian tropical forest

Terrestrial biosphere models (TBMs) include the representation of vertical gradients in leaf traits associated with modeling photosynthesis, respiration, and stomatal conductance. However, model assumptions associated with these gradients have not been tested in complex tropical forest canopies. We compared TBM representation of the vertical gradients of key leaf traits with measurements made in a tropical forest in Panama and then quantified the impact of the observed gradients on simulated canopy-scale CO₂ and water fluxes. Comparison between observed and TBM trait gradients showed divergence that impacted canopy-scale simulations of water vapor and CO₂ exchange. Notably, the ratio between the dark respiration rate and the maximum carboxylation rate was lower near the ground than at the top-of-canopy, leaf-level water-use efficiency was markedly higher at the top-of-canopy, and the decrease in maximum carboxylation rate from the top-of-canopy to the ground was less than TBM assumptions. The representation of the gradients of leaf traits in TBMs is typically derived from measurements made within-individual plants, or, for some traits, assumed constant due to a lack of experimental data. Our work shows that these assumptions are not representative of the trait gradients observed in species-rich, complex tropical forests.

Mechanistic forecasts of species responses to climate change: The promise of biophysical ecology

A core challenge in global change biology is to predict how species will respond to future environmental change and to manage these responses. To make such predictions and management actions robust to novel futures, we need to accurately characterize how organisms experience their environments and the biological mechanisms by which they respond. All organisms are thermodynamically connected to their environments through the exchange of heat and water at fine spatial and temporal scales and this exchange can be captured with biophysical models. Although mechanistic models based on biophysical ecology have a long history of development and application, their use in global change biology remains limited despite their enormous promise and increasingly accessible software. We contend that greater understanding and training in the theory and methods of biophysical ecology is vital to expand their application. Our review shows how biophysical models can be implemented to understand and predict climate change impacts on species' behavior, phenology, survival, distribution, and abundance. It also illustrates the types of outputs that can be generated, and the data inputs required for different implementations. Examples range from simple calculations of body temperature at a particular site and time, to more complex analyses of species' distribution limits based on projected energy and water balances, accounting for behavior and phenology. We outline challenges that currently limit the widespread application of biophysical models relating to data availability, training, and the lack of common software ecosystems. We also discuss progress and future developments that could allow these models to be applied to many species across large spatial extents and timeframes. Finally, we highlight how biophysical models are uniquely suited to solve global change biology problems that involve predicting and interpreting responses to environmental variability and extremes, multiple or shifting constraints, and novel abiotic or biotic environments.

Thermal sensitivity across forest vertical profiles: patterns, mechanisms, and ecological implications

Rising temperatures are influencing forests on many scales, with potentially strong variation vertically across forest strata. Using published research and new analyses, we evaluate how microclimate and leaf temperatures, traits, and gas exchange vary vertically in forests, shaping tree, and ecosystem ecology. In closed-canopy forests, upper canopy leaves are exposed to the highest solar radiation and evaporative demand, which can elevate leaf temperature (T_leaf ), particularly when transpirational cooling is curtailed by limited stomatal conductance. However, foliar traits also vary across height or light gradients, partially mitigating and protecting against the elevation of upper canopy T_leaf . Leaf metabolism generally increases with height across the vertical gradient, yet differences in thermal sensitivity across the gradient appear modest. Scaling from leaves to trees, canopy trees have higher absolute metabolic capacity and growth, yet are more vulnerable to drought and damaging T_leaf than their smaller counterparts, particularly under climate change. By contrast, understory trees experience fewer extreme high T_leaf 's but have fewer cooling mechanisms and thus may be strongly impacted by warming under some conditions, particularly when exposed to a harsher microenvironment through canopy disturbance. As the climate changes, integrating the patterns and mechanisms reviewed here into models will be critical to forecasting forest-climate feedback.

Limited acclimation of leaf traits and leaf temperatures in a subtropical urban heat island

The consequences of rising temperatures for trees will vary between species based on their abilities to acclimate their leaf thermoregulatory traits and photosynthetic thermal tolerances. We tested the hypotheses that adult trees in warmer growing conditions (i) acclimate their thermoregulatory traits to regulate leaf temperatures, (ii) acclimate their thermal tolerances such that tolerances are positively correlated with leaf temperature and (iii) that species with broader thermal niche breadths have greater acclimatory abilities. To test these hypotheses, we measured leaf traits and thermal tolerances of seven focal tree species across steep thermal gradients in Miami's urban heat island. We found that some functional traits varied significantly across air temperatures within species. For example, leaf thickness increased with maximum air temperature in three species, and leaf mass per area and leaf reflectance both increased with air temperature in one species. Only one species was marginally more homeothermic than expected by chance due to acclimation of its thermoregulatory traits, but this acclimation was insufficient to offset elevated air temperatures. Thermal tolerances acclimated to higher maximum air temperatures in two species. As a result of limited acclimation, leaf thermal safety margins (TSMs) were narrower for trees in hotter areas. We found some support for our hypothesis that species with broader thermal niches are better at acclimating to maintain more stable TSMs across the temperature gradients. These findings suggest that trees have limited abilities to acclimate to high temperatures and that thermal niche specialists may be at a heightened risk of thermal stress as global temperatures continue to rise.

Accurate phenology analyses require bud traits and energy budgets

Spring phenology is mainly driven by temperature in extratropical ecosystems. Recent evidence highlighted the key role of micrometeorology and bud temperature on delaying or advancing leaf unfolding. Yet, phenology studies, either using ground-based or remote sensing observations, always substitute plant tissue temperature by air temperature. In fact, temperatures differ substantially between plant tissues and the air because plants absorb and lose energy. Here, we build on recent observations and well-established energy balance theories to discuss how solar radiation, wind and bud traits might affect our interpretation of spring phenology sensitivity to warming. We show that air temperature might be an imprecise and biased predictor of bud temperature. Better characterizing the plants' phenological response to warming will require new observations of bud traits and temperature for accurately quantifying their energy budget. As consistent micrometeorology datasets are still scarce, new approaches coupling energy budget modelling and plant traits could help to improve phenology analyses across scales.

Late-day measurement of excised branches results in uncertainty in the estimation of two stomatal parameters derived from response curves in Populus deltoides Bartr. × Populus nigra L

Many terrestrial biosphere models depend on an understanding of the relationship between stomatal conductance and photosynthesis. However, unlike the measurement of photosynthetic parameters, such as the maximum carboxylation capacity, where standard methods (e.g., CO2 response or ACi curves) are widely accepted, a consensus method for empirically measuring parameters representing stomatal response has not yet emerged. Most models of stomatal response to environment represent stomatal conductance as being bounded by a lower intercept parameter (g0), and linearly scaled based on a multivariate term described by the stomatal slope parameter (g1). Here we employ the widely used Unified Stomatal Optimization model, to test whether g1 and g0 parameters are impacted by the choice of measurement method, either on an intact branch or a cut branch segment stored in water. We measured paired stomatal response curves on intact and excised branches of a hybrid poplar clone (Populus deltoides Bartr. × Populus nigra L. OP367), measured twice over a diurnal period. We found that predawn branch excision did not significantly affect measured g0 and g1 when measured within 4 h of excision. Measurement in the afternoon resulted in significantly higher values of g1 and lower values of g0, with values changing by 55% and 56%, respectively. Excision combined with afternoon measurement resulted in a marked effect on parameter estimates, with g1 increasing 89% from morning to afternoon and a 25% lower g1 for cut branches than those measured in situ. We also show that in hybrid poplar the differences in parameter estimates obtained from plants measured under different conditions can directly impact models of canopy function, reducing modeled transpiration by 18% over a simulated 12.5-h period. Although these results are only for a single isohydric woody species, our findings suggest that stomatal optimality parameters may not remain constant throughout the day.

Tradeoffs between leaf cooling and hydraulic safety in a dominant arid land riparian tree species

Leaf carbon gain optimization in hot environments requires balancing leaf thermoregulation with avoiding excessive water loss via transpiration and hydraulic failure. The tradeoffs between leaf thermoregulation and transpirational water loss can determine the ecological consequences of heat waves that are increasing in frequency and intensity. We evaluated leaf thermoregulation strategies in warm- (>40°C maximum summer temperature) and cool-adapted (<40°C maximum summer temperature) genotypes of the foundation tree species, Populus fremontii, using a common garden near the mid-elevational point of its distribution. We measured leaf temperatures and assessed three modes of leaf thermoregulation: leaf morphology, midday canopy stomatal conductance and stomatal sensitivity to vapour pressure deficit. Data were used to parameterize a leaf energy balance model to estimate contrasts in midday leaf temperature in warm- and cool-adapted genotypes. Warm-adapted genotypes had 39% smaller leaves and 38% higher midday stomatal conductance, reflecting a 3.8°C cooler mean leaf temperature than cool-adapted genotypes. Leaf temperatures modelled over the warmest months were on average 1.1°C cooler in warm- relative to cool-adapted genotypes. Results show that plants adapted to warm environments are predisposed to tightly regulate leaf temperatures during heat waves, potentially at an increased risk of hydraulic failure.

The use of R in photosynthesis research

R is one of the most commonly used analytical tools in the plant sciences. To identify key trends in general reported R use and patterns in photosynthesis research, we explored the frequency of R use in 2966 articles published in the 377 journals with 'photosynthesis' in the title from 2010 to 2019 using the Web of Science search. Solutions provided by each R package cited in the articles or online sources was recorded and classified. The percentage of research articles reporting R use increased linearly from 3.6% in 2010 to 12.5% in 2019. The three main categories of R package solutions were 'general statistical calculations and graph packages' (G); 'photosynthesis special-purpose packages' (S); and 'genetic and evolutionary packages' (E). The top five R packages cited were nlme (G), lme4 (G), multcomp (G), plantecophys (S), and ape (E). The increasing popularity of R use in photosynthesis research is due to its user-friendly and abundant open-source codes online for handling specific issues, particularly in fitting photosynthesis models. These findings are limited by the number of articles and online sources, but they reveal a significant increase in usage in photosynthesis research over the past decade and have a bright prospect in the future.

Principles of resilient coding for plant ecophysiologists

Plant ecophysiology is founded on a rich body of physical and chemical theory, but it is challenging to connect theory with data in unambiguous, analytically rigorous and reproducible ways. Custom scripts written in computer programming languages (coding) enable plant ecophysiologists to model plant processes and fit models to data reproducibly using advanced statistical techniques. Since many ecophysiologists lack formal programming education, we have yet to adopt a unified set of coding principles and standards that could make coding easier to learn, use and modify. We identify eight principles to help in plant ecophysiologists without much programming experience to write resilient code: (i) standardized nomenclature, (ii) consistency in style, (iii) increased modularity/extensibility for easier editing and understanding, (iv) code scalability for application to large data sets, (v) documented contingencies for code maintenance, (vi) documentation to facilitate user understanding; (vii) extensive tutorials and (viii) unit testing and benchmarking. We illustrate these principles using a new R package, {photosynthesis}, which provides a set of analytical and simulation tools for plant ecophysiology. Our goal with these principles is to advance scientific discovery in plant ecophysiology by making it easier to use code for simulation and data analysis, reproduce results and rapidly incorporate new biological understanding and analytical tools.

Seasonal trends in photosynthesis and leaf traits in scarlet oak

Understanding seasonal variation in photosynthesis is important for understanding and modeling plant productivity. Here, we used shotgun sampling to examine physiological, structural and spectral leaf traits of upper canopy, sun-exposed leaves in Quercus coccinea Münchh (scarlet oak) across the growing season in order to understand seasonal trends, explore the mechanisms underpinning physiological change and investigate the impact of extrapolating measurements from a single date to the whole season. We tested the hypothesis that photosynthetic rates and capacities would peak at the summer solstice, i.e., at the time of peak photoperiod. Contrary to expectations, our results reveal a late-season peak in both photosynthetic capacity and rate before the expected sharp decrease at the start of senescence. This late-season maximum occurred after the higher summer temperatures and vapor pressure deficit and was correlated with the recovery of leaf water content and increased stomatal conductance. We modeled photosynthesis at the top of the canopy and found that the simulated results closely tracked the maximum carboxylation capacity of Rubisco. For both photosynthetic capacity and modeled top-of-canopy photosynthesis, the maximum value was therefore not observed at the summer solstice. Rather, in each case, the measurements at and around the solstice were close to the overall seasonal mean, with values later in the season leading to deviations from the mean by up to 41 and 52%, respectively. Overall, we found that the expected Gaussian pattern of photosynthesis was not observed. We conclude that an understanding of species- and environment-specific changes in photosynthesis across the season is essential for correct estimation of seasonal photosynthetic capacity.

Developmental and biophysical determinants of grass leaf size worldwide

One of the most notable ecological trends-described more than 2,300  years ago by Theophrastus-is the association of small leaves with dry and cold climates, which has recently been recognized for eudicotyledonous plants at a global scale^1-3. For eudicotyledons, this pattern has been attributed to the fact that small leaves have a thinner boundary layer that helps to avoid extreme leaf temperatures⁴ and their leaf development results in vein traits that improve water transport under cold or dry climates^5,6. However, the global distribution of leaf size and its adaptive basis have not been tested in the grasses, which represent a diverse lineage that is distinct in leaf morphology and that contributes 33% of terrestrial primary productivity (including the bulk of crop production)⁷. Here we demonstrate that grasses have shorter and narrower leaves under colder and drier climates worldwide. We show that small grass leaves have thermal advantages and vein development that contrast with those of eudicotyledons, but that also explain the abundance of small leaves in cold and dry climates. The worldwide distribution of leaf size in grasses exemplifies how biophysical and developmental processes result in convergence across major lineages in adaptation to climate globally, and highlights the importance of leaf size and venation architecture for grass performance in past, present and future ecosystems.

Leaf temperature impacts canopy water use efficiency independent of changes in leaf level water use efficiency

Canopy water use efficiency (above-ground biomass over lifetime water loss, WUE_canopy) can influence yield in wheat and other crops. Breeding for WUE_canopy is difficult because it is influenced by many component traits. For example, intrinsic water use efficiency (WUE_i), the ratio of net carbon assimilation (A_net) over stomatal conductance, contributes to WUE_canopy and can be estimated from carbon isotope discrimination (Δ). However, Δ is not sensitive to differences in the water vapor pressure deficit between the air and leaf (VPD_leaf). Alternatively, measurements of instantaneous leaf water use efficiency (WUE_leaf) are defined as A_net over transpiration and can be determined with gas exchange, but the dynamic nature of field conditions are not represented. Specifically, fluctuations in canopy temperature lead to changes in VPD_leaf that impact transpiration but not A_net. This alters WUE_leaf and in turn affects WUE_canopy. To test this relationship, WUE_canopy was measured in conjunction with WUE_i, WUE_canopy, and canopy temperature under well-watered and water-limited conditions in two drought-tolerant wheat cultivars that differ in canopy architecture. In this experiment, boundary layer conductance was low and significant changes in leaf temperature occurred between cultivars and treatments that correlated with WUE_canopy likely because of the effect of canopy temperature on VPD_leaf driving T. However, deviations between WUE_i, WUE_leaf, and WUE_canopy were present because measurements made at the leaf level do not account for variations in leaf temperature. This uncoupled the relationship of measured WUE_leaf and WUE_i from WUE_canopy and emphasizes the importance of canopy temperature on carbon uptake and transpired water loss.

Developmental changes in the reflectance spectra of temperate deciduous tree leaves and implications for thermal emissivity and leaf temperature

Leaf optical properties impact leaf energy balance and thus leaf temperature. The effect of leaf development on mid-infrared (MIR) reflectance, and hence thermal emissivity, has not been investigated in detail. We measured a suite of morphological characteristics, as well as directional-hemispherical reflectance from ultraviolet to thermal infrared wavelengths (250 nm to 20 µm) of leaves from five temperate deciduous tree species over the 8 wk following spring leaf emergence. By contrast to reflectance at shorter wavelengths, the shape and magnitude of MIR reflectance spectra changed markedly with development. MIR spectral differences among species became more pronounced and unique as leaves matured. Comparison of reflectance spectra of intact vs dried and ground leaves points to cuticular development - and not internal structural or biochemical changes - as the main driving factor. Accompanying the observed spectral changes was a drop in thermal emissivity from about 0.99 to 0.95 over the 8 wk following leaf emergence. Emissivity changes were not large enough to substantially influence leaf temperature, but they could potentially lead to a bias in radiometrically measured temperatures of up to 3 K. Our results also pointed to the potential for using MIR spectroscopy to better understand species-level differences in cuticular development and composition.

Empirical evidence for resilience of tropical forest photosynthesis in a warmer world

Tropical forests may be vulnerable to climate change^1-3 if photosynthetic carbon uptake currently operates near a high temperature limit^4-6. Predicting tropical forest function requires understanding the relative contributions of two mechanisms of high-temperature photosynthetic declines: stomatal limitation (H1), an indirect response due to temperature-associated changes in atmospheric vapour pressure deficit (VPD)⁷, and biochemical restrictions (H2), a direct temperature response^8,9. Their relative control predicts different outcomes-H1 is expected to diminish with stomatal responses to future co-occurring elevated atmospheric [CO₂], whereas H2 portends declining photosynthesis with increasing temperatures. Distinguishing the two mechanisms at high temperatures is therefore critical, but difficult because VPD is highly correlated with temperature in natural settings. We used a forest mesocosm to quantify the sensitivity of tropical gross ecosystem productivity (GEP) to future temperature regimes while constraining VPD by controlling humidity. We then analytically decoupled temperature and VPD effects under current climate with flux-tower-derived GEP trends in situ from four tropical forest sites. Both approaches showed consistent, negative sensitivity of GEP to VPD but little direct response to temperature. Importantly, in the mesocosm at low VPD, GEP persisted up to 38 °C, a temperature exceeding projections for tropical forests in 2100 (ref. ¹⁰). If elevated [CO₂] mitigates VPD-induced stomatal limitation through enhanced water-use efficiency as hypothesized^9,11, tropical forest photosynthesis may have a margin of resilience to future warming.