Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational seismology”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

A comparison of techniques for inversion of radio-ray phase data in presence of ray bending

Derivations are presented of the straight-line Abel transform and the seismological Herglotz-Wiechert transform (which takes ray bending into account) that are used in the reconstruction of refractivity profiles from radio-wave phase data. Profile inversion utilizing these approaches, performed in computer-simulated experiments, are compared for cases of positive, zero, and negative ray bending. For thin atmospheres and ionospheres, such as the Martian atmosphere and ionosphere, radio wave signals are shown to be inverted accurately with both methods. For dense media, such as the solar corona or the lower Venus atmosphere, the refractive recovered by the seismological Herglotz-Wiechert transform provide a significant improvement compared with the straight-line Abel transform.

Wallio, H. A.↗
DOE OSTI2026

Learning earthquake ground motions via conditional generative modeling

Predicting high-fidelity ground motions for future earthquakes is crucial for seismic hazard assessment and infrastructure resilience. Conventional empirical simulations suffer from sparse sensor distribution and geographically localized earthquake locations, while physics-based methods are computationally intensive and require accurate representations of Earth structures and earthquake sources. We propose an artificial intelligence (AI) spectrogram generator, Conditional Generative Modeling for Ground Motion (CGM-GM). CGM-GM leverages earthquake magnitudes and geographic coordinates of earthquakes and sensors as inputs, when postprocessed with phase information, capturing spatially continuous Fourier amplitude spectra (FAS) as well as properties such as P and S arrivals, and waveform durations, without explicit physics constraints. This is achieved through a probabilistic autoencoder that extracts latent distributions in the time-frequency domain and variational sequential models for prior and posterior distributions. We evaluate the performance of CGM-GM using small-magnitude earthquake records from the San Francisco Bay Area, a region with high seismic risks. Here, we report that CGM-GM demonstrates potential for complementing physics-based simulations and non-ergodic empirical ground motion models, as well as shows promise in seismology and beyond.

geophysics↗
DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗
NASA NTRS1995

The Topographic Torque on a Bounding Surface of a Rotating Gravitating Fluid and the Excitation by Core Motions of Decadal Fluctuations in the Earth's Rotation

General expressions (with potential applications in several areas of geophysical fluid dynamics) are derived for all three components of the contribution made by the geostrophic part of the pressure field associated with flow in a rotating gravitating fluid to the topographic torque exerted by the fluid on a rigid impermeable bounding surface of any shape. When applied to the Earth's liquid metallic core, which is bounded by nearly spherical surfaces and can be divided into two main regions, the "torosphere" and "polosphere," the expressions reduce to formulae given previously by the author, thereby providing further support for his work and that of others on the role of topographic coupling at the core-mantle boundary in the excitation by core motions of Earth rotation fluctuations on decadal time scales. They also show that recent criticisms of that work are vitiated by mathematical and physical errors. Contrary to these criticisms, the author's scheme for exploiting Earth rotation and other geophysical data (either real or simulated in computer models) in quantitative studies of the topography of the core-mantle boundary (CMB) by intercomparing various models of (a) motions in the core based on geomagnetic secular variation data and (b) CMB topography based on seismological and gravity data has a sound theoretical basis. The practical scope of the scheme is of course limited by the accuracy of real data, but this is a matter for investigation, not a priori assessment.

Hide, Raymond↗
DOE OSTI2024

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗
NASA NTRS1982

Astronomy on ice

A solar seismology program based on observations made at the South Pole is discussed. Environmental and observing conditions at the South Pole are reviewed, and the solar telescopes used in the program are described. The history of solar seismology and of observations of solar surface oscillations is briefly summarized. It is noted that the South Pole was chosen as the site for the program to obtain good frequency resolution and that the observations have produced more than 600 million bytes of computer data.

Harvey, J.↗
NASA NTRS1994

Prospects for Jovian seismological observations following the impact of comet Shoemaker-Levy 9

The impact of each fragment of comet SL-9 will produce a downward-propagating pressure wave which will travel at the sound speed through the jovian interior. Since the sound speed increases with depth, most of the energy in the pressure pulse will be strongly refracted and return to the surface, as recently computed by Marley (1994). This wave may in principle be observable as it propagates into the stratosphere, using sufficiently sensitive thermal infrared imaging. If so, it will provide a unique opportunity to constrain models of the jovian interior. This paper extends Marley's calculations to include the effect of the limited spatial resolution which will be characteristic of real observations. The wave pattern on the disk will consist of closely spaced regions of alternating temperature increases and decreases. Spatial averaging will significantly reduce the observed amplitude for resolutions attainable using earth-based telescopes, but the waves should remain above the detection limit.

Deming, Drake↗
NASA NTRS1975

A hybrid computer time code reader used to digitize sounding rocket data

A technique is described for utilizing a hybrid computer as a time code reader and digitizer. The NASA 36-bit time code is an absolute time reference for digitizing analog signals to 1 mV resolution at a 1 kHz sampling rate. The analog computer circuit and the operation of the merged analog and digital computer systems is described. A result of this digitization process applied to analog signals recorded on magnetic tape during an auroral sounding rocket flight is given as an example. This is a general technique that could be applied to other areas such as geology, meteorology, oceanography, and seismology.

Dukich, T.↗
NASA NTRS1992

Sources of uncertainty in direct seismological measurements of the solar helium abundance

The methods by which Dappen et al. (1988, 1990, 1991) and Dziembowski et al. (1990, 1991, 1992) recently obtained discrepant estimates of the helium abundance in the solar convection zone are compared. The aim of the investigation reported in this paper is to identify the main source of the discrepancy. Using as proxy data eigenfrequencies of a set of modes of a theoretical solar model, computed with the same physics as were the frequencies of a reference model with which these data are compared, the two methods yield similar results. Thus we ascertain that the principal source of the discrepancy is not in the inversions themselves, which yield essentially a measure of the variation of the adiabatic exponent gamma of the material in the He II ionization zone. Instead it is in the approximations adopted in the treatment of heavy elements in the equation of state used to relate the variation of gamma to chemical composition. We obtain acceptably consistent results when inverting solar data by two methods using the same equation of state. We attempt to identify the likely residual sources of uncertainty.

Kosovichev, A. G.↗
NASA NTRS2008

Wireless Sensor Network Handles Image Data

To relay data from remote locations for NASA s Earth sciences research, Goddard Space Flight Center contributed to the development of "microservers" (wireless sensor network nodes), which are now used commercially as a quick and affordable means to capture and distribute geographical information, including rich sets of aerial and street-level imagery. NASA began this work out of a necessity for real-time recovery of remote sensor data. These microservers work much like a wireless office network, relaying information between devices. The key difference, however, is that instead of linking workstations within one office, the interconnected microservers operate miles away from one another. This attribute traces back to the technology s original use: The microservers were originally designed for seismology on remote glaciers and ice streams in Alaska, Greenland, and Antarctica-acquiring, storing, and relaying data wirelessly between ground sensors. The microservers boast three key attributes. First, a researcher in the field can establish a "managed network" of microservers and rapidly see the data streams (recovered wirelessly) on a field computer. This rapid feedback permits the researcher to reconfigure the network for different purposes over the course of a field campaign. Second, through careful power management, the microservers can dwell unsupervised in the field for up to 2 years, collecting tremendous amounts of data at a research location. The third attribute is the exciting potential to deploy a microserver network that works in synchrony with robotic explorers (e.g., providing ground truth validation for satellites, supporting rovers as they traverse the local environment). Managed networks of remote microservers that relay data unsupervised for up to 2 years can drastically reduce the costs of field instrumentation and data rec

Source record↗
NASA NTRS1977

Solar seismology. II - The stochastic excitation of the solar p-modes by turbulent convection

We test the hypothesis that the solar p-modes are stabilized by damping due to turbulent viscosity in the convective zone. Starting from the assumption that the modes are stable, we calculate expectation values for the modal energies. We find that the interaction between a p-mode and the turbulent convection is such that the modal energy tends toward equipartition with the kinetic energy of turbulent eddies whose lifetimes are comparable to the modal period. From the calculated values of the modal energies, we compute rms surface velocity amplitudes. Our predicted rms surface velocities range from 0.01 cm/sec for the fundamental radial mode to 0.6 cm/sec for the radial mode whose period is approximately 5 minutes. The predicted surface velocities for the low order p-modes are much smaller than the velocities inferred from recent observations.

Goldreich, P.↗
NASA NTRS1992

Inference in infinite-dimensional inverse problems - Discretization and duality

Many techniques for solving inverse problems involve approximating the unknown model, a function, by a finite-dimensional 'discretization' or parametric representation. The uncertainty in the computed solution is sometimes taken to be the uncertainty within the parametrization; this can result in unwarranted confidence. The theory of conjugate duality can overcome the limitations of discretization within the 'strict bounds' formalism, a technique for constructing confidence intervals for functionals of the unknown model incorporating certain types of prior information. The usual computational approach to strict bounds approximates the 'primal' problem in a way that the resulting confidence intervals are at most long enough to have the nominal coverage probability. There is another approach based on 'dual' optimization problems that gives confidence intervals with at least the nominal coverage probability. The pair of intervals derived by the two approaches bracket a correct confidence interval. The theory is illustrated with gravimetric, seismic, geomagnetic, and helioseismic problems and a numerical example in seismology.

Stark, Philip B.↗
DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗
NASA NTRS2022

Extension and Validation of the Pendulum Model for Longitudinal Solar Prominence Oscillations

Context. Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semicircular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not semicircular. Aims. Our aim is to understand the effects of including the nonuniform solar gravity on longitudinal oscillations and explore the validity of the pendulum model with different flux-tube geometries. Methods. We first derived the equation describing the motion of the plasma along the flux tube including the effects of nonuniform gravity, yielding corrections to the original pendulum model. We also computed the full numerical solutions for the normal modes and compared them with the new pendulum approximation. Results. We find that the nonuniform gravity introduces a significant modification in the pendulum model. We also found a cut-off period; i.e., the longitudinal oscillations cannot have a period longer than 167 min. In addition, considering different tube geometries, the period depends almost exclusively on the radius of curvature at the bottom of the dip. Conclusions. We conclude that nonuniform gravity significantly modifies the pendulum model. These corrections are important for prominence seismology, because the inferred values of the radius of curvature and minimum magnetic-field strength differ substantially from those of the old model. However, we find that the corrected pendulum model is quite robust and is still valid for noncircular dips.

Sun: corona↗
NASA NTRS1994

Investigation of finite element: ABC methods for electromagnetic field simulation

The mechanics of wave propagation in the presence of obstacles is of great interest in many branches of engineering and applied mathematics like electromagnetics, fluid dynamics, geophysics, seismology, etc. Such problems can be broadly classified into two categories: the bounded domain or the closed problem and the unbounded domain or the open problem. Analytical techniques have been derived for the simpler problems; however, the need to model complicated geometrical features, complex material coatings and fillings, and to adapt the model to changing design parameters have inevitably tilted the balance in favor of numerical techniques. The modeling of closed problems presents difficulties primarily in proper meshing of the interior region. However, problems in unbounded domains pose a unique challenge to computation, since the exterior region is inappropriate for direct implementation of numerical techniques. A large number of solutions have been proposed but only a few have stood the test of time and experiment. The goal of this thesis is to develop an efficient and reliable partial differential equation technique to model large three dimensional scattering problems in electromagnetics.

Chatterjee, A.↗
NASA NTRS1997

On the Role of Resonances in Nonradial Pulsators

Resonances or near resonances are ubiquitous among the excited nonradial pulsation modes of variable stars and they must play an important role in determining their pulsational behavior. Here in a first step at nonlinear asteroseismology, we explore some of the consequences of resonances by means of the amplitude equation formalism. We show how parity and angular momentum constraints can be used to eliminate many of the possible nonlinear resonant couplings between modes (and multiplets of modes), and how the amplitude equations can thus be simplified. Even when we may not be able, nor wish, to make an ab initio computation of the values of the coupling coefficients, it is still possible to obtain constraints on the nature of the excited modes if a resonance between observed frequencies can be identified. Resonances can cause nonlinear frequency locking of modes. This means that the observed frequencies appear in exact resonance even though the linear frequencies are only approximately in resonance. The nonlinear frequency lock, when it occurs, it does so over a range of departures from linear resonance, and it is accompanied by constant pulsation amplitudes. The locked, nonlinear frequencies can differ noticeably from their nonresonant counterparts which are usually used in seismology. This is particularly true for multiplets of modes split by rotation. Beyond the regime of the frequency lock, amplitude and frequency modulations can appear in the pulsations. Far from the resonance condition one recovers the regime of steady pulsations with nonresonant frequencies for which the seismological studies, as they are presently carried out, are justified (provided furthermore, of course, that nonlinear frequency shifts are negligible). Success in identifying a resonance in an observed power spectrum depends on the quality of the data. While keeping this limitation in mind, ew discuss the possible existence of peculiar resonances the pulsations specific variable white dwarfs and delta Scuti stars.

Buchler, J. R.↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗