Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Ecological 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 55 records · Page 3

Infection and Genomic Properties of Single- and Double-Stranded DNA Cellulophaga Phages

Bacterial viruses (phages) are abundant and ecologically impactful, but laboratory-based experimental model systems vastly under-represent known phage diversity, particularly for ssDNA phages. Here, we characterize the genomes and infection properties of two unrelated marine flavophages—ssDNA generalist phage phi18:4 (6.5 Kbp) and dsDNA specialist phage phi18:1 (39.2 Kbp)—when infecting the same Cellulophaga baltica strain #18 (Cba18), of the class Flavobacteriia. Phage phi18:4 belongs to a new family of ssDNA phages, has an internal lipid membrane, and its genome encodes primarily structural proteins, as well as a DNA replication protein common to ssDNA phages and a unique lysis protein. Phage phi18:1 is a siphovirus that encodes several virulence genes, despite not having a known temperate lifestyle, a CAZy enzyme likely for regulatory purposes, and four DNA methyltransferases dispersed throughout the genome that suggest both host modulation and phage DNA protection against host restriction. Physiologically, ssDNA phage phi18:4 has a shorter latent period and smaller burst size than dsDNA phage phi18:1, and both phages efficiently infect this host. These results help augment the diversity of characterized environmental phage–host model systems by studying infections of genomically diverse phages (ssDNA vs. dsDNA) on the same host.

Howard-Varona, Cristina (ORCID:0000000241495818)

Discrete global grid system-based flow routing datasets in the Amazon and Yukon basins

Abstract. Discrete global grid systems (DGGS) are emerging spatial data structures widely used to organize geospatial datasets across scales. While DGGS have found applications in various scientific disciplines, including atmospheric science and ecology, their integration into physically based hydrological models and Earth system models (ESMs) has been hindered by the lack of flow routing datasets based on DGGS. In response to this gap, this study pioneers the development of new flow routing datasets using icosahedral Snyder equal-area (ISEA) DGGS and a novel mesh-independent flow direction model. We present flow routing datasets for two large basins, the tropical Amazon River basin and the Arctic Yukon River basin. These datasets (1) facilitate the adoption of DGGS for hydrological models and (2) provide flow routing inputs for evaluation of DGGS-based flow routing in the Amazon and Yukon river basins. The data are available at https://doi.org/10.5281/zenodo.8377765 (Liao, 2023).

54 ENVIRONMENTAL SCIENCES

Microbiome Adaptation Could Amplify Modeled Projections of Global Soil Carbon Loss With Climate Warming

Warming alters soil microbial traits through ecological and evolutionary processes, directly influencing the decomposition of organic matter, which significantly affects global soil carbon emissions. Yet, soil carbon models largely ignore these processes and their implications for global responses to warming. Here, we incorporate eco-evolutionary theory into a mechanistic model describing microbial soil carbon decomposition to address the question of whether such processes could have consequential effects on climate carbon feedbacks globally. We assume that a key trait of microbes, their resource allocation to production of exoenzymes (which facilitate decomposition of organic matter)—is optimized to environmental temperatures by natural selection. We find that eco-evolutionary optimization results in microbes allocating more resources to enzyme production under warming. When applied at the global scale, eco-evolutionary optimization enhances the biological realism of soil carbon models and significantly amplifies global soil carbon loss by 2100. Our results highlight the significant potential of microbial eco-evolutionary responses to influence carbon cycle feedbacks to climate change, and motivate an urgent need for more comprehensive data to accurately quantify the adaptive potential of microbiomes in response to climate change.

Biodiversity & Conservation

Roadmap for the future of extreme wildfire events

Background Extreme wildfire events (EWEs) represent a growing threat globally, posing substantial risks to ecosystems, human communities, and infrastructure. Despite increased recognition of their ecological, social, and economic significance, current definitions of EWEs vary widely, reflecting disciplinary biases and regional contexts. This article emerges from an interdisciplinary workshop convened to reassess and refine the definition of EWEs, examine their impacts across ecological and social dimensions, and identify critical knowledge gaps impeding our understanding of these infrequent but important events. Results Our synthesis highlights significant limitations with existing definitions, particularly their reliance on subjective thresholds and their emphasis on extreme fire behavior alone. EWEs encompass a spectrum of complex, multi-dimensional phenomena that extend beyond immediate biophysical characteristics to include cumulative social, economic, and ecological impacts. These impacts often manifest over extended timeframes and include hazardous environmental contamination, severe geomorphic disturbances, ecosystem transformations, and unintended consequences of post-fire management actions. Current wildfire modeling frameworks inadequately capture these compounding factors, particularly the interactions among social systems, ecological conditions, and extreme fire behavior. To overcome these issues, we advocate for an interdisciplinary and context-sensitive approach to defining and studying EWEs. This revised definition emphasizes wildfires exhibiting anomalies in fire behavior, ecological outcomes, or social impacts relative to historically observed baselines, accommodating variability across different geographic regions and ecological settings. Conclusions Adopting an interdisciplinary framework that integrates biophysical and social sciences will enhance the predictive capability of wildfire models and improve resilience planning and response strategies. Filling identified knowledge gaps—such as limited high-quality empirical fire behavior data and insufficient integration of social dynamics into modeling—will better prepare communities and ecosystems to cope with and adapt to EWEs. This inclusive approach underscores the necessity for collaboration across disciplines and sectors, essential to managing extreme wildfires in an era of increasing climatic and ecological uncertainty.

54 ENVIRONMENTAL SCIENCES

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES

Comparative Study of Enhanced Geothermal System Supply Curves Across CONUS from Two Temperature Models Using the Renewable Energy Potential Model (reV)

Enhanced geothermal systems (EGS) have had recent breakthroughs within the geothermal sector. These breakthroughs are reflected in the National Renewable Energy Laboratory (NREL) 2024 Annual Technology Baseline and will result in updated EGS supply curves (i.e., the available resource capacity relative to cost). Our research uses NREL's Renewable Energy Potential (reV) model to compare EGS supply curves across the conterminous United States (CONUS) for two different temperature models: the Stanford temperature model (STM) and the Southern Methodist University temperature model (SMU). The reV model provides the levelized cost of energy (LCOE) at a consistent resolution across CONUS, taking into consideration transmission costs and constraints as well as technical exclusions pertaining to sensitive cultural, ecological, or infrastructure locations. In addition to a countrywide analysis of both models, we also conducted a regional analysis of Texas. We observed the STM had, on average, lower temperatures across different depths, resulting in slightly higher mean and median LCOEs as compared to the SMU temperature model at the same depths. In the regional analysis for Texas, however, when we compared only the common points between the two temperature models, the STM had lower median and mean LCOEs compared to SMU due to higher temperatures at depths greater than 5 km.

enhanced geothermal systems

Evidence for the existence and ecological relevance of fast-cycling mineral-associated organic matter

Longstanding theories and models classify mineral-associated organic matter as the large ( ~ 60%) but slow-cycling and persistent portion of soil organic matter. Strong physico-chemical interactions and diffusion limitations restrict the turnover of mineral-associated organic matter, allowing carbon and nitrogen bound therein to persist in soil for as long as centuries to millennia. However, mineral-associated organic matter is a chemically and functionally diverse pool with a substantial portion cycling at relatively fast (i.e., minutes to years) timescales. Despite a growing body of evidence for the heterogenous and multi-pool nature of mineral-associated organic matter, we lack consensus on how to conceptualize and directly quantify fast-cycling mineral-associated organic matter and its ecological significance. We demonstrate that the dynamic qualities of fast-cycling mineral-associated organic matter vary based on 1) the chemistry of the mineral particles and organic matter, 2) the complex set of interactions between organic matter and the mineral matrix, and 3) the presence and strength of destabilizing forces that lead to decomposition or loss of mineral-associated organic matter (i.e., plant-microbe interactions, agricultural intensification, and climate change). Finally, we discuss potential implications and research opportunities for how we measure, manage, and model the dynamic subfraction of this otherwise persistent pool of soil organic matter.

Biological and medical sciences

Structured interactions drive abrupt transitions in the spatial organization of microbial communities

Bacteria possess diverse mechanisms to regulate their motility in response to environmental and physiological signals, enabling them to navigate complex habitats and adapt their behavior. Some of these mechanisms are species specific and enable cells to modulate their movement based on the ecological identity of neighboring species. Here, we introduce a model in which bacteria interact via local signals that either enhance or suppress the motility of neighboring cells depending on species type. Through large-scale simulations and a coarse-grained stochastic model, we demonstrate the emergence of a sharp transition driven by nucleation processes: increasing the density of motility-suppressing interactions drives the system from a fully mixed, motile phase to a state characterized by large, stationary bacterial clusters. Remarkably, in systems with a large number of interacting species, this transition can be triggered solely by altering the structure of the motility-regulation interaction matrix while maintaining species and interaction densities constant. In particular, we find that heterogeneous and modular interactions promote the transition more readily than homogeneous random ones. These findings add a dimension to the theory of motility-induced phase separation and contribute to the ongoing effort to understand microbial interactions, suggesting that structured, nonrandom ones may be key to reproducing commonly observed spatial patterns in microbial communities.

bacterial communities

Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics Within Water‐Tagging Enabled Hydrologic Models

Determining the age distribution of water exiting a catchment is important for understanding groundwater storage and mixing. New water-tagging capabilities within models track precipitation events as they move through simulated storages, yet forward modeling of individual events may not systematically capture the full transit time distribution (TTD). Here, we present a “sequential precipitation input tagging” (SPIT) framework to tag all input precipitation at regular intervals during extended model simulations. Monthly tags over 7 years were applied at six National Ecological Observatory Network sites to calculate TTDs and derive mean virtual tracer age, $\overline{T_{V}}$, fractions of young water, F yw , and hydrologic tracer concentrations (water isotopes δ 18 O and δ 2 H) within a tagging enabled version of the Weather Research and Forecast hydrologic model (WRF-Hydro). Throughout seven simulation years, the fraction of simulated discharge derived from tagged events, F tag , increased each year, with the final year's F tag ranging from 66% to 100% and highlights the need to apply SPIT over many years to understand TTDs. When the F tag was >75%, simulated $\overline{T_{V}}$ ranged 179–923 days and F yw 0.6%–23.9%, with daily values exhibiting a power-law relationship with precipitation, discharge, and groundwater. Through implementation of SPIT, we find this hydrologic model configuration performs poorly in estimation of $\overline{T_{V}}$ and F yw (root mean squared error of 469 days and 14.4% respectively), suggesting it misrepresents subsurface mixing. Thus, the SPIT framework provides a reproducible approach to calculate watershed transit times within tagging enabled models and thereby assess and improve representation of hydrologic processes.

fraction of young water

DVMDOSTEM v0.8.3: a terrestrial ecosystem model designed to represent arctic, boreal and permafrost ecosystem dynamics

The impacts of climate change on natural ecosystems are the result of complex physical and ecological processes operating and interacting at a variety of spatio-temporal scales, that can be represented in process-based ecosystem models. DVMDOSTEM is an advanced process-based terrestrial ecosystem model (TEM) designed to study ecosystem responses to climate changes and disturbances. It has a particular focus on permafrost regions (i.e. regions characterized by soils that stay partially frozen all year round for at least two consecutive years), encompassing boreal, arctic, and alpine landscapes. The model couples two previous versions of the Terrestrial Ecosystem Model (TEM) (McGuire et al., 1992): DVMTEM that includes a dynamic vegetation module (DVM) (E. S. Euskirchen et al., 2009), and DOSTEM that includes a dynamic organic soil module (DOS) (H. Genet et al., 2013; Yi et al., 2010). DVMDOSTEM simulates processes at yearly and monthly scales, with some physical processes operating at an even finer temporal resolution. Its versatility allows for site-specific to regional simulations, making it valuable for predicting shifts in permafrost, vegetation, and carbon (C) and nitrogen (N) dynamics. While DVMDOSTEM has been described in the methods sections of many manuscripts, this paper is the first stand alone description of DVMDOSTEM, independent of a particular scientific investigation.

Carman, Tobey B. [Univ. of Alaska, Fairbanks, AK (

Initial Mobility Analysis for ORNL VA-EDH Synthetic Populations

Travel burdens are a major barrier to healthcare access among US Veteran patient populations, particularly those residing in rural areas. Spatial accessibility to points of care for US Veteran populations is commonly assessed in two ways. The first approach uses open data from the US Census to represent collective travel burdens, for example the distance between population-weighted census tract centroids and VHA points of care. The second approach uses restricted-access VHA patient data to measure travel costs (e.g., distance, time) for accessing points of care with respect to geolocated patient addresses and real or approximated transportation networks. While the advantage of the open data approach lies in its reproducibility, it has notable limitations in its tendency to infer individual travel behavior from aggregate population characteristics, a problem known as ecological fallacy. Conversely, while the patient data approach is able to account for individual travel behavior, its ability to account for localized access disparities (e.g., a neighborhood with exceptionally high transportation costs) and patient demographics is limited as protecting individual patient data requires their storage in closed systems with limited capacity for adequately modeling real-world travel patterns or for supplementing patient attributes. Additionally, the patient data approach cannot account for veterans who are not enrolled in the VHA system but who may be eligible for care. These challenges limit the ability to perform “what if” analyses on the effects of place-specific interventions on veteran populations with high access barriers to healthcare. To address these challenges, we explore the application of realistic synthetic populations to examine travel burdens and spatial accessibility issues among veteran patient populations. Synthetic populations provide a virtual, individually-resolved and cross-sectional representation of the veteran patient population that enables investigation of spatial access to points of care in ways in which aggregate data and patient data do not. First, synthetic populations allow one to directly assess how individuals access points of care, from synthesized residential locations to outpatient facilities on real-world transportation networks. Modeling access to points of care at the individual scale addresses the ecological fallacy problem associated with using aggregated census data to represent veteran populations and patterns of movement. Second, synthetic populations provide a means of completely representing an area’s veteran population using only publicly available, anonymized census microdata from the American Community Survey (ACS) to ensure the privacy of real-world individuals. Generating synthetic populations from the ACS also expands descriptive characteristics beyond what patient data typically offers to include socio-demographic, economic, housing, and mobility attributes. More detailed profiles of both VHA patient populations and veterans not enrolled in the VA system will provide a comprehensive picture of groups that may benefit from interventions or outreach. As an initial exercise for using synthetic populations to measure veteran travel burdens to VA care, we apply Oak Ridge National Laboratory’s (ORNL) UrbanPop capability to generate a series of synthetic VHA patient populations for 9 Veterans Integrated Services Networks (VISN) market areas in 9 Census Divisions across the continental United States, which are listed in Table 1. We use UrbanPop to produce synthetic populations for the VISN markets selected for each US Census Division, then assign VA outpatient clinic destinations to synthetic VHA patients based on travel about each VISN market’s road network. To demonstrate using the synthetic populations to evaluate healthcare travel burdens, we compare the time-based impedance between simulated home locations and VA outpatient clinics in each VISN market. We then perform validation exercises on the synthetic populations with respect to neighborhood (block group) demographic composition as well as patient mobility, comparing aggregate origin-destination statistics for the synthetic population to outpatient visits available in restricted patient data from the VA’s Corporate Data Warehouse (CDW) database.

97 MATHEMATICS AND COMPUTING

Synthetic communities as a model for determining interactions between a biofertilizer chassis organism and native microbial consortia

Biofertilizers are critical for sustainable agriculture because they can replace ecologically disruptive chemical fertilizers while improving the trajectory of soil and plant health. However, for improving deployment, the persistence of biofertilizers within native soil consortia must be elucidated and enhanced. In this study we characterized a high-throughput, modular, and automation-friendly in vitro approach to screen for biofertilizer persistence within soil-derived consortia after co-cultivation with stable synthetic soil microbial communities (SynComs) obtained through a top-down cultivation process. Here, we profiled ~1200 SynComs isolated from various soil sources and cultivated in divergent media types, and we detected significant phylogenetic diversity (e.g. Shannon index >4) and richness (observed richness >400) across these communities. We observed high reproducibility in SynCom community structure from common soil and media types, which provided a testbed for assessing biofertilizer persistence within representative native consortia. Furthermore, we demonstrated that the screening method described herein can be coupled with microbial engineering to efficiently identify soil-derived SynComs in which an engineered biofertilizer organism (i.e. Bacillus subtilis) persists. Accordingly, we discovered that B. subtilis persisted in ~10% of SynComs that generally followed the diversity–invasion principle. Additionally, our approach enabled analysis of the ecological impact of B. subtilis inoculation on SynCom structure and profile alterations in community diversity and richness associated with the presence of a genetically modified model bacterium. Ultimately, this work has established a modular pipeline that could be integrated into a variety of microbiology/microbiome-relevant workflows or related applications that would benefit from assessment of the persistence of a specific organism of interest and its interaction with native consortia.

biofertilizers

Ecological connectivity and habitat loss shape patterns of genetic diversity in a threatened salamander

Context The maintenance of genetic diversity is essential for preserving adaptive potential in populations, yet it is increasingly threatened by landscape alteration. The field of landscape genetics offers a framework for assessing how patch-level landscape conditions, modeled at multiple scales, influence genetic diversity. Objectives We sought to assess how local environmental features and connectivity influence genetic diversity across 74 four-toed salamander (Hemidactylium scutatum) breeding wetlands in the southeastern United States. Methods Using next-generation sequencing data and hierarchical Bayesian models, we examined genome-wide heterozygosity in relation to local landscape features and ecological connectivity. We also assessed the scale of effect of landscape features and tested for temporal lag effects. Results Genetic diversity was lower in wetlands with higher levels of historic deforestation and lower connectivity. An interaction between deforestation and connectivity indicated that deforestation had stronger negative effects in isolated wetlands but weaker effects in well-connected wetlands. Accounting for scale of effect and temporal lags was critical for detecting these relationships. Conclusions Our analyses highlight the importance of assessing the spatial scale (scale of effect) and temporal lag of landscape features to detect key drivers of genetic diversity. In line with population genetic theory, our results indicate that the genetic consequences of habitat loss do not affect populations uniformly and are most severe in isolated populations where gene flow cannot buffer against loss of diversity. Altogether, we highlight the importance of considering the interaction of habitat loss and connectivity in conservation genetic management.

Hemidactylium scutatum

Ecological Insights from Transferable Plant Biomass Mapping across the Arctic using High-resolution Structure-from-Motion and LiDAR Data

Warmer temperatures, permafrost thaw, and increased wildfire activity are driving rapid ecological change across the Arctic, significantly altering plant productivity and aboveground biomass (AGB). These rapid changes highlight the urgent need to improve monitoring of vegetation dynamics in the Earth’s northern ecosystems, where high spatiotemporal heterogeneity occurs at scales finer than those captured by traditional satellite observations. The growing use of Unoccupied Aerial Systems (UASs) presents an opportunity to overcome this limitation. Yet, the diversity of UAS platforms, sensors, and data collection and processing workflows presents challenges for developing standardized, generalizable approaches. To address this challenge, we compiled 672 AGB plots co-located with 183 UAS-based Structure-from-Motion (SfM) or Light Detection and Ranging (LiDAR) surveys collected across the Arctic. Here, we: (1) evaluated the generalizability of UAS-derived canopy structure derived from high-resolution SfM and LiDAR for estimating AGB, (2) assessed scaling errors and their sources in two recent satellite-based AGB products derived from Landsat and MODIS, and (3) demonstrated the use of high-resolution AGB maps to quantify biomass variation across tundra plant functional types (PFTs) and to monitor post-fire recovery. Our results show that both SfM and LiDAR accurately captured AGB and its variability across tundra PFTs using a Random Forest (RF) model (overall RMSE: 0.336 kg/m2), with mapping performance varying slightly by region and data source. Using UAS-derived AGB maps as a benchmark, we identified systematic biases in satellite-derived AGB products, largely attributable to the magnitude of AGB and structural heterogeneity within coarse-resolution pixels. Applying our model to repeat UAS surveys following a tundra fire on Seward Peninsula, we observed rapid AGB recovery in non-shrub patches, with biomass recovering to pre-fire levels within 2 years. In contrast, shrub patches recovered more slowly, with AGB gains continuing over 2–4 years through both in-patch growth and lateral expansion (via dispersal) into remaining burned areas. Overall, these findings demonstrate the generalizability of UAS-based SfM and LiDAR data for estimating tundra AGB and highlight the potential of our approach to be broadly applied to generate high-quality AGB data for ecological monitoring and model benchmarking across the Arctic.

Yang, Daryl [ORNL] (ORCID:0000000317057823)

TES: Modelling Microbes to Predict Post-Fire Carbon Cycling in the Boreal Forest across Burn Severities (Exploratory Proposal)

Wildfires in boreal ecosystems represent a globally significant ecological process that is expected to be sensitive to changing climate as well as forest management strategies, but is currently insufficiently understood and represented in models. We sought to determine whether linking belowground microbial community composition, size, and activity to aboveground properties of burn severity and plant community composition (upland jack pine and upland spruce) would allow us to better model post-fire soil CO 2 fluxes using the Carbon, Organisms, Rhizosphere and Protection in the Soil Environment (CORPSE) model. We used an integrated modelling experimental approach, with mechanistic laboratory experiments designed to inform models.

54 ENVIRONMENTAL SCIENCES

Tradeoffs between uniform land protection and biodiversity-specific land protection with <2 °C global warming

Nearly 200 countries have pledged to conserve 30% of terrestrial ecosystems to stop the global biodiversity crisis. However, biodiversity is not uniformly distributed across countries. Adequately addressing this crisis requires a scientific basis for selecting protected land that considers both ecological benefits and impacts to humans. We use the global change analysis model to evaluate land use tradeoffs of four land protection cases under two climate cases. We find that biodiversity-specific land protection up to 39% globally can reduce land use constraints and food prices compared to protecting 30% of land uniformly in each country (’30 × 30’ initiative). Valuing terrestrial carbon for climate change mitigation reduces land conversion pressure and can complement protection strategies. Global impacts to agriculture of additional land protection are small, but regional impacts vary and may be considerable. Overall, biodiversity-specific land protection has greater potential than a uniform target to meet both ecological and human needs.

biodiversity

Wildfire management decisions outweigh mechanical treatment as the keystone to forest landscape adaptation

Modern land management faces unprecedented uncertainty regarding future climates, novel disturbance regimes, and unanticipated ecological feedbacks. Mitigating this uncertainty requires a cohesive landscape management strategy that utilizes multiple methods to optimize benefits while hedging risks amidst uncertain futures. We used a process-based landscape simulation model (LANDIS-II) to forecast forest management, growth, climate effects, and future wildfire dynamics, and we distilled results using a decision support tool allowing us to examine tradeoffs between alternative management strategies. We developed plausible future management scenarios based on factorial combinations of restoration-oriented thinning prescriptions, prescribed fire, and wildland fire use. Results were assessed continuously for a 100-year simulation period, which provided a unique assessment of tradeoffs and benefits among seven primary topics representing social, ecological, and economic aspects of resilience. Projected climatic changes had a substantial impact on modeled wildfire activity. In the Wildfire Only scenario (no treatments, but including active wildfire and climate change), we observed an upwards inflection point in area burned around mid-century (2060) that had detrimental impacts on total landscape carbon storage. While simulated mechanical treatments (~ 3% area per year) reduced the incidence of high-severity fire, it did not eliminate this inflection completely. Scenarios involving wildland fire use resulted in greater reductions in high-severity fire and a more linear trend in cumulative area burned. Mechanical treatments were beneficial for subtopics under the economic topic given their positive financial return on investment, while wildland fire use scenarios were better for ecological subtopics, primarily due to a greater reduction in high-severity fire. Benefits among the social subtopics were mixed, reflecting the inevitability of tradeoffs in landscapes that we rely on for diverse and countervailing ecosystem services. This study provides evidence that optimal future scenarios will involve a mix of active and passive management strategies, allowing different management tactics to coexist within and among ownerships classes. Our results also emphasize the importance of wildfire management decisions as central to building more robust and resilient future landscapes.

54 ENVIRONMENTAL SCIENCES