测高重力异常中央区奇异积分数值求积公式

为简化逆Stokes法和逆Vening-Meinesz法反演中央区重力异常计算过程,提高计算效率,本文采用数值求积公式,分别利用Simpson公式和Cotes公式对逆Stokes法和逆Vening-Meinesz法中的奇异积分问题进行了新的研究,系统地推导出了中央区重力异常普适数值积分计算公式.在大地水准面高和垂线偏差理论模型下的分析表明,此公式可直接利用格网节点处的大地水准面高和垂线偏差计算重力异常值,形式简单,计算效率高,计算精度与解析法计算结果精度相当,可以满足实际应用.研究结果可为高精度卫星测高反演重力异常提供基础理论依据.     

利用垂线偏差计算大地水准面中央区效应的改进方法

为提高利用Molodensky公式反演测高大地水准面中央区效应的精度,视中央区为矩形域,将垂线偏差分量表示成双二次多项式插值形式,引入非奇异变换,推导出了大地水准面的计算公式。垂线偏差理论模型下的分析表明本文导出公式误差为零,而传统公式的误差与纬度以及垂线偏差子午分量与卯酉分量之间的比值有关;以中纬度区域分辨率为2'*2'的垂线偏差数据为背景场进行了实际计算,结果表明在反演计算点本身所在的1个网格对大地水准面的贡献时,传统公式与本文导出公式计算结果差值的最大值达数厘米。本文导出公式可为测高大地水准面的高精度反演提供理论依据。     

基于逆Stokes公式的测高重力反演中央区效应计算

为提高利用逆Stokes公式反演测高重力的精度,将中央区大地水准面高表示成双三次多项式插值形式,引入了非奇异变换,推导出了重力异常的计算公式。大地水准面高理论模型下的分析表明,该公式有较高的精度。以分辨率为2′×2′的大地水准面高数据为背景场进行了实际计算,结果说明中央区对反演重力异常有不容忽视的贡献。本文导出的公式可为高精度重力异常的反演提供理论依据。     

论地形垂线偏差中央区贡献的计算

本文根据Ｎｏｗｔｏｎ‘ｓ引力公式，推导出地形垂线偏差的求积公式，通过奇异积分变量变换，推出了一些新的地形垂线偏差计算公式，这些新的计算公式，有效的处理了中央区的奇异积分问题。     

基于球谐函数的重力异常和垂线偏差误差匹配关系

重力异常和垂线偏差是测高卫星非常重要的产品。二者的精度指标对于未来的测高卫星方案设计至关重要。本文利用球谐函数来对重力异常和垂线偏差的精度指标进行讨论,首先从理论上推导了重力异常和垂线偏差误差的近似匹配关系,然后通过6个超高阶重力场模型验证了有关结论的正确性。数值试验表明:垂线偏差误差和重力异常误差满足近似的比例关系,即若垂线偏差各方位向等精度测量,且假定精度均为1μrad,则所对应的重力异常精度约为1.4mGal;反之,若重力异常的精度为1mGal,则所对应的垂线偏差的精度约为0.7μrad。     

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

从莫洛金斯基(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维卷积严格公式,并做了比较和自校核。     

西太平洋海域卫星测高重力垂直梯度分布

从地球重力场基本理论出发,推导了重力垂直梯度计算的数值积分公式,并利用由多颗测高卫星联合反演的2'×2'海洋重力异常资料,解算了西太平洋海域重力垂直梯度.将计算所得的重力异常垂直梯度与现有资料进行比较,其差值均方根为±10.08E,表明两者精度相当.此外,还对差值大小和空间分布进行了深入分析.     

全球及局部海洋扰动重力反演的快速解析方法

从经典边值问题理论及球谐函数理论出发,在空域推导获得了由大地水准面高以及垂线偏差计算扰动重力的解析计算公式,为利用卫星测高数据反演海洋扰动重力提供了理论基础。针对全球海洋区域和局部海洋区域的扰动重力反演,在前人已有工作基础上,提出了改进的基于一维FFT的精确快速算法,保证了计算结果与原解析方法完全一致,且计算速度提高约20倍。该算法在提高计算效率的同时避免了由于引入FFT而产生的混叠、边缘效应问题,而且对观测数据的序列长度没有硬性要求,使得应用更加灵活。利用EGM2008地球重力场模型分别生成了2.5'分辨率大地水准面高数据和垂线偏差数据,按照本文提出的改进方法(采用全球积分计算)分别反演获得了全球及局部海洋区域的扰动重力。经比较分析,由大地水准面和垂线偏差分别反演获得的扰动重力其差异在0.8×10^-5 m/s²以内,这说明两种反演方法是基本一致的,但在数据包含系统误差的情况下,由垂线偏差反演扰动重力具有一定优势。     

用卫星测高资料估算洋区的垂线偏差和重力异常

本文讨论了利用卫星测高资料估算洋区垂线偏差和重力异常问题。使用两种方法，一是最小二乘拟合法，这时，把高程异常展开为一多项式级数，然后由司托克斯反解公式及相应的微分关系式求出所需之值，另一是最小二乘拟推方法，其中所需的协方差函数由不同类型的洋区资料统计得出。选择了二个洋区并按上述方法进行了计算。最后把计算结果与由重力法所得结果进行了比较。     

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

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

Inverse Vening Meinesz formula and deflection-geoid formula: applications to the predictions of gravity and geoid over the South China Sea

Using the spherical harmonic representations of the earth's disturbing potential and its functionals, we derive the inverse Vening Meinesz formula, which converts deflection of the vertical to gravity anomaly using the gradient of the H function. The deflection-geoid formula is also derived that converts deflection to geoidal undulation using the gradient of the C function. The two formulae are implemented by the 1D FFT and the 2D FFT methods. The innermost zone effect is derived. The inverse Vening Meinesz formula is employed to compute gravity anomalies and geoidal undulations over the South China Sea using deflections from Seasat, Geosat, ERS-1 and TOPEX//POSEIDON satellite altimetry. The 1D FFT yields the best result of 9.9-mgal rms difference with the shipborne gravity anomalies. Using the simulated deflections from EGM96, the deflection-geoid formula yields a 4-cm rms difference with the EGM96-generated geoid. The predicted gravity anomalies and geoidal undulations can be used to study the tectonic structure and the ocean circulations of the South China Sea. Received: 7 April 1997 / Accepted: 7 January 1998     

根据多卫星高度计海面高数据推算南中国海及菲律宾海域重力异常(英文)

t Gravity anomalies on a2.5 ×2.5 arc-minute grid in a non-tidal system were derived over the South China and Philippine Seas from multi-satellite altimetry data. North and east components of deflections of the vertical were computed from altimeter-derived sea surface heights at crossover locations, and gridded onto a 2.5 × 2.5 arc-minute resolution grid. EGM96-derived components of deflections of the vertical and gravity anomalies gridded into 2.5 × 2.5 arc-minute resolutions were then used as reference global geopotential model quantities in a remove-restore procedure to implement the Inverse Vening Meinesz formula via the 1D-FFT technique to predict the gravity anomalies over the South China and Philippine Seas from the gridded altimeter-derived components of deflections of the vertical. Statistical comparisons between the altimeter-derived and the shipboard gravity anomalies showed that there is a root-mean-square agreement of 5.7 mgals between them.     

Analysis of some systematic errors affecting altimeter-derived sea surface gradient with application to geoid determination over Taiwan

This paper analyzes several systematic errors affecting sea surface gradients derived from Seasat, Geosat/ERM, Geosat/GM, ERS-1/35d, ERS-1/GM and TOPEX/POSEIDON altimetry. Considering the data noises, the conclusion is: (1) only Seasat needs to correct for the non-geocentricity induced error, (2) only Seasat and Geosat/GM need to correct for the one cycle per revolution error, (3) only Seasat, ERS-1/GM and Geosat/GM need to correct for the tide model error; over shallow waters it is suggested to use a local tide model not solely from altimetry. The effects of the sea surface topography on gravity and geoid computations from altimetry are significant over areas with major oceanographic phenomena. In conclusion, sea surface gradient is a better data type than sea surface height. Sea surface gradients from altimetry, land gravity anomalies, ship gravity anomalies and elevation data were then used to calculate the geoid over Taiwan by least-squares collocation. The inclusion of sea surface gradients improves the geoid prediction by 27% when comparing the GPS-derived and the predicted geoidal heights, and by 30% when comparing the observed and the geoid-derived deflections of the vertical. The predicted geoid along coastal areas is accurate to 2 cm and can help GPS to do the third-order leveling. Received 22 January 1996; Accepted 13 September 1996     

Relation between geoidal undulation, deflection of the vertical and vertical gravity gradient revisited

The vertical gradients of gravity anomaly and gravity disturbance can be related to horizontal first derivatives of deflection of the vertical or second derivatives of geoidal undulations. These are simplified relations of which different variations have found application in satellite altimetry with the implicit assumption that the neglected terms—using remove-restore—are sufficiently small. In this paper, the different simplified relations are rigorously connected and the neglected terms are made explicit. The main neglected terms are a curvilinear term that accounts for the difference between second derivatives in a Cartesian system and on a spherical surface, and a small circle term that stems from the difference between second derivatives on a great and small circle. The neglected terms were compared with the dynamic ocean topography (DOT) and the requirements on the GOCE gravity gradients. In addition, the signal root-mean-square (RMS) of the neglected terms and vertical gravity gradient were compared, and the effect of a remove-restore procedure was studied. These analyses show that both neglected terms have the same order of magnitude as the DOT gradient signal and may be above the GOCE requirements, and should be accounted for when combining altimetry derived and GOCE measured gradients. The signal RMS of both neglected terms is in general small when compared with the signal RMS of the vertical gravity gradient, but they may introduce gradient errors above the spherical approximation error. Remove-restore with gravity field models reduces the errors in the vertical gravity gradient, but it appears that errors above the spherical approximation error cannot be avoided at individual locations. When computing the vertical gradient of gravity anomaly from satellite altimeter data using deflections of the vertical, the small circle term is readily available and can be included. The direct computation of the vertical gradient of gravity disturbance from satellite altimeter data is more difficult than the computation of the vertical gradient of gravity anomaly because in the former case the curvilinear term is needed, which is not readily available.     

球近似下地球外空间任意类型场元的地形影响

传统的重力归算方法只适用于地球表面上的重力异常，不能用于扰动重力、垂线偏差、重力梯度等其他类型扰动重力场元，不适合处理除地面外其他高度上场元的地形影响问题。当前，地球重力场探测的场元类型越来越丰富，探测的高度也逐渐转向航空和卫星高度，精确处理地球外空间各种类型重力场元的地形影响已成为地球重力场领域面临的重要课题。本文通过直接分解由地形生成的具有调和性质的引力场，从而导出地球外空间任意高度、任意类型扰动重力场元的地形影响，在此基础上给出在球近似下地形影响的严密算法和高精度快速算法。利用本文推荐的地形影响计算方案，可以方便地处理各种类型地面重力、海洋重力、航空重力、卫星重力、卫星测高数据的地形影响，从而丰富重力场数据处理的内涵，改善地球重力场算法的性能。     

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

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

Determination of deflections of the vertical using a combination of spherical harmonics and gravimetric data for the area of Greece

The determination of gravimetric deflections of the vertical for the area of Greece is attempted by combining a spherical hamonics model and gravity nomalies using the method of least squares collocation. The components of deflections of the vertical are estimated on a grid with spacing 15′ in latitude and 20′ in longitude covering only the continental area of Greece, where a sufficient number of point gravity anomalies is available. In order to test the accuracy of the determination, gravimetric deflections of the vertical are computed at stations where astrogeodetic data are available. The results show that in a large region of rugged topography and irregular potential field, the prediction is possible with a standard deviation of 18% ... 28% of the root mean square variation of the observations, without taking into account the topography. Furthermore, the estimation of some systematic differences between observed and computed deflections of the vertical is attempted.