Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial variation”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Constraining Relativistic Bow Shock Properties in Rotation-powered Millisecond Pulsar Binaries

Multiwavelength follow-up of unidentified Fermi sources has vastly expanded the number of known galactic-field “black widow” and “redback” millisecond pulsar binaries. Focusing on their rotation-powered state, we interpret the radio to X-ray phenomenology in a consistent framework. We advocate the existence of two distinct modes differing in their intrabinary shock orientation, distinguished by the phase centering of the double-peaked X-ray orbital modulation originating from mildly relativistic Doppler boosting. By constructing a geometric model for radio eclipses, we constrain the shock geometry as functions of binary inclination and shock standoff R {sub 0}. We develop synthetic X-ray synchrotron orbital light curves and explore the model parameter space allowed by radio eclipse constraints applied on archetypal systems B1957+20 and J1023+0038. For B1957+20, from radio eclipses the standoff is R {sub 0} ∼ 0.15–0.3 fraction of binary separation from the companion center, depending on the orbit inclination. Constructed X-ray light curves for B1957+20 using these values are qualitatively consistent with those observed, and we find occultation of the shock by the companion as a minor influence, demanding significant Doppler factors to yield double peaks. For J1023+0038, radio eclipses imply R {sub 0} ≲ 0.4, while X-ray light curves suggest 0.1 ≲ R {sub 0} ≲ 0.3 (from the pulsar). Degeneracies in the model parameter space encourage further development to include transport considerations. Generically, the spatial variation along the shock of the underlying electron power-law index should yield energy dependence in the shape of light curves, motivating future X-ray phase-resolved spectroscopic studies to probe the unknown physics of pulsar winds and relativistic shock acceleration therein.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Modulation of Galactic Cosmic Rays in the Inner Heliosphere, Comparing with PAMELA Measurements

We develop a numerical model to study the time-dependent modulation of galactic cosmic rays in the inner heliosphere. In the model, a time-delayed modified Parker heliospheric magnetic field (HMF) and a new diffusion coefficient model, NLGCE-F, from Qin and Zhang, are adopted. In addition, the latitudinal dependence of magnetic turbulence magnitude is assumed to be ∼(1+sin{sup 2}θ)/2 from the observations of Ulysses, and the radial dependence is assumed to be ∼r{sup S}, where we choose an expression of S as a function of the heliospheric current sheet tilt angle. We show that the analytical expression used to describe the spatial variation of HMF turbulence magnitude agrees well with the Ulysses, Voyager 1, and Voyager 2 observations. By numerically calculating the modulation code, we get the proton energy spectra as a function of time during the recent solar minimum, it is shown that the modulation results are consistent with the Payload for Antimatter-Matter Exploration and Light-nuclei Astrophysics measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

Abstract We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Parallel Diffusion Coefficient of Energetic Charged Particles in the Inner Heliosphere from the Turbulent Magnetic Fields Measured by Parker Solar Probe

Diffusion coefficients of energetic charged particles in turbulent magnetic fields are a fundamental aspect of diffusive transport theory but remain incompletely understood. In this work, we use quasi-linear theory to evaluate the spatial variation of the parallel diffusion coefficient κ ∥ from the measured magnetic turbulence power spectra in the inner heliosphere. We consider the magnetic field and plasma velocity measurements from Parker Solar Probe made during Orbits 5–13. The parallel diffusion coefficient is calculated as a function of radial distance from 0.062 to 0.8 au, and the particle energy from 100 keV to 1 GeV. We find that κ ∥ increases exponentially with both heliocentric distance and energy of particles. The fluctuations in κ ∥ are related to the episodes of large-scale magnetic structures in the solar wind. By fitting the results, we also provide an empirical formula of κ ∥ = (5.16 ± 1.22) × 10 18 r 1.17 ± 0.08 E 0.71 ± 0.02 (cm 2 s -1 ) in the inner heliosphere, which can be used as a reference in studying the transport and acceleration of solar energetic particles as well as the modulation of cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Testing the LSST Difference Image Analysis Pipeline Using Synthetic Source Injection Analysis

Abstract We evaluate the performance of the Legacy Survey of Space and Time Science Pipelines Difference Image Analysis (DIA) on simulated images. By adding synthetic sources to galaxies on images, we trace the recovery of injected synthetic sources to evaluate the pipeline on images from the Dark Energy Science Collaboration Data Challenge 2. The pipeline performs well, with efficiency and flux accuracy consistent with the signal-to-noise ratio of the input images. We explore different spatial degrees of freedom for the Alard–Lupton polynomial-Gaussian image subtraction kernel and analyze for trade-offs in efficiency versus artifact rate. Increasing the kernel spatial degrees of freedom reduces the artifact rate without loss of efficiency. The flux measurements with different kernel spatial degrees of freedom are consistent. We also here provide a set of DIA flags that substantially filter out artifacts from the DIA source table. We explore the morphology and possible origins of the observed remaining subtraction artifacts and suggest that given the complexity of these artifact origins, a convolution kernel with a set of flexible bases with spatial variation may be needed to yield further improvements.

Liu, S. (ORCID:0000000244612143)↗
DOE OSTI2024

Flat Spectra of Energetic Particles in Interplanetary Shock Precursors

The observed energy spectra of accelerated particles at interplanetary shocks often do not match the diffusive shock acceleration (DSA) theory predictions. In some cases, the particle flux forms a plateau over a wide range of energies, extending upstream of the shock for up to seven flux e-folds before submerging into the background spectrum. Remarkably, at and downstream of the shock we have studied in detail, the flux falls off in energy as ϵ -1 , consistent with the DSA prediction for a strong shock. The upstream plateau suggests a particle transport mechanism different from those traditionally employed in DSA models. We show that a standard (linear) DSA solution based on a widely accepted diffusive particle transport with an underlying resonant wave–particle interaction is inconsistent with the plateau in the particle flux. To resolve this contradiction, we modify the DSA theory in two ways. First, we include a dependence of the particle diffusivity κ on the particle flux F (nonlinear particle transport). Second, we invoke short-scale magnetic perturbations that are self-consistently generated by, but not resonant with, accelerated particles. They lead to the particle diffusivity increasing with the particle energy as ∝ϵ 3/2 that simultaneously decreases with the particle flux as 1/F. The combination of these two trends results in the flat spectrum upstream. We speculate that nonmonotonic spatial variations of the upstream spectrum, apart from being time-dependent, may also result from non-DSA acceleration mechanisms at work upstream, such as stochastic Fermi or magnetic pumping acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Revealing the Nano-Scale Structure and Properties of Pinholes in SiOx Layers for POLO Contacts

In this contribution, we report on the charge collection characteristics of silicon photovoltaic devices with passivated contacts based on the c-Si/SiOx /poly-Si structure. Using electron-beam induced current (EBIC) imaging in plan-view and cross-section orientations, we find that charge collection in a device with a 1.5-nm-thick SiOx layer is fairly uniform in the p-n junction region and does not appear to be influenced by pyramidal surface texture. In contrast, a device with a 2.2-nm-thick oxide layer shows significant spatial variation in the charge-collection signal. The locations where we observe preferential collection of excess carriers is very sensitive to the processing conditions used to grow and breakup the SiOx passivation layer. In some cases the apexes of the pyramids exhibit reduced EBIC signal as the valleys between adjacent pyramids show the strongest collection. Other conditions yield enhanced charge collection at well defined spots on pyramid faces with no observable collection in the valleys.

EBIC↗
DOE OSTI2013

Uncertainty in Modeling Dust Mass Balance and Radiative Forcing from Size Parameterization

This study examines the uncertainties in simulating mass balance and radiative forcing of mineral dust due to biases in the aerosol size parameterization. Simulations are conducted quasi-globally (180oW-180oE and 60oS-70oN) using the WRF24 Chem model with three different approaches to represent aerosol size distribution (8-bin, 4-bin, and 3-mode). The biases in the 3-mode or 4-bin approaches against a relatively more accurate 8-bin approach in simulating dust mass balance and radiative forcing are identified. Compared to the 8-bin approach, the 4-bin approach simulates similar but coarser size distributions of dust particles in the atmosphere, while the 3-mode pproach retains more fine dust particles but fewer coarse dust particles due to its prescribed og of each mode. Although the 3-mode approach yields up to 10 days longer dust mass lifetime over the remote oceanic regions than the 8-bin approach, the three size approaches produce similar dust mass lifetime (3.2 days to 3.5 days) on quasi-global average, reflecting that the global dust mass lifetime is mainly determined by the dust mass lifetime near the dust source regions. With the same global dust emission (~6000 Tg yr-1), the 8-bin approach produces a dust mass loading of 39 Tg, while the 4-bin and 3-mode approaches produce 3% (40.2 Tg) and 25% (49.1 Tg) higher dust mass loading, respectively. The difference in dust mass loading between the 8-bin approach and the 4-bin or 3-mode approaches has large spatial variations, with generally smaller relative difference (<10%) near the surface over the dust source regions. The three size approaches also result in significantly different dry and wet deposition fluxes and number concentrations of dust. The difference in dust aerosol optical depth (AOD) (a factor of 3) among the three size approaches is much larger than their difference (25%) in dust mass loading. Compared to the 8-bin approach, the 4-bin approach yields stronger dust absorptivity, while the 3-mode approach yields weaker dust absorptivity. Overall, on quasi-global average, the three size parameterizations result in a significant difference of a factor of 2~3 in dust surface cooling (-1.02~-2.87 W m-2) and atmospheric warming (0.39~0.96 W m-2) and in a tremendous difference of a factor of ~10 in dust TOA cooling (-0.24~-2.20 W m-2). An uncertainty of a factor of 2 is quantified in dust emission estimation due to the different size parameterizations. This study also highlights the uncertainties in modeling dust mass and number loading, deposition fluxes, and radiative forcing resulting from different size parameterizations, and motivates further investigation of the impact of size parameterizations on modeling dust impacts on air quality, climate, and ecosystem.

balance↗
DOE OSTI2014

What is the importance of climate model bias when projecting the impacts of climate change on land surface processes?

Regional climate change impact (CCI) studies have widely involved downscaling and bias-correcting (BC) Global Climate Model (GCM)-projected climate for driving land surface models. However, BC may cause uncertainties in projecting hydrologic and biogeochemical responses to future climate due to the impaired spatiotemporal covariance of climate variables and a breakdown of physical conservation principles. Here we quantify the impact of BC on simulated climate-driven changes in water variables(evapotranspiration, ET; runoff; snow water equivalent, SWE; and water demand for irrigation), crop yield, biogenic volatile organic compounds (BVOC), nitric oxide (NO) emissions, and dissolved inorganic nitrogen (DIN) export over the Pacific Northwest (PNW) Region. We also quantify the impacts on net primary production (NPP) over a small watershed in the region (HJ Andrews). Simulation results from the coupled ECHAM5/MPI-OM model with A1B emission scenario were firstly dynamically downscaled to 12 km resolutions with WRF model. Then a quantile mapping based statistical downscaling model was used to downscale them into 1/16th degree resolution daily climate data over historical and future periods. Two series climate data were generated according to the option of bias-correction (i.e. with bias-correction (BC) and without bias-correction, NBC). Impact models were then applied to estimate hydrologic and biogeochemical responses to both BC and NBC meteorological datasets. These im20 pact models include a macro-scale hydrologic model (VIC), a coupled cropping system model (VIC-CropSyst), an ecohydrologic model (RHESSys), a biogenic emissions model (MEGAN), and a nutrient export model (Global-NEWS). Results demonstrate that the BC and NBC climate data provide consistent estimates of the climate-driven changes in water fluxes (ET, runoff, and water demand), VOCs (isoprene and monoterpenes) and NO emissions, mean crop yield, and river DIN export over the PNW domain. However, significant differences rise from projected SWE, crop yield from dry lands, and HJ Andrews’s ET between BC and NBC data. Even though BC post-processing has no significant impacts on most of the studied variables when taking PNW as a whole, their effects have large spatial variations and some local areas are substantially influenced. In addition, there are months during which BC and NBC post-processing produces significant differences in projected changes, such as summer runoff. Factor-controlled simulations indicate that BC post-processing of precipitation and temperature both substantially contribute to these differences at region scales. We conclude that there are trade-offs between using BC climate data for offline CCI studies vs. direct modeled climate data. These trade-offs should be considered when designing integrated modeling frameworks for specific applications; e.g., BC may be more important when considering impacts on reservoir operations in mountainous watersheds than when investigating impacts on biogenic emissions and air quality (where VOCs are a primary indicator).

Liu, M. L.↗
DOE OSTI2016

Quantifying soil carbon accumulation in Alaskan terrestrial ecosystems during the last 15 000 years

Abstract. Northern high latitudes contain large amounts of soil organic carbon (SOC), of which Alaskan terrestrial ecosystems account for a substantial proportion. In this study, the SOC accumulation in Alaskan terrestrial ecosystems over the last 15 000 years was simulated using a process-based biogeochemistry model for both peatland and non-peatland ecosystems. Comparable with the previous estimates of 25–70 Pg C in peatland and 13–22 Pg C in non-peatland soils within 1 m depth in Alaska using peat-core data, our model estimated a total SOC of 36–63 Pg C at present, including 27–48 Pg C in peatland soils and 9–15 Pg C in non-peatland soils. Current vegetation stored 2.5–3.7 Pg C in Alaska, with 0.3–0.6 Pg C in peatlands and 2.2–3.1 Pg C in non-peatlands. The simulated average rate of peat C accumulation was 2.3 Tg C yr−1, with a peak value of 5.1 Tg C yr−1 during the Holocene Thermal Maximum (HTM) in the early Holocene, 4-fold higher than the average rate of 1.4 Tg C yr−1 over the rest of the Holocene. The SOC accumulation slowed down, or even ceased, during the neoglacial climate cooling after the mid-Holocene, but increased again in the 20th century. The model-estimated peat depths ranged from 1.1 to 2.7 m, similar to the field-based estimate of 2.29 m for the region. We found that the changes in vegetation and their distributions were the main factors in determining the spatial variations of SOC accumulation during different time periods. Warmer summer temperature and stronger radiation seasonality, along with higher precipitation in the HTM and the 20th century, might have resulted in the extensive peatland expansion and carbon accumulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Biogeochemical cycling at the aquatic–terrestrial interface is linked to parafluvial hyporheic zone inundation history

Abstract. The parafluvial hyporheic zone combines the heightened biogeochemical and microbial interactions indicative of a hyporheic region with direct atmospheric/terrestrial inputs and the effects of wet–dry cycles. Therefore, understanding biogeochemical cycling and microbial interactions in this ecotone is fundamental to understanding biogeochemical cycling at the aquatic–terrestrial interface and to creating robust hydrobiogeochemical models of dynamic river corridors. We aimed to (i) characterize biogeochemical and microbial differences in the parafluvial hyporheic zone across a small spatial domain (6 lateral meters) that spans a breadth of inundation histories and (ii) examine how parafluvial hyporheic sediments respond to laboratory-simulated re-inundation. Surface sediment was collected at four elevations along transects perpendicular to flow of the Columbia River, eastern WA, USA. The sediments were inundated by the river 0, 13, 127, and 398 days prior to sampling. Spatial variation in environmental variables (organic matter, moisture, nitrate, glucose, % C, % N) and microbial communities (16S and internal transcribed spacer (ITS) rRNA gene sequencing, qPCR) were driven by differences in inundation history. Microbial respiration did not differ significantly across inundation histories prior to forced inundation in laboratory incubations. Forced inundation suppressed microbial respiration across all histories, but the degree of suppression was dramatically different between the sediments saturated and unsaturated at the time of sample collection, indicating a binary threshold response to re-inundation. We present a conceptual model in which irregular hydrologic fluctuations facilitate microbial communities adapted to local conditions and a relatively high flux of CO2. Upon rewetting, microbial communities are initially suppressed metabolically, which results in lower CO2 flux rates primarily due to suppression of fungal respiration. Following prolonged inundation, the microbial community adapts to saturation by shifting composition, and the CO2 flux rebounds to prior levels due to the subsequent change in respiration. Our results indicate that the time between inundation events can push the system into alternate states: we suggest (i) that, above some threshold of inundation interval, re-inundation suppresses respiration to a consistent, low rate and (ii) that, below some inundation interval, re-inundation has a minor effect on respiration. Extending reactive transport models to capture processes that govern such dynamics will provide more robust predictions of river corridor biogeochemical function under altered surface water flow regimes in both managed and natural watersheds.

58 GEOSCIENCES↗
DOE OSTI2013

The distribution of soil phosphorus for global biogeochemical modeling

We discuss that phosphorus (P) is a major element required for biological activity in terrestrial ecosystems. Although the total P content in most soils can be large, only a small fraction is available or in an organic form for biological utilization because it is bound either in incompletely weathered mineral particles, adsorbed on mineral surfaces, or, over the time of soil formation, made unavailable by secondary mineral formation (occluded). In order to adequately represent phosphorus availability in global biogeochemistry–climate models, a representation of the amount and form of P in soils globally is required. We develop an approach that builds on existing knowledge of soil P processes and databases of parent material and soil P measurements to provide spatially explicit estimates of different forms of naturally occurring soil P on the global scale. We assembled data on the various forms of phosphorus in soils globally, chronosequence information, and several global spatial databases to develop a map of total soil P and the distribution among mineral bound, labile, organic, occluded, and secondary P forms in soils globally. The amount of P, to 50cm soil depth, in soil labile, organic, occluded, and secondary pools is 3.6 ± 3, 8.6 ± 6, 12.2 ± 8, and 3.2 ± 2 Pg P (Petagrams of P, 1 Pg = 1 × 10 15 g) respectively. The amount in soil mineral particles to the same depth is estimated at 13.0 ± 8 Pg P for a global soil total of 40.6 ± 18 Pg P. The large uncertainty in our estimates reflects our limited understanding of the processes controlling soil P transformations during pedogenesis and a deficiency in the number of soil P measurements. In spite of the large uncertainty, the estimated global spatial variation and distribution of different soil P forms presented in this study will be useful for global biogeochemistry models that include P as a limiting element in biological production by providing initial estimates of the available soil P for plant uptake and microbial utilization.

58 GEOSCIENCES↗
DOE OSTI2023

Heliostat Consortium: Gap Analysis on Wind Load for Achieving a Fully Competitive Heliostat Industry

The Heliostat Consortium (HelioCon) Wind Load Subtask was initiated with the aim of bringing research work pertaining to wind load measurement, characterization, and prediction taking place across several tasks, including Advanced Manufacturing, Components and Controls, and Field Deployment. The cross-cutting wind load subtopic in the HelioCon roadmap report [1] highlighted standardized methods and tools that are needed for a more detailed understanding of the static and dynamic loads on a heliostat. This will enable cost reduction of wind-dependent heliostat components to avoid unnecessarily conservative, overly constrained designs and increase field efficiency/reliability, to reduce the risk of component failures due to high-wind events (>15 m/s). Gaps related to heliostat wind load include site characterization for wind measurements, critical load cases for heliostat design, turbulence impacts on heliostat tracking error, testing of heliostat array configurations, understanding spatial variation of maximum loads across the solar field, and heliostat field layout and operating strategies. The recommended highest-priority pathway as first steps taken by HelioCon to address these gaps are to develop site characterization guidelines for heliostat design and field load measurements.

Emes, Matthew (ORCID:0000000341474387)↗
DOE OSTI2017

An Assessment of the Economic Potential of Offshore Wind in the United States from 2015 to 2030

Output data from an NREL report entitled "An Assessment of the Economic Potential of Offshore Wind in the United States from 2015 to 2030" (NREL/TP-6A20-67675), which analyzes the spatial variation of levelized cost of energy (LCOE) and levelized avoided cost of energy (LACE) to understand the economic potential of fixed-bottom and floating offshore wind technologies across more than 7,000 U.S. coastal sites between 2015 and 2030.

17 WIND ENERGY↗
DOE OSTI2019

Current drive by high intensity, pulsed, electron cyclotron wave packets

The nonlinear interaction of electrons with a high intensity, spatially localized, Gaussian, electromagnetic wave packet, or beam, in the electron cyclotron range of frequencies is described by the relativistic Lorentz equation. There are two distinct sets of electrons that result from wave-particle interactions. One set of electrons is reflected by the ponderomotive force due to the spatial variation of the wave packet. The second set of electrons are energetic enough to traverse across the wave packet. Both sets of electrons can exchange energy and momentum with the wave packet. The trapping of electrons in plane waves, which are constituents of the Gaussian beam, leads to dynamics that is distinctly different from quasilinear modeling of wave-particle interactions. This paper illustrates the changes that occur in the electron motion as a result of the nonlinear interaction. The dynamical differences between electrons interacting with a wave packet composed of ordinary electromagnetic waves and electrons interacting with a wave packet composed of extraordinary waves are exemplified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Spurious Shear in Weak Lensing with LSST

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image {approx} 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to r {approx} 27.5, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, additive systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than {approx} 10{prime} in the single short exposures, which propagates into a spurious shear correlation function at the 10{sup -4}-10{sup -3} level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

A Hydro-mechanical Model and Analytical Solutions for Geomechanical Modeling of Carbon Dioxide Geological Sequestration

We present a hydro-mechanical model for geological sequestration of carbon dioxide. The model considers the poroelastic effects by taking into account the coupling between the geomechanical response and the fluid flow in greater detail. The simplified hydro-mechanical model includes the geomechanical part that relies on the linear elasticity, while the fluid flow is based on the Darcy’s law. Two parts were coupled using the standard linear poroelasticity. Analytical solutions for pressure field were obtained for a typical geological sequestration scenario. The model predicts the temporal and spatial variation of pressure field and effects of permeability and elastic modulus of formation on the fluid pressure distribution.

sequestration, geomechanics, Multi-scale, poro-ela↗
DOE OSTI2016

Characterizing U.S. Heat Demand for Potential Application of Geothermal Direct Use: Preprint

In this paper, we assess the U.S. demand for low-temperature thermal energy at the county resolution for four major end-use sectors: residential buildings, commercial buildings, manufacturing facilities, and agricultural facilities. Existing, publicly available data on the U.S. thermal demand market are characterized by coarse spatial resolution, with assessments typically at the state-level or larger. For many uses, these data are sufficient; however, our research was motivated by an interest in assessing the potential demand for direct use (DU) of low-temperature (30 degrees to 150 degrees C) geothermal heat. The availability and quality of geothermal resources for DU applications are highly spatially heterogeneous; therefore, to assess the potential market for these resources, it is necessary to understand the spatial variation in demand for low-temperature resources at a local resolution. This paper presents the datasets and methods we used to develop county-level estimates of the thermal demand for the residential, commercial, manufacturing, and agricultural sectors. Although this analysis was motivated by an interest in geothermal energy deployment, the results are likely to have broader applications throughout the energy industry. The county-resolution thermal demand data developed in this study for four major U.S. sectors may have far-reaching implications for building technologies, industrial processes, and various distributed renewable energy thermal resources (e.g. biomass, solar).

15 GEOTHERMAL ENERGY↗