Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

数字测图技术在阜新市地理数据库更新中的应用

采用GPS和全站仪联合作业的数字测图方式进行了外业快速数据采集。并利用南方CASS软件对数据进行内业成图工作,最后利用大型的地理信息系统软件ArcGIS进行地理数据的入库,从而实现了阜新市区地理数据库的更新,此方式具有速度快、效率高和精度高等优点,可以借鉴到同类的地区。     

电离层电子浓度总含量梯度对GPS差分定位精度的影响

针对GPS定位测量的信号传播路径误差,分析了电离层电子浓度总含量梯度对差分定位精度的影响,利用L1载波重点讨论了1999年-2000年太阳活动期间低纬度赤道异常地区GPS差分定位精度的问题,同时验证对于长基线采用GPS广域差分技术可以使电离层定位误差得到明显提高。     

中国西部现今地壳形变与地震活动特征

以1999-2004年间中国西部地区近650个GPS站点的观测资料为基础,采用二维“高张力样条”函数内插算法获得了连续地壳形变场。结果表明,相对于稳定欧亚参考框架,中国西部现今地壳运动西强东弱,北向运动从西向东逐渐减弱,东向运动逐渐增强,青藏高原东缘及附近地区是东向运动的消减区带;主应变从西向东由压缩应变转变为拉伸应变;最大剪应变主要位于喜马拉雅地块及其东部地区;中国西部地区地震活动在空间分布上似乎集中发生在印度板块北东向挤压欧亚板块的两条共扼带内,与地壳形变场有一定对应关系,地壳形变的高应变率区为中强地震多发地带。     

Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.     

兰州—白银多目标地球化学调查野外工作方法浅析

多目标区域地球化学调查从调查区表层、深层采样土壤元素分布特征入手建立土壤元素地球化学分配机制,以景观、地貌及地质特征为基础,选择合理的采样位置,正确的采样方法、采样层位,应用掌上机（Pocket PC）及数字填图系统（Rgmap）采集野外样品,本文以兰州—白银黄土覆盖区景观区为例,对表层、深层采样方法和技术要求进行探讨,对全省开展多目标地球化学调查,特别是野外工作方法提供借鉴。     

河西走廊疏勒河流域地质灾害危险性现状评估方法的应用——以甘肃省河西走廊(疏勒河)农业灌溉暨移民安置综合开发建设项目8个移民社区及配套工程建设为例

依据“建设用地地质灾害危险性评估技术要求”、利用野外踏勘及相关资料,对“甘肃省河西走廊（疏勒河）农业灌溉暨移民安置综合开发建设项目8个移民社区及配套工程”进行地质灾害危险性现状评估。结果认为：评估区发育和分布的地质灾害类型有土壤盐渍化、风蚀沙埋和洪水冲蚀。现状评估结果：土壤盐渍化和洪水冲蚀危险性均小;风蚀沙埋有1处危险性中等,其它危险性小。     

斜坡应力分布的系统测试与反演分析

采用简易、快速、经济的洞壁应力恢复法对某工程边坡勘探平硐洞壁应力场作了系统测试, 并采用FLAC2D软件对测试成果作了斜坡初始应力的反演分析。结果表明, 斜坡应力场随深度变化具有明显分带性, 应力量值随硐深曲线具有典型的驼峰式特征, 与理论分析结果相符, 表明该方法能较好地应用于斜坡应力场的系统测试和研究。成果对认识斜坡岩体的变形破裂及卸荷特征等具有重要意义。     

利用参考重力场模型基于能量法确定GRACE加速度计校准参数

利用多个参考重力场模型分别对GRACE一个月的实测加速度计观测数据进行检校.数值计算结果的比较分析表明了利用参考重力场模型确定加速度计校准参数是有效的.