Table of Content

20 June 2020, Volume 42 Issue 3

Previous Issue    Next Issue

FOCAL MECHANISM SOLUTIONS AND STRESS FIELD OF THE 2019 CHANGNING, SICHUAN MAINSHOCK AND ITS MODERATE-STRONG AFTERSHOCKS(M_S≥4.0)

LIANG Shan-shan, XU Zhi-guo, SHENG Shu-zhong, ZHANG Guang-wei, ZHAO-Bo, ZOU Li-ye

2020, 42(3):  547-561.  DOI: 10.3969/j.issn.0253-4967.2020.03.001

Asbtract ( )   HTML   PDF (3702KB) ( )  

References | Related Articles | Metrics

A M_S6.0 earthquake with shallow focal depth of 16km struck Changning County, Yibin City, Sichuan Province at 22:55: 43(Beijing Time)on 17 June 2019. Although the magnitude of the earthquake is moderate, it caused heavy casualties and property losses to Changning County and its surrounding areas. In the following week, a series of aftershocks with M_S≥4.0 occurred in the epicentral area successively. In order to better understand and analyze the seismotectonic structure and generation mechanism of these earthquakes, in this paper, absolute earthquake location by HYPOINVERSE 2000 method is conducted to relocate the main shock of M_S6.0 in Changning using the seismic phase observation data provided by Sichuan Earthquake Administration, and focal mechanism solutions for Changning M_S6.0 main shock and M_S≥4.0 aftershocks are inferred using the gCAP method with the local and regional broadband station waveforms recorded by the regional seismic networks of Sichuan Province, Yunnan Province, Chongqing Municipality, and Guizhou Province. The absolute relocation results show that the epicenter of the main shock is located at 28.35°N, 104.88°E, and it occurred at an unusual shallow depth about only 6.98km, which could be one of the most significant reasons for the heavier damage in the Changning and adjoining areas. The focal plane solution of the Changning M_S6.0 earthquake indicates that the main shock occurred at a thrust fault with a left-lateral strike-slip component. The full moment tensor solution provided by gCAP shows that it contains a certain percentage of non-double couple components. After the occurrence of the main shock, a series of medium and strong aftershocks with M_S≥4.0 occurred continuously along the northwestern direction, the fault plane solutions for those aftershocks show mostly strike-slip and thrust fault-type. It is found that the mode of focal mechanism has an obvious characteristic of segmentation in space, which reflects the complexity of the dislocation process of the seismogenic fault. It also shows that the Changning earthquake sequences occurred in the shallow part of the upper crust. Combining with the results from the seismic sounding profile in Changning anticline, which is the main structure in the focal area, this study finds that the existence of several steep secondary faults in the core of Changning anticline is an important reason for the diversity of focal mechanism of aftershock sequences. The characteristics of regional stress field is estimated using the STRESSINVERSE method by the information of focal mechanism solutions from our study, and the results show that the Changning area is subject to a NEE oriented maximum principal stress field with a very shallow dipping and near-vertical minimum principal stress, which is not associated with the results derived from other stress indicators. Compared with the direction of the maximum principal compressive stress axis in the whole region, the direction of the stress field in the focal area rotates from the NWW direction to the NEE direction. The Changning M_S6.0 earthquake locates in the area with complex geological structure, where there are a large number of small staggered fault zones with unstable geological structure. Combining with the direction of aftershocks distribution in Changning area, we infer that the Changning M_S6.0 earthquake is generated by rupturing of the pre-existing fault in the Changning anticline under the action of the overall large stress field, and the seismogenic fault is a high dip-angle thrust fault with left-lateral strike-slip component, trending NW.

DETERMINATION OF FAULT PLANE PARAMETERS IN THE LONGTAN RESERVOIR BY USING PRECISELY LOCATED SMALL EARTHQUAKE DATA AND REGIONAL STRESS FIELD

YAN Chun-heng, ZHOU Bin, LI Sha, XIANG Wei, GUO Pei-lan

2020, 42(3):  562-580.  DOI: 10.3969/j.issn.0253-4967.2020.03.002

Asbtract ( )   HTML   PDF (3927KB) ( )  

References | Related Articles | Metrics

The Longtan reservoir is located in Tian'e County, Guangxi Zhuang Autonomous Region, southwestern China on the upper reaches of Hongshui River, the main stream of the Pearl River. The dam of the reservoir is 200m high, and the maximum water depth can be up to 194m as the water level reaches 400m. The reservoir storage capacity is 27.3 billion cubic meters, so it is a typical high-dam reservoir with large storage capacity. Terrain of the reservoir is high in the west and low in the east. The reservoir is located at the confluence of the Hongshui River, Buliu River, Nanpan River, Beipan River, Mengjiang River and Caodu River. The construction of Longtan hydropower station officially started in July 2001, and the reservoir impoundment was on September 30, 2006. The power station is equipped with 9 sets of 700 000kW water turbine generator units, with a total installed capacity of 6.3 million kW and an average annual generating capacity of 18.7 billion kW·h. So its storage and hydropower capacity rank third only to the world-famous Three Gorges hydropower project and the ultra-large hydropower project in Xiluodu of Jinsha River in China. Seismicity enhanced rapidly in the reservoir area after the impoundment. More the 5 000 earthquakes have been recorded so far, with the maximum magnitude of M_L4.8, which occurred on September 18, 2010. The earthquakes are mainly concentrated in the deep water area where fault zones run through. Assuming the seismogenic fault can be simulated by a plane and most small earthquakes occur nearby the fault plane, the information of seismogenic fault can be obtained by the hypocenter location parameters of small earthquakes.

THE FINE CRUSTAL STRUCTURE OF THE SOUTHERN MARGIN OF NORTH CHINA BLOCK REVEALED BY DEEP SEISMIC REFLECTION PROFILE

FENG Shao-ying, LIU Bao-jin, LI Qian, YUAN Hong-ke, ZHU Guo-jun, TIAN Yi-ming, WANG Hong-wei, HOU Li-hua, DENG Xiao-juan, TAN Ya-li

2020, 42(3):  581-594.  DOI: 10.3969/j.issn.0253-4967.2020.03.003

Asbtract ( )   HTML   PDF (4756KB) ( )  

References | Related Articles | Metrics

The study area is located at the junction of the northern margin of the Qinling orogenic belt and the southern margin of the North China Block. In order to study the fine crustal structure and the deep-shallow structural features of faults in this area, we conducted deep seismic reflection profiling with the seismic profile of 100km long, directing NE-SW in Zhumadian City, Henan Province, and got clear lithospheric structure images along the profile. As regards the data acquisition, we applied the geometry of 25m group interval, 1000 recording channels and more than 60 folds. Seismic wave exploding applies the 30kg shots of dynamite source with the borehole depth of 25m. The shot interval is 200m. In data processing, we focused on improving the signal-to-noise ratio. Data processing methods mainly include first break removal, tomographic static correction, abnormal amplitude elimination, amplitude compensation, pre-stack denoising, surface consistent deconvolution, velocity analysis, several iterations of the residual static correction, dip moveout, post-stack time migration and post-stack denoising, etc. The profile with high signal-to-noise ratio was obtained. The reflection wave group characteristics is obvious in the crust, which reflects abundant information about geological structure. Along the profile, the crust is characterized by double-layer reflection structure, and the Moho surface is composed of a series of laminated arc-shaped strong reflections. The thickness of the upper crust is about 14.8~20.7km, and the total thickness of the crust is about 32.0~35.1km. The upper crust is dominated by the inclined, densely stratified or arc-shaped reflections. The lower crust is dominated by arc-shaped and inclined reflection, and there is a reflective transparent zone under the Moho surface. The reflection sequences with different directions and shapes in the upper crust constitute the nappe structure in southwest segment and the structural model of two concaves with one uplift in NE segment, which correspond to the north Qinling nappe, Zhumadian-Huaibin depression, Pingyu-Xiping uplift and a secondary depression, respectively. There are abundant arc-shaped reflection sequences in the lower crust, which may represent multi-stage magmatic activities. The deep seismic reflection profile shows that faults in the upper crust are well developed. According to the characteristics of reflected wave field in the profile, four groups of fault structure which contain ten faults with different scales are interpreted. Among them, faults F_P1, F_P2 and F_P3 constitute the thrust fault system in the northern margin of Qinling Mountains, and F_P5 and F_P7 are boundary faults of Zhumadian-Huaibin depression. These faults are all developed within the upper crust. In addition, the Fault F_PM is a large fault that cuts through the lower crust and Moho surface. The deep seismic reflection profile reveals the crustal structure and deep-shallow structural features of faults at the junction of the northern margin of the Qinling orogenic belt and the southern margin of the North China block, which provides seismological evidence for the analysis of structural differences, the deep earth's interior processes and deep-shallow structural relationships between the Qinling-Dabie orogenic belt and the southern margin of the North China block. The lower crust of the study area is divided into two parts by deep faults that dislocate the Moho surface. These two parts have distinct reflective structures, suggesting that the area has experienced intense complex tectonic movements. The faults in the upper crust control the formation of basin-mountain structure and stratigraphic deposition of this area. And deep faults in the crust that disrupt Moho surface create conditions for the upwelling and energy exchange of deep materials. All of these have regulated the composition of material and the distribution of energy in the crust. The deep faults cutting through the lower crust and Moho surface and the south-dipping arc-shaped and inclined strong reflection sequences developed in the lower crust should indicate the large-scale subduction of the southern margin of the North China block towards the south-trending Qinling orogenic belt.

FOCAL MECHANISM SOLUTION AND TECTONIC STRESS FIELD CHARACTERISTICS OF THE MIDDLE TIANSHAN MOUNTAINS, XINJIANG

ZHANG Zhi-bin ZHAO Xiao-cheng REN Lin

2020, 42(3):  595-611.  DOI: 10.3969/j.issn.0253-4967.2020.03.004

Asbtract ( )   HTML   PDF (7799KB) ( )  

References | Related Articles | Metrics

The middle part of the Tianshan Mountains in Xinjiang is located in the north-central part of the Tianshan orogenic belt, between the rigid Tarim Basin and Junggar Basin. It is one of the regions with frequent deformation and strong earthquake activities. In this paper, 492 M_S>2.5 earthquake events recorded by Xinjiang seismograph network from 2009 to 2018 were collected. The M_S3.5 earthquake was taken as the boundary, the focal mechanism solutions of the earthquake events in this region were calculated by CAP method and FOCEMEC method respectively. At the same time the focal mechanism solutions of GCMT recorded historical earthquake events in this region were also collected. According to the global stress map classification standard, the moderate-strong earthquakes in the region are mainly dominated by thrust with a certain slip component, which are distributed near the combined belts of the Tarim Basin, Junggar Basin, Turpan Basin and Yili Basin with Tianshan Mountains. The thrust component decreases from south to north, while the strike-slip component increases. The spatial distribution characteristics of the tectonic stress field in the middle section of the Tianshan Mountains in Xinjiang are obtained by using the damped regional-scale stress field inversion method. The maximum principal compressive stress in axis the study area rotated in a fan shape from west to east, the NW direction in the western section gradually shifted to NE direction, its elevation angle is nearly horizontal, in the state of near horizontal compression. The minimum principal compressive stress axis is nearly EW, and the elevation angle is nearly vertical. Influenced by large fault zones such as Kashi River, Bolhinur, Nalati, Fukang, the southern margin of the Junggar and the north Beiluntai, the local regional stress field presents complex diversity. Under the influence of the northward extrusion of Pamir and Tarim blocks, the whole Tianshan is shortened by compression, but its shortening rate decreases from south to north and from west to east, the stress shape factor increases gradually from west to east, the intermediate principal compressive stress axis exhibits a change in compression to extension. There are some differences in the characteristics of tectonic stress field between the north and south of Tianshan Mountains. The regional maximum principal compressive stress axis is 15° north by east on the south side, while it is nearly NS on the north side. The deformation of the Tianshan Mountains and the two basins on both sides is obviously larger than that in the inside of the mountain. Changes in the crustal shortening rate caused by the rotation of the rigid Tarim block and Junggar block to the relatively soft Tianshan block, as well as the uplifts of Borokonu and Bogda Mountains, the comprehensive influence of the material westward expansion constitute the stress field distribution characteristics of the north and south sides of the middle section of Tianshan Mountains. The recent two M_S6.6 earthquakes in the region caused the regional stress field to rotate counterclockwise. The post-earthquake stress field and the main source focal mechanism solution tend to be consistent. The seismic activity in the study area is week in the south and strong in the north. The focal depth is about 20km. Most strike-slip earthquakes occur near the junction belt of the Tianshan and Junggar Basin.

NEAR-FAULT DISPLACEMENT AND DEFORMATION OBTAINED FROM ONE-KILOMETER-LONG FAULT-CROSSING BASELINE MEASUREMENTS-A PRELIMINARY EXPERIMENT AT 2 SITES ON THE EASTERN BOUNDARY OF THE SICHUAN-YUNNAN BLOCK

CAO Jian-ling, ZHANG Jing, WEN Xue-ze, FENG Wei, SHI Yao-lin

2020, 42(3):  612-627.  DOI: 10.3969/j.issn.0253-4967.2020.03.005

Asbtract ( )   HTML   PDF (3499KB) ( )  

References | Related Articles | Metrics

The current and conventional fault-crossing short baseline measurement has a relatively high precision, but its measurement arrays usually fail to or cannot completely span major active fault zones due to the short length of the baselines, which are only tens to 100 meters. GNSS measurement has relatively low resolution on near-fault deformation and hence is not suitable for monitoring those faults with low motion and deformation rates, due to sparse stations and relatively low accuracy of the GNSS observation. We recently built up two experimental sites on the eastern boundary of the active Sichuan-Yunnan block, one crossing the Daqing section of the Zemuhe Fault and the other crossing the Longshu section of the Zhaotong Fault, aiming to test the measurement of near-fault motion and deformation by using fault-crossing arrays of one-kilometer-long baselines. In this paper, from a three-year-long data set we firstly introduce the selection of the sites and the methods of the measurement. We then calculate and analyze the near-field displacement and strain of the two sites by using three hypothetical models, the rigid body, elastic and composed models, proposed by previous researchers. In the rigid body model, we assume that an observed fault is located between two rigid blocks and the observed variances in baseline lengths result from the relative motion of the blocks. In the elastic model, we assume that a fault deforms uniformly within the fault zone over which a baseline array spans, and in the array baselines in different directions may play roles as strainmeters whose observations allow us to calculate three components of near-fault horizontal strain. In the composed model, we assume that both displacement and strain are accumulated within the fault zone that a baseline array spans, and both contribute to the observed variances in baseline lengths. Our results show that, from the rigid body model, variations in horizontal fault-parallel displacement component of the Zemuhe Fault at the Daqing site fluctuate within 3mm without obvious tendencies. The displacement variation in the fault-normal component keeps dropping in 2015 and 2016 with a cumulative decrease of 6mm, reflecting transverse horizontal compression, and it turns to rise slightly(suggesting extension)in 2017. From the elastic model, the variation in horizontal fault-normal strain component of the fault at Daqing shows mainly compression, with an annual variation close to 10^-5, and variations in the other two strain components are at the order of 10^-6. For the Longshu Fault, the rigid-body displacement of the fault varies totally within a few millimeters, but shows a dextral strike-slip tendency that is consistent with the fault motion known from geological investigation, and the observed dextral-slip rate is about 0.7mm/a on average. The fault-parallel strain component of the Longshu Fault is compressional within 2×10^-6, and the fault-normal strain component is mainly extensional. Restricted by the assumption of rigid-body model, we have to ignore homolateral deformation on either side of an observed fault and attribute such deformation to the fault displacement, resulting in an upper limit estimate of the fault displacement. The elastic model emphasizes more the deformation on an observed fault zone and may give us information about localizations of near-fault strain. The results of the two sites from the composed model suggest that it needs caution when using this model due to that big uncertainty would be introduced in solving relevant equations. Level surveying has also been carried out at the meantime at the two sites. The leveling series of the Daqing site fluctuates within 4mm and shows no tendency, meaning little vertical component of fault motion has been observed at this site; while, from the rigid-body model, the fault-normal motion shows transverse-horizontal compression of up to 6mm, indicating that the motion of the Zemuhe Fault at Daqing is dominantly horizontal. The leveling series of the Longshu site shows a variation with amplitude comparable with that observed from the baseline series here, suggesting a minor component of thrust faulting; while the baseline series of the same site do not present tendencies of fault-normal displacement. Since the steep-dip faults at the two sites are dominantly strike-slip in geological time scale, we ignore probable vertical movement temporarily. In addition, lengths of homolateral baselines on either side of the faults change somewhat over time, and this makes us consider the existence of minor faults on either side of the main faults. These probable minor faults may not reach to the surface and have not been identified through geological mapping; they might result in the observed variances in lengths of homolateral baselines, fortunately such variations are small relative to those in fault-crossing baselines. In summary, the fault-crossing measurement using arrays with one-kilometer-long baselines provides us information about near-fault movement and strain, and has a slightly higher resolution relative to current GNSS observation at similar time and space scales, and therefore this geodetic technology will be used until GNSS networks with dense near-fault stations are available in the future.

SITE EFFECT AND EARTHQUAKE DISASTER CHARACTERISTICS IN GUANGZHOU AREA FROM HORIZONTAL-TO-VERTICAL SPECTRAL RATIO(HVSR)METHOD

ZONG Jian-ye, SUN Xin-lei, ZHANG Peng

2020, 42(3):  628-639.  DOI: 10.3969/j.issn.0253-4967.2020.03.006

Asbtract ( )   HTML   PDF (3191KB) ( )  

References | Related Articles | Metrics

Site effect, which is related to the amplification of seismic waves, is mostly affected by the thickness and softness of sediment layers. The study of site effect in cities is becoming more and more important to the assessment of ground motions, seismic hazard and engineering seismology. However, in highly populated urban cities, traditional seismic surveys cannot be applied extensively due to their destructive consequences and high cost. The ambient noise, including microseisms and microtremor, could be acquired anywhere and anytime, and thus can serve as an effective source for engineering seismology. In order to get the site effect and distribution of sedimentary layers of Guangzhou area, one hundred QS-05A seismographs with frequency bandwidth of 5s-250Hz were deployed in early 2018 for 35 days. The inter-station spacing of these seismometers is approximately 2~5km. Using continuous ambient noise signals, we obtained the resonance frequency and amplification value beneath each station by horizontal-to-vertical spectral ratio(HVSR)method. Then sedimentary layer thicknesses as well as K-values, which are related to the site vulnerability to ground shaking, were calculated. Our results suggest that the resonance frequencies in Guangzhou area are between 1~6.5Hz. The resonance frequencies increase gradually from 1Hz on the north-east side to 6.5Hz on the south-west side of the study area. The sediment thicknesses change from several meters to about 40m, with the maximum thickness at around the estuary of the Pearl River. This distribution is consistent with the topography. The amplifications are mainly between 2~6. The largest amplification is around the Pearl River and the west part of Baiyun District. In general, the K-values are small(<20), less than the dangerous value, suggesting that Guangzhou area is relatively safe in ground shaking. However, there are three small areas beneath Huadu District, Sanshui District and Nanhai District. They all have K-values greater than 20, suggesting those areas are more vulnerable to earthquake destruction, and higher construction standard is needed. The reliability of our results is further supported by its consistency with topography and borehole data in Guangzhou area. Our results provide important information for shallow underground structure in Guangzhou area, and can be referred as guidelines in urban architecture planning and disaster prevention and mitigation.

MAGNETOTELLURIC IMAGING OF MAGMA DISTRIBUTION BENEATH MA'ANLING AND LEIHULING VOLCANOES OF NORTHEASTERN HAINAN, CHINA

SUN Xiang-yu, ZHAN Yan, ZHAO Guo-ze, ZHAO Ling-qiang, DENG Yan, HU Ya-xuan, HU Jiu-chang, XIANG Xiao-juan

2020, 42(3):  640-653.  DOI: 10.3969/j.issn.0253-4967.2020.03.007

Asbtract ( )   HTML   PDF (7170KB) ( )  

References | Related Articles | Metrics

The northeastern Hainan Province is one of the areas subjected to the strongest, most frequent and longest-lasting volcanic activities in China since the Cenozoic era. Under the influence of magma and fault activities, northeastern Hainan Island has experienced many moderate and strong earthquakes in history. The Qiongshan M7.5 earthquake occurred in this region in 1605. The deformation measurement and InSAR data found a subsidence area in the south of the Qiongshan M7.5 earthquake. Small earthquakes frequently occur in this area. It has been inferred by some studies on this subsidence area, namely the Puqianwan-Fengjiawan seismic belt, that the subsidence and frequent seismic activity are related to the development of deep magma systems. Magnetotelluric methods are very sensitive to subsurface fluid, different temperature conditions, and resistivity property of the medium in the molten state. With the development of magnetotelluric three-dimensional inversion technique, using dense array magnetotelluric data in three-dimensional inversion can image the medium resistivity occurrence state and position in the volcanic area. To study the deep structure of the magma system and its relationship with seismic activity, we conducted MT observations on two profiles that cross Leihuling and Ma'anling volcanoes. Phase tensor decomposition was used to analyze the electrical structure. This paper investigates the two MT profiles using three-dimensional electromagnetic imaging technology and obtains the electrical structure of the two profiles. The result reveals the media properties and high conductivity bodies' occurrence range beneath the volcanic area in the northeastern Hainan. There are obvious differences in the electrical structure of the northeastern Hainan. The resistivity values are high in the east and low in the west. In addition, there are two high conductivity bodies in the northeast of Hainan. The high conductivity body C₁ inclines to the west and locates beneath the Chengmai County area in the northwestern Hainan Island(west of the Leihuling-Ma'anling volcanoes). Its resistivity value is less than several Ωm. This low resistive body is 40km long in WE direction and 30km wide in SN direction. Its burial depth is about 2km near the HNL1 profile and 6km near the NHNL1 profile. Its bottom reaches the depth of about 25~30km, which may be close to or through the Moho surface depth of 25~26km in this area. It is speculated that the magma eruption of Leihuling-Ma'anling volcanoes did not migrate vertically from its deep part to the surface. The high conductivity body C₂ locates beneath Longquan. The buried depth of C₂ tends to be shallower from north to south, but there is no exposed surface in the study area, nor is it connected with the shallow low-resistivity layer. It is speculated that the C₂ may be a magmatic sac trapped in the crust, but may have nothing to do with the eruption of Ma'anling-Leihuling volcanoes. The recent volcanic magma in this area comes from the lower crust and upper mantle of the ocean area to the west of Hainan Island. As magma enters the upper and middle crust, it continues to move shallowly and eastward. In this process, it should be blocked by the high resistance structure on the east side of the Changliu-Xiangou Fault and then erupt around this fault, thus forming numerous craters in this area. After the repeated eruption, deep magma channels gradually closed and volcanic activity weakened. The magma in the mid-upper crust cooled consolidated gradually, but the speed was uneven in different areas, resulting in the channels having closed down gradually in some places, and some are in the process of closing. Our results show an uneven rise and fall depth of the low resistivity body in the middle and lower crust. There is no high conductivity body in the deep part of the Puqianwan-Fengjiawan seismic belt and the subsidence area in the northeastern Hainan, which rules out the possibility that the small earthquakes are related to deep magma systems.

THE APPLICATION OF PML BOUNDARY CODITIONS IN THE SIMULATION OF SEISMIC WAVE BY THE HIGH-ORDER STAGGERED-GRID FINITE DIFFERENCES AT THE TRANSVERSELY ISOTROPIC MEDIA

CHEN Jie, ZHU Shou-biao

2020, 42(3):  654-669.  DOI: 10.3969/j.issn.0253-4967.2020.03.008

Asbtract ( )   HTML   PDF (4665KB) ( )  

References | Related Articles | Metrics

In the realm of the numerical simulation, finite difference method and finite element method are more intuitive and effective than other simulation methods. In the process of simulating seismic wave propagation, the finite differences method is widely used because of its high computational efficiency and the advantage of the algorithm is more efficient. With the demand of precision, more and more researchers have proposed more effective methods of finite differences, such as the high-order staggered-grid finite differences method, which can restore the actual process of wave propagation on the premise of ensuring accuracy and improving the efficiency of operation. In the past numerical simulation of seismic wave field, different models of isotropic medium are mostly used, but it is difficult to reflect the true layer situation. With the research demand of natural seismology and seismic exploration, the research on anisotropic media is more and more extensive. Transversely isotropic(TI)media can well simulate the seismic wave propagation in the formation medium, such as gas-bearing sandstone, mudstone, shale et al., the character of TI media is reflected by introducing the Thomsen parameters to reflect its weak anisotropy of vertical direction by using Thomson parameter. Therefore, studying the process of seismic wave propagation in TI media can restore the true information of the formation to the greatest extent, and provide a more reliable simulation basis for the numerical simulation of seismic wave propagation. In the geodynamic simulation and the numerical simulation of the seismic wave field, under the limited influence of the calculation area, if no boundary conditions are added, a strong artificial boundary reflection will be generated, which greatly reduces the validity of the simulation. In order to minimize the influence of model boundaries on the reflection of seismic waves, it is often necessary to introduce absorbing boundary conditions. At present, there are three types of absorption boundary conditions: one-way wave absorption boundary, attenuation absorption boundary, and perfectly matched layer(PML)absorption boundary. In terms of numerical simulation of seismic waves, the boundary absorption effect of PML is stronger than the first two, which is currently the most commonly used method, and it also represents the cutting-edge development direction of absorption boundary technology. The perfectly matched layer absorbing boundary is effectively applied to eliminating the reflective waves from model boundaries, but for transversely isotropic medium, the effect of the absorbing is not very well. For this reason, the elastic dynamic wave equations in transversely isotropic media are derived, and we describe a second-order accurate time, tenth-order accurate space, formulation of the Madariaga-Virieux staggered-grid finite difference methods with the perfectly matched layer(PML)are given. In addition, we have established vertical transversely isotropic(VTI)media and arbitrary inclined tilted transversely isotropic(TTI)media models, using a uniform half-space velocity model and a two-layer velocity model, respectively. By combining the actual geoscience background, we set the corresponding parameters and simulation conditions in order to make our model more research-oriented. When setting model parameters, different PML thickness, incident angle, source frequency and velocity layer models were transformed to verify the inhibition of boundary reflection effect by PML absorption boundary layer. The implementations of this simulation show that the formula is correct and for the transversely isotropic(TI)media of any angular symmetry axis, when the thickness of the PML layer reaches a certain value, the seismic wave reflection effect generated by the artificial boundary can be well suppressed, and the absorption effect of PML is not subject to changes in incident angle and wave frequency. Therefore, the results of our study indicate that our research method can be used to simulate the propagation process of seismic waves in the transversely isotropic(TI)media without being affected by the reflected waves at the model boundary to restore the actual formation information and more valuable geological research.

DEVELOPMENT OF TRANSVERSE DRANAGES AND FORMATION OF WIND GAPS ON ACTIVELY GROWING FOLD: REVIEW AND CASE STUDY

CAO Xi-lin, GENG Hao-peng, PAN Bao-tian, HU Xiao-fei

2020, 42(3):  670-687.  DOI: 10.3969/j.issn.0253-4967.2020.03.009

Asbtract ( )   HTML   PDF (3996KB) ( )  

References | Related Articles | Metrics

The most compelling phenomena for transverse drainage in active fold belt are lateral diversion of channels and development of water/wind gaps. This phenomenon is the result of competition between uplift and erosion, which is controlled by fault vertical/lateral propagation and segment linkage, fault geometry, climate condition and lithology. Previous studies found that the higher the uplift rate is, the greater number of wind gaps form, and the variation of the uplift rate is also critical to the sustainability of transverse rivers. Lateral propagation and linkage of several separate folds in fold-and-thrust belts will lead to defeat of streams and diversion into a trunk drainage; if the trunk is still unable to keep pace with uplift, water gap will be abandoned and left as a wind gap. For lateral propagation of an anticline associated with development of tear faults, the locations of wind/water gaps are likely to coincide with the trace of tear fault and it's not quite clear about the relation between tear faulting and stream deflection. Nonzero dip of the underlying detachment induces a lateral surface slope in the direction of fault propagation, which in turn makes rivers deflection more efficient. Climate and rock erodibility control the water/sediment discharge, and further influence river transport/incision capacity. The changing climate and rock erodibility conditions enable river to abandon the current waterway to create a wind gap unless they could down-cut through a growing fold. However, the role of climate cycle in the formation of wind gap is still controversial. In addition, wind gaps are commonly developed along the divides where parts of longitudinal river have been captured by transverse catchments. Generally, the development of transverse drainages and the formation of wind gaps in nature are result from a combination of tectonic and fluvial process. The wind gap pattern and transverse drainage evolution in fold-and-thrust belts contain plenty of information on fault growth, interaction between tectonic uplift and fluvial erosion, and development of sedimentary basin. Such researches have significant implications in geomorphology, seismic hazard assessment and hydrocarbon exploration. However, there are still many knowledge gaps on the study of transverse river evolution in active fold areas. First, adequate chronology and geomorphic/strata mark to quantify fold growth and erosion is commonly not available, which leads to a poorly constrained rate in both river incision and lateral propagation of growing folds. In addition, more geological and geomorphological processes could influence the evolution of transverse drainages. For examples, (1)during the formation of a young range or anticline, the mechanism of fault-related folding may change over time, e.g. from fault-propagation folding to surface breaking; (2)Besides the knickpoint retreat in downstream, efficient lateral planation and downstream sweep erosion are also important in understanding the erosion of folds by rivers flowing through it. These processes make the development of transverse drainage across folds more complex and should be considered in more comprehensive models. There are lots of rivers originating from the Tibetan plateau and cutting through young surrounding mountains. These surrounding mountains, such as Qilian Mountains, Tianshan Mountains and Longmen Mountains, are ideal areas for the study of transverse river evolution and wind gap formation. In the end, combining with the geological and geomorphological features of the Heli Shan-Jintanan Shan, north of Hexi Corridor, we propose that the Heihe River has experienced deflection, beveling and incision since Mid Pleistocene. These processes have led to 1)the formation of a wind gap on the western Heli Shan, 2)a layer of fluvial gravels from the Qilian Shan preserved on the top surface of the Jintanan Shan, and overlying angular unconformity upon older strata, and 3)the incision of the Heihe River to form the Zhengyi Gorge through the linked structure between Heli Shan and Jintanan Shan. Thus, we propose a general model for the development of transverse drainages in the central Hexi Corridor: deflection-beveling-incision.

THE ACTIVITY FEATURES OF XIADIAN FAULT ZONE REVEALED BY RONGJIABAO TRENCH AND ITS PROBABILISTIC SEISMIC HAZARD EVALUATION

YU Zhong-yuan, PAN Hua, SHEN Jun, LI Jin-chen, ZHANG Meng, DAI Xun-ye

2020, 42(3):  688-702.  DOI: 10.3969/j.issn.0253-4967.2020.03.010

Asbtract ( )   HTML   PDF (3200KB) ( )  

References | Related Articles | Metrics

The time-dependent probabilistic seismic hazard assessment of the active faults based on the quantitative study of seismo-geology has the vital practical significance for the earthquake prevention and disaster management because it describes the seismic risk of active faults by the probability of an earthquake that increases with time and the predicted magnitude. The Poisson model used in the traditional probabilistic method contradicts with the activity characteristics of the fault, so it cannot be used directly to the potential earthquake risk evaluation of the active fault where the time elapsing from the last great earthquake is relatively short. That is to say, the present Poisson model might overestimate the potential earthquake risk of the Xiadian active fault zone in North China because the elapsed time after the historical M8 earthquake that occurred in 1679 is only 341a. Thus, based on paleoearthquake study and geomorphology survey in the field, as well as integrating the data provided by the previous scientists, this paper reveals two paleo-events occurring on the Xiadian active fault zone. The first event E₁ occurred in 1679 with magnitude M8 and ruptured the surface from Sanhe City of Hebei Province to Pinggu District of Beijing at about 341a BP, and the other happened in (4.89±0.68)ka BP(E₂). Our research also found that the average co-seismic displacement is ~(1.4±0.1)m, and the predicted maximum magnitude of the potential earthquake is 8.0. In addition, the probabilistic seismic hazard analysis of great earthquakes for Xiadian active fault zone in the forthcoming 30a is performed based on Poisson model, Brownian time passage model(BPT), stochastic characteristic-slip model(SCS)and NB model to describe time-dependent features of the fault rupture source and its characteristic behavior. The research shows that the probability of strong earthquake in the forthcoming 30a along the Xiadian active fault zone is lower than previously thought, and the seismic hazard level estimated by Poisson model might be overestimated. This result is also helpful for the scientific earthquake potential estimation and earthquake disaster protection of the Xiadian active fault zone, and for the discussion on how to better apply the time-dependent probabilistic methods to the earthquake potential evaluation of active faults in eastern China.

LATE QUATERNARY ACTIVITY OF THE CENTRAL SEGMENT OF THE DARI FAULT AND RESTUDY OF THE SURFACE RUPTURE ZONE OF THE 1947 M73/4 DARI EARTHQUAKE, QINGHAI PROVINCE

LIANG Ming-jian, YANG Yao, DU Fang, GONG Yue, SUN Wei, ZHAO Min, HE Qiang

2020, 42(3):  703-714.  DOI: 10.3969/j.issn.0253-4967.2020.03.011

Asbtract ( )   HTML   PDF (10839KB) ( )  

References | Related Articles | Metrics

Bayan Hara Block is one of the most representative active blocks resulting from the lateral extrusion of Tibet Plateau since the Cenozoic. Its southern and northern boundary faults are characterized by typical strike-slip shear deformation. Its eastern boundary is blocked by the Yangze block and its horizontal movement is transformed into the vertical movement of the Longmen Shan tectonic belt, leading to the uplift of the Longmen Shan Mountains and forming a grand geomorphic barrier on the eastern margin of the Tibet Plateau. A series of large earthquakes occurred along the boundary faults of the Bayan Hara Block in the past twenty years, which have attracted attention of many scholars. At present, the related studies of active tectonics on Bayan Hara Block are mainly concentrated on the boundary faults, such as Yushu-Ganzi-Xianshuihe Fault, East Kunlun Fault and Longmen Shan Fault. However, there are also some large faults inside the block, which not only have late Quaternary activity, but also have tectonic conditions to produce strong earthquake. These faults divide the Bayan Hara Block into some secondary blocks, and may play important roles in the kinematics and dynamics mechanism of the Bayan Hara Block, or even the eastern margin of the Tibet Plateau. The Dari Fault is one of the left-lateral strike-slip faults in the Bayan Hara Block. The Dari Fault starts at the eastern pass of the Kunlun Mountains, extends eastward through the south of Yalazela, Yeniugou and Keshoutan, the fault strike turns to NNE direction at Angcanggou, then turns to NE direction again at Moba town, Qinghai Province, and the fault ends near Nanmuda town, Sichuan Province, with a total length of more than 500km. The fault has been considered to be a late Quaternary active fault and the 1947 M73/4 Dari earthquake was produced by its middle segment. But studies on the late Quaternary activity of the Dari Fault are still weak. The previous research mainly focused on the investigation of the surface rupture and damages of the 1947 M73/4 Dari earthquake. However, there were different opinions about the scale of the M73/4 earthquake surface rupture zone. Dai Hua-guang(1983)thought that the surface rupture of the earthquake was about 150km long, but Qinghai Earthquake Agency(1984)believed that the length of surface rupture zone was only 58km. Based on interpretation of high-resolution images and field investigations, in this paper, we studied the late Quaternary activity of the Dari Fault and the surface rupture zone of the 1947 Dari earthquake. Late Quaternary activity in the central segment of the Dari Fault is particularly significant. A series of linear tectonic landforms, such as fault trough valley, fault scarps, fault springs and gully offsets, etc. are developed along the Dari Fault. And the surface rupture zone of the 1947 Dari earthquake is still relatively well preserved. We conducted a follow-up field investigation for the surface rupture zone of the 1947 Dari earthquake and found that the surface rupture related to the Dari earthquake starts at Longgen village in Moba town, and ends near the northwest of the Yilonggounao in Jianshe town, with a length of about 70km. The surface rupture is primarily characterized by scarps, compressional ridges, pull-apart basins, landslides, cleavage, and the coseismic offset is about 2~4m determined by a series of offset gullies. The surface rupture zone extends to the northwest of Yilonggounao and becomes ambiguous. It is mainly characterized by a series of linear fault springs along the surface rupture zone. Therefore, we suggest that the surface rupture zone of the 1947 Dari earthquake ends at the northwest of Yilonggounao. In summary, the central segment of the Dari Fault can be characterized by strong late Quaternary activity, and the surface rupture zone of the 1947 Dari earthquake is about 70km long.

EXPERIMENTAL STUDY ON THE CHANGES OF ULTRASONIC CODA WAVE AND ACOUSTIC EMISSION DURING ROCK LOADING AND DEFORMATION

YANG Hai-ming, CHEN Shun-yun, LIU Pei-xun, GUO Yan-shuang, ZHUO Yan-qun, QI Wen-bo

2020, 42(3):  715-731.  DOI: 10.3969/j.issn.0253-4967.2020.03.012

Asbtract ( )   HTML   PDF (5285KB) ( )  

References | Related Articles | Metrics

The coda wave propagation path has received extensive attention as it is more sensitive to small changes in the medium than the direct wave. During the process of loading, the wave velocity, medium or source changes may cause the coda wave to change. The physical mechanism of change in the ultrasonic coda wave varies during different deformation stages. Meanwhile, there exist local damages in the rock sample during the deformation, and it will be accompanied by acoustic emission. Combining the ultrasonic coda wave and acoustic emission is beneficial to characterize the coda wave characteristics and damage degree of the sample at different deformation stages. In this paper, three kinds of rocks, including granodiorite, marble and sandstone with the sizes of 50mm×50mm×150mm, are used to carry out observations of ultrasonic coda wave and acoustic emission during the whole process of loading so as to study characteristics of the coda wave at different deformation stages. The major results are given below: 1)There is a good correspondence between the coda wave variation and the acoustic emission evolution process. When the acoustic emission frequency increases, the coda wave changes accordingly. In particular, the coda wave changes in the early stages of increased acoustic emission frequency, which indicates that the early damage information of rock can be obtained by analysis of the coda wave. 2)The physical mechanism of the coda wave change is different in different deformation stages. At the initial stage of loading, there are obvious scatterer changes in the coda wave change; then, in the linear elastic deformation stage, the wave velocity change is dominating; in the late-stage of loading, the scatterer change increases and coexists with the wave velocity change, the scatterer change effect is related with the rock micro-fracture degree, the rock will locally be damaged before rupturing, and the role of the scatterer will be enhanced. 3)With the increase of loading, the amplitude of increase of the wave velocity generally decreases gradually, which is basically consistent with the understanding obtained through the direct wave. The interference of acoustic emission can be eliminated because of the Kaiser effect when analyzing the coda wave. The consistency of the wave velocity change and stress loading and unloading is further verified. 4)The micro-fracture generated during rock deformation will change the physical mechanism of the coda wave change, and the scatterer effect will be significantly enhanced. At the same time, the acoustic emission waveform will cause interference to the ultrasonic coda wave. This means that attention needs to be paid when analyzing rock damage using only coda wave data. In short, the ultrasonic coda wave and acoustic emission can reflect the damage inside the rock, and the change mechanism of the coda wave in different deformation stages is different. The joint observation of the two can play a mutual verification role, which is conducive to improving the reliability of the observation results.

COMPARISON STUDY OF TWO KINDS OF CODES TO MEASURE FAULT-OFFSETS BASED ON MATLAB: A CASE STUDY ON EASTERN ALTYN TAGH FAULT

KANG Wen-jun, XU Xi-wei, YU Gui-hua, LUO Jia-hong

2020, 42(3):  732-747.  DOI: 10.3969/j.issn.0253-4967.2020.03.013

Asbtract ( )   HTML   PDF (6975KB) ( )  

References | Related Articles | Metrics

Geomorphic offsets displaced by coseismic surface rupture can be analyzed to identify earthquake recurrence behavior. Therefore, obtaining a sufficient and precise along-fault offset dataset is vital to identify long-term earthquake recurrence behavior. Furthermore, knowledge of along-fault slip distribution during a single-earthquake or multi-earthquakes is important for other reasons, including a better understanding of the relationship between earthquake size and coseismic displacements, fault kinematics and fault mechanics. A recent flourish of offsets-measuring software and high-resolution topographic data together offer an unprecedented opportunity to measure high-density fault offsets. Here, we introduce and compare two kinds of most popular software, LaDiCaoz and 3D_Fault_Offsets. We describe the workflow and principle of the two codes by taking a fault-offset example on the eastern Altyn Tagh Fault. LaDiCaoz iterates over the channel morphology and position parameters and determines the summed absolute elevation difference Σ[Δ(elevation)] between both transverse profiles. The optimal horizontal offset is defined by the parameter combination that results in the least mismatch between two profiles. Compared with LaDiCaoz, the principle of 3D_Fault_Offsets is more complicated by measuring the offset in three dimensions. It mathematically identifies and represents nine of the most prominent geometric characteristics of common sublinear markers along faults in three dimensions, such as the streambed(minimum elevation), top, free face and base of channel banks or scarps(minimum Laplacian, maximum gradient, and maximum Laplacian), and ridges(maximum elevation). By calculating best fit lines through the nine point clouds on either side of the fault, the code computes the lateral and vertical offsets between the piercing points of these lines onto the fault plane, providing nine lateral and nine vertical offset measures per marker. Through a Monte Carlo approach, the code calculates the total uncertainty on each offset. Although both 3D_Fault_Offsets and LaDiCaoz are developed based on the Matlab platform, there are significant differences in principles, linear marker, software interface, repeatability, input-file types, degree of automation, adaptability, output file types, etc. In this part, we compare and summarize their features, advantages, and disadvantages. Finally, we calculate the correlation of two groups of fault-offset data derived from the two methods along the eastern ATF. By doing this, we try to explore if the two methods can be crosschecked and to study how sinuosity of the linear geomorphic markers affect the measuring results. By discussing and comparing the accuracy of the two measuring methods, we consider that LaDiCaoz is better than 3D_Fault_Offsets in accuracy aspect. In our opinion, there exist some disadvantages in the both software, and higher automation and introduction of artificial intelligence will be the future development direction.

RESEARCH ON EARTHQUAKE EMERGENCY MAPPING METHOD BASED ON TEMPLATE MATCH

XU Jing-hai, ZHOU Hai-jun, NIE Gao-zhong, AN Ji-wen

2020, 42(3):  748-761.  DOI: 10.3969/j.issn.0253-4967.2020.03.014

Asbtract ( )   HTML   PDF (9194KB) ( )  

References | Related Articles | Metrics

After an earthquake, earthquake emergency response and rescue are important ways to mitigate earthquake-induced losses. Various earthquake emergency maps can provide effective references and guidance to those actions. Currently, related studies include the investigation on symbols of emergency maps, remote sensing emergency mapping and GIS-based mapping methods. However, the existing studies overlook the characteristics of rapidity, dynamicity and variety of presentation methods in making earthquake emergency maps. In this paper, a map template matching method is used to quickly make earthquake emergency maps considering their characteristics. We take investigations on the service objects(users)of the earthquake emergency maps to understand the needs of making earthquake emergency maps. The audience theory in mass media field and map information transmission theory are adopted to classify the users of the earthquake emergency maps into four categories: earthquake emergency commanders, technical staffs for decision-making, earthquake emergency rescuers, and the public. The components of different types of users are described and then their diverse demands in earthquake emergency maps are analyzed, such as the needs of on-field disaster information maps, earthquake information maps, physical geography and social economic maps. Following those needs, we introduce the representation methods of the earthquake emergency maps according to their formats(vector or raster)and contents, such as point symbolization method, kilometer grid method, line symbolization method and range method. Then, we study the rapid plotting method of earthquake emergency map based on map template matching method. The core steps of the method include: 1)before earthquake, the templates of different earthquake emergency maps are designed, prepared and connect the earthquake emergency features with their related spatial database. The map layout and map elements are stored in the templates. 2)After earthquake, the earthquake emergency features will be generated from seismic models(such as attenuation model of earthquake magnitude and seismic intensity)or the information obtained from field investigation. 3)Corresponding earthquake emergency map template is selected in accordance with the generated seismic features. And the features are used to update related features inside the selected template. 4)Minor adjustments are made such as to the map scale and some map annotations to finally generate the formal earthquake emergency map. Architecture of template system of the earthquake emergency maps is designed, including map user level, map template level, template layer level and map element level. Regrading to the architecture, the general map template of earthquake emergency is presented which includes four main regions: title region, main picture region, auxiliary region and annotation region. The main picture region is the essential, which lays geographic background maps and earthquake emergency features. Finally, an earthquake emergency mapping system is developed. Based on the system, a case study is presented, which demonstrates making a simulated seismic intensity influence map. From three aspects, the case presents the application of the template-matching method including: generating earthquake emergency features, substituting the features inside the template with the generated features, and revising map annotations. Therefore, the map template matching method is verified so that it can be used to quickly generate various earthquake emergency maps.

INVESTIGATION OF STATUS AND DEMAND ON PUBLIC SERVICE OF EARTHQUAKE DISASTER MITIGATION IN CHINA

DONG Li-na, LIAN Wei-ping, CHEN Wei-tao, MA Hai-jian

2020, 42(3):  762-771.  DOI: 10.3969/j.issn.0253-4967.2020.03.015

Asbtract ( )   HTML   PDF (3243KB) ( )  

References | Related Articles | Metrics

In order to make a comprehensive and systematic evaluation of the status and demands of public service of earthquake disaster mitigation, we launched an investigation with Social Opinion Polling Center of National Bureau of Statistics on the current situation of public service of earthquake disaster mitigation and the demands for it in 2018. The investigation was conducted for 18~75-year-old residents for both urban and rural areas in 31 provinces of China mainland using the Computer Assisted Telephone Interview(CATI)method and 20 078 valid samples were received in total. The questionnaire mainly includes two parts: the current situation assessment and the demand survey. The assessment of public service status consists of five indicators: earthquake emergency drill, knowledge publicity of earthquake prevention and disaster reduction, acquisition of earthquake information, reliability of information released by earthquake agencies, and overall satisfaction. The results indicate that only 27.4% of public have participated the earthquake emergency drill, and 73.5% of public have learned some knowledge about earthquake prevention and disaster reduction in different level. It shows that rural residents have a lower level of scientific knowledge of earthquake disaster mitigation although the rural residents have a higher level of satisfaction. Meanwhile, 82.9% of the public are very concerned about earthquake-related information, and the earthquake information acquisition convenience score is 81.51 points(the full score is 100 points). At present, TV is still the most popular way to obtain the earthquake-related information, besides that, network and the new media app such as Wechat and Weibo play a more and more important role. In terms of the reliability of official information published by the earthquake administration, the score is 83.69 points which indicates that the public tend to believe in official departments. Furthermore, the degree of satisfaction evaluation of public service of earthquake disaster mitigation is marked 76.44 points which shows that there is still much room for improvement. The more educated the public, the less satisfied they are with the public service of earthquake prevention and disaster reduction. The results also show that 51.81% of the public are willing to obtain earthquake warning information at their own expense. As for the demand, the most expectative services for the public are post-earthquake rescue and reconstruction, earthquake warning, and house earthquake resilience test. The results of the investigation may provide some reference for the government to provide better services of earthquake disaster mitigation and to make scientific knowledge popularization in a more scientific way.