红河断裂带闭锁程度和滑动亏损分布特征研究

利用1999~2013年青藏高原东南缘的GPS速度场观测数据,采用DEFNODE程序反演红河断裂带走滑速率、三维闭锁程度和滑动亏损分布。反演结果表明:红河断裂带北、中、南段右旋走滑速率分别约为(5.9±1.2)mm/a,(4.8±0.6)mm/a和(4.3±0.4)mm/a,在地表以下6km的闭锁程度分别为0.43,0.22和0.25,其滑动亏损速率分别为3.3mm/a,2.6mm/a和2.3mm/a,红河断裂带闭锁程度和应变积累程度都比较低,与近年来红河断裂带整体活动微弱的现状相吻合。     

An automatic cycle-slip processing method and its precision analysis

On the basis of analyzing and researching the current algorithms of cycle-slip detection and correction, a new method of cycle-slip detection and correction is put forward in this paper, that is, a reasonable cycle-slip detection condition and algorithm with corresponding program COMPRE (COMpass PRE-processing) to detect and correct cycle-slip automatically, compared with GIPSY and GAMIT software, for example, it is proved that this method is effective and credible to cycle-slip detection and correction in GPS data pre-processing.     

城市水资源GPS网的建立

结合利用CORS虚拟技术,阐述了城市水资源勘测GPS网的布设及数据处理方法,并对精度进行了分析,结果表明在确保成果精度和可靠性的前提下,CORS虚拟技术可以节省成本投入,提高生产效率,从而对丰富CORS服务具有重要意义。     

中性大气掩星标准反演技术

结合JPL和GFZ的中性大气标准反演算法，发展了上海天文台中性大气掩星标准反演算法流程．对CHAMP观测资料进行中性大气反演，并将反演结果与ECMWF(the European Center for Medium Range Weather Forecasting)的预报大气数据进行比较．作为讨论，指出将来的发展方向．     

利用GPS双频数据进行区域电离层TEC提取

目的 从利用GPS提取区域电离层总电子含量(total electron content,TEC)的基本原理出发,解决了伪距观测值优化以及硬件延迟(DCB)处理问题,并将提取的TEC信息与欧洲定轨中心(CODE)计算的全球电离层(GIM)模型内插值应用在单频精密单点定位中,进行电离层延迟改正实验。结果表明,利用本文提取的TEC值进行单频精密单点定位电(PPP)离层延迟改正时,点位精度能提高到0.2～0.4m左右,明显优于利用GIM内插值的改正精度。     

联合GPS与重力资料反演福建地区地壳运动

联合福建省近几年的GPS基本网、基准站观测数据和重力复测数据,反演福建块体的分块运动参数和主要活动断裂运动状态。结果显示,受欧亚板块、菲律宾海板块和太平洋板块相互碰撞的影响,福建地区块体的运动样式大体呈现顺时针方向的运动,造成北东向断裂作左旋走滑运动,而北西向断裂以右旋走滑为主。其中,长乐-诏安断裂带、邵武-河源断裂带和永安-晋江断裂带是块体内最活跃的断裂带,其左旋走滑速率分别为0.7±1.3、0.7±1.7和0.6±1.8mm/a,九龙江下游断裂带具有0.5±1.7mm/a的右旋走滑速率。     

利用GPS水准控制山区高程——以陕西丁家林金矿区为例

应用GPS水准控制山区高程的基本原理是当GPS点布设成一定区域时，可以用数学曲面拟合法求待定点的正常高，根据测区中已知点的平面坐标或大地坐标和高程异常值，用数字拟合法求出该区似大地水准面，得出等求点的正常高。在陕西丁家林金矿区实验结果表明，该方法与传统的三角高程控制法相比，无论在平原或在山区都能获得较好精度。因此，用GPS水准可替代几何水准。     

巴颜喀拉块体东缘形变及九寨沟地震孕震环境数值分析

2017年8月8日四川发生九寨沟M7.0地震,是继2008年汶川M8.0地震后发生在巴颜喀拉块体东部的又一强震.现今GPS速度观测数据显示,2008年汶川地震前后的1999-2007年和2011-2016年两个时间段内巴颜喀拉块体东部地表速度场存在明显的差异.本文以实际GPS速度观测资料为约束,构建三维有限元地球动力学模型,分别计算分析了两个时段内震源区及周边现今地壳形变、弹性应变能和应力积累特征,进一步探讨汶川地震的发生对九寨沟地区变形及应力的影响.数值模拟结果显示,汶川地震之后（2011-2016年）巴颜喀拉块体东部的地壳形变、弹性应变能积累及应力积累速率均明显大于震前,增加量值达1.5~3倍;九寨沟地震发震断裂上库仑应力增长率在1999-2007年约为0.7 kPa·a^-1,2011-2016年间增至1.2 kPa·a^-1.上述结果表明,现今巴颜喀拉块体东部地壳应力积累过程有利于左旋走滑型九寨沟地震的发生,汶川地震的发生调整了区域应力状态,加速了九寨沟地震的孕育过程.

     

GEODETIC AND TELESEISMIC CONSTRAINTS ON SLIP DISTRIBUTION OF 2015 MW6.4 PISHAN EARTHQUAKE

On July 3^rd, 2015, a M_W6.4 earthquake occurred on Pishan County, Xinjiang, located in the front of western Kunlun thrust belt, which is the largest earthquake(M_W6.0~7.0)in the past 40 years in this region. In this study, we collected both the near-filed geodetic coseismic deformation observations including 4 GPS sites and one high-resolution ALOS-2 InSAR imagery, and far-field teleseismic P waveforms from 25 stations provided by IRIS/USGS, to invert the fault parameters(strike and dip)and coseismic rupture model of 2015 M_W6.4 Pishan earthquake. Using the finite fault theory, a non-linear simulated annealing algorithm was employed to resolve our joint inversion problem. The strike (120°~130°) and dip angle(35°~40°)of optimal models are different from that of some previous studies, and the dip change is strongly constrained by combined data than that of strike. In fixing the geometric parameters of optimal fault model, we also considered data weight(5)(geodetic data/teleseismic P waveforms)and constrained weight from moment and smooth factor(2.5). Clearly, our results indicate that the slip distribution mainly concentrates in the depth range from 9 to 16km and a length range of 20km along the strike direction, which is similar to the spatial distribution of the relocated aftershocks. The maximum slip is~95cm. The seismic moment release is 5.45×10¹⁸N·m, corresponding to M_W6.42. Compared with the single data set, geodetic data or teleseismic waveform, our joint inversion model could simultaneously constrain the seismic moment and slip distribution well, thus avoiding effectively a lower-resolution rupture distribution determined by teleseismic-only inversion and a bias released moment estimated by the geodetic-only inversion. Importantly, we should consider both the near-field geodetic data and far-field teleseismic data in retrieving the rupture model for accurately describing the seismogenic structure of active fault in western Kunlun region.     

Motive force of Qinghai-Tibet Plateau moving to east

A calculation formula on spherical pattern of Qinghai-Tibet plateau moving model is established. Tibet massif moves norward by east in speed of 28 mm/a, Ganshu-Qinghai massif moves to northeast in speed of 15 mm/a, Qomolangma Feng moves northward by a few east in speed of 35 ～42 mm/a. The low latitude perimeter is longer than the nigh latitude perimeter. When the Tibet massif moves northward, its latitude perimeter must be contracted and the Tibet massif must move eastward by Coriolis. Coriolis force is inertial in earth rotation. It makes the fall body turning to east and the rising block turning westward. In the Northern Hemisphere, it makes the northward body turning to east and the southward block turning to west.This is the reason why the tectonic zones of western Pacific are different from those of eastern Pacific.