上海VLBI站相对于欧亚板块运动的测定与讨论

利用 1 988～ 1 996年共 8年的上海天文台 VLBI站与国际 VLBI站基线长度变化的精确测定结果 ,重新估计了上海 VL BI站的形变速率 ,证实了基于不同地球参考架系统对估计的上海站的运动有明显差异 ;基于最新的 ITRF96地球参考架和 Sillard等 (1 998)欧亚板块的欧拉矢量 ,估计得到上海站的地壳垂直形变速率为 -1 .95± 1 .0 2 mm/a呈下降趋势 ,水平形变速率为 1 3.5 5± 1 .2 0 mm/a,方位 77.5°± 4.5°;估计得到乌鲁木齐 VL BI站的地壳垂直形变速率为 -6 .0 8± 2 .77mm/a呈下降趋势 ,水平形变速率为 2 1 .87± 2 .5 2 mm/a,方位 45 .5°±6 .6°;进一步分别基于 NNR-NUVEL -1 A地球板块运动模型和 Zhang等 (1 999)的欧亚板块的欧拉矢量 ,得到与上述较一致的结果     

联合卫星重力、卫星测高和海洋资料研究全球海平面变化

利用GRACE、卫星测高和海洋实测温盐数据,探讨了2003~2012年间全球海平面、比容海平面和海水质量等的变化特征,并讨论了南极冰盖和格陵兰冰盖消融对全球海平面变化的影响。全球海平面整体呈上升趋势,上升速度为2.72±0.07 mm/a,且存在显著的空间分布特征。全球海平面、比容海平面和海水质量等的变化还具有显著的季节性特征,其中全球海平面变化的年周期振幅为4.6±0.3 mm。使用经验正交函数分析(EOF)得到全球海平面和比容海平面的季节性变化在南北半球存在显著的差异,但海水质量季节性变化不存在这种差异。南极冰盖和格陵兰冰盖的消融速率分别为-75.7±12.3 Gt/a和-124.1±2.9 Gt/a,对海平面的长期趋势项贡献分别为0.21±0.03 mm/a和0.34±0.01 mm/a,仅占全球海水质量增加速度1.80±0.10 mm/a的12%和19%,总计占31%,因此,两极冰盖质量消融并不是2003-2012年间海水质量增加的最主要因素。     

顾及GRACE季节影响的华北平原水储量变化反演

基于移去-恢复方法,利用GRACE卫星重力数据研究了华北平原2003年1月—2014年6月的陆地水变化,提出了一种顾及季节影响尺度因子计算方法,用于减小GRACE后处理误差。将本文方法所得结果与水文模式和降水模型结合进行比较分析,结果表明：2008年之前华北平原发生比较明显的干旱现象,陆地水和地表水分别以（-7.9±2.4）mm/a和（-7.3±2.8）mm/a的速度下降,但地下水仅以（-0.6±1.4）mm/a的速度减少;2008年之后,陆地水和地表水分别以（4.3±1.3）mm/a和（10.9±2.1）mm/a的速度上升,期间地下水的超采严重,平均以（-6.5±1.2）mm/a的速度下降,但下降速度以0.9 mm/a²的趋势减缓。研究期间内,华北平原的陆地水、地表水和地下水整体分别呈现（-2.0±0.6）、（2.9±0.7）和（-4.8±0.7）mm/a的变化趋势。最后,利用交叉小波谱分析了GRACE滞后于降水的现象。研究结论表明,降水和地下水开采是影响华北平原水储量变化的两大决定性因素。     

基于GPS数据的川滇地区地壳水平运动研究

利用川滇地区2011-2013年年积日206-220 d共计45 d的陆态网络站点的GPS观测数据,得到川滇地区基于ITRF2008的两期水平运动绝对速度场;经分析川滇地区近3年地壳水平运动较为显著,整体呈SE向且有增强趋势。利用单点基准得到相对水平运动速度场,相对于LHAZ站,川滇地区整体上呈现显著SW向运动,顺时针旋转向印度板块。相对于区域内部站点SCSM和YNYM,川滇地区不同区域呈现出相异的运动趋势。     

中国近海1992～1998海平面变化监测与分析

利用Topex/Poseidon和ERS-1卫星测高资料采用“共线”法计算出1992-10～1998-06中国近海海域海平面及其变化。在扣除T/P卫星测高仪的零点漂移影响后,发现与全球海平面上升率(+2.1±1.3)mm/a相比,不同海域的海平面变化趋势大不相同,黄、东、南海的海平面年变化率分别为：(+3.44±0.61)mm/a,(+3.12±0.47)mm/a,(-1.41±0.48)mm/a。从海平面变化异常中可以看出1993,1994,1997～1998年3次El Nino异常对中国近海海域海平面的影响是南海海域最大,东海次之,黄海最小。除了对海平面进行传统的频谱分析外,还进行了多分辨率的小波分析,还发现在上述3个海域中除了年周期变化较为稳定外,半年及季节(100 d)周期项存在着时间漂移。此外,在上述3个海域还存在着明显的两个月(62 d)周期的变化,其激发原因在此做了初步探讨。     

川滇断块的现代地壳运动

一、引言川滇地区现代地壳运动的水平运动一直是个有争议的问题。它地处我国东西两种构造域的分界带——约１００°Ｅ～１０５°Ｅ,西部构造运动通过这一带影响东部,故构造运动比较活跃、复杂,且地震频发。就构造分区而言,川滇断块即扬子地块西部边缘,它不是独立的岩...     

四、海洋测量学

CH20002412 中国近海1992～1998海平面变化检测与分析/王海瑛(中国科学院测量与地球物理研究所)…∥测绘学报.—2000,29(增刊).—32～37利用 Topex/Poseidon 和 ERS-1卫星测高资料采用“共线”法计算出1992.10～1998.06中国近海海域海平面及其变化。在扣除 T/P 卫星测高仪的零点漂移影响后,发现与全球海平面上升率(+2.1±1.3)mm/a 相比,不同海域的海平面变化趋势大不相同,黄、东、南海的海平面年变化率分别为(+3.44±0.61)mm/a.(+3.12±0.47)mm/a,(-1.41±0.48)mm/a。从海平面变化异常中可以看出1993.1994,1997～1998年3次 EINino 异常对中国近海海域海平面的影响是南海海域最大,东海次之,黄海最小。除了对海平面进行传统的频谱分析外,还进行了多分辨率的小波分析,还发现在上述3个海域中除了年周期变化较为稳定外,半年及季节(100d)周期项存在着时间漂移。此外,在上述3个海域还存在着明显的两个月(62d)周期的变化,其激发原因在此做了初步探讨。图3表2参7     

广西区域现今地壳运动与应变的GPS监测分析

利用GAMIT/GLOBK软件分析了1998-2010年间多期GPS观测数据,获得了ITRF2005框架下广西区域现今地壳运动三维速度场和应变场参数,并对区域形变场和现今构造运动动力学机制进行了初步讨论。GPS监测结果显示:广西区域现今地壳运动随华南块体作南东东方向的逆时针构造运动,ITRF2005框架下水平方向平均速率为33.49 mm/a,优势方向为N 109.3° E;垂直方向平均速率为+1.30 mm/a,整体趋势是自桂东向桂西方向垂向速率逐步增加。现今应变状态以N 127.3°±4.8° E挤压为主,最大主压应变率-2.7±0.5×10^-9/a;兼有N 37.3°±4.8° E方向的拉张作用,最大主张应变率1.8±0.2×10^-9/a。本文研究结果为描述广西区域现今地壳运动提供了一个更加精细的三维速度场。     

川滇地块现今地壳运动速度场解算与分析

利用GAMIT10.4和BERNESE5.2软件,分别对川滇地区中国地壳运动观测网络51个GNSS基准站2010—2014年的观测数据进行处理。两者解算所得单日解基线和坐标时间序列离散度优于10mm,且变化趋势一致,吻合度较高。基线重复率水平向优于1+4×10-9 L;垂向优于4+6×10-9 L(单位mm,其中L表示基线长度,单位m)。利用PL+VW噪声模型进行了速度估计,获取了川滇地区相对欧亚板块地壳运动速度场。结果显示:除滇西南地区呈现出由西向东的增大趋势外,川滇块体其它区域地壳运动特征由北往南、由西往东呈逐渐衰减趋势;运动方向持续围绕喜马拉雅东结点作顺时针旋转,且仍然存在着南北向强烈的挤压特征。以安宁河-则木河断裂、小江断裂和红河断裂为界,选用了两侧能反映震间断裂构造特征的GPS测站,分析了断裂两侧的速度差异。并利用平行断层的速度剖面,拟合反正切函数变化趋势,获取了断裂带可能的变形宽度。     

利用多面函数拟合和欧拉矢量法建立福建省GPS速度场模型

通过处理福建省卫星导航定位基准服务系统(FJCORS)2013-2015年观测资料,获得GPS站实测运动速度值,在此基础上,分别利用多面函数拟合法和欧拉矢量法构建了福建省及周边区域GPS速度场模型．内符合精度方面,多面函数拟合的速度精度南北分量为1.06 mm／a、东西分量为2.00 mm／a;欧拉矢量法得到的速度精度南北分量为1.36 mm／a、东西分量为1.22 mm／a．外符合精度方面,两种方法模拟得到的检核点速度值与实测速度值的较差均在±2.28 mm／a以内．在福建省陆域内,两种方法得到的速度场具有高度的一致性,运动速度均为35 mm／a左右,运动方向为SE-E向．      

利用GPS数据反演中国红河断裂带活动特性

为了研究中国红河断裂带近年的活动特性，利用块体模型对红河断裂带周围块体内1999-2014年的GPS水平速度场进行建模，反演得到块体的刚体旋转、内部应变率以及断层的滑动速率与三维闭锁程度。结果表明，考虑内部应变的块体模型可以显著提升速度场的拟合效果，反演得到的红河断裂带滑动速率与二维速度剖面的结果吻合较好。红河断裂带的闭锁程度大致呈现中段闭锁较强、南段与北段闭锁较弱的对称分布，南段、中段与北段在地表以下15 km深度范围内的断层闭锁程度分别约为0.34~0.52、0.67~0.75、0.26~0.53，右旋滑动亏损速率分别为1.7、3.0、1.4 mm/a。利用最小二乘配置方法计算了研究区域的应变场分布，发现红河断裂带从南至北由压缩转换为拉张，中段处于压缩与拉张的交界带，剪应变率低于南北两段，这可能是由于中段闭锁较强、断层错动量较小造成的。     

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

利用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,红河断裂带闭锁程度和应变积累程度都比较低,与近年来红河断裂带整体活动微弱的现状相吻合。     

小江断裂带地壳运动特征及地震危险性研究

利用1999~2013年间川滇地区的GPS速度场观测数据和块体模型研究小江断裂带地壳运动特征及地震危险性。结果表明:小江断裂带闭锁程度较高,应变量积累较大,在地表以下6km的闭锁程度约为0.94,滑动亏损速率约为(9.2±0.4)mm/a,整条断裂带左旋走滑速率约为(10.8~12.1)mm/a。结合小江断裂带上的历史地震数据估算,小江断裂带未来地震最大震级为7.3级,与近年来断裂带活动较为活跃相吻合。     

Plate Motion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements

We analyse geodetically estimated deformation across the Nepal Himalaya in order to determine the geodetic rate of shortening between Southern Tibet and India, previously proposed to range from 12 to 21 mm yr^?1. The dataset includes spirit-levelling data along a road going from the Indian to the Tibetan border across Central Nepal, data from the DORIS station on Everest, which has been analysed since 1993, GPS campaign measurements from surveys carried on between 1995 and 2001, as well as data from continuous GPS stations along a transect at the logitude of Kathmandu operated continuously since 1997. The GPS data were processed in International Terrestrial Reference Frame 2000 (ITRF2000), together with the data from 20 International GNSS Service (IGS) stations and then combined using quasi- observation combination analysis (QOCA). Finally, spatially complementary velocities at stations in Southern Tibet, initially determined in ITRF97, were expressed in ITRF2000. After analysing previous studies by different authors, we determined the pole of rotation of the Indian tectonic plate to be located in ITRF2000 at 51.409±1.560° N and ?10.915±5.556°E, with an angular velocity of 0.483±0.015°. Myr^?1. Internal deformation of India is found to be small, corresponding to less than about 2 mm yr^?1 of baseline change between Southern India and the Himalayan piedmont. Based on an elastic dislocation model of interseismic strain and taking into account the uncertainty on India plate motion, the mean convergence rate across Central and Eastern Nepal is estimated to 19±2.5 mm yr^?1, (at the 67% confidence level). The main himalayan thrust (MHT) fault was found to be locked from the surface to a depth of about 20 km over a width of about 115 km. In these regions, the model parameters are well constrained, thanks to the long and continuous time-series from the permanent GPS as well as DORIS data. Further west, a convergence rate of 13.4±5 mm yr^?1, as well as a fault zone, locked over 150 km, are proposed. The slight discrepancy between the geologically estimated deformation rate of 21±1.5 mm yr^?1 and the 19±2.5 mm yr^?1 geodetic rate in Central and Eastern Nepal, as well as the lower geodetic rate in Western Nepal compared to Eastern Nepal, places bounds on possible temporal variations of the pattern and rate of strain in the period between large earthquakes in this region.     

利用局域复测GPS网研究中国大陆块体现今地壳运动速度场

利用1992年以来我国进行的6项不同规模高精度GPS地壳监测网的资料，采用统一的数据处理策略，重新进行了数据处理，并统一进行了平差，归算到相同的ITRF97坐标框架基准之下，获得了中国前覆最广、测站最多的一个较为完整的、框架统一的GPS地壳运动图像和速度场，并对结果作了初步的地学解释。     

基于GPS观测数据的云南地区地壳形变特征研究

利用1999~2013年间多期GPS观测数据,将云南地区分为5个活动块体,研究区域现今水平运动特征、应变特征以及各块体的运动特征,反演主要活动断裂带的走滑速率。结果表明,区域水平运动具有北强南弱的特征,小江断裂带、红河断裂带对云南地区构造活动具有明显控制作用,改变了断裂带两侧的运动特征;区域应变具有明显的差异,西部以拉张应变为主,东部小江断裂带周边地区拉张和挤压应变均衡且为强应变区,滇东地区应变较小;各块体运动方向与速率存在差异,滇东块体为逆时针旋转,其余块体顺时针旋转;基于块体模型反演的主要断裂带的走滑速率与地质方法获取的走滑速率具有较好的一致性。     

中国大陆GPS速度场的球面小波模型及多尺度特征分析

利用1999—2009年间中国大陆共1068个GPS站点在东方向、北方向的速度值,采用DOG球面小波多尺度分析方法,建立了中国大陆东方向、北方向多尺度速度场。球面小波模型的尺度主要根据观测站点的密度来确定,利用检核点上的已知速度与模型速度之间的均方差来评定模型的精度。利用球面小波模型可以更加清晰地表示速度场的大尺度特征和复杂的局部变化特征。站点稠密区域,模型在东方向、北方向上的精度分别为±0.95mm/a、±0.97mm/a,稀疏区域对应的精度分别为±1.32mm/a和±1.30mm/a。     

An assessment of Bernese GPS software precise point positioning using IGS final products for global site velocities

We assess the use of precise point positioning (PPP) within the Bernese GPS software (BSW) Version 5.0 over the period from 2000 to 2005. In our strategy, we compute a set of daily PPP solutions for international GNSS service (IGS) reference frame (IGb00) sites by fixing IGS final satellite orbits and clock products, followed by a Helmert transformation of these solutions into ITRF2000, forming a set of continuous position time series over the entire time span. We assess BSW PPP by comparing our set of transformation parameters with those produced by the IGS analysis centre coordinator (ACC) and our position time series with those of the Jet Propulsion Laboratory (JPL) and the Scripps Orbit and Permanent Array Centre at the Scripps Institute of Oceanography (SIO). The distributions of the north (N), east (E) and up (U) daily position differences are characterized by means and SD of +2.2 ± 4.8, −0.6 ± 7.9 and +4.8 ± 17.3 mm with respect to JPL, and of +0.1 ± 4.4, −0.1 ± 7.4 and −0.1 ± 11.8 mm with respect to SIO. Similarly, we find sub-millimetre mean velocity differences and SD for the N, E and U components of 0.9, 1.5 and 2.2 mm/year with JPL, and of 1.2, 1.6 and 2.3 mm/year with SIO. A noise analysis using maximum likelihood estimation (MLE) shows that when estimating global site velocities from our position time series, the series need to be on average up to 1.3 times longer than those of JPL and SIO, before an uncertainty of less than 0.5 mm/year is obtained.