《Water》：GIS-Based Identification of Critical Points and Contributing Areas for Concentrated Runoff Hazard Screening
编辑推荐：
由集中地表径流引发的山洪是影响小流域和建成区的主要自然灾害之一。本研究提出了一套更新后的捷克共和国关键点（Critical Points, CPs）及其汇水区识别的国家方法。该方法整合了基于激光雷达（LiDAR）的数字地形模型、基于国家地形数据库ZABAGED
由集中地表径流引发的山洪是影响小流域和建成区的主要自然灾害之一。本研究提出了一套更新后的捷克共和国关键点（Critical Points, CPs）及其汇水区识别的国家方法。该方法整合了基于激光雷达（LiDAR）的数字地形模型、基于国家地形数据库ZABAGED?的建成区更新边界，以及代表两种极限水力连通性条件的水文 conditioned 地形模型。研究人员开发了一套新程序，用于在涵洞化河段及其他排水中断处 delineate 汇水区，从而更精细地表达集中径流路径。与2009年原始数据集的比较表明，由于输入数据更新和方法改进，关键点的数量及空间分布发生了显著变化。两种水力连通性情景在全国尺度上产生了相似的空间格局，但在局部改变了径流路径、关键点位置及汇水区边界。所形成的数据集为山洪易发性筛选和减灾措施的空间定位提供了全国性框架。将其整合到捷克水文气象研究所的业务化山洪指标（Flash Flood Indicator）中，进一步证明了该方法在空间精细化预警中的应用价值。关键点方法可推广至具备合适地理空间数据的地区，但可能需要根据当地条件进行适应性调整。
论文解读：捷克全国尺度集中径流关键点与汇水区识别方法更新
一、研究背景与目标
山洪因其突发性强、强度高且可预报性低，成为最具破坏性的自然灾害之一。山洪的生成不仅受短历时强降雨的控制，还显著受下垫面条件影响，包括地形、土地利用、土壤性质及排水结构等。随着气候变化，与山洪生成相关的短历时极端降水在许多地区呈增强趋势。然而，小流域尤其是无水文观测的流域中，详细观测数据往往稀缺，基于地形和水文地貌学的筛选方法因此成为事件式水文与水动力学模拟的重要补充。捷克共和国原有的关键点方法通过识别集中地表径流抵达建成区的位置，并结合汇水区自然与土地利用特征评估其潜在危险性，已应用于欧盟洪水指令实施、洪灾风险管理、土地整理及土壤侵蚀防治等。但自2009年全国评估以来，土地利用和聚落格局已发生变化，更高精度的地形数据已经可得，小型水利构筑物对径流连通性的影响也需要更明确的表达。因此，本研究旨在提出一套更新的关键点识别方法，以反映当前景观条件，并融入改进的地形表达与水文连通性评估。
二、研究内容与主要结论
研究人员在捷克全国范围内开展了系统性的关键点及汇水区更新识别工作。研究以国家激光雷达高程数据DMR 5G为基础生成5米分辨率数字地形模型（DTM），结合ZABAGED?地形数据库更新建成区边界，并构建了两种极限水力连通性情景：情景A假定所有涵洞、桥梁及暗涵河段完全通畅，径流经这些结构连续输送；情景B假定这些结构在极端条件下被泥沙和杂物堵塞，径流漫过结构沿替代地表路径流动。为实现这两种情景，研究人员生成了三种经水文 conditioning 的DTM变体，并专门开发了跨暗涵河段汇水区 delineation 程序。关键点定义为集中径流路径与建成区边界的交点，其汇水区通过分水岭分析划定，并采用临界条件指标（F）及辅助阈值筛选最终数据集。质量控制流程结合系统几何筛查和人工目视检查，剔除了高度可疑的关键点。最终，情景A保留9919个关键点，情景B保留9997个关键点，其汇水区累计范围分别约占捷克国土面积的27.0%和27.4%。与2009年数据集相比，在统一采用F≥1.85阈值时，更新后的情景B识别出13443个关键点，明显多于2009年的9261个，表明输入数据和方法改进显著影响了关键点的数量与空间分布。两种情景的全国尺度格局高度相似，但局部径流路径、关键点位置和汇水区边界存在差异，凸显了排水基础设施对空间明确径流路径分析的重要性。研究还发现，已提议或实施的防洪和侵蚀控制措施与高F值汇水区的重叠比例随F值增大而升高，说明这些措施已有一定程度的针对性布局。此外，更新后的关键点汇水区已整合进捷克水文气象研究所的业务化山洪指标（FFI），2024年6月27日Blí?ejov山洪案例表明，基于汇水区的FFI比原规则网格版本能更准确地定位高风险区域。该研究发表于《Water》。
三、主要关键技术方法
该研究采用的主要技术方法包括：基于机载激光扫描点云生成5米分辨率数字地形模型（DTM），并利用ArcGIS Pro及Arc Hydro Pro扩展进行水文预处理；通过规则化GIS程序从ZABAGED?数据库提取并泛化建成区多边形；对DTM沿地表河流、涵洞和桥梁进行烧录（burning）以实施水文 conditioning，构建两种极限水力连通性情景；使用D8算法计算流向和汇流累积量，提取集中径流路径；采用专门的工作流处理暗涵河段上方的地表径流路径重建与汇水区分配；基于分水岭分析 delineate 各关键点汇水区；利用核密度估计（KDE）和皮尔逊相关系数等比较数据集空间分布；使用临界条件指标（F）对汇水区进行危险度排序，并通过系统几何筛查和人工目视检查进行质量控制。
四、研究结果分述
1. 情景A和B下的关键点与汇水区
研究人员采用修订后的F≥4.00阈值确定最终数据集。情景A获得9919个关键点，汇水区累计面积21323 km2；情景B获得9997个关键点，累计面积21649 km2。两情景下汇水区的自然地理和土地利用特征非常接近，平均汇水面积分别为324.4 ha和326.9 ha。随着F阈值提高，保留的关键点数量和汇水区累计范围逐步减少，而汇水区内耕地比例则逐步增加，表明F指标与农业用地密切相关。
2. 径流路径情景比较
情景A和B的全国统计差异较小，但空间配置存在局部差异。约90%的累计汇水区范围在两情景中重合，其余区域因水力连通性假设不同而改变了径流路径、关键点位置或汇水区边界。核密度表面显示两情景高度空间对应（r=0.9691，RMSE=0.0151），全国44.7%地区密度变化可忽略，情景B在38.4%地区产生更高或明显更高的关键点密度，在17.0%地区则更低。这表明替代假设主要是在局部重新分配径流路径和相关关键点，而非在全国尺度上系统性改变易发性。
3. 与2009年关键点数据集比较
以10 ha共同最小汇水阈值衡量，更新后的径流路径网密度从2009年的1.94 km/km2增加到2.29 km/km2，增幅约18%。建成区边界变化对潜在径流—聚落交点数量影响更强：固定2009年建成区边界时，更新径流网使交点数量增加约12%；固定径流网时，用当前建成区边界替换2009年边界可使交点数量增加约22–24%。核密度比较显示情景B与2009年数据集空间对应r=0.8884，RMSE=0.0292；情景B在45.3%国土面积上显示更高或明显更高的关键点密度，29.8%面积上更低，山区尤其是捷克北部和东部增幅最大。这些差异应视为地形数据改进、建成区和土地利用数据修订以及方法改进共同作用的结果。
4. 实际应用
在捷克的关键点方法实际应用中，汇水区为结构性、自然型和农业型减灾措施提供了空间定位依据，包括滞洪设施、截流沟、沉沙结构、草带、等高耕作、条带种植以及限制陡坡种植玉米等易蚀作物。对过去17年间已提议或实施的防洪和侵蚀控制措施与汇水区重叠分析表明，提议措施覆盖的汇水区比例随F值增大而升高，从F=4–10时的约3–5%增至F>45时的15–20%以上，最高达27.5%（情景B，F=46–48）。已实施措施的覆盖比例明显低于提议措施，多数F等级下不足10%，表明大部分高易发汇水区仍缺乏已完工的防治工程。此外，更新后的CP数据已纳入捷克水文气象研究所业务化山洪指标（FFI）。Blí?ejov案例中，两场连续对流单体在汇水区产生超过80 mm雷达—雨量计降水，局部可能超过100 mm，前期湿润度较高。更新后的FFI将两个受影响汇水区判定为极高危险，而原3×3 km网格版仅显示低至中等危险，因为固定网格将降水在较大面积上平均化了。
五、讨论与结论总结
讨论部分指出，更新后的关键点数据集是改进的地形表达、修订的建成区与土地利用数据以及方法学改进的综合成果，不应简单视为原始清单的时间更新。水力连通性情景应理解为极限状态表达，而非单次事件预测，其作用在于识别对小型水利设施功能假设敏感的区域。临界条件指标F虽能综合多种相关特征，但森林或山地集水区可能在低F值下仍产生危险径流，因此辅助选择标准用于降低遗漏风险。实践方面，CP汇水区为减灾措施定向和评估提供了框架，且捷克境内独立事件汇编和FFI业务化集成均支持其实际相关性。国际比较显示，CP方法与基于DEM的流动路径法、地貌计量法和复合指数法原理相近，但独特之处在于先识别抵达建成区的集中径流路径，再划定每个关键点对应的水文连通汇水区，并显式表达涵洞、桥梁和暗涵河段的影响。限制方面，该方法本质上是空间筛选和优先级排序工具，不模拟径流随时间演化、洪水波传播或淹没动态，也尚未使用全国一致性的事件清单进行系统验证。未来研究应侧重于将静态水文地貌框架与动态气象水文信息（如降水和前期土壤湿度）结合，并探索将基于CP的空间单元与数据驱动方法融合的混合建模途径。最终结论强调，该国家数据集为识别和优先排序与建成区相连的径流源区提供了实用基础，已应用于防洪和侵蚀措施规划，并通过整合进山洪指标支持业务化预警；方法可向其他地区转移，但需根据当地条件和数据可获得性调整标准和输入数据处理方式。
