共查询到20条相似文献，搜索用时 468 毫秒
1.
With the advent of geodetic satellite missions mapping almost globally the Earth’s gravitational field, new methods and theoretical approaches have been developed and investigated to fully exploit the potential of their new observables. Besides estimating values of numerical coefficients in harmonic series of the gravitational potential, new applications emerged such as data validation and combination. In this contribution, new integral transformations are presented which transform principal components of the terrestrial deflection of the vertical onto disturbing satellite-to-satellite tracking and gradiometric data at altitude. Using spherical approximation, necessary integral kernel functions are derived in both spectral and closed forms. The behaviour of isotropic kernel functions is studied and the new integral transformations are tested in a closed-loop simulation using synthetic terrestrial and satellite data synthesized from a global gravitational model. New integral transformations can be used for data validation and combination purposes. 相似文献
2.
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. 相似文献
3.
New integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients are derived in this article. They provide more options for continuation of gravitational gradient combinations and extend available mathematical apparatus formulated for this purpose up to now. The starting point represents the analytical solution of the spherical gradiometric boundary value problem in the spatial domain. Applying corresponding differential operators on the analytical solution of the spherical gradiometric boundary value problem, a total of 18 integral formulas are provided. Spatial and spectral forms of isotropic kernels are given and their behaviour for parameters of a GOCE-like satellite is investigated. Correctness of the new integral formulas and the isotropic kernels is tested in a closed-loop simulation. The derived integral formulas and the isotropic kernels form a theoretical basis for validation purposes and geophysical applications of satellite gradiometric data as provided currently by the GOCE mission. They also extend the well-known Meissl scheme. 相似文献
4.
Satellite gradiometry using a satellite pair 总被引：1，自引：1，他引：1
The GRACE mission has substantiated the low–low satellite-to-satellite tracking (LL-SST) concept. The LL-SST configuration can be combined with the previously realized high–low SST concept in the CHAMP mission to provide a much higher accuracy. The line of sight (LOS) acceleration difference between the GRACE satellite pair, the simplest form of the combined observable, is mostly used for mapping the global gravity field of the Earth in terms of spherical harmonic coefficients. As an alternative observable, a linear combination of the gravitational gradient tensor components is proposed. Being a one-point function and having a direct relation with the field geometry (curvature of the field at the point) are two noteworthy achievements of the alternative formulation. In addition, using an observation quantity that is related to the second-instead of the first-order derivatives of the gravitational potential amplifies the high-frequency part of the signal. Since the transition from the first- to the second-order derivatives includes the application of a finite-differences scheme, the high-frequency part of the noise is also amplified. Nevertheless, due to the different spectral behaviour of signal and noise, in the end the second-order approach leads to improved gravitational field resolution. Mathematical formulae for the gradiometry approach, for both linear and higher-degree approximations, are derived. The proposed approach is implemented for recovery of the global gravitational field and the results are compared with those of LOS acceleration differences. Moreover, LOS acceleration difference residuals are calculated, which are at the level of a few tenths of mGal. Error analysis shows that the residuals of the estimated degree variances are less than 10–3. Furthermore, the gravity anomaly residuals are less than 2 mGal for most points on the Earth. 相似文献
5.
Vertical and horizontal spheroidal boundary-value problems (BVPs), i.e., determination of the external gravitational potential from the components of the gravitational gradient on the spheroid, are discussed in this article. The gravitational gradient is decomposed into the series of the vertical and horizontal vector spheroidal harmonics, before being orthogonalized in a weighted sense by two different approaches. The vertical and horizontal spheroidal BVPs are then formulated and solved in the spectral and spatial domains. Both orthogonalization methods provide the same analytical solutions for the vertical spheroidal BVP, and give distinct, but equivalent, analytical solutions for the horizontal spheroidal BVP. A closed-loop simulation is performed to test the correctness of the analytical solutions, and we investigate analytical properties of the sub-integral kernels. The systematic treatment of the spheroidal BVPs and the resulting mathematical equations extend the theoretical apparatus of geodesy and of the potential theory. 相似文献
6.
为探究重力场恢复与气候实验(gravity recovery and climate experiment,GRACE)卫星与全球定位系统(global positioning system,GPS)两种独立技术获取的因陆地水储量变化引起的地壳垂向季节性位移的一致性,选取澳大利亚27个GPS站点5～10 a的高程时间序... 相似文献
7.
基于中国大陆构造环境监测网络的连续GPS观测数据,比较分析了中国大陆234个GPS台站和GRACE得到的地表垂直形变。GPS和GRACE垂直形变具有较好的一致性,反映了地表质量变化是引起GPS垂直形变非线性变化的重要因素之一,但二者也存在着一定的差异。为定量分析GPS和GRACE垂直形变的差异,探讨了热膨胀效应对GPS垂直位移的影响及区域地壳结构对GRACE估算地表垂直负荷形变的影响。结果表明,中国大陆50%以上的GPS台站热膨胀垂直形变周年振幅不小于1 mm;对GPS进行热膨胀效应改正后,中国大陆GPS与GRACE垂直形变具有更好的一致性;GPS与GRACE垂直形变周年振幅比值由1.07±0.06变为1.01±0.05;热膨胀效应可以解释6.2%的GPS与GRACE垂直形变的差异,热膨胀效应改正可使GPS和GRACE垂直形变的一致性相对增加11.2%。是否顾及区域地壳结构引起的GRACE估算中国大陆垂直负荷形变的相对差异为2.5%。 相似文献
8.
Local geoid determination from airborne vector gravimetry 总被引：3，自引：2，他引：1
Methods are illustrated to compute the local geoid using the vertical and horizontal components of the gravity disturbance vector derived from an airborne GPS/inertial navigation system. The data were collected by the University of Calgary in a test area of the Canadian Rocky Mountains and consist of multiple parallel tracks and two crossing tracks of accelerometer and gyro measurements, as well as precise GPS positions. Both the boundary-value problem approach (Hotines integral) and the profiling approach (line integral) were applied to compute the disturbing potential at flight altitude. Cross-over adjustments with minimal control were investigated and utilized to remove error biases and trends in the estimated gravity disturbance components. Final estimation of the geoid from the vertical gravity disturbance included downward continuation of the disturbing potential with correction for intervening terrain masses. A comparison of geoid estimates to the Canadian Geoid 2000 (CGG2000) yielded an average standard deviation per track of 14 cm if they were derived from the vertical gravity disturbance (minimally controlled with a cross-over adjustment), and 10 cm if derived from the horizontal components (minimally controlled in part with a simulated cross-over adjustment). Downward continuation improved the estimates slightly by decreasing the average standard deviation by about 0.5 cm. The application of a wave correlation filter to both types of geoid estimates yielded significant improvement by decreasing the average standard deviation per track to 7.6 cm. 相似文献
9.
New spherical integral formulas between components of the second- and third-order gravitational tensors are formulated in this article. First, we review the nomenclature and basic properties of the second- and third-order gravitational tensors. Initial points of mathematical derivations, i.e., the second- and third-order differential operators defined in the spherical local North-oriented reference frame and the analytical solutions of the gradiometric boundary-value problem, are also summarized. Secondly, we apply the third-order differential operators to the analytical solutions of the gradiometric boundary-value problem which gives 30 new integral formulas transforming (1) vertical-vertical, (2) vertical-horizontal and (3) horizontal-horizontal second-order gravitational tensor components onto their third-order counterparts. Using spherical polar coordinates related sub-integral kernels can efficiently be decomposed into azimuthal and isotropic parts. Both spectral and closed forms of the isotropic kernels are provided and their limits are investigated. Thirdly, numerical experiments are performed to test the consistency of the new integral transforms and to investigate properties of the sub-integral kernels. The new mathematical apparatus is valid for any harmonic potential field and may be exploited, e.g., when gravitational/magnetic second- and third-order tensor components become available in the future. The new integral formulas also extend the well-known Meissl diagram and enrich the theoretical apparatus of geodesy. 相似文献
10.
There are two spurious jumps in the atmospheric part of the Gravity Recovery and Climate Experiment-Atmosphere and Ocean De-aliasing level 1B (GRACE-AOD1B) products, which occurred in January-February of the years 2006 and 2010, as a result of the vertical level and horizontal resolution changes in the ECMWFop (European Centre for Medium-Range Weather Forecasts operational analysis). These jumps cause a systematic error in the estimation of mass changes from GRACE time-variable level 2 products, since GRACE-AOD1B mass variations are removed during the computation of GRACE level 2. In this short note, the potential impact of using an improved set of 6-hourly atmospheric de-aliasing products on the computations of linear trends as well as the amplitude of annual and semi-annual mass changes from GRACE is assessed. These improvements result from 1) employing a modified 3D integration approach (ITG3D), and 2) using long-term consistent atmospheric fields from the ECMWF reanalysis (ERA-Interim). The monthly averages of the new ITG3D-ERA-Interim de-aliasing products are then compared to the atmospheric part of GRACE-AOD1B, covering January 2003 to December 2010. These comparisons include the 33 world largest river basins along with Greenland and Antarctica ice sheets. The results indicate a considerable difference in total atmospheric mass derived from the two products over some of the mentioned regions. We suggest that future GRACE studies consider these through updating uncertainty budgets or by applying corrections to estimated trends and amplitudes/phases. 相似文献
11.
Gravitational gradients by tensor analysis with application to spherical coordinates 总被引：1，自引：1，他引：0
This contribution deals with the derivation of explicit expressions of the gradients of first, second and third order of the
gravitational potential. This is accomplished in the framework of tensor analysis which naturally allows to apply general
formulae to the specific coordinate systems in use in geodesy. In particular it is recalled here that when the potential field
is expressed in general coordinates on a 3D manifold, the gradient operation leads to the definition of the covariant derivative
and that the covariant derivative of a tensor can be obtained by application of a simple rule. When applied to the gravitational
potential or to any of its gradients, the rule straightforwardly provides the expressions of the higher-order gradients. It
is also shown that the tensor approach offers a clear distinction between natural and physical components of the gradients.
Two fundamental reference systems—a global, bodycentric system and a local, topocentric system, both body-fixed—are introduced
and transformation rules are derived to convert quantities between the two systems. The results include explicit expressions
for the gradients of the first three orders in both reference systems. 相似文献
12.
E. W. Grafarend 《Journal of Geodesy》2001,75(7-8):363-390
In a comparison of the solution of the spherical horizontal and vertical boundary value problems of physical geodesy it is
aimed to construct downward continuation operators for vertical deflections (surface gradient of the incremental gravitational
potential) and for gravity disturbances (vertical derivative of the incremental gravitational potential) from points on the
Earth's topographic surface or of the three-dimensional (3-D) Euclidean space nearby down to the international reference sphere
(IRS). First the horizontal and vertical components of the gravity vector, namely spherical vertical deflections and spherical
gravity disturbances, are set up. Second, the horizontal and vertical boundary value problem in spherical gravity and geometry
space is considered. The incremental gravity vector is represented in terms of vector spherical harmonics. The solution of
horizontal spherical boundary problem in terms of the horizontal vector-valued Green function converts vertical deflections
given on the IRS to the incremental gravitational potential external in the 3-D Euclidean space. The horizontal Green functions
specialized to evaluation and source points on the IRS coincide with the Stokes kernel for vertical deflections. Third, the
vertical spherical boundary value problem is solved in terms of the vertical scalar-valued Green function. Fourth, the operators
for upward continuation of vertical deflections given on the IRS to vertical deflections in its external 3-D Euclidean space
are constructed. Fifth, the operators for upward continuation of incremental gravity given on the IRS to incremental gravity
to the external 3-D Euclidean space are generated. Finally, Meissl-type diagrams for upward continuation and regularized downward
continuation of horizontal and vertical gravity data, namely vertical deflection and incremental gravity, are produced.
Received: 10 May 2000 / Accepted: 26 February 2001 相似文献
13.
Guillaume Ramillien R. Biancale S. Gratton X. Vasseur S. Bourgogne 《Journal of Geodesy》2011,85(6):313-328
We propose an unconstrained approach to recover regional time-variations of surface mass anomalies using Level-1 Gravity Recovery
and Climate Experiment (GRACE) orbit observations, for reaching spatial resolutions of a few hundreds of kilometers. Potential
differences between the twin GRACE vehicles are determined along short satellite tracks using the energy integral method (i.e.,
integration of orbit parameters vs. time) in a quasi-inertial terrestrial reference frame. Potential differences residuals
corresponding mainly to changes in continental hydrology are then obtained after removing the gravitational effects of the
known geophysical phenomena that are mainly the static part of the Earth’s gravity field and time-varying contributions to
gravity (Sun, Moon, planets, atmosphere, ocean, tides, variations of Earth’s rotation axis) through ad hoc models. Regional
surface mass anomalies are restored from potential difference anomalies of 10 to 30-day orbits onto 1◦ continental grids by regularization techniques based on singular value decomposition. Error budget analysis has been made
by considering the important effects of spectrum truncation, the time length of observation (or spatial coverage of the data
to invert) and for different levels of noise. 相似文献
14.
Cuixian Lu Xingxing Li Florian Zus Robert Heinkelmann Galina Dick Maorong Ge Jens Wickert Harald Schuh 《Journal of Geodesy》2017,91(9):1019-1029
Precise positioning with the current Chinese BeiDou Navigation Satellite System is proven to be of comparable accuracy to the Global Positioning System, which is at centimeter level for the horizontal components and sub-decimeter level for the vertical component. But the BeiDou precise point positioning (PPP) shows its limitation in requiring a relatively long convergence time. In this study, we develop a numerical weather model (NWM) augmented PPP processing algorithm to improve BeiDou precise positioning. Tropospheric delay parameters, i.e., zenith delays, mapping functions, and horizontal delay gradients, derived from short-range forecasts from the Global Forecast System of the National Centers for Environmental Prediction (NCEP) are applied into BeiDou real-time PPP. Observational data from stations that are capable of tracking the BeiDou constellation from the International GNSS Service (IGS) Multi-GNSS Experiments network are processed, with the introduced NWM-augmented PPP and the standard PPP processing. The accuracy of tropospheric delays derived from NCEP is assessed against with the IGS final tropospheric delay products. The positioning results show that an improvement in convergence time up to 60.0 and 66.7% for the east and vertical components, respectively, can be achieved with the NWM-augmented PPP solution compared to the standard PPP solutions, while only slight improvement in the solution convergence can be found for the north component. A positioning accuracy of 5.7 and 5.9 cm for the east component is achieved with the standard PPP that estimates gradients and the one that estimates no gradients, respectively, in comparison to 3.5 cm of the NWM-augmented PPP, showing an improvement of 38.6 and 40.1%. Compared to the accuracy of 3.7 and 4.1 cm for the north component derived from the two standard PPP solutions, the one of the NWM-augmented PPP solution is improved to 2.0 cm, by about 45.9 and 51.2%. The positioning accuracy for the up component improves from 11.4 and 13.2 cm with the two standard PPP solutions to 8.0 cm with the NWM-augmented PPP solution, an improvement of 29.8 and 39.4%, respectively. 相似文献
15.
The static gravity field model DGM-1S from GRACE and GOCE data: computation, validation and an analysis of GOCE mission’s added value 总被引：4，自引：2，他引：2
We present a global static model of the Earth’s gravity field entitled DGM-1S based on GRACE and GOCE data. The collection of used data sets includes nearly 7 years of GRACE KBR data and 10 months of GOCE gravity gradient data. The KBR data are transformed with a 3-point differentiation into quantities that are approximately inter-satellite accelerations. Gravity gradients are processed in the instrumental frame. Noise is handled with a frequency-dependent data weighting. DGM-1S is complete to spherical harmonic degree 250 with a Kaula regularization being applied above degree 179. Its performance is compared with a number of other satellite-only GRACE/GOCE models by confronting them with (i) an independent model of the oceanic mean dynamic topography, and (ii) independent KBR and gravity gradient data. The tests reveal a competitive quality for DGM-1S. Importantly, we study added value of GOCE data by comparing the performance of satellite-only GRACE/GOCE models with models produced without GOCE data: either ITG-Grace2010s or EGM2008 depending on which of the two performs better in a given region. The test executed based on independent gravity gradients quantifies this added value as 25–38 % in the continental areas poorly covered with terrestrial gravimetry data (Equatorial Africa, Himalayas, and South America), 7–17 % in those with a good coverage with these data (Australia, North America, and North Eurasia), and 14 % in the oceans. This added value is shown to be almost entirely related to coefficients below degree 200. It is shown that this gain must be entirely attributed to gravity gradients acquired by the mission. The test executed based on an independent model of the mean dynamic topography suggests that problems still seem to exist in satellite-only GRACE/GOCE models over the Pacific ocean, where noticeable deviations between these models and EGM2008 are detected, too. 相似文献
16.
利用GOCE卫星轨道反演地球重力场模型 总被引：1，自引：1，他引：0
根据积分方程法反演地球重力场的数学模型,利用GOCE卫星2009-11-02～2010-01-02共61d的精密轨道数据反演了几组地球重力场模型。结果表明,GOCE卫星轨道能有效提取地球重力场的长波信息,弥补了GOCE卫星重力梯度带宽的限制,在106阶次的大地水准面误差为±9.6cm,该阶次精度优于EIGEN-CHAMP03S及GRACE卫星两个月轨道反演地球重力场的精度,但由于两极空白,反演的带谐位系数精度偏低。联合GOCE及GRACE卫星轨道反演的模型在106阶次的大地水准面误差为±6.9cm,弥补了GOCE卫星轨道的缺陷。 相似文献
17.
Pavel Ditmar Jo?o Teixeira da Encarna??o Hassan Hashemi Farahani 《Journal of Geodesy》2012,86(6):441-465
Spectral analysis of data noise is performed in the context of gravity field recovery from inter-satellite ranging measurements acquired by the satellite gravimetry mission GRACE. The motivation of the study is two-fold: (i) to promote a further improvement of GRACE data processing techniques and (ii) to assist designing GRACE follow-on missions. The analyzed noise realizations are produced as the difference between the actual GRACE inter-satellite range measurements and the predictions based on state-of-the-art force models. The exploited functional model is based on the so-called “range combinations,” which can be understood as a finite-difference analog of inter-satellite accelerations projected onto the line-of-sight connecting the satellites. It is shown that low-frequency noise is caused by limited accuracy of the computed GRACE orbits. In the first instance, it leads to an inaccurate estimation of the radial component of the inter-satellite velocities. A large impact of this component stems from the fact that it is directly related to centrifugal accelerations, which have to be taken into account when the measured range-accelerations are linked with inter-satellite accelerations. Another effect of orbit inaccuracies is a miscalculation of forces acting on the satellites (particularly, the one described by the zero-degree term of the Earth’s gravitational field). The major contributors to the noise budget at high frequencies (above 9?mHz) are (i) ranging sensor errors and (ii) limited knowledge of the Earth’s static gravity field at high degrees. Importantly, we show that updating the model of the static field on the basis of the available data must be performed with a caution as the result may not be physical due to a non-unique recovery of high-degree coefficients. The source of noise in the range of intermediate frequencies (1–9?mHz), which is particularly critical for an accurate gravity field recovery, is not fully understood yet. We show, however, that it cannot be explained by inaccuracies in background models of time-varying gravity field. It is stressed that most of the obtained results can be treated as sufficiently general (i.e., applicable in the context of a statistically optimal estimation based on any functional model). 相似文献
18.
联合Argo浮标、卫星测高和GRACE数据研究海平面变化 总被引：1，自引：1，他引：0
卫星测高、GRACE、Argo等数据为监测海平面变化提供了丰富的观测数据,利用Argo数据计算的比容海平面变化,可以更加深入地理解卫星测高以及卫星重力获得的海平面变化。利用2004年1月至2010年12月间Argo浮标采集的温度和盐度数据,通过数值积分方法计算了65°S~65°N间的比容海平面异常,并通过最小二乘拟合得到比容海平面变化的长期趋势为0.63±0.45 mm/a,与Llovel得到的结果吻合较好。利用卫星测高数据得到该时间段内海平面变化趋势为2.52±0.71 mm/a,GRACE反演得到的海水质量变化引起的海平面趋势为1.84±0.13mm/a,结果表明海水质量变化成为引起海平面变化的主要因素。最后对联合卫星测高、GRACE得到比容海平面变化与相应Argo浮标数据计算结果的空间分布特征进行了比较。 相似文献
19.
Separation of global time-variable gravity signals into maximally independent components 总被引：6，自引：3，他引：3
The Gravity Recovery and Climate Experiment (GRACE) products provide valuable information about total water storage variations over the whole globe. Since GRACE detects mass variations integrated over vertical columns, it is desirable to separate its total water storage anomalies into their original sources. Among the statistical approaches, the principal component analysis (PCA) method and its extensions have been frequently proposed to decompose the GRACE products into space and time components. However, these methods only search for decorrelated components that on the one hand are not always interpretable and on the other hand often contain a superposition of independent source signals. In contrast, independent component analysis (ICA) represents a technique that separates components based on assumed statistical independence using higher-order statistical information. If one assumes that independent physical processes generate statistically independent signal components added up in the GRACE observations, separating them by ICA is a reliable strategy to identify these processes. In this paper, the performance of the conventional PCA, its rotated extension and ICA are investigated when applied to the GRACE-derived total water storage variations. These analyses have been tested on both a synthetic example and on the real GRACE level-2 monthly solutions derived from GeoForschungsZentrum Potsdam (GFZ RL04) and Bonn University (ITG2010). Within the synthetic example, we can show how imposing statistical independence in the framework of ICA improves the extraction of the ‘original’ signals from a GRACE-type super-position. We are therefore confident that also for the real case the ICA algorithm, without making prior assumptions about the long-term behaviour or on the frequencies contained in the signal, improves over the performance of PCA and its rotated extension in the separation of periodical and long-term components. 相似文献
20.
Topographic–isostatic masses represent an important source of gravity field information, especially in the high-frequency
band, even if the detailed mass-density distribution inside the topographic masses is unknown. If this information is used
within a remove-restore procedure, then the instability problems in downward continuation of gravity observations from aircraft
or satellite altitudes can be reduced. In this article, integral formulae are derived for determination of gravitational effects
of topographic–isostatic masses on the first- and second-order derivatives of the gravitational potential for three topographic–isostatic
models. The application of these formulas is useful for airborne gravimetry/gradiometry and satellite gravity gradiometry.
The formulas are presented in spherical approximation by separating the 3D integration in an analytical integration in the
radial direction and 2D integration over the mean sphere. Therefore, spherical volume elements can be considered as being
approximated by mass-lines located at the centre of the discretization compartments (the mass of the tesseroid is condensed
mathematically along its vertical axis). The errors of this approximation are investigated for the second-order derivatives
of the topographic–isostatic gravitational potential in the vicinity of the Earth’s surface. The formulas are then applied
to various scenarios of airborne gravimetry/gradiometry and satellite gradiometry. The components of the gravitational vector
at aircraft altitudes of 4 and 10 km have been determined, as well as the gravitational tensor components at a satellite altitude
of 250 km envisaged for the forthcoming GOCE (gravity field and steady-state ocean-circulation explorer) mission. The numerical
computations are based on digital elevation models with a 5-arc-minute resolution for satellite gravity gradiometry and 1-arc-minute
resolution for airborne gravity/gradiometry. 相似文献