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.

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA.

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond to (resistance) and recover from (resilience) rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the Ecosystem Demography-2.2 (ED-2.2) model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900 (approximating the clear-cut conditions at the time), and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes such as gross primary production (GPP), and interacted with one another. We then calculated resistance and resilience, two dimensions of ecosystem stability. Modeled GPP exhibited a two-fold decrease in mean resistance across disturbance severities of 45%, 65%, and 85% mortality; conversely, resilience increased by a factor of two with increasing disturbance severity. This pattern held for net primary production and net ecosystem production, indicating a trade-off in which greater initial declines were followed by faster recovery. Notably, however, heterotrophic respiration responded more slowly to disturbance, and it's highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

Historically inconsistent productivity and respiration fluxes in the global terrestrial carbon cycle.

The terrestrial carbon cycle is a major source of uncertainty in climate projections. Its dominant fluxes, gross primary productivity (GPP), and respiration (in particular soil respiration, RS), are typically estimated from independent satellite-driven models and upscaled in situ measurements, respectively. We combine carbon-cycle flux estimates and partitioning coefficients to show that historical estimates of global GPP and RS are irreconcilable. When we estimate GPP based on RS measurements and some assumptions about RS:GPP ratios, we found the resulted global GPP values (bootstrap mean [Formula: see text] Pg C yr-1) are significantly higher than most GPP estimates reported in the literature ([Formula: see text] Pg C yr-1). Similarly, historical GPP estimates imply a soil respiration flux (RsGPP, bootstrap mean of [Formula: see text] Pg C yr-1) statistically inconsistent with most published RS values ([Formula: see text] Pg C yr-1), although recent, higher, GPP estimates are narrowing this gap. Furthermore, global RS:GPP ratios are inconsistent with spatial averages of this ratio calculated from individual sites as well as CMIP6 model results. This discrepancy has implications for our understanding of carbon turnover times and the terrestrial sensitivity to climate change. Future efforts should reconcile the discrepancies associated with calculations for GPP and Rs to improve estimates of the global carbon budget.

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.

Structure and parameter uncertainty in centennial projections of forest community structure and carbon cycling.

Secondary forest regrowth shapes community succession and biogeochemistry for decades, including in the Upper Great Lakes region. Vegetation models encapsulate our understanding of forest function, and whether models can reproduce multi-decadal succession patterns is an indication of our ability to predict forest responses to future change. We test the ability of a vegetation model to simulate C cycling and community composition during 100 years of forest regrowth following stand-replacing disturbance, asking (a) Which processes and parameters are most important to accurately model Upper Midwest forest succession? (b) What is the relative importance of model structure versus parameter values to these predictions? We ran ensembles of the Ecosystem Demography model v2.2 with different representations of processes important to competition for light. We compared the magnitude of structural and parameter uncertainty and assessed which sub-model-parameter combinations best reproduced observed C fluxes and community composition. On average, our simulations underestimated observed net primary productivity (NPP) and leaf area index (LAI) after 100 years and predicted complete dominance by a single plant functional type (PFT). Out of 4,000 simulations, only nine fell within the observed range of both NPP and LAI, but these predicted unrealistically complete dominance by either early hardwood or pine PFTs. A different set of seven simulations were ecologically plausible but under-predicted observed NPP and LAI. Parameter uncertainty was large; NPP and LAI ranged from ~0% to >200% of their mean value, and any PFT could become dominant. The two parameters that contributed most to uncertainty in predicted NPP were plant-soil water conductance and growth respiration, both unobservable empirical coefficients. We conclude that (a) parameter uncertainty is more important than structural uncertainty, at least for ED-2.2 in Upper Midwest forests and (b) simulating both productivity and plant community composition accurately without physically unrealistic parameters remains challenging for demographic vegetation models.

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 [Rd ], maximum RuBisCo carboxylation rate [Vc,max ], maximum electron transport rate [Jmax ]) 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.

Calcium and aluminum cycling in a temperate broadleaved deciduous forest of the eastern USA: relative impacts of tree species, canopy state, and flux type.

Ca/Al molar ratios are commonly used to assess the extent of aluminum stress in forests. This is among the first studies to quantify Ca/Al molar ratios for stemflow. Ca/Al molar ratios in bulk precipitation, throughfall, stemflow, litter leachate, near-trunk soil solution, and soil water were quantified for a deciduous forest in northeastern MD, USA. Data were collected over a 3-year period. The Ca/Al molar ratios in this study were above the threshold for aluminum stress (<1). Fagus grandifolia Ehrh. (American beech) had a median annual stemflow Ca/Al molar ratio of 15.7, with the leafed and leafless values of 12.4 and 19.2, respectively. The corresponding Ca/Al molar ratios for Liriodendron tulipifera L. (yellow poplar) were 11.9 at the annual time scale and 11.9 and 13.6 for leafed and leafless periods, respectively. Bayesian statistical analysis showed no significant effect of canopy state (leafed, leafless) on Ca/Al molar ratios. DOC was consistently an important predictor of calcium, aluminum, and Ca/Al ratios. pH was occasionally an important predictor of calcium and aluminum concentrations, but was not a good predictor of Ca/Al ratio in any of the best-fit models (of >500 examined). This study supplies new data on Ca/Al molar ratios for stemflow from two common deciduous tree species. Future work should examine Ca/Al molar ratios in stemflow of other species and examine both inorganic and organic aluminum species to better gauge the potential for, and understand the dynamics of, aluminum toxicity in the proximal area around tree boles.

Stemflow acid neutralization capacity in a broadleaved deciduous forest: the role of edge effects.

Atmospheric deposition is an important pathway for moisture, nutrient, and pollutant exchange among the atmosphere, forest, and soils. Previous work has shown the importance of proximity to the forest edge to chemical fluxes in throughfall, but far less research has considered stemflow. This study examined the difference in acid neutralization capacity (ANC) of stemflow of nineteen Liriodendron tulipifera L. (yellow poplar) trees between the forest edge and interior in a rural area of northeastern Maryland. We measured ANC directly via potentiometric titration. Stemflow from trees at the forest edge was found to have significantly higher and more variable pH and ANC than in the forest interior (p < 0.01). No mathematical trend between ANC and distance to the forest edge was observed, indicating the importance of individual tree characteristics in stemflow production and chemistry. These results reaffirm the importance of stemflow for acid neutralization by deciduous tree species.