Postseismic Displacements Following the 2007 Noto Peninsula Earthquake Detected by Dense GPS Observation

Total Page:16

File Type:pdf, Size:1020Kb

Postseismic Displacements Following the 2007 Noto Peninsula Earthquake Detected by Dense GPS Observation LETTER Earth Planets Space, 60, 139–144, 2008 Postseismic displacements following the 2007 Noto peninsula earthquake detected by dense GPS observation Manabu Hashimoto1, Hiroaki Takahashi2, Ryosuke Doke3, Minoru Kasahara2, Akira Takeuchi3, Kenusuke Onoue1, Yoshinobu Hoso1, Yo Fukushima1, Kajuro Nakamura1, Fumio Ohya1, Ryo Honda2, Masayoshi Ichiyanagi2, Teruhiro Yamaguchi2, Takahiro Maeda2, and Yoshihiro Hiramatsu4 1Disaster Prevention Research Institute, Kyoto University, Gokasho, Uji, Kyoto 611-0011, Japan 2Institute of Seismology and Volcanology, Graduate School, Hokkaido University, Sapporo 060-0810, Japan 3Graduate School of Science and Engineering, University of Toyama, Gofuku 3190, Toyama, Toyama 930-8555, Japan 4Graduate School of Natural Science and Technology, Kanazawa University, Kakuma, Kanazawa, Ishikawa 920-1192, Japan (Received June 29, 2007; Revised November 3, 2007; Accepted November 13, 2007; Online published February 19, 2008) We have been conducting dense GPS observation in and around the epicentral region of the 2007 Noto peninsula earthquake since March 25, 2007, in order to detect postseismic displacements. Continuous observation has been underway at 12 sites to fill the gap of GEONET. Preliminary analysis of data up to early May shows that initial postseismic displacement rapidly decayed within 20 days after the occurrence of the mainshock. Horizontal displacements do not exceed 20 mm even at sites above the aftershock zone for this period. We also found a maximum uplift of about 20 mm there. Inversion of postseismic displacements with the variable slip model suggests a nearly right-lateral afterslip of less than 5 cm on the shallow portion of the source fault. Fitting a theoretical function to a time series of coordinate changes also suggests that the observed postseismic displacements might have been generated by afterslip. Key words: Noto peninsula earthquake, GPS, crustal deformation, postseismic deformation, afterslip. 1. Introduction a much denser network of GPS observation is required to On March 25, 2007, a Mjma = 6.9 earthquake hit the reveal the detailed distribution of postseismic displacement northwestern part of the Noto peninsula, central Japan, following the mainshock. With this aim, we started carry- claiming one casualty and causing serious structural dam- ing out dense continuous GPS observations in and around age, especially in the town of Monzen, Wajima city. This is the epicentral region following the occurrence of the Noto the fourth shallow inland earthquake of M 6.5 or larger in peninsula earthquake. the Japanese islands during the present decade. Postseismic displacements are often observed following large earthquakes, and these convey important information 2. GPS Observation and Data Processing on frictional properties of the source fault (e.g. Marone et Aftershocks were distributed close to the northwestern al., 1991; Perfettini et al., 2004; Montesi, 2004; Miyazaki coast of Noto peninsula, where only three GEONET sites et al., 2004), viscoelastic structure in the crust and upper are operated. Therefore, a much denser network is essen- mantle (e.g. Thatcher et al., 1980; Pollitz, 2005), pore pres- tial to reveal a detailed distribution of postseismic displace- sure changes (e.g. Peltzer et al., 1998; Jonsson´ et al., 2003), ment. Responding to the occurrence of the Noto peninsula among others. These factors may also control the strain ac- earthquake, three groups consisting of researchers from sev- cumulation and release process associated with earthquake eral universities started continuous GPS observation in and generation. Thus, it is essential to carry out such observa- around the epicentral area. During the period from 25 to tions as soon as possible to reveal the detailed history of the 29 March, 2007, we established 12 observation sites which evolution of postseismic deformations. are equipped with dual-frequency receivers (Fig. 1). We se- In Japan, the Geographical Survey Institute’s Earth Ob- lected apparently stable sites, such as reinforced concrete servation NETwork (hereafter GEONET) is under opera- buildings, basements of monuments etc, and fixed antenna tion and has revealed coseismic and postseismic displace- on the roofs, because of limited time available to select sites ments following large earthquakes (e.g. Sagiya et al., 2000; (Fig. 2). Phase data were sampled at intervals of 30 s, but Sagiya, 2004). However, average spacing is about 20– sampling rate was higher (0.1–1 Hz) at some sites in order 25 km for GEONET, which is as large as the size of af- to conduct a kinematic analysis of data. For the static anal- tershock distribution of the earthquakes of M 7. Therefore, ysis, however, all data were decimated to 30-s sampled data with the utility TEQC (Estey and Meertens, 1999). The el- ◦ Copyright c The Society of Geomagnetism and Earth, Planetary and Space Sci- evation mask was set at 5–10 . ences (SGEPSS); The Seismological Society of Japan; The Volcanological Society For the analysis of GPS data, we used the Bernese of Japan; The Geodetic Society of Japan; The Japanese Society for Planetary Sci- ences; TERRAPUB. software (ver. 5.0) for static positioning of daily co- ordinates (Hugentobler et al., 2005) with IGS final 139 140 M. HASHIMOTO et al.: POSTSEISMIC DISPLACEMENTS DETECTED BY GPS CHHR UIRI UNOY Fig. 2. Photos of representative sites: Unoya (UNOY), Chihara (CHHR) and Uiri (UIRI) sites. Day from Mar.25 -5 0 5 10 15 20 25 30 35 40 45 50 20 20 020969:Mar.20 - May.5,2007 10 10 0 0 Fig. 1. Distribution of sites. Solid diamonds, solid and open squares indi- -10 -10 cate GPS sites deployed by Hokkaido University and Kanazawa Univer- Changes from Mar.21 in mm sity, University of Toyama, and Disaster Prevention Research Institute, Changes from Mar.21 in mm Kyoto University, respectively. Open diamonds show GEONET sites. -20 -20 E-W N-S U-D Large triangle indicates GEONET 020969, which is the reference point Mw6.9 for displacements. Gray circles are epicenters of aftershocks determined -5 0 5 10 15 20 25 30 35 40 45 50 by JMA and MRI (2007). Dark and light gray stars show the epicenter Day from Mar.25 of the mainshock and epicenters of aftershocks larger than M 4.5, re- spectively. Gray lines indicate surface traces of active faults (Research Fig. 3. Time series of changes of coordinates of GEONET site 020969 in Group of Active Faults of Japan, 1991). mm since March 21, 2007, in ITRF2005. Circles, diamonds and trian- gles show eastward, northward and vertical components, respectively. Open and solid symbols are original and spatial-filtered data, respec- tively. Error bars indicate 3-σ. ephemeredes, earth rotation parameters and satellite clock coefficients. All data from the IGS, GEONET and temporal sites were processed simultaneously to Amin (2002). We used the sites far away from the epicenter obtain their coordinates in ITRF2005 (http://itrf.ensg. of the mainshock for the calculation of common mode bias. ign.fr/ITRF solutions/2005/ITRF2005.php). We chose re- The common mode bias that was obtained was subtracted mote IGS sites, such as Tsukuba (TSKB), Daejong (DAEJ), from the original time series. Figure 3 shows the observed Yuzhno-Sakhalinsk (YSSK) and others, as fiducial sites and spatially-filtered time series of GEONET 020969 in and tried to strongly constrain to their predicted coordi- ITRF2005. Horizontal components are rather stable, while nates in ITRF2005 in order to align the site coordinates in vertical ones have as large a scatter as 10 mm. Judging ITRF2005. However, their preliminarily estimated coordi- from the stability of the horizontal components and epicen- nates sometimes largely differed from the ITRF2005 pre- tral distance, the choice of GEONET 020969 as a reference dictions. In such a case, Bernese automatically excludes site is adequate. Due to the large scatter in the vertical com- such sites from the list of fiducial sites and calculates coor- ponent, however, it is hard to discuss deformation smaller dinates. than 10 mm. Time series of coordinates usually have relatively large scatters. We suspect that the scatter may be common noise 3. Observed Postseismic Displacement that can be attributed to errors in ephemeredes and satellite Figures 4–6 show the time series of coordinate changes clock delays and other common effects, such as troposh- relative to 020969. Most GEONET sites do not show any phere/ionosphere disturbances. In order to reduce such dis- significant movements since they are located relatively far tubances, we fix GEONET 020969. Geographical Survey from the epicenter. However, several sites, such as 940053, Institute (2007) (hereafter GSI) selected 020969 as the ref- 020972, AKGM, CHHR, UNOY, WMII and SNDI, show erence point for the calculation of coseismic displacements. a different trend. SNDI in particular rapidly moved west- However, we must check the stability of this site. In the cal- ward during the first 3 days and seemed to decelerate after- culation of time series we often apply the spatial filter tech- ward. It is unfortunate that several data points are missing nique proposed by Wdowinski et al. (1997) and Tabei and due to power failure and we could not know its detailed M. HASHIMOTO et al.: POSTSEISMIC DISPLACEMENTS DETECTED BY GPS 141 Day from Mar.25, 2007 Day from Mar.25, 2007 -4 -2 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 -4 -2 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 55 55 110 110 50 50 100 100 940053 45 45 90 90 950248 940053 80 80 40 40 950253 950248 70 70
Recommended publications
  • Recipe for Predicting Strong Ground Motion from Crustal Earthquake Scenarios
    Pure Appl. Geophys. 168 (2011), 85–104 Ó 2010 The Author(s) This article is published with open access at Springerlink.com DOI 10.1007/s00024-010-0150-9 Pure and Applied Geophysics Recipe for Predicting Strong Ground Motion from Crustal Earthquake Scenarios 1 2 KOJIRO IRIKURA and HIROE MIYAKE Abstract—We developed a recipe for predicting strong ground processes of earthquakes, propagation-path and site- motions based on a characterization of the source model for future amplification effects, and the relation between ground crustal earthquakes. From recent developments of waveform inver- sion of strong motion data used to estimate the rupture process, we motion and damage. Predicting strong ground motion have inferred that strong ground motion is primarily related to the slip is one of the critical factors for mitigating disasters heterogeneity inside the source rather than average slip in the entire caused by future earthquakes. rupture area. Asperities are characterized as regions that have large slip relative to the average slip on the rupture area. The asperity areas, So far, most strong motion predictions in earthquake as well as the total rupture area, scale with seismic moment. We hazard analyses have been made using empirical determined that the areas of strong motion generation approximately attenuation-distance curves for peak ground accelera- coincide with the asperity areas. Based on the scaling relationships, the deductive source model for the prediction of strong ground tion (PGA), peak ground velocity (PGV), and response motions is characterized by three kinds of parameters: outer, inner, spectra, in which the source information is defined and extra fault parameters.
    [Show full text]
  • International Workshop on Earthquake Preparation Process 2017
    International Workshop on Earthquake Preparation Process 2017 - Observation, Validation, Modeling, Forecasting - Program & Abstracts IWEP2017 May 26-27, 2017 Chiba University, Chiba, Japan International Workshop on Earthquake Preparation Process 2017 - Observation, Validation, Modeling, Forecasting - (IWEP4) Sponsor Chiba University, Japan Co-sponsor The Earthquake Prediction Society of Japan integrated Study and Test for Earthquake Precursors (iSTEP-4), Taiwan Supported by National Central University, Taiwan Chapman University, US Tokai University, Japan Date: May 26-27, 2017 Venue:Conference Hall, the 1st floor of Sciences and Technology Building No.1, Nishi-Chiba Campus, Chiba University, Chiba, Japan Contacts: Dr. Katsumi Hattori (Chiba University, Japan) +81-43-290-280 Contents Program……………………………………………….…...……1-8 Oral presentations…………………………………..….…….1-4 Poster presentations…………………………………..………5 Abstracts…………………………..……………….…….….9-57 Oral presentations…………………….………………...……9-40 Poster presentations………………….…………..………….41-57 Program May 26, 2017 Morning 0820-0850 Registration 0850-0900 Opening: Dr. Katsumi Hattori Chair: Dr. Katsumi Hattori 0900-0930 Dimitar Ouzounov et al., Radon activity - the hidden driver behind the atmospheric pre-earthquake anomalies 0930-1000 Ching-Chou Fu et al., Exploring the relationship between soil degassing and seismic activity by continuous radon monitoring in the Longitudinal Valley of eastern Taiwan (Invited) 1000-1030 Ramesh P. Singh et al., Changes in Ground water level and Linkage with the Meteorological/Atmospheric
    [Show full text]
  • Ings Damaged During the Recent Large Earthquakes in Japan
    16th World Conference on Earthquake Engineering, 16WCEE 2017 Santiago Chile, January 9th to 13th 2017 Paper N° 4747 (Abstract ID) Registration Code: S-L1472725668 STUDY ON DAMAGE AND SEISMIC PERFORMANCE OF TIMBER BUILD- INGS DAMAGED DURING THE RECENT LARGE EARTHQUAKES IN JAPAN T. Tsuchimoto(1), N. Kawai(2), T. Nakagawa(3) (1) Chief Research Engineer, Building Research Institute, [email protected] (2) Professor, Kogakuin University, [email protected] (3) Senior Research Officer, National Institute for Land and Infrastructure Management, MLIT, [email protected] Abstract It is said that Japan has come to the term of seismic activity. Once, the damage of 1948 Fukui Earthquake triggered the establishment of Building Standard Law of Japan including the seismic regulation. The regulation for timber building adopted the way to set up enough shear wall by using the brace. After that, the regulations were reformed and revised several times at every damage due to earthquake. The required shear wall length increased at every revision of regulation. Timber buildings in Japan were damaged due to large earthquake since 2000, for example, 2000 West of Tottori Pref., 2004 Chuetsu, 2005 Fukuoka west off the coast, 2007 Noto Peninsula, 2007 offshore of Chuetsu, 2011 off the Pacific coast of Tohoku earthquake and 2014 Nagano Kamishiro Fault. The damage to timber buildings due to every earthquake was surveyed on site. Several wood houses were picked up and investigated in detail, for example, floor plans, shear wall length, damage states and structural specification. The reason and the characteristic of damage were discussed. The results of these surveys and studies were summarized as follows; - The 14 major earthquake suffering heavy damage to wood houses has occurred after the year of 2000.
    [Show full text]
  • Strong Ground Motion Simulation for the 2005 Fukuoka Earthquake (Mw6.6) by Stochastic Finite-Fault Simulations
    th The 14 World Conference on Earthquake Engineering October 12-17, 2008, Beijing, China Strong Ground Motion Simulation for the 2005 Fukuoka Earthquake (Mw6.6) by Stochastic Finite-Fault Simulations 1 1 2 MITSUTAKA OSHIMA , HIROSHI TAKENAKA , and HIROSHI KAWASE 1 Dept. of Earth and Planetary Sciences, Kyushu University,6-10-1 Hakozaki, Fukuoka, Japan 2Disaster Prevention Research Institute, Kyoto University, Gokasho, Uji, Kyoto, Japan Email: [email protected], [email protected], [email protected] ABSTRACT : The theoretical basis of the stochastic finite-fault modeling was proposed by Boore (1983) and this technique has been developed by many researchers (g.e., Beresnev and Atkinson, 1997) and many papers were published applying it to real simulations. Recently, Motazedian and Atkinson (2005) have introduced the concept of “dynamic corner frequency" and “pulsing area" into the stochastic finite-fault modeling. It is possible to simulate strong ground motions from earthquakes in a more realistic way with these new concepts. We made a further modification that enables us to treat stress parameter variable on the fault plane, while the stress parameter is constant on the whole fault plane in the conventional stochastic finite-fault modeling. With this modification, we can introduce stress-variable source models such as asperity model. We apply our modified stochastic finite-fault modeling to the 2005 West Off Fukuoka earthquake (Mw6.6) that occurred in northwest Kyushu, Japan, and check whether we could satisfactorily reproduce the observed strong ground motion data. We adopt an asperity model and assign higher stress parameter to the asperity than that of the background region.
    [Show full text]
  • Aftershock Statistics of the 1999 Chi–Chi, Taiwan Earthquake and the Concept of Omori Times
    Pure Appl. Geophys. 170 (2013), 221–228 Ó 2012 Springer Basel AG DOI 10.1007/s00024-011-0445-5 Pure and Applied Geophysics Aftershock Statistics of the 1999 Chi–Chi, Taiwan Earthquake and the Concept of Omori Times 1,2 1 1,3 2 YA-TING LEE, DONALD L. TURCOTTE, JOHN B. RUNDLE, and CHIEN-CHIH CHEN Abstract—In this paper we consider the statistics of the after- the magnitude of the main shock and the largest shock sequence of the m = 7.65 20 September 1999 Chi–Chi, aftershock. This law states that the difference in Taiwan earthquake. We first consider the frequency-magnitude statistics. We find good agreement with Gutenberg–Richter scaling magnitude, Dm, between the main shock and its but find that the aftershock level is anomalously high. This level is largest aftershock is approximately constant inde- quantified using the difference in magnitude between the main pendent of the magnitude of the main shock shock and the largest inferred aftershock DmÃ: Typically, Dmà is in the range 0.8–1.5, but for the Chi–Chi earthquake the value is Dm ¼ mms À mmax as ð2Þ Dmà = 0.03. We suggest that this may be due to an aseismic slow- earthquake component of rupture. We next consider the decay rate where mms is the magnitude of the main shock and of aftershock activity following the earthquake. The rates are well m is the magnitude of the largest aftershock, approximated by the modified Omori’s law. We show that the max as distribution of interoccurrence times between aftershocks follow a typically Dm & 1.2 (BA˚ TH 1965).
    [Show full text]
  • Try Making a “My Map” Specialists Mrs. Power Fill in Map (My Map
    Try Making a “My Map” Lively dialog and chit chat is welcome! Share your wisdom and knowledge energetically! Information on the Information provided by earthquake local community survivors Gathering Specialists information Mrs. Power Predicting the damage and identifying the causes of the damage (Fieldwork) Survey the field Rediscover your local community Fill in map (My Map) Simulation Applied Science Technical Group, Tohoku Regional Office, The (How to use it) Institution of Professional Engineers, Japan Earthquake Preparedness Working Group Disaster Prevention Committee, Miyagi Professional Engineers Association 1 Brochure created by Applied Science Technical Group It could happen anywhere and at any time. But… Since the beginning of the century, starting with the Great Hanshin earthquake of 1995, earthquakes have caused damage in various regions across Japan one after another- 2000 Western Tottori earthquake, 2001 Geiyo earthquake, 2003 Northern Miyagi earthquake, 2003 Hokkaido earthquake, 2004 Chuetsu earthquake, 2005 Fukuoka earthquake, 2007 Chuetsu-oki Sanriku-oki earthquake, and 2008 Iwate-Miyagi Nairiku earthquake. earthquake It is predicted that, in the next 30 years, a Sanriku-oki earthquake and a Miyagi-oki earthquake will undoubtedly occur (along the Pacific Ocean side of northeastern Japan) and a Tokai earthquake, Tonankai earthquake, and Nankai earthquake (along the Pacific Ocean side of western Japan). There may also be earthquakes in unexpected areas, as was the case with the 2008 Iwate-Miyagi Nairiku earthquake. Miyagi-oki earthquake Even if we are able to predict these earthquakes, and even if we know their mechanism, we cannot prevent the earthquake from occurring. Tokai earthquake Active Tonankai earthquake fault zone Nankai earthquake Therefore, we need to prepare ourselves for calamities caused by these earthquakes.
    [Show full text]
  • Crustal Deformation Model of the Beppu-Shimabara Graben Area, Central Kyushu, Japan, Based on Inversion of Three-Component GNSS
    Mochizuki and Mitsui Earth, Planets and Space (2016) 68:177 DOI 10.1186/s40623-016-0550-x LETTER Open Access Crustal deformation model of the Beppu Shimabara graben area, central Kyushu, Japan, −based on inversion of three‑component GNSS data in 2000–2010 Kazuma Mochizuki and Yuta Mitsui* Abstract The 2016 Kumamoto earthquakes, including an Mw-7 right-lateral earthquake on April 15 (UTC), occurred along faults within the Beppu Shimabara graben in central Kyushu, Japan. Previous studies showed that the graben area was under heterogeneous− stress conditions with north–south T-axes and spreading in a north–south direction. Here, we construct a detailed crustal deformation model using three-component Global Navigation Satellite System data in 2000–2010 and considering the distribution of geological fault traces in this area. Our inversion analysis suggests that the strain accumulation rate for the right-lateral seismic slip segment (corresponding to the Futagawa fault), where the largest of the 2016 Kumamoto earthquakes ruptured, was several times smaller than the other segments in the Beppu Shimabara graben. Furthermore, we observe distinct subsidence along the Beppu Shimabara graben. Our base model− attributes the subsidence to deflation of magma reservoirs beneath volcanoes,− but the observed vertical velocities are poorly fit. In order to improve the fitting results for the vertical deformation, we need more sophisticated volcano-deformation model (such as a sill-like deformation source for Mt. Aso) or graben model. Keywords: The 2016 Kumamoto earthquakes, Kyushu, GNSS, Crustal deformation, Beppu Shimabara graben, Volcano, Block-fault model, Nucleation − Introduction Tada (1984, 1985) investigated crustal deformation The 2016 Kumamoto earthquakes, including an Mw-7 of the Kyushu region using repeated triangulation sur- earthquake on April 15 (UTC), struck a central part of veys for 1892–1982 and found that the Beppu−Shima- the island of Kyushu, southwest Japan.
    [Show full text]
  • Rupture Process of the 2005 West Off Fukuoka Prefecture, Japan, Earthquake
    LETTER Earth Planets Space, 58, 87–92, 2006 Rupture process of the 2005 West Off Fukuoka Prefecture, Japan, earthquake Haruo Horikawa Active Fault Research Center, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8567 Japan (Received August 10, 2005; Revised December 3, 2005; Accepted December 7, 2005; Online published January 27, 2006) The 2005 West off Fukuoka Prefecture, Japan, earthquake (MS 6.7) is a moderate-size crustal event. The rupture process of this earthquake was inferred from strong motion data. It is observed regardless of epicentral distance and azimuth that the initial P-wave with a small amplitude continued for 2.5–4 s, suggesting prominent initial rupture. The duration of the initial rupture was estimated to be 3.3 s from travel time analysis. Inversion analysis shows that the initial rupture continued for 3.5 s and the overall rupture finished within 10 s. An area of large slip with two peaks of slip amount appeared to the southeast of the hypocenter, and the location of the large-slip area is consistent with the origin of the main rupture inferred from travel time analysis. Small amount of slip was found at and around the hypocenter, and the rupture velocity was slow there. The seismic moment 18 of this earthquake was estimated to be 5.7 × 10 Nm (Mw 6.4). The maximum stress drop calculated from the derived slip distribution exceeded 10 MPa. An area of negative stress drop expanded, corresponding to the distribution of small slip. Forward modeling of the observed waveforms suggests that the negative stress drop is not fully attributable to poor resolution of inversion analysis for small amount of slip.
    [Show full text]
  • List of Scientific Papers and Meeting Presentations
    List of scientific papers and meeting presentations • Peer-reviewed papers (ISI Journals) 1. Shimojo, K., Enescu, B., Yagi, Y. et al. Nucleation process of the 2011 northern Nagano earthquake from nearby seismic observations, Scientific Reports, 11, 8143, https://doi.org/10.1038/s41598-021-86837-4, 2021. 2. Borleanu, F., Petrescu, L., Enescu, B., Popa, M., and M. Radulian, The missing craton edge: Crustal structure of the East European Craton beneath the Carpathian Orogen revealed by double-difference tomography, Global and Planetary Change, 197, 103390, doi:10.1016/j.gloplacha.2020.103390, 2021. 3. Li, Y., Wang, D., Xu, S., Fang, L., Cheng, Y., Luo, G., Yan, B., Enescu, B., and J. Mori, Thrust and Conjugate Strike-Slip Faults in the 17 June 2018 MJMA 6.1 (⁠Mw 5.5) Osaka, Japan, Earthquake Sequence, Seismol. Res. Lett., 90(6), 2132–2141. doi:10.1785/0220190122, 2019. 4. Omi, T., Ogata, Y., Shiomi, K., Enescu, B., Sawazaki, K., and K. Aihara, Implementation of a Real-Time System for Automatic Aftershock Forecasting in Japan, Seismol. Res. Lett., 90 (1): 242–250. doi:10.1785/0220180213, 2018. 5. Gulia, L., Rinaldi, A. P., Tormann, T., Vannucci, G., Enescu, B., and S. Wiemer, The effect of a mainshock on the size distribution of the aftershocks, Geophys. Res. Lett., 45, 13,277– 13,287, doi:10.1029/2018GL080619, 2018 (AGU Highlighted paper – Cook, T., How do main shocks affect subsequent earthquakes?, Eos, 100, doi:10.1029/2019EO119129, 2019.) 6. Tamaribuchi, K., Y. Yagi, B. Enescu, and S. Hirano, Characteristics of foreshock activity inferred from the JMA earthquake catalog, Earth, Planets and Space, 70:90, doi:10.1186/s40623-018-0866-92018, 2018.
    [Show full text]