Rupture Characteristics and Bedrock Structural Control of the 2016 Mw 6.0 Intraplate Earthquake in the Petermann Ranges, Australia

2020 ◽  
Vol 110 (3) ◽  
pp. 1037-1045 ◽  
Author(s):  
Januka Attanayake ◽  
Tamarah R. King ◽  
Mark C. Quigley ◽  
Gary Gibson ◽  
Dan Clark ◽  
...  

ABSTRACT The 20 May 2016 surface-rupturing intraplate earthquake in the Petermann Ranges is the largest onshore earthquake to occur in the Australian continent in 19 yr. We use in situ and Interferometric Synthetic Aperture Radar surface observations, aftershock distribution, and the fitting of P-wave source spectra to determine source properties of the Petermann earthquake. Surface observations reveal a 21-km-long surface rupture trace (strike=294°±29°) with heterogeneous vertical displacements (<0.1–0.96  m). Aftershock arrays suggest a triangular-shaped rupture plane (dip  ≈  30°) that intersects the subsurface projection of the major geophysical structure (Woodroffe thrust [WT]) proximal to the preferred location of the mainshock hypocenter, suggesting the mainshock nucleated at a fault junction. Footwall seismicity includes apparent southwest-dipping Riedel-type alignments, including possible activation of the deep segment of the WT. We estimate a moment magnitude (Mw) of 6.0 and a corner frequency (fc) of 0.2 Hz, respectively, from spectral fitting of source spectra in the 0.02–2 Hz frequency band. These translate into a fault area of 124  km2 and an average slip of 0.36 m. The estimated stress drop of 2.2 MPa is low for an intraplate earthquake; we attribute this to low-frictional slip (effective coefficient of friction >0.015) along rupture-parallel phyllosilicate-rich surfaces within the host rock fabric with possible additional contributions from elevated pore-fluid pressures.

2007 ◽  
Vol 4 (3) ◽  
pp. 173-182 ◽  
Author(s):  
Liu Yang ◽  
Zhang Qinghong ◽  
Bao Leiying ◽  
Wei Xiucheng
Keyword(s):  
P Wave ◽  
S Waves ◽  

1976 ◽  
Vol 66 (3) ◽  
pp. 803-825
Author(s):  
William A. Peppin

abstract Some 140 P-wave spectra of explosions, earthquakes, and explosion-induced aftershocks, all within the Nevada Test Site, have been computed from wide-band seismic data at close-in (< 30 km) and near-regional (200 to 300 km) distances. Observed near-regional corner frequencies indicate that source corner frequencies of explosions differ little from those of earthquakes of similar magnitude for 3 < ML < 5. Plots of 0.8 to 1.0 Hz Pg spectral amplitude versus 12-sec Rayleigh-wave amplitude show a linear trend with unit slope over three orders of magnitude for explosions; earthquakes fail to be distinguished from explosions on such a plot. These spectra also indicate similar source spectra for explosions in different media (tuff, alluvium, rhyolite) which corroborates Cherry et al. (1973). Close-in spectra of three large explosions indicate that: (1) source corner frequencies of explosions scale with yield in a way significantly different from previously published scaling laws; (2) explosion source spectra in tuff are flat from 0.2 to 1.0 Hz (no overshoot); (3) the far-field source spectrum decays at least as fast as frequency cubed. Taken together, these data indicate that the following factors are not responsible for Peppin and McEvilly's (1974) near-regional discriminant: (a) source dimension, (b) source rise time, or (c) shape of the source spectrum.


1983 ◽  
Vol 73 (6A) ◽  
pp. 1499-1511
Author(s):  
Paul Silver

Abstract A method is proposed for retrieving source-extent parameters from far-field body-wave data. At low frequency, the normalized P- or S-wave displacement amplitude spectrum can be approximated by |Ω^(r^,ω)| = 1 − τ2(r^)ω2/2 where r^ specifies a point on the focal sphere. For planar dislocation sources, τ2(r^) is linearly related to statistical measures of source dimension, source duration, and directivity. τ2(r^) can be measured as the curvature of |Ω^(r^,ω)| at ω = 0 or the variance of the pulse Ω^(r^,t). The quantity ωc=2τ−1(r^) is contrasted with the traditional corner frequency ω0, defined as the frequency at the intersection of the low- and high-frequency trends of |Ω^(r^,ω)|. For dislocation models without directivity, ωc(P) ≧ ωc(S) for any r^. A mean corner frequency defined by averaging τ2(r^) over the focal sphere, ω¯c=2<τ2(r^)>−1/2, satisfies ωc(P) > ωc(S) for any dislocation source. This behavior is not shared by ω0. It is shown that ω0 is most sensitive to critical times in the rupture history of the source, whereas ωc is determined by the basic parameters of source extent. Evidence is presented that ωc is the corner frequency measured on actual seismograms. Thus, the commonly observed corner frequency shift (P-wave corner greater than the S-wave corner), now viewed as a shift in ωc is simply a result of spatial finiteness and is expected to be a property of any dislocation source. As a result, the shift cannot be used as a criterion for rejecting particular dislocation models.


1990 ◽  
Vol 80 (5) ◽  
pp. 1205-1231
Author(s):  
Jiajun Zhang ◽  
Thorne Lay

Abstract Determination of shallow earthquake source mechanisms by inversion of long-period (150 to 300 sec) Rayleigh waves requires epicentral locations with greater accuracy than that provided by routine source locations of the National Earthquake Information Center (NEIC) and International Seismological Centre (ISC). The effects of epicentral mislocation on such inversions are examined using synthetic calculations as well as actual data for three large Mexican earthquakes. For Rayleigh waves of 150-sec period, an epicentral mislocation of 30 km introduces observed source spectra phase errors of 0.6 radian for stations at opposing azimuths along the source mislocation vector. This is larger than the 0.5-radian azimuthal variation of the phase spectra at the same period for a thrust fault with 15° dip and 24-km depth. The typical landward mislocation of routinely determined epicenters of shallow subduction zone earthquakes causes source moment tensor inversions of long-period Rayleigh waves to predict larger fault dip than indicated by teleseismic P-wave first-motion data. For dip-slip earthquakes, inversions of long-period Rayleigh waves that use an erroneous source location in the down-dip or along-strike directions of a nodal plane, overestimate the strike, dip, and slip of that nodal plane. Inversions of strike-slip earthquakes that utilize an erroneous location along the strike of a nodal plane overestimate the slip of that nodal plane, causing the second nodal plane to dip incorrectly in the direction opposite to the mislocation vector. The effects of epicentral mislocation for earthquakes with 45° dip-slip fault mechanisms are more severe than for events with other fault mechanisms. Existing earth model propagation corrections do not appear to be sufficiently accurate to routinely determine the optimal surface-wave source location without constraints from body-wave information, unless extensive direct path (R1) data are available or empirical path calibrations are performed. However, independent surface-wave and body-wave solutions can be remarkably consistent when the effects of epicentral mislocation are accounted for. This will allow simultaneous unconstrained body-wave and surface-wave inversions to be performed despite the well known difficulties of extracting the complete moment tensor of shallow sources from fundamental modes.


1974 ◽  
Vol 64 (4) ◽  
pp. 1159-1180 ◽  
Author(s):  
F. A. Dahlen

abstract We construct a theoretical three-dimensional kinematical model of shallow-focus earthquake faulting in order to investigate the ratio of the P- and S-wave corner frequencies of the far-field elastic radiation. We attempt to incorporate in this model all of the important gross kinematical features which would arise if ordinary mechanical friction should be the dominant traction resisting fault motion. These features include a self-similar nucleation at a single point, a subsonic spreading of rupture away from that point, and a termination of faulting by smooth deceleration. We show that the ratio of the P-wave corner frequency to the S-wave corner frequency for any model which has these features will be less than unity at all points on the focal sphere.


2021 ◽  
Author(s):  
Janneke de Laat ◽  
Sergei Lebedev ◽  
Bruna Chagas de Melo ◽  
Nicolas Celli ◽  
Raffaele Bonadio

<p>We present a new S-wave velocity tomographic model of the Australian Plate, Aus21.  It is constrained by waveforms of 0.9 million seismograms with both the corresponding sources and stations located within the half-hemisphere centred at the Australian continent. Waveform inversion extracts structural information from surface, S- and multiple S-waves on the seismograms in the form of a set of linear equations. These equations are then combined in a large linear system and inverted jointly to obtain a tomographic model of S- and P-wave speeds and S-wave azimuthal anisotropy of the crust and upper mantle. The model has been validated by resolution tests and, for particular locations in Australia with notable differences with previous models, by independent inter-station measurements of surface-wave phase velocities, which we performed using available array data. </p><p> </p><p>Aus21 offers new insights into the structure and evolution of the Australian Plate and its boundaries. The Australian cratonic lithosphere occupies nearly all of the western and central Australia but shows substantial lateral heterogeneity. It extends up to the northern edge of the plate, where it is colliding with island arcs, without subducting. The rugged eastern boundary of the cratonic lithosphere provides a lithospheric definition of the Tasman Line. The thin, warm lithosphere below the eastern part of the continent, east of the Tasman Line, underlies the Cenozoic volcanism locations in the area. The lithosphere is also thin and warm below much of the Tasman Sea, underlying the Lord Howe hotspot and the submerged part of western Zealandia. A low velocity anomaly that may indicate the single source of the Lord Howe and Tasmanid hotspots is observed in the transition zone offshore the Australian continent, possibly also sourcing the East Australia hotspot. Another potential hotspot source is identified below the Kermadec Trench, causing an apparent slab gap in the overlying slab and possibly related to the Samoa Hotspot to the north. Below a portion of the South East Indian Ridge (the southern boundary of the Australian Plate) a pronounced high velocity anomaly is present in the 200-400 km depth range just east of the Australian-Antarctic Discordance (AAD), probably linked to the evolution of this chaotic ridge system.</p>


Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1519-1527 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan

Reflected P‐to‐P and P‐to‐S converted seismic waves in a two‐component elastic common‐source gather generated with a P‐wave source in a two‐dimensional model can be imaged by two independent scalar reverse‐time depth migrations. The inputs to migration are pure P‐ and S‐waves that are extracted by divergence and curl calculations during (shallow) extrapolation of the elastic data recorded at the earth’s surface. For both P‐to‐P and P‐to‐S converted reflected waves, the imaging time at each point is the P‐wave traveltime from the source to that point. The extracted P‐wave is reverse‐time extrapolated and imaged with a P‐velocity model, using a finite difference solution of the scalar wave equation. The extracted S‐wave is reverse‐time extrapolated and imaged similarly, but with an S‐velocity model. Converted S‐wave data requires a polarity correction prior to migration to ensure constructive interference between data from adjacent sources. Synthetic examples show that the algorithm gives satisfactory results for laterally inhomogeneous models.


1989 ◽  
Vol 60 (3) ◽  
pp. 95-100 ◽  
Author(s):  
S.E. Hough ◽  
K. Jacob ◽  
R. Busby ◽  
P.A. Friberg

Abstract We present analysis of a magnitude 3.5 event which occurred at 9 km epicentral distance from a digital strong motion instrument operated by the National Center for Earthquake Engineering Research. Although the size of this isolated event is such that it can scarcely be considered to be a significant earthquake, a careful analysis of this high quality recording does yield several interesting results: 1) the S-wave spectra can be interpreted in terms of a simple omega-squared source spectrum and frequency-independent attenuation, 2) there is the suggestion of a poorly-resolved resonance in the P-wave spectrum, and perhaps most importantly, 3) the apparently simple S-wave spectra can be fit almost equally well with a surprisingly wide range of seismic corner frequencies, from roughly 5 to 25 Hz. This uncertainty in corner frequency translates into uncertainties in inferred Q values of almost an order of magnitude, and into uncertainties in stress drop of two orders of magnitude. Given the high quality of the data and the short epicentral distance to the station, we consider it likely that resolution of spectral decay and corner frequency will be at least as poor for any other recording of earthquakes with comparable or smaller magnitudes.


2020 ◽  
Author(s):  
Wenzheng Gong ◽  
Xiaofei Chen

<p>Spectra analysis is helpful to understand earthquake rupture processes and estimate source parameters like stress drop. Obtaining real source spectra and source time function isn’t easy, because the station recordings contain path effect and we usually can’t get precise path information. Empirical Green’s function (EGF) method is a popular way to cancel out the path effect, main two of which are the stacking spectra method (Prieto et al, 2006) and the spectral ratio method (Viegas et al, 2010; Imanishi et al, 2006). In our study, we apply the latter with multitaper spectral analysis method (Prieto et al, 2009) to calculate relative source spectra and relative source time function. Target event and EGFs must have similar focal mechanism and be collocated, so we combine correlation coefficient of wave at all stations and focal mechanism similarity to select proper EGFs.</p><p>The Bucaramanga nest has very high seismicity, so it’s suitable to calculate source spectra by using EGF method. We calculate the source spectra and source time function of about 1540 earthquakes (3-5.7ml, 135-160km depth) at Bucaramanga nest in Colombia. Simultaneously we also estimate corner frequency by fitting spectral source model (Brune, 1970; Boatwright, 1980) and stress drop using simple model (Eshelby, 1957) of earthquakes with multiple station recordings or EGFs. We obtain about 30000 events data with 12 stations from National Seismological Network of Colombia (RSNC).</p><p>The result show that the source spectra of most earthquakes fitted well by omega-square model are smooth, and the source spectra of some have obvious ‘holes’ near corner frequency, and the source time function of a few earthquakes appear two separate peeks. The first kind of earthquakes are style of self-arresting ruptures (Xu et al. 2015), which can be autonomously arrested by itself without any outside interference. Abercrombie (2014) and Wen et al. (2018) both researched the second kind of earthquakes and Wen think that this kind of earthquakes are style of the runaway ruptures including subshear and supershear ruptures. The last kind of earthquakes maybe be caused by simultaneous slip on two close rupture zone. Stress drop appear to slightly increase with depth and are very high (assuming rupture velocity/s wave velocity is 0.9). We also investigate the high-frequency falloff n, usually 2, of Brune model and Boatwright model by fitting all spectra, and find that the best value of n for Boatwright model is 2 and for Brune model is 3.5.</p>


Sign in / Sign up

Export Citation Format

Share Document