Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probabilistic Seismic Hazard Assessment for Azerbaijan

Probabilistic Seismic Hazard Assessments (PSHA) underpin the determination of seismic loads in most contemporary seismic provisions of building codes around the world. Modern building codes are migrating towards using the entire uniform hazard spectrum at a range of vibration periods (typically up to 4 s or 10 s) rather than a single peak value such as peak ground acceleration (PGA), thus requiring a larger range of PSHA outputs. The hazard maps in the current version of the building code of Azerbaijan (2011) are in terms of intensity and peak ground acceleration. Recognizing the need for an up-to-date seismic hazard assessment in the country, Seismic Cooperation Program (SCP) under the Lawrence Livermore National Laboratory (LLNL) undertook, in coordination with the Republican Seismic Survey Center of the Azerbaijan National Academy of Sciences and the Azerbaijan Scientific Research Institute of Construction and Architecture, a PSHA study that reflects new seismic data recorded locally and recent research conducted nationally and regionally since the last update to the building code. The PSHA framework for this project was designed to help develop a new earthquake catalogue, to incorporate a novel characterization of ground motions that specifically reflects the attenuation characteristics in the eastern Caucasus, and to generate hazard information in a form that is useful for an update of the current building code or the development of a new building code for Azerbaijan. The project also aimed to provide training and support for the local seismologists and engineers related to the seismic hazard models, probabilistic seismic hazard results and their use towards changes in the building code.

58 GEOSCIENCES↗
NASA NTRS2001

Large Torque Variations in Two Soft Gamma Repeaters

We have monitored the pulse frequencies of the two soft gamma repeaters SGR 1806-20 and SGR 1900+14 through the beginning of year 2001 using primarily Rossi X-ray Timing Explorer Proportional Counter Array observations. In both sources, we observe large changes in the spin-down torque up to a factor of approximately 4, which persist for several months. Using long baseline phase-connected timing solutions as well as the overall frequency histories, we construct torque noise power spectra for each SGR. The power spectrum of each source is very red (power-law slope approximately -3.5). These power spectra are consistent in normalization with some accreting systems, yet much steeper in slope than any known accreting source. To the best of our knowledge, torque noise power spectra with a comparably steep frequency dependence have only been seen in young, glitching radio pulsars (e.g. Vela). The observed changes in spin-down rate do not correlate with burst activity, therefore, the physical mechanisms behind each phenomenon are also likely unrelated. Within the context of the magnetar model, seismic activity cannot account for both the bursts and the long-term torque changes unless the seismically active regions are decoupled from one another.

Woods, Peter M.↗
NASA NTRS2002

Large Torque Variations in Two Soft Gamma Repeaters

We have monitored the pulse frequencies of the two soft gamma repeaters SGR 1806-20 and SGR 1900+14 through the beginning of year 2001 using primarily Rossi X-Ray Timing Explorer Proportional Counter Array observations. In both sources, we observe large changes in the spin-down torque up to a factor of approximately 4, which persist for several months. Using long-baseline phase-connected timing solutions as well as the overall frequency histories, we construct torque noise power spectra for each SGR (Soft Gamma Repeater). The power spectrum of each source is very red (power-law slope is approximately -3.5). The torque noise power levels are consistent with some accreting systems on timescales of approximately 1 yr, yet the full power spectrum is much steeper in frequency than any known accreting source. To the best of our knowledge, torque noise power spectra with a comparably steep frequency dependence have been seen only in young, glitching radio pulsars (e.g., Vela). The observed changes in spin-down rate do not correlate with burst activity; therefore, the physical mechanisms behind each phenomenon are also likely unrelated. Within the context of the magnetar model, seismic activity can not account for both the bursts and the long-term torque changes unless the seismically active regions are decoupled from one another.

Woods, Peter M.↗
DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2022

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗
NASA NTRS2020

Deep Learning Models for Planetary Seismicity Detection

Research in planetary seismology is fundamentally constrained by a lack of data. Seismo-logical science products of future missions can typically only be informed by theoretical signal/noise characteristics of the environment or likely Earth-analogues. Although objectives can be re-assessed after some initial data-collection upon lander arrival, transfer of high-resolution data back to Earth is costly on lander power usage. Over the last several years, development of GPU computing techniques and open-source high-level APIs have led to rapid advances in deep learning within the fields of computer vision, natural language processing, and collaborative filtering. These techniques are actively being adapted in seismology for a variety of tasks, including: earthquake detection, seismic phase discrimination, and ground-motion prediction. Until the recent detection of mars quakes during the Mars InSight mission, the only other measurements of seismicity recorded outside of Earth was on the Moon during the Apollo missions between 1969 to 1977. These unique data sets have been periodically revisited using new seismological methods, including ambient noise interferometry and Hidden Markov Models. Our objective is to develop a deep learning seismic detector and use it to catalog moonquakes from the Apollo 17 Lunar Seismic Profiling Experiment (LSPE) and compare the results with those obtained by other methods. Additionally, we will assess the accuracy tradeoff between using a training set of lunar data and one composed of Earth seismicity. In this document, we present preliminary results using a prototype classifier trained on a small set of earthquakes that was able to obtain detections for LSPE moonquakes with a greater accuracy than a recent study using Hidden Markov Models.

Civilini, F.↗
DOE OSTI2021

3D orthorhombic earth model effects on seismic source characterization

Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).

58 GEOSCIENCES↗
DOE OSTI2024

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗
DOE OSTI2021

Updates to the Regional Seismic Travel Time (RSTT) Model: 2. Path-dependent Travel-time Uncertainty

Abstract The regional seismic travel time (RSTT) model and software were developed to improve travel-time prediction accuracy by accounting for three-dimensional crust and upper mantle structure. Travel-time uncertainty estimates are used in the process of associating seismic phases to events and to accurately calculate location uncertainty bounds (i.e. event location error ellipses). We improve on the current distance-dependent uncertainty parameterization for RSTT using a random effects model to estimate slowness (inverse velocity) uncertainty as a mean squared error for each model parameter. The random effects model separates the error between observed slowness and model predicted slowness into bias and random components. The path-specific travel-time uncertainty is calculated by integrating these mean squared errors along a seismic-phase ray path. We demonstrate that event location error ellipses computed for a 90% coverage ellipse metric (used by the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (IDC)), and using the path-specific travel-time uncertainty approach, are more representative (median 82.5% ellipse percentage) of true location error than error ellipses computed using distance-dependent travel-time uncertainties (median 70.1%). We also demonstrate measurable improvement in location uncertainties using the RSTT method compared to the current station correction approach used at the IDC (median 74.3% coverage ellipse).

58 GEOSCIENCES↗
DOE OSTI2025

Panel 4 Discussion - Rowe

Explore the source record for details and available documents.

SMART Cables, global seismic model, local seismici↗
NASA NTRSDate not supplied

Mantle Plume Magmatism in Elysium Planitia as Constrained by InSight Seismic Observations

Based on cratering statistics, there has been active volcanism in the last 100 Ma and possibly the last 1 to 10 Ma in both central Elysium Planitia and in the Cerberus Fossae region of south-eastern Elysium. Seismic measurements by the InSight mission show a strong concentration of seismicity in Cerberus Fossae which are interpreted as due to magma transport along dike systems. Together, these observations strongly support the existence of geologically recent and possibly of currently active magmatism in Elysium. InSight seismic observations were also used to construct seismic velocity models for the martian mantle. These models are effectively a regional aver-age of the seismic velocity structure for the Elysium region, where both the InSight lander and most of the measured seismic events are located. These seismic velocity models have been interpreted in terms of the lithospheric thickness and mantle potential temperature. Khan et al. found a lithospheric thickness of 400 to 600 km and a mantle potential temperature of 1325 to 1425 °C. Durán et al. found a lithospheric thickness of ~450 km and a mantle potential tempera-ture of 1375-1475 °C. Drilleau et al. found a thermal lithosphere of 420-660 km and a mantle potential temperature of 1380-1560 °C. In this study, we combine these InSight seismic observations with finite element mantle plume magma production models to place new constraints on the properties of the mantle beneath Elysium Planitia.

Walter S Kiefer↗
DOE OSTI2023

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗