利用空间大地测量数据探测地球膨胀效应

地球自转服务局(IERS)采用多种高精度的空间探测技术综合解算得到的国际地球参考框架(ITRF)是国际上公认的精度高、稳定性好的参考框架。为了研究地球的膨胀或收缩效应,本文采用ITRF2000的站坐标和速度,利用Delaunay算法生成的三角网逼近地球形体,计算出了地球的体积变化。     

A review of algebraic constraints in terrestrial reference frame datum definition

 Station coordinates are combined with velocities estimated by space geodesy techniques to produce the International Terrestrial Reference System. The input is sets of coordinates and velocities calculated by International Earth Rotation Service analysis centers using space geodesy techniques. The working reference system of individual analysis centers is generally conventionally defined. However, the implications of such processing can have an effect on the resulting combined set. The problem of datum definition as a function of coordinate combinations is reviewed. In particular, the problem of minimum constraints is clearly emphasized and the reference system effect is defined. The goal is to build a process that could be used generally to remove uncertainties in the underlying coordinate system without disturbing the underlying information with additional unnecessary information. Received: 25 January 1999 / Accepted: 20 September 2000     

Consistent realization of Celestial and Terrestrial Reference Frames

The Celestial Reference System (CRS) is currently realized only by Very Long Baseline Interferometry (VLBI) because it is the space geodetic technique that enables observations in that frame. In contrast, the Terrestrial Reference System (TRS) is realized by means of the combination of four space geodetic techniques: Global Navigation Satellite System (GNSS), VLBI, Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite. The Earth orientation parameters (EOP) are the link between the two types of systems, CRS and TRS. The EOP series of the International Earth Rotation and Reference Systems Service were combined of specifically selected series from various analysis centers. Other EOP series were generated by a simultaneous estimation together with the TRF while the CRF was fixed. Those computation approaches entail inherent inconsistencies between TRF, EOP, and CRF, also because the input data sets are different. A combined normal equation (NEQ) system, which consists of all the parameters, i.e., TRF, EOP, and CRF, would overcome such an inconsistency. In this paper, we simultaneously estimate TRF, EOP, and CRF from an inter-technique combined NEQ using the latest GNSS, VLBI, and SLR data (2005–2015). The results show that the selection of local ties is most critical to the TRF. The combination of pole coordinates is beneficial for the CRF, whereas the combination of \(\varDelta \hbox {UT1}\) results in clear rotations of the estimated CRF. However, the standard deviations of the EOP and the CRF improve by the inter-technique combination which indicates the benefits of a common estimation of all parameters. It became evident that the common determination of TRF, EOP, and CRF systematically influences future ICRF computations at the level of several \(\upmu \)as. Moreover, the CRF is influenced by up to \(50~\upmu \)as if the station coordinates and EOP are dominated by the satellite techniques.     

DORIS and the Determination of the Earth’s Polar Motion

DORIS (Détermination d’Orbite et Radiopositionnement Intégrés par Satellite) is a system used for precise orbit determination (POD) and ground-station positioning. It has been implemented on-board various satellites: the SPOT (Système pour l’Observation de la Terre) remote sensing satellites SPOT-2, SPOT-3, SPOT-4, SPOT-5, TOPEX/Poseidon and more recently on its successors Jason-1 and ENVISAT. DORIS is also a terrestrial positioning system that has found many applications in geophysics and geodesy; in particular, it contributes to the realization of the International Terrestrial Reference Frame, ITRF2000 and the forthcoming ITRF2005. Although not its primary objective, DORIS can bring information on Earth orientation monitoring, mainly polar motion and length of day (LOD) variations that complement other astrogeodetic techniques. In this paper, we have analyzed various recent polar motion solutions derived from independent analysis centers using different software packages and applying various analysis strategies. Comparisons of these solutions to the International Earth Rotation and Reference Systems Service (IERS) C04 solution are performed. Depending on the solutions, the accuracy of DORIS polar components are in the range of 0.5–1 mas corresponding to a few centimeters on the Earth’s surface. This is approximately ten times larger than results derived from GPS, which are typically 0.06 mas in both components. This does not allow DORIS results to be taken into account in the IERS–EOP combinations. A gain in the precision could come from technical improvements to the DORIS system, in addition to improvement of the orbit, tropospheric, ionospheric and Earth gravity field modeling.     

VLBI terrestrial reference frame contributions to ITRF2008

In late 2008, the Product Center for the International Terrestrial Reference Frame (ITRF) of the International Earth Rotation and Reference Systems Service (IERS) issued a call for contributions to the next realization of the International Terrestrial Reference System, ITRF2008. The official contribution of the International VLBI Service for Geodesy and Astrometry (IVS) to ITRF2008 consists of session-wise datum-free normal equations of altogether 4,539 daily Very Long Baseline Interferometry (VLBI) sessions from 1979.7 to 2009.0 including data of 115 different VLBI sites. It is the result of a combination of individual series of session-wise datum-free normal equations provided by seven analysis centers (ACs) of the IVS. All series are completely reprocessed following homogeneous analysis options according to the IERS Conventions 2003 and IVS Analysis Conventions. Altogether, nine IVS ACs analyzed the full history of VLBI observations with four different software packages. Unfortunately, the contributions of two ACs, Institute of Applied Astronomy (IAA) and Geoscience Australia (AUS), had to be excluded from the combination process. This was mostly done because the IAA series exhibits a clear scale offset while the solution computed from normal equations contained in the AUS SINEX files yielded unreliable results. Based on the experience gathered since the combination efforts for ITRF2005, some discrepancies between the individual series were discovered and overcome. Thus, the consistency of the individual VLBI solutions has improved considerably. The agreement in terms of WRMS of the Terrestrial Reference Frame (TRF) horizontal components is 1 mm, of the height component 2 mm. Comparisons between ITRF2005 and the combined TRF solution for ITRF2008 yielded systematic height differences of up to 5 mm with a zonal signature. These differences can be related to a pole tide correction referenced to a zero mean pole used by four of five IVS ACs in the ITRF2005 contribution instead of a linear mean pole path as recommended in the IERS Conventions. Furthermore, these systematics are the reason for an offset in the scale of 0.4 ppb between the IVS’ contribution to ITRF2008 and ITRF2005. The Earth orientation parameters of seven series used as input for the IVS combined series are consistent to a huge amount with about 50 μas WRMS in polar motion and 3 μs in dUT1.     

Quality assessment of GPS reprocessed terrestrial reference frame

The International GNSS Service (IGS) contributes to the construction of the International Terrestrial Reference Frame (ITRF) by submitting time series of station positions and Earth Rotation Parameters (ERP). For the first time, its submission to the ITRF2008 construction is based on a combination of entirely reprocessed GPS solutions delivered by 11 Analysis Centers (ACs). We analyze the IGS submission and four of the individual AC contributions in terms of the GNSS frame origin and scale, station position repeatability and time series seasonal variations. We show here that the GPS Terrestrial Reference Frame (TRF) origin is consistent with Satellite laser Ranging (SLR) at the centimeter level with a drift lower than 1 mm/year. Although the scale drift compared to Very Long baseline Interferometry (VLBI) and SLR mean scale is smaller than 0.4 mm/year, we think that it would be premature to use that information in the ITRF scale definition due to its strong dependence on the GPS satellite and ground antenna phase center variations. The new position time series also show a better repeatability compared to past IGS combined products and their annual variations are shown to be more consistent with loading models. The comparison of GPS station positions and velocities to those of VLBI via local ties in co-located sites demonstrates that the IGS reprocessed solution submitted to the ITRF2008 is more reliable and precise than any of the past submissions. However, we show that some of the remaining inconsistencies between GPS and VLBI positioning may be caused by uncalibrated GNSS radomes.     

顾及测量方式的地球自转参数差异性分析

随着甚长基线干涉测量(VLBI)、卫星激光测距(SLR)、激光测月(LLR)、全球卫星导航系统(GNSS)、多里斯系统(DORIS)等多种空间大地测量手段的使用,地球自转参数(ERP)的测量精度不断提高,为航天器导航、深空探测等诸多领域提供了高精度的国际天球参考系(ICRS)和国际地表参考系统(ITRS)之间的转换参数. 以国际地球自转与参考系服务发布的C04序列为基础序列,选取500天ERP序列,分析不同测量手段得到的ERP数据的误差分布情况,为研究利用不同数据之间的一致性进行精度检核的可行性及精度水平提供数据基础,同时也为ERP预报提供更多的数据选择.      

现代大地测量参考系统

概述现代大地测量参考系统的定义及不同参考系统之间的关系.主要讨论我国当代大地测量界常使用的3种用以表示几何位置的参考系统1980年国家大地坐标系、全球大地测量系统 1984 (WGS 84)、国际地球参考系统(ITRS),和一种用以表示物理位置,高程的参考系统1985年国家高程基准.并讨论大地测量中框架和基准的概念.     

高精度国际地球参考框架的研究进展

针对地球动力学等对毫米级国际地球参考框架的需求,介绍了国际地球参考框架的最新研究进展,探讨了现今顾及测站非线性运动的国际地球参考框架的局限性,并在毫米级国际地球参考框架建立的方法和技术改进方面提出了一些见解,对高精度国际地球参考框架的实现具有一定的参考价值。     

地球参考系的实现

本文简要讨论了组合不同空间技术建立全球最优协议地球参考架（ＣＴＲＦ）的理论和方法；重点分析了影响ＣＴＲＦ稳定的几种主要地球长期形变因素的性质和规律。首次考虑到地球质心运动对ＣＴＲＦ的影响，提出了协议地球质心的概念，给出了更为完整的四维地面点位模型。     

The introduction of the Iau 1980 nutation theory in the computation of the Earth Rotation Parameters by the Bureau International de l'Heure

Summary The status of nutation theory in the computation of the Earth Rotation Parameters from the different observational techniques is reviewed. The impact on the combined solution of the Bureau International de l'Heure is evaluated. In view of the improvement brought by the IAU 1980 Nutation Theory, its introduction in the BIH publications, without waiting for the adoption by the astronomical ephemerides (1984 Jan. 1) is decided.     

Geodesy and relativity

Relativity, or gravitational physics, has widely entered geodetic modelling and parameter determination. This concerns, first of all, the fundamental reference systems used. The Barycentric Celestial Reference System (BCRS) has to be distinguished carefully from the Geocentric Celestial Reference System (GCRS), which is the basic theoretical system for geodetic modelling with a direct link to the International Terrestrial Reference System (ITRS), simply given by a rotation matrix. The relation to the International Celestial Reference System (ICRS) is discussed, as well as various properties and relevance of these systems. Then the representation of the gravitational field is discussed when relativity comes into play. Presently, the so-called post-Newtonian approximation to GRT (general relativity theory) including relativistic effects to lowest order is sufficient for practically all geodetic applications. At the present level of accuracy, space-geodetic techniques like VLBI (Very Long Baseline Interferometry), GPS (Global Positioning System) and SLR/LLR (Satellite/Lunar Laser Ranging) have to be modelled and analysed in the context of a post-Newtonian formalism. In fact, all reference and time frames involved, satellite and planetary orbits, signal propagation and the various observables (frequencies, pulse travel times, phase and travel-time differences) are treated within relativity. This paper reviews to what extent the space-geodetic techniques are affected by such a relativistic treatment and where—vice versa—relativistic parameters can be determined by the analysis of geodetic measurements. At the end, we give a brief outlook on how new or improved measurement techniques (e.g., optical clocks, Galileo) may further push relativistic parameter determination and allow for refined geodetic measurements.     

Precise station positions from VLBI observations to satellites: a simulation study

Very long baseline interferometry (VLBI) tracking of satellites is a topic of increasing interest for the establishment of space ties. This shall strengthen the connection of the various space geodetic techniques that contribute to the International Terrestrial Reference Frame. The concept of observing near-Earth satellites demands research on possible observing strategies. In this paper, we introduce this concept and discuss its possible benefits for improving future realizations of the International Terrestrial Reference System. Using simulated observations, we develop possible observing strategies that allow the determination of radio telescope positions in the satellite system on Earth with accuracies of a few millimeters up to 1–2 cm for weekly station coordinates. This is shown for satellites with orbital heights between 2,000 and 6,000 km, observed by dense regional as well as by global VLBI-networks. The number of observations, as mainly determined by the satellite orbit and the observation interval, is identified as the most critical parameter that affects the expected accuracies. For observations of global positioning system satellites, we propose the combination with classical VLBI to radio sources or a multi-satellite strategy. Both approaches allow station position repeatabilities of a few millimeters for weekly solutions.     

Monitoring Earth orientation using space-geodetic techniques: state-of-the-art and prospective

Earth orientation parameters (EOPs) provide the transformation between the International Terrestrial Reference Frame (ITRF) and the International Celestial Reference Frame (ICRF). The different EOP series computed at the Earth Orientation Centre at the Paris Observatory are obtained from the combination of individual EOP series derived from the various space-geodetic techniques. These individual EOP series contain systematic errors, generally limited to biases and drifts, which introduce inconsistencies between EOPs and the terrestrial and celestial frames. The objectives of this paper are first to present the various combined EOP solutions made available at the EOP Centre for the different users, and second to present analyses concerning the long-term consistency of the EOP system with respect to both terrestrial and celestial reference frames. It appears that the present accuracy in the EOP combined IERS C04 series, which is at the level of 200 as for pole components and 20 s for UT1, does not match its internal precision, respectively 100 as and 5 s, because of propagation errors in the realization of the two reference frames. Rigorous combination methods based on a simultaneous estimation of station coordinates and EOPs, which are now being implemented within the International Earth Rotation Service (IERS), are likely to solve this problem in the future.     

Non-linear station motions in epoch and multi-year reference frames

In the conventions of the International Earth Rotation and Reference Systems Service (e.g. IERS Conventions 2010), it is recommended that the instantaneous station position, which is fixed to the Earth’s crust, is described by a regularized station position and conventional correction models. Current realizations of the International Terrestrial Reference Frame use a station position at a reference epoch and a constant velocity to describe the motion of the regularized station position in time. An advantage of this parameterization is the possibility to provide station coordinates of high accuracy over a long time span. Various publications have shown that residual non-linear station motions can reach a magnitude of a few centimeters due to not considered loading effects. Consistently estimated parameters like the Earth Orientation Parameters (EOP) may be affected if these non-linear station motions are neglected. In this paper, we investigate a new approach, which is based on a frequent (e.g. weekly) estimation of station positions and EOP from a combination of epoch normal equations of the space geodetic techniques Global Positioning System (GPS), Satellite Laser Ranging (SLR) and Very Long Baseline Interferometry (VLBI). The resulting time series of epoch reference frames are studied in detail and are compared with the conventional secular approach. It is shown that both approaches have specific advantages and disadvantages, which are discussed in the paper. A major advantage of the frequently estimated epoch reference frames is that the non-linear station motions are implicitly taken into account, which is a major limiting factor for the accuracy of the secular frames. Various test computations and comparisons between the epoch and secular approach are performed. The authors found that the consistently estimated EOP are systematically affected by the two different combination approaches. The differences between the epoch and secular frames reach magnitudes of $23.6~\upmu \hbox {as}$ (0.73 mm) and $39.8~\upmu \hbox {as}$ (1.23 mm) for the x-pole and y-pole, respectively, in case of the combined solutions. For the SLR-only solutions, significant differences with amplitudes of $77.3~\upmu \hbox {as}$ (2.39 mm) can be found.     

内核地球模型的形变场和地球参考系

讨论了内核地球模型的形变场（ 位移、速度） 和地球参考系的确立, 为内核地球动力学的研究提供了依据。     

地球参考系的基本理论和方法研究进展

地球参考系的研究无论在基础科学,还是在国民经济建设中都具有重要的意义。随着空间大地测量技术(VLBI、LLR、SLR和GPS等)的发展,对地球参考系也提出了越来越高的要求。本文详细地介绍了地球参考系的定义,建立和维持的基本理论与方法,综述了地球参考系目前的发展状况和存在的问题。本文可为从事地球参考系研究的学者提供一定的参考。     

基于IERS 2003规范的坐标系转换实现及其方案应用

主要介绍了IERS2003规范中国际地球参考系(ITRS)到地球质心天球参考系(GCRS)的转换过程,以及与IAU2000决议一致的常规转换公式,并详细分析了以IERS计算程序为工具、“基于CEO”[注]和“基于春分点”的ITRS向GCRS转换的两种等价方案和相应的三种应用方法,最后对新坐标系转换模型的应用进行了探讨。     

Effects of adopting new precession,nutation and equinox corrections on the terrestrial reference frame

First, the paper is devoted to the effects of adopting new definitive precession and equinox corrections on the terrestrial reference frame: The effect on polar motion is a diurnal periodic term with an amplitude increasing linearly in time; on UT1 it is a linear term. Second, general principles are given the use of which can determine the effects of small rotations (such as precession, nutation or equinox corrections) of the frame of a Conventional Inertial Reference System (CIS) on the frame of the Conventional Terrestrial Reference System (CTS). Next, seven CTS options are presented, one of which is necessary to accommodate such rotations (corrections). The last of these options requiring no changes in the origin of terrestrial longitudes and in UT1 is advocated; this option would be maintained by eventually referencing the Greenwich Mean Sidereal Time to a fixed point on the equator, instead of to the mean equinox of date, the current practice. Accommodating possible future changes in the astronomical nutation is discussed in the last section. The Appendix deals with the effects of differences which may exist between the various CTS's and CIS's (inherent in the various observational techniques) on earth rotation parameters (ERP) and how these differences can be determined. It is shown that the CTS differences can be determined from observations made at the same site, while the CIS differences by comparing the ERP's determined by the different techniques during the same time period. Presented at XVIII General Assembly of the International Astronomical Union Patras, Greece, August 17–26, 1982.     

Combination procedure for length-of-day time series according to the noise frequency behavior

 Length-of-day (LOD) estimates from seven Global Positioning System (GPS) and three satellite laser ranging (SLR) analysis centers were combined into an even-spaced time series for a 27-month period (1996–1998). This time series was compared to the multi-technique Earth-orientation-parameter (EOP) combined solution (C04) derived at the Central Bureau of the International Earth Rotation Service (IERS/CB). Due to inhomogeneities in the different series derived from the various techniques (time, length, quality, and spatial resolution), the concept of a combined solution is justified. The noise behavior in LOD for different techniques varies with frequency; the data series were divided into frequency windows after removing both biases and trends. Different weight factors were assigned in each window, discriminating by technique, and produced one-technique combined solutions. Finally, these one-technique combined solutions were combined to obtain the final multi-technique solution. The LOD combined time series obtained by the present method based on the frequency windows combined series (FWCS) is very close to the IERS C04 solution. It could be useful to generate a new LOD reference time series to be used in the study of high-frequency variations of Earth rotation. Received: 28 March 2000 / Accepted: 15 February 2001