Observed variation in soil properties can drive large variation in modelled forest functioning and composition during tropical forest secondary succession

Campoe O, Deshar R, Franklin O, Henriksson N, Littke K, Lutter R, Maier CA, Palmroth S, Rosenvald K, Slesak RA, Tullus A, Oren R. Overlooked branch turnover creates a widespread bias in forest carbon accounting

Most measurements and models of forest carbon cycling neglect the carbon flux associated with the turnover of branch biomass, a physiological process quantified for other organs (fine roots, leaves, and stems). Synthesizing data from boreal, temperate, and tropical forests (184,815 trees), we found that including branch turnover increased empirical estimates of aboveground wood production by 16% (equivalent to 1.9 Pg Cy-1 globally), of similar magnitude to the observed global forest carbon sinks. In addition, reallocating carbon to branch turnover in model simulations reduced stem wood biomass, a long-lasting carbon storage, by 7 to 17%. This prevailing neglect of branch turnover suggests widespread biases in carbon flux estimates across global datasets and model simulations. Branch litterfall, sometimes used as a proxy for branch turnover, ignores carbon lost from attached dead branches, underestimating branch C turnover by 38% in a pine forest. Modifications to field measurement protocols and existing models are needed to allow a more realistic partitioning of wood production and forest carbon storage.

Plant responses to elevated CO2 under competing hypotheses of nitrogen and phosphorus limitations

The future ecosystem carbon cycle has important implications for biosphere-climate feedback. The magnitude of future plant growth and carbon accumulation depends on plant strategies for nutrient uptake under the stresses of nitrogen (N) versus phosphorus (P) limitations. Two archetypal theories have been widely acknowledged in the literature to represent N and P limitations on ecosystem processes: Liebig's Law of the Minimum (LLM) and the Multiple Element Limitation (MEL) approach. LLM states that the more limiting nutrient controls plant growth, and commonly leads to predictions of dramatically dampened ecosystem carbon accumulation over the 21st century. Conversely, the MEL approach recognizes that plants possess multiple pathways to coordinate N and P availability and invest resources to alleviate N or P limitation. We implemented these two contrasting approaches in the E3SM model, and compiled 98 in situ forest N or P fertilization experiments to evaluate how terrestrial ecosystems will respond to N and P limitations. We find that MEL better captured the observed plant responses to nutrient perturbations globally, compared with LLM. Furthermore, LLM and MEL diverged dramatically in responses to elevated CO2 concentrations, leading to a two-fold difference in CO2 fertilization effects on Net Primary Productivity by the end of the 21st century. The larger CO2 fertilization effects indicated by MEL mainly resulted from plant mediation on N and P resource supplies through N2 fixation and phosphatase activities. This analysis provides quantitative evidence of how different N and P limitation strategies can diversely affect future carbon and nutrient dynamics.

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.

G, Soper FM, Gutknecht J, Powers JS. Phosphorus limitation of early growth differs between nitrogen-fixing and nonfixing dry tropical forest tree species

Tropical forests are often characterized by low soil phosphorus (P) availability, suggesting that P limits plant performance. However, how seedlings from different functional types respond to soil P availability is poorly known but important for understanding and modeling forest dynamics under changing environmental conditions. We grew four nitrogen (N)-fixing Fabaceae and seven diverse non-N-fixing tropical dry forest tree species in a shade house under three P fertilization treatments and evaluated carbon (C) allocation responses, P demand, P-use, investment in P acquisition traits, and correlations among P acquisition traits. Nitrogen fixers grew larger with increasing P addition in contrast to non-N fixers, which showed fewer responses in C allocation and P use. Foliar P increased with P addition for both functional types, while P acquisition strategies did not vary among treatments but differed between functional types, with N fixers showing higher root phosphatase activity (RPA) than nonfixers. Growth responses suggest that N fixers are limited by P, but nonfixers may be limited by other resources. However, regardless of limitation, P acquisition traits such as mycorrhizal colonization and RPA were nonplastic across a steep P gradient. Differential limitation among plant functional types has implications for forest succession and earth system models.

Increasing calcium scarcity along Afrotropical forest succession

Secondary forests constitute an increasingly important component of tropical forests worldwide. Although cycling of essential nutrients affects recovery trajectories of secondary forests, the effect of nutrient limitation on forest regrowth is poorly constrained. Here we use three lines of evidence from secondary forest succession sequences in central Africa to identify potential nutrient limitation in regrowing forests. First, we show that atmospheric phosphorus supply exceeds demand along forest succession, whereas forests rely on soil stocks to meet their base cation demands. Second, soil nutrient metrics indicate that available phosphorus increases along the succession, whereas available cations decrease. Finally, fine root, foliar and litter stoichiometry show that tissue calcium concentrations decline relative to those of nitrogen and phosphorus during succession. Taken together, these observations suggest that calcium becomes an increasingly scarce resource in central African forests during secondary succession. Furthermore, ecosystem calcium storage shifts from soil to woody biomass over succession, making it a vulnerable nutrient in the wake of land-use change scenarios that involve woody biomass export. Our results thus call for a broadened focus on elements other than nitrogen and phosphorus regarding tropical forest biogeochemical cycles and identify calcium as a scarce and potentially limiting nutrient in an increasingly disturbed and dynamic tropical forest landscape.

Climate and hydraulic traits interact to set thresholds for liana viability

Lianas, or woody vines, and trees dominate the canopy of tropical forests and comprise the majority of tropical aboveground carbon storage. These growth forms respond differently to contemporary variation in climate and resource availability, but their responses to future climate change are poorly understood because there are very few predictive ecosystem models representing lianas. We compile a database of liana functional traits (846 species) and use it to parameterize a mechanistic model of liana-tree competition. The substantial difference between liana and tree hydraulic conductivity represents a critical source of inter-growth form variation. Here, we show that lianas are many times more sensitive to drying atmospheric conditions than trees as a result of this trait difference. Further, we use our competition model and projections of tropical hydroclimate based on Representative Concentration Pathway 4.5 to show that lianas are more susceptible to reaching a hydraulic threshold for viability by 2100.

A Pantropical Overview of Soils across Tropical Dry Forest Ecoregions

Pantropical variation in soils of the tropical dry forest (TDF) biome is enormously high but has been poorly characterized. To quantify variation in the global distribution of TDF soil physical and chemical properties in relation to climate and geology, we produced a synthesis using 7500 points of data with gridded fields representing lithologic, edaphic, and climatic characteristics. Our analyses reveal that 75 TDF ecoregions across five biogeographic domains (Afrotropical, Australasian, Indo-Malayan, Neotropical, and Oceanian) varied strongly with respect to parent material: sediment (57%), metamorphic (22%), volcanic (13%), and plutonic (7%). TDF ecoregions support remarkably high variability in soil suborders (32), with the Neotropical and Oceanian realms being especially diverse. As a whole, TDF soils trend strongly toward low fertility with strong variation across biogeographic domains. Similarly, the exhibited soil properties marked heterogeneity across biogeographic domains, with soil depth varying by an order of magnitude and total organic C, N, and P pools varying threefold. Organic C and N pool sizes were negatively correlated with mean annual temperature (MAT) and positively correlated with mean annual precipitation (MAP). By contrast, the distribution of soil P pools was positively influenced by both MAT and MAP and likely by soil geochemistry, due to high variations in soil parent material across the biogeographic domains. The results summarized here raise important questions as to how climate and parent material control soil biogeochemical processes in TDFs.

Faunal communities mediate the effects of plant richness, drought, and invasion on ecosystem multifunctional stability

Understanding the stability of ecosystem multifunctionality is imperative for maintaining ecosystem health and sustainability under augmented global change. However it remains unknown whether and how biological communities mediate multifunctional stability in response to biodiversity loss and disturbances. Here, we conducted a 3-year experiment by exposing 270 plant communities of four plant richness levels, i.e., 1, 2, 4, or 8 species, to drought and exotic plant invasion disturbances. Then, the direct effects of plant richness, drought and invasion, and their indirect effects mediated by the stability of plant, litter-faunal, and soil-faunal communities on multifunctional stability were disentangled. We found that plant richness increased, while drought and invasion decreased ecosystem multifunctional stability, which were mediated by plant or faunal community stability. By incorporating the stability of communities into the complex ecological mechanisms, the completeness and goodness of ecological models for explaining and maintaining the stability of ecosystem multifunctionality will be improved.

Nitrification, denitrification, and competition for soil N: Evaluation of two Earth System Models against observations

Earth System Models (ESMs) have implemented nitrogen (N) cycles to account for N limitation on terrestrial carbon uptake. However, representing inputs, losses, and recycling of N in ESMs is challenging. Here, we use global rates and ratios of key soil N fluxes, including nitrification, denitrification, mineralization, leaching, immobilization, and plant uptake (both NH₄ ⁺ and NO₃ ^- ), from the literature to evaluate the N cycles in the land model components of two ESMs. The two land models evaluated here, E3SM Land Model version 1 (ELMv1)-ECA and CLM5.0, originated from a common model but have diverged in their representation of plant-microbe competition for soil N. The models predict similar global rates of gross primary productivity (GPP) but have approximately two-fold to three-fold differences in their underlying global mineralization, immobilization, plant N uptake, nitrification, and denitrification fluxes. Both models dramatically underestimate the immobilization of NO₃ ^- by soil bacteria compared with literature values and predict dominance of plant uptake by a single form of mineral nitrogen (NO₃ ^- for ELM, with regional exceptions, and NH₄ ⁺ for CLM5.0). CLM5.0 strongly underestimates the global ratio of gross nitrification:gross mineralization and both models are likely to substantially underestimate the ratio of nitrification:denitrification. Few experimental data exist to evaluate this last ratio, in part because nitrification and denitrification are quantified using different techniques and because denitrification fluxes are difficult to measure at all. More observational constraints on soil nitrogen fluxes such as nitrification and denitrification, as well as greater scrutiny of the functional impact of introducing separate NH₄ ⁺ and NO₃ ^- pools into ESMs, could help to improve confidence in present and future simulations of N limitation on the carbon cycle.

Plant Diversity and Soil Nutrients in a Tropical Coastal Secondary Forest: Association Ordination and Sampling Year Differences

Studying the patterns of changes in species diversity and soil properties can improve our knowledge of community succession. However, there is still a gap in understanding how soil conditions are related to plant diversity in tropical coastal secondary forests. We sampled plant diversity and soil nutrients spanning two different years (2012 and 2019) to assess the patterns of species diversity and relationships of soil nutrients and species diversity on Hainan Island, southern China. Results showed that the soil pH and total nitrogen (TN) significantly decreased while the soil organic matter (OM) and total phosphorus (TP) significantly increased from 2012 to 2019. Plant species diversity was significantly higher in 2012 than in 2019, and the dominant species significantly changed in two different years. Using multiple regression analysis, we determined that soil TP and TN were significantly related to plant diversity in 2012 and 2019, respectively. Using CCA analysis, TN and OM were the strongest predictors for dominant species in 2012, whereas the soil TP and TN were the strongest predictors for dominant species in 2019. Our findings show a significant change in plant diversity and dominant species after 7 years of development in the tropical coastal secondary forest. The patterns of plant diversity and soil nutrients increase our knowledge of forest restoration in coastal areas.

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.

Patterns and mechanisms of spatial variation in tropical forest productivity, woody residence time, and biomass

Tropical forests vary widely in biomass carbon (C) stocks and fluxes even after controlling for forest age. A mechanistic understanding of this variation is critical to accurately predicting responses to global change. We review empirical studies of spatial variation in tropical forest biomass, productivity and woody residence time, focusing on mature forests. Woody productivity and biomass decrease from wet to dry forests and with elevation. Within lowland forests, productivity and biomass increase with temperature in wet forests, but decrease with temperature where water becomes limiting. Woody productivity increases with soil fertility, whereas residence time decreases, and biomass responses are variable, consistent with an overall unimodal relationship. Areas with higher disturbance rates and intensities have lower woody residence time and biomass. These environmental gradients all involve both direct effects of changing environments on forest C fluxes and shifts in functional composition - including changing abundances of lianas - that substantially mitigate or exacerbate direct effects. Biogeographic realms differ significantly and importantly in productivity and biomass, even after controlling for climate and biogeochemistry, further demonstrating the importance of plant species composition. Capturing these patterns in global vegetation models requires better mechanistic representation of water and nutrient limitation, plant compositional shifts and tree mortality.

Unraveling the relative role of light and water competition between lianas and trees in tropical forests: A vegetation model analysis

Despite their low contribution to forest carbon stocks, lianas (woody vines) play an important role in the carbon dynamics of tropical forests. As structural parasites, they hinder tree survival, growth and fecundity; hence, they negatively impact net ecosystem productivity and long-term carbon sequestration.Competition (for water and light) drives various forest processes and depends on the local abundance of resources over time. However, evaluating the relative role of resource availability on the interactions between lianas and trees from empirical observations is particularly challenging. Previous approaches have used labour-intensive and ecosystem-scale manipulation experiments, which are infeasible in most situations.We propose to circumvent this challenge by evaluating the uncertainty of water and light capture processes of a process-based vegetation model (ED2) including the liana growth form. We further developed the liana plant functional type in ED2 to mechanistically simulate water uptake and transport from roots to leaves, and start the model from prescribed initial conditions. We then used the PEcAn bioinformatics platform to constrain liana parameters and run uncertainty analyses.Baseline runs successfully reproduced ecosystem gas exchange fluxes (gross primary productivity and latent heat) and forest structural features (leaf area index, aboveground biomass) in two sites (Barro Colorado Island, Panama and Paracou, French Guiana) characterized by different rainfall regimes and levels of liana abundance.Model uncertainty analyses revealed that water limitation was the factor driving the competition between trees and lianas at the drier site (BCI), and during the relatively short dry season of the wetter site (Paracou). In young patches, light competition dominated in Paracou but alternated with water competition between the wet and the dry season on BCI according to the model simulations.The modelling workflow also identified key liana traits (photosynthetic quantum efficiency, stomatal regulation parameters, allometric relationships) and processes (water use, respiration, climbing) driving the model uncertainty. They should be considered as priorities for future data acquisition and model development to improve predictions of the carbon dynamics of liana-infested forests. Synthesis. Competition for water plays a larger role in the interaction between lianas and trees than previously hypothesized, as demonstrated by simulations from a process-based vegetation model.

Allometric scaling laws linking biomass and rooting depth vary across ontogeny and functional groups in tropical dry forest lianas and trees

There are two theories about how allocation of metabolic products occurs. The allometric biomass partitioning theory (APT) suggests that all plants follow common allometric scaling rules. The optimal partitioning theory (OPT) predicts that plants allocate more biomass to the organ capturing the most limiting resource. Whole-plant harvests of mature and juvenile tropical deciduous trees, evergreen trees, and lianas and model simulations were used to address the following knowledge gaps: (1) Do mature lianas comply with the APT scaling laws or do they invest less biomass in stems compared to trees? (2) Do juveniles follow the same allocation patterns as mature individuals? (3) Is either leaf phenology or life form a predictor of rooting depth? It was found that: (1) mature lianas followed the same allometric scaling laws as trees; (2) juveniles and mature individuals do not follow the same allocation patterns; and (3) mature lianas had shallowest coarse roots and evergreen trees had the deepest. It was demonstrated that: (1) mature lianas invested proportionally similar biomass to stems as trees and not less, as expected; (2) lianas were not deeper-rooted than trees as had been previously proposed; and (3) evergreen trees had the deepest roots, which is necessary to maintain canopy during simulated dry seasons.

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.

Leveraging plant hydraulics to yield predictive and dynamic plant leaf allocation in vegetation models with climate change

Plant functional traits provide a link in process-based vegetation models between plant-level physiology and ecosystem-level responses. Recent advances in physiological understanding and computational efficiency have allowed for the incorporation of plant hydraulic processes in large-scale vegetation models. However, a more mechanistic representation of water limitation that determines ecosystem responses to plant water stress necessitates a re-evaluation of trait-based constraints for plant carbon allocation, particularly allocation to leaf area. In this review, we examine model representations of plant allocation to leaves, which is often empirically set by plant functional type-specific allometric relationships. We analyze the evolution of the representation of leaf allocation in models of different scales and complexities. We show the impacts of leaf allocation strategy on plant carbon uptake in the context of recent advancements in modeling hydraulic processes. Finally, we posit that deriving allometry from first principles using mechanistic hydraulic processes is possible and should become standard practice, rather than using prescribed allometries. The representation of allocation as an emergent property of scarce resource constraints is likely to be critical to representing how global change processes impact future ecosystem dynamics and carbon fluxes and may reduce the number of poorly constrained parameters in vegetation models.