首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 24 毫秒
1.
 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  相似文献   

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

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

4.
 Two long time series were analysed: the C01 series of the International Earth Rotation Service and the pole series obtained by re-analysis of the classical astronomical observations using the HIPPARCOS reference frame. The linear drift of the pole was determined to be 3.31 ± 0.05 milliarcseconds/year towards 76.1 ± 0.80° west longitude. For the least-squares fit the a priori correlations between simultaneous pole coordinates x p , y p were taken into account, and the weighting function was calculated by estimating empirical variance components. The decadal variations of the pole path were investigated by Fourier and wavelet analysis. Using sliding windows, the periods and amplitudes of the Chandler wobble and annual wobble were determined. Typical periods in the variable Chandler wobble and annual wobble parameters were obtained from wavelet analyses. Received: 21 January 2000 / Accepted: 28 August 2000  相似文献   

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

6.
Recently, effective atmospheric-angular-momentum (AAM) functions as calculated from National Centers for Environmental Prediction (NCEP) (formerly National Meteorological Center, NMC) and National Center for Atmospheric Research (NCAR) Reanalyses have become available for the years 1958 to 1998. Concerning the wind terms, the top level in the atmosphere used here is 10 hPa. Compared with earlier NMC model versions, which incorporate wind fields up to 100 hPa since 1976 and up to 50 hPa since 1981, the reanalyses have produced improved data series over a longer period than before. The axial AAM component χ3 is associated with changes in length of day (LOD). Motivated by better quality and continuity of the series AAM (NCEP) Reanalysis, the problem of the seasonal imbalances in the solid Earth–atmosphere axial angular momentum budget is re-examined. To assess better the estimates of the annual and semiannual oscillations in LOD and AAM and of the residual oscillations derived as difference series between LOD and AAM, the series of LOD data from three analysis centers [International Earth Rotation Service (IERS), GeoForschungZentrum Potsdam (GFZ) and Jet Propulsion Laboratory Pasadena (JPL)] and of AAM data in terms of χ3(W), χ3(P) and χ3(P+IB) from four meteorological centers [NCEP, Japan Meteorological Agency (JMA), European Centre for Medium-Range Weather Forecasts (ECMWF) and the UK Meteorological Office (UKMO)] are used in this study. The main analysis steps were removing gaps, filtering out the seasonal oscillations, calculating optimal estimates of the parameters of the oscillations and calculating the difference series between the LOD and AAM systems as well as the residuals in the axial angular momentum budget in the LOD–AAM systems. The results derived as difference series between the different LOD, AAM and LOD–AAM systems show to what extent the variations reflect systematic differences and significant signals, respectively, which is important for future activities in this field. Received: 2 February 1999 / Accepted: 30 November 1999  相似文献   

7.
The TOPEX/Poseidon (T/P) satellite alti- meter mission marked a new era in determining the geopotential constant W 0. On the basis of T/P data during 1993–2003 (cycles 11–414), long-term variations in W 0 have been investigated. The rounded value W 0 = 62636856.0 ± 0.5) m 2 s −2 has already been adopted by the International Astronomical Union for the definition of the constant L G = W 0/c 2 = 6.969290134 × 10−10 (where c is the speed of light), which is required for the realization of the relativistic atomic time scale. The constant L G , based on the above value of W 0, is also included in the 2003 International Earth Rotation and Reference Frames Service conventions. It has also been suggested that W 0 is used to specify a global vertical reference system (GVRS). W 0 ensures the consistency with the International Terrestrial Reference System, i.e. after adopting W 0, along with the geocentric gravitational constant (GM), the Earth’s rotational velocity (ω) and the second zonal geopotential coefficient (J 2) as primary constants (parameters), then the ellipsoidal parameters (a,α) can be computed and adopted as derived parameters. The scale of the International Terrestrial Reference Frame 2000 (ITRF2000) has also been specified with the use of W 0 to be consistent with the geocentric coordinate time. As an example of using W 0 for a GVRS realization, the geopotential difference between the adopted W 0 and the geopotential at the Rimouski tide-gauge point, specifying the North American Vertical Datum 1988 (NAVD88), has been estimated.  相似文献   

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

9.
 Ten days of GPS data from 1998 were processed to determine how the accuracy of a derived three-dimensional relative position vector between GPS antennas depends on the chord distance (denoted L) between these antennas and on the duration of the GPS observing session (denoted T). It was found that the dependence of accuracy on L is negligibly small when (a) using the `final' GPS satellite orbits disseminated by the International GPS Service, (b) fixing integer ambiguities, (c) estimating appropriate neutral-atmosphere-delay parameters, (d) 26 km ≤ L ≤ 300 km, and (e) 4 h ≤T ≤ 24 h. Under these same conditions, the standard error for the relative position in the north–south dimension (denoted S n and expressed in mm) is adequately approximated by the equation S n =k n /T  0.5 with k n =9.5 ± 2.1 mm · h0.5 and T expressed in hours. Similarly, the standard errors for the relative position in the east–west and in the up-down dimensions are adequately approximated by the equations S e =k e /T  0.5 and S u =k u /T  0.5, respectively, with k e =9.9 ± 3.1 mm · h0.5 and k u =36.5 ± 9.1 mm · h0.5. Received: 5 February 2001 / Accepted: 14 May 2001  相似文献   

10.
 The solutions of the CODE Analysis Center submitted to the IGS, the International Global Position System (GPS) Service for Geodynamics, are based on three days of observation of about 80–100 stations of the IGS network. The Earth rotation parameters (ERPs) are assumed to vary linearly over the three days with respect to an a priori model. Continuity at the day boundaries as well as the continuity of the first derivatives are enforced by constraints. Since early April 1995 CODE has calculated a new ERP series with an increased time resolution of 2 hours. Again continuity is enforced at the 2-hours-interval boundaries. The analysis method is described, particularly how to deal with retrograde diurnal terms in the ERP series which may not be estimated with satellite geodetic methods. The results obtained from the first year of data covered by the time series (time interval from 4 April 1995 to 30 June 1996) are also discussed. The series is relatively homogeneous in the sense of the used orbit model and the a priori model for the ERPs. The largest source of excitation at daily and sub-daily periods is likely to be the effect of the ocean tides. There is good agreement between the present results and Topex/Poseidon ocean tide models, as well as with models based on Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data. Non-oceanic periodic variations are also observed in the series. Their origin is most probably a consequence of the GPS solution strategy; other possible sources are the atmospheric tides. Received: 13 July 1999 / Accepted: 21 March 2000  相似文献   

11.
Starlette was launched in 1975 in order to study temporal variations in the Earth’s gravity field; in particular, tidal and Earth rotation effects. For the period April 1983 to April 1984 over12,700 normal points of laser ranging data to Starlette have been sub-divided into49 near consecutive 5–6 day arcs. Normal equations for each arc as obtained from a least-squares data reduction procedure, were solved for ocean tidal parameters along with other geodetic and geodynamic parameters. The tidal parameters are defined relative to Wahr’s body tides and Wahr’s nutation model and show fair agreement with other satellite derived results and those obtained from spherical harmonic decomposition of global ocean tidal models.  相似文献   

12.
Green's function solution to spherical gradiometric boundary-value problems   总被引：1，自引：1，他引：1  
 Three independent gradiometric boundary-value problems (BVPs) with three types of gradiometric data, {Γ rr }, {Γ r θ,Γ r λ} and {Γθθ−Γλλ,Γθλ}, prescribed on a sphere are solved to determine the gravitational potential on and outside the sphere. The existence and uniqueness conditions on the solutions are formulated showing that the zero- and the first-degree spherical harmonics are to be removed from {Γ r θ,Γ r λ} and {Γθθ−Γλλ,Γθλ}, respectively. The solutions to the gradiometric BVPs are presented in terms of Green's functions, which are expressed in both spectral and closed spatial forms. The logarithmic singularity of the Green's function at the point ψ=0 is investigated for the component Γ rr . The other two Green's functions are finite at this point. Comparisons to the paper by van Gelderen and Rummel [Journal of Geodesy (2001) 75: 1–11] show that the presented solution refines the former solution. Received: 3 October 2001 / Accepted: 4 October 2002  相似文献   

13.
Using the ΔT (integrated variation of the Earth's rotation measured in terestrial time) series (1891.5–1955.5) derived from lunar occultation observations and the UT1–UTC (universal time–coordinated universal time) series (1955.5–1997.5) of the Bureau International de L'Heure/International Earth Rotation Service, a new ΔLOD (variation of the length of day) series in monthly intervals from 1892.0 to 1997.0 is calculated. Using digital filtering, the interannual and decadal components of the ΔLOD series are separated and then compared with those inferred from other geophysical quantities. It is shown that, on the interannual time scale, atmospheric processes can play an important role in exciting astronomical ΔLOD. However, the main oscillation with a mean period of about 5.8 years and peak-to-peak amplitude of about 0.3 ms in the residuals of ΔLOD(Astr) −ΔLOD(Wind) for 1968.0–1997.0 suggests that about half of the amplitude in astronomical ΔLOD must be excited by other geophysical processes, while on the decadal time scale the atmospheric excitation is too small. Geomagnetic core–mantle coupling may be a plausible source of the excitation of ΔLOD on the decadal time scale, but the geomagnetic data are still insufficient and an improved model of core–mantle coupling is required. Received: 3 April 1998 / Accepted: 31 May 1999  相似文献   

14.
The contribution of the International VLBI Service for Geodesy and Astrometry (IVS) to the ITRF2005 (International Terrestrial Reference Frame 2005) has been computed by the IVS Analysis Coordinator’s office at the Geodetic Institute of the University of Bonn, Germany. For this purpose the IVS Analysis Centres (ACs) provided datum-free normal equation matrices in Solution INdependent EXchange (SINEX) format for each 24 h observing session to be combined on a session-by-session basis by a stacking procedure. In this process, common sets of parameters, transformed to identical reference epochs and a prioris, and especially representative relative weights have been taken into account for each session. In order to assess the quality of the combined IVS files, Earth orientation parameters (EOPs) and scaling factors have been derived from the combined normal equation matrices. The agreement of the EOPs of the combined normal equation matrices with those of the individual ACs in terms of weighted root mean square (WRMS) is in the range of 50–60 μas for the two polar motion components and about 3 μs for UT1−UTC. External comparisons with International GNSS Serive (IGS) polar motion components is at the level of 130–170 μas and 21 μs/day for length of day (LOD). The scale of the terrestrial reference frame realized through the IVS SINEX files agrees with ITRF2000 at the level of 0.2 ppb.  相似文献   

15.
Prediction of Earth orientation parameters by artificial neural networks   总被引：3，自引：1，他引：3  
 Earth orientation parameters (EOPs) [polar motion and length of day (LOD), or UT1–UTC] were predicted by artificial neural networks. EOP series from various sources, e.g. the C04 series from the International Earth Rotation Service and the re-analysis optical astrometry series based on the HIPPARCOS frame, served for training the neural network for both short-term and long-term predictions. At first, all effects which can be described by functional models, e.g. effects of the solid Earth tides and the ocean tides or seasonal atmospheric variations of the EOPs, were removed. Only the differences between the modeled and the observed EOPs, i.e. the quasi-periodic and irregular variations, were used for training and prediction. The Stuttgart neural network simulator, which is a very powerful software tool developed at the University of Stuttgart, was applied to construct and to validate different types of neural networks in order to find the optimal topology of the net, the most economical learning algorithm and the best procedure to feed the net with data patterns. The results of the prediction were analyzed and compared with those obtained by other methods. The accuracy of the prediction is equal to or even better than that by other prediction methods. Received: 6 February 2001 / Accepted: 23 October 2001  相似文献   

16.
World Geodetic Datum 2000   总被引：7，自引：1，他引：6  
 Based on the current best estimates of fundamental geodetic parameters {W 0,GM,J 2,Ω} the form parameters of a Somigliana-Pizzetti level ellipsoid, namely the semi-major axis a and semi-minor axis b (or equivalently the linear eccentricity ) are computed and proposed as a new World Geodetic Datum 2000. There are six parameters namely the four fundamental geodetic parameters {W 0,GM,J 2,Ω} and the two form parameters {a,b} or {a,ɛ}, which determine the ellipsoidal reference gravity field of Somigliana-Pizzetti type constraint to two nonlinear condition equations. Their iterative solution leads to best estimates a=(6 378 136.572±0.053)m, b=(6 356 751.920 ± 0.052)m, ɛ=(521 853.580±0.013)m for the tide-free geoide of reference and a=(6 378 136.602±0.053)m, b=(6 356 751.860±0.052)m, ɛ=(521 854.674 ± 0.015)m for the zero-frequency tide geoid of reference. The best estimates of the form parameters of a Somigliana-Pizzetti level ellipsoid, {a,b}, differ significantly by −0.39 m, −0.454 m, respectively, from the data of the Geodetic Reference System 1980. Received: 1 February 1999 / Accepted: 31 August 1999  相似文献   

17.
 Length-of-day (LOD) estimates from seven Global Positioning System (GPS) and three satellite laser ranging (SLR) analysis centers were combined into an even-spaced time series for a 27-month period (1996–1998). This time series was compared to the multi-technique Earth-orientation-parameter (EOP) combined solution (C04) derived at the Central Bureau of the International Earth Rotation Service (IERS/CB). Due to inhomogeneities in the different series derived from the various techniques (time, length, quality, and spatial resolution), the concept of a combined solution is justified. The noise behavior in LOD for different techniques varies with frequency; the data series were divided into frequency windows after removing both biases and trends. Different weight factors were assigned in each window, discriminating by technique, and produced one-technique combined solutions. Finally, these one-technique combined solutions were combined to obtain the final multi-technique solution. The LOD combined time series obtained by the present method based on the frequency windows combined series (FWCS) is very close to the IERS C04 solution. It could be useful to generate a new LOD reference time series to be used in the study of high-frequency variations of Earth rotation. Received: 28 March 2000 / Accepted: 15 February 2001  相似文献   

18.
 Harmonic coefficients of the 2nd degree are separated into the invariant quantitative (the 2nd-degree variance) and the qualitative (the standardized harmonic coefficients) characteristics of the behavior of the potential V 2(t). On this basis the evolution of the Earth's dynamical figure is described as a solution of the time-dependent eigenvalues–eigenvectors problem in the canonical form. Such a canonical quadratic form is defined only by temporal variations of the harmonic coefficients and always remains finite, even within an infinite time interval. An additional condition for the correction or the determination of temporal variations of the 2nd degree is obtained. Temporal variations of the fully normalized sectorial harmonic coefficients are estimated in addition to ˙Cˉ 20, ˙Cˉ 21, and ˙Sˉ 21 of the EGM96 gravity model. In addition, a non-linear hyperbolic model for Cˉ 2m (t), Sˉ 2m (t) is constructed. The trigonometric form of the hyperbolic model leads to the consideration of the potential V 2(ψ) instead of V 2(t) within the closed interval −π/2≤ψ≤+π/2. Thus, it is possible to evaluate the global trend of V 2(t), the Earth's principal axes and the differences of the moments of inertia within the whole infinite time interval. Received: 25 September 1998 / Accepted: 28 June 2000  相似文献   

19.
 The weighted Procrustes algorithm is presented as a very effective tool for solving the three-dimensional datum transformation problem. In particular, the weighted Procrustes algorithm does not require any initial datum parameters for linearization or any iteration procedure. As a closed-form algorithm it only requires the values of Cartesian coordinates in both systems of reference. Where there is some prior information about the variance–covariance matrix of the two sets of Cartesian coordinates, also called pseudo-observations, the weighted Procrustes algorithm is able to incorporate such a quality property of the input data by means of a proper choice of weight matrix. Such a choice is based on a properly designed criterion matrix which is discussed in detail. Thanks to the weighted Procrustes algorithm, the problem of incorporating the stochasticity measures of both systems of coordinates involved in the seven parameter datum transformation problem [conformal group ℂ7(3)] which is free of linearization and any iterative procedure can be considered to be solved. Illustrative examples are given. Received: 7 January 2002 / Accepted: 9 September 2002 Correspondence to: E. W. Grafarend  相似文献   

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

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

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