Search NASA⌕ Search

SEARCH · Search NASA

Results for “extrapolation method”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Practical quantum computation of chemical and nuclear energy levels using quantum imaginary time evolution and Lanczos algorithms

Various methods have been developed for the quantum computation of the ground and excited states of physical and chemical systems, but many of them require either large numbers of ancilla qubits or high-dimensional optimization in the presence of noise. The quantum imaginary-time evolution (QITE) and quantum Lanczos (QLanczos) methods proposed in Motta et al. (2020) eschew the aforementioned issues. In this study, we demonstrate the practical application of these algorithms to challenging quantum computations of relevance for chemistry and nuclear physics, using the deuteron-binding energy and molecular hydrogen binding and excited state energies as examples. With the correct choice of initial and final states, we show that the number of timesteps in QITE and QLanczos can be reduced significantly, which commensurately simplifies the required quantum circuit and improves compatibility with NISQ devices. We have performed these calculations on cloud-accessible IBM Q quantum computers. With the application of readout-error mitigation and Richardson error extrapolation, we have obtained ground and excited state energies that agree well with exact results obtained from diagonalization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Annealing of ion tracks in apatite under pressure characterized in situ by small angle x-ray scattering

Fission track thermochronology is routinely used to investigate the thermal history of sedimentary basins, as well as tectonic uplift and denudation rates. While the effect of temperature on fission track annealing has been studied extensively to calibrate the application of the technique, the effect of pressure during annealing is generally considered to be negligible. However, a previous study suggested elevated pressure results in a significantly different annealing behaviour that was previously unknown. Here, we present a method to study track annealing in situ under high pressure by using synchrotron- based small angle x-ray scattering (SAXS). To simulate fission tracks in a controlled environment, ion tracks were created in apatite from Durango, Mexico using 2 GeV Au or Bi ions provided by an ion accelerator facility. Samples were annealed at 250°C at approximately 1 GPa pressure using diamond anvil cells (DAcs) with heating capabilities. Additional in situ annealing experiments at ambient pressure and temperatures between 320 and 390°C were performed for comparison. At elevated pressure a significantly accelerated annealing rate of the tracks was observed compared with annealing at ambient pressure. However, when extrapolated to geologically relevant temperatures and pressures, the effects become very small. The measurement methodology presented provides a new avenue to study materials behaviour in extreme environments.

58 GEOSCIENCES↗
DOE OSTI2023

Engineering-Scale Integrated Energy System Data Projection Demonstration via the Dynamic Energy Transport and Integration Laboratory

The objective of this study is to demonstrate and validate the Dynamic Energy Transport and Integration Laboratory (DETAIL) preliminary scaling analysis using Modelica language system-code Dymola. The DETAIL preliminary scaling analysis includes a multisystem integral scaling package between thermal-storage and hydrogen-electrolysis systems. To construct the system of scaled equations, dynamical system scaling (DSS) was applied to all governing laws and closure relations associated with the selected integral system. The existing Dymola thermal-energy distribution system (TEDS) facility and high-temperature steam electrolysis (HTSE) facility models in the Idaho National Laboratory HYBRID repository were used to simulate a test case and a corresponding scaled case for integrated system HYBRID demonstration and validation. The DSS projected data based on the test-case simulations and determined scaling ratios were generated and compared with scaled case simulations. The preliminary scaling analysis performance was evaluated, and scaling distortions were investigated based on data magnitude, sequence, and similarity. The results indicated a necessity to change the normalization method for thermal storage generating optimal operating conditions of 261 kW power and mass flow rate of 6.42 kg/s and the possibility of reselecting governing laws for hydrogen electrolysis to improve scaling predictive properties. To enhance system-scaling similarity for TEDS and HTSE, the requirement for scaling validation via physical-facility demonstration was identified.

08 HYDROGEN↗
DOE OSTI2024

Can Error Mitigation Improve Trainability of Noisy Variational Quantum Algorithms?

Variational Quantum Algorithms (VQAs) are often viewed as the best hope for near-term quantum advantage. However, recent studies have shown that noise can severely limit the trainability of VQAs, e.g., by exponentially flattening the cost landscape and suppressing the magnitudes of cost gradients. Error Mitigation (EM) shows promise in reducing the impact of noise on near-term devices. Thus, it is natural to ask whether EM can improve the trainability of VQAs. In this work, we first show that, for a broad class of EM strategies, exponential cost concentration cannot be resolved without committing exponential resources elsewhere. This class of strategies includes as special cases Zero Noise Extrapolation, Virtual Distillation, Probabilistic Error Cancellation, and Clifford Data Regression. Second, we perform analytical and numerical analysis of these EM protocols, and we find that some of them (e.g., Virtual Distillation) can make it harder to resolve cost function values compared to running no EM at all. As a positive result, we do find numerical evidence that Clifford Data Regression (CDR) can aid the training process in certain settings where cost concentration is not too severe. Our results show that care should be taken in applying EM protocols as they can either worsen or not improve trainability. On the other hand, our positive results for CDR highlight the possibility of engineering error mitigation methods to improve trainability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Constraining level densities through quantitative correlations with cross-section data

The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the accurate calculation of cross sections for the different reaction channels. Many common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account in these models: LD at neutron separation energy coming from average spacings of s- and p-wave resonances ( D 0 and D 1 , respectively) whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only in very specific regions of excitation energy, and for specific spins and parities. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for fitting and the determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD as a starting point as the HFB LD provide a more realistic spin and parity distributions than phenomenological models such as Gilbert-Cameron (GC). We then associate variations predicted by the HFB model with the structure observed in double-differential cross sections at low outgoing neutron energy, a region that is dominated by the LD input. We also use (n, p) on 56 Fe , as an example case where angle-integrated cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or values of D 0 or D 1 are unknown. Finally, it also predicts neutron-induced inelastic γ cross sections that in some cases can differ significantly from more standard phenomenological LD models such as GC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Homotypic aggregates contribute to heterogeneity in B cell fates due to an intrinsic gradient of stimulant exposure

Highlights: • CD40-signaled B cells in cultures form three-dimensional homotypic aggregates. • These tight aggregates hinder the free diffusion of large molecule stimulants. • The resulting concentration gradient of stimulants leads to heterogeneous cell fates. • We describe a flow cytometric approach to quantify this positional information. Monocultures of several cell types result in the formation of robust clusters called homotypic aggregates (HAs). How this physical aggregation affects cell fates in immune cell cultures, is poorly understood. We studied anti-CD40-stimulated primary B cell cultures, where cells assembled into large three-dimensional LFA1-driven HAs by 72 h. The dense packing in these aggregates restricts the infiltration of stimulants, such as antibodies, to cells inside the clusters. This creates a concentration gradient of stimulant availability across the cross-section of HAs. We describe a method to retain this positional information even after the disruption of HAs, for analysis by flow cytometry. Comparison of stage-specific cell-surface markers showed that the extent of stimulant-binding affected multiple fates non-uniformly. While germinal center and lineage markers were moderately upregulated, immunoglobulins and markers associated with memory were more than doubled in the peripheral cells binding more anti-CD40. These cells also experienced a strong repression of the plasma cell regulator Prdm1 and an upregulation of the oncogene Myc. Thus, cells at different locations in HAs are subjected to unequal doses of stimulants, leading to a hitherto unreported source of heterogeneity in cell fates. These findings can be extrapolated to understand the dose-dependent effects of stimulants in other three-dimensional cell clusters.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Design Criteria for High-Gradient Radio-Frequency Linacs

This article will review methods used at the SLAC National Accelerator Laboratory and other world accelerator laboratories to design high-gradient normal conducting accelerating structures. A quest for compact radio-frequency linacs fueled decades of studies toward a higher accelerating gradient. A major phenomena limiting the increase of the gradient is vacuum radio-frequency breakdown; therefore, this paper will address the breakdown physics and discuss approaches that reduce the breakdown probability. This discussion will cover both the electrical design and fabrication technology of the accelerating structures to achieve practical operating accelerating gradients in excess of 100 MV/m. Most of the data described here were obtained during the development of 11 GHz linacs for electron–positron linear colliders, so extrapolation of the results to other frequencies should be performed cautiously.

Dolgashev, Valery A. (ORCID:0000000166715402)↗
DOE OSTI2021

Synthetic galaxy clusters and observations based on Dark Energy Survey Year 3 Data

ABSTRACT We develop a novel data-driven method for generating synthetic optical observations of galaxy clusters. In cluster weak lensing, the interplay between analysis choices and systematic effects related to source galaxy selection, shape measurement, and photometric redshift estimation can be best characterized in end-to-end tests going from mock observations to recovered cluster masses. To create such test scenarios, we measure and model the photometric properties of galaxy clusters and their sky environments from the Dark Energy Survey Year 3 (DES Y3) data in two bins of cluster richness $\lambda \in [30; 45)$, $\lambda \in [45; 60)$ and three bins in cluster redshift ($z\in [0.3; 0.35)$, $z\in [0.45; 0.5)$ and $z\in [0.6; 0.65)$. Using deep-field imaging data, we extrapolate galaxy populations beyond the limiting magnitude of DES Y3 and calculate the properties of cluster member galaxies via statistical background subtraction. We construct mock galaxy clusters as random draws from a distribution function, and render mock clusters and line-of-sight catalogues into synthetic images in the same format as actual survey observations. Synthetic galaxy clusters are generated from real observational data, and thus are independent from the assumptions inherent to cosmological simulations. The recipe can be straightforwardly modified to incorporate extra information, and correct for survey incompleteness. New realizations of synthetic clusters can be created at minimal cost, which will allow future analyses to generate the large number of images needed to characterize systematic uncertainties in cluster mass measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Melting and phase relations of Fe-Ni-Si determined by a multi-technique approach

Many studies have suggested silicon as a candidate light element for the cores of Earth and Mercury. However, the effect of silicon on the melting temperatures of core materials and thermal profiles of cores is poorly understood, due to disagreements among melt detection techniques, uncertainties in sample pressure evolution during heating, and sparsity of studies investigating the combined effects of nickel and silicon on the phase diagram of iron. Here, in this study we develop a multi-technique approach for measuring the high-pressure melting and solid phase relations of iron alloys and apply it to Fe 0.8 Ni 0.1 Si 0.1 (Fe-11wt%Ni-5.3wt%Si), a composition compatible with recent estimates for the cores of Earth and Mercury. This approach combines results (20-83 GPa) from two atomic-level techniques: synchrotron Mossbauer spectroscopy (SMS) and synchrotron x-ray diffraction (XRD). Melting is independently detected by the loss of the Mossbauer signal, produced exclusively by solid-bound iron nuclei, and the onset of a liquid diffuse x-ray scattering signal. The use of a burst heating and background updating method for quantifying changes in the reference background during heating facilitates the determination of liquid diffuse signal onsets and leads to strong reproducibility and excellent agreement in melting temperatures determined separately by the two techniques. XRD measurements additionally constrain the hcp-fcc phase boundary and in-situ pressure evolution of the samples during heating. We apply our updated thermal pressure model to published SMS melting data on fcc-Fe and fcc-Fe 0.9 Ni 0.1 to precisely evaluate the effect of silicon on melting temperatures. We find that the addition of 10 mol% Si to Fe 0.9 Ni 0.1 reduces melting temperatures by ~250 K at low pressures (<60 GPa) and flattens the hcp-fcc phase boundary. Extrapolating our results, we constrain the location of the hcp-fcc-liquid quasi-triple point at 147±14 GPa and 3140±90 K, which implies a melting temperature reduction of 500 K compared with Fe 0.9 Ni 0.1 . The results demonstrate the advantages of combining complementary experimental techniques in investigations of melting under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deep-learning-based canopy height model generation from sub-meter resolution panchromatic satellite imagery

Canopy height models (CHMs) with sufficient resolution to distinguish individual trees are useful for a variety of applications. However, standard techniques to acquire such data, such as airborne lidar surveying, are often prohibitively expensive. Deep learning techniques for generating CHMs from high-resolution imagery are an attractive option to reduce costs. To date, success with these methods has been demonstrated using multichannel aerial photography and specialized satellite data products derived from multiple sensors, neither of which is commonly available at temporal resolutions finer than one year. Here we demonstrate a method to generate sub-meter resolution CHMs in three forests in California using a more abundant data source: sub-meter resolution, panchromatic satellite imagery from a single sensor. We show that phenology and species composition play important roles in model transferability; when trained using imagery from a single conifer forest in autumn, the model performs well on autumn imagery from a second conifer forest several hundred kilometers distant with no re-training. With modest additions to the training dataset, the same model generates minimally biased estimates of canopy height in both conifer and deciduous forests during multiple seasons. Because the model operates on satellite data with global coverage and a relatively short return interval, we propose its suitability to extrapolate tree-level canopy height data to remote regions and conduct high-temporal resolution monitoring of forest structure. We furthermore demonstrate the workflow’s applicability to fire modeling by conducting simulations in forests populated by trees measured using both this approach and airborne lidar surveying. We find minimal differences in fire behavior relative to a baseline case in which only statistical distributions of tree height and crown area are known. This result underscores the value of forest structural information derived from our workflow for improving the fidelity of wildland fire simulations, among other ecological applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Scale setting the Möbius domain wall fermion on gradient-flowed HISQ action using the omega baryon mass and the gradient-flow scales 𝑡 0 and 𝑤 0

We report on a subpercent scale determination using the omega baryon mass and gradient-flow methods. The calculations are performed on 22 ensembles of 𝑁 𝑓 =2 +1 +1 highly improved, rooted staggered sea-quark configurations generated by the MILC and CalLat Collaborations. The valence quark action used is Möbius domain wall fermions solved on these configurations after a gradient-flow smearing is applied with a flowtime of 𝑡 gf = 1 in lattice units. The ensembles span four lattice spacings in the range 0.06 ≲ 𝑎 ≲0.15 fm, six pion masses in the range 130 ≲ 𝑚 𝜋 ≲ 400 MeV and multiple lattice volumes. On each ensemble, the gradient-flow scales 𝑡 0 /𝑎 2 and 𝑤 0 /𝑎 and the omega baryon mass 𝑎⁢𝑚 Ω are computed. The dimensionless product of these quantities is then extrapolated to the continuum and infinite volume limits and interpolated to the physical light, strange and charm quark mass point in the isospin limit, resulting in the determination of $\sqrt{t}_0$ = 0.1422⁢(14) fm and 𝑤 0 = 0.1709⁢(11) fm with all sources of statistical and systematic uncertainty accounted for. The dominant uncertainty in both results is the stochastic uncertainty, though for $\sqrt{t}_0$ there are comparable continuum extrapolation uncertainties. For 𝑤 0 , there is a clear path for a few-per-mille uncertainty just through improved stochastic precision, as recently obtained by the Budapest-Marseille-Wuppertal Collaboration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Generating high-resolution total canopy SIF emission from TROPOMI data: Algorithm and application

Solar-induced chlorophyll fluorescence (SIF) is a rapidly advancing front in modeling global terrestrial gross primary production (GPP). Canopy total SIF emissions (SIF total ) are mechanistically linked to the plant photosynthesis, and can be estimated from satellite observed SIF (SIF obs ) through radiative transfer modeling. However, the current satellite SIF obs and thus SIF total are available only at coarse spatial resolutions from several kilometers to tens of kilometers, inhibiting the application at fine spatial scales. Here, in this work, we proposed an algorithm to generate both global high-resolution SIF total (HSIF total ) and high-resolution SIF obs (HSIF obs ) at 1 km from low-resolution SIF obs (LSIF obs ) from the TROPOspheric Monitoring Instrument (TROPOMI), which has a spatial resolution at nadir of 3.5 km by 5.6–7 km. Our statistical method is based on the law of energy conservation and uses satellite derived fraction of absorbed photosynthetically active radiation, fluorescence efficiency, and the escape probability of fluorescence. We evaluated the accuracy of our HSIF total using the Orbiting Carbon Observatory-2 SIF (R 2 = 0.78). We found that the spatial resolution had clear effects on the relationship between HSIF total and GPP. We also compared HSIF total to 8-day averaged tower GPP from 135 flux sites and found that they were better correlated when HSIF total was averaged over a 1-km radius around the tower than when averaged over a larger radius. Our study provided a unique high-resolution HSIF total product, which will advance the estimation of GPP by extrapolating site-level relationships to the global scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

An algebraic convolution formulation for multiple-scattering correction in small-angle neutron scattering

Multiple scattering in small-angle neutron scattering (SANS) redistributes spectral weight and distorts structural interpretation, particularly for thick or strongly scattering samples. We develop a finite-dimensional spectral desmearing framework that corrects multiple scattering without resorting to integral transforms or model-dependent extrapolation. The primary intensity is expanded in an orthonormal basis adapted to the isotropic transverse-momentum measure, under which convolution reduces to a recursive tensor contraction, allowing the Poisson-weighted multiple-scattering series to be evaluated directly in a finite-dimensional basis representation. This formulation yields a stable forward–inverse mapping between apparent and primary spectra. Numerical tests demonstrate convergence under repeated convolution and accurate recovery of the single-scattering intensity. Application to SANS measurements collected at multiple neutron facilities, including the Spallation Neutron Source, the High Flux Isotope Reactor, and the Institut Laue-Langevin, shows the quantitative reconstruction of the underlying primary spectrum across a wide range of transmission conditions, including strongly attenuating samples. Here, the method provides a stable, model-agnostic framework for multiple-scattering correction in SANS and enables consistent structural interpretation across instruments and scattering regimes.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗
DOE OSTI2022

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗