没有合适的资源?快使用搜索试试~ 我知道了~
青藏高原西北缘时序InSAR对流层延迟改正方法评估.docx
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 61 浏览量
2022-11-29
17:43:23
上传
评论
收藏 967KB DOCX 举报
温馨提示
试读
17页
青藏高原西北缘时序InSAR对流层延迟改正方法评估.docx
资源推荐
资源详情
资源评论
时序合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)技术
是测量地表形变的重要手段,可以用来获取大面积缓慢地表形变时变信息,在众多地球科
学研究领域中得到了广泛应用
[1-2]
。然而,时序 InSAR 技术的精度受多种误差因素影响,其
中,大气延迟误差校正是时序 InSAR 技术面临的最大挑战
[3]
。卫星发射的微波信号在穿过
大气介质时,受大气水汽含量、温度和压力等条件的影响会发生不同程度的相位延迟,即
大气延迟效应,通常包括电离层延迟和对流层延迟
[4]
。由于 C 波段波长较短且研究区域位
于中纬度,电离层延迟影响较小
[5]
,本文仅研究对流层延迟。
对流层大气参数(如水汽含量、温度和压力)的时空变化可引起高达几十厘米量级的
InSAR 延迟信号
[6]
,微小尺度的构造运动引起的地表位移可能会被掩盖。已有研究表明,
阿尔金断裂滑动速率自西向东沿断层走向由 15 mm/a 逐渐减小至 1 mm/a
[7]
,绝大部分区
域左旋走滑速率小于 10 mm/a,表示 10 a 间隔的地表位移不足 10 cm,投影到 InSAR 视
线向(line of sight,LOS)则更小,比大气延迟低一个量级
[8]
。
对流层延迟包含垂直分层延迟和湍流混合延迟两部分,大气折射率沿垂直方向的变化
引起垂直分层延迟(本文将大气等效成无数个水平薄层,假设同一层折射率相同),而大
气对流过程会引起湍流混合延迟(大气对流导致大气折射率在不同时刻的空间分布不
同)。目前,国内外很多专家学者对大气对流层延迟改正方法展开研究,所得成果可以分
为 3 类:(1)主要通过叠加、平均或时空滤波等方法去除具有时空随机分布性质的湍流
混合延迟。文献[9]发现相位叠加技术能够有效降低 InSAR 数据的数字高程模型(digital
elevation model,DEM)残差、大气相位及轨道误差的影响。根据大气信号在时间上高
频、空间上低频的特性,时序 InSAR 通过时空滤波即可减少大气延迟的影响
[10]
。(2)利
用气象模型、多光谱观测或与其他外部数据相结合的方法估计对流层延迟。文献[11]利用
动态的气象模型估计和去除大气水汽引起的信号延迟,并与 GPS 的反演结果进行比较,
结果表明两者测量的长波长空间差分结果吻合。搭载在卫星平台的多光谱仪器也可以用于
评估大气对干涉相位的影响,文献[12]验证了中分辨率成像光谱仪(medium resolution
imaging spectrometer,MERIS)数据集校正大气水汽引起相位延迟的能力。此外,文献
[13]利用 MERIS 多光谱数据与中尺度气象模型 MM5 联合估计对流层延迟相位,有效地减
少了大气效应的影响。近年来,欧洲中期天气预报中心(European Centre for Medium-
Range Weather Forecasts,ECMWF)提供的全球中分辨率数值大气再分析产品也被广泛
应用,如文献[14]利用 ECMWF 再分析资料(ECMWF reanalysis,ERA)-Interim(ERA-
I)数据集来估计对流层延迟
[14]
,在后续研究中又利用多种全球大气模型进行验证
[15]
。文献
[16]结合 GPS 观测与高分辨率 ECMWF 资料开发了一种全球覆盖、近实时的通用 InSAR
大气改正在线服务(generic atmospheric correction online service for InSAR,GACOS)
[16]
,目前得到了广泛应用。(3)根据地形与相位间的相关关系进行改正的经验模型改正
法,如线性模型或幂律模型改正。文献[17]通过相位/高程相关分析,对每幅干涉图的对流
层延迟相位进行估计和校正,利用 ERA40 全球大气再分析数据验证了校正结果。文献[18]
提出了一种幂律校正方法,该方法考虑了大气性质的空间变化特性,通过幂律关系描述了
相位延迟随高度的变化。
由于每种改正方法都存在各自的局限性,不同改正方法对垂直分层和湍流混合两个延
迟部分的敏感程度不同(例如,基于气象模型的改正方法通常无法很好地估计湍流效应引
起的延迟分量
[19]
),难以做到在任何区域和时间均能有效削弱对流层延迟的影响。本文以
青藏高原西北缘为研究区域,基于小基线集(small baseline subset,SBAS)原理,利用
哨兵 1 号(Sentinel-1)宽幅影像数据进行时序 InSAR 处理,考虑到研究区域范围广、高
程差异大、气象站点分布稀疏,选择了常用的经验模型线性改正、空间分辨率最高的
GACOS 改正和时间分辨率最高的 ERA5 数据集改正 3 种对流层延迟改正方法,讨论它们
在该研究区域形变速率反演中的效果和适用性。
1. 研究区域与数据源
青藏高原位于亚洲中部,是世界上最高最大的造山带,其北部存在一条规模巨大的左
旋走滑断裂——阿尔金断裂带,该断裂带长达 2 000 km
[20]
,在印度-欧亚板块碰撞中占据
独特的地位,它现今的运动特征对识别青藏高原变形模式具有独特的意义。本文选取了 3
个位于青藏高原西北缘的研究区域(D136、D165、D19),其中 D136 位于
76.5°E~79.7°E,D165 位于 80.5°E~83.7°E,D19 位于 84.9°E~88.1°E。
如图 1 所示,北东东向的阿尔金断裂带将青藏高原与塔里木盆地分隔开来,地形起
伏强烈,高程差异可达 5 000 m。自北向南,随着海拔增高,塔里木暖温带荒漠逐渐过渡
为高山荒漠,在海拔 5 km 以上出现永久冰川。印度-欧亚板块之间的相对运动被青藏高原
内部的地壳缩短和走滑剪切调整吸收,北部的柴达木盆地、祁连山脉和阿尔金山脉吸收
15%~17%
[21]
。尽管阿尔金断裂带构造活动明显,但历史记录到的强震(Mw > 7)的数量
很少,2017—2019 年研究区域内未发生过强震,研究区域及周边的地震分布情况见图
1。
图 1 研究区域与研究数据覆盖范围
Figure 1. Research Area and Data Coverage
下载: 全尺寸图片 幻灯片
本文研究数据采用 Sentinel-1 宽幅(interferometric wide swath,IW)降轨数据,基
于 GAMMA 平台对 D136、D165、D19 轨道近 3 a(2017—2019 年)的单视复数(single
look complex,SLC)影像数据进行差分干涉处理,采用 90 m 分辨率的 DEM 去除差分干
涉图中的地形相位,利用 MintPy(Miami InSAR time-series software in Python)程序包
进行时序分析并去除对流层延迟相位,所使用 Sentinel-1 数据的时空基线如图 2 所示。
图 2 时空基线网络图
Figure 2. Spatial-Temporal Baseline Networks
下载: 全尺寸图片 幻灯片
气象模型数据为 ECMWF 最新发布的全球天气再分析数据集 ERA5 和 GACOS。
ERA5 数据集融合了地面观测、探空气球、卫星等气象观测数据,将地表到 80 km 高度的
大气分为 137 层,网格空间分辨率为 30 km(约 0.25°×0.25°),时间分辨率精确到 1 h,
可以提供 1979 年至今每小时的大气、陆地和海洋气候变化数据。相比于 ECMWF 的前一
代产品 ERA-I 数据集,ERA5 具有更高的时间分辨率和空间分辨率
[22]
,使得气象模型估计
更接近于真实的大气情况。GACOS 提供了 SAR 影像获取时间对应的对流层延迟数据,其
数据源包含最高分辨率的 ECMWF 数值天气模型(水平分辨率 0.1°×0.1°,时间分辨率 6
h)、GNSS 对流层延迟产品、空间分辨率 90 m 的航天飞机雷达地形测绘使命(shuttle
radar topography mission,SRTM)与先进星载热发射和反射辐射仪全球数字高程模型
(advanced spaceborne thermal emission and reflection radiometer global digital
elevation model,ASTER GDEM)地形数据。
2. 研究方法
2.1 时序 InSAR 技术
时序 InSAR 技术的出现大大提高了 InSAR 地表观测的精度,目前在时序 InSAR 技术
上发展了很多软件与算法,如永久散射体差分干涉测量(permanent scatterer InSAR,
PS-InSAR)
[23]
、SBAS
[24]
、斯坦福永久散射体干涉(Stanford method for persistent
scatterers,StaMPS)
[25]
等,不同方法具有各自的优缺点,应根据实际条件(研究区域变
形特点、误差分布情况等)选出最优方法进行应用。由于本文研究区域位于青藏高原西北
缘,地形差异明显,部分区域存在冰雪覆盖,相干性受到影响,因此采用基于 SBAS 原理
的时序处理算法进行 InSAR 时序分析。
基于 SBAS 的时序 InSAR 技术,每个像素点的相位 φ 可以表示为:
ϕ=ϕdis+ϕunw+ϕatm+ϕorb+ϕdem+ϕnoise ]]>
式中,ϕdisϕdis、ϕunwϕunw、ϕatmϕatm、ϕorbϕorb、ϕdemϕdem 和 ϕnoiseϕnoise
分别为地表位移相位、解缠误差相位、大气延迟相位、轨道误差相位、DEM 相位和其他噪
声相位。数据处理时采用统计费用流方法(statistical⁃cost network⁃flow algorithm for
phase unwrapping,SNAPHU)
[26]
进行解缠,在解缠过程中进行了解缠误差的校正;在卫
星数据处理时利用精密轨道数据去除轨道误差;由于 DEM 数据自身带有误差,在数据处
理中不可避免地会将这些误差带入差分干涉图中,假定研究区域的变形为线性变化,常规
的 SBAS 改正方法可获取可靠的高精度 InSAR 时序结果;大气误差主要考虑了对流层延
迟的影响,本文采用地形线性相关的经验方法、GACOS 改正和 ERA5 数据集改正 3 种方
法校正对流层延迟相位。
对 Sentinel-1 数据提供的 SAR 影像差分干涉成图并进行时间序列分析,评估 3 种对
流层延迟校正方法的性能,技术流程如图 3 所示。时序分析基于 MintPy 实现,该软件是
InSAR 时序分析的开源程序包,通过读取解缠后的干涉图信息来生成 LOS 的地表形变信
息
[27]
。
剩余16页未读,继续阅读
资源评论
罗伯特之技术屋
- 粉丝: 3770
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- IMG_9692.PNG
- 基于SpringBoot+Vue 的网上书城项目源码+数据库脚本(毕业设计)
- 基于springboot+vue的网上书店系统源码+数据库+使用文档(高分毕设)
- 前端学习路线推荐+自己做的一个前端项目(旅游线路推广)+通过购买腾讯云服务器上传的个人经历(附带部分关键源码)
- 制造业IT规划咨询汇报方案
- 开车不犯困100首.之1~10,DJ
- Keil 找不到编译器 Missing:Complier Version5 的解决方法
- 00、3500英语词汇(英译汉+汉译英+单词表默写)(Excel表格+PDF双版本).rar
- linux发布第一篇视频解析
- 基于SpringBoot和Vue的网上书店系统源码+数据库(高分毕设)
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功