首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 781 毫秒
1.
The accuracy of the gravity field approximation depends on the amount of the available data and their distribution as well as on the variation of the gravity field. The variation of the gravity field in the Greek mainland, which is the test area in this study, is very high (the variance of point free air gravity anomalies is 3191.5mgal 2). Among well known reductions used to smooth the gravity field, the complete isostatic reduction causes the best possible smoothing, however remain strong local anomalies which disturb the homogeneity of the gravity field in this area. The prediction of free air gravity anomalies using least squares collocation and regional covariance function is obtained within a ±4 ... ±19mgal accuracy depending on the local peculiarities of the free air gravity field. By taking into account the topography and its isostatic compensation with the usual remove-restore technique, the accuracy of the prediction mentioned obove was increased by about a factor of 4 and the prediction results become quite insensitive to the covariance function used (local or regional). But when predicting geoidal heights, in spite of using the smoothed field, the prediction results remain still depend on the covariance function used in such a way that differences up to about 50cm/100km result between relative geoidal heights computed with regional or local covariance functions.  相似文献   

2.
The astrogeodetic—gravimetric method based on the principle of least—squares solution has been used to determine the geocentric Indian geodetic datum making use of the available nongeocentric astrogeodetic data and the gravimetric geocentric geoidal heights in the form of smoothened values. Everett's method of interpolation has been used to obtain the smoothened geoidal heights at the astrogeodetic stations in India from the available generalized values at 1°×1° corners. The values of the geoidal height and deflections of the vertical at the geodetic datum Kalianpur H.S. so obtained have the negligible difference from the values computed earlier by the same method using directly computed gravimetric geoidal heights at the astrogeodetic stations, indicating that the use of the interpolated values in the astrogeodetic—gravimetric method employed would be an economical approach of absolute orientation of a nongeocentric system if the gravimetric geoidal heights are available at 1°×1° corners in the area of interest.  相似文献   

3.
The investigations refer to the compartment method by using mean terrestrial free air anomalies only. Three main error influences of remote areas (distance from the fixed point >9°) on height anomalies and deflections of the vertical are being regarded:
  1. The prediction errors of mean terrestrial free air anomalies have the greatest influence and amount to about ±0″.2 in each component for deflections of the vertical and to ±3 m for height anomalies;
  2. The error of the compartment method, which originates from converting the integral formulas of Stokes and Vening-Meinesz into summation formulas, can be neglected if the anomalies for points and gravity profiles are compiled to 5°×5° mean values.
  3. The influences of the mean gravimetric correction terms of Arnold—estimated for important mountains of the Earth by means of an approximate formula—on height anomalies may amount to 1–2 m and on deflections of the vertical to 0″0.5–0″.1, and, therefore, they have to be taken into account for exact calculations.
The computations of errors are carried out using a global covariance function of point free air anomalies.  相似文献   

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

5.
Terrestrial free-air gravity anomalies form a most essential data source in the framework of gravity field determination. Gravity anomalies depend on the datums of the gravity, vertical, and horizontal networks as well as on the definition of a normal gravity field; thus gravity anomaly data are affected in a systematic way by inconsistencies of the local datums with respect to a global datum, by the use of a simplified free-air reduction procedure and of different kinds of height system. These systematic errors in free-air gravity anomaly data cause systematic effects in gravity field related quantities like e.g. absolute and relative geoidal heights or height anomalies calculated from gravity anomaly data. In detail it is shown that the effects of horizontal datum inconsistencies have been underestimated in the past. The corresponding systematic errors in gravity anomalies are maximum in mid-latitudes and can be as large as the errors induced by gravity and vertical datum and height system inconsistencies. As an example the situation in Australia is evaluated in more detail: The deviations between the national Australian horizontal datum and a global datum produce a systematic error in the free-air gravity anomalies of about −0.10 mgal which value is nearly constant over the continent  相似文献   

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

7.
The evaluation of deflections of the vertical for the area of Greece is attempted using a combination of topographic and astrogeodetic data. Tests carried out in the area bounded by 35°≤ϕ≤42°, 19°≤λ≤27° indicate that an accuracy of ±3″.3 can be obtained in this area for the meridian and prime vertical deflection components when high resolution topographic data in the immediate vicinity of computation points are used, combined with high degree spherical harmonic expansions of the geopotential and isostatic reduction potential. This accuracy is about 25% better than the corresponding topographic-Moho deflection components which are evaluated using topographic and Moho data up to 120 km around each station, without any combination with the spherical harmonic expansion of the geopotential or isostatic reduction potential. The accuracy in both cases is increased to about 2″.6 when the astrogeodetic data available in the area mentioned above are used for the prediction of remaining values. Furthermore the estimation of datum-shift parameters is attempted using least squares collocation.  相似文献   

8.
魏子卿 《测绘学报》2022,51(6):797-803
在空间大地测量时代,GNSS可以测定地面点的大地高,使重力扰动变成了直接观测量,以重力扰动为边界条件的第二边值问题在大地测量中得以实用化。它的解与GNSS组合正在成为一种颇有应用前景的海拔高测量方法。本文原理性地讨论了有两种不同边界面的球近似第二大地边值问题。第一种以地形面为边界面,给出了高程异常与地面垂线偏差的解析延拓解;第二种以参考椭球面为边界面,将其外部地形质量按照Helmert第二压缩法移至参考椭球面,然后将Hotine函数与从地球表面延拓至边界面的Helmert重力扰动进行卷积,并顾及地形间接影响,最后得到大地水准面高、椭球面垂线偏差、高程异常与地面垂线偏差的Helmert解。在讨论部分,进行了第二与第三大地边值问题的比较,提出了现有重力点高程从正高或正常高到大地高的改化方法,并展望了它的应用前景。  相似文献   

9.
A detailed gravimetric geoid in the North Atlantic Ocean, named DGGNA-77, has been computed, based on a satellite and gravimetry derived earth potential model (consisting in spherical harmonic coefficients up to degree and order 30) and mean free air surface gravity anomalies (35180 1°×1° mean values and 245000 4′×4′ mean values). The long wavelength undulations were computed from the spherical harmonics of the reference potential model and the details were obtained by integrating the residual gravity anomalies through the Stokes formula: from 0 to 5° with the 4′×4′ data, and from 5° to 20° with the 1°×1° data. For computer time reasons the final grid was computed with half a degree spacing only. This grid extends from the Gulf of Mexico to the European and African coasts. Comparisons have been made with Geos 3 altimetry derived geoid heights and with the 5′×5′ gravimetric geoid derived byMarsh andChang [8] in the northwestern part of the Atlantic Ocean, which show a good agreement in most places apart from some tilts which porbably come from the satellite orbit recovery.  相似文献   

10.
The method of analytical downward continuation has been used for solving Molodensky’s problem. This method can also be used to reduce the surface free air anomaly to the ellipsoid for the determination of the coefficients of the spherical harmonic expansion of the geopotential. In the reduction of airborne or satellite gradiometry data, if the sea level is chosen as reference surface, we will encounter the problem of the analytical downward continuation of the disturbing potential into the earth, too. The goal of this paper is to find out the topographic effect of solving Stoke’sboundary value problem (determination of the geoid) by using the method of analytical downward continuation. It is shown that the disturbing potential obtained by using the analytical downward continuation is different from the true disturbing potential on the sea level mostly by a −2πGρh 2/p. This correction is important and it is very easy to compute and add to the final results. A terrain effect (effect of the topography from the Bouguer plate) is found to be much smaller than the correction of the Bouguer plate and can be neglected in most cases. It is also shown that the geoid determined by using the Helmert’s second condensation (including the indirect effect) and using the analytical downward continuation procedure (including the topographic effect) are identical. They are different procedures and may be used in different environments, e.g., the analytical downward continuation procedure is also more convenient for processing the aerial gravity gradient data. A numerical test was completed in a rough mountain area, 35°<ϕ<38°, 240°<λ<243°. A digital height model in 30″×30″ point value was used. The test indicated that the terrain effect in the test area has theRMS value ±0.2−0.3 cm for geoid. The topographic effect on the deflections of the vertical is around1 arc second.  相似文献   

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

12.
The AUSGeoid09 model of the Australian Height Datum   总被引：8，自引：6，他引：2  
AUSGeoid09 is the new Australia-wide gravimetric quasigeoid model that has been a posteriori fitted to the Australian Height Datum (AHD) so as to provide a product that is practically useful for the more direct determination of AHD heights from Global Navigation Satellite Systems (GNSS). This approach is necessary because the AHD is predominantly a third-order vertical datum that contains a ~1 m north-south tilt and ~0.5 m regional distortions with respect to the quasigeoid, meaning that GNSS-gravimetric-quasigeoid and AHD heights are inconsistent. Because the AHD remains the official vertical datum in Australia, it is necessary to provide GNSS users with effective means of recovering AHD heights. The gravimetric component of the quasigeoid model was computed using a hybrid of the remove-compute-restore technique with a degree-40 deterministically modified kernel over a one-degree spherical cap, which is superior to the remove-compute-restore technique alone in Australia (with or without a cap). This is because the modified kernel and cap combine to filter long-wavelength errors from the terrestrial gravity anomalies. The zero-tide EGM2008 global gravitational model to degree 2,190 was used as the reference field. Other input data are ~1.4 million land gravity anomalies from Geoscience Australia, 1′ × 1′ DNSC2008GRA altimeter-derived gravity anomalies offshore, the 9′′ × 9′′ GEODATA-DEM9S Australian digital elevation model, and a readjustment of Australian National Levelling Network (ANLN) constrained to the CARS2006 mean dynamic ocean topography model. To determine the numerical integration parameters for the modified kernel, the gravimetric component of AUSGeoid09 was compared with 911 GNSS-observed ellipsoidal heights at benchmarks. The standard deviation of fit to the GNSS-AHD heights is ±222 mm, which dropped to ±134 mm for the readjusted GNSS-ANLN heights showing that careful consideration now needs to be given to the quality of the levelling data used to assess gravimetric quasigeoid models. The publicly released version of AUSGeoid09 also includes a geometric component that models the difference between the gravimetric quasigeoid and the zero surface of the AHD at 6,794 benchmarks. This a posteriori fitting used least-squares collocation (LSC) in cross-validation mode to determine a correlation length of 75 km for the analytical covariance function, whereas the noise was taken from the estimated standard deviation of the GNSS ellipsoidal heights. After this LSC surface fitting, the standard deviation of fit reduced to ±30 mm, one-third of which is attributable to the uncertainty in the GNSS ellipsoidal heights.  相似文献   

13.
A calculation of quasigeoidal heights and plumb-line deflections according to Molodensky formulae was carried out under elimination of the effect of topography from gravity anomalies. After the masses of topography had been removed a smoothed-out surface passing through astronomical and gravity stations was considered as representing the physical surface of the Earth. Thus it has been practically rendered possible to use the first-approximation formulae of Molodensky, and, in many cases, also the “zero-approximation” formulae analogous to the formulae of Stokes and Vening-Meinesz. The effect of the restored masses of topography was then added to the quantities found; the said effect was expressed as the effect of topography condensed on the normal equipotential surface passing through the point under investigation, plus a correction for condensation. Following some transformations, the resulting formulae (13) and (18) were obtained which formulae differ in their “zero-approximation” (15) and (20) from traditional formulas in that they contain terrait reductions added to free-air anomalies. Moreover, in the calculation of plumb-line deflections directly in mountain regions a correction for differing effects of topography before and after its condensation is to be introduced. A tentative expansion of terrain reduction in terms of spherical harmonics up to the third order is given; it can be seen therefrom that the Stokes series in its usual form is subject to a mean arror about 15–20%. It is also shown that the expansion of free-air anomalies in terms of spherical functions contains a first-order harmonic with a mean values about ±0.3 mgl. The said harmonic practically disappears in the expansion of the sum of free-air anomalies and terrain reductions.  相似文献   

14.
Summary Using a data set of 260 000 gravity anomalies it is shown that common characteristics for a local covariance function exist in an area as large as Canada excluding the Rocky Mountains. After eliminating global features by referencing the data to the GEM-10 satellite solution, the shape of the covariance function is remarkably consistent from one sample area to the next. The determination of the essential parameters and the fitting of the covariance function are discussed in detail. To test the reliability of the derived function, deflections of the vertical are estimated at about 230 stations where astrogeodetic data are available. Results show that the standard error obtained from the discrepancies is about1″ for each component and that the error covariance matrix of least-squares collocation reflects this accuracy remarkably well.  相似文献   

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

16.
 Four different implementations of Stokes' formula are employed for the estimation of geoid heights over Sweden: the Vincent and Marsh (1974) model with the high-degree reference gravity field but no kernel modifications; modified Wong and Gore (1969) and Molodenskii et al. (1962) models, which use a high-degree reference gravity field and modification of Stokes' kernel; and a least-squares (LS) spectral weighting proposed by Sj?berg (1991). Classical topographic correction formulae are improved to consider long-wavelength contributions. The effect of a Bouguer shell is also included in the formulae, which is neglected in classical formulae due to planar approximation. The gravimetric geoid is compared with global positioning system (GPS)-levelling-derived geoid heights at 23 Swedish Permanent GPS Network SWEPOS stations distributed over Sweden. The LS method is in best agreement, with a 10.1-cm mean and ±5.5-cm standard deviation in the differences between gravimetric and GPS geoid heights. The gravimetric geoid was also fitted to the GPS-levelling-derived geoid using a four-parameter transformation model. The results after fitting also show the best consistency for the LS method, with the standard deviation of differences reduced to ±1.1 cm. For comparison, the NKG96 geoid yields a 17-cm mean and ±8-cm standard deviation of agreement with the same SWEPOS stations. After four-parameter fitting to the GPS stations, the standard deviation reduces to ±6.1 cm for the NKG96 geoid. It is concluded that the new corrections in this study improve the accuracy of the geoid. The final geoid heights range from 17.22 to 43.62 m with a mean value of 29.01 m. The standard errors of the computed geoid heights, through a simple error propagation of standard errors of mean anomalies, are also computed. They range from ±7.02 to ±13.05 cm. The global root-mean-square error of the LS model is the other estimation of the accuracy of the final geoid, and is computed to be ±28.6 cm. Received: 15 September 1999 / Accepted: 6 November 2000  相似文献   

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

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

19.
Conclusion If we want to compute the height anomalies and the deflections of the vertical for a point where there are height differences exceeding one kilometer in the neighborhood, a very dense net of gravity values must be observed, even when the topographic corrections are used. In smoother regions, the simple Bouguer anomalies with a moderate spacing and with the estimation of the mean heights give reliable results without further reductions.  相似文献   

20.
Any errors in digital elevation models (DEMs) will introduce errors directly in gravity anomalies and geoid models when used in interpolating Bouguer gravity anomalies. Errors are also propagated into the geoid model by the topographic and downward continuation (DWC) corrections in the application of Stokes’s formula. The effects of these errors are assessed by the evaluation of the absolute accuracy of nine independent DEMs for the Iran region. It is shown that the improvement in using the high-resolution Shuttle Radar Topography Mission (SRTM) data versus previously available DEMs in gridding of gravity anomalies, terrain corrections and DWC effects for the geoid model are significant. Based on the Iranian GPS/levelling network data, we estimate the absolute vertical accuracy of the SRTM in Iran to be 6.5 m, which is much better than the estimated global accuracy of the SRTM (say 16 m). Hence, this DEM has a comparable accuracy to a current photogrammetric high-resolution DEM of Iran under development. We also found very large differences between the GLOBE and SRTM models on the range of −750 to 550 m. This difference causes an error in the range of −160 to 140 mGal in interpolating surface gravity anomalies and −60 to 60 mGal in simple Bouguer anomaly correction terms. In the view of geoid heights, we found large differences between the use of GLOBE and SRTM DEMs, in the range of −1.1 to 1 m for the study area. The terrain correction of the geoid model at selected GPS/levelling points only differs by 3 cm for these two DEMs.  相似文献   

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

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