Table of Content

20 February 2019, Volume 41 Issue 1

Previous Issue    Next Issue

RAYLEIGH WAVE TOMOGRAPHY IN THE CRUST AND UPPER MANTLE OF THE DABIE-TANLU OROGENIC ZONE

XIONG Cheng, XIE Zu-jun, ZHENG Yong, XIONG Xiong, AI San-xi, XIE Ren-xian

2019, 41(1):  1-20.  DOI: 10.3969/j.issn.0253-4967.2019.01.001

Asbtract ( )   HTML   PDF (11337KB) ( )  

References | Related Articles | Metrics

Dabie Orogen has a series of special tectonic and geological features which make it important to the study of the tectonic evolution of mainland China and East Asia. The distribution of high pressure/ultra-high pressure metamorphic rocks discovered on the surface, the specific location of a series of deep and shallow sutures in the collisional convergence collage, and the seismogenic environment of shallow earthquakes attract many scientists continuously to study the interesting Dabie Orogen.
In this paper, we used waveform records of 200 broadband seismic stations deployed by China Digital Seismograph Network and vertical component records of 21 mobile seismic stations located in the Dabie-Tanlu orogenic zone and its surrounding areas. Based on seismic ambient noise tomography, we have obtained the phase velocity distributions of Rayleigh surface wave with the periods between 8~40s, with the resolution higher than 50km. The high velocity anomalies are observed on the Hong'an-Tongbo region in the images of 8~16s phase velocity, which decreases with increasing periods. These high velocity anomalies are in consistence with the ultra-high pressure(UHP)metamorphic rocks of the region. It leads us an estimation of the extension of UHP metamorphic rocks at various depths. The distribution of these anomalies found in phase velocity maps of 8s to 16s indicates that the estimated depth is up to~20km. The horizontal distribution forms a heart shape, which is narrower on western side and wider on the eastern side. It is very much consistent with the surface observations. The whole shape is similar to a cone that laterally extends its wings on the southwest. It indicates that the high-pressure/ultra-high pressure metamorphic rocks had experienced quick exhumation after they broke and formed a drag at the tail, and the residual area formed by the fast exhumation was likely to be invaded by magma. We agree that it has experienced complex structural history, such as stretching, magmatic emplacement and tectonic extrusion, resulting in the high-pressure/UHP metamorphic rocks finally exhuming on the surface with the structural pattern of narrower on the western margin and wider on the eastern margin in the Hong'an-Tongbo area.
The significant phase velocity difference from the period of 8s to 35s on both sides of the southern Tanlu fault zone enables us to infer that the Tanlu fault zone is a deep and huge fault, and the entire crust of the eastern zone of Dabie was cut by the Tanlu fault zone. It demonstrates that the Dabie block is separated from the northern Subei Basin and southern Yangtze blocks, which forms a seismogenic environment suitable for the generation of small-to-intermediate earthquakes in this region. Most of earthquakes in Anhui and adjacent provinces are distributed in those areas where the phase velocities changed dramatically, which are in consistence with the small faults of the upper crust in shallow layers of the Dabie-Tanlu orogenic belt. The shallow-source earthquakes mainly occur in velocity contrast regions, as demonstrated by the short period images. Earthquakes distribution and velocity maps show that the possible distribution of tiny faults of the upper crust can be roughly inferred from the geological structure. It helps to understand the seismogenic environment and seismic hazard in the Dabie areas. We conclude that the shallower faults with different velocity on either side of this region are still seismically active. These results have important significance for understanding the tectonic activity of the research areas.More detail work and further discussion are needed on the velocity structure of the Dabie orogen.

THE MAXIMUM LIKELIHOOD ESTIMATION OF b-VALUE IN MAGNITUDE-FREQUENCY RELATION AND ANALYSIS OF ITS INFLUENCING FACTORS

WU Guo, ZHOU Qing, RAN Hong-liu

2019, 41(1):  21-43.  DOI: 10.3969/j.issn.0253-4967.2019.01.002

Asbtract ( )   HTML   PDF (1421KB) ( )  

References | Related Articles | Metrics

b-value in the magnitude-frequency(G-R)relationship plays a vital role in seismicity research and seismic hazard analysis, and the most commonly used techniques to simulate it are least square approach and maximum likelihood method. Least square method is simple and easy to apply, therefore widely used in China. However, many researches show that there exist some limits in least square estimation of b-value. Earthquakes with different magnitudes are not equally weighted in this method, and larger events have higher weights, so b-value is vulnerable to the fluctuation of several big earthquakes; meanwhile, least square method needs to divide magnitude intervals artificially. With a small sample size, data points could be not enough if the magnitude interval is too wide, and events in a magnitude interval may be lacking if it is divided to be too narrow. Especially for incremental G-R relationship, it is possible that N(M_i)equals 0 in an interval with large magnitude, so log(N(M_i))loses meaning and has to be ignored, resulting in a low b-value. Therefore, under certain conditions, maximum likelihood method is recommended as an effective substitution or supplementary for least square estimation of b-value. Among numerous previous researches on maximum likelihood estimation of b-value, lots of equations have been provided, based on varied implicit assumptions and different ways of solution. A brief overview is first presented for these equations, and classification and summary are provided based on whether taking account of the effect of binned magnitude, with finite maximum magnitude, using unequal observation periods for different magnitude intervals, and with analytic solution or not. Following this, a total of 6 influential factors are analyzed, such as binning magnitude, measurement errors of magnitude, sample size, magnitude span, minimum completeness magnitude and fore- and aftershocks. At last, reasonable suggestions are provided for using those equations properly. The equations of Aki(1965), Utsu(1965), Page(1968)and Kijko and Smit(2012)are based on assumption that magnitudes are continuous random variables, and have no corrections for this, so these equations are not recommended here. For simplicity, the equations of Utsu(1966)or Tinti and Mulargia(1987)can be used, but magnitude span should be greater than 2.5 due to without finite maximum magnitude in the formulas. For researchers having capability to write code and calculate numerically, Weichert(1980)or Bender(1983)'s algorithm could be utilized. Especially when it is required to apply data with different observation periods for varied magnitudes, the formula of Weichert(1980)is recommended. This study contributes to more accurately understand and use different formulas of estimating b-value by maximum likelihood technique, which can be used as reference for peers.

MECHANISM OF THE 2016 HUTUBI, XINJIANG, M_S6.2 MAINSHOCK AND RELOCATION OF ITS AFTERSHOCK SEQUENCES

XU Zhi-guo, LIANG Shan-shan, ZOU Li-ye, LIU Jing-guang

2019, 41(1):  44-57.  DOI: 10.3969/j.issn.0253-4967.2019.01.003

Asbtract ( )   HTML   PDF (5884KB) ( )  

References | Related Articles | Metrics

A strong earthquake with magnitude M_S6.2 hit Hutubi, Xinjiang at 13:15:03 on December 8th, 2016(Beijing Time). In order to better understand its mechanism, we performed centroid moment tensor inversion using the broadband waveform data recorded at stations from the Xinjiang regional seismic network by employing gCAP method. The best double couple solution of the M_S6.2 mainshock on December 8th, 2016 estimated from local and near-regional waveforms is strike:271°, dip:64ånd rake:90° for nodal plane I, and strike:91°, dip:26ånd rake:90°for nodal plane Ⅱ; the centroid depth is about 21km and the moment magnitude(M_W)is 5.9. ISO, CLVD and DC, the full moment tensor, of the earthquake accounted for 0.049%, 0.156% and 99.795%, respectively. The share of non-double couple component is merely 0.205%. This indicates that the earthquake is of double-couple fault mode, a typical tectonic earthquake featuring a thrust-type earthquake of squeezing property.The double difference(HypoDD)technique provided good opportunities for a comparative study of spatio-temporal properties and evolution of the aftershock sequences, and the earthquake relocation was done using HypoDD method. 486 aftershocks are relocated accurately and 327 events are obtained, whose residual of the RMS is 0.19, and the standard deviations along the direction of longitude, latitude and depth are 0.57km, 0.6km and 1.07km respectively. The result reveals that the aftershocks sequence is mainly distributed along the southern marginal fault of the Junggar Basin, extending about 35km to the NWW direction as a whole; the focal depths are above 20km for most of earthquakes, while the main shock and the biggest aftershock are deeper than others. The depth profile shows a relatively steep dip angle of the seismogenic fault plane, and the aftershocks dipping northward. Based on the spatial and temporal distribution features of the aftershocks, it is considered that the seismogenic fault plane may be the nodal plane I and the dip angle is about 271°. The structure of the Hutubi earthquake area is extremely complicated. The existing geological structure research results show that the combination zone between the northern Tianshan and the Junggar Basin presents typical intracontinental active tectonic features. There are numerous thrust fold structures, which are characterized by anticlines and reverse faults parallel to the mountains formed during the multi-stage Cenozoic period. The structural deformation shows the deformation characteristics of longitudinal zoning, lateral segmentation and vertical stratification. The ground geological survey and the tectonic interpretation of the seismic data show that the recoil faults are developed near the source area of the Hutubi earthquake, and the recoil faults related to the anticline are all blind thrust faults. The deep reflection seismic profile shows that there are several listric reverse faults dipping southward near the study area, corresponding to the active hidden reverse faults; At the leading edge of the nappe, there are complex fault and fold structures, which, in this area, are the compressional triangular zone, tilted structure and northward bedding backthrust formation. Integrating with geological survey and seismic deep soundings, the seismogenic fault of the M_S6.2 earthquake is classified as a typical blind reverse fault with the opposite direction close to the southern marginal fault of the Junggar Basin, which is caused by the fact that the main fault is reversed by a strong push to the front during the process of thrust slip. Moreover, the Manas earthquake in 1906 also occurred near the southern marginal fault in Junggar, and the seismogenic mechanism was a blind fault. This suggests that there are some hidden thrust fault systems in the piedmont area of the northern Tianshan Mountains. These faults are controlled by active faults in the deep and contain multiple sets of active faults.

FOCAL FAULTS AND STRESS FIELD CHARACTERISTICS OF M7.0 JIUZHAIGOU EARTHQUAKE SEQUENCE IN 2017

LI Jun, WANG Qin-cai, CUI Zi-jian, LIU Geng, ZHOU Lin, LU Zhen, ZHOU Hui

2019, 41(1):  58-71.  DOI: 10.3969/j.issn.0253-4967.2019.01.004

Asbtract ( )   HTML   PDF (5579KB) ( )  

References | Related Articles | Metrics

On August 8, 2017, Beijing time, an earthquake of M7.0 occurred in Jiuzhaigou County, Aba Prefecture, Sichuan Province, with the epicenter located at 33.20°N 103.82°E. The earthquake caused 25 people dead, 525 people injured, 6 people missing and 170000 people affected. Many houses were damaged to various degrees. Up to October 15, 2017, a total of 7679 aftershocks were recorded, including 2099 earthquakes of M ≥ 1.0.
The M7.0 Jiuzhaigou earthquake occurred in the northeastern boundary belt of the Bayan Har block on the Qinghai-Tibet Plateau, where many active faults are developed, including the Tazhong Fault(the eastern segment of the East Kunlun Fault), the Minjiang fault zone, the Xueshan fault zone, the Huya fault zone, the Wenxian fault zone, the Guanggaishan-Daishan Fault, the Bailongjiang Fault, the Longriuba Fault and the Longmenshan Fault. As one of the important passages for the eastward extrusion movement of the Qinghai-Tibet Plateau(Tapponnier et al., 2001), the East Kunlun fault zone has a crucial influence on the tectonic activities of the northeastern boundary belt of Bayan Kala. Meanwhile, the Coulomb stress, fault strain and other research results show that the eastern boundary of the Bayan Har block still has a high risk of strong earthquakes in the future. So the study of the M7.0 Jiuzhaigou earthquake' seismogenic faults and stress fields is of great significance for scientific understanding of the seismogenic environment and geodynamics of the eastern boundary of Bayan Har block.
In this paper, the epicenter of the main shock and its aftershocks were relocated by the double-difference relocation method and the spatial distribution of the aftershock sequence was obtained. Then we determined the focal mechanism solutions of 24 aftershocks(M ≥ 3.0)by using the CAP algorithm with the waveform records of China Digital Seismic Network. After that, we applied the sliding fitting algorithm to invert the stress field of the earthquake area based on the previous results of the mechanism solutions. Combining with the previous research results of seismogeology in this area, we discussed the seismogenic fault structure and dynamic characteristics of the M7.0 Jiuzhaigou earthquake. Our research results indicated that:1)The epicenters of the M7.0 Jiuzhaigou earthquake sequence distribute along NW-SE in a stripe pattern with a long axis of about 35km and a short axis of about 8km, and with high inclination and dipping to the southwest, the focal depths are mainly concentrated in the range of 2~25km, gradually deepening from northwest to southeast along the fault, but the dip angle does not change remarkably on the whole fault. 2)The focal mechanism solution of the M7.0 Jiuzhaigou earthquake is:strike 151°, dip 69° and rake 12° for nodal plane Ⅰ, and 245°, 78° and -158° for nodal plane Ⅱ, the main shock type is pure strike-slip and the centroid depth of the earthquake is about 5km. Most of the focal mechanism of the aftershock sequence is strike-slip type, which is consistent with the main shock's focal mechanism solution; 3)In the earthquake source area, the principal compressive stress and the principal tensile stress are both near horizontal, and the principal compressive stress is near east-west direction, while the principal tensile stress is near north-south direction. The Jiuzhaigou earthquake is a strike-slip event that occurs under the horizontal compressive stress.

STUDY ON RELATIONSHIP BETWEEN THE ACTIVITY OF MIYALUO FAULT AND WENCHUAN EARTHQUAKE SEQUENCE

LI Zhen-yue, WAN Yong-ge, SHENG Shu-zhong

2019, 41(1):  72-83.  DOI: 10.3969/j.issn.0253-4967.2019.01.005

Asbtract ( )   HTML   PDF (2647KB) ( )  

References | Related Articles | Metrics

Under the background of thrusting stress regime, a large number of strike-slip earthquakes occurred on the Miyaluo Fault during the Wenchuan earthquake sequence process, which is in the southern part of the Longmenshan Fault. In order to find the cause of their occurrence, stress tensors in subregions near the Miyaluo Fault are estimated. The result shows that in both north and south side of the Miyaluo Fault, the direction of principal compressive stress is nearly perpendicular to the Longmenshan Fault, and its dip is nearly horizontal, and the direction of tensile stress is nearly vertical. While in the Miyaluo fault zone, the direction of principal compressive stress is SWW-NEE, and its dip is nearly horizontal, the direction of principal tensile stress is NNW-SSE, also its dip is nearly horizontal. It is consistent with sinistral shear stress state in the Miyaluo fault zone. It was referred that the behavior of Miyaluo Fault during the Wenchuan earthquake sequence process was caused by tearing effect generated from unbalanced forces of two sides of the fault. To understand the rupture mode of the aftershocks in subregions as described above, the total seismic moment tensors are estimated by adding the corresponding component separately of the seismic moment tensor of aftershocks in each region. The result shows the similar trend of total seismic moment tensor components in the north and south side of the Miyaluo Fault(indicating the consistency of rupture mode in the north and south side of the Miyaluo Fault), and most seismic moment tensor components in the south side is higher than that in the north side, especially the compression component perpendicular to Longmenshan Fault and expansion component in the vertical direction. It indicates that thrusting component in the southeast direction in the south side is greater than that in the north side, and the thrusting difference causes the sinistral tearing effect of the Miyaluo Fault. We also find that the sinistral tearing component of the Miyaluo Fault is the same order of magnitude with the thrusting difference of its two sides, which indicates that the tearing effect of Miyaluo Fault can be completely explained by thrusting difference of its two sides. According to the analysis, we put forward the dynamic model of the Miyaluo Fault, which can explain the above phenomenon.

THE THREE DIMENSIONAL DENSITY STRUCTURE OF CRUST AND UPPER MANTLE IN THE CENTRAL-SOUTHERN PART OF LONGMENSHAN

XU Zhi-ping, WANG Fu-yun, JIANG Lei, ZHAO Yan-na, YANG Li-pu, TANG Lin

2019, 41(1):  84-98.  DOI: 10.3969/j.issn.0253-4967.2019.01.006

Asbtract ( )   HTML   PDF (4923KB) ( )  

References | Related Articles | Metrics

In recent years, strong earthquakes of M_S8.0 Wenchuan and M_S7.0 Lushan occurred in the central-southern part of Longmenshan fault zone. The distance between the two earthquakes is less than 80 kilometers. So if we can obtain the inner structure of the crust and upper mantle, it will benefit us to understand the mechanism of the two earthquakes. Based on the high resolution dataset of Bouguer gravity anomaly data and the initial model constrained by three-dimensional tomography results of P-wave velocity in Sichuan-Yunnan region, with the help of the preconditioned conjugate gradient(PCG)inversion method, we established the three dimensional density structure model of the crust and upper mantle of the central-southern segment of Longmenshan, the spatial interval of which is 10 kilometers along the horizontal direction and 5 kilometers along the depth which is limited to 0~65km, respectively. This model also provides a new geophysical model for studying the crustal structure of western Sichuan plateau and Sichuan Basin. The results show obvious differences in the crustal density structure on both sides(Songpan-Ganzê block and Sichuan Basin)of Longmenshan fault zone which is a boundary fault and controls the inner crustal structure. In Sichuan Basin, the sedimentary layer is represented as low density structure which is about 10km thick. In contrast, the upper crust of Songpan-Ganzê block shows a thinner sedimentary layer and higher density structure where bedrock is exposed. Furthermore, there is a wide scale low density layer in the middle crust of the Songpan-Ganzê block. Based on this, we inferred that the medium intensity of the Songpan-Ganzê block is significantly lower than that of Sichuan Basin. As a result, the eastward movement of material of the Qinghai-Tibet plateau, blocked by the Sichuan Basin, is inevitably impacted, resulting in compressional deformation and uplift, forming the Longmenshan thrust-nappe tectonic belt at the same time. The result also presents that the crustal structure has a distinct segmental feature along the Longmenshan fault zone, which is characterized by obviously discontinuous changes in crustal density. Moreover, a lot of high- and low-density structures appear around the epicenters of Wenchuan and Lushan earthquakes. Combining with the projection of the precise locating earthquake results, it is found that Longmenshan fault zone in the upper crust shows obvious segmentation, both Wenchuan and Lushan earthquake occurred in the high density side of the density gradient zone. Wenchuan earthquake and its aftershocks are mainly distributed in the west of central Longmenshan fault zone. In the south of Maoxian-Beichuan, its aftershocks occurred in high density area and the majority of them are thrust earthquake. In the north of Maoxian-Beichuan, its aftershocks occurred in the low density area and the majority of them are strike-slip earthquake. The Lushan earthquake and its aftershocks are concentrated near the gradient zone of crustal density and tend to the side of the high density zone. The aftershocks of Lushan earthquake ended at the edge of low-density zone which is in EW direction in the north Baoxing. The leading edge of Sichuan Basin, which has high density in the lower crust, expands toward the Qinghai-Tibet Plateau with the increase of depth, and is close to the west of the Longmenshan fault zone at the top of upper mantle. Our results show that there are a lot of low density bodies in the middle and lower crust of Songpan-Ganzê Block. With the increase of the depth, the low density bodies are moving to the south and its direction changed. This phenomenon shows that the depth and surface structure of Songpan-Ganzê Block are not consistent, suggesting that the crust and upper mantle are decoupled. Although a certain scale of low-density bodies are distributed in the middle and lower crust of Songpan-Ganzê, their connectivity is poor. There are some low-density anomalies in the floor plan. It is hard to give clear evidence to prove whether the lower crust flow exists.

STUDY ON RELATIONSHIP BETWEEN SEISMIC DISTRIBUTIONOF RUSHAN SEQUENCE AND VELOCITY STRUCTURE

QU Jun-hao, WANG Chang-zai, LIU Fang-bin, ZHOU Shao-hui, ZHENG Jian-chang, LI Xin-feng, ZHANG Qin

2019, 41(1):  99-118.  DOI: 10.3969/j.issn.0253-4967.2019.01.007

Asbtract ( )   HTML   PDF (5447KB) ( )  

References | Related Articles | Metrics

Since the earthquake of M_L3.8 occurring on October 1, 2013 in Ruishan, Weihai City, Shandong Province, the sequence has lasted for about 4 years(Aug. 31, 2017). Seismicity is enhanced or weakened and fluctuated continuously. More than 13250 aftershocks have been recorded in Shandong Seismic Network. During this period, the significant earthquake events were magnitude 4.2(M_L4.7)on January 7, 4.0(M_L4.5)on April 4, M3.6(M_L 4.1)on September 16 in 2014 and M4.6(M_L5.0)on May 22, 2015. The earthquake of M_L5.0 was the largest one in the Rushan sequence so far. In order to strengthen the monitoring of aftershocks, 18 temporary stations were set up near the epicenter at the end of April, 2014(official recording began on May 7)by Shandong Earthquake Agency, which constitutes an intensified network in Rushan that surrounds the four quadrants of the small earthquake concentration area together with 12 fixed stations nearby, and provides an effective data foundation for the refinement of Rushan earthquake sequence.
The velocity structure offers important information related to earthquake location and the focal medium, providing an important basis for understanding the background and mechanism of the earthquake. In this paper, double-difference tomography method is used to relocate the seismic events recorded by more than six stations of Rushan array from May 7, 2014 to December 31, 2016, and the inversion on the P-wave velocity structure of the focal area is conducted. The Hyposat positioning method is used to relocate the absolute position. Only the stations with the first wave arrival time less than 0.1 second are involved in the location. A total of 14165 seismic records are obtained, which is much larger than that recorded by Shandong Seismic Network during the same period with 7708 earthquakes and 2048 localizable ones. A total of 1410 earthquakes with M_L ≥ 1.0 were selected to participate in the inversion. Precise relocation of 1376 earthquakes is obtained by using double-difference tomography, in which, there are 14318 absolute traveltime P waves and 63162 relative travel time P waves. The epicenters are located in distribution along NWW-SEE toward SEE and tend to WS, forming a seismic belt with the length about 3km and width about 1km. The focal depths are mainly concentrated between 4km and 9km, occurring mainly at the edge of the high velocity body, and gradually dispersing with time. It has obvious temporal and spatial cluster characteristics. Compared with the precise relocation of Shandong network, the accuracy of the positioning of Rushan array is higher. The main reason is that the epicenter of Rushan earthquake swarm is near the seaside, and the fixed stations of Shandong Seismic Network are located on the one side of the epicenter. The nearest three stations(RSH, HAY, WED)from the epicenter are Rushan station with epicentral distance about 13km, the Haiyang station with epicentral distance about 33km, and Wendeng station with epicentral distance about 42km. The epicentral distance of the rest stations are more than 75km. In addition, the magnitude of most earthquakes in Rushan sequence is small. The accuracy of phase identification is relatively limited due to the slightly larger epicentral distance of the station HAY and station WED in Shandong Seismic Network. Furthermore, the one-dimensional velocity model used in network location is simple with only the depth and velocity of Moho surface and Conrad surface. The epicentral distances of the 18 temporary stations in Rushan are less than 10km, and the initial phase is clear. The island station set up on the southeast side and the Haiyangsuo station on the southwest side form a comprehensive package for the epicenter. Compared with the double-difference algorithm method, the double-difference tomography method used in this paper is more accurate for the velocity structure, thus can obtain the optimal relocation result and velocity structure.
the velocity structure shows that there are three distinct regions with different velocities in the vicinity of the focal area. The earthquakes mainly occur in the intersection of the three regions and on the side of the high velocity body. With the increase of depth, P wave velocity increases gradually and there are two distinct velocity changes. The aftershock activities basically occur near the dividing line to the high velocity side. The south side is low velocity abnormal body and the north side is high velocity abnormal body. High velocity body becomes shallower from south to north, which coincides with the tectonic conditions of Rushan. Considering the spatial relationships between the epicenter distribution and the high-low velocity body and different lithology of geological structure, and other factors, it is inferred that the location of the epicenter should be the boundary of two different rock bodies, and there may be a hidden fault in the transition zone between high velocity abnormal body and low velocity abnormal body. The interface position of the high-low velocity body, the concentrating area of the aftershocks, is often the stress concentration zone, the medium is relatively weak, and the intensity is low. There is almost no earthquake in the high velocity abnormal body, and the energy accumulated in the high velocity body is released at the peripheral positions. It can be seen that the existence of the high-low velocity body has a certain control effect on the distribution of the aftershocks.

VELOCITY CHARACTRISTICS OF SHANXI AND ADJACENT AREA AND ITS TECTONIC SIGNIFICANCE

WANG Xia, SONG Mei-qin, ZHENG Yong, AI San-xi

2019, 41(1):  119-136.  DOI: 10.3969/j.issn.0253-4967.2019.01.008

Asbtract ( )   HTML   PDF (10063KB) ( )  

References | Related Articles | Metrics

In order to acquire a better velocity structure of the crustal and uppermost mantle beneath Shanxi area, we obtain the group and phase velocities of Rayleigh wave of the periods 8s to 50s in Shanxi and adjacent area using ambient seismic noise recorded at 216 broad-band stations. All available vertical-component time series for 2014 have been cross-correlated to yield estimates of empirical Rayleigh wave Green's function. Group and phase velocity dispersion curves for Rayleigh wave are measured for each interstation path by applying frequency-time analysis. It describes finer velocity structure of the crust and upper mantle in Shanxi, which reflects the geological structure characteristics at different depths. The resolution is within 50km and the resolution of part periods can reach 40km.The Rayleigh wave group and phase speed maps at short periods(8~18s and 10~22s)show clear correlations with shallow geological structures. Mountain areas on both sides of Shanxi depression zone show apparent high-velocity anomaly, except for low-velocity anomaly in the Taiyuan Basin, Linfen-Yuncheng Basin and Weihe Basin. Especially, the areas of Youyu County-Pianguan County-Kelan County-Shuozhou City and Jingle County-Lishi District of Lüliang City in Lüliang Mountains, and Yu County-Fuping County-Yi County and Yangcheng County-Licheng County in Taihang Mountains, present higher velocity anomaly. In addition, the velocity is lowest in the Weihe Basin, and the amplitude of low velocity decreases gradually from the south to the north of the basins in Shanxi, which probably is related to the process of gradual stretching and development of the Shanxi rift zone from the southwest to the northeast. The obvious velocity difference across the latitude of 38°N exists at 18~30s period of phase and 24~35s period of group velocity maps, which is probably related to the deep and shallow Moho depth variation in the south and north of Shanxi and the suture zone of ancient blocks including "hard" southern block and "soft" northern block. At the same time, the research result of receiver function reveals that partial melting of the lower crust occurs in the northern Taihang Mountains, while the southern section remains stable(Poisson's ratio is above 0.3 in the northern Taihang Mountains and 0.25~0.26 in the southern section). The phase velocity map at 30~50s period clearly shows NW velocity gradient belt, and the low velocity anomaly in the northeast side may be related to Cenozoic volcanism. Meanwhile, the eastern border of Ordos block is the western faults of central basins in Shanxi depression zone. However, some research results indicate that the above border is Lishi Fault in the surface, inferring that the Ordos block shows a shape of wide in the upper and narrow in the lower part from the surface to deep. The Datong volcanic area at 18~45s period of phase and 24~35s period of group velocity maps shows low velocity of trumpet shape from shallow to deep, related to the upwelling of hot material from lower mantle in the Cenozoic causing a large area of intense magmatic activity. It indicates the more specific upwelling channel of Datong volcanoes simultaneously.

GEODETIC AND TELESEISMIC CONSTRAINTS ON SLIP DISTRIBUTION OF 2015 M_W6.4 PISHAN EARTHQUAKE

CHEN Yun-guo, HE Ping, DING Kai-hua, LI Shui-ping, WANG Qi

2019, 41(1):  137-149.  DOI: 10.3969/j.issn.0253-4967.2019.01.009

Asbtract ( )   HTML   PDF (5974KB) ( )  

References | Related Articles | Metrics

On July 3^rd, 2015, a M_W6.4 earthquake occurred on Pishan County, Xinjiang, located in the front of western Kunlun thrust belt, which is the largest earthquake(M_W6.0~7.0)in the past 40 years in this region. In this study, we collected both the near-filed geodetic coseismic deformation observations including 4 GPS sites and one high-resolution ALOS-2 InSAR imagery, and far-field teleseismic P waveforms from 25 stations provided by IRIS/USGS, to invert the fault parameters(strike and dip)and coseismic rupture model of 2015 M_W6.4 Pishan earthquake. Using the finite fault theory, a non-linear simulated annealing algorithm was employed to resolve our joint inversion problem. The strike (120°~130°) and dip angle(35°~40°)of optimal models are different from that of some previous studies, and the dip change is strongly constrained by combined data than that of strike. In fixing the geometric parameters of optimal fault model, we also considered data weight(5)(geodetic data/teleseismic P waveforms)and constrained weight from moment and smooth factor(2.5). Clearly, our results indicate that the slip distribution mainly concentrates in the depth range from 9 to 16km and a length range of 20km along the strike direction, which is similar to the spatial distribution of the relocated aftershocks. The maximum slip is~95cm. The seismic moment release is 5.45×10¹⁸N·m, corresponding to M_W6.42. Compared with the single data set, geodetic data or teleseismic waveform, our joint inversion model could simultaneously constrain the seismic moment and slip distribution well, thus avoiding effectively a lower-resolution rupture distribution determined by teleseismic-only inversion and a bias released moment estimated by the geodetic-only inversion. Importantly, we should consider both the near-field geodetic data and far-field teleseismic data in retrieving the rupture model for accurately describing the seismogenic structure of active fault in western Kunlun region.

SOURCE PARAMETERS OF THE CANGWU M_S5.4 EARTHQUAKE, 31 JULY, 2016

ZHOU Yi, YAN Chun-heng, XIANG Wei, ZHOU Bin, WEN Xiang

2019, 41(1):  150-161.  DOI: 10.3969/j.issn.0253-4967.2019.01.010

Asbtract ( )   HTML   PDF (3548KB) ( )  

References | Related Articles | Metrics

On July 31^st, 2016, an earthquake of M_S5.4 occurred in Cangwu County, Guangxi Zhuang Autonomous Region, which is the first M_S ≥ 5.0 earthquake in coastal areas of southern China in the past 17a. The moderate earthquake activities have come into a comparatively quiet period in coastal areas of southern China for decades, so the study about the Cangwu M_S5.4 earthquake is very important. However, differernt research institutions and scholars have got different results for the focal depth of the Cangwu M_S5.4 earthquake. For this reason, we further measured the focal depth by using CAP method and sPL phase method.
sPL phase was first put forward by Chong in 2010. It is often observed between P and S wave of continental earthquakes with epicentral distance of about 30km to 50km. The energy of sPL phase is mainly concentrated on the radial component. Arrival time difference between sPL phase and direct P wave is insensitive to epicentral distancs, but increases almost linearly with the increase of focal depth. Based on these characteristics and advantages, sPL phase method is chosen to measure the focal depth of Cangwu M_S5.4 earthquake in the paper.
First of all, we selected the broadband waveform data through seismic stations distributed mainly in Guangxi and adjacent provinces from Data Management Centre of China National Seismic Network and Guangxi Earthquake Networks Center. And an appropriate velocity model of Cangwu area was constructed by the teleseismic receiver function method. Then, the focal mechanism and focal depth of Cangwu M_S5.4 earthquake were determined by using the CAP(Cut and Paste)method. Next, we compared the synthetic waveforms simulated by F-K forward method of different focal depth models with the actual observed waveforms. According to the difference of arrival times between sPL and Pg phases, we finally obtained the focal depth of Cangwu earthquake. The results show that the focal depth is 11km measured by CAP method and 9km by sPL phase method. Based on the focal mechanism solution, isoseismal shapes, aftershocks distributions and investigation on spot, we conclude that the Cangwu M_S5.4 earthquake is a left-lateral strike-slip earthquake which occurred in the upper crust. Our preliminary analysis considers that the seismogenic structure of Cangwu earthquake is a north-northwest branch fault, and the control fault of this earthquake is the Hejie-Xiaying Fault.

APPLICATION OF LOGISTIC REGRESSION MODEL FOR HAZARD ASSESSMENT OF EARTHQUAKE-TRIGGERED LANDSLIDES: A CASE STUDY OF 2017 JIUZHAIGOU(CHINA)M_S7.0 EVENT

MA Si-yuan, XU Chong, TIAN Ying-ying, XU Xi-wei

2019, 41(1):  162-177.  DOI: 10.3969/j.issn.0253-4967.2019.01.011

Asbtract ( )   HTML   PDF (12909KB) ( )  

References | Related Articles | Metrics

The M_S7.0 Jiuzhaigou earthquake in Sichuan Province of 8 August 2017 triggered a large number of landslides. A comprehensive and objective panorama of these landslides is of great significance for understanding the mechanism, intensity, spatial pattern and law of these coseismic landslides, recovery and reconstruction of earthquake affected area, as well as prevention and mitigation of landslide hazard. The main aim of this paper is to present the use of remote sensing images, GIS technology and Logistic Regression(LR)model for earthquake triggered landslide hazard mapping related to the 2017 Jiuzhaigou earthquake. On the basis of a scene post-earthquake Geoeye-1 satellite image(0.5m resolution), we delineated 4834 co-seismic landslides with an area of 9.63km². The ten factors were selected as the influencing factors for earthquake triggered landslide hazard mapping of Jiuzhaigou earthquake, including elevation, slope angle, aspect, horizontal distance to fault, vertical distance to fault, distance to epicenter, distance to roads, distance to rivers, TPI index, and lithology. Both landsliding and non-landsliding samples were needed for LR model. Centroids of the 4834 initial landslide polygons were extracted for landslide samples and the 4832 non-landslide points were randomly selected from the landslide-free area. All samples(4834 landslide sites and 4832 non-landslide sites)were randomly divided into the training set(6767 samples)and validation set(2899 samples). The logistic regression model was used to carry out the landslide hazard assessment of the Jiuzhaigou earthquake and the results show that the landslide hazard assessment map based on LR model is very consistent with the actual landslide distribution. The areas of Wuhuahai-Xiamo, Huohuahai and Inter Continental Hotel of Jiuzhai-Ruyiba are high hazard areas. In order to quantitatively evaluate the prediction results, the trained model calculated with the training set was evaluated by training set and validation set as the input of the model to get the output results of the two sets. The ROC curve was used to evaluate the accuracy of the model. The ROC curve for LR model was drawn and the AUC values were calculated. The evaluation result shows good prediction accuracy. The AUC values for the training and validation data set are 0.91 and 0.89, respectively. On the whole, more than 78.5% of the landslides in the study area are concentrated in the high and extremely high hazard zones. Landslide point density and landslide area density increase very rapidly as the level of hazard increases. This paper provides a scientific reference for earthquake landslides, disaster prevention and mitigation in the earthquake area.

TEXTUAL RESEARCH OF 1568 M7 GAOLING EARTHQUAKE IN SHAANXI AND ANALYSIS OF ITS SEISMOGENIC STRUCTURE

MA Ji, FENG Xi-jie, LI Gao-yang, LI Xiao-ni, SHI Ya-qin

2019, 41(1):  178-188.  DOI: 10.3969/j.issn.0253-4967.2019.01.012

Asbtract ( )   HTML   PDF (6412KB) ( )  

References | Related Articles | Metrics

Study of historical earthquake is one of the important methods to understand the seismic activities and analyze the seismogenic faults. On the May 25^th, 1568 AD, a destructive earthquake occurred to the northeast of the present-day city of Xi'an, Shaanxi Province. Because this earthquake happened shortly after the 1556 M8 earthquake and was regarded as an aftershock, it has received little attention in previous studies. Previous earthquake catalogue agreed in assigning a magnitude 6 3/4 to this earthquake but had different epicentral locations and seismic intensity, and the seismogenic structure remains ambiguous.
Based on textual research of historical earthquake and field investigation, the Jingyang County, Gaoling County, and Xianning County, were the worst hit area by the earthquake, and the areas, including Yongle Town, Gaozhuang Town at southeastern Jingyang County to Gaoling County and its southeastern present-day Jijia and Zhangbu, should be the mesoseismal area of this earthquake. The epicenter intensity of this earthquake is Ⅸ⁺(9~10 degrees), and the magnitude is estimated to be 7. The isoseismal lines were drawn to exhibit the various intensities of the areas damaged during the event, with its major axis directed NWW. Intensities reached Ⅸ⁺ in the zone extending west-northwest parallel to the Weinan-Jingyang Fault. This fault, characterized by a normal fault that developed during the Cenozoic extensional history of the Weihe Basin, dipping to the north at an angle of 60°~80°, is one part of the southern boundary faults in Weihe graben. There are geomorphological and geological evidences of recent activity of the fault during (180±30)a BP to (1 600±30)a BP. At T₁-T₂ fluvial terraces on the north bank of Weihe River, the scarps were faulted during Ming Dynasty, and sandy soil liquefaction, dense structural tensional fissures and faulted strata are noted in stratigraphic profiles and trenches. Thus, we suggest that this fault can reliably be regarded as being active during Holocene, and re-name the earthquake as the Shaanxi Gaoling earthquake.

THE VOLCANIC ACTIVITIES AND HAZARD PREDICTION OF WUDALIANCHI VOLCANIC BELT

ZOU Ying, ZHAO Yong-wei, FAN Qi-cheng

2019, 41(1):  189-207.  DOI: 10.3969/j.issn.0253-4967.2019.01.013

Asbtract ( )   HTML   PDF (9600KB) ( )  

References | Related Articles | Metrics

More than 40 late Cenozoic monogenetic volcanoes formed a volcanic belt striking NNW from Keluo, through Wudalianchi to Erkeshan in NE China. These volcanoes belong to a unified volcano system, namely Wudalianchi volcanic belt(WVB for short). Based on the volcanic evolution history and the nature of monogenetic volcanic system, we estimate that the volcanic system of WVB is still active and has the potential to erupt again. Hence, this paper studied the temporal-spatial distribution and volcanic eruption types to evaluate the possible eruption hazard types and areas of influence in the future.
Volcanic field characteristics and K-Ar radiometric data suggest two episodes of volcanism in the WVB, the Pliocene to early Pleistocene volcanism(4.59~1.00MaBP)and the middle Pleistocene to Holocene volcanism(0.79Ma to now). The early episode volcanoes are distributed only in the north of WVB(mainly in Keluo volcanic field), featured by effusive eruption, and mainly formed monogenetic shield, whose base diameter is large and slope is gentle. However, the late episode eruptions occurred over the entire WVB. The explosive eruption in this stage formed numerous relatively intact scoria cones of explosive origin. Meanwhile the effusive eruption formed widely distributed lava flows.
Both effusive eruption and explosive eruption are common in WVB. The effusive eruption formed monogenetic shields and lava flows. The resulting pahoehoe lava, aa lava and block lava appeared in WVB. There are three end-member types of explosive eruption driven by magmatic volatile. Violent Strombolian eruption has the highest degree of fragmentation and mass flux, characterized by eruption column. Strombolian eruption has the high degree of fragmentation, but low mass flux, featured by pulse eruption. Hawaiian eruption has low degree of fragmentation, but high in mass flux, generating large scoria cones. In addition, this paper for the first time found phreatomagmatic eruption in WVB, which formed tuff cone. Transitional eruptions are also common in WVB, which have certain characteristics among the end-member eruption types. Besides, certain volcanoes displayed multiple explosive eruption types during the whole eruption span.
According to the volcanic temporal-spatial distribution and eruption characteristics in WVB, the potential volcanic hazards in future are constrained. It appears that the violent Strombolian and Strombolian eruption will not have significant impact on aviation safety in the vertical direction. In the radial direction, the ejected volcanic bomb can reach as far as 1km from the vents and the fallout tephra may disperse downwind over a distance ranging from 1~10km. The major hazard of Hawaiian eruption and effusive eruption comes from lava flow, and its migration distance may reach 3.0~13.5km for pahoehoe lava and 2.9~14.9km for aa lava. The base surge in phreatomagmatic eruption can reach a velocity of 200~400m/s, and the migration distance is around 10km. This is a big threat that people should pay more attention to and take precautions in advance. Besides, it is necessary to strengthen the real-time observation of the volcanoes in the WVB, especially those formed in the late episode as well as near the active fault.

A COMPARATIVE STUDY ON THE CHARACTERISTICS OF TWO STAGES OF FALLOUT PUMICES DEPOSITS FROM THE MILLENNIUM ERUPTION OF TIANCHI VOLCANO IN CHANGBAISHAN AREA

WANG Yu-qin, YU Hong-mei, XU Jian-dong, CHEN Zheng-quan, ZHAO Bo

2019, 41(1):  208-224.  DOI: 10.3969/j.issn.0253-4967.2019.01.014

Asbtract ( )   HTML   PDF (8258KB) ( )  

References | Related Articles | Metrics

Tianchi volcano in Changbaishan area is located at the border between China and Democratic People's Republic of Korea, and is one of the most dangerous volcanoes in China. It has experienced several explosive eruptions in late Pleistocene and Holocene, i.e. 50000aBP eruption, 946 AD eruption, 1668 AD eruption, 1702 AD eruption, 1903 AD eruption. Especially, the 946 AD eruption(also known as "Millennium eruption")of this volcano is considered to be one of the largest volcanic eruptions in the world in the past 2000a. The eruption history and strata sequence of Tianchi volcano have long been the focus of attention. The stratigraphic unit division of fallout deposits in the past millennium is controversial, especially for the heterogeneous trachytic pumices(erupted from the Yuanchi stage)above the off-white pumices(erupted from the Chifeng stage).
In this paper, through the detailed field exploration and strata comparation, it was found that there was no depositional interval between the two stage eruptions, or the interval was not long, and thus, it is believed that two stages of fallout pumice should be classified into the Millennium eruption. The off-white fallout pumices in Chifeng stage are relatively homogeneous, with angular shape, normal grading and good sorting. The median size(Md_Φ)and the sorting coefficient(σ_Φ)of Chifeng pumice are in the range of -4.25~-1.3 and 0.93~1.53, respectively. The eruption of Yuanchi stage is in pulsing pattern, and the strata show interbedding of rich khaki pumice layer and rich black pumice layer. The pumices with angular shape show inconspicuous grain grading and good sorting. The median size(Md_Φ)and the sorting coefficient(σ_Φ)of Yuanchi pumice are in the range of -2.55~-0.6 and 1~1.68, respectively. Both the granularities of the pumice particles from two stages are normally distributed and fall into the air-fall field in the median diameter versus sorting diagram. The pumices from 50000aBP and pyroclastic flow of Millennium eruption were also shown in the diagram.
Phenocrysts in pumices are mainly feldspar and pyroxene, but the phenocrysts with obvious resorbed characteristic in Yuanchi black pumice are bigger, and the phenocryst contents are a little higher than those in others. Feldspar content in off-white pumice in Chifeng stage was 0.24%~1.77%, that in khaki pumice in Yuanchi stage was 0.2%~7.5%, and that in black pumice in Yuanchi stage was 3.02%~8.0%. The phenocrysts in Chifeng pumice are broken, which represents more violent explosion. The vesicles inside the pumice also reflect the intensity of the eruption. The Chifeng pumices have large, continuous vesicles and thin vesicle walls. The Yuanchi khaki pumices have continuous vesicles but thicker vesicle wall than the Chifeng pumices. The vesicularity is the lowest and the vesicle walls are the thickest in the black pumices in Yuanchi stage, indicating the eruption strength become weaker from Chifeng stage to Yuanchi stage.
The Chifeng pumices with SiO₂ content of 69.12~72.71wt%, K₂O content of 4.33~4.52wt%, Na₂O content of 5.26~5.39wt%, Al₂O₃ content of 10.32~11.99wt%, CaO content of 0.29~0.95wt%, MgO content of 0.11~0.51wt%, TiO₂ content of 0.23~0.43wt% are comendite in composition. The pumices from 50000aBP eruption are comendite in composition, and their SiO₂ content(65.56~68.28wt%)is slightly lower than Chifeng pumices. The Yuanchi khaki pumices with SiO₂ content of 62.14~63.29wt%, K₂O content of 5.35~5.7wt%, Na₂O content of 5.35~5.62wt%, Al₂O₃ content of 15.00~15.59wt%, CaO content of 1.06~1.61wt%, MgO content of 0.25~0.57wt%, TiO₂ content of 0.4~0.64wt% belong to trachyte in composition, and are close to the composition of the black pumices on the Tianwen Peak. The Yuanchi black pumices are also trachyte in composition, but have obviously lower SiO₂(59.51~60.59wt%), K₂O(4.39~4.84wt%), and Na₂O(4.94~5.08wt%)content, and higher Al₂O₃(15.81~16.42wt%), CaO(2.78~3.66wt%), MgO(1.43~1.9wt%), TiO₂(1.04~1.4wt%)content than the khaki pumices.
The above results show that the eruptive intensity of the Yuanchi stage is weaker than that of the Chifeng stage and the several magmatic compositions of pumices from the Millennium eruption reveal a complex magma system under the Tianchi volcano. The magma layers with different compositions may exist in the magma chamber contemporaneously. At Chifeng stage, only the upper comendite magma erupted, but the magma below erupted in the pulsing pattern at the Yuanchi stage.

TEPHRA RECORD FROM QUANYANG PEAT OF THE CHANGBAISHAN MILLENNIUM ERUPTION

HAN Ling-fei, LIU Jia-qi, YOU Hai-tao, ZHU Ze-yang, SUN Chun-qing

2019, 41(1):  225-236.  DOI: 10.3969/j.issn.0253-4967.2019.01.015

Asbtract ( )   HTML   PDF (2033KB) ( )  

References | Related Articles | Metrics

Tephra, usually produced by explosive eruptions, is deposited rapidly, hence, it can serve as a distinctive and widespread synchronous marker horizon correlating terrestrial, marine and ice core records. The tephra from Changbaishan Millennium eruption, a widely distributed tephra, is an important marker bed across the Japan Sea, Japan Islands and even in the Greenland ice cores 9000km away from volcanic vent. In this study, a discrete tephra was identified in the Quanyang peat~45km northeast to the Changbaishan volcano. Radiocarbon ¹⁴ C dating on the plant remains constrains an age of 886-1013calAD(95.4%)to the tephra layer, which can correspond to the Millennium eruption of Changbaishan in time. In addition, there was no similar volcanic eruption in the surrounding areas except Changbaishan at the same time. This tephra shows rhyolitic glass shards major element compositions similar to those rhyolitic tephra from Millennium eruption. This study illustrates that tephra from Millennium eruption has been transported to Quanyang peat~45km northwest to the Changbaishan volcano. Additionally, the diameter of the pumice lapilli is up to 0.3cm, implying that the tephra must be transported more distal away from Quanyang peat and formed a widely distributed isochronic layer. Glass geochemistry of the Quanyang tephra, different from the distal tephra recorded at Sihailongwan, Japan, and Greenland ice, shows a close affinity to the pyroclastic flow deposits of the Millennium eruption while not from fall deposits. This may indicate that distribution of the Millennium eruption of Changbaishanin in different directions may be controlled by different stages of eruption. This layer with well-defined annual results can be used to optimize the chronological framework of the corresponding sedimentary environment, thus facilitating more accurate discussion of corresponding environmental changes, which can achieve the contrast of the ancient climate records in the whole Northeast China-Japan and arctic regions.

PRECISE AND ACCURATE IN SITU U-Pb DATING OF ZIRCON BY LA-ICP-MS

LI Chao-peng, ZHENG De-wen, WANG Ying, PANG Jian-zhang, XIAO Lin, LI You-juan

2019, 41(1):  237-249.  DOI: 10.3969/j.issn.0253-4967.2019.01.016

Asbtract ( )   HTML   PDF (3660KB) ( )  

References | Related Articles | Metrics

LA-ICP-MS(laser ablation-inductively coupled-mass spectrometry)has been recently used for rapid, accurate and precise U-Pb geochronology on zircon grains. In this paper, we adopted an Agilent 7900 quadrupole ICP-MS coupled with a Resolution M50-LR 193nm excimer laser system to establish integrated measurement procedures. Before analysis, the system is tuned to achieve sensitivities better than 30000 cps/s for ²³⁸ U with a 40μm spot size, at~3.5J/cm². Detailed parameters for laser system and ICP-MS are presented here. Then, we analyzed five reference zircons(91500, GJ-1, Plesovice, FCT, Penglai)with a wide range in age from~1064 to~4.4Ma. Two standard zircons, 91500 and GJ-1, are employed as external reference standards. Generally, second zircon standard is analyzed in an effort to ensure accuracy and evaluate reproducibility. A typical analysis sequence includes one international glass standard(NIST610), two external reference standards, five grains of unknown zircon with every eight ablations. Laser induced time-dependent elemental fractionation is corrected using the intercept method, whereas instrument drift, mass bias and elemental fractional caused by ionization differences are corrected by external reference standard 91500 or GJ-1. Compared with 91500 and GJ-1, common Pb content of Plesovice, FCT, Penglai can't be ignored. Thus, we did common Pb correction for the above three standard zircons. The performance of the established procedure was assessed by analyzing zircon range in age from~1 064 to~4Ma. The results show that the ages of these five references are consistent with the ages of published studies with accuracy for three international references(91500, GJ-1, Plesovice)better than 3% and two young secondary references(FCT, Penglai)lower than 7% at the 2 sigma level, which indicates that our analytical procedure is reliable. For individual laser analysis, the uncertainties are mainly from three sources:Measurement error of isotope ratio, error of correction factors for instrument drift and element fractionation, and error of recommended age of external references. Compared to three international references, there are three extra uncertainties for young reference zircons, including:1)little radioactive isotopes closing to blank level increase the measurement error of isotope ratio; 2) effect of common lead becomes more significant;3) the nonhomogeneous samples couldn't match references well. Therefore, accuracy and precision of measurement depend on absolute age, content of common lead and matching degree between references and samples. In summary, the accuracy and precision obtained using the technique presented in this study are similar to those of other LA-ICP-MS laboratories.