Search NASA⌕ Search

SEARCH · Search NASA

Results for “geophysical methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

FY 2026 Midyear Report: Seismic Monitoring of Underground Vibration Sources Using Distributed Acoustic Sensing and Seismometers

Safeguards-relevant temporal changes in underground facilities can be observed using geophysical monitoring techniques. Seismic waves, in particular, provide valuable insights into subsurface activities and can serve as an important tool for detecting anomalous events that may indicate containment breaches at geological repositories. This midyear report summarizes ongoing efforts to automatically and rapidly detect and locate anomalous vibration signals that could be indicative of potential containment breaches. Previous work during FY25 focused on compiling continuous seismic datasets from two underground sites and developing a database of continuous waveforms and ground-truth event data derived from multiple sensing modalities. Building on this foundation, we are adapting anomaly detection and geolocation algorithms to explore methods for monitoring underground activities using two relatively low-maintenance sensing technologies: a dense surface geophone array deployed at the Pleasant Gap mine in Pennsylvania, and a three-dimensional fiber-optic cable array for distributed acoustic sensing (DAS) installed in the subsurface at the Sanford Underground Research Facility (SURF) in South Dakota. This report summarizes work conducted during the first two quarters of FY26, during which we refined a dynamic power spectral density (PSD)-based detector, applied it independently to each geophone station, and then combined the per‑station detections with density-based spatial clustering of applications with noise (DBSCAN) to cluster events and produce spatial maps over a nine‑day interval. In addition, we outline plans for a field trial at the Waste Isolation Pilot Plant (WIPP) in New Mexico to compare traditional seismic monitoring approaches with DAS techniques and to evaluate the benefits of combined data analysis. Activities during the past two quarters have included the preparation and submission of a Field Test Plan to WIPP for approval, as well as submission to headquarters for review and feedback.

58 GEOSCIENCES↗
NASA NTRS1979

Report of the mid- and low-latitude E and F Region Working Group

The formation and variation of the ionosphere is addressed with regard to the ability to understand, specify, and predict the low and mid latitude E and F regions. A brief survey of prediction methods and techniques is given for long and short term variability in the E and F regions. It is indicated that the usefulness of theoretical models to predict electron density distribution in the low and mid latitude ionospheric E and F regions is limited by the ability to predict the parameters which enter the relevant equations; i.e., neutral atmospheric constituents, neutral and charged particle temperatures, neutral wind, electric fields, and ionizing sources such as solar (E sub uv) radiation and energetic particles. It is recommended that areas for research include improving knowledge of the input parameters and how they respond to changing solar and geophysical conditions.

Rush, C. M.↗
NASA NTRS1988

Kinetics of thermochemical gas-solid reactions important in the Venus sulfur cycle

The thermochemical net reaction CaCO3 + SO2 yields CaSO4 + CO is predicted to be an important sink for incorporation of SO2 into the Venus crust. The reaction rate law was established to understand the dependence of rate on experimental variables such as temperature and partial pressure of SO2, CO2, and O2. The experimental approach was a variant of the thermogravimetric method often employed to study the kinetics of thermochemical gas-solid reactions. Clear calcite crystals were heated at constant temperature in SO2-bearing gas streams for varying time periods. Reaction rate was determined by three independent methods. A weighted linear least squares fit to all rate data yielded a rate equation. Based on the Venera 13, 14 and Vega 2 observations of CaO content of the Venus atmosphere, SO2 at the calculated rate would be removed from the Venus atmosphere in about 1,900,00 years. The most plausible endogenic source of the sulfur needed to replenish atmospheric SO2 is volcanism. The annual amount of erupted material needed for the replenishment depends on sulfur content; three ratios are used to calculate rates ranging from 0.4 to 11 cu km/year. This geochemically derived volcanism rate can be used to test if geophysically derived rates are correct. The work also suggests that Venus is less volcanically active than the Earth.

Fegley, Bruce, Jr.↗
NASA NTRS1993

The effect of dissipative inhomogeneous medium on the statistics of the wave intensity

One of the main theoretical points in the theory of wave propagation in random medium is the derivation of closed form equations to describe the statistics of the propagating waves. In particular, in one dimensional problems, the closed form representation of the multiple scattering effects is important since it contributes in understanding such problems like wave localization, backscattering enhancement, and intensity fluctuations. In this the propagation of plane waves in a layer of one-dimensional dissipative random medium is considered. The medium is modeled by a complex permittivity whose real part is a constant representing the absorption. The one dimensional problem is mathematically equivalent to the analysis of a transmission line with randomly perturbed distributed parameters and a single mode lossy waveguide and the results can be used to study the propagation of radio waves through atmosphere and the remote sensing of geophysical media. It is assumed the scattering medium consists of an ensemble of one-dimensional point scatterers randomly positioned in a layer of thickness L with diffuse boundaries. A Poisson impulse process with density lambda is used to model the position of scatterers in the medium. By employing the Markov properties of this process an exact closed form equation of Kolmogorov-Feller type was obtained for the probability density of the reflection coefficient. This equation was solved by combining two limiting cases: (1) when the density of scatterers is small; and (2) when the medium is weakly dissipative. A two variable perturbation method for small lambda was used to obtain solutions valid for thick layers. These solutions are then asymptotically evaluated for small dissipation. To show the effect of dissipation, the mean and fluctuations of the reflected power are obtained. The results were compared with a lossy homogeneous medium and with a lossless inhomogeneous medium and the regions where the effect of absorption is not essential were discussed.

Saatchi, Sasan S.↗
NASA NTRS1997

Regularization for Atmospheric Temperature Retrieval Problems

Passive remote sensing of the atmosphere is used to determine the atmospheric state. A radiometer measures microwave emissions from earth's atmosphere and surface. The radiance measured by the radiometer is proportional to the brightness temperature. This brightness temperature can be used to estimate atmospheric parameters such as temperature and water vapor content. These quantities are of primary importance for different applications in meteorology, oceanography, and geophysical sciences. Depending on the range in the electromagnetic spectrum being measured by the radiometer and the atmospheric quantities to be estimated, the retrieval or inverse problem of determining atmospheric parameters from brightness temperature might be linear or nonlinear. In most applications, the retrieval problem requires the inversion of a Fredholm integral equation of the first kind making this an ill-posed problem. The numerical solution of the retrieval problem requires the transformation of the continuous problem into a discrete problem. The ill-posedness of the continuous problem translates into ill-conditioning or ill-posedness of the discrete problem. Regularization methods are used to convert the ill-posed problem into a well-posed one. In this paper, we present some results of our work in applying different regularization techniques to atmospheric temperature retrievals using brightness temperatures measured with the SSM/T-1 sensor. Simulation results are presented which show the potential of these techniques to improve temperature retrievals. In particular, no statistical assumptions are needed and the algorithms were capable of correctly estimating the temperature profile corner at the tropopause independent of the initial guess.

Velez-Reyes, Miguel↗
DOE OSTI2024

Perfectly Matched Layers and Characteristic Boundaries in Lattice Boltzmann: Accuracy vs Cost

Artificial boundary conditions (BCs) play a ubiquitous role in numerical simulations of transport phenomena in several diverse fields, such as fluid dynamics, electromagnetism, acoustics, geophysics, and many more. They are essential for accurately capturing the behavior of physical systems whenever the simulation domain is truncated for computational efficiency purposes. Ideally, an artificial BC would allow relevant information to enter or leave the computational domain without introducing artifacts or unphysical effects. Boundary conditions designed to control spurious wave reflections are referred to as nonreflective boundary conditions (NRBCs). Another approach is given by the perfectly matched layers (PMLs), in which the computational domain is extended with multiple dampening layers, where outgoing waves are absorbed exponentially in time. Here, in this work, the definition of PML is revised in the context of the lattice Boltzmann method. The impact of adopting different types of BCs at the edge of the dampening zone is evaluated and compared, in terms of both accuracy and computational costs. It is shown that for sufficiently large buffer zones, PMLs allow stable and accurate simulations even when using a simple zeroth-order extrapolation BC. Moreover, employing PMLs in combination with NRBCs potentially offers significant gains in accuracy at a modest computational overhead, provided the parameters of the BC are properly tuned to match the properties of the underlying fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

EGS Collab Experiment 2: Microseismic Monitoring

This dataset contains continuous seismic waveform data recorded during stimulation and thermal circulation tests for the Enhanced Geothermal Systems (EGS) Collab Experiment #2, conducted from February to September 2022 at the Sanford Underground Research Facility in Lead, South Dakota. This experiment aimed to study and validate models of geothermal systems by injecting high-pressure fluids into rock formations 1200-1500 meters below the surface, inducing microseismic events. The seismic monitoring system included 16 three-component accelerometers and a 24-channel hydrophone array, installed in boreholes surrounding the test area. Data were recorded at high sampling rates using a continuous waveform recording system to monitor seismic activity in real time. The dataset contains the raw data stored in binary format, with files named based on timestamps, and includes calibration certificates for some sensors to facilitate corrections to real units. Users are strongly advised to consult the accompanying detailed report, which outlines the experimental setup, sensor specifications, installation procedures, and data processing methods. The report also describes important nuances, such as the hardware filters on hydrophones, sensor calibration details, and the naming conventions for the recorded data. Proper use of this dataset may require familiarity with seismic data analysis tools, such as the Obspy Python package, and an understanding of the SEED naming conventions used for channel identification.

15 GEOTHERMAL ENERGY↗
NASA NTRS1977

Evaluation of algorithms for geological thermal-inertia mapping

The errors incurred in producing a thermal inertia map are of three general types: measurement, analysis, and model simplification. To emphasize the geophysical relevance of these errors, they were expressed in terms of uncertainty in thermal inertia and compared with the thermal inertia values of geologic materials. Thus the applications and practical limitations of the technique were illustrated. All errors were calculated using the parameter values appropriate to a site at the Raft River, Id. Although these error values serve to illustrate the magnitudes that can be expected from the three general types of errors, extrapolation to other sites should be done using parameter values particular to the area. Three surface temperature algorithms were evaluated: linear Fourier series, finite difference, and Laplace transform. In terms of resulting errors in thermal inertia, the Laplace transform method is the most accurate (260 TIU), the forward finite difference method is intermediate (300 TIU), and the linear Fourier series method the least accurate (460 TIU).

Miller, S. H.↗
NASA NTRS1994

Improvement of the photometric sunspot index and changes of the disk-integrated sunspot contrast with time

The photometric sunspot index (PSI) was developed to study the effects of sunspots on solar irradiance. It is calculated from the sunspot data published in the Solar-Geophysical Data catalog. It has been shown that the former PSI models overestimate the effect of dark sunspots on solar irradiance; furthermore results of direct sunspot photometry indicate that the contrast of spots depends on their area. An improved PSI calculation is presented; it takes into account the area dependence of the contrast and calculates `true' daily means for each observation using the differential rotation of the spots. Moreover, the observations are screened for outliers which improves the homogeneity of the data set substantially, at least for the period after December 1981 when NOAA started to report data from a few instead of one to two stations. A detailed description of the method is provided. The correlation between the newly calculated PSI and total solar irradiance is studied for different phases of the solar cycles 21 and 22 using bi-variate spectral analysis. The results can be used as a `calibration' of PSI in terms of gain, the factor by which PSI has to be multiplied to yield the observed irradiance change. The factor changes with time from about 0.6 in 1980 to 1.1 in 1990. This unexpected result cannot be interpreted by a change of the contrast relative to the quiet Sun (as it is normally defined and determined by direct photometry) but rather as a change of the contrast between the spots and their surrounding as seen in total irradiance (integrated over the solar disk). This may partly be explained by a change in the ratio between the areas of the spots and the surrounding faculae.

Froehlich, Claus↗
DOE OSTI2025

Basin-Scale Structural Features Database

The Basin-Scale Structural Features database provides spatial datasets of faults, fractures, folds, and earthquakes compiled from public, authoritative sources (e.g., U.S. Geological Survey and State Geological Surveys) and aggregated into derivative forms to support subsurface assessments. Recognizing that characterizing basin-scale structural features requires interpreting data that are often ambiguous or lack key information, the source data were evaluated using a knowledge-data framework and geospatial fuzzy logic method (Justman et al., 2020) to represent both measured (observed) and predicted (inferred or potential) structural features as derivative datasets. This workflow employs conceptual models for known structural features and predicted structural features, incorporating geospatial data to estimate potential, even with limited data. The aim is to aid and support an understanding of basin-scale features and identify potential gaps in data and knowledge. As of 4/30/2025, the database includes resources for nine sedimentary basins: Appalachian, Denver, U.S. Gulf Coast, Illinois, Michigan, Permian, Sacramento, San Joquin and Williston. The database is organized by basin and then data category: 1) Faults, fractures, folds, 2) Earthquakes, 3) Topographic, 4) Structural contours and isopachs, 5) Geophysical, and 6) Structural feature density assessment maps.

basin scale↗
NASA NTRS2008

Comment on "Weekly Precipitation Cycles? Lack of Evidence from United States Surface Stations"

There is a good deal of interest lately in whether or not rainfall varies with the day of the week in response to the weekly variations in human activity. The most likely cause of such changes in the U.S. would be from the weekly variations in pollution levels that are known to occur throughout the country. A paper on this topic will soon be published by the Journal of Geophysical Research entitled, "Midweek Increase in U.S. Summer Rain and Storm Heights Suggests Air Pollution Invigorates Rainstorms, by T. L. Bell, D. Rosenfeld, K.-M. Kim, J.-M. Yoo, M.-I. Lee, and M. Hahnenberger (referred to here as "Bell et al."). A paper by D. M. Schultz and co-authors was recently published in Geophysical Research Letters that claimed to contradict some of the results in Bell et al. The paper can be found here: . Our Comment points out that Schultz and co-authors ignored the fact that the results from satellite data obtained by Bell et al. were for a later time period than Schultz et al. examined, and that Bell et al. in fact also analyzed rainfall data for the same time period as Schultz et al. and, like them, also failed to find signs of a weekly cycle in rainfall during this time period. The contradictions claimed by Schultz et al. are non-existent. We point out some other problems with the methods and presentation by Schultz et al.

Bell, Thomas L.↗
NASA NTRS1997

Geophysical Characterization of Some Terranes and the Geophysical Modeling of Candidate Suture Zones

Indian participation in this project was terminated during the last year by a sudden withdrawal of support by the Department of Science and Technology, India, to the Indian Institute of Geomagnetism, Bombay. As a result, significant changes in the project focus had to be undertaken. Much of the work carried out at Southern Illinois University at Carbondale during the first year of the project anticipated the Indian participation and included development of computer programs to be used on gravity and magnetic data from the Indian subcontinent and preparations for fieldwork, tutorials, and workshops in India. Despite these setbacks, which were beyond our control, a number of significant tasks have been accomplished during the project period. These include: (1) Completion of digitization of the regional Bouguer gravity anomaly map of India and the regional ground total intensity magnetic anomaly map of India at an overdetermined spacing of 0.05 degrees. (2) We investigated and assessed the limitations of the Euler method using environmental examples because detailed aeromagnetic maps of parts of India were not available for interpretation by this method. (3) We also undertook an assessment of a suture zone between the Nyaza Craton (Archean) and the Mozambique Belt (Pan African) in the Kenya Rift, Africa, using gravity anomalies and the lithospheric seismological models. (4) We studied Magsat and high-altitude (approx. 4 km) aeromagnetic data over Canada.

Ravat, D.↗
NASA NTRS1994

Modeling of the coupled magnetospheric and neutral wind dynamos

This report summarizes the progress made in the first year of NASA Grant No. NAGW-3508 entitled 'Modeling of the Coupled Magnetospheric and Neutral Wind Dynamos.' The approach taken has been to impose magnetospheric boundary conditions with either pure voltage or current characteristics and solve the neutral wind dynamo equation under these conditions. The imposed boundary conditions determine whether the neutral wind dynamo will contribute to the high-latitude current system or the electric potential. The semi-annual technical report, dated December 15, 1993, provides further detail describing the scientific and numerical approach of the project. The numerical development has progressed and the dynamo solution for the case when the magnetosphere acts as a voltage source has been evaluated completely using spectral techniques. The simulation provides the field-aligned current distribution at high latitudes due to the neutral wind dynamo. A number of geophysical conditions can be simulated to evaluate the importance of the neutral wind dynamo contribution to the field-aligned current system. On average, field-aligned currents generated by the neutral wind dynamo contributed as much as 30 percent to the large-scale field-aligned current system driven by the magnetosphere. A term analysis of the high-latitude neutral wind dynamo equation describing the field aligned current distribution has also been developed to illustrate the important contributing factors involved in the process. The case describing the neutral dynamo response for a magnetosphere acting as a pure current generator requires the existing spectral code to be extended to a pseudo-spectral method and is currently under development.

Thayer, Jeffrey P.↗
DOE OSTI2026

Passive Microwave Radiometry and Active Radar Sounding as Complementary Tools for Geophysical Investigations of Icy Ocean Worlds

Juno Microwave Radiometer (MWR) observations of Europa and Ganymede offer critical insights into the icy shells of these moons ahead of NASA's Europa Clipper and ESA's JUpiter ICy moons Explorer (JUICE) missions. Both missions are equipped with active radar sounders designed to address key unknowns such as ice shell thickness, thermal state, and composition. In this study, we explore how passive microwave radiometry and active radar sounding can constrain ice shell properties, focusing on Europa. Using modeled microwave brightness temperature observations at 0.6 and 1.2 GHz alongside simulated radar attenuation rate observations, we show that each instrument can independently produce robust ice shell thickness constraints under idealized conditions. We then relax these assumptions, quantifying how uncertainties from non-ideal properties—including convective layers, freezing-point depression, and chloride-doped ice—affect thickness estimates. Finally, we demonstrate how combining observations from these complementary techniques breaks degeneracies between ice shell properties, enabling more robust constraints than either method alone. This approach will maximize the science return of Europa Clipper and JUICE, advancing our understanding of the thermophysical structure and habitability of icy ocean worlds.

58 GEOSCIENCES↗
NASA NTRS1992

Phased array feed design technology for Large Aperture Microwave Radiometer (LAMR) Earth observations

An assessment of the potential and limitations of phased array antennas in space-based geophysical precision radiometry is described. Mathematical models exhibiting the dependence of system and scene temperatures and system sensitivity on phased array antenna parameters and components such as phase shifters and low noise amplifiers (LNA) are developed. Emphasis is given to minimum noise temperature designs wherein the LNA's are located at the array level, one per element or subarray. Two types of combiners are considered: array lenses (space feeds) and corporate networks. The result of a survey of suitable components and devices is described. The data obtained from that survey are used in conjunction with the mathematical models to yield an assessment of effective array antenna noise temperature for representative geostationary and low Earth orbit systems. Practical methods of calibrating a space-based, phased array radiometer are briefly addressed as well.

Schuman, H. K.↗
DOE OSTI2024

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2022

A Waveform Skewness Index for Measuring Time Series Nonlinearity and Its Applications to the ENSO-Indian Monsoon Relationship

Many geophysical time series possess nonlinear characteristics that reflect the underlying physics of the phenomena the time series describe. The nonlinear character of times series can change with time, so it is important to quantify time series nonlinearity without assuming stationarity. A common way of quantifying the time evolution of time series nonlinearity is to compute sliding skewness time series, but it is shown here that such an approach can be misleading when time series contain periodicities. To remedy this deficiency of skewness, a new waveform skewness index is proposed for quantifying local nonlinearities embedded in time series. A waveform skewness spectrum is proposed for determining the frequency components that are contributing to time series waveform skewness. The new methods are applied to the El Niño–Southern Oscillation (ENSO) and the Indian monsoon to test a recently proposed hypothesis that states that changes in the ENSO–Indian monsoon relationship are related to ENSO nonlinearity. We show that the ENSO–Indian rainfall relationship weakens during time periods of high ENSO waveform skewness. The results from two different analyses suggest that the breakdown of the ENSO–Indian monsoon relationship during time periods of high ENSO waveform skewness is related to the more frequent occurrence of strong central Pacific El Niño events, supporting arguments that changes in the ENSO–Indian rainfall relationship are not solely related to noise.

Justin Schulte↗
NASA NTRSDate not supplied

New Ways of Facilitating Improved Data Discovery and Access for NASA's Suborbital Earth Science Observations

NASA conducts field research in various Earth Science disciplines utilizing airborne and other non-satellite platforms to acquire in situ and remotely sensed observations indicative of physical processes across a range of scales. Field efforts are key in the development and validation of instruments and satellite algorithm refinements. The heterogeneous data, with a range of file formats, scales, and acquisition methods, support research in several science areas. NASA’s archive process assigns data products to discipline-oriented Distributed Active Archive Centers (DAACs) for stewardship. Over time, individual DAACs have developed tools for data browsing and serving disparate user bases. As science becomes more interdisciplinary, researchers need to incorporate observations from multiple campaigns, and multiple DAACs, into their work. Motivated in part by this shifting paradigm of needs, the Catalog of Archived Suborbital Earth Science Investigations (CASEI) was created. CASEI provides a single starting point to browse, search, and discover airborne and field data. Contextual metadata are organized and inter-linked allowing intuitive, integrated exploration across all NASA DAACs. Campaign science objectives, platform and instrument configurations, geographical details, geophysical concepts, and more are tracked in CASEI’s database, facilitating multi-parameter search, browse, and discovery of relevant data products. Researchers are able to directly access associated data products, via DOI links, regardless of the DAAC where they reside. Significant events, key time periods of high science interest within the longer-duration campaign effort, are also indicated and allow for a more efficient identification of critical data subsets. This presentation describes CASEI’s development, intensive metadata curation process, and demonstrates the web interface experience. Initial content metrics and plans for continued maintenance will also be discussed.

metadata↗