GRACE和GOCE卫星监测的青藏高原重力场特征

文章阐述了对青藏高原重力场进行研究的意义,并进一步利用重力卫星GRACE和GOCE的数据对该区域的重力场特征进行了描述.通过对该区域的重力异常、径向引力梯度的计算和分析,可以得出:在青藏高原的西部,有明显的3条重力异常区,这与当地的地形有关,也与断层的位置有关;引力梯度比重力异常具有更高的空间分辨率;重力变化剧烈的区域与梯度的异常区有一定的对应关系,同时也是地球动力活动变化剧烈的区域.     

Approximating covariance matrices estimated in multivariate models by estimated auto- and cross-covariances

Quantities like tropospheric zenith delays or station coordinates are repeatedly measured at permanent VLBI or GPS stations so that time series for the quantities at each station are obtained. The covariances of these quantities can be estimated in a multivariate linear model. The covariances are needed for computing uncertainties of results derived from these quantities. The covariance matrix for many permanent stations becomes large, the need for simplifying it may therefore arise under the condition that the uncertainties of derived results still agree. This is accomplished by assuming that the different time series of a quantity like the station height for each permanent station can be combined to obtain one time series. The covariance matrix then follows from the estimates of the auto- and cross-covariance functions of the combined time series. A further approximation is found, if compactly supported covariance functions are fitted to an estimated autocovariance function in order to obtain a covariance matrix which is representative of different kinds of measurements. The simplification of a covariance matrix estimated in a multivariate model is investigated here for the coordinates of points of a grid measured repeatedly by a laserscanner. The approximations are checked by determining the uncertainty of the sum of distances to the points of the grid. To obtain a realistic value for this uncertainty, the covariances of the measured coordinates have to be considered. Three different setups of measurements are analyzed and a covariance matrix is found which is representative for all three setups. Covariance matrices for the measurements of laserscanners can therefore be determined in advance without estimating them for each application.     

GRACE和GOCE卫星重力场模型的谱组合研究

首先,介绍了基于不同卫星重力场模型的位系数误差方差谱组合、误差阶方差谱组合模型的算法;其次,根据位系数误差方差定权和位系数误差阶方差定权两种谱组合方法编写计算程序;最后,采用GRACE和GOCE模型进行谱组合计算,并对谱组合计算模型进行内、外符合精度分析与验证,谱组合得到的重力场模型精度和可靠性优于单一重力场模型,验证了谱组合方法的有效性.     

Variations in the accuracy of gravity recovery due to ground track variability: GRACE,CHAMP, and GOCE

Following an earlier recognition of degraded monthly geopotential recovery from GRACE (Gravity Recovery And Climate Experiment) due to prolonged passage through a short repeat (low order resonant) orbit, we extend these insights also to CHAMP (CHAllenging Minisatellite Payload) and GOCE (Gravity field and steady state Ocean Circulation Explorer). We show wide track-density variations over time for these orbits in both latitude and longitude, and estimate that geopotential recovery will be as widely affected as well within all these regimes, with lesser track density leading to poorer recoveries. We then use recent models of atmospheric density to estimate the future orbit of GRACE and warn of degraded performance as other low order resonances are encountered in GRACE’s free fall. Finally implications for the GOCE orbit are discussed.     

对误差传播定律一些应用问题的再讨论

本文结合三角测量的实例对误差传播定律一些应用问题进行了讨论,指出当只有一个独立的观测值时,和函数与倍数函数运用误差传播定律不会出现悖论;采用数学中更为复杂的恒等函数关系式中不同的算式求解相同观测值的函数值,运用误差传播定律也不会出现悖论。如果在测量工作中有多余的直接观测值,就需用平差后的间接观测值按协方差传播律来计算,这样数学中相等的函数关系才能得到同样的函数中误差结果。     

Optimal multi-step collocation: application to the space-wise approach for GOCE data analysis

Collocation is widely used in physical geodesy. Its application requires to solve systems with a dimension equal to the number of observations, causing numerical problems when many observations are available. To overcome this drawback, tailored step-wise techniques are usually applied. An example of these step-wise techniques is the space-wise approach to the GOCE mission data processing. The original idea of this approach was to implement a two-step procedure, which consists of first predicting gridded values at satellite altitude by collocation and then deriving the geo-potential spherical harmonic coefficients by numerical integration. The idea was generalized to a multi-step iterative procedure by introducing a time-wise Wiener filter to reduce the highly correlated observation noise. Recent studies have shown how to optimize the original two-step procedure, while the theoretical optimization of the full multi-step procedure is investigated in this work. An iterative operator is derived so that the final estimated spherical harmonic coefficients are optimal with respect to the Wiener–Kolmogorov principle, as if they were estimated by a direct collocation. The logical scheme used to derive this optimal operator can be applied not only in the case of the space-wise approach but, in general, for any case of step-wise collocation. Several numerical tests based on simulated realistic GOCE data are performed. The results show that adding a pre-processing time-wise filter to the two-step procedure of data gridding and spherical harmonic analysis is useful, in the sense that the accuracy of the estimated geo-potential coefficients is improved. This happens because, in its practical implementation, the gridding is made by collocation over local patches of data, while the observation noise has a time-correlation so long that it cannot be treated inside the patch size. Therefore, the multi-step operator, which is in theory equivalent to the two-step operator and to the direct collocation, is in practice superior thanks to the time-wise filter that reduces the noise correlation before the gridding. The criteria for the choice of this filter are investigated numerically.     

The celestial mechanics approach: application to data of the GRACE mission

The celestial mechanics approach (CMA) has its roots in the Bernese GPS software and was extensively used for determining the orbits of high-orbiting satellites. The CMA was extended to determine the orbits of Low Earth Orbiting satellites (LEOs) equipped with GPS receivers and of constellations of LEOs equipped in addition with inter-satellite links. In recent years the CMA was further developed and used for gravity field determination. The CMA was developed by the Astronomical Institute of the University of Bern (AIUB). The CMA is presented from the theoretical perspective in (Beutler et al. 2010). The key elements of the CMA are illustrated here using data from 50 days of GPS, K-Band, and accelerometer observations gathered by the Gravity Recovery And Climate Experiment (GRACE) mission in 2007. We study in particular the impact of (1) analyzing different observables [Global Positioning System (GPS) observations only, inter-satellite measurements only], (2) analyzing a combination of observations of different types on the level of the normal equation systems (NEQs), (3) using accelerometer data, (4) different orbit parametrizations (short-arc, reduced-dynamic) by imposing different constraints on the stochastic orbit parameters, and (5) using either the inter-satellite ranges or their time derivatives. The so-called GRACE baseline, i.e., the achievable accuracy of the GRACE gravity field for a particular solution strategy, is established for the CMA.     

Instability in spatial error models: an application to the hypothesis of convergence in the European case

This paper focuses on the hypothesis of stability in the mechanisms of spatial dependence that are usually employed in spatial econometric models. We propose a specification strategy for which the first step is to solve a local estimation algorithm, called the Zoom estimation. The aim of this stage is to detect problems of heterogeneity in the parameters and to identify the regimes. Then we resort to a battery of formal Lagrange Multipliers to test the assumption of stability in the processes of spatial dependence. The alternative hypothesis consists of the existence of several regimes in these parameters. A small Monte Carlo serves to confirm the behaviour of this strategy in a context of finite size samples. As an illustration, we solve an application to the case of the hypothesis of convergence for the per capita income in the European regions. Our results reveal the existence of a strong Centre-Periphery dichotomy in which instability extends to all the elements (coefficients of regression as well as parameters of spatial dependence) that intervene in a classical conditional β-convergence model.     

Rigorous transformation of variance–covariance matrices of GPS-derived coordinates and velocities

With the advances in the field of GPS positioning and the global densification of permanent GPS tracking stations, it is now possible to determine at the highest level of accuracy the transformation parameters connecting various international terrestrial reference frame (ITRF) realizations. As a by-product of these refinements, not only the seven usual parameters of the similarity transformations between frames are available, but also their rates, all given at some epoch t _k . This paper introduces rigorous matrix equations to estimate variance–covariance matrices for transformed coordinates at any epoch t based on a stochastic model that takes into consideration all a priori information of the parameters involved at epoch t _k , and the coordinates and velocities at the reference frame initial epoch t ₀. The results of this investigation suggest that in order to attain maximum accuracy, the agencies determining the 14-parameter transformations between reference frames should also publish their full variance–covariance matrix. Electronic Publication     

Validation of GOCE gravity field models by means of orbit residuals and geoid comparisons

Three GOCE-based gravity field solutions have been computed by ESA’s high-level processing facility and were released to the user community. All models are accompanied by variance-covariance information resulting either from the least squares procedure or a Monte-Carlo approach. In order to obtain independent external quality parameters and to assess the current performance of these models, a set of independent tests based on satellite orbit determination and geoid comparisons is applied. Both test methods can be regarded as complementary because they either investigate the performance in the long wavelength spectral domain (orbit determination) or in the spatial domain (geoid comparisons). The test procedure was applied to the three GOCE gravity field solutions and to a number of selected pre-launch models for comparison. Orbit determination results suggest, that a pure GOCE gravity field model does not outperform the multi-year GRACE gravity field solutions. This was expected as GOCE is designed to improve the determination of the medium to high frequencies of the Earth gravity field (in the range of degree and order 50 to 200). Nevertheless, in case of an optimal combination of GOCE and GRACE data, orbit determination results should not deteriorate. So this validation procedure can also be used for testing the optimality of the approach adopted for producing combined GOCE and GRACE models. Results from geoid comparisons indicate that with the 2 months of GOCE data a significant improvement in the determination of the spherical harmonic spectrum of the global gravity field between degree 50 and 200 can be reached. Even though the ultimate mission goal has not yet been reached, especially due to the limited time span of used GOCE data (only 2 months), it was found that existing satellite-only gravity field models, which are based on 7 years of GRACE data, can already be enhanced in terms of spatial resolution. It is expected that with the accumulation of more GOCE data the gravity field model resolution and quality can be further enhanced, and the GOCE mission goal of 1–2 cm geoid accuracy with 100 km spatial resolution can be achieved.     

Optimum data weighting and error calibration for estimation of gravitational parameters

Summary A new approach has been developed for determining consistent satellite-tracking data weights in solutions for the satellite-only gravitational models. The method employs subset least-squares solutions of the satellite data contained within the complete solution and requires that the differences of the parameters of subset solutions and the complete solution to be in agreement with their error estimates by adjusting the data weights. GEM-T2 model was recently computed and adjusted through a direct application of this method. The estimated data weights are markedly smaller than the weights implied by the formal uncertainties of the measurements. Orbital arc tests as well as surface gravity comparisons show significant improvements for solutions when more realistic data weighting is achieved.     

Time transfer using GPS carrier phase: error propagation and results

 A joint time-transfer project between the Astronomical Institute of the University of Berne (AIUB) and the Swiss Federal Office of Metrology and Accreditation (METAS) was initiated to investigate the power of the time transfer using GPS carrier phase observations. Studies carried out in the context of this project are presented. The error propagation for the time-transfer solution using GPS carrier phase observations was investigated. To this purpose a simulation study was performed. Special interest was focussed on errors in the vertical component of the station position, antenna phase-center variations and orbit errors. A constant error in the vertical component introduces a drift in the time-transfer results for long baselines in east–west directions. The simulation study was completed by investigating the profit for time transfer when introducing the integer carrier phase ambiguities from a double-difference solution. This may reduce the drift in the time-transfer results caused by constant vertical error sources. The results from the present time-transfer solution are shown in comparison to results obtained with independent time-transfer techniques. The interpretation of the comparison benefits from the investigations of the error propagation study. Two types of solutions are produced on a regular basis at AIUB: one based on the rapid orbits from CODE, the other on the CODE final orbits. The rapid solution is available the day after the observations and has nearly the same quality as the final solution, which has a latency of about one week. The differences between these two solutions are below the nanosecond level. The differences from independent time-transfer techniques such as TWSTFT (two-way satellite time and frequency transfer) are a few nanoseconds for both products. Received: 15 November 2001 / Accepted: 6 September 2002 Correspondence to:R. Dach     

一种Kalman滤波系统误差及其协方差矩阵的半参数估计方法

提出用半参数估计理论来解决系统误差对Kalman滤波解的影响问题。即用半参数模型中的非参数分量表达观测模型和动力学模型中未知的系统误差,在移动的窗口内,基于观测残差和状态向量预测残差拟合模型系统误差,进而修正相应的观测向量和状态预测向量的协方差矩阵,以消除系统误差对滤波的影响。同时这种方法还有明显的优点,就是在滤波过程中不需要对系统误差做任何假设。文中推导了基于正则核估计来解算导航系统半参数模型的相应公式,并根据一个模拟的算例,证明了算法的有效性。     

联合GRACE/GOCE重力场模型和GPS/水准数据确定我国85高程基准重力位

GRACE、GOCE卫星重力计划的实施,对确定高精度重力场模型具有重要贡献。联合GRACE、GOCE卫星数据建立的重力场模型和我国均匀分布的649个GPS/水准数据可以确定我国高程基准重力位,但我国高程基准对应的参考面为似大地水准面,是非等位面,将似大地水准面转化为大地水准面后确定的大地水准面重力位为62 636 854.395 3m~2s~(-2),为提高高阶项对确定大地水准面的贡献,利用高分辨率重力场模型EGM2008扩展GRACE/GOCE模型至2190阶,同时将重力场模型和GPS/水准数据统一到同一参考框架和潮汐系统,最后利用扩展后的模型确定的我国大地水准面重力位为62 636 852.751 8m~2s~(-2)。其中组合模型TIM_R4+EGM2008确定的我国85高程基准重力位值62 636 852.704 5m~2s~(-2)精度最高。重力场模型截断误差对确定我国大地水准面的影响约16cm,潮汐系统影响约4~6cm。     

非球形引力位中J_3项对轨道的影响及应用

地球的非球形引力摄动是人造卫星环绕地球飞行过程中所受到的摄动中较重要的一项。在人造卫星的精密定轨和卫星位置的预报中,人们往往对J2项比较重视,但是对于J3项的研究比较少。为了解J3项的性质,文中分析了非球形引力位中J3项对于卫星轨道的影响,通过严格的公式推导,得到J3项对轨道偏心率以及轨道倾角的影响,详细地介绍利用它形成冻结轨道的基本原理。通过数值实验对得到的基本理论以及应用进行验证。结果表明J3项是冻结轨道形成的重要因素。     

Generalyzed variance-covariance propagation law formulae and application to explicit least-squares adjustments

Standard formulae overlook the contribution of a number of terms in the derivation of variance-covariance matrices for parameters in nonlinear least squares adjustment. In a large class of nonlinear mathematical models, these terms can contribute to an important error in the estimation of parameter variances. Improved formulae are derived. A numerical example is given and the use of our improved formula in the case of least-squares adjustment in the explicit case (L=F(X)) is fully documented.     

Formation of surface spherical harmonic normal matrices and application to high-degree geopotential modeling

 The structure of normal matrices occurring in the problem of weighted least-squares spherical harmonic analysis of measurements scattered on a sphere with random noises is investigated. Efficient algorithms for the formation of the normal matrices are derived using fundamental relations inherent to the products of two surface spherical harmonic functions. The whole elements of a normal matrix complete to spherical harmonic degree L are recursively obtained from its first row or first column extended to degree 2L with only O(L ⁴) computational operations. Applications of the algorithms to the formation of surface normal matrices from geoid undulations and surface gravity anomalies are discussed in connection with the high-degree geopotential modeling. Received: 22 March 1999 / Accepted: 23 December 1999     

Evaluation of the third- and fourth-generation GOCE Earth gravity field models with Australian terrestrial gravity data in spherical harmonics

In March 2013, the fourth generation of European Space Agency’s (ESA) global gravity field models, DIR4 (Bruinsma et al. in Proceedings of the ESA living planet symposium, 28 June–2 July, Bergen, ESA, Publication SP-686, 2010b) and TIM4 (Migliaccio et al. in Proceedings of the ESA living planet symposium, 28 June–2 July, Bergen, ESA, Publication SP-686, 2010), generated from the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gravity observation satellite was released. We evaluate the models using an independent ground truth data set of gravity anomalies over Australia. Combined with Gravity Recovery and Climate Experiment (GRACE) satellite gravity, a new gravity model is obtained that is used to perform comparisons with GOCE models in spherical harmonics. Over Australia, the new gravity model proves to have significantly higher accuracy in the degrees below 120 as compared to EGM2008 and seems to be at least comparable to the accuracy of this model between degree 150 and degree 260. Comparisons in terms of residual quasi-geoid heights, gravity disturbances, and radial gravity gradients evaluated on the ellipsoid and at approximate GOCE mean satellite altitude ( $h=250$  km) show both fourth generation models to improve significantly w.r.t. their predecessors. Relatively, we find a root-mean-square improvement of 39 % for the DIR4 and 23 % for TIM4 over the respective third release models at a spatial scale of 100 km (degree 200). In terms of absolute errors, TIM4 is found to perform slightly better in the bands from degree 120 up to degree 160 and DIR4 is found to perform slightly better than TIM4 from degree 170 up to degree 250. Our analyses cannot confirm the DIR4 formal error of 1 cm geoid height (0.35 mGal in terms of gravity) at degree 200. The formal errors of TIM4, with 3.2 cm geoid height (0.9 mGal in terms of gravity) at degree 200, seem to be realistic. Due to combination with GRACE and SLR data, the DIR models, at satellite altitude, clearly show lower RMS values compared to TIM models in the long wavelength part of the spectrum (below degree and order 120). Our study shows different spectral sensitivity of different functionals at ground level and at GOCE satellite altitude and establishes the link among these findings and the Meissl scheme (Rummel and van Gelderen in Manusrcipta Geodaetica 20:379–385, 1995).