Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradients”

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 829 records · Page 46

Chemical cartography with LAMOST and Gaia reveal azimuthal and spiral structure in the Galactic disc

Chemical Cartography, or mapping, of our Galaxy has the potential to fully transform our view of its structure and formation. In this work, we use chemical cartography to explore the metallicity distribution of OBAF-type disc stars from the LAMOST survey and a complementary sample of disc giant stars from Gaia DR3. We use these samples to constrain the radial and vertical metallicity gradients across the Galactic disc. We also explore whether there are detectable azimuthal variations in the metallicity distribution on top of the radial gradient. For the OBAF-type star sample from LAMOST, we find a radial metallicity gradient of Δ[Fe/H]/ΔR ~–0.078 ± 0.001 dex kpc –1 in the plane of the disc and a vertical metallicity gradient of Δ[Fe/H]/ΔZ ~–0.15 ± 0.01 dex kpc –1 in the solar neighbourhood. The radial gradient becomes shallower with increasing vertical height, while the vertical gradient becomes shallower with increasing Galactocentric radius, consistent with other studies. We also find detectable spatially dependent azimuthal variations on top of the radial metallicity gradient at the level of ~0.10 dex. Interestingly, the azimuthal variations appear be close to the Galactic spiral arms in one data set (Gaia DR3) but not the other (LAMOST). Furthermore, these results suggest that there is azimuthal structure in the Galactic metallicity distribution and that in some cases it is co-located with spiral arms.

79 ASTRONOMY AND ASTROPHYSICS↗

Molecular dynamics of the early stages of high-temperature corrosion

In this work, we performed molecular dynamics simulations of the high-temperature oxidation of metal alloys composed of Al, Cr, and Fe and compared their behavior with that of pure Fe. The metal alloy elements (Al and Cr) segregated to the surface during oxidation, producing a lower stress gradient at the metal/environment interface compared to pure Fe. We have found that the lowered stress gradients produced in the alloy material appear to play a key role in the development of corrosion. Interfaces with lower stress gradients have reduced rates of H 2 O adsorption, especially for the ferritic (bcc) alloys. The diffusivity of oxygen and hydrogen drops more rapidly for the interfaces with reduced stress gradients. The stress gradient is also diminished when the gas pressure is increased, indicating that the Fe-Cr-Al alloy system is more resistant to oxidation than pure Fe at higher pressures. Therefore, we conclude that the lower stress gradients at the alloy/environment interface reduce the stress concentration and can slow down the rate of the initial oxide scale growth. We also compared bcc and fcc alloys with pure Fe based on our three evaluation criteria (peak stress, stress gradient, and summation of stress in the oxide scale). We found that the alloys have lower values under the three criteria compared to pure Fe. The bcc alloy has the best score under a water rich environment and the fcc alloy is proven to be better for peak stress and summation of stress in the oxide scale under an oxygen rich environment. For surface segregation to occur, we find that a minimum content of Al or Cr in the near-surface region must be achieved. We also learned that the role of Al is more important than that of Cr in terms of corrosion resistant behavior at relatively higher temperatures for the Fe-Cr-Al ternary alloys.

36 MATERIALS SCIENCE↗

The Open Cluster Chemical Abundances and Mapping Survey. IV. Abundances for 128 Open Clusters Using SDSS/APOGEE DR16

The Open Cluster Chemical Abundances and Mapping (OCCAM) survey aims to constrain key Galactic dynamical and chemical evolution parameters by the construction of a large, comprehensive, uniform, infrared-based spectroscopic data set of hundreds of open clusters. This fourth contribution from the OCCAM survey presents analysis using Sloan Digital Sky Survey/APOGEE DR16 of a sample of 128 open clusters, 71 of which we designate to be “high quality” based on the appearance of their color–magnitude diagram. We find the APOGEE DR16 derived [Fe/H] abundances to be in good agreement with previous high-resolution spectroscopic open cluster abundance studies. Using the high-quality sample, we measure Galactic abundance gradients in 16 elements, and find evolution of some of the [X/Fe] gradients as a function of age. We find an overall Galactic [Fe/H] versus R {sub GC} gradient of −0.068 ± 0.001 dex kpc{sup −1} over the range of 6 < R {sub GC} < 13.9 kpc; however, we note that this result is sensitive to the distance catalog used, varying as much as 15%. We formally derive the location of a break in the [Fe/H] abundance gradient as a free parameter in the gradient fit for the first time. We also measure significant Galactic gradients in O, Mg, S, Ca, Mn, Cr, Cu, Na, Al, and K, some of which are measured for the first time. Our large sample allows us to examine four well-populated age bins in order to explore the time evolution of gradients for a large number of elements and comment on possible implications for Galactic chemical evolution and radial migration.

79 ASTRONOMY AND ASTROPHYSICS↗

RCEMIP-II: mock-Walker simulations as phase II of the radiative–convective equilibrium model intercomparison project

Abstract. The radiative–convective equilibrium (RCE) model intercomparison project (RCEMIP) leveraged the simplicity of RCE to focus attention on moist convective processes and their interactions with radiation and circulation across a wide range of model types including cloud-resolving models (CRMs), general circulation models (GCMs), single-column models, global cloud-resolving models, and large-eddy simulations. While several robust results emerged across the spectrum of models that participated in the first phase of RCEMIP (RCEMIP-I), two points that stand out are (1) the strikingly large diversity in simulated climate states and (2) the strong imprint of convective self-aggregation on the climate state. However, the lack of consensus in the structure of self-aggregation and its response to warming is a barrier to understanding. Gaining a deeper understanding of convective aggregation and tropical climate will require reducing the degrees of freedom with which convection can vary. Therefore, we propose phase II of RCEMIP (RCEMIP-II) that utilizes a prescribed sinusoidal sea surface temperature (SST) pattern to provide a constraint on the structure of convection and move one critical step up the model hierarchy. This so-called “mock-Walker” configuration generates features that resemble observed tropical circulations. The specification of the mock-Walker protocol for RCEMIP-II is described, along with example results from one CRM and one GCM. RCEMIP-II will consist of five required simulations: three simulations with the same three mean SSTs as in RCEMIP-I but with an SST gradient and two additional simulations at one of the mean SSTs with different values of the SST gradients. We also test the sensitivity to the imposed SST gradient and the domain size. Under weak SST gradients, unforced self-aggregation emerges across the entire domain, similar to what was found in RCEMIP. As the SST gradient increases, the convective region narrows and is more confined to the warmest SSTs. At warmer mean SSTs and stronger SST gradients, low-frequency variability in the convective aggregation emerges, suggesting that simulations of at least 200 d may be needed to achieve robust equilibrium statistics in this configuration. Simulations with different domain sizes generally have similar mean statistics and convective structures, depending on the value of the SST gradient. The prescribed SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I, which enables comparison between the two; however, we also welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.

Geology↗

Engineering Redox Flow Battery Electrodes with Spatially Varying Porosity Using Non‐Solvent‐Induced Phase Separation

Redox flow batteries (RFBs) are a promising electrochemical platform for efficiently and reliably delivering electricity to the grid. Within the RFB, porous carbonaceous electrodes facilitate electrochemical reactions and distribute the flowing electrolyte. Tailoring electrode microstructure and surface area can improve RFB performance, lowering costs. Electrodes with spatially varying porosity may increase electrode utilization and provide surface area in reaction‐limited zones; however, the efficacy of such designs remains an open area of research. Herein, a non‐solvent‐induced phase‐separation (NIPS) technique that enables the reproducible synthesis of macrovoid‐free electrodes with well‐defined across‐thickness porosity gradients is described. The monotonically varying porosity profile is quantified and the physical properties and surface chemistries of porosity‐gradient electrodes are compared with macrovoid‐containing electrode, also synthesized by NIPS. Then, the electrochemical and fluid dynamic performance of the porosity‐gradient electrodes is evaluated, exploring the effect of changing the direction of the porosity gradient and benchmarking against the macrovoid‐containing electrode. Lastly, the performance is examined in a vanadium RFB, finding that the porosity‐gradient electrode outperforms the macrovoid electrode, is independent of gradient direction, and performs favorably compared to advanced electrodes in the contemporary literature. It is anticipated that the approach motivates further exploration of microstructurally tailored electrodes in electrochemical systems.

25 ENERGY STORAGE↗

Dynamic Excitonic Beam Switching with Atomically‐Thin Binary Blazed Gratings

Abstract Beam steering metasurfaces are ultra‐compact optical coatings that offer on‐demand redirection of optical power to specific diffraction orders. To achieve this, spatial gradients are commonly introduced in the phase of light scattered by plasmon or Mie resonant nanoparticles within the metasurface grating's unit cell. However, these phase gradients are oftentimes difficult to tune post‐fabrication. Recently, excitons in monolayer 2D semiconductors have emerged as a new metasurface building block, due to their strong and electrically‐tunable resonant light‐matter interaction. These 2D excitonic metasurfaces offer the tantalizing prospect of beam switching within a single monolayer. Here, it is demonstrated how the 2D analog of binary blazed gratings enables such beam switching by mere nanopatterning of a large monolayer WS 2 , even though nanoscale ribbons of WS 2 do not support geometrical resonances. By introducing a gradient in the nanoribbon width within the metasurface unit cell, an amplitude gradient combined with a small phase gradient in the scattered fields results in asymmetric diffraction efficiencies. Using a scattered‐field analysis, it is shown that these gradients can be further engineered via interference effects with the substrate reflection. Finally, the electrical tunability of the exciton resonance is leveraged to achieve selective and dynamic beam switching with an atomically‐thin metasurface.

Guarneri, Ludovica [Van der Waals‐Zeeman Institute↗

Degradation Reduces Microbial Richness and Alters Microbial Functions in an Australian Peatland

Abstract Peatland ecosystems cover only 3% of the world’s land area; however, they store one-third of the global soil carbon (C). Microbial communities are the main drivers of C decomposition in peatlands, yet we have limited knowledge of their structure and function. While the microbial communities in the Northern Hemisphere peatlands are well documented, we have limited understanding of microbial community composition and function in the Southern Hemisphere peatlands, especially in Australia. We investigated the vertical stratification of prokaryote and fungal communities from Wellington Plains peatland in the Australian Alps. Within the peatland complex, bog peat was sampled from the intact peatland and dried peat from the degraded peatland along a vertical soil depth gradient (i.e., acrotelm, mesotelm, and catotelm). We analyzed the prokaryote and fungal community structure, predicted functional profiles of prokaryotes using PICRUSt, and assigned soil fungal guilds using FUNGuild. We found that the structure and function of prokaryotes were vertically stratified in the intact bog. Soil carbon, manganese, nitrogen, lead, and sodium content best explained the prokaryote composition. Prokaryote richness was significantly higher in the intact bog acrotelm compared to degraded bog acrotelm. Fungal composition remained similar across the soil depth gradient; however, there was a considerable increase in saprotroph abundance and decrease in endophyte abundance along the vertical soil depth gradient. The abundance of saprotrophs and plant pathogens was two-fold higher in the degraded bog acrotelm. Soil manganese and nitrogen content, electrical conductivity, and water table level (cm) best explained the fungal composition. Our results demonstrate that both fungal and prokaryote communities are shaped by soil abiotic factors and that peatland degradation reduces microbial richness and alters microbial functions. Thus, current and future changes to the environmental conditions in these peatlands may lead to altered microbial community structures and associated functions which may have implications for broader ecosystem function changes in peatlands.

Environmental Sciences & Ecology↗

Interfacial characteristics of P91 steel - Inconel 740H bimetallic structure fabricated using wire-arc additive manufacturing

Although bimetallic structures can offer unique solutions to engineering problems with varying functionality, joining two dissimilar alloys using additive manufacturing involves challenges such as cracking and segregation. In this work, bimetallic structures of Inconel 740H superalloy and P91 steel are processed using wire-arc additive manufacturing (WAAM) for the first time. P91 steel is directly deposited over 740H with a cooling time of 20 minutes in between the two depositions. During microstructure characterization, a large gradient zone with coarse grains spanning for ~2200 μm is observed. Moreover, the microhardness of the gradient zone is the least in comparison with the P91 and 740H deposits. No deleterious phase formation is found; however, intergranular cracks are observed along the gradient zone. In order to deduce the reason for cracking, non-equilibrium solidification simulations are performed for the mixed compositions of P91 and 740H in the gradient zone. Here, the reasons for crack formation in the gradient zone are the sudden change in the volumetric coefficient of thermal expansion due to the formation of MC carbides as well as the development of local strains from the thermal residual stresses. Therefore, altering the deposition sequence and introducing gradient layers with mixed composition between the P91 and 740H deposits will be beneficial in achieving defect-free builds. Through the integration of experiments and simulation, this work gains valuable insights into the processability of dissimilar alloy manufacturing using the WAAM technique.

20 FOSSIL-FUELED POWER PLANTS↗

Cryogenic operation of planar ultrananocrystalline diamond field emission source in SRF injector

Here, we report a demonstration of electron beam generation in a superconducting radio frequency (SRF) Teraelectronvolt Energy Superconducting Linear Accelerator 1.3 GHz gun equipped with a field emission cathode and operated at 2 K. The cathode is a submicrometer film of nitrogen-incorporated ultrananocrystalline diamond [(N)UNCD] deposited atop a Nb RRR300 cathode plug that is 4 mm in diameter. Measurements of the output current showed that it increased exponentially as a function of the gun gradient. Our results demonstrate a feasible path toward simplified fully cryogenic SRF injector technology. At a maximal gradient of 0.9 MV/m, a current density of 22 μA/cm2 was obtained [calculated as a 2.75 μA current over the total cathode surface area of π(2 mm)2]. One important finding is that the electron emitter made of (N)UNCD, a material long been known as a highly efficient field emission material, demonstrated a record low turn-on gradient of 0.6 MV/m. A hypothesis explaining this behavior is proposed.

43 PARTICLE ACCELERATORS↗

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of triggering and evolution of ELM by pellet injection in EAST under BOUT++ framework

A BOUT ++ three-field magnetohydrodynamic model is employed to study the triggering and evolution of edge localized mode (ELM) by Li pellets injected along the outer mid-plane in the EAST configuration. The linear simulation shows that compared with a large deposition on the pedestal top (scenario I), a smaller deposition within the steep-gradient pedestal region (scenario II) can stimulate much larger linear growth rates of all-n peeling-ballooning modes (PBMs). The nonlinear simulation shows that there exists a pellet size threshold for ELM triggering for two deposition locations; the threshold for scenario I predicted in the present study matches the EAST observation well. Comparison of the two scenarios reveals that a smaller deposition is sufficient to trigger an ELM in a much shorter time in scenario II, whose ELM size is comparable to that in scenario I. This conclusion confirms previous DIII-D and ASDEX-Upgrade observations, suggesting that the steep-gradient pedestal region is a favorable deposition location for ELM triggering with minimum pellet size. Simulation analyses also find that the positive radial gradient of the hump-like pressure profile in the outer mid-plane induced by the pellet deposition plays a different role in the two scenarios. In scenario I, the force resulting from the gradient hinders the outflow of core plasmas and in return, the perturbation is suppressed from spreading inwards after ELM crashes. In scenario II, with a sizable deposition, the gradient results in another competitive perturbation growth region during the linear phase, thus dispersing the free energy and reducing the efficiency of destabilizing PBMs by pellet injection. The suppressing effect of saturated zonal flow on other modes, the short ELM fast crash phase, and the restricting transport effect of the positive radial pressure gradient work together to constrain the pedestal energy loss, especially when the pellet deposition amount is high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Natural evolutionary strategies for variational quantum computation

Abstract Natural evolutionary strategies (NES) are a family of gradient-free black-box optimization algorithms. This study illustrates their use for the optimization of randomly initialized parameterized quantum circuits (PQCs) in the region of vanishing gradients. We show that using the NES gradient estimator the exponential decrease in variance can be alleviated. We implement two specific approaches, the exponential and separable NES, for parameter optimization of PQCs and compare them against standard gradient descent. We apply them to two different problems of ground state energy estimation using variational quantum eigensolver and state preparation with circuits of varying depth and length. We also introduce batch optimization for circuits with larger depth to extend the use of ES to a larger number of parameters. We achieve accuracy comparable to state-of-the-art optimization techniques in all the above cases with a lower number of circuit evaluations. Our empirical results indicate that one can use NES as a hybrid tool in tandem with other gradient-based methods for optimization of deep quantum circuits in regions with vanishing gradients.

Anand, Abhinav (ORCID:0000000280812310)↗

On the use of adjoints in the inversion of observed quasi-static deformation

An adjoint-based conjugate gradient algorithm provides an efficient means for imaging sources of deformation within the Earth, such as volume stresses associated with fluid flow in aquifers and reservoirs. For time intervals over which the overburden deforms elastically, one can calculate the gradient elements for a single model update using just two numerical simulations. The first is a forward run that is used to compute the residuals associated with the given iteration. The second simulation is to evaluate the application of the adjoint operator to the residuals. In this adjoint calculation, the residual displacements are applied as sources at the measurement locations, driving the deformation in the simulation. The volume stress on the source grid blocks, in response to the residual displacements, provide the gradient components. We apply this technique to satellite-based interferometric synthetic aperture radar (InSAR) line-of-sight displacements that were observed over an oil reservoir in California’s Central Valley. Furthermore, we find that the adjoint-based gradient estimates, requiring 18 CPU seconds, agree with conventional numerical calculations that take over 3700 CPU seconds to compute. Conjugate gradient algorithms utilizing the conventional approach and adjoint-based gradient computations give roughly the same reductions in misfit and similar final estimates of reservoir volume change.

58 GEOSCIENCES↗

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗

Using AI for Management of Field Emission in SRF Linacs

Field emission control, mitigation, and reduction is critical for reliable operation of high gradient superconducting radio-frequency (SRF) accelerators. With the SRF cavities at high gradients, the field emission of electrons from cavity walls can occur and will impact the operational gradient, radiological environment via activated components, and reliability of CEBAF’s two linacs. A new effort has started to minimize field emission in the CEBAF linacs by re-distributing cavity gradients. To measure radiation levels, newly designed neutron and gamma radiation dose rate monitors have been installed in both linacs. Artificial intelligence (AI) techniques will be used to identify cavities with high levels of field emission based on control system data such as radiation levels, cryogenic readbacks, and vacuum loads. The gradients on the most offending cavities will be reduced and compensated for by increasing the gradients on least offensive cavities. Training data will be collected during this year’s operational program and initial implementation of AI models will be deployed. Preliminary results and future plans are presented.

Carpenter, Adam↗

TRANSFER LEARNING FOR FIELD EMISSION MITIGATION IN CEBAF SRF CAVITIES

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab operates hundreds of super-conducting radio frequency (SRF) cavities in its two linear accelerators (linacs). Field emission (FE) is an ongoing operational challenge in higher gradient SRF cavities. FE generates high levels of neutron and gamma radiation leading to damaged accelerator hardware and a radiation hazard environment. During machine development periods, we performed gradient scans to record data capturing the relationship between cavity gradients and radiation levels measured throughout the linacs. However, the field emission environment at CEBAF varies considerably over time as the configuration of the radio frequency (RF) gradients changes and due to the changing behaviour of field emitters. An artificial intelligence/machine learning (AI/ML) approach with transfer learning could be a valuable tool to mitigate FE and lower the radiation levels. In this work, we mainly focus on leveraging the RF trip data gathered during CEBAF operations. We develop a transfer learning-based surrogate model for radiation detector readings given RF cavity gradients to track the CEBAF?s changing configuration and environment. Then, we could use the developed model as an optimization process for redistributing the RF gradients within a linac to minimize radiation levels.

Ahammed, K.↗

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↗

Gap Flows Along the Columbia River Observed During the WFIP2 Field Campaign

Data collected in the Columbia River Gorge and Basin Area during the second Wind Forecasting Improvement Project (WFIP2) was used to study gap flow events in the region. In particular, the relationship between 200 m zonal wind speed and the zonal mesoscale pressure gradient at Troutdale (West of the gorge), Wasco (east of the gorge), and Boardman (farther east of the gorge) sites was examined. An objective criterion was used to identify 169 (Troutdale), 161 (Wasco), and 113 (Boardman) gap flow days that had weak synoptic forcing. The 200 m zonal wind speed exhibited a linear relationship with the zonal pressure gradient at Wasco during both easterly and westerly gap flow days at daily timescales. Such a relationship was only observed during easterly gap flow days at Troutdale and was not observed at Boardman. The zonal wind speed exhibited a diurnal cycle at Wasco and Boardman, with greater changes in the nighttime winds than the daytime winds due to changes in pressure gradient. The root mean square differences (RMSDs) between the observed pressure gradient at Wasco and that simulated by the High Resolution Rapid Refresh (HRRR) model in 3 km (HRRR-3km) and 750 m (HRRR-750m) horizontal resolution were 3.14 hPa and 0.57 hPa, respectively. However, the RMSDs for the 200 m zonal wind speed were 1.40 m s -1 for HRRR-3km and 1.20 m s -1 for HRRR- 750m. Improvement in simulating the zonal pressure gradient accompanied by negligible changes in simulating zonal wind speed due to higher horizontal resolution points to a complex interplay between different components of the model physics when simulating the meteorology. The HRRR-750m resolution model accurately mimicked the observed covariability (phase diagram) between the zonal pressure gradient and the 200 m zonal wind speed during the nighttime, but not during the daytime, suggesting model errors might be related to the representation of boundary layer processes.

54 ENVIRONMENTAL SCIENCES↗