共查询到20条相似文献，搜索用时 31 毫秒
1.
L. E. Sjöberg 《Journal of Geodesy》2003,77(3-4):139-147
Assuming that the gravity anomaly and disturbing potential are given on a reference ellipsoid, the result of Sjöberg (1988, Bull Geod 62:93–101) is applied to derive the potential coefficients on the bounding sphere of the ellipsoid to order e
2 (i.e. the square of the eccentricity of the ellipsoid). By adding the potential coefficients and continuing the potential downward to the reference ellipsoid, the spherical Stokes formula and its ellipsoidal correction are obtained. The correction is presented in terms of an integral over the unit sphere with the spherical approximation of geoidal height as the argument and only three well-known kernel functions, namely those of Stokes, Vening-Meinesz and the inverse Stokes, lending the correction to practical computations. Finally, the ellipsoidal correction is presented also in terms of spherical harmonic functions. The frequently applied and sometimes questioned approximation of the constant m, a convenient abbreviation in normal gravity field representations, by e
2/2, as introduced by Moritz, is also discussed. It is concluded that this approximation does not significantly affect the ellipsoidal corrections to potential coefficients and Stokes formula. However, whether this standard approach to correct the gravity anomaly agrees with the pure ellipsoidal solution to Stokes formula is still an open question. 相似文献
2.
L. E. Sjöberg 《Journal of Geodesy》2000,74(2):255-268
The topographic potential and the direct topographic effect on the geoid are presented as surface integrals, and the direct
gravity effect is derived as a rigorous surface integral on the unit sphere. By Taylor-expanding the integrals at sea level
with respect to topographic elevation (H) the power series of the effects is derived to arbitrary orders. This study is primarily limited to terms of order H
2. The limitations of the various effects in the frequently used planar approximations are demonstrated. In contrast, it is
shown that the spherical approximation to power H
2 leads to a combined topographic effect on the geoid (direct plus indirect effect) proportional to H˜2 (where terms of degrees 0 and 1 are missing) of the order of several metres, while the combined topographic effect on the
height anomaly vanishes, implying that current frequent efforts to determine the direct effect to this order are not needed.
The last result is in total agreement with Bjerhammar's method in physical geodesy. It is shown that the most frequently applied
remove–restore technique of topographic masses in the application of Stokes' formula suffers from significant errors both
in the terrain correction C (representing the sum of the direct topographic effect on gravity anomaly and the effect of continuing the anomaly to sea
level) and in the term t (mainly representing the indirect effect on the geoidal or quasi-geoidal height).
Received: 18 August 1998 / Accepted: 4 October 1999 相似文献
3.
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 相似文献
4.
H. Nahavandchi 《Journal of Geodesy》2002,76(6-7):345-352
It is suggested that a spherical harmonic representation of the geoidal heights using global Earth gravity models (EGM) might
be accurate enough for many applications, although we know that some short-wavelength signals are missing in a potential coefficient
model. A `direct' method of geoidal height determination from a global Earth gravity model coefficient alone and an `indirect'
approach of geoidal height determination through height anomaly computed from a global gravity model are investigated. In
both methods, suitable correction terms are applied. The results of computations in two test areas show that the direct and
indirect approaches of geoid height determination yield good agreement with the classical gravimetric geoidal heights which
are determined from Stokes' formula. Surprisingly, the results of the indirect method of geoidal height determination yield
better agreement with the global positioning system (GPS)-levelling derived geoid heights, which are used to demonstrate such
improvements, than the results of gravimetric geoid heights at to the same GPS stations. It has been demonstrated that the
application of correction terms in both methods improves the agreement of geoidal heights at GPS-levelling stations. It is
also found that the correction terms in the direct method of geoidal height determination are mostly similar to the correction
terms used for the indirect determination of geoidal heights from height anomalies.
Received: 26 July 2001 / Accepted: 21 February 2002 相似文献
5.
The standard analytical approach which is applied for constructing geopotential models OSU86 and earlier ones, is based on
reducing the boundary value equation to a sphere enveloping the Earth and then solving it directly with respect to the potential
coefficients
n,m
. In an alternative procedure, developed by Jekeli and used for constructing the models OSU91 and EGM96, at first an ellipsoidal
harmonic series is developed for the geopotential and then its coefficients
n,m
e
are transformed to the unknown
n,m
. The second solution is more exact, but much more complicated. The standard procedure is modified and a new simple integral
formula is derived for evaluating the potential coefficients. The efficiency of the standard and new procedures is studied
numerically. In these solutions the same input data are used as for constructing high-degree parts of the EGM96 models. From
two sets of
n,m
(n≤360,|m|≤n), derived by the standard and new approaches, different spectral characteristics of the gravity anomaly and the geoid undulation
are estimated and then compared with similar characteristics evaluated by Jekeli's approach (`etalon' solution). The new solution
appears to be very close to Jekeli's, as opposed to the standard solution. The discrepancies between all the characteristics
of the new and `etalon' solutions are smaller than the corresponding discrepancies between two versions of the final geopotential
model EGM96, one of them (HDM190) constructed by the block-diagonal least squares (LS) adjustment and the other one (V068)
by using Jekeli's approach. On the basis of the derived analytical solution a new simple mathematical model is developed to
apply the LS technique for evaluating geopotential coefficients.
Received: 12 December 2000 / Accepted: 21 June 2001 相似文献
6.
When standard boundary element methods (BEM) are used in order to solve the linearized vector Molodensky problem we are confronted with
two problems: (1) the absence of O(|x|−2) terms in the decay condition is not taken into account, since the single-layer ansatz, which is commonly used as representation
of the disturbing potential, is of the order O(|x|−1) as x→∞. This implies that the standard theory of Galerkin BEM is not applicable since the injectivity of the integral operator
fails; (2) the N×N stiffness matrix is dense, with N typically of the order 105. Without fast algorithms, which provide suitable approximations to the stiffness matrix by a sparse one with O(N(logN)
s
), s≥0, non-zero elements, high-resolution global gravity field recovery is not feasible. Solutions to both problems are proposed.
(1) A proper variational formulation taking the decay condition into account is based on some closed subspace of co-dimension
3 of the space of square integrable functions on the boundary surface. Instead of imposing the constraints directly on the
boundary element trial space, they are incorporated into a variational formulation by penalization with a Lagrange multiplier.
The conforming discretization yields an augmented linear system of equations of dimension N+3×N+3. The penalty term guarantees the well-posedness of the problem, and gives precise information about the incompatibility
of the data. (2) Since the upper left submatrix of dimension N×N of the augmented system is the stiffness matrix of the standard BEM, the approach allows all techniques to be used to generate
sparse approximations to the stiffness matrix, such as wavelets, fast multipole methods, panel clustering etc., without any
modification. A combination of panel clustering and fast multipole method is used in order to solve the augmented linear system
of equations in O(N) operations. The method is based on an approximation of the kernel function of the integral operator by a degenerate kernel
in the far field, which is provided by a multipole expansion of the kernel function. Numerical experiments show that the fast
algorithm is superior to the standard BEM algorithm in terms of CPU time by about three orders of magnitude for N=65 538 unknowns. Similar holds for the storage requirements. About 30 iterations are necessary in order to solve the linear
system of equations using the generalized minimum residual method (GMRES). The number of iterations is almost independent
of the number of unknowns, which indicates good conditioning of the system matrix.
Received: 16 October 1999 / Accepted: 28 February 2001 相似文献
7.
A methodology for precise determination of the fundamental geodetic parameter w
0, the potential value of the Gauss–Listing geoid, as well as its time derivative 0, is presented. The method is based on: (1) ellipsoidal harmonic expansion of the external gravitational field of the Earth
to degree/order 360/360 (130 321 coefficients; http://www.uni-stuttgard.de/gi/research/ index.html projects) with respect
to the International Reference Ellipsoid WGD2000, at the GPS positioned stations; and (2) ellipsoidal free-air gravity reduction
of degree/order 360/360, based on orthometric heights of the GPS-positioned stations. The method has been numerically tested
for the data of three GPS campaigns of the Baltic Sea Level project (epochs 1990.8,1993.4 and 1997.4). New w
0 and 0 values (w
0=62 636 855.75 ± 0.21 m2/s2, 0=−0.0099±0.00079 m2/s2 per year, w
0/&γmacr;=6 379 781.502 m,0/&γmacr;=1.0 mm/year, and &γmacr;= −9.81802523 m2/s2) for the test region (Baltic Sea) were obtained. As by-products of the main study, the following were also determined: (1)
the high-resolution sea surface topography map for the Baltic Sea; (2) the most accurate regional geoid amongst four different
regional Gauss–Listing geoids currently proposed for the Baltic Sea; and (3) the difference between the national height datums
of countries around the Baltic Sea.
Received: 14 August 2000 / Accepted: 19 June 2001 相似文献
8.
A solution to the downward continuation effect on the geoid determined by Stokes' formula 总被引：2，自引：1，他引：2
L.E. Sjöberg 《Journal of Geodesy》2003,77(1-2):94-100
The analytical continuation of the surface gravity anomaly to sea level is a necessary correction in the application of Stokes'
formula for geoid estimation. This process is frequently performed by the inversion of Poisson's integral formula for a sphere.
Unfortunately, this integral equation corresponds to an improperly posed problem, and the solution is both numerically unstable,
unless it is well smoothed, and tedious to compute. A solution that avoids the intermediate step of downward continuation
of the gravity anomaly is presented. Instead the effect on the geoid as provided by Stokes' formula is studied directly. The
practical solution is partly presented in terms of a truncated Taylor series and partly as a truncated series of spherical
harmonics. Some simple numerical estimates show that the solution mostly meets the requests of a 1-cm geoid model, but the
truncation error of the far zone must be studied more precisely for high altitudes of the computation point. In addition,
it should be emphasized that the derived solution is more computer efficient than the detour by Poisson's integral.
Received: 6 February 2002 / Accepted: 18 November 2002
Acknowledgements. Jonas ?gren carried out the numerical calculations and gave some critical and constructive remarks on a draft version of
the paper. This support is cordially acknowledged. Also, the thorough work performed by one unknown reviewer is very much
appreciated. 相似文献
9.
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 相似文献
10.
The structure of normal matrices occurring in the problem of weighted least-squares spherical harmonic analysis of measurements
scattered on a sphere with random noises is investigated. Efficient algorithms for the formation of the normal matrices are
derived using fundamental relations inherent to the products of two surface spherical harmonic functions. The whole elements
of a normal matrix complete to spherical harmonic degree L are recursively obtained from its first row or first column extended to degree 2L with only O(L
4) computational operations. Applications of the algorithms to the formation of surface normal matrices from geoid undulations
and surface gravity anomalies are discussed in connection with the high-degree geopotential modeling.
Received: 22 March 1999 / Accepted: 23 December 1999 相似文献
11.
H. Nahavandchi 《Journal of Geodesy》2000,74(6):488-496
The direct topographical correction is composed of both local effects and long-wavelength contributions. This implies that
the classical integral formula for determining the direct effect may have some numerical problems in representing these different
signals. On the other hand, a representation by a set of harmonic coefficients of the topography to, say, degree and order
360 will omit significant short-wavelength signals. A new formula is derived by combining the classical formula and a set
of spherical harmonics. Finally, the results of this solution are compared with the Moritz topographical correction in a test
area.
Received: 27 July 1998 / Accepted: 29 March 2000 相似文献
12.
13.
Lars E. Sjöberg 《Journal of Geodesy》2006,79(12):675-681
The application of Stokes’s formula to determine the geoid height requires that topographic and atmospheric masses be mathematically removed prior to Stokes integration. This corresponds to the applications of the direct topographic and atmospheric effects. For a proper geoid determination, the external masses must then be restored, yielding the indirect effects. Assuming an ellipsoidal layering of the atmosphere with 15% increase in its density towards the poles, the direct atmospheric effect on the geoid height is estimated to be −5.51 m plus a second-degree zonal harmonic term with an amplitude of 1.1 cm. The indirect effect is +5.50 m and the total geoid correction thus varies between −1.2 cm at the equator to 1.9 cm at the poles. Finally, the correction needed to the atmospheric effect if Stokes’s formula is used in a spherical approximation, rather than an ellipsoidal approximation, of the Earth varies between 0.3 cm and 4.0 cm at the equator and pole, respectively. 相似文献
14.
An operational algorithm for computation of terrain correction (or local gravity field modeling) based on application of closed-form solution of the Newton integral in terms of Cartesian coordinates in multi-cylindrical equal-area map projection of the reference ellipsoid is presented. Multi-cylindrical equal-area map projection of the reference ellipsoid has been derived and is described in detail for the first time. Ellipsoidal mass elements with various sizes on the surface of the reference ellipsoid are selected and the gravitational potential and vector of gravitational intensity (i.e. gravitational acceleration) of the mass elements are computed via numerical solution of the Newton integral in terms of geodetic coordinates {,,h}. Four base- edge points of the ellipsoidal mass elements are transformed into a multi-cylindrical equal-area map projection surface to build Cartesian mass elements by associating the height of the corresponding ellipsoidal mass elements to the transformed area elements. Using the closed-form solution of the Newton integral in terms of Cartesian coordinates, the gravitational potential and vector of gravitational intensity of the transformed Cartesian mass elements are computed and compared with those of the numerical solution of the Newton integral for the ellipsoidal mass elements in terms of geodetic coordinates. Numerical tests indicate that the difference between the two computations, i.e. numerical solution of the Newton integral for ellipsoidal mass elements in terms of geodetic coordinates and closed-form solution of the Newton integral in terms of Cartesian coordinates, in a multi-cylindrical equal-area map projection, is less than 1.6×10–8 m2/s2 for a mass element with a cross section area of 10×10 m and a height of 10,000 m. For a mass element with a cross section area of 1×1 km and a height of 10,000 m the difference is less than 1.5×10–4m2/s2. Since 1.5× 10–4 m2/s2 is equivalent to 1.5×10–5m in the vertical direction, it can be concluded that a method for terrain correction (or local gravity field modeling) based on closed-form solution of the Newton integral in terms of Cartesian coordinates of a multi-cylindrical equal-area map projection of the reference ellipsoid has been developed which has the accuracy of terrain correction (or local gravity field modeling) based on the Newton integral in terms of ellipsoidal coordinates.Acknowledgments. This research has been financially supported by the University of Tehran based on grant number 621/4/859. This support is gratefully acknowledged. The authors are also grateful for the comments and corrections made to the initial version of the paper by Dr. S. Petrovic from GFZ Potsdam and the other two anonymous reviewers. Their comments helped to improve the structure of the paper significantly. 相似文献
15.
Z. Martinec 《Journal of Geodesy》2003,77(1-2):41-49
Three independent gradiometric boundary-value problems (BVPs) with three types of gradiometric data, {Γ
rr
}, {Γ
r
θ,Γ
r
λ} and {Γθθ−Γλλ,Γθλ}, prescribed on a sphere are solved to determine the gravitational potential on and outside the sphere. The existence and
uniqueness conditions on the solutions are formulated showing that the zero- and the first-degree spherical harmonics are
to be removed from {Γ
r
θ,Γ
r
λ} and {Γθθ−Γλλ,Γθλ}, respectively. The solutions to the gradiometric BVPs are presented in terms of Green's functions, which are expressed in
both spectral and closed spatial forms. The logarithmic singularity of the Green's function at the point ψ=0 is investigated
for the component Γ
rr
. The other two Green's functions are finite at this point. Comparisons to the paper by van Gelderen and Rummel [Journal of
Geodesy (2001) 75: 1–11] show that the presented solution refines the former solution.
Received: 3 October 2001 / Accepted: 4 October 2002 相似文献
16.
J. Marshall 《Journal of Geodesy》2002,76(6-7):334-344
Several pre-analysis measures which help to expose the behavior of L
1 -norm minimization solutions are described. The pre-analysis measures are primarily based on familiar elements of the linear
programming solution to L
1-norm minimization, such as slack variables and the reduced-cost vector. By examining certain elements of the linear programming
solution in a probabilistic light, it is possible to derive the cumulative distribution function (CDF) associated with univariate
L
1-norm residuals. Unlike traditional least squares (LS) residual CDFs, it is found that L
1-norm residual CDFs fail to follow the normal distribution in general, and instead are characterized by both discrete and
continuous (i.e. piecewise) segments. It is also found that an L
1 equivalent to LS redundancy numbers exists and that these L
1 equivalents are a byproduct of the univariate L
1 univariate residual CDF. Probing deeper into the linear programming solution, it is found that certain combinations of observations
which are capable of tolerating large-magnitude gross errors can be predicted by comprehensively tabulating the signs of slack
variables associated with the L
1 residuals. The developed techniques are illustrated on a two-dimensional trilateration network.
Received: 6 July 2001 / Accepted: 21 February 2002 相似文献
17.
Traditionally, the evaluation of geoidal height by Stokes formula and the vertical deflection by Vening-Meinesz one, and the estimation of the influence of neglecting the distant zone on computing the geoidal height and the vertical deflection were done by taking the inner zone as a spherical cap. It is not very convenient from the point of view of modern numerical methods such as fast Fourier and Hartley transforms where the inner zone is not a spherical cap, but a spherical trapezoid. So, we generalized the known formulas for evaluating the geoidal height and the vertical deflection for an integration area of arbitrary shape. The corresponding formulas for computing the effects of neglecting the distant zone have been derived. Some issues on computation techniques have been investigated. As an example, the case where the inner zone is modeled as a spherical trapezoid was given special attention, and practical computations were performed. 相似文献
18.
Solutions of the linearized geodetic boundary value problem for an ellipsoidal boundary to order e
3
The geodetic boundary value problem (GBVP) was originally formulated for the topographic surface of the Earth. It degenerates to an ellipsoidal problem, for example when topographic and downward continuation reductions have been applied. Although these ellipsoidal GBVPs possess a simpler structure than the original ones, they cannot be solved analytically, since the boundary condition still contains disturbing terms due to anisotropy, ellipticity and centrifugal components in the reference potential. Solutions of the so-called scalar-free version of the GBVP, upon which most recent practical calculations of geoidal and quasigeoidal heights are based, are considered. Starting at the linearized boundary condition and presupposing a normal field of Somigliana–Pizzetti type, the boundary condition described in spherical coordinates is expanded into a series with respect to the flattening f of the Earth. This series is truncated after the linear terms in f, and first-order solutions of the corresponding GBVP are developed in closed form on the basis of spherical integral formulae, modified by suitable reduction terms. Three alternative representations of the solution are discussed, implying corrections by adding a first-order non-spherical term to the solution, by reducing the boundary data, or by modifying the integration kernel. A numerically efficient procedure for the evaluation of ellipsoidal effects, in the case of the linearized scalar-free version of the GBVP, involving first-order ellipsoidal terms in the boundary condition, is derived, utilizing geopotential models such as EGM96. 相似文献
19.
Iterative vector methods for computing geodetic latitude and height from rectangular coordinates 总被引：4，自引：4，他引：4
J. Pollard 《Journal of Geodesy》2002,76(1):36-40
Two iterative vector methods for computing geodetic coordinates (φ, h) from rectangular coordinates (x, y, z) are presented. The methods are conceptually simple, work without modification at any latitude and are easy to program. Geodetic
latitude and height can be calculated to acceptable precision in one iteration over the height range from −106 to +109 m.
Received: 13 December 2000 / Accepted: 13 July 2001 相似文献
20.
Fast spherical collocation: theory and examples 总被引：2，自引：4，他引：2
It has long been known that a spherical harmonic analysis of gridded (and noisy) data on a sphere (with uniform error for
a fixed latitude) gives rise to simple systems of equations. This idea has been generalized for the method of least-squares
collocation, when using an isotropic covariance function or reproducing kernel. The data only need to be at the same altitude
and of the same kind for each latitude. This permits, for example, the combination of gravity data at the surface of the Earth
and data at satellite altitude, when the orbit is circular. Suppose that data are associated with the points of a grid with
N values in latitude and M values in longitude. The latitudes do not need to be spaced uniformly. Also suppose that it is required to determine the
spherical harmonic coefficients to a maximal degree and order K. Then the method will require that we solve K systems of equations each having a symmetric positive definite matrix of only N × N. Results of simulation studies using the method are described.
Received: 18 October 2001 / Accepted: 4 October 2002
Correspondence to: F. Sansò 相似文献