帕米尔高原现代地壳运动首期GPS观测及处理

帕米尔高原现代地壳运动首期GPS观测及处理摘要：介绍了在帕米尔高原进行全球定位系统（GPS）观测的剖面布设和首期观测情况，详细描述了GAMIT/GLOBK软件的处理流程，对研究区域首期GPS观测数据进行处理，并将结果归算到ITRF08参考框架下，GAMIT计算得到的准观测值均满足结果需要的精度，其中单天解标准化均方差都小于0.3mm，最好站点的单程均方根残差处于35mm，最差站点的单程均方根残差处于79mm，满足GAMIT结算的预期结果。GLOBK在ITRF08参考框架下对GAMIT的单天解联合全球解一起结算，并且对处理结果及其误差进行分析与研究，得到GPS站点的精确坐标。结果表明：所有站点的点位精度水平分量为12mm，高程分量为37mm，达到帕米尔地区GPS观测预期的设计要求。这对于研究帕米尔高原运动学特征及探讨青藏高原的隆升机制具有重要意义。关键词：地壳运动；全球定位系统；GAMIT/GLOBK软件；ITRF08参考框架；单天解；全球解；帕米尔0引言全球定位系统（GPS）利用24颗绕地卫星进行实时精确的测定地面点位，以地面GPS固定站点之间的位置及其基线长度随时间变化来测定局部地区以及全球地壳的相对运动。GPS观测信息丰富，全球覆盖面广，具有高精度、全天候、全球性、实时性的特点，为得到高时空分辨率地壳形变的可靠结果提供了保证12。高精度GPS的相对定位已经达到上千千米长基线的观测中误差可达毫米级精度的水平，这样的精度足以检测出地壳运动与构造变形的细微变化37。近年来，国内外学者先后在帕米尔地区展开了一系列的GPS观测及其研究。王琪等针对天山及塔里木地区，在19921999年间多次利用GPS测量直接观测到跨天山西段（76E）每年约20mm的地壳快速缩短率及天山东段（87E）每年约4mm的汇聚变形率89。杨少敏等利用在19922006年间天山境内外地区近400个测站的GPS观测，获取了天山现今地壳运动速度场1011。牛之俊等利用天山境内外地区在19922005年382个GPS测站的原始观测资料计算了天山地区现今地壳运动速率，并得出天山的汇聚速率由西向东逐步减小的结论12。Reigber等于19941998年间沿着塔吉克盆地北边和东边、北帕米尔高原、塔里木盆地到费尔干纳盆地及其南天山到哈萨克地台，布置了92个GPS观测台站（CATS），得到北天山汇聚速率为每年（123）mm，北北西向到东东南向缩短速率为每年（233）mm13。Banerjee等在喜马拉雅西北地区主冲断带MFT布设了26个GPS台站，测得喀喇昆仑西北部的滑动速率为每年（114）mm，并且得出青藏高原南部区域呈东西向伸展，西北部喜马拉雅区域相对于南加帕尔巴特向西运动14。Mohadjer等在30N44N和60E76E区域布置了一系列的GPS站点，得到阿莱南天山的汇聚速率为每年（122）mm，恰曼伽德兹库纳断裂带滑动速率为每年（181）mm，达瓦孜卡拉库尔断裂带为每年（112）mm，塔拉斯费尔干纳断裂带汇聚速率小于每年2mm15。Zubovich等于19942003年间在中亚地区布置了144个连续GPS观测点和400个左右流动GPS观测点，得到塔里木盆地相对于欧亚板块的汇聚速率为每年（202）mm；GPS观测数据结合地质数据显示，帕米尔北缘俯冲到阿莱谷地至少需要速率为每年10mm，最大达每年15mm，塔拉斯费尔干纳断裂带现今滑动的最大速率为每年2mm16。这些GPS观测点主要分布在帕米尔高原西南侧及其东北侧天山附近，在喀喇昆仑断裂地区及其阿里盆地地区的GPS观测点非常稀少。笔者通过对帕米尔地区GPS测点进行合理设置及补充，分别在南北向及其东西向布置了2条剖面进行观测，认识青藏高原的形成过程和大陆动力学理论，探讨岩石圈变形的行为方式。1研究区概况新生代以来，印度板块与欧亚大陆碰撞和持续汇聚作用在中亚地区造成了强烈的陆内形变，引起帕米尔高原大幅度向北推移旋转以及西昆仑山、天山等古老造山带的复活，使帕米尔及其邻区成为大陆内部构造运动十分强烈的地区1718。在整个青藏高原及其附近地区的GPS速率图上，可以明显看出帕米尔地区的GPS观测相对不足（图1）。现有这些GPS观测资料难以约束帕米尔高原一些大型断裂的运动学特征，为了更加紧密具体地得到帕米尔地区一些主要断裂的滑动速率，更好地约束整个青藏高原及其帕米尔活动区域的应变场11，19，笔者在前人工作基础上，于2011年7月在该地区进行了首次观测，并取得了一些结果（图2）。箭头方向表示运动方向，箭头长短表示运动相对速率2GPS观测剖面布设与首期观测本期GPS测点位于帕米尔高原区域（36N42N、68E78E），南边剖面从塔吉克盆地东北部到塔里木盆地西侧，西边剖面从费尔干纳盆地到帕米尔构造结内弧。野外选取无遮挡的出露基岩布设GPS测点，采用TrimbleNetRs接收机和Zephry天线的GPS观测设备，设置采样率为30s，每个点位进行了不低于36h的连续观测。3GPS数据处理本研究采用GAMIT/GLOBK软件（10.40版）对首期GPS观测数据进行处理。GAMIT/GLOBK软件是一种综合性GPS精密分析软件，由美国麻省理工学院（MIT）、哈佛大学天体物理研究中心（CFA）和斯克里普斯海洋研究所（SIO）共同开发和维护，可用来解算地面观测站的相对位置和卫星轨道及其他参数的数据处理。首期GPS观测数据处理主要分为2步进行：利用GAMIT软件进行单天解解算，估计站位置、卫星轨道参数等；利用GLOBK软件获得时间序列检验坐标的一致性。基线处理时，选择周围的IGS（InternationalGNSSService）站与自己剖面测量的数据一起解算。解算时，主要考虑如下因素及其相关设置2024：电离层折射影响用LC观测值消除；使用8个IGS站（BJFS，BRAO，HYDE，IRKT，POL2，KUNM，SELE，TASH）数据作为参考站约束；使用来自IGS的精密卫星轨道和地球定向参数；同时估计卫星轨道误差和测站坐标，即松弛IGS轨道；截止高度角为15，历元间隔为30s；对流层折射根据标准大气模型用萨斯坦莫宁（Saastamoinen）模型改正；卫星钟差的模型改正（用广播星历中的钟差参数）；接收机钟差的模型改正（用根据伪距观测值计算出的钟差）；卫星和接收机天线相位中心改正，卫星与接收机天线L1、L2相位中心偏差采用GAMIT软件的设定值；测站位置的潮汐改正以及光压模型BERNE（与SIO的全球解一致）。IGS跟踪站（NER测站）坐标约束为0005、0005、0010m；GPS监测站（NER测站）坐标约束为9999、9999、9999m。在完成基线计算之后，进行GLOBK平差计算。计算中，采用ITRF08框架作为网平差的基础，主要约束为：IGS核心站约束，共采用44个测站，每站坐标的约束为0003、0003、0010m；轨道约束，与SIO的全球解H文件联合求解，放松轨道至20m。利用GAMIT解算时，首先需要建立工程目录（project），然后准备GPS数据文件（rinex文件），导航文件（n文件）和星历文件（sp3文件），接下来更新表文件（tables），并建立表文件目录与该工程的连接，最后编辑控制文件包括处理控制文件（process.defaults）、站处理控制文件（site.defaults）、时段控制表（sestbl.）、站控制表（sittbl.）、站信息文件（）和站坐标文件（lfile.）及其准备参考框架文件（itrf08.apr），准备工作完成之后，即可运行GAMIT的批处理命令进行解算。在解算过程中，GAMIT通过makexp和makex准备数据，用fixdrv模块生成批处理文件，用arc模块进行卫星轨道积分，用model模块计算观测值残差和偏导数，用autcln模块自动编辑数据，最后用solve模块进行最小二乘分析计算参数。利用GLOBK解算时，首先须将GAMIT的松弛解H文件由ASCII格式转换成二进制格式，然后准备好先验坐标文件、地球极移文件、计算控制文件等，利用glred和glorg以批处理形式进行解算。4结果分析GAMIT计算得到的单天解标准化均方差（normalizedrms，NRMS）是衡量单天解质量的重要指标之一。据国内外GPS数据处理经验，NRMS值一般应小于03mm，若标准化均方差太大，则说明处理过程中周跳可能没有完全得到修正。本研究数据时间范围为2011年年积日为第199天至第228天，检查q文件后，处理30d数据的标准化均方差基本符合解算要求2528（表1）。在GAMIT得到单天解后，可以得到卫星与测站的单程均方根残差，该指数可以在单天解总结文件（sh_gamit_ddd.summary）中查看，并且在用sh_make_sky_gifs或者批处理命令进行批处理生成每个测站的天空图时，对均方根残差有一个直观的可视化呈现（图3）。其中，0对应子午线。另外，在评价GAMIT解算过程中，可以从GAMIT预处理结算中LC相位残差得到站点的矫正模型，并且得到一些高度角与残差的关系，以KHUF测站为例，均方残差为488mm，基本符合测站的要求（图4、5）。在用GLOBK结算时，也可以对站点的信息进行统计分析，检验GAMIT处理结果，用gdl文件及其GMT进行画图处理。图6左边纵轴表示GPS站点名称，右边纵轴表示站点个数，方形代表单天解的个数（即每个站点的观测时段，一天对应一个观测文件）。利用GLOBK软件对观测期间各测点的点位重复性进行计算，并得到每个测点的时间序列，结果显示测点的点位重复性北向分量精度在3mm以内，东向分量在3mm以内，垂向分量在6mm以内。选取观测时间最长的DASH站点为代表，从中可以看出，该站点水平方向北向分量精度为1.04mm，东向分量为1.01mm，垂向分量为441mm。这体现了该站点在水平方向具有很好的重复精度（图7）。GPS观测的垂直方向敏感度很差，误差较大，因此图（c）中误差棒未标注多天解利用GLOBK软件对观测期间所有测点进行联合解算，并将结果归算到ITRF08参考框架下2934。多天整体解利用GAMIT软件30d（年积日为第199天至第228天）的处理结果H文件，结合SIO的全球解H文件，并且将30d的单天解进行合并，生成合并解GLX文件，利用全球stab_site文件及其core_site文件联合itrf08.apr文件，将其结果归算到ITRF08参考框架中，再用GLOBK命令作综合处理，通过平差计算得到最终点位和基线结果，结果归算到2011年8月16日，点位精度水平分量为23mm，高程分量为37mm（表2）。5结语2011年7月至8月，对帕米尔高原及其邻近地区进行了GPS观测剖面的布设和首期GPS观测，这是2条真正意义上的GPS观测剖面。笔者利用GAMIT/GLOBK软件对首期观测数据进行处理，引入周边8个IGS站点参与联合解算，并将结果归算到ITRF08下，解算得到的准观测值均满足所须精度要求，其中所有站点的标准化均方差均小于03mm。在观测区域单天解的基础上，利用GLOBK合并全球解，并且进行全区域的综合平差，得到GPS站点的精确坐标，结果显示所有站点的点位水平分量为12mm，高程分量为37mm。有必要在未来继续进行23期GPS观测，完善现有的GPS观测资料，帮助认识帕米尔高原的运动学特征和青藏高原隆升机制。QOCA（QuasiobservationCombinationAnalysis）学习班的董大南老师、中国地震局地震预测研究所武艳强副研究员和中国地震局地壳应力研究所田云锋副研究员在GAMIT/GLOBK软件处理中提供了指导，在此一并致谢。KH*2D参考文献：References：KH*2D1党亚民.GPS和地球动力学进展J.测绘科学，2004，29（2）：7779.DANGYamin.ProgressinGPSandItsGeodynamicsJ.ScienceofSurveyingandMapping，2004，29（2）：7779.2宋淑丽，朱文耀，廖新浩.GPS应用于地球动力学研究的进展J.天文学进展，2003，21（2）：95112.SONGShuli，ZHUWenyao，LIAOXinhao.RecentProgressofGPSinGeodynamicsJ.ProgressinAstronomy，2003，21（2）：95112.3王敏，沈正康，董大南.非构造形变对GPS连续站位置时间序列的影响和修正J.地球物理学报，2005，48（5）：10451052.WANGMin，SHENZhengkang，DONGDanan.EffectsofNontectonicCrustalDeformationonContinuousGPSPositionTimeSeriesandCorrectiontoThemJ.ChineseJournalofGeophsysics，2005，48（5）：10451052.4王敏，张祖胜，许明元，等.2000国家GPS大地控制网的数据处理和精度评估J.地球物理学报，2005，48（4）：817823.WANGMin，ZHANGZusheng，XUMingyuan，etal.DataProcessingandAccuracyAnalysisofNational2000GPSGeodeticControlNetworkJ.ChineseJournalofGeophsysics，2005，48（4）：817823.5张培震，王琪，马宗晋.中国大陆现今构造运动的GPS速度场与活动地块J.地学前缘，2002，9（2）：430441.ZHANGPeizhen，WANGQi，MAZongjin.GPSVelocityFieldandActiveCrustalBlocksofContemporaryTectonicDeformationinContinentalChinaJ.EarthScienceFrontiers，2002，9（2）：430441.6张培震，王琪，马宗晋.青藏高原现今构造变形特征与GPS速度场J.地学前缘，2002，9（2）：442450.ZHANGPeizhen，WANGQi，MAZongjin.GPSVelocityFieldandActiveCrustalDeformationinandAroundtheQinghaiTibetPlateauJ.EarthScienceFrontiers，2002，9（2）：442450.7张永志，胡斌，陈宏宾.应用GPS观测青藏高原东北缘应力场变化J.长安大学学报：地球科学版，2003，25（4）：8084.ZHANGYongzhi，HUBin，CHENHongbin.VariationofStressComputedfromGPSMeasurementinNortheasternMarginofQinghaiTibetPlateau，ChinaJ.JournalofChanganUniversity：EarthScienceEdition，2003，25（4）：8084.8王琪，张培震，马宗晋.中国大陆现今构造变形GPS观测数据与速度场J.地学前缘，2002，9（2）：415429.WANGQi，ZHANGPeizhen，MAZongjin.GPSDatabaseandVelocityFieldofContemporaryTectonicDeformationinContinentalChinaJ.EarthScienceFrontiers，2002，9（2）：415429.9王琪，丁国瑜，乔学军，等.天山现今地壳快速缩短与南北地块的相对运动J.科学通报，2000，45（14）：15431547.WANGQi，DINGGuoyu，QIAOXuejun，etal.RapidCrustalShorteningandRelativeMotionofCrustalBlocksintheTianshanJ.ChineseScienceBulletin，2000，45（14）：15431547.10杨少敏，李杰，王琪.GPS研究天山现今变形与断层活动J.中国科学：D辑，2008，38（7）：872880.YANGShaomin，LIJie，WANGQi.TheDeformationPatternandFaultRateintheTianshanMountainsInferredfromGPSObservationsJ.ScienceinChina：SeriesD，2008，38（7）：872880.11杨少敏，王琪，游新兆.中国现今地壳运动GPS速度场的连续变形分析J.地震学报，2005，27（2）：128138.YANGShaomin，WANGQi，YOUXinzhao.NumericalAnalysisofContemporaryHorizontalTectonicDeformationFieldsinChinafromGPSDataJ.ActaSeismologicaSinica，2005，27（2）：128138.12牛之俊，游新兆，杨少敏.利用GPS分析天山现今地壳形变特征J.大地测量与地球动力学，2007，27（2）：19.NIUZhijun，YOUXinzhao，YANGShaomin.AnalysisofContemporaryCrustalDeformationCharacteristicswithGPSDataofTianshanMountainJ.JournalofGeodesyandGeodynamics，2007，27（2）：19.13REIGBERC，MICHELGW，GALASR，etal.NewSpaceGeodeticConstraintsontheDistributionofDeformationinCentralAsiaJ.EarthandPlanetaryScienceLetters，2001，191（1/2）：157165.14BANERJEEP，BURGMANNR.ConvergenceAcrosstheNorthwestHimalayafromGPSMeasurementsJ.GeophysicalResearchLetters，2002，29（13）：1652.15MOHADJERS，BENDICKR，ISCHUKA，etal.PartitioningofIndiaEurasiaConvergenceinthePamirHinduKushfromGPSMeasurementsJ.GeophysicalResearchLetters，2010，37（4）：L04305.16ZUBOVICHAV，WANGXQ，SCHERBAYG，etal.GPSVelocityFieldfortheTianShanandSurroundingRegionsJ.Tectonics，2010，29（6）：TC6014.17陈杰，李涛，李文巧，等.帕米尔构造结及邻区的晚新生代构造与现今变形J.地震地质，2011，33（2）：241259.CHENJie，LITao，LIWenqiao，etal.LateCenozoicandPresentTectonicDeformationinthePamirSalient，NorthwesternChinaJ.SeismologyandGeology，2011，33（2）：241259.18FUBH，NINOMIYAY，GUOJM.SlipPartitioningintheNortheastPamirTianShanConvergenceZoneJ.Tectonophysics，2010，483（3）：344364.19孙文斌，和跃时，常征，等.帕米尔兴都库什地区板块俯冲及其应力状态J.地震地质，2009，31（2）：207217.SUNWenbin，HEYueshi，CHANGZheng，etal.ThePlateSubductionandStressStateinthePamirHinduKushRegionJ.SeismologyandGeology，2009，31（2）：207217.20王敏，张培震，沈正康，等.全球定位系统（GPS）测定的印尼苏门达腊巨震的远场同震地表位移J.科学通报，2006，51（3）：365368.WANGMin，ZHANGPeizhen，SHENZhengkang，etal.FarfieldCoseismicDisplacementsAssociatedwiththeGreatSumatraEarthquakesofDecember26，2004andMarch29，2005ConstrainedbyGlobalPositioningSystemJ.ChineseScienceBulletin，2006，51（3）：365368.21王敏，沈正康，甘卫军.GPS连续监测鲜水河断裂形变场动态演化J.中国科学：D辑，2008，38（5）：575581.WANGMin，SHENZhengkang，GANWeijun.GPSMonitoringofTemporalDeformationoftheXianshuiheFaultJ.ScienceinChina：SereisD，2008，38（5）：575581.22王阎昭，王恩宁，沈正康，等.基于GPS资料约束反演川滇地区主要断裂现今活动速率J.中国科学：D辑，2008，38（5）：582597.WANGYanzhao，WANGEnning，SHENZhengkang，etal.GPSconstrainedInversionofPresentdaySlipRatesAlongMajorFaultsoftheSichuanYunnanRegionJ.ScienceinChina：SereisD，2008，38（5）：582597.23张培震，王敏，甘卫军，等.GPS观测的活动断裂滑动速率及其对现今大陆动力作用的制约J.地学前缘，2003，10（增1）：8192.ZHANGPeizhen，WANGMin，GANWeijun，etal.SlipRatesAlongMajorActiveFaultsfromGPSMeasurementsandConstraintsonContemporaryContinentalTectonicsJ.EarthScienceFrontiers，2003，10（S1）：8192.24张晓亮，张希，王双绪，等.GPS资料揭示的汶川80级地震前地壳动态演化特征J.大地测量与地球动力学，2011，31（4）：15.ZHANGXiaoliang，ZHANGXi，WANGShuangxu，etal.CharacteristicsofCrustalDynamicEvolutionRevealedbyGPSObservationBeforeWenchuan8.0EarthquakeJ.JournalofGeodesyandGeodynamics，2011，31（4）：15.25孟国杰，赵承坤，顾国华，等.中国地壳运动观测网络首期观测数据处理J.地震，2001，21（1）：6368.MENGGuojie，ZHAOChengkun，GUGuohua，etal.DataProcessionoftheFirstEpochMeasurementoftheCrustalMovementObservationNetworkinChinaJ.Earthquake，2001，21（1）：6368.26牛之俊，王敏，孙汉荣，等.中国大陆现今地壳运动速度场的最新观测结果J.科学通报，2005，50（8）：839840.NIUZhijun，WANGMin，SUNHanrong，etal.ContemporaryVelocityFieldofCrustalMovementofChineseMainlandfromGlobalPositioningSystemMeasurementsJ.ChineseScienceBulletin，2005，50（8）：839840.27田云锋，沈正康.GPS坐标时间序列中非构造噪声的剔除方法研究进展J.地震学报，2009，31（1）：6881.TIANYunfeng，SHENZhengkang.ProgressonReductionofNontectonic