Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VOLCANICS”

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.

306 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Preliminary Radiological Risk Assessment Model for Disposition of Remote-Handled Transuranic Wastes at Los Alamos National Laboratory Area G - 20116

The U.S. Department of Energy (DOE) operates a low-level radioactive waste (LLW) disposal site at Material Disposal Area G, in Los Alamos, New Mexico, USA. Area G has been the primary LLW disposal site for Los Alamos National Laboratory (LANL) since the 1960's. In addition to LLW, Area G is host to a variety of other wastes, the disposition of which must be determined before closure of the site. A probabilistic Radiological Risk Assessment (RRA) for Area G is used in order to support decision making regarding some wastes that are not addressed in the extant Area G Performance Assessment (PA) and Composite Analysis (CA). Between 1979 and 1987, 33 special shafts were augered into the Bandelier Tuff at Area G. This volcanic tuff is present across Pajarito Plateau on the eastern slopes of the Jemez Mountains, and varies widely in its consistency, from weakly indurated non-welded layers to welded layers that uphold the mesa cliffs of the Plateau. These mesas are home to LANL, Area G, and the townsites of Los Alamos and White Rock, with residences about 1400 m from Area G. The 33 Shafts were lined with steel casing, and contain remote-handled (RH) transuranic wastes (TRU) resulting from experiments and analysis performed in special glove boxes at the Chemistry and Metallurgy Research (CMR) facility at LANL. Some of these wastes originated as used nuclear fuel. The purpose of the Area G RRA is to evaluate the potential future risk to humans and the environment from the RH TRU in the 33 Shafts in the context of the risk associated with the surrounding wastes at Area G. The analysis is responsive to expectations outlined in DOE Order 458.1, Radiation Protection of the Public and the Environment, and is informed by the Manual and Guidance accompanying DOE O 435.1, Radioactive Waste Management. Because the waste meets the definition of TRU, the regulatory context necessarily takes into consideration the regulation governing the disposal of TRU from the U.S. Environmental Protection Agency (EPA): 40 CFR 191, Environmental Radiation Protection Standards for Management and Disposal of Spent Nuclear Fuel, High-Level and Transuranic Radioactive Wastes. Given the broader regulatory context for the RRA, the analysis is subject to different assumptions from those made in the existing DOE O 435.1 PA and CA, such as allowing for future occupation of the site. The analysis begins with a comprehensive evaluation of features, events, processes, and exposure scenarios (FEPS) for Area G and the wastes it contains. These FEPSs are screened to eliminate from further consideration those of extremely low probability and/or consequence, and a conceptual site model (CSM) is subsequently developed. The scope and structure of the Area G RRA Model is informed by this CSM, and the Area G RRA Model is developed using the GoldSim systems analysis modeling platform. This paper presents the initial version of a defensible, transparent, and reasonably realistic model, which is based on the state of knowledge of the wastes, the site, and the FEPSs that govern contaminant transport from wastes into the environment and subsequent exposures to humans and other biota. Probabilistic model input distributions represent uncertainties inherent in the real and modeled systems. The results of the Area G RRA Model inform decisions regarding the disposition of the RH TRU in the 33 Shafts. (authors)

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

DOE Hanford Tank Side Cesium Removal System (TSCR) Update: Fabrication and Factory Acceptance Testing - 20300

The Tank Side Cesium Removal (TSCR) System is a modular, quickly-deployable, and passively safe system for efficiently removing solids and cesium from alkaline liquid waste raffinate created by fuel processing on the Hanford site. This paper provides an overview of the project with particular emphasis on completion of fabrication and factory acceptance testing (FAT). All TSCR System components are inside three enclosures known as the Process Enclosure, Ancillary Enclosure, and Control Enclosure. The Process Enclosure is seismic design category II (SDC-II) and built to withstand the rigors of transportation and lifting, high winds, volcanic ash fall, and heavy precipitation (snow) over a temperature range of -25 to 115 deg. F. The Process Enclosure and Control Enclosure were fabricated in AVANTech's Columbia, SC facility, while the Ancillary Enclosure was fabricated in AVANTech's Richland, WA facility. The design and preliminary documented safety analysis were completed in December 2019, with all work being performed in accordance with ASME NQA-1 quality assurance requirements. The Process Enclosure, Ion Exchange Columns (IXCs), and several other TSCR components have a nuclear safety function and are designated as safety significant structures, systems, and components (SS-SSCs). The NQA-1 Commercial Grade Dedication (CGD) process was used to identify critical characteristics in order to accept materials and items used for fabrication of several safety significant TSCR components. This paper also addresses unique aspects and challenges of implementing CGD in DOE facilities and systems. After fabrication and preliminary testing, the Process and Control Enclosures were shipped to AVANTech's Richland facility where the complete TSCR System was setup and connected in a represented layout and configuration that will be used on the Hanford site. Test tanks and pumps were put in place and connected to the TSCR Process Enclosure in a manner that mimicked interfacing Double Shell Tanks in the 241-AP Tank Farm. The TSCR System was then tested/operated under conditions equivalent to those planned for onsite operation. Testing was separated into five separate phases as defined below: - Phase 1: Confirmation of TSCR System readiness for the FAT; - Phase 2: Verification of TSCR System operation under standard conditions; - Phase 3: Demonstration of TSCR System operation under non-standard conditions; - Phase 4: Demonstration of ion exchange media removal from an IXC; and - Phase 5: Completion of TSCR System lift plans and site mobilization. The Hanford Tank Operations Contractor, Washington River Protection Solutions (WRPS), participated in all phases of testing. Testing at AVANTech's Richland Facility allowed WRPS personnel to familiarize themselves with TSCR operations in a radiologically clean facility prior to onsite mobilization. Engineering and operations personnel gained valuable experience that allowed them to better refine onsite operational protocols and procedures - including handling and movement of IXCs with a forklift. This paper describes the technologies used by the TSCR System as well as how these technologies are deployed and operated to advance the Hanford mission of safely, efficiently, and effectively treating tank waste. It also provides detailed information on fabrication and FAT, including operational data from simulated tank waste runs. A project schedule for WRPS' planned construction acceptance testing, operational acceptance testing, and actual tank waste treatment is provided. This presentation should be of great interest to parties responsible for the design, testing, and mobilization of modular NQA-1 technologies for the treatment of raffinate and similar liquid wastes containing high concentrations of radioactive cesium. (authors)

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

Exploring Geothermal Potential of Great Basin Sub-Regions: Preprint

The INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems (INGENIOUS) project aims to discover new, economically viable hidden geothermal systems in the Great Basin region by building on previous work in play fairway analysis and machine learning. A key objective of this project is to develop an exploration workflow to reduce geothermal exploration risks for hidden geothermal systems. A single preliminary play fairway workflow was developed from the assessment of the regional INGENIOUS geological, geophysical, and geochemical datasets. This workflow provided new preliminary predictive geothermal fairway maps for the INGENIOUS study area, which encompasses most of Nevada, western Utah, southern Idaho, southeastern Oregon, and easternmost California. However, a recent study (incorporating machine learning techniques) of a portion of Nevada identified four geologic domains and determined that the relative importance of individual datasets or features as indicators of geothermal potential may differ across these domains. The INGENIOUS study area includes a much larger and more geologically diverse region; therefore, additional geologic domains or sub-regions are expected. To assess the sub-regions in the INGENIOUS study area, principal component analysis and k-means clustering were applied. Preliminary results indicate that the INGENIOUS regional data cluster into groups that relate to different geologic domains in the Great Basin region. These include domains such as the Walker Lane, extensional western Great Basin region, broad lower strain region in the eastern Great Basin of western Utah and eastern Nevada, Quaternary volcanic fields, and the area adjacent to the Snake River Plain. These clusters are assessed to determine the key geologic drivers of the identified clusters. Understanding this variability can provide key insights for the exploration and characterization of hidden geothermal systems in the Great Basin region and could indicate the need to develop multiple geothermal conceptual models and play fairway workflows for the INGENIOUS study area.

GEOTHERMAL ENERGY↗
DOE OSTI2025

Exploring Geothermal Potential of Great Basin Sub-Regions

The INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems (INGENIOUS) project aims to discover new, economically viable hidden geothermal systems in the Great Basin region by building on previous work in play fairway analysis and machine learning. A key objective of this project is to develop an exploration workflow to reduce geothermal exploration risks for hidden geothermal systems. A single preliminary play fairway workflow was developed from the assessment of the regional INGENIOUS geological, geophysical, and geochemical datasets. This workflow provided new preliminary predictive geothermal fairway maps for the INGENIOUS study area, which encompasses most of Nevada, western Utah, southern Idaho, southeastern Oregon, and easternmost California. However, a recent study (incorporating machine learning techniques) of a portion of Nevada identified four geologic domains and determined that the relative importance of individual datasets or features as indicators of geothermal potential may differ across these domains. The INGENIOUS study area includes a much larger and more geologically diverse region; therefore, additional geologic domains or sub-regions are expected. To assess the sub-regions in the INGENIOUS study area, principal component analysis and k-means clustering were applied. Preliminary results indicate that the INGENIOUS regional data cluster into groups that relate to different geologic domains in the Great Basin region. These include domains such as the Walker Lane, extensional western Great Basin region, broad lower strain region in the eastern Great Basin of western Utah and eastern Nevada, Quaternary volcanic fields, and the area adjacent to the Snake River Plain. These clusters are assessed to determine the key geologic drivers of the identified clusters. Understanding this variability can provide key insights for the exploration and characterization of hidden geothermal systems in the Great Basin region and could indicate the need to develop multiple geothermal conceptual models and play fairway workflows for the INGENIOUS study area.

exploration↗
DOE OSTI2024

Shallow Geothermal Resources for Cooling Applications at the University of Hawai‘i

Scientists at Berkeley Lab have teamed up with The University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s (DOE’s) Energy Transitions Initiative Partnership Project (ETIPP) to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology for building cooling, energy efficiency, and emissions reduction applications in Hawai’i. The team is assessing the data necessary to model the feasibility of deploying this technology, the actual models that will be used, and what hurdles need to be overcome to install a demonstration case. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project, and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, making it crucial to select sites with large enough values of these properties. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Both closed-loop and open-loop systems are being investigated. Another option being considered is using cool seawater as the source of chill. Currently work is progressing on two fronts. A hydrogeologic model for a closed-loop system is being developed for the Stan Sheriff Center at the UH Manoa campus, where a subsurface karst system immediately downgradient of the borefield may provide efficient removal of heated groundwater. The team will also develop a technoeconomic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air conditioning system. At the state scale, geographic information system (GIS) layers of various attributes relevant for GHE are being combined to develop an overall favorability map for employing GHE in Hawai‘i.

Doughty, C↗
DOE OSTI2025

The CO2 mineralization potential of New Mexico basalts and associated critical element behavior and mobility.

This presentation reports on laboratory experiments assessing the carbon dioxide (CO₂) mineralization potential of New Mexico basalts, including trachybasalt and basalt lithologies. Using simulated groundwater and controlled conditions, the study evaluates the kinetics of CO₂ sequestration, mineral saturation indices, and secondary mineral formation. It also investigates the mobility of critical and trace elements such as Ni, Zn, Cu, and Li during basalt dissolution. Results indicate that New Mexico basalts can effectively sequester CO₂ through carbonate formation, with grain size and mineral composition influencing reaction rates and elemental behavior. These findings support the viability of geologic carbon storage in mafic volcanic terrains and inform strategies for critical element recovery.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Isolated vesicles in submarine pumice: insights from the 2019 Volcano F eruption, Kingdom of Tonga

No reliable method can distinguish between subaerial and submarine pumice, but recent work shows that highly vesicular (total porosity, φ > 65%) submarine pumice from the 2012 Havre and 2021 Fukutoku Oka-no-Ba eruptions contains abundant isolated porosity (connectivity, C < 0.6). This differs from textural measurements of subaerial pumice with similarly high vesicularities, where high connectivities (C > 0.8) are often measured. To investigate the implications of this under-studied texture for bubble nucleation, growth, and coalescence dynamics at submarine eruptions, we analyze the textural characteristics of rafted pumice clasts from the 2019 submarine eruption of Volcano F, located along the Tonga-Kermadec Arc. We examined clasts collected while floating near the vent (50–200 km) and distally on the shores of Fiji (900 km). Using helium pycnometry and X-ray tomography (XRT), we quantified porosity and connectivity in 45 pumice lapilli, classifying them by vesicle macrotexture and edge morphology. Microvesicular clasts, the majority of which have breadcrust and cauliflower textures, exhibit the lowest connectivities (C = 0.61 ± 0.2), while macrovesicular clasts have high connectivities with C ~ 1. We propose that isolated porosity abundance varies according to clast location in the eruption column post-fragmentation. Microvesicular clasts erupted on the edges of the column and experienced high levels of seawater ingestion and subsequent quenching of the entire clast. Macrovesicular clasts, in contrast, erupted in the center of the eruption column and were thus thermally insulated, allowing magmatic bubbles in the molten pumice to nucleate, grow, and coalesce for longer. Our results and interpretation imply that low connectivity, microvesicular clasts preserve pre and syn-eruptive bubble processes at the 2019 Volcano F eruption.

Permeability↗
DOE OSTI2021

DEEPEN Data Catalog for Magmatic Geothermal Systems in the United States

This data catalog contains information related to the Training Site Analysis for the Geothermica project "DE-risking Exploration of geothermal Plays in magmatic ENvironments (DEEPEN)." The DEEPEN project aims to reduce exploration risk for geothermal fluids in magmatic systems by developing improved an improved framework for interpretation of exploration data using the Play Fairway Analysis (PFA) methodology. The Training Site Analysis performed for DEEPEN leverages existing datasets to develop a customized PFA approach to exploration for multiple geothermal resource types in magmatic systems (conventional hydrothermal resources, supercritical fluid and superheated steam resources, and superhot EGS resources). This data catalog contains links to publicly available data files related to 8 training sites in the United States. US training sites are: the Cascades/Aleutians PFA project; the Hawaii PFA project, the Oregon Cascades PFA project, the Snake River Plain, Idaho PFA project, the Washington State PFA project, Newberry Volcano, Coso Geothermal Field, and the Geysers Geothermal field. This database contains an overview of these training sites, data sources, and links to publicly available exploration datasets. For the five PFA projects, details on exploration data related to PFA components (heat, fluid, permeability, sometimes seal) are provided, including a summary of data weighting methodologies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

A collective trigger for widespread planetesimal formation revealed by accretion ages

The formation of planetesimals was an integral part of the cascading series of processes that built the terrestrial planets. To illuminate planetesimal formation, here we develop a refined thermal evolution model to calculate the formation ages of meteorite parent planetesimals. This model includes chemical reactions and phase changes during heating, as well as natural variations in the proportions of the constituent phases of these planetesimals. We find that the parent bodies of non-carbonaceous (NC) and carbonaceous (CC) iron meteorites start forming at very similar times ( ∼ 0.95 Myr after calcium-aluminium-rich inclusion [CAI] formation) and occupy overlapping time windows. NC and CC chondrite parent bodies formed later during non-overlapping periods. We combine these ages with proportions of isotopic end-members we recover from mixing models to construct records of motion throughout the protoplanetary disk. These records argue that NC and CC material traversed the barrier in the disk after ∼ 0.95 Myr after CAI formation. The onset of this motion coincided with planetesimal formation, indicating that the phenomenon that drove motion also triggered planetesimal formation. We argue that this feature also served as the semi-permeable barrier in the disk. Although its identity is uncertain, the effects this phenomenon had on the timing of planetesimal formation and motion through the disk can now serve as constraints on models of disk evolution. Models that reproduce these effects would elucidate the nature and implications of this phenomenon, which is key to unlocking a holistic model of terrestrial planet building.

Astronomy and AstroPhysics↗
DOE OSTI2025

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide↗
DOE OSTI2021

Impacts of Climate Change on Groundwater Availability and Spring Flows: Observations from the Highly Productive Medicine Lake Highlands/Fall River Springs Aquifer System

Medicine Lake Highlands/Fall River Springs Aquifer System, located in northeastern California, is home to some of the largest first-order springs in the United States. This work assesses the likely effects of projected climate change on spring flow. Four anticipated climate futures (GFDL A2, GFDL B1, CCSM4 rcp 8.5, CNRM rcp 8.5) for California, which predict a range of conditions (generally warming and transitioning from snow to rain with variable amounts of total precipitation), are postulated to affect groundwater recharge primarily by changing evapotranspiration. The linkages between climate variables and spring flow are evaluated using a water balance model that represents the physics of evapotranspiration and recharge, the Basin Characterization Model. Three of the four climate scenarios (GFDL A2, GFDL B1, CCSM4 rcp 8.5) project that by the year 2100, groundwater recharge (and consequently decreased spring flow) will decrease by 27%, 21%, and 9%, respectively. The fourth scenario (CNRM rcp 8.5) showed an increase in recharge of 32% due to a significant increase in precipitation (27%). Evapotranspiration increases due to a shift in the type of precipitation and a longer growing season. While the likelihood of each scenario is outside the scope of this work, unless total precipitation increases dramatically in the future, increased temperatures and decreasing precipitation will likely result in reduced spring flows, along with warmer water temperatures in downstream habitats.

42 ENGINEERING↗
DOE OSTI2021

NACHOS, a CubeSat-based high-resolution UV-Visible hyperspectral imager for remote sensing of trace gases: system overview and science objectives

Here, the Nano-satellite Atmospheric Chemistry Hyperspectral Observation System (NACHOS) is a high-throughput (f/2.9), high spectral resolution (~1.3 nm optical resolution, 0.6 nm sampling) Offner-design hyperspectral imager operating in the 300-500 nm spectral region. The 1.5U instrument payload (1U optical system, 0.5U electronics module) is hosted by a 1.5U LANL-designed CubeSat bus to comprise a 3U complete satellite. Spectroscopically similar to NASA’s Ozone Monitoring Instrument (OMI), which provides wide-field global mapping of ozone and other gases at coarse spatial resolution, NACHOS fills the complementary niche of targeted measurements at much higher spatial resolution. With 350 across-track spatial pixels and a 15-degree across-track field of view, NACHOS will provide spectral imaging at roughly 0.4 km per pixel from 500 km altitude. NACHOS incorporates highly streamlined gas-retrieval algorithms for rapid onboard processing, alleviating the need to routinely downlink massive hyperspectral data cubes. We will discuss the instrument design, challenges in achieving mechanical robustness to launch vibration in such a compact instrument, the onboard calibration system, and gas-retrieval data downlink strategy. We will also discuss potential science missions, including monitoring of NO 2 as an easily detected proxy for anthropogenic fossil-fuel greenhouse gases, monitoring lowlevel SO 2 degassing at pre-eruptive volcanoes, H 2 CO from wildfires, and characterization of aerosols. The long-term vision is for a many-satellite constellation that could provide both high spatial resolution and frequent revisits for selected targets of interest. As an initial technology demonstration of this vision, the NACHOS project is currently slated to launch two CubeSats in early 2022.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

The Indian Ocean, its supra-subduction history, and implications for ophiolites

We report Ophiolite complexes represent fragments of ocean crust and mantle formed at spreading centers and emplaced on land. The setting of their origin, whether at mid-ocean ridges, back-arc basins, or forearc basins has been debated. Geochemical classification of many ophiolite extrusive rocks reflect an approach interpreting their tectonic environment as the same as rocks with similar compositions formed in various modern oceanic settings. This approach has pointed to the formation of many ophiolitic extrusive rocks in a supra-subduction zone (SSZ) environment. Paradoxically, structural and stratigraphic evidence suggests that many apparent SSZ-produced ophiolite complexes are more consistent with mid-ocean ridge settings. Compositions of lavas in the southeastern Indian Ocean resemble those of modern SSZ environments and SSZ ophiolites, although Indian Ocean lavas clearly formed in a mid-ocean ridge setting. These facts suggest that an interpretation of the tectonic environment of ophiolite formation based solely on their geochemistry may be unwarranted. New seismic images revealing extensive Mesozoic subduction zones beneath the southern Indian Ocean provide one mechanism to explain this apparent paradox. Cenozoic mid-ocean-ridge–derived ocean floor throughout the southern Indian Ocean apparently formed above former sites of subduction. Compositional remnants of previously subducted mantle in the upper mantle were involved in generation of mid-ocean ridge lavas. The concept of historical contingency may help resolve the ambiguity on understanding the environment of origin of ophiolites. Many ophiolites with “SSZ” compositions may have formed in a mid-ocean ridge setting such as the southeastern Indian Ocean.

58 GEOSCIENCES↗
DOE OSTI2021

Data Arrays for Microearthquake (MEQ) Monitoring using Deep Learning for the Newberry EGS Sites

The 'Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties' project looks to apply machine learning (ML) methods to Microearthquake (MEQ) data for imaging geothermal reservoir properties and forecasting seismic events, in order to advance geothermal exploration and safe geothermal energy production. As part of the project, this submission provides data arrays for 149 microearthquakes between the year 2012 and 2013 at the Newberry EGS Site for use with the Deep Learning Algorithm that has been developed. The data provided includes raw waveform data, location data, normalized waveform data, and processed waveform data. Penn State Geothermal Team has shared the following files from the project: - 149 microearthquakes (MEQs) between 2012 and 2013 at Newberry EGS sites, 'Normalized Waveform Inputs.npz' are normalized waveforms. - labels of 149 MEQs: Processed Waveform Inputs.npz - location labels of 149 MEQs: Location Data.npz Note: .npz is the python file format by NumPy that provides storage of array data.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗
DOE OSTI2022

Estimating Explosion Yields Using Moment Tensor Solutions and Seismic Moment

We report seismic moment, a measurable and well-understood quantity of seismic sources, is used to estimate the yield of explosions. Application of such a method in the past, as in the manner of m b -derived yields, has been complicated by the effect of variations in the explosion working point, depth, and secondary source effects (such as spalling and tectonic release) on the observed moment. We start using the full (six-element) moment tensor solution, which can capture the relevant source physics and, at least in theory, better isolate the primary explosion source. The moment-to-yield ratio is then estimated using an explosion source model which, provided with emplacement conditions, can relate the two parameters. We discuss the major sources of uncertainty associated with the method, and calibrate it with chemical and nuclear explosions at the Nevada National Security Site. We then apply the method to published moment tensor solutions for the six declared North Korean nuclear explosions that occurred between 2006 and 2017. The results are mostly consistent with other yield estimates made using a variety of high-frequency methods. This technique is a new approach to estimating explosive yield and simple to implement, as much of the complexity is captured by the source models.

58 GEOSCIENCES↗
DOE OSTI2023

Gas-mediated trace element incorporation into rhyolite-hosted topaz: A synchrotron microbeam XAS study

Magmatic gas exsolving during late-stage cooling of shallow magmas has been considered an important facilitator of low-pressure alteration and metal transport. However, the chemical properties of such gas, particularly its metal transport mechanisms and capacity, remain elusive. Trace elements in minerals produced by gas-mediated surface reaction or precipitation from gas capture details of gas composition and reaction pathways. However, interpretation of mineral trace element contents is dependent on understanding crystallographic controls on gas/mineral partitioning. Here, this work investigates the structural accommodation of As, Mn, Ga, Ge, Fe, and Ti in vapor-deposited topaz of vesicular topaz rhyolite from the Thomas Range, Utah, through single-crystal synchrotron microbeam X-ray techniques on picogram quantities of those trace elements. X-ray absorption near edge structure (XANES) data indicates that these elements are incorporated into topaz as As 5+ , Fe 3+ , Mn 3+ , Ti 4+ , Ga 3+ , and Ge 4+ . Extended X-ray absorption fine structure (EXAFS) analysis for these trace elements, compared to EXAFS of structural Al and Si, reveals that As 5+ and Ge 4+ are incorporated directly into the tetrahedral site of the topaz structure, with the octahedral site accommodating Mn 3+ , Fe 3+ , Ga 3+ , and Ti 4+ . For As 5+ and Fe 3+ , the structural impact of substitution extends to at least second neighbors (other elements were only resolvable to first neighbors). Further interpretation of the EXAFS results suggests that the substitution of Ti 4+ results in increased distortion of the octahedral site, while the other trace elements induce more uniform expansion correlating in magnitude to their ionic radius. Comparison of quantified X-ray fluorescence (XRF) data for two topaz crystals from this rhyolite reveals variable trace element concentrations for As 5+ , Fe 3+ , Ga 3+ , and Ti 4+ , reflective of a source gas undersaturated in these trace elements changing in concentration over the period of topaz deposition. The identical Ge 4+ content of the two topaz crystals suggests that Ge 4+ in the gas was buffered by the growth of another Ge 4+ -bearing phase, such as quartz. The very low Mn 3+ content in the topaz crystals does not reflect the abundance of Mn 3+ in the gas (saturation of Mn is evidenced by coexisting bixbyite). Instead, it suggests a strong Jahn-Teller inhibitory effect to the substitution of Mn 3+ for Al 3+ in the distorted octahedral site of topaz. It is proposed that exsolution of an HF-enriched gas from cooling rhyolitic magma led to local scouring of Al, Si, and trace metals from the magma. Once topaz crystals nucleated, self-catalyzed reactions that recycle HF led to continued growth of topaz.

36 MATERIALS SCIENCE↗
DOE OSTI2023

RARE EARTH AND CRITICAL ELEMENT RELATIONSHIPS TO COAL DEPOSIT LITHOLOGIES IN SOUTHWEST WYOMING AND NORTHWEST COLORADO

Rare earth elements (REE) and critical minerals (CM) are vital components in a wide variety of applications, not only in the renewable energy sector but also in military, telecommunication, and other civilian uses. The association of REE and CM with coal deposits is widely known but not thoroughly understood. While individual coal samples or beds themselves might not host significant concentrations of REE, recent studies indicate temporal relationships among some REE with coal beds and their immediately adjacent non-coal bounding strata, including interbeds. These strata most commonly include clays, carbonaceous shales, and ash beds, but other related strata are possible in different coal-bearing settings. Variations in REE and CM concentrations and their direct or indirect association with coal beds are not always predictable. This study presents recent findings from studying different associations between coal beds and adjacent rocks to identify possible relationships between REE-CM concentrations and the coals and adjacent rocks. For this study we collected and analyzed several dozen samples from the U.S. Geological Survey’s Core Research Center (CRC). Additionally, a small number of coal mine outcrops and their adjacent non-coal rocks were included. Most samples are of Cretaceous age. Previous work has shown a potential link between coals in the Upper Cretaceous Williams Fork Formation and REE occurrences in the Yampa ash bed in northwest Colorado. While data are uncommon and sparse, this possible link may suggest REE or CM enrichment from leaching or reworking of the ash bed. Trace element geochemistry, whole rock data, and microscopic examinations are used here to explore these relationships and guide future studies.

REE, CM, ash, volcanic ash, coal↗