共查询到20条相似文献,搜索用时 31 毫秒
1.
The standard analytical approach which is applied for constructing geopotential models OSU86 and earlier ones, is based on
reducing the boundary value equation to a sphere enveloping the Earth and then solving it directly with respect to the potential
coefficients
n,m
. In an alternative procedure, developed by Jekeli and used for constructing the models OSU91 and EGM96, at first an ellipsoidal
harmonic series is developed for the geopotential and then its coefficients
n,m
e
are transformed to the unknown
n,m
. The second solution is more exact, but much more complicated. The standard procedure is modified and a new simple integral
formula is derived for evaluating the potential coefficients. The efficiency of the standard and new procedures is studied
numerically. In these solutions the same input data are used as for constructing high-degree parts of the EGM96 models. From
two sets of
n,m
(n≤360,|m|≤n), derived by the standard and new approaches, different spectral characteristics of the gravity anomaly and the geoid undulation
are estimated and then compared with similar characteristics evaluated by Jekeli's approach (`etalon' solution). The new solution
appears to be very close to Jekeli's, as opposed to the standard solution. The discrepancies between all the characteristics
of the new and `etalon' solutions are smaller than the corresponding discrepancies between two versions of the final geopotential
model EGM96, one of them (HDM190) constructed by the block-diagonal least squares (LS) adjustment and the other one (V068)
by using Jekeli's approach. On the basis of the derived analytical solution a new simple mathematical model is developed to
apply the LS technique for evaluating geopotential coefficients.
Received: 12 December 2000 / Accepted: 21 June 2001 相似文献
2.
A methodology for precise determination of the fundamental geodetic parameter w
0, the potential value of the Gauss–Listing geoid, as well as its time derivative 0, is presented. The method is based on: (1) ellipsoidal harmonic expansion of the external gravitational field of the Earth
to degree/order 360/360 (130 321 coefficients; http://www.uni-stuttgard.de/gi/research/ index.html projects) with respect
to the International Reference Ellipsoid WGD2000, at the GPS positioned stations; and (2) ellipsoidal free-air gravity reduction
of degree/order 360/360, based on orthometric heights of the GPS-positioned stations. The method has been numerically tested
for the data of three GPS campaigns of the Baltic Sea Level project (epochs 1990.8,1993.4 and 1997.4). New w
0 and 0 values (w
0=62 636 855.75 ± 0.21 m2/s2, 0=−0.0099±0.00079 m2/s2 per year, w
0/&γmacr;=6 379 781.502 m,0/&γmacr;=1.0 mm/year, and &γmacr;= −9.81802523 m2/s2) for the test region (Baltic Sea) were obtained. As by-products of the main study, the following were also determined: (1)
the high-resolution sea surface topography map for the Baltic Sea; (2) the most accurate regional geoid amongst four different
regional Gauss–Listing geoids currently proposed for the Baltic Sea; and (3) the difference between the national height datums
of countries around the Baltic Sea.
Received: 14 August 2000 / Accepted: 19 June 2001 相似文献
3.
The Somigliana–Pizzetti gravity field (the International gravity formula), namely the gravity field of the level ellipsoid
(the International Reference Ellipsoid), is derived to the sub-nanoGal accuracy level in order to fulfil the demands of modern
gravimetry (absolute gravimeters, super conducting gravimeters, atomic gravimeters). Equations (53), (54) and (59) summarise
Somigliana–Pizzetti gravity Γ(φ,u) as a function of Jacobi spheroidal latitude φ and height u to the order ?(10−10 Gal), and Γ(B,H) as a function of Gauss (surface normal) ellipsoidal latitude B and height H to the order ?(10−10 Gal) as determined by GPS (`global problem solver'). Within the test area of the state of Baden-Württemberg, Somigliana–Pizzetti
gravity disturbances of an average of 25.452 mGal were produced. Computer programs for an operational application of the new
international gravity formula with (L,B,H) or (λ,φ,u) coordinate inputs to a sub-nanoGal level of accuracy are available on the Internet.
Received: 23 June 2000 / Accepted: 2 January 2001 相似文献
4.
The New Hebrides experiment consisted of setting up a pair of DORIS beacons in remote tropical islands in the southwestern
Pacific, between 1993 and 1997. Because of orbitography requirements on TOPEX/Poséidon, the beacons were only transmitting
to SPOT satellites. Root-mean-square (RMS) scatters at the centimeter level on the latitude and vertical components were achieved,
but 2-cm RMS scatters affected the longitude component. Nevertheless, results of relative velocity (123 mm/year N250°) are
very consistent with those obtained using the global positioning system (GPS) (126 mm/yr N246°). The co-seismic step (12 mm
N60°) related to the Walpole event (M
W = 7.7) is consistent with that derived from GPS (10 mm N30°) or from the centroid moment tensor (CMT) of the quake (12 mm
N000°).
Received: 19 November 1999 / Accepted: 17 May 2000 相似文献
5.
Since the beginning of the International Global Navigation Satellite System (GLONASS) Experiment, IGEX, in October 1998,
the Center for Orbit Determination in Europe (CODE) has acted as an analysis center providing precise GLONASS orbits on a
regular basis. In CODE's IGEX routine analysis the Global Positioning System (GPS) orbits and Earth rotation parameters are
introduced as known quantities into the GLONASS processing. A new approach is studied, where data from the IGEX network are
combined with GPS observations from the International GPS Service (IGS) network and all parameters (GPS and GLONASS orbits,
Earth rotation parameters, and site coordinates) are estimated in one processing step. The influence of different solar radiation
pressure parameterizations on the GLONASS orbits is studied using different parameter subsets of the extended CODE orbit model.
Parameterization with three constant terms in the three orthogonal directions, D, Y, and X (D = direction satellite–Sun, Y = direction of the satellite's solar panel axis), and two periodic terms in the X-direction, proves to be adequate for GLONASS satellites. As a result of the processing it is found that the solar radiation
pressure effect for the GLONASS satellites is significantly different in the Y-direction from that for the GPS satellites, and an extensive analysis is carried out to investigate the effect in detail.
SLR observations from the ILRS network are used as an independent check on the quality of the GLONASS orbital solutions. Both
processing aspects, combining the two networks and changing the orbit parameterization, significantly improve the quality
of the determined GLONASS orbits compared to the orbits stemming from CODE's IGEX routine processing.
Received: 10 May 2000 / Accepted: 9 October 2000 相似文献
6.
A potential-type Molodensky telluroid based upon a minimum-distance mapping is derived. With respect to a reference potential
of Somigliana–Pizzetti type which relates to the World Geodetic Datum 2000, it is shown that a point-wise minimum-distance
mapping of the topographical surface of the Earth onto the telluroid surface, constrained to the gauge W(P)=u(p), leads to a system of four nonlinear normal equations. These normal equations are solved by a fast Newton–Raphson iteration.
Received: 7 February 2000 / Accepted: 23 October 2001 相似文献
7.
A 2×2 arc-minute resolution geoid model, CARIB97, has been computed covering the Caribbean Sea. The geoid undulations refer
to the GRS-80 ellipsoid, centered at the ITRF94 (1996.0) origin. The geoid level is defined by adopting the gravity potential
on the geoid as W
0=62 636 856.88 m2/s2 and a gravity-mass constant of GM=3.986 004 418×1014 m3/s2. The geoid model was computed by applying high-frequency corrections to the Earth Gravity Model 1996 global geopotential
model in a remove-compute-restore procedure. The permanent tide system of CARIB97 is non-tidal. Comparison of CARIB97 geoid
heights to 31 GPS/tidal (ITRF94/local) benchmarks shows an average offset (h–H–N) of 51 cm, with an Root Mean Square (RMS) of 62 cm about the average. This represents an improvement over the use of a global
geoid model for the region. However, because the measured orthometric heights (H) refer to many differing tidal datums, these comparisons are biased by localized permanent ocean dynamic topography (PODT).
Therefore, we interpret the 51 cm as partially an estimate of the average PODT in the vicinity of the 31 island benchmarks.
On an island-by-island basis, CARIB97 now offers the ability to analyze local datum problems which were previously unrecognized
due to a lack of high-resolution geoid information in the area.
Received: 2 January 1998 / Accepted: 18 August 1998 相似文献
8.
Considering a GPS satellite and two terrestrial stations, two types of equations are derived relating the heights of the
two stations to the measured data (frequency ratio or clock rate differences) and the coordinates and velocity components
of all three participating objects. The potential possibilities of using such relations for the determination of heights (in
terms of geopotential numbers or orthometric heights) are discussed.
Received: 6 December 2000 / Accepted: 9 July 2001 相似文献
9.
GPS-assisted GLONASS orbit determination 总被引:1,自引:0,他引:1
D. Kuang Y. E. Bar-Sever W. I. Bertiger K. J. Hurst J. F. Zumberge 《Journal of Geodesy》2001,75(11):569-574
Using 1 week of data from a network of GPS/GLONASS dual-tracking receivers, 15-cm accurate GLONASS orbit determination is
demonstrated with an approach that combines GPS and GLONASS data. GPS data are used to define the reference frame, synchronize
receiver clocks and determine troposphere delay for the GLONASS tracking network. GLONASS tracking data are then processed
separately, with the GPS-defined parameters held fixed, to determine the GLONASS orbit. The quality of the GLONASS orbit determination
is currently limited by the size and distribution of the tracking network, and by the unavailability of a sufficiently refined
solar pressure model. Temporal variations in the differential clock bias of the dual-tracking receivers are found to have
secondary impact on the orbit determination accuracy.
Received: 5 January 2000 / Accepted: 15 February 2001 相似文献
10.
Iterative vector methods for computing geodetic latitude and height from rectangular coordinates 总被引:4,自引:4,他引:4
J. Pollard 《Journal of Geodesy》2002,76(1):36-40
Two iterative vector methods for computing geodetic coordinates (φ, h) from rectangular coordinates (x, y, z) are presented. The methods are conceptually simple, work without modification at any latitude and are easy to program. Geodetic
latitude and height can be calculated to acceptable precision in one iteration over the height range from −106 to +109 m.
Received: 13 December 2000 / Accepted: 13 July 2001 相似文献
11.
Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt 总被引:1,自引:0,他引:1
C. C. Tscherning Awar Radwan A. A. Tealeb S. M. Mahmoud M. Abd El-Monum Ramdan Hassan I. El-Syaed K. Saker 《Journal of Geodesy》2001,75(7-8):343-348
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 相似文献
12.
Stochastic assessment of GPS carrier phase measurements for precise static relative positioning 总被引:17,自引:11,他引:17
Global positioning system (GPS) carrier phase measurements are used in all precise static relative positioning applications.
The GPS carrier phase measurements are generally processed using the least-squares method, for which both functional and stochastic
models need to be carefully defined. Whilst the functional model for precise GPS positioning is well documented in the literature,
realistic stochastic modelling for the GPS carrier phase measurements is still both a controversial topic and a difficult
task to accomplish in practice. The common practice of assuming that the raw GPS measurements are statistically independent
in space and time, and have the same accuracy, is certainly not realistic. Any mis-specification in the stochastic model will
inevitably lead to unreliable positioning results. A stochastic assessment procedure has been developed to take into account
the heteroscedastic, space- and time-correlated error structure of the GPS measurements. Test results indicate that the reliability
of the estimated positioning results is improved by applying the developed stochastic assessment procedure. In addition, the
quality of ambiguity resolution can be more realistically evaluated.
Received: 13 February 2001 / Accepted: 3 September 2001 相似文献
13.
The parameter distributions of the integer GPS model 总被引:6,自引:0,他引:6
P. J. G. Teunissen 《Journal of Geodesy》2002,76(1):41-48
A parameter estimation theory is incomplete if no rigorous measures are available for describing the uncertainty of the parameter
estimators. Since the classical theory of linear estimation does not apply to the integer GPS model, rigorous probabilistic
statements cannot be made with reference to the classical results. The fact that integer parameters are involved in the estimation
process forces a reappraisal of the propagation of uncertainty. It is with this purpose in mind that the joint and marginal
distributional properties of both the integer and non-integer parameters of the GPS model are determined. These joint distributions
can also be used to determine the distribution of functions of the parameters. As an important example, the distribution of
the vector of ambiguity residuals is determined.
Received: 30 January 2001 / Accepted: 31 July 2001 相似文献
14.
Ionospheric variation may be considered as a stationary time series under quiet conditions. However, the disturbance of a
stationary random process from stationarity results in the bias of corresponding samples from the stationary observations,
and in the change of statistical model parameters of the process. From a general mathematical aspect, a new method is presented
for monitoring ionospheric variations, based on the characteristic of time-series observation of GPS, and an investigation
of the statistical properties of the estimated auto-covariance of the random ionospheric delay when changing the number of
samples in the time series is carried out. A preliminary scheme for monitoring ionospheric delays is proposed.
Received: 18 August 2000 / Accepted: 12 April 2001 相似文献
15.
GPS measurements of ocean loading and its impact on zenith tropospheric delay estimates: a case study in Brittany, France 总被引:1,自引:0,他引:1
S. Vey E. Calais M. Llubes N. Florsch G. Woppelmann J. Hinderer M. Amalvict M. F. Lalancette B. Simon F. Duquenne J. S. Haase 《Journal of Geodesy》2002,76(8):419-427
The results from a global positioning system (GPS) experiment carried out in Brittany, France, in October 1999, aimed at
measuring crustal displacements caused by ocean loading and quantifying their effects on GPS-derived tropospheric delay estimates,
are presented. The loading effect in the vertical and horizontal position time series is identified, however with significant
disagreement in amplitude compared to ocean loading model predictions. It is shown that these amplitude misfits result from
spatial tropospheric heterogeneities not accounted for in the data processing. The effect of ocean loading on GPS-derived
zenith total delay (ZTD) estimates is investigated and a scaling factor of 4.4 between ZTD and station height for a 10° elevation
cut-off angle is found (i.e. a 4.4-cm station height error would map into a 1-cm ZTD error). Consequently, unmodeled ocean
loading effects map into significant errors in ZTD estimates and ocean loading modeling must be properly implemented when
estimating ZTD parameters from GPS data for meteorological applications. Ocean loading effects must be known with an accuracy
of better than 3 cm in order to meet the accuracy requirements of meteorological and climatological applications of GPS-derived
precipitable water vapor.
Received: 16 July 2001 / Accepted: 25 April 2002
Acknowledgments. The authors are grateful to H.G. Scherneck for fruitful discussions and for his help with the ocean loading calculations.
They thank H. Vedel for making the HIRLAM data available; D. Jerett for helpful discussions; and the city of Rostrenen, the
Laboratoire d'Océanographie of Concarneau, and the Institut de Protection et de S?reté Nucléaire (BERSSIN) for their support
during the GPS measurement campaign. Reviews by C.K. Shum and two anonymous referees significantly improved this paper. This
work was carried out in the framework of the MAGIC project (http://www.acri.fr/magic), funded by the European Commission,
Environment and Climate Program (EC Contract ENV4-CT98–0745).
Correspondence to: E. Calais, Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, IN 47907-1397, USA. e-mail:
ecalais@purdue.edu Tel. : +1-765-496-2915; Fax:+1-765-496-1210 相似文献
16.
This paper generalizes the Stokes formula from the spherical boundary surface to the ellipsoidal boundary surface. The resulting
solution (ellipsoidal geoidal height), consisting of two parts, i.e. the spherical geoidal height N
0 evaluated from Stokes's formula and the ellipsoidal correction N
1, makes the relative geoidal height error decrease from O(e
2) to O(e
4), which can be neglected for most practical purposes. The ellipsoidal correction N
1 is expressed as a sum of an integral about the spherical geoidal height N
0 and a simple analytical function of N
0 and the first three geopotential coefficients. The kernel function in the integral has the same degree of singularity at
the origin as the original Stokes function. A brief comparison among this and other solutions shows that this solution is
more effective than the solutions of Molodensky et al. and Moritz and, when the evaluation of the ellipsoidal correction N
1 is done in an area where the spherical geoidal height N
0 has already been evaluated, it is also more effective than the solution of Martinec and Grafarend.
Received: 27 January 1999 / Accepted: 4 October 1999 相似文献
17.
When standard boundary element methods (BEM) are used in order to solve the linearized vector Molodensky problem we are confronted with
two problems: (1) the absence of O(|x|−2) terms in the decay condition is not taken into account, since the single-layer ansatz, which is commonly used as representation
of the disturbing potential, is of the order O(|x|−1) as x→∞. This implies that the standard theory of Galerkin BEM is not applicable since the injectivity of the integral operator
fails; (2) the N×N stiffness matrix is dense, with N typically of the order 105. Without fast algorithms, which provide suitable approximations to the stiffness matrix by a sparse one with O(N(logN)
s
), s≥0, non-zero elements, high-resolution global gravity field recovery is not feasible. Solutions to both problems are proposed.
(1) A proper variational formulation taking the decay condition into account is based on some closed subspace of co-dimension
3 of the space of square integrable functions on the boundary surface. Instead of imposing the constraints directly on the
boundary element trial space, they are incorporated into a variational formulation by penalization with a Lagrange multiplier.
The conforming discretization yields an augmented linear system of equations of dimension N+3×N+3. The penalty term guarantees the well-posedness of the problem, and gives precise information about the incompatibility
of the data. (2) Since the upper left submatrix of dimension N×N of the augmented system is the stiffness matrix of the standard BEM, the approach allows all techniques to be used to generate
sparse approximations to the stiffness matrix, such as wavelets, fast multipole methods, panel clustering etc., without any
modification. A combination of panel clustering and fast multipole method is used in order to solve the augmented linear system
of equations in O(N) operations. The method is based on an approximation of the kernel function of the integral operator by a degenerate kernel
in the far field, which is provided by a multipole expansion of the kernel function. Numerical experiments show that the fast
algorithm is superior to the standard BEM algorithm in terms of CPU time by about three orders of magnitude for N=65 538 unknowns. Similar holds for the storage requirements. About 30 iterations are necessary in order to solve the linear
system of equations using the generalized minimum residual method (GMRES). The number of iterations is almost independent
of the number of unknowns, which indicates good conditioning of the system matrix.
Received: 16 October 1999 / Accepted: 28 February 2001 相似文献
18.
Until recently, the Global Positioning System (GPS) was the only operational means of distributing time to an arbitrary number
of users and of synchronizing clocks over large distances with a high degree of precision and accuracy. Over the last few
years it has been shown that similar performance can be achieved using the Russian Global Navigation Satellite System (GLONASS).
GLONASS time transfer between continents was initially hampered by the lack of post-processed precise ephemerides. Results
from the International GLONASS Experiment (IGEX) campaign are now available, however, and this paper reports on the first
use of IGEX precise ephemerides for GLONASS P-code intercontinental time links. The results of GLONASS P-code and GPS C/A-code
time transfer are compared under similar conditions.
Received: 31 January 2000 / Accepted: 10 July 2000 相似文献
19.
The problem of phase ambiguity resolution in global positioning system (GPS) theory is considered. The Bayesian approach
is applied to this problem and, using Monte Carlo simulation to search over the integer candidates, a practical expression
for the Bayesian estimator is obtained. The analysis of the integer grid points inside the search ellipsoid and their evolution
with time, while measurements are accumulated, leads to the development of a Bayesian theory based on a mathematical mixture
model for the ambiguity.
Received: 29 March 2001 / Accepted: 3 September 2001 相似文献
20.
C. Vigny J. Chéry T. Duquesnoy F. Jouanne J. Ammann M. Anzidei J.-P. Avouac F. Barlier R. Bayer P. Briole E. Calais F. Cotton F. Duquenne K. L. Feigl G. Ferhat M. Flouzat J.-F. Gamond A. Geiger A. Harmel M. Kasser M. Laplanche M. Le Pape J. Martinod G. Ménard B. Meyer J.-C. Ruegg J.-M. Scheubel O. Scotti G. Vidal 《Journal of Geodesy》2002,76(2):63-76
The Western Alps are among the best studied collisional belts with both detailed structural mapping and also crustal geophysical
investigations such as the ECORS and EGT seismic profile. By contrast, the present-day kinematics of the belt is still largely
unknown due to small relative motions and the insufficient accuracy of the triangulation data. As a consequence, several tectonic
problems still remain to be solved, such as the amount of N–S convergence in the Occidental Alps, the repartition of the deformation
between the Alpine tectonic units, and the relation between deformation and rotation across the Alpine arc. In order to address
these problems, the GPS ALPES group, made up of French, Swiss and Italian research organizations, has achieved the first large-scale
GPS surveys of the Western Alps. More than 60 sites were surveyed in 1993 and 1998 with a minimum observation of 3 days at
each site. GPS data processing has been done by three independent teams using different software. The different solutions
have horizontal repeatabilities (N–E) of 4–7 mm in 1993 and 2–3 mm in 1998 and compare at the 3–5-mm level in position and
2-mm/yr level in velocity. A comparison of 1993 and 1998 coordinates shows that residual velocities of the GPS marks are generally
smaller than 2 mm/yr, precluding a detailed tectonic interpretation of the differential motions. However, these data seem
to suggest that the N–S compression of the Western Alps is quite mild (less than 2 mm/yr) compared to the global convergence
between the African and Eurasian plate (6 mm/yr). This implies that the shortening must be accomodated elsewhere by the deformation
of the Maghrebids and/or by rotations of Mediterranean microplates. Also, E–W velocity components analysis supports the idea
that E–W extension exists, as already suggested by recent structural and seismotectonic data interpretation.
Received: 27 November 2000 / Accepted: 17 September 2001 相似文献