基于重力卫星数据监测地表质量变化的三维点质量模型法

利用卫星重力测量手段监测全球质量变化取得了巨大成功,本文基于牛顿万有引力定律在三维空间直角坐标系中导出利用重力卫星观测数据监测全球质量变化的三维点质量模型法,该方法可直接利用重力卫星的轨道和星间观测数据或时变重力场模型计算全球质量变化,由于利用卫星观测数据计算地表质量变化的向下延拓过程以及观测数据噪声的影响,需要采用合适的空间约束方程或正则化技术对解算结果进行约束或平滑处理.利用合成全球质量变化模型模拟一个月的GRACE双星轨道和星间距离变率数据计算全球质量变化,对三维点质量模型法进行分析验证,采用零阶Tikhonov正则化技术处理病态问题.结果表明,三维点质量模型法可有效用于重力卫星观测数据监测全球质量变化,为利用重力卫星观测数据监测全球质量变化提供一种可选的途径.     

The Earth's gravity field

The Earth's gravity field can be determined from gravity measurements made on the surface of the Earth, and through the analysis of the motion of Earth satellites. Gravity data can be used to solve the boundary value problem of gravimetric geodesy in various ways, from the classical formulation using a geoid to the concept of a reference surface interior to the masses of the Earth to a statistical method. We now have gravity information for 1⁰ data blocks over 46% of the Earth's surface and more than several million point measurements available.Satellite observations such as range, range-rate, and optical data have been analyzed to determine potential coefficients used to describe the Earth's gravitational potential field. Coefficients, in a spherical harmonic expansion to degree 12, can be determined from satellite data alone, and to at least degree 20 when the satellite data is combined with surface gravity material. Recent solutions for potential coefficients agree well to degree 4, but with increasing disagreement at higher degrees.     

Geoid Determination through Ellipsoidal Stokes Boundary-Value Problem

Martinec and Grafarend (1997) have shown how the construction of Green's function in the Stokes boundary-value problem with gravity data distributed on an ellipsoid of revolution is approached in the O(e ₀ ² )-approximation. They have also expressed the ellipsoidal Stokes function describing the effect of ellipticity of the boundary as a finite sum of elementary functions. We present an effective method of avoiding the singularity of spherical and the ellipsoidal Stokes functions, and also an analytical expression for the ellipsoidal Stokes integral around the computational point suitable for numerical solution. We give the numerical results of solving the ellipsoidal Stokes boundary-value problem and their difference with respect to the spherical Stoke boundary-value problem.     

利用余弦变换计算重力异常的向上延拓

利用余弦变换计算重力异常的向上延拓是一种新方法.根据余弦变换的基本性质,推导了二度、三度体异常向上延拓余弦变换谱理论公式,采用离散余弦变换实现了该法的数值计算;研究了无限长水平圆柱体的补偿因子中主频段的特性,给出了二度体的线性补偿方式;补偿后的理论模型异常向上延拓具有较高的计算精度,除边部几个数据因数据的离散和有限截断使误差较大外（最大误差为6.23%）,其余数据的误差均在1%以内,理论值和计算值曲线基本重合.这说明,与Fourier变换相比,离散余弦变换在数值计算中,受非周期性深度因子的影响小,补偿方式易于选择,其计算方法优于Fourier变换.     

积分迭代法的正则性分析及其最优步长的选择

位场积分迭代法是一种解决位场大数据量、大深度向下延拓的有效方法.本文基于Kirsch正则化子理论,推导了积分迭代法对应的正则化滤子函数,并证明积分迭代法为一种求解位场向下延拓不适定反问题的正则化方法.针对积分迭代法迭代步长固定、迭代次数较多影响收敛速度的问题,提出该迭代法最优迭代步长的选择原理.理论模型和实测数据对比分...     

反演地表质量变化的附有空间约束的三维加速度点质量模型法

重力卫星可以在相同误差尺度下对全球质量变化进行连续重复观测,并在近十余年来取得了巨大成功,探索重力卫星数据精化处理方法和相关应用研究具有重要意义.本文基于三维加速度点质量模型法的基本原理,进一步发展建立了时变重力场模型球谐位系数的变化和地面点质量变化的关系,可有效考虑地表质量变化导致的负荷形变的影响;引入等权形式、线性形式、指数形式和高斯形式的空间约束方法处理南北条带噪声和向下延拓导致的病态问题,并与零阶Tikhonov正则化方法进行对比分析.采用模拟数据和一个月的实测GRACE时变重力场模型计算全球质量变化,对三维加速度点质量模型法和几种空间约束方法进行对比分析验证.计算结果表明,对于3°等面积的全球格网质量点,高斯和指数形式空间约束方法的最优相关距离约为500km,等权和线性形式空间约束方法的最优相关距离约为600km,各方法均可有效处理条带噪声的影响,四种空间约束方法的计算效果优于零阶Tikhonov正则化方法,本文的相关方法为进一步利用三维加速度点质量模型法监测全球质量变化提供了借鉴.     

引力场球谐综合的Walsh变换方法

提出直接在序率域内用Walsh变换实现引力场球谐综合的问题。给出球谐函数展开式的Walsh变换及快速算法,讨论了Walsh变换和Walsh-Fourier变换、Fourier变换之间的差异,分析了用地球重力场模型OSU81的位系数作出的Walsh变换和Fourier变换的结果。研究表明:Walsh变换与Walsh-Fourier变换、Fourier变换对应向量在数量方面的差值通常都小于士10~(-5);对于给定的阶数和飞行高度,3种方法求得的球谐综合值总是完全一致的;1°×1°等网格数据和Walsh函数形状相近。在重力场研究中Walsh级数会比Fourier级数收敛得更快;Walsh变换在计算速度、计算准确度、数据储存、收敛速度和方法简单方面都好于Fourier变换。     

重力场向下延拓Milne法

重力场向下延拓能够突出局部和浅部的异常信息,分离叠加的异常特征.但是向下延拓通常具有过程不稳定、下延深度小、结果不准确等问题.针对向下延拓所存在的不足,本文利用重力场及其垂向一阶导数,基于辛普森(Simpson)求积公式,推导出重力场向下延拓米尔尼(Milne)公式.将本文向下延拓方法应用于模型数据,向下延拓模型结果及误差曲线表明,相对于向下延拓快速傅里叶变换(FFT)法和积分迭代法,向下延拓Milne法的深度更大,相对误差更小;相对模型值,向下延拓Milne法能够获得稳定且准确的结果.对加拿大乃查科(Nechako)盆地地区实测航空重力数据进行本文方法向下延拓验证,处理结果表明,相对于实测异常,本文方法向下延拓结果能够很好还原实测数据,并且在进一步向下延拓中反映原始异常的趋势,增强局部和细小异常信息.     

Numerical simulations for the non-linear Molodensky problem

We present a boundary element method to compute numerical approximations to the non-linear Molodensky problem, which reconstructs the surface of the Earth from the gravitational potential and the gravity vector. Our solution procedure solves a sequence of exterior oblique Robin problems and is based on a Nash-Hörmander iteration. We apply smoothing with the heat equation to overcome a loss of derivatives in the surface update. Numerical results show the error between the approximation and the exact solution in a model problem.     

Direct BEM for high-resolution global gravity field modelling

The paper presents a high-resolution global gravity field modelling by the boundary element method (BEM). A direct BEM formulation for the Laplace equation is applied to get a numerical solution of the linearized fixed gravimetric boundary-value problem. The numerical scheme uses the collocation method with linear basis functions. It involves a discretization of the complicated Earth’s surface, which is considered as a fixed boundary. Here 3D positions of collocation points are simulated from the DNSC08 mean sea surface at oceans and from the SRTM30PLUS_V5.0 global topography model added to EGM96 on lands. High-performance computations together with an elimination of the far zones’ interactions allow a very refined integration over the all Earth’s surface with a resolution up to 0.1 deg. Inaccuracy of the approximate coarse solutions used for the elimination of the far zones’ interactions leads to a long-wavelength error surface included in the obtained numerical solution. This paper introduces an iterative procedure how to reduce such long-wavelength error surface. Surface gravity disturbances as oblique derivative boundary conditions are generated from the EGM2008 geopotential model. Numerical experiments demonstrate how the iterative procedure tends to the final numerical solutions that are converging to EGM2008. Finally the input surface gravity disturbances at oceans are replaced by real data obtained from the DNSC08 altimetryderived gravity data. The ITG-GRACE03S satellite geopotential model up to degree 180 is used to eliminate far zones’ interactions. The final high-resolution global gravity field model with the resolution 0.1 deg is compared with EGM2008.     

The iteration method for downward continuation of a potential field from a horizontal plane

This paper introduces an iteration method for the downward continuation of potential‐field data from a horizontal plane, and compares it with the conventional frequency‐domain method (Fourier transform) using 2D and 3D model tests. The paper evaluates the two methods in terms of the results, i.e. downward‐continuation distance and stability. The iteration method proves to be more stable and able to downward continue the potential for a greater distance than the Fourier transform method.     

Gravity potential at the Major Vertical Datum as primary geodetic constant

Summary On the basis of the fundamental relations of the Molodensky's Earth's figure theory (1945), admitting the inequality of the gravity potentials at the Major Vertical Datum W₀ and on the surface of the reference level ellipsoid U₀, and taken into account that potential W₀ enters into equations directly, it is recomended, W₀ should be adopted as a primary geodetic constant. Parameters of the best fitting ellipsoid are not needed for the solution of geodetic problems and for the investigation of the Earth's gravity field. The reason for requiring the reference and actual fields be close is only that the boundary-value problem can be solved in the linear approximation. Dedicated to the Memory of M.S. Molodensky Contribution to the I.A.G. Special Commission SC3 Fundamental Constants (SCFC).     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

重力－重力梯度边值问题的准解

依据超定边值问题的准解理论，推导出关于重力－重力梯度边值问题的准解，其简洁易算的解析形式给联合重力和卫星重力梯度数据精化地球重力场提供了一种新方法．     

Ellipsoidal Neumann geodetic boundary-value problem based on surface gravity disturbances: case study of Iran

An ellipsoidal Neumann type geodetic boundary-value problem (GBVP) for the computation of disturbing potential on the surface of the Earth based on the surface gravity disturbance as the boundary data is formulated. The solution methodology of the GBVP can be algorithmically summarized as follows: (i) using global navigation satellite systems (GNSS) coordinates of the gravity stations, the surface gravity disturbances are generated as the boundary data. (ii) Applying the deflection correction to the gravity disturbances to arrive at the derivative of the surface disturbing potential along the ellipsoidal normal. (iii) Removing the low frequencies part of the gravity field using harmonic expansion to degree and order 110. (iv) Using the short wavelength part of the corrected gravity disturbances derived in the previous section as the boundary data within the constructed GBVP to derive the short wavelength disturbing potential over the Earth surface. (v) The computed shortwave length signals of disturbing potentials are converted to disturbing potential values by restoring the removed effects.     

A Brovar-type solution of the fixed geodetic boundary-value problem

For more than 150 years gravity anomalies have been used for the determination of geoidal heights, height anomalies and the external gravity field. Due to the fact that precise ellipsoidal heights could not be observed directly, traditionally a free geodetic boundary-value problem (GBVP) had to be formulated which after linearisation is related to gravity anomalies. Since nowadays the three-dimensional positions of gravity points can be determined by global navigation satellite systems very precisely, the modern formulation of the GBVP can be based on gravity disturbances which are related to a fixed GBVP using the known topographical surface of the Earth as boundary surface. The paper discusses various approaches into the solution of the fixed GBVP which after linearization corresponds to an oblique-derivative boundary-value problem for the Laplace equation. Among the analytical solution approaches a Brovar-type solution is worked out in detail, showing many similarities with respect to the classical solution of the scalar free GBVP.     

地球外部扰动重力严密改化模型及分析检验

直接积分模型是计算地球外部扰动重力的主要数学工具,将全球积分模型改化为局域积分模型是实现地球外部重力场赋值的前提条件.相比表层积分模型和向上延拓积分模型,Stokes积分模型要求的输入信息种类最少,故在减轻数据保障压力方面具有比较明显的优势.在实施Stokes积分模型工程化应用过程中,我们发现扰动重力径向分量积分模型从...     

Gravity field determination using boundary element methods

The Boundary Element Method (BEM), a numerical technique for solving boundary integral equations, is introduced to determine the earth's gravity field. After a short survey on its main principles, we apply this method to the fixed gravimetric boundary value problem (BVP), i.e. the determination of the earth's gravitational potential from measurements of the intensity of the gravity field in points on the earth's surface. We show how to linearize this nonlinear BVP using an implicit function theorem and how to transform the linearized BVP into a boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation into a linear system of equations. We discuss the major problems of this approach for setting up and solving the linear system. The BVP is numerically solved for a bounded part of the earth's surface using a high resolution reference gravity model, measured gravity values of high density, and a 50 50 m² digital terrain model to describe the earth's surface. We obtain a gravity field resolution of 1 1 km² with an accuracy of the order 10^–3 to 10^–4 in about 1 CPU-hour on a Siemens/Fujitsu SIMD vector pipeline machine using highly sophisticated numerical integration techniques and fast equation solvers. We conclude that BEM is a powerful numerical tool for solving boundary value problems and may be an alternative to classical geodetic techniques.     

Some geophysical and geodetic contributions of satellite-determined gravity results

Satellite orbital data yield reliable values of low degree and order coefficients in the spherical harmonic expansion of the Earth's gravity field. The second degree coefficient yields the shape of the Earth — probably the most important single parameter in geodesy. It is crucial in the numerical evaluation of different forms of the theoretical gravity formula. The new information requires the standardization of gravity anomalies obtained from satellite gravity and terrestrial gravity data in the context of three most commonly used reference figures, e.g.,International Reference Ellipsoid, Reference Ellipsoid 1967, andEquilibrium Reference Ellipsoid. This standardization is important in the comparison and combination of satellite gravity and gravimetric data as well as the integration of surface gravity data, collected with different objectives, in a single reference system.Examination of the nature of satellite gravity anomalies aids in the geophysical and geodetic interpretation of these anomalies in terms of the tectonic features of the Earth and the structure of the Earth's crust and mantle. Satellite results also make it possible to compute the Potsdam correction and Earth's equatorial radius from the satellite-determined geopotential. They enable the decomposition of the total observed gravity anomaly into components of geophysical interest. They also make it possible to study the temporal variations of the geogravity field. In addition, satellite results make significant contributions in the prediction of gravity in unsurveyed areas, as well as in providing a check on marine gravity profiles.On leave from University of Hawaii, Honolulu.     

Spectral Combination of Spherical Gradiometric Boundary-Value Problems: A Theoretical Study

The Earth’s gravity potential can be determined from its second-order partial derivatives using the spherical gradiometric boundary-value problems which have three integral solutions. The problem of merging these solutions by spectral combination is the main subject of this paper. Integral estimators of biased- and unbiased-types are presented for recovering the disturbing gravity potential from gravity gradients. It is shown that only kernels of the biased-type integral estimators are suitable for simultaneous downward continuation and combination of gravity gradients. Numerical results show insignificant practical difference between the biased and unbiased estimators at sea level and the contribution of far-zone gravity gradients remains significant for integration. These contributions depend on the noise level of the gravity gradients at higher levels than sea. In the cases of combining the gravity gradients, contaminated with Gaussian noise, at sea and 250?km levels the errors of the estimated geoid heights are about 10 and 3 times smaller than those obtained by each integral.