Fracture Nucleation due to Occurrence of Cracks in a Seismic Source Region

1992 ◽  
Vol 45 (1) ◽  
pp. 107-117
Author(s):  
Hideo YUKUTAKE
Geophysics ◽  
2002 ◽  
Vol 67 (1) ◽  
pp. 156-165 ◽  
Author(s):  
Graham. J. Hicks

In finite‐difference methods a seismic source can be implemented using either initial wavefield values or body forces. However, body forces can only be specified at finite‐difference nodes, and, if using initial values, a source cannot be located close to a reflecting boundary or interface in the model. Hence, difficulties can exist with these schemes when the region surrounding a source is heterogeneous or when a source either is positioned between nodes or is arbitrarily close to a free surface. A completely general solution to these problems can be obtained by using Kaiser windowed sinc functions to define a small region around the true source location that contains several nodal body forces. Both monopole and dipole point sources can be defined, enabling many source types to be implemented in either acoustic or elastic media. Such a function can also be used to arbitrarily locate receivers. If the number of finite‐difference nodes per wavelength is four or more (and with a source region half‐width of only four nodes) this scheme results in insignificant phase errors and in amplitude errors of no more than 0.1%. Numerical examples for sources located less than one node from either a free surface or an image source demonstrate that the scheme can be used successfully for any surface‐source or multisource configuration.


2021 ◽  
Author(s):  
Gaspard Farge ◽  
Claude Jaupart ◽  
Nikolaï Shapiro

<p>Many subduction zones host intermittent, low-frequency, low-magnitude seismic activity emitted from the vicinity of the plates' interface. For instance, in Guerrero, Mexico, deep (30--50 km) low-frequency earthquakes (LFEs) occur in bursts, and migrate in cascades along the subduction interface. Those patterns are often attributed to episodic pulses of fluid pressure and slow slip that travel within the fault zone. However, the dynamic behavior of the permeable system in which fluid-pressure circulates remains a blindspot in most models of tremor generation, even as geological observations report pervasive imprint of strong, localized fluid pressure and permeability variations in its source region.</p><p>In order to analyze the role of such processes in generating tremor, we design a simple model of how fluid pressure and permeability can interact within the subduction interface, and generate realistic, tremor-like patterns. It is based on seismic source triggering and interaction in a permeable channel. The latter contains a number of low-permeability plugs acting as elementary fault-valves. In a mechanism akin to erosive burst documented in porous media, valve permeability abruptly opens and closes in response to the local fluid pressure. The brutal pressure transient and/or mechanical fracturing associated with valve opening acts as the seismic source of an LFE-like event. The strong fluid pressure transient that it triggers allows valves to interact constructively: as a valve breaks open, neighbor valves are more likely to break. This interaction therefore leads to cascades and migrations of synthetic seismicity along the model fault channel, that can synchronize into larger bursts of activity that migrate more slowly along the channel. In our model, valve activity draws patterns of that closely resemble tremor patterns in Guerrero and other subduction zones.</p><p>The input metamorphic fluid flux at the base of the channel exerts a key control on the occurence of and distribution of synthetic tremor in space and time. A weak input flux will not allow valves to open, conversely a strong flux will not allow them to close. In both cases, no activity will occur. However when the value of the fluid flux is intermediate, permanent regimes of sustained activity arise. Depending on its value, activity can be strongly time-clustered, quasi-periodic or random but constant in time.</p><p>Our model is based on a simple yet powerful and realistic description of the permeability and its dynamics in fault zones. It allows for new interpretations of low-frequency seismicity in terms of effective flux and fault-zone permeability, both for long-term regimes and finer scale, transient dynamics. Eventually, it could lead to deep enhancements of our understanding of fault-zone hydraulic processes and how they are coupled with fault-slip.</p>


1982 ◽  
Vol 72 (6A) ◽  
pp. 1825-1841
Author(s):  
A. McGarr

abstract Two independent arguments indicate an upper bound of about 10 for the ratio ro/ri in the expressions for peak velocity v and peak acceleration a at close hypocentral distances R: v = (βΔτro/μR)[0.10(ro/ri) + 0.15] and a = (Δτ/ρR)[0.30(ro/ri)2 + 0.45], where ri is the radius of the most heavily loaded asperity that fails within an earthquake source region of radius ro, Δτ is the stress drop, β is the shear-wave velocity, μ is the modulus of rigidity, and ρ is the density; these relationships are for ground motion recorded in a whole-space. First, a recently reported data set was augmented by observations for six earthquakes in the magnitude range 4 ≦ ML ≦ 6.6, for which ground motion was recorded at a minimum of five sites at hypocentral distances of the order of 10 km; the new events include the 1979 Coyote Lake and 1979 Imperial Valley shocks. The entire data set of 22 events, spanning a range in seismic moment from 5 × 1016 to over 1026 dyne-cm, is consistent both with the bound ro/ri < 10 and with the previous conclusion that this ratio does not depend systematically on earthquake size. Second, a theoretical argument, using the result of Savage and Wood that the apparent stress acting on the earthquake fault plane is less than half of the stress drop, is made to the effect that ro/ri < 10. In addition, absolute limits, independent of earthquake size, for peak acceleration are related to the state of stress in the crust; for an extensional state of stress a ≦ 0.40 g and for a compressional stress state a ≦ 2.0 g, where a now represents the maximum horizontal acceleration as recorded at the surface directly above the seismic source.


2014 ◽  
Vol 30 (1) ◽  
pp. 1-15 ◽  
Author(s):  
Brendon A. Bradley ◽  
Mark C. Quigley ◽  
Russ J. Van Dissen ◽  
Nicola J. Litchfield

This paper provides an overview of the ground motion and seismic source aspects of the Canterbury earthquake sequence. Common reported attributes among the largest earthquakes in this sequence are complex ruptures, large displacements per unit fault length, and high stress drops. The Darfield earthquake produced an approximately 30 km surface rupture in the Canterbury Plains with dextral surface displacements of several meters, and a subordinate amount of vertical displacement, impacting residential structures, agricultural land, and river channels. The dense set of strong ground motions recorded in the near-source region of all the major events in the sequence provides significant insight into the spatial variability in ground motion characteristics, as well as the significance of directivity, basin-generated surface waves, and nonlinear local site effects. The ground motion amplitudes in the 22 February 2011 earthquake, in particular, produced horizontal ground motion amplitudes in the Central Business District (CBD) well above those specified for the design of conventional structures.


2019 ◽  
Vol 2 (2) ◽  
pp. 192-199
Author(s):  
Tuğba Türker ◽  
Yusuf Bayrak

The aim of this study determined spatial-temporal variations of Gutenberg-Richter b-values using the Classic Method of Moment (CMM) and seismicity in and around Lake Van. We wrote a Matlab program for this method in the Matlab system. The faults for the Van region used as updated in the GIS system given the reference from the different sources by Bayrak and Türker (will be used in the doctoral thesis). Around Lake Van divided into 4 different seismic source regions. For each seismic source region, probabilities ranging from 0 to 1, the annual probability of exceedances and return periods were calculated for the next 25, 50 and 100 years. Van region determined the seismicity b value changes, β-values. We mapped spatial variations of the b-values using the GIS system. As a result, Lake Van region (region 2) occurred 7.1 magnitude an earthquake in the next 100 years with 31 % probability levels. If it occurred in the 50 years with 15 % probability levels, in the 25 years with 7.8 % probability levels. We estimated the low b-values in the Lake Van so, it has been the high seismicity and can be occurred an earthquake as the high magnitude in the next years. This study will be lead to earthquake hazard analysis and seismic hazard studies in and around Lake Van.


2008 ◽  
Vol 14 ◽  
pp. 105-116 ◽  
Author(s):  
A. B. Rabinovich ◽  
L. I. Lobkovsky ◽  
I. V. Fine ◽  
R. E. Thomson ◽  
T. N. Ivelskaya ◽  
...  

Abstract. Two major earthquakes near the Central Kuril Islands (Mw=8.3 on 15 November 2006 and Mw=8.1 on 13 January 2007) generated trans-oceanic tsunamis recorded over the entire Pacific Ocean. The strongest oscillations, exceeding several meters, occurred near the source region of the Kuril Islands. Tide gauge records for both tsunamis have been thoroughly examined and numerical models of the events have been constructed. The models of the 2006 and 2007 events include two important advancements in the simulation of seismically generated tsunamis: (a) the use of the finite failure source models by Ji (2006, 2007) which provide more detailed information than conventional models on spatial displacements in the source areas and which avoid uncertainties in source extent; and (b) the use of the three-dimensional Laplace equation to reconstruct the initial tsunami sea surface elevation (avoiding the usual shallow-water approximation). The close agreement of our simulated results with the observed tsunami waveforms at the open-ocean DART stations support the validity of this approach. Observational and model findings reveal that energy fluxes of the tsunami waves from the source areas were mainly directed southeastward toward the Hawaiian Islands, with relatively little energy propagation into the Sea of Okhotsk. A marked feature of both tsunamis was their high-frequency content, with typical wave periods ranging from 2–3 to 15–20 min. Despite certain similarities, the two tsunamis were essentially different and had opposite polarity: the leading wave of the November 2006 trans-oceanic tsunami was positive, while that for the January 2007 trans-oceanic tsunami was negative. Numerical modeling of both tsunamis indicates that, due to differences in their seismic source properties, the 2006 tsunami was more wide-spread but less focused than the 2007 tsunami.


Author(s):  
Hiroyuki Kumagai ◽  
Pablo Placios ◽  
Mario Ruiz ◽  
Hugo Yepes ◽  
Tomofumi Kozono

Sign in / Sign up

Export Citation Format

Share Document