Search

Stomata play a central role in surface-atmosphere exchange by controlling the flux of water and CO

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.

Tropical leaf phenology—particularly its variability at the tree-crown scale—dominates the seasonality of carbon and water fluxes. However, given enormous species diversity, accurate means of monitoring leaf phenology in tropical forests is still lacking. Time series of the Green Chromatic Coordinate (GCC) metric derived from tower-based red–greenblue (RGB) phenocams have been widely used to monitor leaf phenology in temperate forests, but its application in the tropics remains problematic. To improve monitoring of tropical phenology, we explored the use of a deep learning model (i.e. superpixel-based Residual Networks 50, SP-ResNet50) to automatically differentiate leaves from non-leaves in phenocam images and to derive leaf fraction at the tree-crown scale. To evaluate our model, we used a year of data from six phenocams in two contrasting forests in Panama. We first built a comprehensive library of leaf and non-leaf pixels across various acquisition times, exposure conditions and specific phenocams. We then divided this library into training and testing components. We evaluated the model at three levels: 1) superpixel level with a testing set, 2) crown level by comparing the model-derived leaf fractions with those derived using image-specific supervised classification, and 3) temporally using all daily images to assess the diurnal stability of the model-derived leaf fraction. Finally, we compared the model-derived leaf fraction phenology with leaf phenology derived from GCC. Our results show that: 1) the SP-ResNet50 model accurately differentiates leaves from non-leaves (overall accuracy of 93%) and is robust across all three levels of evaluations; 2) the model accurately quantifies leaf fraction phenology across tree-crowns and forest ecosystems; and 3) the combined use of leaf fraction and GCC helps infer the timing of leaf emergence, maturation and senescence, critical information for modeling photosynthetic seasonality of tropical forests. Collectively, this study offers an improved means for automated tropical phenology monitoring using phenocams.

To investigate the burden of infectious complications following ureteroscopy (URS) for ureteric stones on a national level in England using data from the Hospital Episodes Statistics (HES) data warehouse.A retrospective cohort was identified and followed up in HES during the period April 2013 to March 2020 for all procedure codes relating to ureteroscopic stone treatment (M27.1, M27.2, M27.3). Treatment episodes relating to the first URS ('index ureteroscopy') for each patient were further analysed. All subsequent admissions within 30 days were also captured. The primary outcome was diagnosis of urinary tract infection (UTI; including all codes relating to a UTI/sepsis within the first 30 days of index URS). Secondary outcomes were critical care attendance, attendance at the accident and emergency department (AE) within 30 days, and mortality.A total of 71 305 index ureteroscopies were eligible for analysis. The median age was 55 years, and 81% of procedures were elective and 45% were undertaken as day-cases. At the time of index URS, 16% of patients had diabetes, 0.5% had coexisting neurological disease and 40% had an existing stent/nephrostomy. Overall, 6.8% of the cohort (n = 4822) had a diagnosis of UTI within 30 days of index URS (3.9% immediately after surgery). A total of 339 patients (0.5%) required an unplanned stay in critical care during their index URS admission; 8833 patients (12%) attended AE within 30 days. Overall mortality was 0.18% (60 in-hospital, 65 within 30 days); 40 deaths (0.056%) included infection as a contributing cause of death.We present the largest series evaluating infectious complications after ureteroscopic stone treatment. The procedure is safe, with low inpatient infective complication and critical care admission rates.

Partial least squares regression (PLSR) modelling is a statistical technique for correlating datasets, and involves the fitting of a linear regression between two matrices. One application of PLSR enables leaf traits to be estimated from hyperspectral optical reflectance data, facilitating rapid, high-throughput, non-destructive plant phenotyping. This technique is of interest and importance in a wide range of contexts including crop breeding and ecosystem monitoring. The lack of a consensus in the literature on how to perform PLSR means that interpreting model results can be challenging, applying existing models to novel datasets can be impossible, and unknown or undisclosed assumptions can lead to incorrect or spurious predictions. We address this lack of consensus by proposing best practices for using PLSR to predict plant traits from leaf-level hyperspectral data, including a discussion of when PLSR is applicable, and recommendations for data collection. We provide a tutorial to demonstrate how to develop a PLSR model, in the form of an R script accompanying this manuscript. This practical guide will assist all those interpreting and using PLSR models to predict leaf traits from spectral data, and advocates for a unified approach to using PLSR for predicting traits from spectra in the plant sciences.

Objectives: Urinary biochemistry is used to detect and monitor conditions associated with recurrent kidney stones. There are no predictive machine learning (ML) tools for kidney stone type or recurrence. We therefore aimed to build and validate ML models for these outcomes using age, gender, 24-hour urine biochemistry and stone composition. Materials and Methods: Data from 3 cohorts were used, Southampton, UK (n=3013), Newcastle, UK (n=5984) and Bern, Switzerland (n=794). Of these 3130 had available 24-hour urine biochemistry measurements (calcium, oxalate, urate, pH, volume), and 1684 had clinical data on kidney stone recurrence. Two predictive models were constructed (UK and Swiss) using two ML techniques (Partitioning and Random Forests [RF]) and validated internally with a subset of the same dataset (e.g UK model/UK test set), and externally with the other dataset (UK model/Swiss test set). Results and Limitations: For kidney stone type, on external validation accuracy of UK RF model=0.79 (95% CI: 0.73-0.84), sensitivity: calcium oxalate=0.99 and calcium phosphate/urate=0.00. Specificity: calcium oxalate=0.00 and calcium phosphate/urate=0.99. For the Swiss RF model accuracy=0.87 (95% CI: 0.83-0.89), sensitivity: calcium oxalate=0.99 and calcium phosphate/urate=0.00. Specificity: calcium oxalate=0.00, calcium phosphate=0.00 and urate=1.00. For stone recurrence, on external validation accuracy of UK RF model=0.22 (95% CI: 0.19-0.25), sensitivity=0.93 and specificity=0.09. Swiss RF model accuracy=0.42 (95% CI: 0.39-0.47), sensitivity=0.03 and specificity=0.97. Conclusions: Neither kidney stone type nor kidney stone recurrence can be accurately predicted using modelling tools built using specific 24-hour urinary biochemistry values alone. Further studies to delineate accurate predictive tools should be undertaken using both known and novel risk factors.

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.

The representation of stomatal regulation of transpiration and CO

Stomata play a central role in regulating the exchange of carbon dioxide and water vapor between ecosystems and the atmosphere. Their function is represented in land surface models (LSMs) by conductance models. The Functionally Assembled Terrestrial Ecosystem Simulator (FATES) is a dynamic vegetation demography model that can simulate both detailed plant demographic and physiological dynamics. To evaluate the effect of stomatal conductance model formulation on forest water and carbon fluxes in FATES, we implemented an optimality-based stomatal conductance model – the Medlyn (MED) model – that simulates the relationship between photosynthesis (A) and stomatal conductance to water vapor (gsw) as an alternative to the FATES default Ball–Woodrow–Berry (BWB) model. To evaluate how the behavior of FATES is affected by stomatal model choice, we conducted a model sensitivity analysis to explore the response of gsw to climate forcing, including atmospheric CO2 concentration, air temperature, radiation, and vapor pressure deficit in the air (VPDa). We found that modeled gsw values varied greatly between the BWB and MED formulations due to the different default stomatal slope parameters (g1). After harmonizing g1 and holding the stomatal intercept parameter (g0) constant for both model formulations, we found that the divergence in modeled gsw was limited to conditions when the VPDa exceeded 1.5 kPa. We then evaluated model simulation results against measurements from a wet evergreen forest in Panama. Results showed that both the MED and BWB model formulations were able to capture the magnitude and diurnal changes of measured gsw and A but underestimated both by about 30 % when the soil was predicted to be very dry. Comparison of modeled soil water content from FATES to a reanalysis product showed that FATES captured soil drying well, but translation of drying soil to modeled physiology reduced the models' ability to match observations. Our study suggests that the parameterization of stomatal conductance models and current model response to drought are the critical areas for improving model simulation of CO2 and water fluxes in tropical forests.

Tropical forests absorb large amounts of atmospheric CO2 through photosynthesis, but high surface temperatures suppress this absorption while promoting isoprene emissions. While mechanistic isoprene emission models predict a tight coupling to photosynthetic electron transport (ETR) as a function of temperature, direct field observations of this phenomenon are lacking in the tropics and are necessary to assess the impact of a warming climate on global isoprene emissions. Here we demonstrate that in the early successional species Vismia guianensis in the central Amazon, ETR rates increased with temperature in concert with isoprene emissions, even as stomatal conductance (gs ) and net photosynthetic carbon fixation (Pn ) declined. We observed the highest temperatures of continually increasing isoprene emissions yet reported (50°C). While Pn showed an optimum value of 32.6 ± 0.4°C, isoprene emissions, ETR, and the oxidation state of PSII reaction centers (qL ) increased with leaf temperature with strong linear correlations for ETR (ƿ = 0.98) and qL (ƿ = 0.99) with leaf isoprene emissions. In contrast, other photoprotective mechanisms, such as non-photochemical quenching, were not activated at elevated temperatures. Inhibition of isoprenoid biosynthesis repressed Pn at high temperatures through a mechanism that was independent of stomatal closure. While extreme warming will decrease gs and Pn in tropical species, our observations support a thermal tolerance mechanism where the maintenance of high photosynthetic capacity under extreme warming is assisted by the simultaneous stimulation of ETR and metabolic pathways that consume the direct products of ETR including photorespiration and the biosynthesis of thermoprotective isoprenoids. Our results confirm that models which link isoprene emissions to the rate of ETR hold true in tropical species and provide necessary "ground-truthing" for simulations of the large predicted increases in tropical isoprene emissions with climate warming.

Background: The aim of this study is to investigate outcomes of pre-operative stent dwell time on infectious complications following ureteroscopy and stone treatment to identify a time cut-off. Material and Methods: Three tertiary referral centres in Europe retrospectively collected outcomes of ureteroscopy and laser fragmentation (URSL) for all patients with pre-operative indwelling ureteric stents over a period of up to 5 years. Data was collected on patient details, stone demographics, stent dwell time, complications and stone free rate (SFR). Matching for age, sex, operative time, stone size and post-operative stent insertion. To examine for a threshold effect, monthly cut-offs were used to compare post-ureteroscopic febrile UTIs. Binomial logistic regression was used (SPSS v.24) with a significance level set at 0.0036. The risk ratio (RR) with a 95% confidence interval (CI) and the number needed to harm (NNH) are reported. Results: There were 467 patients with a pre-operative stent for analysis. These patients (n = 315) were matched to non-stented controls after excluding 152 patients to achieve adequate matching. There was a significant difference in rates of post-ureteroscopic febrile UTI between stented vs non-stented patients (RR = 2.67, 95% CI: 1.10–6.48, p = 0.03). On adjustment, a dwell time of more than two months was associated with an increased risk of post-ureteroscopic febrile UTI (RR = 3.94, 95% CI: 1.30–12.01, p = 0.02), this increased risk rose with longer dwell time. At stent time longer than four months was associated with a significantly increased risk of post-ureteroscopic febrile UTI (5% vs. 15%, RR = 3.09, 95% CI: 1.56–6.10, p = 0.001), with the number needed to harm at 10. Conclusions: Overall infectious complication rates from URSL are low. The risk of post-operative UTI after four months of dwell time is nearly tripled compared to less than four months.

Leaf respiration in the dark (Rdark) and light (Rday) is poorly characterized in diverse tropical ecosystems, and little to no information exists on the degree of light suppression in common tree species within the Amazon basin, and their dependences upon plant functional traits and position within the canopy. We quantified Rdark and apparent Rday using the Kok method and measured key leaf traits in 26 tree individuals of different species distributed in three different canopy positions: canopy, lower canopy, and understory. To explore the relationships between the leaf traits we used the standardized major axis (SMA). We found that canopy trees had significantly higher rates of Rdark and Rday than trees in the understory. The difference between Rdark and Rday (the light suppression of respiration) was greatest in the understory (68 ± 9%, 95% CI) and lower canopy (49 ± 9%, 95% CI) when compared to the canopy (37 ± 10%, 95% CI). We also found that Rday was significantly and strongly correlated with Rdark (p < 0.001) for all the canopy positions. Also, leaf mass per area (LMA) and leaf Phosphorus concentration (P) had a significant relationship with Rdark (p < 0.001; p = 0.003), respectively. In addition, a significant relationship was found for LMA in the canopy and lower canopy positions (p = 0.009; p = 0.048) while P was only significant in the canopy (p = 0.044). Finally, no significant relationship was found between Rdark and nitrogen, sugars, and starch. Our results highlight the importance of including representation of the light suppression of leaf respiration in terrestrial biosphere models and also of accounting for vertical gradients within forest canopies and connections with functional traits.

Author(s): Crystal-Ornelas, R; Varadharajan, C; Bond-Lamberty, B; Boye, K; Burrus, M; Cholia, S; Crow, M; Damerow, J; Devarakonda, R; Ely, KS; Goldman, A; Heinz, S; Hendrix, V; Kakalia, Z; Pennington, SC; Robles, E; Rogers, A; Simmonds, M; Velliquette, T; Weierbach, H; Weisenhorn, P; Welch, JN; Agarwal, DA | Abstract: Data standardization combined with descriptive metadata facilitate data reuse, which is the ultimate goal of the Findable, Accessible, Interoperable, and Reusable (FAIR) principles. Community data or metadata standards are increasingly created through an approach that emphasizes collaboration between various stakeholders. Such an approach requires platforms for collaboration on the development process that centers on sharing information and receiving feedback. Our objective in this study was to conduct a systematic review to identify data standards and reporting formats that use version control for developing data standards and to summarize common practices, particularly in earth and environmental sciences. Out of 108 data standards and reporting formats identified in our review, 32 used GitHub as the version control platform, and no other platforms were used. We found no universally accepted methodology for developing and publishing data standards. Many GitHub repositories did not use key features that could help developers to gather user feedback, or to create and revise standards that build on previous work. We provide guidance for community-driven standard development and associated documentation on GitHub based on a systematic review of existing practices.

The Arctic-Boreal Region (ABR) has a large impact on global vegetation-atmosphere interactions and is experiencing markedly greater warming than the rest of the planet, a trend that is projected to continue with anticipated future emissions of CO2 . The ABR is a significant source of uncertainty in estimates of carbon uptake in terrestrial biosphere models such that reducing this uncertainty is critical for more accurately estimating global carbon cycling and understanding the response of the region to global change. Process representation and parameterization associated with gross primary productivity (GPP) drives a large amount of this model uncertainty, particularly within the next 50 years, where the response of existing vegetation to climate change will dominate estimates of GPP for the region. Here we review our current understanding and model representation of GPP in northern latitudes, focusing on vegetation composition, phenology, and physiology, and consider how climate change alters these three components. We highlight challenges in the ABR for predicting GPP, but also focus on the unique opportunities for advancing knowledge and model representation, particularly through the combination of remote sensing and traditional boots-on-the-ground science.

Tropical forests are one of the main carbon sinks on Earth, but the magnitude of CO2 absorbed by tropical vegetation remains uncertain. Terrestrial biosphere models (TBMs) are commonly used to estimate the CO2 absorbed by forests, but their performance is highly sensitive to the parameterization of processes that control leaf-level CO2 exchange. Direct measurements of leaf respiratory and photosynthetic traits that determine vegetation CO2 fluxes are critical, but traditional approaches are time-consuming. Reflectance spectroscopy can be a viable alternative for the estimation of these traits and, because data collection is markedly quicker than traditional gas exchange, the approach can enable the rapid assembly of large datasets. However, the application of spectroscopy to estimate photosynthetic traits across a wide range of tropical species, leaf ages and light environments has not been extensively studied. Here, we used leaf reflectance spectroscopy together with partial least-squares regression (PLSR) modeling to estimate leaf respiration (Rdark25), the maximum rate of carboxylation by the enzyme Rubisco (Vcmax25), the maximum rate of electron transport (Jmax25), and the triose phosphate utilization rate (Tp25), all normalized to 25°C. We collected data from three tropical forest sites and included leaves from fifty-three species sampled at different leaf phenological stages and different leaf light environments. Our resulting spectra-trait models validated on randomly sampled data showed good predictive performance for Vcmax25, Jmax25, Tp25 and Rdark25 (RMSE of 13, 20, 1.5 and 0.3 μmol m-2 s-1, and R2 of 0.74, 0.73, 0.64 and 0.58, respectively). The models showed similar performance when applied to leaves of species not included in the training dataset, illustrating that the approach is robust for capturing the main axes of trait variation in tropical species. We discuss the utility of the spectra-trait and traditional gas exchange approaches for enhancing tropical plant trait studies and improving the parameterization of TBMs.

Leaf mass per area (LMA) is a key plant trait, reflecting tradeoffs between leaf photosynthetic function, longevity, and structural investment. Capturing spatial and temporal variability in LMA has been a long-standing goal of ecological research and is an essential component for advancing Earth system models. Despite the substantial variation in LMA within and across Earth's biomes, an efficient, globally generalizable approach to predict LMA is still lacking. We explored the capacity to predict LMA from leaf spectra across much of the global LMA trait space, with values ranging from 17 to 393 g m-2 . Our dataset contained leaves from a wide range of biomes from the high Arctic to the tropics, included broad- and needleleaf species, and upper- and lower-canopy (i.e. sun and shade) growth environments. Here we demonstrate the capacity to rapidly estimate LMA using only spectral measurements across a wide range of species, leaf age and canopy position from diverse biomes. Our model captures LMA variability with high accuracy and low error (R2 = 0.89; root mean square error (RMSE) = 15.45 g m-2 ). Our finding highlights the fact that the leaf economics spectrum is mirrored by the leaf optical spectrum, paving the way for this technology to predict the diversity of LMA in ecosystems across global biomes.

The triose phosphate utilization (TPU) rate has been identified as one of the processes that can limit terrestrial plant photosynthesis. However, we lack a robust quantitative assessment of TPU limitation of photosynthesis at the global scale. As a result, TPU, and its potential limitation of photosynthesis, is poorly represented in terrestrial biosphere models (TBMs). In this study, we utilized a global data set of photosynthetic CO2 response curves representing 141 species from tropical rainforests to Arctic tundra. We quantified TPU by fitting the standard biochemical model of C3 photosynthesis to measured photosynthetic CO2 response curves and characterized its instantaneous temperature response. Our results demonstrate that TPU does not limit leaf photosynthesis at the current ambient atmospheric CO2 concentration. Furthermore, our results showed that the light-saturated photosynthetic rates of plants growing in cold environments are not more often limited by TPU than those of plants growing in warmer environments. In addition, our study showed that the instantaneous temperature response of TPU is distinct from temperature response of the maximum rate of Rubisco carboxylation. The new formulations of the temperature response of TPU derived in this study may prove useful in quantifying the biochemical limits to terrestrial plant photosynthesis and improve the representation of plant photosynthesis in TBMs.

Understanding the pronounced seasonal and spatial variation in leaf carboxylation capacity (Vc,max ) is critical for determining terrestrial carbon cycling in tropical forests. However, an efficient and scalable approach for predicting Vc,max is still lacking. Here the ability of leaf spectroscopy for rapid estimation of Vc,max was tested. Vc,max was estimated using traditional gas exchange methods, and measured reflectance spectra and leaf age in leaves sampled from tropical forests in Panama and Brazil. These data were used to build a model to predict Vc,max from leaf spectra. The results demonstrated that leaf spectroscopy accurately predicts Vc,max of mature leaves in Panamanian tropical forests (R2 = 0.90). However, this single-age model required recalibration when applied to broader leaf demographic classes (i.e. immature leaves). Combined use of spectroscopy models for Vc,max and leaf age enabled construction of the Vc,max -age relationship solely from leaf spectra, which agreed with field observations. This suggests that the spectroscopy technique can capture the seasonal variability in Vc,max , assuming sufficient sampling across diverse species, leaf ages and canopy environments. This finding will aid development of remote sensing approaches that can be used to characterize Vc,max in moist tropical forests and enable an efficient means to parameterize and evaluate terrestrial biosphere models.

How terrestrial biosphere models (TBMs) represent leaf photosynthesis and its sensitivity to temperature are two critical components of understanding and predicting the response of the Arctic carbon cycle to global change. We measured the effect of temperature on the response of photosynthesis to irradiance in six Arctic plant species and determined the quantum yield of CO2 fixation ( ϕ CO 2 ) and the convexity factor (θ). We also determined leaf absorptance (α) from measured reflectance to calculate ϕ CO 2 on an absorbed light basis ( ϕ CO 2 . a ) and enabled comparison with nine TBMs. The mean ϕ CO 2 . a was 0.045 mol CO2 mol-1 absorbed quanta at 25°C and closely agreed with the mean TBM parameterisation (0.044), but as temperature decreased measured ϕ CO 2 . a diverged from TBMs. At 5°C measured ϕ CO 2 . a was markedly reduced (0.025) and 60% lower than TBM estimates. The θ also showed a significant reduction between 25°C and 5°C. At 5°C θ was 38% lower than the common model parameterisation of 0.7. These data show that TBMs are not accounting for observed reductions in ϕ CO 2 . a and θ that can occur at low temperature. Ignoring these reductions in ϕ CO 2 . a and θ could lead to a marked (45%) overestimation of CO2 assimilation at subsaturating irradiance and low temperature.

Earth system models (ESMs) use photosynthetic capacity, indexed by the maximum Rubisco carboxylation rate (V cmax), to simulate carbon assimilation and typically rely on empirical estimates, including an assumed dependence on leaf nitrogen determined from soil fertility. In contrast, new theory, based on biochemical coordination and co‐optimization of carboxylation and water costs for photosynthesis, suggests that optimal V cmax can be predicted from climate alone, irrespective of soil fertility. Here, we develop this theory and find it captures 64% of observed variability in a global, field‐measured V cmax dataset for C3 plants. Soil fertility indices explained substantially less variation (32%). These results indicate that environmentally regulated biophysical constraints and light availability are the first‐order drivers of global photosynthetic capacity. Through acclimation and adaptation, plants efficiently utilize resources at the leaf level, thus maximizing potential resource use for growth and reproduction. Our theory offers a robust strategy for dynamically predicting photosynthetic capacity in ESMs.

Leaf trait relationships are widely used to predict ecosystem function in terrestrial biosphere models (TBMs), in which leaf maximum carboxylation capacity (Vc,max ), an important trait for modelling photosynthesis, can be inferred from other easier-to-measure traits. However, whether trait-Vc,max relationships are robust across different forest types remains unclear. Here we used measurements of leaf traits, including one morphological trait (leaf mass per area), three biochemical traits (leaf water content, area-based leaf nitrogen content, and leaf chlorophyll content), one physiological trait (Vc,max ), as well as leaf reflectance spectra, and explored their relationships within and across three contrasting forest types in China. We found weak and forest type-specific relationships between Vc,max and the four morphological and biochemical traits (R2 ≤ 0.15), indicated by significantly changing slopes and intercepts across forest types. By contrast, reflectance spectroscopy effectively collapsed the differences in the trait-Vc,max relationships across three forest biomes into a single robust model for Vc,max (R2 = 0.77), and also accurately estimated the four traits (R2 = 0.75-0.94). These findings challenge the traditional use of the empirical trait-Vc,max relationships in TBMs for estimating terrestrial plant photosynthesis, but also highlight spectroscopy as an efficient alternative for characterising Vc,max and multitrait variability, with critical insights into ecosystem modelling and functional trait ecology.