Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data depth”

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 1,513 records · Page 84

Recent Results for AGN Observed by the Rossi X-Ray Timing Explorer

The Rossi X-ray Timing Explorer (RXTE) has produced many excellent observations of active galaxies, providing the best sensitivity in the 10 - 20 keV range so far. This presentation reports selected RTXE data for AGN in the context of the currently popular models. One is the recent result for two Seyfert 1 galaxies, NGC 5548 and IC4329a: both show the "canonical" Seyfert I X-ray spectra, with an underlying power law, plus Gaussian iron K line and Compton reflection. Interestingly, in both cases, the profile of the Fe K line does not extend as far to the red as seen in the famous NCG-6-30-15, and this indicates that the regions where the Fe K lines originate in AGN are diverse. Independently, in both objects we see a strong spectral variability of the primary continua, which soften as the sources brighten. The second result is for the heavily absorbed Seyfert 2 NGC 4945. The RXTE data confirm the strong absorption corresponding to the optical depth to electron scattering of about 2, but also reveal rapid variability of the hard (8-30 keV) X-ray emission on a time scale of a day or less. This suggests that for NGC 4945, the putative parsec-size molecular torus cannot be both geometrically and optically thick, and implies that the Cosmic X-ray Background is unlikely to be made up primarily of AGN with geometry as inferred for this object.

Madejski, G. M.↗

Estimates of Internal Tide Energy Fluxes from Topex/Poseidon Altimetry: Central North Pacific

Energy fluxes for first-mode M(sub 2) internal tides are deduced throughout the central North Pacific Ocean from Topex/Poseidon satellite altimeter data. Temporally coherent internal tide signals in the altimetry, combined with climatological hydrographic data, determine the tidal displacements, pressures, and currents at depth, which yield power transmission rates. For a variety of reasons the deduced rates should be considered lower bounds. Internal tides were found to emanate from several large bathymetric structures, especially the Hawaiian Ridge, where the integrated flux amounts to about six gigawatts. Internal tides are generated at the Aleutian Trench near 172 deg west and propagate southwards nearly 2000 km.

Ray, Richard D.↗

NASA Supercomputer Improves Prospects for Ocean Climate Research

Estimates of ocean circulation constrained by in situ and remotely sensed observations have become routinely available during the past five years, and they are being applied to myriad scientific and operational problems [Stammer et al.,2002]. Under the Global Ocean Data Assimilation Experiment (GODAE), several regional and global estimates have evolved for applications in climate research, seasonal forecasting, naval operations, marine safety, fisheries,the offshore oil industry, coastal management, and other areas. This article reports on recent progress by one effort, the consortium for Estimating the Circulation and Climate of the Ocean (ECCO), toward a next-generation synthesis of ocean and sea-ice data that is global, that covers the full ocean depth, and that permits eddies.

Estimating the Circulation and Climate of the Ocea↗

Simulation Studies of the Effect of Forest Spatial Structure on InSAR Signature

The height of scattering phase retrieved from InSAR data is considered being correlated with the tree height and the spatial structure of the forest stand. Though some researchers have used simple backscattering models to estimate tree height from the height of scattering center, the effect of forest spatial structure on InSAR data is not well understood yet. A three-dimensional coherent radar backscattering model for forest canopies based on realistic three-dimensional scene was used to investigate the effect in this paper. The realistic spatial structure of forest canopies was established either by field measurements (stem map) or through use of forest growth model. Field measurements or a forest growth model parameterized using local environmental parameters provides information of forest species composition and tree sizes in certain growth phases. A fractal tree model (L-system) was used to simulate individual 3- D tree structure of different ages or heights. Trees were positioned in a stand in certain patterns resulting in a 3-D medium of discrete scatterers. The radar coherent backscatter model took the 3-D forest scene as input and simulates the coherent radar backscattering signature. Interferometric SAR images of 3D scenes were simulated and heights of scattering phase centers were estimated from the simulated InSAR data. The effects of tree height, crown cover, crown depth, and the spatial distribution patterns of trees on the scattering phase center were analyzed. The results will be presented in the paper.

Sun, Guoqing↗

Sensitivity of Photolysis Frequencies and Key Tropospheric Oxidants in a Global Model to Cloud Vertical Distributions and Optical Properties

As a follow-up study to our recent assessment of the radiative effects of clouds on tropospheric chemistry, this paper presents an analysis of the sensitivity of such effects to cloud vertical distributions and optical properties in a global 3-D chemical transport model (GEOS4-Chem CTM). GEOS-Chem was driven with a series of meteorological archives (GEOS1-STRAT, GEOS-3 and GEOS-4) generated by the NASA Goddard Earth Observing System data assimilation system, which have significantly different cloud optical depths (CODs) and vertical distributions. Clouds in GEOS1-STRAT and GEOS-3 have more similar vertical distributions while those in GEOS-4 are optically much thinner in the tropical upper troposphere. We find that the radiative impact of clouds on global photolysis frequencies and hydroxyl radical (OH) is more sensitive to the vertical distribution of clouds than to the magnitude of column CODs. Model simulations with each of the three cloud distributions all show that the change in the global burden of O3 due to clouds is less than 5%. Model perturbation experiments with GEOS-3, where the magnitude of 3-D CODs are progressively varied by -100% to 100%, predict only modest changes (<5%) in global mean OH concentrations. J(O1D), J(NO2) and OH concentrations show the strongest sensitivity for small CODs and become insensitive at large CODs due to saturation effects. Caution should be exercised not to use in photochemical models a value for cloud single scattering albedo lower than about 0.999 in order to be consistent with the current knowledge of cloud absorption at the UV wavelength. Our results have important implications for model intercomparisons and climate feedback on tropospheric photochemistry.

Liu, Hongyu↗

A13K-0336: Airborne Multi-Wavelength High Spectral Resolution Lidar for Process Studies and Assessment of Future Satellite Remote Sensing Concepts

NASA Langley recently developed the world's first airborne multi-wavelength high spectral resolution lidar (HSRL). This lidar employs the HSRL technique at 355 and 532 nm to make independent, unambiguous retrievals of aerosol extinction and backscatter. It also employs the standard backscatter technique at 1064 nm and is polarization-sensitive at all three wavelengths. This instrument, dubbed HSRL-2 (the secondgeneration HSRL developed by NASA Langley), is a prototype for the lidar on NASA's planned Aerosols- Clouds-Ecosystems (ACE) mission. HSRL-2 completed its first science mission in July 2012, the Two-Column Aerosol Project (TCAP) conducted by the Department of Energy (DOE) in Hyannis, MA. TCAP presents an excellent opportunity to assess some of the remote sensing concepts planned for ACE: HSRL-2 was deployed on the Langley King Air aircraft with another ACE-relevant instrument, the NASA GISS Research Scanning Polarimeter (RSP), and flights were closely coordinated with the DOE's Gulfstream-1 aircraft, which deployed a variety of in situ aerosol and trace gas instruments and the new Spectrometer for Sky-Scanning, Sun-Tracking Atmospheric Research (4STAR). The DOE also deployed their Atmospheric Radiation Measurement Mobile Facility and their Mobile Aerosol Observing System at a ground site located on the northeastern coast of Cape Cod for this mission. In this presentation we focus on the capabilities, data products, and applications of the new HSRL-2 instrument. Data products include aerosol extinction, backscatter, depolarization, and optical depth; aerosol type identification; mixed layer depth; and rangeresolved aerosol microphysical parameters (e.g., effective radius, index of refraction, single scatter albedo, and concentration). Applications include radiative closure studies, studies of aerosol direct and indirect effects, investigations of aerosol-cloud interactions, assessment of chemical transport models, air quality studies, present (e.g., CALIPSO) and future (e.g., EarthCARE) satellite calibration/validation, and development/assessment of advanced retrieval techniques for future satellite applications (e.g., lidar+polarimeter retrievals of aerosol and cloud properties). We will also discuss the relevance of HSRL-2 measurement capabilities to the ACE remote sensing concept.

Hostetler, Chris A.↗

Video Altimeter and Obstruction Detector for an Aircraft

Video-based altimetric and obstruction detection systems for aircraft have been partially developed. The hardware of a system of this type includes a downward-looking video camera, a video digitizer, a Global Positioning System receiver or other means of measuring the aircraft velocity relative to the ground, a gyroscope based or other attitude-determination subsystem, and a computer running altimetric and/or obstruction-detection software. From the digitized video data, the altimetric software computes the pixel velocity in an appropriate part of the video image and the corresponding angular relative motion of the ground within the field of view of the camera. Then by use of trigonometric relationships among the aircraft velocity, the attitude of the camera, the angular relative motion, and the altitude, the software computes the altitude. The obstruction-detection software performs somewhat similar calculations as part of a larger task in which it uses the pixel velocity data from the entire video image to compute a depth map, which can be correlated with a terrain map, showing locations of potential obstructions. The depth map can be used as real-time hazard display and/or to update an obstruction database.

Delgado, Frank J.↗

The First Brown Dwarf Discovered by the Backyard Worlds: Planet 9 Citizen Science Project

The Wide-field Infrared Survey Explorer (WISE) is a powerful tool for finding nearby brown dwarfs and searching for new planets in the outer solar system, especially with the incorporation of NEOWISE and NEOWISE Reactivation data. However, so far, searches for brown dwarfs in WISE data have yet to take advantage of the full depth of the WISE images. To efficiently search this unexplored space via visual inspection, we have launched anew citizen science project, called "Backyard Worlds: Planet 9," which asks volunteers to examine short animations composed of difference images constructed from time-resolved WISE co adds. We report the first new substellar object discovered by this project, WISEA J110125.95+540052.8, a T5.5 brown dwarf located approximately 34 pc from the Sun with a total proper motion of approx.0. "7/ yr. WISEA J110125.95+540052.8 has a WISE W2 magnitude of W2 = 15.37+/- 0.09; our sensitivity to this source demonstrates the ability of citizen scientists to identify moving objects via visual inspection that are 0.9 mag fainter than the W2 single-exposure sensitivity, a threshold that has limited prior motion-based brown dwarf searches with WISE.

Wide-field Infrared Survey Explorer (WISE)↗

A Qualitative Investigation of Selected Infrared Flow Visualization Image Processing Techniques

The deployment and integration of high-sensitivity infrared cameras in a transonic wind tunnel testenvironment has resulted in a unique capability to image aerodynamic phenomena in real-time. Multi-camera infrared flow visualization data systems are now routinely utilized at the NASA Ames Unitary Plan Wind Tunnel. The small flow-induced temperature gradients on the surface of the wind tunnel test article coupled with the high bit-depth of the infrared camera sensor makes the processing of the image data critically important. An image processing routine must enhance features of interest with minimal artifacts. Additionally, the production wind tunnel test environment demands that these processed images are made available in a real-time, automatic fashion. Therefore, any image processing routine must be computationally economical and enhance the image data with minimal input from a human operator. The following seeks to qualitatively explore selected image processing techniques by assessing their effectiveness to resolve flow features on a wind tunnel test article. A multi-scale contrast enhancement technique is discussed as well as a new implementation of a multi-scale, non-interpolated adaptive histogram equalization. Finally, a novel method is introduced that demonstrates the ability to resolve flow features imaged on bare-steel test articles possessing low emissivity.This method makes use of dynamic mode decomposition and discrete-time filtering to separate the background reflections that dominate low emissivity surfaces from the aerodynamic driven surface temperature gradients.This process will be shown to resolve the onset of boundary layer transition on a bare metal wing as well as identify and resolve hidden features in the image data. While the implementation of this technique is very preliminary it demonstrates the potential to extend the application of infrared flow-visualization within the wind tunnel test environment.

Image Processeing↗

A Qualitative Review of Selected Infrared Flow Visualization Processing Techniques: Contrast Enhancement and Frequency Domain Analysis

The deployment and integration of high-sensitivity infrared cameras in a transonic wind tunnel test environment has resulted in a unique capability to image aerodynamic phenomena in real-time. Multi-camera infrared flow visualization data systems are now routinely utilized at the NASA Ames Unitary Plan Wind Tunnel. The small flow-induced temperature gradients on the surface of the wind tunnel test article coupled with the high bit-depth of the infrared camera sensor makes the processing of the image data critically important. An image processing routine must enhance features of interest with minimal artifacts. Additionally, the production wind tunnel test environment demands that these processed images are made available in a real-time, automatic fashion. Therefore, any image processing routine must be computationally economical and enhance the image data with minimal input from a human operator. The following seeks to qualitatively explore selected image processing techniques by assessing their effectiveness to resolve flow features on a wind tunnel test article. A multi-scale contrast enhancement technique is introduced as well as a new implementation of a multi-scale, non-interpolated adaptive histogram equalization. Finally, a novel method is introduced that demonstrates the ability to resolve flow features imaged on bare-steel test articles possessing low emissivity. This method merges frequency domain analysis with contrast enhancement and has the potential to extend the application of infrared flow-visualization within the wind tunnel test environment.

Infrared Imaging↗

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗

Assimilation of Sea Ice Thickness Derived from Cryosat-2 Along-Track Freeboard Measurements into the Met Office's Forecast Ocean Assimilation Model (FOAM)

The feasibility of assimilating sea ice thickness (SIT) observations derived from CryoSat-2 along-track measurements of sea ice freeboard is successfully demonstrated using a 3D-Var assimilation scheme, NEMOVAR, within the Met Office's global, coupled ocean–sea-ice model, Forecast Ocean Assimilation Model (FOAM). The CryoSat-2 Arctic freeboard measurements are produced by the Centre for Polar Observation and Modelling (CPOM) and are converted to SIT within FOAM using modelled snow depth. This is the first time along-track observations of SIT have been used in this way, with other centres assimilating gridded and temporally averaged observations. The assimilation leads to improvements in the SIT analysis and forecast fields generated by FOAM, particularly in the Canadian Arctic. Arctic-wide observation-minus-background assimilation statistics for 2015–2017 show improvements of 0.75 m mean difference and 0.41 m root-mean-square difference (RMSD) in the freeze-up period and 0.46 m mean difference and 0.33 m RMSD in the ice break-up period. Validation of the SIT analysis against independent springtime in situ SIT observations from NASA Operation IceBridge (OIB) shows improvement in the SIT analysis of 0.61 m mean difference (0.42 m RMSD) compared to a control without SIT assimilation. Similar improvements are seen in the FOAM 5 d SIT forecast. Validation of the SIT assimilation with independent Beaufort Gyre Exploration Project (BGEP) sea ice draft observations does not show an improvement, since the assimilated CryoSat-2 observations compare similarly to the model without assimilation in this region. Comparison with airborne electromagnetic induction (Air-EM) combined measurements of SIT and snow depth shows poorer results for the assimilation compared to the control, despite covering similar locations to the OIB and BGEP datasets. This may be evidence of sampling uncertainty in the matchups with the Air-EM validation dataset, owing to the limited number of observations available over the time period of interest. This may also be evidence of noise in the SIT analysis or uncertainties in the modelled snow depth, in the assimilated SIT observations, or in the data used for validation. The SIT analysis could be improved by upgrading the observation uncertainties used in the assimilation. Despite the lack of CryoSat-2 SIT observations available for assimilation over the summer due to the detrimental effect of melt ponds on retrievals, it is shown that the model is able to retain improvements to the SIT field throughout the summer months due to prior, wintertime SIT assimilation. This also results in regional improvements to the July modelled sea ice concentration (SIC) of 5 % RMSD in the European sector, due to slower melt of the thicker sea ice.

Emma K. Fiedler↗

Interoperable Map Services with Performance Tuning for Earth Science Data through API-Tiles and Dynamic API-Styles

NASA’s Goddard Earth Sciences Data and Information Services Center (GES DISC) provides access to a wide range of global climate data from various satellite missions and models. However, the visualization and analysis of these data can be challenging due to their large volume, complex structure, and diverse formats. This study presents the implementation of interoperable map services (API-Maps) with performance tuning using API-Tiles and dynamic API-Styles. API-Maps is a standard for defining and exposing map services through RESTful (representational state transfer) APIs (application programming interfaces). API-Tiles is a technique for generating and delivering map tiles on demand from any data source. API-Styles is a method for dynamically applying styles to map tiles based on user preferences or data attributes. The use of API-Tiles and dynamic API-Styles enhances the performance and scalability of the map services, allowing for smooth and interactive visualization of large datasets. Two types of Earth Science data sources from the NASA GES DISC are used in the experiment: regularly gridded data, such as Global Precipitation Measurement (GPM) precipitation data, and low processing level data, such as low-level data of atmospheric composite measurements from the TROPOspheric Monitoring Instrument (TROPOMI) mission. Re-gridding of swath data (low level data - e.g. Level 2) of atmospheric composites (e.g. TROPOMI products, such as nitrogen dioxide, ozone and aerosol optical depth) is applied to enable the Web-based, interoperable, tiled, and styled mapping (rendering) services of such data. The results demonstrate the effectiveness of the proposed approach in providing fast and efficient access to Earth science data through interoperable map services.

Geographic Information System↗

Applications of ArcticDEM for measuring volcanic dynamics, landslides, retrogressive thaw slumps, snowdrifts, and vegetation heights

Topographical changes are of fundamental interest to a wide range of Arctic science disciplines faced with the need to anticipate, monitor, and respond to the effects of climate change, including geohazard management, glaciology, hydrology, permafrost, and ecology. This study demonstrates several geomorphological, cryospheric, and biophysical applications of ArcticDEM – a large collection of publicly available, time-dependent digital elevation models (DEMs) of the Arctic. Our study illustrates ArcticDEM's applicability across different disciplines and five orders of magnitude of elevation derivatives, including measuring volcanic lava flows, ice cauldrons, post-failure landslides, retrogressive thaw slumps, snowdrifts, and tundra vegetation heights. We quantified surface elevation changes in different geological settings and conditions using the time series of ArcticDEM. Following the 2014–2015 Bárðarbunga eruption in Iceland, ArcticDEM analysis mapped the lava flow field, and revealed the post-eruptive ice flows and ice cauldron dynamics. The total dense-rock equivalent (DRE) volume of lava flows is estimated to be (1431 ± 2) million m 3 . Then, we present the aftermath of a landslide in Kinnikinnick, Alaska, yielding a total landslide volume of (400 ± 8) × 103 m 3 and a total area of 0.025 km 2 . ArcticDEM is further proven useful for studying retrogressive thaw slumps (RTS). The ArcticDEM-mapped RTS profile is validated by ICESat-2 and drone photogrammetry resulting in a standard deviation of 0.5 m. Volume estimates for lake-side and hillslope RTSs range between 40,000 ± 9000 m 3 and 1,160,000 ± 85,000 m 3 , highlighting applicability across a range of RTS magnitudes. A case study for mapping tundra snow demonstrates ArcticDEM's potential for identifying high-accumulation, late-lying snow areas. The approach proves effective in quantifying relative snow accumulation rather than absolute values (standard deviation of 0.25 m, bias of −0.41 m, and a correlation coefficient of 0.69 with snow depth estimated by unmanned aerial systems photogrammetry). Furthermore, ArcticDEM data show its feasibility for estimating tundra vegetation heights with a standard deviation of 0.3 m (no bias) and a correlation up to 0.8 compared to the light detection and ranging (LiDAR). The demonstrated capabilities of ArcticDEM will pave the way for the broad and pan-Arctic use of this new data source for many disciplines, especially when combined with other imagery products. The wide range of signals embedded in ArcticDEM underscores the potential challenges in deciphering signals in regions affected by various geological processes and environmental influences.

Chunli Dai↗

Summarizing Multiple Aspects of Triple Collocation Analysis in a Single Diagram

With the ongoing expansion of global observation networks, it is expected that we shall routinely analyze records of geophysical variables such as temperature from multiple collocated instruments. Validating datasets in this situation is not a trivial task because every observing system has its own bias and noise. Triple collocation is a general statistical framework to estimate the error characteristics in three or more observational-based datasets. In a triple colocation analysis, several metrics are routinely reported but traditional multiple-panel plots are not the most effective way to display information. A new formula of error variance is derived for connecting the key terms in the triple collocation theory. A diagram based on this formula is devised to facilitate triple collocation analysis of any data from observations, as illustrated using three aerosol optical depth datasets from the recent Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE). An observational-based skill score is also derived to evaluate the quality of three datasets by taking into account both error variance and correlation coefficient. Several applications are discussed and sample plotting routines are provided.

triple collocation↗

ROP model for PDC bits in Geothermal drilling

Geothermal energy is a renewable source of energy, where heat extraction is preferentially balanced with the reservoir's natural heat recharge rate. The objective of this paper is to present and validate a novel rate of penetration (ROP) model for drilling hard and abrasive formations including granite formations for polycrystalline diamond compact (PDC) bits. The ROP model was developed based on a derived relationship of a threshold weight on cutter (WOC) and its corresponding depth of cut (DOC) for a single cutter. Laboratory data was used to scale the derived single cutter relationship to a full-hole ROP model for PDC bits. The ROP model includes a non-linear correlation for Phase I (inefficient drilling due to low WOB values) and a linear Phase II (efficient drilling) ROP response to WOB. The ROP model was verified using measured drilling parameter data from Utah FORGE well# 58-32 and data from Chocolate Mountains well # 17-8 in Southern California. When compared to oil and gas well drilling, geothermal drilling in granitic formations can be more difficult and complicated due to rock hardness and high temperatures. PDC bits can increase ROP and optimize drilling for these types of hard formations. This paper provides novel insight into the ROP response of PDC bits to drilling operational parameters.

15 GEOTHERMAL ENERGY↗

Low Temperature Geothermal Play Fairway Analysis For The Appalachian Basin: Phase 1 Revised Report November 18, 2016

Geothermal energy is an attractive sustainable energy source. Yet project developers need confirmation of the resource base to warrant their time and financial resources. The Geothermal Play Fairway Analysis of the Appalachian Basin evaluated risk metrics that communicate the favorability of potential low-temperature geothermal energy resources in reservoirs more than 1000 m below the surface. This analysis is focused on the direct use of the heat, rather than on electricity production. Four risk factors of concern for direct-use geothermal plays in the Appalachian Basin portions of New York, Pennsylvania, and West Virginia are examined individually, and then in combination: 1) thermal resource quality, 2) natural reservoir quality, 3) induced seismicity, and 4) utilization opportunities. Uncertainty in the risk estimation is quantified. Based on these metrics, geothermal plays in the Appalachian Basin were identified as potentially viable for a variety of direct-use-heat applications. The methodologies developed in this project may be applied in other sedimentary basins as a foundation for low temperature (50-150 °C), direct use geothermal resource, risk, and uncertainty assessment. Three methods with which to combine the four risk factors were used. Among these, the averaging of the individual risk factors indicates the most favorable counties within the study area are the West Virginia counties of Monongalia, Harrison, Lewis (dubbed the Morgantown–Clarksburg play fairway), Putnam, and Kanawha (Charleston play fairway), the New York counties of Chemung and Steuben plus adjacent Bradford county in Pennsylvania (Corning–Ithaca play fairway), and the Pennsylvania counties of Mercer, Crawford, Erie, and Warren, and adjacent Chautauqua county in New York (together, the Meadville–Jamestown play fairway). These higher priority regions are surrounded by broader medium priority zones. Also worthy of additional exploration is a broad region near Pittsburgh Pennsylvania, for which the available geological data are insufficient to fully analyze the geological risks but yet the population is high. First, to assess the spatial variation in the depth to which one would need to drill to obtain geothermal temperatures that are useful to a future project, the project used bottom-hole temperature data from Appalachian Basin oil and gas exploration. These bottom hole temperature data are abundant but of low quality. Second, the project examined the potential for sufficient water flow rates through rocks to harvest heat from a geothermal well field, considering only natural reservoirs. This analysis provides a very incomplete picture of spatial variability of natural reservoirs because the oil and gas reservoir data lack key properties and are spatially biased toward those locations with profitable amounts of hydrocarbons in the rock pore spaces. Third, in light of the fact that earthquake activity has been induced in several states by subsurface work related to the oil and gas industry, this project examined the potential for similar activity in the Appalachian Basin. Acknowledging that data for such a task are insufficient, we utilized what was available: records of seismic activity, regional estimates of the orientations of stress in the rocks, and locations and orientations of zones of lateral change in rock properties at depths down to several kilometers below Earth’s surface. With these data, we created a first approximation of spatially variable risks for induced earthquakes. Because no data existed with which to test the reliability of these methods, the results have a high degree of uncertainty. Fourth, we examined the spatial variability of the above-the-ground factors that contribute to the economical viability of projects to tap low-temperature geothermal resources for direct-use. We worked principally with population density as a regionally known variable that would impact the cost of district heating. The resulting maps omit the costs of producing the hot water from the ground, because the below-ground costs are directly coupled to the thermal resource risk factor and natural reservoir risk factor – later analyses of those costs will be needed. The result of the district heating analysis is highly skewed: few census locations yielded a low estimated surface cost. The team also identified more than 165 prospects for high value direct-use geothermal energy opportunities such as industrial sites, university campuses, and federal facilities, among others. At the closure of this regional analysis, the most significant technical uncertainties are 1) reservoir distribution and capacities; 2) validity of thermal resource maps, and 3) the holistic estimation of Levelized Cost of Heat for favorable geological situations.

15 GEOTHERMAL ENERGY↗

Comparisons of Arctic In-Situ Snow and Ice Data with Airborne Passive Microwave Measurements

As part of the AMSR-E sea ice validation campaign in March 2003, aircraft flights over the Arctic sea ice were coordinated with ground measurements of snow and sea ice properties. The surface-based measurements were in the vicinity of Barrow, AK, and at a Navy ice camp located in the Beaufort Sea. The NASA P-3 aircraft was equipped with the NOAA ETL PSR microwave radiometer that has the same frequencies as the AMSR-E sensor. The goal was to validate the standard AMSR-E products ice temperature and snow depth on sea ice. Ground measurements are the only way to validate these parameters. The higher spatial resolution of the PSR instrument (between 30 and 500 m, depending on altitude) enables a better comparison between ground measurements and microwave data because of the expected smaller spatial variability. Maps of PSR data can then be used for further down-scaling to AMSR-E pixel areas. Initial results show a good qualitative agreement between the in-situ snow depths and the PSR data. Detailed studies are underway and latest results will be presented.

Markus, T.↗