天文测时测纬的一种短时间异常与台站附近强震的对应关系

本文具体处理了自1966年以来,周围曾经发生过强震的北京,天津,东京（日本）,水泽（日本）及贝尔格莱德（南斯拉夫）等五个台站的测时测纬数据的残差,发现大震前的三、五个月内测时残差均出现一对符号相反的异常,而测纬残差则是一个孤立的单峰.这种震前的异常特征在无震的年份和无强震台站的资料（为了对比,又选了世界上观测稳定、质量较好的六架照相天顶筒,三架天顶仪和两架等高仪的资料）中极为罕见.对比关系较为密切,难于用偶然误差和其它影响来解释.各地测时与测纬残差异常特征的差异可能系地质构造不同所致.上述结论对进一步探讨地球自转的不均匀性的原因,提高时纬观测精度,揭示地下物质的运动形式及震源机制的本质等方面提供了新的线索和依据.     

地球自转速度变化在地震预报中的地位初探

本文计算了地球自转速度季节性变化与全球强震季节性频度曲线的相关性,相关系数与震级有关,並发现强震频度的纬度分布与离心力大小有关,从而可能找到地球自转触发地震的某些依据。並用地下水异常、激光锁相应变仪测到的地震前兆,分析了地球自转短周期与地震的关系。     

Earth Rotation Parameters 1899.7–:1992.0 After Reanalysis Within The Hipparcos Frame

Earth rotation parameters (ERP) in the interval 1899.7–1992.0 are obtained from re-analysis of the observed latitude/universal time variations by optical astrometry. Hipparcos Catalogue is used to define the celestial reference frame, within which the ERP are described, with special care devoted to 'problematic' double and/or multiple stars. The terrestrial reference frame is defined by the adopted latitudes/longitudes of participating instruments and their secular motions as given by the NUVEL-1 NNR model of plate motions, and it is chosen to be very close to the International Terrestrial Reference Frame (ITRF). More than four million observations made with 48 different instruments at 31 observatories, located all over the world, are utilized to determine polar motion, celestial pole offsets and (after 1956) universal time UT1, all at 5-day intervals. Along with these parameters, the combinations of Love and Shida numbers, governing the tidal variations of the local verticals at individual observatories, are also determined. The analysis of the results covering almost a century, namely the long-periodic polar motion and length-of-day changes, is presented.     

2004年印尼苏门答腊MW9.0地震前地震活动增强现象及其与地球自转的关系

选取2004年12月26日印尼苏门答腊M_W9.0地震的余震分布空间范围为研究区域,该区域位于缅甸小板块中部,纬度范围为2°N~15°N.根据USGS地震目录,研究了2004年印尼苏门答腊M_W9.0地震前地震频次的时空变化,并分析了地震活动与地球自转的关系.结果显示:1)在余震区域范围内,地震频次自2000年开始出现了明显的上升趋势变化,反映了地震活动的增强过程,临近地震发生前,地震频次有一定幅度的回落;2)通过地震频次空间扫描,发现震前在震中附近区域地震活动增强显著;3)在1990~2004年的15年中,单位时间内发生的地震数在地球自转季节性变化减速时段与加速时段之比值r_da在震前出现明显的上升过程:1998年~2002年r_da缓慢上升,2002年底r_da快速上升达最高值,之后,r_da呈现下降趋势,直至苏门答腊M_W9.0地震发生;4)2000年1月~2004年11月间,在纬度1°N~7°N的空间范围内的震中附近地区, 5.5级以上地震都发生在地球自转季节性变化减速时段内;5)地震活动累积年频次相对于其长期平均值的变化量ΔN与地球自转"十年起伏"变化关系表现为地球自转减速时ΔN上升,地震活动增强,加速时ΔN下降,地震活动减弱.这些现象并非偶然,可能预示着地球自转减速与苏门答腊M_W9.0地震之间存在着某种联系,这些结果对地震孕育过程的认识和地震预测研究具有一定的价值.     

中低纬度电离层对太阳风状态的响应

本文讨论了行星际磁场B₂分量变化时内磁层和中低纬度电离层的响应.指出B₂变化引起的磁层大尺度对流电场的变化在一定条件下有可能透入内磁层,并沿磁力线映射到中低纬度电离层,在那里产生电场和电流体系,从而使S_q电流体系发生畸变,并在地面磁场中反映出来.数值计算表明,当△B₂<0时,S_q电流体系的焦点向东和向高纬移动,地面磁场会观测到数伽马的变化.这就为中低纬地磁观测诊断磁层和太阳风状态提供了一种可能性.此外,本文还用上述物理过程解释了赤道地区一些高空物理现象,如B₂倒转时电离层漂移速度的变化,赤道磁场异常以及赤道q型偶现E层的消失等等.     

Semi-diurnal tidal effects in P.Z.T. observations

Deflexions of the vertical due to changes in luni-solar gravitational attraction can be determined by astronomical observations and the measures are of great importance in the derivation of values for parameters defining elastic properties of the earth. The primary purpose of the observations is the determination of time or latitude, or both, and the systematic deflexions of the vertical are small in comparison with the scatter of the observations so that regular observations made over a number of years are required for the detection of these effects. Eleven years of time and latitude observations made with the Herstmonceux Photographic Zenith Tube have been separately analysed and the two solutions reveal with great clarity the magnitudes and phases of the variations associated with tidal activity. The results are complicated because at Herstmonceux the deflexions arising directly from changes in the gravitational attraction of the moon have associated with them deflexions caused by changes in the attraction of the proximate waters of the English Channel which is itself affected by complex tidal activity.     

Seismic excitation of the polar motion, 1977–1993

The mass redistribution in the earth as a result of an earthquake faulting changes the earth's inertia tensor, and hence its rotation. Using the complete formulae developed byChao andGross (1987) based on the normal mode theory, we calculated the earthquake-induced polar motion excitation for the largest 11,015 earthquakes that occurred during 1977.0–1993.6. The seismic excitations in this period are found to be two orders of magnitude below the detection threshold even with today's high precision earth rotation measurements. However, it was calculated that an earthquake of only one tenth the size of the great 1960 Chile event, if happened today, could be comfortably detected in polar motion observations. Furthermore, collectively these seismic excitations have a strong statistical tendency to nudge the pole towards 140°E, away from the actually observed polar drift direction. This non-random behavior, similarly found in other earthquake-induced changes in earth rotation and low-degree gravitational field byChao andGross (1987), manifests some geodynamic behavior yet to be explored.     

Analyses of zonal atmospheric excitation functions and their correlation with polar motion excitation functions

The atmospheric influence on the Earths, rotation can be described by the effective atmospheric angular momentum (EAAM) functions. In this study we focus on the analysis of short period variations of the equatorial components of the zonal EAAM excitation functions ₁ and ₂ and their influence on similar variations of polar motion. The global objective analysis data of the Japanese Meteorological Agency for the period 1986–1992 were used to compute the EAAM excitation functions in different latitude belts. Time- and latitude-variable amplitude spectra of variations of these functions with periods shorter than 150 days, containing pressure, pressure with the inverted barometric correction, and wind terms were computed. The spectra show distinct latitude and time variations of the prograde and retrograde oscillations which reach their maxima mainly in mid-latitudes. Prograde and retrograde oscillations with periods of about 40–60 days and about 110–120 days are seen in the spectra of pressure terms of the equatorial components of the zonal EAAM excitation functions. Additionally, correlation coefficients and cross-spectra between variations of the geodetic polar motion and equatorial components of the zonal EAAM excitation functions were computed to identify the latitude belts of the globe over which atmospheric circulation changes are correlated mostly with short period variations of the polar motion excitation functions. The correlation coefficients vary in time and latitude and reach maximum values in the northern latitudes from 50°N to 60°N. In the cross-spectra between the polar motion excitation functions and pressure terms of the zonal EAAM excitation functions there are peaks of common prograde oscillations with the periods around 20, 30, 40–50, 60 and 80–150 days and of common retrograde oscillations around 20, 30, 40 and 50–70 days.Paper presented at the IERS Workshop in Paris, March 1994     

地球自转减速对2008年汶川MS8.0地震的作用

本文将地球看成是一个质量分布均匀的球对称弹性球体, 给出了地球以匀角速率旋转情况下产生的应力场的解析表达式, 计算了地球匀速自转时在2008年汶川M_S8.0地震发震断层面上产生的应力, 分析了地球自转速率变化对汶川地震发震断层的影响。 结果表明, 地球自转减速有利于汶川地震发生。 如果地球自转减速对汶川地震的发生起到了触发作用, 那么在震前震区或许会出现中小地震发生受到地球自转减速触发的现象。 分析表明, 从2006年起汶川地震震中附近地区M_L≥2.7地震受到地球自转减速的触发, 这从另一个侧面说明了地球自转减速对汶川地震的发生具有促进作用。     

内核地震波速各向异性的成因

地球内核是轴对称各向异性的,其对称轴与地球的极轴之间有１１°左右的夹角,本 文根据地球内核相对于外部地球有差异转动这一观测结果,利用晶体生长理论,对内核地震波 速度各向异性的成因进行了探讨．当从熔融状态结晶时,晶体的生长速度与晶体和熔融态之 间相对运动的线速度成正比涸此当固态内核在液态外核中生长时,沿赤道方向的生长速度比 两极方向快．在万有引力场的作用下内核始终保持近似球形,生长速度较快的赤道附近的物 质会向两极区域流动,形成轴对称的流变场。这一轴对称的流变场伴随着轴对称的应力场,使 得构成地球内核的ｈｃｐ型铁晶体的ｃ轴沿着内核自转轴的方向排列,导致观测到的地球内核地 震波速度各向异性。作为推论,内核相对于外部地球可能同时存在着进动和章动。     

Secular polar motion,crustal movements,and international latitude service observations

Summary The latitude variations observed by optical astrometry at five International Latitude Service (ILS) stations since the beginning of the century are analyzed to estimate the secular motion of the pole and the correction of the proper motions of the observed stars. The geophysical model of tectonic plate motions NUVEL-1 is used to correct the observed latitude drifts, thus referring the derived secular polar motion to the "no net rotation" terrestrial reference system. It is shown that one of the ILS stations, Ukiah (near the U.S. west coast), very probably exhibits an anomalous drift which differs significantly from the motion of both the Pacific and North American plates.     

Models of solar differential rotation

Abstract

Models of a differentially rotating compressible convection zone are calculated, considering the inertial forces in the poloidal components of the equations of motion. Two driving mechanisms have been considered: latitude dependent heat transport and anisotropic viscosity. In the former case a meridional circulation is induced initially which in turn generates differential rotation, whereas in the latter case differential rotation is directly driven by the anisotropic viscosity, and the meridional circulation is a secondary effect.

In the case of anisotropic viscosity the choice of boundary conditions has a big influence on the results: depending on whether or not the conditions of vanishing pressure perturbation are imposed at the bottom of the convection zone, one obtains differential rotation with a fast (≥ 10 ms^?1) or a slow (～ 1 ms^?1) circulation. In the latter case the rotation law is mainly a function of radius and the rotation rate increases inwards if the viscosity is larger in radial direction than in the horizontal directions.

The models with latitude dependent heat transport exhibit a strong dependence on the Prandtl number. For values of the Prandtl number less than 0.2 the pole-equator temperature difference and the surface velocity of the meridional circulation are compatible with observations. For sufficiently small values of the Prandtl number the convection zone becomes globally unstable like a layer of fluid for which the critical Rayleigh number is exceeded.     

天文观测极移运动周期变化的原因解析

宋贯一（1991,2006）相继发现日-地之间存在一种奇特的能量（动量）相耦合的自然现象:当太阳辐射光压作用于地表之后,地球表面物质的特殊物理性质会自然地把太阳辐射光压立刻分解为两部分,即P₁和P₂.其中P₁为全球各纬度带内单位海洋和陆地表面接收到的等值光压（P₁为随时间和纬度而变化的变量）.在一个回归年的时段内,由P₁在北、南两半球上的不平衡分布激发自转轴摆动,引起转动惯量的变化（该变化是目前所了解到的地球转动惯量的最大变化）并产生极移,该项极移运动的周期为12个月左右.由于P₁在北、南半球上的分布相对赤道是基本对称的和规则的,P₁对自转轴摆动的激发可称谓“规则性激发”;P₂为全球各纬度带内单位陆地和海洋表面接收到的光压值之差（P₂为随时间和纬度而变化的变量）.在一个回归年的时段内,由 P₂在北、南两半球上的不平衡分布激发自转轴摆动而产生的极移运动周期为14个月附近.由于P₂在北、南两半球上的分布相对赤道是不对称和不规则的, P₂对自转轴摆动的激发可称谓“不规则性激发”.这种天然存在的力源,恰恰是近百年来世界各国研究地球自转的地球物理学家渴望寻找的那种既能激发自转轴产生自由章动、又使其摆动含有两个不同周期（12个月和14个月）的激发力源.正是这一奇特自然现象的发现,才使长期以来困扰地球物理学领域内的极移、极移所包含两个周期的涨落变化及由此引起的地球自转速度变化等自然之谜得以破解成为可能. 本文作者仅对天文观测的极移运动周期及其极移运动所包含的两个周期分量在一定的范围内变化的成因作出了详细的解析,并得出如下结论： （1）极移运动主要是由太阳光压P₁ 和P₂共同激发引起的.天文观测的极移摆动周期的涨落变化是太阳光压（P₁＋P₂）激发自转轴摆动过程中,在空间上自转轴的摆动中心相对自转轴中心（地心）移动造成的,涨落范围在395～403±2天之间,即天文观测的极移运动的实际计算周期应在13.0～13.3个月之间变化.（2）极移所含的周年期摆动是由太阳光压P₁激发的.天文观测的周年期摆动周期涨落很小,变化于365.24～365.53天之间.在一个回归年内,由于日-地间距离的变化,使地球表面接受到的太阳辐射光能产生微小差异则是造成观测的周年摆动周期稍有拖长的原因.（3）极移所含的钱德勒摆动周期是由太阳光压P₂激发的.天文观测的钱德勒周期涨落较大,变化于426～437±2天之间,即实际计算周期应在14.0～14.4个月之间变化.观测的钱德勒摆动周期的变化是太阳光压P₂激发自转轴摆动过程中,在空间上自转轴的摆动中心相对自转轴中心（地心）移动造成的.上述的定量解析数据均得到实际观测资料的验证,为极移光压成因理论的正确性提供了具体详实的证据.     

Time variations in the mechanical characteristics of local crustal segments according to seismic observations

The results of the seismic observations made with two different experimental setups are presented. In the first case, the signals produced by underground nuclear explosions at the Semipalatinsk Test Site were measured on a linear profile, which allowed one to definitely outline the areas where the mechanical properties of rocks experienced considerable time variations. In the second case, the waves excited by the open-pit mine blasts recorded at a small-aperture seismic array at the Mikhnevo Geophysical Station (Institute of Geosphere Dynamics, Russian Academy of Sciences) on the East European Platform favored the estimation of variations in the integral characteristics of the seismic path. Measurements in aseismic regions characterized by diverse geological structure and different tectonic conditions revealed similar effects of the strong dependency of seismic parameters on the time of explosions. Here, the variations experienced by the maximum amplitudes of oscillations and irrelevant to seasonal changes or local conditions reached a factor of two. The generic periods of these variations including the distinct annual rhythm are probably the fragments of a lower-frequency process. The obtained results suggest that these variations are due to changes in the stressstrain state of active fault zones, which, in turn, can be associated with the macroscale motion of large blocks triggered by tidal strains, tectonic forces and, possibly, variations in the rate of the Earth’s rotation.     

Complex demodulation in VLBI estimation of high frequency Earth rotation components

The spectrum of high frequency Earth rotation variations contains strong harmonic signal components mainly excited by ocean tides along with much weaker non-harmonic fluctuations driven by irregular processes like the diurnal thermal tides in the atmosphere and oceans. In order to properly investigate non-harmonic phenomena a representation in time domain is inevitable. We present a method, operating in time domain, which is easily applicable within Earth rotation estimation from Very Long Baseline Interferometry (VLBI). It enables the determination of diurnal and subdiurnal variations, and is still effective with merely diurnal parameter sampling. The features of complex demodulation are used in an extended parameterization of polar motion and universal time which was implemented into a dedicated version of the Vienna VLBI Software VieVS. The functionality of the approach was evaluated by comparing amplitudes and phases of harmonic variations at tidal periods (diurnal/semidiurnal), derived from demodulated Earth rotation parameters (ERP), estimated from hourly resolved VLBI ERP time series and taken from a recently published VLBI ERP model to the terms of the conventional model for ocean tidal effects in Earth rotation recommended by the International Earth Rotation and Reference System Service (IERS). The three sets of tidal terms derived from VLBI observations extensively agree among each other within the three-sigma level of the demodulation approach, which is below 6 μas for polar motion and universal time. They also coincide in terms of differences to the IERS model, where significant deviations primarily for several major tidal terms are apparent. An additional spectral analysis of the as well estimated demodulated ERP series of the ter- and quarterdiurnal frequency bands did not reveal any significant signal structure. The complex demodulation applied in VLBI parameter estimation could be demonstrated a suitable procedure for the reliable reproduction of high frequency Earth rotation components and thus represents a qualified tool for future studies of irregular geophysical signals in ERP measured by space geodetic techniques.     

地球自转的天文地球动力学效应研究(Ⅰ)--引力位和重力的变化

详细讨论了地球自转的天文地球动力学效应，地极移动和日长变化导致地球引力位系数产生时变特性并引起重力的摄动，根据理论力学的基本概念，导出了由于地球自转变化引起的地球引力位系数变化、重力摄动、垂线偏差和地球形变的表达式，并定量地研究了极移和日长变化对测站重力观测值和地球形变的影响，建议在高精度的空间大地测量中要顾及到地球自转变化引起的一系列效应。     

根据激光雷达观测结果研究中国武汉地区钠层的分布

通过大量钠层荧光激光雷达观测数据的分析研究，给出我国武汉地区钠层平均分布形态的基本特征. 结合钠层模型分析，指出大气微量元素分布的纬度变化可能是导致不同纬度地区钠层平均分布形态差异的主要原因,同时给出武汉地区随季节变化的钠层平均分布形态. 对武汉地区钠层分布形态时间变化的分析，发现其短期变化与重力波活动及夜间变化与潮汐波活动相关.     

纬度变化和地极运动

远在1687年,牛顿已说到地球自转轴的运动,欧拉在1790年指出:由于地球转动的主要动力矩的方向,并不和地球惯性椭球的短轴方向相重合,将使地球的自转轴发生沿着一个圆锥面的摆动.假定地球为一个刚体时,可以从地球惯性椭球的扁率（或称力学椭率,或动力学扁率）算出地球自转轴运动的周期.     

川滇地区地震活动与地球自转速率变化之间的关系

利用川滇地区1962—2016年间发生的MS≥6.0地震资料,在去除余震的情况下,分析了其发生与地球自转速率变化极值点之间的关系。得到以下结果:川滇地区MS≥6.0地震的发生与地球自转速率变化极值点具有显著关系,对于不同震级范围的地震,其与不同周期的地球自转速率变化有关。对于MS≥7.0的地震,有90%发生在地球自转速率季节性变化极值点前14天至后37天时间段内;对于6.0≤MS≤6.4的地震,有80%发生在地球自转速率短周期变化极值点前68小时至后30小时的时间段内;对于6.5≤MS≤6.9的地震,有87.5%发生在地球自转速率短周期变化极值点前36小时至后64小时的时间段内。地震发生在特定时间段的显著性检验结果表明,川滇地区MS≥6.0地震与地球自转速率变化极值点的关系都可以在α=0.05的显著性水平下通过显著性检验。这个结果表明,在地球自转速率发生转折的期间容易触发地震,对川滇地区地震发生时间预测具有参考意义。