Search

Climatic extreme events are expected to occur more frequently and potentially be stronger in the future, increasing the likelihood of unprecedented climate extremes (UCEs), or record-breaking events such as prolonged droughts, to occur. To prepare for UCEs and their impacts, we need to develop a better understanding of terrestrial ecosystem responses to events such as extreme drought. We know that intense, extreme droughts can substantially affect ecosystem stability and carbon cycling through increased plant mortality and delaying ecosystem recovery. Our ability to predict such effects is limited due to the lack of experiments focusing on climatic excursions beyond the range of historical experience.We explore the response of forest ecosystems to UCEs using two dynamic vegetation demographic models (VDMs), ED2 and LPJ-GUESS, in which the abundances of different plant functional types, as well as tree size- and age-class structure, are emergent properties of resource competition. We investigate the hypothesis that ecosystem responses to UCEs (e.g., unprecedented droughts) cannot be extrapolated from ecosystem responses to milder extremes, as a result of non-linear ecosystem responses (e.g. due to plant plasticity, functional diversity, and trait combinations). We evaluate each model’s mechanisms and state variables prior, during, and after a continuum of drought intensities ultimately reaching very extreme drought scenarios (i.e., 0% to 100% reduction in precipitation for drought durations of 1-year, 2-year, and 4-year scenarios) at two dry forested sites: Palo Verde, Costa Rica (i.e. tropical) and EucFACE, Australia (i.e. temperate). Both models produce nonlinear responses to these UCEs. Due to differences in model structure and process representation, the model’s sensitivity of biomass loss diverged based on either duration or intensity of droughts, as well as different model responses at each site. Biomass losses in ED2 are sensitive to drought duration, while in LPJ-GUESS they are mainly driven by drought intensity. Elevated atmospheric CO2 concentrations alone did not buffer the ecosystems from carbon losses during UCEs in the majority of our simulations. Our findings highlight discrepancies in process formulations and uncertainties in models, notably related to availability in plant carbohydrate storage and the diversity of plant hydraulic schemes. This shows that different hypotheses of plant responses to UCEs exist in two similar models, reflecting knowledge gaps, which should be tested with gap-informed field experiments. This iterative modeling-experiment framework would help improve predictions of terrestrial ecosystem responses and climate feedbacks.

Current climate change scenarios indicate warmer temperatures and the potential for more extreme droughts in the tropics, such that a mechanistic understanding of the water cycle from individual trees to landscapes is needed to adequately predict future changes in forest structure and function. In this study, we contrasted physiological responses of tropical trees during a normal dry season with the extreme dry season due to the 2015–2016 El Niño-Southern Oscillation (ENSO) event. We quantified high resolution temporal dynamics of sap velocity (Vs), stomatal conductance (gs) and leaf water potential (ΨL) of multiple canopy trees, and their correlations with leaf temperature (Tleaf) and environmental conditions [direct solar radiation, air temperature (Tair) and vapor pressure deficit (VPD)]. The experiment leveraged canopy access towers to measure adjacent trees at the ZF2 and Tapajós tropical forest research (near the cities of Manaus and Santarém). The temporal difference between the peak of gs (late morning) and the peak of VPD (early afternoon) is one of the major regulators of sap velocity hysteresis patterns. Sap velocity displayed species-specific diurnal hysteresis patterns reflected by changes in Tleaf. In the morning, Tleaf and sap velocity displayed a sigmoidal relationship. In the afternoon, stomatal conductance declined as Tleaf approached a daily peak, allowing ΨL to begin recovery, while sap velocity declined with an exponential relationship with Tleaf. In Manaus, hysteresis indices of the variables Tleaf-Tair and ΨL-Tleaf were calculated for different species and a significant difference (p < 0.01, α = 0.05) was observed when the 2015 dry season (ENSO period) was compared with the 2017 dry season (“control scenario”). In some days during the 2015 ENSO event, Tleaf approached 40°C for all studied species and the differences between Tleaf and Tair reached as high at 8°C (average difference: 1.65 ± 1.07°C). Generally, Tleaf was higher than Tair during the middle morning to early afternoon, and lower than Tair during the early morning, late afternoon and night. Our results support the hypothesis that partial stomatal closure allows for a recovery in ΨL during the afternoon period giving an observed counterclockwise hysteresis pattern between ΨL and Tleaf. © 2019 Gimenez, Jardine, Higuchi, Negrón-Juárez, Sampaio-Filho, Cobello, Fontes, Dawson, Varadharajan, Christianson, Spanner, Araújo, Warren, Newman, Holm, Koven, McDowell and Chambers.

Uncertainties surrounding vegetation response to increased disturbance rates associated with climate change remains a major global change issue for Amazon forests. Additionally, turnover rates computed as the average of mortality and recruitment rates in the Western Amazon basin are doubled when compared to the Central Amazon, and notable gradients currently exist in specific wood density and aboveground biomass (AGB) between these two regions. This study investigates the extent to which the variation in disturbance regimes contributes to these regional gradients. To address these issues, we evaluated disturbance-recovery processes under two scenarios of increased disturbance rates in a complex Central Amazon forest using first ZELIG-TROP, a dynamic vegetation gap model which we calibrated using long-term inventory data, and second using the Community Land Model (CLM), a global land surface model that is part of the Community Earth System Model (CESM). Upon doubling the mortality rate in the Central Amazon to mirror the natural disturbance regime in the Western Amazon of ∼2% mortality, at steady-state, AGB significantly decreased by 41.9% and there was no significant difference between the modeled AGB of 104 Mg C ha−1 and empirical AGB from the western Amazon datasets of 107 Mg C ha−1. We confirm that increases in natural disturbance rates in the Central Amazon will result in terrestrial carbon loss associated with higher turnover. However, different processes were responsible for the reductions in AGB between the models and empirical datasets. We observed that with increased turnover, the subsequent decrease in wood density drives the reduction in AGB in empirical datasets. However, decrease in stand basal area was the driver of the drop in AGB in ZELIG-TROP, and decreased leaf area index (LAI) was the driver in CLM. Further comparisons found that stem density, specific wood density, and basal area growth rates differed between the two Amazonian regions. This suggests that: (1) the variability between regions cannot be entirely explained by the variability in disturbance regime, but rather potentially sensitive to intrinsic environmental factors; or (2) the models are not accurately simulating all forest characteristics in response to increased disturbances. Last, to help quantify the impacts of increased disturbances on climate and the earth system, we evaluated the fidelity of tree mortality and disturbance in a global land surface model: CLM. For a 100% increase in annual mortality rate, both ZELIG-TROP and CLM were in close agreement with each other and predicted a net carbon loss of 41.9 and 49.9%, respectively, with an insignificant effect on aboveground net primary productivity (ANPP). Likewise, a 20% increase in mortality every 50 years (i.e. periodic disturbance treatment) resulted in a reciprocal biomass loss of 18.3 and 18.7% in ZELIG-TROP and CLM, respectively.

We describe an implementation of the Ecosystem Demography (ED) concept in the Community Land Model. The structure of CLM(ED) and the physiological and structural modifications applied to the CLM are presented. A major motivation of this development is to allow the prediction of biome boundaries directly from plant physiological traits via their competitive interactions. Here we investigate the performance of the model for an example biome boundary in eastern North America. We explore the sensitivity of the predicted biome boundaries and ecosystem properties to the variation of leaf properties using the parameter space defined by the GLOPNET global leaf trait database. Furthermore, we investigate the impact of four sequential alterations to the structural assumptions in the model governing the relative carbon economy of deciduous and evergreen plants. The default assumption is that the costs and benefits of deciduous vs. evergreen leaf strategies, in terms of carbon assimilation and expenditure, can reproduce the geographical structure of biome boundaries and ecosystem functioning. We find some support for this assumption, but only under particular combinations of model traits and structural assumptions. Many questions remain regarding the preferred methods for deployment of plant trait information in land surface models. In some cases, plant traits might best be closely linked to each other, but we also find support for direct linkages to environmental conditions. We advocate intensified study of the costs and benefits of plant life history strategies in different environments and the increased use of parametric and structural ensembles in the development and analysis of complex vegetation models.

No claim to original US government works New Phytologist © 2018 New Phytologist Trust 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.

Disturbance-induced tree mortality is a key factor regulating the carbon balance of a forest, but tree mortality and its subsequent effects are poorly represented processes in terrestrial ecosystem models. It is thus unclear whether models can robustly simulate moderate (non-catastrophic) disturbances, which tend to increase biological and structural complexity and are increasingly common in aging US forests. We tested whether three forest ecosystem models – Biome-BGC (BioGeochemical Cycles), a classic big-leaf model, and the ZELIG and ED (Ecosystem Demography) gap-oriented models – could reproduce the resilience to moderate disturbance observed in an experimentally manipulated forest (the Forest Accelerated Succession Experiment in northern Michigan, USA, in which 38% of canopy dominants were stem girdled and compared to control plots). Each model was parameterized, spun up, and disturbed following similar protocols and run for 5 years post-disturbance. The models replicated observed declines in aboveground biomass well. Biome-BGC captured the timing and rebound of observed leaf area index (LAI), while ZELIG and ED correctly estimated the magnitude of LAI decline. None of the models fully captured the observed post-disturbance C fluxes, in particular gross primary production or net primary production (NPP). Biome-BGC NPP was correctly resilient but for the wrong reasons, and could not match the absolute observational values. ZELIG and ED, in contrast, exhibited large, unobserved drops in NPP and net ecosystem production. The biological mechanisms proposed to explain the observed rapid resilience of the C cycle are typically not incorporated by these or other models. It is thus an open question whether most ecosystem models will simulate correctly the gradual and less extensive tree mortality characteristic of moderate disturbances.

Polymerization of high internal phase emulsions (polyHIPEs) is a relatively new method for the production of high-porosity scaffolds. The tunable architecture of these polyHIPE foams makes them attractive candidates for tissue engineered bone grafts. Previously studied polyHIPE systems require either toxic diluents or high cure temperatures which prohibit their use as an injectable bone graft. In contrast, we have developed an injectable polyHIPE that cures at physiological temperatures to a rigid, high-porosity foam. First, a biodegradable macromer, propylene fumarate dimethacrylate (PFDMA), was synthesized that has appropriate viscosity and hydrophobicity for emulsification. The process of surfactant selection is detailed with particular focus on the key structural features of both polymer (logP values, hydrogen bond acceptor sites) and surfactant (HLB values, hydrogen bond donor sites) that enable stable HIPE formation. Incubation of HIPEs at 37 °C was used to initiate radical cross-linking of the unsaturated double bond of the methacrylate groups to polymerize the continuous phase and lock in the emulsion geometry. The resulting polyHIPEs exhibited ~75% porosity, pore sizes ranging from 4 to 29 μm, and an average compressive modulus and strength of 33 and 5 MPa, respectively. These findings highlight the great potential of these scaffolds as injectable, tissue engineered bone grafts.