共查询到20条相似文献,搜索用时 250 毫秒
1.
A comparison of the tesseroid,prism and point-mass approaches for mass reductions in gravity field modelling 总被引:6,自引:3,他引:6
The calculation of topographic (and iso- static) reductions is one of the most time-consuming operations in gravity field
modelling. For this calculation, the topographic surface of the Earth is often divided with respect to geographical or map-grid
lines, and the topographic heights are averaged over the respective grid elements. The bodies bounded by surfaces of constant
(ellipsoidal) heights and geographical grid lines are denoted as tesseroids. Usually these ellipsoidal (or spherical) tesseroids
are replaced by “equivalent” vertical rectangular prisms of the same mass. This approximation is motivated by the fact that
the volume integrals for the calculation of the potential and its derivatives can be exactly solved for rectangular prisms,
but not for the tesseroids. In this paper, an approximate solution of the spherical tesseroid integrals is provided based
on series expansions including third-order terms. By choosing the geometrical centre of the tesseroid as the Taylor expansion
point, the number of non-vanishing series terms can be greatly reduced. The zero-order term is equivalent to the point-mass
formula. Test computations show the high numerical efficiency of the tesseroid method versus the prism approach, both regarding
computation time and accuracy. Since the approximation errors due to the truncation of the Taylor series decrease very quickly
with increasing distance of the tesseroid from the computation point, only the elements in the direct vicinity of the computation
point have to be separately evaluated, e.g. by the prism formulas. The results are also compared with the point-mass formula.
Further potential refinements of the tesseroid approach, such as considering ellipsoidal tesseroids, are indicated. 相似文献
2.
Topographic–isostatic masses represent an important source of gravity field information, especially in the high-frequency
band, even if the detailed mass-density distribution inside the topographic masses is unknown. If this information is used
within a remove-restore procedure, then the instability problems in downward continuation of gravity observations from aircraft
or satellite altitudes can be reduced. In this article, integral formulae are derived for determination of gravitational effects
of topographic–isostatic masses on the first- and second-order derivatives of the gravitational potential for three topographic–isostatic
models. The application of these formulas is useful for airborne gravimetry/gradiometry and satellite gravity gradiometry.
The formulas are presented in spherical approximation by separating the 3D integration in an analytical integration in the
radial direction and 2D integration over the mean sphere. Therefore, spherical volume elements can be considered as being
approximated by mass-lines located at the centre of the discretization compartments (the mass of the tesseroid is condensed
mathematically along its vertical axis). The errors of this approximation are investigated for the second-order derivatives
of the topographic–isostatic gravitational potential in the vicinity of the Earth’s surface. The formulas are then applied
to various scenarios of airborne gravimetry/gradiometry and satellite gradiometry. The components of the gravitational vector
at aircraft altitudes of 4 and 10 km have been determined, as well as the gravitational tensor components at a satellite altitude
of 250 km envisaged for the forthcoming GOCE (gravity field and steady-state ocean-circulation explorer) mission. The numerical
computations are based on digital elevation models with a 5-arc-minute resolution for satellite gravity gradiometry and 1-arc-minute
resolution for airborne gravity/gradiometry. 相似文献
3.
A new methodology for computing the gravitational effect of a spherical tesseroid has been devised and implemented. The methodology is based on the rotation from the global Earth-Centred Rotational reference frame to the local Earth-Centred P-Rotational reference frame, referred to the computation point P, and it requires knowledge of the height and the angular extension of each topographic column. After rotation, the gravitational effect of the tesseroid is computed via the effect of a sector of the spherical zonal band. In this respect, two possible procedures for handling the rotated tesseroids have been proposed and tested. The results obtained with the devised methodology are in good agreement with those derived by applying other existing methodologies. 相似文献
4.
G. L. Ramillien 《Journal of Geodesy》2017,91(8):881-895
A radial integration of spherical mass elements (i.e. tesseroids) is presented for evaluating the six components of the second-order gravity gradient (i.e. second derivatives of the Newtonian mass integral for the gravitational potential) created by an uneven spherical topography consisting of juxtaposed vertical prisms. The method uses Legendre polynomial series and takes elastic compensation of the topography by the Earth’s surface into account. The speed of computation of the polynomial series increases logically with the observing altitude from the source of anomaly. Such a forward modelling can be easily applied for reduction of observed gravity gradient anomalies by the effects of any spherical interface of density. An iterative least-squares inversion of measured gravity gradient coefficients is also proposed to estimate a regional set of juxtaposed topographic heights. Several tests of recovery have been made by considering simulated gradients created by idealistic conical and irregular Great Meteor seamount topographies, and for varying satellite altitudes and testing different levels of uncertainty. In the case of gravity gradients measured at a GOCE-type altitude of \(\sim \)300 km, the search converges down to a stable but smooth topography after 10–15 iterations, while the final root-mean-square error is \(\sim \)100 m that represents only 2 % of the seamount amplitude. This recovery error decreases with the altitude of the gravity gradient observations by revealing more topographic details in the region of survey. 相似文献
5.
Toshio Fukushima 《Journal of Geodesy》2018,92(12):1371-1386
We developed an accurate method to compute the gravitational field of a tesseroid. The method numerically integrates a surface integral representation of the gravitational potential of the tesseroid by conditionally splitting its line integration intervals and by using the double exponential quadrature rule. Then, it evaluates the gravitational acceleration vector and the gravity gradient tensor by numerically differentiating the numerically integrated potential. The numerical differentiation is conducted by appropriately switching the central and the single-sided second-order difference formulas with a suitable choice of the test argument displacement. If necessary, the new method is extended to the case of a general tesseroid with the variable density profile, the variable surface height functions, and/or the variable intervals in longitude or in latitude. The new method is capable of computing the gravitational field of the tesseroid independently on the location of the evaluation point, namely whether outside, near the surface of, on the surface of, or inside the tesseroid. The achievable precision is 14–15 digits for the potential, 9–11 digits for the acceleration vector, and 6–8 digits for the gradient tensor in the double precision environment. The correct digits are roughly doubled if employing the quadruple precision computation. The new method provides a reliable procedure to compute the topographic gravitational field, especially that near, on, and below the surface. Also, it could potentially serve as a sure reference to complement and elaborate the existing approaches using the Gauss–Legendre quadrature or other standard methods of numerical integration. 相似文献
6.
Far-zone effects for different topographic-compensation models based on a spherical harmonic expansion of the topography 总被引:1,自引:1,他引:0
The determination of the gravimetric geoid is based on the magnitude of gravity observed at the surface of the Earth or at
airborne altitude. To apply the Stokes’s or Hotine’s formulae at the geoid, the potential outside the geoid must be harmonic
and the observed gravity must be reduced to the geoid. For this reason, the topographic (and atmospheric) masses outside the
geoid must be “condensed” or “shifted” inside the geoid so that the disturbing gravity potential T fulfills Laplace’s equation everywhere outside the geoid. The gravitational effects of the topographic-compensation masses
can also be used to subtract these high-frequent gravity signals from the airborne observations and to simplify the downward
continuation procedures. The effects of the topographic-compensation masses can be calculated by numerical integration based
on a digital terrain model or by representing the topographic masses by a spherical harmonic expansion. To reduce the computation
time in the former case, the integration over the Earth can be divided into two parts: a spherical cap around the computation
point, called the near zone, and the rest of the world, called the far zone. The latter one can be also represented by a global
spherical harmonic expansion. This can be performed by a Molodenskii-type spectral approach. This article extends the original
approach derived in Novák et al. (J Geod 75(9–10):491–504, 2001), which is restricted to determine the far-zone effects for
Helmert’s second method of condensation for ground gravimetry. Here formulae for the far-zone effects of the global topography
on gravity and geoidal heights for Helmert’s first method of condensation as well as for the Airy-Heiskanen model are presented
and some improvements given. Furthermore, this approach is generalized for determining the far-zone effects at aeroplane altitudes.
Numerical results for a part of the Canadian Rocky Mountains are presented to illustrate the size and distributions of these
effects. 相似文献
7.
非心摄动引力的快速计算方法研究 总被引:4,自引:0,他引:4
给出了球谐函数不含阶次 (n,m)调制的递推公式 ,推导出非心引力矢量、引力张量的快速计算格式 ,给出了相应的算法。该算法优于传统正常化递推求和算法 ,减少了运算次数 ,使计算速度提高了 5倍。对低轨卫星预报、卫星重力测量反演、动力法定轨等的响应时间都具有重要贡献。 相似文献
8.
Optimized formulas for the gravitational field of a tesseroid 总被引:7,自引:3,他引:4
Various tasks in geodesy, geophysics, and related geosciences require precise information on the impact of mass distributions on gravity field-related quantities, such as the gravitational potential and its partial derivatives. Using forward modeling based on Newton’s integral, mass distributions are generally decomposed into regular elementary bodies. In classical approaches, prisms or point mass approximations are mostly utilized. Considering the effect of the sphericity of the Earth, alternative mass modeling methods based on tesseroid bodies (spherical prisms) should be taken into account, particularly in regional and global applications. Expressions for the gravitational field of a point mass are relatively simple when formulated in Cartesian coordinates. In the case of integrating over a tesseroid volume bounded by geocentric spherical coordinates, it will be shown that it is also beneficial to represent the integral kernel in terms of Cartesian coordinates. This considerably simplifies the determination of the tesseroid’s potential derivatives in comparison with previously published methodologies that make use of integral kernels expressed in spherical coordinates. Based on this idea, optimized formulas for the gravitational potential of a homogeneous tesseroid and its derivatives up to second-order are elaborated in this paper. These new formulas do not suffer from the polar singularity of the spherical coordinate system and can, therefore, be evaluated for any position on the globe. Since integrals over tesseroid volumes cannot be solved analytically, the numerical evaluation is achieved by means of expanding the integral kernel in a Taylor series with fourth-order error in the spatial coordinates of the integration point. As the structure of the Cartesian integral kernel is substantially simplified, Taylor coefficients can be represented in a compact and computationally attractive form. Thus, the use of the optimized tesseroid formulas particularly benefits from a significant decrease in computation time by about 45 % compared to previously used algorithms. In order to show the computational efficiency and to validate the mathematical derivations, the new tesseroid formulas are applied to two realistic numerical experiments and are compared to previously published tesseroid methods and the conventional prism approach. 相似文献
9.
E. Grafarend 《Journal of Geodesy》1973,47(3):237-260
Based on exterior calculus, the G. Frobenius integration theorem, holonomic and anholonomic Riemannian geometry, the typical
geodetic problems are summarized in a unified manner. The E. Cartan pseudotorsion of natural orthogonal coordinates causes
the misclosure of a closed three dimensional traverse. Natural coordinate differences are path dependent, anholonomic, nonintegrable,
nonunique, therefore. The geodetic pseudotorsion form depends only on the components of the A. Marussi tensor of gravity gradients.
A physically defined coordinate system can be found which is pseudotorsion free, whose coordinates are holonomic, integrable,
unique. The G. Frobenius transformation matrix is of rank three, explaining the number of three dimensions of an intrinsic
surface geometry. The matrix elements depend on either the second derivatives of the real gravity potential and the Euclidean
norm of its gravity vector or the second derivatives of the standard gravity potential, the Euclidean norm of its standard
gravity vector and the vertical deflections. Incomplete information of the earth's gravity field leads to the concept of boundary
value problems and satellite geodesy.
相似文献
10.
The topographic effects by Stokes formula are typically considered for a spherical approximation of sea level. For more precise determination of the geoid, sea level is better approximated by an ellipsoid, which justifies the consideration of the ellipsoidal corrections of topographic effects for improved geoid solutions. The aim of this study is to estimate the ellipsoidal effects of the combined topographic correction (direct plus indirect topographic effects) and the downward continuation effect. It is concluded that the ellipsoidal correction to the combined topographic effect on the geoid height is far less than 1 mm. On the contrary, the ellipsoidal correction to the effect of downward continuation of gravity anomaly to sea level may be significant at the 1-cm level in mountainous regions. Nevertheless, if Stokes formula is modified and the integration of gravity anomalies is limited to a cap of a few degrees radius around the computation point, nor this effect is likely to be significant.AcknowledgementsThe author is grateful for constructive remarks by J Ågren and the three reviewers. 相似文献
11.
Georges Blaha 《Journal of Geodesy》1980,54(1):119-135
Second-order derivatives of a general scalar function of position (F) with respect to the length elements along a family of local Cartesian axes are developed in the spheroidal and spherical
coordinate systems. A link between the two kinds of formulations is established when the results in spherical coordinates
are confirmed also indirectly, through a transformation from spheroidal coordinates. IfF becomesW (earth's potential) the six distinct second-order derivatives—which include one vertical and two horizontal gradients of
gravity—relate the symmetric Marussi tensor to the curvature parameters of the field.
The general formulas for the second-order derivatives ofF are specialized to yield the second-order derivatives ofU (standard potential) and ofT (disturbing potential), which allows the latter to be modeled by a suitable set of parameters. The second-order derivatives
ofT in which the property ΔT=0 is explicitly incorporated are also given. According to the required precision, the spherical approximation may or may not
be desirable; both kinds of results are presented. The derived formulas can be used for modeling of the second-order derivatives
ofW orT at the ground level as well as at higher altitudes. They can be further applied in a rotating or a nonrotating field. The
development in this paper is based on the tensor approach to theoretical geodesy, introduced by Marussi [1951] and further
elaborated by Hotine [1969], which can lead to significantly shorter demonstrations when compared to conventional approaches. 相似文献
12.
Spherical cap harmonic analysis: a comment on its proper use for local gravity field representation 总被引:1,自引:0,他引:1
Spherical cap harmonic analysis is the appropriate analytical technique for modelling Laplacian potential and the corresponding
field components over a spherical cap. This paper describes the use of this method by means of a least-squares approach for
local gravity field representation. Formulations for the geoid undulation and the components ξ, η of the deflection of the
vertical are derived, together with some warnings in the application of the technique. Although most of the formulations have
been given by another paper, these were confusing or even incorrect, mainly because of an improper application of the spherical
cap harmonic analysis.
Received: 16 January 1996 / Accepted: 17 March 1997 相似文献
13.
Gravity field terrain effect computations by FFT 总被引:2,自引:2,他引:2
René Forsberg 《Journal of Geodesy》1985,59(4):342-360
The widespread availability of detailed gridded topographic and bathymetric data for many areas of the earth has resulted
in a need for efficient terrain effect computation techniques, especially for applications in gravity field modelling. Compared
to conventional integration techniques, Fourier transform methods provide extremely efficient computations due to the speed
of the Fast Fourier Transform (FFT. The Fourier techniques rely on linearization and series expansions of the basically unlinear terrain effect integrals, typically
involving transformation of the heights/depths and their squares. TheFFT methods will especially be suited for terrain reduction of land gravity data and satellite altimetry geoid data.
In the paper the basic formulas will be outlined, and special emphasis will be put on the practial implementation, where a
special coarse/detailed grid pair formulation must be used in order to minimize the unavoidable edge effects ofFFT, and the special properties ofFFT are utilized to limit the actual number of data transformations needed. Actual results are presented for gravity and geoid
terrain effects in test areas of the USA, Greenland and the North Atlantic. The results are evaluated against a conventional
integration program: thus, e.g., in an area of East Greenland (with terrain corrections up to10 mgal), the accuracy ofFFT-computed terrain corrections in actual gravity stations showed anr.m.s. error of0.25 mgal, using height data from a detailed photogrammetric digital terrain model. Similarly, isostatic ocean geoid effects in the
Faeroe Islands region were found to be computed withr.m.s. errors around0.03 m 相似文献
14.
L. E. Sjöberg 《Journal of Geodesy》2007,81(5):345-350
This study emphasizes that the harmonic downward continuation of an external representation of the Earth’s gravity potential
to sea level through the topographic masses implies a topographic bias. It is shown that the bias is only dependent on the
topographic density along the geocentric radius at the computation point. The bias corresponds to the combined topographic
geoid effect, i.e., the sum of the direct and indirect topographic effects. For a laterally variable topographic density function,
the combined geoid effect is proportional to terms of powers two and three of the topographic height, while all higher order
terms vanish. The result is useful in geoid determination by analytical continuation, e.g., from an Earth gravity model, Stokes’s
formula or a combination thereof. 相似文献
15.
A method is presented with which to verify that the computer software used to compute a gravimetric geoid is capable of producing
the correct results, assuming accurate input data. The Stokes, gravimetric terrain correction and indirect effect formulae
are integrated analytically after applying a transformation to surface spherical coordinates centred on each computation point.
These analytical results can be compared with those from geoid computation software using constant gravity data in order to
verify its integrity. Results of tests conducted with geoid computation software are presented which illustrate the need for
integration weighting factors, especially for those compartments close to the computation point.
Received: 6 February 1996 / Accepted: 19 April 1997 相似文献
16.
17.
Ellipsoidal geoid computation 总被引:1,自引:1,他引:0
Modern geoid computation uses a global gravity model, such as EGM96, as a third component in a remove–restore process. The classical approach uses only two: the reference ellipsoid and a geometrical model representing the topography. The rationale for all three components is reviewed, drawing attention to the much smaller precision now needed when transforming residual gravity anomalies. It is shown that all ellipsoidal effects needed for geoid computation with millimetric accuracy are automatically included provided that the free air anomaly and geoid are calculated correctly from the global model. Both must be consistent with an ellipsoidal Earth and with the treatment of observed gravity data. Further ellipsoidal corrections are then negligible. Precise formulae are developed for the geoid height and the free air anomaly using a global gravity model, given as spherical harmonic coefficients. Although only linear in the anomalous potential, these formulae are otherwise exact for an ellipsoidal reference Earth—they involve closed analytical functions of the eccentricity (and the Earths spin rate), rather than a truncated power series in e2. They are evaluated using EGM96 and give ellipsoidal corrections to the conventional free air anomaly ranging from –0.84 to +1.14 mGal, both extremes occurring in Tibet. The geoid error corresponding to these differences is dominated by longer wavelengths, so extrema occur elsewhere, rising to +766 mm south of India and falling to –594 mm over New Guinea. At short wavelengths, the difference between ellipsoidal corrections based only on EGM96 and those derived from detailed local gravity data for the North Sea geoid GEONZ97 has a standard deviation of only 3.3 mm. However, the long-wavelength components missed by the local computation reach 300 mm and have a significant slope. In Australia, for example, such a slope would amount to a 600-mm rise from Perth to Cairns. 相似文献
18.
19.
D. Arabelos 《Journal of Geodesy》1985,59(2):109-123
The evaluation of deflections of the vertical for the area of Greece is attempted using a combination of topographic and astrogeodetic
data. Tests carried out in the area bounded by 35°≤ϕ≤42°, 19°≤λ≤27° indicate that an accuracy of ±3″.3 can be obtained in
this area for the meridian and prime vertical deflection components when high resolution topographic data in the immediate
vicinity of computation points are used, combined with high degree spherical harmonic expansions of the geopotential and isostatic
reduction potential. This accuracy is about 25% better than the corresponding topographic-Moho deflection components which
are evaluated using topographic and Moho data up to 120 km around each station, without any combination with the spherical
harmonic expansion of the geopotential or isostatic reduction potential. The accuracy in both cases is increased to about
2″.6 when the astrogeodetic data available in the area mentioned above are used for the prediction of remaining values. Furthermore
the estimation of datum-shift parameters is attempted using least squares collocation. 相似文献
20.
A synthetic Earth Gravity Model Designed Specifically for Testing Regional Gravimetric Geoid Determination Algorithms 总被引:1,自引:0,他引:1
I. Baran M. Kuhn S. J. Claessens W. E. Featherstone S. A. Holmes P. Vaníček 《Journal of Geodesy》2006,80(1):1-16
A synthetic [simulated] Earth gravity model (SEGM) of the geoid, gravity and topography has been constructed over Australia specifically for validating regional gravimetric geoid determination theories, techniques and computer software. This regional high-resolution (1-arc-min by 1-arc-min) Australian SEGM (AusSEGM) is a combined source and effect model. The long-wavelength effect part (up to and including spherical harmonic degree and order 360) is taken from an assumed errorless EGM96 global geopotential model. Using forward modelling via numerical Newtonian integration, the short-wavelength source part is computed from a high-resolution (3-arc-sec by 3-arc-sec) synthetic digital elevation model (SDEM), which is a fractal surface based on the GLOBE v1 DEM. All topographic masses are modelled with a constant mass-density of 2,670 kg/m3. Based on these input data, gravity values on the synthetic topography (on a grid and at arbitrarily distributed discrete points) and consistent geoidal heights at regular 1-arc-min geographical grid nodes have been computed. The precision of the synthetic gravity and geoid data (after a first iteration) is estimated to be better than 30 μ Gal and 3 mm, respectively, which reduces to 1 μ Gal and 1 mm after a second iteration. The second iteration accounts for the changes in the geoid due to the superposed synthetic topographic mass distribution. The first iteration of AusSEGM is compared with Australian gravity and GPS-levelling data to verify that it gives a realistic representation of the Earth’s gravity field. As a by-product of this comparison, AusSEGM gives further evidence of the north–south-trending error in the Australian Height Datum. The freely available AusSEGM-derived gravity and SDEM data, included as Electronic Supplementary Material (ESM) with this paper, can be used to compute a geoid model that, if correct, will agree to in 3 mm with the AusSEGM geoidal heights, thus offering independent verification of theories and numerical techniques used for regional geoid modelling.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00190-005-0002-z 相似文献