Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sedimentary”

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

The LAKE model input dataset for three Arctic lakes

This dataset contains meteorological data collected for three Arctic lakes and compiled to satisfy input requirements of the LAKE 2.0 model. The dataset was generated to act as a benchmarking dataset for future model-data inter-comparisons. The LAKE 2.0 model simulates temperatures within the water later and the sedimentary layer of a lake. The LAKE2.0. is an open-source code and available to download via this weblike http://tesla.parallel.ru/Viktor/LAKE/-/wikis/LAKE-model (last visit July 14, 2021). The meteorological data are required to simulate the surface energy balance at the surface of a lake. This dataset includes a compilation of the meteorological data pulled from multiple data streams, including National Oceanic and Atmospheric Administration (NOAA) climate data, Circumarctic Lakes Observation Network (CALON) data, and the United States Geological Survey (USGS) data. The data were compiled for three Arctic lakes: FoxDen (66.55877, -164.45670), Atqasuk (70.452497, -156.951984), and Toolik (68.63150, -149.60740). Each meteorological data is in comma-delimited format (file extension ‘.dat’) and includes eight columns: Temperature [K], Pressure [Pa], longwave downward radiation [W/m2], shortwave downward radiation [W/m2], “U” wind speed [m/s], ”V” wind speed [m/s], humidity [kg/kg], precipitation [m/s]. In addition to the meteorological data file, we included setup and driver files. The Toolik lake is the deepest out of three lakes and has inflowing and outflowing groundwater data. InflowOutflowREADME.txt has more information about inflow and outflow flies. The other two lakes are much shallower and modeled as a closed system (i.e. no water inflow or outflow).

54 ENVIRONMENTAL SCIENCES↗

AmeriFlux US-Jo1 Jornada Experimental Range Bajada Site

This is the AmeriFlux version of the carbon flux data for the site US-Jo1 Jornada Experimental Range Bajada Site. Site Description - The Jornada Basin Experimental Range (JER) covers 783 km2 in the La Jornada del Muerto Plain of the northern Chihuahuan Desert and is located 20km of Las Cruces, NM. Extensive livestock grazing at the JER and througout the US Southwest was coincident with large-scale grassland deterioration and transition to shrubland begining in the 1800s. The JER was established n 1912 to investigate these rangeland changes and has since become a central location for understanding dryland ecology. This flux tower monitors CO2 and H2O dynamics in a representative shrubland on the piedmont slope (bajada) of the San Andreas mountains. The dominant shrubs are evergreen Larrea tridentata (Creosote) and winter-deciduous Prosopis glandulosa (Honey Mesquite). Other cover types include Flourensia cernua (tarbush) and patchy occurences of the grasses Muhlenbergia porteri (Bush Muhly) and Dasyochloa pulchella (Fluff Grass). The site is occasionally visited by stray domestic cattle, free-ranging introduced Oryx, and other native herbivores (Jack Rabbits, Desert Pronghorn). Soils at the site are Ustic Calciargids and parent material consists of limestone, other sedimentary rock, and some igneous rock. Virtual Site Visit: https://youtu.be/v1uJCKuicqs​

Tweedie, Craig↗

AmeriFlux FLUXNET-1F US-Jo1 Jornada Experimental Range Bajada Site

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Jo1 Jornada Experimental Range Bajada Site. This is the FLUXNET version of the carbon flux data for the site US-Jo1 Jornada Experimental Range Bajada Site produced by applying the standard ONEFlux (1F) software. Site Description - The Jornada Basin Experimental Range (JER) covers 783 km2 in the La Jornada del Muerto Plain of the northern Chihuahuan Desert and is located 20km of Las Cruces, NM. Extensive livestock grazing at the JER and througout the US Southwest was coincident with large-scale grassland deterioration and transition to shrubland begining in the 1800s. The JER was established n 1912 to investigate these rangeland changes and has since become a central location for understanding dryland ecology. This flux tower monitors CO2 and H2O dynamics in a representative shrubland on the piedmont slope (bajada) of the San Andreas mountains. The dominant shrubs are evergreen Larrea tridentata (Creosote) and winter-deciduous Prosopis glandulosa (Honey Mesquite). Other cover types include Flourensia cernua (tarbush) and patchy occurences of the grasses Muhlenbergia porteri (Bush Muhly) and Dasyochloa pulchella (Fluff Grass). The site is occasionally visited by stray domestic cattle, free-ranging introduced Oryx, and other native herbivores (Jack Rabbits, Desert Pronghorn). Soils at the site are Ustic Calciargids and parent material consists of limestone, other sedimentary rock, and some igneous rock. Virtual Site Visit: https://youtu.be/v1uJCKuicqs​

Tweedie, Craig↗

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗

Rayleigh Wave Propagation in the Bighorn Mountains Region, Wyoming

Short-period Rayleigh waves, Rg , provide strong constraints on the depth of shallow seismic events and are of interest for monitoring small explosions. Characterizing the seismic sources that generate Rg requires an understanding of how shallow crustal structure affects Rayleigh wave propagation. Here, in support of these efforts, this study utilizes observed waveforms from small shallow explosions recorded on temporary seismic network deployments in the Bighorn region, Wyoming. We study regional near-surface geology by measuring changes in surface-wave amplitude and polarization during propagation through basins, foothills, and mountains. We develop additional insight by carrying out surface-wave eigenfunction analyses and numerical-wave simulations, which together reproduce many characteristics seen in the observed waveforms. Our results show how sedimentary basins in the Bighorn region allow for amplified prograde-polarized higher-mode and retrograde-polarized fundamental-mode Rayleigh waves, whereas adjacent mountains only support retrograde motion. These different modes provide distinct constraints on the Earth structure and source characteristics, potentially enabling targeted inversions in future studies. Finally, our findings provide insight into Rg propagation through complex near-surface geology, improving our understanding of shallow propagation and source effects that are relevant to explosion monitoring efforts.

58 GEOSCIENCES↗

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗

Investigation of Site Amplifications Using Ambient-Noise-Derived Shallow Velocity Structures Under a Dense Array in Oklahoma

The shear-wave velocity (V S ) structure plays an important role in characterizing site amplification. The Large-n Seismic Survey in Oklahoma (LASSO; 1820 stations) revealed large vertical ground-motion variability in a 25 km × 32 km area in northern Oklahoma. The LASSO array has a relatively simple and flat topography, typical in a sedimentary basin environment in the central United States. In this study, we use the dense array to investigate the velocity structure under the LASSO array and how vertical ground motions relate to the shallow-to-deep structures. We extract the fundamental-mode Rayleigh wave by cross-correlating one month of ambient noise (0.7–5 Hz). We use double-beamforming to measure the group and phase velocities and anisotropy. By jointly inverting the group and phase velocities, we obtain the V S structure. Here, we observe correlations between V S at depths of 0.1–1.5 km and vertical ground motions using sites on the stiffer Permian formations. The shallow Quaternary alluvium and terrace deposits can amplify vertical ground motions by a factor of 2–4.5 between 2 and 25 Hz and attenuate signals above 25 Hz. We use 1D V S profiles to simulate the SV-wave transfer functions. An average V S of 250 m/s in the upper 20–40 m may cause the observed amplification between 2 and 40 Hz. V S estimated by topographic slopes cannot predict the relative amplification. Our results highlight the large variability of site-dependent ground motion in a small local region and the importance of characterizing shallow structures to estimate seismic hazards. Small thickness variations of the shallow formation can significantly change the resonance amplitude and frequency, which likely reduces the coherency of the Rayleigh waves extracted from ambient noise.

Chang, Hilary [Massachusetts Inst. of Technology (↗

Distinguishing Unique Earthquakes with Overlapping Signals in Oklahoma

Abstract During routine operations monitoring Oklahoma earthquakes, we found that certain earthquakes occurred closely both in space and time and had overlapping phases at the recording stations. Through further scrutiny and analysis, we determined that rather than being distinctly different earthquakes, some of the earthquakes exhibited multiphase arrivals and longer than expected coda due to unique ray paths that encounter impedance contrasts such as at the sedimentary rock-basement. Of course, some of these events truly were distinct events, which we term overlapping earthquakes, for which perceived coda duration overlaps and obscures the phase arrivals of the second event due to the source proximity in both time and space. We detail our classification scheme to separate the local earthquakes in Oklahoma as single, overlapping earthquakes, or those associated with multiphase arrivals. We forward model seismic wave propagation in a 2D crustal model and develop a methodology that utilizes waveform correlation to distinguish phases from overlapping earthquakes to those from crustal reverberations. Duration analysis shows a more elongated duration, qualitatively similar to the duration produced by overlapping earthquakes, at the sites where multiphase arrivals are observed.

Geochemistry & Geophysics↗

Seismic Monitoring near Ithaca, New York, Reveals Nonuniform Distribution of Microseismicity in an Intraplate Region

Abstract Cornell University intends to use a deep direct-use geothermal system to heat its Ithaca, New York, campus. In preparation for this project, the Cornell Seismic Network has been monitoring the background seismicity in this intraplate region since 2019. From January 2020 to June 2023, 95 events were detected within 20 km of the proposed geothermal well site, with local magnitudes ranging from −1.02 to 0.56. None of these events appear in regional or national catalogs. Events locate in a narrow geographic band, with one-fourth exhibiting multimodal hypocentral probability peaks both near the surface and at 1–4 km depth. We relocate events with a joint hypocenter and 1D velocity model inversion, in addition to a fully nonlinear method, and then compare observations with synthetic waveforms. Together, these approaches provide strong evidence for >95% of events locating at the surface or within the 3-km-thick sedimentary sequence. We explore how anthropogenic activity and regional topographic stress may contribute to frequent surficial events. This information is critical for characterizing the background microseismicity for comparison during future geothermal operations. Ithaca’s geology of Paleozoic sediments overlying Precambrian crystalline basement is typical of many continental interiors, so these results also provide insight into intraplate microseismicity patterns.

Geochemistry & Geophysics↗

AK112: Full Waveform Inversion Tomography of Alaska Improves Waveform Fits While Imaging Crustal, Mantle, and Slab Structure

We report a full waveform inversion tomography model of Alaska and the surrounding regions, inferring radially anisotropic shear and isotropic compressional wavespeeds by fitting complete waveforms from 120 regional earthquakes. Our multiscale approach inverted time–frequency phase misfits (maximum period of 100 s), starting with a minimum period of 40 s and ending at 20 s in 7 stages and 112 total iterations. The model (AK112) was evaluated by computing the misfits for 36 independent validation events. We find that misfit reductions were large and equal (∼55%) for both the inversion and validation data sets, providing confidence in the model. AK112 also provides much better waveform fits compared to other reported models for the region, including an isotropic version of itself, highlighting the importance of anisotropy. The model resolves known crustal, upper mantle, and slab structure to depths of 100 km with new detail: sedimentary basins in the Alaskan Shelf, Cook Inlet, and Colville basins, among others; discontinuous lithospheric structure across major terrane boundaries; and subducting slab geometry and back‐arc volcanic sources. In addition to tectonic interpretations, the model enables full waveform simulations for long‐period earthquake ground motions and source characterization (e.g., moment tensor and finite‐fault inversion).

Rodgers, Arthur [Lawrence Livermore National Labor↗

URC Assessment Method

This tool evaluates the potential occurrence of URC resources using a series of validated heuristics as outlined in "Creason, C.G., Justman, D., Rose, K. et al. A Geo-Data Science Method for Assessing Unconventional Rare-Earth Element Resources in Sedimentary Systems. Nat Resour Res (2023). https://doi.org/10.1007/s11053-023-10163-x"

Critical Minerals↗

Basin-Scale Structural Features Database

The Basin-Scale Structural Features database provides spatial datasets of faults, fractures, folds, and earthquakes compiled from public, authoritative sources (e.g., U.S. Geological Survey and State Geological Surveys) and aggregated into derivative forms to support subsurface assessments. Recognizing that characterizing basin-scale structural features requires interpreting data that are often ambiguous or lack key information, the source data were evaluated using a knowledge-data framework and geospatial fuzzy logic method (Justman et al., 2020) to represent both measured (observed) and predicted (inferred or potential) structural features as derivative datasets. This workflow employs conceptual models for known structural features and predicted structural features, incorporating geospatial data to estimate potential, even with limited data. The aim is to aid and support an understanding of basin-scale features and identify potential gaps in data and knowledge. As of 4/30/2025, the database includes resources for nine sedimentary basins: Appalachian, Denver, U.S. Gulf Coast, Illinois, Michigan, Permian, Sacramento, San Joquin and Williston. The database is organized by basin and then data category: 1) Faults, fractures, folds, 2) Earthquakes, 3) Topographic, 4) Structural contours and isopachs, 5) Geophysical, and 6) Structural feature density assessment maps.

basin scale↗

Discrete Measurements of the Least Horizontal Principal Stress from Core Data: An Application of Viscoelastic Stress Relaxation

Summary The emerging Paradox Oil Play in southeastern Utah is among the most significant unconventional plays in the western USA. The mean total undiscovered oil resources within just the Pennsylvanian Cane Creek interval of the Paradox Basin are believed to exceed 215 million barrels. However, to date, less than 5% (~9 million barrels) of the total Cane Creek resource has been produced from fewer than 40 wells, and only approximately one-half of those are horizontal wells. More than 95% of production is from the central Cane Creek Unit (CCU). Natural fractures are a key feature of many production wells, but stimulation by induced hydraulic fractures is not consistently successful. We hypothesize that more effective production in this play will rely on better fundamental characterization, especially on better quantification of the state of stress. Approximately 110 ft of core, well logs, and a diagnostic fracture injection test (DFIT) were acquired from the State 16-2 well within the CCU. With these data, we applied two methods to constrain and clarify the state of stress. The first technique, the Simpson’s coefficient method, provides lower bounds on the two horizontal principal stresses and relies on only limited data. Alternatively, the viscoelastic stress relaxation (VSR) method is used to estimate the least horizontal principal stress, building on observations that principal stresses become more isotropic as the viscous behavior of a rock is more pronounced. Results of these two methods support the hypothesis that the state of stress in the CCU of the Paradox Basin is nearly lithostatic and isotropic. Other factors consistent with this hypothesis include high formation pore pressure, which tends to reduce the possible stress states by changing the frictional failure equilibrium; lack of induced fractures in the core, which should be present in the case of stress anisotropy; and interbedded halite layers, which given their high degree of ductility, probably lead to greater VSR for the entire sedimentary package.

Energy & Fuels↗

A simple method for obtaining heat capacity coefficients of minerals

Abstract Heat capacity data are unavailable or incomplete for many minerals at geologically relevant temperatures. Despite the availability of entropy and enthalpy values in numerous thermodynamic tables (even sometimes at elevated temperatures), there remains need for extrapolation beyond, or interpolation between, temperatures. This approach inevitably results in estimates for entropy and enthalpy values because the heat capacity coefficients required for optimal thermodynamic treatment are less frequently available. Here we propose a simple method for obtaining heat capacity coefficients of minerals. This method requires only the empirically measured temperature-specific heat capacity for calculation via a matrix algorithm. The system of equations solver is written in the Python computing language and has been made accessible in an online repository. Thermodynamically, the solution to a system of equations represents the heat capacity coefficients that satisfy the mineral-specific polynomial. Direct coefficient calculation will result in more robust thermodynamic data, which are not subject to fitting uncertainties. Using hematite as an example, this method provides results that are comparable to conventional means and is applicable to any solid material. Coefficients vary within the traditional large 950 K temperature interval, indicating that best results should instead utilize a smaller 400 K temperature interval. Examples of large-scale implications include the refinement of geothermal gradient estimation in rapidly subsiding sedimentary basins or metamorphic and hydrothermal evolution.

Geochemistry & Geophysics↗

Hero Carbonsafe Phase 2 Project in the Columbia River Basalt Group: Technical Program Overview

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Hermiston is located near the center of the Columbia River Basalt Group (CRBG), which is one of the largest basalt flows in the US. Basalt CO2 storage has potential advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. associated decreases in pressure and CO2 migration risks, 3. reduced long-term monitoring requirements with respect to plume tracking, 4. widespread geographic distribution and, 5. large storage potential due to thickness, porosity, and CO2 interactions with basalt. For locations such as the Pacific Northwest (PNW), Hawaii, Iceland, India and Japan, whose localities are isolated from large sedimentary basins offering conventional saline storage options, basalt may offer the only feasible option for local CO2 storage. However, mineralization/basalt storage still has many uncertainties, as there are limited field-scale assessments of CO2 storage in basalt. There are significant uncertainties hindering the effective implementation of carbon capture utilization and storage (CCUS) in basalt. These include the lack of proven storage capacities, challenges in methodologies for modeling the area of review in igneous formations, limited understanding of mineralization kinetics and timing, and uncertainties in injectivity. Additionally, the domestic availability of specialized services and drilling expertise is constrained, and existing CCUS permitting and regulatory frameworks, originally developed for conventional saline reservoirs, may not adequately address the unique requirements of basalt systems. HERO CarbonSAFE is designed to address major research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts of the CRBG, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts.

58 GEOSCIENCES↗

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Evaluation of Engineered Barrier Systems (FY19 Report)

This report describes research and development (R&D) activities conducted during fiscal year 2019 (FY19) specifically related to the Engineered Barrier System (EBS) R&D Work Package in the Spent Fuel and Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Eneregy (DOE). The R&D activities focus on understanding EBS component evolution and interactions within the EBS, as well as interactions between the host media and the EBS. A primary goal is to advance the development of process models that can be implemented directly within the Genreric Disposal System Analysis (GDSA) platform or that can contribute to the safety case in some manner such as building confidence, providing further insight into the processes being modeled, establishing better constraints on barrier performance, etc.The FY19 EBS activities involved not only modeling and analysis work, but experimental work as well. The report documents the FY19 progress made in seven different research areas as follows: (1) thermal analysis for the disposal of dual purpose canisters (DPCs) in sedimentary host rock using the semianalytical method, (2) tetravalent uranium solubility and speciation, (3) modeling of high temperature, thermal-hydrologic-mechanical-chemical (THMC) coupled processes, (4) integration of coupled thermalhydrologic- chemical (THC) model with GDSA using a Reduced-Order Model, (5) studying chemical controls on montmorillonite structure and swelling pressure, (6) transmission x-ray microscope for in-situ nanotomography of bentonite and shale, and (7) in-situ electrochemical testing of uranium dioxide under anoxic conditions. The R&D team consisted of subject matter experts from Sandia National Laboratories, Lawrence Berkeley National Laboratory (LBNL), Los Alamos National Laboratory (LANL), Pacific Northwest National Laboratory (PNNL), the Bureau de Recherches Géologiques et Minières (BRGM), the University of California Berkeley, and Mississippi State University. In addition, the EBS R&D work leverages international collaborations to ensure that the DOE program is active and abreast of the latest advances in nuclear waste disposal. For example, the FY19 work on modeling coupled THMC processes at high temperatures relied on the bentonite properties from the Full-scale Engineered Barrier EXperiment (FEBEX) Field Test conducted at the Grimsel Test Site in Switzerland. Overall, significant progress has been made in FY19 towards developing the modeling tools and experimental capabilities needed to investigate the performance of EBS materials and the associated interactions in the drift and the surrounding near-field environment under a variety of conditions including high temperature regimes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗