scholarly journals Self-similar stochastic slip distributions on a non-planar fault for tsunami scenarios for megathrust earthquakes

2020 ◽  
Author(s):  
Masaru Nakano ◽  
Shane Murphy ◽  
Ryoichiro Agata ◽  
Yasuhiko Igarashi ◽  
Masato Okada ◽  
...  

Abstract Megathrust earthquakes that occur repeatedly along the plate interface of subduction zones can cause severe damage due to strong ground motion and the destructive tsunamis they can generate. We developed a set of scenario earthquakes to evaluate tsunami hazards and tsunami early warning systems for such devastating earthquakes. Although it is known that the slip distribution on a fault strongly affects the tsunami height distribution in near-field coastal areas, the slip distribution of future earthquakes cannot be exactly predicted. One way to resolve this difficulty is to create a set of scenario earthquakes in which a set of heterogeneous slip distributions on the source fault is stochastically generated based on a given slip probability density function (SPDF). The slip distributions generated in this manner differ from event to event, but their average over a large ensemble of models converges to a predefined SPDF resembling the long-term average of ruptures on the target fault zone. We created a set of SPDF-based scenario earthquakes for an expected future M w 8.2 Tonankai earthquake in the eastern half of the Nankai trough, off southwest Japan, and computed the ensuing tsunamis. We found that the estimated peak coastal amplitudes among the ensemble of tsunamis along the near-field coast differed by factors of 3 to 9 and the earliest and latest arrivals at each observation site differed by 400 to 700 s. The variations in both peak tsunami amplitude and arrival time at each site were well approximated by a Gaussian distribution. For cases in which the slip distribution is unknown, the average and standard deviation of these scenario datasets can provide first approximations of forecast tsunami height and arrival time and their uncertainties, respectively. At most coastal observation sites, tsunamis modelled similarly but using a uniform slip distribution underpredicted tsunami amplitudes but gave earlier arrival times than those modeled with a heterogeneous slip distribution. Use of these earlier arrival times may be useful for providing conservative early warnings of tsunami arrivals. Therefore, tsunami computations for both heterogeneous and uniform slip distributions are important for tsunami disaster mitigation.

2020 ◽  
Author(s):  
Masaru Nakano ◽  
Shane Murphy ◽  
Ryoichiro Agata ◽  
Yasuhiko Igarashi ◽  
Masato Okada ◽  
...  

Abstract Megathrust earthquakes that occur repeatedly along the plate interface of subduction zones can cause severe damage due to strong ground motion and the destructive tsunamis they can generate. We developed a set of scenario earthquakes to evaluate tsunami hazards and tsunami early warning systems for such devastating earthquakes. Although it is known that the slip distribution on a fault strongly affects the tsunami height distribution in near-field coastal areas, the slip distribution of future earthquakes cannot be exactly predicted. One way to resolve this difficulty is to create a set of scenario earthquakes in which a set of heterogeneous slip distributions on the source fault is stochastically generated based on a given slip probability density function (SPDF). The slip distributions generated in this manner differ from event to event, but their average over a large ensemble of models converges to a predefined SPDF resembling the long-term average of ruptures on the target fault zone. We created a set of SPDF-based scenario earthquakes for an expected future Mw 8.2 Tonankai earthquake in the eastern half of the Nankai trough, off southwest Japan, and computed the ensuing tsunamis. We found that the estimated peak coastal amplitudes among the ensemble of tsunamis along the near-field coast differed by factors of 3 to 9 and the earliest and latest arrivals at each observation site differed by 400 to 700 s. The variations in both peak tsunami amplitude and arrival time at each site were well approximated by a Gaussian distribution. For cases in which the slip distribution is unknown, the average and standard deviation of these scenario datasets can provide first approximations of forecast tsunami height and arrival time and their uncertainties, respectively. At most coastal observation sites, tsunamis modelled similarly but using a uniform slip distribution underpredicted tsunami amplitudes but gave earlier arrival times than those modeled with a heterogeneous slip distribution. Use of these earlier arrival times may be useful for providing conservative early warnings of tsunami arrivals. Therefore, tsunami computations for both heterogeneous and uniform slip distributions are important for tsunami disaster mitigation.


2012 ◽  
Vol 7 (1) ◽  
pp. 19-25 ◽  
Author(s):  
Yuichiro Tanioka ◽  
◽  
Hamzah Latief ◽  
Haris Sunendar ◽  
Aditya Riadi Gusman ◽  
...  

Several large earthquakes have recently occurred along the Sumatra-Java subduction zone, the 2004 great Sumatra-Andaman earthquake, the 2005 great Nias earthquake, the 2006 West Java tsunami earthquake, the 2007 great Bengkulu earthquake, and the 2010Mentawai tsunami earthquakes. Serious tsunami disasters were caused by the great underthrust earthquakes which ruptured the plate interface near the trench such as the 2004 Sumatra-Andaman, 2006West Java, 2010Mentawai earthquakes. At Palabuhanratu, maximum tsunami height distribution and inundation areas were computed from expected fault models located near the Java trench. The results shows that the most populated areas of Palabuhanratu would be severely damaged by the expected tsunami caused by the fault model of Mw 8.5. After discussing tsunami disaster mitigation measures with the local government, the result of tsunami inundation area in this study were used to decide tsunami evacuation areas and evacuation routes. The local government also installed tsunami evacuation sign boards near the coast.


2010 ◽  
Vol 10 (6) ◽  
pp. 1281-1292 ◽  
Author(s):  
B. Poisson ◽  
R. Pedreros

Abstract. Two historical landslide-induced tsunamis that reached the coasts of the French Lesser Antilles are studied. First, the Martinique coast was hit by a tsunami down the western flank of Montagne Pelée at the beginning of the big eruption of May 1902. More recently, the northeastern coast of Guadeloupe was affected by a tsunami that had been generated around Montserrat by pyroclastic flows entering the sea, during the July 2003 eruption of the Soufrière Hills volcano. We use a modified version of the GEOWAVE model to compute numerical simulations of both events. Two source hypotheses are considered for each tsunami. The comparison of the simulation results with reported tsunami height data helps to discriminate between the tested source decriptions. In the Martinique case, we obtain a better fit to data when considering three successive lahars entering the sea, as a simplified single source leads to an overstimation of the tsunami wave heights at the coast. In the Montserrat case, the best model uses a unique source which volume corresponds to published data concerning the peak volume flow. These findings emphasize the importance of an accurate description of the relevant volume as well as the timing sequence of the source event in landslide-generated tsunami modelling. They also show that considering far-field effects in addition to near-field effects may significantly improve tsunami modelling.


2021 ◽  
Vol 8 (1) ◽  
Author(s):  
Hijrah Saputra ◽  
Wahyudi Wahyudi ◽  
Iman Suardi ◽  
Ade Anggraini ◽  
Wiwit Suryanto

AbstractThis study comprehensively investigates the source mechanisms associated with the mainshock and aftershocks of the Mw = 6.3 Yogyakarta earthquake which occurred on May 27, 2006. The process involved using moment tensor inversion to determine the fault plane parameters and joint inversion which were further applied to understand the spatial and temporal slip distributions during the earthquake. Moreover, coseismal slip distribution was overlaid with the relocated aftershock distribution to determine the stress field variations around the tectonic area. Meanwhile, the moment tensor inversion made use of near-field data and its Green’s function was calculated using the extended reflectivity method while the joint inversion used near-field and teleseismic body wave data which were computed using the Kikuchi and Kanamori methods. These data were filtered through a trial-and-error method using a bandpass filter with frequency pairs and velocity models from several previous studies. Furthermore, the Akaike Bayesian Information Criterion (ABIC) method was applied to obtain more stable inversion results and different fault types were discovered. Strike–slip and dip-normal were recorded for the mainshock and similar types were recorded for the 8th aftershock while the 9th and 16th June were strike slips. However, the fault slip distribution from the joint inversion showed two asperities. The maximum slip was 0.78 m with the first asperity observed at 10 km south/north of the mainshock hypocenter. The source parameters discovered include total seismic moment M0 = 0.4311E + 19 (Nm) or Mw = 6.4 with a depth of 12 km and a duration of 28 s. The slip distribution overlaid with the aftershock distribution showed the tendency of the aftershock to occur around the asperities zone while a normal oblique focus mechanism was found using the joint inversion.


2021 ◽  
Author(s):  
Alberto Armigliato ◽  
Martina Zanetti ◽  
Stefano Tinti ◽  
Filippo Zaniboni ◽  
Glauco Gallotti ◽  
...  

<p>It is well known that for earthquake-generated tsunamis impacting near-field coastlines the focal mechanism, the position of the fault with respect to the coastline and the on fault slip distribution are key factors in determining the efficiency of the generation process and the distribution of the maximum run-up and inundation along the nearby coasts. The time needed to obtain the aforementioned information from the analysis of seismic records is usually too long compared to the time required to issue a timely tsunami warning/alert to the nearest coastlines. In the context of tsunami early warning systems, a big challenge is hence to be able to define 1) the relative position of the hypocenter and of the fault and 2) the earthquake focal mechanism, based only on the preliminary earthquake localization and magnitude estimation, which are made available by seismic networks soon after the earthquake occurs.</p><p>In this study, the intrinsic unpredictability of the position of the hypocenter on the fault plane is studied through a probabilistic approach based on the analysis of two finite fault model datasets (SRCMOD and USGS) and by limiting the analysis to moderate-to-large shallow earthquakes (Mw  6 and depth  50 km). After a proper homogenization procedure needed to define a common geometry for all samples in the two datasets, the hypocentral positions are fitted with different probability density functions (PDFs) separately in the along-dip and along-strike directions.</p><p>Regarding the focal mechanism determination, different approaches have been tested: the most successful is restricted to subduction-type earthquakes. It defines average values and uncertainties for strike, dip and rake angles based on a combination of a proper zonation of the main tsunamigenic subduction areas worldwide and of subduction zone geometries available from publicdatabases.</p><p>The general workflow that we propose can be schematically outlined as follows. Once an earthquake occurs and the magnitude and hypocentral solutions are made available by seismic networks, it is possible to assign the focal mechanism by selecting the characteristic values for strike, dip and rake of the zone where the hypocenter falls into. Fault length and width, as well as the slip distribution on the fault plane, are computed through regression laws against magnitude proposed by previous studies. The resulting rectangular fault plane can be discretized into a matrix of subfaults: the position of the center of each subfault can be considered as a “realization” of the hypocenter position, which can then be assigned a probability. In this way, we can define a number of earthquake fault scenarios, each of which is assigned a probability, and we can run tsunami numerical simulations for each scenario to quantify the classical observables, such as water elevation time series in selected offshore/coastal tide-gauges, flow depth, run-up, inundation distance. The final results can be provided as probabilistic distributions of the different observables.</p><p>The general approach, which is still in a proof-of-concept stage, is applied to the 16 September 2015 Illapel (Chile) tsunamigenic earthquake (Mw = 8.2). The comparison with the available tsunami observations is discussed with special attention devoted to the early-warning perspective.</p>


2019 ◽  
pp. 121-127
Author(s):  
Victoria Erofeeva ◽  
Vasilisa Galyamina ◽  
Kseniya Gonta ◽  
Anna Leonova ◽  
Oleg Granichin ◽  
...  

In this paper we consider the problem of ultrasound tomography. Recently, an increased interest in ultrasound tomography has been caused by non-invasiveness of the method and increased detection accuracy (as compared to radiation tomography), and also ultrasound tomography does not put at risk human health. We study possibilities of detection of specific areas and determining their density using ultrasound tomography data. The process of image reconstruction based on ultrasound data is computationally complex and time consuming. It contains the following parts: calculation of the time-of-flight (TOF) of a signal, detection of specific areas, calculation of density of specific areas. The calculation of the arrival time of a signal is a very important part, because the errors in the calculation of quantities strongly influence the total problem solution. We offer ultrasound imaging reconstruction technology that can be easily parallelized. The whole process is described: from extracting the arrival times of signals raw data feeding from physical receivers to obtaining the desired results.


2018 ◽  
Vol 55 (4) ◽  
pp. 1272-1286 ◽  
Author(s):  
Kei Noba ◽  
José-Luis Pérez ◽  
Kazutoshi Yamazaki ◽  
Kouji Yano

Abstract De Finetti’s optimal dividend problem has recently been extended to the case when dividend payments can be made only at Poisson arrival times. In this paper we consider the version with bail-outs where the surplus must be nonnegative uniformly in time. For a general spectrally negative Lévy model, we show the optimality of a Parisian-classical reflection strategy that pays the excess above a given barrier at each Poisson arrival time and also reflects from below at 0 in the classical sense.


Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. KS63-KS73
Author(s):  
Yangyang Ma ◽  
Congcong Yuan ◽  
Jie Zhang

We have applied the cross double-difference (CDD) method to simultaneously determine the microseismic event locations and five Thomsen parameters in vertically layered transversely isotropic media using data from a single vertical monitoring well. Different from the double-difference (DD) method, the CDD method uses the cross-traveltime difference between the S-wave arrival time of one event and the P-wave arrival time of another event. The CDD method can improve the accuracy of the absolute locations and maintain the accuracy of the relative locations because it contains more absolute information than the DD method. We calculate the arrival times of the qP, qSV, and SH waves with a horizontal slowness shooting algorithm. The sensitivities of the arrival times with respect to the five Thomsen parameters are derived using the slowness components. The derivations are analytical, without any weak anisotropic approximation. The input data include the cross-differential traveltimes and absolute arrival times, providing better constraints on the anisotropic parameters and event locations. The synthetic example indicates that the method can produce better event locations and anisotropic velocity model. We apply this method to the field data set acquired from a single vertical monitoring well during a hydraulic fracturing process. We further validate the anisotropic velocity model and microseismic event locations by comparing the modeled and observed waveforms. The observed S-wave splitting also supports the inverted anisotropic results.


Sign in / Sign up

Export Citation Format

Share Document