Stochastic Source Model for Strong Motion Prediction

Author(s):  
S. Sangeetha ◽  
S.T.G. Raghukanth

The article aims at developing a stochastic model which simulates spatial distribution of slip on the fault plane. This is achieved by analysing a large dataset of 303 finite-fault rupture models from 152 past earthquakes with varying fault mechanisms and in the magnitude range of 4.11-9.12. New scaling relations to predict the seismic source parameters such as fault length, fault width, rupture area, mean and standard deviation of slip have been derived for distinct fault mechanisms. The developed methodology models the spatial variability of slip as a two-dimensional von Karman power spectral density function (PSD) and correlation lengths are estimated. The proposed stochastic slip model is validated by comparing the simulated near-field ground response with the recorded data available for the 20th September 1999 Chi-Chi earthquake, Taiwan.

Author(s):  
Hernando Tavera ◽  
Bertrand Delouis ◽  
Arturo Mercado ◽  
David Portugal

Abstract The Loreto earthquake of 26 May 2019 occurred below the extreme northeast part of Peru at a depth of 140 km within the subducting Nazca plate at a distance of 700 km from the trench Peru–Chile. The orientation of the seismic source was obtained from waveform inversion in the near field using velocity and strong-motion data. The rupture occurred in normal faulting corresponding to a tensional process with T axis oriented in east–west direction similar to the direction of convergence between the Nazca and South America plates. The analysis of the strong-motion data shows that the levels of ground shaking are very heterogeneous with values greater than 50 Gal up to distances of 300 km; the maximum recorded acceleration of 122 Gal at a distance of 100 km from the epicenter. The Loreto earthquake is classified as a large extensional event in the descending Nazca slab in the transition from flat-slab geometry to greater dip.


Author(s):  
Basil Margaris ◽  
Emmanuel M. Scordilis ◽  
Jonathan P. Stewart ◽  
David M. Boore ◽  
Nikos Theodoulidis ◽  
...  

Abstract We present a Hellenic database of intensity measures from uniformly processed strong ground motion recordings, together with metadata on earthquake source attributes and recording site conditions. The database consists of information from 471 earthquakes between 1973 and 2015 that produced 2993 usable recordings from 333 sites. A key element of this work is a unified presentation of data from two major data providers that operate in Greece (Institute of Engineering Seismology and Earthquake Engineering and the Institute of Geodynamics, National Observatory of Athens) along with a university-operated local urban array (University of Patras). Consistent procedures were applied to develop source parameters that include hypocenter locations, moment magnitudes (directly estimated or derived using a conversion procedure), fault-plane solutions, and finite-fault parameters (generally, for events with M>6.0). The time-averaged shear-wave velocity in the upper 30 m parameter is provided for all 333 recording sites based on geophysical measurements where available (102) and proxy-based estimates otherwise. Most events are in the magnitude range of 3.8–7, occurred at shallow hypocentral depths (<30  km), and provide data for rupture distances generally between 10 and 300 km. The combined ground motion, seismic source, and site database is anticipated to be useful for engineering applications, including ground-motion model development and time series selection for response-history analyses.


2019 ◽  
Vol 219 (1) ◽  
pp. 514-521 ◽  
Author(s):  
Y Wang ◽  
K Satake ◽  
R Cienfuegos ◽  
M Quiroz ◽  
P Navarrete

SUMMARY The 2015 Illapel earthquake (Mw 8.3) occurred off central Chile on September 16, and generated a tsunami that propagated across the Pacific Ocean. The tsunami was recorded on tide gauges and Deep-ocean Assessment and Reporting of Tsunami (DART) tsunameters in east Pacific. Near-field and far-field tsunami forecasts were issued based on the estimation of seismic source parameters. In this study, we retroactively evaluate the potentiality of forecasting this tsunami in the far field based solely on tsunami data assimilation from DART tsunameters. Since there are limited number of DART buoys, virtual stations are assumed by interpolation to construct a more complete tsunami wavefront for data assimilation. The comparison between forecasted and observed tsunami waveforms suggests that our method accurately forecasts the tsunami amplitudes and arrival time in the east Pacific. This approach could be a complementary method of current tsunami warning systems based on seismic observations.


2019 ◽  
Vol 219 (2) ◽  
pp. 1148-1162
Author(s):  
Jiun-Ting Lin ◽  
Wu-Lung Chang ◽  
Diego Melgar ◽  
Amanda Thomas ◽  
Chi-Yu Chiu

SUMMARY We test the feasibility of GPS-based rapid centroid moment tensor (GPS CMT) methods for Taiwan, one of the most earthquake prone areas in the world. In recent years, Taiwan has become a leading developer of seismometer-based earthquake early warning systems, which have successfully been applied to several large events. The rapid determination of earthquake magnitude and focal mechanism, important for a number of rapid response applications, including tsunami warning, is still challenging because of the limitations of near-field inertial recordings. This instrumental issue can be solved by an entirely different observation system: a GPS network. Taiwan is well posed to take advantage of GPS because in the last decade it has developed a very dense network. Thus, in this research, we explore the suitability of the GPS CMT inversion for Taiwan. We retrospectively investigate six moderate to large (Mw6.0 ∼ 7.0) earthquakes and propose a resolution test for our model, we find that the minimum resolvable earthquake magnitude of this system is ∼Mw5.5 (at 5 km depth). Our tests also suggest that the finite fault complexity, often challenging for the near-field methodology, can be ignored under such good station coverage and thus, can provide a fast and robust solution for large earthquake directly from the near field. Our findings help to understand and quantify how the proposed methodology could be implemented in real time and what its contributions could be to the overall earthquake monitoring system.


2020 ◽  
Author(s):  
Sebastian Heimann ◽  
Marius Isken ◽  
Daniela Kühn ◽  
Hannes Vasyura-Bathke ◽  
Henriette Sudhaus ◽  
...  

<p>Seismic source and moment tensor waveform inversion is often ill-posed or non-unique if station coverage is poor or signals are weak. Three key ingredients can help in these situations: (1) probabilistic inference and global search of the full model space, (2) joint optimisation with datasets yielding complementary information, and (3) robust source parameterisation or additional source constraints. These demands lead to vast technical challenges, on the performance of forward modelling, on the optimisation algorithms, as well as on visualisation, optimisation configuration, and management of the datasets. Implementing a high amount of automation is inevitable.</p><p>To tackle all these challenges, we are developing a sophisticated new seismic source optimisation framework, Grond. With its innovative Bayesian bootstrap optimiser, it is able to efficiently explore large model spaces, the trade-offs and the uncertainties of source parameters. The program is highly flexible with respect to the adaption to specific source problems, the design of objective functions, and the diversity of empirical datasets.</p><p>It uses an integrated, robust waveform data processing, and allows for interactive visual inspection of many aspects of the optimisation problem, including visualisation of the result uncertainties. Grond has been applied to CMT moment tensor and finite-fault optimisations at all scales, to nuclear explosions, to a meteorite atmospheric explosion, and to volcano-tectonic processes during caldera collapse and magma ascent. Hundreds of seismic events can be handled in parallel given a single optimisation setup.</p><p>Grond can be used to optimise simultaneously seismic waveforms, amplitude spectra, waveform features, phase picks, static displacements from InSAR and GNSS, and gravitational signals.</p><p>Grond is developed as an open-source package and community effort. It builds on and integrates with other established open-source packages, like Kite (for InSAR) and Pyrocko (for seismology).</p>


2021 ◽  
Author(s):  
Pınar Büyükakpınar ◽  
Mohammadreza Jamalreyhani ◽  
Mehdi Rezapour ◽  
Stefanie Donner ◽  
Nima Nooshiri ◽  
...  

<p>In May 2020 an earthquake with Mw 5.0 struck at ~40 km east of Tehran metropolis and ~15 km south of the Damavand stratovolcano. It was responsible for 2 casualties and 23 injured. The mainshock was preceded by a foreshock with Ml 2.9 and followed by a significant aftershock sequence, including ten events with Ml 3+. The occurrence of this event raised the question of its relation with volcanic activities and/or concern about the occurrence of larger future earthquakes in the capital of Iran. Tehran megacity is surrounded by several inner-city and adjacent active faults that correspond to high-risk seismic sources in the area. The Mosha fault with ~150 km long is one of the major active faults in central Alborz and east of Tehran. It has hosted several historical earthquakes (i.e. 1665 Mw 6.5 and 1830 Mw 7.1 earthquakes) in the vicinity of the 2020 Mw 5.0 Tehran earthquake’s hypocenter. In this study, we evaluate the seismic sequence of the Tehran earthquake and obtain the full moment tensor inversion of this event and its larger aftershocks, which is a key tool to discriminate between tectonic and volcanic earthquakes. Furthermore, we obtain a robust characterization of the finite fault model of this event applying probabilistic earthquake source inversion framework using near-field strong-motion records and broadband seismograms, with an estimation of the uncertainties of source parameters. Due to the relatively weak magnitude and deeper centroid depth (~12 km), no static surface displacement was observed in the coseismic interferograms, and modeling performed by seismic records. Focal mechanism solution from waveform inversion, with a significant double-couple component, is compatible with the orientation of the sinistral north-dipping Mosha fault at the centroid location. The finite fault model suggests that the mainshock rupture propagated towards the northwest. This directivity enhanced the peak acceleration in the direction of rupture propagation, observed in strong-motion records. The 2020 moderate magnitude earthquake with 2 casualties, highlights the necessity of high-resolution seismic monitoring in the capital of Iran, which is exposed to a risk of destructive earthquakes with more than 10 million population. Our results are important for the hazard and risk assessment, and the forthcoming earthquake early warning system development in Tehran metropolis.</p>


2020 ◽  
Author(s):  
Raghukanth Stg ◽  
Varun K singla

<p><strong> </strong></p><p>In recent times, seismic rotational motions have received significant attention of seismologists as well earthquake engineers. The measurement of rotational motions is useful not only for the complete characterization of the ground motion, but also for estimating the additional risk to civil engineering structures posed by these motions. With the development of state-of-the-art instruments such as the ring-laser and fiber-optic gyroscopes, and mechanical and magneto-hydrodynamic devices, it is now possible to measure these motions over a large frequency bandwidth with accuracy. Nevertheless, nearly all earthquake-prone regions of the world lack the instrumentation to record these motions and the existing database of recorded rotational motions thus remains limited. In such situations analytical methods can provide estimates of rotational ground motions. It is thus common to simulate these motions and most analytical methods simulate rotations as curl of the displacement field. Because this relation (refer to as ‘curl-based’) is based on the ‘small deformations’ assumption, the ‘curl-based’ rotations are in a sense an approximation of the exact rotations, the latter being computed by using the rotational matrix directly extracted from the deformation gradient using the polar decomposition theorem. Despite this, no study has so far discussed the situations in which the ‘curl-based’ relation no longer holds. This paper therefore attempts to shed light on this issue by conducting two numerical studies. In the first study, a simple case of a point kinematic dislocation shear source buried in a homogeneous, elastic half-space is considered to get a qualitative idea about the combinations of the source and medium parameters that can lead to appreciable errors in the ‘curl-based’ rotations. The second study considers the finite-fault source model of Chi-Chi earthquake to illustrate these errors in a realistic scenario. Both these studies indicate that the ‘curl-based’ rotations can be in appreciable error when ‘exact’ rotations are in excess of  20 degrees and that this can happen in the near-field regions of surface-rupturing faults.</p>


1972 ◽  
Vol 62 (2) ◽  
pp. 561-589 ◽  
Author(s):  
Thomas C. Hanks ◽  
Max Wyss

abstract Teleseismic determinations of body-wave (P, S) spectra, interpreted in terms of the Brune (1970) seismic-source model, are used to estimate the parameters seismic moment (Mto) and source dimension (r) for three large, shallow, strike-slip earthquakes occurring on nearly vertical fault planes and for which the same parameters can be determined from field (F) data. These earthquakes are (1) the Borrego Mountain, California, earthquake (April 9, 1968) for which [Mo(P) = 10, Mo(S) = 6.6, and Mo(F) = 3.6] × 1025 dyne-cm and [r(p) = 14, r(S) = 23, and L/2(F) = 17] km; (2) the Mudurnu Valley, Turkey, earthquake (July 22, 1967) for which [-Mo(P) = 9.1, Mo(S) = 8.5, and Mo(F) = 7.4] × 1026 dyne-cm, and [r(P) = 39, r(S) = 48, and L/2(F) = 40] km; and (3) the Dasht-e-Bayāz, Iran, earthquake (August 31, 1968) for which [Mo(P) = 4.8, Mo(S) = 8.6, and Mo(F) = 18] × 1026 dyne-cm, and [r(P) = S1, r(S) = 48, and L/2(F) = 40] km. The Brune (1970) model is well-calibrated with respect to the determination of these parameters for the earthquakes considered. A minimum estimate for the radiated energy can be expressed in terms of Mo and r; this estimate is low by a factor of 10 with respect to the estimate obtained from energy-magnitude relations for these three earthquakes. The stress drops of these events are of the order of 10 bars.


2021 ◽  
Author(s):  
Josip Stipčević ◽  
Valerio Poggi ◽  
Marijan Herak ◽  
Stefano Parolai ◽  
Davorka Herak ◽  
...  

<p>The Department of Geophysics, University of Zagreb and the Italian National Institute of Oceanography and Applied Geophysics (OGS) installed on January 4th 2021, five temporary seismic stations near the town of Petrinja, Croatia, in the aftermath of  the 29 Decembre 2020 Mw 6.4 earthquake. The stations equipped with a seismometer and a strong motion sensor, recorded the aftershock sequence beginning six days after the mainshock allowing to augment the permanent seismic network in the area improving the azimuthal coverage and providing additional near‐field observations.</p><p>In this presentation we summarize the motivation and goals of the deployment; details regarding the station installation, instrumentation, and configurations and observations from the network. The collected data set will be useful for carrying out several seismological studies including the analysis of variability of strong ground motions in near field, the determination of the aftershocks source parameters,  the estimation (if any) of rupture directivity of small events, the clustering of events in space and time, the better imaging of the fault zone, the evolution of crustal properties within and outside of the fault zone.</p>


Sign in / Sign up

Export Citation Format

Share Document