Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geophysical”

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 163 records · Page 9

Temporal and spatial variability in surface roughness and accumulation rate around 88° S from repeat airborne geophysical surveys

We use repeat high-resolution airborne geophysical data consisting of laser altimetry, snow, and Ku-band radar and optical imagery acquired in 2014, 2016, and 2017 to analyze the spatial and temporal variability in surface roughness, slope, wind deposition, and snow accumulation at 88° S, an elevation bias validation site for ICESat-2 and potential validation site for CryoSat-2. We find significant small-scale variability (<10 km) in snow accumulation based on the snow radar subsurface stratigraphy, indicating areas of strong wind redistribution are prevalent at 88° S. In general, highs in snow accumulation rate correspond with topographic lows, resulting in a negative correlation coefficient of r(exp 2)=−0.32 between accumulation rate and MSWD (mean slope in the mean wind direction). This relationship is strongest in areas where the dominant wind direction is parallel to the survey profile, which is expected as the geophysical surveys only capture a two-dimensional cross section of snow redistribution. Variability in snow accumulation appears to correlate with variability in MSWD. The correlation coefficient between the standard deviations of accumulation rate and MSWD is r(exp 2)=0.48, indicating a stronger link between the standard deviations than the actual parameters. Our analysis shows that there is no simple relationship between surface slope, wind direction, and snow accumulation rates for the overall survey area. We find high variability in surface roughness derived from laser altimetry measurements on length scales smaller than 10 km, sometimes with very distinct and sharp transitions. Some areas also show significant temporal variability over the course of the 3 survey years. Ultimately, there is no statistically significant slope-independent relationship between surface roughness and accumulation rates within our survey area. The observed correspondence between the small-scale temporal and spatial variability in surface roughness and backscatter, as evidenced by Ku-band radar signal strength retrievals, will make it difficult to develop elevation bias corrections for radar altimeter retrieval algorithms.

Michael Studinger↗

Geophysical Retrievals and Cloud Analyses from Merged Airborne Radiometer Datasets Covering 10–684 GHz

Airborne microwave radiometers provide insight about numerous aspects of Earth’s atmosphere and yield critical validation datasets for spaceborne radiometers. Three radiometers that are important to NASA’s airborne remote-sensing arsenal include: the Advanced Microwave Precipitation Radiometer (AMPR), covering 10–85 GHz; the Conical Scanning Millimeter-wave Imaging Radiometer (CoSMIR), covering 50–183 GHz; and the Compact Scanning Submillimeter-wave Imaging Radiometer (CoSSIR), covering 170–684 GHz. The NASA field campaigns of interest to this study include: the Integrated Precipitation and Hydrology Experiment (IPHEx) in 2014, the Olympic Mountains Experiment and Radar Definition Experiment (OLYMPEX/RADEX) in 2015–2016, the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) in 2020–2023, and the Airborne Lightning Observatory for FEGS and TGFs (ALOFT) in 2023. To provide a more comprehensive perspective on clouds and precipitation observed during these airborne field campaigns, AMPR data were merged spatiotemporally with CoSMIR data for IPHEx, OLYMPEX/RADEX, and IMPACTS (2020 and 2022), while considering differences in instrument characteristics and operations, providing brightness temperature (Tb) values from 10–183 GHz in a common background grid throughout each flight. AMPR and CoSSIR data were similarly merged for IMPACTS (2023) and ALOFT, providing a common background grid with Tb values covering 10–684 GHz throughout each flight. These merged Tb data were employed in geophysical retrievals using the Community Radiative Transfer Model (CRTM), an Eddington radiative transfer model, and a one-dimensional variational (1DVAR) inversion method. Retrievals of cloud liquid water path were of primary interest. This presentation will include an overview of the methods for the radiometer data mergers, the radiative transfer methods, the geophysical retrievals, and detailed results from examining trends in Tb and cloud liquid water path in clouds, precipitation, and cloud-to-precipitation transition zones.

Corey G Amiot↗

NEAR-SURFACE GEOPHYSICAL CHARACTERIZATION OF A HOLOCENE FAULT CONDUCIVE TO GEOTHERMAL FLOW NEAR PYRAMID LAKE, NEVADA

Linear deposits of calcium carbonate tufa columns mark recent faults that cut 11 ka Lake Lahontan sediments at Astor Pass, north of Pyramid Lake, Nevada. Throughout the Great Basin, faults appear to control the location of geothermal resources by providing pathways for fluid migration. Reservoir-depth (greater than 1 km) seismic imaging at Astor Pass reveals a fault that projects to one of the lines of tufa columns at the surface. The presence of the tufa deposits suggests this fault carried warm geothermal waters through the lakebed clay sediments in recent time. The warm fluids deposited the tufa when they hit cold Lake Lahontan water at the lakebed. Lake Lahontan covered this location 11 ka to a depth of at least 60 m. In collaboration with the Pyramid Lake Paiute Tribe, an Applied Geophysics class at UNR investigated the near-surface geophysical characteristics of this fault. The survey at and near the tufa columns comprises near-surface Pwave seismic reflection and refraction, electrical resistivity tomography, nearsurface refraction microtremor arrays, nine near-surface direct-current resistivity soundings, magnetic surveys, and gravity surveys. The refraction microtremor results show shear velocities near tufa and faults to be marginally lower, compared to Vs away from the faults. Overall, the 30-m depth-averaged shear velocities are low, less than 300 m/s, consistent with the lakebed clay deposits. These results indicate that no seismically fast (> 500 m/s) tufa deposits are present below the surface at or near the tufa columns. Vs30 averages were for example 274 ± 13 m/s on the fault, 287 ± 2 m/s at 150 m east of the fault, and 290 ± 15 m/s at 150 m west of the fault. The P-velocity refraction optimization results similarly indicate a lack of high-velocity tufa buried below the surface in the Lahontan sediments, reinforcing the idea that all tufa was deposited above the lakebed surface. The seismic results provide a negative test of the hypothesis that deposition of the lakebeds in the Quaternary buried and preserved older tufa columns within the section. Near-surface Wenner arrays with a-spacings up to 30 m show a higher resistivity near the faults, and tufa, than away from the faults. Resistivity averages within a few meters of the surface were 33 ± 17 ohm-m on the fault, 13 ± 3 ohm-m east of the fault, and 9 ± 3 ohm-m west of the fault. It is possible that the geothermal waters are fresher, and more resistive, than waters held in the lakebed clays. Water samples from more than 1 km depth in exploration wells have a TDS of 2500 p.p.m., nearly drinking-water quality. The relatively resistive water, perhaps localized by greater permeability along the fault, could explain the higher resistivity measured near the fault. The results show that there is no high-velocity, high-resistivity tufa along the faults below the surface, so we are unable to use buried tufa to locate the faults that may promote geothermal upwelling in this area. We further hypothesize that as sedimentation buried the tufa during the Quaternary, warm geothermal waters re-dissolved it, and re-precipitated it in the cold lake-bottom water.

15 GEOTHERMAL ENERGY↗

Fracture Propagation, Fluid Flow, and Geomechanics of Water-Based Hydraulic Fracturing in Shale Gas Systems and Electromagnetic Geophysical Monitoring of Fluid Migration

We investigate fracture propagation induced by hydraulic fracturing with water injection, using numerical simulation. For rigorous, full 3D modeling, we employ a numerical method that can model failure resulting from tensile and shear stresses, dynamic nonlinear permeability, leak-off in all directions, and thermo-poro-mechanical effects with the double porosity approach. Our numerical results indicate that fracture propagation is not the same as propagation of the water front, because fracturing is governed by geomechanics, whereas water saturation is determined by fluid flow. At early times, the water saturation front is almost identical to the fracture tip, suggesting that the fracture is mostly filled with injected water. However, at late times, advance of the water front is retarded compared to fracture propagation, yielding a significant gap between the water front and the fracture top, which is filled with reservoir gas. We also find considerable leak-off of water to the reservoir. The inconsistency between the fracture volume and the volume of injected water cannot properly calculate the fracture length, when it is estimated based on the simple assumption that the fracture is fully saturated with injected water. As an example of flow-geomechanical responses, we identify pressure fluctuation under constant water injection, because hydraulic fracturing is itself a set of many failure processes, in which pressure consistently drops when failure occurs, but fluctuation decreases as the fracture length grows. We also study application of electromagnetic (EM) geophysical methods, because these methods are highly sensitive to changes in porosity and pore-fluid properties due to water injection into gas reservoirs. Employing a 3D finite-element EM geophysical simulator, we evaluate the sensitivity of the crosswell EM method for monitoring fluid movements in shaly reservoirs. For this sensitivity evaluation, reservoir models are generated through the coupled flow-geomechanical simulator and are transformed via a rock-physics model into electrical conductivity models. It is shown that anomalous conductivity distribution in the resulting models is closely related to injected water saturation, but not closely related to newly created unsaturated fractures. Our numerical modeling experiments demonstrate that the crosswell EM method can be highly sensitive to conductivity changes that directly indicate the migration pathways of the injected fluid. Accordingly, the EM method can serve as an effective monitoring tool for distribution of injected fluids (i.e., migration pathways) during hydraulic fracturing operations

02 PETROLEUM↗

Hiereachical Bayesian Model for Combining Geochemical and Geophysical Data for Environmental Applications Software

Development of a hierarchical Bayesian model to estimate the spatiotemporal distribution of aqueous geochemical parameters associated with in-situ bioremediation using surface spectral induced polarization (SIP) data and borehole geochemical measurements collected during a bioremediation experiment at a uranium-contaminated site near Rifle, Colorado. The SIP data are first inverted for Cole-Cole parameters including chargeability, time constant, resistivity at the DC frequency and dependence factor, at each pixel of two-dimensional grids using a previously developed stochastic method. Correlations between the inverted Cole-Cole parameters and the wellbore-based groundwater chemistry measurements indicative of key metabolic processes within the aquifer (e.g. ferrous iron, sulfate, uranium) were established and used as a basis for petrophysical model development. The developed Bayesian model consists of three levels of statistical sub-models: 1) data model, providing links between geochemical and geophysical attributes, 2) process model, describing the spatial and temporal variability of geochemical properties in the subsurface system, and 3) parameter model, describing prior distributions of various parameters and initial conditions. The unknown parameters are estimated using Markov chain Monte Carlo methods. By combining the temporally distributed geochemical data with the spatially distributed geophysical data, we obtain the spatio-temporal distribution of ferrous iron, sulfate and sulfide, and their associated uncertainity information. The obtained results can be used to assess the efficacy of the bioremediation treatment over space and time and to constrain reactive transport models.

Chen, Jinsong↗

Integrated Geologic and Geophysical Approach for Establishing Geothermal Play Fairways and Discovering Blind Geothermal Systems in the Great Basin Region, Western USA: A Progress Report

We have undertaken an integrated geologic, geochemical, and geophysical study of a broad 240-km-wide, 400-km-long transect stretching from west-central to eastern Nevada in the Great Basin region of the western USA. The main goal of this study is to produce a comprehensive geothermal potential map that incorporates up to 11 parameters and identifies geothermal play fairways that represent potential blind or hidden geothermal systems. Our new geothermal potential map incorporates: 1) heat flow; 2) geochemistry from springs and wells; 3) structural setting; 4) recency of faulting; 5) slip rates on Quaternary faults; 6) regional strain rate; 7) slip and dilation tendency on Quaternary faults; 8) seismologic data; 9) gravity data; 10) magnetotelluric data (where available); and 11) seismic reflection data (primarily from the Carson Sink and Steptoe basins). The transect is respectively anchored on its western and eastern ends by regional 3D modeling of the Carson Sink and Steptoe basins, which will provide more detailed geothermal potential maps of these two promising areas. To date, geological, geochemical, and geophysical data sets have been assembled into an ArcGIS platform and combined into a preliminary predictive geothermal play fairway model using various statistical techniques. The fairway model consists of the following components, each of which are represented in grid-cell format in ArcGIS and combined using specified weights and mathematical operators: 1) structural component of permeability; 2) regional-scale component of permeability; 3) combined permeability, and 4) heat source model. The preliminary model demonstrates that the multiple data sets can be successfully combined into a comprehensive favorability map. An initial evaluation using known geothermal systems as benchmarks to test interpretations indicates that the preliminary modeling has done a good job assigning relative ranks of geothermal potential. However, a major challenge is defining logical relative rankings of each parameter and how best to combine the multiple data sets into the geothermal potential/ permeability map. Ongoing feedback and data analysis are in use to revise the grouping and weighting of some parameters in order to develop a more robust, optimized, final model. The final product will incorporate more parameters into a geothermal potential map than any previous effort in the region and may serve as a prototype to develop comprehensive geothermal potential maps for other regions.

15 GEOTHERMAL ENERGY↗

Integrated Geologic and Geophysical Approach for Establishing Geothermal Play Fairways and Discovering Blind Geothermal Systems in the Great Basin Region, Western USA: A Progress Report

We have undertaken an integrated geologic, geochemical, and geophysical study of a broad 240-km-wide, 400-kmlong transect stretching from west-central to eastern Nevada in the Great Basin region of the western USA. The main goal of this study is to produce a comprehensive geothermal potential map that incorporates up to 11 parameters and identifies geothermal play fairways that represent potential blind or hidden geothermal systems. Our new geothermal potential map incorporates: 1) heat flow; 2) geochemistry from springs and wells; 3) structural setting; 4) recency of faulting; 5) slip rates on Quaternary faults; 6) regional strain rate; 7) slip and dilation tendency on Quaternary faults; 8) seismologic data; 9) gravity data; 10) magnetotelluric data (where available); and 11) seismic reflection data (primarily from the Carson Sink and Steptoe basins). The transect is respectively anchored on its western and eastern ends by regional 3D modeling of the Carson Sink and Steptoe basins, which will provide more detailed geothermal potential maps of these two promising areas. To date, geological, geochemical, and geophysical data sets have been assembled into an ArcGIS platform and combined into a preliminary predictive geothermal play fairway model using various statistical techniques. The fairway model consists of the following components, each of which are represented in grid-cell format in ArcGIS and combined using specified weights and mathematical operators: 1) structural component of permeability; 2) regional-scale component of permeability; 3) combined permeability, and 4) heat source model. The preliminary model demonstrates that the multiple data sets can be successfully combined into a comprehensive favorability map. An initial evaluation using known geothermal systems as benchmarks to test interpretations indicates that the preliminary modeling has done a good job assigning relative ranks of geothermal potential. However, a major challenge is defining logical relative rankings of each parameter and how best to combine the multiple data sets into the geothermal potential/permeability map. Ongoing feedback and data analysis are in use to revise the grouping and weighting of some parameters in order to develop a more robust, optimized, final model. The final product will incorporate more parameters into a geothermal potential map than any previous effort in the region and may serve as a prototype to develop comprehensive geothermal potential maps for other regions.

15 GEOTHERMAL ENERGY↗

Discovering Geothermal Systems in the Great Basin Region: An Integrated Geologic, Geochemical, and Geophysical Approach for Establishing Geothermal Play Fairways

The Great Basin currently hosts ~24 geothermal power plants with more than 600 MW of capacity, but estimates suggest that the region is capable of producing much greater amounts of geothermal energy. Recent studies indicate that a combination of geological, geochemical, and geophysical features characterize the most robust geothermal systems. Considering that most of the geothermal resources in this region are blind, it is imperative that the favorable characteristics for geothermal activity be synthesized and techniques perfected for the discovery of new viable systems. We therefore utilized a multi-disciplinary approach and synthesized 9 geologic, geochemical, and geophysical parameters to produce a new detailed geothermal play fairway map of a large transect across the Great Basin region (96,000 km2) from west-central to eastern Nevada, with the primary objective of facilitating discovery of commercial-grade, blind geothermal fields (i.e. systems with no surface hot springs or fumaroles). Our analysis included: 1) structural settings (i.e., patterns of faulting), 2) age of recent faulting, 3) slip rates on recent faults, 4) regional-scale strain rates, 5) the tendency of faults to slip or dilate based on their orientation in the regional stress field, 6) earthquake density, 7) gravity data, 8) temperature at 3 km depth, and 9) geochemistry from springs and wells. These parameters were grouped into subsets to delineate rankings for local permeability, regional permeability, and heat, which collectively defined the geothermal play fairways (i.e., most likely locations for significant geothermal fluid flow). A major challenge was developing appropriate weightings of individual data types to best predict permeability and overall geothermal potential. Rigorous statistical methods, utilizing 34 benchmarks of known relatively high-temperature (>130°C) geothermal systems in the region, were employed to determine the hierarchal weights of each parameter. The new geothermal potential maps will help to reduce the risks in prospecting for new geothermal systems throughout the region.

15 GEOTHERMAL ENERGY↗

3D Geophysical Inversion Modeling of Ground Magnetic Data at Baker Hot Springs, Washington State, USA

Ground-based magnetic survey data were geophysically modelled in 3D as part of a play fairway exploration program at Baker hot springs. Here, a multi-step inversion method was used to address the issue of remanent magnetization in rocks of the study area. The geophysical modelling effort generated an effective susceptibility model that displays the 3D distribution and orientation of magnetization in the subsurface. Importantly, the magnetic modelling effort found that a strong, dipolar magnetic anomaly observed in the study area can be explained by a remanently magnetized body that is ~1¼ km long x ~200 m wide and extends from ~25-400 m below the surface. Geologically, this body may represent 1) an intrusion of reversely magnetized mafic rocks or 2) a magnetite-bearing fault zone. This rock body coincides spatially with Baker hot springs; the margins of this body or the fault structures within it are likely responsible for providing the permeability that gives rise to the hot springs.

Baker hot springs↗

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗

Geophysical monitoring using active seismic techniques at the Citronelle Alabama CO 2 storage demonstration site

Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.

58 GEOSCIENCES↗

An iterative particle filter approach for coupled hydro-geophysical inversion of a controlled infiltration experiment

The modeling of unsaturated groundwater flow is affected by a high degree of uncertainty related to both measurement and model errors. Geophysical methods such as Electrical Resistivity Tomography (ERT) can provide useful indirect information on the hydrological processes occurring in the vadose zone. In this paper, we propose and test an iterated particle filter method to solve the coupled hydrogeophysical inverse problem. We focus on an infiltration test monitored by time-lapse ERT and modeled using Richards equation. The goal is to identify hydrological model parameters from ERT electrical potential measurements. Traditional uncoupled inversion relies on the solution of two sequential inverse problems, the first one applied to the ERT measurements, the second one to Richards equation. This approach does not ensure an accurate quantitative description of the physical state, typically violating mass balance. To avoid one of these two inversions and incorporate in the process more physical simulation constraints, we cast the problem within the framework of a SIR (Sequential Importance Resampling) data assimilation approach that uses a Richards equation solver to model the hydrological dynamics and a forward ERT simulator combined with Archie's law to serve as measurement model. ERT observations are then used to update the state of the system as well as to estimate the model parameters and their posterior distribution. The limitations of the traditional sequential Bayesian approach are investigated and an innovative iterative approach is proposed to estimate the model parameters with high accuracy. The numerical properties of the developed algorithm are verified on both homogeneous and heterogeneous synthetic test cases based on a real-world field experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗