Analysis Of Ground Motion Research Articles

The strong ground motion and structural response analysis of 06 February 2023 Pazarcık and Elbistan earthquakes (Mw 7.7 and 7.6): a Case study for Malatya-Türkiye, Eastern Anatolia

The strong ground motion and structural response analysis of 06 February 2023 Pazarcık and Elbistan earthquakes (Mw 7.7 and 7.6): a Case study for Malatya-Türkiye, Eastern Anatolia

Seismic Ground Motion Analysis of Saturated Frozen Soil Free Field under S-Wave Incidence

Seismic Ground Motion Analysis of Saturated Frozen Soil Free Field under S-Wave Incidence

Modelling seismic ground motion and its uncertainty in different tectonic contexts: challenges and application to the 2020 European Seismic Hazard Model (ESHM20)

Abstract. Current practice in strong ground motion modelling for probabilistic seismic hazard analysis (PSHA) requires the identification and calibration of empirical models appropriate to the tectonic regimes within the region of application, along with quantification of both their aleatory and epistemic uncertainties. For the development of the 2020 European Seismic Hazard Model (ESHM20) a novel approach for ground motion characterisation was adopted based on the concept of a regionalised scaled-backbone model, wherein a single appropriate ground motion model (GMM) is identified for use in PSHA, to which adjustments or scaling factors are then applied to account for epistemic uncertainty in the underlying seismological properties of the region of interest. While the theory and development of the regionalised scaled-backbone GMM concept have been discussed in earlier publications, implementation in the final ESHM20 required further refinements to the shallow-seismicity GMM in three regions, which were undertaken considering new data and insights gained from the feedback provided by experts in several regions of Europe: France, Portugal and Iceland. Exploration of the geophysical characteristics of these regions and analysis of additional ground motion records prompted recalibrations of the GMM logic tree and/or modifications to the proposed regionalisation. These modifications illustrate how the ESHM20 GMM logic tree can still be refined and adapted to different regions based on new ground motion data and/or expert judgement, without diverging from the proposed regionalised scaled-backbone GMM framework. In addition to the regions of crustal seismicity, the scaled-backbone approach needed to be adapted to earthquakes occurring in Europe's subduction zones and to the Vrancea deep seismogenic source region. Using a novel fuzzy methodology to classify earthquakes according to different seismic regimes within the subduction system, we compare ground motion records from non-crustal earthquakes to existing subduction GMMs and identify a suitable-backbone GMM for application to subduction and deep seismic sources in Europe. The observed ground motion records from moderate- and small-magnitude earthquakes allow us to calibrate the anelastic attenuation of the backbone GMM specifically for the eastern Mediterranean region. Epistemic uncertainty is then calibrated based on the global variability in source and attenuation characteristics of subduction GMMs. With the ESHM20 now completed, we reflect on the lessons learned from implementing this new approach in regional-scale PSHA and highlight where we hope to see new developments and improvements to the characterisation of ground motion in future generations of the European Seismic Hazard Model.

Analysis of impulsive ground motions from February 2023 Kahramanmaraş earthquake sequence

AbstractOn the 6th of February 2023, a large magnitude earthquake (Pazarcık earthquake), $$M_{w}=$$ M w = 7.7, occurred in southeast Türkiye, which caused significant destruction in Türkiye and Syria. Relatively large magnitude aftershocks followed the main shock, and after 9 hours of the main event, another large magnitude earthquake (Elbistan earthquake) occurred, $$M_{w}=$$ M w = 7.6, on a nearby fault. This study analyzes the near-fault seismic signals from earthquakes larger than 5.5 recorded between the main shock and the 31st of March 2023. More than 60 impulsive motions are detected in 3 earthquakes, mostly concentrated in the Pazarcık and Elbistan earthquakes. In the Pazarcık earthquake, many impulsive motions are recorded in near-fault stations with periods of up to 14 s. In contrast, in the Elbistan earthquake, impulsive motions are spatially distributed, with pulse periods of up to 11 s and at distances greater than 150 km. Pulse periods mostly correlate with the magnitude of the earthquake, but pulse probability models do not predict impulsive motions over long distances. The presence of strong impulsive motions in vertical components is also observed. For both earthquakes, peak ground velocities (PGVs) are larger than predicted by ground motion prediction equations. The observation of long-period, large amplitude signals may indicate the presence of a directivity effect for both earthquakes. In some stations, spectral periods exceed the 2018 Turkish building design codes for long periods ($$\ge$$ ≥ 1 s).

A transformable wheel-spoke-paddle hybrid amphibious robot

AbstractThe intricate water-land intermingled nature of wild environments necessitates robots to exhibit multimodal cross-domain mobility capabilities. This paper introduces a novel wheel-spoke-paddle hybrid amphibious robot (WSP-bot) that can operate on flat and rough terrains, water surfaces, and water-land transitional zones. The proposed robot relies on a propulsion mechanism called transformable wheel-spoke-paddle (WSP), which combines the stability of wheeled robots with the obstacle-climbing capability of legged robots, while also providing additional aquatic mobility. The utilization of a crank-slider-based transformation mechanism enables seamless switching between multiple motion modes. An analysis of mode transition and ground motion in spoke mode was conducted, along with an investigation of its obstacle-crossing capability. Simulations were performed for mode transition, ground locomotion, and obstacle-crossing, as well as propulsion of a single WSP module on water. Based on the above work, a prototype robot was manufactured. Prototype tests, including mode transition and mobility tests on land and water surfaces under multimodal states, confirmed the effectiveness of the proposed WSP-bot.

Assessing seismic risk in the built environment of Istanbul: High-resolution hazard mapping and ground motion analysis in the sea of Marmara region

This research scrutinizes the seismic threat looming over Istanbul, potentially subject to a substantial earthquake. We analyze six plausible earthquake scenarios, utilizing six ground motion prediction equations (GMPEs), to forge high-resolution seismic hazard maps. These maps reveal not only peak horizontal ground accelerations but also spectral acceleration values across varying temporal frames, integrating the amplification effects of softer sediments. Our approach delineates that Istanbul’s western shoreline faces heightened risk, with median spectral accelerations at 0.3 s approaching 1 g, signifying intense shaking potential. In contrast, the area encompassing the financial district exhibits lower values, around 0.3 g. The granularity of these findings lays bare the seismic vulnerabilities of the region, offering a window into the risks and potential damages facing this bustling metropolis. This enhanced understanding paves the way for strategic urban planning and risk mitigation efforts aimed at safeguarding Istanbul’s populace and infrastructure. This article succinctly condenses our study’s pivotal conclusions, presenting a clarion call for proactive measures to diminish earthquake impacts on this dynamic urban landscape.

Study on physical explosion equivalent of large diameter and high pressure natural gas pipeline

High pressure natural gas pipelines are prone to material failure and explosion due to various reasons. Over the years, scholars have done in-depth research on the process of chemical explosion. In this paper, the energy components of high-pressure trunk natural gas pipeline under normal operation state were derived. Based on the assumption of adiabatic non isentropic, the pressure distribution in the pipeline within a period of time after the pipeline cracks was derived, and the energy released by the physical explosion of the pipeline was further derived. To verify the theoretical results, an original size pipe explosion experiment was organized, and the pressure distribution obtained from the experiment was found to be in good agreement with the theoretical calculation results. Through the crack propagation time and the duration of ground motion obtained from the experiment, the range of the duration of the pipeline burst crack impact was obtained, and the equivalent range of the pipeline explosion was further calculated. At the same time, through the regression analysis of ground motion, the equivalent of the pipe explosion in this experiment was about 177 kg TNT, which was consistent with the theoretical results. Further, the theoretical model and experimental results were verified by numerical simulation, and deeply studied the influence factors of buried depth on the explosion equivalent, as well as the relationship between the environment of the pipeline and the explosion damage mode and hazard range of the pipeline. The research results of this paper give the calculation method and magnitude range of physical explosion equivalent of high-pressure trunk natural gas pipeline. It provides a reference for the formulation of safety standards for high-pressure trunk natural gas pipelines and inert gas pressure vessels.

Reconnaissance of the 2022 Guanshan and Chihshang earthquakes in eastern Taiwan

Two significant earthquakes struck eastern Taiwan on consecutive days in mid-September of 2022. The first, the Guanshan earthquake, occurred on the night of September 17 (local time) with a local magnitude (ML) of 6.4, and the second, the Chihshang earthquake, occurred on the afternoon of September 18 with a mainshock of ML of 6.8. The strong motion of these two events resulted in a series of ground surface ruptures, 2 collapsed buildings, two collapsed bridges, and more than 100 partially damaged structures along the Chihshang and Yuli faults around the Longitudinal Valley in eastern Taiwan. The observed response spectra at the damage sites of the Chihshang earthquake were significantly larger than the design response spectra of the 1997 Taiwan building code. Near-fault velocity pulses with a maximum value of 132 cm/s were observed along the faults, destructively impacting the damaged structures. The details and causes of the structural damage are presented in this article according to the findings of the on-site reconnaissance and ground motion analysis. Finally, the behavior of selected structures that have a structural health monitoring system at the time of these two destructive earthquakes is also evaluated.

From ambient vibration data analysis to 1D ground-motion prediction of the Mj 5.9 and the Mj 6.5 Kumamoto earthquakes in the Kumamoto alluvial plain, Japan

We present horizontal ground motion predictions at a soft site in the Kumamoto alluvial plain for the Mj 5.9 and Mj 6.5 Kumamoto earthquakes of April 2016, in the framework of an international blind prediction exercise. Such predictions were obtained by leveraging all available information which included: (i) analysis of earthquake ground motions; (ii) processing of ambient vibration data (AMV); and (iii) 1D ground response analysis. Spectral analysis of earthquake ground-motion data were used to obtain empirical estimates of the prediction site amplification function, with evidence of an amplification peak at about 1.2 Hz. Horizontal-to-vertical spectral ratio analysis of AMV confirmed this resonance frequency and pointed out also a low-frequency resonance around 0.3 Hz at the prediction site. AMV were then processed by cross-correlation, modified spatial autocorrelation and high-resolution beamforming methods to retrieve the 1D shear-wave velocity (Vs) structure at the prediction site by joint inversion of surface-wave dispersion and ellipticity curves. The use of low frequency dispersion curve and ellipticity data allowed to retrieve a reference Vs profile down to few thousand meters depth which was then used to perform 1D equivalent-linear simulations of the M 5.9 event, and both equivalent-linear and nonlinear simulations of the M 6.5 event at the target site. Adopting quantitative goodness-of-fit metrics based on time–frequency representation of the signals, we obtained fair-to-good agreement between 1D predictions and observations for the Mj 6.5 earthquake and a poor agreement for the Mj 5.9 earthquake. In terms of acceleration response spectra, while ground-motion overpredictions were obtained for the Mj 5.9 event, simulated ground motions for the Mj 6.5 earthquake severely underestimate the observations, especially those obtained by the nonlinear approach.Graphical

Lumped mass models for use in predicting collapse of an isolated building

Development of fragility functions is an importantstep in the seismic evaluation of isolated buildings. To develop collapse fragility curves, a series of detailed dynamic ground motion analysis with multiple intensity measure levels must be conducted, often requiring considerable computational time. In addition, including more modeling detail to enable the prediction of multiple failure modes can make numerical convergence more challenging. Therefore, the possibility of simplified models to predict the collapse probabilityof isolated buildings is attractive. The collapse fragility curves derived from two lumped mass models with differing degrees of simplification are compared to a full nonlinear model. This is done for both sliding and rubber bearings with varying moat wall or restraining rim (for sliding bearings). The comparison demonstrates that the simplified model proposed in this study can be a potential alternative for estimating fragility functions of isolated buildings, as they can significantly improve the computational efficiency withrelatively smallcompromises on accuracy. This will allow for faster validation of system design when incorporating performance targets under large earthquakes.

Unsupervised detection of InSAR time series patterns based on PCA and K-means clustering

The need for implementing efficient value-adding tools able to optimise Earth Observation data usage, compels the scientific community to find innovative solutions for the downstream of Earth Observation information. In this paper we present an unsupervised and automated approach based on Principal Component Analysis (PCA) and K-means clustering to detect patterns of natural or anthropogenic ground deformation from Interferometric Synthetic Aperture Radar (InSAR) Time Series. For our proof-of-concept, we focus on the Valle d’Aosta region (Northwest Italy) where mass wasting processes frequently occurs, interacting with human activities and infrastructures. The large volumes of Sentinel-1 data produced allows for retrieving horizontal and vertical Time Series from multi-geometry data fusion of Line-of-Sight (LOS) InSAR measurements. The added benefit of combining ascending/descending InSAR data and interpolating displacements in time at different time steps is here explored prior to data dimensionality reduction and feature extraction through PCA. The retrieved principal components serve as a continuous solution for cluster membership indicators in the K-means clustering method, allowing to define spatially and temporally coherent displacement phenomena. The signal of the ground deformation clusters is then deconstructed into the underlying trend and seasonality components to enhance the interpretability of the classified satellite InSAR features. Using InSAR Time series data spanning 2014–2020, the proposed approach detects several slope movements and anthropogenic deformations with both linear and seasonal displacement behaviours. The results demonstrate the potential applicability of our transferable approach to the development of automated ground motion analysis systems.

Quantification of rupture directivity effects on strong ground motion during the 8 January 2022 MS6.9 Menyuan earthquake in Qinghai, China

On 8 January 2022, an MS6.9 earthquake occurred in the Menyuan county, Qinghai province at the northeastern margin of the Tibetan Plateau. A dense observatory network consisting of accelerographs and intensity meters has obtained many acceleration recordings in this region, which previously had only a few strong motion recordings. Using acceleration data recorded by 141 accelerographs and micro-electro mechanical system (MEMS)-based intensity meters within 100 km of the causative fault, we investigated the spatial distribution and attenuation characteristics of the near-source ground motion during this event. The observed response spectra were compared with the prediction obtained using the NGA-West2 models. It was found that, in the long-period section of >1.0 s, the acceleration response spectra in the rupture forward direction (southeast of the epicenter) were greater than the prediction, while observations on the northwest side were lower than the prediction. Through the spatial distribution of the ground motion and quantitative analysis of the residuals, we found that the rupture directivity effects on the ground motions were controlled by the spatial azimuth relative to the fault and the along-strike projected length of the rupture, which mainly affected the long-period portion of the >1.0 s. In the rupture forward direction, the spectral accelerations (SAs) with periods of 1.0–10.0 s were enhanced to 1.18–1.56 times the average. In the backward direction (northwest of the epicenter), the SAs with periods of 1.0–10.0 s were only 0.63–0.86 times the average. The rupture directivity effects during the MS6.9 (MW6.6) earthquake exhibited an obvious period-dependence and narrowband characteristics, and the amplification effects were the strongest at 5.0 s. The amplification factors in the straight forward directions were underestimated by the empirical model. Two typical velocity pulse records were observed in the rupture forward direction, and their response spectra were significantly amplified at long periods of greater than 1.0 s, with peak amplification values of 4.0–6.0 near the characteristic period of the pulses. This amplification was related to the resonance excited by a specific periodic pulse. Our study provides a reference of quantifying near-source effects for both moderate and major shallow crustal earthquakes by using MEMS recordings in this region.

Regression analysis on ground motion parameters for the earthquakes (Mw≥4.0) in NE India with special emphasis on 3 Jan 2016 M6.7, Tamenglong earthquake

In the present study, ground motion parameters (GMPs) are estimated by processing 26 earthquake events containing 125 accelerogram records with 375 components that originated in NE, India, and its vicinity with special emphasis on the Jan 3, 2016, M6.7, Tamenglong Earthquake. Moreover, Ground Motion Prediction Equations were developed through multiple regression analysis on observed data of 8 GMPs namely, Peak Ground Acceleration, Peak Ground Velocity, Arias Intensity, Characteristics Intensity, Housner Intensity, Cumulative Absolute Velocity, Effective Design Acceleration, and Acceleration Spectrum Intensity. The developed equations are representative of statistics on changes in amplitudes of parameters with varying distances and magnitude in connection to Northeast India. Besides, newly developed GMPEs can be applied to 4–6.8 magnitude earthquakes and valid up to 525 km of distance. In the case of January 3, 2016, Mw = 6.7, Tamenglong earthquake, a total of 25 parameters are estimated from ground motion recordings. The highest Peak Ground Acceleration (PGA) of 0.05 (g) has been experienced by Kohima station along with Pseudo Spectral Acceleration of 0.2802 (g) at the frequency of 1.72 Hz. The peak pseudo-spectral accelerations for each recording site of the Tamenglong earthquake are perceived within the periods of 0.2s and 0.5s. These estimations help in the characterization of ground motion with varying distances and magnitude which is a prime input to rate damage potential of ground motion and seismic hazard analysis in Northeast India.

Strong ground motion analysis of Assam using H/V ratio method

Strong ground motion analysis of Assam using H/V ratio method

Optimal stacking of noise cross-correlation functions

SUMMARY Cross-correlations of ambient seismic noise are widely used for seismic velocity imaging, monitoring and ground motion analyses. A typical step in analysing noise cross-correlation functions (NCFs) is stacking short-term NCFs over longer time periods to increase the signal quality. Spurious NCFs could contaminate the stack, degrade its quality and limit its use. Many methods have been developed to improve the stacking of coherent waveforms, including earthquake waveforms, receiver functions and NCFs. This study systematically evaluates and compares the performance of eight stacking methods, including arithmetic mean or linear stacking, robust stacking, selective stacking, cluster stacking, phase-weighted stacking, time–frequency phase-weighted stacking, Nth-root stacking and averaging after applying an adaptive covariance filter. Our results demonstrate that, in most cases, all methods can retrieve clear ballistic or first arrivals. However, they yield significant differences in preserving the phase and amplitude information. This study provides a practical guide for choosing the optimal stacking method for specific research applications in ambient noise seismology. We evaluate the performance using multiple onshore and offshore seismic arrays in the Pacific Northwest region. We compare these stacking methods for NCFs calculated from raw ambient noise (referred to as Raw NCFs) and from ambient noise normalized using a one-bit clipping time normalization method (referred to as One-bit NCFs). We evaluate six metrics, including signal-to-noise ratios, phase dispersion images, convergence rate, temporal changes in the ballistic and coda waves, relative amplitude decays with distance and computational time. We show that robust stacking is the best choice for all applications (velocity tomography, monitoring and attenuation studies) using Raw NCFs. For applications using One-bit NCFs, all methods but phase-weighted and Nth-root stacking are good choices for seismic velocity tomography. Linear, robust and selective stacking methods are all equally appropriate choices when using One-bit NCFs for monitoring applications. For applications relying on accurate relative amplitudes, the linear, robust, selective and cluster stacking methods all perform well with One-bit NCFs. The evaluations in this study can be generalized to a broad range of time-series analysis that utilizes data coherence to perform ensemble stacking. Another contribution of this study is the accompanying open-source software package, StackMaster, which can be used for general purposes of time-series stacking.

Seismic Safety Review of Gravity Dam Based on Finite Element Method

With the change of the seismic parameter zoning map of China (GB 18306-2015), the seismic grade of reservoir dams in some areas had changed. At the same time, the standard for seismic design of hydraulic structures (GB 51247-2018) also put forward new requirements for the seismic calculation of reservoir dams. In order to ensure the safe operation of reservoir dams, it is necessary to review the seismic safety based on the finite element numerical simulation technology. Taking the retaining dam section of a gravity dam as an example, a finite element model of the retaining dam section was established, the corresponding calculation and analysis of ground motion were carried out using the mode decomposition response spectrum method, and the seismic safety evaluation of the retaining dam section was carried out according to the calculation results. The results show the following: (1) after considering the effect of hydrodynamic pressure, the natural frequency of the dam body has been significantly reduced, and the first-order natural frequency has been reduced by about 10%; (2) in addition to the local tensile stress at the dam heel, the rest of the vertical stress on the foundation surface of the retaining dam section is compressive stress, and the length of the tensile stress is less than the distance from the dam heel to the curtain center line, which met the requirements; (3) the antisliding stability safety factor of the retaining dam section is greater than the design value in the specification and meets the safety requirements; and (4) the seismic safety of the retaining dam section meets the standard requirements, and the seismic grade is evaluated as grade A.

Construction and Usefulness Verification of Modeling Method of Subsurface Soil Layers for Numerical Analysis of Urban Area Ground Motion

Estimation of urban seismic damage using numerical simulation needs an automatic modeling method of surface layers and residential buildings. This study focuses on modeling of surface layers and shows a method of constructing models by interpolating boring data. An important property of the modeling method is robustness, that means that the method works for boring data with inconsistent soil layers. To satisfy this, we developed the method using artificial layers. We applied the method to a test site and checked its robustness. This test also showed that the method gave realistic models. Finally, we applied the method to the estimation of urban seismic damage and discussed the usefulness by comparing the result with one obtained by a conventional method.

Modeling response spectrum compatible pulse-like ground motion

The seismic response analysis of near-fault pulse-like ground motions is severely restricted due to the scarcity of pulse-like records. The requirement in regulations that the response spectra of artificial ground motions should be compatible with the target response spectrum makes the relevant studies more difficult. As a result, this study proposes a trigonometric series-based stochastic method to simulate pulse-like ground motions, with the advantage that the corresponding pseudo-spectral acceleration is compatible with the given target response spectrum. This goal is achieved by two parts. (1) The envelope function of pulse-like records obtained by the Hilbert transform is utilized as the amplitude modulation function to ensure that the simulated ground motion contains a pulse. (2) A novel iteration scheme based on random frequency parameters is proposed to guarantee the response spectrum compatibility. The velocity ground motion is first simulated since the pulse usually exists in velocity. The ground-motion acceleration subsequently obtained by differentiating the velocity is adopted to calculate the response spectrum. Two cases are implemented and verified the effectiveness of the proposed method in enriching existing pulse-like databases and generating pulse-like ground motion in areas that lack records. Moreover, the amplitude modulation function and target spectrum, as two key factors in the proposed method, determines the presence of a pulse and the pulse periods, respectively. This property makes the proposed method potentially universal applicability for stochastic pulse-like ground motion simulation in engineering.

Multidirectional Linear Site Response Analysis by Applying the SBSR Technique to KiK-net Data

Numerous works have highlighted the importance of the site response analysis of vertical ground motion in the earthquake engineering field. However, most studies have focused on nonlinear site response analysis. Hence, the knowledge of the linear site response analysis of vertical ground motion is inadequate. In this study, we perform multidirectional linear site response analysis by applying the surface-to-borehole spectral ratio (SBSR) technique to the Kiban–Kyoshin strong-motion observation network (KiK-net) seismic data. We find the hyperbolic relationships between multiple natural frequencies and the corresponding fundamental frequency. We propose empirical formulas, which can reduce the error in the multiple natural frequencies estimated by a single-layer model. Moreover, a brief investigation is conducted to provide an extra evidence for a previous correction method to the horizontal-to-vertical spectral ratio (HVSR) technique. The results of this study provide the empirical tools and can be utilized as a reference for the multidirectional linear site response analysis based on seismic observations.

Ground motion characteristics and their cumulative impacts to burst risks in underground coal mines

Ground motions resulting from seismic waves are one of the main factors to trigger coal burst hazards in underground mines. Ground motion induced dynamic impacts may apply to the highly stressed coal/rock mass and initiate dynamic failure. Although ground motion analysis has been widely used in burst-prone hard rock mines for designing support systems, there are only few studies focusing on the understanding ground motion characteristics in coal mines and their relationship with the coal burst damage. Therefore, based on nine-months seismic monitoring in a case study coal mine, this paper conducted a thorough analysis on quantifying ground motions to roadways in both far-field and near-field zones and correlating strong ground motions to actual coal burst damage. The results showed that most far-field ground motions were insignificant, which were less likely to initiate coal burst damage. However, with the same energy levels and hypocentral distances, the seismic events in coal mines can produce higher far-field ground motions than those in hard rock mines. Compared with the far-field ground motions, the near-field ground motions had much higher intensities. The strong dynamic impacts induced by the near-field ground motions may trigger dynamic failure in coal or rock mass when it is already critically stressed. An index called the number of high ground motions (NHGM) was developed to link ground motions with burst risks in roadways. The roadway zone that had experienced a long history of intensive ground motions, indicated by a higher degree of NHGM, were more likely to expect coal bursts in the future.Article highlightsQuantify ground motions to roadways in both far-field and near-field zones in underground coal mines.Demonstrate the relationship between ground motion intensity and coal burst damage.Develop a method to forecast potential burst damage zones in roadways based on historical ground motions.