汶川地震触发平溪村滑坡特征及成因分析

平溪村滑坡是汶川地震触发的高速岩质滑坡,面积为3.7×104m2,体积约65×104m3。滑坡距离映秀—北川发震断裂带中的石坎断层不足500m,在强大地震力作用下,滑坡滑动面陡峭、粗糙,与重力作用下呈圆弧、平滑的滑面迥异。地震触发的滑坡可分为三个阶段:①震动拉裂阶段;②摩擦阻力降低、"锁固段"剪断阶段;③滑体溃滑,高速流动堆积阶段。然后结合结构面的力学分析,通过数值模拟对地震作用下的应力变化和震裂机制进行了初步分析。     

Preliminary investigation of some large landslides triggered by the 2008 Wenchuan earthquake, Sichuan Province, China

The M _s 8.0 Wenchuan earthquake or “Great Sichuan Earthquake” occurred at 14:28 p.m. local time on 12 May 2008 in Sichuan Province, China. Damage by earthquake-induced landslides was an important part of the total earthquake damage. This report presents preliminary observations on the Hongyan Resort slide located southwest of the main epicenter, shallow mountain surface failures in Xuankou village of Yingxiu Town, the Jiufengchun slide near Longmenshan Town, the Hongsong Hydro-power Station slide near Hongbai Town, the Xiaojiaqiao slide in Chaping Town, two landslides in Beichuan County-town which destroyed a large part of the town, and the Donghekou and Shibangou slides in Qingchuan County which formed the second biggest landslide lake formed in this earthquake. The influences of seismic, topographic, geologic, and hydro-geologic conditions are discussed.     

Statistical analysis of an earthquake-induced landslide distribution — the 1989 Loma Prieta, California event

The 1989 Loma Prieta, California earthquake (moment magnitude, M=6.9) generated landslides throughout an area of about 15,000 km² in central California. Most of these landslides occurred in an area of about 2000 km² in the mountainous terrain around the epicenter, where they were mapped during field investigations immediately following the earthquake. The distribution of these landslides is investigated statistically, using regression and one-way analysis of variance (ANOVA) techniques to determine how the occurrence of landslides correlates with distance from the earthquake source, slope steepness, and rock type. The landslide concentration (defined as the number of landslide sources per unit area) has a strong inverse correlation with distance from the earthquake source and a strong positive correlation with slope steepness. The landslide concentration differs substantially among the various geologic units in the area. The differences correlate to some degree with differences in lithology and degree of induration, but this correlation is less clear, suggesting a more complex relationship between landslide occurrence and rock properties.     

Interpretation of landslide distribution triggered by the 2005 Northern Pakistan earthquake using SPOT 5 imagery

The 2005 northern Pakistan earthquake (magnitude 7.6) of 8 October 2005 occurred in the northwestern part of the Himalayas. We interpreted landslides triggered by the earthquake using black-and-white 2.5-m-resolution System Pour l’Observation de la Terre 5 (SPOT 5) stereo images. As a result, the counts of 2,424 landslides were identified in the study area of 55 by 51 km. About 79% or 1,925 of the landslides were small (less than 0.5 ha in area), whereas 207 of the landslides (about 9%) were large (1 ha and more in area). Judging from our field survey, most of the small landslides are shallow rock falls and slides. However, the resolution and whitish image in the photos prevented interpreting the movement type and geomorphologic features of the landslide sites in detail. It is known that this earthquake took place along preexisting active reverse faults. The landslide distribution was mapped and superimposed on the crustal deformation detected by the environmental satellite/synthetic aperture radar (SAR) data, active faults map, geological map, and shuttle radar topography mission data. The landslide distribution showed the following characteristics: (1) Most of the landslides occurred on the hanging-wall side of the Balakot–Garhi fault; (2) greater than one third of the landslides occurred within 1 km from the active fault; (3) the greatest number of landslides (1,147 counts), landslide density (3.2 counts/km²), and landslide area ratio (2.3 ha/km²) was found within Miocene sandstone and siltstone, Precambrian schist and quartzite, and Eocene and Paleocene limestone and shale, respectively; (4) there was a slight trend that large landslides occurred on vertically convex slopes rather than on concave slopes; furthermore, large landslides occurred on steeper (30° and more) slopes than on gentler slopes; (5) many large landslides occurred on slopes facing S and SW directions, which is consistent with SAR-detected horizontal dominant direction of crustal deformation on the hanging wall.     

Geological and geomorphological characteristics of landslides triggered by the 2004 Mid Niigta prefecture earthquake in Japan

The 2004 Mid Niigta prefecture earthquake (M_JMA 6.8) triggered more than one thousand landslides in the Miocene to Quaternary sedimentary rocks in Japan. The most common landslides were shallow disrupted landslides on steep slopes, which has been common in many previous disastrous earthquakes in the world. The Mid Niigta prefecture earthquake also triggered more than one hundred deep landslides, providing valuable information on the conditions for their occurrence. A field investigation and the interpretation of aerial photographs taken before and after the earthquake suggest that reactivation of existing landslides and undercutting of slopes are the most important factors for deep landslides to be triggered by earthquakes. In addition, planar sliding surfaces seem to be essential for the generation of catastrophic landslides triggered by this earthquake. Planar bedding–parallel sliding surfaces were formed at the boundary between the overlying permeable sandstone and underlying siltstone or along the bedding planes of alternating beds of sandstone and siltstone. Sliding surfaces along the slope-parallel oxidation front were formed in the area of black mudstone. New landslides (rockslide-avalanches) occurred with the sliding surfaces in a several-cm thick tuff interbedded in siltstone. One rockslide-avalanche occurred on a slope where buckling deformation preceded the earthquake. Gentle valley bottom sediments were mobilized in many locations, probably because they were saturated and partial liquefaction had occurred due to the earthquake shaking.     

A comparison of the GIS based landslide susceptibility assessment methods: multivariate versus bivariate

The purpose of this study is to evaluate and to compare the results of multivariate (logical regression) and bivariate (landslide susceptibility) methods in Geographical Information System (GIS) based landslide susceptibility assessment procedures. In order to achieve this goal the Asarsuyu catchment in NW Turkey was selected as a test zone because of its well-known landslide occurrences interfering with the E-5 highway mountain pass.Two methods were applied to the test zone and two separate susceptibility maps were produced. Following this a two-fold comparison scheme was implemented. Both methods were compared by the Seed Cell Area Indexes (SCAI) and by the spatial locations of the resultant susceptibility pixels.It was found that both of the methods converge in 80% of the area; however, the weighting algorithm in the bivariate technique (landslide susceptibility method) had some severe deficiencies, as the resultant hazard classes in overweighed areas did not converge with the factual landslide inventory map. The result of the multivariate technique (logical regression) was more sensitive to the different local features of the test zone and it resulted in more accurate and homogeneous susceptibility maps.     

Geological causes and geomorphological precursors of the Tsaoling landslide triggered by the 1999 Chi-Chi earthquake, Taiwan

The Tsaoling landslide, one of the largest landslide areas in Taiwan, has been affected by catastrophic events triggered by rain or earthquakes six times since 1862. These landslides, including that caused by the 1999 earthquake, have essentially not been reactivated old slides, but were sequential new ones that developed upslope, retrogressively. The landslide area is underlain by Pliocene sandstone and shale to form a dip slope with a bedding plane, dipping uniformly at 14°. The slip surface of the 1999 landslide was smooth and planar, parallel to the bedding plane with a slightly stepped profile; it formed within thinly alternated beds of fine sandstone and shale with ripple lamination or in a shale bed. The shale is weathered by slaking and probably by sulfuric acid, which is inferred to be one of the major causes of the intermittent retrogressive development of the landslides. The weathering was likely accelerated by the removal of overlying beds during earlier landslides in 1941 and 1942. The top margin of the 1999 landslide, in plan view, coincided with a V-shaped scarplet, which can be clearly recognized on aerial photographs taken before the landslide. This geomorphological feature indicates that this landslide had already moved slightly before its 1999 occurrence, providing precursory evidences.     

走滑断裂型地震诱发的滑坡在断裂两盘的空间分布差异

以玉树地震滑坡为实例,选择高程、坡度、坡向、坡位、水系、地层岩性、同震地表破裂、地震动峰值加速度(PGA)8个因子,以地震滑坡面积百分比(LAP)与滑坡点密度(LND)为指标,研究走滑断裂型地震诱发滑坡在断裂两盘的空间分布差异。在分析这些影响因子的断裂两盘的差异的基础上,基于LAP与LND两个指标详细分析断裂两盘每个因子内部级别滑坡发育的情况。结果表明,总体上北盘的滑坡较南盘发育,除了个别因子级别内南北盘滑坡差别较大外,大部分表现为南北盘滑坡分布情况类似。总之,玉树地震滑坡在断裂两盘的空间分布基本类似,只是在某些因子分级内存在一定的差异。     

走滑断裂型地震诱发的滑坡在断裂两盘的空间分布差异

以玉树地震滑坡为实例,选择高程、坡度、坡向、坡位、水系、地层岩性、同震地表破裂、地震动峰值加速度(PGA)8个因子,以地震滑坡面积百分比(LAP)与滑坡点密度(LND)为指标,研究走滑断裂型地震诱发滑坡在断裂两盘的空间分布差异。在分析这些影响因子的断裂两盘的差异的基础上,基于LAP与LND两个指标详细分析断裂两盘每个因子内部级别滑坡发育的情况。结果表明,总体上北盘的滑坡较南盘发育,除了个别因子级别内南北盘滑坡差别较大外,大部分表现为南北盘滑坡分布情况类似。总之,玉树地震滑坡在断裂两盘的空间分布基本类似,只是在某些因子分级内存在一定的差异。     

抚顺西露天矿区滑坡易发性评价与时空特征分析

近年抚顺西露天矿矿区滑坡在规模、活动频率以及空间发育位置等方面均与矿区历史滑坡呈现出较明显的时空差异性,为矿区未来治理提出新的挑战。故基于改进频率比模型,选择人类工程活动（主要是采矿工程扰动）和工程地质条件两类致灾因素,对由矿区滑坡时空差异所导致的易发性变化进行研究。结果表明：影响矿区滑坡灾害的主要因素包括坡体结构、工程岩组和水文地质条件;矿区滑坡易发性较2010年前有所降低,高易发区和较高易发区面积共减少1.294 km2,但北帮中区等地段仍呈高易发态势;滑坡易发区的空间位置与2010年前后对比发生了较大变化,矿坑西北帮滑坡易发性大幅降低,东部环矿区境界50~500 m处滑坡易发性轻微降低,北帮中段与南帮中、西段顶部滑坡易发性增加,其中北帮中段为目前最易发的区段。研究结果可为矿区闭坑停采后边坡整治及后续规划利用提供参考依据。     

基于XGBoost模型的三峡库区燕山乡滑坡易发性评价与区划

滑坡易发性评价是精细化滑坡灾害风险评价的基础。为了提升滑坡易发性评价模型的精度和稳健性,以三峡库区万州区燕山乡为例,选取工程地质岩组、堆积层厚度等九个影响因子构建滑坡易发性评价指标体系,应用信息量模型定量分析滑坡发育与指标之间的关系。在此基础上,随机选取70%/30%的滑坡样本作为训练/验证数据集,应用极致梯度提升模型(extreme gradient boosting, XGBoost)开展易发性评价。随后从模型预测精度和模型稳定性两方面将其与决策树模型(decision tree, DT)和梯度提升树模型(gradient boosting decision tree, GBDT)进行对比。结果表明：研究区堆积层滑坡主要受长江水系、堆积层厚度和工程地质岩组影响。XGBoost模型具有最高的准确率(94.3%)和预测精度(97.3%)。在模型稳定性验证中,平均预测精度最高(97.3%),优于DT(91.3%)和GBDT(95.7%),模型标准差和变异系数均为0.01,低于其余两种模型。XGBoost在区域滑坡易发性评价与制图中得到了可靠的结果,为滑坡灾害空间预测提供了新的技术支撑。     

基于确定性系数方法的龙门山地区地震滑坡影响因子敏感性分析——以汶川地震与芦山地震诱发滑坡为例

在我国西南地区,沿龙门山断裂带分别发生了2008年汶川M_S8.0级地震以及2013年芦山M_S7.0级地震,这两次地震均造成了严重的地表破裂,并诱发了大量的滑坡和崩塌等次生地质灾害。文章选择了位于青藏高原向四川盆地过渡的区域——以龙门山断裂带为中心的30°~34°N,102°~106°E区域作为研究区,借助GIS工具,基于确定性系数（CF）方法,选取了地震、地质构造、自然环境和人类活动4大类因子,包括烈度、震中距、岩性、断裂、高程、坡度、坡向、河流、降雨、公路共10个因子（子集）,对汶川和芦山地震诱发的次生滑坡灾害进行影响因子敏感性分析,基于z值确定该区域内地震滑坡的关键因子类以及基于CF值确定各类集（子集）下各特定因子的具体地震滑坡敏感性。研究结果显示：地震因子具有最高的z值,是龙门山地区地震滑坡产生的关键影响因子,表明地震活动的强弱直接关系到斜坡的稳定性和次生滑坡灾害的面积分布。而对比烈度子集中的具体CF值表明：当烈度小于Ⅷ度时,烈度对次生滑坡发生的影响极低,区域内的主要影响因子则由地震因子转变为震中距、自然因子等其他类别的因子;其次以坡度、高程、与河流的距离为主的自然因子类别以及与断层的距离在地震滑坡过程中也有较高权重,而人类活动对研究区内坡体的稳定性也有着不可忽略的作用。本研究结果可作为该地区后续区域地震滑坡相关研究和发展规划的基础科学依据。

     

Geometry of the Wagner basin,upper Gulf of California based on seismic reflections

The Wagner basin occupies the northernmost spreading centre in the Gulf of California, located along the Pacific‐North America plate boundary. It is filled with sediments from the Colorado River that obscure its bathymetric expression; therefore it is not as well defined as other basins in the central and southern Gulf of California. To define the geometry and extension of the Wagner basin, a 2D multi‐channel seismic reflection database was used. Data were collected by Petroleos Mexicanos (PEMEX) in 1979–1980. The most important regional structural features identified are the Consag and Wagner normal faults and the Cerro Prieto strike‐slip fault. These structures play an important role in the development of the basin. The Consag fault, described for the first time in this paper, marks the western side of the basin. The eastern and northwest limits are bound by the Cerro Prieto and Wagner faults respectively. The Wagner fault intersects the Cerro Prieto fault at an angle of 130°, bending the depocentre in a NW direction, adjacent to the Cerro Prieto fault zone. The northernmost segment of the Consag fault bends 25° in a NE direction and joins the Cerro Prieto fault at an angle of 110°. Greater subsidence (up to 300 m) takes place along the northern trace of the Cerro Prieto fault, with a downthrown displacement of 400 m. The Consag and Wagner breaks obliquely intersect the Cerro Prieto fault, and, inasmuch as both are normal faults, they have small horizontal slip components which generated oblique displacement. This structural pattern is different relative to the pattern of basins located south of Wagner basin, such as the Upper and Lower Delfin basins. The orientations of the normal faults are perpendicular to the master fault (Ballenas transform fault). The relationship between normal and transform faults in the Wagner basin and the observed ‘S’ shape are typical of a basin that has not yet reached maturity. As a result of this study, the previously uncertain area (～1330 km²) and perimeter (158 km) of the Wagner basin were defined.     

基于自组织特征映射网络-随机森林模型的滑坡易发性评价

为深入探讨评价单元和非滑坡样本选取对滑坡易发性预测的影响,构建了一种基于自组织特征映射网络-随机森林模型的滑坡易发性评价模型。该模型针对栅格单元和斜坡单元在滑坡易发性评价中的不足,结合栅格单元和斜坡单元的相互关系,提出了滑坡易发性指数的优化计算方法。在此基础上,基于随机森林Tree Bagger分类器构建滑坡易发性评价模型,通过对比分析自组织特征映射网络和随机方法选取非滑坡样本对评价结果的影响,探讨自组织特征映射网络、随机森林和自组织特征映射网络-随机森林三种评价模型的有效性;将评价模型应用于大余县滑坡易发性评价。结果显示,随机森林模型和自组织特征映射网络-随机森林模型的预测精度较高,分别达到91.19%和94.94%,成功率曲线的AUC值分别为0.822和0.849,表明自组织特征映射网络-随机森林模型具有更高的预测率和成功率, 自组织特征映射网络聚类的预测精度虽然有限,但作为非滑坡样本的选择方法,能够有效提高随机森林模型的评价精度。     

Andean earthquakes triggered by the 2010 Maule,Chile (Mw 8.8) earthquake: Comparisons of geodetic,seismic and geologic constraints

The Maule, Chile, (M_w 8.8) earthquake on 27 February 2010 triggered deformation events over a broad area, allowing investigation of stress redistribution within the upper crust following a mega-thrust subduction event. We explore the role that the Maule earthquake may have played in triggering shallow earthquakes in northwestern Argentina and Chile. We investigate observed ground deformation associated with the M_w 6.2 (GCMT) Salta (1450 km from the Maule hypocenter, 9 h after the Maule earthquake), M_w 5.8 Catamarca (1400 km; nine days), M_w 5.1 Mendoza (350 km; between one to five days) earthquakes, as well as eight additional earthquakes without an observed geodetic signal. We use seismic and Interferometric Synthetic Aperture Radar (InSAR) observations to characterize earthquake location, magnitude and focal mechanism, and characterize how the non-stationary, spatially correlated noise present in the geodetic imagery affects the accuracy of our parameter estimates. The focal mechanisms for the far-field Salta and Catamarca earthquakes are broadly consistent with regional late Cenozoic fault kinematics. We infer that dynamic stresses due to the passage of seismic waves associated with the Maule earthquake likely brought the Salta and Catamarca regions closer to failure but that the involved faults may have already been at a relatively advanced stage of their seismic cycle. The near-field Mendoza earthquake geometry is consistent with triggering related to positive static Coulomb stress changes due to the Maule earthquake but is also aligned with the South America-Nazca shortening direction. None of the earthquakes considered in this study require that the Maule earthquake reactivated faults in a sense that is inconsistent with their long-term behavior.     

基于随机森林-特征递归消除模型的可解释性缓丘岭谷地貌滑坡易发性评价

研究旨在基于随机森林-特征递归消除模型,通过SHAP算法(SHapley Additive exPlanation, SHAP)与部分依赖图(Partial Dependence Plot, PDP)对缓丘岭谷地貌区域进行滑坡易发性评价与内部机制解释,以期为地质灾害防治研究提供参考。利用优化随机森林算法对典型缓丘岭谷地区滑坡易发性进行研究,建立缓丘岭谷滑坡易发性评价模型;利用特征递归消除算法剔除噪声因子,选取地形地貌、地质构造、环境条件、人类活动5个类型16个因子构建重庆合川区滑坡致灾因子数据库;结合合川区754个历史滑坡点,利用随机森林算法对因子重要性进行排序,并根据专家经验法对研究区的滑坡易发性进行划分,将研究区的滑坡易发性分为极低、低、中、高、极高5个等级;应用部分依赖图对合川区滑坡发生影响大的因子进行解释和SHAP算法对个体滑坡进行局部解释。结果表明：与原模型相比,随机森林-特征递归消除模型测试集AUC值提高了0.019,证明了特征递归消除算法的有效性;训练集以及测试集的AUC值分别为0.769、0.755,具有较高的预测精度;缓丘缓坡地区在起伏较大地区滑坡密度较大,历史滑坡多集中于高易发地区;滑坡的空间分布具有不均匀性与复杂性,各致灾因子对滑坡发生的影响有着明显的区域特征与空间异质性,在缓坡丘陵地区多年平均降雨、高程、岩性3个因子对滑坡发生的影响最大;由SHAP算法对合川白塔坪上山公路滑坡事件进行解释,岩性与高程对滑坡起抑制作用,起伏度、坡度、归一化植被指数(NDVI)与POI核密度促进滑坡发生。综上所述,基于随机森林-特征递归消除模型在缓丘岭谷区滑坡易发性评价中具有较高的准确性,通过部分依赖图与SHAP算法对全局滑坡与个体滑坡发生的内在机理进行解释分析,有利于构建与完善不同地貌环境下滑坡易发性评价因子体系并探究滑坡内部决策机理,可为区域滑坡易发性评估与地质灾害防治提供参考。     

基于流固耦合的强震大型滑坡水力激发效应研究

2008年汶川Ms8.0级强震触发了体积近12×108 m3的大光包滑坡。该滑坡发生于古生代碳酸盐岩地层,滑带地质背景为斜坡内部深埋400 m、最大厚度达5 m的先期层间构造错动带。最新调查表明,该错动带是斜坡内部地下水通道,错动带岩体处于饱和状态。为揭示强震过程与地下水相关的大光包滑坡启动机制,提出了一种具有软弱层带的硬质碳酸盐岩边坡简化模型,将层间构造错动带概化为碳酸盐岩硬层内部软弱层带,采用FLAC3D程序中的流固耦合算法模拟了模型的响应特性。研究结果表明:强震过程中软弱层带上下碳酸盐岩硬层的变形响应时间、波型、大小出现明显差异,上硬层相对于下硬层产生了张离、压缩和剪切3种非协调变形模式,由此对软弱层带产生了振动冲压-张拉和振动剪切动力学行为,饱水软弱层带形成了具有瞬间放大和累积增涨特征的超孔隙水压力。这里将上下硬层差异性变形称为非协调变形,认为非协调变形是软弱层带应力放大成因,推测软弱层带应力瞬间放大以及放大应力长持时作用下的岩体致损是超孔隙水压力激发和累积的成因;强震过程软弱层带超孔隙水压力导致其内有效应力快速降低,使得斜坡前部锁固段应力快速集中,而后被突然剪断,滑坡骤然启动,揭示了强震过程中超孔隙水压力是大光包滑坡启动的主要原因。     

Interpretation of earthquake-induced landslides triggered by the 12 May 2008, M7.9 Wenchuan earthquake in the Beichuan area, Sichuan Province, China using satellite imagery and Google Earth

The 12 May 2008 M7.9 Wenchuan earthquake in the People’s Republic of China represented a unique opportunity for the international community to use commonly available GIS (Geographic Information System) tools, like Google Earth (GE), to rapidly evaluate and assess landslide hazards triggered by the destructive earthquake and its aftershocks. In order to map earthquake-triggered landslides, we provide details on the applicability and limitations of publicly available 3-day-post- and pre-earthquake imagery provided by GE from the FORMOSAT-2 (formerly ROCSAT-2; Republic of China Satellite 2). We interpreted landslides on the 8-m-resolution FORMOSAT-2 image by GE; as a result, 257 large landslides were mapped with the highest concentration along the Beichuan fault. An estimated density of 0.3 landslides/km² represents a minimum bound on density given the resolution of available imagery; higher resolution data would have identified more landslides. This is a preliminary study, and further study is needed to understand the landslide characteristics in detail. Although it is best to obtain landslide locations and measurements from satellite imagery having high resolution, it was found that GE is an effective and rapid reconnaissance tool.     

基于危险度区划的县级区域降雨引发滑坡的风险预警方法——以四川省米易县降雨滑坡为例

降雨滑坡预警在世界上开展广泛,但对小区域的降雨滑坡风险进行预警尚不多见。以通用的风险定义为基础,首先对目标区域——四川省米易县的地质灾害进行调查,以典型滑坡进行反演,获得灾害的本底因素,按贡献率权重叠加法进行危险度分区;通过对承灾体的调查确定承灾对象,并将各承灾对象按密度由行政界限向1km×1km的网格单元转化后叠加形成易损度分区。二者相乘完成风险分区。对研究区域近5年逐日降雨数据与172条地质灾害记录进行分析,按前期日降雨量模型得到降雨阈值,根据滑坡空间概率和降雨引发滑坡的时间概率叠加的结果,得到研究区域的降雨滑坡概率,对照风险分区,完成风险预警。