Evaluation of Subsurface Gypsum-rich Rock in the Context of Geologic Carbon Sequestration
A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements for the Degree of Master of Science
SEARCH · Engineering Papers
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.
A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements for the Degree of Master of Science
Here, this study investigates the pore structure evolution of the reaction zones in wellbore cement samples exposed to a CO 2 -rich solution in equilibrium with 17 MPa supercritical CO 2 over 14 days. Through advanced characterization methods of field emission SEM, Quantitative Evaluation of Minerals by Scanning Electron Microscopy (QEMSCAN), and micro-CT, a new mechanism of CO 2 -cement reaction involving filling of nanopores in the interior of cement by gypsum was revealed. Gypsum was formed by the liberation of SO 4 2− from ettringite (AFt) and monosulfate (AFm) caused by a decrease in pH. Based on these experimental observations, a new CO 2 -cement reaction model that incorporates four distinct reaction zones is developed. This model provides a comprehensive framework for understanding the spatial and temporal distribution of minerals in cement due to high pressure CO 2 —cement reactions. This study demonstrates that the major damage induced by high pressure CO 2 alteration occurs in the most exterior region of the cement. The interior region of the cement maintains its integrity due to nanopore filling by gypsum.
Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi-year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage-limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.
Abstract Observing the environment in the vast regions of Earth through remote sensing platforms provides the tools to measure ecological dynamics. The Arctic tundra biome, one of the largest inaccessible terrestrial biomes on Earth, requires remote sensing across multiple spatial and temporal scales, from towers to satellites, particularly those equipped for imaging spectroscopy (IS). We describe a rationale for using IS derived from advances in our understanding of Arctic tundra vegetation communities and their interaction with the environment. To best leverage ongoing and forthcoming IS resources, including National Aeronautics and Space Administration’s Surface Biology and Geology mission, we identify a series of opportunities and challenges based on intrinsic spectral dimensionality analysis and a review of current data and literature that illustrates the unique attributes of the Arctic tundra biome. These opportunities and challenges include thematic vegetation mapping, complicated by low‐stature plants and very fine‐scale surface composition heterogeneity; development of scalable algorithms for retrieval of canopy and leaf traits; nuanced variation in vegetation growth and composition that complicates detection of long‐term trends; and rapid phenological changes across brief growing seasons that may go undetected due to low revisit frequency or be obscured by snow cover and clouds. We recommend improvements to future field campaigns and satellite missions, advocating for research that combines multi‐scale spectroscopy, from lab studies to satellites that enable frequent and continuous long‐term monitoring, to inform statistical and biophysical approaches to model vegetation dynamics.
Between 2017 and 2019, a series of injection tests were performed at pilot vertical well 58-32 at the Utah FORGE site. These injection activities were conducted in three distinct zones, isolated with packers and bridge plugs during the tests. Zone 1 consists of an open-hole section at the toe of the well, whereas Zones 2 and 3 are three-meter long cased and perforated sections. Several injection cycles were conducted in each zone. A detailed description of the geology of the FORGE site, of well 58-32 and of all injection cycles can be found in (Xing et al., 2020). The injection cycles included Diagnostic Fracture Injection Tests (DFITs). DFITs comprise two primary stages: an injection (pump-in) phase, where controlled fluid injection is employed to create a small hydraulic fracture in the surrounding rock, followed by a shut-in phase, where injection ceases, allowing observation of the formation’s pressure response. Utilizing methodologies such as the G-function analysis, DFITs provide valuable insights, including fracture closure pressure and the Instantaneous Shut-in Pressure (ISIP) and can be employed to estimate the minimum horizontal stress.
Abstract Understanding the influence of environmental factors on soil organic carbon (SOC) is critical for quantifying and reducing the uncertainty in carbon climate feedback projections under changing environmental conditions. We explored the effect of climatic variables, land cover types, topographic attributes, soil types and bedrock geology on SOC stocks of top 1 m depth across conterminous United States (US) ecoregions. Using 4559 soil profile observations and high-resolution data of environmental factors, we identified dominant environmental controllers of SOC stocks in 21 US ecoregions using geographically weighted regression. We used projected climatic data of SSP126 and SSP585 scenarios from GFDL-ESM 4 Earth System Model of Coupled Model Intercomparison Project phase 6 to predict SOC stock changes across continental US between 2030 and 2100. Both baseline and predicted changes in SOC stocks were compared with SOC stocks represented in GFDL-ESM4 projections. Among 56 environmental predictors, we found 12 as dominant controllers across all ecoregions. The adjusted geospatial model with the 12 environmental controllers showed an R 2 of 0.48 in testing dataset. Higher precipitation and lower temperatures were associated with higher levels of SOC stocks in majority of ecoregions. Changes in land cover types (vegetation properties) was important in drier ecosystem as North American deserts, whereas soil types and topography were more important in American prairies. Wetlands of the Everglades was highly sensitive to projected temperature changes. The SOC stocks did not change under SSP126 until 2100, however SOC stocks decreased up to 21% under SSP585. Our results, based on environmental controllers of SOC stocks, help to predict impacts of changing environmental conditions on SOC stocks more reliably and may reduce uncertainties found in both, geospatial and Earth System Models. In addition, the description of different environmental controllers for US ecoregions can help to describe the scope and importance of global and local models.
This is a presentation on the Evolution of Permeability and Strength Recovery of Shear Fractures Under Hydrothermal Conditions by United States Geological Survey, presented by Tamara Jeppson. This video slide presentation, by the USGS, discusses the determination of how thermal, hydrologic, mechanical, and chemical (THMC) processes affect the sustainability of fracture networks in geothermal reservoirs. This includes (1) the qualification of rates of change of fracture properties, (2) the parameterization of modes of reaction, (3) the development of micromechanical and empirical fracture models, and (4) extended THMC models for laboratory- and reservoir-scale models. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 15, 2024.
As directed in the Consolidation Appropriations Act, 2016, the National Nuclear Security Administration (NNSA) initiated the preconceptual design and development of a Lifecycle Cost Estimate for the Surplus Plutonium Disposition (SPD) Dilute and Dispose Program. Based on August 2016 Program Requirements Document and subsequent supplemental guidance, LANL prepared the Lifecycle Cost Estimate under key assumptions that meet the program’s requirements. For Dilute and Dispose, the program would disposition surplus Pu by diluting oxide produced at LANL with inhibitor materials, packaging the materials in containers, and shipping the containers to a deep geologic repository for permanent disposal. The base assumption is that LANL would disassemble pits, convert the Pu metal to oxide, and characterize and package the material for shipment to SRS, where it would be diluted prior to geologic disposal at the WIPP site in New Mexico. Another major assumption for the Dilute and Dispose option is that LANL would increase the current oxide production rate (or throughput) to 1500 kg/year, 5 times higher than the maximum annual production of ~300 Kgs executed by the ARIES Oxide Production Program at LANL. Analysis based on the ARIES program’s throughput model revealed that 15 pieces of equipment would need to be installed within PF-4 and certain facility improvements would need to be accomplished in order to meet the desired throughput levels. The additional equipment would be essentially identical to equipment already used within PF-4 for existing operations.
Abstract. In the western United States, prolonged drought, a warming climate, and historical fuel buildup have contributed to larger and more intense wildfires as well as to longer fire seasons. As these costly wildfires become more common, new tools and methods are essential for improving our understanding of the evolution of fires and how extreme weather conditions, including heat waves, windstorms, droughts, and varying levels of active-fire suppression, influence fire spread. Here, we develop the Geostationary Operational Environmental Satellites (GOES)-Observed Fire Event Representation (GOFER) algorithm to derive the hourly fire progression of large wildfires and create a product of hourly fire perimeters, active-fire lines, and fire spread rates. Using GOES-East and GOES-West geostationary satellite detections of active fires, we test the GOFER algorithm on 28 large wildfires in California from 2019 to 2021. The GOFER algorithm includes parameter optimizations for defining the burned-to-unburned boundary and correcting for the parallax effect from elevated terrain. We evaluate GOFER perimeters using 12 h data from the Visible Infrared Imaging Radiometer Suite (VIIRS)-derived Fire Event Data Suite (FEDS) and final fire perimeters from the California's Fire and Resource Assessment Program (FRAP). Although the GOES imagery used to derive GOFER has a coarser resolution (2 km at the Equator), the final fire perimeters from GOFER correspond reasonably well to those obtained from FRAP, with a mean Intersection-over-Union (IoU) of 0.77, in comparison to 0.83 between FEDS and FRAP; the IoU indicates the area of overlap over the area of the union relative to the reference perimeters, in which 0 is no agreement and 1 is perfect agreement. GOFER fills a key temporal gap present in other fire tracking products that rely on low-Earth-orbit imagery, where perimeters are available at intervals of 12 h or longer or at ad hoc intervals from aircraft overflights. This is particularly relevant when a fire spreads rapidly, such as at maximum hourly spread rates of over 5 km h−1. Our GOFER algorithm for deriving the hourly fire progression using GOES can be applied to large wildfires across North and South America and reveals considerable variability in the rates of fire spread on diurnal timescales. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates. The resulting GOFER product has a broad set of potential applications, including the development of predictive models for fire spread and the improvement of atmospheric transport models for surface smoke estimates (https://doi.org/10.5281/zenodo.8327264, Liu et al., 2023).
Abstract. Gas-phase oxidation of isoprene by ozone (O3) and the hydroxyl (OH) and nitrate (NO3) radicals significantly impacts tropospheric oxidant levels and secondary organic aerosol formation. The most comprehensive and up-to-date chemical mechanism for isoprene oxidation consists of several hundred species and over 800 reactions. Therefore, the computational expense of including the entire mechanism in large-scale atmospheric chemical transport models is usually prohibitive, and most models employ reduced isoprene mechanisms ranging in size from ∼ 10 to ∼ 200 species. We have developed a new reduced isoprene oxidation mechanism using a directed-graph path-based automated model reduction approach, with minimal manual adjustment of the output mechanism. The approach takes as inputs a full isoprene oxidation mechanism, the environmental parameter space, and a list of priority species which are protected from elimination during the reduction process. Our reduced mechanism, AMORE-Isoprene (where AMORE stands for Automated Model Reduction), consists of 12 species which are unique to the isoprene mechanism as well as 22 reactions. We demonstrate its performance in a box model in comparison with experimental data from the literature and other current isoprene oxidation mechanisms. AMORE-Isoprene's performance with respect to predicting the time evolution of isoprene oxidation products, including isoprene epoxydiols (IEPOX) and formaldehyde, is favorable compared with other similarly sized mechanisms. When AMORE-Isoprene is included in the Community Regional Atmospheric Chemistry Multiphase Mechanism 1.0 (CRACMM1AMORE) in the Community Multiscale Air Quality Model (CMAQ, v5.3.3), O3 and formaldehyde agreement with Environmental Protection Agency (EPA) Air Quality System observations is improved. O3 bias is reduced by 3.4 ppb under daytime conditions for O3 concentrations over 50 ppb. Formaldehyde bias is reduced by 0.26 ppb on average for all formaldehyde measurements compared with the base CRACMM1. There was no significant change in computation time between CRACMM1AMORE and the base CRACMM. AMORE-Isoprene shows a 35 % improvement in agreement between simulated IEPOX concentrations and chamber data over the base CRACMM1 mechanism when compared in the Framework for 0-D Atmospheric Modeling (F0AM) box model framework. This work demonstrates a new highly reduced isoprene mechanism and shows the potential value of automated model reduction for complex reaction systems.
The primary objective of the EGS Collab Project sponsored by DOE is to increase the understanding needed to efficiently implement enhanced geothermal systems (EGS). One goal of the EGS Collab project is to create a collaborative research environment in which to study stimulation of crystalline rock at the 10 meter scale. Key to this effort is the collection of high quality data to allow comparison to numerical coupled process models in an effort to build confidence in the codes and modeling techniques used. In response to this, the EGS Collab team has created an underground test bed at the Sanford Underground Research Facility (SURF) in Lead SD at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently designing a second test bed aimed at investigating shear stimulation (Experiment 2). At the Experiment 1 location, we have characterized our host rock using laboratory testing and numerous field-based geophysical and geological techniques, and created a well-instrumented test bed to allow us to carefully monitor stimulation events and flow tests. In addition to the installed geophysical sensors, we have used tracer tests, differences in the ambient microbial communities at flow collection locations, and cold water injection to inform us about dynamic flow pathways. In Experiment 1, we have hydraulically stimulated the host rock in a number times at several locations in one well, creating new fractures that connect to existing fractures between the injection and production boreholes. We have performed long-term ambient and chilled water injection tests as an analog to EGS, and have monitored system changes resulting from these water injections through geophysical monitoring, flow and pressure measurements, tracer tests, and microbiology. Here, we summarize the tests performed, issues identified including poroelastic and thermoelastic effects, Joule-Thomson effects, restarting effects, indications of flow channeling, and the primary learnings from Experiment 1.
For geologic carbon storage, the ability to detect secondary CO 2 plumes—defined as those CO 2 plumes accumulating outside the intended storage reservoir—is fundamental to preventing unexpected CO 2 migration into groundwater resources and for risk and liability management. Understanding the sensitivity of various geophysical methods to secondary plumes is crucial for designing cost-effective monitoring schemes. We use several modeling scenarios to demonstrate the process of assessing sensitivities and detection thresholds of three primary geophysical techniques—surface seismic, borehole-to-surface electromagnetic (EM), and surface and borehole gravity—for early detection of secondary CO 2 plumes in the post-injection phase. While seismic reflection methods are often considered in monitoring strategies to track the evolution of CO 2 plumes, they are also the most expensive. Due to cost considerations, especially for long-term post-injection monitoring, other techniques complement seismic monitoring when designing an adaptive monitoring network. Borehole-to-surface EM or surface gravity surveys are feasible for time-lapse monitoring of deep secondary CO 2 plumes. Furthermore, these surveys could be carried at intervals defined by site-specific conditions. If time-lapse EM and/or gravity surveys detect any signal responses beyond the expected change, it would trigger a need for the higher resolution seismic survey.
Remote sensing (RS) scene classification has wide applications in the environmental monitoring and geological survey. In the real-world applications, the RS scene images taken by the satellite might have two scenarios: clear and cloudy environments. However, most of existing methods did not consider these two environments simultaneously. In this paper, we assume that the global and local features are iscriminative in either clear or cloudy environments. Many existing Convolution Neural Networks (CNN) based models have made excellent achievements in the image classification, however they somewhat ignored the global and local features in their network structure. In this paper, we propose a new CNN based network (named GLNet) with the Global Encoder and Local Encoder to extract the discriminative global and local features for the RS scene classification, where the constraints for inter-class dispersion and intra-class compactness are embedded in the GLNet training. The experimental results on two publicized RS scene classification datasets show that the proposed GLNet could achieve better performance based on many existing CNN backbones under both clear and cloudy environments.
The recently released Circularity Gap Report 2023 by Circle Economy Foundation found that the circularity score for the global economy is declining. This indicates that material extraction from virgin sources is climbing over earlier years. As per data released by US Geological Survey, material use in the US economy has increased exponentially with the expanding US economy over last century. EPA tracked municipal solid waste from 1960 to 2018 and found that 50% of the waste is destined for landfills. EPA estimates that US industry is responsible for 2.7 billion ton of solid non-hazardous waste annually in our mostly linear economy model. The circular economy framework aims at decoupling economic value generation from virgin materials extraction from nature. The linear model of material extraction and disposal at end of life is highly unsustainable. Manufacturing companies have realized this and in their commitments to sustainability are adopting ambitious waste reduction targets. Through Better Plants program US Department of Energy has established Waste Reduction Network where it is offering technical assistance to partners to achieve these ambitious waste reduction goals. Basic requirements of establishing a target include identifying a baseline, quantifying waste performance, and measuring progress over time. One problem faced by industry is the unstandardized metrics to quantify waste performance that may not be well suited to demonstrate progress. In this paper, we research methods traditionally used to measure waste performance and highlight advantages, gaps and limitations of each method. Suitability of the methods applicable to different manufacturing circumstances are also examined. Finally, we present a case study of a large manufacturer that faced inconsistencies in their tracked measurement metric. A solution was proposed to alter the methodology to enable more accurate waste performance tracking against a baseline.
Recent classification of El Nino-Southern Oscillation (ENSO) into two types, Eastern (EP) and Central (CP) events, has highlighted the importance of the central Pacific. We show here that the local correlation between ENSO subsurface temperatures (T sub ) in the upper 100-m and thermocline depth anomalies breaks down in the central equatorial Pacific, whereas T sub remains well correlated with sea surface height anomalies. This observed difference in the central equatorial Pacific is simulated by almost all climate models. It arises from a vertically slanted Tsub anomaly structure that is unique to the central equatorial Pacific. We show that this feature is an adiabatic response to wind-driving that is even present in a linear dynamic model, as long as the model has enough baroclinic modes to adequately represent the observed vertical complexity. Our findings have implications for better understanding the different importance of thermocline feedback in CP and EP events.
Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.
Apparent hydrological sensitivity (η a ), the change in the global mean precipitation per degree K of global surface warming, is a key aspect of the climate system's response to increasing CO 2 forcing. To determine whether ηa depends on the forcing amplitude we analyze idealized experiments over a broad range of abrupt CO 2 forcing, from 2× to 8× preindustrial values, with two distinct climate models. We find little change in η a between 2× and 4×CO 2 , and almost no change beyond 5×CO 2 . We validate this finding under transient CO 2 forcing at 1%-per-year, up to 8×CO 2 . We further corroborate this result by analyzing the 1%-per-year output of more than 15 CMIP5/6 models. Lastly, we examine the 1,000-year long LongrunMIP model output, and again find little change in η a . This wealth of results demonstrates that η a is a very weak function of CO 2 forcing.
Abstract After large-magnitude earthquakes, a crucial task for impact assessment is to rapidly and accurately estimate the ground shaking in the affected region. To satisfy real-time constraints, intensity measures are traditionally evaluated with empirical Ground Motion Models that can drastically limit the accuracy of the estimated values. As an alternative, here we present Machine Learning strategies trained on physics-based simulations that require similar evaluation times. We trained and validated the proposed Machine Learning-based Estimator for ground shaking maps with one of the largest existing datasets (<100M simulated seismograms) from CyberShake developed by the Southern California Earthquake Center covering the Los Angeles basin. For a well-tailored synthetic database, our predictions outperform empirical Ground Motion Models provided that the events considered are compatible with the training data. Using the proposed strategy we show significant error reductions not only for synthetic, but also for five real historical earthquakes, relative to empirical Ground Motion Models.