自50-55Ma以来,印度次大陆向北与欧亚大陆碰撞后形成“喜马拉雅-青藏高原”造山带,碰撞导致地壳增厚致使高原大幅隆升,改变了亚洲大陆岩石圈的构造格局,也对东亚地区的气候和环境产生了巨大的影响。阿尔金断裂作为青藏高原北缘的主控边界断裂,其运动学性质在上世纪七十年代受到关注,不同量级的滑动速率引出了“块体运动与东向逃逸”和“连续变形与地壳增厚”两种端元模型。约10^-15 Ma以来,在青藏高原南部与北部出现地堑与裂谷,为高原东西向拉张运动提供了证据,表明青藏高原开始经历地壳减薄过程。青藏高原形成以来形变场经历怎样变化,长时间尺度的地质学构造过程与现今GPS观测是否能够统一?10^-15Ma以来青藏高原地壳减薄过程造成高原高程怎样的变化?青藏高原北缘尤其是跨阿尔金断裂具有怎样的现今三维地壳变形场,地壳应变是如何在北阿尔金断裂、祁漫塔格断裂和阿尔金断裂之间分配的?青藏高原北缘与塔里木盆地具有怎样的力学性质,对跨阿尔金断裂构造形变场造成怎样的影响?最后,GPS观测得到的现今地表形变场能够对青藏高原形变模式的争论作出何种回答?上述科学问题的回答,对于研究青藏高原隆升与变形过程具有十分重要的意义。本文分为两部分研究内容,第一部分是青藏高原北缘三维震间运动场的观测与研究。在青藏高原北缘跨阿尔金断裂中段自建9个GPS连续台站并开展观测,根据区域研究特点设计无人值守的观测台站,具有低成本投入、高质量观测的特点。上述连续GPS台站的建立填补了青藏高原北缘尤其是在阿尔金无人区地壳形变观测研究的空白,积累了宝贵的连续GPS数据;截止2015年7月,共有4年的连续GPS观测。数据分析结果证明设计建站方法行之有效,GPS台站稳定、观测数据质量稳定、数据连续性稳定。结合使用“中国大陆构造环境监测网络”在研究区及邻域GPS连续台站数据作位置时间序列与速度场解算,获得青藏高原北缘地区跨阿尔金断裂中段现今三维形变场。使用三维线弹性后向滑移（backslip）块体运动模型,反演“塔里木块体”、“北阿尔金块体”、“柴达木块体”和“祁漫塔格块体”的三维块体运动。结果表明,北阿尔金山相对于塔里木盆地有1.32±0.2 mm/yr的抬升速率,相对于柴达木盆地具有0.73±0.3 mm/yr的抬升速率,可解释为北阿尔金块体存在显著的造山过程;阿尔金断裂有8.21±0.60 mm/yr的左旋走滑速率、0.66±0.60 mm/yr的缩短速率;祁漫塔格断裂有0.53±0.60 mm/yr的左旋走滑速率、1.53±0.60 mm/yr的缩短速率;北阿尔金断裂有0.87±0.60 mm/yr的左旋速率、0.69±0.60 mm/yr的缩短速率。同时,阿尔金断裂中、西两段滑动速率基本一致,约在8.0¹0.0 mm/yr范围。定量研究结果支持“连续形变与地壳增厚”模型,表明相对塔里木块体,青藏高原北缘地区正在抬升、增厚,以北阿尔金山地区最为明显,抬升速率达¹.3mm/yr。跨青藏高原北缘的阿尔金断裂、北阿尔金断裂和祁漫塔格断裂近200 km的宽泛变形带内南北向地壳缩短并不明显,缩短量仅为².9 mm,且近一半缩短量发生在祁漫塔格山南侧。GPS观测阿尔金断裂车尔臣河段（⁸6oE）剖面（He et al.,2013）表明断裂两侧存在非对称变形特征。本文采用非对称变形模型反演GPS速度剖面数据,获得断裂两侧塔里木盆地和青藏高原北部的地壳介质剪切模量差异。结果显示塔里木盆地地壳介质剪切模量约为青藏高原北部剪切模量1.53倍,相应S波波速比值为1.24,与前人（Yang et al.,2012）得到地壳和上地幔三维VSV模型结果一致。地震学研究结果（Xie et al.,2013）认为青藏高原北部与东部地区在中地壳存在低速层,局部区域可能发生部分熔融;Hacker et al.（2014）进一步确认羌塘地块中地壳到深部地壳存在熔融现象。本文研究运用与地震学完全不同的资料,通过大地测量方法推导青藏高原北部与塔里木盆地的地壳介质力学性质差异,得到与地震学研究得到的S波波速比及其构造物理学解释相当一致的结果。成果为青藏高原力学演化模型提供新的约束。本论文第二部分内容是使用覆盖青藏高原及周边的GPS速度场,计算青藏高原内部应变率场。GPS观测速度场不仅显示了南东东-北西西向的地壳拉张过程,也揭示了青藏高原内部更加重要的地壳减薄过程。结果显示,青藏高原北部和南部的垂向应变率（减薄应变率）分别为8.9±0.8 nanostrain a^-1和7.4±1.2 nanostrain a^-1,青藏高原西南部的垂向应变率为12.0±3.2 nanostrain a^-1,表明青藏高原内部大尺度范围应变率测量结果的一致性。并且青藏高原内部的拉张应变率观测也相当一致,青藏高原北部,沿着N114o±1°E主应变方向的拉张应变率为21.9±0.4 nanostrain a^-1;高原南部沿着N93±1°E主应变方向的拉张应变率为16.9±0.2 nanostrain a^-1;高原西南部沿着N74±3°E主应变方向的拉张应变率为22.2±1.8 nanostrain a^-1。如果地壳减薄开始于10^-15 Ma,并且现今观测得到应变率适用于整个时间跨度,那么地壳累积减薄5.5-8.5 km。应用Airy地壳均衡理论,青藏高原的平均高程将下降近¹ km。青藏高原北部、南部和西南部相似的垂向应变速率也表明在三个区域的地壳拉张、正断裂运动和地壳减薄过程由相同的物理机制所支配。综合上述两部分研究成果,发现青藏高原现今垂向运动在高原内部和边缘地区存在很大差别。高原内部地区正在经历地壳减薄,而高原边缘地区正在经历不同程度的增厚与隆升。青藏高原北缘地区的垂向应变率⁵-20 nanostrain a^-1,如果考虑重力均衡作用,对应的垂向隆升速率在0.04⁰.14 mm/yr左右。但是对于局部地区如北阿尔金块体,其底部受到塔里木盆地南缘下插挠曲板块的支持,在没有重力均衡情况下,垂向隆升速率可能达到¹ mm/yr。喜马拉雅地区呈现不同水平的垂向形变,垂向应变强烈(¹0-80nanostrain a^-1),山脉底部受到印度下插板片的支持,无法通过重力均衡假定由垂向应变率估计隆升速率。但由GPS与水准数据约束的俯冲板片模型推测山脉隆升速率达到⁷mm/yr。而对于祁连山地区,GPS应变率推测得到垂向应变率²0-40 nanostrain a^-1,应用地壳均衡理论,平均隆升速率为0.15⁰.3 mm/yr;而由于逆冲推覆构造与褶皱变形带的存在,中下地壳有可能仍存在弹性变形,不能实现完全重力均衡,实际隆升速率有可能高于这一估计。本文研究给出青藏高原不同地区三维形变场与形变速率的定量估计,是对“连续形变与地壳增厚”形变模型的重要修正。结果并不支持“块体运动与东向逃逸”模型,并认为高原南北双向俯冲模型中的塔里木块体南向俯冲几乎不存在。