Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Volcanism”

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 55 records · Page 3

Constraining the Hydration of Clay Minerals and Abundances of Amorphous Phases in Gale Crater, Mars

Both water and organic matter are required for the development and persistence of life. Phyllosilicates (clay minerals) have high surface areas that easily sorb water and organic matter. The Curiosity rover has investigated several hundred meters of stratigraphy in Gale crater, including where clays were detected from orbit. Previous results have suggested that subsurface hydration is greatest in units with the most abundant clays, suggesting that these minerals may be hydrated. Organics have also been found throughout Gale crater. Smectites are the most common and abundant phyllosilicates in Gale crater samples and can expand and sorb water and organics in interlayer sites. The most common organic sorption processes on Earth typically involve water or hydroxyl, so hydrated phyllosilicates are good candidates for organic preservation. Using newly derived subsurface hydration results with previously published mineralogy and geochemistry, we derived modeled constraints on the abundances of hydrated amorphous phases, “excess” water, and “excess” cations. These “excess” phases are not accounted for by published crystalline phase abundances or by amorphous phases constrained here. We found correlations between smectites and both “excess” water and “excess” cation abundances, indicating that smectites in Gale crater are hydrated and that cation bridging could be a mechanism for sorption of organics. Our results also show the persistence of amorphous sulfates, opal-A, and volcanic or impact glass, which indicate low water-rock interactions. Increased abundances of sulfates and glass in stratigraphically higher samples may indicate lower water availability and environmental aridification during the time these units were being deposited.

58 GEOSCIENCES

Corundum Discovered by SuperCam and the Perseverance Rover at Jezero Crater, Mars

Minerals provide clues to the type of geologic environments in which they form. The primary types of minerals that have been identified on Mars indicate a history dominated by volcanism followed by interaction with water. However, other minerals are sometimes found. Here we present results from the SuperCam instrument on the Perseverance rover in which the mineral corundum (α-Al2O3) has been identified using a technique called time-resolved luminescence spectroscopy. Corundum has a unique light-emitting (luminescence) feature based on two well-known peaks as well as the time it takes for this light to fade after being excited by a laser, which is referred to as the lifetime. Corundum forms in specific environments that are enriched in aluminum and depleted in silicon and have been subjected to high temperature and pressure. Such places include deep subsurface magmas or shallower locations that have been exposed to high temperature through tectonic processes or meteorite impacts. Given their small size, association with a mineral called plagioclase, and location on the heavily impacted crater rim, we propose that these corundum minerals formed through impact processes.

58 GEOSCIENCES

Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints

Subsurface structure investigation on Mars is crucial for understanding its geological evolution and past hydrological conditions. Elysium Planitia (EP), located near the hypothesized ancient ocean shorelines, could contain clues for past water activity and paleoclimate. Here we present better-constrained subsurface models beneath InSight extending to ~800 m depth, obtained from joint inversion of seismic and seismoacoustic coupling data, and use the well-resolved subsurface structure to explore the lithological profile through rock physics models. The derived subsurface lithology agrees well with local geological context and exhibits a shallow 60-m-thick low-rigidity layer consistent with hydrated sedimentary materials. Despite possible contributions of aeolian and volcanic deposits, we favor the interpretation that the low-rigidity layer originated from fluvial activity in EP during the Hesperian or Hesperian-to-Amazonian epoch, as supported by adjacent paleo-shoreline morphology observations. These results hint at a period of warmer paleoclimate at low latitudes, possibly during high-obliquity phases of Mars’ rotational axis.

58 GEOSCIENCES

The formation and structure of iron-dominated planetesimals

Metal-rich asteroids and iron meteorites are considered core remnants of differentiated planetesimals and/or products of oxygen-depleted accretion. Investigating the origins of iron-rich planetesimals could provide key insights into planet formation mechanisms. Using differentiation models, we evaluate the interior structure and composition of representative-sized planetesimals (~200 km diameter), while varying oxygen fugacity and initial bulk meteoritic composition. Under the oxygen-poor conditions that likely existed early in the inner regions of the Solar System and other protoplanetary disks, core fractions remain relatively consistent across a range of bulk compositions (CI, H, EH, and CBa). Some of these cores could incorporate significant amounts of silicon (10–30 weight%) and explain the metal fractions of Fe-rich bodies in the absence of mantle stripping. Conversely, planetesimals forming under more oxidizing conditions, such as beyond snow lines, could exhibit smaller cores, enriched in carbon, sulfur (>1 wt%), and oxides. Sulfur-rich cores, like those formed from EH and H bulk compositions, could remain partly molten, sustain dynamos, and even drive sulfur-rich volcanism. Additionally, bodies with high carbon contents, such as CI compositions, can form graphitic outer layers. These variations highlight the importance of initial formation conditions in shaping planetesimal structures. Future missions, such as NASA’s Psyche mission, offer an opportunity to measure the relative abundances of key elements (Fe, Ni, Si, and S) necessary to distinguish among formation scenarios and structure models for Fe-rich and reduced planetesimals.

meteorites

Seismic imaging of the Ecuadorian forearc and arc from joint ambient noise, local, and teleseismic tomography: catching the Nazca slab in the act of flattening

SUMMARY The Ecuadorian Andes are a complex region characterized by accreted oceanic terranes driven by the ongoing subduction of the oceanic Nazca plate beneath South America. Present-day tectonics in Ecuador are linked to the downgoing plate geometry featuring the subduction of the aseismic, oceanic Carnegie Ridge, which is currently entering the trench. Using seismic tomography, we jointly invert arrival times of P and S waves from local and teleseismic earthquakes with surface wave dispersion curves to image the structure of the forearc and magmatic arc of the Ecuadorian Andes. Our data set includes > 100 000 traveltimes recorded at 294 stations across Ecuador. Our images show the basement of the central forearc is composed of accreted oceanic terranes with high elastic wave speeds. Inboard of the Carnegie Ridge, the westernmost forearc and coastal cordilleras display relatively low Vp and Vs and high Vp/Vs values, which we attribute to the increased hydration and fracturing of the overriding plate due to the subduction of the thick oceanic crust of the Carnegie Ridge. We additionally image across-arc differences in magmatic architecture. The frontal volcanic arc overlies accreted terranes and is characterized by low velocities and high Vp/Vs indicative of partial melt reservoirs which are limited to the upper crust. In contrast, the main arc displays regions of partial melt across a wider range of depths. The Subandean zone of Ecuador has two active volcanoes built on continental crust suggesting the arc is expanding eastwards. The mid to lower crust does not show indications of being modified from the magmatic process. We infer that the slab is in the process of flattening as a consequence of early-stage subduction of the buoyant Carnegie Ridge.

58 GEOSCIENCES

Evaluation of a generalized least squares algorithm for infrasound beamforming with coherent background noise

Infrasonic signals of interest can occur during periods with persistent, coherent, background noise, which may be natural or anthropogenic. For high signal-to-noise (SNR) ratio transient signals, an ‘overprinting’ of the coherent background may occur, and the signal may still be detected. However, this approach fails for low SNR signals of interest, which may be obscured by coherent noise. An infrasound beamforming method based on generalized least squares (GLS) is investigated for detecting transient signals of interest in the presence of coherent and incoherent background noise. This approach relies on an estimate of the noise covariance, captured in a covariance matrix, to effectively null contributions to the array response from noisy directions of arrival. Synthetic array data is used to investigate the performance of the GLS beamformer compared to the Bartlett beamformer when coherent and incoherent backgrounds are present. Additionally, the effects of array element number and relative strength of the interfering signal on the GLS estimates is investigated. GLS empirical area under the curve estimates suggest that the beamformer can recover coherent power for a signal of interest lower in amplitude than the coherent background, but this effectiveness degrades more quickly with SNR for a four element array compared to a six or eight element infrasound array. Finally, infrasound from the Forensic Surface Experiment, a bolide signal observed at IMS array I37NO, and a volcanic signal recorded at the Alaska Volcano Observatory array ADKI are used to evaluate GLS performance on recorded data. A ten minute window was used to capture the background noise, and the coherent background signal was nulled in all three examples.

58 GEOSCIENCES

Arsenic Accumulation in Microbial Biomass and the Interpretation of Signals of Early Arsenic‐Based Metabolisms

Carbonaceous particles that concentrate arsenic in microbialites as old as ~3.5 Ga are similar to As-rich organic globules in modern microbialites. The former particles have been interpreted as tracers of As cycling by early microbial metabolisms. However, it is unclear if arsenic accumulation is a consequence of biological activity or passive postmortem binding of arsenic by organic matter during diagenesis in volcanically influenced, As-rich environments. Here, we address this uncertainty by evaluating the concentrations, speciation, and detectability of As in active or heat-killed biofilms formed by cyanobacteria or anoxygenic photosynthetic microbes exposed to environmentally relevant concentrations of As(III) or As(V) (50 μM to 3 mM). The genomes or metagenomes of these biofilms contain genes involved in detoxifying or energy-yielding As metabolisms. Biomass accumulates As from the solution in a concentration-dependent manner and with a preference for oxidized As(V) over As(III). Autoclaved biomass accumulates As even more strongly than active biomass, likely because living biofilms actively detoxify As. Active biofilms oxidize and reduce As and accumulate both As(III) and As(V), whereas a small fraction of As(V) can be reduced in inactive biofilms that bind As during diagenesis. Arsenic enrichments in the biomass are detectable by X-ray based spectroscopy techniques (XRF, EPMA-WDS) that are commonly used to analyze geological materials. These findings enable the reconstruction of past active and passive interactions of microbial biomass with arsenic in fossilized microbial biofilms and microbialites from the early Earth.

Madrigal‐Trejo, David

Continental rift evolution and drainage reorganization along the Dead Sea rift since the Miocene

The Dead Sea fault is a section of the Arabian-African plate boundary. Widespread field relations indicate that three major drainage systems (stages) occupied the landscape west of the Dead Sea fault since its initiation at ca. 20 Ma. Specifically, (1) an early to middle Miocene drainage system, only minorly reconfigured by the fault (all sediments of this system belong to the Hazeva Formation); (2) a late Miocene to early Pleistocene fault-parallel drainage system named Paran-Neqarot (all sediments of this system belong to the Arava and Zehiha Formations; and (3) the early Pleistocene to present drainage configuration. The temporal and spatial frameworks of drainage stage 1 are generally constrained by radiometric dating of interfingering volcanic units, and the onset and temporal and spatial frameworks of drainage stage 3 are well constrained by cosmogenic 10 Be surface exposure ages. The timing and longevity of the rift-parallel drainage (stage 2) have until now been evasive to direct dating. The overall time gap between stages 1 and 3 is ~12–13 million years. Thus, an early age (within this time gap) of stage 2 would imply an immediate response of drainage reorganization to rift tectonics, while a later age of this drainage system would imply a delayed response. We present 11 10 Be- 26 Al cosmogenic burial ages of alluvial and colluvial units related to the fault-parallel drainage system (stage 2), which collectively constrain the time of deposition of the Arava Formation sediments in the central Negev to ca. 8 Ma. The general lack of stratigraphic order, together with the large dispersion of ages both across and within the sampling sites, attests to significant recycling of sediments from drainage stage 1 into the Arava Formation deposits. The termination of Arava and Zehiha Formation sediment deposition at ca. 1.8 Ma was determined previously using cosmogenic exposure ages of desert pavements that cover the formations. Combining the previously published data with our new data, we established the longevity and character of the Paran-Neqarot drainage system. In conlcusion, this framework highlights the temporal aspect of drainage system build-up and collapse as it responded to transform and extensional plate boundary tectonics during the Neogene.

58 GEOSCIENCES

Tracking the spatial and temporal links between late Cenozoic extension and magmatism in the Death Valley region (California, USA)

The relationship between late Cenozoic magmatism and extension in the central Basin and Range province (western United States) is complex, necessitating high-precision geochronology to understand its spatiotemporal connections. In the Death Valley region (California), the lack of high-precision U-Pb zircon ages has limited our understanding of the timing of pluton formation and its links to regional extension. We present new high-precision chemical abrasion-isotope dilution-thermal ionization mass spectrometry 206 Pb/ 238 U zircon ages and trace element analyses for eight Death Valley plutons. Our findings reveal three distinct phases of intrusive magmatism: (1) emplacement of shallow rapakivi granites at 13.2 Ma, (2) construction of the mid-crustal Black Mountains intrusive complex at 11.3 Ma, and (3) late emplacement of shallow, compositionally diverse intrusions at 8.2 Ma. A gap in zircon crystallization between 10 Ma and 8.2 Ma coincides with exhumation of the Black Mountains and a transition from sill to dike emplacement. The dominance of rapakivi granites in the Death Valley region, which is rare among Cenozoic granitoids, is likely a result of rapid crustal extension that induces adiabatic decompression. A comparison of the timing of volcanism, plutonism, and tectonic events in Death Valley reveals that intrusive magmatism closely tracks the locus of extension, underscoring the plutonic record as a vital link for understanding regional tectonics and changes in plate boundary dynamics during this period.

58 GEOSCIENCES

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

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

Rodgers, Arthur [Lawrence Livermore National Labor

Multiscale Wildfire Simulation Framework and Remote Sensing

Wildfire as one type of climate extreme events causes huge socioeconomic losses and damages. Large wildfires (i.e., generated pyrocumulonimbus (PyroCb)) can inject tremendous amounts of smoke into the stratosphere, where black carbon and organic carbon aerosols can persist months to years and influence climate by imposing a significant reduction in the radiative forcing like that associated with large volcanic eruptions or proposed via climate interventions such as geoengineering. Both observations and numerical modeling results clearly indicate an increasing trend in wildfire frequency and intensity in many regions during the recent decades with climate change. However, current understanding of wildfire remains largely uncertain owing to limitations of modeling capabilities in representing the multiscale wildfire physics and dynamics and a scarcity of observations constraining important wildfire and environmental variables. This study primarily aims to improve the wildfire simulation capabilities in the state-of-the-art climate model by filling in two major gaps: (1) model resolution is typically too coarse to resolve fine scale processes associated with fires, and (2) chemistry and aerosol processes in fire smoke are poorly represented.

54 ENVIRONMENTAL SCIENCES

3D Geologic Framework Modelling of the Los Alamos National Laboratory Site and Pajarito Plateau: Integrating a realistic 3D fault network and modelling subsurface relationships in a sparsely sampled and complex geologic region

The subsurface geology beneath the Pajarito Plateau is critical to understanding the seismic hazard of the Pajarito Fault System, yet our understanding of this geology is relatively poor. While previous 3D geologic framework models of the area have been created for the purposes of understanding hydrogeologic flow, they are inadequate for the purposes of understanding the Pajarito Fault System. The specific challenges of using oil and gas software for this purpose include: (1) the geologic complexities resulting from volcanism and tectonism; (2) a need for a high level of stratigraphic detail over a large area; (3) a near complete lack of seismic data; and (4) sparse wellbore data. Presented here is a workflow that handles these challenges of adapting commercially available software used by the oil and gas industries to this seismic hazard problem.

58 GEOSCIENCES

A Regional Phase Amplitude Model of 2-D Attenuation for North America

We analyzed seismic attenuation patterns across the North America using over 70,000,000 Lg wave amplitudes recorded by the various networks at frequencies of 0.1-32 Hz. Our inversion solved for laterally varying attenuation, site terms, moments, and apparent stress following Phillips et al., (2016). The inversion was anchored by independently constrained: corner frequencies (via coda spectral ratios) to control the attenuation-stress tradeoff and moment measurements of teleseismic (GCMT, USGS) and regional (St. Louis University and UC Berkeley) earthquakes to provided absolute scaling. The quality factor (Q) shows clear regional patterns: low values in coastal, volcanic, and tectonically active regions, and high values in stable areas like the Great Plains and major plateaus throughout North America. These 2-D Q models enable improved regional source characterization, magnitude estimation, and yield determination.

58 GEOSCIENCES

CLimate Impact: Determining Etiology thRough pAthways (CLDERA)

Climate impacts have broad economic, health, political, and national security ramifications. Societally relevant impacts are typically farther downstream, are the product of multiple interacting processes, and can arise over small regions and timeframes because their sources are short-term and localized. Short-term forcings (as can be seen in volcanic eruptions, climatic tipping points (e.g., the collapse of rainforests or the disappearance of sea ice), or in increasingly plausible climate interventions) fundamentally possess low signal-to-noise and could benefit from accounting for the multiple conditional processes through which a downstream impact arises. Under the Grand Challenge LDRD CLDERA (CLimate impacts: Discovering Etiology thRough pAthways), we have developed tools to enable downstream impact attribution from geographically and temporally localized source forcings in the climate. CLDERA developed methods that can distinguish how a localized source drives the climate system to respond with particular impacts. The how is embodied in pathways – the spatio-temporally evolving chain of physical processes that connects a source to a series of increasingly distant impacts. Novel analytic methods in pursuit of downstream impact attribution were developed and demonstrated on simulations and observations of the 1991 eruption of Mt. Pinatubo in the Philippines. As described within this report we have • developed stratospheric expertise and aerosol modeling capabilities in E3SM, • created original methods to detect and model pathways from source-to-impact, and • advanced climate attribution through novel methods, cases, and approaches. Further, CLDERA developed a tiered verification process consisting of controlled datasets to prototype, verify, and refine the original method development. CLDERA increased Sandia’s footprint in the climate analytics community and developed new climate collaborations whilst also creating a cadre of climate analysts at Sandia. The products from CLDERA have been extensive with a total of 9 journal articles published, 12 articles submitted and under review, and an additional 8 articles in preparation. We have produced 1750 simulated years and developed 9 code-bases. This report details these accomplishments and serves as a summary of the work completed during the CLDERA Grand Challenge.

54 ENVIRONMENTAL SCIENCES

University of Hawai‘i, Shallow Geothermal Resources: Energy Technology Innovation Partnership Project (Final Report)

Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. 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, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

15 GEOTHERMAL ENERGY

Computational Prediction of Infrasound Arrival Times and Directions from Stationary and Moving Impulsive Sources

This report addresses the need to predict infrasound signal arrival times and back azimuths at monitoring stations, enabling more focused and efficient searches within recorded waveform data. The primary challenge is estimating expected signal arrival windows for stationary and moving acoustic sources, such as chemical explosions, volcanic eruptions, meteoroids, and spacecraft re-entry events. To address this challenge, a reproducible methodology is described that uses simplified propagation speeds for boundary layer, tropospheric, stratospheric, and thermospheric atmospheric waveguides. While the Python source code itself is not freely available, this document provides detailed, step-by-step instructions, and equations enabling users to replicate and adapt the method independently. The method reliably predicts signal arrival intervals and back azimuths, thereby supporting rapid detection and accurate interpretation of infrasound events. Results demonstrate that this method effectively identifies plausible signal arrival intervals and directions, facilitating faster event detection and more reliable interpretation. This methodology directly supports atmospheric monitoring, planetary defense, and forensic analysis of explosive atmospheric events.

47 OTHER INSTRUMENTATION

INL CRMO XRF Report 2024-01: Provenance Determinations for 135 Haskett Projectile Points from Eastern Idaho

This report presents provenance determinations for 135 obsidian and fine-grained volcanic (FGV) Haskett projectile points from across eastern Idaho. These determinations were made using energy-dispersive X-ray fluorescence (ED-XRF) spectrometry at the Idaho National Laboratory (INL) Cultural Resource Management Office (CRMO) under the direction of Dr. Kyle Freund.

58 GEOSCIENCES

INL CRMO XRF Report 2025-03: Provenance Determinations for 854 Projectile Points from Idaho Falls District Bureau of Land Management Archaeological Collections

This report presents provenance determinations for 854 obsidian and fine-grained volcanic projectile points from Idaho Falls District Bureau of Land Management collections. Provenance determinations were made using energy dispersive X-ray fluorescence (ED-XRF) spectrometry at the Idaho National Laboratory (INL) Cultural Resource Management Office (CRMO) under the direction of Dr. Kyle Freund.

99 - GENERAL AND MISCELLANEOUS