Application of stochastic fractal surface rupture on non-planar faults in tsunami simulation

2020 ◽  
Author(s):  
Shane Murphy ◽  
Andrè Herrero ◽  
Fabrizio Romano ◽  
Stefano Lorito
Keyword(s):  
Fault Plane ◽  
Velocity Model ◽  
Source Model ◽  
Tsunami Source ◽  
Surface Rupture ◽  
Tsunami Simulation ◽  
Model Complex ◽  
Source Models ◽  
Fault Trace ◽  
Subduction Zone Earthquakes

<p>Non-planar faults and surface reached rupture are seldom considered in the source modelling of subduction zone earthquakes. Here we present a new method for accounting for both phenomena in the generation of stochastic slip distribution while still maintaining self similar properties. To do this, we use the composite source model, which involves the placement of numerous circular dislocations on the fault plane. The fault plane is described by an unstructured mesh allowing for a non-planar surface while surface rupture is correctly accounted for by reflecting the slip from circular dislocations that intersect with the fault trace.</p><p>In a case study we demonstrate that the inclusion of rupture at the surface alters the ground or seafloor deformation both in terms of the magnitude (between 60%-20% in 5km zone near the fault trace) and the orientation of the deformation vectors (i.e. by up to 5 degrees). Such changes can have a significant effect on tsunami source and subsequent wave.</p><p>Additionally, with a prescribed rupture velocity model, complex source time functions can also be calculated for each element on the fault plane. Generally, rise time is assumed to be instantaneous in tsunami simulation.</p><p>We will also present preliminary results focused on comparing the tsunami wave height observed along nearby coastlines generated by the different source models (i.e. with/without surface reached rupture and variable source time functions).    </p>

Numerical Simulation of Tsunami Run-Up and Inundation for the 2011 Tōhuku-Oki Tsunami: A Parametric Analysis for Tsunami Run-Up and Wave Height

2014 ◽  
Author(s):  
Debashis Basu ◽  
Robert Sewell ◽  
Kaushik Das ◽  
Ron Janetzke ◽  
Biswajit Dasgupta ◽  
...  
Keyword(s):  
Wave Propagation ◽  
Wave Height ◽  
Seismic Waves ◽  
Source Model ◽  
Tsunami Source ◽  
Computational Framework ◽  
Wave Amplification ◽  
Source Models ◽  
Run Up ◽  
Tsunami Run Up

This paper presents computational results for predicting earthquake-generated tsunami from a developed integrated computational framework. The computational framework encompasses the entire spectrum of modeling the earthquake-generated tsunami source, open-sea wave propagation, and wave run-up including inundation and on-shore effects. The present work develops a simplified source model based on pertinent local geologic and tectonic processes, observed seismic data (i.e., data obtained by inversion of seismic waves from seismographic measurements), and geodetic data (i.e., directly measured seafloor and land deformations). These source models estimated configurations of seafloor deformation used as initial waveforms in tsunami simulations. Together with sufficiently accurate and resolved bathymetric and topographic data, they provided the inputs needed to numerically simulate tsunami wave propagation, inundation and coastal impact. The present work systematically analyzes the effect of the tsunami source model on predicted tsunami behavior and the associated variability for the 2011 Tōhuku-Oki tsunami. Simulations were carried out for the 2011 Tōhuku -Oki Tsunami that took place on March 11, 2011, from an MW 9.1 earthquake. The numerical simulations were performed using the fully nonlinear Boussinesq hydrodynamics code, FUNWAVE-TVD (distributed by the University of Delaware). In addition, a sensitivity analysis was also carried out to study the effect of earthquake magnitude on the predicted wave height. The effect of coastal structure on the wave amplification at the shore is also studied. Simulated tsunami results for wave heights are compared to the available observational data from GPS (Global Positioning System) at the central Miyagi location.

Field surveys and numerical modeling of the December 2018 Anak Krakatau volcanic tsunami

2020 ◽  
Author(s):  
Mohammad Heidarzadeh ◽  
Purna Sulastya Putra ◽  
Abdul Muhari ◽  
Septriono Hari Nugroho
Keyword(s):  
Numerical Modeling ◽  
Source Model ◽  
Tsunami Source ◽  
Tsunami Height ◽  
Tsunami Inundation ◽  
Field Surveys ◽  
Source Models ◽  
Krakatau Volcano ◽  
Source Dimension ◽  
Maximum Tsunami Height

<p>We report results of field surveys and numerical modeling of the tsunami generated by the Anak Krakatau volcano eruption on 22 December 2018. We conducted two sets of field surveys of the coastal areas destroyed by the Anak Krakatau tsunami in 26-30 December 2018 and 4-10 January 2020. Field surveys provided information about the maximum tsunami height as well as the most damaged area. The maximum tsunami height was up to 13 m. Most locations registered a wave height of 3-4 m. Tsunami inundation was limited to approximately 100 m. For modeling, we considered 12 source models and conducted numerical modeling. The scenarios have source dimensions of 1.5–4 km and initial tsunami amplitudes of 10–200 m. By comparing observed and simulated waveforms, we constrained the tsunami source dimension and initial amplitude in the ranges of 1.5–2.5 km and 100–150 m, respectively. The best source model involves potential energy of 7.14 × 10<sup>13</sup>–1.05 × 10<sup>14</sup> J which is equivalent to an earthquake of magnitude 6.0–6.1.</p>

OBSERVATION AND SIMULATION OF TSUNAMIS INDUCED BY THE 2003 TOKACHI-OKI EARTHQUAKE (M 8.0)

2015 ◽  
Vol 2 (3) ◽  
Author(s):  
Tatsuo Ohmachi ◽  
Shusaku Inoue ◽  
Tetsuji Imai
Keyword(s):  
Rayleigh Wave ◽  
Water Pressure ◽  
Near Field ◽  
Source Region ◽  
Tsunami Source ◽  
Band Pass Filter ◽  
Pass Filter ◽  
Tsunami Simulation ◽  
Sea Bed ◽  
Band Pass

The 2003 Tokachi-oki earthquake (MJ 8.0) occurred off the southeastern coast of Tokachi, Japan, and generated a large tsunami which arrived at Tokachi Harbor at 04:56 with a wave height of 4.3 m. Japan Marine Science and Technology Center (JAMSTEC) recovered records of water pressure and sea-bed acceleration at the bottom of the tsunami source region. These records are first introduced with some findings from Fourier analysis and band-pass filter analysis. Water pressure disturbance lasted for over 30 minutes and the duration was longer than those of accelerations. Predominant periods of the pressure looked like those excited by Rayleigh waves. Next, numerical simulation was conducted using the dynamic tsunami simulation technique able to represent generation and propagation of Rayleigh wave and tsunami, with a satisfactory result showing validity and usefulness of this technique. Keywords: Earthquake, Rayleigh wave, tsunami, near-field

A dynamic source model for the San Fernando earthquake

1978 ◽  
Vol 68 (6) ◽  
pp. 1555-1576
Author(s):  
Michel Bouchon
Keyword(s):  
Fault Plane ◽  
High Stress ◽  
Time History ◽  
High Strength ◽  
High Amplitude ◽  
Tectonic Stress ◽  
Source Model ◽  
Dynamic Crack ◽  
San Fernando ◽  
High Level

abstract We model the San Fernando earthquake as a propagating rupture in a half-space, using for the slip-time-history on the fault plane analytical expressions which approximate the slip functions of dynamic crack models obtained by Das and Aki (1977a, b). We synthesize the strong ground motions and accelerations at the Pacoima Dam site and compute the teleseismic signals for different models of cracks. Three major featuras of the data–the strong pulse associated with the beginning of the rupture, the high acceleration phase on the Pacoima Dam records, and the presence of ripples on the teleseismic seismograms–which are not compatible with a smooth rupture process, are well explained by a crack with barriers model where the rupture encounters, along the fault plane, barriers or obstacles of high strength materials which may remain unbroken after the passage of the rupture front. A high-stress drop (400 to 500 bars) is required in the hypocentral area to explain the high-amplitude short-duration first pulse of the teleseismic records. This indicates a high level of tectonic stress in the area. A study of the earthquake series following the main shock shows that the aftershocks which took place in the region where major slipping occurred during the earthquake may represent the release of some of the barriers.

scholarly journals Multilayer modelling of waves generated by explosive submarine volcanism

2021 ◽  
Author(s):  
Matthew W. Hayward ◽  
Colin N. Whittaker ◽  
Emily M. Lane ◽  
William Power ◽  
Stéphane Popinet ◽  
...  
Keyword(s):  
Underwater Explosion ◽  
Theoretical Models ◽  
Mono Lake ◽  
Tsunami Source ◽  
Computationally Efficient ◽  
Relative Significance ◽  
Submarine Volcanism ◽  
Wave Characteristics ◽  
Wide Range ◽  
Source Models

Abstract. Theoretical source models of underwater explosions are often applied in studying tsunami hazards associated with submarine volcanism; however, their use in numerical codes based on the shallow water equations can neglect the significant dispersion of the generated wavefield. A non-hydrostatic multilayer method is validated against a laboratory-scale experiment of wave generation from instantaneous disturbances and at field-scale submarine explosions at Mono Lake, California, utilising the relevant theoretical models. The numerical method accurately reproduces the range of observed wave characteristics for positive disturbances and suggests a previously unreported relationship of extended initial troughs for negative disturbances at low dispersivity and high nonlinearity parameters. Satisfactory amplitudes and phase velocities within the initial wave group are found using underwater explosion models at Mono Lake. The scheme is then applied to modelling tsunamis generated by volcanic explosions at Lake Taupō, New Zealand, for a magnitude range representing ejecta volumes between 0.04–0.4 km3. Waves reach all shores within 15 minutes with maximum incident crest amplitudes around 4 m at shores near the source. This work shows that the multilayer scheme used is computationally efficient and able to capture a wide range of wave characteristics, including dispersive effects, which is necessary when investigating submarine explosions. This research therefore provides the foundation for future studies involving a rigorous probabilistic hazard assessment to quantify the risks and relative significance of this tsunami source mechanism.

scholarly journals Unsupervised Domain Adaptation by Statistics Alignment for Deep Sleep Staging Networks

2021 ◽  
Author(s):  
Jiahao Fan ◽  
Hangyu Zhu ◽  
Xinyu Jiang ◽  
Long Meng ◽  
Cong Fu ◽  
...  
Keyword(s):  
Large Scale ◽  
Domain Adaptation ◽  
Source Model ◽  
Deep Sleep ◽  
Sleep Staging ◽  
Target Domain ◽  
Source Domain ◽  
Unsupervised Domain Adaptation ◽  
Generalization Problem ◽  
Source Models

Deep sleep staging networks have reached top performance on large-scale datasets. However, these models perform poorer when training and testing on small sleep cohorts due to data inefficiency. Transferring well-trained models from large-scale datasets (source domain) to small sleep cohorts (target domain) is a promising solution but still remains challenging due to the domain-shift issue. In this work, an unsupervised domain adaptation approach, domain statistics alignment (DSA), is developed to bridge the gap between the data distribution of source and target domains. DSA adapts the source models on the target domain by modulating the domain-specific statistics of deep features stored in the Batch Normalization (BN) layers. Furthermore, we have extended DSA by introducing cross-domain statistics in each BN layer to perform DSA adaptively (AdaDSA). The proposed methods merely need the well-trained source model without access to the source data, which may be proprietary and inaccessible. DSA and AdaDSA are universally applicable to various deep sleep staging networks that have BN layers. We have validated the proposed methods by extensive experiments on two state-of-the-art deep sleep staging networks, DeepSleepNet+ and U-time. The performance was evaluated by conducting various transfer tasks on six sleep databases, including two large-scale databases, MASS and SHHS, as the source domain, four small sleep databases as the target domain. Thereinto, clinical sleep records acquired in Huashan Hospital, Shanghai, were used. The results show that both DSA and AdaDSA could significantly improve the performance of source models on target domains, providing novel insights into the domain generalization problem in sleep staging tasks.<br>

Average Q and yield estimates from the Pahute Mesa test site

1987 ◽  
Vol 77 (4) ◽  
pp. 1274-1294
Author(s):  
R. W. Burger ◽  
T. Lay ◽  
L. J. Burdick
Keyword(s):  
Time Domain ◽  
Near Field ◽  
Energy Levels ◽  
Source Model ◽  
Spectral Shape ◽  
Field Amplitude ◽  
Far Field ◽  
Source Models ◽  
Attenuation Models ◽  
Shape Data

Abstract Attenuation models, with and without frequency dependence, have been developed through analysis of time-domain amplitude measurements and teleseismic spectral shape data from Pahute Mesa nuclear explosions. The time-domain analysis is based on a near-field to far-field amplitude comparison. The near-field amplitude information is incorporated in two parameterized explosion source models (Mueller-Murphy and Helmberger-Hadley) based on analyses of near-field data. The teleseismic amplitude observations are from a large data set of WWSSN short-period analog recordings. For the narrow-band time-domain data, the various source and attenuation models are indistinguishable. We utilize the spectral shape data in the 0.5- to 4-Hz band as a constraint on the source-attenuation models at higher frequencies, concluding that either source model, when convolved with the appropriate frequency-dependent Q model, can be consistent with both the near-field and far-field time-domain amplitudes and the spectral shape data. Given the trade-off between source and attenuation models and the similarity of the different source models in the 0.5- to 4-Hz band, it is difficult to prefer clearly one source model over the other. The Mueller-Murphy model is more consistent with surface wave amplitude measurements because of larger predicted long-period energy levels. Whether or not frequency dependence is included in the attenuation model, the value of t* near 1 Hz is about 1.0 sec (assuming the Mueller-Murphy source model) or 0.8 sec (assuming the Helmberger-Hadley source model). This 0.2 sec difference results from greater 1-Hz energy levels for the Mueller-Murphy source model. Adopting an average attenuation model, predicted amplitudes and yields are shown to be within the uncertainty of the data for all the events analyzed.

UNCERTAINTY ASSESSMENT OF TSUNAMI HEIGHT FOR FUTURE NANKAI-TONANKAI EARTHQUAKES USING STOCHASTIC TSUNAMI SOURCE MODELS

2021 ◽  
Vol 77 (2) ◽  
pp. I_181-I_186
Author(s):  
Takuya MIYASHITA ◽  
Kazuki KURATA ◽  
Tomohiro YASUDA ◽  
Nobuhito MORI ◽  
Tomoya SHIMURA
Keyword(s):  
Uncertainty Assessment ◽  
Tsunami Source ◽  
Tsunami Height ◽  
Source Models

scholarly journals Surface rupture in stochastic slip models

2020 ◽  
Vol 221 (2) ◽  
pp. 1081-1089 ◽  
Author(s):  
S Murphy ◽  
A Herrero
Keyword(s):  
Probability Density Function ◽  
Probability Density ◽  
Density Function ◽  
Source Model ◽  
Surface Rupture ◽  
Size Frequency Distribution ◽  
Fault Surface ◽  
Planar Fault ◽  
Small Change ◽  
Composite Source Model

SUMMARY As an alternative to spectral methods, stochastic self-similar slip can be produced through a composite source model by placing a power-law scaling size-frequency distribution of circular slip dislocations on a fault surface. However these models do not accurately account for observed surface rupture behaviour. We propose a modification to the composite source model that corrects this issue. The advantage of this technique is that it accommodates the use of fractal slip distributions on non-planar fault surfaces. However to mimic a surface rupture using this technique, releasing the boundary condition at the top of the fault, we observed a systematic decrease in slip at shallow depths. We propose a new strategy whereby the surface is treated like a reflector with the slip being folded back onto the fault. Two different techniques based on this principal are presented: the first is the method of images. It requires a small change to pre-existing codes and works for planar faults. The second involves the use of a multistage trilateration technique. It is applied to non-planar faults described by an unstructured mesh. The reflected slip calculated using the two techniques is near identical on a planar fault, suggesting they are equivalent. Applying this correction, where reflected slip is accounted for in the composite source model, the lack of slip at shallow depths is not observed any more and there is no systematic trend with depth. However, there are other parameters which may affect the spatial distribution of slip across the fault plane. For example, the type of probability density function used in the placement of the subevent is also important. In the case where the location of maximum slip is known to a first order, a Gaussian may be appropriate to describe the probability function. For hazard assessment studies a uniform probability density function is more suitable as it provides no underlying systematic spatial trend.

Sign in / Sign up

Export Citation Format

Share Document