Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibration”

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 127 records · Page 7

Low and High LET Degradation Studies of Metal-Loaded Organic Phase Ligands in the ALSEP Process

Organic solutions comprising the Actinide Lanthanide Separation Process (ALSEP) solvent consisting of 0.5 M 2-ethylhexyl phosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) and 0.05 M N,N,N’,N’-tetra(2-ethylhexyl)diglycolamide (T2EHDGA) in n-dodecane were subjected to low LET and high LET irradiation before and after equilibration with an aqueous phase of 3 M HNO 3 . Degradation dose constants revealed greater ligand degradation due to gamma irradiation than alpha irradiation for both ligands. Furthermore, equilibration with nitric acid did not have a significant impact on ligand degradation for either irradiation source. Identified degradation products were similar for both gamma and alpha irradiation and occurred mostly through the rupture of the N–C carbonyl and C–O ether bonds for T2EHDGA and the C–Oether bond in HEH[EHP]. Acid contact appears to alter the degradation pathway by favoring the formation of higher molecular weight recombination products. Furthermore, mixed T2EHDGA-HEH[EHP]-NO 3 complexes were formed with Nd(III) after extraction from 3 M HNO 3 , and low LET gamma irradiation of the Nd(III) loaded organic solution produced similar degradation products as the organic solution absent of Nd(III). Interestingly, and likely due to the greater radiolytic susceptibility of T2EHDGA than HEH[EHP], a HEH[EHP]-Nd(III) complex appears to form as the T2EHDGA degrades with increasing absorbed dose.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Sensitivity of transverse momentum correlations to early-stage and thermal fluctuations

Transverse momentum correlations were recently measured by the ALICE collaboration at the LHC [Acharya et al (2020 Phys. Lett. B 804 135375)]. A long-range structure in terms of relative pseudorapidity of particle pairs is observed. This may imply some signal of the initial state owing to the sheer spread of the correlation. However, the fluctuations inside a thermally equilibrated medium have to be taken into account, serving as motivation for this paper. Using lattice Quantum Chromodynamics constraints, we predicted the development and spread of balancing correlations caused by energy-momentum conservation. At the same time, we propagated the Gaussian-shaped initial correlation using hydrodynamics to estimate its effects. Our findings suggest that the resulting correlation, known as ‘the ridge,’ is sensitive both to flucutations seeded in the pre-equilibrium stage and to those seeded in the equilibrated medium. This can provide important insight into the early stages of the collision.

ALICE↗

The structure of cluster merger shocks: turbulent width and the electron heating time-scale

We present a new 2MsChandra observation of the cluster merger Abell 2146, which hosts two huge M ~ 2 shock fronts each ~500 kpc across. For the first time, we resolve and measure the width of cluster merger shocks. The best-fitting width for the bow shock is 17 ± 1kpc and for the upstream shock is 10.7 ± 0.3 kpc. A narrow collisionless shock will appear broader in projection if its smooth shape is warped by local gas motions. We show that both shock widths are consistent with collisionless shocks blurred by local gas motions of 290 ± 30 km s –1 . The upstream shock forms later on in the merger than the bow shock and is therefore expected to be significantly narrower. From the electron temperature profile behind the bow shock, we measure the time-scale for the electrons and ions to come back into thermal equilibrium. We rule out rapid thermal equilibration of the electrons with the shock-heated ions at the 6σ level. The observed temperature profile instead favours collisional equilibration. For these cluster merger shocks, which have low sonic Mach numbers and propagate through a high β plasma, we find no evidence for electron heating over that produced by adiabatic compression. Our findings are expected to be valid for collisionless shocks with similar parameters in other environments and support the existing picture from the solar wind and supernova remnants. Furthermore, the upstream shock is consistent with this result but has a more complex structure, including a ~2ke V increase in temperature ~50 kpc ahead of the shock.

79 ASTRONOMY AND ASTROPHYSICS↗

Implications of BBN bounds for cosmic ray upscattered dark matter

We consider the big bang nucleosynthesis (BBN) bounds on light dark matter whose cross section off nucleons is sufficiently large to enable acceleration by scattering off of cosmic rays in the local galaxy. Such accelerated DM could then deposit energy in terrestrial detectors. Since this signal involves DM of mass ~keV-GeV and requires large cross sections ≳10-31 cm2 in a relativistic kinematic regime, we find that the DM population in this scenario is generically equilibrated with Standard Model particles in the early universe. For sufficiently low DM masses ≲10 MeV, corresponding to much of the favored region of many cosmic-ray upscattering studies, this equilibrated DM population adds an additional component to the relativistic energy density around T~few MeV and thereby spoils the successful predictions of BBN. In the remaining ~10 MeV–GeV mass range, the large couplings required in this scenario are either currently excluded or within reach of current or future accelerator-based searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bulk viscosity from Urca processes: npeμ matter in the neutrino-trapped regime

In this work, we extend our previous study of the bulk viscosity of hot and dense npe matter induced by the Urca process in the neutrino trapped regime to npeμ matter by adding the muonic Urca processes as well as the purely leptonic electroweak processes involving electron-muon transition. Here, the nuclear matter is modeled in a relativistic density functional approach with two different parametrizations which predict neutrino dominated matter (DDME2 model) and anti-neutrino dominated matter (NL3 model) at temperatures for which neutrinos/anti-neutrinos are trapped. In the case of neutrino-dominated matter, the main equilibration mechanism is lepton capture, whereas in the case of antineutrino-dominated matter this is due to neutron decay. We find that the equilibration rates of Urca processes are higher than that of the pure leptonic processes, which implies that the Urca bulk viscosity can be computed with the leptonic reactions assumed to be frozen. We find that the bulk viscosity decreases with temperature as ζ~T –2 at moderate temperatures. At high temperatures this scaling breaks down by sharp drops of the bulk viscosity close to the temperature where the proton fraction is density-independent and the matter becomes scale-invariant. This occurs also when the matter undergoes a transition from the antineutrino-dominated regime to the neutrino-dominated regime where the bulk viscosity attains a local maximum. We also estimate the bulk viscous dissipation timescales and find that these are in the range ≳ 1 s for temperatures above the neutrino trapping temperature. These timescales would be relevant only for long-lived objects formed in binary neutron star mergers and hot proto-neutron stars formed in core-collapse supernovas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impact of the equation of state on f - and p - mode oscillations of neutron stars

We investigate the impact of the neutron-star matter equation of state on the f- and p1-mode oscillations of neutron stars obtained within the Cowling approximation and linearized general relativity. The f- and p 1 -mode oscillation frequencies, and their damping times are calculated using representative sets of Skyrme Hartree-Fock and relativistic mean-field models, all of which reproduce nuclear systematics and support 2M ⊙ neutron stars. Our study shows strong correlations between the frequencies of f- and p 1 -modes and their damping times with the pressure of β-equilibrated matter at densities equal to or slightly higher than the nuclear saturation density ρ 0 . Such correlations are found to be almost independent of the composition of the stars. The frequency of the p 1 -mode of 1.4M ⊙ star is strongly correlated with the slope of the symmetry energy L 0 and β-equilibrated pressure at density ρ 0 . Compared to GR calculations, the error in the Cowling approximation for the f-mode is about 30% for neutron stars of low mass, whereas it decreases with increasing mass. Here, the accuracy of the p 1 -mode is better than 15% for neutron stars of maximum mass, and improves for lower masses and higher number of radial nodes. I.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermal decoupling of deuterium and tritium during the ICF shock-convergence phase

We report that a series of thin glass-shell shock-driven DT gas-filled capsule implosions was conducted at the OMEGA laser facility. These experiments generate conditions relevant to the central plasma during the shock-convergence phase of ablatively driven inertial confinement fusion (ICF) implosions. The spectral temperatures inferred from the DTn and DDn spectra are most consistent with a two-ion-temperature plasma, where the initial apparent temperature ratio, $\frac{T_T}{T_P}$, is 1.5. This is an experimental confirmation of the long-standing conjecture that plasma shocks couple energy directly proportional to the species mass in multi-ion plasmas. The apparent temperature ratio trend with equilibration time matches expected thermal equilibration described by hydrodynamic theory. This indicates that deuterium and tritium ions have different energy distributions for the time period surrounding shock convergence in ignition-relevant ICF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bound on hot-spot mix in high-velocity, high-adiabat direct-drive cryogenic implosions based on comparison of absolute x-ray and neutron yields

In laser-driven implosions for laboratory fusion, the comparison of hot-spot x-ray yield to neutron production can serve to infer hot-spot mix. For high-performance direct-drive implosions, this ratio depends sensitively on the degree of equilibration between the ion and electron fluids. A scaling for x-ray yield as a function of neutron yield and characteristic ion and electron hot-spot temperatures is developed on the basis of simulations with varying degrees of equilibration. We apply this model to hot-spot x-ray measurements of direct-drive cryogenic implosions typical of the direct-drive designs with best ignition metrics. Furthermore, the comparison of the measured x-ray and neutron yields indicates that hot-spot mix, if present, is below a sensitivity estimated as ~ 2% by-atom mix of ablator plastic into the hot spot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interface pinning causes the hysteresis of the hydride transformation in binary metal hydrides

Hydriding and dehydriding transitions in bulk and nanocrystalline binary metal hydrides were studied using the Pd-H model system by measuring pressure-composition isotherms with in situ x-ray diffractometry. Nanocrystalline Pd showed a smaller pressure hysteresis, solvus hysteresis, and hysteresis in lattice parameter, compared to bulk Pd. The time-dependence of pressure equilibration was measured after dosing with aliquots of hydrogen, giving equilibration times that were much faster in the single-phase regions than in the two-phase plateaus. In the broad two-phase plateaus, the pressure relaxations were exponential functions of time. An explanation of hysteresis is developed that is based on a dissipative potential barrier that impedes the motion of the interface due to interactions between lattice defects and the two-phase interface. The exponential pressure relaxations and hysteresis are consistent for this mechanism. For a simple model of the pinning potential, the potential barrier maximum is an order of magnitude less than typical grain boundary energies. These pinning effects are substantially different in the nanocrystalline Pd, suggesting differences in the hydriding mechanism.

36 MATERIALS SCIENCE↗

Dynamical control in a prethermalized molecular ultracold plasma: Local dissipation drives global relaxation

Prethermalization occurs as an important phase in the dynamics of isolated many-body systems when coupled degrees of freedom relax in a subspace separated from a state of complete thermodynamic equilibrium by a gap in energy or other conserved quantity. Slow equilibration from a prethermal phase can localize the dynamics of natural and model systems despite high dimensionality and limited disorder. Here, we report the signature of an enduring prethermal regime of arrested relaxation in the molecular ultracold plasma that forms following the avalanche of a state-selected Rydberg gas of nitric oxide. For a wide range of initial conditions, this system enters a critical phase in which a density of NO + and electrons balances a population of Rydberg molecules. Electron collisions mix orbital angular momentum, scattering Rydberg electrons to states of very high ℓ . The rapid predissociation of molecules in states of low- ℓ purifies this high- ℓ ensemble. The angular momentum barrier that separates NO + ions and electrons creates an extraordinary gap between the plasma states of n ≈ ℓ , with measured n > 200 and penetrating states of ℓ = 0 , 1 , and 2. Evolution to a statistically equilibrated state of N and O atoms cannot occur without Rydberg electron penetration, and this gap blocks relaxation for a millisecond or more. A weak radiofrequency (RF) field drives ℓ -mixing electron collisions that erase this gap, causing wholesale dissipation. Remarkably, a local quantum-state transition promoting an exceedingly small fraction of the molecules in the prethermalized ensemble to a predissociative state also acts with a global effect, driving the entire system to a dissipative equilibrium. Using the Lindblad master equation, we illustrate qualitatively similar dynamics for a toy model of an open quantum system that consists of a localized set of spins on which dissipation acts at a single site as a gateway to equilibrium. Published by the American Physical Society 2025

Wang, Ruoxi↗

Impact of neutrino flavor conversions on neutron star merger dynamics, ejecta, nucleosynthesis, and multimessenger signals

We present numerical relativity simulations of binary neutron star mergers incorporating neutrino flavor transformations triggered by fast flavor instability, quantum many-body effects, or potential beyond standard model physics. In both long-lived and short-lived remnant scenarios, neutrino flavor conversions modify species-dependent neutrino luminosities and mean energies, and drive the matter towards more neutron rich conditions. They produce up to more neutron rich ejecta and significantly boost the r-process yields, especially in low-density, near-equatorial outflows. We identify regions unstable to fast flavor instabilities and find that these instabilities persist despite flavor conversions. We further test the sensitivity to the equilibration timescale of the flavor conversions, finding that slower flavor conversions can interact with thermodynamic equilibration, and increase the neutron richness of the ejecta. Flavor conversions may also contribute to stronger gravitational wave and neutrino emissions, pointing to a correlation between neutrino transport and merger dynamics. Furthermore, these results highlight the potential impact of flavor conversions while motivating future work to improve on theoretical understanding of flavor instabilities in global simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical Scaling Reveals Topological Defects and Anomalous Evolution of a Photoinduced Phase Transition

Nonequilibrium states of quantum materials can exhibit exotic properties and enable unprecedented functionality and applications. These transient states are inherently inhomogeneous, characterized by the formation of topologically protected structures, requiring nanometer spatial resolution on femtosecond timescales to resolve their evolution. Using ultrafast total x-ray scattering at a free electron laser and a sophisticated scaling analysis, we gain unique access to the dynamics on the relevant mesoscopic length scales. Our results provide direct evidence that ultrafast excitation of LaTe 3 leads to formation of topological vortex strings of the charge density wave. These dislocations of the charge density wave exhibit anomalous, subdiffusive dynamics, slowing the equilibration process, providing rare insight into the nonequilibrium mesoscopic response in a quantum material. Our findings establish a general framework to investigate properties of topological defects, which are expected to be ubiquitous in nonequilibrium phase transitions and may arrest equilibration and enhance competing orders.

Orenstein, Gal [SLAC National Accelerator Laborato↗

Simulating evaporative wet and dry cycles in Gale crater, Mars using thermochemical modeling techniques

The aim of this work is to provide a model‐backed hypothesis for the formation of evaporites—sulfates, borates—in Gale crater using thermochemical modeling to determine constraints on their formation. We test the hypothesis that primary evaporites required multiple wet–dry cycles to form, akin to how evaporite assemblages form on Earth. Starting with a basalt‐equilibrated Mars fluid, Mars‐relevant concentrations of B and Li were added, and then equilibrated with Gale lacustrine bedrock. We simulated the cycles of evaporation followed by groundwater recharge/dilution to establish an approximate minimum number of wet–dry cycles required to form primary evaporites. We determine that a minimum of 250 wet–dry cycles may be required to start forming primary evaporites that consist of borates and Ca‐sulfates. We estimate that ~14,250 annual cycles (~25.6 k Earth years) of wet and dry periods may form primary borates and Ca‐sulfates in Gale crater. These primary evaporites could have been remobilized during secondary diagenesis to form the veins that the Curiosity rover observes in Gale crater. No Li salts form after 14,250 cycles modeled for the Gale‐relevant scenario (approximately 10 6 cycles would be needed) which implies Li may be leftover in a groundwater brine after the time of the lake. No major deposits of borates are observed to date in Gale crater which also implies that B may be leftover in the subsequent groundwater brine that formed after evaporites were remobilized into Ca‐sulfate veins.

58 GEOSCIENCES↗

Reversible structural transformations in supercooled liquid water from 135 to 245 K

Water has many anomalous properties compared to “simple” liquids, and these anomalies are typically enhanced in supercooled water.1-3 While numerous models have been proposed, including the liquid-liquid critical point,4,5 the singularity-free scenario,6 and the stability limit conjecture,1 a molecular-level understanding remains elusive. The main difficulty in determining which, if any, of these models is correct is the limited amount of data in the relevant temperature and pressure ranges. For water at ambient pressures, which is the focus of this work, data is largely missing from 160 – 232 K (“No Man’s Land”) due to rapid crystallization.2,3 Whether rapid crystallization is just an experimental obstacle, or a fundamental problem signaling the inability of water to thermally equilibrate prior to crystallization is also a major unanswered question.5,7 Here, we investigate the structural transformations of transiently-heated, supercooled water with nanosecond time resolution using infrared vibrational spectroscopy. The experiments demonstrate three key results. First, water’s structure relaxes from its initial configuration to a “steady-state” configuration prior to the onset of crystallization over a wide temperature range. Second, water’s steady-state structure can be reproduced by a linear combination of two, temperature-independent structures that correspond to a “high-temperature liquid” and a “low-temperature liquid.” Third, the observed structural changes are reversible over the full temperature range. Taken together, these results show that supercooled water can equilibrate prior to crystallization for temperatures from the homogeneous nucleation temperature, TH ~232 K,3 down to the glass transition temperature (Tg ~ 136 K). Second, the results provide support for the hypothesis that supercooled water can be described as a mixture of two, structurally-distinct, interconvertible liquids from 135 K to 245 K.5,8-18

Kringle, Loni M.↗

Electrically driven proton transfer promotes Brønsted acid catalysis by orders of magnitude

Electric fields play a key role in enzymatic catalysis and can enhance reaction rates by 100,000-fold, but the same rate enhancements have yet to be achieved in thermochemical heterogeneous catalysis. In this work, we probe the influence of catalyst potential and interfacial electric fields on heterogeneous Brønsted acid catalysis. We observed that variations in applied potential of ~380 mV led to a 100,000-fold rate enhancement for 1-methylcyclopentanol dehydration, which was catalyzed by carbon-supported phosphotungstic acid. Mechanistic studies support a model in which the interfacial electrostatic potential drop drives quasi-equilibrated proton transfer to the adsorbed substrate prior to rate-limiting C–O bond cleavage. Large increases in rate with potential were also observed for the same reaction catalyzed by Ti/TiO y H x and for the Friedel Crafts acylation of anisole with acetic anhydride by carbon-supported phosphotungstic acid.

Science & Technology - Other Topics↗

Estimating time in quantum chaotic systems and black holes

We characterize new universal features of the dynamics of chaotic quantum many-body systems, by considering a hypothetical task of "time estimation". Most macroscopic observables in a chaotic system equilibrate to nearly constant late-time values. Intuitively, it should become increasingly difficult to estimate the precise value of time by making measurements on the state. We use a quantity called the Fisher information from quantum metrology to quantify the minimum uncertainty in estimating time. Due to unitarity, the uncertainty in the time estimate does not grow with time if we have access to optimal measurements on the full system. Restricting the measurements to act on a small subsystem or to have low computational complexity leads to results expected from equilibration, where the time uncertainty becomes large at late times. With optimal measurements on a subsystem larger than half of the system, we regain the ability to estimate the time very precisely, even at late times. Hawking's calculation for the reduced density matrix of the black hole radiation in semiclassical gravity contradicts our general predictions for unitary quantum chaotic systems. Hawking's state always has a large uncertainty for attempts to estimate the time using the radiation, whereas our general results imply that the uncertainty should become small after the Page time. This gives a new version of the black hole information loss paradox in terms of the time estimation task. By restricting to simple measurements on the radiation, the time uncertainty becomes large. This indicates from a new perspective that the observations of computationally bounded agents are consistent with the semiclassical effective description of gravity.

Black holes↗

Decarburization and Carburization Behavior of Grade 91 Steel in Sodium Environments

Two heats of G91 steel have been investigated in sodium at 550-650°C to understand its corrosion behaviour, microstructural evolution, and tensile properties. Sodium exposure tests of subsized tensile and coupon specimens of G91 base metal and weldment were conducted in Argonne’s forced convection sodium materials testing loops. Maximum exposure times of ~39,000 h at 550°C, ~37,000 h at 600°C, and ~20,000 h at 650°C have been achieved, respectively. Thermal aging experiments of G91 steel were conducted in parallel, and the data of the sodium-exposed specimens were compared with the thermal aging data of the same heat to isolate the thermal aging effect from the effect of sodium exposure. It was found that the tensile strength of G91 steel can be reduced by more than 50% after sodium exposures at 650°C. The effect of a sodium environment on tensile properties is largely dependent on the carburization/decarburization behavior of G91 steel exposed to sodium. The carbon concentrations in sodium in Argonne’s SMT-1 and SMT-2 loops were determined by the equilibration method. The carbon activity in sodium was established by equilibrating high-purity nickel in sodium, analyzing the material for carbon, and using the reported data on the carbon solubility for this material. The carbon activity results were combined with the carbon solubility data in sodium to establish the carbon concentration in sodium. The carbon concentrations in the SMT-1 and SMT-2 loop were estimated to be 1.23 and 0.25 ppm, respectively. With this initial success, we will continue to monitor the carbon contents in the two loops to correlate the microstructural and mechanical property changes with the sodium purity. Future experiments will also consider sodium exposure tests in controlled carbon activities. Online monitoring and control of carbon contents in sodium in nuclear reactors would also be of interest. Thermodynamic analysis was conducted to understand the carburization-decarburization behavior of G91 steel exposed to sodium in the SMT-1 and SMT-2 loops. The equilibrium carbon activity in the alloy was evaluated by the thermodynamics of M 23 C 6 carbides. The calculated carburization-decarburization boundary showed that G91 would undergo decarburization at 650°C in the SMT-1 loop and carburization at 600°C in the SMT-2 loop, which is consistent with our findings from the sodium exposure tests. Future work will include thermodynamic analysis of the carbon concentration – activity relationship including the effect of carbide-forming elements, vanadium and niobium in G91 steel. Kinetics of carbon diffusion and the effect on microstructure and mechanical performance under various service conditions will also be investigated. The ultimate goal is to understand the mechanisms and the kinetics of carbon transfer in the sodium-steel system and to provide a basis for predicting the effect of long-term decarburization-carburization process on the integrity of reactor components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗