首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 195 毫秒
1.
利用卫星测高数据反演海洋重力异常研究   总被引:20,自引:2,他引:20  
全面研究了利用卫得测高数据反演海洋重力异常3种主要方法(即Stokes数据解析反解以及逆Vening-Meinesz公式)的技术特点,建立了3种算法的数学模型及其谱计算式,在以1440阶次位模型定义的标准场中完成了3种算法的数值比较和内部检核,通过仿真试验实现了3种算法的可靠性和稳定性检验,最后,本文利用卫得测高实测对南中国海地区的海洋重力异常进行了实际反演,并将反演结果同船测数据进行了比较。  相似文献   

2.
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.  相似文献   

3.
Inverse problems in isostasy will consist in making the isostatic anomalies to be zero under a certain isostatic hypothesis. In the case of the Vening Meinesz isostatic hypothesis, the density contrast is constant, while the Moho depth (depth of the Mohorovičić discontinuity) is variable. Hence, the Vening Meinesz inverse isostatic problem aims to determine a suitable variable Moho depth for a prescribed constant density contrast. The main idea is easy but the theoretical analysis is somewhat difficult. Moreover, the practical determination of the variable Moho depths based on the Vening Meinesz inverse problem is a laborious and time-consuming task. The formulas used for computing the inverse Vening Meinesz Moho depths are derived. The computational tricks essentially needed for computing the inverse Vening Meinesz Moho depths from a set of local and global Bouguer anomalies are described. The Moho depths for a test area are computed based on the inverse Vening Meinesz isostatic problem. These Moho depths fit the Moho depths derived from seismic observations with a good accuracy, in which the parameters used for the fitting agree well with those determined geophysically. Received: 4 February 1999 / Accepted: 4 October 1999  相似文献   

4.
Comparisons of gravimetric and astrogeodetic deflections of the vertical in the Australian region indicate that the former are affected by position dependent systematic errors, even after orientation onto the Australian Geodetic Datum. These are probably due to errors in the predicted mean anomalies for gravimetrically unsurveyed oceanic regions to the east, south and west of the continent. Deflection component residuals (astrogeodetic minus oriented gravimetric) at 83 control stations are made the observables in a set of observation equations, based on the Vening Meinesz equations, from which pseudocorrections to the mean anomalies for a set of arbitrarily selected surface elements are computed. These pseudocorrections compensate for prediction errors in much larger unsurveyed regions. Their effects on individual deflection components are calculated using the Vening Meinesz equations. Statistical tests indicate that pseudocorrections computed for four large offshore elements and six smaller elements in unsurveyed areas produce corrections to the gravimetric deflections which make the ξ and η components in seconds of arc consistent with normally distributed populations N (0.00, 0.702).  相似文献   

5.
A general formula giving Molodenskii coefficientsQ n of the truncation errors for the geoidal height is introduced in this paper. A relation betweenQ n andq n, Cook’s truncation function, is also obtained. Cook (1951) has treated the truncation errors for the deflection of the vertical in the Vening Meinesz integration. Molodenskii et al. (1962) have also derived the truncation error formulas for the deflection of the vertical. It is proved in this paper that these two formulas are equivalent.  相似文献   

6.
1 IntroductionThefastFouriertransform (FFT)techniqueisaverypowerfultoolfortheefficientevaluationofgravityfieldconvolutionintegrals.Thankstothegoodcomputationefficiency ,theFFTtechnique ,inthemid_1 980s ,begantofindwidespreaduseingeoiddetermination ,whencompar…  相似文献   

7.
There exist three types of convolution formulae for the efficient evaluation of gravity field convolution integrals, i.e., the planar 2D convolution, the spherical 2D convolution and the spherical 1D convolution. The largest drawback of both the planar and the spherical 2D FFT methods is that, due to the approximations in the kernel function, only inexact results can be achieved. Apparently, the reason is the meridian convergence at higher latitudes. As the meridians converge, the ??,?λ blocks do not form a rectangular grid, as is assumed in 2D FFT methods. It should be pointed out that the meridian convergence not only leads to an approximation error in the kernel function, but also causes an approximation error during the implementation of 2D FFT in computer. In order to meet the increasing need for precise determination of the vertica deflections, this paper derives a more precise planar 2D FFT formula for the computation of the vertical deflections. After having made a detailed comparison between the planar and the spherical 2D FFT formulae, we find out the main source of errors causing the loss in accuracy by applying the conventional spherical 2D FFT method. And then, a modified spherical 2D FFT formula for the computation of the vertical deflections is developed in this paper. A series of numerical tests have been carried out to illustrate the improvement made upon the old spherical 2D FFT. The second part of this paper is to discuss the influences of the spherical harmonic reference field, the limited capsize, and the singular integral on the computation of the vertical deflections. The results of the vertical deflections over China by applying the spherical 1D FFT formula with different integration radii have been compared to the astro-observed vertical deflections in the South China Sea to obtain a set of optimum deflection computation parameters.  相似文献   

8.
 Global mean sea surface heights (SSHs) and gravity anomalies on a 2×2 grid were determined from Seasat, Geosat (Exact Repeat Mission and Geodetic Mission), ERS-1 (1.5-year mean of 35-day, and GM), TOPEX/POSEIDON (T/P) (5.6-year mean) and ERS-2 (2-year mean) altimeter data over the region 0–360 longitude and –80–80 latitude. To reduce ocean variabilities and data noises, SSHs from non-repeat missions were filtered by Gaussian filters of various wavelengths. A Levitus oceanic dynamic topography was subtracted from the altimeter-derived SSHs, and the resulting heights were used to compute along-track deflection of the vertical (DOV). Geoidal heights and gravity anomalies were then computed from DOV using the deflection-geoid and inverse Vening Meinesz formulae. The Levitus oceanic dynamic topography was added back to the geoidal heights to obtain a preliminary sea surface grid. The difference between the T/P mean sea surface and the preliminary sea surface was computed on a grid by a minimum curvature method and then was added to the preliminary grid. The comparison of the NCTU01 mean sea surface height (MSSH) with the T/P and the ERS-1 MSSH result in overall root-mean-square (RMS) differences of 5.0 and 3.1 cm in SSH, respectively, and 7.1 and 3.2 μrad in SSH gradient, respectively. The RMS differences between the predicted and shipborne gravity anomalies range from 3.0 to 13.4 mGal in 12 areas of the world's oceans. Received: 26 September 2001 / Accepted: 3 April 2002 Correspondence to: C. Hwang Acknowledgements. This research is partly supported by the National Science Council of ROC, under grants NSC89-2611-M-009-003-OP2 and NSC89-2211-E-009-095. This is a contribution to the IAG Special Study Group 3.186. The Geosat and ERS1/2 data are from NOAA and CERSAT/France, respectively. The T/P data were provided by AVISO. The CLS and GSFC00 MSS models were kindly provided by NASA/GSFC and CLS, respectively. Drs. Levitus, Monterey, and Boyer are thanked for providing the SST model. Dr. T. Gruber and two anonymous reviewers provided very detailed reviews that improved the quality of this paper.  相似文献   

9.
A general formula is developed and presented for transformations among geoidal undulation, gravity anomaly, gravity disturbance and other gravimetric quantities. Using a spectral form of the general formula, a criterion has been built in order to classify these transformations into forward and inverse transformations in this paper. Then, the two-dimensional convolution techniques are applied to the general formula to deal with the forward transformation while the two-dimensional deconvolution techniques are employed to treat the inverse transformation and evaluate the inverse general formula. Concepts of convolution and deconvolution are also reviewed in this paper. The stability and edge effect problems related to the deconvolution techniques are investigated using simulated data and numerical tests are done to quantify the stability of the deconvolution techniques for estimated gravity information. Finally, the marine gravity information for the Norwegian-Greenland Sea area has been derived from ERS-1 altimetry data using the deconvolution techniques.  相似文献   

10.
The accumulation of good quality satellite altimetry missions allows us to have a precise geoid with fair resolution and to compute free air gravity anomalies easily by fast Fourier transform (FFT) techniques.In this study we are comparing two methods to get gravity anomalies. The first one is to establish a geoid grid and transform it into anomalies using inverse Stokes formula in the spectral domain via FFT. The second one computes deflection of the vertical grids and transforms them into anomalies.The comparison is made using different data sets: Geosat, ERS-1 and Topex-Poseidon exact repeat misions (ERMs) north of 30°S and Geosat geodetic mission (GM) south of 30°S. The second method which transforms the geoid gradients converted into deflection of the vertical values is much better and the results have been favourably evaluated by comparison with marine gravity data.  相似文献   

11.
This study demonstrates that in mountainous areas the use of residual terrain model (RTM) data significantly improves the accuracy of vertical deflections obtained from high-degree spherical harmonic synthesis. The new Earth gravitational model EGM2008 is used to compute vertical deflections up to a spherical harmonic degree of 2,160. RTM data can be constructed as difference between high-resolution Shuttle Radar Topography Mission (SRTM) elevation data and the terrain model DTM2006.0 (a spherical harmonic terrain model that complements EGM2008) providing the long-wavelength reference surface. Because these RTM elevations imply most of the gravity field signal beyond spherical harmonic degree of 2,160, they can be used to augment EGM2008 vertical deflection predictions in the very high spherical harmonic degrees. In two mountainous test areas—the German and the Swiss Alps—the combined use of EGM2008 and RTM data was successfully tested at 223 stations with high-precision astrogeodetic vertical deflections from recent zenith camera observations (accuracy of about 0.1 arc seconds) available. The comparison of EGM2008 vertical deflections with the ground-truth astrogeodetic observations shows root mean square (RMS) values (from differences) of 3.5 arc seconds for ξ and 3.2 arc seconds for η, respectively. Using a combination of EGM2008 and RTM data for the prediction of vertical deflections considerably reduces the RMS values to the level of 0.8 arc seconds for both vertical deflection components, which is a significant improvement of about 75%. Density anomalies of the real topography with respect to the residual model topography are one factor limiting the accuracy of the approach. The proposed technique for vertical deflection predictions is based on three publicly available data sets: (1) EGM2008, (2) DTM2006.0 and (3) SRTM elevation data. This allows replication of the approach for improving the accuracy of EGM2008 vertical deflection predictions in regions with a rough topography or for improved validation of EGM2008 and future high-degree spherical harmonic models by means of independent ground truth data.  相似文献   

12.
S. Ono 《Journal of Geodesy》1985,59(3):275-288
In order to solve the problems of determining the shape of a part of the earth of national or continental extent, that is, of rigorous constituting and computing of the astrogeodetic network, it is required to determine gravimetric deflections of the vertical with an accuracy of, say, 0″.3. For this it is adequate to carry out additional gravity surveys in the neighborhoods of computation points, in addition to a given uniform gravity survey (normal density gravity survey). The study offers a method to determine the optimal distribution of gravity stations in such a gravity survey, which guarantees a given accuracy of computed gravimetric deflections of the vertical for a given statistical condition which characterizes the variation of the gravity field. The approach used here is based on the concept of the error of representation and the error propagation of Vening Meinesz integrals.  相似文献   

13.
超高阶地球位模型的计算与分析   总被引:8,自引:1,他引:8  
全球重力场模型是当今物理大地测量学最为活跃的研究领域之一。本文基于目前国内外最新的重力场模型理论研究成果 ,提出了利用中国地区细部数据和全球卫星测高 2′× 2′网格重力异常扩展超高阶位模型的计算方法 ,详细讨论了数值解算过程中的稳定性和可靠性问题。以 EGM96和 GPM98CR模型作为参考模型 ,在全球意义上分别解算得到 MOD99a( 360阶 )、MOD99b( 72 0阶 )和 MOD99c/d( 1 80 0阶 ) ,将系列模型 MOD99a/b/c/d同中国地区 72个 GPS水准大地水准面和全球海洋 1 2个地区的卫星测高大地水准面进行了比较 ,并通过功率谱分析方法检验了 4组模型的有效性和可靠性。  相似文献   

14.
高阶地球重力场模型的评价及其优选   总被引:1,自引:0,他引:1  
采用"移去-恢复"技术确定大地水准面,需要一个全球重力场模型作为参考场。本文采用基于模糊集合理论的最大隶属度原则,通过某试验区的33个GPS水准点的实测大地水准面差距与模型大地水准面差距的比较分析,从国际重力场服务提供的EGM96、EIGEN-CG01C、EIGEN-CG03C、GFZ93A、GFZ96、OSU91A、PGM2000A等七个高阶全球重力场模型中,选择OSU91A作为该地区最优的参考重力场模型。  相似文献   

15.
Mean gravity anomalies, deflections of the vertical, and a geopotential model complete to degree and order180 are combined in order to determine geoidal heights in the area bounded by [34°≦ϕ≤42°, 18°≦λ≦28°]. Moreover, employing point gravity anomalies simultaneously with the above data, an attempt is made to predict deflections of the vertical in the same area. The method used in the computations is least squares collocation. Using empirical covariance functions for the data, the suitable errors for the different sources of observations, and the optimum cap radius around each point of evaluation, an accuracy better than±0.60m for geoidal heights and±1″.5 for deflections of the vertical is obtained taking into account existing systematic effects. This accuracy refers to the comparison between observed and predicted values.  相似文献   

16.
 It is suggested that a spherical harmonic representation of the geoidal heights using global Earth gravity models (EGM) might be accurate enough for many applications, although we know that some short-wavelength signals are missing in a potential coefficient model. A `direct' method of geoidal height determination from a global Earth gravity model coefficient alone and an `indirect' approach of geoidal height determination through height anomaly computed from a global gravity model are investigated. In both methods, suitable correction terms are applied. The results of computations in two test areas show that the direct and indirect approaches of geoid height determination yield good agreement with the classical gravimetric geoidal heights which are determined from Stokes' formula. Surprisingly, the results of the indirect method of geoidal height determination yield better agreement with the global positioning system (GPS)-levelling derived geoid heights, which are used to demonstrate such improvements, than the results of gravimetric geoid heights at to the same GPS stations. It has been demonstrated that the application of correction terms in both methods improves the agreement of geoidal heights at GPS-levelling stations. It is also found that the correction terms in the direct method of geoidal height determination are mostly similar to the correction terms used for the indirect determination of geoidal heights from height anomalies. Received: 26 July 2001 / Accepted: 21 February 2002  相似文献   

17.
A new computational procedure for derivation of marine geoid on a 2.5′×2.5′grid in a non-tidal system over the South China Sea and the Philippine Sea from multi-satellite altimeter sea surface heights is discussed. Single-and dual-satellite crossovers were performed, and components of deflections of the vertical were determined at the crossover positions using Sand-well's computational theory, and gridded onto a 2.5′×2.5′resolution grid by employing the Shepard's interpolation procedure. 2.5′×2.5′grid of EGM96-derived components of deflections of the vertical and geoid heights were then used as reference global geopotential model quantities in a remove-restore procedure to implement the Molodensky-like formula via 1D-FFT technique to predict the geoid heights over the South China Sea and the Philippine Sea from the gridded altimeter-derived components of deflec-tions of the vertical. Statistical comparisons between the altimeter-and the EGM96- derived geoid heights showed that there was a root-mean-square agreement of ±0.35 m between them in a region of less tectonically active geological structures. However, over areas of tectonically active structures such as the Philippine trench, differences of about -19.9 m were obtained.  相似文献   

18.
This paper presents a set of efficient formulas to evaluate the deflections of the vertical on the sphere using gridded data. The Vening-Meinesz formula, the topographic indirect effect on the deflections of the vertical as well as the terrain corrections are expressed as both 2D and 1D convolutions on the sphere, and consequently can be evaluated by the 2D and the 1D fast Fourier transform (FFT). When compared with the results obtained from pointwise integration, the use of the 1D FFT gives identical results, and therefore these results were used as control values in this paper. The use of the spherical 2D FFT improves significantly the computational efficiency with little sacrifice of accuracy (0.6 rms difference from the 1D FFT results). The planar 2D FFT, which is as efficient as the spherical 2D FFT, gives worse results (1.2 rms difference from the 1D FFT results) because of the extra approximations. Received: 27 February 1996; Accepted: 24 January 1997  相似文献   

19.
The theoretical differences between the Helmert deflection of the vertical and that computed from a truncated spherical harmonic series of the gravity field, aside from the limited spectral content in the latter, include the curvature of the normal plumb line, the permanent tidal effect, and datum origin and orientation offsets. A numerical comparison between deflections derived from spherical harmonic model EGM96 and astronomic deflections in the conterminous United States (CONUS) shows that correcting these systematic effects reduces the mean differences in some areas. Overall, the mean difference in CONUS is reduced from −0.219 arcsec to −0.058 arcsec for the south–north deflection, and from +0.016 arcsec to +0.004 arcsec for the west–east deflection. Further analysis of the root-mean-square differences indicates that the high-degree spectrum of the EGM96 model has significantly less power than implied by the deflection data. Received: 9 December 1997 / Accepted: 21 August 1998  相似文献   

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

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号