SEARCH · Engineering Papers

Results for “interior”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Flood Estimation under Snowmelt and Rain-on-Snow Processes in Alaska: A Military Installation Perspective

Accurate flood estimation in snow-dominated and high-latitude regions remains challenging due to complex interactions among rainfall, snowmelt, and rain-on-snow (ROS) processes, which are not captured in conventional precipitation-based intensity-duration-frequency (PREC-IDF) curves. This study evaluates the Next-Generation IDF (NG-IDF) framework, an extension of PREC-IDF that incorporates total water available for runoff (precipitation minus changes in snow water equivalent), in two contrasting Alaskan watersheds of Department of Defense (DoD) significance: the Little Chena River Basin (LCRB) in interior Alaska, a tributary of the Chena River that flows through Fort Wainwright and near Eielson Air Force Base, and the Upper Ship Creek Basin (USCB), which drains Joint Base Elmendorf-Richardson near Anchorage. Using long-term SNOTEL observations and event-based rainfall-runoff modeling, NG-IDF and PREC-IDF flood estimates were compared against observation-based flood frequency analyses. Results show that NG-IDF consistently reduces flood-estimation bias by 15–20% relative to PREC-IDF, particularly for snowmelt- and ROS-dominated events. In the interior LCRB, permafrost conditions can substantially amplify flood responses during snowmelt events, an effect not explicitly represented in standard design tools. These findings demonstrate that NG-IDF provides a more physically consistent and transferable framework for flood estimation in cold regions, with potential relevance to mission-critical DoD installation resilience. Projected increases in ROS frequency and permafrost degradation across Alaska further emphasize the need to integrate physics-based hydrologic models that explicitly represent snow and permafrost processes to enhance design resilience and operational readiness for military and other critical infrastructure.

Yan, Hongxiang (ORCID:000000022387403X)↗

DOE OSTI2023

Investigation of Nanoparticle Degradation in Hydrogen Fuel Cell Systems through Automated Electron Microscopy

Proton exchange membrane fuel cells (PEMFC) are promising devices for the deployment of hydrogen-powered heavy-duty vehicles, providing a higher efficiency for similar driving range and fueling time than the existing ones. However, PEMFCs still encounter durability challenges mainly due to catalyst degradation in the cathode. Mitigating these performance losses requires a better understanding of the degradation mechanisms under heavy-duty accelerated stress tests (ASTs) [1]. Scanning transmission electron microscopy (STEM) combined with energy dispersive X-ray spectroscopy (EDS) are key tools for the analysis of Pt and PtCo nanoparticle size, spatial distribution and composition [2]. Electron tomography is also used to determine the rate and type of degradation of catalyst nanoparticles as a function of their position on the carbon support. In this work, automated data acquisition software, paired with a custom Python code, have been used to study the effect of different accelerated stress tests (ASTs) on nanoparticle coarsening [2]. Figure 1 shows high-angle annular dark-field (HAADF)-STEM images and EDS maps comparing the cathodes of membrane electrode assemblies (MEAs) following an electrocatalyst AST performed under H2/N2 with that of the heavy-duty AST performed under H2/air. We will discuss how AST conditions affect considerably the spatial distribution of the nanoparticles across the electrode between the membrane and microporous layer. Although the median particle size increased more in the MEA aged under the heavy-duty AST, as determined using a high-throughput image analysis, the quantitative EDS measurements demonstrate that the electrocatalyst AST resulted in more Pt and Co dissolution from the cathode, which is another important indicator of electrocatalyst degradation. We will further present the impact of the relative humidity (% RH) on the degradation mechanisms demonstrated using the same approach. Electron tomography has been used to distinguish the Pt nanoparticles residing on the carbon support surface (exterior) from those within the pore structure (interior) in order to determine the relative stability of interior and exterior nanoparticles. As shown in Figure 2, we will compare the Pt catalyst particle size at the beginning of test (BOT) and end of test (EOT), and discuss the importance of automating the electron tomography workflow, i.e. acquisition, reconstruction, and visualization, to increase sampling and determine the standard deviation of these measurement. The outlook for utilizing low-dose cryo-tomography for limiting damage to the catalyst, support, and especially proton-conducting ionomer will also be discussed [3].

Amichi, Lynda↗

DOE OSTI2022

Investigating the architecture and internal structure of the TOI-561 system planets with CHEOPS, HARPS-N, and TESS

ABSTRACT We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of 62 HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b (P = 0.45 d, R = 1.42 R⊕, M = 2.0 M⊕), c (P = 10.78 d, R = 2.91 R⊕, M = 5.4 M⊕), d (P = 25.7 d, R = 2.82 R⊕, M = 13.2 M⊕), and e (P = 77 d, R = 2.55 R⊕, M = 12.6 R⊕). Moreover, we identify an additional, long-period signal (>450 d) in the RVs, which could be due to either an external planetary companion or to stellar magnetic activity. The precise masses and radii obtained for the four planets allowed us to conduct interior structure and atmospheric escape modelling. TOI-561 b is confirmed to be the lowest density (ρb = 3.8 ± 0.5 g cm−3) ultra-short period (USP) planet known to date, and the low metallicity of the host star makes it consistent with the general bulk density-stellar metallicity trend. According to our interior structure modelling, planet b has basically no gas envelope, and it could host a certain amount of water. In contrast, TOI-561 c, d, and e likely retained an H/He envelope, in addition to a possibly large water layer. The inferred planetary compositions suggest different atmospheric evolutionary paths, with planets b and c having experienced significant gas loss, and planets d and e showing an atmospheric content consistent with the original one. The uniqueness of the USP planet, the presence of the long-period planet TOI-561 e, and the complex architecture make this system an appealing target for follow-up studies.

Lacedelli, G. (ORCID:0000000241977374)↗

DOE OSTI2022

Discovery of TOI-1260d and the characterization of the multiplanet system

ABSTRACT We report the discovery of a third planet transiting the star TOI-1260, previously known to host two transiting sub-Neptune planets with orbital periods of 3.127 and 7.493 d, respectively. The nature of the third transiting planet with a 16.6-d orbit is supported by ground-based follow-up observations, including time-series photometry, high-angular resolution images, spectroscopy, and archival imagery. Precise photometric monitoring with CHEOPS allows to improve the constraints on the parameters of the system, improving our knowledge on their composition. The improved radii of TOI-1260b and TOI-1260c are $2.36 \pm 0.06 \rm R_{\oplus }$, $2.82 \pm 0.08 \rm R_{\oplus }$, respectively while the newly discovered third planet has a radius of $3.09 \pm 0.09 \rm R_{\oplus }$. The radius uncertainties are in the range of 3 per cent, allowing a precise interpretation of the interior structure of the three planets. Our planet interior composition model suggests that all three planets in the TOI-1260 system contains some fraction of gas. The innermost planet TOI-1260b has most likely lost all of its primordial hydrogen-dominated envelope. Planets c and d were also likely to have experienced significant loss of atmospheric through escape, but to a lesser extent compared to planet b.

Lam, K. W. F. (ORCID:0000000299106088)↗

DOE OSTI2023

TESS giants transiting giants V – two hot Jupiters orbiting red giant hosts

ABSTRACT In this work, we present the discovery and confirmation of two hot Jupiters orbiting red giant stars, TOI-4377 b and TOI-4551 b, observed by Transiting Exoplanet Survey Satellite in the Southern ecliptic hemisphere and later followed-up with radial-velocity (RV) observations. For TOI-4377 b, we report a mass of $0.957^{+0.089}_{-0.087} \ M_\mathrm{J}$ and a inflated radius of 1.348 ± 0.081 RJ orbiting an evolved intermediate-mass star (1.36 M⊙ and 3.52 R⊙; TIC 394918211) on a period of of 4.378 d. For TOI-4551 b, we report a mass of 1.49 ± 0.13 MJ and a radius that is not obviously inflated of $1.058^{+0.110}_{-0.062} \ R_\mathrm{J}$, also orbiting an evolved intermediate-mass star (1.31 M⊙ and 3.55 R⊙; TIC 204650483) on a period of 9.956 d. We place both planets in context of known systems with hot Jupiters orbiting evolved hosts, and note that both planets follow the observed trend of the known stellar incident flux-planetary radius relation observed for these short-period giants. Additionally, we produce planetary interior models to estimate the heating efficiency with which stellar incident flux is deposited in the planet’s interior, estimating values of $1.91 \pm 0.48~{{\ \rm per\ cent}}$ and $2.19 \pm 0.45~{{\ \rm per\ cent}}$ for TOI-4377 b and TOI-4551 b, respectively. These values are in line with the known population of hot Jupiters, including hot Jupiters orbiting main-sequence hosts, which suggests that the radii of our planets have re-inflated in step with their parent star’s brightening as they evolved into the post-main sequence. Finally, we evaluate the potential to observe orbital decay in both systems.

Pereira, Filipe (ORCID:0000000221577146)↗

DOE OSTI2021

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Surface stress tensor and junction conditions on a rotating null horizon

The general form of the surface stress tensor of an infinitesimally thin shell located on a rotating null horizon is derived, when different interior and exterior geometries are joined there. Even though the induced metric on the surface must be the same approached from either side, the first derivatives of the metric need not be. Such discontinuities lead to a Dirac δ-distribution in the Einstein tensor localized on the horizon. For a general stationary axisymmetric geometry the surface stress tensor can be expressed in terms of two geometric invariants that characterize the surface, namely the discontinuities [κ] and [J] of the surface gravity κ and angular momentum density J. The Komar energy and angular momentum are given in coordinates adapted to the Killing symmetries, and the surface contributions to each determined in terms of [κ] and [J] . Guided by these, a simple modification of the original Israel junction conditions is verified directly from the Einstein tensor density to give the correct finite result for the surface stress, when the normal n to the surface is allowed to tend continuously to a null vector. The relation to Israel’s original junction conditions, which fail on null surfaces, is given. The modified junction conditions are suitable to the matching of a rotating “black hole” exterior to any interior geometry joined at the Kerr null horizon surface, even when the surface normal is itself discontinuous and the Barrabès-Israel formalism is also inapplicable. This joining on a rotating null horizon is purely of the matter shell type and does not contain a propagating gravitational shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Constraining solar electron number density via neutrino flavor data at Borexino

Understanding the physics of the deep solar interior, and the more exotic environs of core-collapse supernovae (CCSN) and binary neutron-star (NS) mergers, is of keen interest in many avenues of research. To date, this physics is based largely on simulations via forward integration. While these simulations provide valuable constraints, it could be insightful to adopt the "inverse approach" as a point of comparison. Within this paradigm, parameters of the solar interior are not output based on an assumed model, but rather are inferred based on real data. Here, we take the specific case of solar electron number density, which historically is taken as output from the standard solar model. We show how one may arrive at an independent constraint on that density profile based on available neutrino flavor data from the Earth-based Borexino experiment. The inference technique's ability to offer a unique lens on physics can be extended to other datasets, and to analogous questions for CCSN and NS mergers, albeit with simulated data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Accurate Determination of the Neutron Skin Thickness of Pb 208 through Parity-Violation in Electron Scattering

We report a precision measurement of the parity-violating asymmetry A PV in the elastic scattering of longitudinally polarized electrons from 208 Pb. We measure A PV = 550 ± 16 (stat) ±8 (syst) parts per billion, leading to an extraction of the neutral weak form factor F W (Q 2 = 0.00616 GeV 2 ) = 0.368 ± 0.013. Combined with our previous measurement, the extracted neutron skin thickness is R n -R p = 0.283 ± 0.071 fm. The result also yields the first significant direct measurement of the interior weak density of 208 Pb: ρ$^0_W$ = -0.0796 ± 0.0036(exp) ± 0.0013(theo) fm -3 leading to the interior baryon density ρ$^0_b$ = 0.1480 ± 0.0036(exp) ± 0.0013(theo) fm -3 . Finally, the measurement accurately constrains the density dependence of the symmetry energy of nuclear matter near saturation density, with implications for the size and composition of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Carrier properties of Bi(111) grown on mica and Si(111)

A comparison is presented between properties of Bi(111) films grown by Stranski-Krastanov epitaxy on Si(111) 7 × 7 and by van der Waals epitaxy on mica substrates. Thin film morphologies and electronic transport properties of Bi(111) films of variable thickness are investigated for each growth method. Atomic force micrographs for Bi(111) films on mica reveal clearly defined triangular regions consisting of layered steps with height 0.4 nm, corresponding to the Bi(111) bilayer height. Variable-temperature electronic transport measurements show the existence of a quantum confinement-induced energy gap in the film interiors, resulting in a semimetal-to-semiconductor transition. Magnetotransport analysis in a three-carrier model including metallic electrons in surface states and electrons and holes in the films' interiors provides a detailed study of densities, mobilities, and mean-free paths of the three carrier types. Further, improved electronic transport properties are found in Bi(111) films of higher thickness on mica compared to on Si(111), a likely result of the largely strain-free van der Waals epitaxial growth.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗

DOE OSTI2023

Optimization and Comparison of Modern Offshore Wind Turbine Generators Using GeneratorSE 2.0

As the offshore wind industry keeps growing at a rapid pace, developers are bracing themselves for a huge demand in critical rare earth metals which will threaten an already vulnerable supply chain. The wind energy industry is addressing this problem by investing in modern generator technologies that employ magnets with reduced rare earth content and high-field magnets enabled by rare-earth-free superconductors. In this paper we introduce the National Renewable Energy Laboratory's newly advanced GeneratorSE 2.0, which is a design and optimization tool that was developed to investigate the feasibility of such modern generators. Two direct-drive generator topologies with different magnet materials and mounting arrangements are investigated: an outer-rotor, V-shaped interior permanent magnet generator, and an inner-rotor normally conducting armature, paired with a low-temperature superconducting field with race-track coils. These technologies were evaluated for a range of power ratings between 15 and 25 MW, which represent the next generation of offshore wind turbines for both fixed-bottom and floating applications. The analyses indicate a new trend favoring the low-temperature superconducting technology for the direct-drive system.

direct-drive generators↗

DOE OSTI2021

Planetary normal mode computation: Parallel algorithms, performance and reproducibility

This report is an extension of work entitled “Computing planetary interior normal modes with a highly parallel polynomial filtering eigensolver.” by Shi et al., [1] originally presented at the SC18 conference. A highly parallel polynomial filtered eigensolver was developed and exploited to calculate the planetary normal modes. The proposed method is ideally suited for computing interior eigenpairs for large-scale eigenvalue problems as it greatly enhances memory and computational efficiency. In this article, the second-order finite element method is used to further improve the accuracy as only the first-order finite element method was deployed in the previous work. The parallel algorithm, its parallel performance up to 20k processors, and the great computational accuracy are illustrated. The reproducibility of the previous work was successfully performed on the Student Cluster Competition at the SC19 conference by several participant teams using a completely different Mars-model dataset on different clusters. Both weak and strong scaling performances of the reproducibility by the participant teams were impressive and encouraging. The analysis and reflection of their results are demonstrated and future direction is discussed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Tidal–hydrological dynamics of water temperature across freshwater forested wetlands on the northeastern Pacific coast

Abstract Tidal freshwater forests were once extensive across temperate coastlines, but loss and fragmentation have made estimation of their ecosystem functions challenging. We measured water temperature for 2 years in three Sitka spruce tidal forests, a restoration site, and an adjacent emergent marsh on the Columbia River, Washington, United States. We assessed spatial variability of water temperature within sites including the effects of hydrology, differences among bay and tributary tidal forests, and differences between the tidal forests and the mainstem Columbia, the restoration site, and the emergent marsh. The tidal forests nearest to the bay had lower interior water temperatures than their channel confluences by up to 2.5°C (weekly median temperature) and 2.0°C (weekly maximum temperature), with most cooling occurring during the low‐flow months of July–September. Tributary sites had maximum temperatures up to 1.9°C cooler than bay sites and 4.2°C cooler than the mainstem. Temperatures in the two bay sites decreased by −0.16°C/100 m and −0.07°C/100 m, on average. The restoration site had the smallest within‐site temperature gradient. Differences in maximum temperatures were greatest when tidal range was low, while higher tidal ranges were associated with warmer and more variable site interiors relative to their confluences. These results suggest that water temperatures in these tidal forests can provide temperature refugia for cold water biota including salmon.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Topological Characterization and Uncertainty Visualization of Atmospheric Rivers

Atmospheric rivers (ARs) are long, narrow regions of water vapor in the Earth's atmosphere that transport heat and moisture from the tropics to the mid-latitudes. ARs are often associated with extreme weather events in North America and contribute significantly to water supply and flood risk. However, characterizing ARs has been a major challenge due to the lack of a universal definition and their structural variations. Existing AR detection tools (ARDTs) produce distinct AR boundaries for the same event, making the risk assessment of ARs a difficult task. Understanding these uncertainties is crucial to improving the predictability of AR impacts, including their landfall areas and associated precipitation, which could cause catastrophic flooding and landslides over the coastal regions. In this work, we develop an uncertainty visualization framework that captures boundary and interior uncertainties, i.e., structural variations, of an ensemble of ARs that arise from a set of ARDTs. We first provide a statistical overview of the AR boundaries using the contour boxplots of Whitaker et al. that highlight the structural variations of AR boundaries based on their nesting relationships. We then introduce the topological skeletons of ARs based on Morse complexes that characterize the interior variation of an ensemble of ARs. We propose an uncertainty visualization of these topological skeletons, inspired by MetroSets of Jacobson et al. that emphasizes the agreements and disagreements across the ensemble members. Through case studies and expert feedback, we demonstrate that the two approaches complement each other, and together they could facilitate an effective comparative analysis process and provide a more confident outlook on an AR's shape, area, and onshore impact.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Resolving spatial and energetic distributions of trap states in metal halide perovskite solar cells

We report the profiling of spatial and energetic distributions of trap states in metal halide perovskite single-crystalline and polycrystalline solar cells. The trap densities in single crystals varied by five orders of magnitude, with a lowest value of 2 × 10 11 per cubic centimeter and most of the deep traps located at crystal surfaces. The charge trap densities of all depths of the interfaces of the polycrystalline films were one to two orders of magnitude greater than that of the film interior, and the trap density at the film interior was still two to three orders of magnitude greater than that in high-quality single crystals. Suprisingly, after surface passivation, most deep traps were detected near the interface of perovskites and hole transport layers, where a large density of nanocrystals were embedded, limiting the efficiency of solar cells.

14 SOLAR ENERGY↗

DOE OSTI2021

Rapid retreat of the southwestern Laurentide Ice Sheet during the Bølling-Allerød interval

The timing of Laurentide Ice Sheet deglaciation along its southwestern margin controlled the evolution of large glacial lakes and has implications for human migration into the Americas. Accurate reconstruction of the ice sheet’s retreat also constrains glacial isostatic adjustment models and is important for understanding ice-sheet sensitivity to climate forcing. Despite its significance, retreat of the southwestern Laurentide Ice Sheet (SWLIS) is poorly constrained by minimum-limiting 14 C data. We present 26 new cosmogenic 10 Be exposure ages spanning the western Interior Plains, Canada. Using a Bayesian framework, we combine these data with geomorphic mapping, 10 Be, and high-quality minimum-limiting 14 C ages to provide an updated chronology. This dataset presents an internally consistent retreat record and indicates that the initial detachment of the SWLIS from its convergence with the Cordilleran Ice Sheet began by ca. 15.0 ka, concurrent with or slightly prior to the onset of the Bølling-Allerød interval (14.7–12.9 ka) and retreated >1200 km to its Younger Dryas (YD) position in ~2500 yr. Ice-sheet stabilization at the Cree Lake Moraine facilitated a meltwater drainage route to the Arctic from glacial Lake Agassiz within the YD, but not necessarily at the beginning. Our record of deglaciation and new YD constraints demonstrate deglaciation of the Interior Plains was ~60% faster than suggested by minimum 14 C constraints alone. In conclusion, numerical modeling of this rapid retreat estimates a loss of ~3.7 m of sea-level equivalent from the SWLIS during the Bølling-Allerød interval.

58 GEOSCIENCES↗

DOE OSTI2024

Compositional zoning of the Otowi Member of the Bandelier Tuff, Valles caldera, New Mexico, USA

The Otowi Member of the Bandelier Tuff erupted at ca. 1.60 Ma from the Valles caldera (New Mexico, USA). It consists of as much as 400 km 3 (dense rock equivalent) of strongly differentiated high-silica rhyolite and shows systematic upward variations in crystallinity, mineral chemistry, and trace element concentrations through its thickness, but the major element composition is almost constant and is near the low-pressure granite minimum. Incompatible trace elements in whole pumice fragments and glasses show well-correlated linear covariations. Upward zoning to lower abundances of incompatible trace elements is accompanied by development of overgrowths on quartz and alkali feldspar, although earlier-formed interiors of quartz and feldspar have near-constant compositions throughout the tuff, modified by cation diffusion in the case of feldspar. Melt inclusions in remnant quartz cores show diverse Pb isotope ratios, pointing to a wide range of distinct protoliths that contributed rhyolitic melt to the Otowi magma. Mineral thermometers suggest a modest temperature gradient through the melt body, perhaps of 40 °C, at the time of eruption. Chemical, textural, and mineralogical variations and volume-composition relations through the tuff are consistent with an origin for zoning by melting of a high-crystallinity cumulate layer beneath cognate supernatant liquid to produce denser, remobilized liquid of accumulative composition (i.e., the “modified mush model”). Melting may have occurred in several episodes. The latest of these episodes, probably thousands of years prior to eruption, introduced new rhyolitic liquid into the system and was associated with a thermal excursion, recorded in core compositions of pyroxene, during which much of the earlier crystal mass was dissolved. This left inherited cores and interiors of accumulated quartz and feldspar mantled with new growth having less-evolved compositions (higher Ti, Sr, and Ba). Changing solubility of zircon during cumulate melting produced a reversal of Zr concentrations. There is no clear petrologic evidence of a recharge eruption trigger; nonetheless, compositional zoning resulted mainly from repeated recharge-induced remobilization of quartz-feldspar cumulate. The Otowi system was built, evolved, and modified by several events over the course of a few hundred thousand years.

Wolff, J. A.↗