共查询到20条相似文献,搜索用时 31 毫秒
1.
Low-degree earth deformation from reprocessed GPS observations 总被引:3,自引:1,他引:2
Mathias Fritsche R. Dietrich A. Rülke M. Rothacher P. Steigenberger 《GPS Solutions》2010,14(2):165-175
Surface mass variations of low spherical harmonic degree are derived from residual displacements of continuously tracking
global positioning system (GPS) sites. Reprocessed GPS observations of 14 years are adjusted to obtain surface load coefficients
up to degree n
max = 6 together with station positions and velocities from a rigorous parameter combination. Amplitude and phase estimates of
the degree-1 annual variations are partly in good agreement with previously published results, but also show interannual differences
of up to 2 mm and about 30 days, respectively. The results of this paper reveal significant impacts from different GPS observation
modeling approaches on estimated degree-1 coefficients. We obtain displacements of the center of figure (CF) relative to the
center of mass (CM), Δr
CF–CM, that differ by about 10 mm in maximum when compared to those of the commonly used coordinate residual approach. Neglected
higher-order ionospheric terms are found to induce artificial seasonal and long-term variations especially for the z-component of Δr
CF–CM. Daily degree-1 estimates are examined in the frequency domain to assess alias contributions from model deficiencies with
regard to satellite orbits. Finally, we directly compare our estimated low-degree surface load coefficients with recent results
that involve data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission. 相似文献
2.
Considering present attempts to develop a gradiometer with an accuracy between 10−3
E and 10−4
E, two applications for such a device have been studied: (a) mapping the gravitational field of the Earth, and (b) estimating
the geocentric distance of a satellite carrying the instrument. Given a certain power spectrum for the signal and 10−4
E (rms) of white measurement noise, the results of an error analysis indicate that a six-month mission in polar orbit at a height
of 200 km, with samples taken every three seconds, should provide data for estimating the spherical harmonic potential coefficients
up to degree and order 300 with less than 50% error, and improve the coefficients through degree 30 by up to four orders of
magnitude compared to existing models. A simulation study based on numerical orbit integrations suggests that a simple adjustment
of the initial conditions based on gradiometer data could produce orbits where the geocentric distance is accurate to 10 cm
or better, provided the orbits are 2000 km high and some improvement in the gravity field up to degree 30 is first achieved.
In this sense, the gravity-mapping capability of the gradiometer complements its use in orbit refinement. This idea can be
of use in determining orbits for satellite altimetry. Furthermore, by tracking the gradiometer-carrying spacecraft when it
passes nearly above a terrestrial station, the geocentric distance of this station can also be estimated to about one decimeter
accuracy. This principle could be used in combination with VLBI and other modern methods to set up a world-wide 3-D network
of high accuracy. 相似文献
3.
We can presently construct two independent time series of sea level, each at a precision of a few centimeters, from Geosat
(1985–1988) and TOPEX/Poseidon (1992–1995) collinear altimetry. Both are based on precise satellite orbits computed using
a common geopotential model, JGM-2 (Nerem et al. 1994). We have attempted to connect these series using Geosat-T/P crossover
differences in order to assess long-term ocean changes between these missions. Unfortunately, the observed result are large-scale
sea level differences which appear to be due to a combination of geodetic and geopotential error sources. The most significant
geodetic component seems to be a coordinate system bias for Geosat sea level (relative to T/P) of −7 to −9 cm in the y-axis (towards the Eastern Pacific). The Geosat-T/P sea height differences at crossovers (with JGM-2 orbits) probably also
contain stationary geopotential-orbit error of about the same magnitude which also distort any oceanographic interpretation
of the observed changes. We also found JGM-3 Geosat orbits have not resolved the datum errors evident from the JGM-2 Geosat
-T/P results. We conclude that the direct altimetric approach to accurate determination of sea level change using Geosat and
T/P data still depends on further improvement in the Geosat orbits, including definition of the geocenter.
Received: 11 March 1996; Accepted: 19 September 1996 相似文献
4.
G. Beutler W. Gurtner M. Rothacher T. Schildknecht I. Bauersima 《Journal of Geodesy》1986,60(3):205-220
Summary Carrier phase measurements are potentially the most precise observations available from theGPS satellite system, the formal precision being of the order of one centimeter per observation. If the so called double differences
are used as the basic observable, the analysis is relatively simple, since satellite- and receiver-clocks may be represented
by basic models. We investigate the feasibility of double difference phase observations for orbit determination using the
material of the 1985 High Precision Baseline Test, where the coordinates of the so called fiducial points (Haystack, Ft. Davis
Richmond and Mojave) are held fixed.TI-4100 andAFGL-receiver observations were used in the same orbit determination process.
Although no surface weather data had been available to us, the orbit quality seems to be of the order of0.1 ppm. When we use these orbits to estimate the coordinates of the five “non-fiducial points” Owens Valley, Hat Creek Mammoth Lake,
Austin and Dahlgren we get a repeatability of the order of5 cm for latitude and longitude and10 cm for height, if the observations of the first four days of the campaign are compared to those of the second four days.
If we use our orbits estimated withTI andAFGL observations to process the Mojave—Owens Valley baseline (length245 km) measured by the twoSERIES-X receivers, we obtain day to day repeatabilities of1.6 cm (0.06 ppm) in length,2 cm (0.08 ppm) in latitude,4 cm (0.16 ppm) in longitude and7 cm (0.29 ppm) in height.
Since there are indications that regional networks will be realized in the near future, the results presented here should
encourage the realization of regional high precision orbit determination services. 相似文献
5.
The technique of Evaluating CHAMP satellite orbit with SLR measurements is presented. As an independent evaluation of the orbit solution, SLR data observed from January 1 to 16, 2002 are processed to compute the residuals after fixing the GFZ's post science orbits solutions. The SLR residuals are computed as the differences of the SLR measurements minus the corresponding distances between the SLR station and the GPS-derived orbit positions. On the basis of the SLR residuals analysis, it is found that the accuracy of GFZ' s spost science orbits is better than 10 cm and that there is no systematic error in GFZ's post science orbits. 相似文献
6.
Evaluation of CHAMP satellite orbit with SLR measurements 总被引:1,自引:0,他引:1
QIN Xianping YANG Yuanxi 《地球空间信息科学学报》2005,8(3):209-213
The technique of Evaluating CHAMP satellite orbit with SLR measurements iS presented. As an independent evaluation of the orbit solution, SLR data observed from January 1 to 16, 2002 are processed to compute the residuals after fixing the GFZ's post science orbits solutions. The SLR residuals are computed as the differences of the SLR measurements minus the corresponding distances between the SLR station and the GPS-derived orbit positions. On the basis of the SLR residuals analysis, it is found that the accuracy of GFZ' s post science orbits is better than 10 cm and that there is no systematic error in GFZ's post science orbits. 相似文献
7.
Parametric least squares collocation was used in order to study the detection of systematic errors of satellite gradiometer
data. For this purpose, simulated data sets with a priori known systematic errors were produced using ground gravity data
in the very smooth gravity field of the Canadian plains. Experiments carried out at different satellite altitudes showed that
the recovery of bias parameters from the gradiometer “measurements” is possible with high accuracy, especially in the case
of crossing tracks. The mean value of the differences (original minus estimated bias parameters) was relatively large compared
to the standard deviation of the corresponding second-order derivative component at the corresponding height. This mean value
almost vanished when gravity data at ground level were combined with the second-order derivative data set at satellite altitude.
In the case of simultaneous estimation of bias and tilt parameters from ∂2
T/∂z
2“measurements”, the recovery of both parameters agreed very well with the collocation error estimation.
Received: 10 October 1996 / Accepted 25 May 1998 相似文献
8.
Estimability analysis of variance and covariance components 总被引:1,自引:1,他引:1
Although variance and covariance components have been extensively investigated and a number of elegant formulae to compute
them have been derived, nothing is known, without any ambiguity, about their estimability in the case of a fully unknown variance–covariance
matrix. We prove that variance and covariance components in this case are not estimable, thus clarifying the ambiguity of
the literature on the topic and correcting some erroneous statements in the literature. We also give a new theorem on the
estimability of a linear function of variance and covariance components. Then we propose a new method to estimate the variance–covariance
matrix with special structure, which can presumably be represented by, at most, r(r + 1)/2 independent parameters to guarantee its estimability in such a subspace, by directly implementing the positive definiteness
of the matrix as constraint to the restricted maximum likelihood method, where r is the number of redundant measurements. Therefore, our estimates of the variance and covariance components always reconstruct
a positive definite matrix and are always physically meaningful. 相似文献
9.
The problem of “global height datum unification” is solved in the gravity potential space based on: (1) high-resolution local
gravity field modeling, (2) geocentric coordinates of the reference benchmark, and (3) a known value of the geoid’s potential.
The high-resolution local gravity field model is derived based on a solution of the fixed-free two-boundary-value problem
of the Earth’s gravity field using (a) potential difference values (from precise leveling), (b) modulus of the gravity vector
(from gravimetry), (c) astronomical longitude and latitude (from geodetic astronomy and/or combination of (GNSS) Global Navigation
Satellite System observations with total station measurements), (d) and satellite altimetry. Knowing the height of the reference
benchmark in the national height system and its geocentric GNSS coordinates, and using the derived high-resolution local gravity
field model, the gravity potential value of the zero point of the height system is computed. The difference between the derived
gravity potential value of the zero point of the height system and the geoid’s potential value is computed. This potential
difference gives the offset of the zero point of the height system from geoid in the “potential space”, which is transferred
into “geometry space” using the transformation formula derived in this paper. The method was applied to the computation of
the offset of the zero point of the Iranian height datum from the geoid’s potential value W
0=62636855.8 m2/s2. According to the geometry space computations, the height datum of Iran is 0.09 m below the geoid. 相似文献
10.
The crossover adjustment plays a central role in the processing of satellite altimeter measurements. The usual procedure
is to form sea surface height differences at crossover points, solve for the radial orbit error (with due attention to the
singular nature of the estimation problem) and then to construct altimetric sea-level maps using the mean sea surface heights
at the crossover points. Our approach is very different, to make direct use of measurements at crossover points without differencing
and to estimate simultaneously orbit parameters, mean sea surface height and sea surface height variability in a single, unified
adjustment. The technique is suited for repeat data over an area small enough that adjoining passes may be considered to be
parallel and to permit the solution of a set of linear equations of dimension equal to the number of crossover points. The
size of the numerical problem is almost independent of the number of repeat cycles of the altimeter mission. Explicit recognition
is given to the rank defect of the least-squares estimation problem; we show that, for an orbit model with r parameters, the rank defect of the local crossover problem is exactly r
2. The defect may be overcome by choosing an appropriate set of constraints – either giving a best fit of mean sea surface
heights to a reference surface, or minimising orbit parameters, or a minimum norm solution in which both mean sea surface
heights and orbit parameters are minimised. There is no need to choose a reference pass, all passes are treated equally and
data gaps are easily accommodated. Numerical results are presented for the south-western Indian Ocean, based on the first
2 years of altimeter data from the Geosat Exact Repeat Mission.
Received: 31 May 1996 / Accepted: 19 April 1997 相似文献
11.
Since the beginning of the International Global Navigation Satellite System (GLONASS) Experiment, IGEX, in October 1998,
the Center for Orbit Determination in Europe (CODE) has acted as an analysis center providing precise GLONASS orbits on a
regular basis. In CODE's IGEX routine analysis the Global Positioning System (GPS) orbits and Earth rotation parameters are
introduced as known quantities into the GLONASS processing. A new approach is studied, where data from the IGEX network are
combined with GPS observations from the International GPS Service (IGS) network and all parameters (GPS and GLONASS orbits,
Earth rotation parameters, and site coordinates) are estimated in one processing step. The influence of different solar radiation
pressure parameterizations on the GLONASS orbits is studied using different parameter subsets of the extended CODE orbit model.
Parameterization with three constant terms in the three orthogonal directions, D, Y, and X (D = direction satellite–Sun, Y = direction of the satellite's solar panel axis), and two periodic terms in the X-direction, proves to be adequate for GLONASS satellites. As a result of the processing it is found that the solar radiation
pressure effect for the GLONASS satellites is significantly different in the Y-direction from that for the GPS satellites, and an extensive analysis is carried out to investigate the effect in detail.
SLR observations from the ILRS network are used as an independent check on the quality of the GLONASS orbital solutions. Both
processing aspects, combining the two networks and changing the orbit parameterization, significantly improve the quality
of the determined GLONASS orbits compared to the orbits stemming from CODE's IGEX routine processing.
Received: 10 May 2000 / Accepted: 9 October 2000 相似文献
12.
M. K. Paul 《Journal of Geodesy》1978,52(3):177-190
Recurrence relations for the evaluation of the integrals of associated Legendre functions over an arbitrary interval within
(0°, 90°) have been derived which yield sufficiently accurate results throughout the entire range of their possible applications.
These recurrence relations have been used to compute integrals up to degree 100 and similar computations can be carried out
without any difficulty up to a degree as high as the memory in a computer permits. The computed values have been tested with
independent check formulae, also derived in this work; the corresponding relative errors never exceed 10−23 in magnitude.
Contribution from the Earth Physics Branch No. 719 相似文献
13.
Differential tracking of theGPS satellites in high-earth orbit provides a powerful relative positioning capability, even when a relatively small continental
U.S. fiducial tracking network is used with less than one-third of the fullGPS constellation. To demonstrate this capability, we have determined baselines of up to2000 km in North America by estimating high-accuracyGPS orbits and ground receiver positions simultaneously. The2000 km baselines agree with very long baseline interferometry(VLBI) solutions at the level of1.5 parts in10
8 and showrms daily repeatability of0.3–2 parts in10
8. The orbits determined for the most thoroughly trackedGPS satellites are accurate to better than1 m. GPS orbit accuracy was assessed from orbit predictions, comparisons with independent data sets, and the accuracy of the continental
baselines determined along with the orbits. The bestGPS orbit strategies included data arcs of at least one week, process noise models for tropospheric fluctuations, estimation
ofGPS solar pressure coefficients, and combined processing ofGPS carrier phase and pseudorange data. For data arcs of two weeks, constrained process noise models forGPS dynamic parameters significantly improved the solutions. 相似文献
14.
J. Klokočník Ch. Reigber P. Schwintzer C. A. Wagner J. Kostelecký 《Journal of Geodesy》2002,76(4):189-198
The new GFZ/GRGS gravity field models GRIM5-S1 and GRIM5-C1, currently used as initial models for the CHAMP mission, have
been compared with other recent models (JGM 3, EGM 96) for radial orbit accuracy (by means of latitude lumped coefficients)
in computations on altimetry satellite orbits. The bases for accuracy judgements are multi-year averages of crossover sea
height differences from Geosat and ERS 1/2 missions. This radially sensitive data is fully independent of the data used to
develop these gravity models. There is good agreement between the observed differences in all of the world's oceans and projections
of the same errors from the scaled covariance matrix of their harmonic geopotential coefficients. It was found that the tentative
scale factor of five for the formal standard deviations of the harmonic coefficients of the new GRIM fields is justified,
i.e. the accuracy estimates, provided together with the GRIM geopotential coefficients, are realistic.
Received: 20 February 2001 / Accepted: 24 October 2001 相似文献
15.
R. Dach G. Beutler U. Hugentobler S. Schaer T. Schildknecht T. Springer G. Dudle L. Prost 《Journal of Geodesy》2003,77(1-2):1-14
A joint time-transfer project between the Astronomical Institute of the University of Berne (AIUB) and the Swiss Federal
Office of Metrology and Accreditation (METAS) was initiated to investigate the power of the time transfer using GPS carrier
phase observations. Studies carried out in the context of this project are presented. The error propagation for the time-transfer
solution using GPS carrier phase observations was investigated. To this purpose a simulation study was performed. Special
interest was focussed on errors in the vertical component of the station position, antenna phase-center variations and orbit
errors. A constant error in the vertical component introduces a drift in the time-transfer results for long baselines in east–west
directions. The simulation study was completed by investigating the profit for time transfer when introducing the integer
carrier phase ambiguities from a double-difference solution. This may reduce the drift in the time-transfer results caused
by constant vertical error sources. The results from the present time-transfer solution are shown in comparison to results
obtained with independent time-transfer techniques. The interpretation of the comparison benefits from the investigations
of the error propagation study. Two types of solutions are produced on a regular basis at AIUB: one based on the rapid orbits
from CODE, the other on the CODE final orbits. The rapid solution is available the day after the observations and has nearly
the same quality as the final solution, which has a latency of about one week. The differences between these two solutions
are below the nanosecond level. The differences from independent time-transfer techniques such as TWSTFT (two-way satellite
time and frequency transfer) are a few nanoseconds for both products.
Received: 15 November 2001 / Accepted: 6 September 2002
Correspondence to:R. Dach 相似文献
16.
A new gravity map, a new marine geoid around Japan and the detection of the Kuroshio current 总被引:3,自引:0,他引:3
About half a million marine gravity measurements over a 30∘×30∘ area centered on Japan have been processed and adjusted to produce a new free-air gravity map from a 5′×5′ grid. This map
seems to have a better resolution than those previously published as measured by its correlation with bathymetry. The grid
was used together with a high-degree and -order spherical harmonics geopotential model to compute a detailed geoid with two
methods: Stokes integral and collocation. Comparisons with other available geoidal surfaces derived either from gravity or
from satellite altimetry were made especially to test the ability of this new geoid at showing the sea surface topography
as mapped by the Topex/Poseidon satellite. Over 2 months (6 cycles) the dynamic topography at ascending passes in the region
(23∘47∘N and 123∘147∘E) was mapped to study the variability of the Kuroshio current.
Received: 15 July 1994 / Accepted: 17 February 1997 相似文献
17.
Based on the gravity field models EGM96 and EIGEN-GL04C, the Earth's time-dependent principal moments of inertia A, B, C are obtained, and the variable rotation of the Earth is determined. Numerical results show that A, B, and C have increasing tendencies; the tilt of the rotation axis increases 2.1×10^ 8 mas/yr; the third component of the rotational angular velocity, ω3 , has a decrease of 1.0×10^ 22 rad/s^2, which is around 23% of the present observed value. Studies show in detail that both 0 and ω3 experience complex fluctuations at various time scales due to the variations of A, B and C. 相似文献
18.
The spacetime gravitational field of a deformable body 总被引:3,自引:0,他引:3
The high-resolution analysis of orbit perturbations of terrestrial artificial satellites has documented that the eigengravitation
of a massive body like the Earth changes in time, namely with periodic and aperiodic constituents. For the space-time variation
of the gravitational field the action of internal and external volume as well as surface forces on a deformable massive body
are responsible. Free of any assumption on the symmetry of the constitution of the deformable body we review the incremental
spatial (“Eulerian”) and material (“Lagrangean”) gravitational field equations, in particular the source terms (two constituents:
the divergence of the displacement field as well as the projection of the displacement field onto the gradient of the reference
mass density function) and the `jump conditions' at the boundary surface of the body as well as at internal interfaces both
in linear approximation. A spherical harmonic expansion in terms of multipoles of the incremental Eulerian gravitational potential
is presented. Three types of spherical multipoles are identified, namely the dilatation multipoles, the transport displacement
multipoles and those multipoles which are generated by mass condensation onto the boundary reference surface or internal interfaces.
The degree-one term has been identified as non-zero, thus as a “dipole moment” being responsible for the varying position
of the deformable body's mass centre. Finally, for those deformable bodies which enjoy a spherically symmetric constitution,
emphasis is on the functional relation between Green functions, namely between Fourier-/ Laplace-transformed volume versus
surface Love-Shida functions (h(r),l(r) versus h
′(r),l
′(r)) and Love functions k(r) versus k
′(r). The functional relation is numerically tested for an active tidal force/potential and an active loading force/potential,
proving an excellent agreement with experimental results.
Received: December 1995 / Accepted: 1 February 1997 相似文献
19.
Accuracy of GPS-derived relative positions as a function of interstation distance and observing-session duration 总被引:6,自引:0,他引:6
Ten days of GPS data from 1998 were processed to determine how the accuracy of a derived three-dimensional relative position
vector between GPS antennas depends on the chord distance (denoted L) between these antennas and on the duration of the GPS observing session (denoted T). It was found that the dependence of accuracy on L is negligibly small when (a) using the `final' GPS satellite orbits disseminated by the International GPS Service, (b) fixing
integer ambiguities, (c) estimating appropriate neutral-atmosphere-delay parameters, (d) 26 km ≤ L ≤ 300 km, and (e) 4 h ≤T ≤ 24 h. Under these same conditions, the standard error for the relative position in the north–south dimension (denoted S
n
and expressed in mm) is adequately approximated by the equation S
n
=k
n
/T
0.5 with k
n
=9.5 ± 2.1 mm · h0.5 and T expressed in hours. Similarly, the standard errors for the relative position in the east–west and in the up-down dimensions
are adequately approximated by the equations S
e
=k
e
/T
0.5 and S
u
=k
u
/T
0.5, respectively, with k
e
=9.9 ± 3.1 mm · h0.5 and k
u
=36.5 ± 9.1 mm · h0.5.
Received: 5 February 2001 / Accepted: 14 May 2001 相似文献
20.
The identification of mean semi-major axes (suitably defined) for satellite orbits to satisfy a variety of requirements for
geodesy, geophysics and oceanography, in terms of repeat orbits (with orbital resonances), is investigated. Various options
for the definition of semi-major axis, from the viewpoint of satellite dynamics, are described. Simple simulations of the
expected resonant changes in inclination are presented, and tools for the analysis of orbit resonances to extract certain
lumped harmonic coefficients of the geopotential (e.g. from the very precise CHAMP orbit) are resurrected. Finally, a preliminary
example of the 46th-order resonance analysis possible for CHAMP, based on the mean orbital elements produced by GFZ (GeoForschungs
Zentrum) for ephemeris prediction, is presented.
Received: 10 July 2001 / Accepted: 17 July 2002
Correspondence to: J. Klokočník at Ondřejov Observatory
Acknowledgements. We thank Prof. Dr. Ch. Reigber, Dr. P. Schwintzer, Dr. T. Gruber and Dr. R. K?nig from GFZ Potsdam for various consultations
and discussions, and for the CHAMP two-line mean elements. This investigation was performed under the aegis of CEDR (Center
for Earth's Dynamics Research, Prague-Ondřejov); it has been supported by project LN00A005 (provided by the Ministry of Education
of the Czech Republic) and by grant A 3004 of the Grant Agency of the Academy of Sciences of the Czech Republic. 相似文献