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.     

GOCE卫星重力探测任务

通过分析利用卫星观测技术恢复地球重力场的现状与不足,论述了实施GOCE任务———地球重力场和稳态海洋环流实验的必要性;并阐述了GOCE任务的发展历史及现状、基本特点、组成部分、重力梯度测量原理、数据处理过程、误差特性,以及该任务在固体地球物理学、大地测量学、海洋学、冰河学等地球物理相关科学领域中的应用。     

Scientific objectives of the Solar Mesosphere Explorer mission

The 1981–82 Solar Mesosphere Explorer (SME) mission is described. The SME experiment will provide a comprehensive study of mesospheric ozone and the processes which form and destroy it. Five instruments will be carried on the spinning spacecraft to measure the ozone density and its altitude distribution from 30 to 80 km, monitor the incoming solar ultraviolet radiation, and measure other atmospheric constituent which affect ozone. The polar-orbiting spacecraft will be placed into a 3pm-3 am Sun-synchronous orbit. The atmospheric measurements will scan the Earth's limb and measure: (1) the mesospheric and stratospheric ozone density distribution by inversion of Rayleigh-scattered ultraviolet limb radiance, and the thermal emission from ozone at 9.6 m; (2) the water vapor density distribution by inversion of thermal emission at 6.3 m; (3) the ozone photolysis rate by inversion of the O₂(¹g) 1.27 m limb radiance; (4) the temperature profile by a combination of narrow-band and wide-band measurements of the 15 m thermal emission by CO₂; and, (5) theNO₂ density distribution by inversion of Rayleighscattered limb radiance at 0.439 m. The solar ultraviolet monitor will measure both the 0.2–0.31 m spectral region and the Lyman-alpha (0.1216 m) contribution to the solar irradiance. This combination of measurements will provide a rigorous test of the photochemical equilibrium theory of the mesospheric oxygen-hydrogen system, will determine what changes occur in the ozone distribution as a result of changes in the incoming solar radiation, and will detect changes that may occur as a result of meteorological disturbances.     

Assessment of three numerical solution strategies for gravity field recovery from GOCE satellite gravity gradiometry implemented on a parallel platform

 The recovery of a full set of gravity field parameters from satellite gravity gradiometry (SGG) is a huge numerical and computational task. In practice, parallel computing has to be applied to estimate the more than 90 000 harmonic coefficients parameterizing the Earth's gravity field up to a maximum spherical harmonic degree of 300. Three independent solution strategies (preconditioned conjugate gradient method, semi-analytic approach, and distributed non-approximative adjustment), which are based on different concepts, are assessed and compared both theoretically and on the basis of a realistic-as-possible numerical simulation regarding the accuracy of the results, as well as the computational effort. Special concern is given to the correct treatment of the coloured noise characteristics of the gradiometer. The numerical simulations show that the three methods deliver nearly identical results—even in the case of large data gaps in the observation time series. The newly proposed distributed non-approximative adjustment approach, which is the only one of the three methods that solves the inverse problem in a strict sense, also turns out to be a feasible method for practical applications. Received: 17 December 2001 / Accepted: 17 July 2002 Acknowledgments. We would like to thank Prof. W.-D. Schuh, Institute of Theoretical Geodesy, University of Bonn, for providing us with the serial version of the PCGMA algorithm, which forms the basis for the parallel PCGMA package developed at our institute. This study was partially performed in the course of the GOCE project `From E?tv?s to mGal+', funded by the European Space Agency (ESA) under contract No. 14287/00/NL/DC. Correspondence to: R. Pail     

GOCE�������׼����ܶ���

????PANDA??????GOCE????2009???????306-332??????????????????????????????????GOCE???????????PSO??????????PANDA????????????????????????????????α??????????SLR????????4??????????????????????????????????t????????????????GOCE????????????PSO??￡??侶?????????????????2 cm??     

东北太平洋Explorer Ridge热液羽状流位温浊度异常和物质能量通量估算

深海热液流体与周围海水之间存在明显的物理和化学差异,通过检测海水的位温浊度异常是探测深海热液活动的重要手段之一。本文采用"海底火山带项目（Submarine Ring of Fire 2002）"拖曳式温盐深测量仪数据资料,研究了东北太平洋Explorer Ridge热液场的水文特征及物质能量通量的释放。结果表明Explorer Ridge热液场热液羽状流中性浮力层所在深度范围约为1 600~1 900 m,距离海底的高度约为200 m,最大位温、盐度和浊度异常分别为0.04℃、0.004和0.18 NTU;中性浮力层热液羽状流帽呈椭圆结构,其长轴与洋中脊线重合,羽状流帽总面积约为27 km2;热液羽状流在中性层范围内存在明显的分层现象,通过经验公式计算得到Explorer Ridge热液场观测范围内热液喷口的总的浮力通量为6.19×10-2 m4/s3,平均值为2.063×10-2 m4/s3;总的体积通量为9.884×10-2 m3/s,平均值为3.295×10-2 m3/s;总的热通量为194.9 MW,平均值为64.967 MW。     

基于ArcGIS Explorer建立DOM应用系统的研究

数字正射影像图(DOM)因其获取快捷,内容丰富,直观易读,使用方便等优点,特别适合管理者使用,在城市的规划、建设和管理中发挥着重要作用。现实工作中,数字正射影像图的应用方法很多,本文重点研究基于ESRI的免费软件ArcGIS Explorer建立DOM应用系统的技术、方法和步骤,并简要介绍这一系统的基本功能和实际应用。     

Computation of spherical harmonic coefficients from gravity gradiometry data to be acquired by the GOCE satellite: regularization issues

The issue of optimal regularization is investigated in the context of the processing of satellite gravity gradiometry (SGG) data that will be acquired by the GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) satellite. These data are considered as the input for determination of the Earths gravity field in the form of a series of spherical harmonics. Exploitation of a recently developed fast processing algorithm allowed a very realistic setup of the numerical experiments to be specified, in particular: a non-repeat orbit; 1-s sampling rate; half-year duration of data series; and maximum degree and order set to 300. The first goal of the study is to compare different regularization techniques (regularization matrices). The conclusion is that the first-order Tikhonov regularization matrix (the elements are practically proportional to the degree squared) and the Kaula regularization matrix (the elements are proportional to the fourth power of the degree) are somewhat superior to other regularization techniques. The second goal is to assess the generalized cross-validation method for the selection of the regularization parameter. The inference is that the regularization parameter found this way is very reasonable. The time expenditure required by the generalized cross-validation method remains modest even when a half-year set of SGG data is considered. The numerical study also allows conclusions to be drawn regarding the quality of the Earths gravity field model that can be obtained from the GOCE SGG data. In particular, it is shown that the cumulative geoid height error between degrees 31 and 200 will not exceed 1 cm. AcknowledgmentsThe authors thank Dr. E. Schrama for valuable discussions and for computing the orbit used to generate the long data set. They are also grateful to Prof. Tscherning and two anonymous reviewers for numerous valuable remarks and suggestions. The orbit to generate the short data set was kindly provided by J. van den IJssel. Computing resources were provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.