A Study on the Largest Hydraulic-Fracturing-Induced Earthquake in Canada: Observations and Static Stress-Drop Estimation

ABSTRACT On 17 August 2015, an Mw 4.6 earthquake occurred northwest of Fort St. John, British Columbia, possibly induced by hydraulic fracturing (HF). We use data from eight broadband seismometers located ∼50 km from the hypocenter to detect and estimate source parameters of more than 300 events proximal to the mainshock. Stress-drop values estimated using seismic moment and corner frequency from single-event spectra and spectral ratios range from ∼1 to 35 MPa, within the typical range of tectonic earthquakes. We observe an ∼5-day delay between the onset of fluid injection and the mainshock, a b-value of 0.78 for the sequence, and a maximum earthquake magnitude larger than the prediction based on the total injection volume, suggesting that the Mw 4.6 sequence occurred on a pre-existing fault and that the maximum magnitude is likely controlled by tectonic conditions. Results presented here show that pre-existing fault structures should be taken into consideration to better estimate seismic hazard associated with HF operations and to develop schemes for risk mitigation in close proximity to HF wells.

Geometric controls on megathrust earthquakes

SUMMARY The role of subduction zone geometry in the nucleation and propagation of great-sized earthquake ruptures is an important topic for earthquake hazard, since knowing how big an earthquake can be on a given fault is fundamentally important. Past studies have shown subducting bathymetric features (e.g. ridges, fracture zones, seamount chains) may arrest a propagating rupture. Other studies have correlated the occurrence of great-sized earthquakes with flat megathrusts and homogenous stresses over large distances. It remains unclear, however, how subduction zone geometry and the potential for great-sized earthquakes (M 8+) are quantifiably linked—or indeed whether they can be. Here, we examine the potential role of subduction zone geometry in limiting earthquake rupture by mapping the planarity of seismogenic zones in the Slab2 subduction zone geometry database. We build from the observation that historical great-sized earthquakes have preferentially occurred where the surrounding megathrust is broadly planar, and we use this relationship to search for geometrically similar features elsewhere in subduction zones worldwide. Assuming geometry exerts a primary control on earthquake propagation and termination, we estimate the potential size distribution of large (M 7+) earthquakes and the maximum earthquake magnitude along global subduction faults based on geometrical features alone. Our results suggest that most subduction zones are capable of hosting great-sized earthquakes over much of their area. Many bathymetric features previously identified as barriers are indistinguishable from the surrounding megathrust from the perspective of slab curvature, meaning that they either do not play an important role in arresting earthquake rupture or that their influence on slab geometry at depth is not resolvable at the spatial scale of our subduction zone geometry models.

ESTIMATION OF MAXIMUM EARTHQUAKE MAGNITUDE OF EARTHQUAKE POTENTIALS FOR YOGYAKARTA DEPRESSION AREA, INDONESIA

Maximum magnitudes of earthquake potentials are estimated for Yogyakarta depression area by using the faultlength and earthquake magnitude relations for fault specific seismic sources. For estimation of maximum earthquake magnitude, the fault specific seismic sources are modeled as 18 normal faults and 6 strike-slip faults sources referring the geological map of McDonald, 1984 and Rihardjo et al., 1995. For the present area the subduction zone earthquakes are expected to happen in the offshore region regarding the study on the seismicity of the region and the focal mechanisms of the past earthquakes. So three area sources are also assumed for this region and the possible maximum earthquake magnitudes for these sources are determined by probabilistic approaches.