Solving Three Types of Satellite Gravity Gradient Boundary Value Problems by Least-Squares

The principle and method for solving three types of satellite gravity gradient boundary value problems by least-squares are discussed in detail. Also, kernel function expressions of the least-squares solution of three geodetic boundary value problems with the observations {Γ zz },{Γ xz , Γ yz} and {Γ xx -Γ yy ,2 Γxy}are presented. From the results of recovering gravity field using simulated gravity gradient tensor data, we can draw a conclusion that satellite gravity gradient integral formulas derived from least-squares are valid and rigorous for recovering the gravity field.     

最小二乘法求解三类卫星重力梯度边值问题的研究

研究了最小二乘法求解3类卫星重力梯度边值问题的理论和方法，给出了3类梯度观测值{Гzz}、{Гxz、Гyz}和{Гxx-Гyy,2Гxy}对应边值问题解的核函数严密表达式。模拟试算结果表明，最小二乘法求解的卫星重力梯度积分公式用于恢复地球重力场是有效而严密的。     

联合高低卫-卫跟踪和卫星重力梯度数据恢复地球重力场的谱组合法

GOCE采用的高低卫-卫跟踪和卫星重力梯度测量技术在恢复重力场方面各有所长并互为补充,如何有效利用这两类观测数据最优确定地球重力场是GOCE重力场反演的关键问题。本文研究了联合高低卫-卫跟踪和卫星重力梯度数据恢复地球重力场的最小二乘谱组合法,基于球谐分析方法推导并建立了卫星轨道面扰动位T和径向重力梯度Tzz、以及扰动位T和重力梯度分量组合{Tzz-Txx-Tyy}的谱组合计算模型与误差估计公式。数值模拟结果表明,谱组合计算模型可以有效顾及各类数据的精度和频谱特性进行最优联合求解。采用61天GOCE实测数据反演的两个180阶次地球重力场模型WHU_GOCE_SC01S（扰动位和径向重力梯度数据求解）和WHU_GOCE_SC02S（扰动位和重力梯度分量组合数据求解）,结果显示后者精度优于前者,并且它们的整体精度优于GOCE时域解,而与GOCE空域解的精度接近,验证了谱组合法的可行性与有效性。     

利用重力异常数据构建重力梯度场的方法比较

为实现大范围、高精度基准重力梯度数据库的构建,考虑到重力梯度场对地形质量的敏感效应,一般利用恒密度数字高程模型来求取重力梯度值,从而忽略了地形密度变化以及水准面以下密度异常对重力梯度的影响。根据重力位理论中求解边值问题的数值应用方法,直接利用重力异常数据求取重力梯度场,弥补了密度变化和密度异常在重力梯度上的反映。根据模型算例和实测重力异常数据求取了剖面重力梯度值,结果表明,限于重力数据空间分辨率的影响,利用重力异常数据可恢复中长波段重力梯度场。该方法与地形数据求取重力梯度和卫星重力梯度测量等方法技术相结合,对重力梯度数据库的建设具有实际应用价值。     

卫星重力梯度边值问题的点质量调和分析

研究并建立由全球重力梯度复组合分量及全张量解算全球点质量模型的基本方程,进一步推导得到基于卫星重力梯度的单定边值问题和超定边值问题的点质量调和分析解.通过采用分块循环矩阵分解大型线性方程组的方法,实现点质量调和分析解的稳定解算.最后运用EGM2008模型进行模拟数值计算,验证卫星重力梯度边值问题的点质量调和分析法的有效...     

解混合边值问题直接计算位系数的变分原理

本文提出了利用变分法解混合边值问题直接计算位系数的原理。根据这一原理可解第一、第二和第三边值问题的混合边值问题直接求得位系数。利用这一原理可较简单地联合利用经典重力测量(即重力点的平面位置由天文或三角测量确定,高程由水准或三角高程确定)、卫星重力测量(即利用卫星定位技术确定重力点的平面位置和大地高)以及卫星测高数据研究地球的重力场。     

卫星重力径向梯度数据的最小二乘配置调和分析

本文深入研究了利用卫星重力梯度径向分量确定地球引力场位系数的最小二乘配置(LSC)调和分析方法。首先论述了最小二乘配置法的原理,推导了扰动引力梯度观测量与球谐系数之间的协方差和自协方差矩阵,在扰动引力梯度观测数据为等经差规则网格数据的情况下,引力位与扰动引力梯度之间的协方差矩阵具有分块Toeplitz循环阵的结构,有效的利用FFT变换技术将其降阶;研究利用截断奇异值分解法（TSVD）解决协方差阵的病态性问题;最后得到了引力梯度径向分量的最小二乘配置调和分析的完整计算公式。模拟试算结果表明,基于TSVD的最小二乘配置调和分析方法,能够以较高的精度还原全球重力场,验证了本文算法的有效性和实用性。     

块对角最小二乘方法在确定全球重力场模型中的应用

分析研究了最小二乘求解重力场模型的六种位系数排列方式下的块对角形态及三种不同条件下的块对角近似：BD-1、BD-2、BD-3,探讨了块对角最小二乘方法在联合早期卫星重力场模型和最新GRACE-only模型中的应用,通过实验计算,结果表明,块对角最小二乘方法较之于积分方法,能更好的提高所恢复模型的精度,说明在卫星重力飞速发展、地面重力数据不断完善的今天,块对角最小二乘法在超高阶地球重力场模型构建方面的优势逐渐突出。     

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.     

Le theoreme de conservation de la courbure et de la torsion

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.      

关于卫星重力梯度边值问题的准解的若干注记

研究了卫星重力梯度边值问题的准解的具体计算方法，利用地球重力场模型ＷＤＭ９４模拟的卫星重力梯度数据进行试算，验证了准解模型的有效性，并获得一些重要结论     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

The free versus fixed geodetic boundary value problem for different combinations of geodetic observables

Various formulations of the geodetic fixed and free boundary value problem are presented, depending upon the type of boundary data. For the free problem, boundary data of type astronomical latitude, astronomical longitude and a pair of the triplet potential, zero and first-order vertical gradient of gravity are presupposed. For the fixed problem, either the potential or gravity or the vertical gradient of gravity is assumed to be given on the boundary. The potential and its derivatives on the boundary surface are linearized with respect to a reference potential and a reference surface by Taylor expansion. The Eulerian and Lagrangean concepts of a perturbation theory of the nonlinear geodetic boundary value problem are reviewed. Finally the boundary value problems are solved by Hilbert space techniques leading to new generalized Stokes and Hotine functions. Reduced Stokes and Hotine functions are recommended for numerical reasons. For the case of a boundary surface representing the topography a base representation of the solution is achieved by solving an infinite dimensional system of equations. This system of equations is obtained by means of the product-sum-formula for scalar surface spherical harmonics with Wigner 3j-coefficients.     

The free versus fixed geodetic boundary value problem for different combinations of geodetic observables

Various formulations of the geodetic fixed and free boundary value problem are presented, depending upon the type of boundary data. For the free problem, boundary data of type astronomical latitude, astronomical longitude and a pair of the triplet potential, zero and first-order vertical gradient of gravity are presupposed. For the fixed problem, either the potential or gravity or the vertical gradient of gravity is assumed to be given on the boundary. The potential and its derivatives on the boundary surface are linearized with respect to a reference potential and a reference surface by Taylor expansion. The Eulerian and Lagrangean concepts of a perturbation theory of the nonlinear geodetic boundary value problem are reviewed. Finally the boundary value problems are solved by Hilbert space techniques leading to new generalized Stokes and Hotine functions. Reduced Stokes and Hotine functions are recommended for numerical reasons. For the case of a boundary surface representing the topography a base representation of the solution is achieved by solving an infinite dimensional system of equations. This system of equations is obtained by means of the product-sum-formula for scalar surface spherical harmonics with Wigner 3j-coefficients.     

利用GOCE模拟观测反演重力场的Torus法

在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61dGOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022mm和0.099mm。在沿轨卫星引力梯度模拟数据中加入5mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58cm和1.45cm,累积误差最大值分别为6.37cm和5.55cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。     

Preprocessing of gravity gradients at the GOCE high-level processing facility

One of the products derived from the gravity field and steady-state ocean circulation explorer (GOCE) observations are the gravity gradients. These gravity gradients are provided in the gradiometer reference frame (GRF) and are calibrated in-flight using satellite shaking and star sensor data. To use these gravity gradients for application in Earth scienes and gravity field analysis, additional preprocessing needs to be done, including corrections for temporal gravity field signals to isolate the static gravity field part, screening for outliers, calibration by comparison with existing external gravity field information and error assessment. The temporal gravity gradient corrections consist of tidal and nontidal corrections. These are all generally below the gravity gradient error level, which is predicted to show a 1/f behaviour for low frequencies. In the outlier detection, the 1/f error is compensated for by subtracting a local median from the data, while the data error is assessed using the median absolute deviation. The local median acts as a high-pass filter and it is robust as is the median absolute deviation. Three different methods have been implemented for the calibration of the gravity gradients. All three methods use a high-pass filter to compensate for the 1/f gravity gradient error. The baseline method uses state-of-the-art global gravity field models and the most accurate results are obtained if star sensor misalignments are estimated along with the calibration parameters. A second calibration method uses GOCE GPS data to estimate a low-degree gravity field model as well as gravity gradient scale factors. Both methods allow to estimate gravity gradient scale factors down to the 10⁻³ level. The third calibration method uses high accurate terrestrial gravity data in selected regions to validate the gravity gradient scale factors, focussing on the measurement band. Gravity gradient scale factors may be estimated down to the 10⁻² level with this method.     

地球重力场、重力、重力梯度在三维直角坐标系中的表达式

现代卫星重力测量主要利用星载GPS接收机、加速度计、星载测距仪等来确定重力卫星的轨道 ,削弱非保守力的干扰 ,由此根据卫星的位置、速度及其变率来确定地球重力场。而上述GPS等星载仪器所提供的数据 ,包括卫星轨道坐标及其速率、扰动加速度、星间距离及其变率 ,都是以三维直角坐标 (x ,y ,z)的形式表示的 ,因此 ,地球重力场、重力和重力梯度在三维直角坐标系中的表达式在卫星重力解算中具有实际意义     

国际卫星重力梯度测量计划研究进展

本文首先阐述了重力梯度测量原理、从20世纪初到21世纪初重力梯度仪的研究历程、卫星重力梯度仪(静电悬浮重力梯度仪、超导重力梯度仪和量子重力梯度仪)的技术特征以及卫星重力梯度测量的特点;其次,介绍了基于卫星重力梯度技术恢复250阶GOCE地球重力场以及论证首先开展一维径向重力梯度仪的研制进而恢复高精度和高空间解析度中高频地球重力场可行性方面的研究进展;最后,建议我国尽早开展基于时空域混合法解算中高频地球重力场和卫星重力梯度测量系统误差分析的预先研究。     

Geodetic methods to determine the relativistic redshift at the level of 10$$^{}$$ in the context of international timescales: a review and practical results

The frequency stability and uncertainty of the latest generation of optical atomic clocks is now approaching the one part in \(10^{18}\) level. Comparisons between earthbound clocks at rest must account for the relativistic redshift of the clock frequencies, which is proportional to the corresponding gravity (gravitational plus centrifugal) potential difference. For contributions to international timescales, the relativistic redshift correction must be computed with respect to a conventional zero potential value in order to be consistent with the definition of Terrestrial Time. To benefit fully from the uncertainty of the optical clocks, the gravity potential must be determined with an accuracy of about \(0.1\,\hbox {m}^{2}\,\hbox {s}^{-2}\), equivalent to about 0.01 m in height. This contribution focuses on the static part of the gravity field, assuming that temporal variations are accounted for separately by appropriate reductions. Two geodetic approaches are investigated for the derivation of gravity potential values: geometric levelling and the Global Navigation Satellite Systems (GNSS)/geoid approach. Geometric levelling gives potential differences with millimetre uncertainty over shorter distances (several kilometres), but is susceptible to systematic errors at the decimetre level over large distances. The GNSS/geoid approach gives absolute gravity potential values, but with an uncertainty corresponding to about 2 cm in height. For large distances, the GNSS/geoid approach should therefore be better than geometric levelling. This is demonstrated by the results from practical investigations related to three clock sites in Germany and one in France. The estimated uncertainty for the relativistic redshift correction at each site is about \(2 \times 10^{-18}\).