Abstract

In this paper, in order to confirm the strong activity of the Yaziba Fault and provide information needed for disaster prevention and mitigation and safety of the reservoir, we use field investigations and remote sensing interpretations to obtain the geometry of the Yaziba Fault and the distribution of the associated nearby landslides. Several new events are dated, and the paleo earthquake events in the fault area and their time intervals are analyzed along with previously dated events. Then, a typical landslide case is analyzed, and its instability is determined using the limit equilibrium method. The main results of our study are as follows: (1) the Yaziba Fault extends for 20 km; (2) 164 landslides have occurred in the fault area, the LND reached 0.1 km−2, and the LAP reached 0.07 for the entire research area, and the LND and LAP gradually decrease with increasing distance from the Yaziba Fault; (3) 39 dating dates are obtained and 19 paleo earthquakes are concluded, from which the time interval of earthquakes shortened from 147,000 to 3800 years ago is found; and (4) the Daomakan landslide had a 0.6 g horizontal seismic acceleration, which indicates that it was an earthquake landslide. This comprehensive method combined the geometric characteristics of the fault, paleo earthquake data, the regional landslide distribution, and the geotechnical engineering analysis of typical landslides in the area to indirectly determine the fault’s activity and provide a reference for the study of fault activity in mountainous areas.

1. Introduction

Earthquakes in mountainous areas and the landslides they trigger pose a significant risk to human life and property [14]. For example, (i) in 1786, an M 7.75 earthquake occurred in the Kangding-Luding area in southwestern China, causing a landslide that blocked the Dadu River and formed a barrier lake [5]; (ii) in 1933, an M 7.5 earthquake in Diexi Town, southwestern China, triggered a large number of landslides and formed more than 10 landslide dams and dammed lakes [6]; and (iii) in 2008, the Wenchuan M 7.9 earthquake produced a large number of landslides and barrier lakes [7]. Therefore, the medium-long-term prediction of earthquakes is of great significance to disaster prevention and reduction in mountainous areas.

Fault activity evaluation is the most important basis for the medium-long-term prediction of earthquakes. Most research methods use the deep geophysical field [8] and the geometric characteristics of fault, including the length of the surface rupture, the fault offset [9], the timing of fault activity [10], and surface deformation [11], to qualitatively evaluate the fault. Some scholars have tried to identify quantitative indicators and make quantitative evaluations using probability theory [12, 13]. Two important indicators of fault activity are intensity, which is estimated from the magnitude of the earthquake, and frequency, which is determined from the timing of fault activity. The activity intensity is usually calculated based on historical seismic records [14, 15], but it can also be analyzed from seismic traces, e.g., landslides [16, 17], buildings destroyed by earthquakes, [18], and standing orphans [19]. These seismic traces are used to calculate the seismic peak acceleration, which is used to determine the possible seismic fault events and their magnitudes as well as tectonic environment of the site [18, 2022]. The timing of fault activity is mainly obtained through various dating methods, such as the direct testing of fault gouges in the fault zone [10, 23] or the indirect testing of earthquake artifacts. Frequently used earthquake artifacts include ancient landslides [24, 25] and sediments from barrier lakes [26, 27]. These earthquake artifacts, which can provide information about fault activity, seismic landslides, and barrier lakes, are undoubtedly the most common and reliable ones because large earthquakes in mountain areas often produce a large number of large landslides, which can block rivers and form barrier lakes. These artifacts provide abundant information we can use to understand the related fault activity. Nevertheless, the study of fault activity using seismic landslides is still not well developed, and more case studies are needed to promote the further development of this method.

The Yaziba Fault is located on the Jinsha River between Sichuan and Yunnan only 19 km downstream of the Xiluodu hydropower station. This hydropower station is the largest station on the Jinsha River, has the second largest power generation capacity in China, and is ranked as the third largest in the world. Its dam is 285.5 m high, which makes it the third highest arch dam in China. The Yaziba Fault intersects the Jinsha River several times, so there are several well-developed ancient landslides near the river. Due to the construction of large hydropower stations on the river, there have been many studies on the geology, earthquakes, and landslides in this area, but several questions regarding the seismic hazards associated with this fault remain unsolved. In this paper, geological surveys and tracing and remote sensing interpretations of ancient landslides are used to determine the geometry of the fault and the nearby associated landslides. A variety of geological dating methods, along with previous dates, are employed to infer the distribution of paleo earthquakes around the fault. The results provide valuable information that can increase our understanding of the seismicity, earthquake prediction, and safe operation of hydropower stations in the region. The research methods used can also provide a reference for fault activity research in other similar mountainous areas.

2. Study Area

The lower part of the Jinsha River (Panzhihua-Yibin) in China is 782 km long, with a natural drop in the elevation of 729 m, which provides abundant hydropower resources. Along this section of the river, four hydropower stations have been constructed in a cascading manner: the 8.7 million·kW Wudongde power station, the 12 million·kW Baihetan power station, the 12.6 million·kW Xiluodu power station, and the 6 million·kW Xiangjiaba power station (Figure 1).

The Yaziba Fault is located between the Xiluodu and Xiangjiaba power stations. Tectonically, the fault is located on the slope of the transition from western Sichuan and the northern Yunnan Plateau on the eastern margin of the Tibetan Plateau to the Sichuan Basin in the east. This region ranges in elevation from 300 to 3000 m, making it an intermediate-high mountainous area with a large topographic relief. The Jinsha River flows from south to north. The strong incision of the river has produced steep terrain with many deep narrow valleys with cliffs in many areas. Most of the Jinsha River valleys are narrow V-shaped or U-shaped valleys with bank slope gradients of 30°–60° and maximum slope heights of up to 600 m. The Carboniferous, Devonian, Cretaceous, and Tertiary strata are absent, but all of the other strata from the Upper Proterozoic to the Cenozoic are exposed in the area. The lithology of the area is dominated by limestone, sandstone, dolomite, basalt, and shale.

The Yingjing-Mabian-Yanjin fault zone, which contains the Yaziba Fault, is a complex seismic zone on the eastern side of the boundary between Sichuan and Yunnan, which strikes NNW and extends downward to the Moho interface. Based on its surface expression, it is composed of nine major secondary faults (Figure 2), i.e., the Tianquan-Xingjing, Ebian, Lidian, Zhongdu, Manao, Yaziba, Guancun, Zhongcun, and Dianlanba Faults. In the southern section, the Dianlanba, Manao, Yaziba, Zhongdu, Guancun, and Zhongcun Faults collectively make up the Mabian-Daguan fault zone (Figure 2). Many studies have been conducted on this fault zone. The Yaziba Fault is an important component of this feature (Figure 2), which extends approximately parallel to the Jinsha River and crosses the river twice. Previous dating shows that it has been active since the Quaternary, and a large number of ancient landslides and barrier lakes have been identified on both sides of the fault.

Since 1900, 30 M 5.0 or greater earthquakes have occurred on the Yingjing-Mabian-Yanjin fault zone, including 6 M 6.0–6.9 events and one M 7.1 event, i.e., the 1974 Daguan-Yongshan shock. Major earthquakes have also occurred near the Yaziba Fault, e.g., the 1216 M 7.0 Mahu and 1844 M 7.0 Daguanbei earthquakes (Figure 2).

3. Methods and Data

3.1. Investigation of the Yaziba Fault and Landslides

In this study, we had determined the geometry of the fault and landslides based on field investigations and remote sensing interpretations. We used observations and mapping to investigate the Yaziba Fault along strike and for several cross sections at key sites. For the landslides, their planar ranges were delineated using remote sensing interpretations, which were then corrected based on on-site observations and measurements. In addition, the lithology and structure of the landslides were determined. The landslides were interpreted using the Google Earth platform and images. The defined landslide surfaces were converted to kmz files using the Google Earth platform, and then, they were translated into shape files using the ArcGIS platform. The plotting and analysis were also conducted using the ArcGIS platform. The other basic maps used in this study include a 1 : 50,000 topographic map, its digital terrain DEM, and a 1 : 200,000 Leibo geologic map.

3.2. Acquisition of Geochronologic Data

The geochronologic data were obtained from both historical records of earthquakes and geologic dating. The geologic dating included two methods: the direct method involves the direct analysis of materials inside the fault zone, including fault gouges and fault breccia, while the indirect method involves the analysis of earthquake artifacts, e.g., the liquefaction of sandy soil and the presence of landslides and barrier lakes. The dating techniques used in this study include thermoluminescence (TL), scanning electron microscopy (SEM), electron spin resonance (ESR), radiocarbon dating (14C), optically stimulated luminescence (OSL), and 36Cl cosmogenic nuclide exposure methods. Sixteen of the geochronologic dates are from previous studies, and 23 new dates were obtained. Of these 23 new dates, 10 were obtained using the 36Cl method, 10 were obtained using the ESR method, and 3 were obtained using the OSL method. Several field photos are shown in Figure 3.

The 36Cl samples were tested at the accelerator mass spectrometry (AMS) facility in the Ion Beam Physics Laboratory of Swiss Federal Institute of Technology Zurich (ETH Zurich). The pretreatment and preparation were performed using the isotope dilution method [28]. The sample processing and analysis were conducted using the methods of Ivy-Ochs et al. [24, 29]; and the dating calculations were the same as those described by Ivy-Ochs et al. [24, 29] and Claude et al. [30].

The ESR samples were analyzed at the Division of Radiation Metrology, China Institute of Atomic Energy, using an EMX model spectrometer made by the German BRUKER company. The normal quantifier was an alanine thin film dosimeter with a dose rate of 20 Gy·min−1.

The OSL samples were analyzed at the Groundwater, Mineral water, and Environment Monitoring Center of China. The treated samples were analyzed using a Daybreak 2200 OSL instrument made in the United States. The dose rate of the 90Sr-Y beta radiation source was approximately 0.086273 Gy/sec. All of the samples were regenerated using simple multislice regeneration of the fine particles to obtain equivalent dose values and were fitted using the saturation index method.

4. Results and Analysis

4.1. Geometry of the Yaziba Fault

The available data indicate that the Yaziba Fault strikes approximately NS, dips 70°–80° west, and extend for 30 km (Figure 4). It is a regional reverse fault with a 200–400 m displacement that primarily cuts the Triassic strata, but at its southern end, it extends down into the Sinian strata. From south to north, this fault stretches along the Changping River (a tributary of the Jinsha River) that crosses the Jinsha River and extends along a ridge to the mountaintop on the left bank, turning at the top. Then, it extends along a gully, arriving back at the Jinsha River, runs along the left bank for 10 km, and crosses the Jinsha River for the second time before entering the mountain on the right bank (Figure 4). According to our site survey, the fault does follow this path. However, we found that as claimed by previous studies, at its northern end, the fault continues to extend northward for approximately 20 km (Figure 4). The main evidence for this is the discovery of breccia, fault gouge, and springs in this area (Figure 5). Figure 5(d) shows the river valley landforms of the region where the fault crosses the Jinsha River. To illustrate its structure, a cross section was constructed (Figure 6). This cross section shows that the fault cuts the lower Permian strata and the Permian Emeishan basalt. It is a thrust fault that dips 70° to the west. In addition, the landform to the northwest of the observed fault is a gully, implying that the fault may extend even farther. We speculate that the fault may connect to other faults in the Yingjing-Mabian-Yanjin Fault zone, such as No. 5 Manao Fault (Figure 2). Thus, the total fault length will be more than 50 km. Based on the China Geological Survey’s [31] classification criteria for fault activity, if the length of an active fault is used as the classification index, when the length of the fault is adjusted from 30 km to 50 km, the activity of the Yaziba Fault increases by one level, from weak activity to medium activity.

4.2. Paleo Earthquakes around the Fault

As was previously stated, 39 dates were used in this study (Table 1). Of these dates, 8 samples were for fault gouges and fault breccia, 18 were for landslides, 10 were for liquefied sandy soils and Quaternary overburden beds, and 3 were for barrier lakes. Eight samples were dated using TL, 2 samples were dated using SEM, 13 samples were dated using ESR, 2 samples were dated using 14C, 4 samples were dated using OSL, and 10 samples were dated using 36Cl. To properly consider the errors in each dating method and their large differences in accuracy, the data were analyzed, screened, and merged before it was used to identify paleo earthquakes. As shown in Table 1, samples D1 (80.9 ka) and D5 (81 ka) are adjacent, which indicates that these samples may be from the same earthquake; while sample D3 is Q3, (126–11.7 ka), and thus an earthquake occurred before 81 ka. Samples B1 (113 ka) and B6 (115 ka) represent two landslides that may have been formed by the same earthquake, so a seismic event occurred before 113 ka. Samples S4 and S5 are both from lacustrine sediments near the fault, and the TL and OSL dating results are adjacent, so an earthquake occurred before 43.3–44.0 ka. Sample B16 (29.6 ka) is the average value of three samples with ages of 27.9 ka, 33.1 ka, and 27.7 ka. Sample MH-I (20.8 ka) represents the first event, and sample MH-II (13.2 ka) represents the second event, which formed the Mahu landslide. These two ages were obtained from 10 samples using the 36Cl cosmic nuclide exposure dating method, combined with geomorphological and lithological analysis of the Mahu landslide, so they are deemed reliable. Through these comparative analyses, relatively reliable ages were established for the paleo earthquakes that occurred around the fault since 147,000 years (Table 2). The locations of the dated samples presented in Table 2 are shown in Figure 4. The time intervals listed in Table 2 are plotted in Figure 7. As can be seen from the power function trend line in Figure 7, from 147,000 years to 3800 years ago, the time interval between earthquakes is becoming shorter.

In addition, historical records show that four major earthquakes occurred from 1216 to 1917 AD, i.e., the 1216 Mahu, 1844 Dabei, 1917 Jilipu, and 1974 Daguan earthquakes, which results in an average recurrence interval of 175 years. It should be noted that since 1844, the recurrence interval is only 43 years, indicating an enhancement trend of seismicity in this area.

4.3. Landslides on Both Sides of the Yaziba Fault

The research area is centered on the Yaziba Fault, extends for approximately 17 km to each side of the fault, and covers an area of 1660 km2. In the study area, 164 landslides were identified from field surveys and remote sensing interpretations (Figure 8). Among them, four of the landslides (groups), i.e., the Mahu, Guancun, Damaotan, and Laizigou slope failures, have volumes on the scale of billions of m3. The basic characteristics of the Mahu landslide have been described in [35]. Apparently, many landslides have occurred around the Yaziba Fault. Here, two commonly used statistical indices of landslide distribution, i.e., the landslide area percentage (LAP) and the landslide number density (LND) [1, 36], are used to analyze the relationship between these landslides and the Yaziba Fault:where is the landslides area in each zone, is the number of landslides in each zone, and is the area of each zone.

For the selected 33 km × 50 km rectangular study area, the is 0.07 and the is 0.1 km−2, which is higher than those for most mountainous areas. Taking the Yaziba Fault as the center line, the left and right sides can each be divided into five zones: 0–3 km, 3–6 km, 6–9 km, 9–12 km, and >12 km (Figure 8). Due to the large scale and complexity of the Mahu landslide [35], it was eliminated from the statistical analysis. The LAP and LND of each zone were calculated, and the results are shown in Figure 9. Figure 9(a) shows that on the left side of the fault, the LAP and LND gradually decrease from 0 to 6 km, increase from 6 to 9 km, and reach a peak value at >12 km, which is mainly caused by the Laizigou landslide group (Figure 8). Figure 9(b) shows that on the right side of the fault, the LAP and LND gradually decrease from 0 to 6 km, increase from 6 to 9 km, reach a peak value at 9–12 km, and decrease from >12 km. The peak value at 9–12 km is primarily caused by the landslides in area A in Figure 8. There could be a seismic fault in area A that caused the abundance of landslides. Generally speaking, the landslides are more developed in the 0–6 km zone on both sides of the Yaziba Fault. Based on relevant research on earthquake-triggered landslides [1, 37, 38], these landslides were most likely induced by paleo earthquakes caused by the Yaziba Fault.

One of the landslides previously mentioned, the Daomakan landslide (Figure 10), exhibits characteristics typical of earthquake-induced slope failure, and it is only 2.9 km from the Yaziba Fault. The topography around this landslide is similar to a round-backed armchair, with a current slope of 59° toward the Jinsha River. The upper part of the slope is the back wall of the landslide, the sliding surface is the plane of the limestone bedrock, and the middle of the slope (380 m·a.s.l.) is the landslide deposit accumulation platform. Since the leading edge of the landslide has been eroded by the Jinsha River, the slope in the front is relatively gentle, approximately 15°–20°. The top of the head scarp is 500 m·a.s.l., and the leading edge of the residual deposit is 335 m·a.s.l., 5 m above the water surface of the Jinsha River. The volume of the residual landslide deposit is approximately 150 × 104 m3.

The landslide deposit consists of boulders and granular limestone soils with a low clay content, a relatively high gravel content, and rare large boulders. The bedrock of the slide bed is exposed at the front of the landslide (Figure 11), which is a fractured structure. Although the bedding structure of the original limestone is preserved, the bedding is disordered, indicating that the area was strongly compressed by tectonic movement before the landslide. The limit equilibrium method was used to determine how much seismic force was needed to cause slope instability. The equation of Dai et al. [16] was used:where is the slice index, is the seismic acceleration (), and is the gravitational acceleration. Only the horizontal seismic acceleration was considered, while the vertical seismic acceleration was assumed as zero. is the cohesion on the slide surface; is the internal friction angle on the slide surface; is the bottom length; is the pore water pressure on the bottom surface; and is the gravity, is the mass; and is the angle between the tangent to the center of the base of each slice and the horizontal.

The computational model is shown in Figure 12, in which R is the radius of the circle considered in the slope stability analysis. Referring to the empirical values given by Shi-biao and Zhang [39], is 0, the equivalent of is 60°, and is ignored. The calculation indicates that a horizontal seismic acceleration of 0.6 g is needed to cause landslide instability. In the absence of seismic force, if the bank slope slides, its equivalent should be as low as 38.5°, which is impossible for the limestone rock mass. Therefore, the Daomakan landslide was most likely induced by an earthquake.

5. Discussion

Previous studies of this region have yielded the following results. Liu [40] suggested that possible moderate-large earthquakes in the Mabian region would most likely cause earthquake swarms in the future. Li and Cheng [41] assumed that this fault zone was entering its fourth active period, and M 6 or greater earthquakes may occur in the next few years. Wang and Li [42] pointed out that there may be hidden deep faults beneath the fault zone. Hou et al. [33] analyzed the tectonic stress field of the Mabian-Daguan fault zone and revealed its complicated geologic structure. They concluded that the Yaziba Fault is a Holocene active fault, while most of the other faults were active during the Late and Middle Pleistocene. Ruan et al. [43] concluded that the potential for strong earthquakes exists in the southern part of the Yaziba Fault. Cheng [44] estimated that an M 6.75 earthquake may occur in the southern section of the Mabian-Daguan Fault zone.

Based on our field investigations, the length of the Yaziba Fault is extended from its previously determined 30 km to 50 km. Our analysis of 39 dated events around the fault suggests that 19 paleo earthquakes have occurred since 147,000 years ago. Finally, 164 landslides on both sides of the fault were identified from field investigations and interpretations of remote sensing data. The LAP accounts for 7.0% of the study area, which indicates that landslide susceptibility is very high in this area. As an example, the Daomakan landslide was verified to be triggered by an ancient earthquake. Our results further indicate that the Yaziba Fault has strong activity, and the fault area also has strong seismic activity.

Many paleo earthquakes and several major temblors have been recorded throughout history and many landslides have occurred around the fault, likely caused by seismic ground shaking. All of this information suggests a high level of seismicity in this area. Thus, it is important that we determine if it is possible that a strong earthquake will occur on the Yaziba Fault in the near future. Here, we examine the nearby Longmenshan Fault zone, on which the 2008 Mw 7.9 Wenchuan occurred. The near-surface dip angle of the surface rupture produced by this shock was 70°–80°. This is the first event with such a feature that occurred in the interior of the continent in the historical record. Across the world, mostly M 8.0 or greater thrust earthquakes occur on large fault zones with slip rates of more than 20 mm/a. The slip rate of the Longmenshan Fault is less than 2 mm/a as evidenced by the geologic and GPS data. Zhang et al. [45] considered that the Yangtze quasi-platform, on which the Sichuan Basin is located, is rigid and is not easily deformed. Thus, it obstructs the horizontal southeastward extension of the Tibetan Plateau. As the boundary between the Tibetan Plateau and the Sichuan Basin, the Longmenshan Fault zone has a high frictional strength and a high angle listric thrust structure. Thus, it is also not easily deformed and is capable of accumulating a large amount of stress. When the accumulated stress exceeds the strength of the fault zone, the fault can rupture suddenly, resulting in large earthquakes such as the 2008 Wenchuan earthquake. The Yingjing-Mabian-Yanjin Fault zone and the Longmenshan Fault zone are both located in the transitional region between the Tibetan Plateau and the Yangtze quasi-platform and have similar thrust nappe structures. The Yaziba Fault in the Yingjing-Mabian-Yanjin Fault zone also has a large dip angle of 70°–80°. Furthermore, the slip rate of the Yingjing-Mabian-Yanjin Fault zone is very low, indicating very weak tectonic deformation [42, 46]. In light of the adjacent Longmenshan Fault zone, the small slip rate of the Yingjing-Mabian-Yanjin Fault zone does not mean that it is stable, but rather that it is likely slowly accumulating elastic strain. Therefore, we conclude that there is a high probability of strong earthquakes occurring in the near future, especially on the Yaziba Fault.

Since the northern end of the Yaziba Fault intersects the Jinsha River, once an earthquake occurs on the fault, many landslides may be induced on both sides of the river. These landslides will slide into the Jinsha River, resulting in secondary disasters such as surges and barrier lakes, which will seriously endanger the people in this region and the safety of the two giant hydropower stations. It should be noted that due to the construction of the Xiangjiaba hydropower station, a new migration village called Qingping County, which contains thousands of migrants, has been settled on the northern end of the Yaziba Fault (Figures 2 and 5(d)). Therefore, the risk assessment of earthquakes and geological hazards in the Qingping County area is necessary. Since the rocks in this region are limestone and basalt, their strength reduction is lower after they are saturated, so the influence of reservoir impoundment on slope stability is lower than that of an earthquake. Therefore, more attention should be paid to the study of the Yaziba Fault’s activity. The results of this study provide useful information for risk assessment and disaster prevention.

6. Conclusions

Through our research, four new findings about the Yaziba Fault and the fault area are put forward. (1) Our field investigations provide new evidence for the northward extension of the Yaziba Fault by 20 km along the Jinsha River in southwestern China. (2) Based on our analysis of new and previously published dates, we conclude that the time interval between earthquakes has shortened from 147,000 to 3,800 years ago. In addition, the historical earthquake records from 1216 to the present indicate a high seismicity around the fault in recent years. (3) Based on remote sensing interpretations, a large number of landslides have been identified in the fault area, with the LND reaching 0.1 km−2 and the LAP reaching 0.07, which are higher than those of most mountainous areas. In general, the LND and LAP gradually decrease with the increasing distance from the fault, indirectly reflecting the strong activity of the Yaziba Fault. (4) The limit equilibrium calculation shows that the Daomakan landslide required 0.6 g of horizontal seismic acceleration, which confirms the strong activity of the Yaziba Fault.

Paleo earthquake data are useful for long-term earthquake prediction and the study of fault activity. Based on the results of geologic dating methods and the analysis of ancient landslides associated with seismic ground shaking, it is possible to identify more paleo earthquakes around active faults. These results provide valuable information for the assessment of seismic hazards. This study investigated the area around the Yaziba Fault and identified 19 paleo earthquakes based on 39 geological dating of fault material and nearby landslide artifacts. Considering the fault’s scale, the amount of data obtained is still not enough. In the future, we should collect and analyze more data to identify additional paleo earthquakes in this area.

Data Availability

The chronological data and landslides used to support the findings of this study are all included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Key Research and Development Program of China (2018YFC1505004), the National Natural Science Foundation of China (41807267 and 41661144037), and the University Natural Science Research Project of Anhui Province (501100009558). We thank Fuchu Dai of Beijing University of Technology for his field investigation.