高精度天文水准的布设

天文水准是利用垂线偏差确定高程异常的一种经典方法。过去由于垂线偏差测量的作业效率低而不可能大量布测垂线偏差点，我国的天文水准测线上垂线偏差点间距为２０－５４ｋｍ，因而精度很低（μ＝±０．０７-０．１１ｍ）。当前发展中的垂线偏差快速确定技术为高精度天文水准创造了条件。本文讨论了垂线偏差的代表误差和水准面不平行的改正，提出了高精度（μ＝±０．０１ｍ）天文水准布设方案，用物理大地测量实验区实测数据对结论进行了验证。     

珠峰地区似大地水准面的精确求定

１９９４年和１９９９年利用ＧＰＳ精密星历和ＩＥＲＳ参考站,对１９９２年珠峰地区的ＧＰＳ网观测资料再次进行了计算,精确求定了该地区相对ＷＧＳ８４椭球的高程异常和垂线偏差。与１９９３年计算成果相比,垂线偏差的平均差异在精度范围内,而高程异常的平均差值约为－２５ｍ。     

中国海域大地水准面和重力异常的确定

从莫洛金斯基(Molodensky)等1960年给出的由垂线偏差计算大地水准面空域积分公式出发,导出了其相应谱域1维严密卷积和2维球面及平面卷积公式。由Topex/Poseidon,ERS 1/2及Geosat/GM,ERM测高资料求解的垂线偏差计算了我国海域及其邻区大地水准面,其中计算格网为2.5′×2.5′。为了检核,将测高垂线偏差由逆维宁 迈尼兹(Vening Meinesz)公式反演重力异常,与海上船测重力值进行了外部检核;同时还利用司托克斯(Stokes)公式,由上述反演的重力异常计算大地水准面高,与莫洛金斯基公式直接解得的相应结果进行比较作为内部检核。前者的中误差为±9mGal(1Gal=1cm/s2),后者为±0.025m。本文在积分计算中充分应用了2维平面坐标形式和1维卷积严格公式,并做了比较和自校核。     

广义带限航空矢量重力确定大地水准面的两步积分法EI北大核心CSCD

依据物理大地测量学的广义带限水平边值问题理论,本文提出了基于带限航空矢量重力水平分量确定大地水准面的两步积分法。首先,依据广义带限水平边值问题解将带限航空矢量重力水平分量转化为飞行高度面上的带限扰动位;然后,利用带限Dirichlet边值问题解将飞行高度面上的带限扰动位向下延拓到海平面(海拔高程起算面,也称平均海面)上的带限扰动位;最后,依据Bruns公式将海平面上的扰动位转化为大地水准面高度。利用EGM2008重力位模型开展数值仿真计算试验,结果表明,广义带限水平边值问题解具有较好的低通滤波特性,能有效抑制观测高频噪声的影响;当带限航空矢量重力水平分量观测误差取3×10^(-5)m/s 2、测量飞行高度取6 km时,基于带限Dirichlet边值问题解的带限大地水准面计算精度优于3 cm,初步验证了两步积分法的计算稳定性和有效性。     

最新中国陆地数字高程基准模型:重力似大地水准面CNGG2011

本文回顾了近20年国内外国家局部大地水准面模型研究的概况和发展背景,采用Stokes-Helmert方法,计算了一个新的2′×2′中国重力和1985国家高程基准似大地水准面数值模型（CNGG2011）,采用了1百万余陆地重力数据和SRTM 7″.5×7″.5地形高数据,以及649个B级GPS水准点数据。CNGG2011平均精度为±0.13m,东部地区±0.07m,西部地区±0.14m。各省区局部似大地水准面平均精度为±0.06m,东部为±0.05m,西部为±0.11m。西藏精度为±0.22m。本文还讨论了重力大地水准面与GPS水准的关系,提出了今后进一步精化我国高程基准大地水准面模型的构想。     

高原山区局部大地水准面精化方法研究

局部大地水准面精化的实质是精确计算出大地水准面的起伏变化情况。一般情况下,需要密度足够的重力数据,依重力异常密集计算大地水准面差距或高程异常。但是在大陆西部高原山区重力点密度是不够的,无法达到大地水准面精化的目的。本文从理论上证实了用地形和岩石密度数据进行局部大地水准面精化的可行性。     

ENVISAT测高卫星沿轨大地水准面梯度的海洋垂线偏差法研究

研究了利用沿轨大地水准面梯度数据计算海洋垂线偏差的最小二乘法,首先对ENVISAT测高数据进行各项地球物理改正得到近似测高大地水准面,然后计算沿轨大地水准面的梯度,接着用最小二乘法计算格网垂线偏差东西分量和南北分量的平均值。最后,用该方法计算了南中国海区域及其邻近海域(4°N～25°N,104°E～120°E)的5′×5′垂线偏差南北分量和东西分量,其精度优于7″,并与EGM96模型计算的垂线偏差值进行了比较,证明了该方法的有效性。     

我国陆海统一似大地水准面构建的三维重力矢量法

基于重力场水平分量-垂线偏差对地形信息敏感的特点,根据边值理论由重力与地形数据确定格网垂线偏差模型,在此基础上,首先利用三维重力矢量-格网垂线偏差与格网重力异常,联合格网高程数据求得格网点间高程异常差,然后通过GPS/水准点的控制,构成紧密的几何条件,进行严密平差,从而获得高分辨率、高精度似大地水准面的数值模型。按照本文方法,利用我国6600多个GPS/水准点、1'×1'的格网垂线偏差、格网重力异常、格网高程数据,整体平差计算了我国陆海统一的似大地水准面模型,经GPS/水准点检核,全国似大地水准面的绝对精度达到了4 cm,相对精度优于7 cm。     

Solution of the geodetic boundary value problem

Summary The possibility of improving the convergence of Molodensky’s series is considered. Then new formulas are derived for the solution of the geodetic boundary value problem. One of them presents an expression for the boundary condition which involves a linear combination of Stokes’ constants and surface gravity anomalies. This differs from the usually used relation by the appearance of additional terms dependent on second order terns with respect to the elevations of the earth’s surface. The formulas are derived for Stokes’ constants and the anomalous potential in terms of surface anomalies. As compared to the Taylor’s series of Molodensky, they are presented in the form of surface harmonic series. Due regard is made to the earth’s oblateness, in addition to the terrain topography.     

Variational formulation of the geodetic boundary value problem

A variational principle for the Stokesian boundary value problem is derived using the Euler-Lagrange theory. The resulting variational principle is then transformed into an equation determining the semi-major axis of the best fitting ellipsoid which fulfills the conditionU ₀ =W ₀ . The computations using three different geopotential models yields the semi-major axis of the earth ellipsoid asa=6378145.4 metres for the flatteningf=1/298.2564. The corresponding equatorial gravity and the geopotential number are computed as γ_a=978029.59 mgals andU ₀=W ₀=6.26367371 10⁶ kgalmeters respectively.     

The geodetic boundary value problem and the coordinate choice problem

Summary The geodetic boundary value problem (g.b.v.p.) is a free boundary value problem for the Laplace operator: however, under suitable change of coordinates, it can be transformed into a fixed boundary one. Thus a general coordinate choice problem arises: two particular cases are more closely analyzed, namely the gravity space approach and the intrinsic coordinates (Marussi) approach.     

A bias-free geodetic boundary value problem approach to height datum unification

A geodetic boundary value problem (GBVP) approach has been formulated which can be used for solving the problem of height datum unification. The developed technique is applied to a test area in Southwest Finland with approximate size of 1.5° × 3° and the bias of the corresponding local height datum (local geoid) with respect to the geoid is computed. For this purpose the bias-free potential difference and gravity difference observations of the test area are used and the offset (bias) of the height datum, i.e., Finnish Height Datum 2000 (N2000) fixed to Normaal Amsterdams Peil (NAP) as origin point, with respect to the geoid is computed. The results of this computation show that potential of the origin point of N2000, i.e., NAP, is (62636857.68 ± 0.5) (m²/s²) and as such is (0.191 ± 0.003) (m) under the geoid defined by W ₀ = 62636855.8 (m²/s²). As the validity test of our methodology, the test area is divided into two parts and the corresponding potential difference and gravity difference observations are introduced into our GBVP separately and the bias of height datums of the two parts are computed with respect to the geoid. Obtaining approximately the same bias values for the height datums of the two parts being part of one height datum with one origin point proves the validity of our approach. Besides, the latter test shows the capability of our methodology for patch-wise application.     

On the non-linear geodetic boundary value problem for a fixed boundary surface

Summary The fixed gravimetric boundary value problem of Physical Geodesy is a nonlinear, oblique derivative problem. Expanding the non-linear boundary condition into a Taylor series—based upon some reference potential field approximating the geopotential—it is shown that the numerical convergence of this series is very rapid; only the quadratic term may have some practical impact on the solution. The secondorder solution term can be described by a spherical integral formula involving the deflections of the vertical with respect to the reference field. The influence of nonlinear terms on the figure of the level surfaces (e.g. the geoid) is roughly estimated to have an order of magnitude of some few centimetres, based upon a Somigliana-Pizzetti reference field; if on the other hand some high-degree geopotential model is used as reference then the effects by non-linearity are negligible from a practical point of view.     

An overdetermined geodetic boundary value problem approach to telluroid and quasi-geoid computations

In this paper an overdetermined Geodetic Boundary Value Problem (GBVP) approach for telluroid and quasi-geoid computations is presented. The presented GBVP approach can solve the problem of potential value computation on the surface of the Earth, which when applied to a mapping scheme, e.g., here minimum distance mapping, provides a point-wise approach to telluroid computation. Besides, we have succeeded in reducing the number of equations and unknowns of the minimum distance telluroid mapping by one. The sufficient condition of minimum distance telluroid mapping is also recapitulated. Since the introduced GBVP approach has the advantage of implementing various gravity observables simultaneously as input boundary data, it can be regarded as a data fusion technique that exploits all available gravity data. The developed GBVP is used for the computation of the quasi-geoid within a test area in Southwest Finland.     

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.     

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.     

Singularity free formulations of the geodetic boundary value problem in gravity-space

The idea of transforming the geodetic boundary value problem into a boundary value problem with a fixed boundary dates back to the 1970s of the last century. This transformation was found by F. Sanso and was named as gravity-space transformation. Unfortunately, the advantage of having a fixed boundary for the transformed problem was counterbalanced by the theoretical as well as practical disadvantage of a singularity at the origin. In the present paper two more versions of a gravity-space transformation are investigated, where none of them has a singularity. In both cases the transformed differential equations are nonlinear. Therefore, a special emphasis is laid on the linearized problems and their relationships to the simple Hotine-problem and to the symmetries between both formulations. Finally, in numerical simulation study the accuracy of the solutions of both linearized problems is studied and factors limiting this accuracy are identified.     

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.