An individual-based forest model to jointly simulate carbon and tree diversity in Amazonia: description and applications

Toward a coordinated understanding of hydro-biogeochemical root functions in tropical forests for application in vegetation models

Tropical forest root characteristics and resource acquisition strategies are underrepresented in vegetation and global models, hampering the prediction of forest-climate feedbacks for these carbon-rich ecosystems. Lowland tropical forests often have globally unique combinations of high taxonomic and functional biodiversity, rainfall seasonality, and strongly weathered infertile soils, giving rise to distinct patterns in root traits and functions compared with higher latitude ecosystems. We provide a roadmap for integrating recent advances in our understanding of tropical forest belowground function into vegetation models, focusing on water and nutrient acquisition. We offer comparisons of recent advances in empirical and model understanding of root characteristics that represent important functional processes in tropical forests. We focus on: (1) fine-root strategies for soil resource exploration, (2) coupling and trade-offs in fine-root water vs nutrient acquisition, and (3) aboveground-belowground linkages in plant resource acquisition and use. We suggest avenues for representing these extremely diverse plant communities in computationally manageable and ecologically meaningful groups in models for linked aboveground-belowground hydro-nutrient functions. Tropical forests are undergoing warming, shifting rainfall regimes, and exacerbation of soil nutrient scarcity caused by elevated atmospheric CO2. The accurate model representation of tropical forest functions is crucial for understanding the interactions of this biome with the climate.

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.

Deciphering the fingerprint of disturbance on the three-dimensional structure of the world's forests

Canopy gaps and the processes that generate them play an integral role in shaping the structure and dynamics of forests. However, it is only with recent advances in remote sensing technologies such as airborne laser scanning that studying canopy gaps at scale has become a reality. Consequently, we still lack an understanding of how the size distribution and spatial organization of canopy gaps varies among forests ecosystems, nor have we determined whether these emergent properties can be reconciled with existing theories of forest dynamics. Here, I outline a roadmap for integrating remote sensing with field data and individual-based models to build a comprehensive picture of how environmental constraints and disturbance regimes shape the three-dimensional structure of the world's forests.

Nutrient cycling drives plant community trait assembly and ecosystem functioning in a tropical mountain biodiversity hotspot

Community trait assembly in highly diverse tropical rainforests is still poorly understood. Based on more than a decade of field measurements in a biodiversity hotspot of southern Ecuador, we implemented plant trait variation and improved soil organic matter dynamics in a widely used dynamic vegetation model (the Lund-Potsdam-Jena General Ecosystem Simulator, LPJ-GUESS) to explore the main drivers of community assembly along an elevational gradient. In the model used here (LPJ-GUESS-NTD, where NTD stands for nutrient-trait dynamics), each plant individual can possess different trait combinations, and the community trait composition emerges via ecological sorting. Further model developments include plant growth limitation by phosphorous (P) and mycorrhizal nutrient uptake. The new model version reproduced the main observed community trait shift and related vegetation processes along the elevational gradient, but only if nutrient limitations to plant growth were activated. In turn, when traits were fixed, low productivity communities emerged due to reduced nutrient-use efficiency. Mycorrhizal nutrient uptake, when deactivated, reduced net primary production (NPP) by 61-72% along the gradient. Our results strongly suggest that the elevational temperature gradient drives community assembly and ecosystem functioning indirectly through its effect on soil nutrient dynamics and vegetation traits. This illustrates the importance of considering these processes to yield realistic model predictions.

Hydraulically-vulnerable trees survive on deep-water access during droughts in a tropical forest

Deep-water access is arguably the most effective, but under-studied, mechanism that plants employ to survive during drought. Vulnerability to embolism and hydraulic safety margins can predict mortality risk at given levels of dehydration, but deep-water access may delay plant dehydration. Here, we tested the role of deep-water access in enabling survival within a diverse tropical forest community in Panama using a novel data-model approach. We inversely estimated the effective rooting depth (ERD, as the average depth of water extraction), for 29 canopy species by linking diameter growth dynamics (1990-2015) to vapor pressure deficit, water potentials in the whole-soil column, and leaf hydraulic vulnerability curves. We validated ERD estimates against existing isotopic data of potential water-access depths. Across species, deeper ERD was associated with higher maximum stem hydraulic conductivity, greater vulnerability to xylem embolism, narrower safety margins, and lower mortality rates during extreme droughts over 35 years (1981-2015) among evergreen species. Species exposure to water stress declined with deeper ERD indicating that trees compensate for water stress-related mortality risk through deep-water access. The role of deep-water access in mitigating mortality of hydraulically-vulnerable trees has important implications for our predictive understanding of forest dynamics under current and future climates.

The function-dominance correlation drives the direction and strength of biodiversity-ecosystem functioning relationships

Community composition is a primary determinant of how biodiversity change influences ecosystem functioning and, therefore, the relationship between biodiversity and ecosystem functioning (BEF). We examine the consequences of community composition across six structurally realistic plant community models. We find that a positive correlation between species' functioning in monoculture versus their dominance in mixture with regard to a specific function (the "function-dominance correlation") generates a positive relationship between realised diversity and ecosystem functioning across species richness treatments. However, because realised diversity declines when few species dominate, a positive function-dominance correlation generates a negative relationship between realised diversity and ecosystem functioning within species richness treatments. Removing seed inflow strengthens the link between the function-dominance correlation and BEF relationships across species richness treatments but weakens it within them. These results suggest that changes in species' identities in a local species pool may more strongly affect ecosystem functioning than changes in species richness.

Tackling unresolved questions in forest ecology: The past and future role of simulation models

Understanding the processes that shape forest functioning, structure, and diversity remains challenging, although data on forest systems are being collected at a rapid pace and across scales. Forest models have a long history in bridging data with ecological knowledge and can simulate forest dynamics over spatio-temporal scales unreachable by most empirical investigations.We describe the development that different forest modelling communities have followed to underpin the leverage that simulation models offer for advancing our understanding of forest ecosystems.Using three widely applied but contrasting approaches - species distribution models, individual-based forest models, and dynamic global vegetation models - as examples, we show how scientific and technical advances have led models to transgress their initial objectives and limitations. We provide an overview of recent model applications on current important ecological topics and pinpoint ten key questions that could, and should, be tackled with forest models in the next decade.Synthesis. This overview shows that forest models, due to their complementarity and mutual enrichment, represent an invaluable toolkit to address a wide range of fundamental and applied ecological questions, hence fostering a deeper understanding of forest dynamics in the context of global change.

Continent-wide tree fecundity driven by indirect climate effects

Indirect climate effects on tree fecundity that come through variation in size and growth (climate-condition interactions) are not currently part of models used to predict future forests. Trends in species abundances predicted from meta-analyses and species distribution models will be misleading if they depend on the conditions of individuals. Here we find from a synthesis of tree species in North America that climate-condition interactions dominate responses through two pathways, i) effects of growth that depend on climate, and ii) effects of climate that depend on tree size. Because tree fecundity first increases and then declines with size, climate change that stimulates growth promotes a shift of small trees to more fecund sizes, but the opposite can be true for large sizes. Change the depresses growth also affects fecundity. We find a biogeographic divide, with these interactions reducing fecundity in the West and increasing it in the East. Continental-scale responses of these forests are thus driven largely by indirect effects, recommending management for climate change that considers multiple demographic rates.

An overview of agent-based models in plant biology and ecology

Agent-based modelling (ABM) has become an established methodology in many areas of biology, ranging from the cellular to the ecological population and community levels. In plant science, two different scales have predominated in their use of ABM. One is the scale of populations and communities, through the modelling of collections of agents representing individual plants, interacting with each other and with the environment. The other is the scale of the individual plant, through the modelling, by functional-structural plant models (FSPMs), of agents representing plant building blocks, or metamers, to describe the development of plant architecture and functions within individual plants. The purpose of this review is to show key results and parallels in ABM for growth, mortality, carbon allocation, competition and reproduction across the scales from the plant organ to populations and communities on a range of spatial scales to the whole landscape. Several areas of application of ABMs are reviewed, showing that some issues are addressed by both population-level ABMs and FSPMs. Continued increase in the relevance of ABM to environmental science and management will be helped by greater integration of ABMs across these two scales.

Allometric constraints and competition enable the simulation of size structure and carbon fluxes in a dynamic vegetation model of tropical forests (LM3PPA-TV)

Tropical forests are a key determinant of the functioning of the Earth system, but remain a major source of uncertainty in carbon cycle models and climate change projections. In this study, we present an updated land model (LM3PPA-TV) to improve the representation of tropical forest structure and dynamics in Earth system models (ESMs). The development and parameterization of LM3PPA-TV drew on extensive datasets on tropical tree traits and long-term field censuses from Barro Colorado Island (BCI), Panama. The model defines a new plant functional type (PFT) based on the characteristics of shade-tolerant, tropical tree species, implements a new growth allocation scheme based on realistic tree allometries, incorporates hydraulic constraints on biomass accumulation, and features a new compartment for tree branches and branch fall dynamics. Simulation experiments reproduced observed diurnal and seasonal patterns in stand-level carbon and water fluxes, as well as mean canopy and understory tree growth rates, tree size distributions, and stand-level biomass on BCI. Simulations at multiple sites captured considerable variation in biomass and size structure across the tropical forest biome, including observed responses to precipitation and temperature. Model experiments suggested a major role of water limitation in controlling geographic variation forest biomass and structure. However, the failure to simulate tropical forests under extreme conditions and the systematic underestimation of forest biomass in Paleotropical locations highlighted the need to incorporate variation in hydraulic traits and multiple PFTs that capture the distinct floristic composition across tropical domains. The continued pressure on tropical forests from global change demands models which are able to simulate alternative successional pathways and their pace to recovery. LM3PPA-TV provides a tool to investigate geographic variation in tropical forests and a benchmark to continue improving the representation of tropical forests dynamics and their carbon storage potential in ESMs.

Impacts of Degradation on Water, Energy, and Carbon Cycling of the Amazon Tropical Forests

Selective logging, fragmentation, and understory fires directly degrade forest structure and composition. However, studies addressing the effects of forest degradation on carbon, water, and energy cycles are scarce. Here, we integrate field observations and high-resolution remote sensing from airborne lidar to provide realistic initial conditions to the Ecosystem Demography Model (ED-2.2) and investigate how disturbances from forest degradation affect gross primary production (GPP), evapotranspiration (ET), and sensible heat flux (H). We used forest structural information retrieved from airborne lidar samples (13,500 ha) and calibrated with 817 inventory plots (0.25 ha) across precipitation and degradation gradients in the eastern Amazon as initial conditions to ED-2.2 model. Our results show that the magnitude and seasonality of fluxes were modulated by changes in forest structure caused by degradation. During the dry season and under typical conditions, severely degraded forests (biomass loss ≥66%) experienced water stress with declines in ET (up to 34%) and GPP (up to 35%) and increases of H (up to 43%) and daily mean ground temperatures (up to 6.5°C) relative to intact forests. In contrast, the relative impact of forest degradation on energy, water, and carbon cycles markedly diminishes under extreme, multiyear droughts, as a consequence of severe stress experienced by intact forests. Our results highlight that the water and energy cycles in the Amazon are driven by not only climate and deforestation but also the past disturbance and changes of forest structure from degradation, suggesting a much broader influence of human land use activities on the tropical ecosystems.

Slow rate of secondary forest carbon accumulation in the Guianas compared with the rest of the Neotropics

Secondary forests are a prominent component of tropical landscapes, and they constitute a major atmospheric carbon sink. Rates of carbon accumulation are usually inferred from chronosequence studies, but direct estimates of carbon accumulation based on long-term monitoring of stands are rarely reported. Recent compilations on secondary forest carbon accumulation in the Neotropics are heavily biased geographically as they do not include estimates from the Guiana Shield. We analysed the temporal trajectory of aboveground carbon accumulation and floristic composition at one 25-ha secondary forest site in French Guiana. The site was clear-cut in 1976, abandoned thereafter, and one large plot (6.25 ha) has been monitored continuously since. We used Bayesian modeling to assimilate inventory data and simulate the long-term carbon accumulation trajectory. Canopy change was monitored using two aerial lidar surveys conducted in 2009 and 2017. We compared the dynamics of this site with that of a surrounding old-growth forest. Finally, we compared our results with that from secondary forests in Costa Rica, which is one of the rare long-term monitoring programs reaching a duration comparable to our study. Twenty years after abandonment, aboveground carbon stock was 64.2 (95% credibility interval 46.4, 89.0) Mg C/ha, and this stock increased to 101.3 (78.7, 128.5) Mg C/ha 20 yr later. The time to accumulate one-half of the mean aboveground carbon stored in the nearby old-growth forest (185.6 [155.9, 200.2] Mg C/ha) was estimated at 35.0 [20.9, 55.9] yr. During the first 40 yr, the contribution of the long-lived pioneer species Xylopia nitida, Goupia glabra, and Laetia procera to the aboveground carbon stock increased continuously. Secondary forest mean-canopy height measured by lidar increased by 1.14 m in 8 yr, a canopy-height increase consistent with an aboveground carbon accumulation of 7.1 Mg C/ha (or 0.89 Mg C·ha^-1 ·yr^-1 ) during this period. Long-term AGC accumulation rate in Costa Rica was almost twice as fast as at our site in French Guiana. This may reflect higher fertility of Central American forest communities or a better adaptation of the forest tree community to intense and frequent disturbances. This finding may have important consequences for scaling-up carbon uptake estimates to continental scales.

Improving plant allometry by fusing forest models and remote sensing

Allometry determines how tree shape and function scale with each other, related through size. Allometric relationships help scale processes from the individual to the global scale and constitute a core component of vegetation models. Allometric relationships have been expected to emerge from optimisation theory, yet this does not suitably predict empirical data. Here we argue that the fusion of high-resolution data, such as those derived from airborne laser scanning, with individual-based forest modelling offers insight into how plant size contributes to large-scale biogeochemical processes. We review the challenges in allometric scaling, how they can be tackled by advances in data-model fusion, and how individual-based models can serve as data integrators for dynamic global vegetation models.

A Wood Biology Agenda to Support Global Vegetation Modelling

Realistic forecasting of forest responses to climate change critically depends on key advancements in global vegetation modelling. Compared with traditional 'big-leaf' models that simulate forest stands, 'next-generation' vegetation models aim to track carbon-, light-, water-, and nutrient-limited growth of individual trees. Wood biology can play an important role in delivering the required knowledge at tissue-to-individual levels, at minute-to-century scales and for model parameterization and benchmarking. We propose a wood biology research agenda that contributes to filling six knowledge gaps: sink versus source limitation, drivers of intra-annual growth, drought impacts, functional wood traits, dynamic biomass allocation, and nutrient cycling. Executing this agenda will expedite model development and increase the ability of models to forecast global change impact on forest dynamics.