基于非线性回归和地理探测器的复杂地形系统中山洪暴发的动态及其影响机制

《Journal of Environmental Management》:Flash flood dynamics in complex terrain system and its impact mechanism based on non-linear regression and geodetector

【字体: 时间:2026年01月05日 来源:Journal of Environmental Management 8.4

编辑推荐:

  云南山洪灾害的驱动机制与空间分布研究,基于1955年以来的灾害点数据、地探测器分析及双指数模型与幂律模型模拟,揭示其集中在24°N-27°N、1500-2000米海拔、坡度≤25°及年降水680-1410毫米区域。驱动因素呈现非线性双机制:降水与人类活动符合双指数模型(y=a·ebdm+k),且因子影响力顺序为致灾因子>承灾体>环境背景。研究建立“驱动函数-交互机制”框架,为复杂地形区山洪风险防控提供科学依据。

  
姚高|彭家斌|陈俊旭|吴建婷|朱妍|季轩|赵璐媛|王子尧
中国-老挝-孟加拉-缅甸自然资源遥感监测云南国际联合实验室,云南省教育厅喀斯特过程与水陆资源配置重点实验室,云南大学,昆明,650091,中国

摘要

山洪暴发对生命和财产构成重大威胁,尤其是在地形复杂的地区。本研究旨在利用灾害点调查数据(自1955年以来)、地理探测器分析和函数拟合模拟,揭示中国云南省山洪暴发的空间分布、动态及其驱动因素。研究结果表明,山洪暴发主要集中在北纬24°至27°之间的中部地区,特别是在海拔1500–2000米、坡度≤25°、年降水量在680至1410毫米之间的区域。非线性曲线拟合分析揭示了两种不同的驱动机制:降水量和人为因素最好用双指数模型(y = a * e^b * x + c * e^d * x)来表示,该模型描绘出一个“最大敏感性窗口”,在该窗口内洪水概率在特定阈值附近达到峰值;而地表环境因素则遵循幂律衰减模型(y = l * x^m + k),表明风险随着因素强度的增加而降低。驱动因素的影响顺序为:引发灾害的因素、承载灾害的物体因素以及易受灾的环境因素。这些因素之间的协同效应表现出非线性的增强作用。在北部流域,山洪暴发主要由降水量和地表条件的共同作用驱动;在中部流域,降水量和人类活动起主导作用;而在南部流域,降水量是主要驱动因素。本研究为复杂地形系统中山洪暴发的机制和空间动态提供了新的见解,为山洪风险评估和制定缓解策略以保护云南省及类似地区的脆弱社区提供了宝贵的科学支持。

引言

山地地形在全球环境中起着关键作用,支撑着世界26%的人口,并间接维持了近一半人类的生计(Beniston, 2003)。在这些脆弱的环境中,山洪暴发是特别具有破坏性的自然灾害。它们由山区或陡峭地形的小流域中的短期强降雨引发,其特征是水位突然上升、持续时间短且破坏力强(Barthold et al., 2015; Borga et al., 2011)。山洪暴发对人类生命和财产构成严重威胁,加剧了人与土地之间的冲突,并阻碍了区域经济的有序发展(Yin et al., 2023)。根据EM-DAT-CRED数据库的数据,1900年至2020年间,与洪水相关的死亡人数占全球灾害相关死亡人数的30.34%,其中中国的比例高达60.42%(Guha-Sapir et al., 2015)。历史数据显示,自1949年以来,山洪暴发一直是中国洪水相关死亡的主要原因,特别是在2010年至2015年间,其贡献超过了60%(Liu et al., 2018a, Liu et al., 2018b; He et al., 2018; Liu and Huang, 2020)。与一般洪水灾害不同,山洪暴发发展迅速(Wei et al., 2022),破坏性极强,往往导致更严重的伤亡。在当前气候变化和人类活动的背景下,极端降水事件的频率和强度显著增加(IPCC, 2021; Li et al., 2022)。同时,山区的快速城市化深刻改变了地表条件(Schroeder et al., 2016),从而显著增加了山洪暴发的潜在风险。因此,准确理解山洪暴发的时空分布模式和驱动机制对于有效的灾害预防和缓解至关重要。
为应对这一挑战,全球学术界通过综合研究加深了对山洪暴发的理解。文献计量分析(图S1)显示,当前的研究主要集中在时空变异性(Diakakis et al., 2016; Wang et al., 2021b)、风险评估(Peng et al., 2023; Rosenzweig et al., 2018; Shirani et al., 2023)、早期预警和建模(Wang et al., 2023)以及因果分析(Duan et al., 2022; Hu et al., 2018b; Khosravi et al., 2018; Liu and Huang, 2020)方面。时间分析关注多尺度时间变化(Pereira et al., 2017),研究山洪暴发的时间变化模式,并探讨降水事件的空间变异性与其他因素之间的相互作用(Douinot et al., 2016)。空间差异研究采用标准差椭圆和聚类等技术来评估地理分布模式(Xiong et al., 2019)。在更详细的流域或沟壑尺度上,高分辨率数据提高了灾害敏感性和影响评估的准确性(Lian et al., 2017)。在分析驱动因素时,研究人员经常使用地理探测器(Liu et al., 2018b; Phillips, 2002)、线性/非线性和逻辑回归(Diakakis et al., 2016; Wang et al., 2021b),以及各种机器学习方法(Ahemaitihali and Dong, 2022; Band et al., 2020; Bui et al., 2019)。这些研究共同旨在量化山洪暴发与环境因素(如降水量、径流、地形、土壤、植被和人类活动)之间的复杂关系(Liu and Huang, 2020; Piras et al., 2016; Rosenzweig et al., 2018; Qu et al., 2018; Youssef et al., 2016),认识到在这些条件下(如强降雨、陡坡、薄土层、稀疏植被和人为干扰),这些因素可以协同加速地表径流(Hoang and Liou, 2024)。风险评估研究通常整合风险、脆弱性和暴露维度(Gigovic et al., 2017),使用层次分析过程或熵等方法开发多指标综合评估系统(Chen et al., 2016; Rosenzweig et al., 2018)。
尽管取得了显著的研究进展,但仍存在一些关键空白。大多数研究集中在典型的河流或单个河流流域(Cheng et al., 2023; Li et al., 2021),而对像云南省(Chen et al., 2023, Peng et al., 2026)这样的复杂山区系统的全面比较研究仍然有限。该地区的破碎地形、多个共存的河流流域和强烈的构造活动尚未得到充分探索。此外,虽然触发机制(如降水量(Nie et al., 2025; Wei et al., 2018)以及个别环境因素(如土壤(Posner et al., 2014; Uber et al., 2018)和植被(Ebel, 2020)的影响已经得到了充分研究,但山洪暴发的发生本质上是复杂的,从整合灾害诱发因素、暴露元素和环境易感性的系统视角来看,量化多个因素之间的非线性相互作用的研究仍然有限(Ruin et al., 2008; You et al., 2025)。这些灾害影响因素的相互作用以及灾害的强烈区域异质性给明确的因果归因带来了巨大挑战(Xiong et al., 2020; Zhang et al., 2024)。值得注意的是,尽管使用计算流体动力学(CFD)和先进机器学习方法(Bui et al., 2019, Luo et al., 2025)的精细物理模拟正在受到关注,但在山洪暴发研究中,将地理探测器等统计工具与非线性回归模型结合使用以表征空间异质性和检测因素相互作用的情况仍然很少见。这种方法可以构建一个“驱动函数-相互作用机制”框架,为这些复杂现象提供新的见解。
云南省的地形受到青藏高原和云贵高原抬升的影响。该省南北延伸960公里,海拔从76米急剧上升到6740米,形成了多样化的地形和气候带(Chen et al., 2023)。大多数居民生活在盆地和河谷中,使他们极易受到山洪暴发的影响(Xiong et al., 2019)。过去65年的观测数据显示,云南省中部地区的山洪暴发空间分布发生了显著变化。这种变化与多种因素密切相关,包括归一化植被指数(NDVI)、海拔、地形起伏、人口密度、坡度和降水量(Yang et al., 2024)。然而,这些动态变化与其驱动因素之间的内在关系尚未得到彻底分析,导致在这种复杂地形中山洪暴发的精细空间差异规则和形成机制尚不清楚。
本研究旨在阐明降水量、人类活动、地形和植被覆盖如何影响山洪的形成,并考察这些因素在不同空间尺度上的不同影响。具体目标包括:(1)揭示复杂地形系统中山洪的空间差异模式;(2)量化导致山洪发生的关键环境因素,确定非线性响应阈值、峰值、拐点等影响事件概率的特征;(3)确定云南省省级和流域尺度上的主要驱动机制及其相互作用。在方法论上,本研究通过整合复杂地形、降水量和人类活动等因素,突出了山洪的系统和复杂性,以研究其发生机制。预期这些发现将为分析山洪灾害的驱动因素以及制定区域预防和控制策略提供有价值的见解。

研究区域

研究区域

云南省(北纬21°8′32″–29°15′8″,东经97°31′39″–106°11′47″),位于中国西南边缘(图1)。其地貌受到青藏高原和云贵高原抬升的影响,形成了从西北向东南倾斜的山脉、高地和深谷。海拔在大约960公里的范围内从6740米急剧上升到76米。山区占全省面积的88.6%,而平原、丘陵和台地分别占4.8%、4.9%

材料与方法

本研究利用热点分析和空间自相关性来研究云南省山洪暴发的空间分布模式。随后进行了拟合分析,以评估山洪暴发概率与各种影响因素之间的关系。此外,还应用了地理探测器驱动分析,以探索省级和三级流域尺度上这些因素背后的驱动力。主要内容和方法包括

云南省山洪暴发的空间分布

云南省的山洪暴发频率在0到23次之间,表现出显著的空间异质性(图3a)。较高的频率集中在昭通东部、曲靖北部、昆明东北部、玉溪东部、楚雄北部和保山中部等地。相反,文山、西双版纳、普洱中部和丽江北部的频率较低。为了确定最佳的空间分析单元,进行了全球空间自相关性分析

结论

本研究构建了一个系统的“驱动函数-相互作用机制”框架,以量化云南省山洪暴发的非线性控制因素。通过空间分析、拟合函数和地理探测器,我们研究了云南省山洪暴发的影响因素和动态,重点关注空间分布模式、影响因素、驱动力及其协同效应。本文为研究山洪暴发的影响机制提供了新的思路

CRediT作者贡献声明

姚高:撰写——原始草稿、可视化、软件。彭家斌:撰写——原始草稿、可视化、软件。陈俊旭:撰写——审稿与编辑、监督、资金获取、概念化。吴建婷:撰写——审稿与编辑、调查、数据管理。朱妍:验证、软件。季轩:监督、方法论、调查。赵璐媛:软件、形式分析。王子尧:调查、形式分析、数据管理。

写作过程中生成式AI和AI辅助技术的声明

在准备这项工作时,作者使用了Qwen3工具来改进语言表达。使用该工具后,作者根据需要审查和编辑了内容,并对出版物的内容负全责。

利益冲突声明

作者声明他们没有已知的竞争性财务利益或个人关系可能影响本文所述的工作。

致谢

本手稿得到了云南省基础研究项目(编号202301BF070001-004;202201AT070094)、国家自然科学基金(编号42161006)、云南省高层次人才专项计划(编号C6213001159)、中国地质调查局批准的地质调查工作项目(编号DD20242172)的支持,以及欧洲研究委员会执行机构项目(编号101040939)和云南省科学研究基金的部分支持
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普
  • 急聘职位
  • 高薪职位

知名企业招聘

热点排行

    今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号