尼泊尔M8.1地震震前卫星重力场时变特征

利用多年GRACE卫星重力变化资料,给出尼泊尔及周缘地区的重力场时空变化,以探索2015-04-25尼泊尔M8.1地震震前重力变化特征;并基于球形位错理论,利用断层模型数据模拟同震重力场变化。研究表明,区域重力场变化以2008年为界存在明显的分时段特征,2003~2008年以重力增加为主,而2009~2014年以重力减少为主;且存在明显的分区性,区域重力场在尼泊尔地震震中四周呈正负相间分布;震前局部重力场变化在压缩区重力增加,拉张区重力减少,与理论同震重力变化结果一致。总结研究区域近年来M7.0以上地震震前重力变化特征发现,地震基本发生在重力变化高梯度带的过渡带及重力变化零值线附近。     

Development of community-based landslide early warning system in the earthquake-affected areas of Nepal Himalaya

In the central Nepal Himalaya, landslides form the major natural hazards annually resulting in many casualties and damage. Structural as well as non-structural measures are in place to minimize the risk of landslide hazard. To reduce the landslide risk,a Landslide Early Warning System(LEWS) as a nonstructural measure has been piloted at Sundrawati village(Kalinchowk rural municipality, Dolakha district) to identify its effectiveness. Intensive discussions with stakeholders, aided by landslide susceptibility map, resulted in a better understanding of surface dynamics and the relationship between rainfall and surface movement. This led to the development of a LEWS comprised of extensometers,soil moisture sensors, rain gauge stations, and solar panels as an energy source that blows siren receiving signals via a micro-controller and interfacing circuit.The data generated through the system is transmitted via a Global System for Mobile Communications(GSM) network to responsible organizations in realtime to circulate the warning to local residents. This LEWS is user-friendly and can be easily operated by a community. The successful pilot early warning system has saved 495 people from 117 households in August 2018.However,landslide monitoring and dissemination of warning information remains a complex process where technical and communications skill should work closely together.     

长期静水浸泡对三峡库区滑带土物理-化学-力学性质的影响

三峡水库蓄水使得库岸大量土体长期处于浸泡状态,导致土体软化从而诱发滑坡失稳。为研究长期浸泡对滑坡土体物理-化学-力学性质的影响,以马家沟滑坡原状滑带土为对象开展了浸泡软化试验,通过比较不同浸泡时间滑带土的粒度分布、界限含水率、化学与矿物成分、剪切特性等特征,探讨了滑带土浸泡软化机理。研究结果表明:浸泡过程中滑带土中Ca2+、Mg2+等离子流失较多,但矿物成分无变化;浸泡后滑带土出现阶段性粒度细化现象,液塑限和塑性指数均随黏粒含量增加而增大;随着浸泡时间增加,滑带土应力应变关系在低法向应力下由强软化型变为弱软化型,在高法向应力下由软化型变为硬化型;滑带土抗剪强度参数随着浸泡时间增加呈指数形式降低,黏聚力c降低程度大于内摩擦角φ。研究成果可以为水库滑坡稳定性评价提供理论依据。      

Landslide susceptibility assessment of the region affected by the 25 April 2015 Gorkha earthquake of Nepal

Nepal was hit by a 7.8 magnitude earthquake on 25^th April, 2015. The main shock and many large aftershocks generated a large number of coseismic landslips in central Nepal. We have developed a landslide susceptibility map of the affected region based on the coseismic landslides collected from remotely sensed data and fieldwork, using bivariate statistical model with different landslide causative factors. From the investigation, it is observed that most of the coseismic landslides are independent of previous landslides. Out of 3,716 mapped landslides, we used 80% of them to develop a susceptibility map and the remaining 20% were taken for validating the model. A total of 11 different landslide-influencing parameters were considered. These include slope gradient, slope aspect, plan curvature, elevation, relative relief, Peak Ground Acceleration (PGA), distance from epicenters of the mainshock and major aftershocks, lithology, distance of the landslide from the fault, fold, and drainage line. The success rate of 87.66% and the prediction rate of 86.87% indicate that the model is in good agreement between the developed susceptibility map and the existing landslides data. PGA, lithology, slope angle and elevation have played a major role in triggering the coseismic mass movements. This susceptibility map can be used for relocating the people in the affected regions as well as for future land development.     

Landslide initiation and runout susceptibility modeling in the context of hill cutting and rapid urbanization: a combined approach of weights of evidence and spatial multi-criteria

Rainfall induced landslides are a common threat to the communities living on dangerous hill-slopes in Chittagong Metropolitan Area, Bangladesh. Extreme population pressure, indiscriminate hill cutting, increased precipitation events due to global warming and associated unplanned urbanization in the hills are exaggerating landslide events. The aim of this article is to prepare a scientifically accurate landslide susceptibility map by combining landslide initiation and runout maps. Land cover, slope, soil permeability, surface geology, precipitation, aspect, and distance to hill cut, road cut, drainage and stream network factor maps were selected by conditional independence test. The locations of 56 landslides were collected by field surveying. A weight of evidence (WoE) method was applied to calculate the positive (presence of landslides) and negative (absence of landslides) factor weights. A combination of analytical hierarchical process (AHP) and fuzzy membership standardization (weighs from 0 to 1) was applied for performing a spatial multi-criteria evaluation. Expert opinion guided the decision rule for AHP. The Flow-R tool that allows modeling landslide runout from the initiation sources was applied. The flow direction was calculated using the modified Holmgren’s algorithm. The AHP landslide initiation and runout susceptibility maps were used to prepare a combined landslide susceptibility map. The relative operating characteristic curve was used for model validation purpose. The accuracy of WoE, AHP, and combined susceptibility map was calculated 96%, 97%, and 98%, respectively.     

Landslide Susceptibility Mapping along Bhalubang–Shiwapur Area of Mid-Western Nepal Using Frequency Ratio and Conditional Probability Models

Roads constructed in fragile Siwaliks are prone to large number of instabilities. Bhalubang–Shiwapur section of Mahendra Highway lying in Western Nepal is one of them. To understand the landslide causative factor and to predict future occurrence of the landslides, landslide susceptibility mapping(LSM) of this region was carried out using frequency ratio(FR) and weights-of-evidence(W of E) models. These models are easy to apply and give good results. For this, landslide inventory map of the area was prepared based on the aerial photo interpretation, from previously published/unpublished reposts, and detailed field survey using GPS. About 332 landslides were identified and mapped, among which 226(70%) were randomly selected for model training and the remaining 106(30%) were used for validation purpose. A spatial database was constructed from topographic, geological, and land cover maps. The reclassified maps based on the weight values of frequency ratio and weights-of-evidence were applied to get final susceptibility maps. The resultant landslide susceptibility maps were verified andcompared with the training data, as well as with the validation data. From the analysis, it is seen that both the models were equally capable of predicting landslide susceptibility of the region(W of E model(success rate = 83.39%, prediction rate = 79.59%); FR model(success rate = 83.31%, prediction rate = 78.58%)). In addition, it was observed that the distance from highway and lithology, followed by distance from drainage, slope curvature, and slope gradient played major role in the formation of landsides. The landslide susceptibility maps thus produced can serve as basic tools for planners and engineers to carry out further development works in this landslide prone area.     

基于信息量模型和数据标准化的滑坡易发性评价

本文以北川曲山-擂鼓片区为研究区,将坡度、坡向、高程、地层、距断层的距离、距水系的距离和距道路的距离作为该区域滑坡易发性评价因子。采用信息量模型计算了各项评价因子的信息量值,并运用4种标准化模型对信息量值进行标准化处理。各评价因子的权重由层次分析法（AHP）确定。在GIS中将权重值和各评价因子的标准化信息量值,进行叠加计算得到区域滑坡总信息量值,并基于自然断点法对其进行重分类,将研究区划分为极高易发区、高易发区、中易发区、低易发区和极低易发区5级易发区。将基于4种标准化模型和信息量模型得到的滑坡易发性评价结果进行了对比分析,结果表明：基于最值标准化信息量模型的滑坡易发性评价结果的ROC曲线下面积AUC值为0.807,高于其余模型的AUC值,说明最值标准化信息量模型的滑坡易发性评价效果最好。极高易发区面积占研究区面积的20.03%,离断层和水系较近,主要分布地层为寒武系、志留系和三迭系。研究结果可为区内滑坡风险评价和灾害防治提供参考。     

基于模型试验的动水驱动型顺层岩质滑坡启滑机制初探

动水驱动型顺层岩质滑坡数量多、灾害频发、危害大, 是滑坡地质灾害领域的研究重点, 但目前对于滑坡启滑机制的认识仍不充分, 滑坡的准确预报还面临巨大挑战。鉴于此, 以含软弱夹层的中倾角顺层岩质滑坡为研究对象, 通过构建理想的单层滑带滑坡物理模型, 开展了一系列动水作用下的滑坡模型试验研究。结果表明, 动水作用下顺层岩质滑坡从开始变形至失稳滑动需经历初始变形、缓慢变形、加速变形和失稳破坏4个阶段, 而各个阶段的演化特征与滑面粗糙度和倾角密切相关。滑面倾角越大或粗糙度越小, 滑坡体从开始变形至失稳滑动所需的时间则越短; 相应地, 坡体加速变形阶段越不明显, 滑坡破坏的突发性越强。滑带内的渗流冲蚀作用会使滑带土中的骨料流失, 导致其抗剪强度降低, 进而引发坡体滑动。与此同时, 上覆坡体的压剪作用以及变形演化过程亦将反过来影响冲蚀强度。基于滑带土黏聚力随水力梯度和冲蚀时间的变化关系, 提出了渗流驱动下滑带土黏聚力演化模型, 可较好地描述滑带土黏聚力的退化过程。滑面粗糙度的存在不仅显著影响了滑带的冲蚀劣化规律, 还改变了滑带不同区域的破坏模式。此外, 通过考虑滑面粗糙度对滑带不同区域破坏模式的影响, 开展了动水多效应关联分析, 建立了滑坡地质体力学分析模型, 实现了动水作用下顺层岩质滑坡动态稳定性的有效评估。本研究成果可为实际动水驱动型顺层岩质滑坡的预测和防治提供理论参考。      

Effect of rainfall on a colluvial landslide in a debris flow valley

A colluvial landslide in a debris flow valley is a typical phenomena and is easily influenced by rainfall. The direct destructiveness of this kind of landslide is small, however, if failure occurs the resulting blocking of the channel may lead to a series of magnified secondary hazards. For this reason it is important to investigate the potential response of this type of landslide to rainfall. In the present paper, the Goulingping landslide, one of the colluvial landslides in the Goulingping valley in the middle of the Bailong River catchment in Gansu Province, China, was chosen for the study. Electrical Resistivity Tomography (ERT), Terrestrial Laser Scanning (TLS), together with traditional monitoring methods, were used to monitor changes in water content and the deformation of the landslide caused by rainfall. ERT was used to detect changes in soil water content induced by rainfall. The most significant findings were as follows:(1) the water content in the centralupper part (0~41 m) of the landslide was greater than in the central-front part (41~84 m) and (2) there was a relatively high resistivity zone at depth within the sliding zone. The deformation characteristics at the surface of the landslide were monitored by TLS and the results revealed that rainstorms caused three types of deformation and failure: (1) gully erosion at the slope surface; (2) shallow sliding failure; (3) and slope foot erosion. Subsequent monitoring of continuous changes in pore-water pressure, soil pressure and displacement (using traditional methods) indicated that long duration light rainfall (average 2.22 mm/d) caused the entire landslide to enter a state of creeping deformation at the beginning of the rainy season. Shear-induced dilation occurred for the fast sliding (30.09 mm/d) during the critical failure sub-phase (EF). Pore-water pressure in the sliding zone was affected by rainfall. In addition, the sliding L1 parts of the landslide exerted a discontinuous pressure on the L2 part. Through the monitoring and analysis, we conclude that this kind of landslide may have large deformation at the beginning and the late of the rainy season.     

Numerical modeling and dynamic analysis of the 2017 Xinmo landslide in Maoxian County,China

A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eighty-three people were killed or went missing and more than 103 houses were destroyed. In this paper, the geological conditions of the landslide are analyzed via field investigation and high-resolution imagery. The dynamic process and runout characteristics of the landslide are numerically analyzed using a depth-integrated continuum method and Mac Cormack-TVD finite difference algorithm.Computational results show that the evaluated area of the danger zone matchs well with the results of field investigation. It is worth noting that soil sprayed by the high-speed blast needs to be taken into account for such kind of large high-locality landslide. The maximum velocity is about 55 m/s, which is consistent with most cases. In addition, the potential danger zone of an unstable block is evaluated. The potential risk area evaluated by the efficient depthintegrated continuum method could play a significant role in disaster prevention and secondary hazard avoidance during rescue operations.     

A two-parameter exponential recession model for simulating cropland soil moisture dynamics

To simulate the soil moisture variation in cropland, a two-parameter exponential recession model was derived to depict the recession process of soil moisture in the root zone. The model is based on the assumption that the recession rate of soil water is propor- tional to the potential evapotranspiration rate and the difference of soil water content and steady soil water content. Two parameters in this model are soil texture-dependent recession constant and steady soil water content. The model was calibrated and validated with measured soil water data at two experiment sites in North China with different soil textures and cropping systems. Coefficients of de- termination between measured and model simulated soil water content were all greater than 0.7, indicating that both models gave satis- factory simulation results. Results showed that values of two parameters mentioned above are both larger for finer soil than those for coarser soil. At the same potential evapotranspiration rate and soil water content, the recession rate of finer soil is usually lower than that of coarser soil. The proposed model can be used in irrigation management to predict approximate date for irrigation, as well as be em- bedded into watershed hydrological models to estimate the antecedent precipitation index.     

基于广义回归神经网络的京津冀地区土壤湿度遥感逐日估算研究

土壤湿度是地表水热交换过程和水文循环中的一个关键组成部分,获取高时空分辨率的土壤湿度数据一直是当前研究的热点。SMAP（Soil Moisture Passive and Active）主被动微波土壤湿度产品的精度高,但存在着空间分辨率低和时间分辨率缺失的问题,这限制了其在区域尺度上的应用,为解决这一问题得到更高时空分辨率的土壤湿度产品,本文利用广义回归神经网络模型（GRNN）模拟了MODIS地表温度、反射率、植被指数光学/热红外遥感数据以及高程、坡度、坡向、经纬度数据与SMAP土壤湿度的关系,从而将京津冀地区SMAP L2土壤湿度产品的时间分辨率由不连续（4~20 d）提升至1 d,空间分辨率由3 km提升至1 km,并扩展其在京津冀地区的空间覆盖范围。研究发现:① GRNN模型总体验证结果表明土壤湿度估算值与SMAP原始值的相关性较高（r=0.7392）,均方根误差（RMSE）为0.0757 cm3/cm3;② 不同季节典型日期的GRNN模型估算结果精度相差较大,春季处的相关性相比其他季节最低,精度相对较高（r_春=0.6152,RMSE_春=0.0653cm3/cm3）,秋季和夏季土壤湿度估算精度较为接近（r_夏=0.6957,r_秋=0.7053,RMSE_夏=0.0754cm3/cm3,RMSE_秋=0.0694cm3/cm3）,冬季的估算精度最高（r_冬=0.8214,RMSE_冬=0.0367cm3/cm3）;③ 2016年京津冀夏秋季节的土壤湿度较其他季节要显著提高,空间分布上坝上高原区域较低,而沿海地区的土壤湿度明显较高。本研究对京津冀地区的生态水文、气候预测以及干旱监测等应用领域具有重要价值。     

全球36 km格网土壤水分逐日估算

土壤水分是陆面生态系统和能量循环的核心变量之一,利用微波遥感技术获得的土壤水分产品的时间分辨率一般是2-3 d,因此精确地获得具有较高时间分辨率的土壤水分成了人们关注的焦点。本文尝试将SMAP (the Soil Moisture Passive and Active）土壤水分和MODIS光学数据相结合,利用广义回归神经网络进行全球36 km土壤水分的估算,提升SMAP土壤水分的时间分辨率。结果显示,广义回归神经网络估算土壤水分与SMAP保持了高相关性(r = 0.7528),但其却保留了较高的误差 (rmse = 0.0914 m³/m³)。尽管如此,估算的土壤水分能够很好地保持SMAP土壤水分的整体空间变化,并且提升了土壤水分的时间分辨率（1 d）。此处,本文研究了SMAP土壤水分与MODIS光学数据之间的关系,这对今后利用机器学习进行SMAP土壤水分降尺度研究提供了重要的参考价值。     

基于物理模型试验的多层滑带滑坡变形演化特征

滑坡变形演化特征一直是滑坡灾害预测与防治领域急需解决的关键问题, 但对于多层滑带滑坡的变形演化特征却少有研究。以物理模型试验为手段建立了三层滑带滑坡物理试验模型, 完成了多层滑带滑坡变形演化全过程的模拟。基于PIV技术获取坡表位移数据, 通过柔性测斜仪监测滑坡深部位移, 同时布设土压力盒获取滑坡内部土压力的变化情况, 实现了多层滑带滑坡演化过程多参量数据分析。试验结果表明, 多层滑带滑坡破坏过程可分为初始、等速、加速和破坏4个阶段。不同破坏阶段滑坡的主要变形区域不同, 下层滑体受到上层滑体牵引作用, 在重力和推力作用下滑坡变形逐渐向浅层发展。变形过程中滑坡应力逐渐向滑带集中, 滑坡推力沿埋深方向呈多级梯形分布。加速变形阶段滑带处应力迅速增大, 滑坡体内产生多层应力集中带, 滑带位置推力变化与滑坡位移显著相关。      

呷爬滑坡滑带土蠕变特性及其稳定性

水库滑坡变形破坏受其岩土体蠕变特性及环境因素的影响。当滑坡进入加速变形阶段后,变形骤然增大,失稳概率增加。为了研究滑坡岩土体蠕变特性及其稳定性,选取锦屏一级水电站呷爬滑坡为研究对象,采用坡表位移监测曲线分析与室内三轴蠕变试验相结合的方法,建立了Burgers蠕变模型结合FLAC3D软件进行了滑坡稳定性研究。分析坡表位移-时间曲线发现,坡体变形特征与一般滑坡土体的蠕变特征具有相似性,滑带土室内三轴蠕变试验结果表明,滑带土变形可划分为瞬时蠕变、减速蠕变与稳定蠕变3个阶段,同时其瞬时变形量、稳定蠕变速率均随围压以及应力水平的增大而增大。基于滑带土蠕变特性的Burgers蠕变模型的计算结果,对比了常规强度折减法与考虑蠕变的强度折减法的滑坡稳定性系数,计算结果表明呷爬滑坡目前处于稳定状态,在一个计算周期内考虑蠕变的强度折减法较常规强度折减法的稳定性系数下降了0.04,因此,揭示滑坡土体蠕变特性并在此基础上研究其稳定性具有实际意义。      

‘Initial’ Soil Moisture Effects on the Climate in China——A Regional Climate Model Study

In this study,the effects of ‘initial’ soil moisture(SM) in arid and semi-arid Northwestern China on subsequent climate were investigated with a regional climate model. Besides the control simulations(denoted as CTL),a series of sensitivity experiments were conducted,including the DRY and WET experiments,in which the simulated ‘initial’ SM over the region 30 –50°N,75 –105°E was only 5% and 50%,and up to 150% and 200% of the simulated value in the CTL,respectively. The results show that SM change can modify ...     

Rapid susceptibility mapping of co-seismic landslides triggered by the 2013 Lushan Earthquake using the regression model developed for the 2008 Wenchuan Earthquake

The primary objective of landslide susceptibility mapping is the prediction of potential landslides in landslide-prone areas.The predictive power of a landslide susceptibility mapping model could be tested in an adjacent area of similar geoenvironmental conditions to find out the reliability.Both the 2008 Wenchuan Earthquake and the 2013 Lushan Earthquake occurred in the Longmen Mountain seismic zone,with similar topographical and geological conditions.The two earthquakes are both featured by thrust fault and similar seismic mechanism.This paper adopted the susceptibility mapping model of co-seismic landslides triggered by Wenchuan earthquake to predict the spatial distribution of landslides induced by Lushan earthquake.Six influencing parameters were taken into consideration: distance from the seismic fault,slope gradient,lithology,distance from drainage,elevation and Peak Ground Acceleration(PGA).The preliminary results suggested that the zones with high susceptibility of coseismic landslides were mainly distributed in the mountainous areas of Lushan,Baoxing and Tianquan counties.The co-seismic landslide susceptibility map was completed in two days after the quake and sent to the field investigators to provide guidance for rescue and relief work.The predictive power of the susceptibility map was validated by ROC curve analysis method using 2037 co-seismic landslides in the epicenter area.The AUC value of 0.710 indicated that the susceptibility model derived from Wenchuan Earthquake landslides showed good accuracy in predicting the landslides triggered by Lushan earthquake.     

Regional scale rainfall- and earthquake-triggered landslide susceptibility assessment in Wudu County, China

Wudu County in northwestern China frequently experiences large-scale landslide events.High-magnitude earthquakes and heavy rainfall events are the major triggering factors in the region.The aim of this research is to compare and combine landslide susceptibility assessments of rainfalltriggered and earthquake-triggered landslide events in the study area using Geographical Information System(GIS) and a logistic regression model.Two separate susceptibility maps were produced using inventories reflecting single landslide-triggering events,i.e.,earthquakes and heavy rain storms.Two groups of landslides were utilized: one group containing all landslides triggered by extreme rainfall events between 1995 and 2003 and the other group containing slope failures caused by the 2008 Wenchuan earthquake.Subsequently,the individual maps were combined to illustrate the locations of maximum landslide probability.The use of the resulting three landslide susceptibility maps for landslide forecasting,spatial planning and for developing emergency response actions are discussed.The combined susceptibility map illustrates the total landslide susceptibility in the study area.     

‘Initial’ soil moisture effects on the climate in China

In this study, the effects of ‘initial’ soil moisture (SM) in arid and semi-arid Northwestern China on subsequent climate were investigated with a regional climate model. Besides the control simulations (denoted as CTL), a series of sensitivity experiments were conducted, including the DRY and WET experiments, in which the simulated ‘initial’ SM over the region 30–50°N, 75–105°E was only 5% and 50%, and up to 150% and 200% of the simulated value in the CTL, respectively. The results show that SM change can modify the subsequent climate in not only the SM-change region proper but also the far downstream regions in Eastern and even Northeastern China. The SM-change effects are generally more prominent in the WET than in the DRY experiments. After the SM is initially increased, the SM in the SM-change region is always higher than that in the CTL, the latent (sensible) heat flux there increases (decreases), and the surface air temperature decreases. Spatially, the most prominent changes in the WET experiments are surface air temperature decrease, geopotential height decrease and corresponding abnormal changes of cyclonic wind vectors at the mid-upper troposphere levels. Generally opposite effects exist in the DRY experiments but with much weaker intensity. In addition, the differences between the results obtained from the two sets of sensitivity experiments and those of the CTL are not always consistent with the variation of the initial SM. Being different from the variation of temperature, the rainfall modifications caused by initial SM change are not so distinct and in fact they show some common features in the WET and DRY experiments. This might imply that SM is only one of the factors that impact the subsequent climate, and its effect is involved in complex processes within the atmosphere, which needs further investigation.     

基于DEM的福建省土质滑坡敏感性评价

已有滑坡敏感性研究中对评价指标的选取可以归结为气象、水文、地形、地质、植被、人类活动等方面,这些因子指标来源不一,在缺少数据资料地区难以完整收集。针对这个问题,考虑到目前DEM数据的广泛可获得性及其对滑坡评价的重要性,本文仅利用DEM数据及其派生因子,研究土质滑坡敏感性评价的可行性。研究中把评价因子分为2组：第1组数据仅由DEM派生,包括高程、坡度、坡向、地形起伏度、曲率、水流强度指数（Stream Power Index, SPI）、沉积运输指数（Sediment Transport Index, STI）、地形湿度指数（Topographic Wetness Index, TWI）;第2组数据作为对照组,除了包括上述DEM派生的8个因子外,同时加入植被覆盖度、土地利用、土壤类型、年均降雨量因子。本文分别选取逻辑回归模型和证据权法,基于上述2组评价因子,以德化县为例对比2组因子评价结果,利用第1组和第2组数据进行滑坡敏感性评价,结果精度分别为73%和83%。结果表明,仅利用DEM数据进行土质滑坡敏感性评价方法可行,可以为缺乏资料区滑坡敏感性评价提供借鉴。