Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earthquakes”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

The 2003 Bam (Iran) earthquake: Rupture of a blind strike-slip fault

A magnitude 6.5 earthquake devastated the town of Bam in southeast Iran on 26 December 2003. Surface displacements and decorrelation effects, mapped using Envisat radar data, reveal that over 2 m of slip occurred at depth on a fault that had not previously been identified. It is common for earthquakes to occur on blind faults which, despite their name, usually produce long-term surface effects by which their existence may be recognised. However, in this case there is a complete absence of morphological features associated with the seismogenic fault that destroyed Bam.

earthquakes blind faults Bam Envisat radar remote ↗
NASA NTRS1996

Surface Deformation Due to the May 27, 1995 Sakhalin Earthquake and Related Events Measured by JERS-1 SAR Interferometry

A large (M&subw;=7.0) earthquake on May 27, 1995 completely destroyed the town of Neftegorsk in the northern part of Sakhalin Island and caused more than 2000 human deaths. The shallow, right-lateral, strick-slip earthquake resulted in extensive surface ruptures and up to 7 m of horizontal displacement as reported by field workers. The sourthern part of the mainshock epicenter zone was imaged by the JERS-1 SAR (synthetic aperature radar) one month (April 28) before and two weeks after (June 11) the mainshock. Despite drastically changed surface conditions in the 44 days between the two images, due primarily to spring thaw, we obtained reasonably good interferometric correlation with the L-band (24 cm) SAR pair. The interoferogram records the distribution of deformation reflecting displacement during both the mainshock and aftershocks. The ability to map the deformation pattern can aid the assessment and mitigation of damage.

Sakhalin Island JERS-1↗
NASA NTRS2011

GPS Time Series Analysis of Southern California Associated with the 2010 M7.2 El Mayor/Cucapah Earthquake

The Magnitude 7.2 El-Mayor/Cucapah earthquake the occurred in Mexico on April 4, 2012 was well instrumented with continuous GPS stations in California. Large Offsets were observed at the GPS stations as a result of deformation from the earthquake providing information about the co-seismic fault slip as well as fault slip from large aftershocks. Information can also be obtained from the position time series at each station.

earthquake↗
NASA NTRS2011

Integrating Machine Learning into a Crowdsourced Model for Earthquake-Induced Damage Assessment

On January 12th, 2010, a catastrophic 7.0M earthquake devastated the country of Haiti. In the aftermath of an earthquake, it is important to rapidly assess damaged areas in order to mobilize the appropriate resources. The Haiti damage assessment effort introduced a promising model that uses crowdsourcing to map damaged areas in freely available remotely-sensed data. This paper proposes the application of machine learning methods to improve this model. Specifically, we apply work on learning from multiple, imperfect experts to the assessment of volunteer reliability, and propose the use of image segmentation to automate the detection of damaged areas. We wrap both tasks in an active learning framework in order to shift volunteer effort from mapping a full catalog of images to the generation of high-quality training data. We hypothesize that the integration of machine learning into this model improves its reliability, maintains the speed of damage assessment, and allows the model to scale to higher data volumes.

crowdsourcing↗
NASA NTRS2023

Learnings From Rapid Response Efforts to Remotely Detect Landslides Triggered By the August 2021 Nippes Earthquake and Tropical Storm Grace in Haiti

On August 14, 2021, a Mw 7.2 earthquake struck the Tiburon Peninsula of western Haiti triggering thousands of landslides. Three days after the earthquake on August 17, 2021, Tropical Storm Grace crossed shallow waters offshore of southern Haiti triggering more landslides worsening the situation. In the aftermath of these events, several organizations with disaster response capabilities or programs activated to provide information on the location of landslides to first responders on the ground. Utilizing remote sensing to support rapid response, one organization manually mapped initiation point of landslides and three automatically detected landslides. The 2021 Haiti event also provided a unique opportunity to test different automated landslide detection methods that utilized both SAR and optical data in a rapid response scenario where rapid situational awareness was critical. As the methods used are highly replicable, the main goal of this study is to summarize the landslide rapid response products released by the organizations, detection methods, quantify accuracy and provide guidelines on how some of the shortcomings encountered in this effort might be addressed in the future. To support this validation, a manually mapped polygon-based landslide inventory covering the entire affected area was created and is also released through this effort.

2021 Nippes earthquake↗
DOE OSTI2022

Predicting deviations in ex-core detector response in Krsko NPP during an earthquake with Monte-Carlo neutron transport methods

To study the response of ex-core detectors during an earthquake, an approach using Monte Carlo code MCNP was developed. As a first step to predict the response of the ex-core detectors a detailed MCNP model of the Krsko reactor core was created. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krsko nuclear power plant. Due to the outside the core position of the neutron detectors, the hybrid code ADVANTG is used to generate weight windows to speed up neutron transport outside the reactor core. To be able to use ADVANTG, fixed neutron source had to be reconstructed from the criticality core calculation. It was found that the power range ex-core detector movement has a negligible effect on the value of the thermal neutron flux at the active part of the detector. However, the 5 cm radial movement of the intermediate-range detector leads to 7 % - 8 % change in thermal neutron flux within the active intermediate-range detector region. The analysis continued with the evaluation of the effect of core barrel movement on the ex-core detector response. It was determined that the 2 mm core barrel radial oscillation can lead to ∼ 4% change in thermal neutron flux within the active detector region. It was evaluated that reactor core movement inside the baffle can lead up to ∼ 1.9% and ∼ 1.3% deviation in power-range and intermediate-range detector signal respectively. Analysis showed that the mechanical movement of ex-core neutron detectors could not explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗
DOE OSTI2025

Earthquake-triggered submarine canyon flushing transfers young terrestrial and marine organic carbon into the deep sea

Submarine canyons transfer substantial amounts of sediment and organic carbon (OC) into the deep ocean, nourishing deep-sea ecosystems and contributing to the global carbon cycle through OC burial and sequestration. Tracking lateral OC transport through submarine canyon systems is challenged by the deep-ocean setting, difficulties with constraining episodic depositional events, and the need to assess the composition and age of marine and terrestrial organic matter. We apply innovative parallel ramped pyrolysis oxidation accelerator mass spectrometry and pyrolysis-gas chromatography-mass spectrometry with isotope analyses to track OC age and sources in the 2016 Kaikoura earthquake-triggered, canyon-flushing event that deposited along >1300 km of a submarine canyon-channel system, offshore Aotearoa New Zealand. Specifically, these techniques allow us to determine the ages, sources, and partitioning of OC within the Kaikoura turbidite deposit and test hypotheses of how submarine canyon systems contribute to lateral OC flux and burial. Our results show that, despite considerable canyon floor erosion, substantial amounts of young OC were flushed into the deep sea, with relatively little (~2%) pre-Holocene OC contributions. Even without a direct connection between rivers and submarine canyons, most (~55%) of the OC in the Kaikoura event bed is from terrestrial sources. However, the deposit also contains substantial amounts (~22%) of marine-derived OC and ~23% of the material is of unassignable origin. Particle sorting imparts variability on the age and composition of OC within turbidite deposits and along the turbidity current flow path. Terrestrial-derived OC is preferentially older than marine-derived OC and concentrated in coarser particle sizes found more commonly at the deposit base and in proximal settings. Young, marine-derived OC is concentrated at the surface of the deposits and tends to be enriched in finer particle sizes. Such OC partitioning in turbidites supports the relevance of depositional models for predicting and quantifying distribution of OC in deep-sea deposits. Earthquake-triggered, canyon flushing events and resulting turbidites enhance OC burial efficiency and can sequester OC effectively, contributing an important carbon sink to the sedimentary carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗
DOE OSTI2022

Modeling Sequences of Earthquakes and Aseismic Slip (SEAS) in Elasto-Plastic Fault Zones With a Hybrid Finite Element Spectral Boundary Integral Scheme

We present a coupled finite element spectral boundary integral framework for modeling sequences of earthquakes and aseismic slip on a 2-D planar rate-and-state fault with off-fault visco-plastic response in the plane strain approximation. The model resolves both slow aseismic deformation and inertia effects during rapid slip. As an application, we perform two sets of simulations with different choices of cohesion to explore the co-evolution of fault slip, bulk plasticity and local stress fields. The first set implements a relatively large value of the cohesion parameter, which results in limiting inelastic strain accumulation to dynamic rupture phases. The second set implements a smaller cohesion, allowing for plastic strain to accumulate in both seismic and aseismic phases. For the first model, our results indicate that the extent and distribution of plastic strain depend on the angle of maximum compressive principal stress. At larger angles, inelastic strain accumulates on the extensional side of a dynamically propagating rupture. At smaller angles, the extent of plasticity is limited to the compressional side of the domain. At smaller cohesion values, off-fault plasticity may occur during aseismic slip, which alters the nucleation characteristics and earthquake sequence pattern. Furthermore, our results at lower cohesion values indicate that plastic strain accumulation may occur in both the extensional and compressional sides of the off-fault bulk even at higher angles of maximum compression. In conclusion, this produces damage patterns that deviate from the traditional off-fault fan-like distribution observed in dynamic rupture simulations and emphasizes the significance of long-term deformation in interpreting observations.

58 GEOSCIENCES↗
DOE OSTI2025

Severe Strain‐Induced Olivine‐Ringwoodite Transformation at Room Temperature: Key to Enigmas of Deep‐Focus Earthquake

Deep‐focus earthquakes at 350–660 km are presumably caused by olivine‐spinel phase transformation (PT). This cannot, however, explain the observed high seismic strain rate, which requires PT to complete within seconds, while metastable olivine does not transform for over a million years. Recent theory quantitatively describes how severe plastic deformations (SPD) can solve this dilemma but lacking experimental proof. Here, we introduce dynamic rotational diamond anvil cell with rough diamond anvils to impose SPD on San Carlos olivine. While olivine never transformed to spinel at room temperature, we obtained reversible olivine‐ringwoodite PT under SPD at 15–28 GPa within tens of seconds. The PT pressure reduces with increasing dislocation density, microstrain, plastic strain, and decreasing crystallite size. Results demonstrate a new strain‐induced PT mechanism compared to a pressure/temperature‐induced one. Combined with SPD during olivine subduction, this mechanism can accelerate olivine‐ringwoodite PT from millions of years to timescales relevant to earthquakes.

Lin, F. [Iowa State Univ., Ames, IA (United States↗
DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗
DOE OSTI2025

AI for Earthquake Physics

The core LANL program sponsored by Office of Science, Basic Energy Science, Chemical Sciences, Geosciences, and Biosciences (DOE-BES-CSGB) and led by PI Johnson aims to research earthquake faults to advance fault physics and earthquake hazards. All work completed is required to be made publicly available through publications and open-source codes supporting the published results. All routines are/will-be written in open source python and applied to publicly available data sets. These routines will format data from input into models, develop and test modeling frameworks for the problems addressed, and produce figures applicable to peer-reviewed manuscripts. All work is reviewed for Los Alamos Unlimited Release before submitting to a journal. This summary encompasses recently completed work and work to be complete for the duration of the program.

Johnson, Christopher↗
DOE OSTI2021

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗
DOE OSTI2023

Evaluating the Aftershock Duration of Induced Earthquakes

As the carbon sequestration community prepares to scale up the number and size of commercial operations, the need for tools and methods to assess and mitigate risks associated with these operations becomes increasingly important. One outstanding question is whether aftershocks of induced events decay quickly after injection operations cease or if aftershock activity persists for hundreds of years before returning to background levels more akin to tectonic events (Stein and Liu, 2009). Appropriate estimates of the aftershock duration impact several operational management decisions including mitigation strategies and post-injection monitoring for seismic activity. It is hypothesized that induced earthquake rates may diminish more quickly after injection is stopped, owing to higher stressing rates from injected fluids. Alternatively, it is plausible that only the first event in the sequence is induced by increased fluid overpressures, whereas subsequent events (e.g., aftershocks) respond to the stored tectonic stresses and static and dynamic stress changes due to the mainshock (Keranen et al., 2013). Further, if the aftershock duration can be linked to stressing rates due to injection, then it follows that operational strategies to reduce seismic hazard by reducing injection rates or volumes may be successful. However, if aftershocks of induced events are relieving stored tectonic stresses, then altering injection volumes may not alleviate ongoing seismic activity. Furthermore, knowledge of an aftershock duration could aid in the determination of an appropriate post-injection monitoring period for ongoing seismicity, which is a factor in overall operational costs. In this study, we model induced seismicity sequences in Oklahoma with a coupled Coulomb rate–state earthquake rate model (Dieterich, 1994; Kroll et al., 2017) to estimate aftershocks durations. Results for the current study indicate that elevated rates of aftershock activity following induced mainshocks return to background seismicity rates in less than five years, contrary to the tens to hundreds of years observed for tectonic aftershocks.

58 GEOSCIENCES↗
DOE OSTI2025

Importance of Considering Near-Surface Attenuation in Earthquake Source Parameter Estimation: Insights from Kappa at a Dense Array in Oklahoma

ABSTRACT Separating earthquake source spectra from propagation effects is challenging. The propagation effect contains a site-dependent term related to the high attenuation of shallow sediments. Neglecting the site-dependent attenuation can cause large biases and scattering in the corner-frequency (fc) estimates, resulting in significant stress-drop deviations. In this study, we investigate shallow attenuation at the LArge-n Seismic Survey in Oklahoma (LASSO) and site-related biases and scattering in source parameter measurements due to simplified attenuation models. We measure the high-frequency spectral decay parameter kappa on the vertical acceleration spectra of regional earthquakes (125 km away). The site-dependent kappa (κ0,acc) suggests that attenuation increases rapidly at shallow depth and is highly site-dependent. 10%–75% of the attenuation is site-dependent for S waves and even larger for P waves. The quality factor for S waves (QS) ranges from 10 to 100 in the upper 400 m. QP for P waves is mostly below 10 within the same depth. The Quaternary sediments tend to be more attenuating (QS<30), but the Permian rocks also can have high attenuation. We demonstrate that using a non-site-dependent attenuation model in single-spectra fitting leads to large scattering in fc estimates among stations with apparent good fits. The apparent fc can significantly deviate when the range of site-dependent kappa is large or with a higher assumed source spectral fall-off rate n. The biases in apparent fc depend on site condition and distance; however, the correlation between fc and these factors might not be obvious, depending on model assumptions. An apparent increase of stress drop with magnitude in a previous study for local microearthquakes (1.3

Chang, Hilary↗