Fault modeling and stress drop estimation based on millimeter-scale tsunami records of an M6 earthquake detected by the dense and wide pressure gauge array

Author(s):  
Kubota Tatsuya ◽  
Tatsuhiko Saito ◽  
Wataru Suzuki

<p>Tsunamis observed by offshore ocean-bottom pressure gauges have been used to infer fault models and stress drops for major (M > 7) offshore earthquakes, to understand the earthquake and tsunamigenesis (e.g., Satake et al. 2013). However, it is challenging to observe tsunamis due to moderate (M ~6) earthquakes with reasonable quality by those, previous, few and remote pressure gauge arrays. Recently, a new, dense and wide pressure gauge network, the Seafloor Observation Network for Earthquakes and Tsunamis along the Japan Trench (S-net), was constructed off eastern Japan (Kanazawa et al. 2016). This array observed tsunamis associated with a moderate (M~6) earthquake which occurred inside the array, with amplitudes of less than one cm. We analyzed these millimeter-scale tsunami records to infer the finite fault model and stress drop, and to examine its relationship with other interplate earthquake phenomena.</p><p>We analyzed the pressure data associated with an Mw 6.0 earthquake off Sanriku on August 20, 2016. This earthquake was located at the shallowest part of the plate boundary off Sanriku, Japan, near the northern edge of the rupture area of the 1896 Sanriku tsunami earthquake (Kanamori, 1972). Although the signal-to-noise ratio is not high, the westward tsunami propagation with the velocity of ~0.1 km/s could be recognized when the waveforms were aligned according to the station locations. Using these data, we constrained the rectangular fault model with a uniform slip across the fault. As a result, the fault model was located ~10 km to the west of the Global CMT centroid (a seismic moment M<sub>0</sub>= 1.4 × 10<sup>18</sup> Nm, Mw 6.0, and a stress drop of Δσ = 1.5 MPa). The stress drop seems not so small as expected in tsunami earthquakes such as the 1896 Sanriku tsunami earthquake (≪ ~1 MPa, e.g., Kanamori 1972) even if the uncertainty of the stress drop estimation is considered (Δσ > ~ 0.7 MPa). We also found the rupture area was unlikely to overlap with regions where slow earthquakes are active, such as low-frequency-tremors and very-low-frequency-earthquakes (e.g., Matsuzawa et al. 2015; Nishikawa et al. 2019; Tanaka et al. 2019).</p><p>This result demonstrates that the S-net new dense and wide pressure gauge array dramatically increases the detectability of a millimeter-scale tsunami and the constraints on earthquake source parameters of moderate earthquakes off eastern Japan. It is expected that more tsunamis due to minor-to-moderate offshore earthquakes are recorded by this new array, which will reveal the spatial variation of the stress drops, or mechanical properties, along the plate interface with much higher resolution than previously possible.</p>

2021 ◽  
Vol 40 (9) ◽  
pp. 667-676
Author(s):  
Clay Kurison ◽  
Huseyin S. Kuleli

Microseismic events associated with shale reservoir hydraulic fracturing stimulation (HFS) are interpreted to be reactivations of ubiquitous natural fractures (NFs). Despite adoption of discrete fracture network (DFN) models, accounting for NFs in fluid flow within shale reservoirs has remained a challenge. For an explicit account of NFs, this study introduced the use of seismology-based relations linking seismic moment, moment magnitude, fault rupture area, and stress drop. Microseismic data from HFS monitoring of Marcellus Shale horizontal wells had been used to derive planar hydraulic fracture geometry and source properties. The former was integrated with associated well production data found to exhibit transient linear flow. Analytical solutions led to linear flow parameters (LFPs) and system permeability for scenarios depicting flow through infinite and finite conductivity hydraulic fractures. Published core plug permeability was stress-corrected for in-situ conditions to estimate average matrix permeability. For comparison, the burial and thermal history for the study area was used in 1D Darcy-based modeling of steady and episodic expulsion of petroleum to account for geologic timescale persistence of abnormal pore pressure. Both evaluations resulted in matrix permeability in the same picodarcy (pD) range. Coupled with LFPs, reactivated NF surface area for stochastic DFNs was estimated. Subsequently, the aforementioned seismology-based relations were used for determining average stress drops needed to estimate NF rupture area matching flow-based DFN surface areas. Stress drops, comparable to values for tectonic events, were excluded. One of the determined values matched stress drops for HFS operations in past and recent seismological studies. In addition, calculated changes in pore pressure matched estimates in the aforementioned studies. This study unlocked the full potential of microseismic data beyond extraction of planar geometry attributes and stimulated reservoir volume (SRV). Here, microseismic events were explicitly used in the quantitative account of NFs in fluid flow within shale reservoirs.


1991 ◽  
Vol 81 (2) ◽  
pp. 553-575 ◽  
Author(s):  
Michael Fehler ◽  
W. Scott Phillips

Abstract An inversion that fits spectra of earthquake waveforms and gives robust estimates of corner frequency and low-frequency spectral amplitude has been used to determine source parameters of 223 microearthquakes induced by hydraulic fracturing in granodiorite. Assuming a ω−2 source model, the inversion fits the P-wave spectra of microearthquake waveforms to determine individual values of corner frequency and low-frequency spectral amplitude for each event and one average frequency-independent Q for all source-receiver paths. We also implemented a constraint that stress drops of all microearthquakes be similar but not equal and found that this constraint did not significantly degrade the quality of the fits to the spectra. The waveforms analyzed were recorded by a borehole seismometer. The P-wave Q was found to be 1070. For Q values as low as 600 and as high as 3000, the misfit between model and spectra increased by less than 5 per cent and the average corner frequency changed by less than 15 per cent from those obtained with a Q of 1070. Average stress drop was 3.7 bars. Seismic moments obtained from spectra ranged from 1013 to 1018 dyne-cm. The low stress drops are interpreted to result from underestimation of the actual stress drops because of a nonuniform distribution of stress drop and slip along the fault planes. Spatially varying stress drops and slips result from the strong rock heterogeneity due to the injection of fluid into the rock. Stress drops were found to be larger near the edges of the seismic zone, in regions that had not been seismically active during previous injections. The seismic moments determined from spectra were used to obtain a coda length-to-moment relation. Then, moments were estimated for 1149 events from measurements of coda lengths from events whose moments could not be measured from spectra because of saturation or a low signal-to-noise ratio. The constant of proportionality between cumulative number of events and seismic moment is higher than that found for tectonic regions. The slope is so high that the seismic energy release is dominated by the large number of small events. In the absence of information about the number of events smaller than we studied, we cannot estimate the total seismic energy released by the hydraulic injection.


1984 ◽  
Vol 74 (1) ◽  
pp. 27-40
Author(s):  
M. E. O'Neill

Abstract Source dimensions and stress drops of 30 small Parkfield, California, earthquakes with coda duration magnitudes between 1.2 and 3.9 have been estimated from measurements on short-period velocity-transducer seismograms. Times from the initial onset to the first zero crossing, corrected for attenuation and instrument response, have been interpreted in terms of a circular source model in which rupture expands radially outward from a point until it stops abruptly at radius a. For each earthquake, duration magnitude MD gave an estimate of seismic moment MO and MO and a together gave an estimate of static stress drop. All 30 earthquakes are located on a 6-km-long segment of the San Andreas fault at a depth range of about 8 to 13 km. Source radius systemically increases with magnitude from about 70 m for events near MD 1.4 to about 600 m for an event of MD 3.9. Static stress drop ranges from about 2 to 30 bars and is not strongly correlated with magnitude. Static stress drop does appear to be spatially dependent; the earthquakes with stress drops greater than 20 bars are concentrated in a small region close to the hypocenter of the magnitude 512 1966 Parkfield earthquake.


1979 ◽  
Vol 36 (10) ◽  
pp. 1223-1227
Author(s):  
D. D. Lemon ◽  
P. H. LeBlond ◽  
T. R. Osborn

Seiche motions observed in San Juan Harbour with a bottom-mounted pressure gauge have been Fourier-analyzed and interpreted in terms of a theoretical model of oscillations in a rectangular basin with an exponential depth profile. Two of the observed periods (at 14.6 and 38.5 min) are identified with resonances of the basin; two other significant low frequency peaks (at 21 and 55 min) do not coincide with resonant periods of the basin and must be due to strong external forcing. Higher frequency fluctuations (20–160 s) are attributed to swell and to its subharmonic interactions with edge waves. Key words: water waves, seiches, mathematical model, Juan de Fuca Strait, British Columbia


1995 ◽  
Vol 85 (5) ◽  
pp. 1513-1517
Author(s):  
Z.-M. Yin ◽  
G. C. Rogers

Abstract Earthquake faulting results in stress drop over the rupture area. Because the stress drop is only in the shear stress and there is no or little stress drop in the normal stress on the fault, the principal stress directions must rotate to adapt such a change of the state of stress. Using two constraints, i.e., the normal stress on the fault and the vertical stress (the overburden pressure), which do not change before and after the earthquake, we derive simple expressions for the rotation angle in the σ1 axis. For a dip-slip earthquake, the rotation angle is only a function of the stress-drop ratio (defined as the ratio of the stress drop to the initial shear stress) and the angle between the σ1 axis and the fault plane, but for a strike-slip earthquake the rotation angle is also a function of the stress ratio. Depending on the faulting regimes, the σ1 axis can either rotate toward the direction of fault normal or rotate away from the direction of fault normal. The rotation of the stress field has several important seismological implications. It may play a significant role in the generation of heterogeneous stresses and in the occurrence and distribution of aftershocks. The rotation angle can be used to estimate the stress-drop ratio, which has been a long-lasting topic of debate in seismology.


2020 ◽  
Vol 91 (4) ◽  
pp. 2320-2329 ◽  
Author(s):  
James S. Neely ◽  
Seth Stein ◽  
Bruce D. Spencer

Abstract Earthquake stress drop, the stress change on a fault due to an earthquake, is important for seismic hazard analysis because it controls the level of high-frequency ground motions that damage structures. Numerous studies report that stress drops vary by tectonic environment, providing insight into a region’s seismic hazard. Here, we show that teleseismic stress-drop estimates have large uncertainties that make it challenging to distinguish differences between the stress drops of different earthquakes. We compared stress drops for ∼900 earthquakes derived from two independent studies using teleseismic data and found practically zero correlation. Estimates for the same earthquake can differ by orders of magnitude. Therefore, reported stress-drop differences between earthquakes may not reflect true differences. As a result of these larger uncertainties, some tectonic environment stress-drop patterns that appear in one study do not appear in the other analysis of the same earthquakes. These large uncertainties in teleseismic estimates might lead to erroneous inferences about earthquake hazards. In many applications, it may be more appropriate to assume that earthquakes in different regions have approximately the same average stress drop.


2020 ◽  
Vol 110 (1) ◽  
pp. 241-249
Author(s):  
Kazuhito Hikima ◽  
Akihiro Shimmura

ABSTRACT For crustal earthquakes, the scaling relationship between the seismic moment M0 and rupture area S varies with the size of the earthquake, due to the limited thickness of the seismogenic layer. In those M0–S scaling relations, in most cases, the calculated static stress drop is altered with the size of earthquake, although the change depends on the assumed fault model. However, it is not clear whether the dependence of the stress drop on M0 is physically reasonable. In this study, the scaling relation between M0 and S, which assumes a constant stress drop over a wide M0 range, is discussed based on the analytical stress drop formula of a rectangular strike-slip fault. In the proposed relation, M0 is proportional to S3/2 for small and medium faults and to S1 for long faults. In addition, the relation between M0 and S varies in the intermediate range, depending on the aspect ratio. The scaling relation showed good agreement with past event data when the saturated rupture width was set to around 15–20 km and the stress drop was set to about 3 MPa.


2020 ◽  
Author(s):  
Jonas Folesky ◽  
Joern Kummerow ◽  
Serge A. Shapiro

<p>The Northern Chilean subduction zone has been monitored by the IPOC network for more than ten years. During this time period two very large earthquakes occurred, the 2007 M<sub>W</sub>7.7  Tocopilla earthquake and the 2014 M<sub>W</sub>8.1 Iquique earthquake. Over the entire subduction zone a vast amount of seismic activity has been recorded and a huge catalog was compiled including over 100000 events (Sippl et al. 2018). With this exceptional data base we attempt a systematic analysis of the stress drops of as many events from the catalog as possible. We apply different estimation techniques, namely the spectral ratio type, the spectral stacking approach, and the lower bound method. A goal of our research is a comparison and possibly a combination of the techniques to obtain reliable and well constrained results.</p><p>The data set covers events at the interface, within the subducting plate, crustal events, and intermediate depth events. It therefore bears a great potential to better understand the stress drop distribution within a subduction zone. Also, the long observation interval allows to analyze temporal variations according to pre-, inter-, and post-seismic phases of megathrust earthquakes.   </p><p>We present preliminary results where a subset of 730 events with a magnitude range of M<sub>L</sub>2.7 - M<sub>L</sub>4.8  was used for analysis with the spectral ratio technique. For these events we show maps of spatial stress drop variation, and we analyze the time dependent stress drop variance. </p>


Author(s):  
Seong Ju Jeong ◽  
Brian W. Stump ◽  
Heather R. DeShon ◽  
Louis Quinones

ABSTRACT Earthquakes in the Fort Worth basin (FWB) have been induced by the disposal of recovered wastewater associated with extraction of unconventional gas since 2008. Four of the larger felt earthquakes, each on different faults, prompted deployment of local distance seismic stations and recordings from these four sequences are used to estimate the kinematic source characteristics. Source spectra and the associated source parameters, including corner frequency, seismic moment, and stress drop, are estimated using a modified generalized inversion technique (GIT). As an assessment of the validity of the modified GIT approach, corner frequencies and stress drops from the GIT are compared to estimates using the traditional empirical Green’s function (EGF) method for 14 target events. For these events, corner-frequency residuals (GIT−EGF) have a mean of −0.31 Hz, with a standard deviation of 1.30 Hz. We find consistent mean stress drops using the GIT and EGF methods, 9.56 and 11.50 MPa, respectively, for the common set of target events. The GIT mean stress drop for all 79 earthquakes is 5.33 MPa, similar to estimates for global intraplate earthquakes (1–10 MPa) as well as other estimates for induced earthquakes near the study area (1.7–9.5 MPa). Stress drops exhibit no spatial or temporal correlations or depth dependency. In addition, there are no time or space correlations between estimated FWB stress drops and modeled pore-pressure perturbations. We conclude that induced earthquakes in the FWB occurring on normal faults in the crystalline basement release pre-existing tectonic stresses and that stress drops on the four sequences targeted in this study do not directly reflect perturbations in pore-fluid pressure on the fault.


1983 ◽  
Vol 73 (6A) ◽  
pp. 1735-1751
Author(s):  
J. B. Fletcher ◽  
J. Boatwright ◽  
W. B. Joyner

Abstract Three estimates of stress differences, which include Brune stress drop, stress drop from rms of acceleration (arms), and the apparent stress, have been determined for 13 earthquakes at Monticello, South Carolina, a site of reservoir-induced seismicity. Data for nine of the events come from digitally recorded three-component seismograms at four or five stations that were deployed around the Monticello Reservoir in May and early June 1979. The data from the other four events come from a strong-motion accelerograph located on the dam abutment at the southwest end of the reservoir. Estimates of the seismic moment (Mo) range from 4.6 × 1017 to 3.4 × 1020 dyne-cm (S waves) and radiated energy from about 1011 to 3 × 1016 dyne-cm (S waves). Brune stress drops ranged from 0.5 bars to about 90 bars and show a strong dependence on depth (focal depths range from 0.07 to 1.4 km) and a moderate dependence on Mo. Arms stress drops from the direct S-wave span a similar range of values and also exhibit a strong dependence on depth. Apparent stresses are usually lower than the other estimates of stress differences by a factor of 2 to 4. Seismic stress differences are highest in the topmost 0.2 to 0.3 km, a depth range for which the in situ measurements of stress and pore pressure suggest that the rock is in a state of incipient failure. In this depth range, where the four largest events occurred, the stress drops are of the same order as the ambient shear stress. These data suggest that at Monticello, where pore fluids have a strong influence on the failure process, the largest stresses released seismically are in regions most conducive to failure and that the seismic efficiencies for events at Monticello are larger than have been reported for other tremors in different tectonic settings.


Sign in / Sign up

Export Citation Format

Share Document