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.     

阿尔金地区构造应力场及其对金属矿产分布的控制作用

构造应力场研究可以为区域矿产预测提供依据。主要根据节理、擦痕的测试分析，结合区域构造解析，确定阿尔金山东西向拉配泉－－红柳沟构造带经历了三期不同方向的应力作用：印支期以前的南北向挤压作用，印支－燕山期的北西－南东向的挤压和中生代晚期－新生代北东东－南南西向的挤压；并利用有限元数值模拟前两期构造应力场的演变过程，进而讨论了区域构造应力场对内生金属矿产分布规律的控制作用指出了区域寻找大型内生金属矿产的有利地段。     

河北省峰峰矿区通二井田石炭系—二叠系砂岩显微组构特征分析

对河北省峰峰矿区通二井田煤系地层砂岩中石英光轴的组构特征进行了显微镜下的观察分析。实测的12个岩组图均呈大圆环带状,表现为多组的极密区和次极密区,并以单斜对称组构为主体,具有一个对称面和垂直对称面的1个二次轴,以该轴两两对称的极密区和次极密区有2对或2对以上,表明研究区发生过不容忽视的旋转或剪切运动,它们的构造变形史可能具有多次活动性,并伴随有地温上升的过程。根据水平相岩组图所做的主应力轨迹图显示F56断层是形成张家庄向斜构造区和集贤村构造区的内部边界条件,并表明通二井田最大主压应力方向在北部以北西—南东向为主,向南逐渐转为东西向;最小主压应力方向在集贤村构造区以北东向为主,在中南部张家庄向斜构造区则以近南北向为主。     

江苏高邮地震台地磁变化与地震关系的初探

分析了江苏高邮地磁台1980年至2002年地磁垂直分量日变幅逐日比值的高值异常与江苏及邻区ML5．0以上地震之间的关系，发现它们之间有良好的对应关系，通过分析地磁垂直分量谐波幅度的逐日比，提出了适用于高邮地磁台的地震分析预报参数，并对其内检预报效能进行了评价。     

深圳市1 km高分辨率厘米级高精度大地水准面的确定

利用65个精度优于2 cm的GPS水准数据、5 213个实测重力点数据、100 m分辨率的数字地形模型和WDM94地球重力场模型,采用移去-恢复技术计算了深圳市1 km分辨率的大地水准面模型.将该模型大地水准面高与由29个GPS水准得到的大地水准面高进行比较,其差值的标准差为±1.4 cm.     

Geoid determination using adapted reference field, seismic Moho depths and variable density contrast

 The traditional remove-restore technique for geoid computation suffers from two main drawbacks. The first is the assumption of an isostatic hypothesis to compute the compensation masses. The second is the double consideration of the effect of the topographic–isostatic masses within the data window through removing the reference field and the terrain reduction process. To overcome the first disadvantage, the seismic Moho depths, representing, more or less, the actual compensating masses, have been used with variable density anomalies computed by employing the topographic–isostatic mass balance principle. In order to avoid the double consideration of the effect of the topographic–isostatic masses within the data window, the effect of these masses for the used fixed data window, in terms of potential coefficients, has been subtracted from the reference field, yielding an adapted reference field. This adapted reference field has been used for the remove–restore technique. The necessary harmonic analysis of the topographic–isostatic potential using seismic Moho depths with variable density anomalies is given. A wide comparison among geoids computed by the adapted reference field with both the Airy–Heiskanen isostatic model and seismic Moho depths with variable density anomaly and a geoid computed by the traditional remove–restore technique is made. The results show that using seismic Moho depths with variable density anomaly along with the adapted reference field gives the best relative geoid accuracy compared to the GPS/levelling geoid. Received: 3 October 2001 / Accepted: 20 September 2002 Correspondence to: H.A. Abd-Elmotaal     

一个中尺度流复合体的环境场、云场的特征分析

利用常规天气图、物理量诊断分析、逐时红外云图、单站资料等,对2003年4月12日影响鹰潭市的强对流、暴雨天气过程进行分析,认为该过程是在一定的大尺度环流背景下,由不同层次的天气系统共同作用下产生的;其触发机制是冷空气从近地层楔入暖空气中,在高、低空急流、层结不稳定等诸多有利条件下,由几个小尺度对流云团合并发展成1个中尺度对流复合体(MCC).     

Helmholtz算子逆的多方向分裂

Fourier有限差分(FFD)方法是一个重要的波场深度延拓方法，它需要用到包括Taylor展开、连分式晨开、迭代展开等一系列数学展开，将FFD方法从二维推广到三维时，需要解决二维Helmholtz算子的求逆问题。本文利用多方向分裂方法导出了二维Helmholtz算子逆的近似公式，可以应用于三维波场深度延拓问题的研究。     

The spacetime gravitational field of a deformable body

The high-resolution analysis of orbit perturbations of terrestrial artificial satellites has documented that the eigengravitation of a massive body like the Earth changes in time, namely with periodic and aperiodic constituents. For the space-time variation of the gravitational field the action of internal and external volume as well as surface forces on a deformable massive body are responsible. Free of any assumption on the symmetry of the constitution of the deformable body we review the incremental spatial (“Eulerian”) and material (“Lagrangean”) gravitational field equations, in particular the source terms (two constituents: the divergence of the displacement field as well as the projection of the displacement field onto the gradient of the reference mass density function) and the `jump conditions' at the boundary surface of the body as well as at internal interfaces both in linear approximation. A spherical harmonic expansion in terms of multipoles of the incremental Eulerian gravitational potential is presented. Three types of spherical multipoles are identified, namely the dilatation multipoles, the transport displacement multipoles and those multipoles which are generated by mass condensation onto the boundary reference surface or internal interfaces. The degree-one term has been identified as non-zero, thus as a “dipole moment” being responsible for the varying position of the deformable body's mass centre. Finally, for those deformable bodies which enjoy a spherically symmetric constitution, emphasis is on the functional relation between Green functions, namely between Fourier-/ Laplace-transformed volume versus surface Love-Shida functions (h(r),l(r) versus h ^′(r),l ^′(r)) and Love functions k(r) versus k ^′(r). The functional relation is numerically tested for an active tidal force/potential and an active loading force/potential, proving an excellent agreement with experimental results. Received: December 1995 / Accepted: 1 February 1997     

科学计算可视化技术的研究与应用

本文介绍了科学计算可视化技术的发展和应用情况，论述了这一新技术的特点及其在油气藏描述，工程设计ＣＡＤ等学科的应用和发展趋势。阐述了目前国际国内计算机图形学领域具体的研究方向和方法及图形软件、图形硬件的发展概况。