2020珠峰测量与高程确定

为实现中国和尼泊尔共同宣布珠峰高程,我国于2019—2020年开展了珠峰高程测量工作,并于2020年5月27日完成峰顶测量。首次在珠峰北侧区域实施航空重力测量、开展峰顶地面重力测量,首次联合航空和地面重力等数据确定了基于国际高程参考系统(international height reference system,IHRS)的珠峰区域重力似大地水准面模型和峰顶大地水准面差距。此次珠峰测量,各种先进测量装备尤其是国产测量仪器全面担纲,通过多种技术手段相互验证和严密检核计算,确保了珠峰高程测量成果的精度和可靠性。最后,中尼双方合作开展数据处理,共同确定珠峰峰顶雪面正高(海拔高)为8848.86 m。     

2020与2005珠峰测量与高程确定异同

为了分析2020与2005珠峰高程测量与确定过程中的异同，该文从GNSS数据处理、高程控制网数据处理、峰顶交会数据处理、峰顶大地水准面差距计算4个方面对其进行了异同比较分析，阐述了2020珠峰高程测量的技术进步与创新。分析表明，2020珠峰高程测量从测量装备的国产化、测量手段和数据的丰富性，数据处理的高精度等多个方面，较2005年都有长足的进步。同时，中尼合作开展数据处理，共同确定了基于国际高程参考系统(international height reference system, IHRS)的珠峰正高。     

联合GRACE/GOCE重力场模型和GPS/水准数据确定我国85高程基准重力位

GRACE、GOCE卫星重力计划的实施,对确定高精度重力场模型具有重要贡献。联合GRACE、GOCE卫星数据建立的重力场模型和我国均匀分布的649个GPS/水准数据可以确定我国高程基准重力位,但我国高程基准对应的参考面为似大地水准面,是非等位面,将似大地水准面转化为大地水准面后确定的大地水准面重力位为62 636 854.395 3m~2s~(-2),为提高高阶项对确定大地水准面的贡献,利用高分辨率重力场模型EGM2008扩展GRACE/GOCE模型至2190阶,同时将重力场模型和GPS/水准数据统一到同一参考框架和潮汐系统,最后利用扩展后的模型确定的我国大地水准面重力位为62 636 852.751 8m~2s~(-2)。其中组合模型TIM_R4+EGM2008确定的我国85高程基准重力位值62 636 852.704 5m~2s~(-2)精度最高。重力场模型截断误差对确定我国大地水准面的影响约16cm,潮汐系统影响约4~6cm。     

利用大地边值问题方法统一全球高程基准

为解决世界各国高程基准差异的问题,提出联合卫星重力场模型、地面重力数据、GNSS大地高、局部高程基准的正高或正常高,按大地边值问题法确定局部高程基准重力位差的方法。首先推导了利用传统地面"有偏"重力异常确定高程基准重力位差的方法;接着利用改化Stokes核函数削弱"有偏"重力异常的影响,并联合卫星重力场模型和地面"有偏"重力数据,得到独立于任何局部高程基准的重力水准面,以此来确定局部高程基准重力位差;最后利用GNSS+水准数据和重力大地水准面确定了美国高程基准与全球高程基准W₀的重力位差为-4.82±0.05 m²s^-2。     

全球高程基准研究进展

建立统一的全球高程基准是国际大地测量科学界的核心目标之一,也是全球尺度地球科学研究、跨境工程应用等的必要基础设施。国际大地测量协会（international geodesy association,IAG）2015年发布了国际高程参考系统的定义,并于2019年提出了建立国际高程参考框架的目标。从全球高程参考系统的理论基础和定义出发,对国际高程参考系统与框架的理论、方法和实际问题开展论述与研究,主要包括全球大地水准面重力位$W_{\mathrm{0}}$的确定、基于高阶重力场模型的重力位确定、基于区域重力场建模的重力位确定,并重点论述和分析了IAG组织的科罗拉多大地水准面建模试验和中国2020珠峰高程测量实现国际高程参考系统2项典型案例研究。结果表明,在平坦地区和一般山区,重力大地水准面模型精度能达到1 cm（重力位0.1 m2/s2）,即使在珠穆朗玛峰这样的特大山区,也有望达到2~3 cm精度（重力位0.2~0.3 m2/s2）。综合典型案例研究结果、观测技术、数据资源和区域分布等因素,提出了建立国际高程参考框架的初步策略,包括IHRF参考站布设、重力位确定方法、数据要求、应遵循的标准/约定和预期精度指标等,展望了光学原子钟与相对论大地测量对于全球高程基准统一的潜在贡献。     

2005珠峰高程测定

2005年我国对珠穆朗玛峰高程进行了新的测定,为此在珠峰及其邻近地区开展了大规模的大地测量数据获取和数据处理工作。相对于1975年珠峰测高,2005年在珠峰以北地区的地面控制和珠峰高程测定中采用了GPS技术,采用了雷达探测技术测定珠峰峰顶冰雪覆盖层的深度,利用地球重力场模型、重力和数字地形数据、以及GPS水准等资料,精化珠峰地区的大地水准面,提高了测量珠峰高程和探测峰顶冰雪覆盖层深度的精度和可靠性。由此测得珠峰峰顶雪面正常高为8 846.67 M,珠峰峰顶雪面正高(海拔高)为8 847.93 M,珠峰峰顶岩面正高为8 844.43 M,珠峰峰顶相应点的冰雪层厚度为3.50 M。     

组合法精化胶州市大地水准面

本文利用全球重力位模型、胶州市地面重力观测数据、胶州市GPS水准数据和数字地面模型(DTM),采用组合法应用移去-恢复技术计算剩余大地水准面,并与地球位模型计算的高程异常进行拟合,得到该地区重力似大地水准面,再和布测、计算得到的GPS/水准所构成的几何大地水准面拟合,利用多项式拟合完成系统改正,获得最终的大地水准面结果及相关的精度信息。     

组合法精化胶州市大地水准面

本文利用全球重力位模型、胶州市地面重力观测数据、胶州市GPS水准数据和数字地面模型(DTM),采用组合法应用移去-恢复技术计算剩余大地水准面,并与地球位模型计算的高程异常进行拟合,得到该地区重力似大地水准面,再和布测、计算得到的GPS/水准所构成的几何大地水准面拟合,利用多项式拟合完成系统改正,获得最终的大地水准面结果及相关的精度信息。     

广东地区地球重力场模型的比较与选择

为了优化广东地区的似大地水准面,需要选择最适合该地区的地球重力场模型。本文对该地区6个具有代表性的地球重力场模型及简单组合模型,从高程异常、重力异常和垂线偏差3个方面的精度进行比较分析。结果表明,各模型所代表的似大地水准面与该地区的大地水准面之间均存在着一定的系统偏差,而EGM2008/EIGEN-6C简单组合模型精度则优于其他模型,此模型可作为该地区似大地水准面优化的参考重力场首选模型。     

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

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

利用等效点质量进行(似)大地水准面拟合

基于等效源原理,提出了一种半自由点质量模型,并给出了顾及相邻点空间关系构造虚拟点质量的简单快速迭代算法。实验结果表明,利用该点质量模型对离散GPS/水准观测数据进行拟合是可行的。     

Quasigeoid evaluation with improved levelled height data for Norway

GPS-levelling points are widely used to control gravimetric geoid or quasigeoid models. Direct comparison is often interpreted to reveal the accuracy of the gravimetric model, using GPS-levelling as a reference. However, both GPS and levelled heights contain errors, and in order to achieve a centimeter-accuracy geoid, these should be investigated. The Norwegian Height System NN1954 is known to contain large systematic errors due to postglacial land uplift in the area. In this study, the current height system and two revised versions, corrected for uplift, are applied to compute three sets of control quasigeoid heights in the southern part of Norway. These heights are then compared to various Nordic gravimetric quasigeoid models generated during the last two decades. In contradiction to some earlier studies, the accuracy of gravimetric quasigeoid models for this area are found to improve near-linearly with time. This is in accordance with expectations, since both data coverage and computation methods have progressed during this time. However, this study shows the importance of establishing accurate and error-free control data for geoid comparisons.     

On the geoid–quasigeoid separation in mountain areas

The separation between the reference surfaces for orthometric heights and normal heights—the geoid and the quasigeoid—is typically in the order of a few decimeters but can reach nearly 3 m in extreme cases. The knowledge of the geoid–quasigeoid separation with centimeter accuracy or better, is essential for the realization of national and international height reference frames, and for precision height determination in geodetic engineering. The largest contribution to the geoid–quasigeoid separation is due to the distribution of topographic masses. We develop a compact formulation for the rigorous treatment of topographic masses and apply it to determine the geoid–quasigeoid separation for two test areas in the Alps with very rough topography, using a very fine grid resolution of 100 m. The magnitude of the geoid–quasigeoid separation and its accuracy, its slopes, roughness, and correlation with height are analyzed. Results show that rigorous treatment of topographic masses leads to a rather small geoid–quasigeoid separation—only 30 cm at the highest summit—while results based on approximations are often larger by several decimeters. The accuracy of the topographic contribution to the geoid–quasigeoid separation is estimated to be 2–3 cm for areas with extreme topography. Analysis of roughness of the geoid–quasigeoid separation shows that a resolution of the modeling grid of 200 m or less is required to achieve these accuracies. Gravity and the vertical gravity gradient inside of topographic masses and the mean gravity along the plumbline are modeled which are important intermediate quantities for the determination of the geoid–quasigeoid separation. We conclude that a consistent determination of the geoid and quasigeoid height reference surfaces within an accuracy of few centimeters is feasible even for areas with extreme topography, and that the concepts of orthometric height and normal height can be consistently realized and used within this level of accuracy.     

基于GNSS水准和重力场误差特性的大地水准面精度评估方法

缺乏有效的大地水准面成果精度评估方法,是高程基准现代化及其成果应用面临的关键问题。本文基于GNSS水准高程异常与重力场频域误差特性,研究GNSS水准与重力地面高程异常融合的技术要求,进而提出一种大地水准面成果的误差表达与精度评估方法。经示例测试分析,得出主要结论如下:①实用地面高程异常(即融合后的似大地水准面)精度,应采用随距离非线性变化的高程异常差误差曲线表达;②似大地水准面的精度评估,推荐采用两项误差指标和两条误差曲线共4个要素完整表达,即重力地面高程异常差误差、实用地面高程异常内部误差、实用地面高程异常差误差曲线与GNSS水准高程异常差误差曲线;③当两个GNSS水准点间距离接近或小于所有GNSS水准点平均间距时,GNSS水准高程异常对实用地面高程异常的贡献起主要作用;④较大空间尺度的实用地面高程异常精度主要依靠重力地面高程异常控制。     

Using gravity and topography-implied anomalies to assess data requirements for precise geoid computation

Many regions around the world require improved gravimetric data bases to support very accurate geoid modeling for the modernization of height systems using GPS. We present a simple yet effective method to assess gravity data requirements, particularly the necessary resolution, for a desired precision in geoid computation. The approach is based on simulating high-resolution gravimetry using a topography-correlated model that is adjusted to be consistent with an existing network of gravity data. Analysis of these adjusted, simulated data through Stokes’s integral indicates where existing gravity data must be supplemented by new surveys in order to achieve an acceptable level of omission error in the geoid undulation. The simulated model can equally be used to analyze commission error, as well as model error and data inconsistencies to a limited extent. The proposed method is applied to South Korea and shows clearly where existing gravity data are too scarce for precise geoid computation.     

精化山区大地水准面的一种方法--GPS、测距三角高程与地形的组合

鉴于我国中西部是多山地区,为满足国家统一高程和山区建设的需要,本文提出了用ＧＰＳ和测距三角高程的方法确定高程异常,再用地形数据推求重力异常垂直梯度,由此可以经济和快速地实现精化大地水准面的目标。这时边远山区任意一点的精度可达３０ｃｍ。     

The Bruns formula in three dimensions

The Bruns formula is generalized to three dimensions with the derivation of equations expressing the height anomaly vector or the geoid undulation vector as a function of the disturbing gravity potential and its spatial derivatives. It is shown that the usual scalar Bruns formula provides not the separation along the normal to the reference ellipsoid but the component of the relevant spatial separation along the local direction of normal gravity. The above results which hold for any type of normal potential are specialized for the usual Somigliana-Pizzetti normal field so that the components of the geoid undulation vector are expressed as functions of the parameters of the reference ellipsoid, the disturbing potential and its spatial derivatives with respect to three types of curvilinear coordinates, ellipsoidal, geodetic and spherical. Finally the components of the geoid undulation vector are related to the deflections of the vertical in a spherical approximation.