定位卫星轨道监视系统研究

与卫星测控中心的卫星轨道参数相比,地球同步卫星定位系统中轨道计算子系统所计算出的卫星星历数据具有实时性强、精度高的特点,如何有效利用该星历数据来分析卫星轨道运行情况是该领域的一个热门研究课题。在介绍卫星轨道测定法的基础上,提出了以卫星速度变化率和卫星星下点为主要监视内容的卫星轨道监视系统,通过系统的试运行结果可以看出,该系统能够有效地监视卫星的瞬时和长期运行情况,为卫星的在轨运行管理提供了一种直观、形象的监视方法。     

风云卫星遥感数据高精度地理定位软件系统开发研究

为了适应风云气象卫星遥感数据高精度地理定位要求,在研究和比较分析现有环境气象卫星遥感数据地理定位方法,尤其对其中关键的轨道计算模型进行对比研究之后,我们使用变阶变步长多步卫星轨道数值计算模型(DE/DEABM)进行气象卫星轨道计算,研制开发了新一代卫星轨道计算及遥感数据地理定位软件系统.该软件系统中卫星轨道数值积分计算模型包含了多项摄动因素计算,特别是对低轨卫星影响较大的因素,其中地球的球形引力项使用了高精度高阶EGM-96地球引力场模型,提高了非球形引力摄动计算精度,另外还考虑了太阳、月亮引力项,辐射光压摄动和大气摄动因素,使得轨道计算精度大大提高.在遥感数据定位算法开发工作中,以中国2002年5月发射的风云1号D星10通道扫描辐射计和计划2008年上半年发射的风云3号A卫星中分辨率光谱成像仪等遥感仪器为对象,比较详细地分析和研究了探测器、焦平面、主光学系统和扫描镜等遥感仪器几大关键部件的光学几何关系,提高了坐标转换系统计算精度.经过对风云1号D星多天逐日的轨道计算和定位计算试验,结果表明该软件系统24h风云1号D星轨道计算卫星矢径精度可达到几十厘米至几十米,较原来平根数分析解方法1000m左右精度有显著量级的提高.同时,风云1号D星遥感数据地理定位精度达到星下点1个像元.     

点质模型求定月球重力场及其特征分析

利用Lunar-Prospector扩展任务的视线加速度数据,根据点质模型恢复了月球近区的重力场,将其与LP165月球重力场模型进行了比较和分析,并利用恢复的重力场联合月球地形数据对Mas-con进行分析,总结了月球重力场的主要特征及其研究方向发展趋势,将重力场恢复技术与我国的探月计划-"嫦娥"工程结合起来,为建立高精度的月球重力场和地形模型提供了一种有效而实用的方法。     

随机Dirichlet问题解的调和展开

针对重力学随机Dirichlet问题,通过适当地对边界检验函数的分解,并在随机边界样本空间中提取确定性部分的对偶基,本文将随机Dirichlet问题的一般解展开为一随机系数的调和级数形式。     

贵州中东部剩余重力异常与区域矿产分布关系

重点研究1:20万剩余重力异常,结合航磁、物性、地质等资料,初步圈定区内隐伏、半隐伏岩浆岩体,对研究与重力异常、岩体分布关系密切的矿产的时间、空间分布及其成因,有一定参考价值和指示意义。     

一个新的沉积变质型钛矿床的成矿模式及找矿方法

内蒙古羊蹄子山沉积变质型钛矿床是一个新的钛矿床类型。通过对地质、地球物理前提条件的研究和找矿方法试验,认为高精度磁法是找沉积变质型钛矿床的基本方法,重力为有效的找矿方法。在矿区北带磨石山以金红石为主的富钛矿体上有600 μGal的重力异常,无磁异常;南带的羊蹄子山钛铁矿体上重力异常为500 μGal,有60 nT的弱磁异常。     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals

This paper presents a general modeling strategy for ambiguity resolution (AR) and position estimation (PE) using three or more phase-based ranging signals from a global navigation satellite system (GNSS). The proposed strategy will identify three best “virtual” signals to allow for more reliable AR under certain observational conditions characterized by ionospheric and tropospheric delay variability, level of phase noise and orbit accuracy. The selected virtual signals suffer from minimal or relatively low ionospheric effects, and thus are known as ionosphere-reduced virtual signals. As a result, the ionospheric parameters in the geometry-based observational models can be eliminated for long baselines, typically those of length tens to hundreds of kilometres. The proposed modeling comprises three major steps. Step 1 is the geometry-free determination of the extra-widelane (EWL) formed between the two closest L-band carrier measurements, directly from the two corresponding code measurements. Step 2 forms the second EWL signal and resolves the integer ambiguity with a geometry-based estimator alone or together with the first EWL. This is followed by a procedure to correct for the first-order ionospheric delay using the two ambiguity-fixed widelane (WL) signals derived from the integer-fixed EWL signals. Step 3 finds an independent narrow-lane (NL) signal, which is used together with a refined WL to resolve NL ambiguity with geometry-based integer estimation and search algorithms. As a result, the above two AR processes performed with WL/NL and EWL/WL signals respectively, either in sequence or in parallel, can support real time kinematic (RTK) positioning over baselines of tens to hundreds of kilometres, thus enabling centimetre-to-decimentre positioning at the local, regional and even global scales in the future.     

A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.