多年冻土区铁路路基热状况对工程扰动及气候变化的响应

基于青藏铁路沿线长期地温监测资料,对天然场地及铁路路基下部的浅层地温、多年冻土上限及下伏冻土地温动态变化过程进行对比分析,研究多年冻土区铁路路基热状况对于工程扰动及气候变化的响应过程.监测结果表明,路基修筑后边坡热效应显著,由此导致路基下部多年冻土热状况的不对称分布,必须引起足够的重视.块石路基修筑后,下部多年冻土上限抬升显著,其中阴坡路肩下抬升幅度普遍较阳坡路肩下显著.普通路基修筑后,在年平均地温低于?0.6～?0.7℃的地区下部多年冻土上限有不同程度的抬升,而在年平均地温高于?0.6℃的地区下部冻土上限则出现了一定程度的下降,其中阳坡路肩下降幅显著.受块石层冷却降温作用,低温冻土区块石路基下部浅层冻土地温有明显降温过程,而在高温冻土区这一降温趋势只存在于阴坡路肩下.对于普通路基,多年冻土上限抬升后,浅层冻土地温存在一定的升温过程.对于气候变暖,低温冻土区多年冻土的响应主要集中体现在冻土升温上,而高温冻土区多年冻土的响应则主要表现为冻土上限下降,冻土厚度减小.基于上述监测结果,可将目前青藏铁路路基热状况分为稳定型(低温冻土区块石路基)、亚稳定型(低温冻土区普通路基及高温冻土区块石路基)和不稳定型(高温冻土区普通路基).     

人类工程活动下冻土环境变化评价模型

提出用冻土热稳定性、热融敏感性及地表景观稳定性来评价人类工程活动下冻土环境变化, 并提出冻土环境的定量评价模型. 利用青藏公路沿线28个地温监测资料, 对冻土环境评价模型进行计算. 分析冻土热融敏感性与冻融过程和季节融化深度之间的关系; 冻土热稳定性与多年冻土顶板温度、年平均地温及季节融化深度之间的关系; 地表景观稳定性与冻融灾害之间的关系及产生的可能性. 结果表明, 冻土热稳定性、热融敏感性和地表景观稳定性可用来评价和预报人类工程活动下冻土环境的变化特征.     

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

冻土的冻结和融化的反复交替会造成地质环境与结构的破坏,从而导致房屋和道路等地面工程建筑物的地基破裂或者塌陷,还会引起山体滑坡、洪水暴发以及冰川移动等.因此,监测冻土形变对确保冻土区工程建筑的稳定性和安全性,同时保证冻土区社会经济可持续发展具有重要的意义.目前,在冻土监测方面并没有能大面积监测冻土形变时间演化情况的有效方法,本文提出将InSAR技术中的小基线集方法(SBAS-InSAR)应用于监测冻土来获取其形变时间序列中.考虑到冻土形变呈现明显的季节性特征,本文提出利用周期形变模型来代替传统SBAS方法中的线性形变模型,从而更好地分离出高程残差和大气误差.利用ENVISAT卫星获取的21景ASAR影像图作为实验数据,采用改进的SBAS技术成功获取了青藏高原从羊八井站至当雄站铁路段冻土区的地表形变时间序列图,揭示了该冻土区从2007年到2010年的季节性形变演化情况.通过与研究地区温度变化的联合分析,发现所得到的地表形变结果与冻土的物理变化规律非常吻合,证明了SBAS-InSAR技术在冻土形变监测中具有良好的发展应用前景.     

青藏高原多年冻土区天然气水合物形成条件模拟研究

基于野外气体地球化学调查研究,以及前人有关冻土表层温度、冻土层内地温梯度、冻土层下地温梯度等的资料,对青藏高原多年冻土区天然气水合物的形成条件开展了模拟研究. 结果显示：研究区冻土条件能够满足天然气水合物形成的基本要求;气体组成、冻土特征（如冻土厚度或冻土表层温度、冻土层内地温梯度、冻土层下地温梯度等）是影响研究区天然气水合物稳定带厚度的最重要因素,其在不同点位上的差异性可能导致天然气水合物分布的不均匀性的主要原因;研究区最可能的天然气水合物为甲烷与重烃（乙烷和丙烷）的混合气体型天然气水合物;在天然气水合物分布的区域,其产出的上临界点深度在几十至一百多米间,下临界点深度在几百至近一千米间,厚度可达到几百米. 与Canadian Mallik三角洲多年冻土区相比,青藏高原多年冻土区除了冻土厚度小些外,其他条件,如冻土层内地温梯度、冻土层下地温梯度、气体组成等条件较为相近,具有一定的可比性,预示着良好的天然气水合物潜力.     

气候变化条件下东北地区多年冻土变化预测

东北多年冻土（除非指明是季节冻土,以下将多年冻土简称冻土）是中国第二大冻土分布区,主要发育＂兴安-贝加尔型＂冻土.由于处在欧亚大陆冻土区南缘,冻土的热稳定性差,寒区生态的敏感性强.在气候变暖条件下,冻土已经和正在发生着＂三向＂退化.为预测冻土南界和地温变化,根据47个气象站资料并在SHAW模型对植被影响地表温度修正的基础上,建立了冻土地表温度分布的等效纬度模型.结合非稳态热传导模型的有限元数值计算,以多模型结合的方法,进一步计算和分析了目前、50年和100年后冻土地温分区变化.结果表明,在目前地表温度为1.5℃范围,仍可残留冻土.以0.048℃a-1气温递增速率,在目前地表温度为0.5℃和-0.5℃的区域,50年和100年后各自仍有可能存在冻土;冻土面积将由现在的2.57×105 km2各自减至1.84×105和1.29×105 km2,分别减少28.4%和49.8%,且东部退化幅度大于西部.同时,区域地温升高,冻土厚度减薄;稳定型（年平均地温Tcp≤-1.0℃）冻土面积逐渐减小,将由现在的1.07×105 km2分别减少至8.8×104 km2（50年后）和5.6×104 km2（100年后）.相应地,不稳定型（Tcp〉-1.0℃）多年冻土和季节冻土的面积增加,冻土南界将显著北移.冻土环境的变化,将给东北寒区工程设施和生态环境带来重要影响.减少或避免人为地改变冻土赋存条件,是保护冻土环境较可行的途径.     

不同升温模式下冻土地区装配式基础热稳定性研究

以IPCC SRES A2、A1B、B1三种气候变化模式为基础,利用数值方法研究了青藏直流联网工程冻土区装配式基础的冻融过程以及活动层、融化深度、地温的变化规律.结果表明：工程扰动和气候变暖改变了冻土的热状态,促进了冻土退化,均为影响基础长期稳定性的重要因素,其中混凝土桩基的强化导热作用加剧了冻融过程,气候升温导致活动层厚度增加,土层温度升高;随着深度的增加,冻土响应减弱,冻土温度变化幅度越小;在三种升温模式下50年后融化深度分别达到3.12m、5.07m和6.02m,而同期天然场地活动层厚度为2.07m、4.37 m、5.62 m,说明冻土对不同升温模式的响应程度不同,且中心冻土在气候变暖和工程扰动双重影响下退化更快;从第10年到第50年,这三种模式下桩基中心融化速率分别为1.5 cm·α-1,6.2 cm·α-1,8.6 cm·α-1,即随着升温速率的增加,土层融化深度增加,冻土退化速度加快;低升温率时冻土变化主要受工程作用,而在较高升温模式下冻土退化则主要受气候变暖的影响.     

利用InSAR技术研究青藏高原冻土形变特征

由于自然环境的限制, 青藏高原的大地测量网络十分稀疏, 不能满足区域地壳运动监测的需求。干涉合成孔径雷达(InSAR)是非接触监测地壳运动的一种重要方式, 但在高原上受到冻土的影响。本文基于2014~2018年的Sentinel-1卫星C波段雷达数据, 采用InSAR时序技术分析了冻土形变的时空特征。针对InSAR位移时间序列, 采用空间滤波去除了大气延迟、 地形效应等局部公共误差, 提高了时间序列的信噪比。结果显示, 青藏高原的冻土运动可分为差异较大的两类: 在大部分冻土区域, 与周边高山(基岩)区域相比, 冻土地区显示类似的季节波动或一定的长期沉降; 在部分冻土地区存在异常快速下沉区域, 例如在西藏中部布若错湖西南侧的沉积盆地内, 存在一个直径约2 km的漏斗型沉降区, LOS向沉降速率可达约10±2.1 mm/a。构造运动造就了高原上大量沿断裂线分布的河流、 湖泊, 河床和沉积盆地等广泛分布着冻土, 给准确分析构造形变带来很大挑战, 本文所得的结果可作为区分冻土运动与构造变形的一种有效判据, 也有益于研究高原冻土的物理特性及变形机理。     

青藏公路沿线冻土的地温特征及退化方式

青藏高原多年冻土(以下简称冻土)具有地域分布广、厚度薄及稳定性差等特征. 过去几十年的气候变暖背景下, 冻土广泛退化, 地温升高, 夏季最大融化深度加深, 冬季冻结深度减小. 冻土已经产生下引式、上引式和侧引式退化. 冻土层厚度减薄, 或者在某些地区彻底消失. 冻土退化模式研究在冻土学、寒区工程和寒区环境管理方面具有重要意义. 由南至北穿越560 km冻土区的青藏公路沿线(简称青藏线)冻土在青藏高原腹地具有很好的代表性. 在水平方向上, 冻土退化在多年冻土下界附近的零星冻土分布区、融区边缘和岛状冻土区表现得更为明显. 当最大季节融化深度超过最大季节冻结深度时, 冻土开始下引式退化; 通常形成融化夹层, 造成多年冻土和季节冻结层不衔接. 当多年冻土层中地温梯度减小到小于下伏或周边融土层时, 则产生上引式或侧引式退化. 下引式退化进程可分为4个阶段: (1) 初始退化阶段, (2) 加速退化阶段, (3) 融化夹层阶段, (4) 最终多年冻土彻底融化为季节冻土阶段. 当多年冻土中地温梯度降至下伏融土层地温梯度以下时, 则产生上引式退化. 3种类型冻土温度曲线(稳定型、退化型和相变过渡型)展现了这些退化模式. 虽然存在不同地段和类型的地温特征, 三种退化模式的各种组合最终将使多年冻土消融, 转变成季节冻土. 过去25年来, 青藏线冻土年平均下引式退化速率变化在6~25 cm, 年平均上引式退化速率在12~30 cm, 零星多年冻土区年平均侧引式退化速率为62~94 cm. 这些观测结果超过所报道的过去20年来阿拉斯加亚北极不连续冻土区4 cm的年平均退化速率, 蒙古国不连续冻土区的4~7 cm的年平均退化速率, 以及雅库悌共和国亚北极和阿拉斯加北极稳定性冻土区退化速率.     

基于方差分量估计的多源InSAR数据自适应融合形变测量

近年来,合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar,InSAR)技术在地面沉降监测方面展现了巨大的应用潜力,但受其重访周期和一维形变测量能力的限制,仅利用单一轨道卫星观测数据很难揭示真实的地表形变特征及其演化规律.随着在轨运行的SAR卫星系统不断增加,使得融合相同时间段内覆盖同一区域的多源多轨道InSAR数据成为可能.然而目前普遍采用的多源InSAR数据融合方法均为针对大尺度形变监测设计,或者忽略南北向形变甚至水平形变,容易造成误判.为此,本文对经典小基线集(Small Baseline Subset, SBAS)时序InSAR分析方法进行改进,在其形变反演模型中加入东西向和南北向形变参数,采用方差分量估计方法解算多源观测数据验后方差,通过迭代精化确定权重矩阵,从而获得形变参数的最优估值.使用美国南加州地区的ALOS PALSAR和ENVISAT ASAR数据开展实验,利用南加州综合GPS网(SCIGN)位于研究区域内的9个站点观测数据进行验证,结果表明本文方法得到的融合形变测量结果在垂直向上能够准确反映地表形变波动,周期性与GPS观测比较一致;同时,融合得到的三维形变场显示南加州洛杉矶地区存在不可忽略的水平形变,东西向形变测量精度略高于南北向.因此,基于方差分量估计的多源InSAR融合方法在提高形变测量时间序列连续性的同时,能够更准确地反演研究区域三维形变特征.     

时序InSAR同质样本选取算法研究

分布式雷达目标时序InSAR技术是目前InSAR形变监测领域的主流方向，其中同质样本选取是该技术的基础，其估计精度直接影响SAR影像分辨率与后续参数解算精度.本文在追踪最新研究进展之上，系统回顾了当今统计同质选点算法的优缺点.在参数与非参数两类统计方法的应用中，采用蒙特卡罗方法和真实数据验证定量比较算法差异以及适用场景.根据之前的研究结论，提出一种改进的最优参数统计同质样本选择方法.最后，论文介绍了团队研发的MATLAB开源工具包，涵盖了同质样本提取和时序InSAR协方差矩阵估计两部分内容，为InSAR科研人员和后续数据处理提供高质量、全分辨率的观测源.     

高亚洲冰川质量变化趋势的卫星重力探测

利用高亚洲地区32个Mascon,基于GRACE RL05时变重力场模型频域和空域上的两种计算方法有效分离并提取出高亚洲冰川及其毗邻地区的等效水质量变化,得到2002—2013年期间高亚洲地区更为可靠的Mascon质量变化.高亚洲冰川质量变化的空间特征是:青藏高原内陆地区以正增长为主,边缘地区以负增长为主,在藏东南的最边缘地区冰川质量损失最为严重.天山地区、帕米尔和昆仑山地区、喜马拉雅山和喀喇昆仑山地区、青藏高原内陆地区冰川质量的平均变化趋势分别为-2.8±0.9Gt/a、-3.3±1.5Gt/a、-9.9±2.1Gt/a和5.0±0.8Gt/a,高亚洲冰川质量整体的平均变化趋势为-11.0±2.9Gt/a.印度等北部平原地区地下水平均变化趋势为-35.0±4.2Gt/a,该地区地下水信号泄漏是影响GRACE研究高亚洲冰川质量变化的关键因素,频域法和空域法能有效改正该地区地下水信号泄漏的影响.     

GNSS影像及其时空特征初探

地壳垂向形变在连续空间和时间域内呈现显著特征,探求其时空变化特征有助于理解地球物理过程,为研究地球内部相互作用机制提供支持.本文使用美国西部地区PBO与中国大陆CMONOC两个GNSS网测站的坐标时间序列,通过基于中位数并顾及年际差异的非参数方法(MIDAS方法)估计测站的速度与不确定性;建立空间结构函数(SSF)并确定区域内各测站间的相对关系;以此为基础,构建顾及空间结构的滤波器(MSF)以剔除粗差,增强区域共性;最后,基于MSF与图像处理技术对速度场进行空间加密,生成了研究区域空间内连续的地壳垂向形变图,即区域GNSS影像.随后,针对两个研究区域,分别利用MSF验证实验与棋盘格检测验证了GNSS成像方法的合理性及生成GNSS影像的可靠性;并通过对比使用顾及空间结构滤波前后的各测站速度与不确定性生成的GNSS影像,验证了顾及空间结构的滤波方法在GNSS影像生成中的必要性,并分析了其中存在过度平滑与规则圆弧状突变边缘的问题,讨论了可能的解决方案.最终,将两区域GNSS影像结果与已有的大地测量学及地球动力学结果进行了对比,发现美国西部地区的GNSS影像正确反映出了海岸山脉以峰值速度为2mm·a~(-1),内华达山脉以峰值速度为3mm·a~(-1),以及赫布根湖地区以峰值速度为1.5mm·a~(-1)隆升;洛杉矶地区(峰值速度为-2.5mm·a~(-1)),维多利亚河及其河谷地区(速度为-1.5mm·a~(-1)),以及蛇河平原东部、蒙大拿州西南部(速度为-1mm·a~(-1)左右)的沉降运动;中国大陆的GNSS影像同样反映出喜马拉雅山脉与青藏高原南部(速度呈现为1.0mm·a~(-1)),华北地块与天山地块(速度为1.5mm·a~(-1)与0.3～0.6mm·a~(-1))等典型区域的隆升;长江下游地区以苏锡常地区(速度为-2.1mm·a~(-1))为中心,向外速度逐渐减小的沉降运动,以及华南地区(速度呈现为-0.6～-1.5mm·a~(-1))、东北地区(速度呈现为-0.6～-1.5mm·a~(-1))、塔里木盆地(速度呈现为-1.2mm·a~(-1))等区域的沉降运动.因此,本文认为GNSS影像具有合理性与正确性,有助于地壳垂向形变的整体时空分布特征研究.     

利用GRACE卫星重力数据监测关中地区地下水储量变化

关中地区作为一带一路重要的工农业发达地区之一，开展针对该地区地下水储量变化的监测和分析工作对揭示地下水储量变化特征与经济社会发展具有重要现实意义.本文基于2003—2014年GRACE卫星重力场模型数据，采用组合滤波及单一尺度因子方法反演了关中地区陆地水储量变化，扣除GLDAS地表水平均结果，对关中地区地下水储量变化进行了监测分析.将陆地水储量变化与GLDAS进行相关性分析，将地下水储量变化与WGHM地下水模型及实测地下水位结果进行对比分析.研究结果表明：①关中地区陆地水变化与GLDAS模型结果具有较强的相关性，相关系数多数大于0.7，其中与模型平均结果的相关系数可达0.8.② 2003—2008年关中地区地下水呈正增长趋势，增加速率为0.25 cm·a^-1，与同期实测数据变化趋势一致；但2003—2013年地下水存在长期亏损，亏损速率为-0.37 cm·a^-1等效水高，这与同时期WGHM估算结果-0.35 cm·a^-1十分吻合.③关中地区地下水存在明显的年变化特征，在2003—2014年期间地下水减少速率为-0.44 cm·a^-1，与该地区降雨量有较好的一致性，在降雨偏少的2008、2012和2013年，地下水也显著减少.     

青藏高原色林错流域区冰川消融对湖泊水量变化的影响

青藏高原大部分湖泊近年来持续扩张,湖泊水位和水量明显增加.冰川消融是流域水量平衡和水循环的重要影响因素,直接导致湖泊水量变化.由于缺乏大范围的冰川质量平衡观测结果,青藏高原冰川消融对湖泊水量变化的影响仍存在较大争议.本文选择青藏高原内流区的色林错流域区(水系编号5Z2)作为研究对象,利用SRTM DEM和TanDEM-X双站InSAR数据,精确估算该流域三个主要冰川区(普若岗日、格拉丹东和西念青唐古拉)2000—2012年的冰川质量平衡,依次为:-0.020±0.030、-0.128±0.049、-0.143±0.032m·w.e.·a~(-1).并据此采用面积加权法准确推估出5Z2流域的冰川质量变化为:-0.166±0.021Gt·a~(-1).综合ICESat和Cryosat-2卫星测高数据,计算该流域2003—2012年湖泊水量变化速率(3.006±0.202Gt·a~(-1)),并定量评估冰川质量变化对5Z2流域湖泊水量增加的贡献为:5.52%±1.07%,因此在青藏高原色林错流域区,冰川消融不是导致21世纪初期湖泊水位上升的主要因素.     

Distribution and changes of active layer thickness (ALT) and soil temperature (TTOP) in the source area of the Yellow River using the GIPL model

Active layer thickness (ALT) is critical to the understanding of the surface energy balance, hydrological cycles, plant growth, and cold region engineering projects in permafrost regions. The temperature at the bottom of the active layer, a boundary layer between the equilibrium thermal state (in permafrost below) and transient thermal state (in the atmosphere and surface canopies above), is an important parameter to reflect the existence and thermal stability of permafrost. In this study, the Geophysical Institute Permafrost Model (GIPL) was used to model the spatial distribution of and changes in ALT and soil temperature in the Source Area of the Yellow River (SAYR), where continuous, discontinuous, and sporadic permafrost coexists with seasonally frozen ground. Monthly air temperatures downscaled from the CRU TS3.0 datasets, monthly snow depth derived from the passive microwave remote-sensing data SMMR and SSM/I, and vegetation patterns and soil properties at scale of 1:1000000 were used as input data after modified with GIS techniques. The model validation was carried out carefully with ALT in the SAYR has significantly increased from 1.8 m in 1980 to 2.4 m in 2006 at an average rate of 2.2 cm yr^?1. The mean annual temperature at the bottom of the active layer, or temperature at the top of permafrost (TTOP) rose substantially from ?1.1°C in 1980 to ?0.6°C in 2006 at an average rate of 0.018°C yr^?1. The increasing rate of the ALT and TTOP has accelerated since 2000. Regional warming and degradation of permafrost has also occurred, and the changes in the areal extent of regions with a sub-zero TTOP shrank from 2.4×10⁴ to 2.2×10⁴ km² at an average rate of 74 km² yr^?1. Changes of ALT and temperature have adversely affected the environmental stability in the SAYR.     

鄂尔多斯块体周缘地区现今地壳水平运动与应变

位于青藏块体和华北块体之间的鄂尔多斯块体及其周缘地区是中国大陆构造活动最活跃的地区之一,从1300年至今,在块体周边断陷盆地和西南缘断裂带上发生了五次8级以上的地震.为了了解该地区现今地壳运动、应变状态以及断裂滑动分布,我们收集了中国大陆构造环境监测网络2009—2013年、国家GPS控制网、跨断陷盆地的8个GPS剖面等共527个流动站和32个连续站GPS观测数据,获得了高空间分辨率的地壳水平运动速度场,进一步用均匀弹性模型计算了应变率分布.结果表明,块体内部GPS站点向NEE方向运动,速度变化较小,应变率大多在(-1.0~1.0)×10~(-8)/a之间;山西断陷带构造运动与变形最为强烈,盆地相对于鄂尔多斯块体为拉张变形,应变率为(1.0~3.0)×10~(-8)/a,相对于东部山地则为挤压变形,应变率为(-2.0~-3.0)×10~(-8)/a,盆地西侧断裂(如罗云山断裂、交城断裂)以拉张运动为主,拉张速率为2~3mm·a-1,盆地东侧断裂主要以右旋缩短运动为主,速率为1~3mm·a-1;河套断陷带西部的临河凹陷处于较强的张性应变状态,应变率为(2.0~3.0)×10~(-8)/a;块体西南边缘处于压缩应变状态,应变率为(-1.0~-2.0)×10~(-8)/a,六盘山断裂存在明显的地壳缩短运动,速率约为2.1mm·a-1,速率在断裂附近逐渐减小,反映了断裂处于闭锁状态;相对于鄂尔多斯块体内部渭河断裂带为左旋运动,速率为1.0mm·a-1,盆地处在弱拉张变形状态.     

基于迭代PCA的GPS时间序列震后形变估计方法和应用

GPS时间序列的震后形变分析对于研究区域震后形变机制和岩石圈流变学特性以及维持国际动态地球参考框架具有重要意义.本文在现有参数估计方法的基础上,提出兼顾震后形变衰减特征空间相关性和整体建模的"迭代PCA参数估计方法",并利用模拟数据证实了新方法可以获取更稳健可靠的震后形变、同震形变和震间速度参数.最后,以37个新西兰GPS连续站坐标时间序列为例,利用迭代PCA方法提取了2016年11月13日Kaikoura地震共性的震后形变时间演化过程和各站点的震后形变,并定量分析了震后形变对地表速度的影响.结果表明各站的震后形变在时间域上以衰减常数τ为4天的对数模型持续松弛;空间域上南岛北东部和北岛最南部震后形变较大,其中,最大震后形变点为cmbl站,截至2017年6月10日NEU方向累计的震后形变分别达到107mm,135mm和187mm,地表速度分别达到133.58mm·a~(-1),112.05mm·a~(-1)和175.58mm·a~(-1),仍高于稳定的震间速度.     

利用GRACE时变重力场反演青藏高原的隆升速率

青藏高原隆升对中国、亚洲乃至世界的气候都有着重要影响,研究青藏高原地壳隆升速率具有重大意义.本文利用2004—2015年期间高覆盖度的卫星重力数据,通过去除陆地储水的重力效应获得地壳隆升引起的重力变化速率,基于直立长方体垂直运动与重力变化的关系模型反演了该区域的地壳隆升速率分布.研究结果表明在300 km的空间尺度下青藏高原隆升速率分布具有不均匀的特点,表现为以冈底斯山—唐古拉山—鲜水河断裂带为界线,其两侧的速率差异较大.位于界线以南,沿喜马拉雅推覆构造带的区域平均隆升速率为2.01±0.87 mm·a^-1,其中西侧的印度板块与东侧的缅甸板块隆升速率分别为~2.43 mm·a^-1、~2.89 mm·a^-1,位于两板块之间的区域隆升速率为~0.69 mm·a^-1;位于界线以北,除了天山区域和华北板块的隆升速率为~1 mm·a^-1,其他区域隆升现象不明显,其速率为~0 mm·a^-1.我们发现存在两条均穿过正断裂带区域的隆升速率梯度带,其中一条为从加德满都到塔里木盆地,其恰好穿过青藏高原内部的正断裂带,另一条为从那加山到四川盆地,其恰好穿过大理正断裂带.本文反演的青藏高原隆升速率与以往观测到的GPS结果有很好的一致性,为青藏高原隆升、地壳增厚等科学问题提供理论支持.