GPS观测得到的中国大陆地壳垂直运动

利用中国地壳运动观测网络基准站的GPS连续观测及基本站非连续GPS观测结果，分析了基准站垂直位移的年周期变化特点，讨论了获得中国大陆垂直位移长趋势速率的条件，说明由基本站（连同基准站）多期GPS观测得到的长趋势垂直位移速率是较为可靠的。由于不少GPS观测站有幅度达数cm的年周期变化及大部分区域站观测次数少，由1999年、2001年两期区域站GPS观测难以得到可靠的长趋势垂直构造运动速率。基本站的观测结果表明中国大陆长趋势垂直构造运动主要特点是，速率较低，北升，南降，东强，西弱，西部相对东部略有下降。     

Determination of vertical displacements over the coastal area of Korea due to the ocean tide loading using GPS observations

This paper describes the GPS applicability for detecting the vertical displacements of ground stations caused by ocean tide loading effects. An experiment was carried out using 12 permanent GPS stations located in the coastal area of Korea using data in the period 1 July until 26 August 2003. The relative height differences were calculated from hourly DGPS data processing based on the carrier-phase observation. The power spectra of the M₂ and N₂ constituents of ocean tide loading were derived using the CLEAN algorithm. The differential vertical displacements generated by the ocean tide loading effect are typically 3–25 mm in coastal area of the Korea. We compared the results from GPS with those of the ocean tide models, NAO.99Jb regional model and GOT00.2, FES99 global models. The M₂ (N₂) amplitude differences of vertical displacements between GPS and GOT00.2 is 1.22 ± 3.61 mm (1.01 ± 1.48 mm), and that of the M₂ (N₂) amplitude difference between GPS and FES99 is 0.04 ± 4.64 mm (0.64 ± 1.75 mm), whereas the M₂ (N₂) amplitude difference between GPS and NAO.99Jb are 0.05 ± 1.03 mm (0.86 ± 1.18 mm). The highest vertical displacements at the PALM station are found for 24.5 ± 0.7 mm from GPS observation, and 22.9 mm from the regional model NAO.99Jb and 13.17 and 10.00 mm from the global models GOT00.2 and FES99, respectively. These values show that the vertical displacements derived from GPS are in good agreement with those of the regional model NAO.99Jb around Korea, more than with the global models. This result indicated that GPS is an effective tool to measure the vertical displacement caused by the ocean tide loading effect in the coastal area, and we need to use the NAO.99Jb ocean tide model rather than the global ocean tide models in and around the Korean peninsula for position determination with permanent GPS installations. This work demonstrates that vertical displacement caused by the M₂ and N₂ constituents of ocean tide loading can be measured by carrier-phase DGPS.     

GPS测定的四川九寨沟MW6.6地震同震形变

2017年8月8日四川阿坝州九寨沟发生M_W6.6地震,震源机制解显示该地震为左旋走滑型地震。对震中周围的GPS连续站观测资料进行处理,获得高频GPS动态形变和静态同震水平位移。震中100km范围内四川松潘和甘肃武都站观测到1 Hz动态形变。距离震中约69km的松潘站观测的同震水平位移为7.4mm。根据少量的GPS静态同震位移反演的同震破裂模型显示本次地震的最大滑动量为376mm,地震矩为7.25×1018 N·m,等效矩震级为M_W6.6。正演计算的同震三维形变场显示本次地震的最大水平位移可达4~5cm,垂直位移呈四象限分布,最大可达1.5cm,区域内10个流动GPS站可观测到同震形变。     

长白山天池火山近期形变场演化过程分析

对长白山天池火山区2002～2005年的4期水准测量数据、2000～2005年6期GPS测量数据以及1999年以来的定点地倾斜观测数据进行了系统处理分析,结果表明：（1）2000年9月～2002年6月期间,长白山天池火山地表变形很小,基本上处于较为稳定的状态. 1993～1998年的InSAR观测资料也未发现明显的地表变形;（2）2002～2003年,长白山天池火山出现了明显的地表形变异常变化,其水平形变表现为以天池火山口为中心的放射状膨胀运动,垂直形变则表现为近火山口的快速隆起,显示天池火山下部可能有岩浆压力增大变化;（3）长白山天池火山2003年以来的形变速率有逐年迅速衰减的变化特征,其最大水平运动速率已由2002～2003年期间的38.3mm/a衰减到2004～2005年的15.8mm/a,最大垂直形变速率则由2002～2003年的46.0mm/a 衰减到2004～2005年的5.0mm/a.     

高频GPS双差残差模型监测强震地表运动

根据震时地表震动持续时间短的特点,以及双差残差的相关性,提出了基于双差残差建模的高频GPS强震位移监测方法.首先对震前若干历元各双差残差进行建模,地震发生后,利用模型对各双差残差进行预报,最后利用预报残差实现短时间尺度内的位移解算.利用一条长约1100 km基线的静态数据和El Mayor-Cucapah 7.2级地震的94个测站数据进行试验分析,结果表明,5 min预测时间内,静态数据的动态位移在N、E、U三个方向的定位中误差分别为6 mm、6 mm和13 mm,地震数据的解算结果与实际情况有较好的一致性.     

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

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

Crustal movement in Chinese mainland observed from 1998 to 1999

IntroductionTheCrustalMovementObservationNetworkofChina(CMONOC)isoneofmajorscientificprojectsinChinaduringthe9thFive-YearPlan.ThefirstGPSobservationsof25fiducialstationsand56basicstationsinthenetworkwerecompletedfromAugust26toSeptember6,1998andallthe81stationswereobservedsimultaneously.25fiducialstationshavebeenincontinuousoperationsincelateMarch1999.Inaddition,dataofcontinuousobservationfromtheIGS(InternationalGPSService)stationXIANinXi(anarealsoavailable.Figure1showsthelocation…     

亚太网三期GPS观测数据处理结果及地壳形变特征分析

本文对亚洲太平洋地区地球动力学计划(APRGP,简称亚太网)在1999-2001年采集的三期GPS观测数据,利用GAMIT/GLOBK软件进行了分析处理,得出各期的观测站在ITRF97框架下的三维直角坐标,并根据三期观测数据中的重复观测站进行了速度场计算和地壳形变特征分析。结果表明亚洲大陆板块具有明显的顺时针旋转运动特征;印度板块和澳大利亚板块从西南方向对亚洲大陆板块进行挤压;而西太平洋地区各观测站则向西北方向运动。这些特征与国内GPS区域网观测数据处理得出的结论是一致的,而本文则揭示了更大范围的亚洲太平洋地区地壳形变信息。     

Estimation of the time-dependent crustal movements of the zmit Earthquake

The 17 August 1999 zmit earthquake significantly deformed the earth's crust in the Marmara Region, especially in the Gölcük–Sapanca Zone, Turkey. It broke a 150-km long segment of the northern branch of North Anatolian Fault Zone. The geodetically determined moment magnitude was M_w=7.5. Global Positioning System (GPS) sites, which are a small subset of the Marmara Continuous GPS Network (MAGNET), and survey sites in the region were studied to estimate coseismic and postseismic deformations, using different methodologies with linear, quadratic and exponential kinematic models. Six GPS epochs for these sites, which were carried out before and after the 17 August 1999 zmit earthquake, were used to define the kinematic models. The quadratic deformation model was also applied to determine the time-dependent crustal movement parameters (velocity and acceleration) of the sites, using the Kalman filter technique. In order to show the differences between the models, the estimated deformation fields on the last epoch were compared. In all models, as expected, the faults near the sites show large coseismic displacements with fault parallel direction, whereas the far sites show small coseismic displacements due to the effects of the zmit earthquake. Each kinematic model, fitted to the epochs after the earthquake, shows different behaviour. While the linear model shows insufficient results, the nonlinear models (quadratic and exponential) give the best fitted to the postseismic deformations. As a result of Kalman filter analysis, the fault near-sites shows significant velocities with fault parallel direction, whereas the far sites have insignificant velocities. All stations have insignificant accelerations in the last epoch.     

Pleistocene and recent deglaciation in Svalbard: implications for tide-gauge,GPS and VLBI measurements

Apart from sea-level change caused by ocean-water volume change, tide-gauge measurements are also affected by vertical movement of the earth's surface and by geoid change. To study these phenomena more closely, we consider the tide-gauge station at Ny-Ålesund, Svalbard, where vertical movements are controlled by GPS and VLBI measurements. Whereas the tide-gauge record indicates a relative sea-level fall of about 2.6 mm a⁻¹, the GPS and VLBI measurements suggest a land uplift of about 5.6 mm a⁻¹. We predict the deglaciation-induced vertical movement and geoid change by combining the Pleistocene ice model ICE-3G with the newly developed Recent ice model SVAL. We find that, for particular deglaciation histories and earth models with an asthenosphere, the predicted land uplift matches the measured land uplift rather closely. However, for these combinations of ice and earth models, the predicted sea-level fall is consistently too large by at least 3 mm a⁻¹. Considering the uncertainties and simplifications involved in the study, the discrepancy weakly suggests a sea-level rise due to increased ocean-water volume.     

Loading effects in Metsähovi from the atmosphere and the Baltic Sea

Loading by atmosphere and by the Baltic Sea cause gravity change at Metsähovi, located 15 km from the open sea. Gravity is changed by both the Newtonian attraction of the loading mass and by the crustal deformation. We have performed loading calculations using appropriate Green's function for both gravity and deformation, for both atmospheric and Baltic loading. The loading by atmosphere has been computed using a detailed surface pressure field from high resolution limited area model (HIRLAM) for north Europe up to 10° distances. Baltic Sea level is modelled using tide gauge records. Calculations show that 1 m of uniform layer of water corresponds to 31 nm s⁻² in gravity and −11 mm in height. Modelled loading is compared with observations of the superconducting gravimeter T020 for years 1994–2002. The combination of HIRLAM and a tide gauge record decreases RMS of gravity residuals by 14% compared to single admittance in air pressure corrections without sea level data. Regression of gravity residuals on the tide gauge record at Helsinki (at 30 km distance) gives a gravity effect of 26 nm s⁻² m⁻¹ for Baltic loading.The gravity station is co-located with a permanent GPS station. We have also associated the loading effects of the atmosphere and of the Baltic Sea with temporal height variations. The range of modelled vertical motion due to air pressure was 46 mm and that due to sea level 18 mm. The total range was 38 mm. The effects of the Baltic Sea and of the atmosphere partly cancel each other, since at longer periods the inverse barometer assumption is valid. Regression of the modelled height on local air pressure gives −0.37 mm hPa⁻¹, corresponding approximately to width 6° for pressure system.We have tested the models using one year of daily GPS data. Multilinear regression on local air pressure and sea level in Helsinki gives the coefficient −0.34 mm hPa⁻¹ for pressure, and −11 mm m⁻¹ for sea level. These match model values. Loading by air pressure and Baltic Sea explains nearly 40% of the variance of daily GPS height solutions.     

基于GPS观测的北天山主要断裂现今构造运动特征研究

借助分布在北天山地区最新GPS点位的运动观测资料, 利用GAMIT/GLOBK数据处理软件获取了北天山地区现今地壳的运动位移场． 以该位移场为基础, 利用弹性半空间位错理论, 估算了研究区内博罗科努—阿其克库杜克断裂和准噶尔盆地南缘断裂两条具有代表性的主要断裂的现今活动速率． 结果表明: 博罗科努—阿其克库杜克右旋走滑断裂东、 西两段滑移速率的差异性不明显, 1944年3月10日乌苏南M_S7.2强震发生后, 该断层现今表现为震后微蠕滑运动, 东、 西两段滑动速率均在1—2 mm/a之间; 准噶尔盆地南缘断裂现今滑动速率为(5.6±1.0) mm/a．      

A Reanalysis and Reinterpretation of Geodetic and Geological Evidence of Glacial-Isostatic Adjustment in the Churchill Region, Hudson Bay

We review the historical, geological, tide-gauge, GPS and gravimetric evidence advanced in favour of, or against, continuing land uplift around Hudson Bay, Canada. We also reanalyse the tide-gauge and GPS data for Churchill using longer time series than those available to previous investigators. The dependence of the mean rate of relative sea-level change obtained from the tide-gauge record on the length and mid-epoch of the observation interval considered is investigated by means of a newly developed linear-trend analysis diagram. For studying the shorter-period variability of the tide-gauge record, the wavelet transform is used. The mean rate of land uplift obtained from GPS is based on a new analysis using IGS solutions of GFZ. To include the post-glacial land uplift, sea-level indicators from the Churchill region representing the relative sea-level history during the past 8000 years are also used. Finally, the values of the four observables are jointly inverted in terms of mantle viscosity. The optimum values are ~3.2 × 10²⁰ Pa s and ~1.6 × 10²² Pa s for the upper- and lower-mantle viscosities, respectively.     

Velocity field and tectonic strain in Southern Spain and surrounding areas derived from GPS episodic measurements

The goal of this paper is to study the velocity field and deformation parameters in Southern Spain and surrounding areas (Ibero-Maghrebian region) using GPS episodic measurements. Results are compared to those previously published as well as deformation parameters derived from seismic data. For this purpose, a geodetic GPS network of 12 stations was observed during eight field campaigns from 1998 to 2005 by the San Fernando Naval Observatory (ROA), Spain. Relative GPS velocities in the Gulf of Cadiz with respect to the stable part of Eurasia are ～4.1 mm/yr in a NW–SE to NNW–SSE direction. In the Betics, Alboran Sea and North of Morocco, velocities are ～4.4 mm/yr in a NW–SE direction, and they are ～2.3 mm/yr in a N–S direction in the eastern part of the Iberian Peninsula. These results are in agreement with the anticlockwise rotation of the African plate. GPS strain tensors are determined from the velocity model, to obtain a more realistic crustal deformation model. The Gulf of Cadiz is subjected to uniform horizontal compression in a NNW–SSE direction, with a rotation to N–S in the Alboran Sea and Northern Morocco. An extensional regime in a NW–SE direction, which rotates to W–E, is present in the Internal Betics area. In the Betic, Alboran Sea and North of Morocco regions we compare seismic deformation rates from shallow earthquakes with the determined GPS deformation rates. The comparison indicates a seismic coupling of 27%, while the remaining 73% might be generated in aseismic processes. Deformations measured in the Ibero-Maghrebian region with GPS could be interpreted in terms of either elastic loading or ductile deformation.     

Dynamic mechanisms of the post-seismic deformation following large events: Case study of the 1999 Chi-Chi earthquake in Taiwan of China

The mechanism of postseismic deformation related to strong earthquakes is important in geodynamics, and presumably afterslip or viscoelastic relaxation is responsible for the postsesimic deformation. The 1999 Chi-Chi, Taiwan of China, earthquake occurred in the region where GPS observation station is most densely deployed in the world. The unprecedented GPS data provides a unique opportunity to study the physical processes of postseismic deformation. Here we assume that the interactions of viscoelastic relaxation, afterslip, fault zone collapse, poroelastic rebound, flow of underground fluids, and all these combined contribute to the surface displacements following the main shock. In order to know the essence of the postseismic deformation after the strong event, fault zone collapse, poroelastic rebound, flow of underground fluids, and so on, are represented equivalently by the variations of the focal medium properties. Therefore, the viscoelastic relaxation, afterslip, and the variations of the equivalent focal medium properties are inverted by applying the GPS temporal series measurement data with viscoelastic finite element method. Both the afterslip rate distribution along the fault and the afterslip evolution with time are obtained by means of inversion. Also, the preliminary result suggests that viscosities of the lower crust and the upper mantle in Taiwan region is 2.7×10¹⁸ and 4.2×10²⁰ Pa·s, respectively. Moreover, the inversion results indicate that the afterslip contributing to postseismic deformation of 44.6% in 450 days after the Chi-Chi earthquake, with 34.7% caused by the viscous relaxation and 20.7% by other factors such as fault zone collapse, poroelastic rebound, and the flow of liquids.     

Error analysis for European IGS stations

Each of the GPS time series that describes the changes of topocentric components consists of a deterministic and a stochastic part, whose character influences the errors of the deterministic parameters. As to the uncertainties of reliable velocities of permanent satellite station systems, surveys that estimate and take into account any dependencies that may affect subsequent operational efficiency are very important. For this analysis, we used 42 stations from the IGS (International GNSS Service) network from Europe, processed at the Military University of Technology EUREF Permanent Network Local Analysis Centre (MUT LAC). The deterministic part of the GPS time series was removed using the least squares method. The seasonal periods in topocentric components were determined assuming the existence of the residual Chandler oscillation (1.67 cpy), as well as the annual tropical (1 cpy) and draconitic (1.04 cpy) oscillations with their harmonics up to 4th. We assumed the character of the residue as a combination of white and powerlaw noise. The obtained results show, that in the case of the European sub-network of IGS stations we are dealing with the coloured noise between white and flicker noise with the amplitudes between 3 to 6 mm/year^-k/4 for horizontal components and between 6 to 15 mm/year^-κ/4 for the vertical ones, where κ is a spectral index. Finally, we showed that the amplitudes and spectral indices of noise are reduced after performing a spatio-temporal filtering. All the elicited results referred to the uncertainties of velocities by estimating them before and after filtration and the simulation of their values for different lengths of the time series.     

Test of GPS for permanent ionospheric TEC monitoring at high latitudes

The Global Positioning System (GPS) observables are affected by the ionosphere. The dispersive nature of this effect and the use of two frequencies in the GPS observations make possible to measure the ionospheric total electron content (TEC) from dual frequency GPS data. In this work we test the concept of permanent monitoring of TEC using a network of GPS receivers at high latitudes. We have used GPS data from five permanent receivers in Scandinavia, from 1–30 January 1994, with geographic latitudes ranging from 57.4^°N to 78.9^°N. The results show the capability of the method to monitor the evolution of TEC as a function of time and geographical location. We have detected night-time enhancements almost every night for some of the stations, and we have also been able to produce maps of the instantaneous TEC as a function of both latitude and longitude around the GPS network. We also present some of the current limitations in the use of GPS for estimating TEC at high latitudes such as the difficulties in solving for cycle-slips, and the necessity of reliable values for the receiver and satellite differential instrumental biases.     

呼图壁地下储气库地表盖层变形的GPS研究

利用新疆呼图壁地下储气库地表盖层由13个点位组成的形变监测网的前5期GPS观测资料,研究地下储气库注采过程中地表盖层的变形响应。通过获取地下储气库运行过程中地表盖层形变的三维时间序列,并结合井口压力数据,区分地下储气库在不同过程中的变形信号。研究结果表明,地表盖层在储气库注采过程中水平方向上存在明显的“呼吸效应”,储气库每MPa气井压力变化在注、采周期内对地表变形造成的影响在水平方向上分别达到1.02、1.24mm,垂直方向分别达到-1.11、0.86mm。     

Infiltration-excess caused by Stemflow in a cyclone-prone tropical rainforest

Tropical rainforest canopy trees that have large projected areas of upwardly inclined branches are capable of funnelling large volumes of rainwater down their trunks. During periods of prolonged heavy rainfall on Mount Bellenden Ker in northeast Queensland, Australia, stemflow volumes were found to be as much as two orders of magnitude greater than the volume of incident rainfall expected in a rain gauge occupying an area equal to the trunk basal area. Stemflow totals ranging from 6000 to 70000 litres were generated by individual trees from 7800 mm of rainfall over two successive wet seasons. The combination of high intensity rainfall and the funnelling effect results in significant quantities of infiltration-excess at the ground surface. Stemflow fluxes as high as 31.4 cm³ min^?1 per cm² of basal area (i.e. the equivalent of 314 mm min^?1) were recorded when rainfall intensity was only 2 mm min^?1. The mean infiltration capacity of the topsoil was determined to be 6.2 mm min^?1. The areas over which the stemflow would have had to spread in order to infiltrate were computed to be as much as 3 m² around the bases of individual canopy trees. Approximations of the distances that the infiltration-excess would have travelled away from the tree bases were calculated by assuming that the infiltration area either expands radially outward in the form of an annulus or extends straight downslope from the tree base.     

Systematic errors between VLBI and GPS precipitable water vapor estimations from 5-year co-located measurements

Space geodetic techniques (e.g., Global Positioning System, GPS and very long baseline interference, VLBI) have been widely used to determine the precipitable water vapor (PWV) for meteorology and climatology, which was verified by comparing with co-located independent technique observations. However, most of these comparisons have been conducted using only short-time spanning observations at several stations. The goal of this study is to identify and quantify the systematic errors between VLBI and GPS PWV estimates using a 5-year (2002–2007), PWV data set constructed from co-located measurements and radiosonde data as well. It has found systematic errors between VLBI and GPS PWV estimations from comparisons with long-term co-located GPS measurements. The total mean VLBI PWVs are systematically smaller than GPS estimates with 0.8–2.2 mm for all sites, but can be as much as 15–30%. The subdiurnal PWV variation magnitudes and long-term trends between VLBI and GPS are nearly similar, but the VLBI-derived PWV trends are systematically smaller than GPS estimates with about 0.1±0.02 mm/year. These systematic errors in PWV estimates between VLBI and GPS are probably due to technique own problems, different used elevation angles and co-location separation.