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.

Short-term variation in leaf-level water use efficiency in a tropical forest.

The representation of stomatal regulation of transpiration and CO2 assimilation is key to forecasting terrestrial ecosystem responses to global change. Given its importance in determining the relationship between forest productivity and climate, accurate and mechanistic model representation of the relationship between stomatal conductance (gs ) and assimilation is crucial. We assess possible physiological and mechanistic controls on the estimation of the g1 (stomatal slope, inversely proportional to water use efficiency) and g0 (stomatal intercept) parameters, using diurnal gas exchange surveys and leaf-level response curves of six tropical broadleaf evergreen tree species. g1 estimated from ex situ response curves averaged 50% less than g1 estimated from survey data. While g0 and g1 varied between leaves of different phenological stages, the trend was not consistent among species. We identified a diurnal trend associated with g1 and g0 that significantly improved model projections of diurnal trends in transpiration. The accuracy of modeled gs can be improved by accounting for variation in stomatal behavior across diurnal periods, and between measurement approaches, rather than focusing on phenological variation in stomatal behavior. Additional investigation into the primary mechanisms responsible for diurnal variation in g1 will be required to account for this phenomenon in land-surface models.

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 CO2 and water fluxes. Comparison between observed and TBM trait gradients showed divergence that impacted canopy-scale simulations of water vapor and CO2 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.

Integrating plant physiology into simulation of fire behavior and effects.

Wildfires are a global crisis, but current fire models fail to capture vegetation response to changing climate. With drought and elevated temperature increasing the importance of vegetation dynamics to fire behavior, and the advent of next generation models capable of capturing increasingly complex physical processes, we provide a renewed focus on representation of woody vegetation in fire models. Currently, the most advanced representations of fire behavior and biophysical fire effects are found in distinct classes of fine-scale models and do not capture variation in live fuel (i.e. living plant) properties. We demonstrate that plant water and carbon dynamics, which influence combustion and heat transfer into the plant and often dictate plant survival, provide the mechanistic linkage between fire behavior and effects. Our conceptual framework linking remotely sensed estimates of plant water and carbon to fine-scale models of fire behavior and effects could be a critical first step toward improving the fidelity of the coarse scale models that are now relied upon for global fire forecasting. This process-based approach will be essential to capturing the influence of physiological responses to drought and warming on live fuel conditions, strengthening the science needed to guide fire managers in an uncertain future.

Leaves as bottlenecks: The contribution of tree leaves to hydraulic resistance within the soil-plant-atmosphere continuum.

Within vascular plants, the partitioning of hydraulic resistance along the soil-to-leaf continuum affects transpiration and its response to environmental conditions. In trees, the fractional contribution of leaf hydraulic resistance (Rleaf ) to total soil-to-leaf hydraulic resistance (Rtotal ), or fRleaf (=Rleaf /Rtotal ), is thought to be large, but this has not been tested comprehensively. We compiled a multibiome data set of fRleaf using new and previously published measurements of pressure differences within trees in situ. Across 80 samples, fRleaf averaged 0.51 (95% confidence interval [CI] = 0.46-0.57) and it declined with tree height. We also used the allometric relationship between field-based measurements of soil-to-leaf hydraulic conductance and laboratory-based measurements of leaf hydraulic conductance to compute the average fRleaf for 19 tree samples, which was 0.40 (95% CI = 0.29-0.56). The in situ technique produces a more accurate descriptor of fRleaf because it accounts for dynamic leaf hydraulic conductance. Both approaches demonstrate the outsized role of leaves in controlling tree hydrodynamics. A larger fRleaf may help stems from loss of hydraulic conductance. Thus, the decline in fRleaf with tree height would contribute to greater drought vulnerability in taller trees and potentially to their observed disproportionate drought mortality.

Reducing model uncertainty of climate change impacts on high latitude carbon assimilation.

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.

An improved representation of the relationship between photosynthesis and stomatal conductance leads to more stable estimation of conductance parameters and improves the goodness-of-fit across diverse data sets.

Stomata play a central role in surface-atmosphere exchange by controlling the flux of water and CO2 between the leaf and the atmosphere. Representation of stomatal conductance (gsw ) is therefore an essential component of models that seek to simulate water and CO2 exchange in plants and ecosystems. For given environmental conditions at the leaf surface (CO2 concentration and vapor pressure deficit or relative humidity), models typically assume a linear relationship between gsw and photosynthetic CO2 assimilation (A). However, measurement of leaf-level gsw response curves to changes in A are rare, particularly in the tropics, resulting in only limited data to evaluate this key assumption. Here, we measured the response of gsw and A to irradiance in six tropical species at different leaf phenological stages. We showed that the relationship between gsw and A was not linear, challenging the key assumption upon which optimality theory is based-that the marginal cost of water gain is constant. Our data showed that increasing A resulted in a small increase in gsw at low irradiance, but a much larger increase at high irradiance. We reformulated the popular Unified Stomatal Optimization (USO) model to account for this phenomenon and to enable consistent estimation of the key conductance parameters g0 and g1 . Our modification of the USO model improved the goodness-of-fit and reduced bias, enabling robust estimation of conductance parameters at any irradiance. In addition, our modification revealed previously undetectable relationships between the stomatal slope parameter g1 and other leaf traits. We also observed nonlinear behavior between A and gsw in independent data sets that included data collected from attached and detached leaves, and from plants grown at elevated CO2 concentration. We propose that this empirical modification of the USO model can improve the measurement of gsw parameters and the estimation of plant and ecosystem-scale water and CO2  fluxes.

Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities.

As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiagvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.

Spectroscopy outperforms leaf trait relationships for predicting photosynthetic capacity across different forest types.

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.

Source:sink imbalance detected with leaf- and canopy-level spectroscopy in a field-grown crop.

The finely tuned balance between sources and sinks determines plant resource partitioning and regulates growth and development. Understanding and measuring metabolic indicators of source or sink limitation forms a vital part of global efforts to increase crop yield for future food security. We measured metabolic profiles of Cucurbita pepo (zucchini) grown in the field under carbon sink limitation and control conditions. We demonstrate that these profiles can be measured non-destructively using hyperspectral reflectance at both leaf and canopy scales. Total non-structural carbohydrates (TNC) increased 82% in sink-limited plants; leaf mass per unit area (LMA) increased 38% and free amino acids increased 22%. Partial least-squares regression (PLSR) models link these measured functional traits with reflectance data, enabling high-throughput estimation of traits comprising the sink limitation response. Leaf- and canopy-scale models for TNC had R2 values of 0.93 and 0.64 and %RMSE of 13 and 38%, respectively. For LMA, R2 values were 0.91 and 0.60 and %RMSE 7 and 14%; for free amino acids, R2 was 0.53 and 0.21 with %RMSE 20 and 26%. Remote sensing can enable accurate, rapid detection of sink limitation in the field at the leaf and canopy scale, greatly expanding our ability to understand and measure metabolic responses to stress.

Detection of the metabolic response to drought stress using hyperspectral reflectance.

Drought is the most important limitation on crop yield. Understanding and detecting drought stress in crops is vital for improving water use efficiency through effective breeding and management. Leaf reflectance spectroscopy offers a rapid, non-destructive alternative to traditional techniques for measuring plant traits involved in a drought response. We measured drought stress in six glasshouse-grown agronomic species using physiological, biochemical, and spectral data. In contrast to physiological traits, leaf metabolite concentrations revealed drought stress before it was visible to the naked eye. We used full-spectrum leaf reflectance data to predict metabolite concentrations using partial least-squares regression, with validation R2 values of 0.49-0.87. We show for the first time that spectroscopy may be used for the quantitative estimation of proline and abscisic acid, demonstrating the first use of hyperspectral data to detect a phytohormone. We used linear discriminant analysis and partial least squares discriminant analysis to differentiate between watered plants and those subjected to drought based on measured traits (accuracy: 71%) and raw spectral data (66%). Finally, we validated our glasshouse-developed models in an independent field trial. We demonstrate that spectroscopy can detect drought stress via underlying biochemical changes, before visual differences occur, representing a powerful advance for measuring limitations on yield.

A best-practice guide to predicting plant traits from leaf-level hyperspectral data using partial least squares regression.

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.

Multi-hypothesis comparison of Farquhar and Collatz photosynthesis models reveals the unexpected influence of empirical assumptions at leaf and global scales.

Mechanistic photosynthesis models are at the heart of terrestrial biosphere models (TBMs) simulating the daily, monthly, annual and decadal rhythms of carbon assimilation (A). These models are founded on robust mathematical hypotheses that describe how A responds to changes in light and atmospheric CO2 concentration. Two predominant photosynthesis models are in common usage: Farquhar (FvCB) and Collatz (CBGB). However, a detailed quantitative comparison of these two models has never been undertaken. In this study, we unify the FvCB and CBGB models to a common parameter set and use novel multi-hypothesis methods (that account for both hypothesis and parameter variability) for process-level sensitivity analysis. These models represent three key biological processes: carboxylation, electron transport, triose phosphate use (TPU) and an additional model process: limiting-rate selection. Each of the four processes comprises 1-3 alternative hypotheses giving 12 possible individual models with a total of 14 parameters. To broaden inference, TBM simulations were run and novel, high-resolution photosynthesis measurements were made. We show that parameters associated with carboxylation are the most influential parameters but also reveal the surprising and marked dominance of the limiting-rate selection process (accounting for 57% of the variation in A vs. 22% for carboxylation). The limiting-rate selection assumption proposed by CBGB smooths the transition between limiting rates and always reduces A below the minimum of all potentially limiting rates, by up to 25%, effectively imposing a fourth limitation on A. Evaluation of the CBGB smoothing function in three TBMs demonstrated a reduction in global A by 4%-10%, equivalent to 50%-160% of current annual fossil fuel emissions. This analysis reveals a surprising and previously unquantified influence of a process that has been integral to many TBMs for decades, highlighting the value of multi-hypothesis methods.

Beyond ecosystem modeling: A roadmap to community cyberinfrastructure for ecological data-model integration.

In an era of rapid global change, our ability to understand and predict Earth's natural systems is lagging behind our ability to monitor and measure changes in the biosphere. Bottlenecks to informing models with observations have reduced our capacity to fully exploit the growing volume and variety of available data. Here, we take a critical look at the information infrastructure that connects ecosystem modeling and measurement efforts, and propose a roadmap to community cyberinfrastructure development that can reduce the divisions between empirical research and modeling and accelerate the pace of discovery. A new era of data-model integration requires investment in accessible, scalable, and transparent tools that integrate the expertise of the whole community, including both modelers and empiricists. This roadmap focuses on five key opportunities for community tools: the underlying foundations of community cyberinfrastructure; data ingest; calibration of models to data; model-data benchmarking; and data assimilation and ecological forecasting. This community-driven approach is a key to meeting the pressing needs of science and society in the 21st century.

Leaf traits and canopy structure together explain canopy functional diversity: an airborne remote sensing approach.

Plant functional diversity is strongly connected to photosynthetic carbon assimilation in terrestrial ecosystems. However, many of the plant functional traits that regulate photosynthetic capacity, including foliar nitrogen concentration and leaf mass per area, vary significantly between and within plant functional types and vertically through forest canopies, resulting in considerable landscape-scale heterogeneity in three dimensions. Hyperspectral imagery has been used extensively to quantify functional traits across a range of ecosystems but is generally limited to providing information for top of canopy leaves only. On the other hand, lidar data can be used to retrieve the vertical structure of forest canopies. Because these data are rarely collected at the same time, there are unanswered questions about the effect of forest structure on the three -dimensional spatial patterns of functional traits across ecosystems. In the United States, the National Ecological Observatory Network's Airborne Observation Platform (NEON AOP) provides an opportunity to address this structure-function relationship by collecting lidar and hyperspectral data together across a variety of ecoregions. With a fusion of hyperspectral and lidar data from the NEON AOP and field-collected foliar trait data, we assessed the impacts of forest structure on spatial patterns of N. In addition, we examine the influence of abiotic gradients and management regimes on top-of-canopy percent N and total canopy N (i.e., the total amount of N [g/m2 ] within a forest canopy) at a NEON site consisting of a mosaic of open longleaf pine and dense broadleaf deciduous forests. Our resulting maps suggest that, in contrast to top of canopy values, total canopy N variation is dampened across this landscape resulting in relatively homogeneous spatial patterns. At the same time, we found that leaf functional diversity and canopy structural diversity showed distinct dendritic patterns related to the spatial distribution of plant functional types.

The response of stomatal conductance to seasonal drought in tropical forests.

Stomata regulate CO2 uptake for photosynthesis and water loss through transpiration. The approaches used to represent stomatal conductance (gs ) in models vary. In particular, current understanding of drivers of the variation in a key parameter in those models, the slope parameter (i.e. a measure of intrinsic plant water-use-efficiency), is still limited, particularly in the tropics. Here we collected diurnal measurements of leaf gas exchange and leaf water potential (Ψleaf ), and a suite of plant traits from the upper canopy of 15 tropical trees in two contrasting Panamanian forests throughout the dry season of the 2016 El Niño. The plant traits included wood density, leaf-mass-per-area (LMA), leaf carboxylation capacity (Vc,max ), leaf water content, the degree of isohydry, and predawn Ψleaf . We first investigated how the choice of four commonly used leaf-level gs models with and without the inclusion of Ψleaf as an additional predictor variable influence the ability to predict gs , and then explored the abiotic (i.e. month, site-month interaction) and biotic (i.e. tree-species-specific characteristics) drivers of slope parameter variation. Our results show that the inclusion of Ψleaf did not improve model performance and that the models that represent the response of gs to vapor pressure deficit performed better than corresponding models that respond to relative humidity. Within each gs model, we found large variation in the slope parameter, and this variation was attributable to the biotic driver, rather than abiotic drivers. We further investigated potential relationships between the slope parameter and the six available plant traits mentioned above, and found that only one trait, LMA, had a significant correlation with the slope parameter (R2  = 0.66, n = 15), highlighting a potential path towards improved model parameterization. This study advances understanding of gs dynamics over seasonal drought, and identifies a practical, trait-based approach to improve modeling of carbon and water exchange in tropical forests.

Global photosynthetic capacity is optimized to the environment.

Earth system models (ESMs) use photosynthetic capacity, indexed by the maximum Rubisco carboxylation rate (Vcmax ), 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 Vcmax 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 Vcmax 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.

Terrestrial biosphere models may overestimate Arctic CO2 assimilation if they do not account for decreased quantum yield and convexity at low temperature.

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 ( ϕCO2 ) and the convexity factor (θ). We also determined leaf absorptance (α) from measured reflectance to calculate ϕCO2 on an absorbed light basis ( ϕCO2.a ) and enabled comparison with nine TBMs. The mean ϕCO2.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 ϕCO2.a diverged from TBMs. At 5°C measured ϕCO2.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 ϕCO2.a and θ that can occur at low temperature. Ignoring these reductions in ϕCO2.a and θ could lead to a marked (45%) overestimation of CO2 assimilation at subsaturating irradiance and low temperature.

From the Arctic to the tropics: multibiome prediction of leaf mass per area using leaf reflectance.

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.

Leaf reflectance spectroscopy captures variation in carboxylation capacity across species, canopy environment and leaf age in lowland moist tropical forests.

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.