共查询到20条相似文献,搜索用时 15 毫秒
1.
A technique for the analysis of low–low intersatellite range-rate data in a gravity mapping mission is explored. The technique
is based on standard tracking data analysis for orbit determination but uses a spherical coordinate representation of the
12 epoch state parameters describing the baseline between the two satellites. This representation of the state parameters
is exploited to allow the intersatellite range-rate analysis to benefit from information provided by other tracking data types
without large simultaneous multiple-data-type solutions. The technique appears especially valuable for estimating gravity
from short arcs (e.g. less than 15 minutes) of data. Gravity recovery simulations which use short arcs are compared with those
using arcs a day in length. For a high-inclination orbit, the short-arc analysis recovers low-order gravity coefficients remarkably
well, although higher-order terms, especially sectorial terms, are less accurate. Simulations suggest that either long or
short arcs of the Gravity Recovery and Climate Experiment (GRACE) data are likely to improve parts of the geopotential spectrum
by orders of magnitude.
Received: 26 June 2001 / Accepted: 21 January 2002 相似文献
2.
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 相似文献
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.
Low-low satellite-to-satellite tracking: a comparison between analytical linear orbit perturbation theory and numerical integration 总被引:1,自引:0,他引:1
P.N.A.M. Visser 《Journal of Geodesy》2005,79(1-3):160-166
Low-low satellite-to-satellite tracking (ll-SST) range-rate observations have been predicted by two methods: one based on a linear perturbation theory in combination with the Hill equations, and one based on solving the equations of motion of two low-flying satellites by numerical integration. The two methods produce almost equivalent Fourier spectra of the range-rate observations after properly taking into account a few resonant terms. For a typical GRACE-type configuration, where the two satellites trail each other at a distance of 300 km at an altitude of 460 km, and in the presence of the EGM96 gravity field model, complete to spherical harmonic degree and order 70, the agreement between the Fourier spectra is about 1 mm/s compared to a root-mean-square (RMS) value of more than 220 mm/s for the range-rate signal. The discrepancy of 1 mm/s can be reduced significantly when not taking into account perturbations caused by the J2 term. Excluding the J2 term, the agreement between the two methods improves to 0.4 mm/s compared to a RMS value of 6 mm/s for the range-rate signal. These values are 0.01 and 2.3 mm/s when ignoring the spectrum for frequencies below two cycles per orbital revolution, reducing the discrepancy even further to about 0.5% of the signal. The selected linear perturbation theory is thus capable of modeling gravity field induced range-rate observations with very high precision for a large part of the spectrum. 相似文献
5.
P. Schwintzer C. Reigber A. Bode Z. Kang S. Y. Zhu F.-H. Massmann J. C. Raimondo R. Biancale G. Balmino J. M. Lemoine B. Moynot J. C. Marty F. Barlier Y. Boudon 《Journal of Geodesy》1997,71(4):189-208
Summary. GFZ Potsdam and GRGS Toulouse/Grasse jointly developed a new pair of global models of the Earth's gravity field to satisfy
the requirements of the recent and future geodetic and altimeter satellite missions. A precise gravity model is a prerequisite
for precise satellite orbit restitution, tracking station positioning and altimeter data reduction. According to different
applications envisaged, the new model exists in two parallel versions: the first one being derived exclusively from satellite
tracking data acquired on 34 satellites, the second one further incorporating satellite altimeter data over the oceans and
terrestrial gravity data. The most recent “satellite-only” gravity model is labelled GRIM4-S4 and the “combined” gravity model
GRIM4-C4. The models are solutions in spherical harmonics and have a resolution up to degree and order 60 plus a few resonance
terms in the case of GRIM4-S4, and up to degree/order 72 in the case of GRIM4-C4, corresponding to a spatial resolution of
555 km at the Earth's surface. The gravitational coefficients were estimated in a rigorous least squares adjustment simultaneously
with ocean tidal terms and tracking station position parameters, so that each gravity model is associated with a consistent
ocean tide model and a terrestrial reference frame built up by over 300 optical, laser and Doppler tracking stations. Comprehensive
quality tests with external data and models, and test arc computations over a wide range of satellites have demonstrated the
state-of-the-art capabilities of both solutions in long-wavelength geoid representation and in precise orbit computation.
Received 1 February 1996; Accepted 17 July 1996 相似文献
6.
In October 1998 the IGEX field campaign, the first coordinated international effort to monitor GLONASS satellites on global
basis, was started. Currently about 40 institutions worldwide support this effort either by providing GLONASS tracking data
or in operating related data and analysis centers. The increasing quality and consistency of the calculated GLONASS orbits
(about 25 cm early in 2000), even after the end of the official IGEX field campaign, are shown. Particular attention is drawn
to the combination of precise ephemerides in order to generate a robust, reliable and complete IGEX orbits product. Some problems
in modeling the effect of solar radiation pressure on GLONASS satellites are demonstrated. Finally, the expected benefits
and prospects of the upcoming International GLOnass Service-Pilot Project (IGLOS-PP) of the International GPS Service (IGS)
are discussed in more detail.
Received: 17 August 2000 / Accepted: 12 April 2001 相似文献
7.
Global mean sea surface heights (SSHs) and gravity anomalies on a 2′×2′ grid were determined from Seasat, Geosat (Exact Repeat Mission and Geodetic Mission), ERS-1 (1.5-year mean of 35-day, and
GM), TOPEX/POSEIDON (T/P) (5.6-year mean) and ERS-2 (2-year mean) altimeter data over the region 0∘–360∘ longitude and –80∘–80∘ latitude. To reduce ocean variabilities and data noises, SSHs from non-repeat missions were filtered by Gaussian filters
of various wavelengths. A Levitus oceanic dynamic topography was subtracted from the altimeter-derived SSHs, and the resulting
heights were used to compute along-track deflection of the vertical (DOV). Geoidal heights and gravity anomalies were then
computed from DOV using the deflection-geoid and inverse Vening Meinesz formulae. The Levitus oceanic dynamic topography was
added back to the geoidal heights to obtain a preliminary sea surface grid. The difference between the T/P mean sea surface
and the preliminary sea surface was computed on a grid by a minimum curvature method and then was added to the preliminary
grid. The comparison of the NCTU01 mean sea surface height (MSSH) with the T/P and the ERS-1 MSSH result in overall root-mean-square
(RMS) differences of 5.0 and 3.1 cm in SSH, respectively, and 7.1 and 3.2 μrad in SSH gradient, respectively. The RMS differences
between the predicted and shipborne gravity anomalies range from 3.0 to 13.4 mGal in 12 areas of the world's oceans.
Received: 26 September 2001 / Accepted: 3 April 2002
Correspondence to: C. Hwang
Acknowledgements. This research is partly supported by the National Science Council of ROC, under grants NSC89-2611-M-009-003-OP2 and NSC89-2211-E-009-095.
This is a contribution to the IAG Special Study Group 3.186. The Geosat and ERS1/2 data are from NOAA and CERSAT/France, respectively.
The T/P data were provided by AVISO. The CLS and GSFC00 MSS models were kindly provided by NASA/GSFC and CLS, respectively.
Drs. Levitus, Monterey, and Boyer are thanked for providing the SST model. Dr. T. Gruber and two anonymous reviewers provided
very detailed reviews that improved the quality of this paper. 相似文献
8.
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 相似文献
9.
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 相似文献
10.
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 相似文献
11.
Since its launch in 2002, the Gravity Recovery and Climate Experiment (GRACE) mission has been providing measurements of the
time-varying Earth gravity field. The GRACE mission architecture includes two satellites in near-circular, near-polar orbits
separated in the along-track direction by approximately 220 km (e.g. collinear). A microwave ranging instrument measures changes
in the distance between the spacecraft, while accelerometers on each spacecraft are used to measure changes in distance due
to non-gravitational forces. The fact that the satellites are in near-polar orbits coupled with the fact that the inter-satellite
range measurements are directed in the along-track direction, contributes to longitudinal striping in the estimated gravity
fields. This paper examines four candidate mission architectures for a future gravity recovery satellite mission to assess
their potential in measuring the gravity field more accurately than GRACE. All satellites were assumed to have an improved
measurement system, with an inter-satellite laser ranging instrument and a drag-free system for removal of non-gravitational
accelerations. Four formations were studied: a two-satellite collinear pair similar to GRACE; a four-satellite architecture
with two collinear pairs; a two-satellite cartwheel formation; and a four-satellite cartwheel formation. A cartwheel formation
consists of satellites performing in-plane, relative elliptical motion about their geometric center, so that inter-satellite
measurements are, at times, directed radially (e.g. parallel to the direction towards the center of the Earth) rather than
along-track. Radial measurements, unlike along-track measurements, have equal sensitivity to mass distribution in all directions
along the Earth’s surface and can lead to higher spatial resolution in the derived gravity field. The ability of each architecture
to recover the gravity field was evaluated using numerical simulations performed with JPL’s GIPSY-OASIS software package.
Thirty days of data were used to estimate gravity fields complete to degree and order 60. Evaluations were done for 250 and
400 km nominal orbit altitudes. The sensitivity of the recovered gravity field to under-sampled effects was assessed using
simulated errors in atmospheric/ocean dealiasing (AOD) models. Results showed the gravity field errors associated with the
four-satellite cartwheel formation were approximately one order of magnitude lower than the collinear satellite pair when
only measurement system errors were included. When short-period AOD model errors were introduced, the gravity field errors
for each formation were approximately the same. The cartwheel formations eliminated most of the longitudinal striping seen
in the gravity field errors. A covariance analysis showed the error spectrum of the cartwheel formations to be lower and more
isotropic than that of the collinear formations. 相似文献
12.
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 相似文献
13.
GPS measurements of ocean loading and its impact on zenith tropospheric delay estimates: a case study in Brittany, France 总被引:1,自引:0,他引:1
S. Vey E. Calais M. Llubes N. Florsch G. Woppelmann J. Hinderer M. Amalvict M. F. Lalancette B. Simon F. Duquenne J. S. Haase 《Journal of Geodesy》2002,76(8):419-427
The results from a global positioning system (GPS) experiment carried out in Brittany, France, in October 1999, aimed at
measuring crustal displacements caused by ocean loading and quantifying their effects on GPS-derived tropospheric delay estimates,
are presented. The loading effect in the vertical and horizontal position time series is identified, however with significant
disagreement in amplitude compared to ocean loading model predictions. It is shown that these amplitude misfits result from
spatial tropospheric heterogeneities not accounted for in the data processing. The effect of ocean loading on GPS-derived
zenith total delay (ZTD) estimates is investigated and a scaling factor of 4.4 between ZTD and station height for a 10° elevation
cut-off angle is found (i.e. a 4.4-cm station height error would map into a 1-cm ZTD error). Consequently, unmodeled ocean
loading effects map into significant errors in ZTD estimates and ocean loading modeling must be properly implemented when
estimating ZTD parameters from GPS data for meteorological applications. Ocean loading effects must be known with an accuracy
of better than 3 cm in order to meet the accuracy requirements of meteorological and climatological applications of GPS-derived
precipitable water vapor.
Received: 16 July 2001 / Accepted: 25 April 2002
Acknowledgments. The authors are grateful to H.G. Scherneck for fruitful discussions and for his help with the ocean loading calculations.
They thank H. Vedel for making the HIRLAM data available; D. Jerett for helpful discussions; and the city of Rostrenen, the
Laboratoire d'Océanographie of Concarneau, and the Institut de Protection et de S?reté Nucléaire (BERSSIN) for their support
during the GPS measurement campaign. Reviews by C.K. Shum and two anonymous referees significantly improved this paper. This
work was carried out in the framework of the MAGIC project (http://www.acri.fr/magic), funded by the European Commission,
Environment and Climate Program (EC Contract ENV4-CT98–0745).
Correspondence to: E. Calais, Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, IN 47907-1397, USA. e-mail:
ecalais@purdue.edu Tel. : +1-765-496-2915; Fax:+1-765-496-1210 相似文献
14.
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 相似文献
15.
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 相似文献
16.
Starting from the analytical theory of perturbed␣circular motions presented in Celestial Mechanics (Bois 1994) and from specific
extended formulations of the perturbations in a uniformly rotating plane of constant inclination, this paper presents an extended
formulation of the solution. The actual gain made through this extension is the establishment of a first-order predictive
theory written in spherical coordinates and thus free of singularities, whose perturbations are directly expressed in the
local orbital frame generally used in satellite geodesy. This new formulation improves the generality, the precision and the
field of applications of the theory. It is particularly devoted to the analysis of satellite position perturbations for satellites
in low eccentricity orbits usually used for many Earth observation applications. An application to the TOPEX/Poseidon (T/P)
orbit is performed. In particular, contour maps are provided which show the geographical location of orbit differences coming
from geopotential coefficient differences of two recent gravity field models. Comparison of predicted radial and along-track
orbit differences with respect to numerical results provided by the French group (CNES, in Toulouse) in charge of the T/P
orbit are convincing.
Received 22 January 1996; Accepted 19 September 1996 相似文献
17.
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 相似文献
18.
卫星自主定轨是提高全球卫星导航系统(GNSS)可靠性、稳健性、完整性和生存能力的重要保证。新一代的北斗卫星已可以进行星间链路测距,从而达到提高卫星全球跟踪能力以及实现整个卫星导航系统的自主定轨。然而由于卫星运行会受到多种摄动力的影响,如果不能对这些摄动力进行精密的改正,在没有地面或其他天体提供绝对约束的条件下,导航系统会随着自主定轨时间的延长出现星座整体旋转。卫星所受摄动力分为保守力和非保守力两部分:对于保守力,如地球非球形摄动、潮汐摄动、太阳月球和其他三体引力,现在已有的力学模型可以很精确地进行改正;而非保守力(如太阳光压摄动),则难以用精确的模型进行改正,因此成为影响卫星定轨精度的主要因素。星载加速度计可以高精度地测量非保守力,并已成功应用于重力卫星(CHAMP、GRACE、GOCE)的重力场反演与大气研究中。本文研究主要探讨采用星上加速度计提高北斗卫星自主定轨精度和延长自主定轨时长的可行性。利用模拟的卫星轨道和星间链路数据,以及现有的星载加速度计误差模型,对北斗卫星系统分别使用星间链路数据和星间链路与加速度计组合数据,进行自主定轨与精度评定。计算结果表明,使用星间链路与星载加速度计数据进行自主定轨,较单纯使用星间链路数据精度具有明显改进。在模拟的星间测距观测数据具有0.33m随机噪声以及分米级系统误差,自主定轨两个月的情况下,联合使用加速度计数据的自主定轨IGSO和MEO卫星精度为分米级,而仅使用星间链路数据的定轨精度约为3~6m,比使用加速度计精度低一个量级。 相似文献
19.
G. Ramillien 《Journal of Geodesy》2002,76(3):139-149
A fast spherical harmonic approach enables the computation of gravitational or magnetic potential created by a non-uniform
shell of material bounded by uneven topographies. The resulting field can be evaluated outside or inside the sphere, assuming
that density of the shell varies with latitude, longitude, and radial distance. To simplify, the density (or magnetization)
source inside the sphere is assumed to be the product of a surface function and a power series expansion of the radial distance.
This formalism is applied to compute the gravity signal of a steady, dry atmosphere. It provides geoid/gravity maps at sea
level as well as satellite altitude. Results of this application agree closely with those of earlier studies, where the atmosphere
contribution to the Earth's gravity field was determined using more time-consuming methods.
Received: 14 August 2000 / Accepted: 19 March 2001 相似文献
20.
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 相似文献