利用时空表面模型和鲁棒的自适应卡尔曼滤波器来表征局部地表动态

《Geomorphology》:Characterizing local land-surface dynamics with spatio-temporal surface models and robust adaptive Kalman filters

【字体: 时间:2026年02月28日 来源:Geomorphology 3.3

编辑推荐:

  在地形变化检测中,本研究通过改进的自适应卡尔曼滤波方法,实现了对流域地形几何特征的空间分布分析,并建立了量化不确定性的概率阈值模型,有效区分了真实地表位移与噪声干扰。

  
丹尼尔·R·纽曼(Daniel R. Newman)| 浅川裕一(Yuichi S. Hayakawa)
日本学术振兴会(JSPS)国际研究学者

摘要

过程与形态之间的关系是地貌学研究的基础。低成本遥感平台和传感器提高了收集密集时空地形观测数据的能力,从而可以实证观察地表变化。在过去的几十年中,人们开发出了地貌变化检测(Geomorphic Change Detection,简称GCD)方法,以解决由于分析高程时间序列而产生的不确定性问题。本研究旨在通过创建一组指数来扩展GCD方法,这些指数能够表征地表几何变化,改进局部不确定性估计,并将不确定性纳入指数中。为此,应用了二维自适应卡尔曼滤波器(robust adaptive Kalman filters)来估计从共震滑坡中恢复的流域的地表及其变化系数的空间分布。研究结果表明,侵蚀作用具有不对称性,更倾向于北坡;同时,在一个子流域中观察到了沟壑源头的迁移现象。与简单的减法相比,基于概率加权的垂直位移空间分布估计显示出了更好的信噪比,符号上的不一致性也更小。此外,还利用模型不确定性估计来评估模型性能。总体而言,这些方法上的进步扩展了对地表变化的分析范围,使其类似于对地形表面的地貌测量分析,并且更重要的是,提高了所有衍生结果的可靠性。

引言

作用于地球表面的地貌过程会改变地表的位置并定义其几何形态。研究这些地表形态模式是地貌学的主要关注点,因为地形是潜在过程的可观测反映(Urban, 2013)。数字数据的日益普及以及各种计算机硬件和软件的发展,直接提升了进行地貌学研究的能力(Otto et al., 2018)。事实上,这些数字技术催生了多个与地貌学相关的交叉研究领域。例如,地貌测量学(geomorphometry)是一门定量地表分析的科学,通常使用数字地形模型(Pike et al., 2009)。这导致了多种定义明确的地形指数,这些指数能够通过分析数字高程模型(Digital Elevation Models,简称DEM)来表征局部地表几何形态或地表相互作用,如日照或水流方向(Florinsky, 2017)。
几十年来,全球研究工作的重点一直是能够以越来越精细的分辨率对大面积区域进行建模和分析(例如,Miliaresis and Argialas, 1999, Meadows et al., 2024, Guth et al., 2024)。最近,无人机系统(Unmanned Aerial Systems,简称UAS)平台及多光谱相机、光探测与测距(Light Detection and Ranging,简称LiDAR)和干涉合成孔径雷达(Interferometric Synthetic Aperture Radar,简称InSAR)等传感器的普及,使得能够在较短的时间间隔和中等空间尺度上进行成本效益高的地形测量(Nex and Remondino, 2014, Cook, 2017)。这些技术非常适合监测动态过程,如海岸过程(Gon?alves and Henriques, 2015)、河流过程(Cook, 2017)或自然灾害(Giordan et al., 2018)。然而,无论是地形测量还是DEM都不是时间模型,因此通常通过减去不同时间收集的两个DEM来纳入时间维度(James et al., 2012)。由此得到的差异DEM(Difference DEM,简称DoD)量化了时间间隔内地表位置的逐像素垂直变化(Williams, 2012)。这种操作需要谨慎处理,因为两种来源的误差都会传递到结果中(Wheaton et al., 2010),从而实际上使不确定性加倍。这种累积的不确定性削弱了对结果的信心,因为真实的地表变化与噪声难以区分。
最常见的解决方案是通过定义最低细节水平来进行概率阈值处理(Wheaton et al., 2010),这种方法有时被称为地貌变化检测(Geomorphic Change Detection,简称GCD)方法(Balaguer-Puig et al., 2018)。GCD方法需要估计不确定性场,并设定一个置信水平,低于该水平的差异将被视为无法与噪声区分开来。然而,这引入了一个非平凡的挑战,即如何估计不确定性场,而不确定性场会因地形、传感器类型、伪影的存在、系统误差和遗留数据问题而有所不同(Wheaton et al., 2010, Schaffrath et al., 2015, James et al., 2020)。一些最近的方法试图通过使用方差-协方差矩阵来改进不确定性估计(Winiwarter et al., 2021, Newman and Hayakawa, 2025b)。实际上,局部不确定性模型允许使用局部自适应的鲁棒函数来处理伪影和误差,当不确定性被低估时可以相应增加(Newman and Hayakawa, 2025b)。
然而,目前的GCD方法仅限于分析垂直位移。因此,本研究的目的是扩展GCD的概念,以分析地表变化的局部几何特性,并量化可观测地表动态的时空特征。这通过以下方式实现:(1)扩展常见的地表变量以模拟局部位移场的几何特性;(2)改进局部不确定性估计和处理;(3)将概率阈值方法的原则应用于这些几何特性。第一个目标是通过基于局部地表运动学开发一个时空地表模型来实现的。然后,开发了几种指数来利用局部时空关系表征局部地表变化。Newman和Hayakawa(2025b)开发的自适应卡尔曼滤波器被应用于空间分布阵列中,以递归估计局部时空地表模型的系数并量化不确定性。最后,开发了一个具有灵活置信水平的概率加权函数,以纳入所有地表模型系数的不确定性。同时考虑了对多尺度分析的影响。

研究地点

2018年东部伊布里地震(Eastern Iburi Earthquake)是一次震级为Mw 6.6的地震,发生于2018年9月6日(YYYY-MM-DD),地点位于日本北海道(震中坐标为46.691°N, 142.007°E,深度为37.0公里),导致440平方公里区域内发生了6117次滑坡(Kasai and Yamada, 2019, Osanai et al., 2019)。虽然有些滑坡规模较大且深度较深,但大多数滑坡发生在坡度在20°至40°之间的丘陵地带(Osanai et al., 2019)。一些低矮的山脊将小流域分隔开来。

时空地形表面模型

简化的地形表面被定义为嵌入三维欧几里得空间中的二维封闭定向流形(Florinsky, 2017)。此外,地形表面由一个单值连续双变量函数定义,该函数位于一个重力场中,其范围小于平均行星半径的0.1倍(Evans, 1980, Florinsky, 2017, Shary, 1995)。因此,地表的高度是方向坐标位置的函数。

结果

图6展示了2022年8月25日至2022年11月8日期间时空梯度、方向和梯度曲率的空间分布,其中右侧列显示了一个小流域内侵蚀发生的详细图像。图6A1和A2显示了梯度的方向,梯度的大小不同颜色表示:明亮的黄色表示快速变化,洋红色表示朝北,品红色表示向东,青色表示向南,绿色表示向西。

讨论

图6和图8中呈现的时空梯度、方向和曲率信息通过量化基本几何特性增强了垂直位移信息的解释能力。图6A2中的梯度和方向分布显示,朝北的沟壑边缘变化率高于朝南的沟壑边缘。虽然高程剖面没有显示出明显的侵蚀事件,但它们显示出了持续的变化趋势。

结论

地球表面不断变化,反映了各种过程。检测和建模这些地表变化的能力为研究这些过程提供了信息支持。本研究旨在提高在不确定性存在的情况下估计局部不确定性的能力,估计地表模型,并分析地表模型以表征与观测到的地表变化相关的时空和几何特性。

CRediT作者贡献声明

丹尼尔·R·纽曼(Daniel R. Newman):负责撰写初稿、可视化处理、验证、软件开发、项目管理、方法论设计、数据收集、正式分析以及概念构思。浅川裕一(Yuichi S. Hayakawa):负责审稿与编辑、监督工作、资源协调、项目管理、资金筹集以及数据整理。

利益冲突声明

作者声明他们没有已知的财务利益或个人关系可能影响本文的研究结果。

致谢

本研究得到了日本学术振兴会(JSPS KAKENHI)的资助,资助编号包括:JP23KF0180JP23K20541JP23K23639
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普

知名企业招聘

热点排行

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

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号