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.     

GGM02 – An improved Earth gravity field model from GRACE

A new generation of Earth gravity field models called GGM02 are derived using approximately 14 months of data spanning from April 2002 to December 2003 from the Gravity Recovery And Climate Experiment (GRACE). Relative to the preceding generation, GGM01, there have been improvements to the data products, the gravity estimation methods and the background models. Based on the calibrated covariances, GGM02 (both the GRACE-only model GGM02S and the combination model GGM02C) represents an improvement greater than a factor of two over the previous GGM01 models. Error estimates indicate a cumulative error less than 1 cm geoid height to spherical harmonic degree 70, which can be said to have met the GRACE minimum mission goals. Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Investigating the propagation mechanism of unmodelled systematic errors on coordinate time series estimated using least squares

The propagation of unmodelled systematic errors into coordinate time series computed using least squares is investigated, to improve the understanding of unexplained signals and apparent noise in geodetic (especially GPS) coordinate time series. Such coordinate time series are invariably based on a functional model linearised using only zero and first-order terms of a (Taylor) series expansion about the approximate coordinates of the unknown point. The effect of such truncation errors is investigated through the derivation of a generalised systematic error model for the simple case of range observations from a single known reference point to a point which is assumed to be at rest by the least squares model but is in fact in motion. The systematic error function for a one pseudo-satellite two-dimensional case, designed to be as simple but as analogous to GPS positioning as possible, is quantified. It is shown that the combination of a moving reference point and unmodelled periodic displacement at the unknown point of interest, due to ocean tide loading, for example, results in an output coordinate time series containing many periodic terms when only zero and first-order expansion terms are used in the linearisation of the functional model. The amplitude, phase and period of these terms is dependent on the input amplitude, the locations of the unknown point and reference point, and the period of the reference point's motion. The dominant output signals that arise due to truncation errors match those found in coordinate time series obtained from both simulated data and real three-dimensional GPS data.     

Digital topography models for Martian surfaces

We propose to use an unsupervised automated classification of topographic features on Mars in order to speed up geomorphic and geologic mapping of the planet. We construct a digital topography model (DTM), a multilayer grid that stores various kinds of topographical information for every pixel in a site. The method uses a probabilistic clustering algorithm to assign topographically meaningful labels to all pixels in the DTM. The results are displayed as a thematic map of topography. Resultant topographical features are characterized and compared using statistics of their constituent pixels. We demonstrate the usage of our method by classifying and characterizing the topography of a landscape in the Tisia Valles region on Mars. We discuss extensions and further applications of our method.     

Analysis of geophysical networks derived from multiscale digital elevation models: a morphological approach

We provide a simple and elegant framework based on morphological transformations to generate multiscale digital elevation models (DEMs) and to extract topologically significant multiscale geophysical networks. These terrain features at multiple scales are collectively useful in deriving scaling laws, which exhibit several significant terrain characteristics. We present results derived from a part of Cameron Highlands DEM.     

斜交型不稳定谱点分布侧半圆定理及其增长率上界的估计

本对斜交型扰动不稳定谱点的分布做了理论分析，得到了该谱点分布的半圆定理一该谱点分布在复一面上以原点为圆心以R0为半径的上半平面上，同时还对该不稳定增长率的上界作了估计。发现水平永度越小，模式顶越高则该估计值越大；垂直风切变的增大和纬度的增高对该增长率的增大有正贡献；当层结稳定度减小时，最大增长率随相对最大增长率得增大而减小。     

青藏高原近地面层微气象学特征

利用1998年5月-7月在改则、当雄和昌都三测站获得的近地面层气象要素变化的观测资料，分析了青藏高原近地面层风速、温度和湿度日变化特征及廓线规律，发现高原近地面层微气象学特征具有自己的特点；同时还讨论了高原近地面层白天出现的逆湿现象。     

一类广义Hamilton系统的辛格式

通过引人泊松括号，分析了无限维Hamilton的性质，并将其推广到广义Hamilton系统，且从理论和实用角度讨论了这类广义Hamilton系统的辛格式构造问题，从而为辛几何算法在一般的时间发展方程的数值求解提供新的具体途径。