海洋重力场特征统计模型计算与分析

海洋重力场特征参数在地球重力场逼近计算和海上测量优化设计中具有重要的应用价值。基于卫星测高重力在海域具有覆盖范围广且分布均匀的独特优势，提出了利用最新卫星测高重力数据集开展海洋重力场特征统计模型计算和分析的研究方案，给出了代表误差和协方差函数模型参数的计算公式，定义并研究了海洋广义布格重力异常的变化特征，提出了等精度和非等精度拟合经验协方差函数的计算模型。利用中国近海及西太平洋海区超过50万个5'×5'方块的1'×1'网格卫星测高重力异常数据，首次计算得到一组有代表性的中国周边海域重力场特征统计模型参数，较好地揭示了海洋重力场有别于陆地重力场的变化特征，利用海面船测重力数据对计算结果进行了可靠性检核，提出了相应的模型参数修正方案和使用建议。     

CHAMP重力卫星结果在中国海及邻近海域的初步分析

利用世界上第一个采用高低卫_卫跟踪技术的CHAMP重力卫星计划导出的全新的高精度全球长波重力场模型EIGEN_1S结果 ,根据二维高斯滤波原理 ,基于相同空域尺度对卫星重力结果和由卫星测高解算的几种版本的海洋重力异常数据进行了长波部分的分析研究。结果表明 ,在中国海及其邻近海域卫星重力数据与卫星测高解算的海洋重力异常数据之间存在明显的偏差。     

联合卫星重力和卫星测高数据确定稳态海洋动力地形

新一代卫星重力探测任务GRACE大大提高了地球重力场模型中长波分量的精度,使得联合卫星测高平均海面分离更精细稳态海洋动力地形成为可能。本文利用T/P(1994年~2003年)和JASON-1(2003年~2005年)卫星测高数据确定了全球30′×30′平均海面高;基于重力场模型WHU-GM-05,计算得到对应于海面高的GRACE海洋大地水准面格网值;利用“移去-恢复法”和高斯滤波求得全球稳态海面地形。与EGM96、R io05、ECCO和GGM02模型进行比较,检验结果表明GRACE任务有效的改善了海洋大地水准面的精度,使得稳态海洋动力地形能够呈现更多细部。     

垂线偏差-大地水准面与逆范宁-梅尼兹公式的应用——联合多种卫星测高数据求解全球海平面与海洋重力异常

卫星测高开辟了地球科学研究的新纪元。平均海面高(MSSH)和海洋重力异常是卫星测高在大地测量和地球物理应用领域的重要成果。平均海面高对于建立全球潮汐模型、海平面变化研究以及海洋测深等方面有重要意义。关于海洋重力异常的应用,Sandweel与Smith(1997)已进行了说明。根据卫星测高数据计算垂线偏差(DOV)是恢复重力异常的一个新途径。对利用改进的卫星测高数据进一步提高全球平均海面高和重力异     

卫星重力与卫星测高的研究进展

该文介绍了GRACE和GOCE重力卫星在地球重力场模型以及局部变化等方面的研究成果,介绍了卫星测高的检核与校验、全球和局部海平面变化监测、内陆湖和水库水位变化监测等方面的研究进展,对综合利用卫星重力、卫星测高以及合成孔径雷达差分干涉测量(DInSAR)等技术方法监测冰川和冰盖变化的成果进行了介绍,重点介绍了南极冰盖变化监测的研究成果。对综合卫星重力、卫星测高、DInSAR、GNSS/水准、GNSS等多种技术在冰川冰盖、局部形变监测等方面的应用进行了展望。     

融合卫星测高和船测重力数据的非均匀密度点质量方法

融合多源海洋重力数据是提升极地、近岸和海底地形复杂区域卫星测高反演海洋重力场模型精度的重要途径.利用国际公开的船测重力数据与测高重力场模型,以琉球群岛附近海域为例,本文分析了点质量密度和船测重力数据分布对点质量模型融合精度的影响,结合船测重力数据分布,构建了基于迭代移去恢复的非均匀密度点质量模型.在相同船测重力约束下,非均匀密度点质量模型融合精度较均匀密度的点质量模型提升可达12％.     

联合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。     

卫星测高误差对海洋重力场反演的影响

本文研究了卫星测高误差对海洋重力场探测的影响,导出了卫星测高随机误差在海洋重力异常推估过程中的误差传播公式,实验分析了卫星测高随机误差在多种测量分辨率下对重力异常推估值影响的大小,可为我国未来的测高卫星系统论证设计提供技术参考。     

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

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

全国及部分省市地区高精度、高分辨率大地水准面的研究及其实施

结合我国重力和地形资料及国内外较优的重力场模型，研制了适合我国重力场特征的360阶重力场模型WDM94，建立了中国新一代包括全部陆海国土的dm级（似）大地水准面CQG2000，建立了中国以GPS／水准为基础的高程异常控制网，利用海洋卫星测高数据进行我国海洋大地水准面的计算、我国陆地重力（似）大地水准面的研制厦我国陆海（似）大地水准面的拼接；研制了江苏省、海南省、深圳市、大连市、南京市及南水北调西线工程具有cm级精度的省市地区（似）大地水准面模型；结合GPS技术和高精度（似）大地水准面模型，研制了GPS测图软硬件一体化系统。     

A strategy for determining the regional geoid by combining limited ground data with satellite-based global geopotential and topographical models: a case study of Iran

The computation of regional gravimetric geoid models with reasonable accuracy, in developing countries, with sparse data is a difficult task that needs great care. Here we investigate the procedure for gathering, evaluating and combining different data for the determination of a gravimetric geoid model for Iran, where limited ground gravity data are available. Heterogeneous data, including gravity anomalies, the high-resolution Shuttle Radar Topography Mission global digital terrain model and different global geopotential models including recently published Gravity Recovery and Climate Experiment models, are combined through least-squares modification of the Stokes formula. The new gravimetric geoid model, IRG04, agrees considerably better with GPS/levelling than any of the other recent local geoid model in the area. Its RMS fit with GPS/levelling is 0.27 m and 3.8 ppm in the absolute and relative view, respectively. The relative accuracy of IRG04 is four times better than the most recently published global and regional geoid models available in this area. This progress shows the practical potential of the method of least-squares modification of Stokes’s formula in combination with heterogeneous data for regional geoid determination     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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

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

On the combined adjustment of a gravimetrically determined geoid and GPS levelling stations

Summary A new basic geodetic network using the GPS technique is now being set up in France. There will be altogether 1000 benchmarks connected to the French levelling network. Obviously, the GPS levelling points are not dense enough to produce a national levelling reference surface. A gravimetrically determined geoid has therefore been proposed to be used for the interpolation between the GPS levelling points. However, because of long-wavelength errors, we consider that a gravimetric geoid does not have sufficient accuracy. A regression by fitting the gravimetrically determined geoid to the GPS levelling points is generally proposed. Unfortunately, this country-wide geoid fitting work cannot eliminate local deformations in the geoid, which happen in areas where there are errors or shortages of gravity or DTM data. This paper proposes and discusses a combined adjustment method. The principle is to divide up the geoid into small pieces and then to adjust them to the GPS levelling points locally with constraint conditions for the common points of the adjacent pieces. In order to benefit from the advantages of the high resolution and high relative accuracy of the gravimetric geoid, as well as the high absolute accuracy of the GPS levelling points, we establish respectively a relative observation equation for the difference of the gravimetric geoid undulation and an absolute observation equation for the GPS levelling points. Finally, we adjust the observation equations as a whole. Several global and local systematic errors are also taken into account and some special cases, such as adjustment in groups and blunder detection, are also discussed.     

Earth gravity field recovered from CHAMP science orbit and accelerometer data

The earth gravity field model CDS01S of degree and order 36 has been recovered from the post processed Science Orbits and on-board accelerometer data of GFZ's CHAMP satellite. The model resolves the geoid with an accuracy of better than 4 cm at a resolution of 700 km half-wavelength. By using the degree difference variances of geopotential coefficients to compare the model CDS01S with EIGEN3P, EIGEN1S and EGM96, the result indicates that the coefficients of CDS01S are most close to those of EIGEN3P. The result of the comparison between the accuracies of geopotential coefficients in the above models, indicates that the accuracy of coefficients in CDS01S is higher than that in EGM96. The geoid undulations of CDS01S and GGM01C up to 30 degrees are calculated and the standard deviation is 4. 7 cm between them.     

Earth gravity field recovered from CHAMP science orbit and accelerometer data

IntroductionSince the launch of man-made satellite early in1957 ,the research for satellite gravity has beentaken a wide attentioninfield of geodesy .Early ,the ground-based satellite tracking has providedan observational data set which has been used tode…     

On the adjustment of combined GPS/levelling/geoid networks

A detailed treatment of adjustment problems in combined global positioning system (GPS)/levelling/geoid networks is given. The two main types of `unknowns' in this kind of multi-data 1D networks are usually the gravimetric geoid accuracy and a 2D spatial field that describes all the datum/systematic distortions among the available height data sets. An accurate knowledge of the latter becomes especially important when we consider employing GPS techniques for levelling purposes with respect to a local vertical datum. Two modelling alternatives for the correction field are presented, namely a pure deterministic parametric model, and a hybrid deterministic and stochastic model. The concept of variance component estimation is also proposed as an important statistical tool for assessing the actual gravimetric geoid noise level and/or testing a priori determined geoid error models. Finally, conclusions are drawn and recommendations for further study are suggested. Received: 9 September 1998 / Accepted: 8 June 1999     

用重力异常逐级余差计算重力大地水准面

本文将计算重力大地水准面的频域方法推广至空域，提出了一种新的用重力数据和重力模型位系数联合确定大地水准面的方法。利用重力异常的逐级余差实施积分，使得通常的Ｓｔｏｋｅｓ积分方法具有明确的频域分析含义，可按精度要求确定出使用重力异常余差的块形大小及积分半径ψｏ。     

星载原子干涉技术用于地球重力场测量及其精度评估

重力梯度卫星GOCE通过搭载静电式重力梯度仪,将全球静态重力场恢复至200阶以上。目前GOCE卫星已结束寿命,亟须发展下一代更高分辨率的卫星重力梯度测量来完善200~360阶的全球静态重力场模型。原子干涉型的重力梯度测量在空间微重力环境下可获得较长的干涉时间,因此具有很高的星载测量精度,是下一代卫星重力梯度测量的候选技术之一。本文针对未来更高分辨率全球重力场测量的科学需求,提出了一种适用于空间微重力环境下的原子干涉重力梯度测量方案,其梯度测量噪声可低至0.85mE/Hz1/2。文中对不同类型的卫星重力梯度测量方案进行了重力场反演精度的对比评估,仿真结果表明,相比于现有静电式卫星重力梯度测量,原子干涉型的卫星重力梯度测量有望将重力场的恢复阶数提升至252~290阶,对应的累积大地水准面误差7~8cm,累积重力异常误差3×10-5 m/s2。     

Numerical study on the mixed model in the GOCE polar gap problem

Gravity gradients acquired by the Gravity field and steady-state Ocean Circulation Explorer(GOCE) do not cover the entire earth because of its sun-synchronous orbit leaving data gaps with a radius of about 6.5° in the polar regions.Previous studies showed that the loss of data in the polar regions deteriorates the accuracy of the low order(or near zonal) coefficients of the earth gravity model,which is the so-called polar gap problem in geodesy.In order to find a stable solution for the earth gravity model from the GOCE gravity gradients,three models,i.e.the Gauss-Markov model,light constraint model and the mixed model,are compared and evaluated numerically with the gravity gradient simulated with the EGM2008.The comparison shows that the Best Linear Uniformly Unbiased Estimation(BLUUE) estimator of the mixed model can solve the polar gap problem as effectively as the light constraint model;furthermore,the mixed model is more rigorous in dealing with the supplementary information and leads to a better accuracy in determining the global geoid.