Vegetation demographics in Earth System Models: A review of progress and priorities

Preserving isohydricity: vertical environmental variability explains Amazon forest water-use strategies

Increases in hydrological extremes, including drought, are expected for Amazon forests. A fundamental challenge for predicting forest responses lies in identifying ecological strategies which underlie such responses. Characterization of species-specific hydraulic strategies for regulating water-use, thought to be arrayed along an 'isohydric-anisohydric' spectrum, is a widely used approach. However, recent studies have questioned the usefulness of this classification scheme, because its metrics are strongly influenced by environments, and hence can lead to divergent classifications even within the same species. Here, we propose an alternative approach positing that individual hydraulic regulation strategies emerge from the interaction of environments with traits. Specifically, we hypothesize that the vertical forest profile represents a key gradient in drought-related environments (atmospheric vapor pressure deficit, soil water availability) that drives divergent tree water-use strategies for coordinated regulation of stomatal conductance (gs) and leaf water potentials (ΨL) with tree rooting depth, a proxy for water availability. Testing this hypothesis in a seasonal eastern Amazon forest in Brazil, we found that hydraulic strategies indeed depend on height-associated environments. Upper canopy trees, experiencing high vapor pressure deficit (VPD), but stable soil water access through deep rooting, exhibited isohydric strategies, defined by little seasonal change in the diurnal pattern of gs and steady seasonal minimum ΨL. In contrast, understory trees, exposed to less variable VPD but highly variable soil water availability, exhibited anisohydric strategies, with fluctuations in diurnal gs that increased in the dry season along with increasing variation in ΨL. Our finding that canopy height structures the coordination between drought-related environmental stressors and hydraulic traits provides a basis for preserving the applicability of the isohydric-to-anisohydric spectrum, which we show here may consistently emerge from environmental context. Our work highlights the importance of understanding how environmental heterogeneity structures forest responses to climate change, providing a mechanistic basis for improving models of tropical ecosystems.

Remotely sensed crown nutrient concentrations modulate forest reproduction across the contiguous United States

Global forests are increasingly lost to climate change, disturbance, and human management. Evaluating forests' capacities to regenerate and colonize new habitats has to start with the seed production of individual trees and how it depends on nutrient access. Studies on the linkage between reproduction and foliar nutrients are limited to a few locations and few species, due to the large investment needed for field measurements on both variables. We synthesized tree fecundity estimates from the Masting Inference and Forecasting (MASTIF) network with foliar nutrient concentrations from hyperspectral remote sensing at the National Ecological Observatory Network (NEON) across the contiguous United States. We evaluated the relationships between seed production and foliar nutrients for 56,544 tree-years from 26 species at individual and community scales. We found a prevalent association between high foliar phosphorous (P) concentration and low individual seed production (ISP) across the continent. Within-species coefficients to nitrogen (N), potassium (K), calcium (Ca), and magnesium (Mg) are related to species differences in nutrient demand, with distinct biogeographic patterns. Community seed production (CSP) decreased four orders of magnitude from the lowest to the highest foliar P. This first continental-scale study sheds light on the relationship between seed production and foliar nutrients, highlighting the potential of using combined Light Detection And Ranging (LiDAR) and hyperspectral remote sensing to evaluate forest regeneration. The fact that both ISP and CSP decline in the presence of high foliar P levels has immediate application in improving forest demographic and regeneration models by providing more realistic nutrient effects at multiple scales.

Community assembly influences plant trait economic spectra and functional trade-offs at ecosystem scales

Plant functional traits hold the potential to greatly improve the understanding and prediction of climate impacts on ecosystems and carbon cycle feedback to climate change. Traits are commonly used to place species along a global conservative-acquisitive trade-off, yet how and if functional traits and conservative-acquisitive trade-offs scale up to mediate community and ecosystem fluxes is largely unknown. Here, we combine functional trait datasets and multibiome datasets of forest water and carbon fluxes at the species, community, and ecosystem-levels to quantify the scaling of the tradeoff between maximum flux and sensitivity to vapor pressure deficit. We find a strong conservative-acquisitive trade-off at the species scale, which weakens modestly at the community scale and largely disappears at the ecosystem scale. Functional traits, particularly plant water transport (hydraulic) traits, are strongly associated with the key dimensions of the conservative-acquisitive trade-off at community and ecosystem scales, highlighting that trait composition appears to influence community and ecosystem flux dynamics. Our findings provide a foundation for improving carbon cycle models by revealing i) that plant hydraulic traits are most strongly associated with community- and ecosystem scale flux dynamics and ii) community assembly dynamics likely need to be considered explicitly, as they give rise to ecosystem-level flux dynamics that differ substantially from trade-offs identified at the species-level.

Bomb radiocarbon evidence for strong global carbon uptake and turnover in terrestrial vegetation

Vegetation and soils are taking up approximately 30% of anthropogenic carbon dioxide emissions because of small imbalances in large gross carbon exchanges from productivity and turnover that are poorly constrained. We combined a new budget of radiocarbon produced by nuclear bomb testing in the 1960s with model simulations to evaluate carbon cycling in terrestrial vegetation. We found that most state-of-the-art vegetation models used in the Coupled Model Intercomparison Project underestimated the radiocarbon accumulation in vegetation biomass. Our findings, combined with constraints on vegetation carbon stocks and productivity trends, imply that net primary productivity is likely at least 80 petagrams of carbon per year presently, compared with the 43 to 76 petagrams per year predicted by current models. Storage of anthropogenic carbon in terrestrial vegetation is likely more short-lived and vulnerable than previously predicted.

Positive feedbacks and alternative stable states in forest leaf types

The emergence of alternative stable states in forest systems has significant implications for the functioning and structure of the terrestrial biosphere, yet empirical evidence remains scarce. Here, we combine global forest biodiversity observations and simulations to test for alternative stable states in the presence of evergreen and deciduous forest types. We reveal a bimodal distribution of forest leaf types across temperate regions of the Northern Hemisphere that cannot be explained by the environment alone, suggesting signatures of alternative forest states. Moreover, we empirically demonstrate the existence of positive feedbacks in tree growth, recruitment and mortality, with trees having 4-43% higher growth rates, 14-17% higher survival rates and 4-7 times higher recruitment rates when they are surrounded by trees of their own leaf type. Simulations show that the observed positive feedbacks are necessary and sufficient to generate alternative forest states, which also lead to dependency on history (hysteresis) during ecosystem transition from evergreen to deciduous forests and vice versa. We identify hotspots of bistable forest types in evergreen-deciduous ecotones, which are likely driven by soil-related positive feedbacks. These findings are integral to predicting the distribution of forest biomes, and aid to our understanding of biodiversity, carbon turnover, and terrestrial climate feedbacks.

Developing a predictive science of the biosphere requires the integration of scientific cultures

Increasing the speed of scientific progress is urgently needed to address the many challenges associated with the biosphere in the Anthropocene. Consequently, the critical question becomes: How can science most rapidly progress to address large, complex global problems? We suggest that the lag in the development of a more predictive science of the biosphere is not only because the biosphere is so much more complex, or because we do not have enough data, or are not doing enough experiments, but, in large part, because of unresolved tension between the three dominant scientific cultures that pervade the research community. We introduce and explain the concept of the three scientific cultures and present a novel analysis of their characteristics, supported by examples and a formal mathematical definition/representation of what this means and implies. The three cultures operate, to varying degrees, across all of science. However, within the biosciences, and in contrast to some of the other sciences, they remain relatively more separated, and their lack of integration has hindered their potential power and insight. Our solution to accelerating a broader, predictive science of the biosphere is to enhance integration of scientific cultures. The process of integration-Scientific Transculturalism-recognizes that the push for interdisciplinary research, in general, is just not enough. Unless these cultures of science are formally appreciated and their thinking iteratively integrated into scientific discovery and advancement, there will continue to be numerous significant challenges that will increasingly limit forecasting and prediction efforts.

Trait-environment relationships are timescale dependent

This article is a Commentary on Famiglietti et al. (2024), 241: 2423–2434.

Geographic range of plants drives long-term climate change

Long computation times in vegetation and climate models hamper our ability to evaluate the potentially powerful role of plants on weathering and carbon sequestration over the Phanerozoic Eon. Simulated vegetation over deep time is often homogenous, and disregards the spatial distribution of plants and the impact of local climatic variables on plant function. Here we couple a fast vegetation model (FLORA) to a spatially-resolved long-term climate-biogeochemical model (SCION), to assess links between plant geographical range, the long-term carbon cycle and climate. Model results show lower rates of carbon fixation and up to double the previously predicted atmospheric CO2 concentration due to a limited plant geographical range over the arid Pangea supercontinent. The Mesozoic dispersion of the continents increases modelled plant geographical range from 65% to > 90%, amplifying global CO2 removal, consistent with geological data. We demonstrate that plant geographical range likely exerted a major, under-explored control on long-term climate change.

Climate regulation processes are linked to the functional composition of plant communities in European forests, shrublands, and grasslands

Terrestrial ecosystems affect climate by reflecting solar irradiation, evaporative cooling, and carbon sequestration. Yet very little is known about how plant traits affect climate regulation processes (CRPs) in different habitat types. Here, we used linear and random forest models to relate the community-weighted mean and variance values of 19 plant traits (summarized into eight trait axes) to the climate-adjusted proportion of reflected solar irradiation, evapotranspiration, and net primary productivity across 36,630 grid cells at the European extent, classified into 10 types of forest, shrubland, and grassland habitats. We found that these trait axes were more tightly linked to log evapotranspiration (with an average of 6.2% explained variation) and the proportion of reflected solar irradiation (6.1%) than to net primary productivity (4.9%). The highest variation in CRPs was explained in forest and temperate shrubland habitats. Yet, the strength and direction of these relationships were strongly habitat-dependent. We conclude that any spatial upscaling of the effects of plant communities on CRPs must consider the relative contribution of different habitat types.

Tree size distribution as the stationary limit of an evolutionary master equation

The diameter distribution of a given species of deciduous trees is well approximated by a Gamma distribution. Here we give new experimental evidence for this conjecture by analyzing deciduous tree size data in mature semi-natural forest and ancient, traditionally managed wood-pasture from Central Europe. These distribution functions collapse on a universal shape if the tree sizes are normalized to the mean value in the considered sample. A new evolutionary master equation is used to model the observed distribution. The model incorporates four ecological processes: tree growth, mortality, recruitment, and diversification. Utilizing simple and realistic kernel functions describing the first three, along with an assumed multiplicative dilution due to diversification, the stationary solution of the master equation yields the experimentally observed Gamma distribution. The model as it is formulated allows an analytically compact solution and has only two fitting parameters whose values are consistent with the experimental data related to these processes. We found that the equilibrium size distribution of tree species with different ecology, originating from two contrastingly different semi-natural ecosystem types can be accurately described by a single dynamical mean-field model.

Constraining long-term model predictions for woody growth using tropical tree rings

The strength and persistence of the tropical carbon sink hinges on the long-term responses of woody growth to climatic variations and increasing CO2 . However, the sensitivity of tropical woody growth to these environmental changes is poorly understood, leading to large uncertainties in growth predictions. Here, we used tree ring records from a Southeast Asian tropical forest to constrain ED2.2-hydro, a terrestrial biosphere model with explicit vegetation demography. Specifically, we assessed individual-level woody growth responses to historical climate variability and increases in atmospheric CO2 (Ca ). When forced with historical Ca , ED2.2-hydro reproduced the magnitude of increases in intercellular CO2 concentration (a major determinant of photosynthesis) estimated from tree ring carbon isotope records. In contrast, simulated growth trends were considerably larger than those obtained from tree rings, suggesting that woody biomass production efficiency (WBPE = woody biomass production:gross primary productivity) was overestimated by the model. The estimated WBPE decline under increasing Ca based on model-data discrepancy was comparable to or stronger than (depending on tree species and size) the observed WBPE changes from a multi-year mature-forest CO2 fertilization experiment. In addition, we found that ED2.2-hydro generally overestimated climatic sensitivity of woody growth, especially for late-successional plant functional types. The model-data discrepancy in growth sensitivity to climate was likely caused by underestimating WBPE in hot and dry years due to commonly used model assumptions on carbon use efficiency and allocation. To our knowledge, this is the first study to constrain model predictions of individual tree-level growth sensitivity to Ca and climate against tropical tree-ring data. Our results suggest that improving model processes related to WBPE is crucial to obtain better predictions of tropical forest responses to droughts and increasing Ca . More accurate parameterization of WBPE will likely reduce the stimulation of woody growth by Ca rise predicted by biosphere models.

The anthropogenic imprint on temperate and boreal forest demography and carbon turnover

Aim

The sweeping transformation of the biosphere by humans over the last millennia leaves only limited windows into its natural state. Much of the forests that dominated temperate and southern boreal regions have been lost and those that remain typically bear a strong imprint of forestry activities and past land-use change, which have changed forest age structure and composition. Here, we ask how would the dynamics, structure and function of temperate and boreal forests differ in the absence of forestry and the legacies of land-use change?

Location

Global.

Time Period

2001-2014, integrating over the legacy of disturbance events from 1875 to 2014.

Major Taxa Studied

Trees.

Methods

We constructed an empirical model of natural disturbance probability as a function of community traits and climate, based on observed disturbance rate and form across 77 protected forest landscapes distributed across three continents. Coupling this within a dynamic vegetation model simulating forest composition and structure, we generated estimates of stand-replacing disturbance return intervals in the absence of forestry for northern hemisphere temperate and boreal forests. We then applied this model to calculate forest stand age structure and carbon turnover rates.

Results

Comparison with observed disturbance rates revealed human activities to have almost halved the median return interval of stand-replacing disturbances across temperate forest, with more moderate changes in the boreal region. The resulting forests are typically much younger, especially in northern Europe and south-eastern North America, resulting in a 32% reduction in vegetation carbon turnover time across temperate forests and a 7% reduction for boreal forests.

Conclusions

The current northern hemisphere temperate forest age structure is dramatically out of equilibrium with its natural disturbance regimes. Shifts towards more nature-based approaches to forest policy and management should more explicitly consider the current disturbance surplus, as it substantially impacts carbon dynamics and litter (including deadwood) stocks.

Accounting for herbaceous communities in process-based models will advance our understanding of "grassy" ecosystems

Grassland and other herbaceous communities cover significant portions of Earth's terrestrial surface and provide many critical services, such as carbon sequestration, wildlife habitat, and food production. Forecasts of global change impacts on these services will require predictive tools, such as process-based dynamic vegetation models. Yet, model representation of herbaceous communities and ecosystems lags substantially behind that of tree communities and forests. The limited representation of herbaceous communities within models arises from two important knowledge gaps: first, our empirical understanding of the principles governing herbaceous vegetation dynamics is either incomplete or does not provide mechanistic information necessary to drive herbaceous community processes with models; second, current model structure and parameterization of grass and other herbaceous plant functional types limits the ability of models to predict outcomes of competition and growth for herbaceous vegetation. In this review, we provide direction for addressing these gaps by: (1) presenting a brief history of how vegetation dynamics have been developed and incorporated into earth system models, (2) reporting on a model simulation activity to evaluate current model capability to represent herbaceous vegetation dynamics and ecosystem function, and (3) detailing several ecological properties and phenomena that should be a focus for both empiricists and modelers to improve representation of herbaceous vegetation in models. Together, empiricists and modelers can improve representation of herbaceous ecosystem processes within models. In so doing, we will greatly enhance our ability to forecast future states of the earth system, which is of high importance given the rapid rate of environmental change on our planet.

Wood-density has no effect on stomatal control of leaf-level water use efficiency in an Amazonian forest

Forest disturbances increase the proportion of fast-growing tree species compared to slow-growing ones. To understand their relative capacity for carbon uptake and their vulnerability to climate change, and to represent those differences in Earth system models, it is necessary to characterise the physiological differences in their leaf-level control of water use efficiency and carbon assimilation. We used wood density as a proxy for the fast-slow growth spectrum and tested the assumption that trees with a low wood density (LWD) have a lower water-use efficiency than trees with a high wood density (HWD). We selected 5 LWD tree species and 5 HWD tree species growing in the same location in an Amazonian tropical forest and measured in situ steady-state gas exchange on top-of-canopy leaves with parallel sampling and measurement of leaf mass area and leaf nitrogen content. We found that LWD species invested more nitrogen in photosynthetic capacity than HWD species, had higher photosynthetic rates and higher stomatal conductance. However, contrary to expectations, we showed that the stomatal control of the balance between transpiration and carbon assimilation was similar in LWD and HWD species and that they had the same dark respiration rates.

A new approach to derive productivity of tropical forests using radar remote sensing measurements

Deriving gross & net primary productivity (GPP & NPP) and carbon turnover time of forests from remote sensing remains challenging. This study presents a novel approach to estimate forest productivity by combining radar remote sensing measurements, machine learning and an individual-based forest model. In this study, we analyse the role of different spatial resolutions on predictions in the context of the Radar BIOMASS mission (by ESA). In our analysis, we use the forest gap model FORMIND in combination with a boosted regression tree (BRT) to explore how spatial biomass distributions can be used to predict GPP, NPP and carbon turnover time (τ) at different resolutions. We simulate different spatial biomass resolutions (4 ha, 1 ha and 0.04 ha) in combination with different vertical resolutions (20, 10 and 2 m). Additionally, we analysed the robustness of this approach and applied it to disturbed and mature forests. Disturbed forests have a strong influence on the predictions which leads to high correlations (R2 > 0.8) at the spatial scale of 4 ha and 1 ha. Increased vertical resolution leads generally to better predictions for productivity (GPP & NPP). Increasing spatial resolution leads to better predictions for mature forests and lower correlations for disturbed forests. Our results emphasize the value of the forthcoming BIOMASS satellite mission and highlight the potential of deriving estimates for forest productivity from information on forest structure. If applied to more and larger areas, the approach might ultimately contribute to a better understanding of forest ecosystems.

A theory of demographic optimality in forests

Carbon uptake by the land is a key determinant of future climate change. Unfortunately, Dynamic Global Vegetation Models have many unknown internal parameters which leads to significant uncertainty in projections of the future land carbon sink. By contrast, observed forest inventories in both Amazonia and the USA show strikingly common tree-size distributions, pointing to a simpler modelling paradigm. The curvature of these size-distributions is related to the ratio of mortality to growth in Demographic Equilibrium Theory (DET). We extend DET to include recruitment limited by competitive exclusion from existing trees. From this, we find simultaneous maxima of tree density and biomass in terms of respectively the ratio of mortality to growth and the proportion of primary productivity allocated to reproduction, an idea we call Demographic Optimality (DO). Combining DO with the ratio of mortality to growth common to the US and Amazon forests, results in the prediction that about an eighth of productivity should be allocated to reproduction, which is broadly consistent with observations. Another prediction of the model is that seed mortality should decrease with increasing seed size, such that the advantage of having many small seeds is nullified by the higher seed mortality. Demographic Optimality is therefore consistent with the common shape of tree-size distributions seen in very different forests, and an allocation to reproduction that is independent of seed size.

Accounting for trait variability and coordination in predictions of drought-induced range shifts in woody plants

Functional traits offer a promising avenue to improve predictions of species range shifts under climate change, which will entail warmer and often drier conditions. Although the conceptual foundation linking traits with plant performance and range shifts appears solid, the predictive ability of individual traits remains generally low. In this review, we address this apparent paradox, emphasizing examples of woody plants and traits associated with drought responses at the species' rear edge. Low predictive ability reflects the fact not only that range dynamics tend to be complex and multifactorial, as well as uncertainty in the identification of relevant traits and limited data availability, but also that trait effects are scale- and context-dependent. The latter results from the complex interactions among traits (e.g. compensatory effects) and between them and the environment (e.g. exposure), which ultimately determine persistence and colonization capacity. To confront this complexity, a more balanced coverage of the main functional dimensions involved (stress tolerance, resource use, regeneration and dispersal) is needed, and modelling approaches must be developed that explicitly account for: trait coordination in a hierarchical context; trait variability in space and time and its relationship with exposure; and the effect of biotic interactions in an ecological community context.

Models of underlying autotrophic biomass dynamics fit to daily river ecosystem productivity estimates improve understanding of ecosystem disturbance and resilience

Directly observing autotrophic biomass at ecologically relevant frequencies is difficult in many ecosystems, hampering our ability to predict productivity through time. Since disturbances can impart distinct reductions in river productivity through time by modifying underlying standing stocks of biomass, mechanistic models fit to productivity time series can infer underlying biomass dynamics. We incorporated biomass dynamics into a river ecosystem productivity model for six rivers to identify disturbance flow thresholds and understand the resilience of primary producers. The magnitude of flood necessary to disturb biomass and thereby reduce ecosystem productivity was consistently lower than the more commonly used disturbance flow threshold of the flood magnitude necessary to mobilize river bed sediment. The estimated daily maximum percent increase in biomass (a proxy for resilience) ranged from 5% to 42% across rivers. Our latent biomass model improves understanding of disturbance thresholds and recovery patterns of autotrophic biomass within river ecosystems.

Modelling the impact of forest management and CO2-fertilisation on growth and demography in a Sitka spruce plantation

Afforestation and reforestation to meet 'Net Zero' emissions targets are considered a necessary policy by many countries. Their potential benefits are usually assessed through forest carbon and growth models. The implementation of vegetation demography gives scope to represent forest management and other size-dependent processes within land surface models (LSMs). In this paper, we evaluate the impact of including management within an LSM that represents demography, using both in-situ and reanalysis climate drivers at a mature, upland Sitka spruce plantation in Northumberland, UK. We compare historical simulations with fixed and variable CO2 concentrations, and with and without tree thinning implemented. Simulations are evaluated against the observed vegetation structure and carbon fluxes. Including thinning and the impact of increasing CO2 concentration ('CO2 fertilisation') gave more realistic estimates of stand-structure and physical characteristics. Historical CO2 fertilisation had a noticeable effect on the Gross Primary Productivity seasonal-diurnal cycle and contributed to approximately 7% higher stand biomass by 2018. The net effect of both processes resulted in a decrease of tree density and biomass, but an increase in tree height and leaf area index.

Climate-ecosystem modelling made easy: The Land Sites Platform

Dynamic Global Vegetation Models (DGVMs) provide a state-of-the-art process-based approach to study the complex interplay between vegetation and its physical environment. For example, they help to predict how terrestrial plants interact with climate, soils, disturbance and competition for resources. We argue that there is untapped potential for the use of DGVMs in ecological and ecophysiological research. One fundamental barrier to realize this potential is that many researchers with relevant expertize (ecology, plant physiology, soil science, etc.) lack access to the technical resources or awareness of the research potential of DGVMs. Here we present the Land Sites Platform (LSP): new software that facilitates single-site simulations with the Functionally Assembled Terrestrial Ecosystem Simulator, an advanced DGVM coupled with the Community Land Model. The LSP includes a Graphical User Interface and an Application Programming Interface, which improve the user experience and lower the technical thresholds for installing these model architectures and setting up model experiments. The software is distributed via version-controlled containers; researchers and students can run simulations directly on their personal computers or servers, with relatively low hardware requirements, and on different operating systems. Version 1.0 of the LSP supports site-level simulations. We provide input data for 20 established geo-ecological observation sites in Norway and workflows to add generic sites from public global datasets. The LSP makes standard model experiments with default data easily achievable (e.g., for educational or introductory purposes) while retaining flexibility for more advanced scientific uses. We further provide tools to visualize the model input and output, including simple examples to relate predictions to local observations. The LSP improves access to land surface and DGVM modelling as a building block of community cyberinfrastructure that may inspire new avenues for mechanistic ecosystem research across disciplines.

Plant traits poorly predict winner and loser shrub species in a warming tundra biome

Climate change is leading to species redistributions. In the tundra biome, shrubs are generally expanding, but not all tundra shrub species will benefit from warming. Winner and loser species, and the characteristics that may determine success or failure, have not yet been fully identified. Here, we investigate whether past abundance changes, current range sizes and projected range shifts derived from species distribution models are related to plant trait values and intraspecific trait variation. We combined 17,921 trait records with observed past and modelled future distributions from 62 tundra shrub species across three continents. We found that species with greater variation in seed mass and specific leaf area had larger projected range shifts, and projected winner species had greater seed mass values. However, trait values and variation were not consistently related to current and projected ranges, nor to past abundance change. Overall, our findings indicate that abundance change and range shifts will not lead to directional modifications in shrub trait composition, since winner and loser species share relatively similar trait spaces.

Modeling the mechanisms of conifer mortality under seawater exposure

Relative sea level rise (SLR) increasingly impacts coastal ecosystems through the formation of ghost forests. To predict the future of coastal ecosystems under SLR and changing climate, it is important to understand the physiological mechanisms underlying coastal tree mortality and to integrate this knowledge into dynamic vegetation models. We incorporate the physiological effect of salinity and hypoxia in a dynamic vegetation model in the Earth system land model, and used the model to investigate the mechanisms of mortality of conifer forests on the west and east coast sites of USA, where trees experience different form of sea water exposure. Simulations suggest similar physiological mechanisms can result in different mortality patterns. At the east coast site that experienced severe increases in seawater exposure, trees loose photosynthetic capacity and roots rapidly, and both storage carbon and hydraulic conductance decrease significantly within a year. Over time, further consumption of storage carbon that leads to carbon starvation dominates mortality. At the west coast site that gradually exposed to seawater through SLR, hydraulic failure dominates mortality because root loss impacts on conductance are greater than the degree of storage carbon depletion. Measurements and modeling focused on understanding the physiological mechanisms of mortality is critical to reducing predictive uncertainty.

Spatial heterogeneity of global forest aboveground carbon stocks and fluxes constrained by spaceborne lidar data and mechanistic modeling

Forest carbon is a large and uncertain component of the global carbon cycle. An important source of complexity is the spatial heterogeneity of vegetation vertical structure and extent, which results from variations in climate, soils, and disturbances and influences both contemporary carbon stocks and fluxes. Recent advances in remote sensing and ecosystem modeling have the potential to significantly improve the characterization of vegetation structure and its resulting influence on carbon. Here, we used novel remote sensing observations of tree canopy height collected by two NASA spaceborne lidar missions, Global Ecosystem Dynamics Investigation and ICE, Cloud, and Land Elevation Satellite 2, together with a newly developed global Ecosystem Demography model (v3.0) to characterize the spatial heterogeneity of global forest structure and quantify the corresponding implications for forest carbon stocks and fluxes. Multiple-scale evaluations suggested favorable results relative to other estimates including field inventory, remote sensing-based products, and national statistics. However, this approach utilized several orders of magnitude more data (3.77 billion lidar samples) on vegetation structure than used previously and enabled a qualitative increase in the spatial resolution of model estimates achievable (0.25° to 0.01°). At this resolution, process-based models are now able to capture detailed spatial patterns of forest structure previously unattainable, including patterns of natural and anthropogenic disturbance and recovery. Through the novel integration of new remote sensing data and ecosystem modeling, this study bridges the gap between existing empirically based remote sensing approaches and process-based modeling approaches. This study more generally demonstrates the promising value of spaceborne lidar observations for advancing carbon modeling at a global scale.

Plant water use theory should incorporate hypotheses about extreme environments, population ecology, and community ecology

Plant water use theory has largely been developed within a plant-performance paradigm that conceptualizes water use in terms of value for carbon gain and that sits within a neoclassical economic framework. This theory works very well in many contexts but does not consider other values of water to plants that could impact their fitness. Here, we survey a range of alternative hypotheses for drivers of water use and stomatal regulation. These hypotheses are organized around relevance to extreme environments, population ecology, and community ecology. Most of these hypotheses are not yet empirically tested and some are controversial (e.g. requiring more agency and behavior than is commonly believed possible for plants). Some hypotheses, especially those focused around using water to avoid thermal stress, using water to promote reproduction instead of growth, and using water to hoard it, may be useful to incorporate into theory or to implement in Earth System Models.

The effect of the vertical gradients of photosynthetic parameters on the CO2 assimilation and transpiration of a Panamanian tropical forest

Terrestrial biosphere models (TBMs) include the representation of vertical gradients in leaf traits associated with modeling photosynthesis, respiration, and stomatal conductance. However, model assumptions associated with these gradients have not been tested in complex tropical forest canopies. We compared TBM representation of the vertical gradients of key leaf traits with measurements made in a tropical forest in Panama and then quantified the impact of the observed gradients on simulated canopy-scale CO₂ and water fluxes. Comparison between observed and TBM trait gradients showed divergence that impacted canopy-scale simulations of water vapor and CO₂ exchange. Notably, the ratio between the dark respiration rate and the maximum carboxylation rate was lower near the ground than at the top-of-canopy, leaf-level water-use efficiency was markedly higher at the top-of-canopy, and the decrease in maximum carboxylation rate from the top-of-canopy to the ground was less than TBM assumptions. The representation of the gradients of leaf traits in TBMs is typically derived from measurements made within-individual plants, or, for some traits, assumed constant due to a lack of experimental data. Our work shows that these assumptions are not representative of the trait gradients observed in species-rich, complex tropical forests.

Embracing fine-root system complexity in terrestrial ecosystem modeling

Projecting the dynamics and functioning of the biosphere requires a holistic consideration of whole-ecosystem processes. However, biases toward leaf, canopy, and soil modeling since the 1970s have constantly left fine-root systems being rudimentarily treated. As accelerated empirical advances in the last two decades establish clearly functional differentiation conferred by the hierarchical structure of fine-root orders and associations with mycorrhizal fungi, a need emerges to embrace this complexity to bridge the data-model gap in still extremely uncertain models. Here, we propose a three-pool structure comprising transport and absorptive fine roots with mycorrhizal fungi (TAM) to model vertically resolved fine-root systems across organizational and spatial-temporal scales. Emerging from a conceptual shift away from arbitrary homogenization, TAM builds upon theoretical and empirical foundations as an effective and efficient approximation that balances realism and simplicity. A proof-of-concept demonstration of TAM in a big-leaf model both conservatively and radically shows robust impacts of differentiation within fine-root systems on simulating carbon cycling in temperate forests. Theoretical and quantitative support warrants exploiting its rich potentials across ecosystems and models to confront uncertainties and challenges for a predictive understanding of the biosphere. Echoing a broad trend of embracing ecological complexity in integrative ecosystem modeling, TAM may offer a consistent framework where modelers and empiricists can work together toward this grand goal.

Shower thoughts: why scientists should spend more time in the rain

Stormwater is a vital resource and dynamic driver of terrestrial ecosystem processes. However, processes controlling interactions during and shortly after storms are often poorly seen and poorly sensed when direct observations are substituted with technological ones. We discuss how human observations complement technological ones and the benefits of scientists spending more time in the storm. Human observation can reveal ephemeral storm-related phenomena such as biogeochemical hot moments, organismal responses, and sedimentary processes that can then be explored in greater resolution using sensors and virtual experiments. Storm-related phenomena trigger lasting, oversized impacts on hydrologic and biogeochemical processes, organismal traits or functions, and ecosystem services at all scales. We provide examples of phenomena in forests, across disciplines and scales, that have been overlooked in past research to inspire mindful, holistic observation of ecosystems during storms. We conclude that technological observations alone are insufficient to trace the process complexity and unpredictability of fleeting biogeochemical or ecological events without the shower thoughts produced by scientists' human sensory and cognitive systems during storms.

Modeling strategies and data needs for representing coastal wetland vegetation in land surface models

Vegetated coastal ecosystems sequester carbon rapidly relative to terrestrial ecosystems. Coastal wetlands are poorly represented in land surface models, but work is underway to improve process-based, predictive modeling of these ecosystems. Here, we identify guiding questions, potential simulations, and data needs to make progress in improving representation of vegetation in terrestrial-aquatic interfaces, with a focus on coastal and estuarine ecosystems. We synthesize relevant plant traits and environmental controls on vegetation that influence carbon cycling in coastal ecosystems. We propose that models include separate plant functional types (PFTs) for mangroves, graminoid salt marshes, and succulent salt marshes to adequately represent the variation in aboveground and belowground productivity between common coastal wetland vegetation types. We also discuss the drivers and carbon storage consequences of shifts in dominant PFTs. We suggest several potential approaches to represent the diversity in vegetation tolerance and adaptations to fluctuations in salinity and water level, which drive key gradients in coastal wetland ecosystems. Finally, we discuss data needs for parameterizing and evaluating model implementations of coastal wetland vegetation types and function.

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

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

Amazon drought resilience - emerging results point to new empirical needs

This article is a Commentary on Costa et al., 237: 714–733.

Linking remotely sensed ecosystem resilience with forest mortality across the continental United States

Episodes of forest mortality have been observed worldwide associated with climate change, impacting species composition and ecosystem services such as water resources and carbon sequestration. Yet our ability to predict forest mortality remains limited, especially across large scales. Time series of satellite imagery has been used to document ecosystem resilience globally, but it is not clear how well remotely sensed resilience can inform the prediction of forest mortality across continental, multi-biome scales. Here, we leverage forest inventories across the continental United States to systematically assess the potential of ecosystem resilience derived using different data sets and methods to predict forest mortality. We found high resilience was associated with low mortality in eastern forests but was associated with high mortality in western regions. The unexpected resilience-mortality relation in western United States may be due to several factors including plant trait acclimation, insect population dynamics, or resource competition. Overall, our results not only supported the opportunity to use remotely sensed ecosystem resilience to predict forest mortality but also highlighted that ecological factors may have crucial influences because they can reverse the sign of the resilience-mortality relationships.

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.

Multiscale predictors of small tree survival across a heterogeneous tropical landscape

Uncertainties about controls on tree mortality make forest responses to land-use and climate change difficult to predict. We tracked biomass of tree functional groups in tropical forest inventories across Puerto Rico and the U.S. Virgin Islands, and with random forests we ranked 86 potential predictors of small tree survival (young or mature stems 2.5-12.6 cm diameter at breast height). Forests span dry to cloud forests, range in age, geology and past land use and experienced severe drought and storms. When excluding species as a predictor, top predictors are tree crown ratio and height, two to three species traits and stand to regional factors reflecting local disturbance and the system state (widespread recovery, drought, hurricanes). Native species, and species with denser wood, taller maximum height, or medium typical height survive longer, but short trees and species survive hurricanes better. Trees survive longer in older stands and with less disturbed canopies, harsher geoclimates (dry, edaphically dry, e.g., serpentine substrates, and highest-elevation cloud forest), or in intervals removed from hurricanes. Satellite image phenology and bands, even from past decades, are top predictors, being sensitive to vegetation type and disturbance. Covariation between stand-level species traits and geoclimate, disturbance and neighboring species types may explain why most neighbor variables, including introduced vs. native species, had low or no importance, despite univariate correlations with survival. As forests recovered from a hurricane in 1998 and earlier deforestation, small trees of introduced species, which on average have lighter wood, died at twice the rate of natives. After hurricanes in 2017, the total biomass of trees ≥12.7 cm dbh of the introduced species Spathodea campanulata spiked, suggesting that more frequent hurricanes might perpetuate this light-wooded species commonness. If hurricane recovery favors light-wooded species while drought favors others, climate change influences on forest composition and ecosystem services may depend on the frequency and severity of extreme climate events.

Long-term hydrological response emerges from forest self-thinning behaviour and tree sapwood allometry

Fires in forested catchments are of great concern to catchment managers due to their potential effect on water yield. Among other factors such as meteorological conditions and topography, dominant vegetation and its regeneration traits can play a key role in controlling the variability in the type and recovery-time of the hydrological response between forested catchments after stand-replacing fires. In temperate South-Eastern Australia, a long-term reduction in streamflow from catchments dominated by regenerating tall-wet Eucalyptus obligate seeder forests was observed, which has substantial implications for Melbourne's water supply. While the unusual hydrological response has been attributed to the higher water-use of the regrowth forests, the dominant underlying mechanism has not yet been identified. Here we show analytically and with a closed-form solution that this streamflow pattern can emerge from forest dynamics, namely the combination of growth and tree mortality as constrained by the self-thinning line (STL) and sapwood allometry of the dominant overstory tree species under non-limiting rainfall regimes. A sensitivity analysis shows that observed variations in the relative streamflow anomaly trend can be explained by parameters controlling: (i) the shape of the STL; (ii) regeneration success; (iii) radial tree growth rate; and (iv) fire severity. We conclude that the observed variation in long-term post-disturbance streamflow behaviour might have resulted from different trajectories of forest dynamics and suggest that to minimize uncertainty in future water-balance predictions, eco-hydrological models for even aged forests include a mechanistic representation of stand demography processes that are constrained by forest inventory data.

Open science priorities for rigorous nature-based climate solutions

Nature-based climate solutions (NbCS) hold promise, but must be based on the best available science to be successful. We outline key ingredients of open data and science crucial for robust and scalable nature-based climate solutions efforts, as an urgent call to action for academic researchers, nongovernmental organizations, government agencies, and private companies.

Using the atmospheric CO2 growth rate to constrain the CO2 flux from land use and land cover change since 1900

We explore the ability of the atmospheric CO₂ record since 1900 to constrain the source of CO₂ from land use and land cover change (hereafter "land use"), taking account of uncertainties in other terms in the global carbon budget. We find that the atmospheric constraint favors land use CO₂ flux estimates with lower decadal variability and can identify potentially erroneous features, such as emission peaks around 1960 and after 2000, in some published estimates. Furthermore, we resolve an offset in the global carbon budget that is most plausibly attributed to the land use flux. This correction shifts the mean land use flux since 1900 across 20 published estimates down by 0.35 PgC year^-1 to 1.04 ± 0.57 PgC year^-1 , which is within the range but at the low end of these estimates. We show that the atmospheric CO₂ record can provide insights into the time history of the land use flux that may reduce uncertainty in this term and improve current understanding and projections of the global carbon cycle.

The evolution, complexity and diversity of models of long-term forest dynamics

To assess the impacts of climate change on vegetation from stand to global scales, models of forest dynamics that include tree demography are needed. Such models are now available for 50 years, but the currently existing diversity of model formulations and its evolution over time are poorly documented. This hampers systematic assessments of structural uncertainties in model-based studies.We conducted a meta-analysis of 28 models, focusing on models that were used in the past five years for climate change studies. We defined 52 model attributes in five groups (basic assumptions, growth, regeneration, mortality and soil moisture) and characterized each model according to these attributes. Analyses of model complexity and diversity included hierarchical cluster analysis and redundancy analysis.Model complexity evolved considerably over the past 50 years. Increases in complexity were largest for growth processes, while complexity of modelled establishment processes increased only moderately. Model diversity was lowest at the global scale, and highest at the landscape scale. We identified five distinct clusters of models, ranging from very simple models to models where specific attribute groups are rendered in a complex manner and models that feature high complexity across all attributes.Most models in use today are not balanced in the level of complexity with which they represent different processes. This is the result of different model purposes, but also reflects legacies in model code, modelers' preferences, and the 'prevailing spirit of the epoch'. The lack of firm theories, laws and 'first principles' in ecology provides high degrees of freedom in model development, but also results in high responsibilities for model developers and the need for rigorous model evaluation. Synthesis. The currently available model diversity is beneficial: convergence in simulations of structurally different models indicates robust projections, while convergence of similar models may convey a false sense of certainty. The existing model diversity-with the exception of global models-can be exploited for improved projections based on multiple models. We strongly recommend balanced further developments of forest models that should particularly focus on establishment and mortality processes, in order to provide robust information for decisions in ecosystem management and policymaking.

Dynamic Energy Budget models: fertile ground for understanding resource allocation in plants in a changing world

Climate change is having dramatic effects on the diversity and distribution of species. Many of these effects are mediated by how an organism's physiological patterns of resource allocation translate into fitness through effects on growth, survival and reproduction. Empirically, resource allocation is challenging to measure directly and so has often been approached using mathematical models, such as Dynamic Energy Budget (DEB) models. The fact that all plants require a very similar set of exogenous resources, namely light, water and nutrients, integrates well with the DEB framework in which a small number of variables and processes linked through pathways represent an organism's state as it changes through time. Most DEB theory has been developed in reference to animals and microorganisms. However, terrestrial vascular plants differ from these organisms in fundamental ways that make resource allocation, and the trade-offs and feedbacks arising from it, particularly fundamental to their life histories, but also challenging to represent using existing DEB theory. Here, we describe key features of the anatomy, morphology, physiology, biochemistry, and ecology of terrestrial vascular plants that should be considered in the development of a generic DEB model for plants. We then describe possible approaches to doing so using existing DEB theory and point out features that may require significant development for DEB theory to accommodate them. We end by presenting a generic DEB model for plants that accounts for many of these key features and describing gaps that would need to be addressed for DEB theory to predict the responses of plants to climate change. DEB models offer a powerful and generalizable framework for modelling resource allocation in terrestrial vascular plants, and our review contributes a framework for expansion and development of DEB theory to address how plants respond to anthropogenic change.

A climate risk analysis of Earth's forests in the 21st century

Earth's forests harbor extensive biodiversity and are currently a major carbon sink. Forest conservation and restoration can help mitigate climate change; however, climate change could fundamentally imperil forests in many regions and undermine their ability to provide such mitigation. The extent of climate risks facing forests has not been synthesized globally nor have different approaches to quantifying forest climate risks been systematically compared. We combine outputs from multiple mechanistic and empirical approaches to modeling carbon, biodiversity, and disturbance risks to conduct a synthetic climate risk analysis for Earth's forests in the 21st century. Despite large uncertainty in most regions we find that some forests are consistently at higher risk, including southern boreal forests and those in western North America and parts of the Amazon.

Tree crown damage and its effects on forest carbon cycling in a tropical forest

Crown damage can account for over 23% of canopy biomass turnover in tropical forests and is a strong predictor of tree mortality; yet, it is not typically represented in vegetation models. We incorporate crown damage into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to evaluate how lags between damage and tree recovery or death alter demographic rates and patterns of carbon turnover. We represent crown damage as a reduction in a tree's crown area and leaf and branch biomass, and allow associated variation in the ratio of aboveground to belowground plant tissue. We compare simulations with crown damage to simulations with equivalent instant increases in mortality and benchmark results against data from Barro Colorado Island (BCI), Panama. In FATES, crown damage causes decreases in growth rates that match observations from BCI. Crown damage leads to increases in carbon starvation mortality in FATES, but only in configurations with high root respiration and decreases in carbon storage following damage. Crown damage also alters competitive dynamics, as plant functional types that can recover from crown damage outcompete those that cannot. This is a first exploration of the trade-off between the additional complexity of the novel crown damage module and improved predictive capabilities. At BCI, a tropical forest that does not experience high levels of disturbance, both the crown damage simulations and simulations with equivalent increases in mortality does a reasonable job of capturing observations. The crown damage module provides functionality for exploring dynamics in forests with more extreme disturbances such as cyclones and for capturing the synergistic effects of disturbances that overlap in space and time.

A Conceptual Framework to Integrate Biodiversity, Ecosystem Function, and Ecosystem Service Models

Global biodiversity and ecosystem service models typically operate independently. Ecosystem service projections may therefore be overly optimistic because they do not always account for the role of biodiversity in maintaining ecological functions. We review models used in recent global model intercomparison projects and develop a novel model integration framework to more fully account for the role of biodiversity in ecosystem function, a key gap for linking biodiversity changes to ecosystem services. We propose two integration pathways. The first uses empirical data on biodiversity-ecosystem function relationships to bridge biodiversity and ecosystem function models and could currently be implemented globally for systems and taxa with sufficient data. We also propose a trait-based approach involving greater incorporation of biodiversity into ecosystem function models. Pursuing both approaches will provide greater insight into biodiversity and ecosystem services projections. Integrating biodiversity, ecosystem function, and ecosystem service modeling will enhance policy development to meet global sustainability goals.

Summer temperature—but not growing season length—influences radial growth of Salix arctica in coastal Arctic tundra

Arctic climate change is leading to an advance of plant phenology (the timing of life history events) with uncertain impacts on tundra ecosystems. Although the lengthening of the growing season is thought to lead to increased plant growth, we have few studies of how plant phenology change is altering tundra plant productivity. Here, we test the correspondence between 14 years of Salix arctica phenology data and radial growth on Qikiqtaruk–Herschel Island, Yukon Territory, Canada. We analysed stems from 28 individuals using dendroecology and linear mixed-effect models to test the statistical power of growing season length and climate variables to individually predict radial growth. We found that summer temperature best explained annual variation in radial growth. We found no strong evidence that leaf emergence date, earlier leaf senescence date, or total growing season length had any direct or lagged effects on radial growth. Radial growth was also not explained by interannual variation in precipitation, MODIS surface greenness (NDVI), or sea ice concentration. Our results demonstrate that at this site, for the widely distributed species S. arctica, temperature—but not growing season length—influences radial growth. These findings challenge the assumption that advancing phenology and longer growing seasons will increase the productivity of all plant species in Arctic tundra ecosystems.

Forest regeneration within Earth system models: current process representations and ways forward

Earth system models must predict forest responses to global change in order to simulate future global climate, hydrology, and ecosystem dynamics. These models are increasingly adopting vegetation demographic approaches that explicitly represent tree growth, mortality, and recruitment, enabling advances in the projection of forest vulnerability and resilience, as well as evaluation with field data. To date, simulation of regeneration processes has received far less attention than simulation of processes that affect growth and mortality, in spite of their critical role maintaining forest structure, facilitating turnover in forest composition over space and time, enabling recovery from disturbance, and regulating climate-driven range shifts. Our critical review of regeneration process representations within current Earth system vegetation demographic models reveals the need to improve parameter values and algorithms for reproductive allocation, dispersal, seed survival and germination, environmental filtering in the seedling layer, and tree regeneration strategies adapted to wind, fire, and anthropogenic disturbance regimes. These improvements require synthesis of existing data, specific field data-collection protocols, and novel model algorithms compatible with global-scale simulations. Vegetation demographic models offer the opportunity to more fully integrate ecological understanding into Earth system prediction; regeneration processes need to be a critical part of the effort.

Simulating environmentally-sensitive tree recruitment in vegetation demographic models

Vegetation demographic models (VDMs) endeavor to predict how global forests will respond to climate change. This requires simulating which trees, if any, are able to recruit under changing environmental conditions. We present a new recruitment scheme for VDMs in which functional-type-specific recruitment rates are sensitive to light, soil moisture and the productivity of reproductive trees. We evaluate the scheme by predicting tree recruitment for four tropical tree functional types under varying meteorology and canopy structure at Barro Colorado Island, Panama. We compare predictions to those of a current VDM, quantitative observations and ecological expectations. We find that the scheme improves the magnitude and rank order of recruitment rates among functional types and captures recruitment limitations in response to variable understory light, soil moisture and precipitation regimes. Our results indicate that adopting this framework will improve VDM capacity to predict functional-type-specific tree recruitment in response to climate change, thereby improving predictions of future forest distribution, composition and function.