SEARCH · Engineering Papers

Results for “variational”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Agent-based simulations of shared automated vehicle operations: reflecting travel-party size, season and day-of-week demand variations

Here, this paper explores the effects of day of week and season of year demand variations for shared rides, along with realistic travel party sizes, on shared autonomous vehicle (SAV) services across the Austin, Texas region. Using the agent-based POLARIS program, synthetic person-trips that reflect travel-party size (from one to four persons) and demand variations over days and months, as evident in the National Household Travel Survey data were simulated in each scenario over a 24 h travel day. Results show that realistic party sizes can bring considerable changes to SAV fleet performance, including up to 8.5% higher service rates (number of requests accepted within 15 min), 5 min shorter journey times (wait time + travel time), 28% higher vehicle occupancies on weekends, and roughly 4% lower empty fleet VMT. Weekend travel is most impacted by season of year, with weekday travel patterns looking more uniform (thanks to work and school trips). Various performance metrics for the Austin network, like total and empty VMT, change by up to 30% when considering realistic variations in party size and time of year. This paper underscores the value of recognizing day-to-day and month-to-month variations in travel demand, and the importance of agent-based model equations to reflect travel-party size. Such realism can help quantify SAV seat occupancies more accurately, highlighting the importance of shared mobility. However, it also creates demand and supply issues for operators that now need more information on party size to manage dynamic ride-sharing, or those that may wish to shift their fleet vehicles to other regions for special events to protect profits while offering reasonable wait times to customers throughout the year.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Significance of slight ambient temperature variation on the electrocatalyst performance toward oxygen evolution reaction

Hundreds of oxygen evolution reaction (OER) electrocatalysts have been developed over the past few decades, and their performances are evaluated and compared at ambient temperature. However, the effect of ambient temperature variation on OER electrocatalyst performance has received less attention, which may play a remarkable role in the electrocatalytic activity. In this work, we systematically investigated the influence of ambient temperature variation on electrocatalyst performance toward OER. The results show that the slight ambient temperature variation has a significant effect on OER catalyst performance based on the changes of overpotential (10mAcm -2 ) and Tafel slope. Both remarkable chances are observed on transition metal (Ni) and noble metal (IrO 2 ) electrocatalysts, and the overpotentials decrease around 81 mV with a temperature increase by 20°C (from 10° to 30°C) for both Ni and IrO 2 electrocatalysts with the Tafel slope drops of 36.9 and 29.5 mV dec -1 , respectively. Further, a similar trend is also found in the electrochemically active surface (ECSA) normalized performance and the charge transfer resistance. This study demonstrates that reporting the actual operating temperature for OER is not only recommended but also necessary to evaluate and compare electrocatalyst activities from different materials systems properly, and neglecting the ambient temperature variation effect can highly mislead conclusions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Carbonated apatite lattice parameter variation across incremental growth lines in teeth

Here, results of microbeam mapping (x-ray excited x-ray fluorescence and transmission x-ray diffraction) of Beluga whale cementum and of beaver dentin reveal growth bands are defined by variation in mineral content (carbonated hydroxyapatite, cAp), by variation in Zn content and by variation in cAp lattice parameter a. For both tissues, diffraction and fluorescence data show that Zn signal and a are anti-correlated, that is high Zn signal is co-located with low a, and vice versa. In cementum and bone, treatment with acid removes Ca, Zn and cAp. The data are strong evidence that Zn is incorporated into the cAp lattice in these tissues and that the amount of Zn incorporated varies across two very different types of growth bands, annual bands in Beluga cementum and daily incremental bands in beaver dentin. To the best of the authors’ knowledge, this is the first observation of spatial variation of lattice parameter within growth bands.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Using SIMS to decode noisy stratigraphic δ 13 C variations in Ediacaran carbonates

Carbonate carbon isotope (δ 13 C carb ) chemostratigraphy is a valuable tool in Precambrian stratigraphic correlation. The effectiveness of this tool rests on the assumption that δ 13 C carb data record global seawater signals. However, in some cases δ 13 C carb data exhibit rapid and noisy stratigraphic variations that may have been influenced by authigenic or diagenetic carbonate minerals. To assess the contribution of non-primary minerals towards bulk carbonate carbon isotope values, we acquired SIMS (secondary ion mass spectrometry) δ 13 C carb data, electron microscopic data, and CL (cathodoluminescence) microscopic data from two Ediacaran successions—the lower Member II of the Doushantuo Formation in South China and the Mooifontein Member of the Zaris Formation in southern Namibia. The Doushantuo samples came from a stratigraphic interval with noisy meter-scale δ 13 C carb variations of up to 10‰, whereas the Mooifontein Member is characterized by consistent δ13Ccarb values with limited meter-scale variations of less than ~4‰. Our data show that the meter-scale stratigraphic variations in δ 13 C carb are also mirrored in the SIMS data at μm-mm scales in both Doushantuo and Mooifontein samples. In the Doushantuo samples, SIMS δ 13 C carb values of authigenic calcite vary by up to 10‰ over μm-mm scales and can be either higher or lower than those of the co-existing dolomite matrix, which also appears to be affected by authigenic carbonate on the basis of petrographic observation. Bulk-sample δ 13 C carb values measured on powders microdrilled from the same SIMS specimens are within the range of SIMS δ 13 C carb values. Thus, we infer that bulk-sample δ 13 C carb values of Doushantuo samples represent mixtures of different carbonate components. The Doushantuo and Mooifontein SIMS δ 13 C carb data validate the traditional method of using chemostratigraphic consistency to evaluate authigenic/diagenetic alteration, but given that secondary phases may be 13 C-enriched relative to matrix carbonate, caution is warranted when taking maximum values to approximate primary chemostratigraphic trends.

58 GEOSCIENCES↗

DOE OSTI2023

Daily and Seasonal Variations of Shallow Seismic Velocities in Southern California From Joint Analysis of H/V Ratios and Autocorrelations of Seismic Waveforms

Establishing a baseline of ongoing secular velocity variations at the subsurface can improve the accuracy of detecting and interpreting short-term velocity changes, and advance the understanding of observed seismic motions and the behavior of subsurface materials. Toward these goals, we develop and apply a deconvolved autocorrelation (DA) method to estimate regional daily and seasonal changes of seismic velocities in southern California. The DA method combines advantages of traditional autocorrelation and Horizontal-to-Vertical Spectral Ratio, and is used to analyze over 10 years of data recorded by 50 stations. The results indicate widespread daily and seasonal changes of up to 10% and 4%, respectively, in the top tens of meters of the crust. The thickness of the surface layer, distance from the coast, and topographic variations are important factors controlling the amplitudes of the resolved velocity variations. The results suggest that changes of soil moisture and thermoelastic strain are likely dominant factors affecting the daily and seasonal variations, respectively. The developed DA method can improve the accuracy and robustness of estimated changes of subsurface materials at other locations.

58 GEOSCIENCES↗

DOE OSTI2024

Sub-Diurnal Methane Variations on Mars Driven by Barometric Pumping and Planetary Boundary Layer Evolution

In recent years, the Tunable Laser Spectrometer within the Sample Analysis at Mars (TLS-SAM) instrument on board the Mars Science Laboratory (MSL) Curiosity rover has detected methane variations in the atmosphere at Gale crater. Methane concentrations appear to fluctuate seasonally as well as sub-diurnally, which is difficult to reconcile with an as-yet-unknown transport mechanism delivering the gas from underground to the atmosphere. To potentially explain the fluctuations, we consider barometrically induced transport of methane from an underground source to the surface, modulated by temperature-dependent adsorption. The subsurface fractured-rock seepage model is coupled to a simplified 1-D atmospheric mixing model to provide insights on the pattern of atmospheric methane concentrations in response to transient surface methane emissions, as well as to predict sub-diurnal variation in methane abundance for the northern summer period, which is a candidate time frame for a MSL Curiosity sampling campaign. Our analysis suggests that there is a lower limit to the subsurface fracture density that can produce the observed methane patterns, below which the atmospheric methane variations would be out of phase with the observations. The best-performing model scenarios indicate a significant, short-lived methane pulse just prior to sunrise, the detection of which by TLS-SAM would be a potential indicator of the contribution of barometric pumping to Mars' atmospheric methane variations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Genetic variation at transcription factor binding sites largely explains phenotypic heritability in maize

Abstract Comprehensive maps of functional variation at transcription factor (TF) binding sites (cis-elements) are crucial for elucidating how genotype shapes phenotype. Here, we report the construction of a pan-cistrome of the maize leaf under well-watered and drought conditions. We quantified haplotype-specific TF footprints across a pan-genome of 25 maize hybrids and mapped over 200,000 variants, genetic, epigenetic, or both (termed binding quantitative trait loci (bQTL)), linked tocis-element occupancy. Three lines of evidence support the functional significance of bQTL: (1) coincidence with causative loci that regulate traits, includingvgt1,ZmTRE1and the MITE transposon nearZmNAC111under drought; (2) bQTL allelic bias is shared between inbred parents and matches chromatin immunoprecipitation sequencing results; and (3) partitioning genetic variation across genomic regions demonstrates that bQTL capture the majority of heritable trait variation across ~72% of 143 phenotypes. Our study provides an auspicious approach to make functionalcis-variation accessible at scale for genetic studies and targeted engineering of complex traits.

Genetics & Heredity↗

DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗

DOE OSTI2020

Convergence acceleration of Monte Carlo many-body perturbation methods by using many control variates

Herein, the use of many control variates is proposed as a method to accelerate the second- and third-order Monte Carlo (MC) many-body perturbation (MC-MP2 and MC-MP3) calculations. A control variate is an exactly integrable function that is strongly correlated or anti-correlated with the target function to be integrated by the MC method. Evaluating both integrals as well as their covariances in the same MC run, one can effect a mutual cancellation of the statistical uncertainties and biases in the MC integrations, thereby accelerating its convergence considerably. Six and thirty-six control variates, whose integrals are known a priori, are generated for MC-MP2 and MC-MP3, respectively, by systematically replacing one or more two-electron-integral vertexes of certain configurations by zero-valued overlap-integral vertexes in their Goldstone diagrams. The variances and co- variances of these control variates are computed at a marginal cost, enhancing the overall efficiency of the MC-MP2 and MC-MP3 calculations by a factor of up to 14 and 20, respectively.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Efficient Berry phase calculation via adaptive variational quantum computing approach

We present an adaptive variational quantum algorithm to estimate the Berry phase accumulated by a nondegenerate ground state under cyclic, adiabatic evolution of a time-dependent Hamiltonian. Our method leverages cyclic adiabatic evolution of the Hamiltonian and employs adaptive variational quantum algorithms for state preparation and evolution, optimizing circuit efficiency while maintaining high accuracy. We benchmark our approach on dimerized Fermi–Hubbard chains with four sites, demonstrating precise Berry phase simulations in both noninteracting and interacting regimes. Our results show that circuit depths reach up to 106 layers for noninteracting systems and increase to 279 layers for interacting systems due to added complexity. In addition, we demonstrate the robustness of our scheme across a wide range of parameters governing adiabatic evolution and variational algorithms. These findings highlight the potential of adaptive variational quantum algorithms for advancing quantum simulations of topological materials and computing geometric phases in strongly correlated systems.

Mootz, Martin [Ames Laboratory (AMES), Ames, IA (U↗

DOE OSTI2025

Environmental drivers of spatial variation in tropical forest canopy height: Insights from NASA’s GEDI spaceborne LiDAR

Forest canopy height is a fundamental ecosystem property—influencing patterns of forest carbon storage and forest ecosystem responses to climate variability and change. Previous studies have analyzed environmental drivers influencing spatial variation in canopy height at landscape-to-regional scales; however, far less is known about the environmental determinants underlying regional and global scale variation in forest canopy height. Using the canopy height metrics products from Global Ecosystem Dynamics Investigation (GEDI), a space-borne Light Detection and Ranging (LiDAR) instrument specifically designed to characterize forest structure, we analyze the environmental correlates of spatial variation of global tropical forest canopy height. Our study demonstrates that climate, topography, and soil properties account for 75% of the variation in tropical forest canopy height. Elevation, dry season length, and solar radiation are the most important drivers in determining canopy height both locally and regionally. These results emphasize the vulnerability of tropical forest structure to ongoing changes in the earth’s climate and provide a valuable empirical baseline for tropical forest management.

Science & Technology - Other Topics↗

DOE OSTI2021

Variations in pharmacokinetic-pharmacodynamic target values across MICs and their potential impact on determination of susceptibility test interpretive criteria

Abstract Background An antibacterial drug’s susceptibility test interpretive criteria (STIC) are determined by integrating clinical, microbiological and pharmacokinetic-pharmacodynamic (PK-PD) data. PTA analysis plays a pivotal or supportive role in STIC determination and is heavily dependent on the PK-PD target values determined from animal PK-PD studies. Therefore, variations in PK-PD target values may impact STIC determination. Factors contributing to variation in the PK-PD target values include the number of and MICs for bacterial isolates used in animal PK-PD studies. Objectives To analyse the relationship between PK-PD target values and MICs, describe the variations in PK-PD target values of isolates and evaluate whether the proposed/target STICs were within the ranges of the MICs for isolates used in animal PK-PD studies. Methods A database was compiled for this research by screening animal PK-PD study reports submitted to the FDA from 10 new drug applications (NDAs). Results A relationship evaluation between PK-PD target values and MICs for tested isolates for seven drugs (that used AUC/MIC ratio as the PK-PD index) showed that, generally, the AUC/MIC values decreased with an increase in MIC. These target values were highly variable, with the percentage coefficient of variation ranging between 1% and 132% for isolates having the same MIC. For 16/27 (59%) drug/bacteria combinations from all 10 drugs, the proposed/target STICs were higher than the highest MIC for bacteria isolates evaluated, while 6/27 (22.5%) were lower. Conclusions This research suggests that careful considerations related to selection of bacterial isolates for animal PK-PD studies could strengthen the STIC determination process.

Waack, Ursula↗

DOE OSTI2020

Impacts of Soil Microbiome Variations on Root Colonization by Fungi and Bacteria and on the Metabolome of Populus tremula × alba

Trees depend on beneficial interactions between roots and soil microbes for their nutrition and protection against stresses. The soil microbiome provides the main reservoir of microbes for root colonization and is subject to natural variations that can affect its composition. It is not clear whether the tree’s root system is able to buffer the natural variations occurring in the soil microbiome to capture a stable and effective microbiome or whether these variations affect its microbiome to impact its physiology. To address this question, we planted cuttings of Gray Poplar (Populus tremula × alba clone 717-1B4) in natural soil taken from a poplar stand under the same tree over two consecutive years and grew them in a greenhouse. We analyzed the soil and root microbiomes by high throughput Illumina MiSeq sequencing of fungal rDNA internal transcribed spacer and bacterial 16S rRNA amplicons and we characterized the root metabolome by gas chromatography-mass spectrometry. Soil and root microbial communities significantly shifted over the 2 years. A modification of the balance between endophytes, saprophytes, and mycorrhizal fungi occurred in the roots as well as a replacement of some dominant operational taxonomic units by others. These modifications were correlated with a significant alteration of the levels of about 10% of primary and secondary metabolites, suggesting that natural fluctuations in soil microbial communities can have a profound impact on tree root metabolism and physiology. Tree roots functioning may thus be indirectly strongly affected by the effects of future extreme climatic variations on the soil microbiome.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Application of the variational autoencoder to detect the critical points of the anisotropic Ising model

We generalize the previous study on the application of variational autoencoders to the two-dimensional Ising model to a system with anisotropy. Due to the self-duality property of the system, the critical points can be located exactly for the entire range of anisotropic coupling. This presents an excellent test bed for the validity of using a variational autoencoder to characterize an anisotropic classical model. Furthermore, we reproduce the phase diagram for a wide range of anisotropic couplings and temperatures via a variational autoencoder without the explicit construction of an order parameter. Considering that the partition function of ($d$ + 1)-dimensional anisotropic models can be mapped to that of the $d$-dimensional quantum spin models, the present study provides numerical evidence that a variational autoencoder can be applied to analyze quantum systems via the quantum Monte Carlo method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Genome‐wide association studies reveal genetic basis of ionomic variation in cassava

SUMMARY As one of the most important food crops, cassava ( Manihot esculenta ) is the main dietary source of micronutrients for about 1 billion people. However, the ionomic variation in cassava and the underlying genetic mechanisms remain unclear so far. Herein, genome‐wide association studies were performed to reveal the specific single nucleotide polymorphisms (SNPs) that affect the ionomic variation in cassava. We identified 164 SNPs with P‐ values lower than the threshold located in 88 loci associated with divergent ionomic variations. Among them, 13 SNPs are related to both calcium (Ca) and magnesium (Mg), and many loci for different ionomic traits seem to be clustered on specific chromosome regions. Moreover, we identified the peak SNPs in the promoter regions of Sc10g003170 (encoding methionyl‐tRNA synthetase [MetRS]) and Sc18g015190 (encoding the transcriptional regulatory protein AlgP) for nitrogen (N) and phosphorus (P) accumulation, respectively. Notably, these two SNPs (chr10_32807962 and chr18_31343738) were directly correlated with the transcript levels of Sc10g003170 ( MetRS ) and Sc18g015190 ( AlgP ), which positively modulated N accumulation and P concentration in cassava, respectively. Taken together, this study provides important insight into the genetic basis of cassava natural ionomic variation, which will promote genetic breeding to improve nutrient use and accumulation of elements in cassava.

Yin, Hongyan↗

DOE OSTI2025

Cold-air outbreaks in the continental US: Connections with stratospheric variations

Mid-latitude Northern Hemisphere extreme cold events continue to occur despite overall winter warming trends. These events have been linked to weakened stratospheric polar vortex (SPV) states. In this study, we analyze both the upper and lower polar stratosphere for links to extreme winter cold and snow in the continental US, finding two SPV variations of interest. The first features an upper-level vortex displaced toward western Canada and linked to northwestern US severe winter weather. The second features a weakened upper-level vortex displaced toward the North Atlantic and linked to central-eastern US severe winter weather. Both variations feature lower-level stretched vortices and stratospheric wave reflection. Since 2015, a northwestward shift in severe winter weather across the US is concurrent with an increase in the frequency of the westward-focused variation relative to the eastward-focused variation and a shift to more negative phases of the El Niño–Southern Oscillation.

Science & Technology - Other Topics↗

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↗