Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magma”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Geothermal potential of Harrat Rahat, Northern Arabian Shield: geological constraints

Arabian Shield offers a strong potential for geothermal energy. However, only a limited amount of exploration effort has taken place. Harrat Rahat is a recent volcanic area located along the Makkah-Madinah-Nafud volcanic line in western Saudi Arabia. This study is oriented to identify the geothermal potential of Harrat Rahat. The geologic data was collected and processed. Here, we focus on evaluating the potential sources of subsurface heat and groundwater. Subsurface heat is based on the age of erupted volcanoes and the expected volume of subsurface magma chambers. An analogy with an existing geothermal field in a similar setting as Mexico suggests that the faults and the coincident seismicity may be associated with higher permeability. It is concluded that Harrat Rahat is characterized by a good geothermal potential for future geothermal energy source production.

58 GEOSCIENCES↗
DOE OSTI2025

The structure of high-Mg alkali-bearing aluminosilicate glasses investigated in situ at ambient and high pressure by multi-angle energy dispersive X-ray diffraction and infrared microspectroscopy

The structural properties of synthetic high-Mg alkali-bearing aluminosilicate glasses analogues of natural picritic-to-komatiitic magmas were investigated in situ by multiangle energy dispersive X-ray diffraction at 2.1 GPa and ambient pressure and by Fourier Transform infrared spectroscopy up to 5.4 GPa in a cycle of compression and decompression experiments. Our results show that the intermediate range ordering of the glass structure at 2.1 GPa is 3.14 Å, increasing to 3.19 Å when decompressed. The local structure shows T-O lengths of 1.66 Å (2.1 GPa) and 1.65 Å (ambient pressure), T-T distances of 3.19 Å at high pressure, which lengthen to 3.21 Å at ambient pressure, causing the T-O-T angle of 147° determined at 2.1 GPa to widen to 154° upon decompression. The deconvoluted infrared spectra result in the presence of Q 1 , Q 2 , Q 3 populations in the aluminosilicate spectral region, whose proportions remain relatively unchanged up to 5.4 GPa. The structural response of the investigated glasses to cold-compression does not involve changes in polymerization, but rather a shrinking and compaction of the structure as evidenced by the Qn species shifting to higher wavenumbers as a function of pressure. The structural properties determined from X-ray diffraction for this glass composition are discussed together with those of glasses emerging from previous studies to highlight a compositional dependence mainly dictated by the amount of SiO 2 and Al 2 O 3 .

glass structure↗
DOE OSTI2019

Krypton storage capacity of the Earth's lower mantle

Noble gases are important geochemical tracers allowing reconstructing global volatile cycles in Earth's reservoirs. To constrain these fundamental processes, precise data on their partitioning behavior at deep Earth conditions are needed. Such data are only available at moderate pressures up to 25 GPa due to experimental challenges. We have investigated the possibility of noble gas storage in the Earth's lower mantle up to 115 GPa. We studied the incorporation of krypton in the second most abundant lower mantle mineral (Mg 1-x ,Fe x )O (ferropericlase) as well as in liquid metal-alloys by performing experiments up to 115 GPa and 3700 K using the laser-heated diamond anvil cell coupled to post-mortem EMPA analysis and X-ray absorption spectroscopy. The results reveal that, at these extreme conditions, up to 3 wt.% of krypton can be stored in (Mg 1-x ,Fe x )O and 3000 ppm in the Fe-rich liquid metal. For both phases the storage capacities increase with pressure (between 40 GPa and 60 GPa) at a constant high temperature of 2300 K. Fpc has never been considered as a NG host, despite being the second most abundant mineral in the Earth's LM. Using recent accurate compressibility data, we demonstrate that a substitution of krypton into the anion site of (Mg 1-x ,Fe x )O in form of neutral oxygen Schottky defects at diluted lower mantle conditions is possible. This noble gas incorporation mechanism is in agreement with a previous study on bridgmanite. We show that (Mg 1-x ,Fe x )O exhibits higher noble gas storage capacities than bridgmanite through the lower mantle using lattice strain modeling and including experimental solubility and thermoelastic data for neon, argon, krypton and xenon. We also demonstrate that both phases exhibit the highest solubilities for argon and krypton. We used the solubility data from lattice strain modeling to predict noble gas abundances stored in the solid lower mantle after magma ocean crystallization. The modeled abundances show apparent similarities with estimates for the deep noble gas reservoir that are based on either 3He abundances in ocean island basalts or radiogenic 40Ar abundances in the bulk Earth. This strongly indicates that the crystalline lower mantle may play an important role as deep noble gas storage reservoir. We propose, based on considerations on noble gas replenishment from the lower mantle to the atmosphere, that the lower mantle can only contribute to a small fraction of the present-day atmospheric noble gases. This suggests that the lower mantle is an un-degassed reservoir.

58 GEOSCIENCES↗
DOE OSTI2022

Boron isotopes in blue diamond record seawater-derived fluids in the lower mantle

Boron is a quintessential crustal element but its conspicuous presence in diamond – a mantle mineral – raises questions about potential subduction pathways for boron and other volatiles. It has been a long-standing goal to characterize the isotopic composition of boron in blue, boron-bearing (Type IIb) diamonds to reveal its origin. Mineral inclusions indicate that Type IIb diamonds crystallize at transition zone to lower mantle depths, meaning that if the boron is subducted it would trace a pathway of volatile elements into the deep mantle. Here, using off-line laser ablation sampling, we present the first boron isotope compositions, along with trace element contents and carbon isotope compositions, of a suite of blue diamonds mainly sourced from the Cullinan diamond mine in South Africa. The ten analyzed blue diamonds have a wide range in δ 11 B, between -9.2 ± 2‰ and -0.5 ± 2‰, compared to the more restricted range for mid-ocean ridge basalts (-7.1 ± 0.9‰). Carbon isotope values for the blue diamonds range between -20.6 and -1.8‰, with a mode at -17‰, significantly more negative than the main mantle mode at -5‰. Combined, the boron and carbon isotope compositions require fluid input from subducted oceanic lithosphere down into deep mantle source regions of blue diamonds. This finding highlights a deep subduction pathway for volatiles to stoke the source regions of deeply-derived magmas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Gallium isotopic constraints for the origin of the Earth-Moon system

Comparatively heavy isotopic compositions of moderately volatile elements (MVE) in lunar rocks have been advocated to reflect the loss of light isotopes during devolatilization processes from the Moon. In this study we present new gallium (Ga) isotope data for lunar highland rocks, with a focus on the Ferroan Anorthosite Suite (FAS). These are commonly thought to be direct crystallization products from the late lunar magma ocean (LMO) and should contain the majority of the Ga inventory of the Moon. As such, FAS rocks are crucial for identifying the processes that drove Ga isotope fractionation as well as for inferring the Ga isotopic composition of the bulk Moon. Our data reveal that FAS samples have a range in δ 71 Ga from -0.27 to 0.22‰ and are generally isotopically light in Ga compared to other lunar lithologies, but straddle values typical of terrestrial rocks. Although Ga is defined as an MVE, these Ga isotope variations do not correspond with concentrations of more volatile elements, indicating that Ga isotope variations in the FAS are not primarily controlled by devolatilization processes. Instead, the Ga isotopic compositions of bulk FAS rocks broadly correlate with the composition of plagioclase, with the calcium content of plagioclase decreasing as Ga becomes isotopically heavier. This suggests that fractionation of Ga isotopes in FAS rocks was caused by the preferential incorporation of isotopically light Ga into plagioclase during the later solidification stages of the LMO. The progressive crystallization and extraction of plagioclase forces the residual melt towards increasingly heavier Ga isotope ratios, corroborating similar conclusions derived from correlations between δ 71 and Eu* in the mare basalt suite. Using Ga isotope partitioning calculations, we demonstrate that an isotope fractionation coefficient between plagioclase and coexisting melt of -0.3 to -0.4‰ could explain the observed range of δ 71 values in FAS, mare basalt suite rocks, and KREEP. These calculations allow for a first order estimate of the Ga isotopic composition of the bulk silicate Moon prior to plagioclase fractionation and suggest it was close to the composition of the bulk silicate Earth. This would imply that the Moon did not lose a substantial fraction of its Ga inventory during accretion, consistent with new constraints from Rb-Sr isotope systematics that indicate the Moon's volatile deficit was primarily inherited from Theia. In conjunction with the overlap in non-mass-dependent isotope ratios, these collective observations could be reconciled if Theia and the proto-Earth formed in similar regions of the inner Solar System that were already volatile-depleted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Citation for the 2019 Victor M. Goldschmidt award to Donald J. DePaolo

Perhaps the comments most commonly heard when a colleague receives professional recognition are “well-deserved” and “long-overdue”. As such, it is my pleasure to introduce Donald DePaolo as the recipient of this year’s Goldschmidt Award as it is both long-overdue and so well-deserved. While Don has received many honors over the course of his career, the Goldschmidt best recognizes his broad impact on geochemistry. There are few who have contributed as much to geochemistry as Don DePaolo. Finally, his scientific contributions are voluminous and broad, ranging from pioneering efforts to establish the Sm/Nd isotopic system, and its applications to mantle evolution and basalt petrogenesis, to studies of magma mixing, assimilation and crystal fractionation, to evolution of the lower crust and granite petrogenesis, to chemical oceanography, diagenesis and seafloor alteration, to reactive transport accompanying rock-water interaction, and to kinetic isotopic fractionation associated with crystal growth, diffusion and evaporation.

58 GEOSCIENCES↗
DOE OSTI2020

Petrogenesis, alteration, and shock history of intermediate shergottite Northwest Africa 7042: Evidence for hydrous magmatism on Mars?

We report Northwest Africa (NWA) 7042 is an intermediate, permafic shergottite consisting of two generations of zoned olivine (early olivine Fo 41-76 , and late-stage olivine Fo 46-56 ), complexly zoned pyroxene (En 35-64 Fs 22-46 Wo 5-34 ), maskelynitized feldspar and shock melt (An 5-30 Ab 16-61 Or 1-47 ), and accessory merrillite, apatite, ilmenite, titanomagnetite, Fe-Cr-Ti spinels, pyrrhotite, and baddeleyite. The olivine grains have been pervasively altered resulting in orange-brown stained cores with dark-brown margins, surrounded by significantly less altered overgrowth rims. This textural relationship suggests that the alteration of the cores took place at magmatic temperatures prior to crystallization of the rims on Mars. While this type of olivine alteration is rare in shergottites, it is similar to deuterically altered olivine in basalts and gabbros on Earth, caused by residual water-rich magmatic fluids. Launch-generated shock veins in NWA 7042 crosscut and displace several of the altered olivine grains indicating that the alteration occurred before ejection of the meteorite from Mars. However, NWA 7042 has also been subjected to extensive terrestrial weathering, with terrestrial calcite/barite veins occurring throughout the sample. Transmission electron microscopy analysis of the olivine alteration also reveals an absence of high-temperature phases typically associated with hydrous deuteric alteration on Earth, complicating our interpretation. The potential presence of deuterically altered olivine in NWA 7042 has significant implications, as it is the third martian meteorite where deuteric alteration of olivine has been observed (the others being NWA 10416, ALH 77005). The different mantle sources for the parental melts of these three meteorites would suggest many, if not all martian mantle reservoirs have the potential to produce water-rich magmas.

58 GEOSCIENCES↗
DOE OSTI2023

Activity model for 36 elements in Fe-Ni-Si-S-C liquids with application to terrestrial planet accretion and mantle geochemistry: New data for Ru, Re, Pt, Os, Ti, Nb, and Ta

Understanding siderophile element partitioning between metal and silicate melts under diverse conditions can be used to place important constraints on the materials and conditions of planetary accretion and core formation, as well as post core formation processes. However, the effects of Si on the partitioning and activity coefficients for these elements are not well known, despite Si likely being one of the dominant light elements in Earth’s core. To address this gap in understanding, we have undertaken a systematic study of the highly siderophile elements Re, Pt, Os, and Ru, and the refractory lithophile elements Nb, Ta and Ti at 1600 °C and 1 GPa, to derive epsilon interaction parameters for these elements in FeSi metallic liquids. Positive epsilon interaction parameters were measured for Nb, Ta, Ti, Ru, Re, Pt, and Os, indicating that dissolved Si in Fe liquids causes a decrease in their metal/silicate partition coefficients (or ‘silicophobic’ behavior). Furthermore, ε$_{Re, Os, or Ru}^{Si}$ > ε$_{Re, Os, or Ru}^{S}$ which means Si causes a larger decrease in D(metal/silicate) than S, and the chalcophile behavior expected from some elements will be completely masked by the presence of Si in a metallic liquid. The new parameters are used to update an activity model that now includes 36 siderophile elements in Fe-Ni-Si-S-C liquids (27 trace elements considered here). Systematic assessment of these 27 elements shows which have the strongest affinity for Si, C, and S, and also how activity coefficients for these elements would vary during accretion and core formation in Earth, Mars, and Mercury of widely differing fO 2 and core compositional conditions. The activity model is combined with new partitioning expressions for Mo, W, Cr, Re Ru, Pt, and Os and applied to aspects of post core formation mantle geochemistry of Earth, Mars, and Mercury. Our updated expressions show that the BSE Mo/W ratio can easily be achieved with metal/silicate partitioning during growth of the Earth, whereas Re, Os and Ru become lower than and highly fractionated compared with BSE values during core formation and accretion, and thus nearly 99% of their BSE abundances are likely contributed by late accretion. Ru isotopes should be a very good indicator of the source material for the late accretion. The high Pt/Os and Re/Os developed in a deepening magma ocean during the growth of the Earth, indicates 186 Os and 187 Os isotopes could be coupled if this ancient material remained isolated and subsequently became entrained in mantle plumes and measured in surficial lavas. The extent to which this occurred will be limited by the low Os content of this ancient material, thus requiring mixing as a major component in plume sources. Martian mantle Hf/W ratio stays low during accretion and core formation modelling, suggesting that W isotope anomalies are more likely due to solid/liquid silicate fractionation than to core formation. Finally, Ti contents measured by MESSENGER at Mercury’s surface can be explained by segregation of either a metallic core (IW-6 to -8) or metallic core + sulfide (IW-4 to -7.5) followed by mantle melting.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Chronology of ungrouped achondrites Northwest Africa 11119 and 7325/8486: implications for early evolution of a heterogeneous crust on a differentiated planetesimal

The chronology of ungrouped achondrites provides key insights into the timeline of igneous activity in the early Solar System. Two ungrouped achondrites – Northwest Africa (NWA) 11119 (andesite-dacite) and NWA 8486 and its pair NWA 7325 (olivine gabbro) – are the focus of this study, where their chronologies are investigated using multiple high-resolution techniques. Here we report the lead-lead (207Pb-206Pb) and manganese-chromium (53Mn-53Cr) systematics of NWA 11119, as well as the 207Pb-206Pb systematics of NWA 8486 alongside aluminum-magnesium (26Al-26Mg) systematics for NWA 7325. The U-corrected 207Pb-206Pb ages of NWA 11119 and the combined ages of NWA 7325/8486 are 4566.4 ± 0.8 Ma and 4563.8 ± 1.9 Ma, respectively. Additionally, we report the 53Mn-53Cr age of NWA 11119 to be 4564.4 ± 2.5 Ma and the 26Al-26Mg age of NWA 7325/8486 to be 4563.1 ± 0.3 Ma. The formation of NWA 11119 requires partial melting of a (likely chondritic) source reservoir leading to eruption of Si rich, alkali depleted magmas, while NWA 7325/8486 likely formed from a chemically fractionated reservoir with superchondritic Al/Mg. The clear geochemical and isotopic similarities of these achondrites, combined with the chronology reported here, is suggestive of formation of these two ungrouped achondrites on a common parent body which likely formed in the inner Solar System and experienced early differentiation under reducing conditions. If these achondrites did share a parent body, it would suggest that primary asteroids in the early Solar System commonly produced mineralogically and geochemically heterogeneous crusts. While mineralogical and geochemical heterogeneity is known to exist in the (mostly mafic) crust of asteroid Vesta, our findings show that even more significant crustal heterogeneity (representing felsic and mafic compositions) may exist on other asteroids.

Cartwright, Julia [University of Leicester, UK]↗
DOE OSTI2024

Time-dependent thermal degradation of lost circulation materials in geothermal systems

Treatment of lost circulation can represent anywhere from 5 to 25 % of the cost in drilling geothermal wells. The cost of the materials used for lost circulation treatment is less important than their effectiveness at reducing fluid losses. In geothermal systems, the high temperatures (>90 °C) are expected to degrade many commonly used lost circulation materials over time. This degradation could compromise different materials ability to mitigate fluid loss, creating more non-productive time as multiple treatments are needed, but may result in recovering desired permeability zones within the reservoir section over time. This research aimed to study how thermal degradation of eight different lost circulation materials affected their properties relevant to sealing loss zones in geothermal wells. Here, mass loss experiments were conducted with each material at temperatures of 90–250 °C for 1–42 days to measure the breakdown of the material at geothermal conditions, collecting gases during several experiments to determine the waste produced during degradation. Compaction experiments were conducted with the degraded materials to show how temperatures reduced the rigidity and increased packing of the materials. Viscosity tests were conducted to show the impact of different materials on drilling fluid rheology. Microscope observations were conducted to characterize the alterations to each material due to thermal degradation. Organic materials tend to degrade more than inorganic materials, with organics like microcellulose, cotton seed hulls and sawdust losing 30–50 % of their mass after 1 day of heating at 200 °C, while inorganics like magma fiber only lose ~5–10 % of its mass after one day of heating at 200 °C. Granular materials are the strongest when compacted despite any mass loss, while fibrous and flaky materials are fairly weak and breakdown easily under stress. The materials do not generally affect fluid rheology unless they have a viscosifying agent as part of the mixture. Microscopic analysis showed that more rigid materials like microcellulose and cedar fiber degrade in brittle manners with splitting and fracturing, while others like cotton seed hulls degrade in more ductile manners forming meshes or clumps of material. The thermal breakdown of lost circulation materials tested suggests that each material should also be classified by its degree of thermal degradability, as at certain temperatures the materials can lose the capability to bridge loss zones around the wellbore.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Magmatic, hydrothermal and ore element transfer processes of the southeastern Archean Superior Province implied from electrical resistivity structure

Although magmatic- and metamorphic- derived fluids are widely recognized mineralizing agents, the role of crustal architecture in defining source and sink zones within the middle to lower crust and upper mantle of ancient orogens remains enigmatic. The globally largest and best-preserved Archean greenstone belts lie in the Superior Province, Canada. They provide an ideal location to investigate the influence of igneous construction and subsequent syn-deformational plutonism and metamorphism on the localization of metal-rich melts and fluids throughout the crustal column. Integration of three-dimensional magnetotelluric modelling and seismic reflection sections across the Abitibi subprovince reveals details of a 'whole-of-crust' magmatic and hydrothermal system. East-west low resistivity structures broadly underlie the surface traces of the major deformation zones that are host to significant gold endowment (>200 Moz), while mid-crustal cross trends suggest mineralized fluids flowed along 'pipes' within fault planes. Most low resistivity structures are inferred to represent domains containing interconnected zones of graphite and/or sulfide. Furthermore, these delineate relict mantle source/transit domains and crustal pathways enriched by the flow of magmas or metamorphic fluids genetically related to a late-stage pulse of ore-bearing magmatism, possibly as a result of slab break-off or delamination. Using the combined MT and seismic data, we develop a 3-D crustal-scale model which highlights how evolving orogenic architecture-controlled mass transfer and metallogenic processes developed.

58 GEOSCIENCES↗
DOE OSTI2022

Development towards stable chlorine isotope measurements of astromaterials using the modified Middleton source of an accelerator mass spectrometer

The volatile element Cl can be lost during the formation and evolution of planetary bodies, leading to fractionation of its two stable isotopes 35 Cl and 37 Cl. Chlorine isotope variations (reported as δ 37 Cl in parts per thousand (‰) relative to Standard Mean Ocean Chloride, SMOC) are documented to exceed 80‰ between different lunar rock samples and have been variably interpreted as the fingerprint of degassing during accretion, magma ocean, or volcanic portions of lunar history. The large intersample and intrasample variations observed by both bulk isotope ratio mass spectrometry (IRMS) and in-situ secondary ion mass spectrometry (SIMS) methods are difficult to interpret in part because of a paucity of bulk Cl isotope measurements. This lack of high-precision bulk data is due to the relative rarity of IRMS laboratories capable of making these high precision measurements on small samples of precious planetary materials such as those returned by human or robotic exploration. Here we present a new method for performing high precision δ 37 Cl measurements using the modified Middleton ion source of an existing accelerator mass spectrometer. For samples with as little as 1 μg Cl–the equivalent of 2–4 mm 3 of a typical lunar rock sample–the average cathode accuracy is ~1‰. Finally, cathode reproducibility is typically ~1‰ (2σ) for samples with at least 10 μg of Cl, increasing to ~3–6‰ for aliquots with ~1–2 μg Cl, similar to published SIMS results and sufficient to study astromaterials from the Moon, Mars, or 4 Vesta, which have tens of ‰ observed variations.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

In-situ investigation of the vibrational properties of H 2 O CO 2 -bearing and dry K-rich basaltic glasses at high pressure by mid infrared spectroscopy

Here, the vibrational properties of CO 2 –H 2 O-bearing and dry synthetic K-rich basaltic glasses were investigated at room temperature and pressures between 0.0001–5.5 GPa using the diamond anvil cell combined with in situ reflectance and transmittance Fourier Transform infrared micro-spectroscopy. The absorption coefficient a(v) calculated from the Kramers-Kronig relation shows that glasses are dominated by the Q 2 aluminosilicate unit followed by Q 1 , Q 3 . The variation in Q n concentration upon compression suggests that glasses undergo polymerization from ambient pressure to 2 GPa followed by less marked structural changes up to 4 GPa, above which the structure is further polymerized. Once decompressed, glasses retained a polymerized structure. Our results show that 1.03 wt% CO 2 -1.42 wt% H 2 O in a glass with 3.10 wt% K 2 O and 46.77 wt% SiO 2 prevent the formation of fully polymerized connections when cold-compressed. Our results can be used to explain the effect of volatiles and alkali on the rheology of natural basaltic magmas.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Recent seismicity at Ceboruco Volcano (Mexico)

Magma and related hydrothermal fluid movement, and their interaction with solid rock, in active volcanic regions, generate a wide variety of seismic waves whose characterization can mitigate the risk of a potential eruption. Located in the western region of the Trans-Mexican Volcanic Belt, Ceboruco Volcano, whose last eruptive period was 1870–1875, is considered to be one of the most hazardous volcanoes in Mexico. We have conducted a detailed study of the seismicity in the surroundings of Ceboruco's volcanic edifice to assess the current state of this volcano. A dense temporary seismic network with 25 seismic stations in an area of 16 km × 16 km was deployed between November 2016 and July 2017, as part of the P-24 project of the CeMIEGeo consortium; this effort has allowed the detection of 81 earthquakes concentrated beneath the crater with depths between 4 and 8 km. In this study, we observe that the recorded seismicity occurs in swarms, and we specifically identify four sequences that we characterize in detail via the first focal mechanisms available for this volcano. Our results suggest a change in the local seismicity distribution compared to earlier observations, which reported seismic activity near the volcano edifice associated with fluid migration along zones of weakness related to the extensional stresses of the Tepic-Zacoalco rift. The changes in seismic patterns and obtained focal mechanisms are consistent with observed fluid effects at many geothermal sites worldwide, but also could suggest resumption of activity at this currently dormant volcano.

58 GEOSCIENCES↗
DOE OSTI2021

GPU algorithms for Efficient Exascale Discretizations

In this paper we describe the research and development activities in the Center for Efficient Exascale Discretization within the US Exascale Computing Project, targeting state-of-the-art high-order finite-element algorithms for high-order applications on GPU-accelerated platforms. Furthermore, we discuss the GPU developments in several components of the CEED software stack, including the libCEED, MAGMA, MFEM, libParanumal, and Nek projects. We report performance and capability improvements in several CEED-enabled applications on both NVIDIA and AMD GPU systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High-resolution lithospheric structure of continental China from joint inversion of surface wave and gravity data

The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.

58 GEOSCIENCES↗
DOE OSTI2021

Uplift and Seismicity Driven by Magmatic Inflation at Sierra Negra Volcano, Galápagos Islands

Abstract Although episodes of surface uplift and elevated seismicity precede many volcanic eruptions, their temporal evolution is often complex, and apparently in contradiction to simple trends predicted by mechanical deformation models. Here, we use continuous global positioning system and seismic data recorded at Sierra Negra volcano, Galápagos Islands, to show how the edifice responded to stress changes driven by magma accumulation in a shallow sill. The rate of uplift varied during the 13 years and 6.5 m of inflation before the 2018 eruption. The number of earthquakes per unit of uplift increased exponentially with total uplift as the differential stress increased. Accordingly, the temporal seismicity rate varied in time as a function of both the total uplift and the uplift rate. The Gutenberg‐Richter b ‐value decreased as a function of total uplift. In the final six months before the eruption, a sequence of large ( M > 4) earthquakes regulated the state of stress on the fault, each being followed by 2–3 days of postseismic quiescence, and retarding the increase in seismicity rate. These earthquakes did not affect the overall uplift rate. Subsidence of 8.5 m accompanied the 2‐month eruption. On resumption of uplift, the number of earthquakes per unit of uplift was very low, and the b ‐value high, reflecting the relaxed stress state of the fault system. These observations show that crustal deformation becomes increasingly brittle at higher stress states, and supports theoretical models based on elastic‐brittle mechanics. They suggest that joint interpretation of deformation and seismicity is key for forecasting future eruptions in similar volcanic settings.

Bell, Andrew Forbes↗
DOE OSTI2023

Lightning Rings and Gravity Waves: Insights Into the Giant Eruption Plume From Tonga's Hunga Volcano on 15 January 2022

Abstract On 15 January 2022, Hunga Volcano in Tonga produced the most violent eruption in the modern satellite era, sending a water‐rich plume at least 58 km high. Using a combination of satellite‐ and ground‐based sensors, we investigate the astonishing rate of volcanic lightning (>2,600 flashes min −1 ) and what it reveals about the dynamics of the submarine eruption. In map view, lightning locations form radially expanding rings. We show that the initial lightning ring is co‐located with an internal gravity wave traveling >80 m s −1 in the stratospheric umbrella cloud. Buoyant oscillations of the plume's overshooting top generated the gravity waves, which enhanced turbulent particle interactions and triggered high‐current electrical discharges at unusually high altitudes. Our analysis attributes the intense lightning activity to an exceptional mass eruption rate (>5 × 10 9 kg s −1 ), rapidly expanding umbrella cloud, and entrainment of abundant seawater vaporized from magma‐water interaction at the submarine vent.

58 GEOSCIENCES↗