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.

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.

Reduced ecosystem resilience quantifies fine-scale heterogeneity in tropical forest mortality responses to drought.

Sensitivity of forest mortality to drought in carbon-dense tropical forests remains fraught with uncertainty, while extreme droughts are predicted to be more frequent and intense. Here, the potential of temporal autocorrelation of high-frequency variability in Landsat Enhanced Vegetation Index (EVI), an indicator of ecosystem resilience, to predict spatial and temporal variations of forest biomass mortality is evaluated against in situ census observations for 64 site-year combinations in Costa Rican tropical dry forests during the 2015 ENSO drought. Temporal autocorrelation, within the optimal moving window of 24 months, demonstrated robust predictive power for in situ mortality (leave-one-out cross-validation R2  = 0.54), which allows for estimates of annual biomass mortality patterns at 30 m resolution. Subsequent spatial analysis showed substantial fine-scale heterogeneity of forest mortality patterns, largely driven by drought intensity and ecosystem properties related to plant water use such as forest deciduousness and topography. Highly deciduous forest patches demonstrated much lower mortality sensitivity to drought stress than less deciduous forest patches after elevation was controlled. Our results highlight the potential of high-resolution remote sensing to "fingerprint" forest mortality and the significant role of ecosystem heterogeneity in forest biomass resistance to drought.

Beyond leaf habit: generalities in plant function across 97 tropical dry forest tree species.

Leaf habit has been hypothesized to define a linkage between the slow-fast plant economic spectrum and the drought resistance-avoidance trade-off in tropical forests ('slow-safe vs fast-risky'). However, variation in hydraulic traits as a function of leaf habit has rarely been explored for a large number of species. We sampled leaf and branch functional traits of 97 tropical dry forest tree species from four sites to investigate whether patterns of trait variation varied consistently in relation to leaf habit along the 'slow-safe vs fast-risky' trade-off. Leaf habit explained from 0% to 43.69% of individual trait variation. We found that evergreen and semi-deciduous species differed in their location along the multivariate trait ordination when compared to deciduous species. While deciduous species showed consistent trait values, evergreen species trait values varied as a function of the site. Last, trait values varied in relation to the proportion of deciduous species in the plant community. We found that leaf habit describes the strategies that define drought avoidance and plant economics in tropical trees. However, leaf habit alone does not explain patterns of trait variation, which suggests quantifying site-specific or species-specific uncertainty in trait variation as the way forward.

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.

Biomass increases attributed to both faster tree growth and altered allometric relationships under long-term carbon dioxide enrichment at a temperate forest.

Increases in atmospheric carbon dioxide (CO2 ) concentrations are expected to lead to increases in the rate of tree biomass accumulation, at least temporarily. On the one hand, trees may simply grow faster under higher CO2 concentrations, preserving the allometric relations that prevailed under lower CO2 concentrations. Alternatively, the allometric relations themselves may change. In this study, the effects of elevated CO2 (eCO2 ) on tree biomass and allometric relations were jointly assessed. Over 100 trees, grown at Duke Forest, NC, USA, were harvested from eight plots. Half of the plots had been subjected to CO2 enrichment from 1996 to 2010. Several subplots had also been subjected to nitrogen fertilization from 2005 to 2010. Allometric equations were developed to predict tree height, stem volume, and aboveground biomass components for loblolly pine (Pinus taeda L.), the dominant tree species, and broad-leaved species. Using the same diameter-based allometric equations for biomass, it was estimated that plots with eCO2 contained 21% more aboveground biomass, consistent with previous studies. However, eCO2 significantly affected allometry, and these changes had an additional effect on biomass. In particular, P. taeda trees at a given diameter were observed to be taller under eCO2 than under ambient CO2 due to changes in both the allometric scaling exponent and intercept. Accounting for allometric change increased the treatment effect of eCO2 on aboveground biomass from a 21% to a 27% increase. No allometric changes for the nondominant broad-leaved species were identified, nor were allometric changes associated with nitrogen fertilization. For P. taeda, it is concluded that eCO2 affects allometries, and that knowledge of allometry changes is necessary to accurately compute biomass under eCO2 . Further observations are needed to determine whether this assessment holds for other taxa.

A catastrophic tropical drought kills hydraulically vulnerable tree species.

Drought-related tree mortality is now a widespread phenomenon predicted to increase in magnitude with climate change. However, the patterns of which species and trees are most vulnerable to drought, and the underlying mechanisms have remained elusive, in part due to the lack of relevant data and difficulty of predicting the location of catastrophic drought years in advance. We used long-term demographic records and extensive databases of functional traits and distribution patterns to understand the responses of 20-53 species to an extreme drought in a seasonally dry tropical forest in Costa Rica, which occurred during the 2015 El Niño Southern Oscillation event. Overall, species-specific mortality rates during the drought ranged from 0% to 34%, and varied little as a function of tree size. By contrast, hydraulic safety margins correlated well with probability of mortality among species, while morphological or leaf economics spectrum traits did not. This firmly suggests hydraulic traits as targets for future research.

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

Censuses of tropical forest plots reveal large variation in biomass and plant composition. This paper evaluates whether such variation can emerge solely from realistic variation in a set of commonly measured soil chemical and physical properties. Controlled simulations were performed using a mechanistic model that includes forest dynamics, microbe-mediated biogeochemistry, and competition for nitrogen and phosphorus. Observations from 18 forest inventory plots in Guanacaste, Costa Rica were used to determine realistic variation in soil properties. In simulations of secondary succession, the across-plot range in plant biomass reached 30% of the mean and was attributable primarily to nutrient limitation and secondarily to soil texture differences that affected water availability. The contributions of different plant functional types to total biomass varied widely across plots and depended on soil nutrient status. In Central America, soil-induced variation in plant biomass increased with mean annual precipitation because of changes in nutrient limitation. In Central America, large variation in plant biomass and ecosystem composition arises mechanistically from realistic variation in soil properties. The degree of biomass and compositional variation is climate sensitive. In general, model predictions can be improved through better representation of soil nutrient processes, including their spatial variation.

Differential declines in Alaskan boreal forest vitality related to climate and competition.

Rapid warming and changes in water availability at high latitudes alter resource abundance, tree competition, and disturbance regimes. While these changes are expected to disrupt the functioning of boreal forests, their ultimate implications for forest composition are uncertain. In particular, recent site-level studies of the Alaskan boreal forest have reported both increases and decreases in productivity over the past few decades. Here, we test the idea that variations in Alaskan forest growth and mortality rates are contingent on species composition. Using forest inventory measurements and climate data from plots located throughout interior and south-central Alaska, we show significant growth and mortality responses associated with competition, midsummer vapor pressure deficit, and increased growing season length. The governing climate and competition processes differed substantially across species. Surprisingly, the most dramatic climate response occurred in the drought tolerant angiosperm species, trembling aspen, and linked high midsummer vapor pressure deficits to decreased growth and increased insect-related mortality. Given that species composition in the Alaskan and western Canadian boreal forests is projected to shift toward early-successional angiosperm species due to fire regime, these results underscore the potential for a reduction in boreal productivity stemming from increases in midsummer evaporative demand.

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

Numerous current efforts seek to improve the representation of ecosystem ecology and vegetation demographic processes within Earth System Models (ESMs). These developments are widely viewed as an important step in developing greater realism in predictions of future ecosystem states and fluxes. Increased realism, however, leads to increased model complexity, with new features raising a suite of ecological questions that require empirical constraints. Here, we review the developments that permit the representation of plant demographics in ESMs, and identify issues raised by these developments that highlight important gaps in ecological understanding. These issues inevitably translate into uncertainty in model projections but also allow models to be applied to new processes and questions concerning the dynamics of real-world ecosystems. We argue that stronger and more innovative connections to data, across the range of scales considered, are required to address these gaps in understanding. The development of first-generation land surface models as a unifying framework for ecophysiological understanding stimulated much research into plant physiological traits and gas exchange. Constraining predictions at ecologically relevant spatial and temporal scales will require a similar investment of effort and intensified inter-disciplinary communication.

Relation between rainfall intensity and savanna tree abundance explained by water use strategies.

Tree abundance in tropical savannas exhibits large and unexplained spatial variability. Here, we propose that differentiated tree and grass water use strategies can explain the observed negative relation between maximum tree abundance and rainfall intensity (defined as the characteristic rainfall depth on rainy days), and we present a biophysical tree-grass competition model to test this idea. The model is founded on a premise that has been well established in empirical studies, namely, that the relative growth rate of grasses is much higher compared with trees in wet conditions but that grasses are more susceptible to water stress and lose biomass more quickly in dry conditions. The model is coupled with a stochastic rainfall generator and then calibrated and tested using field observations from several African savanna sites. We show that the observed negative relation between maximum tree abundance and rainfall intensity can be explained only when differentiated water use strategies are accounted for. Numerical experiments reveal that this effect is more significant than the effect of root niche separation. Our results emphasize the importance of vegetation physiology in determining the responses of tree abundance to climate variations in tropical savannas and suggest that projected increases in rainfall intensity may lead to an increase in grass in this biome.

Variations of leaf longevity in tropical moist forests predicted by a trait-driven carbon optimality model.

Leaf longevity (LL) varies more than 20-fold in tropical evergreen forests, but it remains unclear how to capture these variations using predictive models. Current theories of LL that are based on carbon optimisation principles are challenging to quantitatively assess because of uncertainty across species in the 'ageing rate:' the rate at which leaf photosynthetic capacity declines with age. Here, we present a meta-analysis of 49 species across temperate and tropical biomes, demonstrating that the ageing rate of photosynthetic capacity is positively correlated with the mass-based carboxylation rate of mature leaves. We assess an improved trait-driven carbon optimality model with in situLL data for 105 species in two Panamanian forests. We show that our model explains over 40% of the cross-species variation in LL under contrasting light environment. Collectively, our results reveal how variation in LL emerges from carbon optimisation constrained by both leaf structural traits and abiotic environment.

Diversity in plant hydraulic traits explains seasonal and inter-annual variations of vegetation dynamics in seasonally dry tropical forests.

We assessed whether diversity in plant hydraulic traits can explain the observed diversity in plant responses to water stress in seasonally dry tropical forests (SDTFs). The Ecosystem Demography model 2 (ED2) was updated with a trait-driven mechanistic plant hydraulic module, as well as novel drought-phenology and plant water stress schemes. Four plant functional types were parameterized on the basis of meta-analysis of plant hydraulic traits. Simulations from both the original and the updated ED2 were evaluated against 5 yr of field data from a Costa Rican SDTF site and remote-sensing data over Central America. The updated model generated realistic plant hydraulic dynamics, such as leaf water potential and stem sap flow. Compared with the original ED2, predictions from our novel trait-driven model matched better with observed growth, phenology and their variations among functional groups. Most notably, the original ED2 produced unrealistically small leaf area index (LAI) and underestimated cumulative leaf litter. Both of these biases were corrected by the updated model. The updated model was also better able to simulate spatial patterns of LAI dynamics in Central America. Plant hydraulic traits are intercorrelated in SDTFs. Mechanistic incorporation of plant hydraulic traits is necessary for the simulation of spatiotemporal patterns of vegetation dynamics in SDTFs in vegetation models.

Dynamically downscaling predictions for deciduous tree leaf emergence in California under current and future climate.

Models that predict the timing of deciduous tree leaf emergence are typically very sensitive to temperature. However, many temperature data products, including those from climate models, have been developed at a very coarse spatial resolution. Such coarse-resolution temperature products can lead to highly biased predictions of leaf emergence. This study investigates how dynamical downscaling of climate models impacts simulations of deciduous tree leaf emergence in California. Models for leaf emergence are forced with temperatures simulated by a general circulation model (GCM) at ~200-km resolution for 1981-2000 and 2031-2050 conditions. GCM simulations are then dynamically downscaled to 32- and 8-km resolution, and leaf emergence is again simulated. For 1981-2000, the regional average leaf emergence date is 30.8 days earlier in 32-km simulations than in ~200-km simulations. Differences between the 32 and 8 km simulations are small and mostly local. The impact of downscaling from 200 to 8 km is ~15 % smaller in 2031-2050 than in 1981-2000, indicating that the impacts of downscaling are unlikely to be stationary.

Responses of terrestrial ecosystems and carbon budgets to current and future environmental variability.

We assess the significance of high-frequency variability of environmental parameters (sunlight, precipitation, temperature) for the structure and function of terrestrial ecosystems under current and future climate. We examine the influence of hourly, daily, and monthly variance using the Ecosystem Demography model version 2 in conjunction with the long-term record of carbon fluxes measured at Harvard Forest. We find that fluctuations of sunlight and precipitation are strongly and nonlinearly coupled to ecosystem function, with effects that accumulate through annual and decadal timescales. Increasing variability in sunlight and precipitation leads to lower rates of carbon sequestration and favors broad-leaved deciduous trees over conifers. Temperature variability has only minor impacts by comparison. We also find that projected changes in sunlight and precipitation variability have important implications for carbon storage and ecosystem structure and composition. Based on Intergovernmental Panel on Climate Change model estimates for changes in high-frequency meteorological variability over the next 100 years, we expect that terrestrial ecosystems will be affected by changes in variability almost as much as by changes in mean climate. We conclude that terrestrial ecosystems are highly sensitive to high-frequency meteorological variability, and that accurate knowledge of the statistics of this variability is essential for realistic predictions of ecosystem structure and functioning.

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.

Tropical carbon sink accelerated by symbiotic dinitrogen fixation.

A major uncertainty in the land carbon cycle is whether symbiotic nitrogen fixation acts to enhance the tropical forest carbon sink. Nitrogen-fixing trees can supply vital quantities of the growth-limiting nutrient nitrogen, but the extent to which the resulting carbon-nitrogen feedback safeguards ecosystem carbon sequestration remains unclear. We combine (i) field observations from 112 plots spanning 300 years of succession in Panamanian tropical forests, and (ii) a new model that resolves nitrogen and light competition at the scale of individual trees. Fixation doubled carbon accumulation in early succession and enhanced total carbon in mature forests by ~10% (~12MgC ha-1) through two mechanisms: (i) a direct fixation effect on tree growth, and (ii) an indirect effect on the successional sequence of non-fixing trees. We estimate that including nitrogen-fixing trees in Neotropical reforestation projects could safeguard the sequestration of 6.7 Gt CO2 over the next 20 years. Our results highlight the connection between functional diversity of plant communities and the critical ecosystem service of carbon sequestration for mitigating climate change.

Accelerating rates of Arctic carbon cycling revealed by long-term atmospheric CO2 measurements.

The contemporary Arctic carbon balance is uncertain, and the potential for a permafrost carbon feedback of anywhere from 50 to 200 petagrams of carbon (Schuur et al., 2015) compromises accurate 21st-century global climate system projections. The 42-year record of atmospheric CO2 measurements at Barrow, Alaska (71.29 N, 156.79 W), reveals significant trends in regional land-surface CO2 anomalies (ΔCO2), indicating long-term changes in seasonal carbon uptake and respiration. Using a carbon balance model constrained by ΔCO2, we find a 13.4% decrease in mean carbon residence time (50% confidence range = 9.2 to 17.6%) in North Slope tundra ecosystems during the past four decades, suggesting a transition toward a boreal carbon cycling regime. Temperature dependencies of respiration and carbon uptake suggest that increases in cold season Arctic labile carbon release will likely continue to exceed increases in net growing season carbon uptake under continued warming trends.

Predicting ecosystem dynamics at regional scales: an evaluation of a terrestrial biosphere model for the forests of northeastern North America.

Terrestrial biosphere models are important tools for diagnosing both the current state of the terrestrial carbon cycle and forecasting terrestrial ecosystem responses to global change. While there are a number of ongoing assessments of the short-term predictive capabilities of terrestrial biosphere models using flux-tower measurements, to date there have been relatively few assessments of their ability to predict longer term, decadal-scale biomass dynamics. Here, we present the results of a regional-scale evaluation of the Ecosystem Demography version 2 (ED2)-structured terrestrial biosphere model, evaluating the model's predictions against forest inventory measurements for the northeast USA and Quebec from 1985 to 1995. Simulations were conducted using a default parametrization, which used parameter values from the literature, and a constrained model parametrization, which had been developed by constraining the model's predictions against 2 years of measurements from a single site, Harvard Forest (42.5° N, 72.1° W). The analysis shows that the constrained model parametrization offered marked improvements over the default model formulation, capturing large-scale variation in patterns of biomass dynamics despite marked differences in climate forcing, land-use history and species-composition across the region. These results imply that data-constrained parametrizations of structured biosphere models such as ED2 can be successfully used for regional-scale ecosystem prediction and forecasting. We also assess the model's ability to capture sub-grid scale heterogeneity in the dynamics of biomass growth and mortality of different sizes and types of trees, and then discuss the implications of these analyses for further reducing the remaining biases in the model's predictions.