Missing early aftershocks following relatively large or moderate earthquakes can cause significant bias in the analysis of seismic catalogs.In this paper,we systematically address the aftershock missing problem for fi...Missing early aftershocks following relatively large or moderate earthquakes can cause significant bias in the analysis of seismic catalogs.In this paper,we systematically address the aftershock missing problem for five earthquake sequences associated with moderate-size events that occurred inland Japan,by using a stochastic replenishing method.The method is based on the notion that if a point process(e.g.,earthquake sequence)with timeindependent marks(e.g.,magnitudes)is completely observed,it can be transformed into a homogeneous Poisson process by a bi-scale empirical transformation.We use the Japan Meteorological Agency(JMA)earthquake catalog to select the aftershock data and replenish the missing early events using the later complete part of each aftershock sequence.The time windows for each sequence span from 6 months before the mainshock to three months after.The semi-automatic spatial selection uses a clustering method for the epicentral selection of earthquakes.The results obtained for the original JMA catalog and replenished datasets are compared to get insight into the biases that the missing early aftershocks may cause on the Omori-Utsu law parameters’estimation,characterizing the aftershock decay with time from the mainshock.We have also compared the Omori-Utsu law parameter estimates for two datasets following the same mainshock;the first dataset is the replenished sequence,while the second dataset has been obtained by waveform-based analysis to detect early aftershocks that are not recorded in the JMA catalog.Our results demonstrate that the Omori-Utsu law parameters estimated for the replenished datasets are robust with respect to the threshold magnitude used for the analyzed datasets.Even when using aftershock time windows as short as three days,the replenished datasets provide stable Omori-Utsu law parameter estimations.The p-values for all the analyzed sequences are about 1.1 and c-values are significantly smaller compared to those of original datasets.Our findings prove that the replenishment method is a fast,reliable approach to address the missing aftershock problem.展开更多
On April 25, 2015, Nepal was struck by the MW7.8 Gorkha earthquake followed by an intense aftershock sequence. It was one of the most destructive earthquakes in the Himalayan arc, causing more than 8900 fatalities. In...On April 25, 2015, Nepal was struck by the MW7.8 Gorkha earthquake followed by an intense aftershock sequence. It was one of the most destructive earthquakes in the Himalayan arc, causing more than 8900 fatalities. In this study, we analyzed the dataset (429 events, magnitude of completeness (Mc) ≥ 4.2 local magnitude) of the first 45 days after the Gorkha earthquake to estimate the seismicity parameters b-value, D-value, and p-value. We used the maximum likelihood method to estimate the b-value and Omori-Utsu parameters, whereas the correlation integral method was applied to estimate the fractal dimension (D-value). The analysis was carried out using running and sliding window techniques. The lowest b-value (0.57 ± 0.04) and the highest D-value (1.65 ± 0.02) were computed at the time of the Gorkha earthquake, after which the b-value significantly increased to a maximum of 1.57. It again dropped to 0.93 at the time of the major aftershock on May 12, 2015. The D-value showed an initial quick drop and then decreased in a wavy pattern until the end of the study period, indicating the clustering and scattering of earthquakes in a fault region. The b-value contour map identified the eastern part of the study area as a high stress region (b = ~0.8), implying that the stress shifted to that region. The D-value contour map reveals that the seismogenic structure shifted from linear to planar in the region. The rate of aftershock decay (p = 0.86 ± 0.04) for a short period reflects that the level of stress decreased rapidly. This study helps to understand the level of stress and seismicity pattern of a region, which could be useful for aftershock studies.展开更多
基金Bogdan Enescu is grateful to the Executive Agency for Higher Education,Research,Development and Innovation Funding(UEFISCDI),Romania,through the project PNIII-P4-ID-PCE-2020-1361,119 PCE/2021(AFROS)for support.Jiancang Zhuang was supported by MEXT Project for Seismology toward Research Innovation with Data of Earthquake(STAR-E)Grant Number JPJ010217.
文摘Missing early aftershocks following relatively large or moderate earthquakes can cause significant bias in the analysis of seismic catalogs.In this paper,we systematically address the aftershock missing problem for five earthquake sequences associated with moderate-size events that occurred inland Japan,by using a stochastic replenishing method.The method is based on the notion that if a point process(e.g.,earthquake sequence)with timeindependent marks(e.g.,magnitudes)is completely observed,it can be transformed into a homogeneous Poisson process by a bi-scale empirical transformation.We use the Japan Meteorological Agency(JMA)earthquake catalog to select the aftershock data and replenish the missing early events using the later complete part of each aftershock sequence.The time windows for each sequence span from 6 months before the mainshock to three months after.The semi-automatic spatial selection uses a clustering method for the epicentral selection of earthquakes.The results obtained for the original JMA catalog and replenished datasets are compared to get insight into the biases that the missing early aftershocks may cause on the Omori-Utsu law parameters’estimation,characterizing the aftershock decay with time from the mainshock.We have also compared the Omori-Utsu law parameter estimates for two datasets following the same mainshock;the first dataset is the replenished sequence,while the second dataset has been obtained by waveform-based analysis to detect early aftershocks that are not recorded in the JMA catalog.Our results demonstrate that the Omori-Utsu law parameters estimated for the replenished datasets are robust with respect to the threshold magnitude used for the analyzed datasets.Even when using aftershock time windows as short as three days,the replenished datasets provide stable Omori-Utsu law parameter estimations.The p-values for all the analyzed sequences are about 1.1 and c-values are significantly smaller compared to those of original datasets.Our findings prove that the replenishment method is a fast,reliable approach to address the missing aftershock problem.
文摘On April 25, 2015, Nepal was struck by the MW7.8 Gorkha earthquake followed by an intense aftershock sequence. It was one of the most destructive earthquakes in the Himalayan arc, causing more than 8900 fatalities. In this study, we analyzed the dataset (429 events, magnitude of completeness (Mc) ≥ 4.2 local magnitude) of the first 45 days after the Gorkha earthquake to estimate the seismicity parameters b-value, D-value, and p-value. We used the maximum likelihood method to estimate the b-value and Omori-Utsu parameters, whereas the correlation integral method was applied to estimate the fractal dimension (D-value). The analysis was carried out using running and sliding window techniques. The lowest b-value (0.57 ± 0.04) and the highest D-value (1.65 ± 0.02) were computed at the time of the Gorkha earthquake, after which the b-value significantly increased to a maximum of 1.57. It again dropped to 0.93 at the time of the major aftershock on May 12, 2015. The D-value showed an initial quick drop and then decreased in a wavy pattern until the end of the study period, indicating the clustering and scattering of earthquakes in a fault region. The b-value contour map identified the eastern part of the study area as a high stress region (b = ~0.8), implying that the stress shifted to that region. The D-value contour map reveals that the seismogenic structure shifted from linear to planar in the region. The rate of aftershock decay (p = 0.86 ± 0.04) for a short period reflects that the level of stress decreased rapidly. This study helps to understand the level of stress and seismicity pattern of a region, which could be useful for aftershock studies.