利用GPS与GRACE监测陆地水负荷导致的季节性水平形变：以喜马拉雅山地区为例

地表陆地水负荷变化是引起重力场和地壳形变呈现季节性特征的主要因素,并且能够利用地表及空间大地测量技术对其进行有效的监测.本文通过对质量负荷形变效应的理论模拟,描述了水平分量的形变指向以及垂直与水平分量的幅值比可以提高对负荷区域的辨别程度,并且联合GPS坐标时间序列及GRACE模型对喜马拉雅山地区的季节性负荷形变进行了详细对比分析,研究结果显示两者垂直分量的季节性变化具有较好的一致性,且GPS周年项幅值要大于GRACE.而由GRACE解算得到的水平分量结果表明该地区季节性形变主要受东南亚及印度东北部地区的陆地水负荷控制,位于喜马拉雅山地区多数GPS台站的垂直分量及北向分量的初相位与GRACE模型解算结果相近,而部分GPS台站的东向分量与GRACE模型存在明显不同,由此导致GPS与GRACE监测到的形变指向存在差异.通过对GRACE估算精度以及GPS垂直与水平分量幅值比的深入分析,发现GPS对局部周边地区的河流、谷地及农田灌溉等负荷变化造成的形变效应较为敏感,而GRACE由于截断阶次及平滑滤波等影响因素,不仅造成在水平分量上的分辨率远低于垂直分量,而且整体估算精度要低于GPS观测得到的形变信息.     

Repeated Geodetic Measurements in the West Bohemia Seismoactive Region: Period 1993 - 1996

The western part of the Bohemian Massif is characterized by repeated occurrences of intraplate earthquake swarms. To study surface deformations of this anomalous region, a network covering about 2000 square kilometres for repeated geodetic measurements was established in 1993 - 1994. The positions of the individual points of the network were carefully picked with respect to local tectonic structure and earthquake foci distribution. GPS and precise levelling measurements were performed 1-2 times a year. The GPS data were processed by Bernese GPS software. No tendency to any displacement - either horizontal or vertical - of geological blocks was derived from the geodetic data for the period 1993 - 96. Only displacements of less than 5 mm/year in average could occur in the whole region in that period; larger displacements would have been revealed by our measurements.     

GPS与PS-InSAR联网监测的台湾屏东地区三维地表形变场

联合高精度的GPS水平位移观测和高密度的PS-InSAR雷达视线位移测量,实现地表三维形变的精确反演.本文在准确计算卫星轨道方位角基础上,使用GPS观测位移与星载雷达LOS方向形变的投影转换模型,将雷达LOS方向形变转换为垂直方向位移,并基于地面GPS与SAR影像PS目标联合构建形变监测网,采用参数平差算法估计区域地表形变场.以地质构造活动极其活跃的台湾岛及其西南屏东高雄地区为例,联合屏东地区48个GPS监测台站与雷达PS目标,监测该地区从1995-1999年间由于板块构造挤压运动和地下水抽取导致的三维地表形变.结果表明,该地区年均水平位移量为向西30~50 mm/a,高雄沿海地区发生明显的逆时针西偏南的逐渐增大的水平位移;垂直位移为屏东平原南部呈现-10 mm/a~-15 mm/a的地面沉降,而平原北部和高雄地区呈现约+5 mm/a~+10 mm/a的地面抬升.     

日本M9.0级巨震对山东地区地壳活动的影响研究

以山东地壳运动GPS观测网为基础,结合周边IGS提供的全球GPS观测资料,计算了2011年3月11日发生的日本9.0级巨震对山东所处地块产生的影响,并给出了该地区站点的真实运动轨迹和运动方式.结果显示:山东各GPS基准站受地震波影响产生的震时最大振幅超过0.2 m,持续时间约160 s;同震位移方向为东东南向,位移量最大超过10 mm.日本9.0级地震对山东所在地块的影响主要是在水平方向,对垂直方向影响不大,它对山东主要断裂带-沂沭断裂带起到了拉张的作用,且N段幅度明显大于S段.根据地震前后站点时间序列的发展趋势来看,近期山东所处地块处于向西回调运动状态,推测震后3～5个月左右将恢复到震前的轨迹,为本地区的震情判定提供了可靠依据.     

利用GPS观测结果对我国地壳水平形变强度的分析

介绍利用中国地壳运动观测网络区域网(包括基本网与基准网)1999年与2001年2期GPS观测所获得的中国大陆地壳水平位移速率结果, 研究块体的划分、 块体位移与块体变形, 提出了用位移(速率)离散度分析块体或局部水平变形强度的方法, 分析了中国大陆现今地壳运动的基本特征及2001年11月14日昆仑山口西8.1级大地震前的地壳水平运动。     

印尼8.7级大震前后GPS观测站的地壳水平与垂直位移时间序列结果

利用IGS及中国地壳运动观测网络GPS连续观测站6年多时间的观测结果, 以中国大陆东部的稳定点组作为水平位移与垂直位移解的参考基准, 得到了2004年12月26日印尼苏门答腊8.7级与2005年3月28日8.5级大震前后GPS观测站的地壳水平及垂直位移时间序列的结果。 尽管所采用的GPS站点少, 且分布范围大, 但仍可看到在如此大的地震震前、 同震与震后地壳运动的特征, 为今后观测研究提供了十分有意义的结果。 文中所采用稳定点组基准有效消除水平位移场的平移与旋转, 而局部椭球面的不平性对计算区域地壳运动结果, 特别是对垂直位移的影响很小, 说明了位移解所描述的大范围区域位移场是合理的。 这两次地震的同震水平位移及震后的垂直位移影响量级达数厘米, 范围达远离8.7级地震震中4500 km之外, 甚至更远。     

Horizontal crustal movement in Chinese mainland from 1999 to 2001

The paper introduces the horizontal crustal movement obtained from GPS observations in the regional networks (including the basic network and the fiducial network) of the Crustal Movement Observation Network of China (CMONOC) carried out in 1999 and 2001. This paper is characterized by the acquisition of the horizontal displacement velocities during the period from 1999 to 2001 at the observation stations in the regional networks with datum definition of a group of stable stations with small mutual displacements in east China. Based on the most detailed map of horizontal crustal movement in Chinese mainland, the division of blocks, their displacements and deformations are studied. An approach to analysis of the intensity of the horizontal crustal deformation is proposed. The general characteristics of the recent horizontal crustal movement in Chinese mainland and that before the Kunlunshan earthquake of M=8.1 on November 14, 2001 are analyzed. Foundation item: The National Development and Programming Project for Key Basic Research (95-13-03-07).     

2016年8月24日意大利佩鲁贾MW 6.2级地震同震位移场估计和震源滑动模型反演

利用改进的自动经验基线校正方法SMBLOC,对2016年8月24日意大利佩鲁贾MW 6.2级地震震中周围约60 km内的近场强震记录进行基线校正并尝试给出同震位移场,与GPS观测结果进行对比分析,分别独立和联合两种资料反演震源滑动模型,并根据震源模型进一步给出全空间预测位移场分布.研究结果表明:(1)两种不同的资料给出的水平位移场幅值均为cm级,且均表明断层的错动以正断为主.(2)两种同震位移场分别独立和联合反演所得的震源静态滑动范围基本一致,最大滑动均发生在震中东北侧,强震模型表现出明显的双事件特征,较大滑动分布在震中东北侧和东南侧,GPS模型在震中东南侧的滑动相对较小,其双事件特征不明显.两种模型的最大滑动量分别为0.96 m和0.86 m,较为一致,反演的矩震级均在MW 6.3左右.(3)根据震源滑动模型计算所得的佩鲁贾地震全空间预测的水平同震位移场中最大位移分布区域与震后报告中受灾严重的地区基本一致.表明在一定的条件下,利用SMBLOC方法解算震级较小的MW 6.0左右地震强震记录的同震位移场,并反演震源滑动模型具有一定的可行性,且其同震位移场和滑动模型可为震后灾害快速评估、救援力量分配、余震趋势判定等快速应急响应工作提供参考依据.     

Study on the Horizontal Deformation Strain Field in the Central and Northern Parts of Yunnan Province

Based on the horizontal deformation field and the strain field derived from the GPS data over the period of 1999～2001 in the Yunnan area, the characteristics of deformation and strain in the northern part of Yunnan Province have been studied. The results indicate that the central part of the studied area is rather stable with little crustal displacement, while the western and eastern parts are active with larger displacement. The strain field reveals that the orientations of the principal compressive strain axis of the crust and the sub-blocks in the area are NW-SE, while the orientations of the principal tensile strain axis is NE-SW. In the studied area, the tensile strain is predominatly in the northern part and the compressive strain is predominatly in the central and southern parts. The stretching direction of the shear-strain contour is basically consistent with the strike of the active fault. The strain and stress fields of the fault activity are related to the structure where the fault is located, while the activity properties of the faults are different.     

地形效应和局部地质构造对计算同震形变的影响——以2011年日本东北大地震（Mw9.0）为例

本文利用二维有限元数值模拟方法,以2011年日本东北M_w9.0级大地震为例,探讨了地形和局部地质构造等物理因素对计算同震形变的影响.计算中使用了双节点技术构成断层滑动面,同时考虑了均质、层状、地形以及俯冲板片等不同构造介质模型.对不同介质模型分别计算并进行比较,以便考察不同物理因素的影响.数值计算结果表明：海沟地形对计算同震形变的影响非常大,局部地质构造的影响也很明显,而同时考虑地形和局部地质构造所产生的影响要比仅考虑地形效应更大,其影响在水平方向和垂直方向上分别达到-1.78~0.8 m和-1.4~0.64 m,相应的百分比分别为34%和92%.这些结果表明起伏较大的地形（特别是海山）及局部地质构造对计算同震形变的影响都不容忽略,它们均能被现代大地测量技术（如GPS、InSAR等）观测到,在计算同震形变或进行断层反演时应该加以考虑.     

拟准检定法用于划分不同运动趋势块体

利用位移观测数据区分不同形变趋势块体一般分三步进行:基准检测、块体范围区分和各块体内异常形变检测．本文提出用拟准检定法来实现这三步．前两步初选拟准观测是关键．对于基准检测，按形变位移量大小；对于块体范围区分，则按形变方向．最后，通过一个模拟算例，简单介绍了每一步的实施步骤，并比较了拟准检定法与抗差估计(Huber)法、假设检验法的检测效果．结果表明，拟准检定法可以将上述三步成功实现，而抗差估计法只能检测出块体内的异常形变，对稳定基准检测和块体范围划分是失败的；假设检验法在该三步检测中均为失败．从而验证了拟准检定法具有粗差定位准确、崩溃污染率高的优越性．      

中国大陆构造块体的现今活动和变形

在重新对 1998年和 2 0 0 0年的中国地壳运动观测网络基准站和基本站的 2期观测资料进行预处理的基础上 ,得到了ITRF97坐标框架下 ,参考时刻分别为 1998年 9月 5日和 2 0 0 0年 6月 8日 ,分布在全国各主要构造块体上的 79个GPS站的坐标和协方差矩阵。分别以中国岩石圈动力学地图集 (马杏垣 ,1989)中的中国大陆主要构造单元 (称之为亚板块 )和张培震等 ( 2 0 0 2 )给出的中国主要活动块体为格架 ,用笔者提出的 1种推广了的QUAD方法对中国大陆的 2 0个主要构造块体逐个进行判别检验。那些现今无明显相对运动的相邻块体则被归并起来 ,从而确定了活动块体和它们的边界。采用刚体运动 +块体均匀应变 +局部变形的模型作为描述中国大陆构造块体的现今活动和变形的模型。求出了有明显相对运动块体的欧拉运动矢量和块体的整体均匀变形参数、各块体内部的不均匀局部变形以及活动边界的活动方式和强度。在此基础上 ,除了一般地指出中国大陆地壳运动西强东弱的特征之外 ,还对西部主要活动块体和边界活动强弱给出了定量比较结果 ,从而为强震危险区的判别提供了形变背景依据     

巴颜喀拉块体北东地区现今水平运动与变形

本文利用GPS数据研究了巴颜喀拉块体北东地区现今水平运动与变形特征。 在球坐标系中解算了各应变分量, 分析了应变率场的空间分布特征, 并与地球物理学和地震地质学研究结果进行了综合对比分析。 最新的GPS速度场结果表明, 巴颜喀拉块体北东地区与高原整体运动性质一样具有顺时针向南东方向旋转的特征, 自西向东和北东方向测站水平运动速度呈现明显的衰减特征。 应变场结果显示, 研究区以北东向的主压应变为主, 伴随着近北西向的张性应变。 应变较强的区域主要分布在活动块体的边界断裂东昆仑断裂带的东段塔藏段和龙门山断裂带上。 东昆仑断裂带东段塔藏段的主压应变明显, 结合地震地质和活动构造资料, 认为东昆仑断裂带东段塔藏段的运动性质自西向东发生了改变, 水平滑动速率逐渐减小, 垂向运动逐渐增强。 研究区GPS速度场和应变场的这一变形特征表明, 青藏高原内部的块体运动特征较为明显, 变形主要集中在作为活动块体边界的活动断裂带上, 边界断裂带的运动特征在调节活动块体间的相互运动中起着重要作用。     

Contemporary horizontal movements and seismicity of the south Baikal Basin (Baikal rift system)

The contemporary horizontal movements and deformations in the central and southern parts of the Baikal depression are analyzed, and their relationship with contemporary seismicity is studied. Based on the long-term measurements by the Baikal geodynamical GPS monitoring network, the refined estimate is obtained for the velocity of the divergence of the Siberian and Transbaikalian blocks, which is found to occur in the southeastward direction (130°) at 3.4 ± 0.7 mm per annum. This agrees with the parameters of the long-term extension component estimated from the geological data and with the direction of extension determined from the seismic data. The distribution of the displacement velocity across the strike of the rift, which gradually increases from one block to another, suggests a nonrigid behavior of the continental lithospheric plates at the divergent boundary. About 30% (1.0–1.5 mm per annum) of the total increase in the velocity is accommodated by the Baikal Basin. The strain rate within the trough reaches 3.1 × 10^?8 yr^?1 and decreases on either side across the structure. The character of distribution of the horizontal displacement velocities on the Baikal divergent boundary between the Eurasian and Amurian plates favors the model of passive rifting. The zones of highly contrasting topography and increased seismicity are localized within the area of contemporary deformations, and the seismic moment release rate directly depends on the strain rate. Here, the rate of the seismic moment release rate makes up a few percent of the geodetic moment accumulation rate calculated by the approach suggested by Anderson (1979). Based on the coherence between the graphs of the rates of geodetic moment accumulation and seismic moment release rate by the earthquakes with M ≥ 5.0 during the historical and instrumental observation periods, the contemporary seismic hazard for the South Baikal Basin could be assessed at a level of seismic event with M = 7.5–7.6.     

Dynamic mechanism for horizontal deformation in part of North China area

ｎｔｒｏｄｕｃｔｉｏｎＴｈｅｈｉｇｈｐｒｅｃｉｓｉｏｎＧＰＳｍｅａｓｕｒｅｍｅｎｔｈａｓｐｒｏｖｉｄｅｄａｎｅｆｅｃｔｉｖｅｔｅｃｈｎｉｃａｌｍｅｔｈｏｄｆｏｒｔｈｅｒｅｓｅａｒｃｈｏｆｃｒｕｓｔａｌｈｏｒｉｚｏｎｔａｌｍｏｖｅｍｅｎｔ．Ｉｎ１９...     

Surface Deformation and Seismic Rebound: Implications and Applications

An earthquake process includes pre-seismic stress accumulation, co-seismic rock rupture and post-seismic elastic and/or viscoelastic rebound. Although co-seismic and post-seismic deformations have been readily observed using the global positioning system (GPS), detecting pre-seismic stress accumulation hidden in time-series data remains challenging. This study applies the Hilbert–Huang transform to extract non-linear and non-stationary pre-earthquake deformation data from GPS records for central Taiwan. By converting the derived surface deformation into horizontal azimuths, the randomly oriented GPS-azimuths are reoriented in a similar direction several days before and after earthquakes due to loading and rebound stress, respectively. Analytical results demonstrate that the stress accumulation and release along the entire course of an earthquake process provide significant evidence supporting the seismic rebound theory. This finding would be applicable to areas with dense GPS networks and active plate interactions. Surface deformations detected by the proposed analytical technique have encouraging potential for mitigating future seismic hazards.     

A new approach for distinguishing different deformation trend blocks with displacement observations

Introduction The accumulation of high precision GPS data has great significance for understanding and explaining crustal movement and establishing the geodynamical model correctly. However, it is an important basic task to study block movement with the displacement data which covers the datum influence information, as well as the deformation information of the block and abnormal deforma-tion information in the block such as fault or zone of fracture and so on. Firstly steady or com-parativel…     

青藏块体东北缘近期水平运动与变形

利用青藏块体东北缘地区13、1年GPS观测资料，给出了本区地壳水平运动速度场及视应变场分布图，提出了由位移观测值直接求解块体旋转和变形参数的方法，初步研究了本区构造块体运动与变形特征.结果表明：①本区存在整体性向东-东南方的运动(速率约mm/a)；②南部的甘肃-青海块体的运动较快，而北部的阿拉善块体的运动较慢，二者运动速率相差近6mm/a，祁连-海原断裂带左旋走滑运动显著.③自西向东存在北北东-北东东向压性运动；④阿拉善块体、甘肃-青海块体内部存在北西西向张性变形，阿拉善块体的整体张性变形更显著，鄂尔多斯块体西侧的块体交接地带为压性运动.     

Experiment of Dividing the Active Blocks Based on Horizontal Crustal Movement Observed by GPS

In this paper, the discrepancy between the movements of intraplate blocks and plates isdiscussed, and the method to divide the intraplate active blocks is presented by selecting Bursaformula as the kinematic model for the intraplate blocks. Based on the data of three GPScampaigns in North-China network in l995, 1996 and l999, respectively, the northern area inNorth China is divided into eight small blocks with the mathematical model and methodpresented in this paper. The divided blocks based on tbe horizontal and vertical crustalmovements in the paper are consistent or approximate with each other in the area as a whole.The divided blocks in the paper is also basically accordant with the neotectonic movement,which indicates that the current movement of active blocks in this area is the succession anddevelopment of neotectonic movement. Moreover, some new activity characteristics in the areahave been revealed by the tectonic units divided with the horizontal crustal movement.     

The influence of geologic structures on deformation due to ground water withdrawal

A 62 day controlled aquifer test was conducted in thick alluvial deposits at Mesquite, Nevada, for the purpose of monitoring horizontal and vertical surface deformations using a high-precision global positioning system (GPS) network. Initial analysis of the data indicated an anisotropic aquifer system on the basis of the observed radial and tangential deformations. However, new InSAR data seem to indicate that the site may be bounded by an oblique normal fault as the subsidence bowl is both truncated to the northwest and offset from the pumping well to the south. A finite-element numerical model was developed using ABAQUS to evaluate the potential location and hydromechanical properties of the fault based on the observed horizontal deformations. Simulation results indicate that for the magnitude and direction of motion at the pumping well and at other GPS stations, which is toward the southeast (away from the inferred fault), the fault zone (5 m wide) must possess a very high permeability and storage coefficient and cross the study area in a northeast-southwest direction. Simulated horizontal and vertical displacements that include the fault zone closely match observed displacements and indicate the likelihood of the presence of the inferred fault. This analysis shows how monitoring horizontal displacements can provide valuable information about faults, and boundary conditions in general, in evaluating aquifer systems during an aquifer test.