Why Do Aftershocks Occur?

Total Page:16

File Type:pdf, Size:1020Kb

Why Do Aftershocks Occur? Yukutake and Iio Earth, Planets and Space (2017) 69:68 DOI 10.1186/s40623-017-0650-2 FULL PAPER Open Access Why do aftershocks occur? Relationship between mainshock rupture and aftershock sequence based on highly resolved hypocenter and focal mechanism distributions Yohei Yukutake1* and Yoshihisa Iio2 Abstract In order to clarify the origin of aftershocks, we precisely analyze the hypocenters and focal mechanisms of the after- shocks following the 2000 Western Tottori Earthquake, which occurred in the western part of Japan, using data from dense seismic observations. We investigate whether aftershocks occur on the mainshock fault plane on which coseis- mic slip occurred or they represent the rupture of fractures surrounding the mainshock fault plane. Based on the hypocenter distribution of the aftershocks, the subsurface fault structure of the mainshock is estimated using princi- pal component analysis. As a result, we can obtain the detail fault structure composed of 8 best-ft planes. We demon- strate that the aftershocks around the mainshock fault are distributed within zones of 1.0–1.5 km in thicknesses, and their focal mechanisms are signifcantly diverse. This result suggests that most of the aftershocks represent the rupture of fractures surrounding the mainshock fault rather than the rerupture of the mainshock fault. The aftershocks have a much wider zone compared with the exhumed fault zone in feld observations, suggesting that many aftershocks occur outside the fault damage zone. We fnd that most aftershocks except in and around the large-slip region are well explained by coseismic stress changes. These results suggest that the thickness of the aftershock distribution may be controlled by the stress changes caused by the heterogeneous slip distribution during the mainshock. The aftershock is also distributed within a much wider zone than the hypocenter distribution observed in swarm activity in the geothermal region, which is thought to be caused by the migration of hydrothermal fuid. This result implies a diference in generation processes: Stress changes due to the mainshock contribute primarily to the occurrence of aftershocks, whereas earthquake swarms in the geothermal region are caused by fuid migration within the localized zone. Keywords: Aftershock, Fault plane of mainshock, Fault damage zone, Hypocenter distribution, Focal mechanism Background of the mainshock slip (e.g., Das and Henry 2003) and to Understanding why aftershocks occur is one of the major be triggered by strength reductions due to high pore fuid unsolved problems in seismology, and several genera- pressure caused by the mainshock slip (e.g., Miller et al. tion mechanisms of aftershocks have been proposed. For 2004). Although all of these mechanisms have reasonable example, aftershocks have been thought to be reruptures physical bases, which mechanism is actually at work has of remnants on the fault surface that were not broken by rarely been confrmed based on precise aftershock data. the mainshock, to be generated by stress concentrations Te hypocenter distribution of aftershocks has been recognized to refect the structure of the fault plane that was ruptured during the mainshock. Tus, a number of *Correspondence: [email protected] previous studies estimated the subsurface structure of the 1 Hot Springs Research Institute of Kanagawa Prefectural Government, Odawara, Japan mainshock fault plane based on the aftershock distribu- Full list of author information is available at the end of the article tion (e.g., Shibutani et al. 2005a; Kato 2005). On the other © The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Yukutake and Iio Earth, Planets and Space (2017) 69:68 Page 2 of 15 hand, feld observation of an exhumed outcrop has dem- was released on the fault plane in the southern part of onstrated that a fault damage zone is developed around the aftershock region, whereas afterslip occurred on the a fault core that is the localized slip zone and contains fault planes in the northern part of the aftershock region. high-strain products, such as gouges, cataclasites, and Fukuyama et al. (2003) estimated the fault model of the ultracataclasites (Mitchell and Faulkner 2009). Te fault Western Tottori Earthquake using the detailed distribu- damage zone has also been identifed using geophysical tion of the aftershock hypocenters and their centroid data, such as seismic travel times or trapped waves (e.g., moment tensor (CMT) solutions. However, they used Cochran et al. 2009; Ben-Zion et al. 2003). We have not only seismic data obtained by permanent stations in clarifed what an aftershock distribution represents, spe- and around the aftershock region with a station inter- cifcally, whether aftershocks refect the rerupture of the val of approximately 20 km. Terefore, the locations of mainshock fault plane or the rupture of the damage zone the aftershock hypocenters in their study contain large surrounding the mainshock fault plane or whether after- uncertainty. shocks occur outside the damage zone. Tis is important to clarifying the generation mechanism of aftershocks. Data and methods Liu et al. (2003) investigated the hypocenter catalog Determination of hypocenters and focal mechanisms of the aftershock sequence following the 1993 Mw 7.3 A total of 48 temporary stations were installed in and Landers earthquake and concluded that only 5% of after- around the aftershock region during the period from 13 shocks occurred on the mainshock fault plane. However, October to early December 2000 (Shibutani et al. 2005b). since the location errors of the aftershock hypocenters in Together with data from seven permanent stations, the their study were up to 1 km, discussion based on more dense seismic observation network was developed with precisely determined hypocenters is essential. Moreo- an average station spacing of 4–5 km near the center of ver, a focal mechanism solution of aftershocks provides the aftershock region (Fig. 1). We manually picked the additional information associated with the orientation of arrival time of the P- and S-waves and the polarities of the fault plane. We analyzed the aftershock sequence of the P-waves for 4770 local events that occurred during the 2000 Western Tottori Earthquake that occurred in an the period from October 15, to November 31, 2000. Te intraplate region. A dense seismic observation network local magnitudes of these events ranged from 0.4 to 3.9. composed of 48 temporary stations was installed imme- In order to obtain the initial hypocenter locations, we diately after the mainshock. Te high-quality observation selected 4190 events for which the P- or S- wave arrival data provide an excellent opportunity to clarify the above times were obtained at eight or more stations. Te initial issue. hypocenter locations were obtained by applying the one- In the present study, we estimate the hypocenter loca- dimensional velocity model and the station corrections tions and focal mechanisms of the aftershocks of the estimated by Shibutani et al. (2005b) to the method of 2000 Western Tottori Earthquake using dense seismic hypocenter determination developed by Kawanishi et al. observations. We examine the fault model of the main- (2009), which is a modifed version of the Hypomh pro- shock based on precisely determined hypocenter loca- gram code (Hirata and Matsu’ura 1987). tions. Moreover, we discuss the relationship between the Te double-diference relocation (DD) method (Wald- subsurface fault structure of the mainshock and the after- hauser and Ellsworth 2000) was applied to the ini- shock distribution. tial hypocenter data. Te diferential arrival time data for the manually picked P- and S-waves consisted of Overview of the 2000 Western Tottori Earthquake 597,744 and 270,359 pairs, respectively. We also used Te 2000 Western Tottori Earthquake occurred on Octo- the diferential arrival times obtained by cross-corre- ber 6, 2000, in the western part of Tottori Prefecture in lation analysis. Te correlation analysis was conducted southwestern Honshu, Japan (Fig. 1). Te moment mag- using a velocity waveform with a time window of 0.5 s, nitude was determined to be 6.6 (Fukuyama et al. 2003). including the manually picked arrival times, applying a Iwata and Sekiguchi (2002) estimated the slip distribution band-pass flter between 3 and 20 Hz. We used only the of the mainshock using near-feld strong motion wave- diferential arrival time data with normalized cross-cor- forms and geodetic data and found that coseismic slip relation coefcients greater than or equal to 0.8. Using occurred on the southern part of the aftershock region, this threshold, we obtained cross-correlation data con- while there was little moment release on the north- sisting of 307,447 pairs for the P-wave and 177,677 pairs ern part of the aftershock region (Fig. 2b). Sagiya et al. for the S-wave. After application of the DD method, (2002) also estimated the slip distribution using global 4089 events were successfully relocated (Fig. 2). Te navigation satellite system (GNSS) data and reported root-mean-square (RMS) of the double-diference time that most of the seismic moment during the mainshock residual decreased from 95 to 30 ms for the manually Yukutake and Iio Earth, Planets and Space (2017) 69:68 Page 3 of 15 0 km 30 35.6˚N 35.4˚ 35.2˚ Mainshock Temporary station a First motion b CMT Permanent station Active fault 133.0˚E 133.2˚ 133.4˚ 133.6˚ Fig. 1 Map showing the spatial distribution of seismic stations (blue: temporary stations, white: permanent stations). The red star indicates the start- ing point of rupture during the 2000 Western Tottori Earthquake.
Recommended publications
  • Coulomb Stresses Imparted by the 25 March
    LETTER Earth Planets Space, 60, 1041–1046, 2008 Coulomb stresses imparted by the 25 March 2007 Mw=6.6 Noto-Hanto, Japan, earthquake explain its ‘butterfly’ distribution of aftershocks and suggest a heightened seismic hazard Shinji Toda Active Fault Research Center, Geological Survey of Japan, National Institute of Advanced Industrial Science and Technology (AIST), site 7, 1-1-1 Higashi Tsukuba, Ibaraki 305-8567, Japan (Received June 26, 2007; Revised November 17, 2007; Accepted November 22, 2007; Online published November 7, 2008) The well-recorded aftershocks and well-determined source model of the Noto Hanto earthquake provide an excellent opportunity to examine earthquake triggering associated with a blind thrust event. The aftershock zone rapidly expanded into a ‘butterfly pattern’ predicted by static Coulomb stress transfer associated with thrust faulting. We found that abundant aftershocks occurred where the static Coulomb stress increased by more than 0.5 bars, while few shocks occurred in the stress shadow calculated to extend northwest and southeast of the Noto Hanto rupture. To explore the three-dimensional distribution of the observed aftershocks and the calculated stress imparted by the mainshock, we further resolved Coulomb stress changes on the nodal planes of all aftershocks for which focal mechanisms are available. About 75% of the possible faults associated with the moderate-sized aftershocks were calculated to have been brought closer to failure by the mainshock, with the correlation best for low apparent fault friction. Our interpretation is that most of the aftershocks struck on the steeply dipping source fault and on a conjugate northwest-dipping reverse fault contiguous with the source fault.
    [Show full text]
  • Intraplate Earthquakes in North China
    5 Intraplate earthquakes in North China mian liu, hui wang, jiyang ye, and cheng jia Abstract North China, or geologically the North China Block (NCB), is one of the most active intracontinental seismic regions in the world. More than 100 large (M > 6) earthquakes have occurred here since 23 BC, including the 1556 Huax- ian earthquake (M 8.3), the deadliest one in human history with a death toll of 830,000, and the 1976 Tangshan earthquake (M 7.8) which killed 250,000 people. The cause of active crustal deformation and earthquakes in North China remains uncertain. The NCB is part of the Archean Sino-Korean craton; ther- mal rejuvenation of the craton during the Mesozoic and early Cenozoic caused widespread extension and volcanism in the eastern part of the NCB. Today, this region is characterized by a thin lithosphere, low seismic velocity in the upper mantle, and a low and flat topography. The western part of the NCB consists of the Ordos Plateau, a relic of the craton with a thick lithosphere and little inter- nal deformation and seismicity, and the surrounding rift zones of concentrated earthquakes. The spatial pattern of the present-day crustal strain rates based on GPS data is comparable to that of the total seismic moment release over the past 2,000 years, but the comparison breaks down when using shorter time windows for seismic moment release. The Chinese catalog shows long-distance roaming of large earthquakes between widespread fault systems, such that no M ࣙ 7.0 events ruptured twice on the same fault segment during the past 2,000 years.
    [Show full text]
  • Seismic Rate Variations Prior to the 2010 Maule, Chile MW 8.8 Giant Megathrust Earthquake
    www.nature.com/scientificreports OPEN Seismic rate variations prior to the 2010 Maule, Chile MW 8.8 giant megathrust earthquake Benoit Derode1*, Raúl Madariaga1,2 & Jaime Campos1 The MW 8.8 Maule earthquake is the largest well-recorded megathrust earthquake reported in South America. It is known to have had very few foreshocks due to its locking degree, and a strong aftershock activity. We analyze seismic activity in the area of the 27 February 2010, MW 8.8 Maule earthquake at diferent time scales from 2000 to 2019. We diferentiate the seismicity located inside the coseismic rupture zone of the main shock from that located in the areas surrounding the rupture zone. Using an original spatial and temporal method of seismic comparison, we fnd that after a period of seismic activity, the rupture zone at the plate interface experienced a long-term seismic quiescence before the main shock. Furthermore, a few days before the main shock, a set of seismic bursts of foreshocks located within the highest coseismic displacement area is observed. We show that after the main shock, the seismic rate decelerates during a period of 3 years, until reaching its initial interseismic value. We conclude that this megathrust earthquake is the consequence of various preparation stages increasing the locking degree at the plate interface and following an irregular pattern of seismic activity at large and short time scales. Giant subduction earthquakes are the result of a long-term stress localization due to the relative movement of two adjacent plates. Before a large earthquake, the interface between plates is locked and concentrates the exter- nal forces, until the rock strength becomes insufcient, initiating the sudden rupture along the plate interface.
    [Show full text]
  • The Growing Wealth of Aseismic Deformation Data: What's a Modeler to Model?
    The Growing Wealth of Aseismic Deformation Data: What's a Modeler to Model? Evelyn Roeloffs U.S. Geological Survey Earthquake Hazards Team Vancouver, WA Topics • Pre-earthquake deformation-rate changes – Some credible examples • High-resolution crustal deformation observations – borehole strain – fluid pressure data • Aseismic processes linking mainshocks to aftershocks • Observations possibly related to dynamic triggering The Scientific Method • The “hypothesis-testing” stage is a bottleneck for earthquake research because data are hard to obtain • Modeling has a role in the hypothesis-building stage http://www.indiana.edu/~geol116/ Modeling needs to lead data collection • Compared to acquiring high resolution deformation data in the near field of large earthquakes, modeling is fast and inexpensive • So modeling should perhaps get ahead of reproducing observations • Or modeling could look harder at observations that are significant but controversial, and could explore a wider range of hypotheses Earthquakes can happen without detectable pre- earthquake changes e.g.Parkfield M6 2004 Deformation-Rate Changes before the Mw 6.6 Chuetsu earthquake, 23 October 2004 Ogata, JGR 2007 • Intraplate thrust earthquake, depth 11 km • GPS-detected rate changes about 3 years earlier – Moment of pre-slip approximately Mw 6.0 (1 div=1 cm) – deviations mostly in direction of coseismic displacement – not all consistent with pre-slip on the rupture plane Great Subduction Earthquakes with Evidence for Pre-Earthquake Aseismic Deformation-Rate Changes • Chile 1960, Mw9.2 – 20-30 m of slow interplate slip over a rupture zone 920+/-100 km long, starting 20 minutes prior to mainshock [Kanamori & Cipar (1974); Kanamori & Anderson (1975); Cifuentes & Silver (1989) ] – 33-hour foreshock sequence north of the mainshock, propagating toward the mainshock hypocenter at 86 km day-1 (Cifuentes, 1989) • Alaska 1964, Mw9.2 • Cascadia 1700, M9 Microfossils => Sea level rise before1964 Alaska M9.2 • 0.12± 0.13 m sea level rise at 4 sites between 1952 and 1964 Hamilton et al.
    [Show full text]
  • Foreshock Sequences and Short-Term Earthquake Predictability on East Pacific Rise Transform Faults
    NATURE 3377—9/3/2005—VBICKNELL—137936 articles Foreshock sequences and short-term earthquake predictability on East Pacific Rise transform faults Jeffrey J. McGuire1, Margaret S. Boettcher2 & Thomas H. Jordan3 1Department of Geology and Geophysics, Woods Hole Oceanographic Institution, and 2MIT-Woods Hole Oceanographic Institution Joint Program, Woods Hole, Massachusetts 02543-1541, USA 3Department of Earth Sciences, University of Southern California, Los Angeles, California 90089-7042, USA ........................................................................................................................................................................................................................... East Pacific Rise transform faults are characterized by high slip rates (more than ten centimetres a year), predominately aseismic slip and maximum earthquake magnitudes of about 6.5. Using recordings from a hydroacoustic array deployed by the National Oceanic and Atmospheric Administration, we show here that East Pacific Rise transform faults also have a low number of aftershocks and high foreshock rates compared to continental strike-slip faults. The high ratio of foreshocks to aftershocks implies that such transform-fault seismicity cannot be explained by seismic triggering models in which there is no fundamental distinction between foreshocks, mainshocks and aftershocks. The foreshock sequences on East Pacific Rise transform faults can be used to predict (retrospectively) earthquakes of magnitude 5.4 or greater, in narrow spatial and temporal windows and with a high probability gain. The predictability of such transform earthquakes is consistent with a model in which slow slip transients trigger earthquakes, enrich their low-frequency radiation and accommodate much of the aseismic plate motion. On average, before large earthquakes occur, local seismicity rates support the inference of slow slip transients, but the subject remains show a significant increase1. In continental regions, where dense controversial23.
    [Show full text]
  • Earthquake Measurements
    EARTHQUAKE MEASUREMENTS The vibrations produced by earthquakes are detected, recorded, and measured by instruments call seismographs1. The zig-zag line made by a seismograph, called a "seismogram," reflects the changing intensity of the vibrations by responding to the motion of the ground surface beneath the instrument. From the data expressed in seismograms, scientists can determine the time, the epicenter, the focal depth, and the type of faulting of an earthquake and can estimate how much energy was released. Seismograph/Seismometer Earthquake recording instrument, seismograph has a base that sets firmly in the ground, and a heavy weight that hangs free2. When an earthquake causes the ground to shake, the base of the seismograph shakes too, but the hanging weight does not. Instead the spring or string that it is hanging from absorbs all the movement. The difference in position between the shaking part of the seismograph and the motionless part is Seismograph what is recorded. Measuring Size of Earthquakes The size of an earthquake depends on the size of the fault and the amount of slip on the fault, but that’s not something scientists can simply measure with a measuring tape since faults are many kilometers deep beneath the earth’s surface. They use the seismogram recordings made on the seismographs at the surface of the earth to determine how large the earthquake was. A short wiggly line that doesn’t wiggle very much means a small earthquake, and a long wiggly line that wiggles a lot means a large earthquake2. The length of the wiggle depends on the size of the fault, and the size of the wiggle depends on the amount of slip.
    [Show full text]
  • Energy and Magnitude: a Historical Perspective
    Pure Appl. Geophys. 176 (2019), 3815–3849 Ó 2018 Springer Nature Switzerland AG https://doi.org/10.1007/s00024-018-1994-7 Pure and Applied Geophysics Energy and Magnitude: A Historical Perspective 1 EMILE A. OKAL Abstract—We present a detailed historical review of early referred to as ‘‘Gutenberg [and Richter]’s energy– attempts to quantify seismic sources through a measure of the magnitude relation’’ features a slope of 1.5 which is energy radiated into seismic waves, in connection with the parallel development of the concept of magnitude. In particular, we explore not predicted a priori by simple physical arguments. the derivation of the widely quoted ‘‘Gutenberg–Richter energy– We will use Gutenberg and Richter’s (1956a) nota- magnitude relationship’’ tion, Q [their Eq. (16) p. 133], for the slope of log10 E versus magnitude [1.5 in (1)]. log10 E ¼ 1:5Ms þ 11:8 ð1Þ We are motivated by the fact that Eq. (1)istobe (E in ergs), and especially the origin of the value 1.5 for the slope. found nowhere in this exact form in any of the tra- By examining all of the relevant papers by Gutenberg and Richter, we note that estimates of this slope kept decreasing for more than ditional references in its support, which incidentally 20 years before Gutenberg’s sudden death, and that the value 1.5 were most probably copied from one referring pub- was obtained through the complex computation of an estimate of lication to the next. They consist of Gutenberg and the energy flux above the hypocenter, based on a number of assumptions and models lacking robustness in the context of Richter (1954)(Seismicity of the Earth), Gutenberg modern seismological theory.
    [Show full text]
  • What Is an Earthquake?
    A Violent Pulse: Earthquakes Chapter 8 part 2 Earthquakes and the Earth’s Interior What is an Earthquake? Seismicity • ‘Earth shaking caused by a rapid release of energy.’ • Seismicity (‘quake or shake) cause by… – Energy buildup due tectonic – Motion along a newly formed crustal fracture (or, stresses. fault). – Cause rocks to break. – Motion on an existing fault. – Energy moves outward as an expanding sphere of – A sudden change in mineral structure. waves. – Inflation of a – This waveform energy can magma chamber. be measured around the globe. – Volcanic eruption. • Earthquakes destroy – Giant landslides. buildings and kill people. – Meteorite impacts. – 3.5 million deaths in the last 2000 years. – Nuclear detonations. • Earthquakes are common. Faults and Earthquakes Earthquake Concepts • Focus (or Hypocenter) - The place within Earth where • Most earthquakes occur along faults. earthquake waves originate. – Faults are breaks or fractures in the crust… – Usually occurs on a fault surface. – Across which motion has occurred. – Earthquake waves expand outward from the • Over geologic time, faulting produces much change. hypocenter. • The amount of movement is termed displacement. • Epicenter – Land surface above the focus pocenter. • Displacement is also called… – Offset, or – Slip • Markers may reveal the amount of offset. Fence separated by fault 1 Faults and Fault Motion Fault Types • Faults are like planar breaks in blocks of crust. • Fault type based on relative block motion. • Most faults slope (although some are vertical). – Normal fault • On a sloping fault, crustal blocks are classified as: • Hanging wall moves down. – Footwall (block • Result from extension (stretching). below the fault). – Hanging wall – Reverse fault (block above • Hanging wall moves up.
    [Show full text]
  • Hypocenter and Focal Mechanism Determination of the August 23, 2011 Virginia Earthquake Aftershock Sequence: Collaborative Research with VA Tech and Boston College
    Final Technical Report Award Numbers G13AP00044, G13AP00043 Hypocenter and Focal Mechanism Determination of the August 23, 2011 Virginia Earthquake Aftershock Sequence: Collaborative Research with VA Tech and Boston College Martin Chapman, John Ebel, Qimin Wu and Stephen Hilfiker Department of Geosciences Virginia Polytechnic Institute and State University 4044 Derring Hall Blacksburg, Virginia, 24061 (MC, QW) Department of Earth and Environmental Sciences Boston College Devlin Hall 213 140 Commonwealth Avenue Chestnut Hill, Massachusetts 02467 (JE, SH) Phone (Chapman): (540) 231-5036 Fax (Chapman): (540) 231-3386 Phone (Ebel): (617) 552-8300 Fax (Ebel): (617) 552-8388 Email: [email protected] (Chapman), [email protected] (Ebel), [email protected] (Wu), [email protected] (Hilfiker) Project Period: July 2013 - December, 2014 1 Abstract The aftershocks of the Mw 5.7, August 23, 2011 Mineral, Virginia, earthquake were recorded by 36 temporary stations installed by several institutions. We located 3,960 aftershocks from August 25, 2011 through December 31, 2011. A subset of 1,666 aftershocks resolves details of the hypocenter distribution. We determined 393 focal mechanism solutions. Aftershocks near the mainshock define a previously recognized tabular cluster with orientation similar to a mainshock nodal plane; other aftershocks occurred 10-20 kilometers to the northeast. Detailed relocation of events in the main tabular cluster, and hundreds of focal mechanisms, indicate that it is not a single extensive fault, but instead is comprised of at least three and probably many more faults with variable orientation. A large percentage of the aftershocks occurred in regions of positive Coulomb static stress change and approximately 80% of the focal mechanism nodal planes were brought closer to failure.
    [Show full text]
  • Earthquake Location Accuracy
    Theme IV - Understanding Seismicity Catalogs and their Problems Earthquake Location Accuracy 1 2 Stephan Husen • Jeanne L. Hardebeck 1. Swiss Seismological Service, ETH Zurich 2. United States Geological Survey How to cite this article: Husen, S., and J.L. Hardebeck (2010), Earthquake location accuracy, Community Online Resource for Statistical Seismicity Analysis, doi:10.5078/corssa-55815573. Available at http://www.corssa.org. Document Information: Issue date: 1 September 2010 Version: 1.0 2 www.corssa.org Contents 1 Motivation .................................................................................................................................................. 3 2 Location Techniques ................................................................................................................................... 4 3 Uncertainty and Artifacts .......................................................................................................................... 9 4 Choosing a Catalog, and What to Expect ................................................................................................ 25 5 Summary, Further Reading, Next Steps ................................................................................................... 30 Earthquake Location Accuracy 3 Abstract Earthquake location catalogs are not an exact representation of the true earthquake locations. They contain random error, for example from errors in the arrival time picks, as well as systematic biases. The most important source of systematic
    [Show full text]
  • Qt88c3k67v.Pdf
    UC Berkeley UC Berkeley Previously Published Works Title Early aftershocks and afterslip surrounding the 2015 Mw 8.4 Illapel rupture Permalink https://escholarship.org/uc/item/88c3k67v Authors Huang, H Xu, W Meng, L et al. Publication Date 2017 DOI 10.1016/j.epsl.2016.09.055 Peer reviewed eScholarship.org Powered by the California Digital Library University of California Early aftershocks and afterslip surrounding the 2015 Mw 8.4 Illapel rupture Author links open overlay panel HuiHuang a WenbinXu bc LingsenMeng a RolandBürgmann b Juan CarlosBaez d Show more https://doi.org/10.1016/j.epsl.2016.09.055 Get rights and content Highlights • Missing early aftershocks and repeaters are recovered by the matched- filtermethod. • Differential southward and northward expansion of early aftershocks are observed. • Repeaters and geodetic data reveal afterslip around the Illapel mainshock rupture. Abstract On 16 September 2015, the Mw 8.4 Illapel earthquake ruptured a section of the subduction thrust on the west coast of central Chile. The mainshock was followed by numerous aftershocks including some normal-faulting events near the trench. We apply a template matching approach to improve the completeness of early aftershocks within one month of the mainshock. To constrain the distribution of afterslip, we utilize repeating earthquakes among the aftershocks and perform a joint slip inversion of postseismic GPS and InSAR data. The results show that the aftershock zone abruptly expands to the south ∼14 h after the mainshock while growing relatively continuously to the north within the first day. The repeating earthquakes accompanying the early expansion suggest that aseismic afterslip on the subduction thrust surrounding the coseismic rupture is an important triggering mechanism of aftershocks in addition to stress transfer or poroelastic effects.
    [Show full text]
  • Evaluation of Automatic Hypocenter Determination in the JMA Unified
    Tamaribuchi Earth, Planets and Space (2018) 70:141 https://doi.org/10.1186/s40623-018-0915-4 FULL PAPER Open Access Evaluation of automatic hypocenter determination in the JMA unifed catalog Koji Tamaribuchi* Abstract The Japan Meteorological Agency (JMA) unifed seismic catalog has been widely used for research and disaster pre- vention purposes for more than 20 years. Since the introduction in April 2016 of an improved method of automatic hypocenter determinations (PF method), the number of detected earthquakes has almost doubled due to a decrease in the completeness magnitude around the Tohoku region, where seismicity has been very active in the aftermath of the 2011 Tohoku earthquake. Automatically processed hypocenters of small events, accepted without manual modifcation, now make up approximately 70% of new events in the JMA unifed catalog. In this paper, we show that the introduction of automated processing did not systematically bias the quality of the JMA unifed catalog. Approxi- mately 90% of automatically processed hypocenters were less than 1 km from their manually reviewed locations in inland and shallow areas. We also considered the use of automated event characterization in real-time monitoring of earthquake sequences using the example of the April 2016 Kumamoto earthquake sequence, when the PF method could have supplied the catalog with about 70,000 events in real time over the course of 2 months. We show that the PF method is capable of monitoring the migration or expansion of the hypocentral distribution and can support sta- tistical analyses such as variations of the b-value distribution. Further improvements in automatic hypocenter deter- mination will contribute to a better understanding of seismicity as well as rapid risk assessment, especially in cases of swarms and aftershocks.
    [Show full text]