用极移数据直接测算钱德勒摆动Q值

运用已经发展起来的共振激发和参数共振模型,在频率调制的基础上,对ＩＥＲＳ极移数据测算出钱德勒摆动的Ｑ值,得到平均值为６３,并以每百年０．８的速率在增大。这一结果与Ｑ的滞弹性ＰＲＥＭ模型理论值６９十分接近,表明共振激发模型和参数共振模型完全与钱德勒摆动滞弹理论相容。     

用极移数据直接测算钱德勒摆动Q值

运用已经发展起来的共振激发和参数共振模型,在频率调制的基础上,对IERS极移数据测算出钱德勒摆动的Q值,得到平均值为63,并以每百年0.8的速率在增大。这一结果与Q的滞弹性PREM模型理论值69十分接近,表明共振激发模型和参数共振模型完全与钱德勒摆动滞弹理论相容。     

地球自转形变的动力学效应

本文详细讨论了地球自转形变的动力学效应。地球自转形变导数地球引力位系数产生时变摄动并引起测站位置的变化。导出了由于地球自转引起的地球引力位系数变化的表达式，并给出了表示地球自转角速度矢量的三个参数m1、m2、m3与地球自转参数xp、yp、D间的关系。建议在高精度空间大地测量中要顾及由极潮引起的测站位置的变化。     

升降轨道ASAR雷达干涉揭示的巴姆地震（Mw6．5）3D同震形变场

通过对覆盖同一地区的升降轨道ASAR雷达数据进行两路差分干涉（D-INSAR）处理，得到了2003年12月26日伊朗巴姆（Bam）地震（Mw6．5）在两种成像几何状态下的视线向（LOS）同震形变场。利用Okada弹性形变模型对形变场进行正演计算，得到地震断层的几何参数及形变的南北向分量。对得到的干涉条纹图进行相位解绕处理，综合分析两种不同的雷达成像几何关系，并利用模型估计得到的南北向分量，计算了该地震造成的3D同震形变场，结果显示了典型的单断层右旋走滑活动特征。研究表明，合并升降轨道雷达数据，能够为形变模型提供有效约束，用一条简化的断层模型就可以解释升降两条轨能的干涉形变场，更为复杂的双断层模型是不必要的。     

汶川地震震后GNSS形变分析

采用国际通用高精度数据处理软件对1999—2015年汶川地震震后龙门山断裂两侧分布的109个GNSS连续站和区域站资料进行数据处理与分析,得到各测站点水平向坐标时间序列。考虑2013年芦山地震同震位移影响,通过时间序列分析,获取震后该区域109个测站在2010—2015年间累计震后形变场,断裂带上盘靠近震中的近场区域,水平向最大震后位移达到5~7 cm。基于震后黏弹性松弛模型,通过格网搜索,根据3种不同破裂模型分别反演龙门山断裂上盘中上地壳最佳弹性层厚度估值与下地壳/上地幔黏弹性层最佳黏滞系数,获取3种不同破裂模型在2010—2015年间震后形变观测值与模拟值分布。根据模型二反演得到的模型参数,推估汶川地震在未来几十年内对周边区域地壳形变影响,其中仅2018—2058年40年尺度累计震后位移最大能达到19 cm。     

板块内部层状负位错模型及其反演

提出了一种层状负位错模型，用于描述下层地壳拖带上层地壳的耦合作用所造成的上伏地壳变形，并通过模拟计算，给出了这种层状耦合模型的形变场特征，以及利用模拟形变数据反演层状负位错模型参数的结果。     

重力场长周期变化的观测与理论结果比较

利用Chandler摆动周期作为约束，估计了地幔滞弹性对地球二阶带谐响应系数κ、带谐位系数J2和卫星轨道升交点Ω的影响，理论预测的长周期潮汐参数被用来与人卫激光测距（SLR）观测结果进行比较。结果表明，SLR确定的长周期潮汐解已能检测到地幔滞弹性的影响。考虑地幔滞弹性和非平衡海洋潮汐效应后，理论预测的18．6a潮汐参数与SLR潮汐解基本相符。     

力模型误差对近地卫星定轨的影响分析

本文针对近地卫星所受摄动力的特点，详细分析了大气阻尼系数、大气模式、地球重力场模型、系统性误差等对定轨的影响。利用待估大气阻尼系数及其变率，可有效地消除大气模型差的影响，提高定轨精度。在得克萨斯大学空间研究中心编制的精密定轨软件ＵＴＯＰＩＡ８６的基础上，重新编制和调试出能有效处理近地卫星定轨的实用软件。     

GPS反射信号用于海冰识别与海面测高

针对全球定位系统反射信号干涉测量(GPS-IR)反演海面参数技术模型简单且应用单一的问题,该文提出引入表征信号振幅变化的阻尼系数来修正粗糙面对反射信号的影响,以非线性最小二乘算法对修正模型的参数进行识别.实验表明:在海冰识别上,阻尼系数的变化可清晰反映测站附近海冰状态的变化,证明了岸基GPS-IR识别海冰的可行性;在海面测高上,反演值与验潮站的实测值的相关系数达到0.98以上,相比经典频谱分析理论,均方根误差和平均绝对误差降低了约50％.     

由PREM地球模型计算的洛夫数

本文选取迄今为止最完善的地球模型ＰＲＥＭ,利用它所定义的地球内部基本间断面及所给物质密度的弹性等参数分布,通过求解弹性地幔运动方程,得到形变位移矢量,由此计算出二阶洛夫数。本文的理论值与观测结果符合良好。     

Nonlinear mechanism for polar motion with period of 7 months

 The 7 month-period (sub-Chandler) wobble is considered with respect to the nonlinear dynamical equation of polar motion. Starting with the frequency modulation of Chandler wobble (CW) in the model developed by introduction of damping from perturbed visco-elastic deformation, the rotation equation of the CW becomes a resonance model with a time-dependent parameter. According to evolution calculation, the parameter resonance model is essentially identical to the reality of CW observations. If the frequency of CW is modulated about 3 by visco-elastic deformation, then the amplitude of CW will be modulated by greater than 70. On the other hand, bifurcation may occur according to the nonlinear dynamical system of the parameter resonance model, i.e. a pair of solutions splitting from the main CW. One is the 7-month-period wobble and the other is a motion with a period of about 28 months. Although the latter is very weak, the 7-month-period wobble will be observed as the stability condition is satisfied. The maximum amplitude is about 22.89 mas and the average 12.65 mas. This is identical to what is observed in reality. Received: 7 April 2000 / Accepted: 5 March 2001     

Comparison of polar motion with oceanic and atmospheric angular momentum time series for 2-day to Chandler periods

During a 4-year period starting in July 1996 and using intervals ranging from 3 days to 4 years, four precise polar motion (PM) series have been compared to excitation by atmospheric angular momentum (AAM) augmented with oceanic angular momentum (OAM) data. The first three series (C03, C04 and Bulletin A) are multi-technique combinations generated by the International Earth Rotation and Reference Systems Service (IERS) and the fourth combined series (IGS00P02) is produced by the International GPS Service (IGS) using only GPS data. The IGS PM compared the best with the combined excitations of atmosphere and oceans (AAM+OAM) at all intervals, showing high overall correlation of 0.8–0.9. Even for the interval of only three days, the IGS PM gave a significant correlation of about 0.6. Moreover, during the interval of February 1999 – July 2000, which should be representative of the current precision of the IGS PM, a significant correlation (>0.4) extended to periods as short as 2.2 days and 2.5 days for the x_p and y_p PM components, respectively. When using the IERS Bulletin B (C04) PM and an interval of almost 6 years, starting in November 1994, the combined OAM+AAM accounted for practically all the annual, semi-annual and Chandler wobble (CW) PM signals. When only AAM was used, either the US National Centers for Environment Prediction reanalysis data, which were used throughout this study, or the Japanese Meteorological Agency data, two large and well-resolved amplitude peaks of about 0.1 mas/day, remained at the retrograde annual and CW periods.     

对钱德勒参数进行时变修正的CLS+AR模型在极移预测中的应用

根据LS模型参数中钱德勒项的时变性质,提出了基于钱德勒参数时变修正的CLS模型,并用CLS+AR对极移进行预测。实验结果表明,CLS+AR模型在极移预测精度上较LS+AR模型有较大改善。     

Atmospheric excitation of polar motion

The polar motion excited by the fluctuation of global atmospheric angular momentum (AAM) is investigated. Based on the global AAM data, numerical results demonstrate that the fluctuation of AAM can excite the seasonal wobbles (e.g., the 18-month wobble) and the Chandler wobble, which agree well with previous studies. In addition, by filtering the dominant low frequency components, some distinct polar wobbles corresponding to some great diurnal and semi-diurnal atmospheric tides are found.     

Atmospheric excitation of polar motion

The polar motion excited by the fluctuation of global atmospheric angular momentum (AAM) is investigated. Based on the global AAM data, numerical results demonstrate that the fluctuation of AAM can excite the seasonal wobbles (e.g., the 18-month wobble) and the Chandler wobble, which agree well with previous studies. In addition, by filtering the dominant low frequency components, some distinct polar wobbles corresponding to some great diurnal and semi-diurnal atmospheric tides are found.     

Polar motion modeling,analysis, and prediction with time dependent harmonic coefficients

A time dependent amplitude model was proposed for the analysis and prediction of polar motion time series. The formulation was implemented to analyze part of the new combined solution, EOP (IERS) C 04, daily polar motion time series of 14 years length using a statistical model with first order autoregressive disturbances. A new solution approach, where the serial correlations of the disturbances are eliminated by sequentially differencing the measurements, was used to estimate the model parameters using weighted least squares. The new model parsimoniously represents the 14-year time series with 0.5 mas rms fit, close to the reported 0.1 mas observed pole position precisions for the x and y components. The model can also predict 6 months into the future with less than 4 mas rms prediction error for both polar motion components, and down to sub mas for one-step ahead prediction as validated using a set of daily time series data that are not used in the estimation. This study is dedicated to the memory of Prof. Urho Uotila (1923–2006) whose teaching of “Adjustment Computations” over the years influenced so much, so many of us who had the privilege of being his students.     

利用LS_AR模型对极移参数的中长期预报

为了满足深空探测器自主导航定位对极移参数中长期预报的需求,阐述了基于LS_AR模型的极移参数中长期预报和精度评定的原理,提出了4种改进方案对LS_AR模型的构建进行优化,并利用IERS提供的1990～1996年的极移参数的时间序列检验4种优化方案,得到了最优的预报模型,在400 d跨度上,其预报结果的平均绝对误差比未优化的模型小3 mas左右。     

A study of gravity variations caused by polar motion using superconducting gravimeter data from the GGP network

Long-term continuous gravity observations, recorded at five superconducting gravimeter (SG) stations in the Global Geodynamic Project (GGP) network, as well as data on orientation variations in the Earths rotation axis (i.e. polar motion), have been used to investigate the characteristics of gravity variations on the Earths surface caused by polar motion. All the SG gravity data sets were pre-processed using identical techniques to remove the luni-solar gravity tides, the long-term trends of the instrumental drift, and the effects of atmospheric pressure. The analysis indicates that the spectral peaks, related to the Chandler and annual wobbles, were identified in both the power and product spectral density estimates. The magnitude of gravity variations, as well as the gravimetric amplitude factor associated with the Chandler wobble, changed significantly at different SG stations and during different observation periods. However, when all the SG observations at these five sites were combined, the gravimetric parameters of the Chandler wobble were retrieved accurately: 1.1613 ± 0.0737 for the amplitude factor and –1°.30 ± 1°.33 for the phase difference. The value of the estimated amplitude factor is in agreement with that predicted theoretically for the zonal tides of an elastic Earth model.     

New Solution for the Earth’s Free Wobble and Its Geophysical Implications

In this paper, the theory of the free wobble of the triaxial Earth is developed and new conclusions are drawn: the Euler period should be actually expressed by the first kind of complete elliptic integral; the trace of the free polar motion is elliptic and the orientations of its semi-minor and major axes are approximately parallel to the Earth’s principal axes A and B, respectively. In addition, the present theory shows that there is a mechanism of frequency-amplitude modulation in the Chandler wobble, whi...