共查询到20条相似文献,搜索用时 109 毫秒
1.
Lars E. Sj?berg 《Journal of Geodesy》1988,62(2):93-101
The spherical harmonic coefficients of the Earth’s gravitational potential are conveniently determined by integration of gravity
data or potential data (derived from satellite altimetry) over a sphere. The major problem of such a method is that the data,
given on the non-spherical surface of the Earth, must be reduced to the sphere.
A new integral formula over the surface of the Earth is derived. With this formula improved first order topographic corrections
to the spherical formulas are obtained. 相似文献
2.
The solution of the linear Molodensky problem by analytical continuation to point level is numerically the most convenient
of all the theoretically equivalent solutions. It is obtained by successively applying the same integral operator and it does
not depend explicitly on the terrain inclination. However, its dependence on the computation point restricts somehow the computational
efficiency. The use of the Fourier transform for the evaluation of the integral operator in planar approximation lessens significantly
the burden of computations. Using this spectral approach, the problem has been reformulated and solved in the frequency domain.
Moreover, it is shown that the solution can be easily split into two steps: (a) “downward” continuation to sea level, which
is independent of the computation point, and (b) “upward” continuation from sea to point level, using the values computed
at sea level. Such a treatment not only simplifies the formulas and increases the numerical efficiency but also clarifies
the physical interpretation and the theoretical equivalence of the continuation solution with respect to the other solution
types. Numerical tests have been performed to investigate which terms in the Molodensky series are of significance for geoid
and deflection computations. The practical difficulty of differences in the grid spacings of gravity and height data has been
overcome by frequency domain interpolation.
Presented at theXIX IUGG General Assembly, Vancouver, B.C., August 9–22, 1987. 相似文献
3.
The classical integral formula for determining the indirect effect in connection with the Stokes–Helmert method is related
to a planar approximation of the sea level. A strict integral formula, as well as some approximations to it, are derived.
It is concluded that the cap- size truncated integral formulas will suffer from the omission of some long-wavelength contributions,
of the order of 50 cm in high mountains for the classical formula. This long-wavelength information can be represented by
a set of spherical harmonic coefficients of the topography to, say, degree and order 360. Hence, for practical use, a combination
of the classical formula and a set of spherical harmonics is recommended.
Received: 10 March 1998 / Accepted: 16 November 1998 相似文献
4.
Laser-based validation of GLONASS orbits by short-arc technique 总被引:1,自引:0,他引:1
F. Barlier C. Berger P. Bonnefond P. Exertier O. Laurain J. F. Mangin J. M. Torre 《Journal of Geodesy》2001,75(11):600-612
The International GLONASS Experiment (IGEX-98) was carried out between 19 October 1998 and 19 April 1999. Among several objectives
was the precise orbit determination of GPS and GLONASS satellites and its validation by laser ranging observations. Local
laser-based orbit corrections (radial, tangential and normal components in a rotating orbital local reference frame) are computed
using a geometrical short-arc technique. The order of magnitude of these corrections is at the level of few decimeters, depending
on the considered components. The orbit corrections are analyzed as a function of several parameters (date, orbital plane,
geographical area). The mean corrections are at the level of several centimeters. However, when averaging over the entire
campaign and for all the satellites, no mean radial, tangential and normal orbit corrections are found. The origin of the
observed corrections is considered (errors due to the geocentric gravitational constant, the non-gravitational forces, the
thermal equilibrium of on-board equipment, the reference systems, the location and the signature of the retroreflector array,
and the precision of the satellite laser ranges). Some features are also due to errors in the radio-tracking GLONASS orbits.
Further investigations will be needed to better understand the origin of various biases.
Received: 17 February 2000 / Accepted: 31 January 2001 相似文献
5.
Y. M. Wang 《Journal of Geodesy》1989,63(4):359-370
The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of
the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies
have to be analytically downward continued from the earth's surface to a sphere—at least to the ellipsoid. The goal of this
paper is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models.
The basic method for the downward continuation is the gradient solution (theg
1 term). The terrain correction has also been computed because of the role it can play as a correction term when calculating
harmonic coefficients from surface gravity data.
Theg
1 term and the terrain correction were expanded into the spherical harmonics up to180
th
order. The corrections (theg
1 term and the terrain correction) have the order of about 2% of theRMS value of degree variance of the disturbing potential per degree. The influences of theg
1 term and the terrain correction on the geoid take the order of 1 meter (RMS value of corrections of the geoid undulation) and on the deflections of the vertical is of the order 0.1″ (RMS value of correction of the deflections of the vertical). 相似文献
6.
An inverse Poisson integral technique has been used to determine a gravity field on the geoid which, when continued by analytic
free space methods to the topographic surface, agrees with the observed field. The computation is performed in three stages,
each stage refining the previous solution using data at progressively increasing resolution (1o×1o, 5′×5′, 5/8′×5/8′) from a decreasing area of integration. Reduction corrections are computed at 5/8′×5/8′ granularity by
differencing the geoidal and surface values, smoothed by low-pass filtering and sub-sampled at 5′ intervals. This paper discusses
1o×1o averages of the reduction corrections thus obtained for 172 1o×1o squares in western North America.
The 1o×1o mean reduction corrections are predominantly positive, varying from −3 to +15mgal, with values in excess of 5mgal for 26 squares. Their mean andrms values are +2.4 and 3.6mgal respectively and they correlate well with the mean terrain corrections as predicted byPellinen in 1962. The mean andrms contributions from the three stages of computation are: 1o×1o stage +0.15 and 0.7mgal; 5′×5′ stage +1.0 and 1.6mgal; and 5/8′×5/8′ stage +1.3 and 1.8mgal. These results reflect a tendency for the contributions to become larger and more systematically positive as the wavelengths
involved become shorter. The results are discussed in terms of two mechanisms; the first is a tendency for the absolute values
of both positive and negative anomalies to become larger when continued downwards and, the second, a non-linear rectification,
due to the correlation between gravity anomaly and topographic height, which results in the values continued to a level surface
being systematically more positive than those on the topography. 相似文献
7.
Prediction of surface horizontal displacements, and gravity and tilt changes caused by filling the Three Gorges Reservoir 总被引:11,自引:0,他引:11
Horizontal displacements, and gravity and tilt changes induced by filling the Three Gorges Reservoir are modeled using elastic
loading Green functions. When the water surface reaches its highest level, the effects become maximum on the reservoir banks.
The longitudinal and latitudinal components of the horizontal displacements reach −8.2 and 7.7 mm respectively, gravity is
increased by up to 3.4 mGal, and the prime vertical and meridian components of the tilt changes are −7.8 and −17.5 arcseconds
respectively. Accordingly, the filling of the reservoir will influence values observed from global positioning system (GPS),
gravimetry and tilt measurements in the area. The results given can be used to provide important corrections for extracting
earthquake-related signals from observed data.
Received: 19 January 2001 / Accepted: 3 September 2001 相似文献
8.
Four different implementations of Stokes' formula are employed for the estimation of geoid heights over Sweden: the Vincent
and Marsh (1974) model with the high-degree reference gravity field but no kernel modifications; modified Wong and Gore (1969)
and Molodenskii et al. (1962) models, which use a high-degree reference gravity field and modification of Stokes' kernel;
and a least-squares (LS) spectral weighting proposed by Sj?berg (1991). Classical topographic correction formulae are improved
to consider long-wavelength contributions. The effect of a Bouguer shell is also included in the formulae, which is neglected
in classical formulae due to planar approximation. The gravimetric geoid is compared with global positioning system (GPS)-levelling-derived
geoid heights at 23 Swedish Permanent GPS Network SWEPOS stations distributed over Sweden. The LS method is in best agreement,
with a 10.1-cm mean and ±5.5-cm standard deviation in the differences between gravimetric and GPS geoid heights. The gravimetric
geoid was also fitted to the GPS-levelling-derived geoid using a four-parameter transformation model. The results after fitting
also show the best consistency for the LS method, with the standard deviation of differences reduced to ±1.1 cm. For comparison,
the NKG96 geoid yields a 17-cm mean and ±8-cm standard deviation of agreement with the same SWEPOS stations. After four-parameter
fitting to the GPS stations, the standard deviation reduces to ±6.1 cm for the NKG96 geoid. It is concluded that the new corrections
in this study improve the accuracy of the geoid. The final geoid heights range from 17.22 to 43.62 m with a mean value of
29.01 m. The standard errors of the computed geoid heights, through a simple error propagation of standard errors of mean
anomalies, are also computed. They range from ±7.02 to ±13.05 cm. The global root-mean-square error of the LS model is the
other estimation of the accuracy of the final geoid, and is computed to be ±28.6 cm.
Received: 15 September 1999 / Accepted: 6 November 2000 相似文献
9.
Evaluation of deflections of the vertical on the sphere and the plane: a comparison of FFT techniques 总被引:1,自引:0,他引:1
This paper presents a set of efficient formulas to evaluate the deflections of the vertical on the sphere using gridded data.
The Vening-Meinesz formula, the topographic indirect effect on the deflections of the vertical as well as the terrain corrections
are expressed as both 2D and 1D convolutions on the sphere, and consequently can be evaluated by the 2D and the 1D fast Fourier
transform (FFT). When compared with the results obtained from pointwise integration, the use of the 1D FFT gives identical
results, and therefore these results were used as control values in this paper. The use of the spherical 2D FFT improves significantly
the computational efficiency with little sacrifice of accuracy (0.6″ rms difference from the 1D FFT results). The planar 2D FFT, which is as efficient as the spherical 2D FFT, gives worse results
(1.2″ rms difference from the 1D FFT results) because of the extra approximations.
Received: 27 February 1996; Accepted: 24 January 1997 相似文献
10.
Lars E. Sjöberg 《Journal of Geodesy》1989,63(2):213-221
Integral formulas are derived for the determination of geopotential coefficients from gravity anomalies and gravity disturbances
over the surface of the Earth. First order topographic corrections to spherical formulas are presented. In addition new integral
formulas are derived for the determination of the external gravity field from surface gravity.
Taking advantage of modern satellite positioning techniques, it is suggested that, in general, the external gravity field
as well as individual coefficients are better determined from gravity disturbances than from gravity anomalies. 相似文献
11.
A 2×2 arc-minute resolution geoid model, CARIB97, has been computed covering the Caribbean Sea. The geoid undulations refer
to the GRS-80 ellipsoid, centered at the ITRF94 (1996.0) origin. The geoid level is defined by adopting the gravity potential
on the geoid as W
0=62 636 856.88 m2/s2 and a gravity-mass constant of GM=3.986 004 418×1014 m3/s2. The geoid model was computed by applying high-frequency corrections to the Earth Gravity Model 1996 global geopotential
model in a remove-compute-restore procedure. The permanent tide system of CARIB97 is non-tidal. Comparison of CARIB97 geoid
heights to 31 GPS/tidal (ITRF94/local) benchmarks shows an average offset (h–H–N) of 51 cm, with an Root Mean Square (RMS) of 62 cm about the average. This represents an improvement over the use of a global
geoid model for the region. However, because the measured orthometric heights (H) refer to many differing tidal datums, these comparisons are biased by localized permanent ocean dynamic topography (PODT).
Therefore, we interpret the 51 cm as partially an estimate of the average PODT in the vicinity of the 31 island benchmarks.
On an island-by-island basis, CARIB97 now offers the ability to analyze local datum problems which were previously unrecognized
due to a lack of high-resolution geoid information in the area.
Received: 2 January 1998 / Accepted: 18 August 1998 相似文献
12.
Gravity recovery using COSMIC GPS data: application of orbital perturbation theory 总被引:14,自引:0,他引:14
C. Hwang 《Journal of Geodesy》2001,75(2-3):117-136
COSMIC is a joint Taiwan–US mission to study the atmosphere using the Global Positioning System (GPS) occultation technique.
Improved formulas are developed for the radial, along-track, and cross-track perturbations, which are more accurate than the
commonly used order-zero formulas. The formulas are used to simulate gravity recovery using the geodetic GPS data of COSMIC
in the operational phase. Results show that the EGM96 model can be improved up to degree 26 using 1 year of COSMIC data. TOPEX/POSEIDON
altimeter data are used to derive a temporal gravity variation. COSMIC cannot reproduce this gravity variation perfectly because
of data noise and orbital configuration, but the recovered field clearly shows the gravity signature due to mass movement
in an El Ni?o.
Received: 3 March 2000 / Accepted: 10 November 2000 相似文献
13.
F. L. Clarke 《Journal of Geodesy》1981,55(1):1-16
Comparisons of gravimetric and astrogeodetic deflections of the vertical in the Australian region indicate that the former
are affected by position dependent systematic errors, even after orientation onto the Australian Geodetic Datum. These are
probably due to errors in the predicted mean anomalies for gravimetrically unsurveyed oceanic regions to the east, south and
west of the continent. Deflection component residuals (astrogeodetic minus oriented gravimetric) at 83 control stations are
made the observables in a set of observation equations, based on the Vening Meinesz equations, from which pseudocorrections
to the mean anomalies for a set of arbitrarily selected surface elements are computed. These pseudocorrections compensate
for prediction errors in much larger unsurveyed regions. Their effects on individual deflection components are calculated
using the Vening Meinesz equations. Statistical tests indicate that pseudocorrections computed for four large offshore elements
and six smaller elements in unsurveyed areas produce corrections to the gravimetric deflections which make the ξ and η components
in seconds of arc consistent with normally distributed populations N (0.00, 0.702). 相似文献
14.
Demosthenes C. Christodoulidis 《Journal of Geodesy》1979,53(1):61-77
Seasonal and latitude dependent corrections to the gravity and height anomalies are developed in order to account for the
neglect of the atmospheric masses outside the geold, when using Stokes’ equation. It is shown that the atmospheric correction
to gravity at sea level is almost constant, equal to0.871 mgals with a variation of2 μ gals whereas the height anomaly correction varies between −0.1 cm and −1.3 cm. Further, when the combined latitudinal/seasonal dependence is neglected in the atmospheric corrections, the maximum error
introduced is of the order of40 μ gals for the gravity corrections and0.7 cm for the height anomaly corrections. 相似文献
15.
Z. Martinec 《Journal of Geodesy》1998,72(7-8):460-472
Green's function for the boundary-value problem of Stokes's type with ellipsoidal corrections in the boundary condition for
anomalous gravity is constructed in a closed form. The `spherical-ellipsoidal' Stokes function describing the effect of two
ellipsoidal correcting terms occurring in the boundary condition for anomalous gravity is expressed in O(e
2
0)-approximation as a finite sum of elementary functions analytically representing the behaviour of the integration kernel
at the singular point ψ=0. We show that the `spherical-ellipsoidal' Stokes function has only a logarithmic singularity in
the vicinity of its singular point. The constructed Green function enables us to avoid applying an iterative approach to solve
Stokes's boundary-value problem with ellipsoidal correction terms involved in the boundary condition for anomalous gravity.
A new Green-function approach is more convenient from the numerical point of view since the solution of the boundary-value
problem is determined in one step by computing a Stokes-type integral. The question of the convergence of an iterative scheme
recommended so far to solve this boundary-value problem is thus irrelevant.
Received: 5 June 1997 / Accepted: 20 February 1998 相似文献
16.
The upward-downward continuation of a harmonic function like the gravitational potential is conventionally based on the direct-inverse
Abel-Poisson integral with respect to a sphere of reference. Here we aim at an error estimation of the “planar approximation”
of the Abel-Poisson kernel, which is often used due to its convolution form. Such a convolution form is a prerequisite to
applying fast Fourier transformation techniques. By means of an oblique azimuthal map projection / projection onto the local
tangent plane at an evaluation point of the reference sphere of type “equiareal” we arrive at a rigorous transformation of
the Abel-Poisson kernel/Abel-Poisson integral in a convolution form. As soon as we expand the “equiareal” Abel-Poisson kernel/Abel-Poisson
integral we gain the “planar approximation”. The differences between the exact Abel-Poisson kernel of type “equiareal” and
the “planar approximation” are plotted and tabulated. Six configurations are studied in detail in order to document the error
budget, which varies from 0.1% for points at a spherical height H=10km above the terrestrial reference sphere up to 98% for points at a spherical height H = 6.3×106km.
Received: 18 March 1997 / Accepted: 19 January 1998 相似文献
17.
P. J. G. Teunissen 《Journal of Geodesy》2001,75(7-8):399-407
Carrier phase ambiguity resolution is the key to fast and high-precision GNSS (Global Navigation Satellite System) kinematic
positioning. Critical in the application of ambiguity resolution is the quality of the computed integer ambiguities. Unsuccessful
ambiguity resolution, when passed unnoticed, will too often lead to unacceptable errors in the positioning results. Very high
success rates are therefore required for ambiguity resolution to be reliable. Biases which are unaccounted for will lower
the success rate and thus increase the chance of unsuccessful ambiguity resolution. The performance of integer ambiguity estimation
in the presence of such biases is studied. Particular attention is given to integer rounding, integer bootstrapping and integer
least squares. Lower and upper bounds, as well as an exact and easy-to-compute formula for the bias-affected success rate,
are presented. These results will enable the evaluation of the bias robustness of ambiguity resolution.
Received: 28 September 2000 / Accepted: 29 March 2001 相似文献
18.
H. Abd-Elmotaal 《Journal of Geodesy》2000,74(5):390-398
Inverse problems in isostasy will consist in making the isostatic anomalies to be zero under a certain isostatic hypothesis.
In the case of the Vening Meinesz isostatic hypothesis, the density contrast is constant, while the Moho depth (depth of the
Mohorovičić discontinuity) is variable. Hence, the Vening Meinesz inverse isostatic problem aims to determine a suitable variable
Moho depth for a prescribed constant density contrast. The main idea is easy but the theoretical analysis is somewhat difficult.
Moreover, the practical determination of the variable Moho depths based on the Vening Meinesz inverse problem is a laborious
and time-consuming task. The formulas used for computing the inverse Vening Meinesz Moho depths are derived. The computational
tricks essentially needed for computing the inverse Vening Meinesz Moho depths from a set of local and global Bouguer anomalies
are described. The Moho depths for a test area are computed based on the inverse Vening Meinesz isostatic problem. These Moho
depths fit the Moho depths derived from seismic observations with a good accuracy, in which the parameters used for the fitting
agree well with those determined geophysically.
Received: 4 February 1999 / Accepted: 4 October 1999 相似文献
19.
XU Xinyu LI Jiancheng ZOU Xiancai CHU Yonghai 《地球空间信息科学学报》2007,10(3):168-172
The principle and method for solving three types of satellite gravity gradient boundary value problems by least-squares are discussed in detail. Also, kernel function expressions of the least-squares solution of three geodetic boundary value problems with the observations {Γ zz },{Γ xz , Γ yz} and {Γ xx -Γ yy ,2 Γxy}are presented. From the results of recovering gravity field using simulated gravity gradient tensor data, we can draw a conclusion that satellite gravity gradient integral formulas derived from least-squares are valid and rigorous for recovering the gravity field. 相似文献
20.
New integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients are derived in this article. They provide more options for continuation of gravitational gradient combinations and extend available mathematical apparatus formulated for this purpose up to now. The starting point represents the analytical solution of the spherical gradiometric boundary value problem in the spatial domain. Applying corresponding differential operators on the analytical solution of the spherical gradiometric boundary value problem, a total of 18 integral formulas are provided. Spatial and spectral forms of isotropic kernels are given and their behaviour for parameters of a GOCE-like satellite is investigated. Correctness of the new integral formulas and the isotropic kernels is tested in a closed-loop simulation. The derived integral formulas and the isotropic kernels form a theoretical basis for validation purposes and geophysical applications of satellite gradiometric data as provided currently by the GOCE mission. They also extend the well-known Meissl scheme. 相似文献