首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 765 毫秒
1.
One of the main objectives of ESA’s Gravity Field and Steady-State Ocean Circulation mission GOCE (Gravity field and steady-state ocean circulation mission, 1999) is to allow global unification of height systems by directly providing potential differences between benchmarks in different height datum zones. In other words, GOCE provides a globally consistent and unbiased geoid. If this information is combined with ellipsoidal (derived from geodetic space techniques) and physical heights (derived from leveling/gravimetry) at the same benchmarks, datum offsets between the datum zones can be determined and all zones unified. The expected accuracy of GOCE is around 2–3 cm up to spherical harmonic degree n max ≈ 200. The omission error above this degree amounts to about 30 cm which cannot be neglected. Therefore, terrestrial residual gravity anomalies are necessary to evaluate the medium and short wavelengths of the geoid, i.e. one has to solve the Geodetic Boundary Value Problem (GBVP). The theory of height unification by the GBVP approach is well developed, see e.g. Colombo (A World Vertical Network. Report 296, Department of Geodetic Science and Surveying, 1980) or Rummel and Teunissen (Bull Geod 62:477–498, 1988). Thereby, it must be considered that terrestrial gravity anomalies referring to different datum zones are biased due to the respective datum offsets. Consequently, the height reference surface of a specific datum zone deviates from the unbiased geoid not only due to its own datum offset (direct bias term) but is also indirectly affected by the integration of biased gravity anomalies. The latter effect is called the indirect bias term and it considerably complicates the adjustment model for global height unification. If no satellite based gravity model is employed, this error amounts to about the same size as the datum offsets, i.e. 1–2 m globally. We show that this value decreases if a satellite-only gravity model is used. Specifically for GOCE with n max ≈ 200, the error can be expected not to exceed the level of 1 cm, allowing the effect to be neglected in practical height unification. The results are supported by recent findings by Gatti et al. (J Geod, 2012).  相似文献   

2.
We present new insights on the time-averaged surface velocities, convergence and extension rates along arc-normal transects in Kumaon, Garhwal and Kashmir–Himachal regions in the Indian Himalaya from 13 years of high-precision Global Positioning System (GPS) time series (1995–2008) derived from GPS data at 14 GPS permanent and 42 campaign stations between $29.5{-}35^{\circ }\hbox {N}$ and $76{-}81^{\circ }\hbox {E}$ . The GPS surface horizontal velocities vary significantly from the Higher to Lesser Himalaya and are of the order of 30 to 48 mm/year NE in ITRF 2005 reference frame, and 17 to 2 mm/year SW in an India fixed reference frame indicating that this region is accommodating less than 2 cm/year of the India–Eurasia plate motion ( ${\sim }4~\hbox {cm/year}$ ). The total arc-normal shortening varies between ${\sim }10{-}14~\hbox {mm/year}$ along the different transects of the northwest Himalayan wedge, between the Indo-Tsangpo suture to the north and the Indo-Gangetic foreland to the south indicating high strain accumulation in the Himalayan wedge. This convergence is being accommodated differentially along the arc-normal transects; ${\sim } 5{-}10~\hbox {mm/year}$ in Lesser Himalaya and 3–4 mm/year in Higher Himalaya south of South Tibetan Detachment. Most of the convergence in the Lesser Himalaya of Garhwal and Kumaon is being accommodated just south of the Main Central Thrust fault trace, indicating high strain accumulation in this region which is also consistent with the high seismic activity in this region. In addition, for the first time an arc-normal extension of ${\sim }6~\hbox {mm/year}$ has also been observed in the Tethyan Himalaya of Kumaon. Inverse modeling of GPS-derived surface deformation rates in Garhwal and Kumaon Himalaya using a single dislocation indicate that the Main Himalayan Thrust is locked from the surface to a depth of ${\sim }15{-}20~\hbox {km}$ over a width of 110 km with associated slip rate of ${\sim }16{-}18~\hbox {mm/year}$ . These results indicate that the arc-normal rates in the Northwest Himalaya have a complex deformation pattern involving both convergence and extension, and rigorous seismo-tectonic models in the Himalaya are necessary to account for this pattern. In addition, the results also gave an estimate of co-seismic and post-seismic motion associated with the 1999 Chamoli earthquake, which is modeled to derive the slip and geometry of the rupture plane.  相似文献   

3.
A terrestrial survey, called the Geoid Slope Validation Survey of 2011 (GSVS11), encompassing leveling, GPS, astrogeodetic deflections of the vertical (DOV) and surface gravity was performed in the United States. The general purpose of that survey was to evaluate the current accuracy of gravimetric geoid models, and also to determine the impact of introducing new airborne gravity data from the ‘Gravity for the Redefinition of the American Vertical Datum’ (GRAV-D) project. More specifically, the GSVS11 survey was performed to determine whether or not the GRAV-D airborne gravimetry, flown at 11 km altitude, can reduce differential geoid error to below 1 cm in a low, flat gravimetrically uncomplicated region. GSVS11 comprises a 325 km traverse from Austin to Rockport in Southern Texas, and includes 218 GPS stations ( $\sigma _{\Delta h }= 0.4$ cm over any distance from 0.4 to 325 km) co-located with first-order spirit leveled orthometric heights ( $\sigma _{\Delta H }= 1.3$ cm end-to-end), including new surface gravimetry, and 216 astronomically determined vertical deflections $(\sigma _{\mathrm{DOV}}= 0.1^{\prime \prime })$ . The terrestrial survey data were compared in various ways to specific geoid models, including analysis of RMS residuals between all pairs of points on the line, direct comparison of DOVs to geoid slopes, and a harmonic analysis of the differences between the terrestrial data and various geoid models. These comparisons of the terrestrial survey data with specific geoid models showed conclusively that, in this type of region (low, flat) the geoid models computed using existing terrestrial gravity, combined with digital elevation models (DEMs) and GRACE and GOCE data, differential geoid accuracy of 1 to 3 cm (1 $\sigma )$ over distances from 0.4 to 325 km were currently being achieved. However, the addition of a contemporaneous airborne gravity data set, flown at 11 km altitude, brought the estimated differential geoid accuracy down to 1 cm over nearly all distances from 0.4 to 325 km.  相似文献   

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

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

6.
We present results from a new vertical deflection (VD) traverse observed in Perth, Western Australia, which is the first of its kind in the Southern Hemisphere. A digital astrogeodetic QDaedalus instrument was deployed to measure VDs with \({\sim }\)0.2\(''\) precision at 39 benchmarks with a \({{\sim }}1~\hbox {km}\) spacing. For the conversion of VDs to quasigeoid height differences, the method of astronomical–topographical levelling was applied, based on topographical information from the Shuttle Radar Topography Mission. The astronomical quasigeoid heights are in 20–30 mm (RMS) agreement with three independent gravimetric quasigeoid models, and the astrogeodetic VDs agree to 0.2–0.3\(''\) (north–south) and 0.6–0.9\(''\) (east–west) RMS. Tilt-like biases of \({\sim }1\,\,\hbox {mm}\) over \({\sim }1\,\,\hbox {km}\) are present for all quasigeoid models within \({\sim }20\,\,\hbox {km}\) of the coastline, suggesting inconsistencies in the coastal zone gravity data. The VD campaign in Perth was designed as a low-cost effort, possibly allowing replication in other Southern Hemisphere countries (e.g., Asia, Africa, South America and Antarctica), where VD data are particularly scarce.  相似文献   

7.
Continental hydrology loading observed by VLBI measurements   总被引:1,自引:1,他引:0  
Variations in continental water storage lead to loading deformation of the crust with typical peak-to-peak variations at very long baseline interferometry (VLBI) sites of 3–15 mm in the vertical component and 1–2 mm in the horizontal component. The hydrology signal at VLBI sites has annual and semi-annual components and clear interannual variations. We have calculated the hydrology loading series using mass loading distributions derived from the global land data assimilation system (GLDAS) hydrology model and alternatively from a global grid of equal-area gravity recovery and climate experiment (GRACE) mascons. In the analysis of the two weekly VLBI 24-h R1 and R4 network sessions from 2003 to 2010 the baseline length repeatabilities are reduced in 79 % (80 %) of baselines when GLDAS (GRACE) loading corrections are applied. Site vertical coordinate repeatabilities are reduced in about 80 % of the sites when either GLDAS or GRACE loading is used. In the horizontal components, reduction occurs in 70–80 % of the sites. Estimates of the annual site vertical amplitudes were reduced for 16 out of 18 sites if either loading series was applied. We estimated loading admittance factors for each site and found that the average admittances were 1.01 \(\pm \) 0.05 for GRACE and 1.39 \(\pm \) 0.07 for GLDAS. The standard deviations of the GRACE admittances and GLDAS admittances were 0.31 and 0.68, respectively. For sites that have been observed in a set of sufficiently temporally dense daily sessions, the average correlation between VLBI vertical monthly averaged series and GLDAS or GRACE loading series was 0.47 and 0.43, respectively.  相似文献   

8.
Water vapor tomography has been developed as a powerful tool to model spatial and temporal distribution of atmospheric water vapor. Global navigation satellite systems (GNSS) water vapor tomography refers to the 3D structural construction of tropospheric water vapor using a large number of GNSS signals that penetrate the tomographic modeling area from different positions. The modeling area is usually discretized into a number of voxels. A major issue involved is that some voxels are not crossed by any GNSS signal rays, resulting in an undetermined solution to the tomographic system. To alleviate this problem, the number of voxels crossed by GNSS signal rays should be as large as possible. An important way to achieve this is to optimize the geographic distribution of tomographic voxels. We propose an approach to optimize voxel distribution in both vertical and horizontal domains. In the vertical domain, water vapor profiles derived from radiosonde data are exploited to identify the maximum height of tomography and the optimal vertical resolution. In the horizontal domain, the optimal horizontal distribution of voxels is obtained by searching the maximum number of ray-crossing voxels in both latitude and longitude directions. The water vapor tomography optimization procedures are implemented using GPS water vapor data from the Hong Kong Satellite Positioning Reference Station Network. The tomographic water vapor fields solved from the optimized tomographic voxels are evaluated using radiosonde data and a numerical weather prediction non-hydrostatic model (NHM) obtained for the Hong Kong station. The comparisons of tomographic integrated water vapor (IWV) with the radiosonde and NHM IWV show that RMS errors of their differences are 1.41 and 3.09 mm, respectively. Moreover, the tomographic water vapor density results are compared with those of radiosonde and NHM. The RMS error of the density differences between tomography and radiosonde data is 1.05  \(\mathrm{g/m}^{3}\) . For the comparison between tomography and NHM, an overall RMS error of \(1.43\,\mathrm{g/m^{3}}\) is achieved.  相似文献   

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

10.
It is difficult to calculate the accurate ground movement due to deep underground mining because of the complexity of the geotechnical environment. Guan-Zhuang iron mine is a pillarless sublevel caving mine operated by Luzhong Metallurgical Mining Company, south-east of Jinan, PR China. It mines the Zhangjiawa Seam at a depth of approximately 520 m. Although the towers are outside the conventional ‘angle of draw’ subsidence influence criteria, and have seen only negligible vertical displacement as a result of deep mining, there has been widespread evidence of regional horizontal displacement of the land surface, large distances away from the mining area. Possible explanations of these displacements include one or a combination of mechanisms such as pre-mining stress relaxation, regional joint patterns, soft rock strata, displacement toward active goaf areas. Luzhong Metallurgical Mining Company have been making precise measurements of distances near the shaft towers in the Guan-Zhuang iron mine since 2003. The results show horizontal displacements of up to 96 mm occur even when underground mining is about 0.8 km from the survey displacements. From an analysis of these and other survey results it is concluded that mining effects extend a long way from deep mining. The results also show that ground horizontal displacements are typically at least as great as the vertical component, that the maximum horizontal displacement occurs soon after undermining.  相似文献   

11.
Well credited and widely used ionospheric models, such as the International Reference Ionosphere or NeQuick, describe the variation of the electron density with height by means of a piecewise profile tied to the F2-peak parameters: the electron density, $N_m \mathrm{F2}$ N m F 2 , and the height, $h_m \mathrm{F2}$ h m F 2 . Accurate values of these parameters are crucial for retrieving reliable electron density estimations from those models. When direct measurements of these parameters are not available, the models compute the parameters using the so-called ITU-R database, which was established in the early 1960s. This paper presents a technique aimed at routinely updating the ITU-R database using radio occultation electron density profiles derived from GPS measurements gathered from low Earth orbit satellites. Before being used, these radio occultation profiles are validated by fitting to them an electron density model. A re-weighted Least Squares algorithm is used for down-weighting unreliable measurements (occasionally, entire profiles) and to retrieve $N_m \mathrm{F2}$ N m F 2 and $h_m \mathrm{F2}$ h m F 2 values—together with their error estimates—from the profiles. These values are used to monthly update the database, which consists of two sets of ITU-R-like coefficients that could easily be implemented in the IRI or NeQuick models. The technique was tested with radio occultation electron density profiles that are delivered to the community by the COSMIC/FORMOSAT-3 mission team. Tests were performed for solstices and equinoxes seasons in high and low-solar activity conditions. The global mean error of the resulting maps—estimated by the Least Squares technique—is between $0.5\times 10^{10}$ 0.5 × 10 10 and $3.6\times 10^{10}$ 3.6 × 10 10 elec/m $^{-3}$ ? 3 for the F2-peak electron density (which is equivalent to 7 % of the value of the estimated parameter) and from 2.0 to 5.6 km for the height ( $\sim $ 2 %).  相似文献   

12.
Error analysis of the NGS’ surface gravity database   总被引:1,自引:1,他引:0  
Are the National Geodetic Survey’s surface gravity data sufficient for supporting the computation of a 1 cm-accurate geoid? This paper attempts to answer this question by deriving a few measures of accuracy for this data and estimating their effects on the US geoid. We use a data set which comprises ${\sim }1.4$ million gravity observations collected in 1,489 surveys. Comparisons to GRACE-derived gravity and geoid are made to estimate the long-wavelength errors. Crossover analysis and $K$ -nearest neighbor predictions are used for estimating local gravity biases and high-frequency gravity errors, and the corresponding geoid biases and high-frequency geoid errors are evaluated. Results indicate that 244 of all 1,489 surface gravity surveys have significant biases ${>}2$  mGal, with geoid implications that reach 20 cm. Some of the biased surveys are large enough in horizontal extent to be reliably corrected by satellite-derived gravity models, but many others are not. In addition, the results suggest that the data are contaminated by high-frequency errors with an RMS of ${\sim }2.2$  mGal. This causes high-frequency geoid errors of a few centimeters in and to the west of the Rocky Mountains and in the Appalachians and a few millimeters or less everywhere else. Finally, long-wavelength ( ${>}3^{\circ }$ ) surface gravity errors on the sub-mGal level but with large horizontal extent are found. All of the south and southeast of the USA is biased by +0.3 to +0.8 mGal and the Rocky Mountains by $-0.1$ to $-0.3$  mGal. These small but extensive gravity errors lead to long-wavelength geoid errors that reach 60 cm in the interior of the USA.  相似文献   

13.
We show that the current levels of accuracy being achieved for the precise orbit determination (POD) of low-Earth orbiters demonstrate the need for the self-consistent treatment of tidal variations in the geocenter. Our study uses as an example the POD of the OSTM/Jason-2 satellite altimeter mission based upon Global Positioning System (GPS) tracking data. Current GPS-based POD solutions are demonstrating root-mean-square (RMS) radial orbit accuracy and precision of \({<}1\)  cm and 1 mm, respectively. Meanwhile, we show that the RMS of three-dimensional tidal geocenter variations is \({<}6\)  mm, but can be as large as 15 mm, with the largest component along the Earth’s spin axis. Our results demonstrate that GPS-based POD of Earth orbiters is best performed using GPS satellite orbit positions that are defined in a reference frame whose origin is at the center of mass of the entire Earth system, including the ocean tides. Errors in the GPS-based POD solutions for OSTM/Jason-2 of \({<}4\)  mm (3D RMS) and \({<}2\)  mm (radial RMS) are introduced when tidal geocenter variations are not treated consistently. Nevertheless, inconsistent treatment is measurable in the OSTM/Jason-2 POD solutions and manifests through degraded post-fit tracking data residuals, orbit precision, and relative orbit accuracy. For the latter metric, sea surface height crossover variance is higher by \(6~\hbox {mm}^{2}\) when tidal geocenter variations are treated inconsistently.  相似文献   

14.
Increasing population and natural disasters like drought, flood, cyclone etc., has impacted global agriculture area and hence continuously modifying cropping pattern and associated statistics. The present study analysed agriculture dynamics over one of the densely populated and disaster prone state (Bihar) in India and derived vital statistics (single, double and triple cropping area, and monthly, seasonal, annual and long term status at the state and district level) for the years 2001–2012. The study used time-series MODIS vegetation index (EVI; MOD13A2, 1 km, 16 day, 2001–2012), MODIS annual Land Cover product (MCD12Q1, 500 m, 2001–2012) and Global Land Cover map (Scasia_V4, 1 km, 2000; Globcover_V2.2, 300 m, 2005/2006 and V2.3, 2009, 300 m), and extracted horizontal (i.e., area change) and vertical (i.e., cropping intensification) agriculture change pattern. The results were inter-compared, and validated using government reports as well as with high spatial resolution data (IRS-LISS III 23.5 m). From 2001–2006 to 2007–2012, the net horizontal and vertical change in agriculture area is +145.24 and +907.82 km2, respectively, and net change in seasonal crop area (winter, summer and monsoon) is +959.21, +1009.84 and ?1061.64 km2, respectively. The districts which are located along the eastern part of Ganges experienced maximum positive changes and the districts along Gandak river in the north-western part of the study area experienced maximum negative changes. Overall, the study has quantified and revealed interesting space–time agriculture change patterns over 12 years including impacts caused by droughts and floods in the study area.  相似文献   

15.
The proper identification and removal of outliers in the combination of rates of vertical displacements derived from GPS, tide gauges/satellite altimetry, and GRACE observations is presented. Outlier detection is a necessary pre-screening procedure in order to ensure reliable estimates of stochastic properties of the observations in the combined least-squares adjustment (via rescaling of covariance matrices) and to ensure that the final vertical motion model is not corrupted and/or distorted by erroneous data. Results from this study indicate that typical data snooping methods are inadequate in dealing with these heterogeneous data sets and their stochastic properties. Using simulated vertical displacement rates, it is demonstrated that a large variety of outliers (random scattered and adjacent, as well as jointly influential) can be dealt with if an iterative re-weighting least-squares adjustment is combined with a robust median estimator. Moreover, robust estimators are efficient in areas weakly constrained by the data, where even high quality observations may appear to be erroneous if their estimates are largely influenced by outliers. Four combined models for the vertical motion in the region of the Great Lakes are presented. The computed vertical displacements vary between  − 2 mm/year (subsidence) along the southern shores and 3 mm/year (uplift) along the northern shores. The derived models provide reliable empirical constraints and error bounds for postglacial rebound models in the region.  相似文献   

16.
A radial integration of spherical mass elements (i.e. tesseroids) is presented for evaluating the six components of the second-order gravity gradient (i.e. second derivatives of the Newtonian mass integral for the gravitational potential) created by an uneven spherical topography consisting of juxtaposed vertical prisms. The method uses Legendre polynomial series and takes elastic compensation of the topography by the Earth’s surface into account. The speed of computation of the polynomial series increases logically with the observing altitude from the source of anomaly. Such a forward modelling can be easily applied for reduction of observed gravity gradient anomalies by the effects of any spherical interface of density. An iterative least-squares inversion of measured gravity gradient coefficients is also proposed to estimate a regional set of juxtaposed topographic heights. Several tests of recovery have been made by considering simulated gradients created by idealistic conical and irregular Great Meteor seamount topographies, and for varying satellite altitudes and testing different levels of uncertainty. In the case of gravity gradients measured at a GOCE-type altitude of \(\sim \)300 km, the search converges down to a stable but smooth topography after 10–15 iterations, while the final root-mean-square error is \(\sim \)100 m that represents only 2 % of the seamount amplitude. This recovery error decreases with the altitude of the gravity gradient observations by revealing more topographic details in the region of survey.  相似文献   

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

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

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

20.
重力卫星GRACE(gravity recovery and climate experiment)监测斯堪的纳维亚半岛陆地水储量变化会受到冰川均衡调整(GIA)信号的严重影响。首先根据该地区绝对重力和GPS并址观测数据计算了GIA重力和垂直位移的实测线性比值,利用该比值和GPS网观测的垂直位移速度场得到了GIA重力。然后,对GRACE观测的重力变化速率进行GIA重力改正,进而可分离陆地水储量变化趋势,避免了使用GIA模型所带来的巨大不确定性,并根据观测数据完整估计了所得结果的不确定性。最后与水文模型作对比分析。结果表明,实测的GIA重力-垂直位移线性比值为0.148±0.020μGal/mm(1Gal=10-2 m/s2),该结果检验了Wahr的理论近似值且与北美实测的结果非常接近。2003年1月至2011年3月期间,斯堪的纳维亚半岛陆地水储量存在明显的增加趋势,信号的主体位于半岛南端的维纳恩湖附近,总的水量增加速率为4.6±2.1km3/a,数据观测期间的累积增加水量为38±17km3。研究结果与WGHM水文模型的结果有较好的一致性,相关系数达到0.69,而与GLDAS水文模型的相关性略小。  相似文献   

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

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