Search

The temporal distribution of earthquakes provides important basis for earthquake prediction and seismic hazard analysis. The relatively limited records of strong earthquakes have often made it difficult to study the temporal distribution models of regional strong earthquakes. However, there are hundreds of years of complete strong earthquake records in the North China Seismic Zone, providing abundant basic data for studying temporal distribution models. Using the data of M ≥ 6.5 earthquakes in North China as inputs, this paper estimates the model parameters using the maximum likelihood method with Poisson, Gamma, Weibull, Lognormal and Brownian passage time (BPT) distributions as target models. The optimal model for describing the temporal distribution of earthquakes is determined according to Akaike information criterion (AIC),and Kolmogorov–Smirnov test (K–S test). The results show that Lognormal and BPT models perform better in describing the temporal distribution of strong earthquakes in North China. The mean recurrence periods of strong earthquakes (M ≥ 6.5) calculated based on these two models are 8.1 years and 13.2 years, respectively. In addition, we used the likelihood profile method to estimate the uncertainty of model parameters. For the BPT model, the mean and 95% confidence interval of recurrence interval μ is 13.2 (8.9–19.1) years, and the mean and 95% confidence interval of α is 1.29 (1.0–1.78). For the Lognormal model, the mean value and 95% confidence interval of v is 2.09 (1.68–2.49), the mean value exp (v) corresponding to earthquake recurrence interval is 8.1 (5.4–12.1) years. In this study, we also calculated the occurrence probability of M ≥ 6.5 earthquakes in the North China Seismic Zone in the future, and found that the probability and 95% confidence interval in the next 10 years based on the BPT model is 35.3% (26.8%-44.9%); the mean value and 95% confidence interval of earthquake occurrence probability based on the Lognormal distribution is 35.4% (22.9%-49.7%); the mean probability and 95% confidence interval based on the Poisson model is 53.1% (41.1%-64%). The results of this study may provide important reference for temporal distribution model selection and earthquake recurrence period calculation in future seismic hazard analysis in North China. [ABSTRACT FROM AUTHOR]

Seismically triggered landslides can cause great damage to the road construction in mountainous areas. The permanent displacement analysis based on Newmark sliding-block model can evaluate risk of these landslides from the perspective of deformation damage and overall failure probability of slopes. However, the sliding-block model does not consider the attenuation effect of the shear strength on the sliding surface during earthquake, causing the calculated value of Jibson method to be less than the actual value. Therefore, the Newmark model was modified by adding attenuation coefficients to the effective internal friction angle and the effective cohesion of geologic units. The landslide areal density was proposed for hazard zoning with the Wenchuan earthquake data. The results showed that the predicted values agreed well with the real distribution of the landslides triggered by the Lushan earthquake. The proposed hazard zoning method in this paper can predict the severity of seismic landslides in consideration of the environmental changes in mountainous regions after the earthquake and provide support for the site selection in highly seismic areas. [ABSTRACT FROM AUTHOR]

In the Three Gorges reservoir region, central China, seismic activity increased substantially after the reservoir impoundment in 2003 which continues till date. Previous studies show that these are reservoir-triggered earthquakes and various factors are responsible for the increase in seismic activity after the reservoir impoundment. However, these studies do not provide a comprehensive assessment of influence of reservoir water level variations on spatiotemporal distribution of earthquakes. In this study, we statistically analyze the influence of the water level variations on the increased seismic activity in the reservoir region for the period from May 2003 to April 2020, using the power spectrum and singular spectrum techniques. Our statistical analyses confirm the influence of long-term variations in the water level time series on the occurrence of earthquakes after the reservoir impoundment. The analysis also indicates a positive role of annual reservoir water level fluctuations in the total seismicity of the region. Depending on the cluster patterns and relationship with faults, the earthquakes of the Three Gorges reservoir region are divided into three seismic zones (A, B, and C). For zone C, both the power spectrum and singular spectrum analyses confirm the strong periodic influence of reservoir water level variations on the earthquakes. Increase in seismicity of zone B is only in the initial period but not in the later stages of water impoundment and our statistical analyses indicate that the seismicity of this zone is not directly related to the annual reservoir water level variations. This confirms the conjecture in the earlier studies that the seismicity of this zone is related to the collapse of coal mines present in the area in the initial stages of reservoir impoundment. For zone A, our statistical analyses do not show strong influence of the annual reservoir water level variations on the occurrence of earthquakes. We suggest that this is due to the contribution of various other factors along with reservoir impoundment in the occurrence of earthquakes in this zone, as also opined in some earlier studies. [ABSTRACT FROM AUTHOR]

This paper analyzes the seismicity in Bohai Sea, introducing a shape factor K to characterize the seismic risk distribution in sub-regions of the sea. Based on the seismic design ground motions for 46 platforms located in the Bohai Sea, a statistical analysis was performed for different peak ground acceleration (PGA) ratios at two different probability levels. In accordance with the two-stage design method, a scheme of two seismic design levels is proposed, and two seismic design objectives are established respectively for the strength level earthquake and the ductility level earthquake. By analogy with and comparison to the Chinese seismic design code for buildings, it is proposed that the probability level for the strength level earthquake and ductility level earthquake have a return period of 200 and 1000–2500 years, respectively. The validity of these proposed values is discussed. Finally, the PGAs corresponding to these two probability levels are calculated for different sub-regions of the Bohai Sea. [ABSTRACT FROM AUTHOR]

This paper aims at implementing and introducing the use of conditional mean spectrum (CMS) and conditional spectrum (CS) as the main input parameters in the practice of seismic safety evaluation (SSE) in China, instead of the currently used uniform hazard spectrum (UHS). For this purpose, a procedure for M-R-epsilon seismic hazard deaggregation in China was first developed. For illustration purposes, two different typical sites in China, with one to two dominant seismic zones, were considered as examples to carry out seismic hazard deaggregation and illustrate the construction of CMS/CS. Two types of correlation coefficients were used to generate CMS and the results were compared over a vibration period range of interest. Ground motion records were selected from the NSMONS (2007-2015) and PEER NGA-West2 databases to correspond to the target CMS and CS. Hazard consistency of the spectral accelerations of the selected ground motion records was evaluated and validated by computing the annual exceedance probability rate of the response spectra and comparing the results to the hazard curve corresponding to each site of concern at different periods. The tools developed in this work and their illustrative application to specific case studies in China are a first step towards the adoption of CMS and CS into the practice of seismic safety evaluation in this country. [ABSTRACT FROM AUTHOR]

Detailed and accurate historical earthquake disaster data are an important basis for understanding the characteristics of earthquake disaster losses and for conducting research on earthquake disaster prevention strategies and post-earthquake emergency rescue work. Due to different project scopes and discrepancies in data accuracy, the currently available databases may contain errors such as missing data or data with low accuracies. In this study, we developed the Mainland China Composite Earthquake causing Casualties Database (MCCEC-DAT), which contains 252 historical earthquake records that span from 1966 to 2021. The database contains earthquake attribute information (earthquake time, location, epicenter latitude and longitude, focal depth, epicentral intensity, etc.), loss information (number of deaths and injuries, building damage ratio, economic loss, cause of death, etc.), and socioeconomic data (population, population density, area of each zone of seismic intensity, mortality rate, etc.). This database also contains the lethality levels generated by this study. These levels are the result of the quantitative evaluation of the overall characteristic levels of the earthquake area. Furthermore, this metric permits the direct comparison of different spatial locations. Based on the results of the statistical analysis, this work shows that the distribution of historical earthquakes in mainland China has obvious spatial and temporal characteristics. Events are mainly distributed throughout the southwest region, and the magnitudes of those events are relatively diverse. The mortality rate is significantly correlated with the building damage ratio and the lethality level. The lethality level can not only facilitate an improved understanding of the comprehensive level of historical earthquake regions, but this metric also provides a quantitative foundation for the comparative analysis of earthquakes between different regions. Furthermore, the database can also provide a foundation for future research, such as research into assessing the loss from earthquake disasters. [ABSTRACT FROM AUTHOR]

As the largest water conservancy and hydropower project in China, the Three Gorges Reservoir is a weak seismic activity area before impoundment, but the frequency of earthquakes increases significantly after impoundment. The spatial density scanning method was used to obtain the characteristics of spatio-temporal earthquake distribution in the reservoir area during loading and unloading processes. The results show that the frequencies of earthquakes during the loading and unloading processes were higher than that during the low-water-level operation period, which is well explained by the acoustic emission test results. The seismic b-value, fractal dimension D, and spatial correlation length SCL can be used together to indicate stress criticality. To analyze the impacts of reservoir water loading and unloading on seismicity in the reservoir area, time-scan analyses were performed on the b-value, D-value, and SCL of earthquakes near the Zigui segment and the Badong segment. Previous studies argued that the time-varying characteristics of b-values do not hold predictive significance for earthquakes in the M4.0–6.2 range. However, our study found that the time-varying characteristics of b-values are of predictive significance for earthquakes around M4.0. These seismic parameters decrease significantly before moderate earthquakes but at different rates in different regions. [ABSTRACT FROM AUTHOR]

Based on the existing materials of fault segmentation, characteristic earthquakes, and their empirical relationships, we calculated the parameters of the fault segments, such as length, width, magnitudes of characteristic earthquakes, etc. Constrained by GPS velocity field, the slip rates of these fault segments in depth were inversed using the 3-D half-space elastic dislocation model. As not all of the recurrence periods and co-seismic displacements of characteristic earthquakes are known, we selected the fault segments with these two parameters known and calculated the accumulation rate of average co-seismic displacement, which shows the faults' slip rate in seismogenic layer. Then, the slip rate in depth was compared with that in seismogenic layer, the relationship between them was obtained, and this relationship was used to get the recurrence periods and co-seismic displacements of all fault segments. After the studies above, we calculated the co-seismic deformation field of all the earthquakes larger than M6.8 from AD 1700 one by one and inversed the potential displacement in the co-seismic deformation field. Then, we divided the potential displacement by the slip rate from GPS inversion to get the influences of these fault segments, added the influences into the elapsed time of the characteristic earthquakes, and obtained the earthquake hazard degree of all the segments we studied in the form of the ratio of elapsed time to recurrence period; so, we name the ratio as the Impending Earthquake Risk (IER). Historical earthquake cases show that the fault segment is in safety when the IER is less than 1 but in danger after the IER becomes larger than 1. In 2009, the IER is larger than 1 on the following segments, 1.35 on the Tagong segment of Xianshuihe fault, 1 on the Menggu-Dongchuan segment, 1.04 on the Dongchuan-Xundian segment, and 1.09 on the Yiliang-Chengjiang segment of Xiaojiang fault. [ABSTRACT FROM AUTHOR]

The East Kunlun Fault (EKF) is one of the largest and most active sinistral strike-slip faults in the Qinghai–Tibet Plateau. Since 1900, many earthquakes of magnitude > 7 have occurred in this region. The Maqin–Maqu Segment (MMS), a seismic gap underlying a densely populated zone, is particularly important. While the area has been extensively researched to assess seismic risk, refined asperity geodetic imaging is deficient. We combined multiple sets of GPS and leveling observations to obtain spatial–temporal distributions of locking along the EKF. We delineated the probable asperity ranges according to the fault plane locking and the distribution of small earthquakes. Long-term observations on GPS velocity (1999–2013) and leveling (1970–2013) indicate that the Dulan–Maduo segment (DMS) of the EKF has a 140-km-long strong coupling section, and both the MMS and the eastern Ruoergai segments have a 40-km-long locking range. The 2013–2015 GPS velocity observations show that from Maqu to Maqin, the locked DMS and MMS lengths increased to 200 km and 240 km, respectively, while the locked eastern Ruoergai segment decreased and migrated 60 km to the west. This change in locking over time implies that the post-seismic effects of surrounding earthquakes have had significant impacts on the EKF. The locking area derived from 2015 to 2017 GPS velocity shows a growing trend and indicates the connectedness between the two speculative asperities. Combining the distribution of historical earthquakes and the precision positioning of small earthquakes, we conclude that the geodetic imaging of the EKF indicates that the MMS and DMS are high-risk zones for destructive earthquakes. [ABSTRACT FROM AUTHOR]

The conducted paleoseismological and archaeoseismological studies reveal three zones of concentration of the ancient and historical earthquakes in Gorny Altai which are related to the Kurai Fault zone, Katun, and South Terekta faults. The surface ruptures are detected within the Kurai Fault zone, which were formed in the epicentral zones of the paleoearthquakes that occurred 6500, 5800, 3200, and 1300 years ago and had magnitudes Mw = 6.7–7.6. The recurrence period of the paleoearthquakes is 700 to 2600 years. The detected secondary seismogenic deformations indicate that an epicentral zone of the paleoearthquake with an age of less than 12.5 ka (Mw = 7.2–7.6, intensity I = 10–11), the traces of earthquakes and their clusters with M ≥ 5–5.5 and I ≥ 6–7, which occurred about 150 and 90 ka ago, in the intervals of 38–19 ka ago (with a recurrence period of about 2 ka), and 19–12.5 ka ago are related to the southern part of the Katun Fault. The earthquake of I ≥ 5–6 which damaged the constructions of the Chultukov Log 1 burial mound in the period from IV century B.C. to the beginning of I century A.D. is associated with the northern part of the Katun Fault. In the zone of the South Terekhta Fault, the seismogenic displacements that occurred in VII–VIII centuries A.D. (Mw = 7.4–7.7, I = 9–11) and about 16 ka ago (M ≥ 7, I = 9–10) are revealed. The latter triggered the formation of a landslide-dammed lake which was destroyed by the earthquake about 6 ka ago (M ≥ 7, I = 9–10). Secondary paleoseismic deformations of the ancient earthquakes (M ≥ 5–5.5, I ≥ 6–7) are recorded in the sediments of the Uimon Basin with an age of 100–90 ka and about 77 ka. These results should be taken into account in designing a gas pipeline in the People's Republic of China, building infrastructure for tourism, and elaborating the seismic zoning maps for the territory of the Russian Federation. [ABSTRACT FROM AUTHOR]

We observed abnormal changes of the water temperature in Yushu (YSWT) well, China, most of which were followed by earthquakes. This study statistically analyzes the correlation between the magnitude and duration of the anomalies in YSWT and earthquakes in the Tibetan block and its margins. The effectiveness of using observed YSWT data to predict earthquakes was quantitatively examined by the Molchan error diagram method. The results show that (1) the YSWT underwent several abnormal changes marked by “V”-shaped patterns, which might be related to several earthquakes that occurred in the Tibetan block and its margins. The extent and duration of the abnormal changes in the YSWT were linearly related to the magnitude of the earthquake; i.e., the higher the magnitude, the greater the change in the YSWT, and the shorter the duration. (2) Abnormal changes in the YSWT are somewhat predictive of earthquakes with magnitudes ≥5.5 (≥M5.5) within 800 km of the Yushu well and ≥M6.5 earthquakes in the Tibetan block and its margins. The prediction has a probability gain of approximately 2, and the most likely time period for an earthquake to occur is within approximately 3 months after the occurrence of an YSWT anomaly. Most of the anomalies in YSWT appeared before earthquakes in the thrust of block margins. Notably the larger strains from the earthquake did not produce any response. We speculate that the preseismic responses reflect the regional tectonics, such as the motion of the Indian plate, straining sub-blocks of the Tibetan block. [ABSTRACT FROM AUTHOR]

Coulomb stress accumulation and releasing history and its relationship with the occurrence of strong historical earthquakes could deepen our understanding of the occurrence pattern of strong earthquakes and hence its seismic potential in future. The sinistral strike-slip Xianshuihe- Xiaojiang fault zone (XXFS) is one of the most dangerous fault zones in China, extending 1 500-km-long from the central Tibetan Plateau to the Red River fault zone. There are 35 M≥6.5 historical earthquakes occurred since 1327, hence it is an ideal site for studying the Coulomb stress evolution history and its relationship with the occurrences of strong earthquakes. In this study, we evaluated the Coulomb stress change history along the XXFS by synthesizing fault geometry, GPS data and historical earthquakes. Coulomb stress change history also revealed different patterns of historical earthquakes on different segments of the XXFS, such as characteristic recurrence intervals along the Salaha-Moxi fault and super-cycles along the Xianshuihe fault. Based on the occurrence pattern of past historical earthquakes and current Coulomb stress field obtained in this study, we suggest positive ΔCFS and hence high seismic potential along the Salaha-Moxi fault and the Anninghe fault. [ABSTRACT FROM AUTHOR]

On August 8, 2017, a M7.0 earthquake occurred in Jiuzhaigou County, Sichuan Province, China, resulting in significant casualties and property damage. Therefore, it is critical to identify the areas of potential aftershocks before reconstruction and re-settling people to avoid future disasters. Based on the elastic dislocation theory and a multi-layered lithospheric model, we calculate the Coulomb failure stress changes caused by the Wenchuan and Jiuzhaigou earthquakes, discuss the relationship between the M 7.9 Wenchuan and M7.0 Jiuzhaigou earthquakes, and analyze the influence of the aftershock distribution and stress changes on the major faults in this region caused by the Jiuzhaigou earthquake. The co- and post-seismic stress changes caused by the Wenchuan earthquake significantly increased the stress accumulation at the hypocenter of the Jiuzhaigou earthquake. Therefore, the occurrence of the Jiuzhaigou earthquake was probably stimulated by the Wenchuan earthquake. The aftershock distribution is well explained by the co-seismic stress changes of the Jiuzhaigou earthquake. The stress accumulation and corresponding seismic hazard on the Maqu-Heye segment of the East Kunlun fault and the northern extremity of the Huya fault has been further increased by the Jiuzhaigou earthquake. [ABSTRACT FROM AUTHOR]

From the 1960s to 1970s, North China has been hit by a series of large earthquakes. During the past half century, geophysicists have carried out numerous surveys of the crustal and upper mantle structure, and associated studies in North China. They have made significant progress on several key issues in the geosciences, such as the crustal and upper mantle structure and the seismogenic environment of strong earthquakes. Deep seismic profiling results indicate a complex tectonic setting in the strong earthquake areas of North China, where a listric normal fault and a low-angle detachment in the upper crust coexist with a high-angle deep fault passing through the lower crust to the Moho beneath the hypocenter. Seismic tomography images reveal that most of the large earthquakes occurred in the transition between the high- and low-velocity zones, and the Tangshan earthquake area is characterized by a low-velocity anomaly in the middle-lower crust. Comprehensive analysis of geophysical data identified that the deep seismogenic environment in the North China extensional tectonic region is generally characterized by a low-velocity anomalous belt beneath the hypocenter, inconsistency of the deep and shallow structures in the crust, a steep crustalal-scale fault, relative lower velocities in the uppermost mantle, and local Moho uplift, etc. This indicates that the lithospheric structure of North China has strong heterogeneities. Geologically, the North China region had been a stable craton named the North China Craton or in brief the NCC, containing crustal rocks as old as ~3.8 Ga. The present-day strong seismic activity and the lower velocity of the lower crust in the NCC are much different from typical stable cratons around the world. These findings provide significant evidence for the destruction of the NCC. Although deep seismic profiling and seismic tomography have greatly enhanced knowledge about the deep-seated structure and seismogenic environment, some fundamental issues still remain and require further work. [ABSTRACT FROM AUTHOR]

Instrumentally recorded seismicity from 1970/01/01 to 2014/01/01 of the central China north-south seismic belt (21.0°-41.5°N, 97.5°-107.5°E) was analyzed, emphasizing the applicability of the predictive algorithms based on the assumptions of meta-stable equilibrium. The seismicity in this region was shown to exhibit ergodicity from 1980 to the present, with sub-region dependence, and interrupted by the 2008 Wenchuan earthquake. pattern informatics algorithm, a statistical physics-based predictive model for five-year time scale, is put to forward forecast test for the period 2014/01/01 to 2019/01/01. [ABSTRACT FROM AUTHOR]

After the 20 April 2013 Lushan M 6.6 earthquake occurred, investigation and identification of the seismogenic fault for this event have become a focused and debatable issue. This work prepared an initial landslide inventory map related to the Lushan earthquake based on field investigations and visual interpretation of high-resolution aerial photographs and provided evidence for solving the issue aforementioned. The analysis of three landslide-density profiles perpendicular to strike direction of the probable seismogenic fault shows that many landslides occurred on the footwall of the Shuangshi-Dachuan fault (SDF), without sudden change of landslide density near the fault. Very few landslides were detected near the Dayi fault (DF) and also no change of landslide density there. While obvious sudden change of landslide density appeared about 1-2 km from the northwest to the western Shangli fault (WSF), and the landslide density on the hanging wall of the fault is obviously higher than that of on the footwall. Therefore, we infer that the seismogenic fault for the Lushan earthquake is neither the SDF nor the DF, rather probably the WSF located between these two faults, which is an evident linear trace on the earth surface. Meanwhile, the coseismic slip did not propagate upward to the ground, implying the Lushan earthquake was spawned by a blind-thrust-fault beneath the WSF. [ABSTRACT FROM AUTHOR]

Following the Lushan M7.0 earthquake on 20 April 2013, a topic of much concern is whether events of M7 or greater could occur again on the southern segment of the Longmenshan fault zone. In providing evidence to answer this question, this work analyzes the tectonic relationship between the Lushan event and the 2008 Wenchuan earthquake and the rupture history of the southern segment of the Longmenshan fault zone, through field investigations of active tectonics and paleoearthquake research, and our preliminary conclusions are as follows. The activity of the southern segment of the Longmenshan fault zone is much different to that of its central section, and the late Quaternary activity has propagated forward to the basin in the east. The seismogenic structure of the 2008 Wenchuan earthquake is the central-fore-range fault system, whereas that of the 2013 Lushan event is attributed to the fore-range-range-front fault system, rather than the central fault. The southern segment of the Longmenshan fault zone becomes wider towards the south with an increasing number of secondary faults, of which the individual faults exhibit much weaker surface activity. Therefore, this section is not as capable of generating a major earthquake as is the central segment. It is most likely that the 2013 earthquake fills the seismic gap around Lushan on the southern segment of the Longmenshan fault zone. [ABSTRACT FROM AUTHOR]

Using global positioning system (GPS) technology, significant postseismic surface displacements were observed within the first 4 months after the 2001 Mw 7.8 Kunlun earthquake which occurred in China. In this study, we investigated the mechanisms that may have possibly contributed to the postseismic deformations that have been observed. Based on the modeling results, we find that an afterslip model can interpret postseismic displacements in the near field even when the fault plane is extended to the bottom of the crust (~70 km). Models based on the viscoelastic relaxation theory showed a large discrepancy in the spatial pattern of the deformation compared with what has been observed. Thus, we infer that both mechanisms cannot interpret the observed postseismic deformation independently. A combination of afterslip and viscoelastic relaxation can further improve the data fit, especially at sites far from the fault. With maximum afterslip of ~0.4 m occurring at a depth of 10 km in the central section, the combined model shows that the estimated afterslip occurred mostly on and below the coseismic rupture plane, as well as on its eastern extension. The estimated moment released by the afterslip in the first 4 months is almost 40% of that released by the coseismic slip. The best-fitting viscoelastic relaxation model shows a 'weak' upper mantle with a viscosity of ~1.0 × 10 Pa s. The combined model also suggests the existence of a lower crust with viscosity larger than 1.0 × 10 Pa s, although it cannot be constrained accurately. [ABSTRACT FROM AUTHOR]

— For the realistic modeling of the seismic ground motion in lateral heterogeneous anelastic media, the database of 3-D geophysical structures for Beijing City has been built up to model the seismic ground motion in the City, caused by the 1976 Tangshan and the 1998 Zhangbei earthquakes. The hybrid method, which combines the modal summation and the finite-difference algorithms, is used in the simulation. The modeling of the seismic ground motion, for both the Tangshan and the Zhangbei earthquakes, shows that the thick Quaternary sedimentary cover amplifies the peak values and increases the duration of the seismic ground motion in the northwestern part of the City. Therefore the thickness of the Quaternary sediments in Beijing City is the key factor controling the local ground effects. Four zones are defined on the base of the different thickness of the Quaternary sediments. The response spectra for each zone are computed, indicating that peak spectral values as high as 0.1 g are compatible with past seismicity and can be well exceeded if an event similar to the 1697 Sanhe-Pinggu occurs. [ABSTRACT FROM AUTHOR]

On 3 July 2015, the Mw 6.5 Pishan earthquake occurred at the junction of the southwestern margin of the Tarim Basin and the northwestern margin of the Tibetan Plateau. To understand the seismogenic mechanism and the post-seismic deformation behavior, we investigated the characteristics of the post-seismic deformation fields in the seismic area, using 9 Sentinel-1A TOPS synthetic aperture radar (SAR) images acquired from 18 July 2015 to 22 September 2016 with the Small Baseline Subset Interferometric SAR (SBAS-InSAR) technique. Postseismic LOS deformation displayed logarithmic behavior, and the temporal evolution of the post-seismic deformation is consistent with the aftershock sequence. The main driving mechanism of near-field post-seismic displacement was most likely to be afterslip on the fault and the entire creep process consists of three creeping stages. Afterward, we used the steepest descent method to invert the afterslip evolution process and analyzed the relationship between post-seismic afterslip and co-seismic slip. The results witness that 447 days after the mainshock (22 September 2016), the afterslip was concentrated within one principal slip center. It was located 5–25 km along the fault strike, 0–10 km along with the fault dip, with a cumulative peak slip of 0.18 m. The 447 days afterslip seismic moment was approximately 2.65 × 1017 N m, accounting for approximately 4.1% of the co-seismic geodetic moment. The deep afterslip revealed that a creeping process from steady-state "secondary" creeping to accelerating "tertiary" creep in the deep of fault. The future seismic hazard deserves further attention and research. [ABSTRACT FROM AUTHOR]