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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quality Assurance Project Plan for the Salish Sea Model – Continuing Development of New Capabilities and Applications: European Green Crab (EGC) Larval Transport, Coupling to VELMA and Atlantis, Performance Improvements and Diagnostic Applications

This quality assurance project plan (QAPP) is for proposed FY24-FY27 efforts to (a) develop/expand a predictive European Green Crab (EGC) larval dispersal model for the Salish Sea, (b) add linkage between VELMA (Visualizing Ecosystem Land Management Assessments) watershed model and Atlantis model (ecosystem dynamics and food web) through SSM, and (c) continue SSM capability development through the addition of metals, harmful algal blooms (HABs), and FVCOM-LTRANS (Lagrangian Larval Transport) modules.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Comprehensive Framework for Valuation of Ecosystem Services in the Context of Hydropower

While models to value ecosystem services exist, they suffer from a lack of “portability” in the sense that a single modelling framework developed for one site cannot be imported and applied to another. Additionally, benefit estimations are often fragmented, incomplete, incomparable, and in heterogenous metrics, thus preventing aggregation to arrive at one comprehensive value. In the context of hydropower, the methods, as well as the metrics used to value the associated ecosystem services, vary between stakeholders (e.g., competing water users, agencies, etc.), detracting from the ability to assess the total value of a hydropower project, water management schemes, or comparing value between competing water users. Further, current methods to value water neglect non-market ecosystem services’ values by failing to account for externalities, resulting in misallocation of costs (e.g., obligatory stakeholder payments for water use) and misinterpretation of hydropower benefits. The benefit-cost ratios in hydropower development are often incomplete because non-market benefits are excluded (e.g., benefits, such as fish habitat resiliency and agricultural water reliability, derived from hydropower development). An approach to enable consistent, standardized hydropower benefit-cost analyses does not currently exist. This paper conceptualizes a methodology to standardize ecosystem services valuation in the context of hydropower projects through the integration of the U.S. Environmental Protection Agency’s existing Final Ecosystem Goods and Services Classification System framework and economic valuation techniques, and demonstrates the applicability of this methodology through a limited proof-of-concept application (the New Waddell Plant, a pumped hydroelectric energy storage facility in Phoenix, AZ). The beauty of our proposed framework lies in its ability to be imported and applied to any other pumped storage facility after accounting for site-specific features.

13 HYDRO ENERGY↗
DOE OSTI2020

Moisture availability mediates the relationship between terrestrial gross primary production and solar-induced chlorophyll fluorescence: Insights from global-scale variations

Effective use of solar-induced chlorophyll fluorescence (SIF) to estimate and monitor gross primary production (GPP) in terrestrial ecosystems requires a comprehensive understanding and quantification of the relationship between SIF and GPP. To date, this understanding is incomplete and somewhat controversial in the literature. Here we derived the GPP/SIF ratio from multiple data sources as a diagnostic metric to explore its global-scale patterns of spatial variation and potential climatic dependence. We found that the growing season GPP/SIF ratio varied substantially across global land surfaces, with the highest ratios consistently found in boreal regions. Spatial variation in GPP/SIF was strongly modulated by climate variables. The most striking pattern was a consistent decrease in GPP/SIF from cold-and-wet climates to hot-and-dry climates. We propose that the reduction in GPP/SIF with decreasing moisture availability may be related to stomatal responses to aridity. Furthermore, we show that GPP/SIF can be empirically modeled from climate variables using a machine learning (random forest) framework, which can improve the modeling of ecosystem production and quantify its uncertainty in global terrestrial biosphere models. Finally, our results point to the need for targeted field and experimental studies to better understand the patterns observed and to improve the modeling of the relationship between SIF and GPP over broad scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Monthly averages of ED2 model simulations initialized with airborne lidar structure, Jan 1981-Dec 2018, Brazilian Amazon

Deforestation and forest degradation (selective logging, fires, fragmentation) have impacted nearly 40% of the original extent of the Brazilian Amazon, and have markedly impacted forest structure across the region. To date, few studies analysed how shifts in forest structure from degradation influence the forest sensitivity to climate extremes, because of the complex interactions between forest structure and micro-environmental conditions. To address this knowledge gap, we carried out a series of simulations across the Brazilian Amazon using the Ecosystem Demography Model (ED2), using observed forest structure derived from 541 airborne lidar transects (375 ha each) and two scenarios representing forest recovery and expansion of degradation to investigate how shifts in forest structure impact ecosystem function under near-average and extreme climate conditions, as part of the manuscript Longo et al 2025 "Degradation and Deforestation Increase the Sensitivity of the Amazon Forest to Climate Extremes". This dataset provides the output results from the ED2 model simulations for the three simulations at monthly time scales, in NetCDF format. For all simulations, we used bias-corrected hourly reanalyses (WFDE5) for most meteorological drivers, except for precipitation, which was obtained from CHIRPS. The meteorological drivers used in the study span 38 years (Jan 1981–Dec 2018). The output results correspond to the last 38 years of simulation (one full cycle of meteorological drivers), in which ED2 simulations used static stand structure (i.e., the forest structure was held constant). The following files are provided:ED2_emean_Global_R004_BrAmaz_s1c0t0l0f0.nc. This corresponds to the Control simulation. The forest structure was obtained from the airborne lidar.ED2_emean_Global_R005_BrAmaz_s1c0t1l1f0.nc. This corresponds to the Degraded simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that expanded deforestation and selective logging across the Amazon.ED2_emean_Global_R006_BrAmaz_s1c0t1l0f0.nc. This corresponds to the Recovery simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that completely halted deforestation and degradation, allowing degraded forests to recover for 38 years.We also provide file ED2_zones_R004_BrAmaz_s1c0t0l0f0.nc, which classifies each grid cell into zones used in the reference manuscript: 1: Southeast. 2: South. 3: West. 4: Central. 5: Northeast. 6: North. 7: Northwest". Index 0 corresponds to grid cells excluded from sub-region analyses because they were dominated by flooded forests, deforestation, and naturally non-forest vegetation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Across‐model spread and shrinking in predicting peatland carbon dynamics under global change

Abstract Large across‐model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across‐model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO 2 levels and five temperature levels. We developed a matrix‐based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter‐model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter‐model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO 2 (eCO 2 ) concentrations in the measurements, simulated NEP responded positively to eCO 2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across‐model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.

Hou, Enqing↗
DOE OSTI2020

DayCent data and results for "Robust paths to net greenhouse gas mitigation and negative emissions via advanced biofuels"

DayCent data and results for:J. L. Field, T. L. Richard, E. A. Smithwick, H. Cai, M. S. Laser, D. S. LeBauer, S. P. Long, K. Paustian, Z. Qin, J. J. Sheehan, P. Smith, M. Q. Wang, L. R. Lynd, Robust paths to net greenhouse gas mitigation and negative emissions via advanced biofuels. Proceedings of the National Academy of Sciences (2020). https://doi.org/10.1073/pnas.1920877117 This zip file contains a UNIX-format DayCent model executable, input files, automation code, and associated directory structure necessary to re-produce the DayCent analysis underlying the manuscript. The main script 'autodaycent.py' (written for Python 2.7) opens an interactive command line routine that facilitates:* Calibrating the DayCent pine growth model.* Initializing DayCent for a set of case studies sites.* Executing an ensemble of model runs representing case study site reforestation, grassland restoration, or conversion to switchgrass cultivation. * Results analysis & generation of manuscript Fig. 3. Note that the interactive analysis code requires that all input files to be contained in the directory structure as uploaded, without modification. Executable versions of the DayCent model (https://www.nrel.colostate.edu/projects/daycent/) compatible with other operating systems are available upon request.Please send questions/comments to John.L.Field@gmail.com

Agro-ecosystem function and prediction↗
DOE OSTI2026

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 capture and storage↗
DOE OSTI2022

Confronting the water potential information gap

Water potential directly controls the function of leaves, roots and microbes, and gradients in water potential drive water flows throughout the soil–plant–atmosphere continuum. Notwithstanding its clear relevance for many ecosystem processes, soil water potential is rarely measured in situ, and plant water potential observations are generally discrete, sparse, and not yet aggregated into accessible databases. These gaps limit our conceptual understanding of biophysical responses to moisture stress and inject large uncertainty into hydrologic and land-surface models. Here, we outline the conceptual and predictive gains that could be made with more continuous and discoverable observations of water potential in soils and plants. We discuss improvements to sensor technologies that facilitate in situ characterization of water potential, as well as strategies for building new networks that aggregate water potential data across sites. Here, we end by highlighting novel opportunities for linking more representative site-level observations of water potential to remotely sensed proxies. Together, these considerations offer a road map for clearer links between ecohydrological processes and the water potential gradients that have the ‘potential’ to substantially reduce conceptual and modelling uncertainties.

58 GEOSCIENCES↗
DOE OSTI2019

Plant roots stimulate the decomposition of complex, but not simple, soil carbon

Roots release carbon into soil and can alleviate energy limitation of microbial organic matter decomposition. We know little about the effects of roots on microbial decomposition of different organic matter substrates, despite the importance for soil carbon stocks and turnover. Through implementing root-microbe interactions, the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model was previously shown to represent dynamics of total soil carbon in temperate forest field experiments. However, the model permits alternative hypotheses concerning microbial-substrate affinity. We investigated how root inputs affect decomposition of soil organic carbon (SOC) with variable decomposability. We simulated SOC stocks in CORPSE and compared microbial degradation of two substrates types with varying root-microbe interactions under two alternative hypotheses that varied in microbial-substrate affinity. We compared our modeled hypotheses to a forest field experiment where we quantified decomposition of isotopically-labeled starch and leaf tissues in soils with manipulated root access to microbes. Here, we tested the hypothesis that decomposition of leaves would be more sensitive to root inputs than decomposition of starch, corresponding to the alternative model hypothesis. In the field study, leaf decomposition increased with root density while starch decomposition was unchanged by root density. Microbial biomass and enzyme activity consistently increased with root inputs in CORPSE and the field study. Our field experiment supported the CORPSE simulations with high microbial-substrate affinity. Roots stimulated microbial growth and enzyme production, which increased degradation of more complex substrates such as leaf tissues. Substrates that were easily decomposed, such as starch, may already be degrading at a maximum rate in the absence of rhizosphere influence because their decomposition rate was unchanged by root inputs. We found that the degree to which roots stimulate microbial decomposition depends on the substrate being decomposed, and that root-microbe interactions influenced SOC stocks in both our model and field experiment. Environmental changes that alter root-microbe interactions could, therefore, alter soil C stocks and biogeochemical cycling, and models of these interactions should incorporate differential influence of rhizosphere inputs on different substrates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The circular nutrient economy of terrestrial ecosystems and the consequences for rock weathering

Earth’s biosphere is thought to exert a substantial influence on regolith evolution and chemical weathering rates. However, ecosystems are also highly efficient at retaining and recycling nutrients. Thus, when the ecological demand for rock-derived nutrients (e.g., P, Ca, K) exceeds the rates of regolith supply, ecological retention and recycling strategies can minimize nutrient limitations. To evaluate the balance between nutrient recycling and new nutrient input, we combined a plant model that drives growth according to foliar P levels with a weathering model that includes regolith rejuvenation via erosion and export via chemical weathering according to water flow, regolith thickness, mineral dissolution rates, secondary minerals, and nutrient storage in organic and mineral phases. We find that plant growth is strongly dependent on the total regolith nutrient inventory, resulting in a strong correlation between plant productivity and erosion. Increased water export or decreased regolith thickness diminish the total inventory of nutrient corresponding to lower rates of recycling and lower plant growth. In contrast, purported biogenic drivers of weathering, such as enhanced mineral dissolution, only support higher growth rates at high erosion rates. At erosion rates typical of the global land surface, more rapid mineral dissolution combined with enhanced formation of secondary minerals, depletes the inventory of mineral P, resulting in no benefit for plant growth. We also find that the increased chemical weathering export does not scale directly with plant growth. For example, accelerated mineral weathering does increase chemical weathering export but not potential plant growth. Conversely, thicker regolith is associated with a small increase in weathering export, but a large increase in potential plant growth. Collectively, when plant growth is coupled to regolith weathering our calculations suggest that plant productivity is not directly correlated with silicate weathering fluxes, and that biotic drivers of silicate weathering may only be effective at high erosion rates not typical at the Earth’s surface.

58 GEOSCIENCES↗
DOE OSTI2022

Threat by marine heatwaves to adaptive large marine ecosystems in an eddy-resolving model

Marine heatwaves (MHWs), episodic periods of abnormally high sea surface temperature, severely affect marine ecosystems. Large marine ecosystems (LMEs) cover ~22% of the global ocean but account for 95% of global fisheries catches. Yet how climate change affects MHWs over LMEs remains unknown because such LMEs are confined to the coast where low-resolution climate models are known to have biases. Here, using a high-resolution Earth system model and applying a ‘future threshold’ that considers MHWs as anomalous warming above the long-term mean warming of sea surface temperatures, we find that future intensity and annual days of MHWs over the majority of the LMEs remain higher than in the present-day climate. Better resolution of ocean mesoscale eddies enables simulation of more realistic MHWs than low-resolution models. These increases in MHWs under global warming pose a serious threat to LMEs, even if resident organisms could adapt fully to the long-term mean warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Multispectral UAV imagery of experimental freshwater wetlands under 5 ppt saltwater intrusion, Louisiana, 2023 and 2024

Multispectral imagery was collected using an unmanned aerial vehicle (UAV) to evaluate how freshwater vegetation responds to short-term simulated saltwater intrusion events. The purpose of this data collection was to understand how plant health changes in response to acute salinity exposure, which is increasingly relevant in coastal wetland ecosystems facing sea level rise and storm surge events, such as in coastal Louisiana. Three experimental saltwater intrusions were conducted at a salinity of approximately 5 parts per thousand (ppt) for durations of 6-days, 10-days, and 17-days. UAV flights occurred both before and after each treatment. The resulting imagery was processed using Pix4DMapper software to georeference the images and generate orthomosaics. The multispectral sensor used in this study captures reflectance in five bands: blue, green, red, red-edge, and near-infrared. The uploaded data consist of georeferenced .tif orthomosaics for each spectral band, which are compatible with GIS software for vegetation analysis. This imagery can be utilized in investigations into vegetation stress, remote sensing of freshwater wetland ecosystems, and modeling of plant response to environmental changes.

EARTH SCIENCE > BIOSPHERE > ECOSYSTEMS↗
DOE OSTI2021

Mercury methylation by metabolically versatile and cosmopolitan marine bacteria

Microbes transform aqueous mercury (Hg) into methylmercury (MeHg), a potent neurotoxin that accumulates in terrestrial and marine food webs, with potential impacts on human health. This process requires the gene pair hgcAB, which encodes for proteins that actuate Hg methylation, and has been well described for anoxic environments. However, recent studies report potential MeHg formation in suboxic seawater, although the microorganisms involved remain poorly understood. In this study, we conducted large-scale multi-omic analyses to search for putative microbial Hg methylators along defined redox gradients in Saanich Inlet, British Columbia, a model natural ecosystem with previously measured Hg and MeHg concentration profiles. Analysis of gene expression profiles along the redoxcline identified several putative Hg methylating microbial groups, including Calditrichaeota, SAR324 and Marinimicrobia, with the last the most active based on hgc transcription levels. Marinimicrobia hgc genes were identified from multiple publicly available marine metagenomes, consistent with a potential key role in marine Hg methylation. Computational homology modelling predicts that Marinimicrobia HgcAB proteins contain the highly conserved amino acid sites and folding structures required for functional Hg methylation. Furthermore, a number of terminal oxidases from aerobic respiratory chains were associated with several putative novel Hg methylators. Our findings thus reveal potential novel marine Hg-methylating microorganisms with a greater oxygen tolerance and broader habitat range than previously recognized.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Multi-Physics Coupled Seismic Safety Analysis of Molten-Salt Reactors

Pool-type Molten Salt Reactors (MSRs) are a promising advanced nuclear reactor concept relying on passive physical behavior to provide increased safety characteristics. The high heat capacity, high boiling point, unpressurized liquid salt contains the nuclear fuel and as it passes through the vessel achieves criticality and produces nuclear power. The salt then carries away the heat produced along with the delayed neutron precursors towards heat exchangers and primary pumps before entering the core. This forms a highly coupled multiphysics problem between neutronics and fluid flow, which makes modeling these reactors challenging. The safety/licensing case heavily relies on modeling and simulation to prove their safety, and significant effort has been expanded by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to develop the appropriate simulation tools. Their safety/licensing will require a thorough analysis of their behavior during earthquake-based transients. Earthquakes are a transient condition commonly experienced in many regions, and are safely survived by dozens of reactors every year. Nonetheless, significant conservatisms were introduced with limited predictive modeling and simulation available. The first-of-a-kind capability developed by this research will allow for engineers to analyze ranges of design-basis and beyond-design basis accident scenarios to assess the integrity of the reactor vessel and surrounding system, possibly allowing for limiting the costly conservatisms in the design. This projects uses MASTODON, Griffin, and Pronghorn in the NEAMS ecosystem to model the coupled multiphysics problem posed by an earthquake in a molten salt reactor. It marches through a number of number of coupling schemes, from uncoupled simulations to tightly two-way coupled simulations. The appropriate level of coupling for these simulations will then be determined based on the accuracy in the quantities of interest (QoI) and computational effort involved in each scheme. The QoIs chosen for this paper include the mass flow rate across the heat exchanger (for Pronghorn simulations), power output of the reactor (for Griffin simulations), and the maximum Von Mises stress (for MASTODON simulations). This summary paper reports on the current progress of this project, with 2D tightly coupled multiphysics results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond (resistance) and recover from (resilience) to rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the ED-2.2 model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900, and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes, and interacted with one another. We then calculated metrics of functional ecosystem stability. Modeled gross primary production exhibited a mean resistance of -0.13, -0.23, and -0.42 (unitless) in response to increasing disturbance severity of 45%, 65%, and 85% mortality, respectively; and resilience of 0.002, 0.004, 0.008 (yr-1) respectively. This pattern held true for other carbon fluxes, meaning that in general simulated GPP declined more initially but recovered more quickly from higher mortality. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Atmospheric methane consumption in arid ecosystems acts as a reverse chimney and is accelerated by plant-methanotroph biomes

Drylands cover one-third of the Earth’s surface and are one of the largest terrestrial sinks for methane. Understanding the structure–function interplay between members of arid biomes can provide critical insights into mechanisms of resilience toward anthropogenic and climate-change-driven environmental stressors—water scarcity, heatwaves, and increased atmospheric greenhouse gases. This study integrates in situ measurements with culture-independent and enrichment-based investigations of methane-consuming microbiomes inhabiting soil in the Anza-Borrego Desert, a model arid ecosystem in Southern California, United States. The atmospheric methane consumption ranged between 2.26 and 12.73 μmol m 2 h −1 , peaking during the daytime at vegetated sites. Metagenomic studies revealed similar soil-microbiome compositions at vegetated and unvegetated sites, with Methylocaldum being the major methanotrophic clade. Eighty-four metagenome-assembled genomes were recovered, six represented by methanotrophic bacteria (three Methylocaldum , two Methylobacter , and uncultivated Methylococcaceae ). The prevalence of copper-containing methane monooxygenases in metagenomic datasets suggests a diverse potential for methane oxidation in canonical methanotrophs and uncultivated Gammaproteobacteria. Five pure cultures of methanotrophic bacteria were obtained, including four Methylocaldum . Genomic analysis of Methylocaldum isolates and metagenome-assembled genomes revealed the presence of multiple stand-alone methane monooxygenase subunit C paralogs, which may have functions beyond methane oxidation. Furthermore, these methanotrophs have genetic signatures typically linked to symbiotic interactions with plants, including tryptophan synthesis and indole-3-acetic acid production. Based on in situ fluxes and soil microbiome compositions, we propose the existence of arid-soil reverse chimneys, an empowered methane sink represented by yet-to-be-defined cooperation between desert vegetation and methane-consuming microbiomes.

59 BASIC BIOLOGICAL SCIENCES↗