北京山区暴雨泥石流激发雨量条件

涂剑，马超，杨海龙

(1.北京林业大学 水土保持国家林业局重点实验室，100083，北京；2.北京林业大学 北京市水土保持工程技术研究中心，100083，北京； 3.北京林业大学水土保持学院，100083，北京)

北京山区暴雨泥石流激发雨量条件

涂剑，马超†，杨海龙

(1.北京林业大学 水土保持国家林业局重点实验室，100083，北京；2.北京林业大学 北京市水土保持工程技术研究中心，100083，北京； 3.北京林业大学水土保持学院，100083，北京)

泥石流是北京山区常见的一种山地灾害，也是威胁山区生态环境的头号安全问题。分析该区的泥石流激发雨量条件对于该区山洪泥石流预测预报具有指导意义。笔者搜集18场降雨数据，分析该区泥石流的激发雨量特征。结果表明：激发降雨的峰值降雨时段介于3～9 h，峰值时段平均降雨强度达22.6～50.0 mm/h，累积雨量达91.9～350.0 mm，峰值降雨时段雨量占累积雨量的63.9%～100.0%，这说明激发泥石流的降雨过程具有历时短、雨量集中、突发性强的特点。用修正后的李氏法进行雨场分割，得到该区泥石流激发的降雨强度-历时和累积雨量-历时关系。最后，以密云区北部山区为例分析3期植被盖度和泥石流激发雨量条件。归一化植被指数和像元线性分解模型分析表明，1989—1991、2004—2005和2011—2012年植被覆盖度分为64.48%、68.57%和72.48%；相应的泥石流雨量激发条件也随着植被覆盖度的增加而增大，可见植被恢复盖度对泥石流的形成具有一定影响。本文统计分析得出的泥石流激发雨量条件可为该区山洪泥石预测预报和监测预警指标选取提供参考。

暴雨； 泥石流； 降雨特征； 雨量激发条件； 植被盖度； 北京山区

泥石流预警预报是防灾减灾的主要手段之一，国内外针对泥石流的激发雨量条件开展了较多研究。大体上，泥石流激发雨量条件获取手段可分为2种：一种是从泥石流形成机制出发，根据降雨诱发浅表层滑坡模型和沟床物质起动的水力学模型反算得到雨量指标[1-2]，这种方法适用于缺少长系列降雨资料地区[3]；另一种是统计学方法，根据前期有效降雨、降雨强度等各类组合指标之间的关系绘制雨量阈值曲线[4-5]，该方法需要长序列资料和灾害资料，相应的雨量指标也较多，如10 min、1 h降雨量，累积降雨量，前期降雨量等[6-8]。其中，平均降雨强度-历时(I-D)表征短历时、高强度降雨激发条件，适用于泥石流暴发时段与峰值降雨时段重合的情况[9]。

泥石流是北京山区常见的一种山地灾害，也是威胁山区生态环境的头号安全问题。近年来，因极端暴雨天气以及人类工程活动的叠加影响，泥石流事件逐渐增多。据统计：1971—1980，1981—1990，1991—2000和2001—2010年间，分别有3、7、9、10场泥石流暴发[10-12]。2011年“7·24”和2012年“7·21”2场特大暴雨引发的山洪泥石流导致密云区龙潭沟、房山区十渡镇受损严重，造成巨大经济损失。此外，自2000年以来，北京市对北京山区开展了植被恢复、小流域治理等工作。截至2009年，植被覆盖度从59%上升至72%[13]。植被在很大程度上抑制浅表层滑坡和崩塌，有效地控制了区域水土流失，对山区生态环境恢复起到了积极作用。随着泥石流孕育环境发生变化[14]，泥石流的临界雨量激发条件也会发生变化。目前，北京地区现有泥石流雨量激发条件研究大多为植被恢复前，缺乏植被恢复后的激发条件研究。鉴于此，笔者通过收集激发泥石流灾害的雨量过程，分析激发雨量特征，采用统计方法开展该区泥石流的激发雨量条件研究，以密云区为例，分析不同植被盖度下泥石流激发雨量条件变化，以期为该区山洪泥石流预测预报、监测预警提供参考。

1 材料与方法

1.1 北京山区概况及雨量数据获取

北京山区面积达1万417.5 km2，占北京市总面积的62%，是我国暴雨泥石流多发区之一。大部分泥石流分布在白草畔与桃山连成的一条1 000 m高程的障碍线附近(图1)，包括北部的云蒙山至海沱山一带及西部的百花山至笔架山一带[10,15]。目前，北京市气象局、区(县)防汛办以及国家气象局在北京山区共计有300多个雨量站点，这些站点大多位于山洪泥石流集中暴发区(图1)。结合泥石流灾害点及雨量站分布，共搜集了1989—2012年间18次泥石流事件的雨量资料(表1)。其中：1989、1991年雨量及灾害资料获取自《密云县统计年鉴》《密云县水利志》；2002—2012年雨量资料购买自北京市气象局、国家气象局，同时结合中央气象数据网获得的《中国暴雨洪涝灾害数据集》整理而成；泥石流灾害点资料则获取自北京各区统计年鉴及实地调查。泥石流事件的日期、地点以及附近雨量站点见表1。

图1 气象站与泥石流灾害点的空间位置关系Fig.1 Spatial distribution of the meteorological stations and the debris flow sites

1.2 降雨参数和阈值方法选取

北京山区土层薄，黏粒、粉粒等细颗粒物质含量较少，物源体主要以粗颗粒为主。泥石流主要是由于高强度的降雨转化为地表径流掀揭沟床物质形成[9,15]。此外，区域泥石流多集中在6—8月份暴发，此段时间的降雨占全年降雨74.9%，其他时段降雨较少[15]；因此，该区的泥石流受前期降雨的影响较小，与降雨过程中峰值降雨时段的高强度降雨密切相关。平均降雨强度-历时(I-D)关系表征短历时、高强度降雨激发条件，适用于北京山区泥石流暴发时段与峰值降雨时段重合的情况[12]；选取降雨强度-降雨历时(I-D)、累积雨量-历时(C-D)2种雨量特征参数，采用单线法确定该地区泥石流的雨量激发条件。

2 北京山区降雨特征

2.1 激发泥石流的降雨特征

图2为1989—2012年间18次泥石流事件累积降雨过程。累积雨量曲线的斜率代表此时段的降雨强度，斜率越大，降雨强度越大，斜率最高处为该场降雨的峰值时段。大多数降雨过程曲线都有1个斜率较高的时段(峰值降雨时段)。如2012年“7·21”特大暴雨中：房山(当天11:00—13:00)及河北镇(当天08:00—10:00)记录的3 h内最大雨量分别达到178.0和247.9 mm，远远超过中央气象台制定的暴雨红色预警标准(3 h内降雨达到100 mm以上)。这18场暴雨泥石流激发的降雨过程中(表2)：3 h最大降雨量超过100 mm(红色预警)的场次达11场，最高为2012年7月21日的河北镇，达247.9 mm；其余7场降雨中，最大3 h雨量最低的为1991年6月10日密云四合堂村，为79.9 mm，但四合堂村整个降雨历时7 h内均维持了较高水平，总累积雨量达161 mm。2011年6月23日门头沟的累积雨量曲线趋于平缓，但仍然激发了泥石流。这是因为该场降雨突发性强，降雨过程中第1小时(当天09:00)雨量达到57.2 mm，3 h内累积雨量达99.0 mm。可见，北京山区激发泥石流的降雨过程具有历时短、雨量集中、突发性强的特点。

2.2 激发泥石流的雨量条件

对18次暴雨泥石流的降雨过程进行处理，从雨量站获取的原始数据不能直接作为分析泥石流发生的降雨资料，需要用雨场分割方法将其分割成不同时间段的雨量数据。目前，用于连续降雨的雨场分割方法有6种[16]，但这些方法均不能体现历时短、雨量集中、突发性强特征。如李氏法是以1场降雨结束为连续3 h雨量低于4 mm对雨场进行分割，采用该方法会延长降雨历时使平均降雨强度变小；因此，本研究采用中央气象台制定的大中雨分界值8 mm/h作为标准，即在一场连续降雨过程中以时降雨量<8 mm处为有效降雨的结束(称为修正法)，以能更好体现北京山区峰值时段的降雨特征。

用修正雨场分割法对图2中18次暴雨泥石流的降雨过程进行处理(表2)。可以看出，北京山区暴雨泥石流峰值时段持续仅3～9 h，峰值时段平均降雨强度达22.6～50.0 mm/h，峰值时段累积雨量达91.9～350.0 mm，雨量大且集中，占该次降雨过程的63.9%～100%。

将表2中获得的平均降雨强度/历时数据绘于双对数坐标系中(图3)，将累积雨量及其降雨历时的数据点都绘于双对数坐标系中(图4)，根据各个时段的下限值拟合出1条临界雨量线。泥石流发生的雨量数据点(平均降雨强度及累积雨量)均应在临界雨量线的上方，即当泥石流流域内雨量数据点(平均降雨强度及累积雨量)及其对应的降雨历时的相关点在临界雨量线上方时，该沟暴发泥石流的可能性较大。

对图3和图4中的暴雨泥石流平均降雨强度及累积雨量的下限值进行拟合，得出以下关系式：

I-D类型I=35.4D-0.18。

(1)

C-D类型C=35.4D0.82。

(2)

式中：I为平均降雨强度，mm/h；C为累积雨量，mm；D为降雨历时，h。

表1泥石流事件的日期、地点以及附近雨量站点

Tab.1Occurrence dates and locations of debris flow events, as well as the locations of rainfall stations

编号No.日期Data区域(气象站点)Location(Weatherstation)泥石流灾害点Debrisflowsite泥石流沟与雨量站位置的距离Distancebetweenrainfallstationanddebrisflowsite/km11989-07-21密云番字牌村(番字牌)PanzipaiinMiyun(Panzipai)番字牌西沟、小西天等PanzipaiWestGullly,Xiaoxitian,etc.0.51～5.12密云冯家峪镇(西白莲裕)FengjiayuinMiyun(XibailianyuWatershed)大烂碴沟、王虎沟等DalanzhaGullly,WanghuGullly,etc.0.72～3.4521991-06-10密云四合堂村(白河四合堂)SihetanginMiyun(SihetangofBaiRiver)水泉沟、骡子沟等ShuiquanGullly,LuoziGullly,etc.0.33～3.37密云张家坟村(捧河岩村)ZhangjiafeninMiyun(Pengheyan)黄土板地村、小南沟等Huangtubandi,XiaonanGullly,etc.0.77～6.6532002-06-24门头沟区(达摩庄村)Mengtougou(Damozhuang)北山沟BeishanGullly1.8542002-08-01密云石城镇(石城)ShichenginMiyun(Shicheng)北峪沟BeiyuGully5.2052004-07-22密云北庄镇(北庄)BeizhuanginMiyun(Beizhuang)葡萄峪、达峪沟等Putaoyu,DayuGully,etc.3.11～4.2762005-08-12密云张家坟村(捧河岩村)ZhangjiafeninMiyun(Pengheyan)黄土梁沟、柳棵峪村等HuangtuliangGully,Liukeyu,etc.0.42～1.9172006-07-09门头沟区(王平)Mentougou(Wangping)草甸村、平原村等Caodian,Pingyuan,etc.1.15～2.1782011-06-23门头沟区(斋堂)Mentougou(Zhaitang)川底下村东EastofChuandixia3.3792011-07-24密云龙潭沟村(太师屯)LongtangouinMiyun(Taishitun)龙潭沟北沟、黑龙潭景区入口处等LongtangouNorthGully,theEntranceofHeilongtanScenic,etc.0.57～1.27102012-07-21门头沟区(门头沟国家站)Mentougou(MentougouNationalStation)桑峪村东、灵水村西等EastofSangyu,WestofLingshui,etc.2.01～7.17门头沟区(门头沟)Mentougou(Mentougou)黄塔村东、塔河村西EastofHuangta,WestofTahe0.45～5.14房山霞云岭乡(霞云岭国家站)XiayunlinginFangshan(XiayunlingNationalStation)栗树港沟、十八湖沟内等LishugangGully,ShibahuGully,etc.0.78～2.21房山区(房山国家站)Fangshan(FangshanNationalStation)河东村村道、三福村村道等TheroadsofHedongandSanfu,etc.1.14～2.23密云区(密云国家站)Miyun(MiyunNationalStation)区域内各大小沟均发生Allaroundthedistrict0～10.00房山河北镇(河北镇)HebeiTowninFangshan(HebeiTown)寺上村及宝地洼村村道TheroadsofSishangandBaodiwa1.34～4.10密云大城子镇(大城子)DachengziinMiyun(Dachengzi)庄户峪水库河沟、蔡峪南沟等ZhuanghuyuReservoir,CaiyuSouthGully,etc.1.50～4.50

图2 18次泥石流事件的累积降雨过程Fig.2 Cumulative rainfall process of 18 triggering debris flow events

日期Date区域Location3h累积雨量Accumulativerainfallin3hours/mm3h平均降雨强度Meanrainfallintensityin3hours/(mm·h-1)降雨峰值时长Peakrainfallduration/h峰值时段平均降雨强度Meanrainfallintensityatpeakrainfallintensity/(mm·h-1)峰值时段雨量Cumulativepeakrainfall/mm总累积雨量Totalcumulativerainfall/mm峰值时段雨量占总雨量比例Percentageofpeakrainfallintotalrainfall/%1989-07-21密云冯家峪镇FengjiayuinMiyun91.930.6330.691.9143.963.9密云番字牌村PanzipaiinMiyun116.138.7927.9251.2353.871.01991-06-10密云张家坟村ZhangjiafeninMiyun101.533.8823.2185.9219.084.9密云四合堂村SihetanginMiyun79.926.6723.0161.0161.0100.02002-06-24门头沟区Mengtougou114.638.2432.0128.1140.091.52002-08-01密云石城镇ShichenginMiyun190.663.5646.6279.3280.199.72004-07-22密云半北庄镇BeizhuanginMiyun86.528.8522.6113.0117.096.62005-08-12密云张家坟村ZhangjiafeninMiyun87.629.2524.7123.7145.385.12006-07-09门头沟区Mengtougou82.827.6429.8119.0122.597.12011-06-23门头沟区Mengtougou99.033.0333.099.0125.878.72011-07-24密云龙潭沟村LongtangouinMiyun85.628.5625.9155.6163.395.3门头沟国家站MentougouNationalSta-tion131.944.0437.5149.9200.674.7门头沟区Mengtougou132.344.1830.7245.8305.780.42012-07-21房山霞云岭乡站XiayunlinginFangshan177.359.1931.6284.4337.584.3房山区Fangshan178.059.3636.2217.0248.587.3密云区Miyun118.139.4627.7166.4201.882.5房山河北镇HebeiTowninFangshan247.982.6750.0350.0535.065.4密云大城子镇DachengziinMiyun10735.7627.7166185.589.5

图3 平均降雨强度与降雨历时Fig.3 Relationship between mean rainfall intensity and rainfall duration

图4 累积雨量与降雨历时Fig.4 Relationship between cumulative rainfall and rainfall duration

3 密云区激发泥石流的雨量条件变化

Ir: The maximum value of mean rainfall intensity determined with successive n hours. Ir1, Ir2, Ir3 indicates 2011—2012, 2004—2005, 1989—1991, respectively. Dr: Rainfall duration determined with successive n hours. The same below. 图5 密云区暴雨泥石流发生所需的平均降雨强度-历时关系下限Fig.5 Lower limit of the mean rainfall intensity and duration for debris flows occurring in Miyun between 1989 and 2012

为探讨该区植被恢复和激发泥石流发生的降雨特征值变化的关系，选取密云区1989—1991年(植被恢复前)4次、2004—2005年(植被恢复中)2次、2011—2012年(植被恢复后)3次，共9次降雨激发泥石流事件进行分析。由于激发泥石流的往往是因为峰值时段的强降雨而形成径流冲刷沟床沉积物形成的。采用获取连续n(n>1，n为自然数)小时最大平均降雨强度法，对图2中密云区降雨过程进行处理，将获得的平均降雨强度/降雨历时的数据点都绘于双对数坐标系中(图5)，分时段求出暴雨泥石流发生所需的平均降雨强度/历时关系。从中国科学院地理空间数据云获取密云区1989、2004和2010年3期Landsat数据。通过ENVY5.1、ArcMap10.3进行数据提取及归一化植被指数(Normalized Difference Vegetation Index，N)计算，将植被盖度分为0～30%，30%～60%，60%～80%和80%～100%共4个等级，以获取密云区这3个年份植被覆盖的时空动态(图6)。通过像元线性分解模型(式3)，得到植被覆盖各等级变化(表3)：

fc=(N-Nsoil)/(Nveg-Nsoil)。

(3)

式中：fc为植被覆盖度，N为观测像元的归一化植被指数，Nsoil为裸土的归一化植被指数，Nveg为100%覆盖条件下的归一化植被指数。

图5分别对2011—2012、2004—2005和1989—1991年中的暴雨泥石流平均降雨强度的下限值进行拟合，可以得出密云区这3个时间段泥石流发生的降雨条件：

(4)

(5)

(6)

式中：Ir1、Ir2和Ir3分别为采用连续1，2和3 h最大降雨强度法所获取的平均降雨强度，mm/h；Dr为连续n小时最大降雨强度的降雨历时，h。式(4)～(6)代表暴雨泥石流发生所需的平均降雨强度-历时关系下限。由于图5中的雨量数据点主要来源于密云北部山区，因而代表该区域中的泥石流激发雨量条件变化。可以看出：在1989—1991、2004—2005和2011—2012这3个时间段，泥石流激发雨量条件有升高的趋势。

如表3所示：植被覆盖度在30%～60%之间的山区面积占整个区域面积比例从1989年的24.82%降低到2010年的8.43%，而植被覆盖度>80%的山区面积比例却增加了。植被覆盖度<30%和60%～80%之间的比例相对持平。整体上，这3个时期的植被覆盖度逐渐增加。这说明，植被恢复后，相当部分的低覆盖度地区植被盖度增加了。进一步对植被覆盖度变化较大地区的分析表明，这些地方集中在密云北部山区，该地区白马关河流域、张家坟流域等是泥石流的多发区。比如：1989和1991年白马关河流域内番字牌、冯家峪镇、四合堂村以及小西天流域等爆发了大规模泥石流灾害，这一时期的植被覆盖度最低(64.48%)；2002—2005年该区只有石城镇、张家坟流域内爆发了泥石流，该时期植被覆盖度较1989—1991要稍高(68.57%)；2010—2012年植被覆盖度最高，只有零星的泥石流灾害暴发，规模相对于1989和1991年的群发性泥石流以及2002—2005年石城镇、张家坟泥石流而言，相对要小。

图6 密云区1989、2004和2010植被覆盖时空动态(Ⅰ：高覆盖度；Ⅱ：低覆盖度)Fig.6 Vegetation variation in 1989, 2004 and 2010 in Miyun (Ⅰ: High vegetation coverage; Ⅱ: Low vegetation coverage)

年份Year植被覆盖度等级Vegetationcoveragelevel/%植被覆盖条件Vegetationcoveragecondition激发雨量条件Rainfallconditionoftriggeringdebrisflows0～3030～6060～8080～100均值VMean/%归一化值Normalization表达式Expression归一化值Normalization2011—20124.558.4351.5835.44V1=72.481Ir1=37.9D-0.26r12004—20054.5416.9951.9226.55V2=68.570.41Ir2=34.6D-0.38r0.28～0.311989—19915.7824.8250.0919.31V3=64.480Ir3=33.1D-0.41r0

由于在该区的泥石流孕育环境中，植被是变化最为明显的因素。植被能通过固结地表土体、削减径流、截阻泥沙、削减水动力条件等抑制泥石流的形成。为定量评价植被覆盖度和泥石流雨量激发条件二者间的关系，通过式(7)将植被覆盖度及雨量激发条件进行归一化处理：

(7)

式中：x归为归一化值，取值为[0，1]；x为样本值；xmax为样本最大值，即Ir1或V1；xmin为样本最小值，即Ir3或V3。计算后得到激发雨量条件及植被覆盖条件归一化值(表3)：Ir2归一化值为0.28～0.31，小于V2的归一化值0.41。这说明，泥石流激发的雨量条件随着植被覆盖度的升高而增高。

4 结论与讨论

1)北京山区泥石流主要是由于高强度的降雨转化为地表径流掀揭沟床物质形成，激发泥石流的降雨过程具有历时短、降雨集中、突发性强的特点。笔者研究在分析北京山区降雨峰值特征的同时，提出了修正法分割雨场对雨量数据进行处理。采用平均降雨强度-降雨历时(I-D)、累积雨量-降雨历时(C-D)组合指标，用单线法拟合了泥石流发生的雨量激发条件：I=35.4D-0.18和C=35.4D0.82，该结果可为新建地质灾害预警系统提供参考。

2)通过归一化植被指数和像元线性分解模型分析了密云区泥石流激发雨量条件和植被盖度的关系。1989—1991、2004—2005和2011—2012年植被覆盖度分为64.48%、68.57%和72.48%。密云北部山区的覆盖度变化最为明显，而这些区域的泥石流雨量激发条件也增高了。由于在该区的泥石流孕育环境中，植被是变化最为明显的因素。植被能通过固结地表土体、削减径流、截阻泥沙、削减水动力条件等抑制泥石流的形成，在一定程度上能改变泥石流的激发雨量条件。

3)因本文采取统计分析的方法分析该区泥石流激发的雨量条件，选取的雨量站点虽然多位于泥石流多发区，但这些站点多位于泥石流沟沟口；而山区降雨条件具有非常明显的时空不均匀性，泥石流源区的降雨量、强度多大于中游和下游：因此，实际泥石流的激发雨量条件可能要比本研究中的要大。

4)笔者只收集了激发泥石流的降雨过程，为准确界定可能激发泥石流的雨量条件，应收集这些雨量站点中未激发泥石流的雨量过程，进而获取较为准确的泥石流雨量激发条件，并结合该区泥石流形成特征，建立基于泥石流起动机理的预警预报模型，为该区泥石流灾害预测预报提供支撑。

[1] JAKOB M, WEATHERLY H. A hydroclimatic threshold for landslide initiation on the North Shore Mountains of Vancouver, British Columbia [J]. Geomorphology, 2003, 54(3/4):137.

[2] BURTON A, BATHURST J C. Physically based modelling of shallow landslide sediment yield at a catchment scale [J]. Water Resources Research, 1998, 35(2):89.

[3] 潘华利,欧国强,黄江成,等. 缺资料地区泥石流预警雨量阈值研究[J]. 岩土力学,2012,33(7):2122.

PAN Huali, OU Guoqiang, HUANG Jiangcheng, et al. Study of rainfall threshold of debris flow forewarning in data lack areas [J]. Rock and Soil Mechanics, 2012, 33(7):2122.

[4] 白利平,孙佳丽,南赟. 北京地区泥石流灾害临界雨量阈值分析[J]. 地质通报,2008,27(5):674.

BAI Liping, SUN Jiali, NAN Yun. Analysis of the critical rainfall thresholds for mudflow in Beijing, China [J]. Geological Bulletin of China, 2008, 27(5):674.

[5] 庄建琦,崔鹏,葛永刚,等. 降雨特征与泥石流总量的关系分析[J]. 北京林业大学学报,2009,31(4):77.

ZHUANG Jianqi, CUI Peng, GE Yonggang, et al. Relationship between rainfall characteristics and total amount of debris flow [J] .Journal of Beijing Forestry University, 2009, 31(4):77.

[6] CHLEBORAD A F, BAUM R I, GODT J W. Rainfall thresholds for forecasting landslides in the Seattle, Washington, area-exceedance and probability [R]. U. S. Geological Survey Open-File Report 2006-1064, 2006.

[7] FIORILLO F, WILSON R C. Rainfall induced debris flows in pyroclastic deposits, Campania (southern Italy) [J]. Engineering Geology, 2004, 75(3/4):263.

[8] GUZZETTIL F, PERUCCACCIL S, ROSSIL M, et al. Rainfall thresholds for the initiation of landslides in central and southern Europe [J]. Meteorology and Atmospheric Physics, 2007, 98(3/4):239.

[9] MA Chao, WANG Yujie, Du Cui, et al. Variation in initiation condition of debris flows in the mountain regions surrounding Beijing [J]. Geomorphology, 2016, 273:323.

[10] 何媛,袁爱萍,胡影,等. 北京泥石流沟道现状调查及特征分析[J]. 中国水利,2013,(21):49.

HE Yuan, YUAN Aiping, HU Ying, et al. Investigation and analysis of the debris flow gullies in Beijing [J]. China Water Resources, 2013, (21):49.

[11] 吴斌,高甲荣,甘敬,等. 北京山区泥石流灾害现状及防治对策[J]. 中国水土保持科学,2008,6(4):1.

WU Bin, GAO Jiarong, GAN Jing, et al. Situation and countermeasures of debris flow disaster in Beijing mountainous area [J]. Science of Soil and Water Conservation, 2008, 6(4):1.

[12] 谢洪,钟敦伦,韦方强,等. 北京山区泥石流的分类与类型[J]. 山地学报,2004,22(2):212.

XIE Hong, ZHONG Dunlun, WEI Fangqiang, et al. Classification of debris flow in the mountains of Beijing [J]. Journal of Mountain Science, 2004, 22(2):212.

[13] 胡勇,刘良云,贾建华. 北京山区植被动态及生态恢复的遥感监测[J]. 应用生态学报,2010,21(11):2876.

HU Yong, LIU Liangyun, JIA Jianhua. Remote sensing based monitoring of vegetation dynamics and ecological restoration in Beijing mountainous area [J].Chinese Journal of Applied Ecology, 2010, 21(11):2876.

[14] 陈晓清,崔鹏,韦方强. 良好植被区泥石流防治初探[J]. 山地学报,2006(3):333.

CHEN Xiaoqing, CUI Peng, WEI Fangqiang. Study of control debris flow in high-covered vegetation region [J]. Journal of Mountain Science, 2006(3):333.

[15] 张嫱,马超,杨海龙,等. 北京山区典型低频泥石流特征及危险性研究 [J]. 北京林业大学学报,2015,37(12):92.

ZHANG Qiang, MA Chao, YANG Hailong. The risk assessment research of typical low frequency debris flow in Beijing mountainous region [J]. Journal of Beijing Forestry University, 2015, 37(12):92.

[16] 周伟,唐川. 汶川震区暴雨泥石流发生的降雨阈值[J]. 水科学进展,2013,24(6):786.

ZHOU Wei, TANG Chuan. Rainfall thresholds for debris flows occurrence in the Wenchuan earthquake area [J]. Advances in Water Science, 2013, 24(6):786.

RainfallconditionoftriggeringdebrisflowsinBeijingmountainregions

TU Jian，MA Chao，YANG Hailong

(1.Key Laboratory of State Forestry Administration on Soil and Water Conservation, Beijing Forestry University, 100083, Beijing, China;2.Beijing Engineering Research Center of Soil and Water Conservation, Beijing Forestry University, 100083, Beijing, China;3.School of Soil and Water Conservation, Beijing Forestry University, 100083, Beijing, China)

BackgroundDebris flows are one of the mountain hazards in mountain regions near Beijing, and cause serious threatening to local ecological safety. In recent years, the debris flow events increased a lot than last few decades as a result of the double influence of increasing extreme rainstorm events and human being activities. Researches on the triggering rainfall condition are the basis of early-warning and prediction of debris flows.MethodsUsing the distribution of debris flow sites and rainfall stations, 18 rainfall records were selected to analyze the rainfall characteristics. A rainfall event is defined when the hourly rainfall is greater than 8 mm/h in the beginning and ends in less than 8 mm afterwards. The rainfall records and rainfall events were used to identify the rainfall conditions of triggering debris flows in study area. Two rainfall thresholds, including the mean rainfall intensity-duration (I-D) and the accumulative rainfall-duration (C-D), were developed. The formula of rainfall thresholds with single-line method after analysis of theI-DandC-Dwere built. The Landsat data in year of 1989, 2004 and 2010 in Miyun District from geographic spatial data cloud of Chinese Academy of Sciences were analyzed. Through the ENVY5.1 and ArcMap10.3, the normalized difference vegetation index was calculated. Vegetation coverages in the three periods were obtained with pixel linear decomposition model. The relationship between triggering rainfall condition and vegetation coverage were built by standardized calculation.ResultsThe intensified rainfall duration lasted merely 3 h～9 h, during which the mean rainfall intensity ranged from 22.6 mm/h to 50.0 mm/h. Accumulative rainfall ranged from 91.9 mm to 350.0 mm, accounting for 63.9%～100% of total rainfall. Two functions with expressions asI-DandC-Din Beijing mountain regions were identified. Rainfall conditions of triggering debris flow in Miyun at 2009-2011, 2004-2005 and 1989-1991 were identified respectively. Vegetation coverage in Miyun at 2009, 2004 and 1989 were 72.48%, 68.57% and 64.48% respectively.ConclusionsThe results show that rainstorm in Beijing mountain regions characterizes in short duration, high intensity and high rainfall amount. The formula of rainfall thresholds with single-line method after analysis of theI-DandC-Dcan be used as a reference in forecasting of debris flows. Then, relationship between vegetation cover and the rainfall conditions of triggering debris flows in three stages of Miyun was studied. Results found that the vegetation cover degree in Miyun increased as time went on, corresponding rainfall triggering also increased. This study may provide a suggestion for the debris flow early-warning in the study area.

rainstorm; debris flow; rainfall characteristics; triggering rainfall condition; vegetation coverage; Beijing mountain regions

P642.23

A

2096-2673(2017)05-0103-08

10.16843/j.sswc.2017.05.013

2017-04-07

2017-08-02

项目名称： 国家自然科学基金“北方土石山区低黏度泥石流起动机理研究”(41702369)

涂剑(1991—)，男，硕士研究生。主要研究方向：山地灾害预警。E-mail:tj19911104@126.com

†

马超(1986—)，男，博士，讲师。主要研究方向：山地灾害预测预报。E-mai:sanguoxumei@163.com