多源高程模型不确定性向水文模拟与泥沙连通性传播的评估——以黄土区小流域为例

《Water》：Propagation of Uncertainty in Multi-Source Elevation Models into Hydrological Modelling and Sediment Connectivity in a Loess Catchment

　　

数字高程模型（Digital Elevation Model, DEM）的精度与分辨率决定了水文模拟所用输出的可靠性。然而，对其适用性的评估不应仅限于垂直精度。特别是在小型且高度破碎的流域中，高程表达的不确定性可根据局部地形结构传播至后续派生产品，并由此影响对径流组织和泥沙输移的解读。本研究旨在确定多源高程模型之间的差异如何从基础高程表达，经坡度和径流组织，传播至排水网络和泥沙连通性空间格局。研究人员在卢布林高地（波兰东南部）的Stocki Creek流域比较了ASTER GDEM、SRTM、EU-DEM和局部TIN模型，并以机载激光扫描（Airborne Laser Scanning, ALS）数据作为高分辨率参考来验证高程及派生输出。采用均方根误差（Root Mean Square Error, RMSE）和归一化中位绝对偏差（Normalized Median Absolute Deviation, NMAD）评估垂直精度，随后比较了地形起伏、坡度、水文路径输出、排水网络和泥沙连通性指数（Sediment Connectivity Index, SCI）的表达。与ALS相比，TIN的偏差最小（RMSE 1.48 m；NMAD 0.41 m），ASTER的偏差最大（6.78 m；2.37 m）。森林区域误差最高。所有模型均较好地再现了流域总体高程格局，而在依赖局部坡度和地形连续性的产品中，差异则更为显著。卫星模型平滑了沟壑的陡坡，概括了短径流线，并导致泥沙连通性空间格局的简化。结果表明，高程模型在水文应用中的适用性不仅应依据其垂直精度评估，还应考虑不确定性向后续处理产品的传播。因此，基于ALS的验证可评估多源对地观测（Earth Observation, EO）产品在整个分析链中的可靠性：从高程，经局部形态测量和径流路径，到排水网络和泥沙连通性。

研究背景与科学问题：数字高程模型（Digital Elevation Model, DEM）和数字地形模型（Digital Terrain Model, DTM）是水文模拟与地貌分析的主要地形信息来源，用于提取坡度、径流方向与累积量、流域边界、排水网络及侵蚀输沙指数。DEM质量对水文分析结果的影响虽已被长期认识，但以往评估多集中于垂直精度和分辨率，较少系统追踪同一分析链中多源高程模型差异的逐级传播。在黄土区高度破碎的小流域，沟壑密集、短径流路径众多，高程表达的不确定性可能通过局部地形结构影响径流组织与泥沙输移的解读。因此，研究人员以波兰东南部卢布林高地Na??czów高原Stocki Creek流域（约17 km²，沟壑密度局部可达10 km/km²）为研究对象，比较ASTER GDEM、SRTM、EU-DEM和局部TIN模型，并以ISOK项目机载激光扫描（Airborne Laser Scanning, ALS）数据作为高分辨率参考，检验差异如何沿“高程→坡度→径流路径→排水网络→泥沙连通性”链条传播。该研究不追求简单排序垂直精度，而是评估多源对地观测（Earth Observation, EO）产品在整套分析链中的可靠性。

关键技术与方法：研究人员保留各DEM原始分辨率，统一投影至PL-1992（EPSG:2180）并裁剪至共同范围；对每个模型执行相同的填洼（Fill）、流向（Flow Direction）、汇流累积（Flow Accumulation）流程。以ALS/ISOK生成的5 m和1 m DTM为参考，选取350个控制点，用RMSE和NMAD评估高程一致性；随后计算坡度、坡向、流域边界、排水网络（1000个汇流单元阈值）和基于Borselli等与Cavalli等方法的泥沙连通性指数（Sediment Connectivity Index, SCI），其中C因子采用ESDAC/JRC的泛欧洲覆盖管理因子栅格并统一重采样至各DEM网格。该研究为地理空间数据比较，不涉及人群队列。

研究结果：3.1 与ALS参考的高程一致性：TIN的RMSE和NMAD最小（1.48 m；0.41 m），ASTER最大（6.78 m；2.37 m），SRTM和EU-DEM居中（分别为5.03 m；1.48 m和4.32 m；1.44 m）。森林区误差最大，ASTER的RMSE达9.64 m，TIN为2.67 m；ASTER高程范围（132–248 m）比ALS（约136–225 m）更宽。3.2 地貌参数：各模型的高程面积曲线总体趋势一致，但卫星模型平滑了沟谷陡坡并压缩高坡度值范围，TIN保留更多细节但存在局部伪影；坡向图在沟壑网络内亦出现与分辨率相关的概括化。3.3 水文参数：卫星模型提取的排水网络总长仅约11–13 km，而ALS 5 m模型为65–68 km，说明25–30 m分辨率无法可靠表达短小排水线和沟壑轴线。自动流域边界面积分别为ASTER 11.94 km²、SRTM 12.70 km²、EU-DEM 11.98 km²、TIN 17.60 km²、ALS 13.30 km²，而参考流域约17.0 km²；多数DEM低估，TIN略高估。因采用固定1000单元阈值，网络长度差异既反映地形概括，也反映阈值对应的实际面积不同。3.4 泥沙连通性指数：所有模型均显示谷底和沟坡连通性较高，但卫星模型概括了局部地形和排水网络，TIN与ALS参考的格局最接近。由于各DEM共用同一C因子，SCI差异主要反映地形表达的影响，而非土地覆盖变化。

讨论与结论：讨论表明，DEM质量的影响在分析链各层次并不均一：整体高程分布较稳定，而依赖局部梯度和地形连续性的产品（坡度、径流组织、排水网络）对数据质量和分辨率更敏感。误差在森林区最大，说明植被冠层影响全球高程产品；但该研究未评估FABDEM等植被校正产品。黄土沟壑系统中，短历时水文事件常集中在小尺度侵蚀地貌内，局部地形信息的平滑或缺失可能不成比例地影响径流路径和坡面—沟道连接的表达。研究的主要价值在于跨平台比较，证明仅验证高程不足以判断DEM的水文适用性。研究结论可译为：对地观测产品在水文学中的适用性评估不应止步于原始数据验证；必须对与所模拟过程直接相关的派生产品进行验证，以判断输入数据的不确定性经过后续建模环节后是否仍处于可接受范围。