How vertical patterns in leaf traits shift seasonally and the implications for modeling canopy photosynthesis in a temperate deciduous forest

Dynamic Simulation of the Leaf Mass per Area (LMA) in Multilayer Crowns of Young Larix principis-rupprechtii

Leaf mass per area (LMA) is a key structural parameter that reflects the functional traits of leaves and plays a vital role in simulating the material and energy cycles of plant ecosystems. In this study, vertical whorl-by-whorl sampling of LMA was conducted in a young Larix principis-rupprechtii plantation during the growing season at the Saihanba Forest Farm. The vertical and seasonal variations in LMA were analysed. Subsequently, a predictive model of LMA was constructed. The results revealed that the LMA varied significantly between different crown whorls and growing periods. In the vertical direction of the crown, the LMA decreased with increasing crown depth, but the range of LMA values from the tree top to the bottom was, on average, 30.4 g/m2, which was approximately 2.5 times greater in the fully expanded phase than in the early leaf-expanding phase. During different growing periods, the LMA exhibited an allometric growth trend that increased during the leaf-expanding phase and then tended to stabilize. However, the range of LMA values throughout the growing period was, on average, 40.4 g/m2. Among the univariate models, the leaf dry matter content (LDMC) performed well (adjusted determination coefficient (Ra2) = 0.45, root mean square error (RMSE) = 13.48 g/m2) in estimating the LMA. The correlation between LMA and LDMC significantly differed at different growth stages and at different vertical crown whorls. The dynamic predictive model of LMA constructed with the relative depth in the crown (RDINC) and date of the year (DOY) as independent variables was reliable in both the assessments (Ra2 = 0.68, RMSE = 10.25 g/m2) and the validation (absolute mean error (MAE) = 8.05 g/m2, fit index (FI) = 0.682). Dynamic simulations of crown LMA provide a basis for elucidating the mechanism of crown development and laying the foundation for the construction of an ecological process model.

Intraspecific plasticity and co-variation of leaf traits facilitate Ficus tinctoria to acclimate hemiepiphytic and terrestrial habitats

Despite intensive studies on plant functional traits, the intraspecific variation and their co-variation at the multi-scale remains poorly studied, which holds the potential to unveil plant responses to changing environmental conditions. In this study, intraspecific variations of 16 leaf functional traits of a common fig species, Ficus tinctoria G. Frost., were investigated in relation to different scales: habitat types (hemiepiphytic and terrestrial), growth stages (small, medium and large) and tree crown positions (upper, middle and lower) in Xishuangbanna, Southwest China. Remarkable intraspecific variation was observed in leaf functional traits, which was mainly influenced by tree crown position, growth stage and their interaction. Stable nitrogen isotope (δ15N) and leaf area (LA) showed large variations, while stable carbon isotope (δ13C), stomata width and leaf water content showed relatively small variations, suggesting that light- and nitrogen-use strategies of F. tinctoria were plastic, while the water-use strategies have relatively low plasticity. The crown layers are formed with the growth of figs, and leaves in the lower crown increase their chlorophyll concentration and LA to improve the light energy conversion efficiency and the ability to capture weak light. Meanwhile, leaves in the upper crown increase the water-use efficiency to maintain their carbon assimilation. Moreover, hemiepiphytic medium (transitional stage) and large (free-standing stage) figs exhibited more significant trait differentiation (chlorophyll concentration, δ13C, stomata density, etc.) within the crown positions, and stronger trait co-variation compared with their terrestrial counterparts. This pattern demonstrates their acclimation to the changing microhabitats formed by their hemiepiphytic life history. Our study emphasizes the importance of multi-scaled intraspecific variation and co-variation in trait-based strategies of hemiepiphyte and terrestrial F. tinctoria, which facilitate them to cope with different environmental conditions.

Variation in leaf photosynthetic capacity within plant canopies: optimization, structural, and physiological constraints and inefficiencies

Leaf photosynthetic capacity (light-saturated net assimilation rate, AA) increases from bottom to top of plant canopies as the most prominent acclimation response to the conspicuous within-canopy gradients in light availability. Light-dependent variation in AA through plant canopies is associated with changes in key leaf structural (leaf dry mass per unit leaf area), chemical (nitrogen (N) content per area and dry mass, N partitioning between components of photosynthetic machinery), and physiological (stomatal and mesophyll conductance) traits, whereas the contribution of different traits to within-canopy AA gradients varies across sites, species, and plant functional types. Optimality models maximizing canopy carbon gain for a given total canopy N content predict that AA should be proportionally related to canopy light availability. However, comparison of model expectations with experimental data of within-canopy photosynthetic trait variations in representative plant functional types indicates that such proportionality is not observed in real canopies, and AA vs. canopy light relationships are curvilinear. The factors responsible for deviations from full optimality include stronger stomatal and mesophyll diffusion limitations at higher light, reflecting greater water limitations and more robust foliage in higher light. In addition, limits on efficient packing of photosynthetic machinery within leaf structural scaffolding, high costs of N redistribution among leaves, and limited plasticity of N partitioning among components of photosynthesis machinery constrain AA plasticity. Overall, this review highlights that the variation of AA through plant canopies reflects a complex interplay between adjustments of leaf structure and function to multiple environmental drivers, and that AA plasticity is limited by inherent constraints on and trade-offs between structural, chemical, and physiological traits. I conclude that models trying to simulate photosynthesis gradients in plant canopies should consider co-variations among environmental drivers, and the limitation of functional trait variation by physical constraints and include the key trade-offs between structural, chemical, and physiological leaf characteristics.

Patterns of physical, chemical, and metabolic characteristics of sugar maple leaves with depth in the crown and in response to nitrogen and phosphorus addition

Few previous studies have described patterns of leaf characteristics in response to nutrient availability and depth in the crown. Sugar maple has been studied for both sensitivity to light, as a shade-tolerant species, and sensitivity to soil nutrient availability, as a species in decline due to acid rain. To explore leaf characteristics from the top to bottom of the canopy, we collected leaves along a vertical gradient within mature sugar maple crowns in a full-factorial nitrogen by phosphorus addition experiment in three forest stands in central New Hampshire, USA. Thirty-two of the 44 leaf characteristics had significant relationships with depth in the crown, with the effect of depth in the crown strongest for leaf area, photosynthetic pigments, and polyamines. Nitrogen addition had a strong impact on the concentration of foliar N, chlorophyll, carotenoids, alanine, and glutamate. For several other elements and amino acids, N addition changed patterns with depth in the crown. Phosphorus addition increased foliar P and B; it also caused a steeper increase of P and B with depth in the crown. Since most of these leaf characteristics play a direct or indirect role in photosynthesis, metabolic regulation, or cell division, studies that ignore the vertical gradient may not accurately represent whole-canopy performance.

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

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

Long-term leaf C:N ratio change under elevated CO2 and nitrogen deposition in China: Evidence from observations and process-based modeling

Climate change, elevating atmosphere CO₂ (eCO₂) and increased nitrogen deposition (iNDEP) are altering the biogeochemical interactions between plants, microbes and soils, which further modify plant leaf carbon‑nitrogen (C:N) stoichiometry and their carbon assimilation capability. Many field experiments have observed large sensitivity of leaf C:N ratio to eCO₂ and iNDEP. However, the large-scale pattern of this sensitivity is still unclear, because eCO₂ and iNDEP drive leaf C:N ratio toward opposite directions, which are further compounded by the complex processes of nitrogen acquisition and plant-and-microbial nitrogen competition. Here, we attempt to map the leaf C:N ratio spatial variation in the past 5 decades in China with a combination of data-driven model and process-based modeling. These two approaches showed consistent results. Over different regions, we found that leaf C:N ratio had significant but uneven changes between 2 time periods (1960-1989 and 1990-2015): a 5% ± 8% increase for temperate grasslands in northern China, a 3% ± 6% increase for boreal grasslands in western China, and by contrast, a 7% ± 6% decrease for temperate forests in southern China, and a 3% ± 5% decrease for boreal forests in northeastern China. Additionally, the structural equation models indicated that the leaf C:N change was sensitive to ΔNDEP, ΔCO₂ and ΔMAT rather than ΔMAP and ecosystem types. Process-based modeling suggested that iNDEP was the main source of soil mineral nitrogen change, dominating leaf C:N ratio change in most areas in China, while eCO₂ led to leaf C:N ratio increase in low iNDEP area. This study also indicates that the long-term leaf C:N ratio acclimation was dominated by climate constraint, especially temperature, but was constrained by soil N availability over decade scale.

Determination of the most effective design for the measurement of photosynthetic light-response curves for planted Larix olgensis trees

A photosynthetic light-response (PLR) curve is a mathematical description of a single biochemical process and has been widely applied in many eco-physiological models. To date, many PLR measurement designs have been suggested, although their differences have rarely been explored, and the most effective design has not been determined. In this study, we measured three types of PLR curves (High, Middle and Low) from planted Larix olgensis trees by setting 31 photosynthetically active radiation (PAR) gradients. More than 530 million designs with different combinations of PAR gradients from 5 to 30 measured points were conducted to fit each of the three types of PLR curves. The influence of different PLR measurement designs on the goodness of fit of the PLR curves and the accuracy of the estimated photosynthetic indicators were analysed, and the optimal design was determined. The results showed that the measurement designs with fewer PAR gradients generally resulted in worse predicted accuracy for the photosynthetic indicators. However, the accuracy increased and remained stable when more than ten measurement points were used for the PAR gradients. The mean percent error (M%E) of the estimated maximum net photosynthetic rate (Pmax) and dark respiratory rate (Rd) for the designs with less than ten measurement points were, on average, 16.4 times and 20.1 times greater than those for the designs with more than ten measurement points. For a single tree, a unique PLR curve design generally reduced the accuracy of the predicted photosynthetic indicators. Thus, three optimal measurement designs were provided for the three PLR curve types, in which the root mean square error (RMSE) values reduced by an average of 8.3% and the coefficient of determination (R2) values increased by 0.3%. The optimal design for the High PLR curve type should shift more towards high-intensity PAR values, which is in contrast to the optimal design for the Low PLR curve type, which should shift more towards low-intensity PAR values.

Dynamic Simulation of the Crown Net Photosynthetic Rate for Young Larix olgensis Henry Trees

Numerical integration of the instantaneous net photosynthetic rate (An) is a common method for calculating the long-term CO2 uptake of trees, and accurate dynamic simulation of the crown An has been receiving substantial attention. Tree characteristics are challenging to assess given their aerodynamically coarse crown properties, spatiotemporal variation in leaf functional traits and microenvironments. Therefore, the variables associated with the dynamic variations in the crown An must be identified. The relationships of leaf temperature (Tleaf), the vapor pressure deficit (VPD), leaf mass per area (LMA) and the relative depth into the crown (RDINC) with the parameters of the photosynthetic light-response (PLR) model of Larix olgensis Henry were analyzed. The LMA, RDINC and VPD were highly correlated with the maximum net photosynthetic rate (Amax). The VPD was the key variable that mainly determined the variation in the apparent quantum yield (AQY). Tleaf exhibited a significant exponential correlation with the dark respiration rate (Rd). According to the above correlations, the crown PLR model of L. olgensis trees was constructed by linking VPD, LMA and RDINC to the original PLR equation. The model performed well, with a high coefficient of determination (R2) value (0.883) and low root mean square error (RMSE) value (1.440 μmol m−2 s−1). The extinction coefficient (k) of different pseudowhorls within a crown was calculated by the Beer–Lambert equation based on the observed photosynthetically active radiation (PAR) distribution. The results showed that k was not a constant value but varied with the RDINC, solar elevation angle (ψ) and cumulative leaf area of the whole crown (CLA). Thus, we constructed a k model by reparameterizing the power function of RDINC with the ψ and CLA, and the PAR distribution within a crown was therefore well estimated (R2 = 0.698 and RMSE = 174.4 μmol m−2 s−1). Dynamic simulation of the crown An for L. olgensis trees was achieved by combining the crown PLR model and dynamic PAR distribution model. Although the models showed some weakened physiological biochemical processes during photosynthesis, they enabled the estimation of long-term CO2 uptake for an L. olgensis plantation, and the results could be easily fitted to gas-exchange measurements.

Vertical patterns of photosynthesis and related leaf traits in two contrasting agricultural crops

To include within-canopy leaf acclimation responses to light and other resource gradients in photosynthesis modelling, it is imperative to understand the variation of leaf structural, biochemical and physiological traits from canopy top to bottom. In the present study, leaf photosynthetic traits for top and bottom canopy leaves, canopy structure and light profiles, were measured over one growing season for two contrasting crop types, winter barley (Hordeum vulgare L.) and rape seed (Brassica napus L.). With the exception of quantum yield, other traits such as maximum photosynthetic capacity (Amax), dark respiration, leaf nitrogen and chlorophyll contents, and leaf mass per area, showed consistently higher (P<0.05) values for top leaves throughout the growing season and for both crop types. Even though Amax was higher for top leaves, the bottom half of the canopy intercepted more light and thus contributed the most to total canopy photosynthesis up until senescence set in. Incorporating this knowledge into a simple top/bottom-leaf upscaling scheme, separating top and bottom leaves, resulted in a better match between estimated and measured total canopy photosynthesis, compared with a one-leaf upscaling scheme. Moreover, aggregating to daily and weekly temporal resolutions progressively increased the linearity of the leaf photosynthetic responses to light for top leaves.

Spatial and Seasonal Variations of Standardized Photosynthetic Parameters under Different Environmental Conditions for Young Planted Larix olgensis Henry Trees

Understanding the spatial and seasonal variations in leaf physiology is critical for accurately modeling the carbon uptake, physiological processes and growth of entire canopies and stands. For a 17-year-old Larix olgensis Henry plantation, vertical whorl-by-whorl sampling and analyses of seasonally repeated measurements of major photosynthetic parameters were conducted, and the correlations between photosynthetic parameters and environmental conditions, leaf morphological traits and spatial position within the crown were analyzed. According to the correlations, the photosynthetic parameters were standardized based on the environmental conditions to avoid the influence of the changing environment on the patterns of spatial and seasonal variations of photosynthetic parameters. The results showed that the standardized light-saturated net photosynthetic rate (SPmax), standardized dark respiration (SRd) and standardized stomatal conductance under saturated light (Sgs-sat) were all negatively related to the relative depth into the crown (RDINC) throughout the growing season. However, their vertical patterns were different during the development of the phenological phase. In addition, different gradients of environmental conditions also influenced the values and the range of the vertical variation in photosynthesis. High temperature and low humidity usually resulted in smaller values and weaker vertical variations of SPmax and Sgs-sat, but larger values and more obvious vertical variations in SRd. SPmax and Sgs-sat usually exhibited a parabolic seasonal pattern in different vertical positions within the crown; however, SRd generally followed a concave pattern. These seasonal patterns were all weaker with increasing RDINC. Different environments also exhibited a significant influence on the seasonal patterns of photosynthesis. We suggested that standardization is necessary before analyzing spatial and seasonal variations. A single environmental condition could not represent the spatial and seasonal patterns under all gradients of the environment. Spatial and seasonal variations should be simultaneously analyzed because they are related to each other.

Effects of Forest Canopy Vertical Stratification on the Estimation of Gross Primary Production by Remote Sensing

Gross primary production (GPP) in forests is the most important carbon flux in terrestrial ecosystems. Forest ecosystems with high leaf area index (LAI) values have diverse species or complex forest structures with vertical stratifications that influence the carbon–water–energy cycles. In this study, we used three light use efficiency (LUE) GPP models and site-level experiment data to analyze the effects of the vertical stratification of dense forest vegetation on the estimates of remotely sensed GPP during the growing season of two forest sites in East Asia: Dinghushan (DHS) and Tomakomai (TMK). The results showed that different controlling environmental factors of the vertical layers, such as temperature and vapor pressure deficit (VPD), produce different responses for the same LUE value in the different sub-ecosystems (defined as the tree, shrub, and grass layers), which influences the GPP estimation. Air temperature and VPD play important roles in the effects of vertical stratification on the GPP estimates in dense forests, which led to differences in GPP uncertainties from −50% to 30% because of the distinct temperature responses in TMK. The unequal vertical LAI distributions in the different sub-ecosystems led to GPP variations of 1–2 gC/m2/day with uncertainties of approximately −30% to 20% because sub-ecosystems have unique absorbed fractions of photosynthetically active radiation (APAR) and LUE. A comparison with the flux tower-based GPP data indicated that the GPP estimations from the LUE and APAR values from separate vertical layers exhibited better model performance than those calculated using the single-layer method, with 10% less bias in DHS and more than 70% less bias in TMK. The precision of the estimated GPP in regions with thick understory vegetation could be effectively improved by considering the vertical variations in environmental parameters and the LAI values of different sub-ecosystems as separate factors when calculating the GPP of different components. Our results provide useful insight that can be used to improve the accuracy of remote sensing GPP estimations by considering vertical stratification parameters along with the LAI of sub-ecosystems in dense forests.

Vertical leaf mass per area gradient of mature sugar maple reflects both height-driven increases in vascular tissue and light-driven increases in palisade layer thickness

A key trait used in canopy and ecosystem function modeling, leaf mass per area (LMA), is influenced by changes in both leaf thickness and leaf density (LMA = Thickness × Density). In tall trees, LMA is understood to increase with height through two primary mechanisms: (i) increasing palisade layer thickness (and thus leaf thickness) in response to light and/or (ii) reduced cell expansion and intercellular air space in response to hydrostatic constraints, leading to increased leaf density. Our objective was to investigate within-canopy gradients in leaf anatomical traits in order to understand environmental factors that influence leaf morphology in a sugar maple (Acer saccharum Marshall) forest canopy. We teased apart the effects of light and height on anatomical traits by sampling at exposed and closed canopies that had different light conditions at similar heights. As expected, palisade layer thickness responded strongly to cumulative light exposure. Mesophyll porosity, however, was weakly and negatively correlated with light and height (i.e., hydrostatic gradients). Reduced mesophyll porosity was not likely caused by limitations on cell expansion; in fact, epidermal cell width increased with height. Palisade layer thickness was better related to LMA, leaf density and leaf thickness than was mesophyll porosity. Vein diameter and fraction of vascular tissue also increased with height and LMA, density and thickness, revealing that greater investment in vascular and support tissue may be a third mechanism for increased LMA with height. Overall, decreasing mesophyll porosity with height was likely due to palisade cells expanding into the available air space and also greater investments in vascular and support tissue, rather than a reduction of cell expansion due to hydrostatic constraints. Our results provide evidence that light influences both palisade layer thickness and mesophyll porosity and indicate that hydrostatic gradients influence leaf vascular and support tissues in mature Acer saccharum trees.

Canopy gradients in leaf functional traits for species that differ in growth strategies and shade tolerance

In temperate deciduous forests, vertical gradients in leaf mass per area (LMA) and area-based leaf nitrogen (Narea) are strongly controlled by gradients in light availability. While there is evidence that hydrostatic constraints on leaf development may diminish LMA and Narea responses to light, inherent differences among tree species may also influence leaf developmental and morphological response to light. We investigated vertical gradients in LMA, Narea and leaf carbon isotope composition (δ13C) for three temperate deciduous species (Carpinus caroliniana Walter, Fagus grandifolia Ehrh., Liriodendron tulipifera L.) that differed in growth strategy (e.g., indeterminate and determinate growth), shade tolerance and leaf area to sapwood ratio (Al:As). Leaves were sampled across a broad range of light conditions within three vertical layers of tree crowns to maximize variation in light availability at each height and to minimize collinearity between light and height. All species displayed similar responses to light with respect to Narea and δ13C, but not for LMA. Light was more important for gradients in LMA for the shade-tolerant (C. caroliniana) and -intolerant (L. tulipifera) species with indeterminate growth, and height (e.g., hydrostatic gradients) and light were equally important for the shade-tolerant (F. grandifolia) species with determinate growth. Fagus grandifolia had a higher morphological plasticity in response to light, which may offer a competitive advantage in occupying a broader range of light conditions throughout the canopy. Differences in responses to light and height for the taller tree species, L. tulipifera and F. grandifolia, may be attributed to differences in growth strategy or Al:As, which may alter morphological and functional responses to light availability. While height was important in F. grandifolia, height was no more robust in predicting LMA than light in any of the species, confirming the strong role of light availability in determining LMA for temperate deciduous species.