Search NASA⌕ Search

SEARCH · Search NASA

Results for “Deep Space Transport”

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.

67 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment (BIRDIE) (Field Campaign Report)

The Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment or BIRDIE was a precipitation collection campaign designed to measure cosmogenic isotopes of beryllium in rainfall and other forms of precipitation across the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) atmospheric observatory. Precipitation samples were collected over three weeks (May 21-June 10, 2022) at 11 SGP sites across a 120-km 2 domain in northern and central Oklahoma. The field campaign targeted summertime convection rain events. We aimed to capture any spatial and/or temporal variability in the concentrations of beryllium-10 (10Be) and beryllium-7 (7Be) across the ARM SGP site during these rainfall events. The ratio of 10Be to 7Be (10/7Be) has been shown to be a powerful tracer for transport processes in the stratosphere-troposphere exchange (STE) system (e.g., Koch and Rind 1998). Any observed variability is analyzed to see if it could be linked to STE processes, including overshooting of deep convection clouds that penetrate the lower stratosphere. Six wet weather system events were collected within our three-week span. Five out of six of these rainfall events are classified as moderately to strongly convective. Those events are studied in more detail. The five systems also follow each other over a short period of 11 days. In our rainfall samples we observed spatial and event variability in the concentrations of beryllium-10 and -7 that resulted in variability of 10/7Be ratios. Each weather event appeared to exhibit a unique 10/7Be pattern across the domain that was not directly linked to storm total precipitation amount. To better understand the causes for this spatial variability, cloud characteristics in the domain were looked at to assess cloud type, cloud height, and depth information. We also ran Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model back trajectories to identify possible air mass and moisture origins across our events. Detailed analysis is ongoing. The purpose of this technical report is to provide background on using cosmogenic isotopes as atmospheric tracers, describe the experimental design, and present the beryllium concentration and ratio data. Lastly, we note that BIRDIE coincided with two other deep convection field campaigns in the area: The National Aeronautics and Space Administration (NASA)’s Dynamics and Chemistry of the Summer Stratosphere (DCOTSS) experiment (https://dcotss.org/) centered in Salina, Kansas, and the Targeted Observation by Radars and (uncrewed aerial systems) UAS of Supercells (TORUS) experiment (https://catalog.eol.ucar.edu/torus_2022) across the U.S. Central Plains. Additionally, a Lawrence Livermore National Laboratory ( LLNL) team at DCOTSS attempted to measure 10Be/7Be ratios on aerosols collected in situ in the lower and mid-stratosphere over Kansas using high-altitude balloons. Precipitation samples were also collected on the ground in Salina, Kansas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Processes in Salt Repositories for Radioactive Waste Disposal

This document summarizes the key processes (thermal, hydrological, mechanical, and chemical; THMC) impacting the features of a deep geological repository for radioactive waste in salt. Some processes are natural and on-going whether the repository is there or not, and other processes are driven by the perturbation associated with the repository. The features considered here include both engineered and natural components of the repository system. The engineered barrier system (EBS) in a salt repository is quite different from those implemented for a repository in clay or crystalline rocks, because it is comprised mostly of granular salt and salt-compatible cements, rather than bentonite. When compared to other rocks (i.e., silicates), salt has unique properties that make it an excellent potential host rock. Openings and fractures in salt creep closed readily. Salt has high thermal conductivity, which can reduce peak temperatures. Additionally, far away from the excavations the porosity of salt is unconnected, which leads to essentially zero advective or diffusive transport. The small amount of hypersaline brine occurring in salt minimizes microbial activity, reduces colloid-assisted transport, and eliminates in-package criticality (i.e., chloride is a neutron poison). At the end of the report, we present a brief outline for a potential salt repository, including considerations avoided in previous repository disposal concepts. We propose considering higher-temperature processes in future disposal concepts, rather than trying to minimize the thermal perturbation of the repository. Since hot salt is drier, a dry repository would limit corrosion, gas generation, and solute transport. Openings and fractures creep shut faster in hot salt. Therefore, higher temperatures could be seen as beneficial, rather than something to minimize, through increased spacing between waste packages (increasing repository costs).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Anatomy of Continuous Mars SEIS and Pressure Data from Unsupervised Learning

The seismic noise recorded by the Interior Exploration using Seismic Investigations, Geodesy, and Heat Transport (InSight) seismometer (Seismic Experiment for Interior Structure [SEIS]) has a strong daily quasi-periodicity and numerous transient microevents, associated mostly with an active Martian environment with wind bursts, pressure drops, in addition to thermally induced lander and instrument cracks. That noise is far from the Earth’s microseismic noise. Quantifying the importance of nonstochasticity and identifying these microevents is mandatory for improving continuous data quality and noise analysis techniques, including autocorrelation. Cataloging these events has so far been made with specific algorithms and operator’s visual inspection. We investigate here the continuous data with an unsupervised deep-learning approach built on a deep scattering network. This leads to the successful detection and clustering of these microevents as well as better determination of daily cycles associated with changes in the intensity and color of the background noise. We first provide a description of our approach, and then present the learned clusters followed by a study of their origin and associated physical phenomena. We show that the clustering is robust over several Martian days, showing distinct types of glitches that repeat at a rate of several tens per sol with stable time differences. We show that the clustering and detection efficiency for pressure drops and glitches is comparable to or better than manual or targeted detection techniques proposed to date, noticeably with an unsupervised approach. Finally, here we discuss the origin of other clusters found, especially glitch sequences with stable time offsets that might generate artifacts in autocorrelation analyses. We conclude with presenting the potential of unsupervised learning for long-term space mission operations, in particular, for geophysical and environmental observatories.

58 GEOSCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The ATTO-Campina site: A new observatory for tropical convection and gas-aerosol-cloud-precipitation interactions in the Amazon

We present results from ATTO-Campina, a new permanent observational site in central Amazon, about 4 km from the ATTO towers. Operational since 2020, ATTO-Campina characterizes atmospheric, cloud and rainfall properties through remote sensing. The goal is to provide continuous, complementary measurements to the ATTO towers, addressing the rainforest’s complex gas-aerosol-cloud-precipitation dynamics. Using a 3.5-year dataset, we classified convective clouds into three types: shallow cumulus (ShCu), congestus (Con) or (Deep) clouds. The shallow-to-deep transition takes about three hours, starting with ShCu formation at 11:00 local time. The accumulated rainfall peak follows at about 16:00. Only weak downdrafts are present in the upper troposphere where previous studies indicate new particle formation (NPF) occurrence. Strong downdrafts are mostly limited to heights below 5 km. Con and Deep convective days have higher concentrations of ultrafine aerosol and lower concentrations of accumulation-mode particles compared to ShCu. Convective clouds also significantly modify gas mixing ratios. Deep convective clouds are associated with high near-surface O3, consistent with downward transport from the midtroposphere. Our results showcase the added detail achieved by integrating data from the ATTO towers and ATTO-Campina sites. Together, these sites support better understanding of interconnected gas-aerosol-cloud-precipitation processes in the Amazon and their evolution under climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Comparative Study of Deep Learning Models for Fracture and Pore Space Segmentation in Synthetic Fractured Digital Rocks

This study focuses on the comparative study of deep learning (DL) models for pore space and discrete fracture networks (DFNs) segmentation in synthetic fractured digital rocks, specifically targeting low-permeability rock formations, such as shale and tight sandstones. Accurate characterization of pore space and DFNs is critical for subsequent property analysis and fluid flow modeling. Four DL models, SegNet, U-Net, U-Net-wide, and nested U-Net (i.e., U-Net++), were trained, validated, and tested using synthetic datasets, including input and label image pairs with varying properties. The model performance was assessed regarding pixel-wise metrics, including the F1 score and pixel-wise difference maps. In addition, the physics-based metrics were considered for further analysis, including sample porosity and absolute permeability. Particularly, We first simulated the permeability of porous media containing only pore space and then simulated the permeability of porous media with DFNs added. The difference between these two values is used to quantify the connectivity of segmented DFNs, which is an important parameter for low-permeability rocks. The pixel-wise metrics showed that the nested U-Net model outperformed the rest of the DL models in pore space and DFNs segmentation, with the SegNet model exhibiting the second-best performance. Particularly, nested U-Net enhanced segmentation accuracy for challenging boundary pixels affected by partial volume effects. The U-Net-wide model achieved improved accuracy compared to the U-Net model, which indicated the influence of parameter numbers. Similarly, nested U-Net has the closest match to the ground truth of physics-based metrics, including the porosity of pore space and DFNs, and the permeability difference quantifying the connectivity of DFNs. The findings highlight the effectiveness of DL models, especially the U-Net++ model with nested architecture and redesigned skip connections, in accurately segmenting pore spaces and DFNs, which are crucial for pore-scale fluid flow and transport simulation in low-permeability rocks.

Wang, Hongsheng↗
DOE OSTI2026

Short‐Term Hourly Weather Forecasting Using PredRNN With Image Preprocessing

Global weather forecast models are vital tools with numerous applications, including public safety, agriculture, and transportation. Recent advancements in artificial intelligence (AI) and deep learning (DL) have shown the potential to enhance weather forecasting accuracy and speed. In this study, we developed a short-term hourly weather forecast framework with a wavelet transform function for data preprocessing and a spatiotemporal DL model, PredRNN, for predicting five surface atmospheric variables, including wind speed and direction, mean sea level pressure (MSLP), temperature, and precipitation. The framework demonstrated promising results. It produces global forecasts at 0.25° (∼25 km) with a 1-day lead time RMSE of 1.8 m/s for wind components, 180 Pa for MSLP, and 1.8 K for temperature. Although our model does not surpass state-of-the-art AI weather forecast models across all metrics, it outperforms these models in precipitation forecasting and wind prediction at short lead times and achieves comparable accuracy for MSLP. Its native hourly forecasting capability, together with training on widely accessible GPU hardware, contributes meaningfully to the advancement of accessible DL weather forecasting methods. Our work highlights the importance of integrating temporal components and data transformation techniques to improve the predictability and accuracy of weather forecasts.

Tran, Hoang [Pacific Northwest National Laboratory↗
DOE OSTI2023

Converging towards atomic and nuclear pressures (Final Report)

his report details the development of transformational tools and techniques, allowing the accurate characterization of matter at the most extreme conditions yet studied in the high energy density (HED) domain. Before this work, most experiments quantitatively mapping the nature of matter in the HED realm were performed using planar geometry, which allows for the isolation and measurement of variables in a single thermodynamic state. Such measurements include equation of state variables, optical conductivity, heat transport and more. However, due to a variety of plasma processes, such experiments are limited to below ~10 TPa (100 Mbar) pressures. To achieve higher pressures, convergent experiments are needed. Historically, convergent experiments have not been used for benchmark data because they are integral measurements. That is each part of a convergent target undergoes a time dependent wide range of states, making it difficult to accurately isolate any particular quantity for a given thermodynamic state. This effort developed innovative convergent HED platforms and techniques enabling the exploration of matter from atomic to nuclear scale pressures. This effort also performed pioneering experiments that yielded the first rigorous benchmark data at these extreme conditions. This funding award began with the overarching goal to understand the behavior of matter at extreme (atomic-to-nuclear scale) pressures. The motivation for these goals lie in the fact that every time scientists explore matter beyond the threshold of an atomic unit, there is a fundamental shift in science. The atomic unit for energy, mass, charge, length, and time have each been explored, and each time such a threshold was crossed, a new sequence of discoveries was made resulting in significant awards such as the Nobel Prize. The only unexplored atomic unit is pressure, and this effort set the course for exploring matter at and beyond such pressures. That most of the recently discovered extrasolar planets and stars, as well as matter in the late implosion stages of inertial fusion targets, have deep internal pressures at and beyond atomic pressures amplifies the importance of this effort. Much of the initial work in developing techniques to create these atomic-to-nuclear pressures already existed within the HED community, predominantly at large scale laser facilities such as the National Ignition Facility (NIF) and the Omega60 laser the University of Rochester and although these experiments were routinely performed the ability to extract information about the underlying physical states and processes remained elusive due to the complexity of the experiments, extreme scales in both time and space, and the integrated nature of all the measurements. This work built a rigorous framework so such measurements can routinely be made. This was done in part through the introduction of Bayesian inference techniques into the field of HED science. These techniques allow for the self-consistent extraction of the relevant variables and their explicit and implicit correlations, so as to make use of integrated experimental data to constrain physical models and provide accurate uncertainty bounds for the data. The techniques developed here are fully transparent and proved immediately useful providing quantitative rigorous benchmarks for physical states and processes at some of the most extreme conditions yet explored on Earth. This funding is directly or partly responsible for 7 publications in major peer-reviewed scientific journals, a doctoral thesis, 3 invited talks at major conferences, and the training of a post-doc, 2 graduate students, and 1 undergraduate student. Beyond the effort supported by this award launched the introduction of Bayesian inference into the HED physics community and helped push the usage if modern data-science techniques within the physics community at large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗
DOE OSTI2020

Neutronics Analysis of Cold Critical KRUSTY Experiments using MCNP and Serpent

Kilowatt Reactor Using Stirling Technology (KRUSTY) is a prototype designed as a proof of concept for NASA’s Kilo power program, which was funded in order to create a small reactor for various space applications such as providing power for a colony on Mars or the moon and eventually as a possible option for powering rockets on deep space missions. KRUSTY was tested at the National Criticality Experiments Research Center (NCERC) located at the Nevada National Security Site (NNSS) from November 2017, through March 2018. The experimental campaign was conducted in four phases: com ponent critical experiments, cold critical experiments, warm critical runs, and high temperature demonstration. The component critical and cold critical experiments were used to determine the worth of the beryllium oxide (BeO) reflector rings and the boron carbide (B4C) control rod disks under different conditions. These conditions were altered by the addition and subtraction of other components in the assembly. The component critical phase was performed at atmospheric pressure and had four solid support rods spanning the center core; whereas the cold critical phase had the center core placed in a vacuum chamber and eight sodium-filled heat pipes spanning the center core region. This summary discusses the static neutronic analysis completed on the cold critical phase of the campaign using MCNP ® and Serpent. These two codes were used to assess any potential differences or biases in the simulation results. Understanding these potential biases is a key step toward the final goal of being able to fully analyze the transient experiments performed during phases 3 and 4 of the campaign.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Explicit physics-informed neural networks for nonlinear closure: The case of transport in tissues

In upscaling methods, closures for nonlinear problems present a well-known challenge. While a number of theoretical methods have been proposed for handling such closures, nonlinearities still remain a significant obstacle for many problems. In this work, we use a combination of formal upscaling and data-driven machine learning for explicitly closing a nonlinear transport and reaction process in multiscale tissues. The classical effectiveness factor model is used to formulate the macroscale reaction kinetics. We train a multilayer perceptron network using training data generated by direct numerical simulations over microscale examples. Once trained, the network is used in an algorithm for numerically solving the upscaled (coarse-grained) differential equation describing mass transport and reaction in two example tissues. The network is described as being explicit in the sense that the network is trained using macroscale concentrations and gradients of concentration as components of the feature space rather than incorporating them as part of a constraint in the optimization process. Network training and solutions to the macroscale transport equations were computed for two different tissues. The two tissue types (brain and liver) exhibit markedly different geometrical complexity and spatial scale (cell size and sample size). The upscaled solutions for the average concentration are compared with numerical solutions derived from the microscale concentration fields by a posteriori averaging. There are three outcomes of this work of particular note. 1) Our overall approach results in an upscaled nonlinear PDE. The PDE is closed using a neural network, and our approach results in the definition of the classical effectiveness factor for effecting closure. 2) We identify particular source terms for the closure problem that are important for representing the structure of the closure. These source terms involve macroscale concentrations and their gradients. We adopt these source terms to use as explicit features in the learning algorithm. We find the trained networks that include the macroscale source terms generate models that are able to predict the correction factor with increased fidelity over those that do not. 3) We find that the trained network exhibits good generalizability, and it is able to predict the effectiveness factor with high fidelity for realistically-structured tissues despite the significantly different scale and geometrical complexity of the two example tissue types. This latter result emphasizes our purposeful connection between conventional averaging methods with the use of machine learning for closure; this contrasts with some machine learning methods for upscaling where the exact form of the macroscale equation remains unknown.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

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↗