首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 906 毫秒
1.
Error analyses of CHAMP data for recovery of the Earth’s gravity field   总被引：1，自引：0，他引：1  
A preliminary commission error analysis whereby orbit perturbation theory and other techniques are used to assess and predict the recovery of the Earths gravity field from the challenging microsatellite payload (CHAMP) mission is developed and implemented. With CHAMP launched in July 2000, accumulated evidence is now available to quantify the errors in the recovery procedure including the orbital precision from GPS, attitude errors, accelerometer noise and thruster mismatch/misalignment. For the latter, numerical integrations using a variable length single-step Runge–Kutta integrator and a fixed length multi-step method are compared to assess the error associated with assuming that the thruster misalignment can be spread uniformly across a step interval. Error degree variances from simulated studies are compared to results from a recently released CHAMP-based gravity field, EIGEN-1S. It is seen that the orbital positioning, as derived from the onboard GPS receiver, is critical, with accelerometer noise contributing at a lower level. Attitude error, at currently quoted accuracy, is not significant as an error source. AcknowledgementsThe authors would like to thank the UK Natural Environment Research Council (Grant No. NER/A/0000/00612) for financing this study and GFZ for supplying the data and technical support.  相似文献   

2.
Theory of integer equivariant estimation with application to GNSS   总被引：4，自引：4，他引：0  
Carrier phase ambiguity resolution is the key to high-precision global navigation satellite system (GNSS) positioning and navigation. It applies to a great variety of current and future models of GPS, modernized GPS and Galileo. The so-called fixed baseline estimator is known to be superior to its float counterpart in the sense that its probability of being close to the unknown but true baseline is larger than that of the float baseline, provided that the ambiguity success rate is sufficiently close to its maximum value of one. Although this is a strong result, the necessary condition on the success rate does not make it hold for all measurement scenarios. It is discussed whether or not it is possible to take advantage of the integer nature of the ambiguities so as to come up with a baseline estimator that is always superior to both its float and its fixed counterparts. It is shown that this is indeed possible, be it that the result comes at the price of having to use a weaker performance criterion. The main result of this work is a Gauss–Markov-like theorem which introduces a new minimum variance unbiased estimator that is always superior to the well-known best linear unbiased (BLU) estimator of the Gauss–Markov theorem. This result is made possible by introducing a new class of estimators. This class of integer equivariant estimators obeys the integer remove–restore principle and is shown to be larger than the class of integer estimators as well as larger than the class of linear unbiased estimators. The minimum variance unbiased estimator within this larger class is referred to as the best integer equivariant (BIE) estimator. The theory presented applies to any model of observation equations having both integer and real-valued parameters, as well as for any probability density function the data might have. AcknowledgementsThis contribution was finalized during the authors stay, as a Tan Chin Tuan Professor, at the Nanyang Technological Universitys GPS Centre (GPSC) in Singapore. The hospitality of the GPSCs director Prof Law Choi Look and his colleagues is greatly appreciated.  相似文献   

3.
Time variations in the Earths gravity field at periods longer than 1 year, for degree-two spherical harmonics, C21, S21, and C20, are estimated from accurately measured Earth rotational variations. These are compared with predictions of atmospheric, oceanic, and hydrologic models, and with independent satellite laser ranging (SLR) results. There is remarkably good agreement between Earth rotation and model predictions of C21 and S21 over a 22-year period. After decadal signals are removed, Earth-rotation-derived interannual C20 variations are dominated by a strong oscillation of period about 5.6 years, probably due to uncertainties in wind and ocean current estimates. The model-predicted C20 agrees reasonably well with SLR observations during the 22-year period, with the exception of the recent anomaly since 1997/1998.  相似文献   

4.
Today the combination of Stokes formula and an Earth gravity model (EGM) for geoid determination is a standard procedure. However, the method of modifying Stokes formula varies from author to author, and numerous methods of modification exist. Most methods modify Stokes kernel, but the most widely applied method, the remove compute restore technique, removes the EGM from the gravity anomaly to attain a residual gravity anomaly under Stokes integral, and at least one known method modifies both Stokes kernel and the gravity anomaly. A general model for modifying Stokes formula is presented; it includes most of the well-known techniques of modification as special cases. By assuming that the error spectra of the gravity anomalies and the EGM are known, the optimum model of modification is derived based on the least-squares principle. This solution minimizes the expected mean square error (MSE) of all possible solutions of the general geoid model. A practical formula for estimating the MSE is also presented. The power of the optimum method is demonstrated in two special cases. AcknowledgementsThis paper was partly written whilst the author was a visiting scientist at The University of New South Wales, Sydney, Australia. He is indebted to Professor W. Kearsley and his colleagues, and their hospitality is acknowledged.  相似文献   

5.
This study makes an initial comparison of three GPS-like constellations. Starting with a simplified constellation of 25 GPS satellites as a reference, GPS(25), we determine what kinematic positioning improvements would result from a constellation comprising a Hi component of 16 GPS satellites (at roughly 16.8 earth radii) coupled with a Lo component of 49 GPS satellites (at roughly 2.1 earth radii). We also include a GPS constellation of 49 GPS satellites, GPS(49), which comprises orbits like the GPS(25) constellation. The GPS(49) and the Hi(16)/Lo(49) constellations have semi-major axes selected so that they have exactly the same average number of satellites above 7.5 degrees elevation (averaged over 24 hours). What motivated this study was a need to measure the benefits, to precision differential kinematic positioning methods (i.e., RTK), which result from the higher Doppler shifts (hence speedier integrated Doppler) generated by the Lo component. Quicker initial convergence was anticipated, of course.  相似文献   

6.
Productivity convergence and spatial dependence among Spanish regions   总被引：3，自引：0，他引：3  
This paper estimates the evolution of labor productivity disparities among 48 Spanish regions over 1980–1996 according to the concepts of - and -convergence. The results of -convergence emphasize the importance of including the impact of neighboring locations productivity and a disaggregate analysis at a sectoral level. In order to measure the narrowing of inequalities, we examine -convergence and reveal that convergence occurs in aggregate labor productivity but not in productivities per sector. The reason comes from a transfer of resources from agriculture towards more productive sectors that has been more pronounced in the poor regions than in the rich ones.The author would like to thank Julie Le Gallo, an anonymous referee, and the participants of the 50th North American Meetings of the RSAI and of the 43rd Annual Meeting of the WRSA for their valuable comments. This paper won the first place the 2004 Tiebout Prize competition, which was awarded at the WRSA meeting, Hawaii, USA, February 26–28.  相似文献   

7.
  相似文献   

8.
A simple statistical approach has been applied to the repeated electro-optical distance measurements (EDM) of 1,358 lines in the Tohoku district of Japan to obtain knowledge about the precision of EDM and the possible accumulation of strain. The average time interval between measurements is about seven or eight years. It is shown that the whole data of the difference between distance measurements repeated over a given lineD are interpreted in terms of EDM errors comprising distance proportional systematic errors and standard errors expressed by the usual form . The rate of horizontal deformation must therefore be much smaller than the strain rates of about 0.7 0.8 ppm over 7 to 8 years which have been hitherto expected.  相似文献   

9.
In a modern application of Stokes formula for geoid determination, regional terrestrial gravity is combined with long-wavelength gravity information supplied by an Earth gravity model. Usually, several corrections must be added to gravity to be consistent with Stokes formula. In contrast, here all such corrections are applied directly to the approximate geoid height determined from the surface gravity anomalies. In this way, a more efficient workload is obtained. As an example, in applications of the direct and first and second indirect topographic effects significant long-wavelength contributions must be considered, all of which are time consuming to compute. By adding all three effects to produce a combined geoid effect, these long-wavelength features largely cancel. The computational scheme, including two least squares modifications of Stokes formula, is outlined, and the specific advantages of this technique, compared to traditional gravity reduction prior to Stokes integration, are summarised in the conclusions and final remarks. AcknowledgementsThis paper was written whilst the author was a visiting scientist at Curtin University of Technology, Perth, Australia. The hospitality and fruitful discussions with Professor W. Featherstone and his colleagues are gratefully acknowledged.  相似文献   

10.
On Helmert’s methods of condensation   总被引：2，自引：0，他引：2  
B. Heck 《Journal of Geodesy》2003,77(3-4):155-170
Helmerts first and second method of condensation are reviewed and generalized in two respects: First, the point at which the effects of topographical and condensation masses are calculated may be situated on or outside the topographical surface; second, the depth of the condensation layer below the geoid is arbitrary. While the first extension permits the application of the generalized model to the evaluation of airborne and satellite data, the second one gives an additional degree of freedom which can be used to provide a smooth gravity field after reducing the observation data. The respective formulae are derived for the generalized condensation model in both planar and spherical approximation. A comparison of the planar and the spherical model shows some structural differences, which are primarily visible in the out-of-integral terms. Considering the respective formulae for the combined topographic–condensation reduction on the background of the density structure of the Earths lithosphere, the consequences for the residual gravity field are investigated; it is shown that the residual field after applying Helmerts second model of reduction is very rough, making this procedure unfavourable for downward continuation. Further considerations refer to the question of which sets of formulae should be used in geoid and quasigeoid determination. It is concluded that for high-precision applications the generalized spherical model, involving a depth of the condensation layer of between 20 and 30 km, should be superior to Helmerts second model of condensation, although it requires the direct calculation of the indirect effect, which is larger than in the case of Helmerts second method of condensation.  相似文献   

11.
Integrated adjustment of CHAMP, GRACE, and GPS data   总被引：16，自引：3，他引：13  
Various types of observations, such as space-borne Global positioning system (GPS) code and phase data, accelerometer data, K-band range and range-rate data, and ground-based satellite laser ranging data of the CHAllenging Minisatellite Payload (CHAMP) and GRAvity Climate Experiment (GRACE) satellite missions, are used together with ground-based GPS code and phase data in a rigorous adjustment to eventually solve for the ephemerides of the CHAMP, GRACE, and GPS satellites, geocenter variations, and low-degree gravity field parameters. It turns out that this integrated adjustment considerably improves the accuracy of the ephemerides for the high and low satellites, geocenter variations, and gravity field parameters, compared to the case when the adjustment is carried out stepwise or in individual satellite solutions.Acknowledgments. This study has been supported by the German Ministry of Education and Research through the Geotechnologies Programme grants 03F0333A/CHAMP and 03F0326A/GRACE.  相似文献   

12.
Mean 5 × 5 heights and depths from ETOPO5U (Earth Topography at 5 spacing Updated) Digital Terrain Model (DTM) were compared with corresponding quantities of a local DTM in the test area [38° 40°, 21° 24°]. From this comparison a shift of ETOPO5U with respect to the local DTM in the longitudinal direction equal to 5 min was found after applying an efficient fast Fourier transform (FFT) technique. Furthermore, sparse mean height differences larger than 1,000 m were observed between ETOPO5U and the local DTM due rather to errors of ETOPO5U. The effect of these errors on gravity and height anomalies was computed in a subregion of the area under consideration.  相似文献   

13.
W. Sun 《Journal of Geodesy》2003,77(7-8):381-387
An asymptotic theory is presented for calculating co-seismic potential and geoid changes, as an approximation of the dislocation theory for a spherical Earth. This theory is given by a closed-form mathematical expression, so that it is mathematically simple and can be applied easily. Moreover, since the asymptotic theory includes sphericity and vertical structure effects, it is physically more reasonable than the flat-Earth theory. A comparison between results calculated by three dislocation theories (the flat-Earth theory, the theory for a spherical Earth and its asymptotic solution) shows that the true co-seismic geoid changes are approximated better by the asymptotic results than by those of a flat Earth. Numerical results indicate that the sphericity effect is obvious large, especially for a tensile source on a vertical fault plane. AcknowledgementsThe author is grateful to Dr S. Okubo for his helpful suggestions and discussions. Comments by anonymous reviewers are also greatly acknowledged. This research was financially supported by JSPS research grants (C13640420) and Basic design and feasibility studies for the future missions for monitoring Earths environment.  相似文献   

14.
Harmonic maps     
Harmonic maps are generated as a certain class of optimal map projections. For instance, if the distortion energy over a meridian strip of the International Reference Ellipsoid is minimized, we are led to the Laplace–Beltrami vector-valued partial differential equation. Harmonic functions x(L,B), y(L,B) given as functions of ellipsoidal surface parameters of Gauss ellipsoidal longitude L and Gauss ellipsoidal latitude B, as well as x(,q), y(,q) given as functions of relative isometric longitude =L–L0 and relative isometric latitude q=Q–Q0 gauged to a vector-valued boundary condition of special symmetry are constructed. The easting and northing {x(b,),y(b,)} of the new harmonic map is then given. Distortion energy analysis of the new harmonic map is presented, as well as case studies for (1) B[–40°,+40°], L[–31°,+49°], B0= ±30°, L0=9° and (2) B[46°,56°], L{[4.5°, 7.5°]; [7.5°, 10.5°]; [10.5°,13.5°]; [13.5°,16.5°]}, B0= 51°, L0 {6°,9°,12°,15°}.  相似文献   

15.
Marine gravity surveying line system adjustment   总被引：6，自引：0，他引：6  
The general theories and methods of marine surveying line system adjustment were introduced in (1979) and Tang (1991) . According to the characteristics of marine gravity measurement, this paper presents a new method of combined adjustment which takes into account both direct and indirect influence of position errors. The method is particularly suitable to be used in the post- processing of marine gravity observation data. With some practical applications, it is proved to be effective in improving the quality of marine gravity data.  相似文献   

16.
A new theory for high-resolution regional geoid computation without applying Stokess formula is presented. Operationally, it uses various types of gravity functionals, namely data of type gravity potential (gravimetric leveling), vertical derivatives of the gravity potential (modulus of gravity intensity from gravimetric surveys), horizontal derivatives of the gravity potential (vertical deflections from astrogeodetic observations) or higher-order derivatives such as gravity gradients. Its algorithmic version can be described as follows: (1) Remove the effect of a very high degree/order potential reference field at the point of measurement (POM), in particular GPS positioned, either on the Earths surface or in its external space. (2) Remove the centrifugal potential and its higher-order derivatives at the POM. (3) Remove the gravitational field of topographic masses (terrain effect) in a zone of influence of radius r. A proper choice of such a radius of influence is 2r=4×104 km/n, where n is the highest degree of the harmonic expansion. (cf. Nyquist frequency). This third remove step aims at generating a harmonic gravitational field outside a reference ellipsoid, which is an equipotential surface of a reference potential field. (4) The residual gravitational functionals are downward continued to the reference ellipsoid by means of the inverse solution of the ellipsoidal Dirichlet boundary-value problem based upon the ellipsoidal Abel–Poisson kernel. As a discretized integral equation of the first kind, downward continuation is Phillips–Tikhonov regularized by an optimal choice of the regularization factor. (5) Restore the effect of a very high degree/order potential reference field at the corresponding point to the POM on the reference ellipsoid. (6) Restore the centrifugal potential and its higher-order derivatives at the ellipsoidal corresponding point to the POM. (7) Restore the gravitational field of topographic masses ( terrain effect) at the ellipsoidal corresponding point to the POM. (8) Convert the gravitational potential on the reference ellipsoid to geoidal undulations by means of the ellipsoidal Bruns formula. A large-scale application of the new concept of geoid computation is made for the Iran geoid. According to the numerical investigations based on the applied methodology, a new geoid solution for Iran with an accuracy of a few centimeters is achieved.Acknowledgments. The project of high-resolution geoid computation of Iran has been support by National Cartographic Center (NCC) of Iran. The University of Tehran, via grant number 621/3/602, supported the computation of a global geoid solution for Iran. Their support is gratefully acknowledged. A. Ardalan would like to thank Mr. Y. Hatam, and Mr. K. Ghazavi from NCC and Mr. M. Sharifi, Mr. A. Safari, and Mr. M. Motagh from the University of Tehran for their support in data gathering and computations. The authors would like to thank the comments and corrections made by the four reviewers and the editor of the paper, Professor Will Featherstone. Their comments helped us to correct the mistakes and improve the paper.  相似文献   

17.
Summary The least-squares collocation method has been used for the computation of a geoid solution in central Spain, combining a geopotential model complete to degree and order 360, gravity anomalies and topographic information. The area has been divided in two 1°× 1° blocks and predictions have been done in each block with gravity data spacing about 5 × 5 within each block, extended 1/2°. Topographic effects have been calculated from 6 × 9 heights using an RTM reduction with a reference terrain model of 30 × 30 mean heights.  相似文献   

18.
The three-dimensional (3-D) resection problem is usually solved by first obtaining the distances connecting the unknown point P{X,Y,Z} to the known points Pi{Xi,Yi,Zi}i=1,2,3 through the solution of the three nonlinear Grunert equations and then using the obtained distances to determine the position {X,Y,Z} and the 3-D orientation parameters {,, }. Starting from the work of the German J. A. Grunert (1841), the Grunert equations have been solved in several substitutional steps and the desire as evidenced by several publications has been to reduce these number of steps. Similarly, the 3-D ranging step for position determination which follows the distance determination step involves the solution of three nonlinear ranging (`Bogenschnitt') equations solved in several substitution steps. It is illustrated how the algebraic technique of Groebner basis solves explicitly the nonlinear Grunert distance equations and the nonlinear 3-D ranging (`Bogenschnitt') equations in a single step once the equations have been converted into algebraic (polynomial) form. In particular, the algebraic tool of the Groebner basis provides symbolic solutions to the problem of 3-D resection. The various forward and backward substitution steps inherent in the classical closed-form solutions of the problem are avoided. Similar to the Gauss elimination technique in linear systems of equations, the Groebner basis eliminates several variables in a multivariate system of nonlinear equations in such a manner that the end product normally consists of a univariate polynomial whose roots can be determined by existing programs e.g. by using the roots command in Matlab.Acknowledgments.The first author wishes to acknowledge the support of JSPS (Japan Society of Promotion of Science) for the financial support that enabled the completion of the write-up of the paper at Kyoto University, Japan. The author is further grateful for the warm welcome and the good working atmosphere provided by his hosts Professors S. Takemoto and Y. Fukuda of the Department of Geophysics, Graduate School of Science, Kyoto University, Japan.  相似文献   

19.
Geoid determination in Turkey (TG-91)   总被引：1，自引：0，他引：1  
It is considered that precise geoid determination is one of the main current geodetic problems in Turkey since GPS defined coordinates require geoidal heights in practice. In order to determine the geoid by least squares collocation (LSC) the area covering Turkey was divided into 114 blocks of size 1° × 1°. LSC approximation to the geoid based upon the tailored geopotential model GPM2-T1 is constructed within each block. The model GPM2-T1 complete to degree and order 200 has been developed by tailoring of the model GPM2 to mean free-air anomalies and mean heights of one degree blocks in Turkey. Terrain effect reduced point gravity data spaced 5 × 5 within each block which the sides extended 0°.5 were used in LSC. Residual terrain model (RTM) depends on point heights at 15×20 griding and 5×5 and 15×15 mean heights has been carried out in terrain effect reduction. Indirect effect of RTM on geoid is also taken into account. The geoid, called Turkish Geoid 1991 (TG-91), referenced to GRS-80 ellipsoid has been computed at 3 × 3 griding nodes within each block. The quality of the TG-91 is also evaluated by comparing computed and GPS derived geoidal height differences, and 2.1 – 2.6 ppm accuracy for average baseline lenght of 45 km is obtained.  相似文献   

20.
Summary Riemann polar/normal coordinates are the constituents to generate the oblique azimuthal projection of geodesic type, here applied to the reference ellipsoid of revolution (biaxial ellipsoid).Firstly we constitute a minimal atlas of the biaxial ellipsoid built on {ellipsoidal longitude, ellipsoidal latitude} and {metalongitude, metalatitude}. TheDarboux equations of a 1-dimensional submanifold (curve) in a 2-dimensional manifold (biaxial ellipsoid) are reviewed, in particular to represent geodetic curvature, geodetic torsion and normal curvature in terms of elements of the first and second fundamental form as well as theChristoffel symbols. The notion of ageodesic anda geodesic circle is given and illustrated by two examples. The system of twosecond order ordinary differential equations of ageodesic (Lagrange portrait) is presented in contrast to the system of twothird order ordinary differential equations of ageodesic circle (Proofs are collected inAppendix A andB). A precise definition of theRiemann mapping/mapping of geodesics into the local tangent space/tangent plane has been found.Secondly we computeRiemann polar/normal coordinates for the biaxial ellipsoid, both in theLagrange portrait (Legendre series) and in theHamilton portrait (Lie series).Thirdly we have succeeded in a detailed deformation analysis/Tissot distortion analysis of theRiemann mapping. The eigenvalues — the eigenvectors of the Cauchy-Green deformation tensor by means of ageneral eigenvalue-eigenvector problem have been computed inTable 3.1 andTable 3.2 (1, 2 = 1) illustrated inFigures 3.1, 3.2 and3.3. Table 3.3 contains the representation ofmaximum angular distortion of theRiemann mapping. Fourthly an elaborate global distortion analysis with respect toconformal Gau-Krüger, parallel Soldner andgeodesic Riemann coordinates based upon theAiry total deformation (energy) measure is presented in a corollary and numerically tested inTable 4.1. In a local strip [-l E,l E] = [-2°, +2°], [b S,b N] = [-2°, +2°]Riemann normal coordinates generate the smallest distortion, next are theparallel Soldner coordinates; the largest distortion by far is met by theconformal Gau-Krüger coordinates. Thus it can be concluded that for mapping of local areas of the biaxial ellipsoid surface the oblique azimuthal projection of geodesic type/Riemann polar/normal coordinates has to be favored with respect to others.  相似文献   

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

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