青藏高原中部的东西向扩张构造运动

系统分析了1933~2003年间青藏高原及其周缘发生的745个中、强地震的震源机制解,研究了高原地壳构造运动及其动力学特征。结果表明,大量正断层型地震集中发生在青藏高原中部海拔4000m以上的地区,其中许多地震是纯正断层型地震。震源机制结果显示,该区正断层型地震的断层走向多为南北方向,断层位错矢量的水平分量均位于近东西方向,这表明青藏高原高海拔地区存在着近东西方向的扩张构造运动。地震震源应力场的研究结果表明,在高原中部高海拔地区,E-W向或WNW-ESE向的水平扩张作用控制着该区的地壳应力场。青藏高原高海拔地区近东西方向的扩张构造运动是该区引张应力场的作用结果,其动力学原因可能与持续隆升的高原自重增大引起的重力崩塌及其周边区域构造应力状况有关。而青藏高原周缘地区,除了东部边缘外,南部的喜马拉雅山前沿以及青藏高原的北部、西部边缘所发生的绝大部分地震都是逆断层型或走滑逆断层型地震。在青藏高原周缘地区,北东或者北北东方向水平挤压的构造应力场为优势应力场。在中国西部的大范围内,主压应力P轴水平分量位于NE-SW方向,形成了一个广域的NE-SW方向的挤压应力场。青藏高原及其周缘应力场特征表明,印度板块的北上运动以及它与欧亚板块之间的碰撞所形成的挤压应力场是高原强烈隆起的直接原因。在青藏高原中南部形成了近东西向引张应力场为主的区域,并以东西向扩张构造运动部分释放其应力积累。研究高原高海拔地区的引张应力场和近东西向扩张构造运动的特征,对于认识青藏高原强烈隆起的地球动力学过程与机制,有着重要的理论意义。     

邯峰矿区古构造应力场与构造演化

采用共轭剪节理应力反演方法,恢复了邯郸-峰峰矿区晚古生代以来的3期古构造应力场,进而探讨了煤田构造的演化历史,将其分为4大阶段：①中生代早期近NS向挤压,煤系后期改造初动期;②中生代晚期SE-NW向挤压,奠定煤田构造格架的基础;③中生代末至古近纪NW-SE向拉张,煤田构造格架定型;④新近纪以来近东西向拉张,煤田构造的现代活动.     

内蒙古中部重、磁场特征与地壳密度结构

中亚造山带中的古生代——中生代花岗岩普遍具有正εNd值，在世界上是十分独特的。美国西部加利福尼亚中生代-新生代花岗岩同样具有正εNd值，并且其地壳的速度分布特征与内蒙古中部十分相似。本文通过与美国加利福尼亚的地球物理研究成果进行对比，来研究内蒙古中部的地壳密度结构，特别是下地壳的组成，试图探讨产生正￡。值花岗岩的深部地质原因。研究结果表明，在大兴安岭-内蒙古造山带下地壳中可能存在与美国加里福尼亚类似的洋壳物质。此外，还对重、磁异常进行了处理，以确定蛇绿岩带的延展情况。由于在西拉木伦河附近存在切割至莫霍面的深断裂，结合地表出现的蛇绿岩带，故提出温都尔庙-西拉木伦河一线可能是内蒙古中部最重要的地质构造界线。     

三峡深水高土石围堰工程的渗流研究

为了对三峡工程二期高土石围堰防渗设施的布置方案及其阻渗效果进行比较,采用有限元法对二期高土石围堰在不利的运行工况下的渗流场进行数值分析.比较了双排混凝土防渗墙方案、单排塑性混凝土低防渗墙加土工膜斜墙和单排厚塑性混凝土防渗墙等3种方案,分别采用恒定与非恒定模型计算,对立面二维和三维绕渗及防渗墙局部开裂等不利工况分别进行了数值模拟.结果表明:三方案均可有效抑制渗流场,双排混凝土防渗墙的防渗效果最好;墙体的局部开裂仅对局部区域的流场有影响.非恒定数值分析表明,堰体、基础不均匀沙石料及基坑抽水速度对渗流场影响极大,为保证堰体稳定,应限制基坑水位降落速度小于2 m/d.     

"嫦娥"卫星绕月飞行轨道的激光测定法

依据地月激光测距的成功实践和对卫星激光定轨的基础研究,提出了用地面对嫦娥卫星作激光测距的方法,高精度地测定嫦娥卫星绕月飞行时的实时在轨位置,论述了多站激光定轨和单站激光定轨的解算数模。     

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.     

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     

Rapid orbit determination of LEO satellites using IGS clock and ephemeris products

Different types of GPS clock and orbit data provided by the International GPS Service (IGS) have been used to assess the accuracy of rapid orbit determination for satellites in low Earth orbit (LEO) using spaceborne GPS measurements. To avoid the need for reference measurements from ground-based reference receivers, the analysis is based on an undifferenced processing of GPS code and carrier-phase measurements. Special attention is therefore given to the quality of GPS clock data that directly affects the resulting orbit determination accuracy. Interpolation of clock data from the available 15 min grid points is identified as a limiting factor in the use of IGS ultra-rapid ephemerides. Despite this restriction, a 10-cm orbit determination accuracy can be obtained with these products data as demonstrated for the GRACE-B spacecraft during selected data arcs between 2002 and 2004. This performance may be compared with a 5-cm orbit determination accuracy achievable with IGS rapid and final products using 5 min clock samples. For improved accuracy, high-rate (30 s) clock solutions are recommended that are presently only available from individual IGS centers. Likewise, a reduced latency and more frequent updates of IGS ultra-rapid ephemerides are desirable to meet the requirements of upcoming satellite missions for near real-time and precise orbit determination.     

Spaceborne gravitational field determination by means of locally supported wavelets

In space-borne gravitational field determination, two challenges are inherent. First, the continuation of the data down to the surface of the Earth is an ill-posed problem, requiring therefore regularization techniques. Second huge data sets result requiring efficient numerical methods. In this paper, we show how locally supported wavelets on the sphere can be developed by means of a spherical version of the so-called up function. By construction, the corresponding scaling functions and wavelets are infinitely smooth, so that they can be used for regularization purposes. In particular, we show how the ill-posed pseudo-differential equations coming from satellite missions can be regularized by efficient numerical schemes using locally supported wavelets. These methods seem in particular to be interesting for regional gravity field modelling.