Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient projection”

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 307 records · Page 17

Examining transport and integrated modeling predictive capabilities for negative-triangularity scenarios

This paper investigates the predictive capabilities of TGYRO and TGLF models in assessing the performance of negative triangularity (NT) plasmas compared to positive triangularity (PT) plasmas in fusion devices. TGYRO predicts kinetic profiles, while TGLF analyzes turbulent transport. The study reveals that TGYRO reasonably predicts NT profiles similar to PT, although it overpredicts the high-power scenarios where there is increased experimental MHD activity. TGLF analysis finds reduced linear growth rates in NT and altered flux spectra relative to PT. Additionally, the TGLF SAT0 saturation model is observed to predict high-k transport and a reduction of particle transport with the electron temperature gradient. These findings are further corroborated by core-pedestal modeling using the Stability Transport Equilibrium Pedestal workflow, showing stronger confinement improvements in NT, particularly at higher power densities for the SAT0 saturation model. Furthermore, the study underscores the importance of accurately capturing turbulence saturation mechanisms for NT in order to project its performance accurately in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning and deep learning for mineralogy interpretation and CO 2 saturation estimation in geological carbon Storage: A case study in the Illinois Basin

Carbon capture and storage (CCS) is a promising approach to simultaneously maintaining energy security and reducing carbon dioxide (CO 2 ) emissions under the current energy portfolio that is dominated by fossil fuel energy. Pre-injection formation characterization and post-injection CO 2 monitoring are two critical tasks to guarantee storage efficiency in CCS. The CCS projects in the Illinois Basin, the first large-scale CO 2 injection into saline aquifers in the United States, employed conventional and the latest pulsed neutron logging (PNL) tools for mineralogy interpretation and CO 2 saturation estimation, which provide valuable references for future CCS projects. Because of the inherent fuzziness of petrophysical measurements and complex subsurface heterogeneity, interpreting well-logging data is time-consuming, and its accuracy can be user-biased. In recent years, data-driven methods have been widely used to capture the non-linear patterns between input features and interpretation results. This work applied and evaluated four commonly used machine learning (ML) models, including ridge regression (RR), random forest (RF), gradient boosting regression (GBR), support vector regression (SVR), and one deep learning (DL) model, the artificial neural network (ANN). We optimized the hyperparameters of the four ML models and the DL model using the simulated annealing algorithm and the grid search strategy, respectively. The input features of the mineralogy interpretation models were eleven conventional well-logging parameters, and the label data (i.e., ground truth) were the porosity and volumetric fractions of six minerals, including quartz, feldspar, dolomite, calcite, clay, and iron minerals. The results demonstrated that the GBR and RF models were superior in predicting volumetric fractions of minerals and porosity; label data with low coefficient of variation (CV) values tended to yield better performance. For CO 2 saturation estimation, the RF was the best-performing model, followed by SVR, ANN, GBR, and RR. Furthermore, we conducted feature importance ranking using the permutation importance algorithm and found that the formation sigma and well pressure were the most important features in this study. In conclusion, the study of CCS projects in the Illinois Basin bridges the gap between the limited knowledge and understanding of geological carbon storage and the increasing demand for reliable, cost-effective, and sustainable energy solutions.

58 GEOSCIENCES↗

Entanglement, trace anomaly, and confinement in QCD

We formulate confinement in quantum chromodynamics (QCD) as an entropic surface phenomenon. Quark and gluon quantum information is localized on a transverse, entangling two-sphere of radius 𝑅 𝐸⁢𝐸 ; at this radius the QCD vacuum—partitioned by a hadron into interior and exterior regions—reaches its maximal entanglement entropy. Lattice-QCD determinations of the scalar (trace) gravitational form factors fix both 𝑅 𝐸⁢𝐸 and the transverse trace-anomaly density 𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ), yielding a parameter-free slope 𝑐 ℎ = 8⁢𝜋 2 ⁢𝑅𝑆$^2_{𝐸⁢𝐸}$𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ) and a mechanical entropy 𝑆 𝐸⁢𝐸 ⁡(𝑦) = 𝑐 ℎ ⁢𝑦 that grows linearly with rapidity 𝑦. The entropy gradient ∂ 𝑅 𝑆 𝐸⁢𝐸 changes sign at 𝑅 𝐸⁢𝐸 : it pushes colored degrees of freedom outward for 𝑟 <𝑅 𝐸⁢𝐸 and pulls them inward for 𝑟 >𝑅 𝐸⁢𝐸 , thereby localizing them on the codimension-2 entangling two-sphere Σ ⊥ = 𝑆$^2_{𝑅_{𝐸⁢𝐸}}$ (which, in the infinite-momentum frame (IMF), projects onto the transverse plane)—the “information wall.” This provides a high-energy (large-𝑦) entropic confinement diagnostic that complements—rather than replaces—Wilson’s area-law criterion, which probes long-distance dynamics near the rest frame (𝑦 → 0). Imposing unitarity on an entropic ansatz for the amplitude yields 𝜎⁡(𝑠) ∝ 𝑦 𝛿 . World data favor 𝛿 = 2 for elastic 𝑝⁡𝑝⁡($𝑝\bar{⁡𝑝}$) scattering and heavy-quark photoproduction, whereas 𝜙 photoproduction favors a softer 𝛿 = 0.387. All extracted cross sections remain well below the Froissart-Martin bound. These results provide a confinement criterion quantified directly from nonperturbative QCD inputs, unifying the trace anomaly, entanglement entropy, and high-energy scattering within a single quantitative framework.

Color confinement↗

Maximizing the efficiency of plasma-based lepton accelerators

The next generation of particle accelerators will likely rely on compact accelerating structures to achieve very high gradients; one such candidate is accelerators using plasmas for the accelerating structures. Because of the small length scales of these structures, they can also be prone to the beam break-up instability, which would spoil the quality of the beam and leave it unsuitable for many applications. The goal of this project is to develop a framework for understanding the instabilities in plasma accelerators using language from conventional linear accelerator instabilities, such as wake fields and wake functions. During the course of this Phase I, we benchmarked a pseudo-spectral particle-in-cell code to determine its suitability for computing wake functions, and developed a procedure for extracting wake functions from a sequence of simulations using this code.

43 PARTICLE ACCELERATORS↗

Physical vapor deposition of organic glasses (Final Report)

As of 2009, when this project began, it had been established that physical vapor deposition (PVD) could produce high density and high stability (“ultrastable”) glasses of two organic molecules. With this project, we established that ultrastability was a quite general result for PVD glasses of organic molecules and, in particular, for the organic semiconductors used to produce organic light emitting diodes (OLEDs). We established that ultrastable glasses have enhanced photochemical and chemical stability, and increased resistance to crystallization. We developed a temperature-gradient deposition stage to more efficiently screen the properties of PVD glasses. We established that PVD glasses are usually anisotropic and that this anisotropy can be controlled. We characterized the unusual transformation mechanism of PVD glasses, which superficially is more similar to crystals than liquid-cooled glasses. In the most recent grant cycle, we extended these results to co-deposition of two organic semiconductors, showing that mixed PVD glasses of these systems usually showed ultrastability and high density; in two cases, co-deposition results in component segregation, and deposition conditions can be adjusted to control domain size. In our work since 2009, we established that the surface equilibration mechanism can successfully explain almost all of the properties of PVD glasses, and that this mechanism can provide useful predictions regarding stability and anisotropy.

36 MATERIALS SCIENCE↗

Montane forest productivity across a semiarid climatic gradient

Abstract High‐elevation montane forests are disproportionately important to carbon sequestration in semiarid climates where low elevations are dry and characterized by low carbon density ecosystems. However, these ecosystems are increasingly threatened by climate change with seasonal implications for photosynthesis and forest growth. As a result, we leveraged eddy covariance data from six evergreen conifer forest sites in the semiarid western United States to extrapolate the status of carbon sequestration within a framework of projected warming and drying. At colder locations, the seasonal evolution of gross primary productivity (GPP) was characterized by a single broad maximum during the summer that corresponded to snow melt‐derived moisture and a transition from winter dormancy to spring activity. Conversely, winter dormancy was transient at warmer locations, and GPP was responsive to both winter and summer precipitation such that two distinct GPP maxima were separated by a period of foresummer drought. This resulted in a predictable sequence of primary limiting factors to GPP beginning with air temperature in winter and proceeding to moisture and leaf area during the summer. Due to counteracting winter (positive) and summer (negative) GPP responses to warming, leaf area index and moisture availability were the best predictors of annual GPP differences across sites. Overall, mean annual GPP was greatest at the warmest site due to persistent vegetation photosynthetic activity throughout the winter. These results indicate that the trajectory of this region's carbon sequestration will be sensitive to reduced or delayed summer precipitation, especially if coupled to snow drought and earlier soil moisture recession, but summer precipitation changes remain highly uncertain. Given the demonstrated potential for seasonally offsetting responses to warming, we project that decadal semiarid montane forest carbon sequestration will remain relatively stable in the absence of severe disturbance.

Knowles, John F.↗

Improved Performance of Cu(InGa)(SeS) 2 PV Modules Using the Reaction of Metal Precursors. Final Report

This project “Improved Performance of Cu(InGa)(SeS) 2 PV Modules using the Reaction of Metal Precursors” was a partnership led by the Institute of Energy Conversion (IEC) at the University of Delaware with Columbia University and the Molecular Foundry at the Lawrence Berkeley National Laboratory. The aim was to develop pathways to improve Cu(InGa)(SeS) 2 (CIGSS) thin film photovoltaic modules using processes compatible with low manufacturing cost. The CIGSS approach investigated was a two-step process including deposition of metal precursor films following by reaction in hydride gases utilizing IEC’s novel reactor. The process was similar to that under commercial development by the project’s industry partner Stion. When Stion went out of business mid-project the focus changed to a rapid thermal process considered more commercially viable. Approaches to improve the performance of solar cells using the reacted films focused on two material innovations. First, the overall Ga content was increased to increase the operating voltage, which is desirable for scale-up to commercial modules. Second, the processing and performance advantages arising from Ag alloying were investigated. Advanced characterization guided process and material development including control of relative composition gradients. Research on the formation of Cu-Ga-In metal precursors utilized sputtering deposition which is normally used in commercial applications. The work resulted in processes for deposition of precursor stacks with increased relative Ga content and effects of deposition parameters on morphology and phase composition were established. It was shown that the metal precursor films have comparable phase composition and morphology so subsequent reaction follows from the same starting point. The addition of Ag to the metal precursors gave more uniform morphology and improved adhesion of reacted films which enable higher reaction temperature for faster processing. A novel outcome was the discovery of a previously undocumented material phase in sputter-deposited and evaporated Ag-Cu-In-Ga thin films. Hydride gas reaction processes including time-temperature-concentration profiles were developed for different precursor compositions. This enables control of composition profiles to engineer through-film gradients for solar cell optimization with characterization and simulations used to correlate measured film composition profiles to measurements of devices. In particular, the gradient of sulfur at the front of the CIGSS film was found to be critical. The simulations guided process development leading to improved reproducibility of devices improved performance with higher Ga content and higher voltage. With Ag-alloyed precursors, the reaction pathways leading were determined. A significant finding was that Ag-alloying increases the reaction rate to completely convert precursor films to the final chalcopyrite which could enable reduced reaction time to benefit manufacturability. To maintain potential commercial viability, the process under investigation was refocused to a rapid thermal process that could potentially be incorporated into an in-line process for manufacturing. Precursors with different composition were capped with an extra selenium layer and reacted in hydrogen sulfide 5-15 minutes, compared to typically 2 hours in the previous multi-step batch process. Critical RTP parameters were identified to control the reaction. Further optimization would be needed for high efficiency solar cells but pathways to high quality devices with further optimization and improved heating uniformity were developed. The project also developed new optoelectronic characterization approaches with a focus on development and application of spatial- and time-resolved photoluminescence and a custom mapping photoluminescence microscope built. It was shown how critical electronic transport properties strongly depend on the chemical composition of the material and that a wide range of samples show inhomogeneity on a length scale larger than the grains in the films. Additionally, two-photon excitation capability was developed to distinguish bulk vs surface losses. The project advances the state-of-the -art for precursor reaction processes in several ways that could impact manufacturing. This includes validation of approaches to increase voltage and establishment of model-guided control to form optimal composition profiles. The application of process control approaches with knowledge of phase formation and reaction pathways can be critically valuable in designing a large-scale process.

14 SOLAR ENERGY↗

The Low-Beta Quadrupoles for the High-Luminosity Large Hadron Collider: enabling higher luminosity by using Nb_3Sn conductor

The large-aperture, high-gradient Low Beta Quadrupoles for the High Luminosity Large Hadron Collider at CERN are the cornerstones of this collider. They are first-of-the-kind accelerator magnets bringing Nb_3Sn conductor to High Energy Physics colliders. In the US they are being fabricated and tested by the HL-LHC Accelerator Upgrade Project (AUP), and in Europe by CERN. The design was developed by AUP and CERN in close collaboration. The Q1/Q3 cryo-assemblies contain two 4.2 m long magnets (MQXFA) fabricated and by AUP. They are tested vertically before use in a cry-assembly. In this paper we are going to present the main features of the MQXFA design, results of vertical tests, and main lessons learned. We will also discuss lessons learned for future use of Nb_3Sn magnets in HEP colliders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Techno-Economic Simulation Results Using dGeo for EGS-Based District Heating in the Northeastern United States

This dataset presents the results of techno-economic simulations performed using the Distributed Geothermal Market Demand Model (dGeo) to evaluate the feasibility of Enhanced Geothermal Systems (EGS)-based district heating in the Northeastern United States. Developed by the National Renewable Energy Laboratory (NREL), dGeo is a geospatially resolved, bottom-up modeling framework designed to explore the deployment potential of geothermal distributed energy resources. The dataset, created as part of the Cornell EGS Ground-Truthing Project, provides census tract-level data that includes inputs and outputs such as thermal demand, road length, energy prices, geothermal system sizing, annual energy contributions from geothermal and natural gas peaking boilers, system capital costs (CAPEX), operation and maintenance costs (OPEX), and the levelized cost of heat (LCOH). Key simulation parameters include geothermal gradients, measured well depths, production temperatures, and district heating piping lengths based on S1400 neighborhood road lengths. The simulations assume a target bottom hole temperature of 80C and the development of new district heating networks in each census tract.

15 GEOTHERMAL ENERGY↗

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

RLAs WITH FFA ARCS FOR PROTONS AND ELECTRONS

Recirculating Linear Accelerators (RLAs) provide an efficient way of producing high-power, high-quality, continuous-wave hadron and lepton beams. However, their attractiveness had been limited by the cumbersomeness of multiple recirculating arcs and by the complexity of the spreader and recombiner regions. The latter problem sets one of the practical limitations on the maximum number of recirculations. We present an RLA design concept where the problem of multiple arcs is solved using the Fixed-Field Alternating gradient (FFA) design as in CBETA. The spreader/recombiner design is greatly simplified using an adiabatic matching approach. It allows for the spreader/recombiner function to be accomplished by a single beam line. The concept is applied to the designs of a high-power hadron accelerator being considered at ORNL and a CEBAF electron energy doubling project, FFA@CEBAF, being developed at Jefferson lab.

Morozov, V. S.↗

TEAM project review, year 1

This note summarizes our research activities within the TEAM project between October 2019 and November 2020, funded by the ASCR Advanced Research in Quantum Computing program. Three research efforts have begun during the review period. The first effort, resulting in the JuQBox software, implemented a novel algorithm for optimal control in closed quantum systems based on symplectic time integration and exact calculations of the gradient of the discretized objective function. In our second effort, resulting in a code called Quandary, we consider the optimal control problem for open quantum systems. Here the evolution of the quantum system is modeled by the Lindblad master equation. In this code, we are considering both the realization of quantum gates as well as the optimal reset problem in which a general mixed quantum state is driven towards its ground state. In our third research effort we have begun the development of computational software for what we call Lindblad learning, i.e. the process of characterizing the equations governing our quantum computing hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Feasibility Evaluation of a Solid Phase Scalable HEA Cladding Manufacturing Route

First of a kind development result on two low-energy solid-phase processes applied on an irradiation-resistant alloy, NiCoFeCrCu 0.12 , are achieved and demonstrate moderate feasibility of successful tube fabrication using shear assisted processing and extrusion (ShAPE™) and friction stir layer deposition as a bulk manufacturing process. The scope of the work is performed in four phases: 1) direct tube manufacturing of the irradiation-resistant high-entropy alloy (HEA) composite with increased strength, 2) co-shear lining manufacturing process for the increased strength and corrosion-resistant, irradiation-tolerant HEAs, 3) ShAPE of the radially gradient corrosion resistance alloy, and 4) alloy development and fabrication enabled through friction stir additive manufacturing processes among others. This report describes the development activities from April to December 2023 to manufacture a direct customizable thin-walled tubular product from irradiation-tolerant composite high-entropy alloys (C-HEAs) while the overall project is continuing in 2024.

36 MATERIALS SCIENCE↗

A Measurement of Atmospheric Circular Polarization with POLARBEAR

At millimeter wavelengths, the atmospheric emission is circularly polarized owing to the Zeeman splitting of molecular oxygen by the Earth's magnetic field. We report a measurement of the signal in the 150 GHz band using 3 yr of observational data with the POLARBEAR project. Nonidealities of a continuously rotating half-wave plate (HWP) partially convert circularly polarized light to linearly polarized light. While POLARBEAR detectors are sensitive to linear polarization, this effect makes them sensitive to circular polarization. Although this was not the intended use, we utilized this conversion to measure circular polarization. We reconstruct the azimuthal gradient of the circular polarization signal and measure its dependency from the scanning direction and the detector bandpass. We compare the signal with a simulation based on atmospheric emission theory, the detector bandpass, and the HWP leakage spectrum model. We find the ratio of the observed azimuthal slope to the simulated slope is 0.92 ± 0.01(stat) ± 0.07(sys). This ratio corresponds to a brightness temperature of 3.8 mK at the effective band center of 121.8 GHz and bandwidth of 3.5 GHz estimated from representative detector bandpass and the spectrum of Zeeman emission. This result validates our understanding of the instrument and reinforces the feasibility of measuring the circular polarization using the imperfection of continuously rotating HWP. Continuously rotating HWP is popular in ongoing and future cosmic microwave background experiments to modulate the polarized signal. This work shows a method for signal extraction and leakage subtraction that can help measure circular polarization in such experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Data and scripts associated with a manuscript analyzing ELM-FATES parameter sensitivity under pre-fire and postfire scenarios using machine learning

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).

Aboveground biomass↗

Aero-Optics of Hypersonic Turbulent Boundary Layers

Aero-optics refers to optical distortions due to index-of-refraction gradients that are induced by aerodynamic density gradients. At hypersonic flow conditions, the bulk velocity is many times the speed of sound and density gradients may originate from shock waves, compressible turbulent structures, acoustic waves, thermal variations, etc. Due to the combination of these factors, aero-optic distortions are expected to differ from those common to sub-sonic and lower super-sonic speeds. This report summarizes the results from a 2019-2022 Laboratory Directed Research and Development (LDRD) project led by Sandia National Laboratories in collaboration with the University of Notre Dame, New Mexico State University, and the Georgia Institute of Technology. Efforts extended experimental and simulation methodologies for the study of turbulent hypersonic boundary layers. Notable experimental advancements include development of spectral de-aliasing techniques for highspeed wavefront measurements, a Spatially Selective Wavefront Sensor (SSWFS) technique, new experimental data at Mach 8 and 14, a Quadrature Fringe Imaging Interferometer (QFII) technique for time-resolved index-of-refraction measures, and application of QFII to shock-heated air. At the same time, model advancements include aero-optic analysis of several Direct Numerical Simulation (DNS) datasets from Mach 0.5 to 14 and development of wall-modeled Large Eddy Simulation (LES) techniques for aero-optic predictions. At Mach 8 measured and predicted root mean square Optical Path Differences agree within confidence bounds but are higher than semi-empirical trends extrapolated from lower Mach conditions. Overall, results show that aero-optic effects in the hypersonic flow regime are not simple extensions from prior knowledge at lower speeds and instead reflect the added complexity of compressible hypersonic flow physics.

36 MATERIALS SCIENCE↗

NGEE Arctic Integrated Modeling (IM2): Improved subgrid hillslope hydrologic connectivity

This data product represents the integration of new code capability for arctic tundra hillslope hydrologic processes into the Energy Exascale Earth System Model (E3SM), through the E3SM Land Model (ELM) component. This code integration is the result of collaborative effort between the NGEE Arctic project and the E3SM project. The current ELM represents water movement primarily through vertical processes, such as precipitation, canopy interception, evaporation, infiltration, and soil water movement. Lateral water movement—such as surface runoff, subsurface flow, and river transport—plays a significant role in the hydrological cycle, especially in regions with varied topography. While E3SM includes a runoff routing component representing water transport in the river network, the lateral transport of water at the subgrid scale within the land model has previously not been taken into account. With the recent development of topographic units within the ELM subgrid data structure, there is an opportunity to simulate hillslope hydrologic connectivity by introducing water transport along topographic gradients. We expect that more realistic representation of hillslope hydrologic processes will lead to improved predictions of both soil water content and river network flows. Lateral transport of water at and near the surface is represented as a sub-grid process in this new code development. Water is tracked as it moves from higher to lower elevations within a gridcell. This capability uses the nested hierarchical sub-grid scheme within ELM to connect water fluxes from sub-grid elements with higher elevation to those with lower elevation. This data record consists of a single document (pdf format) that describes the theoretical basis for the hillslope hydrology processes added to ELM, and describes the modifications made to the ELM code. The Methods section of this metadata record includes a link to the public E3SM code repository where the exact code modifications as integrated in E3SM can be accessed. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

Thornton, Peter E [ORNL] (ORCID:0000000247595158)↗

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↗