局部大地水准面精化系统软件的设计与实现

给出了局部大地水准面精化系统的总体结构以及各分系统的详细设计与实现方法,开发了旨在对相关重力测量数据实现自动化处理的局部大地水准面精化系统工具软件,并应用于某区域实际大地水准面的精确计算。实践证明,该系统能有效提高相关数据的处理效率,有助于局部大地水准面精化过程的规范化。     

我国海域大地水准面与大陆大地水准面的拼接研究

由重力和GPS水准数据确定的陆地大地水准面和主要由卫星测高数据确定的海洋大地水准面 ,两者之间一般都存在以系统误差为主的拼接差 ,分析了产生这一现象的主要原因。顾及这一现象 ,并结合我国在陆海大地水准面拼接区重力资料稀疏的实际 ,提出了扩展拼接技术。以国家GPS水准网确定的 (似 )大地水准面作控制 ,对陆海重力大地水准面作拟合校正 ,得到我国校正的陆海统一的重力大地水准面     

海洋重力似大地水准面与区域测高似大地水准面的拟合问题

研究了将陆地重力似大地水准面与GPS水；住似大地水准面拟合的处理方法推广到海洋的问题．首先从理论上证明了当存在海面地形．则海洋大地水准面与似大地水准面不重合．导出了在海洋上大地水；住面差距与高程异常之间差值的公式．由此给出了求定平均海面相对于区域高程基准的正常高以及测高似大地水准面的计算公式。由于测高平均海面与GPS大地高有相近的精度．提出了将海洋重力似大地水准面与区域测高似大地水准面拟合的处理方法．并利用当前最新的海面地形模型和测高平均海面模型做了数值估计。     

论(N-ζ)公式的内涵及推求N的精度

介绍了大地水准面(N)与似大地水准面(ξ)之差公式的推导方法，揭示了其物理内涵，论述了其在精化大地水准面中的重要作用。当似大地水准面精度达到cm级时，大地水准面也将达到同样的精度。     

数据源对大地水准面精度的影响分析

组合法是目前我国省市进行大地水准面精化普遍采用的方法。通过对大地水准面精化时所用到的各种数据源进行分析,探讨各数据源对大地水准面精度的影响,并提出提高大地水准面精度的办法。     

无锡市厘米级似大地水准面的研究

利用大地水准面逼近的严密理论M olodenskii级数解和重力归算的地形均衡理论,结合全球定位系统(G PS)、水准和重力资料在无锡确定了精度达厘米级的似大地水准面,并讨论了无锡G PS控制网的布设及数据处理。通过检核,所确定的2.5'×2.5'似大地水准面精度为2.2cm。     

论相对论重力位及相对论大地水准面

讨论了相对论意义下的重力位及大地水准面，指出了等时率大地水准面的缺陷，建议今后采用等频面及等频大地水准面的概念，给出了等频大地水准面与经典大地水准面的差异，同时给出了等频大地水准面的近似表达式。     

大庆市区域似大地水准面精化

精化大地水准面是现代重力场和大地测量学工作的重要任务之一,本文介绍大庆市区域似大地水准面项目所采用的技术路线和方法。论述了GPS和水准路线的布测精度,GPS、水准数据处理和似大地水准面确定所采用的方案,及似大地水准面所获取的精度。     

青藏高原地区大地水准面与似大地水准面的转换

研究了联合重力场模型和地形系数模型直接计算大地水准面与似大地水准面差值的方法;分析了影响计算结果精度的误差源。提出利用实际地形数据SRTM3以及密度数据CRUST2.0改善计算结果的建议,并分别进行了比较。青藏高原地区的实际计算结果表明,该地区大地水准面与似大地水准面的差值均值约为0.8 m。     

《区域大地水准面精化外业作业指导书》的编制

对《区域大地水准面精化外业作业指导书》编制情况做了介绍,阐述了外业作业的流程、方法、内容、技术要求等,有利于规范区域大地水准面精化外业生产,有助于提高产品质量。     

卫星测高与卫星重力对洋流的研究

卫星测高与卫星重力的发展为进一步研究洋流提供了前所未有的机遇。本文从EGM 96 ,GGM0 1与未来的GOCE任务获取的高精度高分辨率海洋大地水准面的角度对洋流的研究方法与可行性进行了分析     

似大地水准面的误差分析与抑制技术

大地水准面误差分析与精度评定是局部重力场逼近技术的重要组成部分,是大地水准面精化工程外业方案优化、算法设计和工程质量评价的基本依据。本文分别从地面重力数据误差和局部重力场算法两个方面,分析cm级大地水准面误差的影响特性,提出重力数据误差与大地水准面精度之间普遍适用的规律,推荐一种GPS水准和地面重力数据联合平差的精度评定方法,结合实例和模拟计算分析,介绍大地水准面误差分析与误差抑制方法。     

Geoid determination using one-step integration

A residual (high-frequency) gravimetric geoid is usually computed from geographically limited ground, sea and/or airborne gravimetric data. The mathematical model for its determination from ground gravity is based on the transformation of observed discrete values of gravity into gravity potential related to either the international ellipsoid or the geoid. The two reference surfaces are used depending on height information that accompanies ground gravity data: traditionally orthometric heights determined by geodetic levelling were used while GPS positioning nowadays allows for estimation of geodetic (ellipsoidal) heights. This transformation is usually performed in two steps: (1) observed values of gravity are downward continued to the ellipsoid or the geoid, and (2) gravity at the ellipsoid or the geoid is transformed into the corresponding potential. Each of these two steps represents the solution of one geodetic boundary-value problem of potential theory, namely the first and second or third problem. Thus two different geodetic boundary-value problems must be formulated and solved, which requires numerical evaluation of two surface integrals. In this contribution, a mathematical model in the form of a single Fredholm integral equation of the first kind is presented and numerically investigated. This model combines the solution of the first and second/third boundary-value problems and transforms ground gravity disturbances or anomalies into the harmonically downward continued disturbing potential at the ellipsoid or the geoid directly. Numerical tests show that the new approach offers an efficient and stable solution for the determination of the residual geoid from ground gravity data.     

我国现代高程测定关键技术若干问题的研究及进展

阐述了大地水准面在大地测量学及相关地球科学中的科学意义、作用及其应用,提出了利用GPS技术结合高精度似大地水准面代替我国高程测定的思路,分析了目前地球重力场长波、中波和短波在计算cm级大地水准面中具备的条件,论述了确定1 cm精度城市和5 cm省级似大地水准面的若干关键技术问题及其必要性和可行性,提出了我国cm级大地水准面确定中需解决的主要理论问题和若干关键技术,评述了我国局部大地水准面的主要进展和取得的成果。     

Geoid and high resolution sea surface topography modelling in the mediterranean from gravimetry,altimetry and GOCE data: evaluation by simulation

The determination of local geoid models has traditionally been carried out on land and at sea using gravity anomaly and satellite altimetry data, while it will be aided by the data expected from satellite missions such as those from the Gravity field and steady-state ocean circulation explorer (GOCE). To assess the performance of heterogeneous data combination to local geoid determination, simulated data for the central Mediterranean Sea are analyzed. These data include marine and land gravity anomalies, altimetric sea surface heights, and GOCE observations processed with the space-wise approach. A spectral analysis of the aforementioned data shows their complementary character. GOCE data cover long wavelengths and account for the lack of such information from gravity anomalies. This is exploited for the estimation of local covariance function models, where it is seen that models computed with GOCE data and gravity anomaly empirical covariance functions perform better than models computed without GOCE data. The geoid is estimated by different data combinations and the results show that GOCE data improve the solutions for areas covered poorly with other data types, while also accounting for any long wavelength errors of the adopted reference model that exist even when the ground gravity data are dense. At sea, the altimetric data provide the dominant geoid information. However, the geoid accuracy is sensitive to orbit calibration errors and unmodeled sea surface topography (SST) effects. If such effects are present, the combination of GOCE and gravity anomaly data can improve the geoid accuracy. The present work also presents results from simulations for the recovery of the stationary SST, which show that the combination of geoid heights obtained from a spherical harmonic geopotential model derived from GOCE with satellite altimetry data can provide SST models with some centimeters of error. However, combining data from GOCE with gravity anomalies in a collocation approach can result in the estimation of a higher resolution geoid, more suitable for high resolution mean dynamic SST modeling. Such simulations can be performed toward the development and evaluation of SST recovery methods.     

EGM96,WDM94和GPM98CR高阶地球重力场模型表示深圳局部重力场的比较与评价

为计算深圳精密重力大地水准面，利用62个高精度GPS水准点和4871个实测重力点数据对EGM96，WDM94和GPM98CR全球重力场模型表示深圳局部重力场进行了比较与评价。结果表明，由上述3个重力场模型计算的大地水准面高和重力异常与实测值之间存在明显的系统偏差，当采用GPS水准数据尽可能消除系统偏差以后，大地水准面高的精度得到显著提高，若应用移去-恢复技术确定深圳高精度大地水准面，则WDM94应该是首选的参考重力场模型。     

On the combination of gravity anomalies and gravity disturbances for geoid determination in Western Australia

The geoid gradient over the Darling Fault in Western Australia is extremely high, rising by as much as 38 cm over only 2 km. This poses problems for gravimetric-only geoid models of the area, whose frequency content is limited by the spatial distribution of the gravity data. The gravimetric-only version of AUSGeoid98, for instance, is only able to resolve 46% of the gradient across the fault. Hence, the ability of GPS surveys to obtain accurate orthometric heights is reduced. It is described how further gravity data were collected over the Darling Fault, augmenting the existing gravity observations at key locations so as to obtain a more representative geoid gradient. As many of the gravity observations were collected at stations with a well-known GRS80 ellipsoidal height, the opportunity arose to compute a geoid model via both the Stokes and the Hotine approaches. A scheme was devised to convert free-air anomaly data to gravity disturbances using existing geoid models, followed by a Hotine integration to geoid heights. Interestingly, these results depended very weakly upon the choice of input geoid model. The extra gravity data did indeed improve the fit of the computed geoid to local GPS/Australian Height Datum (AHD) observations by 58% over the gravimetric-only AUSGeoid98. While the conventional Stokesian approach to geoid determination proved to be slightly better than the Hotine method, the latter still improved upon the gravimetric-only AUSGeoid98 solution, supporting the viability of conducting gravity surveys with GPS control for the purposes of geoid determination. AcknowledgementsThe author would like to thank Will Featherstone, Ron Gower, Ron Hackney, Linda Morgan, Geoscience Australia, Scripps Oceanographic Institute and the three anonymous reviewers of this paper. This research was funded by the Australian Research Council.     

A local geoid determination based on ellipsoid approximation in Western China

A new methodology for precise geoid determination with finest local details based on ellipsoidal approximation is presented.This methodology is formulated through the “fixed-free two-boundary value problem“ based on the observable of the type modulus of gravity intensity,gravity acceleration and gravity potential at the GPS positioned stations,with support of the known geoid‘s potential value,W0.     

A new isostatic model of the lithosphere and gravity field

Based on the analysis of various factors controlling isostatic gravity anomalies and geoid undulations, it is concluded that it is essential to model the lithospheric density structure as accurately as possible. Otherwise, if computed in the classical way (i.e. based on the surface topography and the simple Airy compensation scheme), isostatic anomalies mostly reflect differences of the real lithosphere structure from the simplified compensation model, and not necessarily the deviations from isostatic equilibrium. Starting with global gravity, topography and crustal density models, isostatic gravity anomalies and geoid undulations have been determined. The initial crust and upper-mantle density structure has been corrected in a least squares adjustment using gravity. To model the long-wavelength (>2000 km) features in the gravity field, the isostatic condition (i.e. equal mass for all columns above the compensation level) is applied in the adjustment to uncover the signals from the deep-Earth interior, including dynamic deformations of the Earths surface. The isostatic gravity anomalies and geoid undulations, rather than the observed fields, then represent the signals from mantle convection and deep density inhomogeneities including remnants of subducted slabs. The long-wavelength non-isostatic (i.e. the dynamic) topography was estimated to range from –0.4 to 0.5 km. For shorter wavelengths (<2000 km), the isostatic condition is not applied in the adjustment in order to obtain the non-isostatic topography due to regional deviations from classical Airy isostasy. The maximum deviations from Airy isostasy (–1.5 to 1 km) occur at currently active plate boundaries. As another result, a new global model of the lithosphere density distribution is generated. The most pronounced negative density anomalies in the upper mantle are found near large plume provinces, such as Iceland and East Africa, and in the vicinity of the mid-ocean ridge axes. Positive density anomalies in the upper mantle under the continents are not correlated with the cold and thick lithosphere of cratons, indicating a compensation mechanism due to thermal and compositional density.     

我国精化大地水准面工作中若干问题的讨论

对我国正进行的精化大地水准面工作中若干问题作了讨论,包括:1)推估内插(或拟合)模型;2)利用莫霍面起伏资料代替均衡假说;3)GPS水准布网的优化设计;4)利用GPS水准及重力资料精化大地水准面的精度估计。