Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic characterization”

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 595 records · Page 33

Linking Dissolved Organic Matter Composition to Landscape Properties in Wetlands Across the United States of America

Abstract Wetlands are integral to the global carbon cycle, serving as both a source and a sink for organic carbon. Their potential for carbon storage will likely change in the coming decades in response to higher temperatures and variable precipitation patterns. We characterized the dissolved organic carbon (DOC) and dissolved organic matter (DOM) composition from 12 different wetland sites across the USA spanning gradients in climate, landcover, sampling depth, and hydroperiod for comparison to DOM in other inland waters. Using absorption spectroscopy, parallel factor analysis modeling, and ultra‐high resolution mass spectroscopy, we identified differences in DOM sourcing and processing by geographic site. Wetland DOM composition was driven primarily by differences in landcover where forested sites contained greater aromatic and oxygenated DOM content compared to grassland/herbaceous sites which were more aliphatic and enriched in N and S molecular formulae. Furthermore, surface and porewater DOM was also influenced by properties such as soil type, organic matter content, and precipitation. Surface water DOM was relatively enriched in oxygenated higher molecular weight formulae representing HUP High O/C compounds than porewaters, whose DOM composition suggests abiotic sulfurization from dissolved inorganic sulfide. Finally, we identified a group of persistent molecular formulae (3,489) present across all sites and sampling depths (i.e., the signature of wetland DOM) that are likely important for riverine‐to‐coastal DOM transport. As anthropogenic disturbances continue to impact temperate wetlands, this study highlights drivers of DOM composition fundamental for understanding how wetland organic carbon will change, and thus its role in biogeochemical cycling.

Environmental Sciences & Ecology↗

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗

The root-zone soil moisture spectrum in a mediterranean ecosystem

Storage of water within soil pores of the root zone introduce memory effects in the dynamics of soil moisture that are considerably longer than the integral timescale of many atmospheric processes. Thus, hydro-climatic states can be “sustained” through land-surface heat and water vapor fluxes primarily because they can “feed off” on this long-term soil moisture memory. Root-zone soil moisture memory is only but one feature characterizing the spectrum of soil moisture dynamics, which is analyzed here using a combination of long-term measurements and models. In particular, the spectrum of root-zone soil moisture content in a Mediterranean ecosystem is examined using 14-years of half-hourly measurements. A distinguishing hydro-climatic feature in such ecosystems is that sources (mainly rainfall) and sinks (mainly evapotranspiration) of soil moisture are roughly out of phase with each other. For over 4 decades of time scales and 7 decades of energy, the canonical shape of the measured soil moisture spectrum is shown to be approximately Lorentzian determined by the soil moisture variance and its memory but with two exceptions: the occurrences of a peak at diurnal-to-daily time scales and a weaker peak at near annual time scales. Model calculations and spectral analysis demonstrate that diurnal and seasonal variations in hydroclimate forcing responsible for variability in evapotranspiration had minor impact on the normalized shape of the soil moisture spectrum. However, their impact was captured by adjustments in the temporal variance. These findings indicate that precipitation and not evapotranspiration variability dominates the multi-scaling properties of soil moisture variability consistent with prior climate model simulations. Furthermore, the soil moisture memory inferred by the annual peak of soil moisture (340 d) is consistent with climate model simulations, while the memory evaluated from the loss function of a linearized mass balance approach leads to a smaller value (50 d), highlighting the effect of weak non-stationarity on soil moisture variability.

54 ENVIRONMENTAL SCIENCES↗

Efficient marine environmental characterisation to support monitoring of geological CO 2 storage

Carbon capture and storage is key for mitigating greenhouse gas emissions, and offshore geological formations provide vast CO 2 storage potential. Monitoring of sub-seabed CO 2 storage sites requires that anomalies signifying a loss of containment be detected, and if attributed to storage, quantified and their impact assessed. However, monitoring at or above the seabed is only useful if one can reliably differentiate abnormal signals from natural variability. Baseline acquisition is the default option for describing the natural state, however we argue that a comprehensive baseline assessment is likely expensive and time-bound, given the multi-decadal nature of CCS operations and the dynamic heterogeneity of the marine environment. We present an outline of the elements comprising an efficient marine environmental baseline to support offshore monitoring. We demonstrate that many of these elements can be derived from pre-existing and ongoing sources, not necessarily related to CCS project development. We argue that a sufficient baseline can be achieved by identifying key emergent properties of the system rather than assembling an extensive description of the physical, chemical and biological states. Further, that contemporary comparisons between impacted and non-impacted sites are likely to be as valuable as before and after comparisons. Furthermore, as these emergent properties may be nuanced between sites and seasons and comparative studies need to be validated by the careful choice of reference site, a site-specific understanding of the scales of heterogeneity will be an invaluable component of a baseline.

54 ENVIRONMENTAL SCIENCES↗

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗

Hydrological and Thermal Dynamics of a Supra‐Permafrost Subterranean Estuary

Subterranean estuaries (STEs), where groundwater interacts with seawater, influence surface and subsurface coastal ecology and biogeochemistry. In Arctic-STEs overlying permafrost, groundwater flow and heat transport determine the fate of organic matter. Yet, direct observations of groundwater flow and heat and solute transport processes in Arctic STEs remain limited. This study characterized groundwater flow paths and fluxes and heat transport within an Arctic-STE along Alaska's Beaufort Sea coast during thawing, summer, and freeze-up. Intertidal seabed temperature-depth profiles collected along a 10-m transect captured the active groundwater flow period, from thaw and flow onset in mid-June to freeze-up in late-September. During this period, aquifer geometry evolved non-uniformly due to spatially varying thaw rates across the STE (mean (m) thaw depths–beach: 0.25 to 0.55–0.6 m on 20 June, 25 July, 1 October; seabed: 0.6–0.9 m from 25 July to 1 October). Groundwater and surface water levels, salinity, and subsurface temperature profiles measured over tidal time scales were interpreted alongside groundwater flow-heat transport numerical simulations. Fresh groundwater discharge was sporadic during thawing (m: 0.32 m 3 /day/m), abundant in summer (m: 0.45 m 3 /day/m), and was largely absent during freeze-up. During freeze-up, groundwater flow was driven exclusively by seawater recirculation via tidal pumping (from thawing to summer to freeze-up: 0.00025–0.15–0.5 m3/day/m) and convection. Heat advection dominated near aquatic interfaces (shaping intertidal ice), and conduction controlled vertical temperature gradients in low-flow and unsaturated sediments. These findings will help predict how prolonged summers will alter Arctic-STE cryo-hydrology and biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗

Formation of wide-blocky calcite veins by extreme growth competition

Liassic limestones on the coast of Somerset in the UK contain dense arrays of calcite microveins with a common, but poorly understood microstructure, characterized by laterally wide crystals that form bridges across the vein. We investigated the mechanisms of formation and evolution of these ‘wide-blocky’ vein microstructures using a combination of high-resolution analytical methods, including virtual petrography, optical cathodoluminescence and scanning electron microscopy techniques (e.g. energy-dispersive X-ray spectrometry, back-scattered electron imaging, cathodoluminescence and electron back-scattered diffraction), laboratory experiments and multiphase field modelling. Our results indicate that the studied veins formed in open, fluid-filled fractures, each in a single opening and sealing episode. As shown by the optical and electron back-scattered diffraction images, the vein crystals grew epitaxially on grains of the wall rock and we hypothesize that their growth rates differed depending on whether the crystals were on a wall rock grain substrate that fractured intergranularly (slow growth rates) or transgranularly (rapid growth rates). Our multiphase field models support this hypothesis, showing that wide, blocky crystals only form where there are significant differences in the growth rate and are dependent on the type of seed grain. These results provide strong evidence for extreme growth competition, a process that we propose controls vein-filling in many micritic carbonate reservoirs, as well as demonstrate that the characteristics of the fracture wall can affect the filling processes in syntaxial veins. Supplementary material: The description and images of the studied thin sections are available at https://doi.org/10.6084/m9.figshare.c.5172371 . High-resolution optical microscopy mosaics (under plane-polarized- and crossed polarized light) of the thin section collection in PetroScan file format are available on request from the authors.

Geology↗

Transmission electron microscopy with in-situ ion irradiation: Facilities and community

Whilst there is a clear scientific and technological need for the technical capabilities of transmission electron microscopes with in-situ ion irradiation, it also requires a collaborative community of international researchers to support such facilities in successfully meeting this demand. Instruments of this type serve to provide fundamental understanding of the mechanisms which drive changes in materials important to nuclear fission and fusion energy, the semiconductor industry, quantum information systems, space travel, astronomy, geology and many more applications. As these areas continue to evolve and the instrumentation possibilities expand, the capacity of in-situ ion irradiation facilities must also develop hand-in-hand with the user community to deliver an ever-greater diversity of high-fidelity extreme-environment experimentation. Future directions for the field, such as miniaturization from MEMS/microfluidic devices and advanced controls with ML-based analysis, continuously emerge to advance both the hardware and software which support the coupling of TEMs with ion beams. This review sets out to provide up-to-date insights into the community and advancement of current, and development of future, facilities which have the potential to further unlock access to the nanoscale exploration of coupled extreme environments crucial to many of the important science and engineering challenges we face today.

In-situ irradiation↗

Characterizing the Atmospheric Mn Cycle and Its Impact on Terrestrial Biogeochemistry

Abstract The role of manganese (Mn) in ecosystem carbon (C) biogeochemical cycling is gaining increasing attention. While soil Mn is mainly derived from bedrock, atmospheric deposition could be a major source of Mn to surface soils, with implications for soil C cycling. However, quantification of the atmospheric Mn cycle, which comprises emissions from natural (desert dust, sea salts, volcanoes, primary biogenic particles, and wildfires) and anthropogenic sources (e.g., industrialization and land‐use change due to agriculture), transport, and deposition, remains uncertain. Here, we use compiled emission data sets for each identified source to model and quantify the atmospheric Mn cycle by combining an atmospheric model and in situ atmospheric concentration measurements. We estimated global emissions of atmospheric Mn in aerosols (<10 μm in aerodynamic diameter) to be 1,400 Gg Mn year −1 . Approximately 31% of the emissions come from anthropogenic sources. Deposition of the anthropogenic Mn shortened Mn “pseudo” turnover times in 1‐m‐thick surface soils (ranging from 1,000 to over 10,000,000 years) by 1–2 orders of magnitude in industrialized regions. Such anthropogenic Mn inputs boosted the Mn‐to‐N ratio of the atmospheric deposition in non‐desert dominated regions (between 5 × 10 −5 and 0.02) across industrialized areas, but that was still lower than soil Mn‐to‐N ratio by 1–3 orders of magnitude. Correlation analysis revealed a negative relationship between Mn deposition and topsoil C density across temperate and (sub)tropical forests, consisting with atmospheric Mn deposition enhancing carbon respiration as seen in in situ biogeochemical studies.

Environmental Sciences & Ecology↗

Controls on Stable Methane Isotope Values in Northern Peatlands and Potential Shifts in Values Under Permafrost Thaw Scenarios

Abstract Northern peatlands are a globally significant source of methane (CH 4 ), and emissions are projected to increase due to warming and permafrost loss. Understanding the microbial mechanisms behind patterns in CH 4 production in peatlands will be key to predicting annual emissions changes, with stable carbon isotopes (δ 13 C‐CH 4 ) being a powerful tool for characterizing these drivers. Given that δ 13 C‐CH 4 is used in top‐down atmospheric inversion models to partition sources, our ability to model CH 4 production pathways and associated δ 13 C‐CH 4 values is critical. We sought to characterize the role of environmental conditions, including hydrologic and vegetation patterns associated with permafrost thaw, on δ 13 C‐CH 4 values from high‐latitude peatlands. We measured porewater and emitted CH 4 stable isotopes, pH, and vegetation composition from five boreal‐Arctic peatlands. Porewater δ 13 C‐CH 4 was strongly associated with peatland type, with δ 13 C enriched values obtained from more minerotrophic fens (−61.2 ± 9.1‰) compared to permafrost‐free bogs (−74.1 ± 9.4‰) and raised permafrost bogs (−81.6 ± 11.5‰). Variation in porewater δ 13 C‐CH 4 was best explained by sedge cover, CH 4 concentration, and the interactive effect of peatland type and pH (r 2 = 0.50,p < 0.001). Emitted δ 13 C‐CH 4 varied greatly but was positively correlated with porewater δ 13 C‐CH 4 . We calculated a mixed atmospheric δ 13 C‐CH 4 value for northern peatlands of −65.3 ± 7‰ and show that this value is more sensitive to landscape drying than wetting under permafrost thaw scenarios. Our results suggest northern peatland δ 13 C‐CH 4 values are likely to shift in the future which has important implications for source partitioning in atmospheric inversion models.

Environmental Sciences & Ecology↗

Synthesis and Single Crystals of Refractory Oxides of Lanthanides and Thorium

At the completion of this program, we can report that we developed a considerable degree of technical improvements in our ability to perform hydrothermal reactions at high temperatures and pressures. We can now routinely perform reactions at 700-750°C and 200 MPa. Currently we are in the process of exploiting this new technology synthesizing a range of exotic new materials investigating relatively poorly understood materials. Our initial efforts focused on the chemistry of rare earth oxides with tetravalent and pentavalent oxides. We recently published a study of the lanthanides with Nb 5+ and Ta 5+ ions, where we grew oxides such as RENdO 4 and RETaO 4 as high quality single crystals. These compounds were targeted as potential hosts for luminescent and scintillation materials, particularly given that they are among the densest oxide hosts and hence have good potential as absorbers for high energy radiation like X-rays and gamma rays. We also isolated a range of unusual new rare earth tantalates with very complex structures. indicating that the chemistry is very sensitive to conditions. We performed some fairly comprehensive examinations of the solid-state chemistry of rare earth ions with various tetravalent metal ions especially Si 4+ , Ge 4+ , Sn 4+ and Ti 4+ . Given the potential role of rare earth silicates in immobilizing radioactive waste elements in long-term storage, and the similarity of our hydrothermal fluids with known geological conditions, this chemistry continues to be relevant. We prepared an extensive series of new lanthanide germanates (e.g. RE 13 Ge 6 O 31 (OH), BaRE 10 (GeO 4 ) 4 O 8 ). and found that there is there is almost no overlap between the chemistry of the rare earth silicates. Stannic oxide (SnO 2 ) is much more refractory and requires higher temperatures and of mineralizer concentrations. One significant result is the growth of RE 2 Sn 2 O 7 pyrochlore single crystals. These are of interest because the rare earth stannate pyrochlores are known to display a wide range of magnetic frustration such as spin ice behavior. We grew high quality single crystals of rare earth germanate and stannate pyrochlores and this led to a collaboration with Professor Kate Ross at Colorado State. Preliminary measurements, indicate that the crystals contain no detectable defects or site disorder. Initial neutron diffraction on single crystals was performed at Oak Ridge, and more detailed experiments involving the Ross group are underway at both NIST and ORNL. This particular chemistry has turned out to be the most potentially significant work on this project and the collaborative effort with Prof. Ross is the topic of a DoE renewal project on quantum materials. Our initial foray into the hydrothermal chemistry of rare earth titanates has also been very promising and a range of cubic and polar ferroic phases of the light rare earths RE 2 Ti 2 O 7 (RE = La - Pr) in the P2 1 phase. We also discovered an interesting new phase Ce 2 Ti 4 O 11 that can have implications in heavy metal immobilization and storage. along with a series of new rare earth titanates (La 5 Ti 4 O 15 (OH) Sm 3 TiO 5 (OH) 3 and Lu 5 Ti 2 O 11 (OH) with exceptionally complex structures. One interesting sidelight has been high temperature hydrothermal chemistry terbium, including the growth of large crystals of TbO(OH). This is not a new compound but it is the first time it has been grown as large single crystals. The Tb atom density is almost as high as that in Tb 2 O 3 and has a very high Verdet constant (ca. 70), making it a very attractive candidate as a Faraday rotator. Unfortunately it is not in a cubic structure but he material is hard, stable, pure and inexpensive, so should still be an attractive Faraday oscillator. We recently received a patent on this material. We also synthesized K 2 Tb(Ge 2 O 7 ) containing stable octahedral Tb 4+ ions, which appears to be the first example of a well-characterized Tb 4+ complex. Given that Tb 4+ has been proposed as a benign surrogate for more treacherous tetravalent ions such as Cf 4+ and Bk 4+ , we think that Tb 4+ silicates can be a particularly useful study for actinide immobilization and related work. We also began reaction studies with rare earths and both ReO 2 and RuO 2 . These resulted in large single crystals of species like RE 5 Ru 2 O 12 , RE 4 Re 2 O 11 , REReO 4 and RE 2 ReO 5 . Several of these samples have already been sent to ORNL for magnetic and neutron diffraction studies.

36 MATERIALS SCIENCE↗

Robust In-Situ Strain Measurements to Monitor CO 2 Storage

The goal of this project was to develop and demonstrate robust instrumentation to monitor the in-situ strain tensor in order to improve the reliability and security of CO 2 storage in geologic formations. We met the original goals of the project and the major overarching accomplishment is the advancement of strain tensor monitoring from an intriguing concept to a commercially available technology with a solid foundation of novel instruments supported by theoretical analyses and validation experiments. The main accomplishments of the project are summarized below. We designed, built and evaluated nine new optical fiber strainmeters and tiltmeters using Michelson interferometers to measure deformation with ultra-high resolution at both shallow and deep point locations in the subsurface. These are the robust strainmeters that motivated the title of the project. We designed, built and evaluated a novel method of measuring distributed strain in optical fibers with nanostrain resolution, and cm-scale location, and sampling into the seismic band. The new method is called Coherence-length-gated Microwave Photonics Interfereometry (CMPI). CMPI technology has advantages over existing commercial DAS and DSS methods. We developed and demonstrated capabilities to deploy instruments in the field and used them to measure strain caused by ambient signals like barometric pressure and tides, as well as induced signals like surface loading and pore pressure changes from pumping tests. We deployed a working strainmeter at 1,700 ft depth, slightly above an active reservoir. This is to our knowledge the greatest depth a strainmeter has been deployed and the techniques we used can readily be extended to greater depths. Optical fiber borehole tensor strainmeter techology was advanced from a TRL 4 at the start, to a TRL of 7 at the conclusion of the project. The project included advances in simulations and theoretical analyses. We developed and demonstrated a computational workflow that uses machine learning to reduce the computational requirements and make it practical to use Bayesian inversion to solve large numerical poroelastic analyses needed to interpret strain tensor field data. We evaluated the strain tensor fields and time series that would be caused by leaks of CO 2 or other fluids from reservoirs. These simulations demonstrated that signals from leaks could be measured with instruments developed for the project, opening a potentially new method for ensuring storage security. We showed that strains in caprock can be used to estimate pressure in a reservoir. This avoids the need to drill monitoring wells into the reservoir, and it expands the capabilities of monitoring in the caprock. The project includes a derivation and application of a novel analytical solution to the strains in the vicinity of a pressurized poroelastic inclusion. This solution explains field data measured during injeciton tests at the North Avant Field, and it will simplify future interpretation of strain tensor data. The project included a broad range of experiments, and of the most significant is the characterization of the strain tensor at an array three strainmeters during six injection tests at the North Avant Field, Oklahoma. This demonstrated repeatability of the strain signal measured by the new instruments developed for the project, and it showed similarities between the strain signal at shallow depths and pressure in the underlying reservoir. We also demonstrated that useful strain data can be measured at reservoir depths. This confirms that strain tensor data can be measured throughout the caprock over a reservoir. The project demonstrated the feasibility of using the strain tensor and distributed strain measured in caprock during a variety of different well tests where the pumping rate was constant, sinusoidal and positive, or a periodic square wave with zero net rate. This further strengthens the validity of using strain data to characterize reservoirs and aquifers. We also demonstrated that strain caused be fluctuations of air pressure and water pressure in the vadose zone can be measured and interpreted, suggesting that high resolution distributed strain measurements hold promise for monitoring the vadose zone. The project partially supported nine graduate students in the Environmental Engineering, Hydrogeology, Electrical Engineering programs at Clemson University. The research was described in nine journal papers, 23 talks and conference abstracts. Additional journal papers are in preparation. A new company called Tensora was started to provide strainmeter technology for commercial applications.

01 COAL, LIGNITE, AND PEAT↗

Dynamics of streamflow permanence in a headwater network: Insights from catchment-scale model simulations

The hillslope and channel dynamics that govern streamflow permanence in headwater systems have important implications for ecosystem functioning and downstream water quality. Recent advancements in process-based, semi-distributed hydrologic models that build upon empirical studies of streamflow permanence in well-monitored headwater catchments show promise for characterizing the dynamics of streamflow permanence in headwater systems. However, few process-based models consider the continuum of hillslope-stream network connectivity as a control on streamflow permanence in headwater systems. The objective of this study was to expand a process-based, catchment-scale hydrologic model to better understand the spatiotemporal dynamics of headwater streamflow permanence and to identify controls of streamflow expansion and contraction in a headwater network. Further, we aimed to develop an approach that enhanced the fidelity of model simulations, yet required little additional data, with the intent that the model might be later transferred to catchments with limited long-term and spatially explicit measurements. This approach facilitated network-scale estimates of the controls of streamflow expansion and contraction, albeit with higher degrees of uncertainty in individual reaches due to data constraints. Our model simulated that streamflow permanence was highly dynamic in first-order reaches with steep slopes and variable contributing areas. The simulated stream network length ranged from nearly 98±2% of the geomorphic channel extent during wet periods to nearly 50±10% during dry periods. The model identified a discharge threshold of approximately 1 mm d –1 , above which the rate of streamflow expansion decreases by nearly an order of magnitude, indicating a lack of sensitivity of streamflow expansion to hydrologic forcing during high-flow periods. Overall, we demonstrate that process-based, catchment-scale models offer important insights on the controls of streamflow permanence, despite uncertainties and limitations of the model. We encourage researchers to increase data collection efforts and develop benchmarks to better evaluate such models.

54 ENVIRONMENTAL SCIENCES↗

Subtask 1.5 – CO2 Injection Monitoring with an Optimized Scalable, Automated, Semipermanent Seismic Array

The scalable, automated, semipermanent seismic array (SASSA) method is a flexible and relatively cost-effective surface geophysical method for regular time-lapse monitoring of the movement of injected carbon dioxide (CO2) in a reservoir for CO2 enhanced oil recovery (EOR) or geologic CO2 storage operations. It has the advantages of a low-environmental-footprint while monitoring regions of a reservoir from the surface without the need for a regular grid distribution of receivers. Automated data collection is possible. As only time-lapse amplitude changes at the reservoir level due to CO2 movement within the reservoir are monitored, the turnaround time to deliver results from the SASSA method can be short, without the need for long, time-consuming data-processing workflows. As data is collected and processed, incremental information can be provided to the field operator. The Energy & Environmental Research Center (EERC) conducted a SASSA field test from September 2018 to November 2020 in a portion of the Bell Creek Field in Montana, which implemented new CO2 EOR field activities during the study period. Lessons learned from a proof-of-concept study were incorporated to improve the data quality of the SASSA method and demonstrate the viability of the technology. The EERC implemented several enhancements to improve data quality, including 1) an iterative survey design, which allowed placing the receivers in strategic locations where the movement of the CO2 in the reservoir could be tracked with minimum interference by the cultural noise in the study area; 2) the use of powerful seismic sources in the form of surface orbital vibrators, and 3) data acquisition during optimal periods. History-matched reservoir simulation was performed to predict gas saturation and pressure response induced by CO2 injection in the study area. The results were compared with the SASSA-measured responses to CO2 injection as a partial validation technique. A match between the two methods was observed for most of the SASSA points predicted to have intersected a CO2 saturation change. The validated results provide confidence that the SASSA method can be used independently as a CO2 saturation monitoring technique. As data are collected and processed, incremental information can be provided to the field operator. The critical components of the SASSA workflow for a successful application of the method are the following: Iterative survey design with information about CO2 injection activities from the oilfield operator. A detailed CO2 injection plan is the key driver to select the strategic monitoring location of the SASSA sensors. After this information is incorporated in the initial distribution of sources and receivers in the study area, high-resolution satellite images are used to identify ground locations not affected by cultural noise sources, such as power lines, pipelines/flow lines, or roadways. In the next iteration of the survey design, a scouting trip to the study area is needed to understand more details of the noise sources identified in the previous step and the intensity of the field activities that can also generate noise during the monitoring. Integrating the information from the scouting trip into the survey design to select the optimum source and receiver locations is the final step. Noise attenuation. The variety of noise types during seismic monitoring of an oil field is enormous. Tailored noise characterization and processing at a node-by-node level can enhance the performance and sensitivity of the SASSA technique. Future advancements that could improve the efficiency and application of the SASSA technology include: Gaining a better understanding of the noise field produced by the seismic source to aid the choice of receiver location. Surface noise from the source can overwhelm the small signal changes due to CO2 that the SASSA method measures. Improved data-processing workflow to automatically analyze and adapt to dynamic noise conditions associated with industrial settings. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE- FE0024233.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

High-resolution analysis of clay minerals and amorphous materials in martian analog environments

Mudrocks and mud-sized sediments (i.e., silt to clay) dominate the surface of Earth and Mars. These fine-grained sediments preserve a rich history of sedimentary processes from source to sink and shed light on ancient climates. However, both the physical and chemical nature of these materials make them difficult to fully characterize with traditional laboratory techniques. Here, we explore a cross-disciplinary and high-resolution approach using synchrotron radiation for X-ray diffraction, pair distribution function analysis, and submicrometer-scale X-ray fluorescence, combined with transmission electron microscopy, to better understand the nanostructure and composition of mud-sized sediments from a glacio-fluvial watershed in southwest Iceland. Our results demonstrate that sediments in the cold and wet climate of Iceland are more altered than previously thought, as evidenced by the identification of kaolinite and mixed-layer kaolinite-smectite. Additionally, sediments are enriched in amorphous materials and nanocrystalline phases, as determined from grain morphologies and compositions consistent with allophane, hisingerite, ferrihydrite, and halloysite. These alteration products are present as intimate mixtures that vary across depositional sites, demonstrating the dynamic nature of the secondary assemblage from source to sink. This work has implications for Mars, where, for example, basalt-sourced sedimentary rocks from Gale crater are abundant in clay minerals and amorphous materials. Finally, this work underpins the importance of using high-resolution techniques, a coordinated methodology, and developing innovative approaches for future planetary sample return missions (e.g., Mars sample return).

36 MATERIALS SCIENCE↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

Topical Report – Findings on Subtask 2.7 – Wet ESP and Aerosol Testing at Coal Creek Station

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion carbon capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery. The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the unique components present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants by contributing to aerosol formation. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and flue gas particulate provides nucleation sites where amine vapors can condense as aerosols. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, particulate components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to fully characterize aerosol behavior with various control technologies installed to better optimize aerosol mitigation technology for CO 2 capture. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on mitigating formation of problematic aerosols at Great River Energy’s Coal Creek Station, upstream of the PCCC system. Determine the effectiveness of the Mitsubishi Heavy Industries (MHI) proprietary amine emission reduction unit (AERU) as a postcapture solvent recovery system for reducing aerosol emissions and extending solvent life downstream of the PCCC system. Determine the impact of aerosols on the efficiency and degradation products of both commercial and advanced solvents within the PCCC system. Work was conducted at Coal Creek Station Unit 1 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then optionally routed through a WESP before passing to the CO 2 absorber columns. The MHI solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . Flue gas exiting the absorber column was routed to MHI’s proprietary AERU to recover entrained solvent. The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Particulate was measured, collected, and analyzed from multiple locations throughout the pilot-scale system. Unlike the performance observed in prior work, the inlet FGD and DCC did not remove significant particulate matter from the flue gas. This appears to be due to a difference in the nature of the particulate. The DCC seemed to increase particulate size and count, most likely owing to water condensing onto the surfaces of fly ash particles. When the WESP was operated, it achieved >95% particulate capture. Very little particulate matter or indications of solvent were detected at the AERU outlet. When operating with advanced KS 21 solvent, the particulate material at the AERU outlet was even further decreased. Solvent analysis showed that some species derived from flue gas and ash were slowly concentrating in the solvent over the duration of the test. The levels observed were reported to be within expected ranges and were not of concern to MHI. A high-level techno-economic assessment of installing CO 2 capture at Coal Creek Station suggested that, when using a standard monoethanolamine (MEA)-based solution with simple heat integration, the energy penalty to net generation would be 34%. The bulk of this was due to steam losses for regenerating solvent, followed by parasitic electrical demand for CO 2 compression and then by increased parasitic load for pushing flue gas through the absorber column. These demands could be decreased with a more advanced solvent that exhibits lower heat of regeneration and lower pressure drop than does a simple MEA solution. Further energy could be saved with more thorough heat integration to recover useful energy from the steam used for solvent regeneration. Installing a WESP was predicted to increase the cost of electricity by nearly $5/MWh. This would become cost-effective if solvent losses were roughly 10 times the baseline estimate when not using a WESP but could be returned to baseline by installing the WESP. Piping CO 2 for storage in more favorable geology could help with carbon capture and storage economics. Although storing off-site would necessitate construction of a CO 2 transport pipeline, the cost of this pipeline might be more than offset by reducing the number of wells required, the depths of the wells required, and the electrical demand for the CO 2 compressor. Additional factors that favor off-site storage costs include smaller expected CO 2 plume sizes, which translates to less monitoring and fewer landowner agreements. More detailed assessment of the specific geology in the region under and around Coal Creek Station would be needed to accurately assess the costs and benefits of different storage site options. This subtask was cofunded through the Energy & Environmental Research Center–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission.

Strege, Joshua↗