首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
For science applications of the gravity recovery and climate experiment (GRACE) monthly solutions, the GRACE estimates of \(C_{20}\) (or \(J_{2}\)) are typically replaced by the value determined from satellite laser ranging (SLR) due to an unexpectedly strong, clearly non-geophysical, variation at a period of \(\sim \)160 days. This signal has sometimes been referred to as a tide-like variation since the period is close to the perturbation period on the GRACE orbits due to the spherical harmonic coefficient pair \(C_{22}/S_{22}\) of S2 ocean tide. Errors in the S2 tide model used in GRACE data processing could produce a significant perturbation to the GRACE orbits, but it cannot contribute to the \(\sim \)160-day signal appearing in \(C_{20}\). Since the dominant contribution to the GRACE estimate of \(C_{20}\) is from the global positioning system tracking data, a time series of 138 monthly solutions up to degree and order 10 (\(10\times 10\)) were derived along with estimates of ocean tide parameters up to degree 6 for eight major tides. The results show that the \(\sim \)160-day signal remains in the \(C_{20}\) time series. Consequently, the anomalous signal in GRACE \(C_{20}\) cannot be attributed to aliasing from the errors in the S2 tide. A preliminary analysis of the cross-track forces acting on GRACE and the cross-track component of the accelerometer data suggests that a temperature-dependent systematic error in the accelerometer data could be a cause. Because a wide variety of science applications relies on the replacement values for \(C_{20}\), it is essential that the SLR estimates are as reliable as possible. An ongoing concern has been the influence of higher degree even zonal terms on the SLR estimates of \(C_{20}\), since only \(C_{20}\) and \(C_{40}\) are currently estimated. To investigate whether a better separation between \(C_{20}\) and the higher-degree terms could be achieved, several combinations of additional SLR satellites were investigated. In addition, a series of monthly gravity field solutions (\(60\times 60\)) were estimated from a combination of GRACE and SLR data. The results indicate that the combination of GRACE and SLR data might benefit the resonant orders in the GRACE-derived gravity fields, but it appears to degrade the recovery of the \(C_{20}\) variations. In fact, the results suggest that the poorer recovery of \(C_{40}\) by GRACE, where the annual variation is significantly underestimated, may be affecting the estimates of \(C_{20}\). Consequently, it appears appropriate to continue using the SLR-based estimates of \(C_{20}\), and possibly also \(C_{40}\), to augment the existing GRACE mission.  相似文献   

2.
Large-scale mass redistribution in the terrestrial water storage (TWS) leads to changes in the low-degree spherical harmonic coefficients of the Earth’s surface mass density field. Studying these low-degree fluctuations is an important task that contributes to our understanding of continental hydrology. In this study, we use global GNSS measurements of vertical and horizontal crustal displacements that we correct for atmospheric and oceanic effects, and use a set of modified basis functions similar to Clarke et al. (Geophys J Int 171:1–10, 2007) to perform an inversion of the corrected measurements in order to recover changes in the coefficients of degree-0 (hydrological mass change), degree-1 (centre of mass shift) and degree-2 (flattening of the Earth) caused by variations in the TWS over the period January 2003–January 2015. We infer from the GNSS-derived degree-0 estimate an annual variation in total continental water mass with an amplitude of \((3.49 \pm 0.19) \times 10^{3}\) Gt and a phase of \(70^{\circ } \pm 3^{\circ }\) (implying a peak in early March), in excellent agreement with corresponding values derived from the Global Land Data Assimilation System (GLDAS) water storage model that amount to \((3.39 \pm 0.10) \times 10^{3}\) Gt and \(71^{\circ } \pm 2^{\circ }\), respectively. The degree-1 coefficients we recover from GNSS predict annual geocentre motion (i.e. the offset change between the centre of common mass and the centre of figure) caused by changes in TWS with amplitudes of \(0.69 \pm 0.07\) mm for GX, \(1.31 \pm 0.08\) mm for GY and \(2.60 \pm 0.13\) mm for GZ. These values agree with GLDAS and estimates obtained from the combination of GRACE and the output of an ocean model using the approach of Swenson et al. (J Geophys Res 113(B8), 2008) at the level of about 0.5, 0.3 and 0.9 mm for GX, GY and GZ, respectively. Corresponding degree-1 coefficients from SLR, however, generally show higher variability and predict larger amplitudes for GX and GZ. The results we obtain for the degree-2 coefficients from GNSS are slightly mixed, and the level of agreement with the other sources heavily depends on the individual coefficient being investigated. The best agreement is observed for \(T_{20}^C\) and \(T_{22}^S\), which contain the most prominent annual signals among the degree-2 coefficients, with amplitudes amounting to \((5.47 \pm 0.44) \times 10^{-3}\) and \((4.52 \pm 0.31) \times 10^{-3}\) m of equivalent water height (EWH), respectively, as inferred from GNSS. Corresponding agreement with values from SLR and GRACE is at the level of or better than \(0.4 \times 10^{-3}\) and \(0.9 \times 10^{-3}\) m of EWH for \(T_{20}^C\) and \(T_{22}^S\), respectively, while for both coefficients, GLDAS predicts smaller amplitudes. Somewhat lower agreement is obtained for the order-1 coefficients, \(T_{21}^C\) and \(T_{21}^S\), while our GNSS inversion seems unable to reliably recover \(T_{22}^C\). For all the coefficients we consider, the GNSS-derived estimates from the modified inversion approach are more consistent with the solutions from the other sources than corresponding estimates obtained from an unconstrained standard inversion.  相似文献   

3.
The geodetic VLBI community began using VLBA antennas in 1989 for geodesy and astrometry. We examine how usage of the VLBA has improved the celestial reference frame, the terrestrial reference frame, and Earth orientation parameters. Without the VLBA, ICRF2 would have had only 1011 sources instead of 3414. ICRF3 will contain at least 4121 sources, with approximately 70 % or more coming exclusively from VLBA astrometry and geodesy sessions. The terrestrial reference frame is also more stable and precise due to VLBA geodesy sessions. Approximately two dozen geodesy stations that have participated in VLBA sessions show average position formal errors that are \(\sim \)13–14 % better in the horizontal components and \(\sim \)5 % better in the vertical component than would be expected solely from the increased number of observations. Also the Earth orientation parameters obtained from the RDV sessions represent the most accurate EOP series of any of the long-term VLBI session types.  相似文献   

4.
We analyze the high-resolution dilatation data for the October 2013 \(M_w\) 6.2 Ruisui, Taiwan, earthquake, which occurred at a distance of 15–20 km away from a Sacks–Evertson dilatometer network. Based on well-constrained source parameters (\(\hbox {strike}=217^\circ \), \(\hbox {dip}=48^\circ \), \(\hbox {rake}=49^\circ \)), we propose a simple rupture model that explains the permanent static deformation and the dynamic vibrations at short period (\(\sim \)3.5–4.5 s) for most of the four sites with less than 20 % of discrepancies. This study represents a first attempt of modeling simultaneously the dynamic and static crustal strain using dilatation data. The results illustrate the potential for strain recordings of high-frequency seismic waves in the near-field of an earthquake to add constraints on the properties of seismic sources.  相似文献   

5.
The consistent estimation of terrestrial reference frames (TRF), celestial reference frames (CRF) and Earth orientation parameters (EOP) is still an open subject and offers a large field of investigations. Until now, source positions resulting from Very Long Baseline Interferometry (VLBI) observations are not routinely combined on the level of normal equations in the same way as it is a common process for station coordinates and EOPs. The combination of source positions based on VLBI observations is now integrated in the IVS combination process. We present the studies carried out to evaluate the benefit of the combination compared to individual solutions. On the level of source time series, improved statistics regarding weighted root mean square have been found for the combination in comparison with the individual contributions. In total, 67 stations and 907 sources (including 291 ICRF2 defining sources) are included in the consistently generated CRF and TRF covering 30 years of VLBI contributions. The rotation angles \(A_1\), \(A_2\) and \(A_3\) relative to ICRF2 are ?12.7, 51.7 and 1.8 \({\upmu }\) as, the drifts \(D_\alpha \) and \(D_\delta \) are ?67.2 and 19.1 \(\upmu \) as/rad and the bias \(B_\delta \) is 26.1 \(\upmu \) as. The comparison of the TRF solution with the IVS routinely combined quarterly TRF solution shows no significant impact on the TRF, when the CRF is estimated consistently with the TRF. The root mean square value of the post-fit station coordinate residuals is 0.9 cm.  相似文献   

6.
We develop a slope correction model to improve the accuracy of mean sea surface topography models as well as marine gravity models. The correction is greatest above ocean trenches and large seamounts where the slope of the geoid exceeds 100  \(\upmu \) rad. In extreme cases, the correction to the mean sea surface height is 40 mm and the correction to the along-track altimeter slope is 1–2  \(\upmu \) rad which maps into a 1–2 mGal gravity error. Both corrections are easily applied using existing grids of sea surface slope from satellite altimetry.  相似文献   

7.
Three combined celestial pole offset (CPO) series computed at the Paris Observatory (C04), the United States Naval Observatory (USNO), and the International VLBI Service for Geodesy and Astrometry (IVS), as well as six free core nutation (FCN) models, were compared from different perspectives, such as stochastic and systematic differences, and FCN amplitude and phase variations. The differences between the C04 and IVS CPO series were mostly stochastic, whereas a low-frequency bias at the level of several tens of \(\upmu \)as was found between the C04 and USNO CPO series. The stochastic differences between the C04 and USNO series became considerably smaller when computed at the IVS epochs, which can indicate possible problems with the interpolation of the IVS data at the midnight epochs during the computation of the C04 and USNO series. The comparison of the FCN series showed that the series computed with similar window widths of 1.1–1.2 years were close to one another at a level of 10–20 \(\upmu \)as, whereas the differences between these series and the series computed with a larger window width of 4 and 7 years reached 100 \(\upmu \)as. The dependence of the FCN model on the underlying CPO series was investigated. The RMS differences between the FCN models derived from the C04, USNO, and IVS CPO series were at a level of approximately 15 \(\upmu \)as, which was considerably smaller than the differences among the CPO series. The analysis of the differences between the IVS, C04, and USNO CPO series suggested that the IVS series would be preferable for both precession-nutation and FCN-related studies.  相似文献   

8.
This paper describes the historical sea level data that we have rescued from a tide gauge, especially devised originally for geodesy. This gauge was installed in Marseille in 1884 with the primary objective of defining the origin of the height system in France. Hourly values for 1885–1988 have been digitized from the original tidal charts. They are supplemented by hourly values from an older tide gauge record (1849–1851) that was rediscovered during a survey in 2009. Both recovered data sets have been critically edited for errors and their reliability assessed. The hourly values are thoroughly analysed for the first time after their original recording. A consistent high-frequency time series is reported, increasing notably the length of one of the few European sea level records in the Mediterranean Sea spanning more than one hundred years. Changes in sea levels are examined, and previous results revisited with the extended time series. The rate of relative sea level change for the period 1849–2012 is estimated to have been \(1.08\pm 0.04\)  mm/year at Marseille, a value that is slightly lower but in close agreement with the longest time series of Brest over the common period ( \(1.26\pm 0.04\)  mm/year). The data from a permanent global positioning system station installed on the roof of the solid tide gauge building suggests a remarkable stability of the ground ( \(-0.04\pm 0.25\)  mm/year) since 1998, confirming the choice made by our predecessor geodesists in the nineteenth century regarding this site selection.  相似文献   

9.
Homogeneous reprocessing of GPS,GLONASS and SLR observations   总被引：3，自引：2，他引：1  
The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from \(-35\) and \(-38\)  mm to \(-12\) and \(-13\)  mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.  相似文献   

10.
Impact of seasonal station motions on VLBI UT1 intensives results   总被引：1，自引：1，他引：0  
UT1 estimates obtained from the very long baseline interferometry (VLBI) Intensives data depend on the station displacement model used during processing. In particular, because of seasonal variations, the instantaneous station position during the specific intensive session differs from the position predicted by the linear model generally used. This can cause systematic errors in UT1 Intensives results. In this paper, we first investigated the seasonal signal in the station displacements for the 5 VLBI antennas participating in UT1 Intensives observing programs, along with the 8 collocated GPS stations. It was found that a significant annual term is present in the time series for most stations, and its amplitude can reach 8 mm in the height component, and 2 mm in horizontal components. However, the annual signals found in the displacements of the collocated VLBI and GPS stations at some sites differ substantially in amplitude and phase. The semiannual harmonics are relatively small and unstable, and for most stations no prevailing signal was found in the corresponding frequency band. Then two UT1 Intensives series were computed with and without including the seasonal term found in the previous step in the station movement model. Comparison of these series has shown that neglecting the seasonal station position variations can cause a systematic error in UT1 estimates, which can exceed 1  $\upmu $ s, depending on the observing program.  相似文献   

11.
Deformations of radio telescopes used in geodetic and astrometric very long baseline interferometry (VLBI) observations belong to the class of systematic error sources which require correction in data analysis. In this paper we present a model for all path length variations in the geometrical optics of radio telescopes which are due to gravitational deformation. The Effelsberg 100 m radio telescope of the Max Planck Institute for Radio Astronomy, Bonn, Germany, has been surveyed by various terrestrial methods. Thus, all necessary information that is needed to model the path length variations is available. Additionally, a ray tracing program has been developed which uses as input the parameters of the measured deformations to produce an independent check of the theoretical model. In this program as well as in the theoretical model, the illumination function plays an important role because it serves as the weighting function for the individual path lengths depending on the distance from the optical axis. For the Effelsberg telescope, the biggest contribution to the total path length variations is the bending of the main beam located along the elevation axis which partly carries the weight of the paraboloid at its vertex. The difference in total path length is almost \(-\) 100 mm when comparing observations at 90 \(^\circ \) and at 0 \(^\circ \) elevation angle. The impact of the path length corrections is validated in a global VLBI analysis. The application of the correction model leads to a change in the vertical position of \(+120\)  mm. This is more than the maximum path length, but the effect can be explained by the shape of the correction function.  相似文献   

12.
Earth orientation parameters estimated from VLBI during the CONT11 campaign   总被引：1，自引：1，他引：0  
In this paper we investigate the accuracy of the earth orientation parameters (EOP) estimated from the continuous VLBI campaign CONT11. We first estimated EOP with daily resolution and compared these to EOP estimated from GNSS data. We find that the WRMS differences are about 31  $\upmu $ as for polar motion and 7  $\upmu $ s for length of day. This is about the precision we could expect, based on Monte Carlo simulations and the results of the previous CONT campaigns. We also estimated EOP with hourly resolution to study the sub-diurnal variations. The results confirm the results of previous studies, showing that the current IERS model for high-frequency EOP variations does not explain all the sub-diurnal variations seen in the estimated time series. We then compared our results to various empirical high-frequency EOP models. However, we did not find that any of these gave any unambiguous improvement. Several simulations testing the impact of various aspects of, e.g. the observing network were also made. For example, we made simulations assuming that all CONT11 stations were equipped with fast VLBI2010 antennas. We found that the WRMS error decreased by about a factor five compared to the current VLBI system. Furthermore, the simulations showed that it is very important to have a homogenous global distribution of the stations for achieving the highest precision for the EOP.  相似文献   

13.
Observations of spatio-temporal variations in the geopotential using the GRACE satellites have been used to estimate recent mass fluxes from polar ice sheets and glaciers. However, these estimates have not considered the potential bias associated with the migration of water that accompanies the ice melt. This migration is driven by the diminished gravitational attraction of the melting ice reservoir, and this migration, as well as the crustal loading it induces, will contribute to the observed geopotential anomaly. The extent to which this contribution contaminates the ice mass flux estimates depends on how far the smoothing filters applied to the GRACE data extend beyond the ice margins into the ocean. Using the Antarctic Peninsula as a case study, we estimate the magnitude of this bias for a range of melt areas and Gaussian smoothing filter radii. We conclude that GRACE estimates of ice mass loss over the Antarctic Peninsula are systematically overestimating the loss by up to 10  $\%$ for filter radii of less than 500 km.  相似文献   

14.
M-estimation with probabilistic models of geodetic observations   总被引：1，自引：1，他引：0  
The paper concerns \(M\) -estimation with probabilistic models of geodetic observations that is called \(M_{\mathcal {P}}\) estimation. The special attention is paid to \(M_{\mathcal {P}}\) estimation that includes the asymmetry and the excess kurtosis, which are basic anomalies of empiric distributions of errors of geodetic or astrometric observations (in comparison to the Gaussian errors). It is assumed that the influence function of \(M_{\mathcal {P}}\) estimation is equal to the differential equation that defines the system of the Pearson distributions. The central moments \(\mu _{k},\, k=2,3,4\) , are the parameters of that system and thus, they are also the parameters of the chosen influence function. The \(M_{\mathcal {P}}\) estimation that includes the Pearson type IV and VII distributions ( \(M_{\mathrm{PD(l)}}\) method) is analyzed in great detail from a theoretical point of view as well as by applying numerical tests. The chosen distributions are leptokurtic with asymmetry which refers to the general characteristic of empirical distributions. Considering \(M\) -estimation with probabilistic models, the Gram–Charlier series are also applied to approximate the models in question ( \(M_{\mathrm{G-C}}\) method). The paper shows that \(M_{\mathcal {P}}\) estimation with the application of probabilistic models belongs to the class of robust estimations; \(M_{\mathrm{PD(l)}}\) method is especially effective in that case. It is suggested that even in the absence of significant anomalies the method in question should be regarded as robust against gross errors while its robustness is controlled by the pseudo-kurtosis.  相似文献   

15.
This paper presents deformation analysis of Lake Urmia causeway (LUC) embankments in northwest Iran using observations from interferometry synthetic aperture radar (InSAR) and finite element model (FEM) simulation. 58 SAR images including 10 ALOS, 30 Envisat and 18 TerraSAR-X are used to assess settlement of the embankments during 2003–2013. The interferometric dataset includes 140 differential interferograms which are processed using InSAR time series technique of small baseline subset approach. The results show a clear indication of large deformation on the embankments with peak amplitude of \(>\) 50 mm/year in 2003–2010, increasing to \(>\!\!80\)  mm/year in 2012–2013 in the line of sight (LOS) direction from ground to the satellite. 2D decomposition of InSAR observations from Envisat and ALOS satellites that overlap in the years 2007–2010 shows that the rate of the vertical settlement and horizontal motion is not uniform along the embankments; Both eastern and western embankments show significant vertical motion, while horizontal motion plays a more significant role in eastern embankment than western embankment. The InSAR results are then used to simulate deformation using FEM at two cross-sections at the distance of 4 and 9 km from the most western edge of the LUC for which detailed stratigraphy data are available. Results suggest that consolidation due to dissipation of excess pore pressure in embankments can satisfactory predict settlement of the LUC embankments. Our numerical modeling indicates that nearly half of the consolidation since the construction time of the causeway 30 years ago has been done.  相似文献   

16.
The recent vertical displacement field of the Upper Rhine Graben (URG) located in the tri-national region between Germany, France and Switzerland is investigated using repeatedly measured leveling data. We estimate vertical displacement rates at leveling benchmarks by applying a kinematic network adjustment on more than 40,000 height differences measured by German, French and Swiss surveying agencies. Focusing on an optimal solution for the adjusted rates in the URG area also historical data (measured before 1900) are used, significantly increasing the time span of available measurements and the number of transnational connections between the three countries. To account for inhomogeneities apparent in the database, we apply an iterative variance component estimation within the adjustment procedure, particularly revealing more realistic information on the accuracy of the estimated rates. A special focus within our analysis is put on the statistical testing of gross errors in the observations and model-related errors at benchmarks with non-linear movement. As some of the estimated vertical rates behave significantly different compared to the vertical rates of adjacent benchmarks, a filtering of outliers is applied after the adjustment procedure. The resulting map of linear height changes in an area of 280 km in N–S and 230 km in E–W direction provides detailed insight into the recent vertical displacements of the URG and neighboring regions. In the German part of the study area, it was possible for the first time to consistently constrain an average subsidence rate of 0.5 mm/a ( \(\pm \) 0.2 mm/a) of the Graben interior w.r.t. the Black Forest. In addition to the tectonic displacements, some man-induced surface movements, e.g., caused by oil and groundwater extraction, are observed and discussed.  相似文献   

17.
We examine the relationship between source position stability and astrophysical properties of radio-loud quasars making up the International Celestial Reference Frame (ICRF2). Understanding this relationship is important for improving quasar selection and analysis strategies, and therefore reference frame stability. We construct flux density time series, known as light curves, for 95 of the most frequently observed ICRF2 quasars at both the 2.3 and 8.4 GHz geodetic very long baseline interferometry (VLBI) observing bands. Because the appearance of new quasar components corresponds to an increase in quasar flux density, these light curves alert us about potential changes in source structure before they appear in VLBI images. We test how source position stability depends on three astrophysical parameters: (1) flux density variability at X band; (2) time lag between flares in S and X bands; (3) spectral index root-mean-square (rms), defined as the variability in the ratio between S and X band flux densities. We find that the time lag between S and X band light curves provides a good indicator of position stability: sources with time lags $<$ 0.06 years are significantly more stable ( $>$ 20 % improvement in weighted rms) than sources with larger time lags. A similar improvement is obtained by observing sources with low $(<$ 0.12) spectral index variability. On the other hand, there is no strong dependence of source position stability on flux density variability in a single frequency band. These findings can be understood by interpreting the time lag between S and X band light curves as a measure of the size of the source structure. Monitoring of source flux density at multiple frequencies therefore appears to provide a useful probe of quasar structure on scales important to geodesy. The observed astrometric position of the brightest quasar component (the core) is known to depend on observing frequency. We show how multi-frequency flux density monitoring may allow the dependence on frequency of the relative core positions along the jet to be elucidated. Knowledge of the position–frequency relation has important implications for current and future geodetic VLBI programs, as well as the alignment between the radio and optical celestial reference frames.  相似文献   

18.
Gravimetric quantities are commonly represented in terms of high degree surface or solid spherical harmonics. After EGM2008, such expansions routinely extend to spherical harmonic degree 2190, which makes the computation of gravimetric quantities at a large number of arbitrarily scattered points in space using harmonic synthesis, a very computationally demanding process. We present here the development of an algorithm and its associated software for the efficient and precise evaluation of gravimetric quantities, represented in high degree solid spherical harmonics, at arbitrarily scattered points in the space exterior to the surface of the Earth. The new algorithm is based on representation of the quantities of interest in solid ellipsoidal harmonics and application of the tensor product trigonometric needlets. A FORTRAN implementation of this algorithm has been developed and extensively tested. The capabilities of the code are demonstrated using as examples the disturbing potential T, height anomaly \(\zeta \), gravity anomaly \(\Delta g\), gravity disturbance \(\delta g\), north–south deflection of the vertical \(\xi \), east–west deflection of the vertical \(\eta \), and the second radial derivative \(T_{rr}\) of the disturbing potential. After a pre-computational step that takes between 1 and 2 h per quantity, the current version of the software is capable of computing on a standard PC each of these quantities in the range from the surface of the Earth up to 544 km above that surface at speeds between 20,000 and 40,000 point evaluations per second, depending on the gravimetric quantity being evaluated, while the relative error does not exceed \(10^{-6}\) and the memory (RAM) use is 9.3 GB.  相似文献   

19.
Canadian gravimetric geoid model 2010   总被引：4，自引：1，他引：3  
A new gravimetric geoid model, Canadian Gravimetric Geoid 2010 (CGG2010), has been developed to upgrade the previous geoid model CGG2005. CGG2010 represents the separation between the reference ellipsoid of GRS80 and the Earth’s equipotential surface of $W_0=62{,}636{,}855.69~\mathrm{m}^2\mathrm{s}^{-2}$ W 0 = 62 , 636 , 855.69 m 2 s ? 2 . The Stokes–Helmert method has been re-formulated for the determination of CGG2010 by a new Stokes kernel modification. It reduces the effect of the systematic error in the Canadian terrestrial gravity data on the geoid to the level below 2 cm from about 20 cm using other existing modification techniques, and renders a smooth spectral combination of the satellite and terrestrial gravity data. The long wavelength components of CGG2010 include the GOCE contribution contained in a combined GRACE and GOCE geopotential model: GOCO01S, which ranges from $-20.1$ ? 20.1 to 16.7 cm with an RMS of 2.9 cm. Improvement has been also achieved through the refinement of geoid modelling procedure and the use of new data. (1) The downward continuation effect has been accounted accurately ranging from $-22.1$ ? 22.1 to 16.5 cm with an RMS of 0.9 cm. (2) The geoid residual from the Stokes integral is reduced to 4 cm in RMS by the use of an ultra-high degree spherical harmonic representation of global elevation model for deriving the reference Helmert field in conjunction with a derived global geopotential model. (3) The Canadian gravimetric geoid model is published for the first time with associated error estimates. In addition, CGG2010 includes the new marine gravity data, ArcGP gravity grids, and the new Canadian Digital Elevation Data (CDED) 1:50K. CGG2010 is compared to GPS-levelling data in Canada. The standard deviations are estimated to vary from 2 to 10 cm with the largest error in the mountainous areas of western Canada. We demonstrate its improvement over the previous models CGG2005 and EGM2008.  相似文献   

20.
Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.  相似文献   

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

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