共查询到20条相似文献,搜索用时 15 毫秒
1.
本文结合椭球域的特点导出了三类椭球域边值问题的准格林函数解。这些解的主项分别是椭球域泊松积分,纽曼积分和司托克斯积分;次项是O(e^2)量级的非谐和级数改正项,可由位系数模型按相应的公式算得。 相似文献
2.
3.
Gravity data observed on or reduced to the ellipsoid are preferably represented using ellipsoidal harmonics instead of spherical harmonics. Ellipsoidal harmonics, however, are difficult to use in practice because the computation of the associated Legendre functions of the second kind that occur in the ellipsoidal harmonic expansions is not straightforward. Jekeli’s renormalization simplifies the computation of the associated Legendre functions. We extended the direct computation of these functions—as well as that of their ratio—up to the second derivatives and minimized the number of required recurrences by a suitable hypergeometric transformation. Compared with the original Jekeli’s renormalization the associated Legendre differential equation is fulfilled up to much higher degrees and orders for our optimized recurrences. The derived functions were tested by comparing functionals of the gravitational potential computed with both ellipsoidal and spherical harmonic syntheses. As an input, the high resolution global gravity field model EGM2008 was used. The relative agreement we found between the results of ellipsoidal and spherical syntheses is 10?14, 10?12 and 10?8 for the potential and its first and second derivatives, respectively. Using the original renormalization, this agreement is 10?12, 10?8 and 10?5, respectively. In addition, our optimized recurrences require less computation time as the number of required terms for the hypergeometric functions is less. 相似文献
4.
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. 相似文献
5.
Comparisons between high-degree models of the Earth’s topographic and gravitational potential may give insight into the quality and resolution of the source data sets, provide feedback on the modelling techniques and help to better understand the gravity field composition. Degree correlations (cross-correlation coefficients) or reduction rates (quantifying the amount of topographic signal contained in the gravitational potential) are indicators used in a number of contemporary studies. However, depending on the modelling techniques and underlying levels of approximation, the correlation at high degrees may vary significantly, as do the conclusions drawn. The present paper addresses this problem by attempting to provide a guide on global correlation measures with particular emphasis on approximation effects and variants of topographic potential modelling. We investigate and discuss the impact of different effects (e.g., truncation of series expansions of the topographic potential, mass compression, ellipsoidal versus spherical approximation, ellipsoidal harmonic coefficient versus spherical harmonic coefficient (SHC) representation) on correlation measures. Our study demonstrates that the correlation coefficients are realistic only when the model’s harmonic coefficients of a given degree are largely independent of the coefficients of other degrees, permitting degree-wise evaluations. This is the case, e.g., when both models are represented in terms of SHCs and spherical approximation (i.e. spherical arrangement of field-generating masses). Alternatively, a representation in ellipsoidal harmonics can be combined with ellipsoidal approximation. The usual ellipsoidal approximation level (i.e. ellipsoidal mass arrangement) is shown to bias correlation coefficients when SHCs are used. Importantly, gravity models from the International Centre for Global Earth Models (ICGEM) are inherently based on this approximation level. A transformation is presented that enables a transformation of ICGEM geopotential models from ellipsoidal to spherical approximation. The transformation is applied to generate a spherical transform of EGM2008 (sphEGM2008) that can meaningfully be correlated degree-wise with the topographic potential. We exploit this new technique and compare a number of models of topographic potential constituents (e.g., potential implied by land topography, ocean water masses) based on the Earth2014 global relief model and a mass-layer forward modelling technique with sphEGM2008. Different to previous findings, our results show very significant short-scale correlation between Earth’s gravitational potential and the potential generated by Earth’s land topography (correlation +0.92, and 60% of EGM2008 signals are delivered through the forward modelling). Our tests reveal that the potential generated by Earth’s oceans water masses is largely unrelated to the geopotential at short scales, suggesting that altimetry-derived gravity and/or bathymetric data sets are significantly underpowered at 5 arc-min scales. We further decompose the topographic potential into the Bouguer shell and terrain correction and show that they are responsible for about 20 and 25% of EGM2008 short-scale signals, respectively. As a general conclusion, the paper shows the importance of using compatible models in topographic/gravitational potential comparisons and recommends the use of SHCs together with spherical approximation or EHCs with ellipsoidal approximation in order to avoid biases in the correlation measures. 相似文献
6.
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 相似文献
7.
On the explicit determination of stability constants for linearized geodetic boundary value problems
The theory of GBVPs provide the basis to the approximate methods used to compute global gravity models. A standard approximation
procedure is least squares, which implicitly assumes that data, e.g. gravity disturbance and gravity anomaly, are given functions
in L
2(S). We know that solutions in these cases exist, but uniqueness (and coerciveness which implies stability of the numerical
solutions) is the real problem. Conditions of uniqueness for the linearized fixed boundary and Molodensky problems are studied
in detail. They depend on the geometry of the boundary; however, the case of linearized fixed boundary BVP puts practically
no constraint on the surface S, while the linearized Molodensky BVP requires the previous knowledge of very low harmonics,
for instance up to degree 25, if we want the telluroid to be free to have inclinations up to 60°. 相似文献
8.
This paper generalizes the Stokes formula from the spherical boundary surface to the ellipsoidal boundary surface. The resulting
solution (ellipsoidal geoidal height), consisting of two parts, i.e. the spherical geoidal height N
0 evaluated from Stokes's formula and the ellipsoidal correction N
1, makes the relative geoidal height error decrease from O(e
2) to O(e
4), which can be neglected for most practical purposes. The ellipsoidal correction N
1 is expressed as a sum of an integral about the spherical geoidal height N
0 and a simple analytical function of N
0 and the first three geopotential coefficients. The kernel function in the integral has the same degree of singularity at
the origin as the original Stokes function. A brief comparison among this and other solutions shows that this solution is
more effective than the solutions of Molodensky et al. and Moritz and, when the evaluation of the ellipsoidal correction N
1 is done in an area where the spherical geoidal height N
0 has already been evaluated, it is also more effective than the solution of Martinec and Grafarend.
Received: 27 January 1999 / Accepted: 4 October 1999 相似文献
9.
We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based
on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations.
The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity
field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal
energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite.
Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential
at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy
spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation
by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges
to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite
altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical
experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n
max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting
spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical
evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting
the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute
the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics
of satellite coverage, which could not be done using a global spherical harmonic representation.
The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577. 相似文献
10.
Construction of Green's function to the external Dirichlet boundary-value problem for the Laplace equation on an ellipsoid of revolution 总被引:1,自引:0,他引:1
Green's function to the external Dirichlet boundary-value problem for the Laplace equation with data distributed on an ellipsoid
of revolution has been constructed in a closed form. The ellipsoidal Poisson kernel describing the effect of the ellipticity
of the boundary on the solution of the investigated boundary-value problem has been expressed as a finite sum of elementary
functions which describe analytically the behaviour of the ellipsoidal Poisson kernel at the singular point ψ = 0. We have
shown that the degree of singularity of the ellipsoidal Poisson kernel in the vicinity of its singular point is of the same
degree as that of the original spherical Poisson kernel.
Received: 4 June 1996 / Accepted: 7 April 1997 相似文献
11.
This paper is devoted to the spherical and spheroidal harmonic expansion of the gravitational potential of the topographic masses in the most rigorous way. Such an expansion can be used to compute gravimetric topographic effects for geodetic and geophysical applications. It can also be used to augment a global gravity model to a much higher resolution of the gravitational potential of the topography. A formulation for a spherical harmonic expansion is developed without the spherical approximation. Then, formulas for the spheroidal harmonic expansion are derived. For the latter, Legendre’s functions of the first and second kinds with imaginary variable are expanded in Laurent series. They are then scaled into two real power series of the second eccentricity of the reference ellipsoid. Using these series, formulas for computing the spheroidal harmonic coefficients are reduced to surface harmonic analysis. Two numerical examples are presented. The first is a spherical harmonic expansion to degree and order 2700 by taking advantage of existing software. It demonstrates that rigorous spherical harmonic expansion is possible, but the computed potential on the geoid shows noticeable error pattern at Polar Regions due to the downward continuation from the bounding sphere to the geoid. The second numerical example is the spheroidal expansion to degree and order 180 for the exterior space. The power series of the second eccentricity of the reference ellipsoid is truncated at the eighth order leading to omission errors of 25 nm (RMS) for land areas, with extreme values around 0.5 mm to geoid height. The results show that the ellipsoidal correction is 1.65 m (RMS) over land areas, with maximum value of 13.19 m in the Andes. It shows also that the correction resembles the topography closely, implying that the ellipsoidal correction is rich in all frequencies of the gravity field and not only long wavelength as it is commonly assumed. 相似文献
12.
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. 相似文献
13.
Willi Freeden 《Journal of Geodesy》1980,54(1):1-20
Summary Let S be the (regular) boundary-surface of an exterior regionE
e
in Euclidean space ℜ3 (for instance: sphere, ellipsoid, geoid, earth's surface). Denote by {φn} a countable, linearly independent system of trial functions (e.g., solid spherical harmonics or certain singularity functions)
which are harmonic in some domain containingE
e
∪ S. It is the purpose of this paper to show that the restrictions {ϕn} of the functions {φn} onS form a closed system in the spaceC (S), i.e. any functionf, defined and continuous onS, can be approximated uniformly by a linear combination of the functions ϕn.
Consequences of this result are versions of Runge and Keldysh-Lavrentiev theorems adapted to the chosen system {φn} and the mathematical justification of the use of trial functions in numerical (especially: collocational) procedures. 相似文献
14.
15.
S. Ritter 《Journal of Geodesy》1998,72(2):101-106
The ellipsoidal Stokes problem is one of the basic boundary-value problems for the Laplace equation which arises in physical
geodesy. Up to now, geodecists have treated this and related problems with high-order series expansions of spherical and spheroidal
(ellipsoidal) harmonics. In view of increasing computational power and modern numerical techniques, boundary element methods
have become more and more popular in the last decade. This article demonstrates and investigates the nullfield method for
a class of Robin boundary-value problems. The ellipsoidal Stokes problem belongs to this class. An integral equation formulation
is achieved, and existence and uniqueness conditions are attained in view of the Fredholm alternative. Explicit expressions
for the eigenvalues and eigenfunctions for the boundary integral operator are provided.
Received: 22 October 1996 / Accepted: 4 August 1997 相似文献
16.
Integer least-squares theory for the GNSS compass 总被引:7,自引:2,他引:5
P. J. G. Teunissen 《Journal of Geodesy》2010,84(7):433-447
Global navigation satellite system (GNSS) carrier phase integer ambiguity resolution is the key to high-precision positioning
and attitude determination. In this contribution, we develop new integer least-squares (ILS) theory for the GNSS compass model,
together with efficient integer search strategies. It extends current unconstrained ILS theory to the nonlinearly constrained
case, an extension that is particularly suited for precise attitude determination. As opposed to current practice, our method
does proper justice to the a priori given information. The nonlinear baseline constraint is fully integrated into the ambiguity
objective function, thereby receiving a proper weighting in its minimization and providing guidance for the integer search.
Different search strategies are developed to compute exact and approximate solutions of the nonlinear constrained ILS problem.
Their applicability depends on the strength of the GNSS model and on the length of the baseline. Two of the presented search
strategies, a global and a local one, are based on the use of an ellipsoidal search space. This has the advantage that standard
methods can be applied. The global ellipsoidal search strategy is applicable to GNSS models of sufficient strength, while
the local ellipsoidal search strategy is applicable to models for which the baseline lengths are not too small. We also develop
search strategies for the most challenging case, namely when the curvature of the non-ellipsoidal ambiguity search space needs
to be taken into account. Two such strategies are presented, an approximate one and a rigorous, somewhat more complex, one.
The approximate one is applicable when the fixed baseline variance matrix is close to diagonal. Both methods make use of a
search and shrink strategy. The rigorous solution is efficiently obtained by means of a search and shrink strategy that uses
non-quadratic, but easy-to-evaluate, bounding functions of the ambiguity objective function. The theory presented is generally
valid and it is not restricted to any particular GNSS or combination of GNSSs. Its general applicability also applies to the
measurement scenarios (e.g. single-epoch vs. multi-epoch, or single-frequency vs. multi-frequency). In particular it is applicable
to the most challenging case of unaided, single frequency, single epoch GNSS attitude determination. The success rate performance
of the different methods is also illustrated. 相似文献
17.
The Cartesian moments of the mass density of a gravitating body and the spherical harmonic coefficients of its gravitational
field are related in a peculiar way. In particular, the products of inertia can be expressed by the spherical harmonic coefficients
of the gravitational potential as was derived by MacCullagh for a rigid body. Here the MacCullagh formulae are extended to
a deformable body which is restricted to radial symmetry in order to apply the Love–Shida hypothesis. The mass conservation
law allows a representation of the incremental mass density by the respective excitation function. A representation of an
arbitrary Cartesian monome is always possible by sums of solid spherical harmonics multiplied by powers of the radius. Introducing
these representations into the definition of the Cartesian moments, an extension of the MacCullagh formulae is obtained. In
particular, for excitation functions with a vanishing harmonic coefficient of degree zero, the (diagonal) incremental moments
of inertia also can be represented by the excitation coefficients. Four types of excitation functions are considered, namely:
(1) tidal excitation; (2) loading potential; (3) centrifugal potential; and (4) transverse surface stress. One application
of the results could be model computation of the length-of-day variations and polar motion, which depend on the moments of
inertia.
Received: 27 July 1999 / Accepted: 24 May 2000 相似文献
18.
An effort is made at developing a theory of map readability, defined as the process of the user's representation of the information of the map in his/her own mind. This can be estimated quantitatively by surrogate measures, which include the speed of map comprehension and accuracy of map interpretation. Levels of psychological representation of map information are incorporated to determine at what stage particular aspects of map knowledge are understood. Thus what readability entails, and the approaches used to measure it (what types of things are comprehended, how fast, and how accurately) will depend upon the level of representation. Translated from: Vestnik Leningradskogo Universiteta, seriya 7 [geologiya, geografiya], 1988, No. 1, pp. 32-37. 相似文献
19.
The products of Wuhan University with 5-min sampling are used to analyze the characteristics of BeiDou satellite clocks. Two nanoseconds root-mean-square (RMS) variations are obtained for 1-day quadratic fits in the sub-daily region. The relativistic effects of BDS clocks are also studied. General relativity predicts that linear variation of the semimajor axes of geostationary and inclined geosynchronous satellites causes a quadratic clock drift with a magnitude at the 10?16/day level. The observed drift is higher than what general relativity theory would produce. Several periodic terms are found in the satellite clock variations through spectrum analysis. In order to identify the origin of the BDS clock harmonics, a correlation analysis between the period or amplitude of the harmonics and properties of the satellite orbits is performed. It is found that the period of the harmonics is not exactly equal to the orbit period, but rather the ratio of the orbit period to clock period is almost the same as that of a sidereal day to solar day. The BDS clocks obey white frequency noise statistics for intervals from 300 s to several thousands seconds. For intervals greater than 10,000 s, all the BDS satellites display more complex, non-power-law behavior due to the effects of periodic clock variations. 相似文献
20.
Existing research on DEM vertical accuracy assessment uses mainly statistical methods, in particular variance and RMSE which are both based on the error propagation theory in statistics. This article demonstrates that error propagation theory is not applicable because the critical assumption behind it cannot be satisfied. In fact, the non‐random, non‐normal, and non‐stationary nature of DEM error makes it very challenging to apply statistical methods. This article presents approximation theory as a new methodology and illustrates its application to DEMs created by linear interpolation using contour lines as the source data. Applying approximation theory, a DEM's accuracy is determined by the largest error of any point (not samples) in the entire study area. The error at a point is bounded by max(|δnode|+M2h2/8) where |δnode| is the error in the source data used to interpolate the point, M2 is the maximum norm of the second‐order derivative which can be interpreted as curvature, and h is the length of the line on which linear interpolation is conducted. The article explains how to compute each term and illustrates how this new methodology based on approximation theory effectively facilitates DEM accuracy assessment and quality control. 相似文献