冰川冻土, 2021, 43(4): 964-975 doi: 10.7522/j.issn.1000-0240.2021.0033

第二次青藏高原综合科学考察研究

InSAR技术在多年冻土区形变监测的应用

刘世博,1,3, 赵林,2, 汪凌霄2, 邹德富1, 周华云1,3, 谢昌卫1, 乔永平1, 岳广阳1, 史健宗1

1.中国科学院 西北生态环境资源研究院 冰冻圈科学国家重点实验室 藏北高原冰冻圈特殊环境与灾害国家野外科学观测研究站,甘肃 兰州 730000

2.南京信息工程大学 地理科学学院,江苏 南京 210044

3.中国科学院大学,北京 100049

Application of InSAR technology to monitor deformation in permafrost areas

LIU Shibo,1,3, ZHAO Lin,2, WANG Lingxiao2, ZOU Defu1, ZHOU Huayun1,3, XIE Changwei1, QIAO Yongping1, YUE Guangyang1, SHI Jianzong1

1.Cryosphere Research Station on Qinghai-Xizang Plateau,State Key Laboratory of Cryospheric Science,Northwest Institute of Eco-Environment and Resources,Chinese Academy of Sciences,Lanzhou 730000,China

2.School of Geographical Sciences,Nanjing University of Information Science & Technology,Nanjing 210044,China

3.University of Chinese Academy of Sciences,Beijing 100049,China

通讯作者: 赵林,教授,主要从事多年冻土与环境研究. E-mail: linzhao@lzb.ac.cn

编委: 武俊杰

收稿日期: 2020-05-13   修回日期: 2020-11-12  

基金资助: 第二次青藏高原综合科学考察研究项目.  2019QZKK0201
国家自然科学基金项目.  41931180

Received: 2020-05-13   Revised: 2020-11-12  

作者简介 About authors

刘世博,博士研究生,主要从事多年冻土形变研究.E-mail:liushibo18@mails.ucas.ac.cn , E-mail:liushibo18@mails.ucas.ac.cn

摘要

作为大地测量的一种新兴空间技术,合成孔径雷达干涉(synthetic aperture radar interferometry, InSAR)具有全天时、高精度、大范围和速度快的优点,逐渐被应用于多年冻土区地表形变监测中。通过综述多年冻土形变原理及InSAR监测多年冻土形变的应用实例,研究表明:在气候变暖的背景下,多年冻土区地表年际形变以下沉为主,多年冻土上限附近地下冰含量的大小是影响年际形变量的主要因素;活动层内土壤含水量影响着地表季节形变量的大小,不同类型多年冻土区的地表年际形变量和季节形变量存在着较大的差异。研究还表明,不同波长的SAR产品在不同类型多年冻土区的适用性不同,下垫面特征对利用InSAR获取地表形变量有较大影响,L波段的SAR数据在植被覆盖度较好的区域有更好的效果。由于InSAR的失相干问题,加之目前还缺少长时间、多类型、高频率的实测形变结果作为验证和标校数据,获取准确且连续的大范围形变数据较为困难。针对目前寒区研究需求,布设野外长期观测站点,建立适用于不同多年冻土区的地表形变反演算法,构建具有较高精度和较高时空分辨率的地表形变数据集具有重要的实践和科学意义。

关键词: InSAR ; 多年冻土 ; 地表形变 ; 冻胀 ; 融沉

Abstract

As an emerging space technology in geodesy, synthetic aperture radar interferometry (InSAR) has the advantages of all-day, high-precision, large range and fast speed. It has a great potential for application in surface deformation measurement, and is increasingly used in the monitoring of permafrost areas. In this paper, the principle of permafrost deformation and the application of InSAR to monitor permafrost deformation are reviewed. The results show that the interannual deformation of the permafrost is mainly settlement in the context of climate warming, and the amount of underground ice near the upper limit of permafrost is the main factor affecting the interannual variation. Soil moisture content within the active layer affects the magnitude of the seasonal deformation variables. There are large differences in interannual and seasonal surface deformation variables for different types of permafrost areas. It was also shown that the applicability of SAR products at different wavelengths varies for different types of permafrost zones, and that subsurface features have a greater influence on the magnitude of deformation obtained by InSAR data, while the L-band SAR data are more effective in the areas with better vegetation cover. Owing to the decoherence problem of InSAR processing and the lack of validation and calibration data for long-term, multi-type, high-frequency permafrost deformation monitoring, it is of great practical and scientific significance to establish surface deformation inversion algorithms for different permafrost areas and to construct surface deformation data sets with high precision spatial and temporal resolution.

Keywords: InSAR ; permafrost ; surface deformation ; frost heave ; thaw settlement

PDF (1944KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘世博, 赵林, 汪凌霄, 邹德富, 周华云, 谢昌卫, 乔永平, 岳广阳, 史健宗. InSAR技术在多年冻土区形变监测的应用[J]. 冰川冻土, 2021, 43(4): 964-975 doi:10.7522/j.issn.1000-0240.2021.0033

LIU Shibo, ZHAO Lin, WANG Lingxiao, ZOU Defu, ZHOU Huayun, XIE Changwei, QIAO Yongping, YUE Guangyang, SHI Jianzong. Application of InSAR technology to monitor deformation in permafrost areas[J]. Journal of Glaciology and Geocryology, 2021, 43(4): 964-975 doi:10.7522/j.issn.1000-0240.2021.0033

0 引言

多年冻土主要分布于环北极的高纬度地区和中低纬度的高海拔地区,面积约为23×106 km2,占北半球陆地面积的25%1-3。岩土层中冰的存在是多年冻土的主要特征,由于冰的密度约比水小9%4-5,故随着活动层中水分的季节冻结和融化,地表会发生周期性冻胀和融沉。气候变暖会导致活动层融化深度增大,多年冻土上限附近地下冰融化,致使地表出现下沉的趋势6-8。地表的冻胀融沉对寒区工程建筑的影响一直是冻土学的主要研究方向之一。在过去几十年,由于人为活动的干扰和气候变化的共同作用,地表形变已经导致青藏公路、东北铁路的破坏率达到30%以上9-12,严重危害了寒区的交通运输安全。因此长期且连续的地表形变监测不仅可以获取多年冻土对气候变化的响应,还可以对多年冻土区工程的建设与维护提供宝贵的信息。

地表形变的传统监测手段主要以埋设仪器、水准测量、差分GPS(global positioning system)为主13-16。在兰州的马衔山17,通过多次水准观测并结合温度与土壤水分数据,发现该地区地表10月开始冻胀,4月开始融沉,形变的较大值一般出现在冻土发育良好的小湖滩附近。水分迁移是导致该地区出现冻胀的主要因素。水准测量是监测多年冻土区道路形变的主要手段,在我国的青藏公路和哈大高速铁路都对不同结构的路基进行长期水准监测18-19,获得的形变信息为道路维护提供帮助。差分GPS在监测多年冻土地表形变也具有极大的潜力,在阿拉斯加中北部地区利用DGPS技术获取了该地区2001—2002年多年冻土地表形变状况。相比较于水准测量,内置的高精度大地坐标系统和快速连续的测量是其主要优势16。在西伯利亚勒拿河流域,将树脂玻璃板用金属管水平固定在地表(可由地表垂直移动),通过反复测量深入多年冻土内部的钢管顶部与地面上的有机玻璃板之间的距离来测量地表形变,作为InSAR结果的验证数据20。上述传统方法可以获取连续准确的形变信息,但不能满足大面积、大范围、长时间序列的形变监测需求21

随着InSAR技术的不断发展22-25,利用其获取多年冻土区地表形变逐渐成为研究热点。其中D-InSAR技术在加拿大北部26、阿拉斯加27和青藏高原28等多年冻土区都得到了较为准确的多年冻土地表形变信息。但D-InSAR技术的监测精度受到时空失相干和大气延迟的影响29,且无法获取研究区域地表形变的时间序列。为了弥补D-InSAR的不足,Berardino等30于2002年提出时间序列InSAR的小基线技术,可以获取长时间连续准确的多年冻土地表形变信息。InSAR时间序列方法目前主要分为三种:①永久散射体技术(persistent scatter, PS)31-32,此方法在生成差分干涉相位后仅对高相干的目标进行建模和形变解算;②小基线技术(small baseline subset, SBAS)3033,通过识别单元内观测目标的相位和强度均值作为此像素观测量,处理多个干涉对来获取各影像获取时刻相对于初始时刻的形变量;③临时相干点技术(temporal coherent point, TCP)34,相比于PS技术,此方法增加了观测目标密度,不只是选取高相干目标。还有将上述三种方法相互组合进行形变计算,例如StaMPS算法结合PS和SBAS技术,将PS点作为观测目标,计算各目标的形变时序解,在形变监测中也有着较多的应用35

近年来,许多学者利用InSAR技术实现了对多年冻土区地表形变监测的研究,获取了较为可靠的地表形变结果,证实了该技术在多年冻土区探测地表形变的能力36。Liu等37利用InSAR技术获取了阿拉斯加多年冻土区北坡地表的形变信息,结果显示在1992—2000年间由于多年冻土上限地下冰融化而导致该地区发生了1~4 cm的沉降。这一结果也很好的解释了该地区多年冻土温度升高,活动层厚度却没有增厚的原因。在青藏高原多年冻土区,利用InSAR技术监测多年冻土的形变量为公路铁路的路基维护提供了快捷准确的信息2838,节省了大量的人力物力。利用InSAR技术得到的准确、连续、大范围的地表形变信息不仅可以使人们更好地了解气候变暖下多年冻土的变化,而且随着研究不断的深入,根据形变结果模拟活动层厚度变化、地下冰融化状态等方面有着很好的应用39-41

1 InSAR及多年冻土形变原理

1.1 多年冻土形变原理

多年冻土地表季节形变主要由活动层中水分相变引起地表冻胀融沉所致。如图1所示,活动层的冻融过程可分为四个阶段,分别为夏季融化、秋季冻结、冬季降温和春季升温9。夏季融化阶段活动层处于吸热状态,热量传输由上向下,融化锋面逐渐向下迁移,其附近的自由水在重力作用下迁移至0 ℃附近,由此导致的融水再冻结使得在融化初期地表有可能出现轻微的抬升。随着活动层不断吸收热量,地表的融沉量逐渐增大。在秋季冻结阶段,活动层融化至最大深度后由底部向上开始冻结,而此时活动层上部白天吸收热量,晚上放出热量,此时地表处于较为稳定的状态。随着温度的降低,活动层进入双向冻结过程,也就是进入“零幕层”阶段,此时地表形成稳定冻结层,融化层上部冻结锋面快速向下移动,水分不断向冻结锋面迁移和冻结,此时地表多为抬升的趋势。随着冻结锋面向下迁移速度变缓,活动层进入相对稳定的冻结阶段。而冬季降温过程是冰分凝过程与冷缩过程共同作用,地表以冻胀为主。到了春季升温阶段,活动层温度升高,地表会出现轻微的白天融化下沉,夜晚冻结抬升的情况。

图1

图1   活动层季节冻融过程

Fig.1   Seasonal freeze-thaw process of active layer


图2所示,多年冻土地表的年际形变主要是由其上限附近的地下冰融化导致。研究表明,多年冻土上限附近的地下冰是由于重复分凝作用造成的,具体地说,它是由正冻土中的成冰作用、已冻土中的成冰作用、正融土中的成冰作用、未冻水的不等量迁移规律、冰的自净作用和冰的共生生长等多种作用的综合并重复发生的结果42

图2

图2   多年冻土年际形变原理及趋势

Fig.2   Principles and trends of interannual deformation for permafrost


当活动层中存在温度梯度,就会产生自由能梯度,从而引起水分向温度低的方向移动。当多年冻土发生自上而下冻结时,水分向多年冻土上限聚集,形成地下冰。这种水分迁移不仅在正冻土中出现,在已冻土和正融土中同样存在。多年冻土中还存在未冻水,当温度梯度为正值时,未冻水将发生向上的迁移。由于这种迁移发生在非剧烈相变区,故其迁移量较小;而当温度梯度为负值时,水分将向下迁移,且发生在剧烈相变区,故其迁移量较大。夏季由负梯度造成的向下的迁移水量显然要大于冬季由正梯度造成的向上迁移的未冻水量,因此水分的向下迁移并在多年冻土上部积聚,使其含冰量增加43。同时冰的自净作用使未冻水向下积聚而土颗粒向上积聚,这种作用的逐年重复使含冰量继续提高。最后,冰的共生生长作用致使多年冻土上限附近含冰量进一步增加。

在气候变暖的背景下多年冻土上部活动层增厚,温度升高44,其上限厚层地下冰逐渐融化,地表在长期状态下多呈现下沉趋势。多年冻土地表形变是季节形变和年际形变共同作用的结果,活动层土壤含水量和多年冻土上限地下冰融化速度决定着地表的形变状态。多年冻土的自上而下融化可能引起地表大面积的沉降10,致使多年冻土土壤基底承载强度降低,对多年冻土区工程的稳定性造成较为严重的影响45-46

1.2 InSAR技术原理

干涉测量的本质就是依靠地面目标两次后向散射回波信号的高相干性,求取两次信号在传播过程中的相位差ϕ。如图3所示,假设卫星两次通过地面点P,而卫星过境时的位置分别为A1A2,目标点P到卫星传感器的斜距分别为RR+∆R。图中z为P相对于参考椭球的高程。A1A2之间的距离B被称为基线,α为基线与水平面的夹角。由于卫星的高度在运算时已知,假设A1到参考椭球的高程为H。由图中几何关系可得,如果求解P的高程,首先要解决卫星在A1位置时的入射角θ。信号在两次传播过程中的相位差ϕ可以表示为

图3

图3   InSAR获取地表形变原理

Fig.3   Principle of InSAR to obtain the surface deformation


ϕ=ϕ1-ϕ2=- 4πλR-R+R=4πλR

式中:λ为雷达发射信号的波长。

图3中的几何关系可得

sinθ-α=R2-(R+R)2+B22RB=- R(2R+R)2RB+B2R

由于RBR≫∆R,将式(1)代入式(2)得

sinθ-α- RB=- λϕ4πB

如果已知ϕαBλ以及H,即可计算P点的高程z

θ=α-arcsin(λϕ4πB)
z=H-Rcosθ

式(1)可知波长是计算相位差的重要参数,而SAR影像是卫星发射的电磁波与地物目标相互作用的结果,不同波长的电磁波与地物目标相互作用会有不同的表现。SAR数据的波长与卫星雷达天线的设计密切相关,从图4可以看出,星载SAR卫星主要以C、L和X波段为主。不同波长的数据对于地表的穿透能力也不同,因此在多年冻土区的适用性也有所差异。

图4

图4   不同波段SAR产品

Fig.4   SAR products for different bands


在加拿大育空地区的赫歇尔岛,Short等47针对TerraSAR-X(X波段),RADARSAT-2(C波段)和ALOS-PALSAR(L波段)三种不同波长的SAR产品在连续多年冻土区监测地表形变的能力进行比较。TerraSAR-X较高的时空分辨率,在数据处理过程中降低了失相干现象出现的可能,可以更好地监测一些突发事件,如滑塌造成的地表形变的动态变化。而较高的空间分辨率在监测道路形变等较小目标物上有着很好地应用。但较短的波长更容易受到空气中水汽的干扰,较差的穿透力在植被覆盖度高的区域往往会受到干扰。TerraSAR-X数据目前处于收费阶段,尽管数据时间较长,但获取仍较为困难。而波长较长的L波段数据具有更强的穿透力,在植被覆盖高、土壤含水量大和地势较为复杂的地区有着更好的应用,可以提供更为完整准确的形变结果,也具有长时间监测多年冻土形变的优势。但由于ALOS数据基线较长,误差相较于其他数据可能会更高,所以对DEM数据的准确性要求更精确。在高纬地区电离层对L波段的SAR数据影响较大,会对信号的传播带来较大的干扰。同时较长的重访周期在长时间形变监测中会漏掉很多有用的信息,尤其是在冻融初始阶段,对短期的季节形变和突发性形变监测有着较大的影响。而C波段波长介于L波段和X波段之间,兼顾了二者的优势与不足,是目前使用最为广泛的波段。目前有ERS1/2、Sentinel-1、RADARSAT等卫星,在多年冻土区形变监测方面都取得了较好的结果,在该区域的研究中RADARSAT数据的表现也是最为稳定的。

2 InSAR在多年冻土区形变监测的应用

2.1 季节形变

由于多年冻土所处的纬度、海拔、局地气候条件不同,其季节形变也有着不同的特征。图5显示目前通过InSAR技术获取的多年冻土区季节形变特征表现为,青藏高原北麓河地区48、黑龙江北黑公路沿线49和俄罗斯勒拿河三角洲地区50为高值区,地表的最大季节形变量可达100 mm。北麓河地区于3月地表达到最大的抬升量,而北黑公路较该区融化时间更早,且冻胀量相比沉降量较小,高分辨率的TerraSAR-X数据在公路形变研究中有着较好的应用,可以获得更高的形变信息。Zhang等48发现在北麓河的草甸区域发育了许多冻融湖,剧烈的冻融作用致使该地区形变量高值多出现在热融湖塘周围,Strozzi等50在勒拿河三角洲地区也发现了同样的现象。季节形变量较小的区域主要出现在青藏高原北部地区51和格陵兰岛东北部52,干旱的气候致使这些地区土壤未冻水含量较低,加之较短的暖季和较少的人类活动造成了这些地区地表季节形变量分别只有12 mm和7.4 mm。

图5

图5   多年冻土区的季节形变量

Fig.5   Seasonal deformation in permafrost areas


在青藏高原北部,形变量较大的地区主要分布在东西走向山脉之间的狭长盆地,地势平坦地区的季节形变量往往大于陡坡地区。季节形变量较大区域融化开始时间往往要比形变量较小区域早1~2个月,Daout等51认为是由于土壤水分不饱和导致水分主要集中在活动层底部,对于气温变化响应有所延迟。在南极的利文斯顿岛,InSAR结果显示在2016年12月至2017年4月期间发生了20~30 mm的冻胀50,且表现出了明显的空间异质性。Strozzi等50认为这是由于该区域靠近海洋,土壤中盐分含量差异造成了形变的不同。同时土壤沉积物颗粒差异也造成了冻胀量的不同,细颗粒地表往往具有更大的冻胀量。在加拿大的马更些河山谷,根据RADARSAT的数据结果显示在2006年7月至10月该地区发生过滑坡或火烧现象的区域的地表融沉量可达13~18 mm,而没有出现火烧现象的区域其融沉量只有3~7 mm,火烧现象明显增加了该区域的融沉量,且发生火烧现象的区域地表沉降的结束时间也要晚于其他地区,Singhroy等53认为是由于燃烧后产生的黑色腐殖质使土壤吸收更多的热量,更深的融化深度造成了更大的下沉量。综上所述,多年冻土区的地形和土壤质地通过控制水分和热量影响着季节形变量的大小。

Envisat数据显示,青藏铁路拉萨—当雄试验路段2007年6—8月期间的最大沉降量达到30 mm,从9月至次年的2月最大冻胀量可达40 mm54,从山顶至山底形变振幅总体呈增大趋势。青藏公路沿线附近区域的季节形变量更大55,在该区域的公路沿线融化季地表下沉量接近80 mm。青藏公路沿线季节形变量也有着明显的空间异质性,根据ALOS数据显示两道河地区季节形变量最大值只有20 mm,相对于拉萨当雄段公路形变量较小56。在阿拉斯加的Prudhoe Bay地区,根据ERS1/2数据37和全球差分定位系统(DGPS)16监测结果显示,该区域季节融沉量在10~40 mm之间,港口和石油开采区是该地区形变的高值区域。在格陵兰岛西部,由于较多的人类活动使其季节形变量要远远大于东部,最大值可达30 mm50。从上述几个例子可以看出,人类活动加剧了多年冻土地表形变的发生。

通过上述多年冻土区形变研究发现,多年冻土地表季节形变受到地形、土壤质地、未冻水含量和人类活动的影响,表现出明显的空间差异性,地表融沉开始和结束时间也有所不同。在欧洲的格陵兰岛东北部,该区域1995—1999年间融化季地表形变量在0.4~6.1 mm之间,而在2006—2009年间则增加到2.3~7.4 mm之间,该地区多年冻土的季节融化深度也有着明显的增加。在全球气候变暖的背景下,多年冻土区的季节形变量也有增加的趋势,对气候变暖有着很好的反馈52

2.2 年际形变

随着SAR数据的不断增加和时序InSAR技术的发展,可以获取较长时间的多年冻土形变信息。由表1可知,多年冻土年际形变的高值区主要集中在青藏高原的北麓河57-58、五道梁59地区和俄罗斯的勒拿河三角洲地区60。多年冻土上限附近的地下冰融化速率较快导致了较高的年间沉降量,五道梁地区年际沉降量可达10.28 mm。地表年际形变量较小的地区主要分布在青藏高原的西北部5161、格陵兰岛东北部52和西藏当雄县附近54,这些地区多年冻土上限附近地下冰含量往往较低,地表年际形变量多在3 mm以下51-5254。多年冻土的年际形变也有明显的空间差异性,2003—2010年45景Envisat数据显示北麓河高山地区多年冻土地表状态较为稳定,年平均形变量只有-0.75 mm。而地势较低的稀疏植被区多年冻土地表形变较大,年形变量可达-16~0 mm之间。

表1   多年冻土区的年际形变量

Table 1  Interannual deformation in permafrost areas

区域研究区年份形变量/mmSAR数据
青藏高原北麓河2007—2010-5~5ALOS
2003—2010-1.78Envisat
五道梁2014—2017-10.28Sentinel-1
青藏公路2007—2010-0.52Envisat
青藏铁路2007—2010-0.14Envisat
高原西北部2003—2011-2~3Envisat
2007—2009-2~0ALOS
北美洲阿拉斯加北部1992—2000-4~-1(每年)ERS1/2
欧洲格陵兰岛东北部1995—1999-2.4~-0.3(每年)ERS1/2
2006—2009-2.7~-0.8(每年)Envisat
西伯利亚勒拿河三角洲2016—201710Sentinel-1

新窗口打开| 下载CSV


北麓河地区铁路沿线地表下沉量为每年1.01 mm,公路附近区域年沉降量为1.47 mm,更好的保护措施和更高的通车标准使得铁路沿线年际形变量要低于公路。青藏高原的其他区域也很好地显示了出这一特征,根据2007—2010年Envisat数据显示,青藏铁路西藏当雄县沿线多年冻土年平均形变量为-0.14 mm,而公路年平均形变量为-0.52 mm,公路附近的形变量要远大于铁路。多年冻土区的交通运输对地表年际形变也会产生较大影响,研究人员认为交通运输线的铺设改变了多年冻土区下垫面性质,人类活动进一步加剧了地表的形变,因此运输沿线地区年际形变一般要高于其他地区54

在格陵兰岛东北部地区,利用ERS1/2及Envisat数据监测了该地区两个时间段的地表形变状况,结果显示在1995—1999年间该区域多年冻土区年平均形变量数值为-2.4~-0.3 mm,而2006—2009年间地表下沉量增长到0.8~2.7 mm·a-1,气候变暖导致多年冻土上限地下冰融化速度加快,年际形变量增大。阿拉斯加北部地表的长期下沉也很好地解释了该地区尽管气温升高冻土退化但活动层却没有明显增加的原因,地表形变信息相比于活动层厚度数据可以更好地反映该地区多年冻土的退化趋势。但是在气候变暖及多年冻土区地表年际普遍发生下沉的趋势下,还是有一些区域会出现一些异常的情况。例如在勒拿河三角洲地区2017年融化季地表形变相对于2016年发生了10 mm左右的上升,Chen等60根据气温统计发现这很可能是由于2017年夏季融化开始时间延迟所导致,2017年7月的平均温度要比2016年7月低0.6 ℃。多年冻土年际形变对于气候的长期变化趋势有着很好的反馈,对于气候变化研究是极有意义的。

2.3 不同下垫面的形变

下垫面情况的差异造成了多年冻土区地表形变的异质性。青藏高原现存多年冻土的总面积约为106万km2[62,由于其土壤含水量、气温、海拔、植被覆盖度的不同,其地表形变特征也有所也差异,以青藏高原唐古拉和两道河这两个观测场为例。

唐古拉观测场位于青藏高原多年冻土区腹地,由表2可知该地区多年冻土未冻水含量和植被覆盖度较低。而两道河地区位于多年冻土区南界边缘,土壤含水量和植被覆盖度均高于唐古拉地区。通过InSAR反演和水准实测都显示这两个地区出现了缓慢的长期沉降趋势,且年间形变量差距不大。而由于不同土层含水量的不同,两道河地区季节形变量明显要大于唐古拉地区。通过表2可知,在4月和7月两道河地区各层土壤质量含水量均远远高于唐古拉地区,较高的未冻水含量导致了较大的季节形变量。唐古拉的降水主要集中在活动层融化阶段,孔隙度较大的砾质砂壤土会导致降水或融水下渗到活动层融化锋面以下,再次冻结导致活动层中冰体积的增加而使得地表高程轻微增加。而同样土壤孔隙度较大的两道河地区却没有出现类似的现象,很可能由于该区域活动层土壤在融化期间其含水量较高,导致降水不会继续下渗而是转化为地表流走或者蒸发56。地形的起伏通过影响土壤含水量也对季节形变产生了影响,通过对两道河地区的观测发现在个别地势较低且平缓的观测点,两个月的形变量可达153 mm,远远高于其他地区,造成了该地区内部地表形变的异质性。

表2   唐古拉和两道河观测场植被覆盖度及土壤含水量

Table 2  Vegetation coverage and soil water content in Tanggula and Liangdaohe study sites

观测站点年平均气温/℃植被覆盖度/%深度/cm土壤孔隙度/%土壤质量含水量/%
4月7月
0~1052.52.78.6
10~2043.74.111.8
唐古拉-4.920~3020~3046.89.112.9
30~4033.16.012.8
40~6042.94.711.5
0~1071.374.777.3
10~2062.0188.6104.0
两道河-1.3>8020~3057.7164.0130.9
30~4052.193.083.9
40~6048.756.756.1

新窗口打开| 下载CSV


多年冻土区地表季节形变主要由土壤含水量决定。随着高原降雨量的增加,降水或者融水下渗对多年冻土区水文过程、感热和潜热通量变化有很大影响,对地表形变状态也会产生很大的不确定性。较高的未冻水含量往往伴随着较好的植被覆盖度,结果显示C波段数据在两道河地区结果与水准实测数据及冻融规律不太符合,较高的未冻水含量导致复杂的相位变化,会造成较大的误差。植被生长过程中的变化也会影响地表形变的反演结果,而波长较长的L波段数据在植被覆盖较高的地区较C波段数据效果更好。

2.4 活动层厚度模拟

活动层是指覆盖在多年冻土之上的夏季融化,冬季冻结的土层63。活动层是多年冻土区地层水热交换最为频繁的区域64,其厚度的变化是地气能量交换的直接结果,是影响寒区生态环境最活跃的因素。同时活动层厚度变化也直接影响多年冻土区的水文特征和植物生态环境,在多年冻土研究中有重要意义65

活动层厚度的获取方法主要以实测和模型反演为主66。实测是通过在野外埋设相关仪器记录地下土层温度信息或坑探的方式来获取多年冻土的融化和冻结深度。野外实测数据精度高,作为模型模拟的精度验证数据极为重要,但其空间分辨率受到限制,只能得到小范围内的信息。而模型模拟方法综合考虑积雪、植被、土壤含水量等多种因素,可以获得大范围的活动层厚度信息。但模型反演法原理复杂,需要考虑众多的影响因子,获取准确的活动层厚度信息较为困难。

随着InSAR技术的发展,人们已经不限于利用其获得多年冻土区地表形变信息,还根据多年冻土季节形变与活动层厚度之间的关系,计算出了较为可靠活动层厚度数据。Liu等39假设多年冻土地表季节沉降完全是由活动层水分相变所致,通过1992—2000年ERS1/2数据获取的融化季沉降信息,结合土壤含水量、孔隙度、饱和含水率和土壤水分相变导致的体积差,获取了阿拉斯加北坡8 000 km2的多年冻土最大融化深度,也就是活动层厚度。可表示为

H=ρiPS(ρw-ρi)d

式中:H为融化深度;P为土壤孔隙度;S为饱和土壤含水率;ρw为水的密度;ρi为冰的密度;d为地表沉降量。文献[39]默认该地区活动层土壤含水量饱和,即S=1。

将模拟结果与该地区CALM观测的活动层厚度数据进行比较,结果较为一致,证明了InSAR监测多年冻土活动层厚度的潜力。随后,李珊珊54、Jia等61等利用类似的反演模型并对土壤孔隙度和土壤饱和度进行了简化,反演了青藏高原地区的活动层厚度,均取得了较好的结果。但上述方法没有考虑土壤湿度的影响,或假设土壤水分处于饱和状态,这种简单的处理会引起较大的误差。张正加67在进行青藏高原活动层厚度反演时发现不同地形地貌的土壤含水量具有差异性,同时地下土壤水分随深度而变化,如果不充分考虑土壤含水量的话,会导致活动层厚度的估计结果存在误差。因此根据不同土壤类型和土壤含水量的情况,建立了适合青藏高原多年冻土区的活动层厚度反演模型,并使用两个实测点数据对北麓河地区反演的结果进行验证,误差分别为0.3 m和0.5 m。模拟结果均取得了较好的结果,弥补了气象站点监测空间性不足的缺陷。

3 讨论

InSAR技术在多年冻土区地表形变监测方面取得了较多的成果,但目前还是存在一些问题亟待解决,其中失相干现象是限制InSAR技术在多年冻土区应用的主要问题。在差分干涉测量时,由于两影像获取的卫星天线中心与观测目标的连线并不完全重合,这种几何结构差异就引起几何失相干。而两幅SAR数据时间间隔较长进而导致地面目标的介电常数变化较大时,会引起时间失相干。同时大气条件、基线长度以及轨道轻微不平行、变形运动过快等情况都会导致相位的失相干,尤其是冬季积雪会对形变反演带来极大的困难。失相干现象降低了相位结果的相干性,导致后续解缠的困难和失败,影响研究区信息的准确性与完整性。尽管近年来发展了永久散射体技术、小基线子集等方法大大缓解了失相干现象的发生,例如在2019年迈阿密大学推出最新的时序算法Mintpy68开源软件包,可以读取ISCE、ARIA、FRInGE、SNAP、GAMMA和ROI_PAC等主流InSAR软件格式数据,可以生成三维的地表形变结果。该算法通过加权最小二乘法进行时间序列反演,利用最小树桥接、相位闭合减小干涉图的积分模糊度和基于相干的网络修正等方法减小了相位解缠结果地误差且内置的PYAPS模块可以对结果进行大气校正,减小大气对时序结果造成的误差。但目前失相干情况还是较为普遍的存在,尤其是波长较短的SAR数据在含水量较高,植被覆盖的较好的地区更为明显。此外,相位解缠也是数据处理中的难点,由于地面状况的复杂性和干涉像对本身的质量差异,使得相位成功解缠的难度很大。目前提出了如瞬时频率算法(IFA)、卡尔曼滤波算法、自动抑制全局扩散误差的算法69以及SNAPHU70算法大大提高了解缠成功率,为多年冻土区地表形变监测提供了更为准确完整的信息。

实测数据的缺失对遥感反演的验证造成了很多困扰。目前很多形变反演结果并没有相应的实测数据验证,只能根据是否符合冻胀融沉规律来判断结果的准确性。但在青藏高原多年冻土区的长期监测结果发现地表形变并不一定完全符合融化季下沉,冻结季抬升的规律。例如在唐古拉地区2015年5月至7月中旬和2016年6月至10月地表出现轻微抬升,正常情况下处于融化季多年冻土区地表应处于融沉的状态。这一反常现象可能是由于活动层融化阶段降水或融水下渗到冻结锋面导致再次冻结,从而引起地表出现轻微抬升而导致。这些现象已经在实验室或者野外实验中证实,Cheng71在青藏高原多年冻土区、Mackay72在北极西海岸多年冻土区均有发现。因此,单纯以符合冻融规律作为验证标准并不一定能反映结果的准确性。

尽管InSAR技术在多年冻土区的应用仍存在一定的问题,但该技术有着极为广阔的前景,越来越多的人将目光投到了该技术对多年冻土的监测中。目前处理SAR数据软件和算法多样,生成的数据格式也不同,得到的结果数据格式和文件类型等都有较大差异,在进行时序运算时会由于格式和内容的不同造成很多麻烦。因此,数据输出格式的统一和各软件间相互支持是日后不断发展的趋势。同时,大范围、长时间、高分辨率的SAR产品在自然条件较为恶劣的多年冻土区有着较大的应用需求,地表形变产品的推出会加速多年冻土形变的研究。最后,随着研究的不断深入,InSAR不仅仅局限于对多年冻土形变监测的应用,还要不断挖掘其潜力,应用于更多的研究领域中。地表形变对于多年冻土退化73-74、全球变暖75、碳源碳汇转变76研究都有着重要的意义,值得不断深入的研究与探索。

4 结论

由于SAR数据种类的不断增多和算法的不断进步,InSAR技术在多年冻土区有着广阔的应用前景。因此,对于多年冻土区大范围连续的地表形变监测及形变原理研究迫在眉睫。通过总结前人的研究,发现:

(1) InSAR技术在多年冻土区地表形变监测方面有着较好的应用,在欧洲、中国青藏高原、俄罗斯、美国阿拉斯加、加拿大等地都获取了大范围的多年冻土形变结果,且具有明显的空间异质性。多年冻土地表季节形变量大的区域主要集中在土壤含水量大或人类活动较多的公路铁路沿线地区,而气候较为寒冷干燥土壤含水量低的地区季节形变量较小。年际形变量较大的地区主要分布在富冰多年冻土区,多年冻土上限附近的地下冰在气候变暖背景下不断融化造成了地表出现年际下沉的趋势,例如阿拉斯加北部及青藏高原北麓河地区。随着SAR卫星产品的不断增加及InSAR技术的不断进步,人们可以获取更高的时空分辨率及更加准确的多年冻土区地表形变信息。

(2) 由于多年冻土区下垫面条件有所差异,因此不同波段的SAR数据在不同区域的适用效果是不同的。波长较短的SAR产品在土壤干燥、植被覆盖度较低的地区效果较好。而波长较长的SAR产品在含水量高、植被条件较好的地区准确性更高,在地形起伏较大或融冻泥流造成较大形变的情况下会有更好的表现,但对DEM数据的精度要求更高,同时高纬度地区电离层对L波段的SAR数据影响较大。因此选取合适的SAR产品在形变研究中极为重要。

(3) 随着卫星产品的不断增加和处理技术的不断进步,InSAR在多年冻土区有着更好的应用前景,但也存在着许多问题亟待解决。失相干、相位解缠错误及大气对信号的干扰等问题都需要科学家对算法不断改进,寻求更优的方案。同时多年冻土区形变研究也需要更好的数据获取渠道及下载速度,由于一些产品处于收费的状态,影响了数据在形变监测中的应用。最重要的是,利用InSAR技术在多年冻土区监测地表形变需要掌握必要的多年冻土理论知识,对于多年冻土形变及活动层变化的机理及过程有着较为清楚地认识,这样才能不断推动这项技术在多年冻土区应用、进步和发展。

参考文献

Zhao LinWu QingbaiMarchenko S Set al.

Thermal state of permafrost and active layer in central Asia during the international polar year

[J]. Permafrost and Periglacial Processes, 2010212): 198-207.

[本文引用: 1]

Cheng GuodongZhao Lin.

The problems associated with permafrost in the development of the Qinghai-Xizang Plateau

[J]. Quaternary Sciences, 2000206): 521-531.

程国栋赵林.

青藏高原开发中的冻土问题

[J]. 第四纪研究, 2000206): 521-531.

Bockheim J G.

Permafrost distribution in the southern circumpolar region and its relation to the environment: a review and recommendations for further research

[J]. Permafrost and Periglacial Processes, 199561): 27-45.

[本文引用: 1]

IPCC. Climate change 2013: the physical science basis: contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change[M]. New YorkCambridge University Press2013.

[本文引用: 1]

Cheng G DWu T H.

Responses of permafrost to climate change and their environmental significance, Qinghai-Tibet Plateau

[J]. Journal of Geophysical Research: Earth Surface, 2007112F2): 93-104.

[本文引用: 1]

Hinzman L DKane D LGieck R Eet al.

Hydrologic and thermal-properties of the active layer in the Alaskan Arctic

[J]. Cold Regions Science and Technology, 1991192): 95-110.

[本文引用: 1]

Kane D LHinzman L DZarling J P.

Thermal response of the active layer to climatic warming in a permafrost environment

[J]. Cold Regions Science and Technology, 1991192): 111-122.

Yershov E D

Principles of permafrost[M]. Liu Jingren, trans. LanzhouLanzhou University Press2015. [

[本文引用: 1]

叶尔绍夫Э Д. 冻土学原理[M]. 刘经仁, 译. 兰州兰州大学出版社2015.]

[本文引用: 1]

Zhao LinCheng GuodongLi Shuxunet al.

Thawing and freezing processes of active layer in Wudaoliang region of Tibetan Plateau

[J]. Chinese Science Bulletin, 20004511): 1205-1211.

[本文引用: 2]

赵林程国栋李述训.

青藏高原五道梁附近多年冻土活动层冻结和融化过程

[J]. 科学通报, 20004511): 1205-1211.

[本文引用: 2]

Shiklomanov N IStreletskiy D ALittle J Det al.

Isotropic thaw subsidence in undisturbed permafrost landscapes

[J]. Geophysical Research Letters, 20134024): 6356-6361.

[本文引用: 1]

Wu QingbaiTong Changjiang.

Permafrost change and stability of Qinghai-Tibet Highway

[J]. Journal of Glaciology and Geocryology, 1995174): 350-355.

吴青柏童长江.

冻土变化与青藏公路的稳定性问题

[J]. 冰川冻土, 1995174): 350-355.

Wang ShuangjieLi ZhulongWu Jingmin.

The state-of-the-art and new subjects of highway research in permafrost regions of China

[J]. Journal of Glaciology and Geocryology, 2003254): 471-476.

[本文引用: 1]

汪双杰李祝龙武憼民.

多年冻土地区公路筑路技术研究现状与新课题

[J]. 冰川冻土, 2003254): 471-476.

[本文引用: 1]

Wu QingbaiZhang TingjunLiu Yongzhi.

Permafrost temperatures and thickness on the Qinghai-Tibet Plateau

[J]. Global and Planetary Change, 2010721/2): 32-38.

[本文引用: 1]

Hinkel K MDoolittle J ABockheim J Get al.

Detection of subsurface permafrost features with ground-penetrating radar, Barrow, Alaska

[J]. Permafrost and Periglacial Processes, 2001122): 179-190.

Zorigt MKwadijk JVan Beek Eet al.

Estimating thawing depths and mean annual ground temperatures in the Khuvsgul region of Mongolia

[J]. Environmental Earth Sciences, 20167510): 897.

Little J DSandall HWalegur M Tet al.

Application of differential global positioning systems to monitor frost heave and thaw settlement in tundra environments

[J]. Permafrost and Periglacial Processes, 2003144): 349-357.

[本文引用: 3]

Zhu Rui.

Study on ground deformation of permafrost in Maxian Mountain

[D]. LanzhouLanzhou Jiaotong University2017.

[本文引用: 1]

朱瑞.

马啣山多年冻土地面变形研究

[D]. 兰州兰州交通大学2017.

[本文引用: 1]

Ma WeiLiu DuanWu Qingbai.

Monitoring and analysis of embankment deformation in permafrost regions of Qinghai-Tibet Railway

[J]. Rock and Soil Mechanics, 2008293): 571-579.

[本文引用: 1]

马巍刘端吴青柏.

青藏铁路冻土路基变形监测与分析

[J]. 岩土力学, 2008293): 571-579.

[本文引用: 1]

Shi GangqiangZhao ShiyunLi Xianminget al.

The frost heaving deformation of high-speed railway subgrades in cold regions: monitoring and analyzing

[J]. Journal of Glaciology and Geocryology, 2014362): 360-368.

[本文引用: 1]

石刚强赵世运李先明.

严寒地区高速铁路路基冻胀变形监测分析

[J]. 冰川冻土, 2014362): 360-368.

[本文引用: 1]

Sofia AHenriette STazio Set al.

Thaw subsidence of a yedoma landscape in northern Siberia, measured in situ and estimated from TerraSAR-X interferometry

[J]. Remote Sensing, 2018104): 494.

[本文引用: 1]

Li ShanshanLi ZhiweiHu Junet al.

Investigation of the seasonal oscillation of the permafrost over Qinghai-Tibet Plateau with SBAS-InSAR algorithm

[J]. Chinese Journal of Geophysics, 2013565): 1476-1486.

[本文引用: 1]

李珊珊李志伟胡俊.

SBAS-InSAR技术监测青藏高原季节性冻土形变

[J]. 地球物理学报, 2013565): 1476-1486.

[本文引用: 1]

Graham L C.

Synthetic interferometer radar for topographic mapping

[J]. Proceedings of the IEEE, 1974626): 763-768.

[本文引用: 1]

Evans D LAlpers WCazenave Aet al.

Seasat: a 25-year legacy of success

[J]. Remote Sensing of Environment, 2005943): 384-404.

Gabriel A KZebker H AGoldstein R M.

Mapping small elevation changes over large areas: differential radar interferometry

[J]. Journal of Geophysical Research, 198994B7): 9183-9191.

Massonnet DRossi MCarmona Cet al.

The displacement field of the Landers earthquake mapped by radar interferometry

[J]. Nature, 19933646433): 138-142.

[本文引用: 1]

Wang ZhijunLi Shusun.

Detection of winter frost heaving of the active layer of Arctic permafrost using SAR differential interferograms

[C]// Proceedings of IEEE International Geoscience and Remote Sensing Symposium 1999: vol 4. New YorkIEEE, 1999: 1946-1948.

[本文引用: 1]

Hoffmann JZebker H AGalloway Det al.

Seasonal subsidence and rebound in Las Vegas Valley, Nevada, observed by synthetic aperture radar interferometry

[J]. Water Resources Research, 2001376): 1551-1566.

[本文引用: 1]

Li ZhenLi XinwuLiu Yongzhiet al.

Detecting the displacement field of thaw settlement by means of SAR interferometry

[J]. Journal of Glaciology and Geocryology, 2004264): 389-396.

[本文引用: 2]

李震李新武刘永智.

差分干涉SAR冻土形变检测方法研究

[J]. 冰川冻土, 2004264): 389-396.

[本文引用: 2]

Zebker H AVillasenor J.

Decorrelation in interferometric radar echoes

[J]. IEEE Transactions on Geoscience and Remote Sensing, 1992305): 950-959.

[本文引用: 1]

Berardino PFornaro GLanari Ret al.

A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms

[J]. IEEE Transactions on Geoscience and Remote Sensing, 20024011): 2375-2383.

[本文引用: 2]

Ferretti APrati CRocca F.

Multibaseline InSAR DEM reconstruction: the wavelet approach

[J]. IEEE Transactions on Geoscience and Remote Sensing, 1999372): 705-715.

[本文引用: 1]

Ferretti APrati CRocca F.

Permanent scatters in SAR interferometry

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001391): 8-19.

[本文引用: 1]

Casu FManzo MLanari R.

A quantitative assessment of the SBAS algorithm performance for surface deformation retrieval from DInSAR data

[J]. Remote Sensing of Environment, 20061023/4): 195-210.

[本文引用: 1]

Zhang YonghongWu Hong’anSun Guangtong.

Deformation model of time series interferometric SAR techniques

[J]. Acta Geodaetica et Cartographica Sinica, 2012416): 864-869.

[本文引用: 1]

张永红吴宏安孙广通.

时间序列InSAR技术中的形变模型研究

[J]. 测绘学报, 2012416): 864-869.

[本文引用: 1]

Hooper ABekaert DSpaans Ket al.

Recent advances in SAR interferometry time series analysis for measuring crustal deformation

[J]. Tectonophysics, 2012514/515/516/5171-13.

[本文引用: 1]

Lanari RMora OManunta Met al.

A small-baseline approach for investigating deformations on full-resolution differential SAR interferograms

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004427): 1377-1386.

[本文引用: 1]

Liu LinZhang TingjunWahr J.

InSAR measurements of surface deformation over permafrost on the North Slope of Alaska

[J/OL]. Journal of Geophysical Research: Earth Surface, 2010115(F3) [2021-08-01]. .

[本文引用: 2]

Xie ChouLi ZhenLi Xinwu.

A study of deformation in permafrost regions of Qinghai-Tibet Plateau based on ALOS/PALSAR D-InSAR interferometry

[J]. Remote Sensing for Land and Resources, 2008203): 15-19.

[本文引用: 1]

谢酬李震李新武.

基于PALSAR数据的青藏高原冻土形变检测方法研究

[J]. 国土资源遥感, 2008203): 15-19.

[本文引用: 1]

Liu LSchaefer KZhang Tet al.

Estimating

1992-2000 average active layer thickness on the Alaskan North Slope from remotely sensed surface subsidence[J/OL]. Journal of Geophysical Research: Earth Surface, 2012117(F1) [2021-08-01]. .

[本文引用: 3]

Liu LSchaefer KGusmeroli Aet al.

Seasonal thaw settlement at drained thermokarst lake basins, Arctic Alaska

[J]. The Cryosphere, 201483): 815-826.

Liu LJafarov E ESchaefer K Met al.

InSAR detects increase in surface subsidence caused by an Arctic tundra fire

[J]. Geophysical Research Letters, 20144111): 3906-3913.

[本文引用: 1]

Cheng Guodong.

The formation process of thick subsurface ice

[J]. Science in China: Series B, 1982123): 281-288.

[本文引用: 1]

程国栋.

厚层地下冰的形成过程

[J]. 中国科学: B辑, 1982123): 281-288.

[本文引用: 1]

Cheng Guodong.

Unidirectional accumulation effect of unfrozen water in seasonal freezing and thawing layers

[J]. Chinese Science Bulletin, 198123): 42-45.

[本文引用: 1]

程国栋.

季节冻结和融化层中未冻水的单向积聚效应

[J]. 科学通报, 198123): 42-45.

[本文引用: 1]

Zhao LinHu GuojieZou Defuet al.

Permafrost changes and its effects on hydrological processes on Qinghai-Tibet Plateau

[J]. Bulletin of Chinese Academy of Sciences, 20193411): 1199-1212.

[本文引用: 1]

赵林胡国杰邹德富.

青藏高原多年冻土变化对水文过程的影响

[J]. 中国科学院院刊, 20193411): 1199-1212.

[本文引用: 1]

Xu XiaozuWang JiachengZhang Lixin. Physics of frozen soil[M]. 2nd ed. BeijingScience Press2010.

[本文引用: 1]

徐敩祖王家澄张立新. 冻土物理学[M]. 2版. 北京科学出版社2010.

[本文引用: 1]

Zhou YouwuGuo DongxinQiu Guoqinget al. Geocryology in China[M]. BeijingScience Press2000.

[本文引用: 1]

周幼吾郭东信邱国庆. 中国冻土[M]. 北京科学出版社2000.

[本文引用: 1]

Short NBrisco BCouture Net al.

A comparison of TerraSAR-X, RADARSAT-2 and ALOS-PALSAR interferometry for monitoring permafrost environments, case study from Herschel Island, Canada

[J]. Remote Sensing of Environment, 201111512): 3491-3506.

[本文引用: 1]

Zhang ZhengjiaWang ChaoTang Yixianet al.

Surface deformation monitoring using time series TerraSAR-X images over permafrost of Qinghai-Tibet Plateau, China

[C]// 2016 IEEE International Geoscience and Remote Sensing Symposium. New YorkIEEE20165963-5966.

[本文引用: 2]

Wang Chunjiao.

Land surface deformation research of permafrost degradation area in Northeast China based on D-InSAR

[D]. HarbinNortheast Forestry University2015.

[本文引用: 1]

王春娇.

基于D-InSAR的东北多年冻土退化区地表形变研究

[D]. 哈尔滨东北林业大学2015.

[本文引用: 1]

Strozzi TAntonova SGünther Fet al.

Sentinel-1 SAR interferometry for surface deformation monitoring in low-land permafrost areas

[J]. Remote Sensing, 2018101360.

[本文引用: 5]

Daout SDoin M PPeltzer Get al.

Large-scale InSAR monitoring of permafrost freeze-thaw cycles on the Tibetan Plateau

[J]. Geophysical Research Letters, 2017442): 901-909.

[本文引用: 4]

Yuan Y.

Measuring surface deformation caused by permafrost thawing using radar interferometry, case study: Zackenberg, NE Greenland

[D]. Delft, the NetherlandsDelft University of Technology2011.

[本文引用: 4]

Singhroy VCouture RAlasset P Jet al.

InSAR monitoring of landslides on permafrost terrain in Canada

[C]// 2007 IEEE International Geoscience and Remote Sensing Symposium. New YorkIEEE20072451-2454.

[本文引用: 1]

Li Shanshan.

The study of using SBAS to monitor the motion of frozen soil along Qinghai-Tibet Railway

[D]. ChangshaCentral South University2012.

[本文引用: 5]

李珊珊.

基于SBAS技术的青藏铁路区冻土形变监测研究

[D]. 长沙中南大学2012.

[本文引用: 5]

Cong Chong.

Study on the surface deformation of permafrost along transportation corridor in plateau based on Sentinel-1A satellite interferometric data

[D]. BeijingBeijing Jiaotong University2018.

[本文引用: 1]

丛充.

基于Sentinel-1A卫星干涉数据的高原冻土交通走廊地表形变研究

[D]. 北京北京交通大学2018.

[本文引用: 1]

Wu Z MZhao LLiu Let al.

Surface-deformation monitoring in the permafrost regions over the Tibetan Plateau, using Sentinel-1 data

[J]. Sciences in Cold and Arid Regions, 2018102): 114-125.

[本文引用: 2]

Chen FulongLin HuiZhou Weiet al.

Surface deformation detected by ALOS PALSAR small baseline SAR interferometry over permafrost environment of Beiluhe section, Tibet Plateau, China

[J]. Remote Sensing of Environment, 201313810-18.

[本文引用: 1]

Li ZhenTang PanpanZhou Jianminet al.

Permafrost environment monitoring on the Qinghai-Tibet Plateau using time series ASAR images

[J]. International Journal of Digital Earth, 2015810): 840-860.

[本文引用: 1]

Zhou HuayunZhao LinTian Liminget al.

Monitoring and analysis of surface deformation in the permafrost area of Wudaoliang on the Tibetan Plateau based on Sentinel-1 data

[J]. Journal of Glaciology and Geocryology, 2019413): 525-536.

[本文引用: 1]

周华云赵林田黎明.

基于Sentinel-1数据对青藏高原五道梁多年冻土区地面形变的监测与分析

[J]. 冰川冻土, 2019413): 525-536.

[本文引用: 1]

Chen JieGünther FGrosse Get al.

Sentinel-1 InSAR measurements of elevation changes over yedoma uplands on Sobo-Sise Island, Lena Delta

[J]. Remote Sensing, 2018107): 1152.

[本文引用: 2]

Jia YuanyuanKim J WShum C Ket al.

Characterization of active layer thickening rate over the northern Qinghai-Tibetan Plateau permafrost region using ALOS interferometric synthetic aperture radar data, 2007-2009

[J]. Remote Sensing, 2017984.

[本文引用: 2]

Zou DefuZhao LinSheng Yuet al.

A new map of permafrost distribution on the Tibetan Plateau

[J]. The Cryosphere, 2017116): 2527-2542.

[本文引用: 1]

Qin Dahe. Glossary of cryosphere science[M]. BeijingChina Meteorological Press2014.

[本文引用: 1]

秦大河. 冰冻圈科学辞典[M]. 北京气象出版社2014.

[本文引用: 1]

Pang QiangqiangLi ShuxunWu Tonghuaet al.

Simulated distribution of active layer depths in the frozen ground regions of Tibetan Plateau

[J]. Journal of Glaciology and Geocryology, 2006283): 390-395.

[本文引用: 1]

庞强强李述训吴通华.

青藏高原冻土区活动层厚度分布模拟

[J]. 冰川冻土, 2006283): 390-395.

[本文引用: 1]

Zhang ZhongqiongWu Qingbai.

Predicting changes of active layer thickness on the Qinghai-Tibet Plateau as climate warming

[J]. Journal of Glaciology and Geocryology, 2012343): 505-511.

[本文引用: 1]

张中琼吴青柏.

气候变化情景下青藏高原多年冻土活动层厚度变化预测

[J]. 冰川冻土, 2012343): 505-511.

[本文引用: 1]

Wang ChenghaiJin ShuanglongWu Zhongyuanet al.

Evaluation and application of the estimation methods of frozen (thawing) depth over China

[J]. Advances in Earth Science, 2009242): 132-140.

[本文引用: 1]

王澄海靳双龙吴忠元.

估算冻结(融化)深度方法的比较及在中国地区的修正和应用

[J]. 地球科学进展, 2009242): 132-140.

[本文引用: 1]

Zhang Zhengjia.

Research on Qinghai-Tibet permafrost environment and engineering using high resolution SAR images

[D]. BeijingChinese Academy of Sciences2017.

[本文引用: 1]

张正加.

高分辨率SAR数据青藏高原冻土环境与工程应用研究

[D]. 北京中国科学院大学2017.

[本文引用: 1]

Zhang YunjunFattahi HAmelung F.

Small baseline InSAR time series analysis: unwrapping error correction and noise reduction

[J]. Computers & Geosciences, 2019133104331.

[本文引用: 1]

Li PingxiangYang Jie. Principle and application of radar measurement[M]. BeijingSurveying and Mapping Press2006.

[本文引用: 1]

李平湘杨杰. 雷达测量原理与应用[M]. 北京测绘出版社2006.

[本文引用: 1]

Chen C WZebker H A.

Phase unwrapping for large SAR interferograms: statistical segmentation and generalized network models

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002408): 1709-1719.

[本文引用: 1]

Cheng Guodong.

The forming process of thick layered ground ice

[J]. Science in China: Series B, 19827): 100-111.

[本文引用: 1]

Mackay J R.

Downward water movement into frozen ground, western arctic coast, Canada

[J]. Canadian Journal of Earth Sciences, 1983201): 120-134.

[本文引用: 1]

Nelson F EAnisimov O AShiklomanov N I.

Subsidence risk from thawing permafrost

[J]. Nature, 20014106831): 889-890.

[本文引用: 1]

Jorgenson M TRacine C HWalters J Cet al.

Permafrost degradation and ecological changes associated with a warming climate in central Alaska

[J]. Climatic Change, 2001484): 551-579.

[本文引用: 1]

Osterkamp T.

The recent warming of permafrost in Alaska

[J]. Global and Planetary Change, 2005493/4): 187-202.

[本文引用: 1]

Tarnocai CCanadell J GSchuur E A Get al.

Soil organic carbon pools in the northern circumpolar permafrost region

[J/OL]. Global Biogeochemical Cycles, 2009232) [2021-08-01]. .

[本文引用: 1]

/