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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

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↗

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.↗

DOE OSTI2025

Analysis and optimization of seismic monitoring networks with Bayesian optimal experimental design

SUMMARY Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network’s ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. This framework requires five elements: (i) A likelihood function that describes the distribution of detection and traveltime data from the sensor network, (ii) A prior distribution that describes a priori belief about seismic events, (iii) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, (iv) An algorithm to compute EIG about seismic events over a data set of hypothetical prior events, (v) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number and locations influence uncertainty; and how prior models and constraints influence sensor placement.

58 GEOSCIENCES↗

NASA NTRS1993

New progress of ranging technology at Wuhan Satellite Laser Ranging Station

A satellite laser ranging system with an accuracy of the level of centimeter has been successfully developed at the Institute of Seismology, State Seismological Bureau with the cooperation of the Institute of Geodesy and Geophysics, Chinese Academy of Science. With significant improvements on the base of the second generation SLR system developed in 1985, ranging accuracy of the new system has been upgraded from 15 cm to 3-4 cm. Measuring range has also been expanded, so that the ETALON satellite with an orbit height of 20,000 km launched by the former U.S.S.R. can now be tracked. Compared with the 2nd generation SLR system, the newly developed system has the following improvements. A Q modulated laser is replaced by a mode-locked YAG laser. The new device has a pulse width of 150 ps and a repetition rate of 1-4 pps. A quick response photomultiplier has been adopted as the receiver for echo; for example, the adoption of the MCP tube has obviously reduced the jitter error of the transit time and has improved the ranging accuracy. The whole system is controlled by an IBM PC/XT Computer to guide automatic tracking and measurement. It can carry out these functions for satellite orbit calculation, real-time tracking and adjusting, data acquisition and the preprocessed of observing data, etc. The automatization level and reliability of the observation have obviously improved.

Xia, Zhiz-Hong↗

DOE OSTI2024

Integration of Waveform Simulation Methods

The generation of synthetic seismograms through simulation is a fundamental tool of seismology required to run quantitative hypothesis tests. A variety of approaches have been developed throughout the seismological community and each has their own specific user interface based on their implementation. This causes a challenge to researchers who will need to learn new interfaces with each new software they wish to use and create substantial challenges when attempting to compare results from different tools. Here we provide a unified interface that facilitates interoperability amongst several simulation tools through a modern containerized Python package. Further, this package includes post-processing analysis modules designed to facilitate end-to-end analysis of synthetic seismograms. In this report we present the conceptual guidance and an example implementation of the new Waveform Simulation Framework.

58 GEOSCIENCES↗

DOE OSTI2022

Real-time and post-hoc compression for data from Distributed Acoustic Sensing

Distributed Acoustic Sensing (DAS) is an emerging sensing technology that records the strain-rate along fiber optic cables at high spatial and temporal resolution. This technique is becoming a popular tool in seismology, hydrology, and other subsurface monitoring applications. However, due to the large coverage (10’s of km) and high density of measurements (1m spacing at 100’s of Hz), a DAS installation could produce terabytes of data records per day. Because many DAS instruments are deployed in remote locations, this large data size poses significant challenges to its transfer and storage. In this paper, we explore lossless compression methods to reduce the storage requirement in both real-time and post-hoc scenarios. Here we propose a two-stage compression method to improve the compression ratio and compression speed. This two-stage compression method could reduce the storage requirement by 40%, which is 20% more than other lossless methods, such as ZSTD. We demonstrate that the compression method could complete its operation well before the DAS instrument needs to output the next file, making it suitable for real-time DAS acquisition. We also implement a parallel compression method for a post-hoc scenario and demonstrate that our method could effectively utilize a parallel computer. With 256 CPU cores, our parallel compression method achieves the speed of 26GB/second.

58 GEOSCIENCES↗

DOE OSTI2021

Pycheron: A Python-Based Seismic Waveform Data Quality Control Software Package

Supplementing an existing high-quality seismic monitoring network with openly available station data could improve coverage and decrease magnitudes of completeness; however, this can present challenges when varying levels of data quality exist. Without discerning the quality of openly available data, using it poses significant data management, analysis, and interpretation issues. Incorporating additional stations without properly identifying and mitigating data quality problems can degrade overall monitoring capability. If openly available stations are to be used routinely, a robust, automated data quality assessment for a wide range of quality control (QC) issues is essential. To meet this need, we developed Pycheron, a Python-based library for QC of seismic waveform data. Pycheron was initially based on the Incorporated Research Institutions for Seismology’s Modular Utility for STAtistical kNowledge Gathering but has been expanded to include more functionality. Pycheron can be implemented at the beginning of a data processing pipeline or can process stand-alone data sets. Its objectives are to (1) identify specific QC issues; (2) automatically assess data quality and instrumentation health; (3) serve as a basic service that all data processing builds on by alerting downstream processing algorithms to any quality degradation; and (4) improve our ability to process orders of magnitudes more data through performance optimizations. This article provides an overview of Pycheron, its features, basic workflow, and an example application using a synthetic QC data set.

58 GEOSCIENCES↗

DOE OSTI2025

The Dynamic Networks Experiments: Virtual Experiments to Quantify Gains in Nuclear Explosion Monitoring

We describe an ongoing series of virtual experiments conducted collaboratively by four United States National Laboratories: Sandia National Laboratories, Los Alamos National Laboratory, Lawrence Livermore National Laboratory, and Pacific Northwest National Laboratory. These Dynamic Network Experiments (DNEs) provide an experimental framework to evaluate the potential impact of new research tools on nuclear explosion monitoring. The second DNE (DNE2), completed in 2024, exploited waveform data (seismic, infrasound, and electromagnetic) that was recorded by multi-modal sensors within and near the Nevada National Security Site and synthetic radionuclide signatures over multiple time periods. During the execution of DNE2, we processed and analyzed data through a multi-stage event processing pipeline that ingested raw data, performed quality control, detected signals, built events from these signals, located these events, and characterized the events’ source types and sizes. For each stage and over the entire event processing pipeline, we evaluated performance changes by comparing the performance of new data processing methods, models, and algorithms against a baseline. We also performed an additional execution phase to assess event processing pipeline function, speed, and efficiency against that of an expert analyst, including computational and manual efforts. Finally, we assessed the impact and effort of modern computing infrastructure on the monitoring pipeline. This paper describes key elements of the DNEs, from formulation through execution, as demonstrated in DNE2. The DNEs introduce several novel concepts to quantitatively measure the potential impact of new methods on explosion monitoring, including the collaborative design of multi-modal datasets, performance and logistical metrics, and integrated analyses.

42 ENGINEERING↗

DOE OSTI2025

Advances in geophysical forensic event monitoring

Forensic analysis of man-made, non-nuclear events (such as industrial accidents, explosion experiments and mine collapses) has become more frequent and detailed owing to advancements in geophysical monitoring. Here, in this Technical Review, we demonstrate how geophysical forensic monitoring using seismic, infrasound and hydroacoustic recordings provides insights on events in the solid earth, atmosphere and underwater. Advanced techniques, including machine-learning-based models, have been developed to detect, identify and investigate these events, providing information on location, subevents, sources and explosive yield. The increase in data availability, application of advanced methods and computation and the growth of multitechnology approaches have increased the accuracy of forensic event analysis and enabled more realistic characterization of uncertainties. For example, the 2020 Beirut explosion in Lebanon demonstrated that various seismic, acoustic and other methods could be used to estimate explosive yield (and yield uncertainties) of about 1 ktonne, providing confidence in the application of these methods to smaller events where data are available. However, forensic investigations remain largely limited to known events with identified sources. Increased access to data, sophisticated analysis methods and high-resolution earth models will improve forensic event analysis further, enabling civil and scientific applications, such as localization in the search for the lost ARA San Juan submarine.

geophysics↗

DOE OSTI2021

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗