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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Magnetotelluric Investigations of the Kilauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. Data was acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kilauea. Melt fractions between 0.096 and 0.117 are calculated for the Kilauea and Pu'u'ō'ō low resistivity anomalies, respectively.

58 GEOSCIENCES↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2024

The Global Biogeochemical Cycle of Rhenium

Here, this paper is the first comprehensive synthesis of what is currently known about the different natural and anthropogenic fluxes of rhenium (Re) on Earth's surface. We highlight the significant role of anthropogenic mobilization of Re, which is an important consideration in utilizing Re in the context of a biogeochemical tracer or proxy. The largest natural flux of Re derives from chemical weathering and riverine transport to the ocean (dissolved = 62 × 10 6 g yr -1 and particulate = 5 × 10 6 g yr -1 ). This review reports a new global average [Re] of 16 ± 2 pmol L -1 , or 10 ± 1 pmol L -1 for the inferred pre-anthropogenic concentration without human impact, for rivers draining to the ocean. Human activity via mining (including secondary mobilization), coal combustion, and petroleum combustion mobilize approximately 560 × 10 6 g yr -1 Re, which is more than any natural flux of Re. There are several poorly constrained fluxes of Re that merit further research, including: submarine groundwater discharge, precipitation (terrestrial and oceanic), magma degassing, and hydrothermal activity. The mechanisms and the main host phases responsible for releasing (sources) or sequestrating (sinks) these fluxes remain poorly understood. This study also highlights the use of dissolved [Re] concentrations as a tracer of oxidation of petrogenic organic carbon, and stable Re isotopes as proxies for changes in global redox conditions.

58 GEOSCIENCES↗
DOE OSTI2024

Oxidized Sulfur Species in Slab Fluids as a Source of Enriched Sulfur Isotope Signatures in Arcs

Recycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO 2 ) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ 34 S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H 2 O-NaCl solutions to calibrate a high-pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at log fO 2 from FMQ-2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite at fO 2 conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur-rich and high δ 34 S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke 34 S-rich subducted lithologies.

58 GEOSCIENCES↗
DOE OSTI2025

Sound Velocity and Grüneisen Parameter in Shock‐Melted Silica at Deep Earth Pressures

Silica is a primary component of rocky planet interiors and its melt properties are important for understanding planetary formation and differentiation, magma oceans, and the deep mantle. Although well understood in the solid state, the high-pressure behavior of liquid silica is poorly constrained at lower mantle pressures. Using laser interferometry to measure shock wave profiles, we report measured stress-density states and longitudinal sound speeds in shock-synthesized stishovite, from fused silica staring material, and across the solid-liquid phase boundary up to 154 GPa. Our results constrain completion of melt at 80 GPa and show that at pressures relevant to the deep mantles of Earth-sized, rocky planets, the Grüneisen parameter for liquid silica increases with compression. This finding is consistent with a continuous increase in Si-O coordination above six for liquid silica at core-mantle boundary relevant pressures and temperatures.

Ocampo, I. K. [Princeton University, NJ (United St↗
DOE OSTI2026

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↗
DOE OSTI2024

Full-waveform tomography reveals iron spin crossover in Earth’s lower mantle

Abstract Three-dimensional models of Earth’s seismic structure can be used to identify temperature-dependent phenomena, including mineralogical phase and spin transformations, that are obscured in 1-D spherical averages. Full-waveform tomography maps seismic wave-speeds inside the Earth in three dimensions, at a higher resolution than classical methods. By providing absolute wave speeds (rather than perturbations) and simultaneously constraining bulk and shear wave speeds over the same frequency range, it becomes feasible to distinguish variations in temperature from changes in composition or spin state. We present a quantitative joint interpretation of bulk and shear wave speeds in the lower mantle, using a recently published full-waveform tomography model. At all depths the diversity of wave speeds cannot be explained by an isochemical mantle. Between 1000 and 2500 km depth, hypothetical mantle models containing an electronic spin crossover in ferropericlase provide a significantly better fit to the wave-speed distributions, as well as more realistic temperatures and silica contents, than models without a spin crossover. Below 2500 km, wave speed distributions are explained by an enrichment in silica towards the core-mantle boundary. This silica enrichment may represent the fractionated remains of an ancient basal magma ocean.

58 GEOSCIENCES↗
DOE OSTI2023

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

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

A fast, dense Chebyshev solver for electronic structure on GPUs

Matrix diagonalization is almost always involved in computing the density matrix needed in quantum chemistry calculations. In the case of modest matrix sizes (≲4000), performance of traditional dense diagonalization algorithms on modern GPUs is underwhelming compared to the peak performance of these devices. This motivates the exploration of alternative algorithms better suited to these types of architectures. We newly derive, and present in detail, an existing Chebyshev expansion algorithm whose number of required matrix multiplications scales with the square root of the number of terms in the expansion. Focusing on dense matrices of modest size, our implementation on GPUs results in large speed ups when compared to diagonalization. Additionally, we improve upon this existing method by capitalizing on the inherent task parallelism and concurrency in the algorithm. Furthermore, this improvement is implemented on GPUs by using CUDA and HIP streams via the MAGMA library and leads to a significant speed up over the serial-only approach for smaller (≲1000) matrix sizes. Finally, we apply our technique to a model system with a high density of states around the Fermi level, which typically presents significant challenges.

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

The carbon content of Earth and its core

Earth’s core is likely the largest reservoir of carbon (C) in the planet, but its C abundance has been poorly constrained because measurements of carbon’s preference for core versus mantle materials at the pressures and temperatures of core formation are lacking. Using metal–silicate partitioning experiments in a laser-heated diamond anvil cell, we show that carbon becomes significantly less siderophile as pressures and temperatures increase to those expected in a deep magma ocean during formation of Earth’s core. Based on a multistage model of core formation, the core likely contains a maximum of 0.09(4) to 0.20(10) wt% C, making carbon a negligible contributor to the core’s composition and density. However, this accounts for ~80 to 90% of Earth’s overall carbon inventory, which totals 370(150) to 740(370) ppm. The bulk Earth’s carbon/sulfur ratio is best explained by the delivery of most of Earth’s volatiles from carbonaceous chondrite-like precursors.

Geosciences↗
DOE OSTI2020

In situ X-ray diffraction of silicate liquids and glasses under dynamic and static compression to megabar pressures

Properties of liquid silicates under high-pressure and high-temperature conditions are critical for modeling the dynamics and solidification mechanisms of the magma ocean in the early Earth, as well as for constraining entrainment of melts in the mantle and in the present-day core–mantle boundary. In this paper, we present in situ structural measurements by X-ray diffraction of selected amorphous silicates compressed statically in diamond anvil cells (up to 157 GPa at room temperature) or dynamically by laser-generated shock compression (up to 130 GPa and 6,000 K along the MgSiO 3 glass Hugoniot). The X-ray diffraction patterns of silicate glasses and liquids reveal similar characteristics over a wide pressure and temperature range. Beyond the increase in Si coordination observed at 20 GPa, we find no evidence for major structural changes occurring in the silicate melts studied up to pressures and temperatures exceeding Earth’s core mantle boundary conditions. This result is supported by molecular dynamics calculations. Our findings reinforce the widely used assumption that the silicate glasses studies are appropriate structural analogs for understanding the atomic arrangement of silicate liquids at these high pressures.

58 GEOSCIENCES↗
DOE OSTI2022

A dry ancient plume mantle from noble gas isotopes

Primordial volatiles were delivered to terrestrial reservoirs during Earth’s accretion, and the mantle plume source is thought to have retained a greater proportion of primordial volatiles compared with the upper mantle. This study shows that mantle He, Ne, and Xe isotopes require that the plume mantle had low concentrations of volatiles like Xe and H 2 O at the end of accretion compared with the upper mantle. A lower extent of mantle processing alone is not sufficient to explain plume noble gas signatures. Ratios of primordial isotopes are used to determine proportions of solar, chondritic, and regassed atmospheric volatiles in the plume mantle and upper mantle. The regassed Ne flux exceeds the regassed Xe flux but has a small impact on the mantle Ne budget. Pairing primordial isotopes with radiogenic systems gives an absolute concentration of 130 Xe in the plume source of ∼1.5 × 10 7 atoms 130 Xe/g at the end of accretion, ∼4 times less than that determined for the ancient upper mantle. A record of limited accretion of volatile-rich solids thus survives in the He-Ne-Xe signatures of mantle rocks today. A primordial viscosity contrast originating from a factor of ∼4 to ∼250 times lower H 2 O concentration in the plume mantle compared with the upper mantle may explain (a) why giant impacts that triggered whole mantle magma oceans did not homogenize the growing planet, (b) why the plume mantle has experienced less processing by partial melting over Earth’s history, and (c) how early-formed isotopic heterogeneities may have survived ∼4.5 Gy of solid-state mantle convection.

58 GEOSCIENCES↗
DOE OSTI2025

A tale of two planets: Disparate evolutionary models for Mars inferred from radiogenic isotope compositions of Martian meteorites

The radiogenic isotopic compositions of basaltic Martian meteorites (shergottites) and clinopyroxene/olivine cumulate meteorites (nakhlite/chassignites) are used to define the global evolution of Mars. However, the two main groups of meteorites demonstrate that their sources underwent divergent styles of magmatic evolution. The shergottites portray a planet that differentiated ~4.52 billion years ago via solidification of a magma ocean, producing incompatible element-depleted and -enriched reservoirs that remained isolated until melt production. In contrast, the reservoir from which the nakhlite/chassignites derive may have formed earlier, produced melts that fractionated Sm/Nd and Hf/W differently, was compositionally less variable, and experienced a significantly more complex history following primordial differentiation than the shergottite sources. The disparate histories recorded by these two groups of meteorites elucidate important questions that could be addressed by acquiring additional samples. Obtaining samples that shared the isotopic systematics of the shergottites would provide confidence that extrapolating the primordial differentiation history of Mars from shergottite radiogenic isotope systematics is reasonable. Returned samples from Mars will also constrain the physical locations of the meteorite source regions, providing insights into the general structure of the Martian mantle. In addition, they will help constrain the phases present in the martian mantle during melting and the conditions under which they are stable. Finally, identifying an evolved lithology that satisfies the geochemical and isotopic constraints placed on the incompatible element-enriched endmember observed in the shergottites would define the nature of magmatic evolution on Mars and whether it is more akin to processes on the Earth or the Moon.

Mars↗
DOE OSTI2022

Crustal and uppermost mantle structure of Cape Verde from ambient noise tomography

SUMMARY We present a seismic ambient noise tomography of the Cape Verde archipelago, located in the Atlantic Ocean, approximately 600 km west of Senegal. We used 38 seismic broad-band stations that continuously recorded for 10 months, in order to construct the first 3-D model of Sv-wave velocities for the crust and uppermost mantle beneath the Cape Verde region. We started by computing phase cross-correlations for vertical component recordings using all possible inter-island station pairs. Next, a time–frequency phase-weighted stack was applied to obtain robust Rayleigh-wave group-velocity dispersion curves in the period band between 10 and 24 s. Group-velocity maps at different periods are obtained by inverting the dispersion curves. We then inverted the group-velocity maps to obtain the 3-D shear wave velocity structure of the crust and uppermost mantle beneath Cape Verde. The final 3-D model extends from 8 km down to 23 km and has a lateral resolution of about 50 km. The crust in the southwestern sector, encompassing Fogo, presents lower S-wave velocities that may be caused by the presence of melt pockets and/or hydrothermal fluids circulation. The uppermost mantle beneath the northwestern sector is characterized by higher S-wave velocities in agreement with previous results obtained from Ps and Sp receiver functions. Those high-velocity anomalies can reflect non-altered crust or remnants of magma chambers or solidified basaltic intrusions, which fed the volcanism in these islands. Our maps revealed the presence of crustal underplating across the entire archipelago, yet stronger beneath the groups Santo Antão—São Vicente—São Nicolau and Fogo—Santiago—Maio.

Carvalho, J. (ORCID:0000000310778994)↗
DOE OSTI2023

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

Pinzón, Juan I. (ORCID:0000000338078068)↗
DOE OSTI2021

Screening of Lost Circulation Materials for Geothermal Applications: Experimental Study at High Temperature

Abstract This study presents a laboratory experimental research to determine the characteristics of lost circulation materials (LCMs) capable of addressing thermal degradation, providing bridging and sealing in geothermal conditions. Eleven different materials were tested: Walnut Fine, Walnut Medium, Sawdust, Altavert, Graphite Blend, Bentonite Chips, Micronized Cellulose (MICRO-C), Magma Fiber Fine, diatomaceous earth/amorphous silica powder (DEASP), Cotton Seed Hulls, and a Calcium Carbonate Blend. The filtration and sealing pressure of the LCMs were measured with HPHT equipment up to 149 °C (300 °F). Besides, the particle size distribution (PSD) of fine granular materials was measured. The results show that the performance of some LCM materials commonly used in geothermal operations is affected by high temperature. Characteristics such as shape and size made some materials more prone to thermal degradation. Also, it was found that the PSD of LCMs is a key factor in the effectiveness of bridging and sealing fractures. The results suggest that granular materials with a wide particle size distribution PSD are suitable for geothermal applications.

Energy & Fuels↗
DOE OSTI2025

Elasticity of davemaoite as a primary contributor to lower-mantle heterogeneities

Geophysical detection of subducted mid–ocean ridge basalt (MORB) in the lower mantle is hindered by uncertainties in the elasticity of Fe,Al,Mg,Ti–bearing davemaoite, a key MORB component. Using Brillouin spectroscopy and x-ray diffraction, we determined the elasticity of a Ca 0.906(1) Fe 2+ 0.027(1) Fe 3+ 0.042(1) Mg 0.033(1) Al 0.072(1) Ti 0.020(1) Si 0.912(1) O 3 davemaoite up to 113 gigapascals and 2294 K. Here, we found that it exhibited a shear wave velocity 10 to 20% slower than end-member davemaoite, making it the slowest phase among major lower-mantle minerals. Our models show that MORB, containing 20 to 25 volume percent davemaoite, potentially contributes to large low-shear-velocity provinces (LLSVPs), whereas a cumulate layer enriched in davemaoite crystallized from basal magma ocean may comprise ultralow-velocity zones (ULVZs). Davemaoite’s ability to host incompatible and heat-producing elements possibly links LLSVPs and ULVZs to mantle plume initiation and geochemical signatures of ocean island basalts.

Zhou, Wen-Yi [Texas A & M Univ., College Station, ↗
DOE OSTI2024

Timescales of mafic magmatic fractionation documented by paleosecular variation in basaltic drill core, Snake River Plain volcanic province, Idaho, USA

Abstract The timescales over which fractional crystallization and recharge work in mafic volcano-plutonic provinces is subject to great uncertainty. Currently modeled processes are subject to the scale of measurement: monogenetic basaltic fields accumulate over hundreds of thousands of years, consistent with U-Th-Ra isotopic variations that imply 50% crystallization of basic magmas on timescales of 100,000 years or more, whereas crystal diffusion modeling implies phenocryst residence times of ~1–1000 years. Monogenetic basalts of the Snake River Plain in southern Idaho, USA, are up to 2 km thick and postdate passage over the Yellowstone–Snake River Plain hotspot. Detailed lithologic and geophysical logging of core from deep drill holes, along with chemical stratigraphy and high-resolution paleomagnetic inclination measurements, document individual eruptive units, compound lava flows, and basaltic flow groups that accumulated over 1–6 m.y. Hiatuses are commonly marked by loess or fluvial interbeds that vary from ~0.1 m thick to 20 m thick. Radiometric (40Ar-39Ar, detrital zircon U-Pb) and paleomagnetic timescale ages show that the deepest hole (Kimama drill hole, 1912 m total depth) accumulated over ~6 m.y. Cycles of fractional crystallization and recharge are recognized in the chemical stratigraphy as up-section shifts in major and trace elements; these fractionation cycles commonly represent 40%–50% fractionation. Individual fractionation cycles may comprise 20–40 eruptive units (8–17 lava flows) with little to no change in paleomagnetic inclination (0°–1°), whereas adjacent cycles may differ by several degrees from one another or reflect changes in polarity. Rates of paleosecular variation in Holocene lavas and sediments dated using 14C document significant shifts in magnetic inclination over short timescales, ranging from ~0.05° to 2°/decade, with an average of ~0.5°/decade and a minimum rate of 0.05°/decade. This implies that fractionation cycles with ≤1° variation in magnetic inclination formed on timescales of a few decades up to a few centuries (20–200 years). Thus, the lavas collectively represent only a few thousand years of eruptive activity, with major flow groups separated in time by tens to hundreds of thousands of years. We suggest that the rates defined by paleosecular variation capture the timescales of magmatic chamber evolution (fractionation/recharge) in the seismically imaged mid-crustal sill complex; in contrast, we suggest that crystal diffusion modeling captures the residence times in shallow subvolcanic magmatic chambers that underlie individual monogenetic volcanoes.

Geology↗