The harmonic representation of the disturbing potential

Summary The alternative harmonic representations of the disturbing potential, correct to the order of the flattening, are examined and an example is given where the incorrect use of a spherical harmonic expansion can give rise to fallacious results. The correct usage of the spherical harmonic expansion for the disturbing potential is given in the solution of the general surface integral to define the indirect effect in the case of the non-regularised geoid.     

不同扰动位泛函间的积分变换广义核函数

基于物理大地测量边值问题的解,利用一阶边界算子定义,推导重力异常Δg、单层密度μ、大地水准面高N,垂线偏差ε、扰动重力δg等扰动场元的解。利用球谐函数的正交特性,通过对核函数的算子运算,可以得到上述扰动场元的有关逆变换公式。相对经典物理大地测量公式应用的边界面条件,笔者将含有因子r的对应扰动场元反演关系的公式称为广义积分公式。针对常用的重力异常Δg、大地水准面高N,垂线偏差ε、扰动重力δg计算,重点分析它们之间的变换关系,给出利用某个选定扰动场元计算其他扰动场元的广义积分公式。同时,通过对积分边界面的讨论,分析经典公式与广义积分公式的差异和联系。最后,给出所有外部扰动场元与核函数映射的关系表。     

A formula for computing the gravity disturbance from the second radial derivative of the disturbing potential

 A formula for computing the gravity disturbance and gravity anomaly from the second radial derivative of the disturbing potential is derived in detail using the basic differential equation with spherical approximation in physical geodesy and the modified Poisson integral formula. The derived integral in the space domain, expressed by a spherical geometric quantity, is then converted to a convolution form in the local planar rectangular coordinate system tangent to the geoid at the computing point, and the corresponding spectral formulae of 1-D FFT and 2-D FFT are presented for numerical computation. Received: 27 December 2000 / Accepted: 3 September 2001     

Local relationships between the disturbing density, the disturbing potential and the disturbing gravity of the Earth's gravity field

A function having some properties of a wavelet and being harmonic around a given point in R ³ is defined, and three models showing the local relationships between the disturbing density, the disturbing potential and the disturbing gravity are established by using the function as the kernel function of the integrals in the models. The local relationship has two meanings. One is that we can evaluate with a high accuracy the integrals in the models by using mainly high-accuracy and high-resolution data in a local area. The other is that we can obtain a stable solution with high resolution when we invert the integrals in the models because of the rapid decrease of the kernel function of the integrals. As a result, with these models we evaluate one quantity with high resolution, in a band limited by the maximum degree of a set of geopotential coefficients or by the resolution (spacing) of the local data, from another quantity (or quantities) in a local area, and the resulting solution is stable. Received: 6 April 1998 / Accepted: 16 June 1999     

Accuracy Analysis of Geopotential Coefficients Recovered from In-situ Disturbing Potential by Energy Conservation Method

The characteristics of the normal equation created in recovering the Earth gravity model (EGM) by least-squares(LS) adjustment from the in-situ disturbing potential is discussed in detail. It can be co...     

A new modelling method for evaluating external disturbing potential based on theory of unified representation of gravitational field

1　ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｖｅｎｔｉｏｎａｌｅｖａｌｕａｔｉｎｇｍｏｄｅｌｓｏｆｅｘｔｅｒｎａｌｇｒａｖｉ ｔａｔｉｏｎａｌｆｉｅｌｄｏｆｔｈｅｅａｒｔｈａｒｅｎｏｔｈｉｎｇｅｌｓｅｔｈａｎａｋｉｎｄｏｆｒｅａｌｉｚａｔｉｏｎｓｏｆｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｓｔｏｇｅｏｄｅ ｔｉｃｂｏｕｎｄａｒｙｖａｌｕｅ ｐｒｏｂｌｅｍ .Ｉｎｔｈｅｓｅｎｓｅｏｆｌｉｎ ｅａｒｉｚａｔｉｏｎ ,ｇｅｏｄｅｔｉｃｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍｓａｒｅｒｅ ｄｕｃｅｄｔｏｅｘｔｅｒｎａｌｂｏｕｎｄ…     

A new modelling method for evaluating external disturbing potential based on theory of unified representation of gravitational field

For a special use a new modelling method of evaluating external disturbing potential is presented in this paper. Being different from classical methods in physical geodesy this method is grounded upon the theory of unified representation of gravitational field. The models created in this way are particularly satisfactory for a high-speed computation of gravitational field in low altitude because they take account of topographic effects and have their kernel functions with simple structure and weak singularity.     

四种改进积分法的低空扰动引力计算

针对Stokes积分方法计算扰动引力中计算点从空中趋近地面时存在积分奇异和不连续的问题,该文提出了去中央奇异点法、奇异点积分值修正法、中央格网加密算法和改进积分式法4种改进Stokes积分的计算公式,并进行了实验计算。计算结果表明:近地空间范围内,4种改进算法都能在一定程度上改进原始积分的奇异性问题;相同条件下,奇异点积分值修正法和改进积分式法计算精度最高,适宜于低空计算;改进积分式法通过理论推导,得到了从球外部到球面统一、连续且无奇异的改进Stokes积分公式,理论严谨。     

Poisson重力边值问题

提出了Poisson重力边值问题,即关于扰动位的Poisson方程的Stokes问题和Neumann问题。作为导引,先研究Poisson方程的Dirichlet问题．再分别引入一种辅助函数,将Stokes问题和Neumann问题改化为Dirichlet问题,从而立即得到它们的积分解。最终解式表现为两部分叠加：一部分仅与边界观测相关,另一部分为对地形测量的响应,为研究地形测量对外部重力场和大地水准面的精化提供新的途径。     

Helmert扰动位及其积分核函数的椭球实用公式

借助以地心参考椭球面为边界面的第二大地边值问题的理论,基于Helmert空间的Neumann边值条件,给定Helmert扰动位的椭球解表达式,并详细推导第二类勒让德函数及其导数的递推关系、Helmert扰动位函数的椭球积分解以及类椭球Hotine积分核函数的实用计算公式,便于后续椭球域第二大地边值问题的实际研究。     

Adaptive filtering of GOCE-derived gravity gradients of the disturbing potential in the context of the space-wise approach

Filtering and signal processing techniques have been widely used in the processing of satellite gravity observations to reduce measurement noise and correlation errors. The parameters and types of filters used depend on the statistical and spectral properties of the signal under investigation. Filtering is usually applied in a non-real-time environment. The present work focuses on the implementation of an adaptive filtering technique to process satellite gravity gradiometry data for gravity field modeling. Adaptive filtering algorithms are commonly used in communication systems, noise and echo cancellation, and biomedical applications. Two independent studies have been performed to introduce adaptive signal processing techniques and test the performance of the least mean-squared (LMS) adaptive algorithm for filtering satellite measurements obtained by the gravity field and steady-state ocean circulation explorer (GOCE) mission. In the first study, a Monte Carlo simulation is performed in order to gain insights about the implementation of the LMS algorithm on data with spectral behavior close to that of real GOCE data. In the second study, the LMS algorithm is implemented on real GOCE data. Experiments are also performed to determine suitable filtering parameters. Only the four accurate components of the full GOCE gravity gradient tensor of the disturbing potential are used. The characteristics of the filtered gravity gradients are examined in the time and spectral domain. The obtained filtered GOCE gravity gradients show an agreement of 63–84 mEötvös (depending on the gravity gradient component), in terms of RMS error, when compared to the gravity gradients derived from the EGM2008 geopotential model. Spectral-domain analysis of the filtered gradients shows that the adaptive filters slightly suppress frequencies in the bandwidth of approximately 10–30 mHz. The limitations of the adaptive LMS algorithm are also discussed. The tested filtering algorithm can be connected to and employed in the first computational steps of the space-wise approach, where a time-wise Wiener filter is applied at the first stage of GOCE gravity gradient filtering. The results of this work can be extended to using other adaptive filtering algorithms, such as the recursive least-squares and recursive least-squares lattice filters.     

快速确定扰动引力的广域多项式逼近方法的模拟实验

飞行器在飞行过程中时刻受到地球重力场的作用,地球外部空间扰动引力会对飞行器的运行轨迹和姿态等方面产生影响,因此,研究地球重力场信息并实时快速的计算出飞行器运行处的扰动引力,对飞行器轨迹和姿态等方面的调节具有重要的作用。利用模拟数据,本文基于广域多项式逼近方法来快速逼近扰动引力,该法在单元边界点上满足插值条件,而又利用插值单元周围节点信息,实现拟合过程。试验表明,这一方法在高度方向的步长为30km的条件下可达到1mgal的计算精度。     

应用格林积分直接以地面边值确定外部扰动重力场

以地面为边界的Molodensky问题通常得到的是级数形式的解,高阶项体现了地面边值到某一光滑面上的改正,在应用中不仅遇到计算的复杂性、稳定性问题,也存在对数据密集要求的困难。本文从推求外部扰动重力场的应用出发,将格林公式用于数字地形面,在忽略水平分量的积分影响的情况下,得到以地面重力异常和高程异常差为边值确定外部扰动位的表达式,其核函数分别为距离倒数和Poisson核。该方法不需要对地面数据进行延拓处理,且核函数形式简洁,适于外部扰动重力场的随机计算。     

Integral transformations of deflections of the vertical onto satellite-to-satellite tracking and gradiometric data

With the advent of geodetic satellite missions mapping almost globally the Earth’s gravitational field, new methods and theoretical approaches have been developed and investigated to fully exploit the potential of their new observables. Besides estimating values of numerical coefficients in harmonic series of the gravitational potential, new applications emerged such as data validation and combination. In this contribution, new integral transformations are presented which transform principal components of the terrestrial deflection of the vertical onto disturbing satellite-to-satellite tracking and gradiometric data at altitude. Using spherical approximation, necessary integral kernel functions are derived in both spectral and closed forms. The behaviour of isotropic kernel functions is studied and the new integral transformations are tested in a closed-loop simulation using synthetic terrestrial and satellite data synthesized from a global gravitational model. New integral transformations can be used for data validation and combination purposes.     

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.     

Truncation error formulae for the disturbing gravity vector

Recurrence relations have been derived for truncation error coefficients of the extended Stokes' function and its partial derivatives required in the computation of the disturbing gravity vector at any elevation above the earth's surface. The corresponding formulae, the example of values of the truncation error coefficients for H=30.1 km and ψ₀=30^∘ and the estimations of truncation error are given in this article. Received: 26 January 1996 / Accepted: 11 June 1997     

A comparison of methods for the inversion of airborne gravity data

Four integral-based methods for the inversion of gravity disturbances, derived from airborne gravity measurements, into the disturbing potential on the Bjerhammar sphere and the Earths surface are investigated and compared with least-squares (LS) collocation. The performance of the methods is numerically investigated using noise-free and noisy observations, which have been generated using a synthetic gravity field model. It is found that advanced interpolation of gravity disturbances at the nodes of higher-order numerical integration formulas significantly improves the performance of the integral-based methods. This is preferable to the commonly used one-point composed Newton–Cotes integration formulas, which intrinsically imply a piecewise constant interpolation over a patch centered at the observation point. It is shown that the investigated methods behave similarly for noise-free observations, but differently for noisy observations. The best results in terms of root-mean-square (RMS) height-anomaly errors are obtained when the gravity disturbances are first downward continued (inverse Poisson integral) and then transformed into potential values (Hotine integral). The latter has a strong smoothing effect, which damps high-frequency errors inherent in the downward-continued gravity disturbances. An integral method based on the single-layer representation of the disturbing potential shows a similar performance. This representation has the advantage that it can be used directly on surfaces with non-spherical geometry, whereas classical integral-based methods require an additional step if gravity field functionals have to be computed on non-spherical geometries. It is shown that defining the single-layer density on the Bjerhammar sphere gives results with the same quality as obtained when using the Earths topography as support for the single-layer density. A comparison of the four integral-based methods with LS collocation shows that the latter method performs slightly better in terms of RMS height-anomaly errors.     

球外部空间的Neumann逆问题及逆Stokes公式

基于修改的Poisson积分,首先给出了球面扰动位向上延拓的积分表达式。在此基础上,由微分原理得出了球外部空间Neumann逆问题的解式,利用物理大地测量学的基本微分方程,导出了球外部空间的逆Stokes公式,并对这两类积分公式的核函数进行了讨论。     

远区重力数据对扰动位径向二阶梯度的影响

提出了利用地面重力异常数据计算地面扰动位径向二阶梯度时，远距离带区影响的误差公式，结果以级数形式给出，以其数学模型为基础，利用全球重力位模型GPM98，计算确定了不同平均异常网格大小情况下截断半径取一系列值时的远区影响，在将卫星重力梯度和航空重力梯度观测数据与地面重力数据进行联合处理时，文中提出了远区影响的数量关系，对于确定合理的截断半径大小具有一定的参考价值。     

The geopotential from gravity measurements,levelling data and satellite results

The geodetic boundary value problem is formulated which uses as boundary values the differences between the geopotential of points at the surface of the continents and the potential of the geoid. These differences are computed by gravity measurements and levelling data. In addition, the shape of the geoid over the oceans is assumed to be known from satellite altimetry and the shape of the continents from satellite results together with three-dimensional triangulation. The boundary value problem thus formulated is equivalent to Dirichlet's exterior problem except for the unknown potential of the geoid. This constant is determined by an integral equation for the normal derivative of the gravitational potential which results from the first derivative of Green's fundamental formula. The general solution, which exists, of the integral equation gives besides the potential of the geoid the solution of the geodetic boundary value problem. In addition approximate solutions for a spherical surface of the earth are derived.