Latitudinal patterns in stabilizing density dependence of forest communities.

Numerous studies have shown reduced performance in plants that are surrounded by neighbours of the same species1,2, a phenomenon known as conspecific negative density dependence (CNDD)3. A long-held ecological hypothesis posits that CNDD is more pronounced in tropical than in temperate forests4,5, which increases community stabilization, species coexistence and the diversity of local tree species6,7. Previous analyses supporting such a latitudinal gradient in CNDD8,9 have suffered from methodological limitations related to the use of static data10-12. Here we present a comprehensive assessment of latitudinal CNDD patterns using dynamic mortality data to estimate species-site-specific CNDD across 23 sites. Averaged across species, we found that stabilizing CNDD was present at all except one site, but that average stabilizing CNDD was not stronger toward the tropics. However, in tropical tree communities, rare and intermediate abundant species experienced stronger stabilizing CNDD than did common species. This pattern was absent in temperate forests, which suggests that CNDD influences species abundances more strongly in tropical forests than it does in temperate ones13. We also found that interspecific variation in CNDD, which might attenuate its stabilizing effect on species diversity14,15, was high but not significantly different across latitudes. Although the consequences of these patterns for latitudinal diversity gradients are difficult to evaluate, we speculate that a more effective regulation of population abundances could translate into greater stabilization of tropical tree communities and thus contribute to the high local diversity of tropical forests.

Leaf gene expression trajectories during the growing season are consistent between sites and years in American beech.

Transcriptomics provides a versatile tool for ecological monitoring. Here, through genome-guided profiling of transcripts mapping to 33 042 gene models, expression differences can be discerned among multi-year and seasonal leaf samples collected from American beech trees at two latitudinally separated sites. Despite a bottleneck due to post-Columbian deforestation, the single nucleotide polymorphism-based population genetic background analysis has yielded sufficient variation to account for differences between populations and among individuals. Our expression analyses during spring-summer and summer-autumn transitions for two consecutive years involved 4197 differentially expressed protein coding genes. Using Populus orthologues we reconstructed a protein-protein interactome representing leaf physiological states of trees during the seasonal transitions. Gene set enrichment analysis revealed gene ontology terms that highlight molecular functions and biological processes possibly influenced by abiotic forcings such as recovery from drought and response to excess precipitation. Further, based on 324 co-regulated transcripts, we focused on a subset of GO terms that could be putatively attributed to late spring phenological shifts. Our conservative results indicate that extended transcriptome-based monitoring of forests can capture diverse ranges of responses including air quality, chronic disease, as well as herbivore outbreaks that require activation and/or downregulation of genes collectively tuning reaction norms maintaining the survival of long living trees such as the American beech.

Toward a forest biomass reference measurement system for remote sensing applications.

Forests contribute to climate change mitigation through carbon storage and uptake, but the extent to which this carbon pool varies in space and time is still poorly known. Several Earth Observation missions have been specifically designed to address this issue, for example, NASA's GEDI, NASA-ISRO's NISAR and ESA's BIOMASS. Yet, all these missions' products require independent and consistent validation. A permanent, global, in situ, site-based forest biomass reference measurement system relying on ground data of the highest possible quality is therefore needed. Here, we have assembled a list of almost 200 high-quality sites through an in-depth review of the literature and expert knowledge. In this study, we explore how representative these sites are in terms of their coverage of environmental conditions, geographical space and biomass-related forest structure, compared to those experienced by forests worldwide. This work also aims at identifying which sites are the most representative, and where to invest to improve the representativeness of the proposed system. We show that the environmental coverage of the system does not seem to improve after at least the 175 most representative sites are included, but geographical and structural coverages continue to improve as more sites are added. We highlight the areas of poor environmental, geographical, or structural coverage, including, but not limited to, Canada, the western half of the USA, Mexico, Patagonia, Angola, Zambia, eastern Russia, and tropical and subtropical highlands (e.g. in Colombia, the Himalayas, Borneo, Papua). For the proposed system to succeed, we stress that (1) data must be collected and processed applying the same standards across all countries and continents; (2) system establishment and management must be inclusive and equitable, with careful consideration of working conditions; and (3) training and site partner involvement in downstream activities should be mandatory.

Damage to living trees contributes to almost half of the biomass losses in tropical forests.

Accurate estimates of forest biomass stocks and fluxes are needed to quantify global carbon budgets and assess the response of forests to climate change. However, most forest inventories consider tree mortality as the only aboveground biomass (AGB) loss without accounting for losses via damage to living trees: branchfall, trunk breakage, and wood decay. Here, we use ~151,000 annual records of tree survival and structural completeness to compare AGB loss via damage to living trees to total AGB loss (mortality + damage) in seven tropical forests widely distributed across environmental conditions. We find that 42% (3.62 Mg ha-1 year-1 ; 95% confidence interval [CI] 2.36-5.25) of total AGB loss (8.72 Mg ha-1 year-1 ; CI 5.57-12.86) is due to damage to living trees. Total AGB loss was highly variable among forests, but these differences were mainly caused by site variability in damage-related AGB losses rather than by mortality-related AGB losses. We show that conventional forest inventories overestimate stand-level AGB stocks by 4% (1%-17% range across forests) because assume structurally complete trees, underestimate total AGB loss by 29% (6%-57% range across forests) due to overlooked damage-related AGB losses, and overestimate AGB loss via mortality by 22% (7%-80% range across forests) because of the assumption that trees are undamaged before dying. Our results indicate that forest carbon fluxes are higher than previously thought. Damage on living trees is an underappreciated component of the forest carbon cycle that is likely to become even more important as the frequency and severity of forest disturbances increase.

Individual tree damage dominates mortality risk factors across six tropical forests.

The relative importance of tree mortality risk factors remains unknown, especially in diverse tropical forests where species may vary widely in their responses to particular conditions. We present a new framework for quantifying the importance of mortality risk factors and apply it to compare 19 risks on 31 203 trees (1977 species) in 14 one-year periods in six tropical forests. We defined a condition as a risk factor for a species if it was associated with at least a doubling of mortality rate in univariate analyses. For each risk, we estimated prevalence (frequency), lethality (difference in mortality between trees with and without the risk) and impact ('excess mortality' associated with the risk, relative to stand-level mortality). The most impactful risk factors were light limitation and crown/trunk loss; the most prevalent were light limitation and small size; the most lethal were leaf damage and wounds. Modes of death (standing, broken and uprooted) had limited links with previous conditions and mortality risk factors. We provide the first ranking of importance of tree-level mortality risk factors in tropical forests. Future research should focus on the links between these risks, their climatic drivers and the physiological processes to enable mechanistic predictions of future tree mortality.

Distribution of biomass dynamics in relation to tree size in forests across the world.

Tree size shapes forest carbon dynamics and determines how trees interact with their environment, including a changing climate. Here, we conduct the first global analysis of among-site differences in how aboveground biomass stocks and fluxes are distributed with tree size. We analyzed repeat tree censuses from 25 large-scale (4-52 ha) forest plots spanning a broad climatic range over five continents to characterize how aboveground biomass, woody productivity, and woody mortality vary with tree diameter. We examined how the median, dispersion, and skewness of these size-related distributions vary with mean annual temperature and precipitation. In warmer forests, aboveground biomass, woody productivity, and woody mortality were more broadly distributed with respect to tree size. In warmer and wetter forests, aboveground biomass and woody productivity were more right skewed, with a long tail towards large trees. Small trees (1-10 cm diameter) contributed more to productivity and mortality than to biomass, highlighting the importance of including these trees in analyses of forest dynamics. Our findings provide an improved characterization of climate-driven forest differences in the size structure of aboveground biomass and dynamics of that biomass, as well as refined benchmarks for capturing climate influences in vegetation demographic models.

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.

Joint effects of climate, tree size, and year on annual tree growth derived from tree-ring records of ten globally distributed forests.

Tree rings provide an invaluable long-term record for understanding how climate and other drivers shape tree growth and forest productivity. However, conventional tree-ring analysis methods were not designed to simultaneously test effects of climate, tree size, and other drivers on individual growth. This has limited the potential to test ecologically relevant hypotheses on tree growth sensitivity to environmental drivers and their interactions with tree size. Here, we develop and apply a new method to simultaneously model nonlinear effects of primary climate drivers, reconstructed tree diameter at breast height (DBH), and calendar year in generalized least squares models that account for the temporal autocorrelation inherent to each individual tree's growth. We analyze data from 3811 trees representing 40 species at 10 globally distributed sites, showing that precipitation, temperature, DBH, and calendar year have additively, and often interactively, influenced annual growth over the past 120 years. Growth responses were predominantly positive to precipitation (usually over ≥3-month seasonal windows) and negative to temperature (usually maximum temperature, over ≤3-month seasonal windows), with concave-down responses in 63% of relationships. Climate sensitivity commonly varied with DBH (45% of cases tested), with larger trees usually more sensitive. Trends in ring width at small DBH were linked to the light environment under which trees established, but basal area or biomass increments consistently reached maxima at intermediate DBH. Accounting for climate and DBH, growth rate declined over time for 92% of species in secondary or disturbed stands, whereas growth trends were mixed in older forests. These trends were largely attributable to stand dynamics as cohorts and stands age, which remain challenging to disentangle from global change drivers. By providing a parsimonious approach for characterizing multiple interacting drivers of tree growth, our method reveals a more complete picture of the factors influencing growth than has previously been possible.

Demographic composition, not demographic diversity, predicts biomass and turnover across temperate and tropical forests.

The growth and survival of individual trees determine the physical structure of a forest with important consequences for forest function. However, given the diversity of tree species and forest biomes, quantifying the multitude of demographic strategies within and across forests and the way that they translate into forest structure and function remains a significant challenge. Here, we quantify the demographic rates of 1961 tree species from temperate and tropical forests and evaluate how demographic diversity (DD) and demographic composition (DC) differ across forests, and how these differences in demography relate to species richness, aboveground biomass (AGB), and carbon residence time. We find wide variation in DD and DC across forest plots, patterns that are not explained by species richness or climate variables alone. There is no evidence that DD has an effect on either AGB or carbon residence time. Rather, the DC of forests, specifically the relative abundance of large statured species, predicted both biomass and carbon residence time. Our results demonstrate the distinct DCs of globally distributed forests, reflecting biogeography, recent history, and current plot conditions. Linking the DC of forests to resilience or vulnerability to climate change, will improve the precision and accuracy of predictions of future forest composition, structure, and function.

Interactions between all pairs of neighboring trees in 16 forests worldwide reveal details of unique ecological processes in each forest, and provide windows into their evolutionary histories.

When Darwin visited the Galapagos archipelago, he observed that, in spite of the islands' physical similarity, members of species that had dispersed to them recently were beginning to diverge from each other. He postulated that these divergences must have resulted primarily from interactions with sets of other species that had also diverged across these otherwise similar islands. By extrapolation, if Darwin is correct, such complex interactions must be driving species divergences across all ecosystems. However, many current general ecological theories that predict observed distributions of species in ecosystems do not take the details of between-species interactions into account. Here we quantify, in sixteen forest diversity plots (FDPs) worldwide, highly significant negative density-dependent (NDD) components of both conspecific and heterospecific between-tree interactions that affect the trees' distributions, growth, recruitment, and mortality. These interactions decline smoothly in significance with increasing physical distance between trees. They also tend to decline in significance with increasing phylogenetic distance between the trees, but each FDP exhibits its own unique pattern of exceptions to this overall decline. Unique patterns of between-species interactions in ecosystems, of the general type that Darwin postulated, are likely to have contributed to the exceptions. We test the power of our null-model method by using a deliberately modified data set, and show that the method easily identifies the modifications. We examine how some of the exceptions, at the Wind River (USA) FDP, reveal new details of a known allelopathic effect of one of the Wind River gymnosperm species. Finally, we explore how similar analyses can be used to investigate details of many types of interactions in these complex ecosystems, and can provide clues to the evolution of these interactions.

Wildlife disturbances as a source of conspecific negative density-dependent mortality in tropical trees.

Large vertebrates are rarely considered important drivers of conspecific negative density-dependent mortality (CNDD) in plants because they are generalist consumers. However, disturbances like trampling and nesting also cause plant mortality, and their impact on plant diversity depends on the spatial overlap between wildlife habitat preferences and plant species composition. We studied the impact of native wildlife on a hyperdiverse tree community in Malaysia. Pigs (Sus scrofa) are abnormally abundant at the site due to food subsidies in nearby farmland and they construct birthing nests using hundreds of tree saplings. We tagged 34 950 tree saplings in a 25 ha plot during an initial census and assessed the source mortality by recovering tree tags from pig nests (n = 1672 pig-induced deaths). At the stand scale, pigs nested in flat dry habitats, and at the local neighbourhood scale, they nested within clumps of saplings, both of which are intuitive for safe and efficient nest building. At the stand scale, flat dry habitats contained higher sapling densities and higher proportions of common species, so pig nesting increased the weighted average species evenness across habitats. At the neighbourhood scale, pig-induced sapling mortality was associated with higher heterospecific and especially conspecific sapling densities. Tree species have clumped distributions due to dispersal limitation and habitat filtering, so pig disturbances in sapling clumps indirectly caused CNDD. As a result, Pielou species evenness in 400 m2 quadrats increased 105% more in areas with pig-induced deaths than areas without disturbances. Wildlife induced CNDD and this supported tree species evenness, but they also drove a 62% decline in sapling densities from 1996 to 2010, which is unsustainable. We suspect pig nesting is an important feature shaping tree composition throughout the region.

Leaf turgor loss point shapes local and regional distributions of evergreen but not deciduous tropical trees.

The effects of climate change on tropical forests will depend on how diverse tropical tree species respond to drought. Current distributions of evergreen and deciduous tree species across local and regional moisture gradients reflect their ability to tolerate drought stress, and might be explained by functional traits. We measured leaf water potential at turgor loss (i.e. 'wilting point'; πtlp ), wood density (WD) and leaf mass per area (LMA) on 50 of the most abundant tree species in central Panama. We then tested their ability to explain distributions of evergreen and deciduous species within a 50 ha plot on Barro Colorado Island and across a 70 km rainfall gradient spanning the Isthmus of Panama. Among evergreen trees, species with lower πtlp were associated with drier habitats, with πtlp explaining 28% and 32% of habitat association on local and regional scales, respectively, greatly exceeding the predictive power of WD and LMA. In contrast, πtlp did not predict habitat associations among deciduous species. Across spatial scales, πtlp is a useful indicator of habitat preference for tropical tree species that retain their leaves during periods of water stress, and holds the potential to predict vegetation responses to climate change.

Managing access to advanced therapy medicinal products: Challenges for NHS Wales.

Advanced Therapy Medicinal Products (ATMPs), which include gene, somatic cell therapies and tissue-engineered medicines, have the potential to transform current care pathways by offering durable and potentially curative outcomes. However, they are exceptionally expensive, with prices exceeding £1m per patient in some cases. With an expectation that a large number of ATMPs will soon gain marketing authorisation (global market is estimated to reach £9bn to £14bn by 2025), healthcare payers and providers face a number of challenges to facilitate patient access to this new category of medicines. This viewpoint reflects on the experience of introducing ATMPs into the National Health Service in Wales where £1 in every £200 spent on medicines (2019/2020) is expected to be on ATMPs for just 20 patients. Evidence to date makes it apparent that decisions regarding clinical and cost-effectiveness and the scale of the budget impact of implementing ATMPs create both financial and health service risks. Consequently, there are significant policy implications. A critical examination is made of the approaches taken for the health technology assessment and appraisal of ATMPs, the methods of payment and service impacts of these medicines, and the approach taken to horizon scanning and subsequent modelling of the financial impact over the next 10 years.

Temporal population variability in local forest communities has mixed effects on tree species richness across a latitudinal gradient.

Among the local processes that determine species diversity in ecological communities, fluctuation-dependent mechanisms that are mediated by temporal variability in the abundances of species populations have received significant attention. Higher temporal variability in the abundances of species populations can increase the strength of temporal niche partitioning but can also increase the risk of species extinctions, such that the net effect on species coexistence is not clear. We quantified this temporal population variability for tree species in 21 large forest plots and found much greater variability for higher latitude plots with fewer tree species. A fitted mechanistic model showed that among the forest plots, the net effect of temporal population variability on tree species coexistence was usually negative, but sometimes positive or negligible. Therefore, our results suggest that temporal variability in the abundances of species populations has no clear negative or positive contribution to the latitudinal gradient in tree species richness.

Lithological constraints on resource economies shape the mycorrhizal composition of a Bornean rain forest.

Arbuscular mycorrhizal fungi (AMF) and ectomycorrhizal fungi (EMF) produce contrasting plant-soil feedbacks, but how these feedbacks are constrained by lithology is poorly understood. We investigated the hypothesis that lithological drivers of soil fertility filter plant resource economic strategies in ways that influence the relative fitness of trees with AMF or EMF symbioses in a Bornean rain forest containing species with both mycorrhizal strategies. Using forest inventory data on 1245 tree species, we found that although AMF-hosting trees had greater relative dominance on all soil types, with declining lithological soil fertility EMF-hosting trees became more dominant. Data on 13 leaf traits and wood density for a total of 150 species showed that variation was almost always associated with soil type, whereas for six leaf traits (structural properties; carbon, nitrogen, phosphorus ratios, nitrogen isotopes), variation was also associated with mycorrhizal strategy. EMF-hosting species had slower leaf economics than AMF-hosts, demonstrating the central role of mycorrhizal symbiosis in plant resource economies. At the global scale, climate has been shown to shape forest mycorrhizal composition, but here we show that in communities it depends on soil lithology, suggesting scale-dependent abiotic factors influence feedbacks underlying the relative fitness of different mycorrhizal strategies.

Testing for changes in biomass dynamics in large-scale forest datasets.

Tropical forest responses to climate and atmospheric change are critical to the future of the global carbon budget. Recent studies have reported increases in estimated above-ground biomass (EAGB) stocks, productivity, and mortality in old-growth tropical forests. These increases could reflect a shift in forest functioning due to global change and/or long-lasting recovery from past disturbance. We introduce a novel approach to disentangle the relative contributions of these mechanisms by decomposing changes in whole-plot biomass fluxes into contributions from changes in the distribution of gap-successional stages and changes in fluxes for a given stage. Using 30 years of forest dynamic data at Barro Colorado Island, Panama, we investigated temporal variation in EAGB fluxes as a function of initial EAGB (EAGBi ) in 10 × 10 m quadrats. Productivity and mortality fluxes both increased strongly with initial quadrat EAGB. The distribution of EAGB (and thus EAGBi ) across quadrats hardly varied over 30 years (and seven censuses). EAGB fluxes as a function of EAGBi varied largely and significantly among census intervals, with notably higher productivity in 1985-1990 associated with recovery from the 1982-1983 El Niño event. Variation in whole-plot fluxes among census intervals was explained overwhelmingly by variation in fluxes as a function of EAGBi , with essentially no contribution from changes in EAGBi distributions. The high observed temporal variation in productivity and mortality suggests that this forest is very sensitive to climate variability. There was no consistent long-term trend in productivity, mortality, or biomass in this forest over 30 years, although the temporal variability in productivity and mortality was so strong that it could well mask a substantial trend. Accurate prediction of future tropical forest carbon budgets will require accounting for disturbance-recovery dynamics and understanding temporal variability in productivity and mortality.

Crown damage and the mortality of tropical trees.

What causes individual tree death in tropical forests remains a major gap in our understanding of the biology of tropical trees and leads to significant uncertainty in predicting global carbon cycle dynamics. We measured individual characteristics (diameter at breast height, wood density, growth rate, crown illumination and crown form) and environmental conditions (soil fertility and habitat suitability) for 26 425 trees ≥ 10 cm diameter at breast height belonging to 416 species in a 52-ha plot in Lambir Hills National Park, Malaysia. We used structural equation models to investigate the relationships among the different factors and tree mortality. Crown form (a proxy for mechanical damage and other stresses) and prior growth were the two most important factors related to mortality. The effect of all variables on mortality (except habitat suitability) was substantially greater than expected by chance. Tree death is the result of interactions between factors, including direct and indirect effects. Crown form/damage and prior growth mediated most of the effect of tree size, wood density, fertility and habitat suitability on mortality. Large-scale assessment of crown form or status may result in improved prediction of individual tree death at the landscape scale.

Drivers and mechanisms of tree mortality in moist tropical forests.

Tree mortality rates appear to be increasing in moist tropical forests (MTFs) with significant carbon cycle consequences. Here, we review the state of knowledge regarding MTF tree mortality, create a conceptual framework with testable hypotheses regarding the drivers, mechanisms and interactions that may underlie increasing MTF mortality rates, and identify the next steps for improved understanding and reduced prediction. Increasing mortality rates are associated with rising temperature and vapor pressure deficit, liana abundance, drought, wind events, fire and, possibly, CO2 fertilization-induced increases in stand thinning or acceleration of trees reaching larger, more vulnerable heights. The majority of these mortality drivers may kill trees in part through carbon starvation and hydraulic failure. The relative importance of each driver is unknown. High species diversity may buffer MTFs against large-scale mortality events, but recent and expected trends in mortality drivers give reason for concern regarding increasing mortality within MTFs. Models of tropical tree mortality are advancing the representation of hydraulics, carbon and demography, but require more empirical knowledge regarding the most common drivers and their subsequent mechanisms. We outline critical datasets and model developments required to test hypotheses regarding the underlying causes of increasing MTF mortality rates, and improve prediction of future mortality under climate change.

Soil drivers of local-scale tree growth in a lowland tropical forest.

Soil nutrients influence the distribution of tree species in lowland tropical forests, but their effect on productivity, especially at local scales, remains unclear. We used tree census, canopy occupancy, and soil data from the Barro Colorado Island (BCI; Panama) 50-ha forest dynamics plot to investigate the influence of soil nutrients and potential toxins on aboveground tree productivity. Growth was calculated as the increase in diameter of 150,000 individual stems ≥1 cm diameter at breast height, representing 207 species. The effects of soil variables and other strong predictors of growth (e.g., light) were estimated using hierarchical, linear, mixed-effects models. Growth was weakly positively associated with phosphorus (P), particularly for understory tree species that are typically considered to be limited by light. In contrast, growth was strongly negatively related to manganese (Mn) and aluminum (Al), although the latter effect was confounded by strong correlations between Al and other soil variables. The negative response to increasing Mn (and Al) suggests a toxicity effect due to solubilization and uptake of amorphous pools of metal oxides in the soil. These results show that P limits tropical tree growth at local scale on BCI, but that toxic metals represent an even greater constraint on productivity.

Drought-induced mortality patterns and rapid biomass recovery in a terra firme forest in the Colombian Amazon.

Extreme climatic events affecting the Amazon region are expected to become more frequent under ongoing climate change. In this study, we assessed the responses to the 2010 drought of over 14,000 trees ≥10 cm dbh in a 25 ha lowland forest plot in the Colombian Amazon and how these responses varied among topographically defined habitats, with tree size, and with species wood density. Tree mortality was significantly higher during the 2010-2013 period immediately after the drought than in 2007-2010. The post-drought increase in mortality was stronger for trees located in valleys (+243%) than for those located on slopes (+67%) and ridges (+57%). Tree-based generalized linear mixed models showed a significant negative effect of species wood density on mortality and no effect of tree size. Despite the elevated post-drought mortality, aboveground biomass increased from 2007 to 2013 by 1.62 Mg ha-1  yr-1 (95% CI 0.80-2.43 Mg ha-1  yr-1 ). Biomass change varied among habitats, with no significant increase on the slopes (1.05, 95% CI -0.76 to 2.85 Mg ha-1  yr-1 ), a significant increase in the valleys (1.33, 95% CI 0.37-2.34 Mg ha-1  yr-1 ), and a strong increase on the ridges (2.79, 95% CI 1.20-4.21 Mg ha-1  yr-1 ). These results indicate a high carbon resilience of this forest to the 2010 drought due to habitat-associated and interspecific heterogeneity in responses including directional changes in functional composition driven by enhanced performance of drought-tolerant species that inhabit the drier ridges.