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.     

GPS observations of the ionospheric F2-layer behavior during the 20th November 2003 geomagnetic storm over South Korea

The ionospheric F2-layer peak density (NmF2) and its height (hmF2) are of great influence on the shape of the ionospheric electron density profile Ne (h) and may be indicative of other physical processes within the ionosphere, especially those due to geomagnetic storms. Such parameters are often estimated using models such as the semiempirical international reference ionosphere (IRI) models or are measured using moderately priced to expensive instrumentation, such as ionosondes or incoherent scatter radars. Global positioning system (GPS) observations have become a powerful tool for mapping high-resolution ionospheric structures, which can be used to study the ionospheric response to geomagnetic storms. In this paper, we describe how 3-D ionospheric electron density profiles were produced from data of the dense permanent Korean GPS network using the tomography reconstruction technique. These profiles are verified by independent ionosonde data. The responses of GPS-derived parameters at the ionospheric F2-layer to the 20th November 2003 geomagnetic storm over South Korea are investigated. A fairly large increase in the electron density at the F2-layer peak (the NmF2) (positive storm) has been observed during this storm, which is accompanied by a significant uplift in the height of the F2 layer peak (the hmF2). This is confirmed by independent ionosonde observations. We suggest that the F2-layer peak height uplift and NmF2 increase are mainly associated with a strong eastward electric field, and are not associated with the increase of the O/N₂ ratio obtained from the GUVI instruments aboard the TIMED satellite. It is also inferred that the increase in NmF2 is not caused by the changes in neutral composition, but is related to other nonchemical effects, such as dynamical changes of vertical ion motions induced by winds and E × B drifts, tides and waves in the mesosphere/lower thermosphere region, which can be dynamically coupled upward to generate ionospheric perturbations and oscillations.     

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.     

非量测数码相机（哈苏H1D型）在南极项目生产中的应用

简要分析了非量测数码相机哈苏H1D在南极中山站及拉斯曼丘陵地区所获取的影像数据的几何精度,阐明了利用该影像数据进行空中三角测量,DLG,DOM制作的技术流程及作业过程中出现问题的处理方法。     

模拟退火算法在控制网平差中的应用

线性最小二乘估计在对非线性函数进行线性近似的过程中会产生模型误差,而一些非线性参数估计方法可能因为函数复杂而难以求导,法方程系数矩阵秩亏或呈病态矩阵时难以求解,非线性迭代解法有时对初始值的选择存在依赖性,不恰当的初始值会导致迭代无法收敛。针对这些问题,引入了模拟退火算法,介绍了该算法的基本原理、计算步骤和收敛性,并以3个控制网平差应用为例,说明该算法具有无需求导求逆,简洁实用,易于编程等优势,并能实现全局优化,获得高精度的平差结果。     

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

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

二三维联动的虚拟校园设计与实现

传统的2维地理信息系统技术成熟但缺乏3维可视化和分析的能力,而目前的3维地理信息系统发展可拥有3维可视化分析能力,然而还缺少2维地理信息系统的便捷。以二三维联动为开发思路的地理信息系统能够解决这对矛盾。本文基于二三维联动GIS的设计方案,以虚拟校园设计与实现为应用实例,利用ArcEngine（AE）实现二三维联动GIS系统的开发设计,对GIS系统的开发做出了一种新的、有益的尝试。     

红枫湖夏季分层期间pCO2分布规律的研究

在夏季分层期间对红枫湖南、北湖湖心的水样进行分层采集,同时测定了分层水样的温度、pH、HCO₃-浓度、溶解氧(DO)、叶绿素a(Chl-a)及铵根离子(NH₄+)、硝酸根离子(NO₃-)、磷酸根离子(PO₄3-)的浓度,水体中CO₂的分压(pCO₂)由化学平衡及亨利定律求得。研究结果表明:光合作用、有机质降解及水体热分层是影响红枫湖夏季pCO₂分布的主要因素。其中,温水层CO₂欠饱和是光合作用吸收CO₂引起的,温跃层中pCO₂的急剧增加是光合产物降解产生CO₂引起的。静水层沉积物附近pCO₂最高并且还有持续增加的趋势,说明沉积物中有机质降解是静水层中CO₂增加的主要原因,夏季湖底水温较高加快了沉积物中有机质的降解。分层现象使pCO₂在水体中的分布差别明显,并且使静水层中CO₂得到积累。此外,夏季红枫湖水体中pCO₂的变化与NH₄+、PO₄3-的变化密不可分,表现为温水层中光合作用消耗NH₄+、PO₄3-,有机质降解过程伴随NH₄+、PO₄3-的释放。     

我国灌区地下水调控研究现状及其发展态势

现代灌溉引水工程的发展,改善了灌区的农业用水条件,但同时改变了灌区上的自然生态环境、地质条件与水文地质条件等。分阶段总结了地下水的相关研究与应用;剖析了地下水调控的必要性与重要性;分析了地下水调控研究现状;最后通过总结给出了进一步研究的建议。     

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.