Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Practical Research of IKONOS-2 Geopositioning and Its Accuracy in Tibet of China

This paper presents the studies of the refining of IKONOS-2 RPC, the transform of the datum, the mode of the control point distribution and the method of IKONOS stereo triangulation, so that IKONOS imagery can be used to collect the precise geospatial data and produce the large scale map. The transform between the IKONOS-2 image space and the national coordinate system based on the RPC have been developed, and the results of block adjustment with various control schemes in a practical project near Himalayas have been examined and analysed. The encouraging results of high positioning accuracy have been obtained.     

藏南地区金和锑矿床(点)类型及其时空分布特征

本文首次对藏南地区金和锑矿床(点)地质特征、成因类型和空间分布特点进行了总结，对金和锑成矿作用与中新生代构造一岩浆活动的关系进行了论证，对区域地壳演化过程中金和锑成矿动力学机制进行了讨论。该区的金和锑矿床(点)可按容矿围岩划分为：①变质岩为主要容矿围岩矿床(点)，包括金和金一锑矿床(点)；②沉积岩为主要容矿围岩矿床(点)，包括金、锑、金一锑和锑多金属矿床(点)；③火山岩为主要容矿围岩矿床(点)，以金一银矿床为代表。根据金和锑矿床(点)空间分布特点和地质特点，将本区划分为4个矿化集中区：拉昂错一马攸木一帕羊金矿化集中区(A)，然勒金和锑矿化集中区(B)，洛扎一措美锑矿化集中区(C)和邛多江金一锑矿化集中区(D)。各矿化集中区内大多数金和锑矿床(点)与燕山晚期～喜马拉雅早期富碱火成岩具密切时空分布关系，它们是古板块对接碰撞期和碰撞期后大规模构造一岩浆活动的产物。近东西向挤压性和南北向张性断裂交汇部位以及富碱火成岩发育区是寻找金和锑矿床(点)的有利场所。     

基于MPI的机群并行计算系统平台构建

在高性能计算机领域,机群并行系统已成为一种重要的系统结构。这里介绍了机群并行系统的特点和发展现状;给出了将PC机或工作站通过高速以太网连接,使用TCP/IP作为标准的通信协议,利用MPI作为分布式的并行计算软件环境,在Windows平台和Unix平台上搭建用户自己的PC机群系统的两种方法;还给出了MPI安装中需要注意的关键步骤。实践表明,PC机群系统具有高性能、高可用性和极高的性价比。     

不等电极距情况下高密度电阻率法的反演结果分析

通过高密度电阻率法的模型实验,得出不等电极距情况对直立高阻薄板和低阻球体反演结果的影响:缺失电极对它正下方的ρ_s影响最大,而异常形态主要在水平方向有明显的变化,异常体随电极缺失的位置有趋势变化,电极缺失的数目越多,异常形态变化越大,但是从整体来看反演剖面图基本相似。     

长江中下游夏季高温灾害机理及预测

利用我国1961-2003年夏季（6—8月）高温资料，建立长江中下游地区主要城市强高温及高温过程较完整的时间序列，探讨了该地区主要城市高温气候特征。分析该地区南京、杭州、南昌等城市夏季高温灾害机理，东亚副热带高压是造成长江中下游地区城市夏季高温的主要影响系统。在此基础上用均生函数-最佳子回归集构造预测模型，预测夏季月高温出现日数，通过42a高温资料预报检验，有较好的预测效果，值得在业务中应用。     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

A parametric study on the impact of satellite attitude errors on GOCE gravity field recovery

In the recent design of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, the gravity gradients are defined in the gradiometer reference frame (GRF), which deviates from the actual flight direction (local orbit reference frame, LORF) by up to 3–4°. The main objective of this paper is to investigate the effect of uncertainties in the knowledge of the gradiometer orientation due to attitude reconstitution errors on the gravity field solution. In the framework of several numerical simulations, which are based on a realistic mission configuration, different scenarios are investigated, to provide the accuracy requirements of the orientation information. It turns out that orientation errors have to be seriously considered, because they may represent a significant error component of the gravity field solution. While in a realistic mission scenario (colored gradiometer noise) the gravity field solutions are quite insensitive to small orientation biases, random noise applied to the attitude information can have a considerable impact on the accuracy of the resolved gravity field models.