共查询到20条相似文献,搜索用时 30 毫秒
1.
The performance of the L-curve criterion and of the generalized cross-validation (GCV) method for the Tikhonov regularization
of the ill-conditioned normal equations associated with the determination of the gravity field from satellite gravity gradiometry
is investigated. Special attention is devoted to the computation of the corner point of the L-curve, to the numerically efficient
computation of the trace term in the GCV target function, and to the choice of the norm of the residuals, which is important
for the Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) in the presence of colored observation noise. The
trace term in the GCV target function is estimated using an unbiased minimum-variance stochastic estimator. The performance
analysis is based on a simulation of gravity gradients along a 60-day repeat circular orbit and a gravity field recovery complete
up to degree and order 300. Randomized GCV yields the optimal regularization parameter in all the simulations if the colored
noise is properly taken into account. Moreover, it seems to be quite robust against the choice of the norm of the residuals.
It performs much better than the L-curve criterion, which always yields over-smooth solutions. The numerical costs for randomized
GCV are limited provided that a reasonable first guess of the regularization parameter can be found.
Received: 17 May 2001 / Accepted: 17 January 2002 相似文献
2.
The recovery of a full set of gravity field parameters from satellite gravity gradiometry (SGG) is a huge numerical and computational
task. In practice, parallel computing has to be applied to estimate the more than 90 000 harmonic coefficients parameterizing
the Earth's gravity field up to a maximum spherical harmonic degree of 300. Three independent solution strategies (preconditioned
conjugate gradient method, semi-analytic approach, and distributed non-approximative adjustment), which are based on different
concepts, are assessed and compared both theoretically and on the basis of a realistic-as-possible numerical simulation regarding
the accuracy of the results, as well as the computational effort. Special concern is given to the correct treatment of the
coloured noise characteristics of the gradiometer. The numerical simulations show that the three methods deliver nearly identical
results—even in the case of large data gaps in the observation time series. The newly proposed distributed non-approximative
adjustment approach, which is the only one of the three methods that solves the inverse problem in a strict sense, also turns
out to be a feasible method for practical applications.
Received: 17 December 2001 / Accepted: 17 July 2002
Acknowledgments. We would like to thank Prof. W.-D. Schuh, Institute of Theoretical Geodesy, University of Bonn, for providing us with the
serial version of the PCGMA algorithm, which forms the basis for the parallel PCGMA package developed at our institute. This
study was partially performed in the course of the GOCE project `From E?tv?s to mGal+', funded by the European Space Agency
(ESA) under contract No. 14287/00/NL/DC.
Correspondence to: R. Pail 相似文献
3.
GPS-based precise orbit determination of the very low Earth-orbiting gravity mission GOCE 总被引:5,自引:0,他引:5
A prerequisite for the success of future gravity missions like the European Gravity field and steady-state Ocean Circulation
Explorer (GOCE) is a precise orbit determination (POD). A detailed simulation study has been carried out to assess the achievable
orbit accuracy based on satellite-to-satellite tracking (SST) by the US global positioning system (GPS) and in conjunction
the implications for gravity field determination. An orbit accuracy at the few centimeter level seems possible, sufficient
to support the GOCE gravity mission and in particular its gravity gradiometer.
Received: 21 January 2000 / Accepted: 4 July 2000 相似文献
4.
Exploring gravity field determination from orbit perturbations of the European Gravity Mission GOCE 总被引:5,自引:0,他引:5
A comparison was made between two methods for gravity field recovery from orbit perturbations that can be derived from global
positioning system satellite-to-satellite tracking observations of the future European gravity field mission GOCE (Gravity
Field and Steady-State Ocean Circulation Explorer). The first method is based on the analytical linear orbit perturbation
theory that leads under certain conditions to a block-diagonal normal matrix for the gravity unknowns, significantly reducing
the required computation time. The second method makes use of numerical integration to derive the observation equations, leading
to a full set of normal equations requiring powerful computer facilities. Simulations were carried out for gravity field recovery
experiments up to spherical harmonic degree and order 80 from 10 days of observation. It was found that the first method leads
to large approximation errors as soon as the maximum degree surpasses the first resonance orders and great care has to be
taken with modeling resonance orbit perturbations, thereby loosing the block-diagonal structure. The second method proved
to be successful, provided a proper division of the data period into orbital arcs that are not too long.
Received: 28 April 2000 / Accepted: 6 November 2000 相似文献
5.
One of the aims of the Earth Explorer Gravity Field and Steady-State Ocean Circulation (GOCE) mission is to provide global
and regional models of the Earth's gravity field and of the geoid with high spatial resolution and accuracy. Using the GOCE
error model, simulation studies were performed in order to estimate the accuracy of datum transfer in different areas of the
Earth. The results showed that with the GOCE error model, the standard deviation of the height anomaly differences is about
one order of magnitude better than the corresponding value with the EGM96 error model. As an example, the accuracy of the
vertical datum transfer from the tide gauge of Amsterdam to New York was estimated equal to 57 cm when the EGM96 error model
was used, while in the case of GOCE error model this accuracy was increased to 6 cm. The geoid undulation difference between
the two places is about 76.5 m. Scaling the GOCE errors to the local gravity variance, the estimated accuracy varied between
3 and 7 cm, depending on the scaling model.
Received: 1 March 2000 / Accepted: 21 February 2001 相似文献
6.
Regularization of geopotential determination from satellite data by variance components 总被引:11,自引:18,他引:11
Different types of present or future satellite data have to be combined by applying appropriate weighting for the determination
of the gravity field of the Earth, for instance GPS observations for CHAMP with satellite to satellite tracking for the coming
mission GRACE as well as gradiometer measurements for GOCE. In addition, the estimate of the geopotential has to be smoothed
or regularized because of the inversion problem. It is proposed to solve these two tasks by Bayesian inference on variance
components. The estimates of the variance components are computed by a stochastic estimator of the traces of matrices connected
with the inverse of the matrix of normal equations, thus leading to a new method for determining variance components for large
linear systems. The posterior density function for the variance components, weighting factors and regularization parameters
are given in order to compute the confidence intervals for these quantities. Test computations with simulated gradiometer
observations for GOCE and satellite to satellite tracking for GRACE show the validity of the approach.
Received: 5 June 2001 / Accepted: 28 November 2001 相似文献
7.
Downward continuation and geoid determination based on band-limited airborne gravity data 总被引:4,自引:3,他引:4
The downward continuation of the harmonic disturbing gravity potential, derived at flight level from discrete observations
of airborne gravity by the spherical Hotine integral, to the geoid is discussed. The initial-boundary-value approach, based
on both the direct and inverse solution to Dirichlet's problem of potential theory, is used. Evaluation of the discretized
Fredholm integral equation of the first kind and its inverse is numerically tested using synthetic airborne gravity data.
Characteristics of the synthetic gravity data correspond to typical airborne data used for geoid determination today and in
the foreseeable future: discrete gravity observations at a mean flight height of 2 to 6 km above mean sea level with minimum
spatial resolution of 2.5 arcmin and a noise level of 1.5 mGal. Numerical results for both approaches are presented and discussed.
The direct approach can successfully be used for the downward continuation of airborne potential without any numerical instabilities
associated with the inverse approach. In addition to these two-step approaches, a one-step procedure is also discussed. This
procedure is based on a direct relationship between gravity disturbances at flight level and the disturbing gravity potential
at sea level. This procedure provided the best results in terms of accuracy, stability and numerical efficiency. As a general
result, numerically stable downward continuation of airborne gravity data can be seen as another advantage of airborne gravimetry
in the field of geoid determination.
Received: 6 June 2001 / Accepted: 3 January 2002 相似文献
8.
Adaptive filtering of GOCE-derived gravity gradients of the disturbing potential in the context of the space-wise approach 总被引:1,自引:0,他引:1
Filtering and signal processing techniques have been widely used in the processing of satellite gravity observations to reduce measurement noise and correlation errors. The parameters and types of filters used depend on the statistical and spectral properties of the signal under investigation. Filtering is usually applied in a non-real-time environment. The present work focuses on the implementation of an adaptive filtering technique to process satellite gravity gradiometry data for gravity field modeling. Adaptive filtering algorithms are commonly used in communication systems, noise and echo cancellation, and biomedical applications. Two independent studies have been performed to introduce adaptive signal processing techniques and test the performance of the least mean-squared (LMS) adaptive algorithm for filtering satellite measurements obtained by the gravity field and steady-state ocean circulation explorer (GOCE) mission. In the first study, a Monte Carlo simulation is performed in order to gain insights about the implementation of the LMS algorithm on data with spectral behavior close to that of real GOCE data. In the second study, the LMS algorithm is implemented on real GOCE data. Experiments are also performed to determine suitable filtering parameters. Only the four accurate components of the full GOCE gravity gradient tensor of the disturbing potential are used. The characteristics of the filtered gravity gradients are examined in the time and spectral domain. The obtained filtered GOCE gravity gradients show an agreement of 63–84 mEötvös (depending on the gravity gradient component), in terms of RMS error, when compared to the gravity gradients derived from the EGM2008 geopotential model. Spectral-domain analysis of the filtered gradients shows that the adaptive filters slightly suppress frequencies in the bandwidth of approximately 10–30 mHz. The limitations of the adaptive LMS algorithm are also discussed. The tested filtering algorithm can be connected to and employed in the first computational steps of the space-wise approach, where a time-wise Wiener filter is applied at the first stage of GOCE gravity gradient filtering. The results of this work can be extended to using other adaptive filtering algorithms, such as the recursive least-squares and recursive least-squares lattice filters. 相似文献
9.
One of the most basic and important tools in optimal spectral gravity field modelling is the method of Wiener filtering.
Originally developed for applications in analogue signal analysis and communication engineering, Wiener filtering has become
a standard linear estimation technique of modern operational geodesy, either as an independent practical tool for data de-noising
in the frequency domain or as an integral component of a more general signal estimation methodology (input–output systems
theory). Its theoretical framework is based on the Wiener–Kolmogorov linear prediction theory for stationary random fields
in the presence of additive external noise, and thus it is closely related to the (more familiar to geodesists) method of
least-squares collocation with random observation errors. The main drawback of Wiener filtering that makes its use in many
geodetic applications problematic stems from the stationarity assumption for both the signal and the noise involved in the
approximation problem. A modified Wiener-type linear estimation filter is introduced that can be used with noisy data obtained
from an arbitrary deterministic field under the masking of non-stationary random observation errors. In addition, the sampling
resolution of the input data is explicitly taken into account within the estimation algorithm, resulting in a resolution-dependent
optimal noise filter. This provides a more insightful approach to spectral filtering techniques for noise reduction, since
the data resolution parameter has not been directly incorporated in previous formulations of frequency-domain estimation problems
for gravity field signals with discrete noisy data.
Received: 1 November 2000 / Accepted: 19 June 2001 相似文献
10.
Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt 总被引:1,自引:0,他引:1
C. C. Tscherning Awar Radwan A. A. Tealeb S. M. Mahmoud M. Abd El-Monum Ramdan Hassan I. El-Syaed K. Saker 《Journal of Geodesy》2001,75(7-8):343-348
The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and
the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity
data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not
using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with
gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general
gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan
Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using
only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was
in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used,
as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave
difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the
agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances
and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account
for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads.
Received: 14 August 2000 / Accepted: 28 February 2001 相似文献
11.
The Somigliana–Pizzetti gravity field (the International gravity formula), namely the gravity field of the level ellipsoid
(the International Reference Ellipsoid), is derived to the sub-nanoGal accuracy level in order to fulfil the demands of modern
gravimetry (absolute gravimeters, super conducting gravimeters, atomic gravimeters). Equations (53), (54) and (59) summarise
Somigliana–Pizzetti gravity Γ(φ,u) as a function of Jacobi spheroidal latitude φ and height u to the order ?(10−10 Gal), and Γ(B,H) as a function of Gauss (surface normal) ellipsoidal latitude B and height H to the order ?(10−10 Gal) as determined by GPS (`global problem solver'). Within the test area of the state of Baden-Württemberg, Somigliana–Pizzetti
gravity disturbances of an average of 25.452 mGal were produced. Computer programs for an operational application of the new
international gravity formula with (L,B,H) or (λ,φ,u) coordinate inputs to a sub-nanoGal level of accuracy are available on the Internet.
Received: 23 June 2000 / Accepted: 2 January 2001 相似文献
12.
不同于当前广泛使用的空域法、时域法、直接解法,本文尝试采用Torus方法处理GOCE实测数据,利用71 d的GOCE卫星引力梯度数据反演了200阶次GOCE地球重力场模型,实现了对参考模型的精化。首先,采用Butterworth零相移滤波方法加移去—恢复技术,处理引力梯度观测值中的有色噪声,并利用泰勒级数展开和Kriging方法对GOCE卫星引力梯度数据进行归算和格网化,计算得到了名义轨道上格网点处的引力梯度数据。然后,利用2D-FFT技术和块对角最小二乘方法处理名义轨道上数据,获得了200阶次的GOCE地球重力场模型GOCE_Torus。利用中国和美国的GPS/水准数据进行外部检核结果说明,GOCE_Torus与ESA发布的同期模型的精度相当;GOCE_Torus模型与200阶次的EGM2008模型相比,在美国区域精度相当,但在中国区域精度提高了4.6 cm,这充分体现了GOCE卫星观测数据对地面重力稀疏区的贡献。Torus方法拥有快速高精度反演卫星重力场模型的优势,可以在重力梯度卫星的设计、误差分析及在轨快速评估等方面得到充分应用。 相似文献
13.
Evaluation of the first GOCE static gravity field models using terrestrial gravity,vertical deflections and EGM2008 quasigeoid heights 总被引:1,自引:1,他引:0
Recently, four global geopotential models (GGMs) were computed and released based on the first 2 months of data collected
by the Gravity field and steady-state Ocean Circulation Explorer (GOCE) dedicated satellite gravity field mission. Given that
GOCE is a technologically complex mission and different processing strategies were applied to real space-collected GOCE data
for the first time, evaluation of the new models is an important aspect. As a first assessment strategy, we use terrestrial
gravity data over Switzerland and Australia and astrogeodetic vertical deflections over Europe and Australia as ground-truth
data sets for GOCE model evaluation. We apply a spectral enhancement method (SEM) to the truncated GOCE GGMs to make their
spectral content more comparable with the terrestrial data. The SEM utilises the high-degree bands of EGM2008 and residual
terrain model data as a data source to widely bridge the spectral gap between the satellite and terrestrial data. Analysis
of root mean square (RMS) errors is carried out as a function of (i) the GOCE GGM expansion degree and (ii) the four different
GOCE GGMs. The RMS curves are also compared against those from EGM2008 and GRACE-based GGMs. As a second assessment strategy,
we compare global grids of GOCE GGM and EGM2008 quasigeoid heights. In connection with EGM2008 error estimates, this allows
location of regions where GOCE is likely to deliver improved knowledge on the Earth’s gravity field. Our ground truth data
sets, together with the EGM2008 quasigeoid comparisons, signal clear improvements in the spectral band ~160–165 to ~180–185
in terms of spherical harmonic degrees for the GOCE-based GGMs, fairly independently of the individual GOCE model used. The
results from both assessments together provide strong evidence that the first 2 months of GOCE observations improve the knowledge
of the Earth’s static gravity field at spatial scales between ~125 and ~110 km, particularly over parts of Asia, Africa, South
America and Antarctica, in comparison with the pre-GOCE-era. 相似文献
14.
Johannes Bouman 《Journal of Geodesy》2012,86(4):287-304
The vertical gradients of gravity anomaly and gravity disturbance can be related to horizontal first derivatives of deflection
of the vertical or second derivatives of geoidal undulations. These are simplified relations of which different variations
have found application in satellite altimetry with the implicit assumption that the neglected terms—using remove-restore—are
sufficiently small. In this paper, the different simplified relations are rigorously connected and the neglected terms are
made explicit. The main neglected terms are a curvilinear term that accounts for the difference between second derivatives
in a Cartesian system and on a spherical surface, and a small circle term that stems from the difference between second derivatives
on a great and small circle. The neglected terms were compared with the dynamic ocean topography (DOT) and the requirements
on the GOCE gravity gradients. In addition, the signal root-mean-square (RMS) of the neglected terms and vertical gravity
gradient were compared, and the effect of a remove-restore procedure was studied. These analyses show that both neglected
terms have the same order of magnitude as the DOT gradient signal and may be above the GOCE requirements, and should be accounted
for when combining altimetry derived and GOCE measured gradients. The signal RMS of both neglected terms is in general small
when compared with the signal RMS of the vertical gravity gradient, but they may introduce gradient errors above the spherical
approximation error. Remove-restore with gravity field models reduces the errors in the vertical gravity gradient, but it
appears that errors above the spherical approximation error cannot be avoided at individual locations. When computing the
vertical gradient of gravity anomaly from satellite altimeter data using deflections of the vertical, the small circle term
is readily available and can be included. The direct computation of the vertical gradient of gravity disturbance from satellite
altimeter data is more difficult than the computation of the vertical gradient of gravity anomaly because in the former case
the curvilinear term is needed, which is not readily available. 相似文献
15.
R. Pail 《Journal of Geodesy》2005,79(4-5):231-241
In the recent design of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, the gravity gradients are defined in the gradiometer reference frame (GRF), which deviates from the actual flight direction (local orbit reference frame, LORF) by up to 3–4°. The main objective of this paper is to investigate the effect of uncertainties in the knowledge of the gradiometer orientation due to attitude reconstitution errors on the gravity field solution. In the framework of several numerical simulations, which are based on a realistic mission configuration, different scenarios are investigated, to provide the accuracy requirements of the orientation information. It turns out that orientation errors have to be seriously considered, because they may represent a significant error component of the gravity field solution. While in a realistic mission scenario (colored gradiometer noise) the gravity field solutions are quite insensitive to small orientation biases, random noise applied to the attitude information can have a considerable impact on the accuracy of the resolved gravity field models. 相似文献
16.
Biao Lu Zhicai Luo Bo Zhong Hao Zhou Frank Flechtner Christoph Förste Franz Barthelmes Rui Zhou 《Journal of Geodesy》2018,92(5):561-572
Based on tensor theory, three invariants of the gravitational gradient tensor (IGGT) are independent of the gradiometer reference frame (GRF). Compared to traditional methods for calculation of gravity field models based on the gravity field and steady-state ocean circulation explorer (GOCE) data, which are affected by errors in the attitude indicator, using IGGT and least squares method avoids the problem of inaccurate rotation matrices. The IGGT approach as studied in this paper is a quadratic function of the gravity field model’s spherical harmonic coefficients. The linearized observation equations for the least squares method are obtained using a Taylor expansion, and the weighting equation is derived using the law of error propagation. We also investigate the linearization errors using existing gravity field models and find that this error can be ignored since the used a-priori model EIGEN-5C is sufficiently accurate. One problem when using this approach is that it needs all six independent gravitational gradients (GGs), but the components \(V_{xy}\) and \(V_{yz}\) of GOCE are worse due to the non-sensitive axes of the GOCE gradiometer. Therefore, we use synthetic GGs for both inaccurate gravitational gradient components derived from the a-priori gravity field model EIGEN-5C. Another problem is that the GOCE GGs are measured in a band-limited manner. Therefore, a forward and backward finite impulse response band-pass filter is applied to the data, which can also eliminate filter caused phase change. The spherical cap regularization approach (SCRA) and the Kaula rule are then applied to solve the polar gap problem caused by GOCE’s inclination of \(96.7^{\circ }\). With the techniques described above, a degree/order 240 gravity field model called IGGT_R1 is computed. Since the synthetic components of \(V_{xy}\) and \(V_{yz}\) are not band-pass filtered, the signals outside the measurement bandwidth are replaced by the a-priori model EIGEN-5C. Therefore, this model is practically a combined gravity field model which contains GOCE GGs signals and long wavelength signals from the a-priori model EIGEN-5C. Finally, IGGT_R1’s accuracy is evaluated by comparison with other gravity field models in terms of difference degree amplitudes, the geostrophic velocity in the Agulhas current area, gravity anomaly differences as well as by comparison to GNSS/leveling data. 相似文献
17.
New solutions for the geodetic coordinate transformation 总被引:5,自引:2,他引:5
G. C. Jones 《Journal of Geodesy》2002,76(8):437-446
The Cartesian-to-geodetic-coordinate transformation is approached from a new perspective. Existence and uniqueness of geodetic
representation are presented, along with a clear geometric picture of the problem and the role of the ellipse evolute. A new
solution is found with a Newton-method iteration in the reduced latitude; this solution is proved to work for all points in
space. Care is given to error propagation when calculating the geodetic latitude and height.
Received: 9 August 2001 / Accepted: 27 March 2002
Acknowledgments. The author would like to thank the Clifford W.␣Tompson scholarship fund, Dr. Brian DeFacio, the University of Missouri College
of Arts &Sciences, and the United States Air Force. He also thanks a reviewer for suggesting and providing a prototype MATLAB
code. A MATLAB program for the iterative sequence is presented at the end of the paper (Appendix A). 相似文献
18.
A potential-type Molodensky telluroid based upon a minimum-distance mapping is derived. With respect to a reference potential
of Somigliana–Pizzetti type which relates to the World Geodetic Datum 2000, it is shown that a point-wise minimum-distance
mapping of the topographical surface of the Earth onto the telluroid surface, constrained to the gauge W(P)=u(p), leads to a system of four nonlinear normal equations. These normal equations are solved by a fast Newton–Raphson iteration.
Received: 7 February 2000 / Accepted: 23 October 2001 相似文献
19.
P. Moore 《Journal of Geodesy》2001,75(5-6):241-254
Dual satellite crossovers (DXO) between the two European Remote Sensing satellites ERS-1 and ERS-2 and TOPEX/Poseidon are
used to (1) refine the Earth's gravity field and (2) extend the study of the ERS-2 altimetric range stability to cover the
first four years of its operation. The enhanced gravity field model, AGM-98, is validated by several methodologies and will
be shown to provide, in particular, low geographically correlated orbital error for ERS-2. For the ERS-2 altimetric range
study, TOPEX/Poseidon is first calibrated through comparison against in situ tide gauge data. A time series of the ERS-2 altimeter
bias has been recovered along with other geophysical correction terms using tables for bias jumps in the range measurements
at the single point target response (SPTR) events. On utilising the original version of the SPTR tables the overall bias drift
is seen to be 2.6±1.0 mm/yr with an RMS of fit of 12.2 mm but with discontinuities at the centimetre level at the SPTR events.
On utilising the recently released revised tables, SPTR2000, the drift is better defined at 2.4±0.6 mm/yr with the RMS of
fit reduced to 3.7 mm. Investigations identify the sea-state bias as a source of error with corrections affecting the overall
drift by close to 1.2 mm/yr.
Received: 25 May 2000 / Accepted: 24 January 2001 相似文献
20.
Regional height systems do not refer to a common equipotential surface, such as the geoid. They are usually referred to the mean sea level at a reference tide gauge. As mean sea level varies (by ±1 to 2 m) from place to place and from continent to continent each tide gauge has an unknown bias with respect to a common reference surface, whose determination is what the height datum problem is concerned with. This paper deals with this problem, in connection to the availability of satellite gravity missions data. Since biased heights enter into the computation of terrestrial gravity anomalies, which in turn are used for geoid determination, the biases enter as secondary or indirect effect also in such a geoid model. In contrast to terrestrial gravity anomalies, gravity and geoid models derived from satellite gravity missions, and in particular GRACE and GOCE, do not suffer from those inconsistencies. Those models can be regarded as unbiased. After a review of the mathematical formulation of the problem, the paper examines two alternative approaches to its solution. The first one compares the gravity potential coefficients in the range of degrees from 100 to 200 of an unbiased gravity field from GOCE with those of the combined model EGM2008, that in this range is affected by the height biases. This first proposal yields a solution too inaccurate to be useful. The second approach compares height anomalies derived from GNSS ellipsoidal heights and biased normal heights, with anomalies derived from an anomalous potential which combines a satellite-only model up to degree 200 and a high-resolution global model above 200. The point is to show that in this last combination the indirect effects of the height biases are negligible. To this aim, an error budget analysis is performed. The biases of the high frequency part are proved to be irrelevant, so that an accuracy of 5 cm per individual GNSS station is found. This seems to be a promising practical method to solve the problem. 相似文献