青藏高原隆升的非线性动态有限元仿真研究

根据青藏高原的地质特征建立分析模型，采用3维动态有限元方法，在计算仿真板块速度场的基础上，计算在青藏高原的隆升过程中该地区地壳岩石的等效应力和位移随时间的变化，计算仿真得到的速度场与1998年GPS观测的速度场吻合良好；与过去一贯的假设相反，计算结果反映出地壳应力场不是静态的，而是此起彼伏，不断变化的，应力值最大且变化最剧烈的地区在克什米尔地区、鄂尔多斯地区和鲜水河－小江断裂带，与地震多发区域吻合。     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

GPS/LEO掩星观测的变分同化技术

在简单介绍GPS/LEO掩星探测大气的发展历史和科学意义之后,详细阐述了反演的基本原理;分析了标准反演中存在的问题,并说明一维变分同化(1DVAR)在反演方法中的重要性;给出了一维变分同化中价值函数的求解,以及各种同化因子;简单介绍了对当前气象学中普遍使用的四维变分同化(4DVAR);重点讨论了各种同化方法,以及使用各种同化因子的优缺点。最后,通过CHAMP卫星的观测实例分析,验证了GPS数据在数值天气预报(NWP)中的作用,以及相对于标准反演法一维变分对气象要素的改进。     

小波分析理论在GPS技术中的应用

在GPS数据处理中 ,存在着误差影响、影响波的干扰、周跳和数据量大等问题。误差影响和影响波的干扰实质是在接收卫星信号时受到其它因素的影响 ;周跳是由于卫星信号的失锁而造成信号的不连续 ;数据量大是因为GPS观测需要采样间隔小又连续观测所致。由于小波理论具有时频分析、波形分解、特征提取和快速小波变换等特性 ,应用小波变换和波形分解可以解决误差影响和影响波的干扰的问题 ;应用特征提取可以解决周跳检测问题 ;应用快速小波变换可进行数据压缩     

差分法周跳探测与修复方法改进

在目前常用的周跳探测与修复方法基础上 ,提出了首先将观测资料按照观测历元不连续分成若干小弧段 ,然后利用差分法进行周跳探测 ,根据差分后周跳放大的特性判断周跳和野值 ,并确定其位置利用宽带组合和电离层组合的方法解算周跳大小。通过实例验证了其有效性。     

应用GPS观测青藏高原东北缘应力场变化

采用各向同性弹性地球模型推导了地面位移场速率与地壳内任意点应力场变化的边界积分关系,同时利用青藏高原东北缘1999～2001年观测的GPS资料对观测区地壳深度为5 km和25 km的主应力和最大剪应力进行了计算分析.结果表明,青藏高原东北缘的主应力变化主要集中在祁连山断裂、海原断裂等,在1920～1954年间历史上发生过多次震级为7.0～8.5级强震的断裂附近,并具有主应力变化沿断层走向分布、最大剪应力沿断层走向交替变化等特征.     

The SuperCOSMOS Sky Survey – III. Astrometry

In this, the third in a series of three papers concerning the SuperCOSMOS Sky Survey, we describe the astrometric properties of the data base. We describe the algorithms employed in the derivation of the astrometric parameters of the data, and demonstrate their accuracies by comparison with external data sets using the first release of data, the South Galactic Cap survey. We show that the celestial coordinates, which are tied to the International Celestial Reference Frame via the Tycho–2 reference catalogue, are accurate to better than ±0.2 arcsec at J , R ∼19,18 , rising to ±0.3 arcsec at J , R ∼22,21 , with positional-dependent systematic effects from bright to faint magnitudes at the ∼0.1-arcsec level. The proper motion measurements are shown to be accurate to typically ±10 mas yr⁻¹ at J , R ∼19,18 , rising to ±50 mas yr⁻¹ at J , R ∼22,21 , and are tied to zero using the extragalactic reference frame. We show that the zero-point errors in the proper motions are ≤1 mas yr⁻¹ for R >17 , and are no larger than ∼10 mas yr⁻¹ for R <17 mas yr⁻¹ .     

Parallax distortion by the weak microlensing effect

Parallax measurements allow distances to celestial objects to be determined. Together with measurements of their position on the celestial sphere, they give a full three-dimensional picture of the location of the objects relative to the observer. The distortion of the parallax value of a distant source affected by weak microlensing is considered. This means that the weak microlensing leads to distortion of the distance scale. The gravitational deflection causes a retrograde apparent motion of the image, which is revealed as a negative parallax. It is shown that the distortions may significantly change the parallax values when they amount to several microseconds of arc. In particular, at this level many measured values of parallaxes should be negative.     

Evidence of climatic warming in the southern Urals region derived from borehole temperatures and meteorological data

Thirty borehole temperature–depth profiles in the central and southern Urals, Russia were scrutinized for evidence of ground surface temperature histories. We explored two inversion schemes: a simple ramp inversion in which solutions are parameterized in terms of an onset time and magnitude of change and a more sophisticated functional space inverse algorithm in which the functional form of the solution is left unspecified. To enhance and potentially identify latitudinal differences in the ground surface temperature signal, we subdivided the data into three groups based on geographic proximity and simultaneously inverted the borehole temperature–depth logs. The simultaneous inversions highlighted 13 temperature–depth logs that could not both fit a common ground surface temperature history and a priori models within reasonable bounds. Our results confirm that this is an effective way to reduce site-specific noise from an ensemble of boreholes. Each inversion scheme gives comparable results indicating locally variable warming on the order of 1°C starting between 1800 and 1900 AD. Similarly surface air temperature records from 12 nearby meteorological stations exhibit locally variable warming also on the order of 1°C of warming during the 20th century. To explore the degree to which borehole temperatures and surface air temperature (SAT) time series are responding to the same signal, we average the SAT data into the same three groups and used these averages as a forcing function at the Earth's surface to generate synthetic transient temperature profiles. Root mean square (RMS) misfits between these synthetic temperature profiles and averaged temperature–depth profiles are low, suggesting that first-order curvature in borehole temperatures and variations in SAT records are correlated.     

Report of the International Astronomical Union Division I Working Group on Precession and the Ecliptic

The IAU Working Group on Precession and the Equinox looked at several solutions for replacing the precession part of the IAU 2000A precession–nutation model, which is not consistent with dynamical theory. These comparisons show that the (Capitaine et al., Astron. Astrophys., 412, 2003a) precession theory, P03, is both consistent with dynamical theory and the solution most compatible with the IAU 2000A nutation model. Thus, the working group recommends the adoption of the P03 precession theory for use with the IAU 2000A nutation. The two greatest sources of uncertainty in the precession theory are the rate of change of the Earth’s dynamical flattening, ΔJ₂, and the precession rates (i.e. the constants of integration used in deriving the precession). The combined uncertainties limit the accuracy in the precession theory to approximately 2 mas cent⁻². Given that there are difficulties with the traditional angles used to parameterize the precession, z_A, ζ_A, and θ_A, the working group has decided that the choice of parameters should be left to the user. We provide a consistent set of parameters that may be used with either the traditional rotation matrix, or those rotation matrices described in (Capitaine et al., Astron. Astrophys., 412, 2003a) and (Fukushima Astron. J., 126, 2003). We recommend that the ecliptic pole be explicitly defined by the mean orbital angular momentum vector of the Earth–Moon barycenter in the Barycentric Celestial Reference System (BCRS), and explicitly state that this definition is being used to avoid confusion with previous definitions of the ecliptic. Finally, we recommend that the terms precession of the equator and precession of the ecliptic replace the terms lunisolar precession and planetary precession, respectively.