The Eryuan area is located in the central part of the northwest Yunnan region on the southeastern edge of the Qinghai-Xizang Plateau. The geological structure in the area is complex, including the Weixi-Qiaohou-Weishan Fault and Honghe Fault, as well as the Longpan-Qiaohou Fault and Heqing-Eryuan Fault, which intersect in an “X” shape. The geothermal activity in the area is also active and exhibits strong fault structural characteristics. Earthquake activity is also frequent in the Eryuan area. Since 2013, multiple earthquakes with MS≥5.0 have occurred on the west side of the Weixi-Qiaohou-Weishan Fault. Given the complex geological and geothermal background and seismic activity in the Eryuan area, we use the seismic travel time data recorded by the Yunnan Regional Seismic Network and the Northwest Yunnan Dense Array from January 1, 2008, to July 20, 2023, VP/VS model consistency-constrained DD tomography method to obtain the 3D VP, VS, and VP/VS models and relocation results in the Eryuan and its surrounding areas. The research results indicate that: 1)based on the relocation results, the dense small seismic cluster develops at the intersection of the Weixi-Qiaohou-Weishan Fault, Honghe Fault, and Heqing-eryuan Fault deserves special attention. Four earthquakes with MS≥5.0 that have occurred since 2013 are mainly distributed on the west side of the Weixi-Qiaohou-Weishan Fault, with an NNW-SSE trend distribution. The fault system at the western boundary of the Sichuan-Yunnan block is very complex, and its seismic risk deserves our special attention. 2)Based on the tomography results, the widely distributed low-velocity anomalies in the upper crust of the study area may be related to the crustal material migration pathway caused by the escape of the Sichuan-Yunnan diamond block towards SSE. 3)Combining imaging results with geochemical research results for inference, the high VP/VS below 10km depth beneath the small seismic cluster at the intersection of the Weixi-Qiaohou-Weishan Fault, Honghe Fault, and Heqing-Eryuan Fault may correspond to hot spring geothermal fluids. Moreover, the circulation process of hot spring fluids in complex fault systems may have penetrated to a depth of 7~10km below the seismic cluster, accompanied by some small earthquakes. However, there is no imaging evidence of fluid within a depth of 5~7km where some earthquakes are densely distributed. The occurrence of these earthquakes may be related to the distribution of brittle rocks, but it cannot be ruled out that fluid may have a potential infiltration effect on them in the future. Based on the results of velocity structure, they are combined with pre-existing seismology research results. It is found that the Eryuan earthquake sequences on March 3 and April 17, 2013 were mainly located within low VP, low VS, and low VP/VS anomaly. Generally, if there is fluid present, VS decreases faster than VP, resulting in high VP/VS. So low VP/VS indicates that low VS is not caused by fluids, which may be due to lithology. Therefore, the imaging evidence indicates that there is no fluid in the region where the sequence is located, thus inferring that the occurrence of the sequence is not directly related to the fluid. The MS5.1 Yangbi earthquake sequence on March 27, 2017, is spatially close to the Eryuan sequence and also has the same velocity structure characteristics, so it may not be directly related to fluids. The mainshock and some aftershocks of the MS5.1 Yunlong earthquake sequence on May 18, 2016, were mainly located in high VP, high VS, and relatively high VP/VS regions. High VS indicates that it is not the fluid that causes relatively high VP/VS. It is speculated that there may not be fluid in the region where the mainshock is located, and the fluid did not directly participate in the occurrence process of the sequence.
Earthquake catalog is the foundational data for analyzing seismic activity, assessing seismic hazard, and studying earthquake prediction. The majority of historical earthquake records are sourced from historical documents, with a significant portion of these records found in local gazetteers. Compiling historical literature is an essential way in analyzing seismic activity because historical accounts of earthquakes often provide more detailed and accurate information than geological data. Among these sources, analyzing relevant content in local gazetteers, such as the historical development of local governance, military garrison, official records, and descriptions of disasters and auspicious events, plays a crucial role in seismic activity research. This article aims to acquire historical earthquake records by consulting local gazetteers, folk books, and other historical sources containing natural, social, and political records. These records serve as historical foundations for analyzing the completeness of seismic data records.
The border region between Sichuan, Yunnan, and Tibet is located in the northwest secondary block of the Sichuan-Yunnan block, which is one of the areas with frequent strong earthquakes in China. The Xianshuihe fault zone and the Jinshajiang fault zone are the northeastern and northwestern boundary faults of the Sichuan-Yunnan block, respectively. They are large-scale and highly active fault zones formed due to the eastward escape of the Tibetan plateau caused by the relative movement between the Indian and Eurasian plates. Previous studies on active tectonics have shown that major earthquakes with magnitudes of 8 and above, as well as over 80% of strong earthquakes with magnitudes of 7, mainly occur in the boundary zones of active blocks with intense structural deformation and high stress accumulation. Moreover, the known active faults in the study area, such as the Batang fault and Litang fault, are also major faults that significantly have influence on the occurrence strong earthquakes. The Sichuan-Yunnan-Tibet adjacent region is home to significant infrastructure, including the Sichuan-Tibet railway and hydropower stations. Analyzing the completeness of earthquake data in the border region of Sichuan, Yunnan, and Tibet can contribute to the assessment of fault hazards and the analysis of regional seismic activity trends. This, in turn, can help minimize the damage caused by earthquakes to critical infrastructure and further enhance the safety and security of people’s lives and properties.
This study reviewed the local gazetteers of 44 counties in the border region between Sichuan, Yunnan, and Tibet, and summarized the establishment and historical evolution of each county. Based on the analysis of the road evolution from Sichuan to Tibet and from Yunnan to Tibet, we examined the significant roles of important transportation hubs and nodes, such as stations, pond flood, and grain platforms, in regarding of recording earthquakes. Combining various historical sources and previous research on the completeness of earthquake data in the region, we conducted a comprehensive analysis to determine the probable starting years for the availability of seismic records of magnitude 7 and above in the Xianshuihe area and the three parallel rivers area. Additionally, based on the data of the length and short axis of isoseismal lines from 88 earthquakes, an elliptical model was used to derive the seismic intensity attenuation relationship for the Sichuan-Yunnan block. By placing the fitted isoseismal lines of magnitude 6 and 7 earthquakes in the study area, we analyzed their impact range, providing a spatial dimension basis for the completeness analysis of seismic data.
This article provides a comprehensive analysis and demonstration of the complete starting years of seismic data in the border region between Sichuan, Yunnan, and Tibet from both temporal and spatial perspectives. The results indicate that due to the establishment of grain stations and Tangxun along the Sichuan-Tibet road, as well as the appointment of officials, several counties in the Xianshuihe area, including Kangding, Luhuo, Garzê, Litang, and Yajiang, were developed between 1719 and 1736. At the same time, there are relatively abundant historical documents related to earthquakes in the Xianshuihe area. Local chronicles, reports from governors and resident ministers, written records in Tibetan temples, and accounts from lamas have documented earthquake surveys, disaster assessments, and relief efforts. By combining these historical sources with the analysis of intensity attenuation relationships in the Sichuan-Yunnan block, the affected areas of earthquakes with magnitudes 6 and 7 can be determined that the period from 1719 to 1736 marks the starting years with complete M≥7 earthquake data in the Xianshuihe area. The towns of Batang, Mangkang, and Changdu in the three parallel rivers area are also significant nodes and hubs along the road to Tibet. They were established with administrative institutions and granaries between 1719 and 1728, and the road network extensively covered Tangxun in the region. In considering the seismic records and historical sources in the three parallel rivers area, as well as referencing the recording capabilities of granaries, administrative institutions, and Tangxun in the Xianshuihe area, and estimating the potential recorded seismic magnitudes based on the intensity attenuation relationships of the Sichuan-Yunnan block, it can be suggested that the period from 1719 to 1728 is a possible starting point for complete earthquake data with M≥7 in the three parallel rivers area. In areas farther away from the road to Tibet, such as Jiangda, Gongjue, Baiyu, Xinlong, and the northern regions of Batang and Litang, as well as the large contiguous regions of Derong, Xiangcheng, Daocheng, and Jiulong, the eastern boundary is the Xianshuihe fault zone, while the area between the two zones is divided by the northeast-oriented Batang fault. Previous seismic geological investigations have found that within the aforementioned regions, the influence of the Jinshajiang fault zone extends along the Batang-Derong-Benzilan line. In remote areas away from the road and with sparse population, the possibility of individual earthquakes with magnitudes above 7 occurring but being missed cannot be ruled out. However, in other areas not located on active fault zones, it can be considered unlikely to experience earthquakes with magnitudes above 7. Based on the analysis of the data, the starting years of earthquakes with a magnitude of 6 and above should be the same as those of earthquakes with a magnitude of 7 and above. However, according to the analysis of the average occurrence rate of earthquakes per year, there is a significant lack of records for earthquakes of magnitude 6 and above. This may be due to the sparsely populated and vast nature of the Tibetan region during historical times, limited administrative capabilities of officials, and lack of earthquake historical records and documents. Therefore, it is not possible to determine the exact starting year for complete data on earthquakes of magnitude 6, which would be the same as for earthquakes of magnitude 7 and above.
14 survey lines, with a total of 167 measuring points, were laid out in the northern section of the Red River Fault, the Longpan-Qiaohou Fault, the Heqing-Eryuan Fault, the eastern piedmont fault of Yulong Mountains, and the Lijiang-Jianchuan Fault in the northwest of Yunnan Province, China. Cross-fault soil gas radon, hydrogen, and carbon dioxide have been measured on the above-mentioned faults. The concentration intensity and distribution characteristics of soil gas in the study area were calculated and analyzed. The results show that: (1)The concentrations and distribution patterns of soil gas radon and hydrogen vary greatly in different faults. The concentrations of radon vary from 6.18Bq/L to 168.32Bq/L, while that of hydrogen are between 7.72ppm to 429ppm, and carbon dioxide are from 0.73% to 4.04%. (2)The average results of soil gas measurement show that the concentrations of radon are higher than 40kBq/m3 in the sampling sites of Yinjie, Niujie, Gantangzi, while the concentrations of radon in En’nu and Tiger Leaping Gorge measuring lines are smaller; The concentrations of hydrogen are higher than 60ppm in the sampling sites of Yangwang village, Houqing, Dawa, Yangcaoqing and Tiger Leaping Gorge, while the concentrations in Gantangzi and Niujie measuring lines are smaller. (3)The spatial distribution characteristics of soil gas concentration in faults in northwest Yunnan are obvious, and the intensity of radon and hydrogen concentrations in different active fault zones vary greatly. The intensities of radon and hydrogen concentration are higher and have good consistency in Yinjie and Yangwang village measuring lines located in the northern section of the Red River Fault, the Houqing survey line located in Longpan-Qiaohou Fault, Dawa survey line in the Lijiang-Jianchuan Fault and Yangcaoqing in the south of Chenghai Fault. The soil gas concentration in such sample sites is high and the degassing ability is strong, indicating the different activity characteristics of different segments of the above faults to some extent. Under the action of tectonic stress, the fault will slip and the rock properties and material structure of the fault will change, thus causing changes of underground material, gas transport channel and transport mode, which is characterized by the change of the concentration and distribution characteristics of escaped soil gas. Combined with the active characteristics of faults, slip rate and geomorphological features, the characteristics of concentration and spatial distribution of two soil gases(radon and hydrogen)are discussed, and the following conclusions are obtained. (1)There is a large difference in the concentration of escaped soil gas from different faults in the study area, indicating that the content of soil gas is controlled by regional geochemical background values, and there are certain differences in the gas concentration of different sections of the same fault, indicating that the local concentration/flux change is greatly affected by the transport. (2)The concentration of fault soil gas is related to fault activity, and for different faults, the higher the degree of fault activity is, the higher the concentration of fault gas will be. From the point of fault gas concentration characteristics, the concentrations in the survey lines in the northern section of the Red River Fault and Heqing-Eryuan Fault in the study area vary greatly, suggesting that the fault segmentation is obvious. Compared with other faults, the northern section of the Red River Fault has a higher concentration of soil gas, indicating that the fault is more active. However, there is no simple linear relationship between the soil gas concentration and the fault slip rate, and it may also depend on its material source, transport channel structure, etc. (3)The concentration of soil hydrogen at the outcrop of faults(especially normal faults and strike-slip faults)is generally higher, which shows that hydrogen has better indicative significance in revealing the location of fault rupture, and the distribution pattern of soil gas radon concentration is a good indicator for analyzing the characteristics of fault movement.
Ocean earthquakes pose a serious threat to the security of the economic construction in coastal areas and the marine resource exploitation of China, so building proper seismicity models for China seas and adjacent regions is one of the focuses of the next generation seismic zoning map of China. Different from mainland China, there is lack of multidisciplinary data for the sea areas, such as seismic geology and geophysical exploration, thus the earthquake catalogue is the most important basic data for analyzing its seismic activity characteristics. Recently, a unified catalogue was compiled by researchers for China seas and neighboring regions, which provides a basis for further work. The Smoothed Seismicity Model (SSM) is a classic model based on earthquake catalogue and has become a basic model for the United States National Seismic Hazard Map. The Adaptively Smoothed Seismicity Model (ASSM) is an improved version of SSM. Compared with SSM, ASSM has optimized the value algorithm of the kernel function, thereby improving the model’s capability on mid- and long-term earthquake forecast. Due to the outstanding performance of ASSM in the CSEP project, it has become a research hotspot for seismologists. Based on a further improved ASSM algorithm and the newly compiled catalogue, this study established for the first time an adaptively smoothed seismicity model for China seas and adjacent areas. First, similar to most studies on mid-to-long-term seismicity models, we removed the foreshocks and aftershocks from the catalogue. Then we used seismic zones as the unit to estimate the start and end time of completely recorded earthquakes in different magnitude intervals. Furthermore, the maximum likelihood method was used to estimate the seismic activity parameters such as a-value and b-value for every seismic zone. On this basis, an improved adaptively smoothed algorithm was used to build the model. The function of probability gain per earthquake was applied to evaluate the performance of models under different parameter settings. Finally, our model was compared with the traditional SSM, and the advantages and limitations of our model were analyzed. Results show that: Compared with SSM, our model has a better performance on the probability gain function, and this advantage is not affected by the minimum magnitude(Mmin)of the input earthquakes. When Mmin is set as M4.0, our model achieves its best performance. The performance of the model does not necessarily improve as Mmin decreases or the number of earthquakes increases. This indicates that we need to comprehensively consider the distribution of earthquakes in the study area and the integrity of the earthquake catalogue when selecting parameters. The algorithm used in this study can make full use of seismic data with varying record level over time and space, and has a certain application value in seismic hazard analysis and mid- and long-term earthquake forecast. Meanwhile, our model can be used as one of the basic models to analyze the seismic hazard for China’s maritime areas, and provide technical support for the compilation of seismic zoning maps in China’s sea areas. In addition, the completeness analysis results and seismic activity parameters obtained by this study can also provide references for other peer research. Since ASSM is only based on historical and instrumental earthquake catalogues, it has certain limitations. For example, it cannot calculate the upper limit of the magnitude and reflect the distribution of faults, nor can it consider the time-dependence of the recurrence of large earthquakes. Therefore, this model can be used alone to describe the occurrence probability of small to moderately strong earthquakes, and it can also be used as one of the important factors to determine the spatial distribution functions for potential seismic source zones. However, when conducting seismic hazard analysis, it is recommended to combine the results of other disciplines such as geology and geodesy to form a hybrid model, so as to further improve the applicability and effectiveness of the model.
The seismogenic structure of the Lushan earthquake has remained in suspensed until now. Several faults or tectonics, including basal slipping zone, unknown blind thrust fault and piedmont buried fault, etc, are all considered as the possible seismogenic structure. This paper tries to make some new insights into this unsolved problem. Firstly, based on the data collected from the dynamic seismic stations located on the southern segment of the Longmenshan fault deployed by the Institute of Earthquake Science from 2008 to 2009 and the result of the aftershock relocation and the location of the known faults on the surface, we analyze and interpret the deep structures. Secondly, based on the terrace deformation across the main earthquake zone obtained from the dirrerential GPS meaturement of topography along the Qingyijiang River, combining with the geological interpretation of the high resolution remote sensing image and the regional geological data, we analyze the surface tectonic deformation. Furthermore, we combined the data of the deep structure and the surface deformation above to construct tectonic deformation model and research the seismogenic structure of the Lushan earthquake. Preliminarily, we think that the deformation model of the Lushan earthquake is different from that of the northern thrust segment ruptured in the Wenchuan earthquake due to the dip angle of the fault plane. On the southern segment, the main deformation is the compression of the footwall due to the nearly vertical fault plane of the frontal fault, and the new active thrust faults formed in the footwall. While on the northern segment, the main deformation is the thrusting of the hanging wall due to the less steep fault plane of the central fault. An active anticline formed on the hanging wall of the new active thrust fault, and the terrace surface on this anticline have deformed evidently since the Quaterary, and the latest activity of this anticline caused the Lushan earthquake, so the newly formed active thrust fault is probably the seismogenic structure of the Lushan earthquake. Huge displacement or tectonic deformation has been accumulated on the fault segment curved towards southeast from the Daxi country to the Taiping town during a long time, and the release of the strain and the tectonic movement all concentrate on this fault segment. The Lushan earthquake is just one event during the whole process of tectonic evolution, and the newly formed active thrust faults in the footwall may still cause similar earthquake in the future.
Focal mechanism solutions, distribution of aftershocks and field investigation suggest that the 3 August 2014 Ludian, Yunnan MS6.5 earthquake was spawned by the NW-trending steep left-slip fault that stretches along Baogunao-Xiaohe. The landsides, dominated by avalanches, induced by the shock are distributed in a rectangular area(15km×12km)with a NW-orientated long axis. The bedrock avalanches at many sites indicate that the primary direction of ground motion changes from SE in the north to SN in the south. The distribution of landslides triggered by the shock can be explained by the following two models of seismogenic structure: 1)Left-lateral strike-slip faulting occurred on a generally NW-trending arc-like fault, which caused seismic ground motion with the direction changing gradually from SE in the north to near-SN in the south; and 2)in addition to the left-lateral slip of the NW-trending fault, the NE-trending fault has also contributed to the hazard, which thrust from NW toward SE. So the Ludian MS6.5 event was the result of the joint action of both the NW and NE directed faults, of which the NW left-lateral slip fault is the dominant and the NE thrust fault is the secondary. Another evidence is that most aftershocks are distributed linearly in NW direction, meanwhile some are clustered in nearly NE around the epicenter, which also implies the possibility that both faults moved simultaneously during the event. Besides, focal mechanism solutions of a lot of earthquakes in the Yongshan-Zhaotong area, northeastern Yunnan Province, where the epicenter of the 2014 Ludian MS6.5 event is located, indicate that the seismogenic faults are primarily dominated by thrust faulting, associated with strike-slip faulting.
The Zhaotong-Ludian Fault zone, composed mainly of three right-step en echelon faults, namely, the Zhaotong-Ludian Fault, the Sayuhe Fault and the Longshu Fault, strikes 40°~60° on the whole, with the Sayuhe Fault and the Longshu Fault dipping SE and the Zhaotong-Ludian Fault dipping NW, and they all together constitute a complicated thrust fault system. Based on years of field investigation results of geology and geomorphography, we elaborate the late Quaternary active features, the geological and geomorphic evidences of the latest activity of the Zhaotong-Ludian Faults. Our observation shows that: the late Cenozoic basins along the Zhaotong-Luian Fault zone are obviously dominated by the fault; there are many neo-active fault landforms, such as, flat and straight fault troughs, directional aligned fault facets and fault scarps, and the upper Pleistocene to Holocene strata are offset by the fault. The fault zone has been active since the late Quaternary. For example, the fault at Daqiaobian dislocated a set of strata of the Pliocene, and middle to upper Pleistocene, with an apparently reverse character. The fault trending NE is developed in the Holocene diluvium with oblique striation on the fault plane at Guangming Village. Deposits with an OSL age of(23.4±1.8)ka BP on T2 terrace of a small river near Beizha town was offset by the fault. There is a fault scarp trending NE 40°, 0.5~2.0m in height, on the first terrace of the Longshu River near the Longshu Village. Several Quaternary faults are revealed by the trench which offset the late Pleistocene to Holocene strata and there are three poleo-earthquake events discovered in the trench. At Yanjiao Village the gravel layer has risen steeply and is aligned in a line because of squeezing effect of the fault; the rivers and ridges nearby are synchronously offset dextrally up to 30~40m. The fault zone is dominated by reverse faulting with a small amount of right-lateral motion. Besides, there are some NW-trending faults interweaving with the NE-trending fault zone, some of which are active since late Quaternary as well, and they are the conjugate structures with the NE-trending faults. Surface deformation, such as NE- and NW-trending ground fissures and reverse scarp landforms, has been generated in the epicenter area of the 2014 Ludian M6.5 earthquake, the distribution of which is in consistence with the NE- and NW-trending faults. Because of far-field deformation response and energy exchange and transfer between blocks, the Liangshan active sub-block formed on the east of the Sichuan-Yunnan block, and the Zhaotong-Ludian Fault zone lies in the forefront of the SE movement of this sub-block. On account of its distinct location and its complicated geometric structure, the Zhaotong-Ludian Fault zone is one of main carriers of the tectonic deformation of the Liangshan active sub-block to absorb and accommodate the strains produced by the block's SE movement, and is the southern boundary of the Liangshan sub-block. From the point of view of the regional tectonic positions and the kinematic characteristics, the relation of Zhaotong-Ludian Fault zone to the Liangshan active sub-block is exactly as the relation of the Longmanshan Faults to Bayan Har block. Consequently, the Zhaotong-Ludian Fault zone has an important significance in the division of active block boundaries and the regional tectonic framework, and meanwhile, it is also an important seismogenic structure in the northeastern Yunnan.
The June 24th 2012 Ninglang-Yanyuan MS 5.7 earthquake happened at 30km northwest of the Lijiang-Xiaojinhe Fault, a region seismically active in history and prone to earthquake in northwestern Yunnan. Tectonics in the earthquake region is complex,where two groups of faults are developed,trending NW and NE,respectively,and distributed in a chessboard pattern. Field survey results reveal that there are the NW-trending Yongning Fault and the NE-trending Rigulu-yanwa Fault developed near the epicenter,both are active in late Pleistocene.The Yongning Fault,composed of Wenquan Fault,Yongning Fault,and Alaao Fault,shows obvious fault landforms and clear linear features on satellite imagery.The fault has played an obvious control role in the development of the Yongning Quaternary Basin and Lugu Lake Basin,with several hot springs developed along the fault. Several tributaries of Qiansuo River run along the fault,and there are dextral displacements observed in many parts of the rivers along the fault,such as between Baqi and Haiyijiao,Shancuo village east of Rigulu. Near Alaao,the fault offset the late Pleistocene deposit on the T2 terrace,and the latest TL age of the offset stratum is (21.19±1.80)ka,indicating it is a normal with dextral strike-slip,late Pleistocene active fault.The Rigulu-Wayan Fault has played a noticeable control role in the development of the Tertiary Basins such as Wayan,Rigulu,and Lijiazui and the Quaternary Basin of Yongning. It offset the mid Pleistocene and the upper Pleistocene strata. Between Zhongwadu and Lijiazui,several streams were synchronously displaced left-laterally. There are obvious signs showing the fault was active in the late Pleistocene,dominated by sinistral strike-slip. According to the focal mechanism solutions,the Ninglang-Yanyuan MS 5.7 earthquake is of normal faulting with dextral strike-slip,the attitude of the NW nodal plane is basically consistent to the Yongning Fault,and the seismic rupture pattern is identical to the kinematical characteristic of Yongning Fault.The major axis of the isoseismals,the linear distribution of intensity Ⅷ anomaly sites and the direction of tectonic ground fissures are all consistent to the strike of Yongning Fault. Through analysis,it is believed that Yongning Fault is the seismogenic fault of Ninglang-Yanyuan MS 5.7 earthquake. Furthermore,the 1996 Lijiang M7.0 earthquake,the 1976 Zhongdian M5.5 earthquake and this M5.7 earthquake all have apparent normal dip-slip component. These earthquakes are located in the periphery of the neo-tectonic uplift of Haba Snow Mountain and Yulong Snow Mountain. Based on analyses of regional topography,the normal faulting in this area is most likely related to the gravitational potential energy resulting from the big topography contrast.