共查询到20条相似文献,搜索用时 12 毫秒
1.
Downward continuation and geoid determination based on band-limited airborne gravity data 总被引:4,自引:3,他引:4
The downward continuation of the harmonic disturbing gravity potential, derived at flight level from discrete observations
of airborne gravity by the spherical Hotine integral, to the geoid is discussed. The initial-boundary-value approach, based
on both the direct and inverse solution to Dirichlet's problem of potential theory, is used. Evaluation of the discretized
Fredholm integral equation of the first kind and its inverse is numerically tested using synthetic airborne gravity data.
Characteristics of the synthetic gravity data correspond to typical airborne data used for geoid determination today and in
the foreseeable future: discrete gravity observations at a mean flight height of 2 to 6 km above mean sea level with minimum
spatial resolution of 2.5 arcmin and a noise level of 1.5 mGal. Numerical results for both approaches are presented and discussed.
The direct approach can successfully be used for the downward continuation of airborne potential without any numerical instabilities
associated with the inverse approach. In addition to these two-step approaches, a one-step procedure is also discussed. This
procedure is based on a direct relationship between gravity disturbances at flight level and the disturbing gravity potential
at sea level. This procedure provided the best results in terms of accuracy, stability and numerical efficiency. As a general
result, numerically stable downward continuation of airborne gravity data can be seen as another advantage of airborne gravimetry
in the field of geoid determination.
Received: 6 June 2001 / Accepted: 3 January 2002 相似文献
2.
The long-wavelength geoid errors on large-scale geoid solutions, and the use of modified kernels to mitigate these effects,
are studied. The geoid around the Nordic area, from Greenland to the Ural mountains, is considered. The effect of including
additional gravity data around the Nordic/Baltic land area, originating from both marine, satellite and ground-based measurements,
is studied. It is found that additional data appear to increase the noise level in computations, indicating the presence of
systematic errors. Therefore, the Wong–Gore modification to the Stokes kernel is applied. This method of removing lower-order
terms in the Stokes kernel appears to improve the geoid. The best fit to the global positioning system (GPS) leveling points
is obtained with a degree of modification of approximately 30. In addition to the study of modification errors, the results
of different methods of combining satellite altimetry gravity and other gravimetry are presented. They all gave comparable
results, at the 6-cm level, when evaluated for the Nordic GPS networks. One dimensional (1-D) and 2-D fast Fourier transform
(FFT) methods are also compared. It is shown that even though methods differ by up to 6 cm, the fit to the GPS is essentially
the same. A surprising conclusion is that the addition of more data does not always produce a better geoid, illustrating the
danger of systematic errors in data.
Received: 4 July 2001 / Accepted: 21 February 2002 相似文献
3.
Three different methods of handling topography in geoid determination were investigated. The first two methods employ the
residual terrain model (RTM) remove–restore technique, yielding the quasigeoid, whereas the third method uses the classical
Helmert condensation method, yielding the geoid. All three methods were used with the geopotential model Earth Gravity Model
(1996) (EGM96) as a reference, and the results were compared to precise global positioning system (GPS) levelling networks
in Scandinavia. An investigation of the Helmert method, focusing on the different types of indirect effects and their effects
on the geoid, was also carried out. The three different methods used produce almost identical results at the 5-cm level, when
compared to the GPS levelling networks. However, small systematic differences existed.
Received: 18 March 1999 / Accepted: 21 March 2000 相似文献
4.
Improved convergence rates for the truncation error in gravimetric geoid determination 总被引:2,自引:2,他引:0
When Stokes's integral is used over a spherical cap to compute a gravimetric estimate of the geoid, a truncation error results
due to the neglect of gravity data over the remainder of the Earth. Associated with the truncation error is an error kernel
defined over these two complementary regions. An important observation is that the rate of decay of the coefficients of the
series expansion for the truncation error in terms of Legendre polynomials is determined by the smoothness properties of the
error kernel. Previously published deterministic modifications of Stokes's integration kernel involve either a discontinuity
in the error kernel or its first derivative at the spherical cap radius. These kernels are generalised and extended by constructing
error kernels whose derivatives at the spherical cap radius are continuous up to an arbitrary order. This construction is
achieved by smoothly continuing the error kernel function into the spherical cap using a suitable degree polynomial. Accordingly,
an improved rate of convergence of the spectral series representation of the truncation error is obtained.
Received: 21 April 1998 / Accepted: 4 October 1999 相似文献
5.
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 相似文献
6.
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 相似文献
7.
W. E. Featherstone J. F. Kirby A. H. W. Kearsley J. R. Gilliland G. M. Johnston J. Steed R. Forsberg M. G. Sideris 《Journal of Geodesy》2001,75(5-6):313-330
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.
J. F. Kirby 《Journal of Geodesy》2003,77(7-8):433-439
The geoid gradient over the Darling Fault in Western Australia is extremely high, rising by as much as 38 cm over only 2 km. This poses problems for gravimetric-only geoid models of the area, whose frequency content is limited by the spatial distribution of the gravity data. The gravimetric-only version of AUSGeoid98, for instance, is only able to resolve 46% of the gradient across the fault. Hence, the ability of GPS surveys to obtain accurate orthometric heights is reduced. It is described how further gravity data were collected over the Darling Fault, augmenting the existing gravity observations at key locations so as to obtain a more representative geoid gradient. As many of the gravity observations were collected at stations with a well-known GRS80 ellipsoidal height, the opportunity arose to compute a geoid model via both the Stokes and the Hotine approaches. A scheme was devised to convert free-air anomaly data to gravity disturbances using existing geoid models, followed by a Hotine integration to geoid heights. Interestingly, these results depended very weakly upon the choice of input geoid model. The extra gravity data did indeed improve the fit of the computed geoid to local GPS/Australian Height Datum (AHD) observations by 58% over the gravimetric-only AUSGeoid98. While the conventional Stokesian approach to geoid determination proved to be slightly better than the Hotine method, the latter still improved upon the gravimetric-only AUSGeoid98 solution, supporting the viability of conducting gravity surveys with GPS control for the purposes of geoid determination.
AcknowledgementsThe author would like to thank Will Featherstone, Ron Gower, Ron Hackney, Linda Morgan, Geoscience Australia, Scripps Oceanographic Institute and the three anonymous reviewers of this paper. This research was funded by the Australian Research Council. 相似文献
9.
L. E. Sjöberg 《Journal of Geodesy》2001,75(5-6):283-290
The topographic and atmospheric effects of gravimetric geoid determination by the modified Stokes formula, which combines
terrestrial gravity and a global geopotential model, are presented. Special emphasis is given to the zero- and first-degree
effects. The normal potential is defined in the traditional way, such that the disturbing potential in the exterior of the
masses contains no zero- and first-degree harmonics. In contrast, it is shown that, as a result of the topographic masses,
the gravimetric geoid includes such harmonics of the order of several centimetres. In addition, the atmosphere contributes
with a zero-degree harmonic of magnitude within 1 cm.
Received: 5 November 1999 / Accepted: 22 January 2001 相似文献
10.
GPS-based precise orbit determination of the very low Earth-orbiting gravity mission GOCE 总被引:5,自引:0,他引:5
A prerequisite for the success of future gravity missions like the European Gravity field and steady-state Ocean Circulation
Explorer (GOCE) is a precise orbit determination (POD). A detailed simulation study has been carried out to assess the achievable
orbit accuracy based on satellite-to-satellite tracking (SST) by the US global positioning system (GPS) and in conjunction
the implications for gravity field determination. An orbit accuracy at the few centimeter level seems possible, sufficient
to support the GOCE gravity mission and in particular its gravity gradiometer.
Received: 21 January 2000 / Accepted: 4 July 2000 相似文献
11.
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 相似文献
12.
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 相似文献
13.
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 相似文献
14.
P. Moore 《Journal of Geodesy》2001,75(5-6):241-254
Dual satellite crossovers (DXO) between the two European Remote Sensing satellites ERS-1 and ERS-2 and TOPEX/Poseidon are
used to (1) refine the Earth's gravity field and (2) extend the study of the ERS-2 altimetric range stability to cover the
first four years of its operation. The enhanced gravity field model, AGM-98, is validated by several methodologies and will
be shown to provide, in particular, low geographically correlated orbital error for ERS-2. For the ERS-2 altimetric range
study, TOPEX/Poseidon is first calibrated through comparison against in situ tide gauge data. A time series of the ERS-2 altimeter
bias has been recovered along with other geophysical correction terms using tables for bias jumps in the range measurements
at the single point target response (SPTR) events. On utilising the original version of the SPTR tables the overall bias drift
is seen to be 2.6±1.0 mm/yr with an RMS of fit of 12.2 mm but with discontinuities at the centimetre level at the SPTR events.
On utilising the recently released revised tables, SPTR2000, the drift is better defined at 2.4±0.6 mm/yr with the RMS of
fit reduced to 3.7 mm. Investigations identify the sea-state bias as a source of error with corrections affecting the overall
drift by close to 1.2 mm/yr.
Received: 25 May 2000 / Accepted: 24 January 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.
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 相似文献
17.
Rapid GPS ambiguity resolution for short and long baselines 总被引:3,自引:0,他引:3
A method of quick initial carrier cycle ambiguity resolution is described. The method applies to high-quality dual-band global
positioning system observations. Code measurements on both frequencies must be available. The rapidity of the method is achieved
through smoothing pseudoranges by phase observables and forming linear combinations between the phase observables. Two cases
are investigated. Case 1: ionospheric bias is neglected (short distances); and case 2: the bias is taken into account (longer
distances, more than, say, 10 km). The method was tested on six baselines, from 1 to 31 km long. In most cases, single-epoch
ambiguity resolution was achieved.
Received: 6 October 1999 / Accepted: 4 March 2002 相似文献
18.
The latest gravimetric geoid model for Japan, JGEOID2000, was successfully combined with the nationwide net of GPS at benchmarks,
yielding a new hybrid geoid model for Japan, GSIGEO2000. The least-squares collocation (LSC) method was applied as an interpolation
for fitting JGEOID2000 to the GPS/leveling geoid undulations. The GPS/leveling geoid undulation data were reanalyzed in advance,
in terms of three-dimensional positions from GPS and orthometric heights from leveling. The new hybrid geoid model is, therefore,
compatible with the new Japanese geodetic reference frame. GSIGEO2000 was evaluated internally and independently and the precision
was estimated at 4 cm throughout nearly the whole region.
Received: 15 October 2001 / Accepted: 27 March 2002
Acknowledgments. Messrs. Toshio Kunimi and Tadashi Saito at the Third Geodetic Division of the Geographical Survey Institute (GSI) mainly
carried out the computations of most of the updated leveled heights. With regard to the reanalysis of GPS data, the discussions
with Messrs. Yuki Hatanaka and Shoichi Matsumura of GSI were of great help in building the analysis strategy. Messrs. Kazuyuki
Tanaka and Hiromi Shigematsu collaborated in the preparatory stages of GPS data computation. The authors' thanks are extended
to these colleagues. Some plots were made by GMT software (Wessel and Smith 1991).
Correspondence to: Y. Kuroishi 相似文献
19.
R. Dach G. Beutler U. Hugentobler S. Schaer T. Schildknecht T. Springer G. Dudle L. Prost 《Journal of Geodesy》2003,77(1-2):1-14
A joint time-transfer project between the Astronomical Institute of the University of Berne (AIUB) and the Swiss Federal
Office of Metrology and Accreditation (METAS) was initiated to investigate the power of the time transfer using GPS carrier
phase observations. Studies carried out in the context of this project are presented. The error propagation for the time-transfer
solution using GPS carrier phase observations was investigated. To this purpose a simulation study was performed. Special
interest was focussed on errors in the vertical component of the station position, antenna phase-center variations and orbit
errors. A constant error in the vertical component introduces a drift in the time-transfer results for long baselines in east–west
directions. The simulation study was completed by investigating the profit for time transfer when introducing the integer
carrier phase ambiguities from a double-difference solution. This may reduce the drift in the time-transfer results caused
by constant vertical error sources. The results from the present time-transfer solution are shown in comparison to results
obtained with independent time-transfer techniques. The interpretation of the comparison benefits from the investigations
of the error propagation study. Two types of solutions are produced on a regular basis at AIUB: one based on the rapid orbits
from CODE, the other on the CODE final orbits. The rapid solution is available the day after the observations and has nearly
the same quality as the final solution, which has a latency of about one week. The differences between these two solutions
are below the nanosecond level. The differences from independent time-transfer techniques such as TWSTFT (two-way satellite
time and frequency transfer) are a few nanoseconds for both products.
Received: 15 November 2001 / Accepted: 6 September 2002
Correspondence to:R. Dach 相似文献
20.
Random simulation and GPS decorrelation 总被引:13,自引:1,他引:13
Peiliang Xu 《Journal of Geodesy》2001,75(7-8):408-423
(i) A random simulation approach is proposed, which is at the centre of a numerical comparison of the performances of different
GPS decorrelation methods. The most significant advantage of the approach is that it does not depend on nor favour any particular
satellite–receiver geometry and weighting system. (ii) An inverse integer Cholesky decorrelation method is proposed, which
will be shown to out-perform the integer Gaussian decorrelation and the Lenstra, Lenstra and Lovász (LLL) algorithm, and thus
indicates that the integer Gaussian decorrelation is not the best decorrelation technique and that further improvement is
possible. (iii) The performance study of the LLL algorithm is the first of its kind and the results have shown that the algorithm
can indeed be used for decorrelation, but that it performs worse than the integer Gaussian decorrelation and the inverse integer
Cholesky decorrelation. (iv) Simulations have also shown that no decorrelation techniques available to date can guarantee
a smaller condition number, especially in the case of high dimension, although reducing the condition number is the goal of
decorrelation.
Received: 26 April 2000 / Accepted: 5 March 2001 相似文献