Testing of the IERS2000 Sub-Daily Earth Rotation Parameter Model

The differences between the new International Earth Rotation Service (IERS) 2000 and the previous IERS1996 sub-daily Earth rotation parameters (ERP) models can reach 0.1 mas (0.001 arc sec) and 0.1 mas/day. The largest differences are seen for the aliasing periods of 14.2 and 360 days, which correspond to the diurnal tidal waves of O1 and (K1, P1), respectively. Precise independent polar motion (PM) rate solutions effectively doubles the sampling rate and allows for effective testing of sub-daily ERP models and other periodical effects at the diurnal and semi-diurnal frequency bands. Since November 12, 2000, when the Jet Propulsion Laboratory (JPL) Analysis Center of International GPS Service (IGS) has switched to the conventional IERS1996 sub-daily ERP model, from the older model of Herring and Dog (1994), the JPL daily PM rate solutions show no, or greatly reduced 14.2 day amplitude (O1) peaks. This confirmed that the anomalistic amplitudes at 14.2 day period seen for JPL PM solutions prior November 12, 2000 was largely due to the effects of the older sub-daily ERP model on independent PM rate solutions. As indicated by the latest IGS PM rate solutions, which were corrected for the IERS1996 and 2000 model differences, the new IERS2000 sub-daily ERP model is expected to perform equally well as the conventional IERS1996 model.     

A Possible Detection of the 26 December 2004 Great Sumatra-Andaman Islands Earthquake with Solution Products of the International GNSS Service

The main goal of this work is to critically review the IGS solution products and Precise Point Positioning (PPP) in order to demonstrate their potential to contribute to studies of large earthquakes such as the one that devastated Southeast Asia on December 26th, 2004. In view of a possible detection of the Mw 9.0 Sumatra-Andaman Islands Earthquake of December 26, 2004, position solutions, ranging from intervals of years to one second, of four International GNSS Service (IGS) stations within 3000 km of the epicenter were examined. The IGS combined, cumulative solution product (IGS04P51), consisting of epoch and station velocity solutions and based on data spans of several years prior to the earthquake, was used as a reference. Four IGS combined weekly position solutions (igs04P1301-4), two weeks before and after the earthquake, were utilized for the weekly solution resolution. PPP static and kinematic solutions with IGS Final combined orbits and clocks were used for the mean daily and instantaneous 5-min and 1-sec epoch solutions, respectively. The most significant changes, detected by both weekly and daily solutions occurred in longitude. The nearest IGS station ntus, about 1000 km east of the epicenter, moved westward about 15 mm, while the more distant Indian station iisc (∼ 2300 km NW from the epicenter), shifted about 15 mm eastward. In spite of position errors caused by interpolation of the 5-min IGS clocks, the 1-sec solutions, based on separate data sets, available only for two stations (iisc, dgar), still showed seismic surface waves, in particular at the Indian station iisc. Precise daily IGS combined polar motion and length-of-day products, after correcting for the atmospheric effects, also likely detected, statistically significant, anomalistic excitations on December 26, 2004 that could be caused by this great earthquake.     

地球参考框架的发展现状和未来展望

地球参考框架是国家重要的空间基础设施,是地球上人类所有活动的空间参考基准。本文首先阐述了国际地球参考框架（International Terrestrial Reference Frame,ITRF）的发展现状,重点评述了ITRF的建立与维持,针对ITRF的发展现状提出了存在的问题;其次,以ITRF与2000国家大地坐标系（China Geodetic Coordinate System 2000,CGCS2000）的关系及现状为切入点,探讨了我国建立北斗坐标系的必要性,介绍了建立北斗坐标系的基本思路以及初始实现;最后,对地球参考框架的未来发展进行了展望。     

Measuring Seismic Waves Induced by Large Earthquakes with GPS

Independent GPS position solutions at 1-sec interval, derived from the International GPS Service (IGS) data and orbit÷clock products, clearly show seismic waves generated by the magnitude 7.9 Denali Fault, Alaska earthquake of November 3, 2002. Surface seismic waves with periods of about 20 sec and amplitudes of up to 20 cm were detected up to 4,000 km from the epicenter. This confirms the previous findings reported by Larson et al. (2003); we use additional station data along with different processing software and strategies. The seismic waves from the May 26, 2003 magnitude 7.0 Japanese earthquake were also observed in the 1-sec position solution series at station MIZU, about 80 km from the epicenter. This earthquake, however, could not be detected by GPS at station USUD, about 410 km away. Similarly, the Algerian May 21, 2003 earthquake of magnitude 6.8 could not be detected by GPS at the nearest IGS station located approximately 800 km from the epicenter.     

利用中国地壳运动观测网络研究中国大陆相对于ITRF97板块模型形变

在建立全球ITRF97板块运动模型的基础上,利用"中国地壳运动观测网络"79个GPS基本站的数据,建立我国新的地壳运动方向图和块体运动模型.通过与NNR-NUVEL1A地质模型比较认为,ITRF97板块运动模型反映了现今十几年跨度的地壳运动,在研究我国现今几年到十几年时间跨度的地壳形变时,地壳运动背景场应采用基于ITRF97实测速度场建立欧亚板块运动模型.     

Terrestrial reference frame requirements within GGOS perspective

One of the main objectives of the promising and challenging IAG project GGOS (Global Geodetic Observing System) is the availability of a global and accurate Terrestrial Reference Frame for Earth Science applications, particularly Earth Rotation, Gravity Field and geophysics. With the experience gained within the activities related to the International Terrestrial Reference System (ITRS) and its realization, the International Terrestrial Reference Frame (ITRF), the combination method proved its efficiency to establish a global frame benefiting from the strengths of the various space geodetic techniques and, in the same time, underlining their biases and weaknesses. In this paper we focus on the limitation factors inherent to each individual technique and to the combination, such as the current status of the observing networks, distribution of the co-location sites and their quality and accuracy of the combined frame parameters. Results of some TRF and EOP simultaneous combinations using CATREF software will be used to illustrate the current achievement and to help drawing up future goals and improvements in the GGOS framework. Beyond these technical aspects, the overall visibility and acceptance of ITRS/ITRF as international standard for science and applications is also discussed.     

Terrestrial reference frame requirements within GGOS perspective

One of the main objectives of the promising and challenging IAG project GGOS (Global Geodetic Observing System) is the availability of a global and accurate Terrestrial Reference Frame for Earth Science applications, particularly Earth Rotation, Gravity Field and geophysics. With the experience gained within the activities related to the International Terrestrial Reference System (ITRS) and its realization, the International Terrestrial Reference Frame (ITRF), the combination method proved its efficiency to establish a global frame benefiting from the strengths of the various space geodetic techniques and, in the same time, underlining their biases and weaknesses. In this paper we focus on the limitation factors inherent to each individual technique and to the combination, such as the current status of the observing networks, distribution of the co-location sites and their quality and accuracy of the combined frame parameters. Results of some TRF and EOP simultaneous combinations using CATREF software will be used to illustrate the current achievement and to help drawing up future goals and improvements in the GGOS framework. Beyond these technical aspects, the overall visibility and acceptance of ITRS/ITRF as international standard for science and applications is also discussed.     

Earth Rotation Parameters 1899.7–:1992.0 After Reanalysis Within The Hipparcos Frame

Earth rotation parameters (ERP) in the interval 1899.7–1992.0 are obtained from re-analysis of the observed latitude/universal time variations by optical astrometry. Hipparcos Catalogue is used to define the celestial reference frame, within which the ERP are described, with special care devoted to 'problematic' double and/or multiple stars. The terrestrial reference frame is defined by the adopted latitudes/longitudes of participating instruments and their secular motions as given by the NUVEL-1 NNR model of plate motions, and it is chosen to be very close to the International Terrestrial Reference Frame (ITRF). More than four million observations made with 48 different instruments at 31 observatories, located all over the world, are utilized to determine polar motion, celestial pole offsets and (after 1956) universal time UT1, all at 5-day intervals. Along with these parameters, the combinations of Love and Shida numbers, governing the tidal variations of the local verticals at individual observatories, are also determined. The analysis of the results covering almost a century, namely the long-periodic polar motion and length-of-day changes, is presented.     

The International Laser Ranging Service and its support for IGGOS

The International Laser Ranging Service (ILRS) was established in September 1998 as a service within the IAG to support programs in geodetic, geophysical, and lunar research activities and to provide data products to the International Earth Rotation Service (IERS) in support of its prime objectives. Now in operation for 5 years, the ILRS develops: (1) the standards and specifications necessary for product consistency and (2) the priorities and tracking strategies required to maximize network efficiency. The service collects, merges, analyzes, archives and distributes satellite and lunar laser ranging data to satisfy a variety of scientific, engineering, and operational needs and encourages the application of new technologies to enhance the quality, quantity, and cost effectiveness of its data products. The ILRS works with: (1) the global network to improve station performance; (2) new satellite missions in the design and building of retroreflector targets to maximize data quality and quantity and (3) science programs to optimize scientific data yield. The ILRS Central Bureau maintains a comprehensive web site as the primary vehicle for the distribution of information within the ILRS community. The site, which can be accessed at: http://ilrs.gsfc.nasa.gov is also available at mirrored sites at the Communications Research Laboratory (CRL) in Tokyo and the European Data Center (EDC) in Munich.During the last 2 years, the ILRS has addressed very important challenges: (1) data from the field stations are now submitted hourly and made available immediately through the data centers for access by the user community; (2) tracking on low satellites has been significantly improved through the sub-daily issue of predictions, drag functions, and the real-time exchange of time biases; (3) analysis products are now submitted in SINEX format for compatibility with the other space geodesy techniques; (4) the Analysis Working Group is heavily engaged in Pilot Projects as it works toward an ILRS “standard” global solution and (5) SLR has significantly increased its participation in the International Terrestrial Reference Frame (ITRF) activity, which is important to the success of IGGOS.     

基于服务系统的实时精密单点定位技术及应用研究

讨论了实时精密单点定位的若干关键技术,着重研究了实时精密钟差估计算法;利用全球40个IGS跟踪站连续、实时的观测数据流进行实时精密卫星钟差估计,解算出的实时精密钟差与CODE事后精密钟差具有较好的一致性,二者互差的RMS优于0.2 ns.采用实时估计的卫星钟差和IGS发布的超快速精密卫星轨道进行实时精密单点定位模拟实验,结果表明：经过大约15~30 min初始化后,实时精密单点定位滤波收敛之后水平方向的定位精度优于5 cm,高程方向的精度优于10 cm,能满足在该精度级别的实时用户的定位需求.处理汶川地震期间震区的GPS观测数据,结果表明实时精密单点定位能够探测地震发生期间的地面同震位移.     

基于多颗低轨卫星的SLR台站坐标几何解算方法

卫星激光测距（SLR）技术作为卫星精密定轨手段和轨道检核重要方法,激光反射器已经成为重力卫星和测高卫星等低轨卫星的基本载荷.经典的SLR台站坐标是使用动力学方法计算的,本文根据多颗低轨卫星（LEO）多历元的激光观测数据,采用几何方法开展地面SLR测站坐标计算.通过组建低轨卫星群实现对全球激光站的动态观测,为了合理配置不同低轨卫星间观测值权重,削弱低轨卫星群可能存在的系统性偏差,提出采用方差分量估计组合的最小二乘法进行解算.实测结果显示,解算出SLR台站坐标框架解与SLRF2014差异平均值在25.1 mm,外符合精度达到1~2 cm.该方法避免了复杂动力学模型,SLR台站坐标的几何计算方法既可以作为激光测站框架解算手段之一,同时将LEO卫星群作为空间并址站实现不同技术地球参考框架间的融合.     

地表负载及GPS测站分布对参考框架转换的影响分析

GPS数据处理通常采用Helmert七参数转换将瞬时站坐标转换到指定的框架下,但瞬时站坐标中尚未模型化的季节性地表负载会影响平移参数(地心运动)和尺度参数的估值,进而影响测站坐标;不均匀的测站分布会加剧这一影响.本文利用GRACE重力场系数仿真GPS地表负载的实验表明,基于网平移法采用实际的IGS站至少能够恢复90%的地心运动信号.地表负载及GPS实际测站的不均匀分布可以解释大约30%的GPS尺度的周年变化.相对于IGb08的所有框架站,目前采用91个全球均匀分布的核心站作为框架转换的基准是合理的.采用IGb08的所有框架站进行转换会导致U方向误差增加,特别是对框架站密集的欧洲区域(误差均值约为1mm).因此框架转换时,应尽量选取均匀分布的测站,同时不估计尺度参数.     

关于ITRF96参考架整体旋转性的探讨

利用国际地球自转服务IERS发布的国际地球参考架ITRF96的速度场 ,建立了一个基于现代空间大地测量实测结果的全球板块运动模型ITRF96VEL ,该模型与NNR -NUVEL1A模型有着较好的一致性 .利用该模型 ,本文对ITRF96参考架相对无整体旋转参考架是否存在一个整体旋转的问题进行了初步的探讨 ,结果表明ITRF96参考架相对无整体旋转参考架有一个整体旋转 ,旋转角速度为 0 1 61°/Ma ,旋转极指向南纬 50 5° ,东经 65 5°,这与要求ITRF96无整体旋转的定义不相符 ,这个整体旋转将会对国际地球参考架的高精度应用和长期维持 ,特别是对地球自转参数长期变化的研究产生一定的影响     

Tectonic plates parameters estimated in International Terrestrial Reference Frame ITRF2008 based on DORIS stations

This paper concerns an analysis of the accuracy of the estimated parameters Ω (Φ, Λ, ω) which define the tectonic plate motions. The study is based on the velocities of station positions in the IERS (International Earth Rotation and Reference Systems Service) which has published new realization of the International Terrestrial Reference System—ITRF2008 for Doppler Orbitography by Radiopositioning Integrated on Satellite DORIS technique. Eurasian, African, Australian, North American, Australian, Pacific, Antarctic and South American plates were used in the analysis. The influence of the number and localization of stations on the plate surface on the estimation accuracy of the tectonic plate motion parameters were discussed. The results were compared with the APKIM 2005 IGN model and our earlier estimation for the SLR technique. In general, a remarkable concurrent agreement between the present and the APKIM 2005 solutions was found.     

固定非差整数模糊度的PPP快速精密定位定轨

从GPS基本观测模型出发,给出并推导了分离相位小数偏差求解非差整数模糊度的精密单点定位数学模型和算法.利用少量IGS跟踪站组成服务端观测网计算未检校的相位小数偏差改正信息,用于改正用户端接收机的相位观测值,实现了固定非差整数模糊度的快速精密单点定位与定轨.试验结果表明: 利用30 min的地面动态或静态观测数据进行精密单点定位,非差模糊度固定成整数后,其定位结果较PPP浮点解有明显改善,水平方向提高了近一个数量级,可达到1 cm甚至毫米级的精度;高程方向与对流层延迟解算精度也改善了20%~60%.与浮点解相比,固定解能显著改善PPP的定轨精度,仅用15 min的短弧段观测数据,切向与法向的定轨精度可达到1 cm左右;径向方向为3~5 cm左右,较浮点解定轨精度改善了50%~70%.因此,固定非差整数模糊度后的PPP能够满足毫米至厘米级的快速精密定位和定轨的要求,这在GPS(准)实时应用与服务中具有很好的应用前景.     

基于GNSS网络的实时精密单点定位及精度分析

基于局域GNSS网络,以历元间、星间差分技术实时估计了GPS卫星相对钟差的历元间差值;针对所估计实时精密钟差的特征,推导了实时精密单点定位的估计模型.对所估计得到的相对钟差的历元间差值、实时定位结果分别进行了比较、分析.结果表明,相对钟差的历元间差值与IGS的最终星历相比其精度可以达到0.08 ns;每小时观测的实时静态定位结果在N、E、U三个方向的精度分别为1.47、3.62、4.09 cm.动态模式,实时结果在N、E、U三个方向的精度分别为2.63、3.82、5.20 cm.与采用IGS最终轨道和钟差解算的结果相比较,实时计算结果优于采用精密轨道和精密钟差计算结果.     

Differences in GPS coordinate time series resulting from the use of individual instead of type-mean antenna phase center calibration model

It is well-known that the phase center of a Global Navigation Satellite System (GNSS) antenna is not a stable point coinciding with a mechanical reference. The phase center position depends on the direction of the received signal, and is antenna-and signaldependent. Phase center corrections (PCC) models of GNSS antennas have been available for several years. The first method to create antenna PCC models was the relative field calibration procedure. Currently only absolute calibration models are generally recommended for use. In this study we investigate the differences between position estimates obtained using individual and type-mean absolute antenna calibrations in order to better understand how receiver antenna calibration models contribute to the Global Positioning System (GPS) positioning error budget. The station positions were estimated with two absolute calibration models: the igs08.atx model, which contains typemean calibration results, and individual antenna calibration models. Continuous GPS observations from selected Polish European Permanent Network (EPN) stations were used for these studies. The position time series were derived from the precise point positioning (PPP) technique using the NAPEOS scientific GNSS software package. The results show that the differences in the calibrations models propagate directly into the position domain, affecting daily as well sub-daily results. In daily solutions, the position offsets, resulting from the use of individual calibrations instead of type-mean igs08.atx calibrations, can reach up to 5 mm in the Up component, while in the horizontal one they generally stay below 1 mm. It was found that increasing the frequency of sub-daily coordinate solutions amplifies the effects of type-mean vs individual PCC-dependent differences, and also gives visible periodic variations in time series of GPS position differences.     

Assessment of GPS data for meteorological applications over Africa: Study of error sources and analysis of positioning accuracy

The aim of this study is to assess the availability and quality of data from the International GNSS Service (IGS) Global Positioning System (GPS) network in Africa, especially for retrieving zenith tropospheric delay (ZTD), from which precipitable water vapour (PWV) can be derived, in view of application to the African Monsoon Multidisciplinary Analysis (AMMA) project. Three major error sources for the GPS data analysis evaluating PWV in Africa are the accuracy of the satellite orbits, the correction for the radio delay induced by the ionosphere and the vertical site displacements due to ocean loading. The first part of this study examines these error sources and the validity of GPS data for meteorological applications in Africa in dedicated analyses spanning the year 2001. These analyses were performed using the IGS precise orbits. Weak degradation of baseline precision with increasing baseline lengths suggests that the average orbital error is not limiting the GPS analysis in Africa. The impact of the ionosphere has been evaluated during a maximum of solar activity in 2001. The loss of L₂ data has actually been observed. It amounts to 2% on average for 2001, with maxima of 8% during magnetic storm events. A slight decrease in formal accuracy of ZTD seems to be related to the loss of L₂ data at the end of the day. This indicates that scintillation effects are present in the GPS observations but however are not a major limitation. The impact of ocean loading is found to be significant on ZTD estimates (up to ±2 mm in equivalent PWV). The use of a proper ocean loading model eliminates this effect.The second aspect of this study concerns the IGS analysis quality for the African stations. The accuracy has been assessed through position dispersion between individual solutions and the most recent version of the IGS combined solution IGb00, and residuals from the transformation of the IGS combined solution into the International Terrestrial Reference Frame 2005. The positioning performance of the IGS analysis is consistent with an accuracy in ZTD of ±6 mm (±1 mm in PWV), as requested for meteorological applications such as planned in AMMA.     

Cause of scale inconsistencies in DORIS time series

In this paper we analyze the scale of the DORIS (Doppler Orbitography and Radiopositioning Integrated by Satellite) solutions with respect to DORIS extension of the International Terrestrial Reference Frame (ITRF) for Precise Orbit Determination DPOD2014. The main goal is to explain the scale inconsistencies and to find the optimal solution reaching low-biased and consistent scale time series. Our analysis profits from 4 different strategies based only on the Geodetic Observatory Pecný analysis center solution, using DORIS exchange format data 2.2. A difference in the sequence of the solutions directly corresponds to one of the changes in the solution settings: data elevation dependent weighting, application of data validity indicators and application of phase center - reference point correction. We process multi-satellite and single-satellite solutions for the time period 2011.0–2017.0. Our analysis examines scale inconsistency issues in 2011/2012 and in 2015. The scale increment in 2011/2012 is explained as a result of the concurrence of changes in satellite constellation and change in the provider data validity standards for Cryosat-2 and Jason-2 satellites. The scale increment in 2015 is explained as the effect of change in the standards for phase center - reference center corrections for Saral, Jason-2 and Cryosat-2 satellites. Moreover, comparing the solutions with and without elevation dependent data downweighting using the same elevation cutoff (10°), we found a significant reduction of scale bias and scale variation applying the data downweighting. The data downweighting improved also the station positioning repeatability. We demonstrate that the solution, which is completely free from the additional data associated with observations in DORIS exchange format 2.2 and includes the data downweighting law, eventuates in a consistent scale time series with the lowest offset with respect to DPOD2014 (version 1.0) (12.7 ± 2.3 mm for 2011.0–2017.0). The only remaining scale issue is the part of 2011/2012 increment of around 5 mm, explained by a change in the DORIS satellite constellation.     

利用精密单点定位求解电离层延迟

近年来,高时空分辨率的全球导航卫星系统(GNSS)观测信号已成为电离层研究的重要资源.利用GNSS研究电离层,需首先将观测资料转换成包含电离层信息的可观测量(Ionospheric 0bservables,称之为"电离层观测值").目前,最常用的电离层观测值一般采用联合无几何影响组合的码和相位观测,利用相位平滑伪距方法...