Convergence in relationships between leaf traits, spectra and age across diverse canopy environments and two contrasting tropical forests

Unveiling the transferability of PLSR models for leaf trait estimation: lessons from a comprehensive analysis with a novel global dataset

Leaf traits are essential for understanding many physiological and ecological processes. Partial least squares regression (PLSR) models with leaf spectroscopy are widely applied for trait estimation, but their transferability across space, time, and plant functional types (PFTs) remains unclear. We compiled a novel dataset of paired leaf traits and spectra, with 47 393 records for > 700 species and eight PFTs at 101 globally distributed locations across multiple seasons. Using this dataset, we conducted an unprecedented comprehensive analysis to assess the transferability of PLSR models in estimating leaf traits. While PLSR models demonstrate commendable performance in predicting chlorophyll content, carotenoid, leaf water, and leaf mass per area prediction within their training data space, their efficacy diminishes when extrapolating to new contexts. Specifically, extrapolating to locations, seasons, and PFTs beyond the training data leads to reduced R2 (0.12-0.49, 0.15-0.42, and 0.25-0.56) and increased NRMSE (3.58-18.24%, 6.27-11.55%, and 7.0-33.12%) compared with nonspatial random cross-validation. The results underscore the importance of incorporating greater spectral diversity in model training to boost its transferability. These findings highlight potential errors in estimating leaf traits across large spatial domains, diverse PFTs, and time due to biased validation schemes, and provide guidance for future field sampling strategies and remote sensing applications.

Foliar spectra accurately distinguish most temperate tree species and show strong phylogenetic signal

PREMISE

Spectroscopy is a powerful remote sensing tool for monitoring plant biodiversity over broad geographic areas. Increasing evidence suggests that foliar spectral reflectance can be used to identify trees at the species level. However, most studies have focused on only a limited number of species at a time, and few studies have explored the underlying phylogenetic structure of leaf spectra. Accurate species identifications are important for reliable estimations of biodiversity from spectral data.

METHODS

Using over 3500 leaf-level spectral measurements, we evaluated whether foliar reflectance spectra (400-2400 nm) can accurately differentiate most tree species from a regional species pool in eastern North America. We explored relationships between spectral, phylogenetic, and leaf functional trait variation as well as their influence on species classification using a hurdle regression model.

RESULTS

Spectral reflectance accurately differentiated tree species (κ = 0.736, ±0.005). Foliar spectra showed strong phylogenetic signal, and classification errors from foliar spectra, although present at higher taxonomic levels, were found predominantly between closely related species, often of the same genus. In addition, we find functional and phylogenetic distance broadly control the occurrence and frequency of spectral classification mistakes among species.

CONCLUSIONS

Our results further support the link between leaf spectral diversity, taxonomic hierarchy, and phylogenetic and functional diversity, and highlight the potential of spectroscopy to remotely sense plant biodiversity and vegetation response to global change.

Innovations in plant science from integrative remote sensing research: an introduction to a Virtual Issue

This article is an Editorial to the Virtual issue on ‘Remote sensing’ that includes the following papers Chavana‐Bryant et al. (2017), Coupel‐Ledru et al. (2022), Cushman & Machado (2020), Disney (2019), D'Odorico et al. (2020), Dong et al. (2022), Fischer et al. (2019), Gamon et al. (2023), Gu et al. (2019), Guillemot et al. (2020), Jucker (2021), Koh et al. (2022), Konings et al. (2019), Kothari et al. (2023), Martini et al. (2022), Richardson (2019), Santini et al. (2021), Schimel et al. (2019), Serbin et al. (2019), Smith et al. (2019, 2020), Still et al. (2021), Stovall et al. (2021), Wang et al. (2020), Wong et al. (2020), Wu et al. (2021), Wu et al. (2017), Wu et al. (2018), Wu et al. (2019), Xu et al. (2021), Yan et al. (2021). Access the Virtual Issue at www.newphytologist.com/virtualissues.

Limiting resource and leaf functional traits jointly determine distribution patterns of leaf intrinsic water use efficiency along aridity gradients

Intrinsic water use efficiency (iWUE) is a critical eco-physiological function allowing plants to adapt to water- and nutrient-limited habitats in arid and semi-arid regions. However, the distribution of iWUE in coexisting species along aridity gradients and its controlling factors are unknown. We established two transects along an aridity gradient in the grasslands of Losses Plateau (LP) and Inner Mongolia Plateau (MP) to elucidate the patterns and underlying mechanisms of iWUE distribution in coexisting species along aridity gradient. We determined leaf carbon (δ¹³C) and oxygen (δ¹⁸O) stable isotopes, functional traits related to carbon fixation, and limiting resources. Bulk leaf δ¹³C and δ¹⁸O were used as proxies for time-integrated iWUE and stomatal conductance (gs) during the growing season. Our results showed that variability in iWUE within transect was primarily controlled by species, sampling sites and an interactive effect between species and sampling sites. Mean values of iWUE (iWUE_Mean) increased and coefficient of variation (CV) in iWUE (iWUE_CV) decreased with an increase in aridity, demonstrating that increases in aridity lead to conservative and convergent water use strategies. Patterns of iWUE_Mean and iWUE_CV were controlled primarily by the ratio of soil organic carbon to total nitrogen in LP and soil moisture in MP. This revealed that the most limited resource drove the distribution patterns of iWUE along aridity gradients. Interspecific variation in iWUE within transect was positively correlated with Δ¹⁸O, indicating that interspecific variation in iWUE was primarily regulated by gs. Furthermore, relationship between iWUE and multi-dimensional functional trait spectrum indicated that species evolved species-specific strategies to adapt to a harsh habitat by partitioning limiting resources. Overall, these findings highlighted the interactive effects of limiting resources and leaf functional traits on plant adaptation strategies for iWUE, and emphasized the importance of considering biological processes in dissecting the underlying mechanisms of plant adaptation strategies at large regional scales.

Identification of Guiboutia species by NIR-HSI spectroscopy

Near infrared hyperspectral imaging (NIR-HSI) spectroscopy can be a rapid, precise, low-cost and non-destructive way for wood identification. In this study, samples of five Guiboutia species were analyzed by means of NIR-HSI. Partial least squares discriminant analysis (PLS-DA) and support vector machine (SVM) were used after different data treatment in order to improve the performance of models. Transverse, radial, and tangential section were analyzed separately to select the best sample section for wood identification. The results obtained demonstrated that NIR-HSI combined with successive projections algorithm (SPA) and SVM can achieve high prediction accuracy and low computing cost. Pre-processing methods of SNV and Normalize can increase the prediction accuracy slightly, however, high modelling accuracy can still be achieved by raw pre-processing. Both models for the classification of G. conjugate, G. ehie and G. demeusei perform nearly 100% accuracy. Prediction for G. coleosperma and G. tessmannii were more difficult when using PLS-DA model. It is evidently clear from the findings that the transverse section of wood is more suitable for wood identification. NIR-HSI spectroscopy technique has great potential for Guiboutia species analysis.

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 CO₂ between the leaf and the atmosphere. Representation of stomatal conductance (g_sw ) is therefore an essential component of models that seek to simulate water and CO₂ exchange in plants and ecosystems. For given environmental conditions at the leaf surface (CO₂ concentration and vapor pressure deficit or relative humidity), models typically assume a linear relationship between g_sw and photosynthetic CO₂ assimilation (A). However, measurement of leaf-level g_sw 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 g_sw and A to irradiance in six tropical species at different leaf phenological stages. We showed that the relationship between g_sw 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 g_sw 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 g₀ and g₁ . 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 g₁ and other leaf traits. We also observed nonlinear behavior between A and g_sw in independent data sets that included data collected from attached and detached leaves, and from plants grown at elevated CO₂ concentration. We propose that this empirical modification of the USO model can improve the measurement of g_sw parameters and the estimation of plant and ecosystem-scale water and CO₂  fluxes.

A Perspective on Plant Phenomics: Coupling Deep Learning and Near-Infrared Spectroscopy

The trait-based approach in plant ecology aims at understanding and classifying the diversity of ecological strategies by comparing plant morphology and physiology across organisms. The major drawback of the approach is that the time and financial cost of measuring the traits on many individuals and environments can be prohibitive. We show that combining near-infrared spectroscopy (NIRS) with deep learning resolves this limitation by quickly, non-destructively, and accurately measuring a suite of traits, including plant morphology, chemistry, and metabolism. Such an approach also allows to position plants within the well-known CSR triangle that depicts the diversity of plant ecological strategies. The processing of NIRS through deep learning identifies the effect of growth conditions on trait values, an issue that plagues traditional statistical approaches. Together, the coupling of NIRS and deep learning is a promising high-throughput approach to capture a range of ecological information on plant diversity and functioning and can accelerate the creation of extensive trait databases.

Changes in leaf functional traits with leaf age: when do leaves decrease their photosynthetic capacity in Amazonian trees?

Most leaf functional trait studies in the Amazon basin do not consider ontogenetic variations (leaf age), which may influence ecosystem productivity throughout the year. When leaf age is taken into account, it is generally considered discontinuous, and leaves are classified into age categories based on qualitative observations. Here, we quantified age-dependent changes in leaf functional traits such as the maximum carboxylation rate of ribulose-1,5-biphosphate carboxylase/oxygenase (Rubisco) (Vcmax), stomatal control (Cgs%), leaf dry mass per area and leaf macronutrient concentrations for nine naturally growing Amazon tropical trees with variable phenological strategies. Leaf ages were assessed by monthly censuses of branch-level leaf demography; we also performed leaf trait measurements accounting for leaf chronological age based on days elapsed since the first inclusion in the leaf demography, not predetermined age classes. At the tree community scale, a nonlinear relationship between Vcmax and leaf age existed: young, developing leaves showed the lowest mean photosynthetic capacity, increasing to a maximum at 45 days and then decreasing gradually with age in both continuous and categorical age group analyses. Maturation times among species and phenological habits differed substantially, from 8 ± 30 to 238 ± 30 days, and the rate of decline of Vcmax varied from -0.003 to -0.065 μmol CO2 m-2 s-1 day-1. Stomatal control increased significantly in young leaves but remained constant after peaking. Mass-based phosphorus and potassium concentrations displayed negative relationships with leaf age, whereas nitrogen did not vary temporally. Differences in life strategies, leaf nutrient concentrations and phenological types, not the leaf age effect alone, may thus be important factors for understanding observed photosynthesis seasonality in Amazonian forests. Furthermore, assigning leaf age categories in diverse tree communities may not be the recommended method for studying carbon uptake seasonality in the Amazon, since the relationship between Vcmax and leaf age could not be confirmed for all trees.

Linking Land Use and Plant Functional Diversity Patterns in Sabah, Borneo, through Large-Scale Spatially Continuous Sentinel-2 Inference

Global biodiversity losses erode the functioning of our vital ecosystems. Functional diversity is increasingly recognized as a critical link between biodiversity and ecosystem functioning. Satellite earth observation was proposed to address the current absence of information on large-scale continuous patterns of plant functional diversity. This study demonstrates the inference and spatial mapping of functional diversity metrics through satellite remote sensing over a large key biodiversity region (Sabah, Malaysian Borneo, ~53,000 km2) and compares the derived estimates across a land-use gradient as an initial qualitative assessment to test the potential merits of the approach. Functional traits (leaf water content, chlorophyll-a and -b, and leaf area index) were estimated from Sentinel-2 spectral reflectance using a pre-trained neural network on radiative transfer modeling simulations. Multivariate functional diversity metrics were calculated, including functional richness, divergence, and evenness. Spatial patterns of functional diversity were related to land-use data distinguishing intact forest, logged forest, and oil palm plantations. Spatial patterns of satellite remotely sensed functional diversity are significantly related to differences in land use. Intact forests, as well as logged forests, featured consistently higher functional diversity compared to oil palm plantations. Differences were profound for functional divergence, whereas functional richness exhibited relatively large variances within land-use classes. By linking large-scale patterns of functional diversity as derived from satellite remote sensing to land-use information, this study indicated initial responsiveness to broad human disturbance gradients over large geographical and spatially contiguous extents. Despite uncertainties about the accuracy of the spatial patterns, this study provides a coherent early application of satellite-derived functional diversity toward further validation of its responsiveness across ecological gradients.

Rapid estimation of photosynthetic leaf traits of tropical plants in diverse environmental conditions using reflectance spectroscopy

Tropical forests are one of the main carbon sinks on Earth, but the magnitude of CO₂ absorbed by tropical vegetation remains uncertain. Terrestrial biosphere models (TBMs) are commonly used to estimate the CO₂ absorbed by forests, but their performance is highly sensitive to the parameterization of processes that control leaf-level CO₂ exchange. Direct measurements of leaf respiratory and photosynthetic traits that determine vegetation CO₂ 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 (R_dark25), the maximum rate of carboxylation by the enzyme Rubisco (V_cmax25), the maximum rate of electron transport (J_max25), and the triose phosphate utilization rate (T_p25), 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 V_cmax25, J_max25, T_p25 and R_dark25 (RMSE of 13, 20, 1.5 and 0.3 μmol m^-2 s^-1, and R² 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.

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.

Dynamic iridescent signals of male copperwing damselflies coupled with wing-clapping displays: the perspective of different receivers

Dynamic signals are a widespread phenomenon in several taxa, usually associated with intraspecific communication. In contrast, dynamic iridescent signals are detectable only at specific angles of illumination; hence, the animal can hide the signal to avoid detection when necessary. This structural coloration is mostly dependent on the illumination, the contrast against the background and the vision of the receiver. Complex behavioural displays can be coupled with structural coloration to create dynamic visual signals that enhance these functions. Here, we address whether iridescence of the males of a damselfly that inhabits dark rainforests, Chalcopteryx scintillans, can be considered a dynamic visual signal. We analyse whether coloration is perceived by conspecifics, while reducing detectability to eavesdroppers against three types of backgrounds. Our results suggest that the visual background affects the detectability of male hindwings by different receivers, mostly predators and prey. We discuss whether these results and the angle dependence of colour could indicate a mechanism to avoid unwanted intraspecific interactions or even to lure both predators and prey. We conclude that the main functions of the dynamic iridescent signal are to communicate with conspecifics while hindering the signal for prey, adding evidence of the multifunctionality of structural coloration coupled with behavioural displays in animals.

Leaf Age Matters in Remote Sensing: Taking Ground Truth for Spectroscopic Studies in Hemiboreal Deciduous Trees with Continuous Leaf Formation

We examined the seasonal changes in biophysical, anatomical, and optical traits of young leaves, formed throughout the vegetative season due to sylleptic growth, and mature leaves formed by proleptic growth in spring. Leaf developmental categories contribute to the top-of-canopy reflectance and should be considered when taking ground truth for remote sensing studies (RS). Deciduous tree species, Betula pendula, Populus tremula, and Alnus incana, were sampled from May to October 2018 in an Estonian hemiboreal forest. Chlorophyll and carotenoid content were detected biochemically; leaf anatomical traits (leaf, palisade, and spongy mesophyll thickness) were measured on leaf cross-sections; leaf reflectance was measured by a spectroradiometer with an integrating sphere (350–2500 nm). Biophysical and anatomical leaf traits were related to 64 vegetation indices (VIs). Linear models based on VIs for all tested leaf traits were more robust if both juvenile and mature leaves were included. This study provides information on which VIs are interchangeable or independent. Pigment and leaf thickness sensitive indices formed PC1; water and structural trait related VIs formed an independent group associated with PC3. Type of growth and leaf age could affect the validation of biophysical and anatomical leaf trait retrieval from the optical signal. It is, therefore, necessary to sample both leaf developmental categories—young and mature—in RS, especially if sampling is only once within the vegetation season.

Evaluating the effect of stand properties and site conditions on the forest reflectance from Sentinel-2 time series

Forest stand reflectance at the canopy level results from various factors, such as vegetation chemical properties, leaf morphology, canopy structure, and tree sizes. These factors are dependent on the species, age, and health statuses of trees, as well as the site conditions. Sentinel-2 imagery with the high spatial, spectral, and temporal resolution, has enabled analysis of the relationships between vegetation properties and their spectral responses at large spatial scales. A comprehensive study of these relationships is needed to understand the drivers of vegetation spectral patterns and is essential from the point of view of remote sensing data interpretation. Our study aimed to quantify the site and forest parameters affecting forest stands reflectance. The analysis was conducted for common beech-, silver fir- and Scots pine-dominated stands in a mountainous area of the Polish Carpathians. The effect of stands and site properties on reflectance in different parts of the growing season was captured using the dense time series provided by Sentinel-2 from 2018-2019. The results indicate that the reflectance of common beech stands is mainly influenced by elevation, particularly during spring and autumn. Other factors influencing beech stand reflectance include the share of the broadleaved understory, aspect, and, during summer, the age of stands. The reflectance of coniferous species, i.e., Scots pine and silver fir, is mainly influenced by the age and stand properties, namely the crown closure and stand density. The age is a primary driver for silver fir stands reflectance changes, while the stand properties have a large impact on Scots pine stands reflectance. Also, the understory influences Scots pine stands reflectance, while there appears to be no impact on silver fir stands. The influence of the abovementioned factors is highly diverse, depending on the used band and time of the season.

Exploring the Influence of Biological Traits and Environmental Drivers on Water Use Variations across Contrasting Forests

Understanding species-specific water use patterns across contrasting sites and how sensitivity of responses to environmental variables changes for different species is critical for evaluating potential forest dynamics and land use changes under global change. To quantify water use patterns and the sensitivity of tree transpiration to environmental drivers among sites and species, sap flow and meteorological data sets from three contrasting climatic zones were combined and compared in this analysis. Agathis australis from NZHP site, Schima wallichii Choisy (native) and Acacia mangium Willd (exotic) from CHS site, Liquidamber formosana Hance, Quercus variabilis Blume and Quercus acutissima Carruth from CJGS site were the dominant trees chosen as our study species. Biological traits were collected to explain the underlying physiological mechanisms for water use variation. Results showed that the strongest environmental drivers of sap flow were photosynthetically active radiation (PAR), vapor pressure deficit (VPD) and temperature across sites, indicating that the response of water use to abiotic drivers converged across sites. Water use magnitude was site specific, which was controlled by site characteristics, species composition and local weather conditions. The species with higher sap flow density (Fd) generally had greater stomatal conductance. Native deciduous broadleaved species had a higher Fd and faster response to stomatal regulation than that of native evergreen broadleaved species (S. wallichii) and conifer species A. australis. The analysis also showed that exotic species (A. mangium) consumed more water than native species (S. wallichii). Trees with diffuse porous and lower wood density had relatively higher Fd for angiosperms, suggesting that water use was regulated by physiological differences. Water use characteristics across sites are controlled by both external factors such as site-specific characteristics (local environmental conditions and species composition) and internal factors such as biological traits (xylem anatomy, root biomass and leaf area), which highlights the complexity of quantifying land water budgets for areas covered by different species.

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.

Leaf reflectance spectra capture the evolutionary history of seed plants

Leaf reflectance spectra have been increasingly used to assess plant diversity. However, we do not yet understand how spectra vary across the tree of life or how the evolution of leaf traits affects the differentiation of spectra among species and lineages. Here we describe a framework that integrates spectra with phylogenies and apply it to a global dataset of over 16 000 leaf-level spectra (400-2400 nm) for 544 seed plant species. We test for phylogenetic signal in spectra, evaluate their ability to classify lineages, and characterize their evolutionary dynamics. We show that phylogenetic signal is present in leaf spectra but that the spectral regions most strongly associated with the phylogeny vary among lineages. Despite among-lineage heterogeneity, broad plant groups, orders, and families can be identified from reflectance spectra. Evolutionary models also reveal that different spectral regions evolve at different rates and under different constraint levels, mirroring the evolution of their underlying traits. Leaf spectra capture the phylogenetic history of seed plants and the evolutionary dynamics of leaf chemistry and structure. Consequently, spectra have the potential to provide breakthrough assessments of leaf evolution and plant phylogenetic diversity at global scales.

Hyperspectral and Thermal Sensing of Stomatal Conductance, Transpiration, and Photosynthesis for Soybean and Maize under Drought

During water stress, crops undertake adjustments in functional, structural, and biochemical traits. Hyperspectral data and machine learning techniques (PLS-R) can be used to assess water stress responses in plant physiology. In this study, we investigated the potential of hyperspectral optical (VNIR) measurements supplemented with thermal remote sensing and canopy height (hc) to detect changes in leaf physiology of soybean (C3) and maize (C4) plants under three levels of soil moisture in controlled environmental conditions. We measured canopy evapotranspiration (ET), leaf transpiration (Tr), leaf stomatal conductance (gs), leaf photosynthesis (A), leaf chlorophyll content and morphological properties (hc and LAI), as well as vegetation cover reflectance and radiometric temperature (TL,Rad). Our results showed that water stress caused significant ET decreases in both crops. This reduction was linked to tighter stomatal control for soybean plants, whereas LAI changes were the primary control on maize ET. Spectral vegetation indices (VIs) and TL,Rad were able to track these different responses to drought, but only after controlling for confounding changes in phenology. PLS-R modeling of gs, Tr, and A using hyperspectral data was more accurate when pooling data from both crops together rather than individually. Nonetheless, separated PLS-R crop models are useful to identify the most relevant variables in each crop such as TL,Rad for soybean and hc for maize under our experimental conditions. Interestingly, the most important spectral bands sensitive to drought, derived from PLS-R analysis, were not exactly centered at the same wavelengths of the studied VIs sensitive to drought, highlighting the benefit of having contiguous narrow spectral bands to predict leaf physiology and suggesting different wavelength combinations based on crop type. Our results are only a first but a promising step towards larger scale remote sensing applications (e.g., airborne and satellite). PLS-R estimates of leaf physiology could help to parameterize canopy level GPP or ET models and to identify different photosynthetic paths or the degree of stomatal closure in response to drought.

Does the leaf economic spectrum hold within plant functional types? A Bayesian multivariate trait meta-analysis

The leaf economic spectrum is a widely studied axis of plant trait variability that defines a trade-off between leaf longevity and productivity. While this has been investigated at the global scale, where it is robust, and at local scales, where deviations from it are common, it has received less attention at the intermediate scale of plant functional types (PFTs). We investigated whether global leaf economic relationships are also present within the scale of plant functional types (PFTs) commonly used by Earth System models, and the extent to which this global-PFT hierarchy can be used to constrain trait estimates. We developed a hierarchical multivariate Bayesian model that assumes separate means and covariance structures within and across PFTs and fit this model to seven leaf traits from the TRY database related to leaf longevity, morphology, biochemistry, and photosynthetic metabolism. Although patterns of trait covariation were generally consistent with the leaf economic spectrum, we found three approximate tiers to this consistency. Relationships among morphological and biochemical traits (specific leaf area [SLA], N, P) were the most robust within and across PFTs, suggesting that covariation in these traits is driven by universal leaf construction trade-offs and stoichiometry. Relationships among metabolic traits (dark respiration [R_d ], maximum RuBisCo carboxylation rate [V_c,max ], maximum electron transport rate [J_max ]) were slightly less consistent, reflecting in part their much sparser sampling (especially for high-latitude PFTs), but also pointing to more flexible plasticity in plant metabolistm. Finally, relationships involving leaf lifespan were the least consistent, indicating that leaf economic relationships related to leaf lifespan are dominated by across-PFT differences and that within-PFT variation in leaf lifespan is more complex and idiosyncratic. Across all traits, this covariance was an important source of information, as evidenced by the improved imputation accuracy and reduced predictive uncertainty in multivariate models compared to univariate models. Ultimately, our study reaffirms the value of studying not just individual traits but the multivariate trait space and the utility of hierarchical modeling for studying the scale dependence of trait relationships.

Estimation of leaf area index using PROSAIL based LUT inversion, MLRA-GPR and empirical models: Case study of tropical deciduous forest plantation, North India

Forests play a vital role in biological cycles and environmental regulation. To understand the key processes of forest canopies (e.g., photosynthesis, respiration and transpiration), reliable and accurate information on spatial variability of Leaf Area Index (LAI), and its seasonal dynamics is essential. In the present study, we assessed the performance of biophysical parameter (LAI) retrieval methods viz. Look-Up Table (LUT)-inversion, MLRA-GPR (Machine Learning Regression Algorithm-Gaussian Processes Regression) and empirical models, for estimating the LAI of tropical deciduous plantation using ARTMO (Automated Radiative Transfer Models Operator) tool and Sentinel-2 satellite images. The study was conducted in Central Tarai Forest Division, Haldwani, located in the Uttarakhand state, India. A total of 49 ESUs (Elementary Sampling Unit) of 30m×30m size were established based on variability in composition and age of plantation stands. In-situ LAI was recorded using plant canopy imager during the leaf growing, peak and senescence seasons. The PROSAIL model was calibrated with site-specific biophysical and biochemical parameters before used to the predicted LAI. The plantation LAI was also predicted by an empirical approach using optimally chosen Sentinel-2 vegetation indices. In addition, Sentinel-2 and MODIS LAI products were evaluated with respect to LAI measurements. MLRA-GPR offered best results for predicting LAI of leaf growing (R2 = 0.9, RMSE = 0.14), peak (R2 = 0.87, RMSE = 0.21) and senescence (R2 = 0.86, RMSE = 0.31) seasons while LUT inverted model outperformed VI's based parametric regression model. Vegetation indices (VIs) derived from 740 nm, 783 nm and 2190 nm band combinations of Sentinel-2 offered the best prediction of LAI.

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 (R²  = 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.

Cryptic phenology in plants: Case studies, implications, and recommendations

Plant phenology-the timing of cyclic or recurrent biological events in plants-offers insight into the ecology, evolution, and seasonality of plant-mediated ecosystem processes. Traditionally studied phenologies are readily apparent, such as flowering events, germination timing, and season-initiating budbreak. However, a broad range of phenologies that are fundamental to the ecology and evolution of plants, and to global biogeochemical cycles and climate change predictions, have been neglected because they are "cryptic"-that is, hidden from view (e.g., root production) or difficult to distinguish and interpret based on common measurements at typical scales of examination (e.g., leaf turnover in evergreen forests). We illustrate how capturing cryptic phenology can advance scientific understanding with two case studies: wood phenology in a deciduous forest of the northeastern USA and leaf phenology in tropical evergreen forests of Amazonia. Drawing on these case studies and other literature, we argue that conceptualizing and characterizing cryptic plant phenology is needed for understanding and accurate prediction at many scales from organisms to ecosystems. We recommend avenues of empirical and modeling research to accelerate discovery of cryptic phenological patterns, to understand their causes and consequences, and to represent these processes in terrestrial biosphere models.

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 (V_c,max ) is critical for determining terrestrial carbon cycling in tropical forests. However, an efficient and scalable approach for predicting V_c,max is still lacking. Here the ability of leaf spectroscopy for rapid estimation of V_c,max was tested. V_c,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 V_c,max from leaf spectra. The results demonstrated that leaf spectroscopy accurately predicts V_c,max of mature leaves in Panamanian tropical forests (R²  = 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 V_c,max and leaf age enabled construction of the V_c,max -age relationship solely from leaf spectra, which agreed with field observations. This suggests that the spectroscopy technique can capture the seasonal variability in V_c,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 V_c,max in moist tropical forests and enable an efficient means to parameterize and evaluate terrestrial biosphere models.

High-throughput field phenotyping using hyperspectral reflectance and partial least squares regression (PLSR) reveals genetic modifications to photosynthetic capacity

Spectroscopy is becoming an increasingly powerful tool to alleviate the challenges of traditional measurements of key plant traits at the leaf, canopy, and ecosystem scales. Spectroscopic methods often rely on statistical approaches to reduce data redundancy and enhance useful prediction of physiological traits. Given the mechanistic uncertainty of spectroscopic techniques, genetic modification of plant biochemical pathways may affect reflectance spectra causing predictive models to lose power. The objectives of this research were to assess over two separate years, whether a predictive model can represent natural and imposed variation in leaf photosynthetic potential for different crop cultivars and genetically modified plants, to assess the interannual capabilities of a partial least square regression (PLSR) model, and to determine whether leaf N is a dominant driver of photosynthesis in PLSR models. In 2016, a PLSR analysis of reflectance spectra coupled with gas exchange data was used to build predictive models for photosynthetic parameters including maximum carboxylation rate of Rubisco (V _c,max ), maximum electron transport rate (J _max ) and percentage leaf nitrogen ([N]). The model was developed for wild type and genetically modified plants that represent a wide range of photosynthetic capacities. Results show that hyperspectral reflectance accurately predicted V _c,max, J _max and [N] for all plants measured in 2016. Applying these PLSR models to plants grown in 2017 resulted in a strong predictive ability relative to gas exchange measurements for V _c,max, but not for J _max, and not for genotypes unique to 2017. Building a new model including data collected in 2017 resulted in more robust predictions, with R² increases of 17% for V _c,max . and 13% J _max . Plants generally have a positive correlation between leaf nitrogen and photosynthesis, however, tobacco with reduced Rubisco (SSuD) had significantly higher [N] despite much lower V _c,max. The PLSR model was able to accurately predict both lower V _c,max and higher leaf [N] for this genotype suggesting that the spectral based estimates of V _c,max and leaf nitrogen [N] are independent. These results suggest that the PLSR model can be applied across years, but only to genotypes used to build the model and that the actual mechanism measured with the PLSR technique is not directly related to leaf [N]. The success of the leaf-scale analysis suggests that similar approaches may be successful at the canopy and ecosystem scales but to use these methods across years and between genotypes at any scale, application of accurately populated physical based models based on radiative transfer principles may be required.

Quantifying Leaf Phenology of Individual Trees and Species in a Tropical Forest Using Unmanned Aerial Vehicle (UAV) Images

Tropical forests exhibit complex but poorly understood patterns of leaf phenology. Understanding species- and individual-level phenological patterns in tropical forests requires datasets covering large numbers of trees, which can be provided by Unmanned Aerial Vehicles (UAVs). In this paper, we test a workflow combining high-resolution RGB images (7 cm/pixel) acquired from UAVs with a machine learning algorithm to monitor tree and species leaf phenology in a tropical forest in Panama. We acquired images for 34 flight dates over a 12-month period. Crown boundaries were digitized in images and linked with forest inventory data to identify species. We evaluated predictions of leaf cover from different models that included up to 14 image features extracted for each crown on each date. The models were trained and tested with visual estimates of leaf cover from 2422 images from 85 crowns belonging to eight species spanning a range of phenological patterns. The best-performing model included both standard color metrics, as well as texture metrics that quantify within-crown variation, with r2 of 0.84 and mean absolute error (MAE) of 7.8% in 10-fold cross-validation. In contrast, the model based only on the widely-used Green Chromatic Coordinate (GCC) index performed relatively poorly (r2 = 0.52, MAE = 13.6%). These results highlight the utility of texture features for image analysis of tropical forest canopies, where illumination changes may diminish the utility of color indices, such as GCC. The algorithm successfully predicted both individual-tree and species patterns, with mean r2 of 0.82 and 0.89 and mean MAE of 8.1% and 6.0% for individual- and species-level analyses, respectively. Our study is the first to develop and test methods for landscape-scale UAV monitoring of individual trees and species in diverse tropical forests. Our analyses revealed undescribed patterns of high intraspecific variation and complex leaf cover changes for some species.

Leaf age effects on the spectral predictability of leaf traits in Amazonian canopy trees

Recent work has shown that leaf traits and spectral properties change through time and/or seasonally as leaves age. Current field and hyperspectral methods used to estimate canopy leaf traits could, therefore, be significantly biased by variation in leaf age. To explore the magnitude of this effect, we used a phenological dataset comprised of leaves of different leaf age groups -developmental, mature, senescent and mixed-age- from canopy and emergent tropical trees in southern Peru. We tested the performance of partial least squares regression models developed from these different age groups when predicting traits for leaves of different ages on both a mass and area basis. Overall, area-based models outperformed mass-based models with a striking improvement in prediction observed for area-based leaf carbon (C_area) estimates. We observed trait-specific age effects in all mass-based models while area-based models displayed age effects in mixed-age leaf groups for P_area and N_area. Spectral coefficients and variable importance in projection (VIPs) also reflected age effects. Both mass- and area-based models for all five leaf traits displayed age/temporal sensitivity when we tested their ability to predict the traits of leaves of other age groups. Importantly, mass-based mature models displayed the worst overall performance when predicting the traits of leaves from other age groups. These results indicate that the widely adopted approach of using fully expanded mature leaves to calibrate models that estimate remotely-sensed tree canopy traits introduces error that can bias results depending on the phenological stage of canopy leaves. To achieve temporally stable models, spectroscopic studies should consider producing area-based estimates as well as calibrating models with leaves of different age groups as they present themselves through the growing season. We discuss the implications of this for surveys of canopies with synchronised and unsynchronised leaf phenology.

Mapping the Leaf Economic Spectrum across West African Tropical Forests Using UAV-Acquired Hyperspectral Imagery

The leaf economic spectrum (LES) describes a set of universal trade-offs between leaf mass per area (LMA), leaf nitrogen (N), leaf phosphorus (P) and leaf photosynthesis that influence patterns of primary productivity and nutrient cycling. Many questions regarding vegetation-climate feedbacks can be addressed with a better understanding of LES traits and their controls. Remote sensing offers enormous potential for generating large-scale LES trait data. Yet so far, canopy studies have been limited to imaging spectrometers onboard aircraft, which are rare, expensive to deploy and lack fine-scale resolution. In this study, we measured VNIR (visible-near infrared (400–1050 nm)) reflectance of individual sun and shade leaves in 7 one-ha tropical forest plots located along a 1200–2000 mm precipitation gradient in West Africa. We collected hyperspectral imaging data from 3 of the 7 plots, using an octocopter-based unmanned aerial vehicle (UAV), mounted with a hyperspectral mapping system (450–950 nm, 9 nm FWHM). Using partial least squares regression (PLSR), we found that the spectra of individual sun leaves demonstrated significant (p < 0.01) correlations with LMA and leaf chemical traits: r2 = 0.42 (LMA), r2 = 0.43 (N), r2 = 0.21 (P), r2 = 0.20 (leaf potassium (K)), r2 = 0.23 (leaf calcium (Ca)) and r2 = 0.14 (leaf magnesium (Mg)). Shade leaf spectra displayed stronger relationships with all leaf traits. At the airborne level, four of the six leaf traits demonstrated weak (p < 0.10) correlations with the UAV-collected spectra of 58 tree crowns: r2 = 0.25 (LMA), r2 = 0.22 (N), r2 = 0.22 (P), and r2 = 0.25 (Ca). From the airborne imaging data, we used LMA, N and P values to map the LES across the three plots, revealing precipitation and substrate as co-dominant drivers of trait distributions and relationships. Positive N-P correlations and LMA-P anticorrelations followed typical LES theory, but we found no classic trade-offs between LMA and N. Overall, this study demonstrates the application of UAVs to generating LES information and advancing the study and monitoring tropical forest functional diversity.

Age-dependent leaf physiology and consequences for crown-scale carbon uptake during the dry season in an Amazon evergreen forest

Satellite and tower-based metrics of forest-scale photosynthesis generally increase with dry season progression across central Amazônia, but the underlying mechanisms lack consensus. We conducted demographic surveys of leaf age composition, and measured the age dependence of leaf physiology in broadleaf canopy trees of abundant species at a central eastern Amazon site. Using a novel leaf-to-branch scaling approach, we used these data to independently test the much-debated hypothesis - arising from satellite and tower-based observations - that leaf phenology could explain the forest-scale pattern of dry season photosynthesis. Stomatal conductance and biochemical parameters of photosynthesis were higher for recently mature leaves than for old leaves. Most branches had multiple leaf age categories simultaneously present, and the number of recently mature leaves increased as the dry season progressed because old leaves were exchanged for new leaves. These findings provide the first direct field evidence that branch-scale photosynthetic capacity increases during the dry season, with a magnitude consistent with increases in ecosystem-scale photosynthetic capacity derived from flux towers. Interactions between leaf age-dependent physiology and shifting leaf age-demographic composition are sufficient to explain the dry season photosynthetic capacity pattern at this site, and should be considered in vegetation models of tropical evergreen forests.

Near-infrared spectrometry allows fast and extensive predictions of functional traits from dry leaves and branches

The linking of individual functional traits to ecosystem processes is the basis for making generalizations in ecology, but the measurement of individual values is laborious and time consuming, preventing large-scale trait mapping. Also, in hyper-diverse systems, errors occur because identification is difficult, and species level values ignore intra-specific variation. To allow extensive trait mapping at the individual level, we evaluated the potential of Fourrier-Transformed Near Infra-Red Spectrometry (FT-NIR) to adequately describe 14 traits that are key for plant carbon, water, and nutrient balance. FT-NIR absorption spectra (1,000-2,500 nm) were obtained from dry leaves and branches of 1,324 trees of 432 species from a hyper-diverse Amazonian forest. FT-NIR spectra were related to measured traits for the same plants using partial least squares regressions. A further 80 plants were collected from a different site to evaluate model applicability across sites. Relative prediction error (RMSE_rel ) was calculated as the percentage of the trait value range represented by the final model RMSE. The key traits used in most functional trait studies; specific leaf area, leaf dry matter content, wood density and wood dry matter content can be well predicted by the model (R²  = 0.69-0.78, RMSE_rel  = 9-11%), while leaf density, xylem proportion, bark density and bark dry matter content can be moderately well predicted (R²  = 0.53-0.61, RMSE_rel  = 14-17%). Community-weighted means of all traits were well estimated with NIR, as did the shape of the frequency distribution of the community values for the above key traits. The model developed at the core site provided good estimations of the key traits of a different site. An evaluation of the sampling effort indicated that 400 or less individuals may be sufficient for establishing a good local model. We conclude that FT-NIR is an easy, fast and cheap method for the large-scale estimation of individual plant traits that was previously impossible. The ability to use dry intact leaves and branches unlocks the potential for using herbarium material to estimate functional traits; thus advancing our knowledge of community and ecosystem functioning from local to global scales.

Biological processes dominate seasonality of remotely sensed canopy greenness in an Amazon evergreen forest

Satellite observations of Amazon forests show seasonal and interannual variations, but the underlying biological processes remain debated. Here we combined radiative transfer models (RTMs) with field observations of Amazon forest leaf and canopy characteristics to test three hypotheses for satellite-observed canopy reflectance seasonality: seasonal changes in leaf area index, in canopy-surface leafless crown fraction and/or in leaf demography. Canopy RTMs (PROSAIL and FLiES), driven by these three factors combined, simulated satellite-observed seasonal patterns well, explaining c. 70% of the variability in a key reflectance-based vegetation index (MAIAC EVI, which removes artifacts that would otherwise arise from clouds/aerosols and sun-sensor geometry). Leaf area index, leafless crown fraction and leaf demography independently accounted for 1, 33 and 66% of FLiES-simulated EVI seasonality, respectively. These factors also strongly influenced modeled near-infrared (NIR) reflectance, explaining why both modeled and observed EVI, which is especially sensitive to NIR, captures canopy seasonal dynamics well. Our improved analysis of canopy-scale biophysics rules out satellite artifacts as significant causes of satellite-observed seasonal patterns at this site, implying that aggregated phenology explains the larger scale remotely observed patterns. This work significantly reconciles current controversies about satellite-detected Amazon phenology, and improves our use of satellite observations to study climate-phenology relationships in the tropics.

The phenology of leaf quality and its within-canopy variation is essential for accurate modeling of photosynthesis in tropical evergreen forests

Leaf quantity (i.e., canopy leaf area index, LAI), quality (i.e., per-area photosynthetic capacity), and longevity all influence the photosynthetic seasonality of tropical evergreen forests. However, these components of tropical leaf phenology are poorly represented in most terrestrial biosphere models (TBMs). Here, we explored alternative options for the representation of leaf phenology effects in TBMs that employ the Farquahar, von Caemmerer & Berry (FvCB) representation of CO₂ assimilation. We developed a two-fraction leaf (sun and shade), two-layer canopy (upper and lower) photosynthesis model to evaluate different modeling approaches and assessed three components of phenological variations (i.e., leaf quantity, quality, and within-canopy variation in leaf longevity). Our model was driven by the prescribed seasonality of leaf quantity and quality derived from ground-based measurements within an Amazonian evergreen forest. Modeled photosynthetic seasonality was not sensitive to leaf quantity, but was highly sensitive to leaf quality and its vertical distribution within the canopy, with markedly more sensitivity to upper canopy leaf quality. This is because light absorption in tropical canopies is near maximal for the entire year, implying that seasonal changes in LAI have little impact on total canopy light absorption; and because leaf quality has a greater effect on photosynthesis of sunlit leaves than light limited, shade leaves and sunlit foliage are more abundant in the upper canopy. Our two-fraction leaf, two-layer canopy model, which accounted for all three phenological components, was able to simulate photosynthetic seasonality, explaining ~90% of the average seasonal variation in eddy covariance-derived CO₂ assimilation. This work identifies a parsimonious approach for representing tropical evergreen forest photosynthetic seasonality in TBMs that utilize the FvCB model of CO₂ assimilation and highlights the importance of incorporating more realistic phenological mechanisms in models that seek to improve the projection of future carbon dynamics in tropical evergreen forests.

A roadmap for improving the representation of photosynthesis in Earth system models

Accurate representation of photosynthesis in terrestrial biosphere models (TBMs) is essential for robust projections of global change. However, current representations vary markedly between TBMs, contributing uncertainty to projections of global carbon fluxes. Here we compared the representation of photosynthesis in seven TBMs by examining leaf and canopy level responses of photosynthetic CO₂ assimilation (A) to key environmental variables: light, temperature, CO₂ concentration, vapor pressure deficit and soil water content. We identified research areas where limited process knowledge prevents inclusion of physiological phenomena in current TBMs and research areas where data are urgently needed for model parameterization or evaluation. We provide a roadmap for new science needed to improve the representation of photosynthesis in the next generation of terrestrial biosphere and Earth system models.