Combined Earth orientation parameters based on homogeneous and continuous VLBI and GPS data

The CONT02 campaign is of great interest for studies combining very long baseline interferometry (VLBI) with other space-geodetic techniques, because of the continuously available VLBI observations over 2 weeks in October 2002 from a homogeneous network. Especially, the combination with the Global Positioning System (GPS) offers a broad spectrum of common parameters. We combined station coordinates, Earth orientation parameters (EOPs) and troposphere parameters consistently in one solution using technique- specific datum-free normal equation systems. In this paper, we focus on the analyses concerning the EOPs, whereas the comparison and combination of the troposphere parameters and station coordinates are covered in a companion paper in Journal of Geodesy. In order to demonstrate the potential of the VLBI and GPS space-geodetic techniques, we chose a sub-daily resolution for polar motion (PM) and universal time (UT). A consequence of this solution set-up is the presence of a one-to-one correlation between the nutation angles and a retrograde diurnal signal in PM. The Bernese GPS Software used for the combination provides a constraining approach to handle this singularity. Simulation studies involving both nutation offsets and rates helped to get a deeper understanding of this singularity. With a rigorous combination of UT1–UTC and length of day (LOD) from VLBI and GPS, we showed that such a combination works very well and does not suffer from the systematic effects present in the GPS-derived LOD values. By means of wavelet analyses and the formal errors of the estimates, we explain this important result. The same holds for the combination of nutation offsets and rates. The local geodetic ties between GPS and VLBI antennas play an essential role within the inter-technique combination. Several studies already revealed non-negligible discrepancies between the terrestrial measurements and the space-geodetic solutions. We demonstrate to what extent these discrepancies propagate into the combined EOP solution.     

Error analysis of exterior orientation elements on geolocation for a Moon-based Earth observation optical sensor

ABSTRACT

The Moon is a potential new platform for Earth observation. The advantages of its large-scale observational scope, long temporal duration, and multi-layer detecting of the Earth will undoubtedly advance our understanding of the Earth system. To carry out the observations from a Moon-based optical sensor, the geolocation error caused by exterior orientation elements need to be investigated. This paper analyses the error effects of exterior orientation elements on geolocation for an optical sensor. To estimate the error, we present a geometric image model and utilise some parameters to measure the image offsets. Through a large number of numerical simulations, the results demonstrate that the image offsets are not obvious influenced by the distance and observation angle at mid-high latitude of the Moon and have linear correlation with the increasing errors of the exterior orientation elements. Further, the relationship between the spatial resolution and errors of exterior orientation elements are revealed. Finally, the error characteristics for Moon-based Earth observation are discussed. It is expected that the conclusion drawn in this paper could support the study of a Moon-based Earth observation optical sensor.     

Dependency of geodynamic parameters on the GNSS constellation

Significant differences in time series of geodynamic parameters determined with different Global Navigation Satellite Systems (GNSS) exist and are only partially explained. We study whether the different number of orbital planes within a particular GNSS contributes to the observed differences by analyzing time series of geocenter coordinates (GCCs) and pole coordinates estimated from several real and virtual GNSS constellations: GPS, GLONASS, a combined GPS/GLONASS constellation, and two virtual GPS sub-systems, which are obtained by splitting up the original GPS constellation into two groups of three orbital planes each. The computed constellation-specific GCCs and pole coordinates are analyzed for systematic differences, and their spectral behavior and formal errors are inspected. We show that the number of orbital planes barely influences the geocenter estimates. GLONASS’ larger inclination and formal errors of the orbits seem to be the main reason for the initially observed differences. A smaller number of orbital planes may lead, however, to degradations in the estimates of the pole coordinates. A clear signal at three cycles per year is visible in the spectra of the differences between our estimates of the pole coordinates and the corresponding IERS 08 C04 values. Combinations of two 3-plane systems, even with similar ascending nodes, reduce this signal. The understanding of the relation between the satellite constellations and the resulting geodynamic parameters is important, because the GNSS currently under development, such as the European Galileo and the medium Earth orbit constellation of the Chinese BeiDou system, also consist of only three orbital planes.     

Recent Improvements to IERS Bulletin A Combination and Prediction

Driven by a need for increased accuracy in real-time Earth orientation parameters (EOPs), the Bulletin A (Rapid Servce and Predictions) of the International Earth Rotation Service (IERS) has recently made several major changes to its combination and prediction procedures. Changes to the process ob combining multi-technique results include creation of a daily Bulletin A updata, inclusion of several new data sets, and use of polar motion rantes for the latest epoch. Notably, the contributions from GPS observations have grown steadily in significance, both for polar motion and Universal Time (UT1). The prediction procedure has, in turn, benefited from these changes as well as improvements to the polar motion prediction model. As a result, demanding real-time applications, such as for satellite orbit extrapolations should observe a major improvement in the accuracy of our real-time EOP products. All results, together with supporting and diagnostic information, are available at the website http://maia.usno.navy.mil. The maximum EOP errors (root-mean-squared sense) that a real-time user would experience using the latest available update of Bulletin A are currently estimated to be ∼0.9 milliarcseconds (mas) for polar motion and ∼0.15 milliseconds (ms) for UT1-UTC. The data latency (the lag since the most recent observations) for EOP predictions need not exceed ∼41 hours for users who avail themselves of the daily updates. Over the past four years, the accuracy for real-time applications has improved by nearly a factor of 4 in polar motion and a factor of 10 in UT1. This is primarily due to the large reduction in data latency, which in turn is mostly possible due to the Rapid product delivery of the International GPS Service (IGS) (see Mireault et al, 1999). ? 2001 John Wiley & Sons, Inc.     

Global GPS data analysis at the National Geodetic Survey

NOAA’s National Geodetic Survey (NGS) has been one of the Analysis Centers (ACs) of the International GNSS Service (IGS) since its inception in 1994. Solutions for daily GPS orbits and Earth orientation parameters are regularly contributed to the IGS Rapid and Final products, as well as solutions of weekly station positions. These solutions are combined with those of the other ACs and then the resultant IGS products are distributed to users. To perform these tasks, NGS has developed and refined the Program for the Adjustment of GPS EphemerideS (PAGES) software. Although PAGES has continuously evolved over the past 15 years, recent efforts have focused mostly on updating models and procedures to conform more closely to IGS and the International Earth Rotation Service (IERS) conventions. Details of our processing updates and demonstrations of the improvements will be provided.     

The international global navigation satellite systems service (IGS): development and achievements

Since 21 June 1992 the International GPS Service (IGS), renamed International GNSS Service in 2005, produces and makes available uninterrupted time series of its products, in particular GPS observations from the IGS Global Network, GPS orbits, Earth orientation parameters (components x and y of polar motion, length of day) with daily time resolution, satellite and receiver clock information for each day with different latencies and accuracies, and station coordinates and velocities in weekly batches for further analysis by the IERS (International Earth Rotation and Reference Systems Service). At a later stage the IGS started exploiting its network for atmosphere monitoring, in particular for ionosphere mapping, for troposphere monitoring, and time and frequency transfer. This is why new IGS products encompass ionosphere maps and tropospheric zenith delays. This development became even more important when more and more space-missions carrying space-borne GPS for various purposes were launched. This article offers an overview for the broader scientific community of the development of the IGS and of the spectrum of topics addressed today with IGS data and products.     

Generation of a consistent absolute phase-center correction model for GPS receiver and satellite antennas

The development and numerical values of the new absolute phase-center correction model for GPS receiver and satellite antennas, as adopted by the International GNSS (global navigation satellite systems) Service, are presented. Fixing absolute receiver antenna phase-center corrections to robot-based calibrations, the GeoForschungsZentrum Potsdam (GFZ) and the Technische Universität München reprocessed more than 10 years of GPS data in order to generate a consistent set of nadir-dependent phase-center variations (PCVs) and offsets in the z-direction pointing toward the Earth for all GPS satellites in orbit during that period. The agreement between the two solutions estimated by independent software packages is better than 1 mm for the PCVs and about 4 cm for the z-offsets. In addition, the long time-series facilitates the study of correlations of the satellite antenna corrections with several other parameters such as the global terrestrial scale or the orientation of the orbital planes with respect to the Sun. Finally, completely reprocessed GPS solutions using different phase-center correction models demonstrate the benefits from switching from relative to absolute antenna phase-center corrections. For example, tropospheric zenith delay biases between GPS and very long baseline interferometry (VLBI), as well as the drift of the terrestrial scale, are reduced and the GPS orbit consistency is improved.     

IGS Near Real-Time Products and Their Applications

The primary IGS products, including precise GPS orbits, Earth orientation parameters, and estimated and predicted GPS satellite clocks, are no longer used exclusively for essential geodetic support of scientic research. They are increasingly being used by a wide range of non-academic applications. In these applications, timeliness is extremely critical. To address the timeliness issue, the strengths and weaknesses of current IGS production processes are discussed, new ways to improve the timeliness and quality are explored, and recommendations are proposed to fulfill the application requirements. ? 2001 John Wiley & Sons, Inc.     

Geocenter variations derived from a combined processing of LEO- and ground-based GPS observations

GNSS observations provided by the global tracking network of the International GNSS Service (IGS, Dow et al. in J Geod 83(3):191–198, 2009) play an important role in the realization of a unique terrestrial reference frame that is accurate enough to allow a detailed monitoring of the Earth’s system. Combining these ground-based data with GPS observations tracked by high-quality dual-frequency receivers on-board low earth orbiters (LEOs) is a promising way to further improve the realization of the terrestrial reference frame and the estimation of geocenter coordinates, GPS satellite orbits and Earth rotation parameters. To assess the scope of the improvement on the geocenter coordinates, we processed a network of 53 globally distributed and stable IGS stations together with four LEOs (GRACE-A, GRACE-B, OSTM/Jason-2 and GOCE) over a time interval of 3 years (2010–2012). To ensure fully consistent solutions, the zero-difference phase observations of the ground stations and LEOs were processed in a common least-squares adjustment, estimating all the relevant parameters such as GPS and LEO orbits, station coordinates, Earth rotation parameters and geocenter motion. We present the significant impact of the individual LEO and a combination of all four LEOs on the geocenter coordinates. The formal errors are reduced by around 20% due to the inclusion of one LEO into the ground-only solution, while in a solution with four LEOs LEO-specific characteristics are significantly reduced. We compare the derived geocenter coordinates w.r.t. LAGEOS results and external solutions based on GPS and SLR data. We found good agreement in the amplitudes of all components; however, the phases in x- and z-direction do not agree well.     

A Discussion of IGS Solutions and Their Impact on Geodetic and Geophysical Applications

The International Association of Geodesy officially established the International GPS Service (IGS) on Janaury 1, 1994. Its prime objective is to provide support and a rerefence system for a wide variety of scientific and practical applications involving GPS. To fulfill its role the IGS also generates, in addition to its fundamental products (orbital/staion positions and consistent Earth orientation parameters), additional reference-system products providing the necessary infrastructure, standards, and means of calibrations for timing and various atmospheric applications of GPS. The generation and efficient application of IGS products and their impact on a number of positioning and atmospheric applications, including low earth orbit satellites, is reviewed and discussed. @ 1998 John Wiley & Sons, Inc.     

GEO-REFERENCING WITHOUT GROUND CONTROL

A geo-reference is a global or regional geographical or geodetic coor-dinate system to which sensors or spatial object data are related.Hence,geo-ref-erencing is close to the well known photogrammetric concept of exterior or abso-lute orientation,the common execution of which is indirect,via the use of groundcontrol point.GPS and INS technologies have changed the situation,permittingdirect measurement of position and attitude parameters and making exterior orien-tation feasible without ground control at all,in principle.The analysis of accuracyand reliability performance discloses,however,that especially INS does not yetmeet the high demands of photogrammetry.Moreover,control of systematic er-rors,the problem of datum transformation,and reliability conditions make theuse of some ground control points still mandatory,at least for any high perfor-mance geo-referencing.     

利用Allan方差分析GPS非差随机模型特性

非差GPS定位中,通常采用以观测值服从高斯白噪声条件的估计准则进行参数求解。研究表明卫星端误差、传播路径误差、测站环境误差等会破坏观测值的白噪声特性,并且未模型化误差同样具有不利影响。这不仅破坏了估计准则的假设条件,而且部分非白噪声有可能被状态参数吸收,影响估计的准确性。本文将观测值白噪声、有色噪声和未模型化误差一同纳入GPS非差随机模型,以验后残差来表征GPS数据的随机特性,进行Allan方差分析,研究噪声成分及其参数。结果表明,GPS非差噪声组合主要为WN+GM,相位白噪声为2.392mm,GM过程噪声为4.450mm/s,相关时间为52.074s,伪距白噪声为0.936m,GM过程噪声为0.833m/s,相关时间为14.737s,相位的GM过程噪声与卫星相关性较大,而其余噪声则与测站相关性较大,大量分析结果表明GPS非差随机模型并不服从高斯白噪声假设,有待精化。     

Characteristics of GPS positioning error with non-uniform pseudorange error

We study the characteristics of the random GPS positioning errors when the pseudorange errors differ for each satellite. A concise, explicit, analytical formula is derived for the covariance of the positioning error by using singular value decomposition. It is composed of a uniform error covariance together with additional contributions from those satellites with larger pseudorange errors. The eigenvectors of the uniform error covariance define the principal directions of the 4-dimensional error ellipsoid, and the eigenvalues are the squares of the semi-axes. The additional part from individual satellites has only one eigenvector and one eigenvalue. This makes the error ellipsoid enlarge mainly along a direction related to both the overall satellite geometry and the position of the specific satellites. The theory is validated by simulating the GPS constellation and pseudorange measurements. The random positioning error is examined while any one or more pseudorange errors are increased. Horizontal positioning error distributions are presented to demonstrate the variations of the orientation and size of the error ellipses with the pseudorange error of a specific satellite. The results show that the analytical formula describes the positioning error accurately.     

A strip adjustment procedure to mitigate the impact of inaccurate mounting parameters in parallel lidar strips

Lidar (laser scanning) technology has been proven as a prominent technique for the acquisition of high-density and accurate topographic information. Because of systematic errors in the lidar measurements (drifts in the position and orientation information and biases in the mirror angles and ranges) and/or in the parameters relating the system components (mounting parameters), adjacent lidar strips may exhibit discrepancies. Although position and orientation drifts can have a more significant impact, these errors and their impact do not come as a surprise if the quality of the GPS/INS integration process is carefully examined. Therefore, the mounting errors are singled out in this work. The ideal solution for improving the compatibility of neighbouring strips in the presence of errors in the mounting parameters is the implementation of a rigorous calibration procedure. However, such a calibration requires the original observations, which may not be usually available. In this paper, a strip adjustment procedure to improve the compatibility between parallel lidar strips with moderate flight dynamics (for example, acquired by a fixed-wing aircraft) over an area with moderately varying elevation is proposed. The proposed method is similar to the photogrammetric block adjustment of independent models. Instead of point features, planar patches and linear features, which are represented by sets of non-conjugate points, are used for the strip adjustment. The feasibility and the performance of the proposed procedure together with its impact on subsequent activities are illustrated using experimental results from real data.     

Strategies to mitigate aliasing of loading signals while estimating GPS frame parameters

Although GNSS techniques are theoretically sensitive to the Earth center of mass, it is often preferable to remove intrinsic origin and scale information from the estimated station positions since they are known to be affected by systematic errors. This is usually done by estimating the parameters of a linearized similarity transformation which relates the quasi-instantaneous frames to a long-term frame such as the International Terrestrial Reference Frame (ITRF). It is well known that non-linear station motions can partially alias into these parameters. We discuss in this paper some procedures that may allow reducing these aliasing effects in the case of the GPS techniques. The options include the use of well-distributed sub-networks for the frame transformation estimation, the use of site loading corrections, a modification of the stochastic model by downweighting heights, or the joint estimation of the low degrees of the deformation field. We confirm that the standard approach consisting of estimating the transformation over the whole network is particularly harmful for the loading signals if the network is not well distributed. Downweighting the height component, using a uniform sub-network, or estimating the deformation field perform similarly in drastically reducing the amplitude of the aliasing effect. The application of these methods to reprocessed GPS terrestrial frames permits an assessment of the level of agreement between GPS and our loading model, which is found to be about 1.5 mm WRMS in height and 0.8 mm WRMS in the horizontal at the annual frequency. Aliased loading signals are not the main source of discrepancies between loading displacement models and GPS position time series.     

Impact of the VLBA on reference frames and earth orientation studies

The geodetic VLBI community began using VLBA antennas in 1989 for geodesy and astrometry. We examine how usage of the VLBA has improved the celestial reference frame, the terrestrial reference frame, and Earth orientation parameters. Without the VLBA, ICRF2 would have had only 1011 sources instead of 3414. ICRF3 will contain at least 4121 sources, with approximately 70 % or more coming exclusively from VLBA astrometry and geodesy sessions. The terrestrial reference frame is also more stable and precise due to VLBA geodesy sessions. Approximately two dozen geodesy stations that have participated in VLBA sessions show average position formal errors that are \(\sim \)13–14 % better in the horizontal components and \(\sim \)5 % better in the vertical component than would be expected solely from the increased number of observations. Also the Earth orientation parameters obtained from the RDV sessions represent the most accurate EOP series of any of the long-term VLBI session types.     

GGOS-D: homogeneous reprocessing and rigorous combination of space geodetic observations

In preparation of activities planned for the realization of the Global Geodetic Observing System (GGOS), a group of German scientists has carried out a study under the acronym GGOS-D which closely resembles the ideas behind the GGOS initiative. The objective of the GGOS-D project was the investigation of the methodological and information-technological realization of a global geodetic-geophysical observing system and especially the integration and combination of the space geodetic observations. In the course of this project, highly consistent time series of GPS, VLBI, and SLR results were generated based on common state-of-the-art standards for modeling and parameterization. These series were then combined to consistently and accurately compute a Terrestrial Reference Frame (TRF). This TRF was subsequently used as the basis to produce time series of station coordinates, Earth orientation, and troposphere parameters. In this publication, we present results of processing algorithms and strategies for the integration of the space-geodetic observations which had been developed in the project GGOS-D serving as a prototype or a small and limited version of the data handling and processing part of a global geodetic observing system. From a comparison of the GGOS-D terrestrial reference frame results and the ITRF2005, the accuracy of the datum parameters is about 5?C7?mm for the positions and 1.0?C1.5?mm/year for the rates. The residuals of the station positions are about 3?mm and between 0.5 and 1.0?mm/year for the station velocities. Applying the GGOS-D TRF, the offset of the polar motion time series from GPS and VLBI is reduced to 50 ??as (equivalent to 1.5?mm at the Earth??s surface). With respect to troposphere parameter time series, the offset of the estimates of total zenith delays from co-located VLBI and GPS observations for most stations in this study is smaller than 1.5?mm. The combined polar motion components show a significantly better WRMS agreement with the IERS 05C04 series (96.0/96.0???as) than VLBI (109.0/100.7???as) or GPS (98.0/99.5???as) alone. The time series of the estimated parameters have not yet been combined and exploited to the extent that would be possible. However, the results presented here demonstrate that the experiences made by the GGOS-D project are very valuable for similar developments on an international level as part of the GGOS development.     

LEO GPS attitude determination algorithm for a micro-satellite using boom-arm deployed antennas

This paper describes a low earth orbiter micro-satellite attitude determination algorithm using GPS phase and pseudorange data as the only observables. It is designed to run in real-time, at a rate of 10 Hz, on-board the spacecraft, using minimal chip and memory resources. The spacecraft design includes four GPS antennas deployed on boom arms to improve the antenna separations. The boom arms feature smart sensors, from which time-varying deformation data are used to calculate changes in the body-fixed system (BFS) co-ordinates of the attitude antennas. These data are used as input to the attitude algorithm to improve the accuracy of the output. The conventional double-difference phase observation equations have been re-arranged so that the only unknown parameters in the functions (once the ambiguities have been determined) are the spacecraft Euler angles. This greatly increases the redundancy in the mathematical model, and is exploited to enhance the algorithm's ability to trap observations contaminated by unmodelled multipath. This approach has been shown to be successful in identifying phase outliers at the 5–10 mm level. Speed of execution of the program is improved by utilising numerical differentiation of the model equations in the linearisation process. Furthermore, as the number of solve-for parameters is reduced to three by the chosen mathematical model, matrix inversion requirements are minimised. A novel approach to ambiguity resolution and determination of initial estimates of the attitude parameters has been developed utilising a heuristic technique and the known, and time varying, BFS co-ordinates of the antenna array. Algorithm testing is based on a simulation of the micro-satellite trajectory combined with variations in attitude derived from spin-stabilisation and periodic roll and pitch parameters. The trajectory of the spacecraft centre of mass was calculated by numerical integration of a force model using Earth gravity field parameters, third body effects due to the Sun and the Moon, dynamic Earth tide effects (solar and lunar), and a solar radiation pressure model. Frame transformations between J2000 and ITRF97 used the IERS conventions. A similar approach was used to calculate the trajectories of all available GPS satellites during the same period, using initial conditions of position and velocity from IGS precise orbits. RMS differences between the published precise orbit and the integrated satellite positions were at the 5-mm level. Phase observables are derived from these trajectories, biased by simulation of receiver and satellite clock errors, cycle slips, random or systematic noise and initial integer ambiguities. In the actual simulation of the attitude determination process in orbit, GPS satellite positions are calculated using broadcast ephemerides. The results show that the aim of 0.05° (two sigma) attitude precision can be met provided that the phase noise can be reduced to the level of 1–2 mm. Attitude precision was found to vary strongly with constellation geometry, which can change quite rapidly depending on the variations in spacecraft attitude. The redundancy in the mathematical model was found to be very effective in trapping and isolating cycle slips to the double difference observations that are contaminated. This allows for the possibility of correcting for cycle slips without full recourse to the ambiguity resolution algorithm. Electronic Publication     

Impact of Earth radiation pressure on GPS position estimates

GPS satellite orbits available from the International GNSS Service (IGS) show a consistent radial bias of up to several cm and a particular pattern in the Satellite Laser Ranging (SLR) residuals, which are suggested to be related to radiation pressure mismodeling. In addition, orbit-related frequencies were identified in geodetic time series such as apparent geocenter motion and station displacements derived from GPS tracking data. A potential solution to these discrepancies is the inclusion of Earth radiation pressure (visible and infrared) modeling in the orbit determination process. This is currently not yet considered by all analysis centers contributing to the IGS final orbits. The acceleration, accounting for Earth radiation and satellite models, is introduced in this paper in the computation of a global GPS network (around 200 IGS sites) adopting the analysis strategies from the Center for Orbit Determination in Europe (CODE). Two solutions covering 9 years (2000–2008) with and without Earth radiation pressure were computed and form the basis for this study. In previous studies, it has been shown that Earth radiation pressure has a non-negligible effect on the GPS orbits, mainly in the radial component. In this paper, the effect on the along-track and cross-track components is studied in more detail. Also in this paper, it is shown that Earth radiation pressure leads to a change in the estimates of GPS ground station positions, which is systematic over large regions of the Earth. This observed “deformation” of the Earth is towards North–South and with large scale patterns that repeat six times per GPS draconitic year (350 days), reaching a magnitude of up to 1 mm. The impact of Earth radiation pressure on the geocenter and length of day estimates was also investigated, but the effect is found to be less significant as compared to the orbits and position estimates.     

三天线GPS姿态解算误差分析

介绍了一种通过三个天线的GPS测量定位值来求解载体姿态参数的方法,通过坐标系的变换得到了载体三个姿态角的求解公式;用误差传递的分析方法得到了载体三个姿态角测量误差的表达式;通过对GPS定位误差的简化处理,得到了载体三个姿态角测量误差的简化表达式;给出了一定条件下GPS定位误差的三个姿态参数测量误差的仿真结果。