极端暴雨下典型小流域重力侵蚀的分布及影响因素.pdf

第3 7卷第4期2 0 2 3年8月水土保持学报J o u r n a l o fS o i l a n dW a t e rC o n s e r v a t i o nV o l.3 7N o.4A u g.,2 0 2 3 收稿日期:2 0 2 2-1 1-0 7 资助项目:国家自然科学基金项目(4 2 1 7 7 3 2 7);黄河水科学研究联合基金重点项目(U 2 2 4 3 2 1 2)第一作者:喻涵(2 0 0 0),女,硕士研究生,主要从事土壤侵蚀研究。E-m a i l:y u h a n 0 61 6 3.c o m 通信作者:韩剑桥(1 9 8 7),男,博士,副研究员,主要从事流域洪水泥沙灾害研究。E-m a i l:h j q 1 31 6 3.c o m张曼(1 9 8 9),女,博士,助理研究员,主要从事农田碳氮循环研究。E-m a i l:m a n m a n 0 2 0 31 2 6.c o m极端暴雨下典型小流域重力侵蚀的分布及影响因素喻 涵1,2,周子渊3,王 一4,韩剑桥1,2,张 曼1,2,党维勤3(1.西北农林科技大学水土保持研究所,陕西 杨凌7 1 2 1 0 0;2.中国科学院水利部水土保持研究所,陕西 杨凌7 1 2 1 0 0;3.黄河水利委员会绥德水土保持科学试验站,陕西 榆林7 1 9 0 0 0;4.山东省莱阳市沐浴水库管理中心,山东 烟台2 6 4 0 0 0)摘要:为探究极端暴雨下重力侵蚀分布规律及其环境条件,以2 0 1 9年“利奇马”台风引发的特大暴雨事件为对象,采用无人机航拍、室内解译等方法,研究了位于暴雨中心的山东省临朐县曾家沟小流域内不同类型重力侵蚀的数量、空间分布特征及其影响因素。结果表明:(1)小流域内发生的重力侵蚀共7 0处,主要类型为滑坡(4 9处),其次为泥石流(1 4处)和崩塌(7处)。(2)重力侵蚀面积、侵蚀量分别与坡度变率、坡度、坡位、高程、坡向和坡向变率呈正相关,与距梯田距离、距道路距离和土地利用类型呈负相关。(3)坡度是影响重力侵蚀的最重要因素,侵蚀易发生在坡度为3 0 5 0 的斜坡。(4)人类活动对重力侵蚀的影响较大,特大暴雨下的严重重力侵蚀多发生在梯田田坎、道路边坡等位置,需将重力侵蚀防治的重心放在加强道路两侧边坡、梯田田坎和坡度3 0 的陡坡边坡防护上。研究结果可为应对和预防重力侵蚀灾害提供科学依据。关键词:台风“利奇马”;极端暴雨;重力侵蚀;坡度;防灾减灾中图分类号:S 1 5 7.1 文献标识码:A 文章编号:1 0 0 9-2 2 4 2(2 0 2 3)0 4-0 0 6 9-0 6D O I:1 0.1 3 8 7 0/j.c n k i.s t b c x b.2 0 2 3.0 4.0 1 0D i s t r i b u t i o na n dI n f l u e n c i n gF a c t o r so fG r a v i t a t i o n a lE r o s i o ni nT y p i c a l S m a l lW a t e r s h e dU n d e rE x t r e m eR a i n s t o r mYU H a n1,2,Z HOUZ i y u a n3,WANGY i4,HANJ i a n q i a o1,2,Z HANG M a n1,2,D ANG W e i q i n3(1.I n s t i t u t eo fS o i la n dW a t e rC o n s e r v a t i o n,N o r t h w e s tA&FU n i v e r s i t y,Y a n g l i n g,S h a a n x i7 1 2 1 0 0;2.I n s t i t u t eo fS o i la n dW a t e rC o n s e r v a t i o n,C h i n e s eA c a d e m ya n dS c i e n c e sa n dM i n i s t r yo fW a t e rR e s o u r c e s,Y a n g l i n g,S h a a n x i7 1 2 1 0 0;3.S u i d eS o i la n dW a t e rC o n s e r v a t i o nS c i e n t i f i cE x p e r i m e n t a lS t a t i o no ft h eY e l l o wR i v e rC o n s e r v a n c yC o mm i s s i o n,Y u l i n,S h a a n x i7 1 9 0 0 0;4.M a n a g e m e n tC e n t e ro fM u y uR e s e r v o i r,Y a n t a i,S h a n d o n g2 6 4 0 0 0)A b s t r a c t:T oe x p l o r e t h ed i s t r i b u t i o no f g r a v i t a t i o n a l e r o s i o na n d i t s e n v i r o n m e n t a l c o n d i t i o n su n d e r e x t r e m er a i n s t o r m.T h ee x t r e m e l yr a i n s t o r me v e n tc a u s e db yT y p h o o n“L e k i m a”i n2 0 1 9w a sc h o s e na st h et a r g e t.U s i n gU A Va e r i a l i m a g e r ya n dr e m o t es e n s i n gi n t e r p r e t a t i o n,t h eq u a n t i t i e s,s p a t i a ld i s t r i b u t i o nc h a r a c t e r i s t i ca n di n f l u e n c i n gf a c t o r so fv a r i o u sg r a v i t a t i o n a le r o s i o ni nZ e n g j i a g o us m a l lw a t e r s h e do fL i n q uC o u n t y,S h a n d o n gP r o v i n c e,w h e r et h ec e n t e ro fe x t r e m e l yr a i n s t o r m w a sl o c a t e d.T h er e s u l t ss h o w e dt h a t:(1)T h e r ew e r ea t o t a l o f 7 0g r a v i t a t i o n a l e r o s i o n s i nZ e n g j i a g o u.Am o n gw h i c h,t h ep r e d o m i n a n t t y p ew a sL a n d s l i d e(4 9s i t e s),f o l l o w e db yd e b r i sf l o w(1 4s i t e s)a n dc o l l a p s e(7s i t e s).(2)G r a v i t a t i o n a le r o s i o na r e aa n de r o s i o na m o u n tw e r ep o s i t i v e l yc o r r e l a t e dw i t hs l o p ev a r i a t i o n,s l o p e,s l o p ep o s i t i o n,e l e v a t i o n,s l o p ed i r e c t i o na n ds l o p ed i r e c t i o nv a r i a t i o n,r e s p e c t i v e l y,a n dn e g a t i v e l yc o r r e l a t e dw i t hd i s t a n c e t o t e r r a c e,d i s t a n c e t or o a da n d l a n du s et y p e.(3)S l o p ew a st h em o s t i m p o r t a n t f a c t o ra f f e c t i n gg r a v i t a t i o n a l e r o s i o n,w h i c hw a sp r o n e t oo c c u ro ns l o p e s r a n g i n gf r o m3 0 t o5 0 .(4)G r a v i t a t i o n a le r o s i o nh a db e e np r o f o u n d l yi n f l u e n c e db yh u m a na c t i v i t i e s.U n d e rh e a v yr a i n s t o r m,s e v e r eg r a v i t ye r o s i o nm o s t l yo c c u r r e di nt e r r a c e df i e l da n d r o a ds l o p e s.T h eg r a v i t a t i o n a l e r o s i o np r e v e n t i o na n dc o n t r o l s h o u l db e f o c u s e do ns t r e n g t h e n i n g t h ep r o t e c t i o no nb o t hs i d e so f t h es l o p er o a d,t e r r a c e df i e l d,a n ds t e e ps l o p e sw i t hs l o p e3 0 .T h e r e s u l t sc a np r o v i d es c i e n t i f i cb a s i s f o rc o p i n gw i t ha n dp r e v e n t i n gg r a v i t a t i o n a l e r o s i o nd i s a s t e r s.K e y w o r d s:t y p h o o n“L e k i m a”;e x t r e m er a i n s t o r m;g r a v i t a t i o n a le r o s i o n;s l o p e;d i s a s t e rp r e v e n t i o na n dr e d u c t i o n 重力侵蚀是山区普遍的土壤侵蚀类型之一1,其形式包括崩塌、滑坡、泥石流等,危害巨大2。降雨是重力侵蚀的主要触发因素3,大多数重力侵蚀往往发生在极端暴雨期间或之后1,4。在全球气候变化的背景下,气候变暖导致水文循环发生变化,极端暴雨事件增加5,造成重力侵蚀灾害频发。因此,探寻极端暴雨下重力侵蚀的发生特点和影响因素,对重力侵蚀防灾减灾具有重要的理论和实践意义。除降雨外,地形地貌和植被也是重力侵蚀形成、演变、分布的重要影响因素6-7。陡峭地形增加了发生重力侵蚀的可能性,而重力侵蚀对地表的长期作用导致破碎地貌的形成,同时削弱了沟坡的稳定性,即两者相互影响、互为因果8。坡度对重力侵蚀总量、滑坡侵蚀量和崩塌侵蚀量的影响较大6。坡向和坡位的不同导致坡面土壤含水量等因素之间存在差异,间接影响植被生长状况,进而造成重力侵蚀的不同9。而植被对沟坡重力侵蚀的影响则具有两面性,植被根系固结土壤,使重力侵蚀量减少3,但植被促进地表径流入渗,又加剧重力侵蚀的发生1 0。由于重力侵蚀成因复杂,且区域地形、土地利用等环境因子均有差异,极端暴雨下重力侵蚀问题应因地制宜分析。受2 0 1 9年第9号台风“利奇马”影响,山东局部地区发生特大暴雨。其中,位于暴雨中心的临朐县次降雨量超过5 7 0mm,占临朐县多年平均降雨量的8 2%,达到百年一遇标准,引发了严重的土壤侵蚀灾害1 1,为研究极端暴雨下的重力侵蚀灾害提供了典型案例。因此,本文选择山东临朐县曾家沟小流域为研究对象,基于无人机航拍和室内解译,结合不同类型重力侵蚀的数量和分布特征,采用冗余分析研究坡向、坡度、高程等9个因素对重力侵蚀的影响程度及差异,以期为应对和预防重力侵蚀灾害提供科学依据。1 研究区概况曾家沟小流域位于3 6 1 2 5 7 3 5 1 3 3 5 N,1 1 8 3 4 0 8 1 1 8 3 5 2 7 E,总面积为1.6 1k m2,隶属山东省潍坊市临朐县九山镇(图1)。海拔3 6 8 6 6 4m,地貌以低山丘陵为主。气候属亚湿润大陆性季风气候,雨热同期,年平均气温1 2.8,年最大降水量为12 0 0mm。流域内多为棕壤土,植被以刺槐、黑松、蜀桧、荆条等为主。曾家沟小流域为北方土石山区,属于国家级水土流失重点治理区。在台风“利奇马”影响下,小流域受灾严重,台风路径和影响范围见图1。图1 研究区概况2 研究方法2.1 数据来源在特大暴雨发生5天后(2 0 1 9年8月1 5-2 2日),对曾家沟小流域进行了现场调查和测量,运用D J IP h a n t o m 4R T K无人机对曾家沟流域进行正射影像测量。采用P i x 4 Dm a p p e r建模,完成相机参数的解算与图像拼接,生成0.0 7m分辨率的曾家沟小流域无人机航拍正射影像(D OM)与0.0 9m分辨率的数字地表模型(D S M)。2.2 重力侵蚀解译基于曾家沟小流域正射影像,结合现场调查情况进行目视解析,解译标志见表1。2.3 重力侵蚀量计算重力侵蚀量计算公式为:M=V(1)式中:M为重力侵蚀量(t);V为重力侵 蚀土方量(m3);为土方容重(t/m3)。通过A r c G I S1 0.7软件中填挖方工具得到各侵蚀图斑的重力侵蚀土方数1 2。根据已有研究1 3,研究区土壤容重设为1.3 6t/m3。2.4 环境因子的提取重力侵蚀是各种因素综合作用的结果,参考已有研07水土保持学报 第3 7卷究1 4,本研究从地形地貌和人类活动等方面选取9个环境因子进行分析,分别是坡向、坡度、高程、坡位、坡向变率、坡度变率、土地利用,以及侵蚀点距梯田和道路的距离。其中,坡向和坡度基于D S M数据由栅格表面工具提取。变率由表面分析模块提取,坡度变率指地面坡度在微分空间的变化率1 5,坡向变率指坡向在水平方向上的变率,可以反映出地表所有的山脊线、山谷线,影响地表径流收敛和发散1 6。坡位在坡向数据的基础上经栅格计算得到,按距坡顶水平距离三等分,由上至下为上坡、中坡和下坡。道路距离、梯田距离由距离分析模块提取。土地利用数据通过正射影像经目视解译获得,结合研究区水土保持措施,将小流域土地利用方式划分为5类:梯田、坡耕地、林地、草地和其他用地。以上所有数据的提取均使用A r c G I S1 0.7软件完成。表1 重力侵蚀类型解译标志重力侵蚀类型影像特征现场判别方法崩塌影像色调:浅白色或深灰色纹理:有凹凸感,粗糙感或呈花斑状,既有大块石,也有小的砾石形状:边缘清晰,下部堆积物散乱,呈锥状岩土体在较陡坡上崩落和滚动,堆积在下部,崩塌体各部分相对位置完全打乱,大小混杂,形成大块滚动较远小块较近的堆积体滑坡影像色调:既有白灰色,也有滑塌物的绿色纹理:不清晰,明显凹凸不平形状:不规则,下部存在堆积物岩土体沿一定滑坡面顺坡向下滑动到下部。滑坡体多数在滑坡面上有残留,也有相应的植物生长。表面有一定纵横裂缝泥石流影像色调:较为一致,灰白色纹理:有凹凸感,不规则形状:山坡陡峻,岩石破碎强烈,有沟槽形状,下部有堆积物存在,无植被生长泥石流发生后有明显的形成区、流动的路径和下部有大量散乱的堆积物2.5 数据分析本文基于冗余分析(r e d u n d a n c ya n a l y s i s,R D A)确定环境因子对重力侵蚀的贡献1 7。分别采用去趋势对应分析法(d e t r e n d e dc o r r e s p o n d e n c ea n a l y s i s,D C A)和蒙特卡洛(m o n t ec a r l o)置换检验法检验冗余分析的可行性和显著性1 8。蒙特卡洛置换检验p0.0 5时在统计学上通常被认为是可靠的1 9。将重力侵蚀量和侵蚀面积作为响应变量,坡向、坡度等9个环境因子作为解释变量。本研究的D C A结果前4轴的长度均小于3,因此重力侵蚀与坡度等9个环境因子之间的关系可以通过冗余分析反映。同时,冗余分析结果也通过了蒙特卡洛检验(F=4.3 4 3,p=0.0 0 2 5 0(1 5.7%),侵蚀面积分别为40 6 0.2,77 7 8.4,89 8 7.6m2,侵蚀量分别为58 1 7.0,78 5 0.1,53 3 9.4t。重力侵蚀面积和侵蚀量随着坡 度变率的增 加表现为上 升趋势(图2 b)。当重力侵蚀发生在坡度变率8 5 9 0 时,侵蚀面积最广、水土流失最严重,侵蚀面积为1 14 9 5.3m2,占侵蚀总面积的5 1.5%,侵蚀量为90 1 3.2t,占侵蚀总量的3 9.6%。从坡向(图3)来看,重力侵蚀发生在北坡的可能性最高(3 0.0%),其次为东南(1 2.9%)和东北坡(1 2.9%)。北坡侵蚀面积为1 45 8 0.5m2,侵蚀量为1 06 2 8.6t。随着坡向变率的增加,重力侵蚀面积和侵蚀量均呈现17第4期 喻涵等:极端暴雨下典型小流域重力侵蚀的分布及影响因素上升趋势(图4)。重力侵蚀集中发 生在坡向变 率8 5 9 0 的斜坡上(6 8.6%),侵蚀面积为1 78 6 7.6m2,占侵蚀总面积的8 0.1%,侵蚀量为1 71 4 9.6t,占侵蚀总量的7 5.3%。图2 不同坡度和坡度变率侵蚀特征图3 不同坡向侵蚀特征及频数分布图4 不同坡向变率侵蚀特征研究区的高程在3 6 86 6 3m,重力侵蚀主要发生在高程4 5 05 0 0m的范围内(图5),大多处于小流域中坡位(图6),共4 9处重力侵蚀发生在中坡位,侵蚀面积为1 41 6 0.8m2,占侵蚀总面积的6 7.9%,侵蚀量为1 58 2 1.0t,占侵蚀总量的6 7.9%,造成大范围、严重的水土流失,其次是下坡位(1 2处),侵蚀面积为17 2 7.0m2,侵蚀量为32 6 7.3t。发生在上坡位的重力侵蚀,单位侵蚀面积和单位侵蚀量最大,分别为6 0 3.1m2/处,4 1 0.5t/处。3.2.2 不同土地利用类型的重力侵蚀特征 通过解译小流域重力侵蚀地块上的土地利用,表明发生重力侵蚀的土地利用类型主要包括梯田、林地、草地和其他土地(图7)。在该次暴雨中,重力侵蚀集中发生在梯田(2 8处),侵蚀面积为1 08 0 1.1m2,侵蚀量为1 04 7 6.0t。而发生在草地的重力侵蚀仅1 0处,少于林地,但其侵蚀面积(47 1 2.3m2)和侵蚀量(38 3 7.5t)均超过林地(25 2 1.4m2,33 7 7.7t)。图5 不同高程侵蚀特征图6 不同坡位侵蚀特征27水土保持学报 第3 7卷图7 不同土地利用类型侵蚀特征3.2.3 人类活动影响下的重力侵蚀特征 由图8可知,重力侵蚀集中分布在距梯田4 0m范围内(9 4.3%)。其中,发生在梯田四周01 0m内的侵蚀最为严重,侵蚀面积为1 30 2 4.7m2,侵蚀量为1 35 2 9.9t。其次是距梯田1 02 0m范围内,侵蚀面积为48 0 1.4m2,侵蚀量为4 8 1 5.7t。随着距道路距离的增加,侵蚀面积和侵蚀量呈减少的趋势(图9)。6 7.1%的重力侵蚀分布在距离道路4 0m范围内,侵蚀面积为1 61 0 9.9m2,占侵蚀总面积的7 2.2%;侵蚀量为1 49 4 2.3t,占侵蚀总量的6 5.6%。5.7 1%的侵蚀分布在距离道路8 0m范围外,侵蚀面积为11 0 3.9m2,侵蚀量为13 0 6.4t。3.3 环境因子对重力侵蚀的贡献重力侵蚀与环境因子的冗余分析结果(图1 0)显示,第1轴和第2轴分别解释了3 8.1%和1.4%的变异。第1轴中贡献最大的是坡度,第2轴中距道路和梯田距离与典范轴的相关系数分别为0.5 4 8,0.5 1 5。环境因子中坡度对重力侵蚀的影响最显著(p=0.0 0 1),其次是坡度变率(p=0.0 1 3)和坡向(p=0.0 3 2)。重力侵蚀面积、侵蚀量均与坡度变率、坡度、坡位、高程、坡向和坡向变率表现为正相关,与距梯田距离、距道路距离和土地利用呈负相关。重力侵蚀量与坡度、坡度变率的相关性较强,与高程、坡位的相关性较弱;侵蚀面积与坡度间夹角小,相关性较强,与高程夹角接近9 0,呈较弱的负相关关系。图8 距梯田不同距离重力侵蚀频数图9 距道路不同距离侵蚀特征 注:蒙特卡洛检验F=4.3 4 3,p=0.0 0 2;实线箭头为解释变量;虚线箭头为响应变量。图1 0 重力侵蚀与环境因子的R D A二维排序图4 讨 论坡度是重力侵蚀最重要的限制性因素。重力侵蚀主要发生在3 0 5 0 的斜坡上,随着坡度的升高,侵蚀量和侵蚀面积持续增加。但当坡度5 0 时,侵蚀发生概率却呈下降趋势。当坡度5 0 时,重力分量的作用会阻止土层停留在5 0 的斜坡上,导致重力侵蚀发生的概率降低1 4,2 0。坡度变率反映地面坡度的变化率,其对重力侵蚀面积的影响大于侵蚀量。坡度变率较高的地区在降雨后易发生侵蚀2 1,且坡度的剧烈突变增加了侵蚀力2 2,造成严重的水土流失。由于小流域的北坡光照时间短,且在“利奇马”台风的影响下,降雨量超过5 7 0mm,强降雨无法完全下渗,导致地表径流增加。同时,北坡的土石道路和耕地汇集的雨水来不及随着排水沟流走,最终致使小37第4期 喻涵等:极端暴雨下典型小流域重力侵蚀的分布及影响因素流域北坡梯田田坎、道路两侧滑塌的概率增加。在全球气候变化的背景下,极端暴雨事件加重区域水土流失防治的难度2 3。目前,人们难以直接干预重力侵蚀发生的内外营力,因此改变人类活动应成为重力侵蚀防治的重心2 4。基于本研究结果,建议在道路两侧、坡度3 0 的陡坡采用防护网等边坡防护措施,增加坡面植被覆盖,改善土壤结构,在边坡上部设置截、排水沟等;可为现有梯田设置暗管排水系统,排出土壤多余水分,并采用种植苗木等措施防护梯田田坎。5 结 论(1)小流域内重力侵蚀受灾严重。侵蚀共7 0处,主要类型为滑坡(4 9处),其造成的水土流失危害和影响范围最大,其次是泥石流(1 4处)和崩塌(7处)。(2)重力侵蚀的分布受环境因素影响。从坡度来看,重力侵蚀集中分布于3 0 5 0 的斜坡上。从坡向和坡位来看,重力侵蚀发生在北坡和中坡位的概率最高。侵蚀位置多集中在高程4 5 05 0 0m、距道路4 0m的范围内。当坡度变率和坡向变率为8 5 9 0 时,发生重力侵蚀的概率更大。在不同的土地利用类型中,梯田最易发生重力侵蚀。(3)坡度是影响重力侵蚀的最重要因素。重力侵蚀面积、侵蚀量均分别与坡度变率、坡度、坡位、高程、坡向和坡向变率呈正相关,与距梯田距离、距道路距离和土地利用类型呈负相关。(4)人类活动导致发生重力侵蚀的概率增大。特大暴雨下,严重的重力侵蚀多发生在梯田田坎、道路边坡等位置。因此,建议加强坡度3 0 的陡坡边坡防护,增加 坡面植被覆 盖,在边坡上 部设置截、排水沟;可为现有梯田设置暗管排水系统,加强梯田田坎防护。参考文献:1 X uXZ,L i uZY,X i a oPQ,e ta l.G r a v i t ye r o s i o no nt h es t e e pl o e s ss l o p e:B e h a v i o r,t r i g g e ra n ds e n s i t i v i t yJ.C a t e n a,2 0 1 5,1 3 5:2 3 1-2 3 9.2 高航,徐向舟,肖培青,等.暴雨条件下植被对沟坡重力侵蚀速率和规模的影响J.水土保持研究,2 0 2 1,2 8(6):1 7-2 4.3 胡胜,邱海军,王宁练,等.地形对黄土高原滑坡的影响J.地理学报,2 0 2 1,7 6(1 1):2 6 9 7-2 7 0 9.4 G u oWZ,C h e nZX,W a n gWL,e t a l.T e l l i n gad i f f e r-e n t s t o r y:T h ep r o m o t er o l eo fv e g e t a t i o ni nt h ei n i t i a-t i o no f s h a l l o wl a n d s l i d e sd u r i n gr a i n f a l lo nt h eC h i n e s eL o e s sP l a t e a uJ.G e o m o r p h o l o g y,2 0 2 0,3 5 0:e 1 0 6 8 7 9.5 舒章康,李文鑫,张建云,等.中国极端降水和高温历史变化及未来趋势J.中国工程科学,2 0 2 2,2 4(5):1 1 6-1 2 5.6 余璐,徐向舟,张茂省,等.黄土高原重力侵蚀对地貌因素的敏感性分析J.水土保持学报,2 0 1 9,3 3(4):1 1 9-1 2 5.7 高健健,艾琦森,韩立钦,等.2 0 1 72 0 1 8年黄土丘陵沟壑区第一副区辛店沟流域重力侵蚀观测数据集J.中国科学数据(中英文网络版),2 0 2 1,6(3):1 1 3-1 2 0.8 薛海,何沛华,王文成.黄河中游地形河网分形与重力侵蚀关系研究J.人民黄河,2 0 0 8,3 0(9):6 7-6 8.9 李琦.黄土丘陵地区重力侵蚀区域植被分布特征研究D.陕西 杨凌:西北农林科技大学,2 0 1 9.1 0 赵兴阳,徐向舟,蒋云钟,等.暴雨条件下植被对黄土沟坡重力侵蚀的影响J.水土保持学报,2 0 2 0,3 4(1):5 8-6 3.1 1 林祚顶,刘宝元,丛佩娟,等.山东临朐2 0 1 9年“81 0”特大暴雨水土保持调查J.水土保持学报,2 0 2 1,3 5(1):1 4 9-1 5 3.1 2 李小光,郭少武,宋星杰,等.基于A r c G I S的露天矿山土石方量测算方法研究J.金属矿山,2 0 2 0(1 2):2 1 0-2 1 5.1 3 从辰宇,韩剑桥,焦菊英,等.台风“利奇马”暴雨引发的土壤侵蚀调查研究:以山东省临朐县为例J.水土保持通报,2 0 1 9,3 9(5):3 3 7-3 4 4.1 4 C h e nX,C h e n W.G I S-b a s e dl a n d s l i d es u s c e p t i b i l i t ya s s e s s m e n tu s i n go p t i m i z e dh y b r i d m a c h i n el e a r n i n gm e t h o d sJ.C a t e n a,2 0 2 1,1 9 6:e 1 0 4 8 3 3.1 5 李强,张景发,罗毅,等.2 0 1 7年“88”九寨沟地震滑坡自动识别与空间分布特征J.遥感学报,2 0 1 9,2 3(4):7 8 5-7 9 5.1 6 李星,杨赛,李远耀,等.面向区域滑坡易发性精细化评价的改进斜坡单元法J.地质科技通报,2 0 2 3,4 2(3):8 1-9 2.1 7 龙健,廖洪凯,李娟,等.基于冗余分析的典型喀斯特山区土壤-石漠化关系研究J.环境科学,2 0 1 2,3 3(6):2 1 3 1-2 1 3 8.1 8 G u o M M,Y a n gB,W a n g W L,e ta l.D i s t r i b u t i o n,m o r p h o l o g ya n di n f l u e n c i n gf a c t o r so fr i l l su n d e re x t r e m er a i n f a l l c o n d i t i o n s i nm a i n l a n du s e s o n t h eL o e s sP l a t e a uo fC h i n aJ.G e o m o r p h o l o g y,2 0 1 9,3 4 5:e 1 0 6 8 4 7.1 9 刘振坤,吴华勇,刘峰,等.中国焦化场地近2 0年时空演变特征及驱动因素J.生态环境学报,2 0 2 1,3 0(3):6 0 4-6 1 3.2 0 田媛,巨能攀,解明礼,等.滑坡编录表达模式对易发性评价结果的影响J.成都理工大学学报(自然科学版),2 0 2 2,4 9(5):6 0 6-6 1 5.2 1 南富森,李宗省,张小平,等.基于MGWR的黄河流域土壤有机碳密度和影响因素分析J.环境科学,2 0 2 2,4 4(2):9 1 2-9 2 3.2 2 W a n gZG,Z h a n gG H,W a n gCS,e t a l.A s s e s s m e n to ft h eg u l l ye r o s i o ns u s c e p t i b i l i t yu s i n gt h r e eh y b r i dm o d e l si no n es m a l lw a t e r s h e do nt h eL o e s sP l a t e a uJ.S o i l a n dT i l l a g eR e s e a r c h,2 0 2 2,2 2 3:e 1 0 5 4 8 1.2 3 黄萱.植被覆盖和降雨变化对流域水沙过程的影响机制D.武汉:华中农业大学,2 0 1 9.2 4 张金良.黄河泥沙入黄的机理及过程探讨J.人民黄河,2 0 1 7,3 9(9):8-1 2.47水土保持学报 第3 7卷