首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
 The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used, as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads. Received: 14 August 2000 / Accepted: 28 February 2001  相似文献   

2.
 Two numerical techniques are used in recent regional high-frequency geoid computations in Canada: discrete numerical integration and fast Fourier transform. These two techniques have been tested for their numerical accuracy using a synthetic gravity field. The synthetic field was generated by artificially extending the EGM96 spherical harmonic coefficients to degree 2160, which is commensurate with the regular 5′ geographical grid used in Canada. This field was used to generate self-consistent sets of synthetic gravity anomalies and synthetic geoid heights with different degree variance spectra, which were used as control on the numerical geoid computation techniques. Both the discrete integration and the fast Fourier transform were applied within a 6∘ spherical cap centered at each computation point. The effect of the gravity data outside the spherical cap was computed using the spheroidal Molodenskij approach. Comparisons of these geoid solutions with the synthetic geoid heights over western Canada indicate that the high-frequency geoid can be computed with an accuracy of approximately 1 cm using the modified Stokes technique, with discrete numerical integration giving a slightly, though not significantly, better result than fast Fourier transform. Received: 2 November 1999 / Accepted: 11 July 2000  相似文献   

3.
A synthetic [simulated] Earth gravity model (SEGM) of the geoid, gravity and topography has been constructed over Australia specifically for validating regional gravimetric geoid determination theories, techniques and computer software. This regional high-resolution (1-arc-min by 1-arc-min) Australian SEGM (AusSEGM) is a combined source and effect model. The long-wavelength effect part (up to and including spherical harmonic degree and order 360) is taken from an assumed errorless EGM96 global geopotential model. Using forward modelling via numerical Newtonian integration, the short-wavelength source part is computed from a high-resolution (3-arc-sec by 3-arc-sec) synthetic digital elevation model (SDEM), which is a fractal surface based on the GLOBE v1 DEM. All topographic masses are modelled with a constant mass-density of 2,670 kg/m3. Based on these input data, gravity values on the synthetic topography (on a grid and at arbitrarily distributed discrete points) and consistent geoidal heights at regular 1-arc-min geographical grid nodes have been computed. The precision of the synthetic gravity and geoid data (after a first iteration) is estimated to be better than 30 μ Gal and 3 mm, respectively, which reduces to 1 μ Gal and 1 mm after a second iteration. The second iteration accounts for the changes in the geoid due to the superposed synthetic topographic mass distribution. The first iteration of AusSEGM is compared with Australian gravity and GPS-levelling data to verify that it gives a realistic representation of the Earth’s gravity field. As a by-product of this comparison, AusSEGM gives further evidence of the north–south-trending error in the Australian Height Datum. The freely available AusSEGM-derived gravity and SDEM data, included as Electronic Supplementary Material (ESM) with this paper, can be used to compute a geoid model that, if correct, will agree to in 3 mm with the AusSEGM geoidal heights, thus offering independent verification of theories and numerical techniques used for regional geoid modelling.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00190-005-0002-z  相似文献   

4.
When planning a satellite gravity gradiometer (SGG) mission, it is important to know the quality of the quantities to be recovered at ground level as a function of e.g. satellite altitude, data type and sampling rate, and signal variance and noise. This kind of knowledge may be provided either using the formal error estimates of wanted quantities using least-squares collocation (LSC) or by comparing simulated data at ground level with results computed by methods like LSC or Fast Fourier Transform (FFT). Results of a regional gravity field recovery in a 10o×20o area surrounding the Alps using LSC and FFT are reported. Data used as observations in satellite altitude (202 or161 km) and for comparison at ground level were generated using theOSU86F coefficient set, complete to degree 360. These observations are referred to points across simulated orbits. The simulated quantities were computed for a 45 days mission period and 4 s sampling. A covariance function which also included terms above degree 360 was used for prediction and error estimation. This had the effect that the formal error standard deviation for gravity anomalies were considerably larger than the standard deviations of predicted minus simulated quantities. This shows the importance of using data with frequency content above degree 360 in simulation studies. Using data at202 km altitude the standard deviation of the predicted minus simulated data was equal to8.3 mgal for gravity and0.33 m for geoid heights.  相似文献   

5.
Spherical harmonic synthesis (SHS) of gravity field functionals at the Earth’s surface requires the use of heights. The present study investigates the gradient approach as an efficient yet accurate strategy to incorporate height information in SHS at densely spaced multiple points. Taylor series expansions of commonly used functionals quasigeoid heights, gravity disturbances and vertical deflections are formulated, and expressions of their radial derivatives are presented to arbitrary order. Numerical tests show that first-order gradients, as introduced by Rapp (J Geod 71(5):282–289, 1997) for degree 360 models, produce cm- to dm-level RMS approximation errors over rugged terrain when applied with EGM2008 to degree 2190. Instead, higher-order Taylor expansions are recommended that are capable of reducing approximation errors to insignificance for practical applications. Because the height information is separated from the actual synthesis, the gradient approach can be applied along with existing highly efficient SHS routines to compute surface functionals at arbitrarily dense grid points. This confers considerable computational savings (above or well above one order of magnitude) over conventional point-by-point SHS. As an application example, an ultra-high resolution model of surface gravity functionals (EurAlpGM2011) is constructed over the entire European Alps that incorporates height information in the SHS at 12,000,000 surface points. Based on EGM2008 and residual topography data, quasigeoid heights, gravity disturbances and vertical deflections are estimated at ~200m resolution. As a conclusion, the gradient approach is efficient and accurate for high-degree SHS at multiple points at the Earth’s surface.  相似文献   

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 AUSGeoid98 gravimetric geoid model of Australia has been computed using data from the EGM96 global geopotential model, the 1996 release of the Australian gravity database, a nationwide digital elevation model, and satellite altimeter-derived marine gravity anomalies. The geoid heights are on a 2 by 2 arc-minute grid with respect to the GRS80 ellipsoid, and residual geoid heights were computed using the 1-D fast Fourier transform technique. This has been adapted to include a deterministically modified kernel over a spherical cap of limited spatial extent in the generalised Stokes scheme. Comparisons of AUSGeoid98 with GPS and Australian Height Datum (AHD) heights across the continent give an RMS agreement of ±0.364 m, although this apparently large value is attributed partly to distortions in the AHD. Received: 10 March 2000 / Accepted: 21 February 2001  相似文献   

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

9.
Ellipsoidal geoid computation   总被引：1，自引：1，他引：0  
Modern geoid computation uses a global gravity model, such as EGM96, as a third component in a remove–restore process. The classical approach uses only two: the reference ellipsoid and a geometrical model representing the topography. The rationale for all three components is reviewed, drawing attention to the much smaller precision now needed when transforming residual gravity anomalies. It is shown that all ellipsoidal effects needed for geoid computation with millimetric accuracy are automatically included provided that the free air anomaly and geoid are calculated correctly from the global model. Both must be consistent with an ellipsoidal Earth and with the treatment of observed gravity data. Further ellipsoidal corrections are then negligible. Precise formulae are developed for the geoid height and the free air anomaly using a global gravity model, given as spherical harmonic coefficients. Although only linear in the anomalous potential, these formulae are otherwise exact for an ellipsoidal reference Earth—they involve closed analytical functions of the eccentricity (and the Earths spin rate), rather than a truncated power series in e2. They are evaluated using EGM96 and give ellipsoidal corrections to the conventional free air anomaly ranging from –0.84 to +1.14 mGal, both extremes occurring in Tibet. The geoid error corresponding to these differences is dominated by longer wavelengths, so extrema occur elsewhere, rising to +766 mm south of India and falling to –594 mm over New Guinea. At short wavelengths, the difference between ellipsoidal corrections based only on EGM96 and those derived from detailed local gravity data for the North Sea geoid GEONZ97 has a standard deviation of only 3.3 mm. However, the long-wavelength components missed by the local computation reach 300 mm and have a significant slope. In Australia, for example, such a slope would amount to a 600-mm rise from Perth to Cairns.  相似文献   

10.
LSC法（最小二乘配置法）因能融合不同种类重力观测数据确定大地水准面的特性而受到广泛关注,但由于协方差矩阵存在病态性,微小的观测误差将被协方差矩阵的小奇异值放大,导致计算的配置结果不稳定且精度偏低。本文提出Tikhonov_LSC法,即在LSC法中引入Tikhonov正则化算法,基于GCV法选择协方差矩阵的正则化参数,利用正则化参数修正协方差矩阵的小奇异值,以抑制其对观测误差的放大影响。基于Tikhonov_LSC法计算大地水准面,能有效提高其稳定性和精度。通过以EGM2008重力场模型分别计算山区、丘陵和海域重力异常作为基础数据确定相应区域大地水准面的实验,验证了该方法的有效性。  相似文献   

11.
1　ＩｎｔｒｏｄｕｃｔｉｏｎＤｉｆｆｅｒｅｎｔｇｅｏｉｄｓｏｌｕｔｉｏｎｓｗｅｒｅｃａｒｒｉｅｄｏｕｔｆｏｒＥ ｇｙｐｔｕｓｉｎｇｈｅｔｅｒｏｇｅｎｅｏｕｓｄａｔａａｎｄｄｉｆｆｅｒｅｎｔｍｅｔｈｏｄｏｌｏｇｉｅｓ (Ｅｌ＿Ｔｏｋｈｅｙ ,1 993) .ＴｈｅｍａｉｎｇｏａｌｏｆｔｈｉｓｐａｐｅｒｉｓｔｏｄｅｔｅｒｍｉｎｅａｍｏｓｔａｃｃｕｒａｔｅｎｅｗｇｅｏｉｄｆｏｒＥｇｙｐｔｔａｋｉｎｇａｄｖａｎｔａｇｅｏｆａｎｅｗｕｐｄａｔｅｄｇｒａｖｉｔｙｄａｔａｂａｓｅ,ｔｈｅｉｎｆｏｒｍａｔｉｏｎｇｉｖｅｎｂｙ…  相似文献   

12.
Precise terrain elevation information is required in various remote sensing and Engineering projects. There are many technologies to derive the terrain elevation information like GPS, ground surveys, LiDAR, Photogrammetry. GPS is the most widely used technology to obtain information due to its ease of operation. However the usage of ellipsoidal heights, i.e. with respect to WGS84 has limited usage in hydrological applications. GPS heights must be converted into orthometric heights for use in hydrological applications, and this requires geoid undulation information. These geoid undulations can be deduced from earth gravity models. There are various earth gravity models available for ready usage like EGM96, EGM2008, GFZ96 in the public domain. This paper discusses the improvements observed in deriving orthometric heights using EGM2008 over its predecessor model EGM96. The utilization of the new model in topographical mapping projects are also presented.  相似文献   

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

14.
A global geopotential model, like EGM2008, is not capable of representing the high-frequency components of Earth’s gravity field. This is known as the omission error. In mountainous terrain, omission errors in EGM2008, even when expanded to degree 2,190, may reach amplitudes of 10 cm and more for height anomalies. The present paper proposes the utilisation of high-resolution residual terrain model (RTM) data for computing estimates of the omission error in rugged terrain. RTM elevations may be constructed as the difference between the SRTM (Shuttle Radar Topography Mission) elevation model and the DTM2006.0 spherical harmonic topographic expansion. Numerical tests, carried out in the German Alps with a precise gravimetric quasigeoid model (GCG05) and GPS/levelling data as references, demonstrate that RTM-based omission error estimates improve EGM2008 height anomaly differences by 10 cm in many cases. The comparisons of EGM2008-only height anomalies and the GCG05 model showed 3.7 cm standard deviation after a bias-fit. Applying RTM omission error estimates to EGM2008 reduces the standard deviation to 1.9 cm which equates to a significant improvement rate of 47%. Using GPS/levelling data strongly corroborates these findings with an improvement rate of 49%. The proposed RTM approach may be of practical value to improve quasigeoid determination in mountainous areas without sufficient regional gravity data coverage, e.g., in parts of Asia, South America or Africa. As a further application, RTM omission error estimates will allow refined validation of global gravity field models like EGM2008 from GPS/levelling data.  相似文献   

15.
由CHAMP星载GPS相位双差数据解算地球引力场模型   总被引：1，自引：0，他引：1  
利用7d的CHAMP星载GPS相位观测数据和48个IGS跟踪站的观测数据，构造星地双差相位观测量，进行GPS数据预处理；利用Cowell Ⅱ数值法进行轨道积分和分块Bayes最小二乘参数估计，解算了地球引力场位系数。该模型与EGM96相比（70阶次），大地水准面起伏差异最大为2．872m，差弄精度为0．522m，平均差异为-0．003m，这说明本文解算的地球重力场模型与EGM96没有系统性差异。  相似文献   

16.
The height datum problem and the role of satellite gravity models   总被引：1，自引：0，他引：1  
Regional height systems do not refer to a common equipotential surface, such as the geoid. They are usually referred to the mean sea level at a reference tide gauge. As mean sea level varies (by ±1 to 2 m) from place to place and from continent to continent each tide gauge has an unknown bias with respect to a common reference surface, whose determination is what the height datum problem is concerned with. This paper deals with this problem, in connection to the availability of satellite gravity missions data. Since biased heights enter into the computation of terrestrial gravity anomalies, which in turn are used for geoid determination, the biases enter as secondary or indirect effect also in such a geoid model. In contrast to terrestrial gravity anomalies, gravity and geoid models derived from satellite gravity missions, and in particular GRACE and GOCE, do not suffer from those inconsistencies. Those models can be regarded as unbiased. After a review of the mathematical formulation of the problem, the paper examines two alternative approaches to its solution. The first one compares the gravity potential coefficients in the range of degrees from 100 to 200 of an unbiased gravity field from GOCE with those of the combined model EGM2008, that in this range is affected by the height biases. This first proposal yields a solution too inaccurate to be useful. The second approach compares height anomalies derived from GNSS ellipsoidal heights and biased normal heights, with anomalies derived from an anomalous potential which combines a satellite-only model up to degree 200 and a high-resolution global model above 200. The point is to show that in this last combination the indirect effects of the height biases are negligible. To this aim, an error budget analysis is performed. The biases of the high frequency part are proved to be irrelevant, so that an accuracy of 5 cm per individual GNSS station is found. This seems to be a promising practical method to solve the problem.  相似文献   

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

18.
The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere—at least to the ellipsoid. The goal of this paper is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (theg 1 term). The terrain correction has also been computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. Theg 1 term and the terrain correction were expanded into the spherical harmonics up to180 th order. The corrections (theg 1 term and the terrain correction) have the order of about 2% of theRMS value of degree variance of the disturbing potential per degree. The influences of theg 1 term and the terrain correction on the geoid take the order of 1 meter (RMS value of corrections of the geoid undulation) and on the deflections of the vertical is of the order 0.1″ (RMS value of correction of the deflections of the vertical).  相似文献   

19.
Regional geopotential model improvement for the Iranian geoid determination   总被引：1，自引：0，他引：1  
Spherical harmonic expansions of the geopotential are frequently used for modelling the earth’s gravity field. Degree and order of recently available models go up to 360, corresponding to a resolution of about50 km. Thus, the high degree potential coefficients can be verified nowadays even by locally distributed sets of terrestrial gravity anomalies. These verifications are important when combining the short wavelength model impact, e.g. for regional geoid determinations by means of collocation solutions. A method based on integral formulae is presented, enabling the improvement of geopotential models with respect to non-global distributed gravity anomalies. To illustrate the foregoing, geoid computations are carried out for the area of Iran, introducing theGPM2 geopotential model in combination with available regional gravity data. The accuracy of the geoid determination is estimated from a comparison with Doppler and levelling data to ±1.4m.  相似文献   

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

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

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