Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Population Dynamics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

DNA viruses from different stages of a wastewater treatment plant in southwest Ohio

Wastewater treatment plants (WWTPs) house diverse populations of microbial communities that are dynamic across treatment stages. Although viruses, especially non-pathogenic ones are relatively less studied compared to their bacterial counterparts, they are equally important players in WWTP microbiomes. In this study, we sampled influent, activated sludge, and effluent stages of a WWTP in southwest Ohio to sequence their viral population using multiple displacement amplification (MDA) and metagenomic approaches. We recovered 1003 viral genomes, the majority of which were ssDNA (single-stranded DNA) viruses that formed distinct clusters representing novel species and phylogenetic groups not closely related to known human gut or terrestrial viruses. Additionally, tracking viruses across treatment stages showed several persisted across multiple treatment stages. These results suggest that further studies are needed to understand the persistence of viral populations, particularly non-pathogenic ones, across different stages of WWTPs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Too dim, too bright, and just right: Systems analysis of the Chlamydomonas diurnal program under limiting and excess light

Photosynthetic organisms coordinate their metabolism and growth with diurnal light, which can range in intensity from limiting to excessive. Little is known about how light intensity impacts the diurnal program in Chlamydomonas reinhardtii, or how diurnal rhythms in gene expression and metabolism shape photoprotective responses at different times of day. To address these questions, we performed a systems analysis of synchronized Chlamydomonas populations acclimated to low, moderate, and high diurnal light. Transcriptomic and proteomic data revealed that the Chlamydomonas rhythmic gene expression program is resilient to limiting and excess light: genome-wide, waves of transcripts, and proteins peak at the same times in populations acclimated to stressful light intensities as in populations acclimated to moderate light. Yet, diurnal photoacclimation gives rise to hundreds of gene expression changes, even at night. Time course measurements of photosynthetic efficiency and pigments responsive to excess light showed that high light-acclimated cells partially overcome photodamage in the latter half of the day prior to cell division. Although gene expression and photodamage are dynamic over the diurnal cycle, Chlamydomonas populations acclimated to low and high diurnal light maintain altered photosystem abundance, thylakoid architecture, and non-photochemical quenching capacity through the night phase. This suggests that cells remember or anticipate the light intensities that they have typically encountered during the day. The integrated data constitute an excellent resource for understanding photoacclimation in eukaryotes under environmentally relevant conditions.

Dupuis, Sunnyjoy↗
DOE OSTI2026

Spatio-temporal dynamics of Hendra virus in Australia reveal stable maintenance of diverse viral clades among Pteropus bats

Hendra virus (HeV) was discovered in 1994 in Australia. Limited genomic data have hindered comprehensive understanding of HeV’s evolutionary dynamics. Here, in this work, we recovered 48 HeV genomes from bats and 9 from horses from Australia between 2016 and 2020, revealing four distinct clades. Each clade was distributed over a large spatial area with multiple clades co-circulating within a single bat roost on the same day and over consecutive years. The diversity and temporal stability of co-circulating clades suggest that viral dynamics are driven by episodic shedding of existing lineages maintained at the population level, rather than immune-driven strain-replacement dynamics. HeV isolates of different clades displayed variation in phenotypic properties but minimal antigenic differences. We provide an overview of evolutionary dynamics, phenotypic properties and assessment of countermeasures for HeV, and provide insights into the processes that maintain virus diversity in bats and influence the potential for viral emergence.

genetic variation↗
DOE OSTI2026

Field-driven ion pairing dynamics in concentrated electrolytes

Here, we investigate ion pairing dynamics in electrolytes driven far from equilibrium using molecular simulations and nonequilibrium rate theory. Focusing on 0.5M LiPF 6 in water and acetonitrile under uniform electric fields, we compute transition path theory observables, including reactive fluxes and mean first-passage times of ion pairing. Moreover, we introduce a dynamical proxy of free-ion population, where its field-induced change is strongly correlated with the nonlinear enhancement of conductivity, yielding an increase of 40% at 50 mV/Å in acetonitrile, compared to that of less than 10% in aqueous electrolytes. Further kinetic analysis elucidates that Onsager’s classical theory substantially overestimates field-induced enhancement of ion pair dissociation in molecular electrolytes. This discrepancy arises from solvent-mediated dynamical pathways and field-induced dielectric decrement that suppress ion pair dissociation within explicit solvents, highlighting that a faithful description of molecular details is essential. Our results provide a molecular interpretation of nonlinear electrolyte transport beyond continuum theories and establish a general framework for quantifying nonequilibrium reaction kinetics in condensed phase systems.

Moon, Seokjin [University of California, Berkeley,↗
DOE OSTI2026

A Multi-Region SEIR Model Incorporating Inter-County Mobility and Time-Dependent Transmission Dynamics: Application to COVID-19 Disease Outbreak Data in North Carolina.

Classical infectious disease compartmental models typically do not incorporate spatial heterogeneity or mobility. We develop a multi-region susceptible-exposed-infected-recovered (SEIR) model in which disease dynamics are coupled via inter-region mobility and the transmission rate is both region and time dependent. We calibrate the model using rolling averages of daily COVID-19 data in all 100 North Carolina counties. Mobility parameters are prescribed using daily inter-county commuter data. The number of transmission rate parameters is substantially reduced by hypothesizing that the dynamics correlate with county-level population density. Parameter estimation is carried out using several objective functions with error terms at different scales. An additive combination of least squares error at the county-level and the state-level, along with a quadratic transmission rate polynomial, yields the lowest overall error at both spatial scales. The calibrated model is used to simulate regional effects of perturbing disease transmission rates in adjacent counties and to illustrate effects of the state’s mobility infrastructure on disease dynamics and spread for a new disease outbreak.

COVID-19 modeling↗
DOE OSTI2024

Tracking dissociation pathways of nitrobenzene via mega-electron-volt ultrafast electron diffraction

As the simplest nitroaromatic compound, nitrobenzene is an interesting model system to explore the rich photochemistry of nitroaromatic compounds. Previous investigations of nitrobenzene's photochemical dynamics have probed structural and electronic properties. These investigations paint, at times, a convoluted and sometimes contradictory description of the photochemical landscape. Here, we investigate the ultrafast dynamics of nitrobenzene triggered by photoexcitation at 267 nm for the first time using a structural probe with femtosecond time resolution. Our probe complements previous measurements of nitrobenzene's electronic structure evolution and aids in determining the photochemical dynamics with less ambiguity. We employ megaelectronvolt ultrafast electron diffraction to follow nitrobenzene's structural evolution within the first 5 ps after photoexcitation. We observe ground state recovery within 160 ± 60 fs through nonadiabatic dynamics. Based on comparisons of the experimental signal with molecular dynamics simulations, we exclude a significant population of the triplet manifold. Furthermore, we do not observe fragmentation of nitrobenzene within the investigated time window, which indicates that previously observed photofragmentation reactions take place in the vibrationally 'hot' ground state on timescales considerably beyond 5 ps.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Vortex-Controlled Quasiparticle Multiplication and Self-Growth Dynamics in Superconducting Resonators

Even in the quantum limit, non-equilibrium quasiparticle (QP) populations induce QP poisoning that irreversibly relaxes the quantum state and significantly degrades the coherence of transmon qubits. A particularly detrimental yet previously unexplored mechanism arises from QP multiplication facilitated by vortex trapping in superconducting quantum circuits, where a high-energy QP relaxes by breaking additional Cooper pairs and amplifying the QP population due to the locally reduced excitation gap and enhanced quantum confinement within the vortex core. Here we directly resolve this elusive QP multiplication process by revealing vortex-controlled QP self-generation in a highly nonequilibrium regime preceding the phonon bottleneck of QP relaxation. At sufficiently low fluence, femtosecond-resolved magneto-reflection spectroscopy directly reveals a continuously increasing QP population that is strongly dependent on magnetic-field-tuned vortex density and absent at higher excitation fluences. Quantitative analysis of the emergent QP pre-bottleneck dynamics further reveals that, although the phonon population saturates within $\simeq$10~ps, both free and trapped QPs continue to grow in a self-sustained manner--hallmarks of the long-anticipated QP-vortex interactions in nonequilibrium superconductivity. We estimate a substantial increase of $\sim$34% in QP density at vortex densities of $\sim$ 100 magnetic flux quanta per $\mathrm{μm^{2}}$. Our findings establish a powerful spectroscopic tool for uncovering QP multiplication and reveal vortex-assisted QP relaxation as a critical materials bottleneck whose mitigation will be essential for resolving QP poisoning and enhancing coherence in superconducting qubits.

Park, Joong M. [Ames Lab]↗
DOE OSTI2026

Surface Hydroxyls of Imogolite Nanotubes Drive Distinct Structures and Mobility Differentiation of Nanoconfined Water

Abstract Nanoconfined fluids, particularly water, govern subsurface geochemistry, yet the molecular-level mechanisms by which mineral surfaces dictate confined water structure and mobility remain poorly resolved. Here, we combine solution- and solid-state proton (1H) NMR spectroscopy, NMR relaxometry, modulated-gradient spin–echo (MGSE) NMR diffusometry, infrared spectroscopy, and molecular dynamics simulations to demonstrate that surface hydroxyls drive the structural and dynamic differentiation of water in imogolite nanotubes. In saturated suspensions, we observe the coexistence of distinct water 1H environments, each exhibiting significantly reduced mobility, which we attribute to strong interactions with the imogolite surfaces. As more mobile water is removed, long-range water diffusivity in the fibrous solid samples slows to 1.6 × 10–10 m2·s–1 while local fluctuations increase to 3.4 × 10–8 m2·s–1, indicating a significant increase in molecular restriction through surface interactions. Together, our experiments reveal three coexisting populations with varied structures and mobilities, all distinct from liquid bulk water. We resolve a persistent solid-like interfacial layer strongly bound to surface hydroxyls, characterized by an unusually short T2 (<1 ms), and a very close effective 1H–1H distance of ∼1.55 Å between water and the inner-surface silanol group. An inner-core population occupying the inner cavity exhibits an intermediate dynamic regime with restricted axial diffusion, while outer-surface water associated with aluminum hydroxyls retains relatively higher mobility. These results experimentally substantiate a hierarchy of water populations in imogolite and show how surface chemistry dictates the structures and dynamics of confined water.

Fleming, Xander B. [Earth and Environmental Scienc↗
DOE OSTI2025

DELVE-ing into the Milky Way’s Globular Clusters: Assessing Extratidal Features in NGC 5897, NGC 7492, and Testing Detectability with Deeper Photometry

Extratidal features around globular clusters (GCs) are tracers of their disruption, stellar stream formation, and their host’s gravitational potential. However, these features remain challenging to detect due to their low surface brightness. We conduct a systematic search for such features around 19 GCs in the DECam Local Volume Exploration (DELVE) survey Data Release 2, discovering a new extra-tidal envelope around NGC 5897 and find tentative evidence for an extended envelope surrounding NGC 7492. Through a combination of dynamical modeling and analyzing synthetic stellar populations, we demonstrate these envelopes may have formed through tidal disruption. We use these models to explore the detectability of these features in the upcoming Legacy Survey of Space and Time (LSST), finding that while LSST’s deeper photometry will enhance detection significance, additional methods for foreground removal like proper motions or metallicities may be important for robust stream detection. Our results both add to the sample of globular clusters with extratidal features and provide insights on interpreting similar features in current and upcoming data.

Chiti, A. [Univ. of Chicago, IL (United States); S↗
DOE OSTI2026

Synthesizing land use and demographic change in Southeast Asia’s smaller urbanized areas from 2000–2015

The majority of the human population now reside in urban areas today. The United Nations estimates that nearly half of all urban dwellers currently live in cities smaller than 500 000 persons and the majority of future urban growth will take place in Asia and Africa, likely in these smaller urban areas, not mega cities. Thus, understanding the factors that influence urban demographic trajectories in small urban areas is critical to address sustainable and equitable policy initiatives related to food security, changing climate hazard exposure, and economic opportunities. Here we focus on Southeast Asia—a region historically characterized by lower urban population proportions, yet with a rapidly shifting dynamic demographic—to examine correlates of demographic change among smaller cities. We combine two open-source satellite-informed datasets: GHS urban center database (2015) and age-sex gridded data from WorldPop to calculate socio-demographic characteristics to model drivers of change in annualized urban population growth from 2000–2015 for 505 urbanized places. We find a general pattern of decreasing dependency ratios as city-size increases for most urban areas in Southeast Asia. Higher rates of growth and more variation is observed for smaller cities—those with fewer than 300 000 persons, the lowest population limit for UN data on urbanization. When examining covariates of urban population growth, we find significant statistical associations of population change in smaller urbanized areas with climatic, economic, and land cover/land use variables, but with country-specific variations. Characterizing a continuum of urban population development in the context of changing environmental, economic and climate conditions has been an important sustainable development and equity issue for decades, but newer analysis of city-level drivers allows for systematic inquiry thus moving beyond total population counts for policy-relevant insight.

Southeast Asia synthesis↗
DOE OSTI2025

Integrated multi-omic characterizations of the synapse reveal RNA processing factors and ubiquitin ligases associated with neurodevelopmental disorders

The molecular composition of the excitatory synapse is incompletely defined due to its dynamic nature across developmental stages and neuronal populations. To address this gap, we apply proteomic mass spectrometry to characterize the synapse in multiple biological models including the fetal human brain and hiPSC-derived neurons. To prioritize the identified proteins, we develop an orthogonal multi-omic screen of genomic, transcriptomic, interactomic, and structural data. This data-driven framework identifies proteins with key molecular features intrinsic to the synapse, including characteristic patterns of biophysical interactions and cross-tissue expression. The multi-omic analysis captures synaptic proteins across developmental stages and experimental systems, including 493 synaptic candidates supported by proteomics. We further investigate three such proteins that are associated with neurodevelopmental disorders – the CUL3 E3 ubiquitin ligase, the DDX3X and YBX1 nucleic-acid binding proteins – by mapping their networks of physically interacting synapse proteins or transcripts. Our study demonstrates the potential of an integrated multi-omic approach to systematically and more comprehensively resolve the synaptic architecture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Polar primary aerosols across the ocean-sea ice-snow-atmosphere interface: From sources to impacts

Primary aerosols play a critical role in polar climate systems, influencing cloud formation, precipitation, radiative balance, and surface energy budgets. This paper provides a comprehensive synthesis of primary aerosol sources, transformation and removal processes, and broader atmospheric impacts in polar regions, emphasizing their links to ocean and sea ice biogeochemistry. These aerosols (including sea salt, primary organic aerosol, and primary biological aerosol particles) originate from marine and cryospheric environments and are emitted through physical processes, such as wave breaking, bubble bursting, and blowing snow. Emission sources include seawater, sea ice, snow, and freshwater from river discharge and glacial runoff. Once airborne, these particles can serve as a chemical reservoir, influencing atmospheric composition and reactivity, and as seeds for cloud droplet and ice crystal formation, influencing cloud microphysics and polar climate. Despite their importance, many of the processes governing primary aerosol emissions and transformations remain poorly constrained. The most pressing knowledge gaps pertain to emission processes, limited spatiotemporal observational coverage, instrumentation constraints, parameterization development, and the integration of interdisciplinary expertise. To improve our understanding of primary aerosol drivers and their response to climate, future research efforts should prioritize strategically coordinated and cross-disciplinary process studies, advancements in measurement technologies and coverage, and close collaboration between modelers and observational scientists to inform and refine model parameterizations. As polar regions continue to undergo profound changes marked by increased precipitation, reduced sea and land ice, freshening oceans, and shifting ecosystem dynamics, characterizing present-day primary aerosol populations is vital. Improved understanding will be essential for anticipating future changes in aerosol-radiation and aerosol-cloud interactions and their implications for polar and global climate systems.

Aerosol-cloud↗
DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Quantitative stable isotope probing (qSIP)-informed metagenomics identifies viruses infecting chemoautotrophs

Aquatic environments absorb ~2.5 gigatonnes of atmospheric carbon each year1, more than the carbon stored in the atmosphere, soils, and all biomass combined. Primary producers transform this dissolved inorganic carbon into biomass that can subsequently flow into other trophic levels, or be released back into the environment through viral lysis. While there is substantial knowledge about the diversity and activity of viruses infecting photoautotrophic primary producers and the ecosystem impact, little is known about viruses infecting chemoautotrophs, representing a gap in our understanding of key processes driving microbial carbon cycling. Here, we combine metagenomics with quantitative 12/13C stable isotopic probing (qSIP) mesocosm experiments in a marine-derived meromictic pond to quantify population-specific isotopic enrichment, identify key chemoautotrophic primary producers, and virus-host dynamics. Isotopically enriched carbon is tracked from the genomes of chemoautotrophs to putative viruses, showing that active populations of hydrogen/sulfur-oxidizing chemoautotrophs (Thiomicrorhabdus, Hydrogenovibrio, Sulfurimonas, Sulfurovum) are targeted by viruses. This work provides the foundation for revealing the diversity and activity of viruses infecting globally-widespread chemoautotrophs. Our study sheds light on trophic interactions that impact microbial carbon cycling in aphotic environments and builds toward biogeochemical models that incorporate viral impacts on chemoautotrophic microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

SAXS of murine amelogenin identifies a persistent dimeric species from pH 5.0 to 8.0

Amelogenin is an intrinsically disordered protein essential to tooth enamel formation in mammals. Here, using advanced small angle X-ray scattering (SAXS) capabilities at synchrotrons and computational models, we revisited measuring the quaternary structure of murine amelogenin as a function of pH and phosphorylation at serine-16. The SAXS data shows that at the pH extremes, amelogenin exists as an extended monomer at pH 3.0 (R g = 38.4 Å) and nanospheres at pH 8.0 (R g = 84.0 Å), consistent with multiple previous observations. At pH 5.0 and above there was no evidence for a significant population of monomeric species. Instead, at pH 5.0 ~ 80% of the population is a heterogenous dimeric species that increases to ~ 100% at pH 5.5. The dimer population was observed at all pH > 5 conditions in dynamic equilibrium with a species in the pentamer range at pH < 6.5 and nanospheres at pH 8.0. At pH 8.0 ~ 40% of the amelogenin still remained in the dimeric state. In general, serine-16 phosphorylation of amelogenin appears to modestly stabilize the population of the dimeric species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Scalable Generation of High-fidelity Synthetic Population Ensembles

Used within social simulations, synthetic population ensembles enable uncertainty quantification (UQ) methods for obtaining more robust model inference and prediction. A synthetic population ensemble is a series of plausible virtual reconstructions of an area’s population at the granularity of people and residences, generated stochastically to preserve privacy of the source population survey’s respondents. In this paper, we demonstrate the production of large synthetic population ensembles for the U.S. via Oak Ridge National Laboratory’s UrbanPop framework to support modeling of high spatial resolution energy affordability metrics from nationwide social surveys in collaboration with the fusionACS project. The study involves two scenarios: creating ensembles for (1) 17 U.S. metropolitan areas in 2019 and (2) full U.S. Census Divisions in 2023, with each scenario consisting of 41 population instances (a base realization and 40 replicates). To accomplish this task at scale, we configured an integrated system within a research cloud, comprised of virtual containerizations, GPU-enhanced functionality, and orchestrated deployments of UrbanPop’s maturing Likeness Python ecosystem. Results demonstrate we maintained high-fidelity approximations of residential totals by areas of interest and the demographic characteristics of neighborhoods while reducing manual workflow burdens. Finally, we discuss plans to fine-tune and further develop our automated workflows for truly distributed job orchestration to increase computational efficiency, as well as provide an outlook for broadening applications of the ensembles.

Cluster computing↗
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗