Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ecosystem modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Assessing Impacts of Plant Stoichiometric Traits on Terrestrial Ecosystem Carbon Accumulation Using the E3SM Land Model

Carbon (C) enters into the terrestrial ecosystems via photosynthesis and cycles through the system together with other essential nutrients (i.e., nitrogen [N] and phosphorus [P]). Such a strong coupling of C, N, and P leads to the theoretical prediction that limited nutrient availability will limit photosynthesis rate, plant growth, and future terrestrial C dynamics. However, the lack of reliable information about plant tissue stoichiometric constraints remains a challenge for quantifying nutrient limitations on projected global C cycling. In this study, we harmonized observed plant tissue C:N:P stoichiometry from more than 6,000 plant species with the commonly used plant functional type framework in global land models. Using observed C:N:P stoichiometry and the flexibility of these ratios as emergent plant traits, we show that observationally constrained fixed plant stoichiometry does not improve model estimates of present-day C dynamics compared with unconstrained stoichiometry. However, adopting stoichiometric flexibility significantly improves model predictions of C fluxes and stocks. The 21st century simulations with RCP8.5 CO 2 concentrations show that stoichiometric flexibility, rather than baseline stoichiometric ratios, is the dominant controller of plant productivity and ecosystem C accumulation in modeled responses to CO 2 fertilization. The enhanced nutrient limitations and plant P use efficiency mainly explain this result. This study is consistent with the previous consensus that nutrient availability will limit xfuture land carbon sequestration but challenges the idea that imbalances between C and nutrient supplies and fixed stoichiometry limit future land C sinks. We show here that it is necessary to represent nutrient stoichiometric flexibility in models to accurately project future terrestrial ecosystem carbon sequestration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Pyomo.GDP: an ecosystem for logic based modeling and optimization development

We present three core principles for engineering-oriented integrated modeling and optimization tool sets—intuitive modeling contexts, systematic computer-aided reformulations, and flexible solution strategies—and describe how new developments in Pyomo.GDP for Generalized Disjunctive Programming (GDP) advance this vision. We describe a new logical expression system implementation for Pyomo.GDP allowing for a more intuitive description of logical propositions. The logical expression system supports automated reformulation of these logical constraints to linear constraints. We also describe two new logic-based global optimization solver implementations built on Pyomo.GDP that exploit logical structure to avoid “zero-flow” numerical difficulties that arise in nonlinear network design problems when nodes or streams disappear. These new solvers also demonstrate the capability to link to external libraries for expanded functionality within an integrated implementation. We present these new solvers in the context of a flexible array of solution paths available to GDP models. Finally, we present results on a new library of GDP models demonstrating the value of multiple solution approaches.

42 ENGINEERING↗
DOE OSTI2022

Pathways to a Sustainable Aviation Ecosystem: Electrified Aviation Demand Modeling

With leading energy experts, high-performance computing, and powerful analysis and modeling capabilities, the National Renewable Energy Laboratory (NREL) empowers aviation stakeholders with tools and data to project future energy demand and explore supply options. As a result, NREL can help electric utilities and airport operations identify aviation electrification opportunities and prepare for future electricity needs.

aerial vehicles↗
DOE OSTI2020

Landsat near-infrared (NIR) band and ELM-FATES sensitivity to forest disturbances and regrowth in the Central Amazon

Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data for The Value of Reversible Carbon Storage in a Zero-Emissions World

Atmospheric carbon dioxide removal (CDR) is required to stabilize global temperature. CDR can be achieved via ecosystem-based approaches that are cost-effective but reversible (e.g., soil and forest management) or by more durable but expensive approaches (e.g., direct air capture coupled with geologic storage). Here, we examine trade-offs between these approaches, focusing on timing, climate impacts, and cost. We simulated reversible carbon accrual for a range of CDR contract structures using a general minimalist model of ecosystem carbon cycling, and parameterized it to simulate US agricultural soil management─specifically cover cropping─as a case study. We then quantified the resulting impact on atmospheric carbon and global temperature using a climate model emulator. We find that maintaining a patchwork of reversible CDR projects by replacing lapsed projects with new projects can reduce warming by 22–195 μ°C in 2100 and that the magnitude of this cooling effect depends on how effectively the patchwork is maintained. Long-term maintenance of reversible CDR projects requires institutional stability that cannot be guaranteed over multiple decades. Consequently, effective CDR ultimately requires replacing reversible projects with durable projects. To address this problem, we modeled the cost of replacing reversible agricultural soil CDR with geologic CDR. We found that using reversible CDR as a bridge to durable CDR is potentially more cost-effective as a global cooling strategy (0.20–0.81 billion USD per μ°C avoided) than perpetual maintenance of reversible CDR (0.32–1.31 billion USD per μ°C avoided) or an immediate transition to durable CDR (1.37–2.19 billion USD per μ°C avoided). However, we emphasize that institutional commitments to maintain reversible CDR projects cannot be guaranteed. Reliance on reversible CDR as a bridge to durable CDR therefore carries an unknown amount of risk and will only function if efforts to maintain reversible CDR are robust.

Carbon↗
DOE OSTI2022

Beyond carbon flux partitioning: Carbon allocation and nonstructural carbon dynamics inferred from continuous fluxes

Carbon (C) allocation and nonstructural carbon (NSC) dynamics play essential roles in plant growth and survival under stress and disturbance. However, quantitative understanding of these processes remains limited. Here, in this work, we propose a framework where we connect commonly measured carbon cycle components (eddy covariance fluxes of canopy CO 2 exchange, soil CO 2 efflux, and allometry-based biomass and net primary production) by a simple mass balance model to derive ecosystem-level NSC dynamics (NSC i ), C translocation (dC i ), and the biomass production efficiency (BPE i ) in above- and belowground plant (i = agp and bgp) compartments. We applied this framework to two long-term monitored loblolly pine (Pinus taeda) plantations of different ages in North Carolina and characterized the variations of NSC and allocation in years under normal and drought conditions. The results indicated that the young stand did not have net NSC flux at the annual scale, whereas the mature stand stored a near-constant proportion of new assimilates as NSC every year under normal conditions, which was comparable in magnitude to new structural growth. Roots consumed NSC in drought and stored a significant amount of NSC post drought. The above- and belowground dC i and BPE i varied more from year to year in the young stand and approached a relatively stable pattern in the mature stand. The belowground BPE bgp differed the most between the young and mature stands and was most responsive to drought. With the internal C dynamics quantified, this framework may also improve biomass production estimation, which reveals the variations resulting from droughts. Overall, these quantified ecosystem-scale dynamics were consistent with existing evidence from tree-based manipulative experiments and measurements and demonstrated that combining the continuous fluxes as proposed here can provide additional information about plant internal C dynamics. Given that it is based on broadly available flux data, the proposed framework is promising to improve the allocation algorithms in ecosystem C cycle models and offers new insights into observed variability in soil–plant–climate interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Kokkos EcoSystem: Comprehensive Performance Portability For High Performance Computing

State of the art Engineering and Science codes have grown in complexity dramatically over the last two decades. As a consequence application teams have adopted more sophisticated development strategies, leveraging third party libraries, deploying comprehensive testing and using advanced debugging and profiling tools. In todays environment of diverse hardware platforms, these applications also desire performance portability - avoiding the need to duplicate work for various platforms - which makes it necessary that these tools and libraries also work across the various systems. The Kokkos EcoSystem provides that portable software stack. Based on the Kokkos Core Programming Model, the EcoSystem provides math libraries, interoperability capabilities with Python and Fortran, and Tools for analysing, debugging, and optimizing applications. In this paper we will provide an overview of the components, discuss some specific use cases, and highlight how co-designing these components enables a more developer friendly experience.

42 ENGINEERING↗
DOE OSTI2021

Plant-microbe interactions: from genes to ecosystems using Populus as a model system

Plant-microbe symbioses span a continuum from pathogenic to mutualistic with functional consequences for both organisms in the symbiosis. In order to increase sustainable food and fuel production in the future, it is imperative that we harness these symbioses. The tree genus Populus is an excellent model system for studies examining plant-microbe interactions due to the wealth of genomic information available and the molecular tools that have been developed to manipulate Populus-microbe symbioses. In this review, we highlight how Populus can serve as a model system to explore plant-microbe interactions. Specifically, highlighting research linking Populus-microbe interactions from the gene to the ecosystem level. We explore why Populus is an excellent model for perennial plant systems, the molecular underpinnings of Populus-microbe interactions, how host genetics influence microbial community composition, and how microbial communities vary at fine spatial scales and between Populus species. Further, we explore how the patterns of the microbiome may affect ecosystem level functions in managed and natural ecosystems. Understanding and manipulating these interactions in Populus has the potential to improve plant health and impact ecosystem sustainability and processes as Populus trees function as foundational species in many natural ecosystems and are also deployed in managed ecosystems for various agroforestry applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year -1 (1Tg = 10 12 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH 4 year -1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of -34 to -46 Tg CH 4 year -1 , we estimated that the global net land methane emissions are 149–176 Tg CH 4 year -1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

58 GEOSCIENCES↗
DOE OSTI2025

Degradation and deforestation increase the sensitivity of the Amazon Forest to climate extremes

About 40% of the Brazilian Amazon has been deforested or suffered changes in forest structure through degradation (selective logging, fires, and fragmentation). The impact of forest degradation on the forest’s sensitivity to climate extremes has not been fully explored because of a lack of data and the complex interplay of forest structure and climate drivers. Here, we combined forest structure data from 545 airborne lidar transects (375 ha each) across the Brazilian Amazon with the Ecosystem Demography Model (ED2). We explore the forest’s functional response to near-present (1981–2019) climate extremes under observed forest structure from lidar ( Control ) and two forest structure change scenarios: (1) forest recovery by excluding all future deforestation and degradation ( Recovery ) and (2) expansion of selective logging and deforestation ( Degradation ). Using the Control simulation, we found a close and positive association between local forest aboveground biomass and the predicted gross primary productivity (GPP) and evapotranspiration (ET). Moreover, both GPP and ET respond negatively to extremes in vapor pressure deficit and downwelling shortwave irradiance in degraded forests in Eastern and Southern Amazon, indicating high sensitivity to droughts. Locally high-biomass forest patches showed little or no negative response of GPP and ET to extreme drought conditions whereas low-biomass forest patches in the same locations—typically degraded forest canopies—responded negatively to higher moisture stress. The results from the Recovery scenario showed similar results to simulations with observed structure; however, under the Degradation scenario, low-biomass forest patches became more abundant, resulting in more regions where GPP and ET are negatively impacted by hot drought conditions according to the ED2 model. Our results suggest that local forest structure is a critical determinant of an ecosystem’s response to climate variability, and that the loss of canopy trees in the Amazon through forest degradation could increase and expand forest vulnerability to droughts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Cutting out the middleman: calibrating and validating a dynamic vegetation model (ED2-PROSPECT5) using remotely sensed surface reflectance

Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant processes as well as understanding land–atmosphere feedbacks. In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty. Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified. Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state. We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Trade-offs among restored ecosystem functions are context-dependent in Mediterranean-type regions

Global biodiversity hotspots, including Mediterranean-type ecosystems worldwide, are highly threatened by global change that alters biodiversity, ecosystem functions, and services. Some restoration activities enhance ecosystem functions by reintroducing plant species based on known relationships between plant traits and ecosystem processes. Achieving multiple functions across different site conditions, however, requires understanding how abiotic factors like climate and soil, along with plant assemblages, influence ecosystem functions, including their trade-offs and synergies. We used the ModEST ecosystem simulation model, which integrates carbon, water, and nutrient processes with plant traits, to assess the relationships between restored plant assemblages and ecosystem functions in Mediterranean-type climates and soils. We investigated whether maximised carbon increment, water use efficiency, and nitrogen use efficiency, along with their trade-offs and synergies, varied across different abiotic contexts. Further, we asked whether assemblages that maximised functions varied across environments and among these functions. We found that maximised ecosystem carbon increment and nitrogen use efficiency occurred under moist, warm conditions, while water use efficiency peaked under drier conditions. Generally, the assemblage that maximised one function differed from those for other maximised functions. Synergies were rare, except between water and nitrogen use efficiencies in loam soils across most climates. Trade-offs among maximised functions were common, varying in strength with abiotic context and plant assemblages, and were more pronounced in sandy loam soils compared to clay-rich soils. Our findings suggest that due to variation in abiotic conditions within and across Mediterranean-type regions at the global scale, site-specific plant assemblages are required to maximise ecosystem functions. Thus, lessons from a single site cannot be transferred to another site, even where the same plant functional types are available for restoration. Our simulation results offer valuable insights into potential ecosystem performance under specific abiotic conditions following restoration with particular plant functional types, thereby informing local restoration efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

COMPASS-FME Synoptic Site Characterization

This dataset contains soil biogeochemical and physicochemical characterization data for the COMPASS-FME synoptic sites.This dataset also contains data for the paper Patel et al. 2025 "Transition zones at the changing coastal terrestrial-aquatic interface", https://doi.org/10.1029/2025JG008978.Coastal soils are a significant but highly uncertain component of global biogeochemical cycles. These systems experience unique spatial and temporal variability in biogeochemical processes, driven by wetland-to-upland gradients and hydrological fluctuations. We studied drivers of coastal soil variability (a) at regional scales and (b) across transects from upland forest to wetland, in two contrasting regions — Lake Erie, a freshwater lacustrine system, and Chesapeake Bay, a saltwater estuarine system. Salinity-related analytes were a key driver of soil variability, not just in the saltwater system, but surprisingly, also in the freshwater system. We had hypothesized linear trends in biogeochemical parameters along the TAI – however, contrary to expectations, transition soils were not consistently intermediate between upland and wetland endmembers; the non-monotonic trends of carbon, phosphorus, iron along our transects suggest that these are key analytes to study in our regions. Rapidly changing soil factors across coastal gradients provide insights into which soil processes may act as precursors to ecosystem shifts. Our comprehensive soil characterization across the coastal transects provides essential data for mechanistic modeling of ecosystem dynamics.The data are provided as processed, csv files. Raw data and processing scripts can be accessed on GitHub (https://github.com/COMPASS-DOE/cmps-soil_characterization).A note on the nomenclature: the experimental design represents three points along the coastal gradient -- upland, transition, and wetland. "wetland" is referred to as "marsh" in the corresponding paper. The two terms can be used interchangeably for the sites in this study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Consequences of Altered Root Nutrient Uptake for Soil Carbon Stabilization (Final Report)

Objectives: The primary objectives of this project are to improve our understanding of tropical forest belowground processes, and to increase representation of this biome in ecosystem-scale and global C cycle models. Belowground dynamics present a major source of uncertainty inhibiting our ability to predict C cycle responses to climate change. Representation of tropical forests in these models is particularly lacking, even though these ecosystems hold >25% of terrestrial C stocks. While the majority of humid tropical forests exist on relatively infertile soils (i.e. poor in rock-derived nutrients), enormous gradients in soil nutrient availability exist at landscape-scales due to shifts in geology and soil order. Soil fertility is likely to greatly influence how tropical soil C storage will respond to the declines in precipitation predicted for the tropics. In particular, root dynamics vary across soil fertility gradients, with roots representing the major input of C to subsurface soils. Root characteristics related to soil fertility include biomass, turnover, C exudates, tissue chemistry, and nutrient uptake rates, with each of these likely sensitive to changes in moisture. The proposed project will fully integrate field research with an existing ecosystem model to assess potential effects of drying on belowground tropical C dynamics across a range of soil fertilities. Project Description:. This project will identify key linkages among soil fertility, root nutrient uptake, root dynamics, and soil C storage. The hypothesized framework for understanding these linkages is: Tropical rainforest roots in relatively fertile soils have lesser root biomass, faster turnover, fewer exudates, and improved root tissue quality relative to infertile sites, all resulting from increased root nutrient uptake. Ultimately, these root characteristics in fertile soils promote proportionally greater long-term soil C storage in organo-mineral associations. Specifically, lower root exudates and improved root tissue quality in fertile soils reduce microbial respiration of extant soil C (i.e. “priming) and increase microbial C use efficiency, leading to sorption of protein-rich microbial and plant organic matter on soil mineral surfaces. However, the larger pools of mineral-associated C in fertile tropical soils are more vulnerable to decreased rainfall than smaller C pools in infertile soils, in part because of the lesser root biomass, faster turnover, and softer, more labile root tissues. Also, the shallower lateral rooting structures typical in fertile soils experience greater desiccation and death than deeper roots in infertile sites, with little short-term adaptation to drying.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Non-Destructive, Three-Dimensional Imaging of Processes in the Rhizosphere Utilizing High Energy Photons

Soil structure, which can be described as the aggregation and distribution of pore spaces, regulates carbon, nutrient, and water cycling across the Earth system. Yet the inability to make quantitative, dynamic, in situ measurements of soil structure and rhizosphere carbon flow has prevented meaningful incorporation of soil structural processes into Earth System Models (ESMs). The key limitations are: (i) Lack of scale integration between micron-scale soil structure and ecosystem-scale models, (ii) Poor functional linkage between soil structural properties and biogeochemical processes, and (iii) Absence of dynamic 3D measurements of rhizosphere structural changes and carbon transformations. To address these challenges, we developed a new integrated positron emission tomography (PET) - microcomputed tomography (CT) imaging platform that enables the first 4D (spatial 3D and time), non-invasive, quantitative imaging of carbon allocation and rhizosphere structural dynamics in living plants and intact soils. The system combines: • Rhizo-PET (R-PET): a high-resolution positron emission tomography scanner optimized for 11 CO 2 tracing in plant roots. • a-Se micro-CT: a high-contrast, high-spatial resolution CT system based on amorphous selenium (a-Se) direct conversion technology, enabling micron-scale visualization of soil structure and root–soil interfaces. Together, these advances allow us to quantify how carbon exudates move, transform, and stabilize within the rhizosphere, directly informing missing processes in BER-relevant carbon cycle models.

42 ENGINEERING↗
DOE OSTI2020

Robust paths to net greenhouse gas mitigation and negative emissions via advanced biofuels

Biofuel and bioenergy systems are integral to most climate stabilization scenarios for displacement of transport sector fossil fuel use and for producing negative emissions via carbon capture and storage (CCS). However, the net greenhouse gas mitigation benefit of such pathways is controversial due to concerns around ecosystem carbon losses from land use change and foregone sequestration benefits from alternative land uses. Here, we couple bottom-up ecosystem simulation with models of cellulosic biofuel production and CCS in order to track ecosystem and supply chain carbon flows for current and future biofuel systems, with comparison to competing land-based biological mitigation schemes. Analyzing three contrasting US case study sites, we show that on land transitioning out of crops or pasture, switchgrass cultivation for cellulosic ethanol production has per-hectare mitigation potential comparable to reforestation and severalfold greater than grassland restoration. In contrast, harvesting and converting existing secondary forest at those sites incurs large initial carbon debt requiring long payback periods. We also highlight how plausible future improvements in energy crop yields and biorefining technology together with CCS would achieve mitigation potential 4 and 15 times greater than forest and grassland restoration, respectively. Finally, we show that recent estimates of induced land use change are small relative to the opportunities for improving system performance that we quantify here. While climate and other ecosystem service benefits cannot be taken for granted from cellulosic biofuel deployment, our scenarios illustrate how conventional and carbon-negative biofuel systems could make a near-term, robust, and distinctive contribution to the climate challenge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The impact of hurricane disturbances on a tropical forest: implementing a palm plant functional type and hurricane disturbance module in ED2-HuDi V1.0

Hurricanes commonly disturb and damage tropical forests. Hurricane frequency and intensity are predicted to change under the changing climate. The short-term impacts of hurricane disturbances to tropical forests have been widely studied, but the long-term impacts are rarely investigated. Modeling is critical to investigate the potential response of forests to future disturbances, particularly if the nature of the disturbances is changing with climate. Unfortunately, existing models of forest dynamics are not presently able to account for hurricane disturbances. Therefore, we implement the Hurricane Disturbance in the Ecosystem Demography model (ED2) (ED2-HuDi). The hurricane disturbance includes hurricane-induced immediate mortality and subsequent recovery modules. The parameterizations are based on observations at the Bisley Experimental Watersheds (BEW) in the Luquillo Experimental Forest in Puerto Rico. We add one new plant functional type (PFT) to the model – Palm, as palms cannot be categorized into one of the current existing PFTs and are known to be an abundant component of tropical forests worldwide. The model is calibrated with observations at BEW using the generalized likelihood uncertainty estimation (GLUE) approach. The optimal simulation obtained from GLUE has a mean relative error of –21 %, –12 %, and –15 % for stem density, basal area, and aboveground biomass, respectively. The optimal simulation also agrees well with the observation in terms of PFT composition (+1 %, –8 %, –2 %, and +9 % differences in the percentages of “Early”, “Mid”, “Late”, and “Palm” PFTs, respectively) and size structure of the forest (+0.8 % differences in the percentage of large stems). Lastly, using the optimal parameter set, we study the impact of forest initial condition on the recovery of the forest from a single hurricane disturbance. The results indicate that, compared to a no-hurricane scenario, a single hurricane disturbance has little impact on forest structure (+1 % change in the percentage of large stems) and composition (<1 % change in the percentage of each of the four PFTs) but leads to 5 % higher aboveground biomass after 80 years of succession. The assumption of a less severe hurricane disturbance leads to a 4 % increase in aboveground biomass.

54 ENVIRONMENTAL SCIENCES↗