用Hamilton方法计算滞弹地球自转的潮汐变化

在日月引潮力势作用下地球产生弹性形变．地幔粘滞性子致这个形变对于引潮力滞后,成为引起地球自转长期减慢的原因之一．地幔滞弹性也使有效洛夫数k增加,并使自转变化的周期项位相滞后,即产生反常位相项．本文首先用Hamilton方法计算了地球的形变．然后考虑到地幔的滞弹性,计算了在日月引潮力作用下的地球自转长期减慢和滞弹性对周期（带谐）变化的影响．     

IERS1996规范在参考系方面的改进

对ＩＥＲＳ１９９６规范中有关参考系方面作了简单而系统的介绍，重点叙述了与ＩＥＲＳ１９９２标准相比ＩＥＲＳ１９９６规范在参考系方面的主要改进：天球参考架中的基本源从５７颗增加到２３６颗；动力学参考架采用ＪＰＬ ＤＥ４０３／ＬＥ４０３历表代替ＤＥ２００／ＬＥ２００历表；采用ＮＵＶＥＬ ＮＮＲ－１Ａ板块运动模型代替ＮＵＶＥＬ ＮＮＲ－１模型；变更了９个基本常数值；给出了天极坐标的观测和理论间差的经验模型     

IERS天球参考架的指向维持

本文分析了河外射电源坐标系不统一问题对ＩＥＲＳ天球参考架指向维持的影响。由于１９８８年至１９９１年各年度的ＩＥＲＳ河外射电源综合表中基本源与非基本源坐标系不统一的程度各不相同，所以不能在每年的指向上加入相同的改正量。     

星族综合模型与星系的颜色演化

本文利用Ｂｅｒｔｅｌｉ等人新算出的等龄线，以及Ｇｕｎｎ和Ｓｔｒｙｋｅｒ观测的恒星光谱（３１３０－１０８００）构造了一个星族综合模型．该模型能预计金属丰度在Ｚ＝０．０００４－０．０５，年龄在１０７－２×１０１０年范围内的简单星族（ＳＳＰ）的光谱能量分布（ＳＥＤ）．利用ＳＳＰ，通过对不同时刻恒星形成率的卷积能够得到星系的ＳＥＤ的演化．我们计算出ＳＳＰ以及星系在ＵＢＶＲＩ系统和ＢＡＴＣ１５色中带测光系统中的颜色．计算出的从Ｅ／ＳＯ到Ｓｄｍ各类星系的ＵＢＶＲＩ颜色在年龄为１２Ｇｙｒ时都和观测符合较好．因此，我们的模型为正在进行的关于星系ＳＥＤ的多色测光研究提供了理论上的探索．     

1988年至1995年IERS天球参考架的指向维持

本文分析了１９９２年至１９９５年期间国际地球自转服务（ＩＥＲＳ）天球参考架的指向维持情况，结合１９８８年至１９９２年的指向维持情况指出，该参考架在１９８８年至１９９５年期间历年的指向均存在显著变化，赤经零点和天极的最大年度漂移约为０．４ｍａｓ。至１９９５年，ＩＥＲＳ天球参考架的赤经零点和天极与各自预期位置的偏差不于０．４ｍａｓ。１９９４年和１９９５年的ＩＥＲＳ天球参考架基本上维持了１９９３年的ＩＥＲＳ天球参考指向，但是其实际维持精度约为０．０５ｍａｓ，并非０．００５ｍａｓ。０．００５ｍａｓ的维持精度只是数学上的平均效果，相当一部分基本源在相邻年度的ＩＥＲＳ天球参考架中的坐标差大于０．５ｍａｓ，这说明只有采用恰当的消除局部相对形变的方法，才能将天球参考架的指向真正维持在较高水平。     

由PREM模型参数计算地球自转的周期变化

弹性地球在日月引潮力势作用下的形变引起其转动惯量的改变，从而导致地球自转速率的变化．本文利用PREM地球模型所给的物质密度和弹性等参数分布．计算日月引潮力势产生的地球形变附加势，进而计算转动惯量的变化．最后得到一系列包含周期同引潮势带谐项、振幅大于1微秒的自转速率周期变化系数．     

1988至1992各年度IERS综合河外射电源表的重新归算

鉴于ＩＥＲＳ综合河外射电源表在建立ＩＡＵ协议天球参考架方面的重要意义和它们的不足之处，本文重新组建了１９８８至１９９２各年度的综合河外射电源表，称之为Ｃａｔ系列表。并依据综合表的组建原理，对重新计算得到的Ｃａｔ系列表进行了检验，结果表明Ｃａｔ系列表克服了ＩＥＲＳ系列表中坐标系不统一的缺点，指向的定义及维持亦较之ＩＥＲＳ系列表有较大程度地改观。比较结果显示，ＩＥＲＳ系列表中部分源的坐标存在毫角秒量级的偏差，各年度综合表的指向与ＩＥＲＳ预期指向亦有不同程度的偏离，个别轴的最大指向偏差至０．６ｍａｓ。同时指出，１９９３年ＩＥＲＳ综合表的指向亦有待讨论。     

1988至1992各年度IERS综合河外射电源表的重新归算

鉴于ＩＥＲＳ综合河外射电源表在建立ＩＡＵ协议天球参考架方面的重要意义和它们的不足之处，本文重新组建了１９８８至１９９２各年度的综合河外射电源表，称之为Ｃａｔ称列表，并依据综合表的组建原理，对重新计算得到的Ｃａｔ系列表进行了检验，结果表明Ｃａｔ系列表克服了ＩＥＲＳ系列表中坐标系不统一的缺点，指向的定义及维持亦较之ＩＥＲＳ系列表有较大程度地改观。比较结果显示，ＩＥＲＳ系列表中部分源的坐标存在毫角秒量级     

Starlette卫星精密定轨中大气模型效度的比较

为满足大气模型比较的需要，把ＣＩＲＡ－１９８６大气模式（简称ＣＩＲＡ８６）引入大气阻力选择模块，对上海天文台人造卫星精密定轨软件ＳＨＯＲＤＥ进行大气模式的扩充。利用多种大气模式分别对Ｓｔａｒｌｅｔｔｅ卫星的观测资料进行解算，结果表明：ＳＨＯＲＤＥ中指数大气模式更适合对Ｓｔａｒｌｅｔｔｅ卫星进行精密定轨研究。     

活动星系核黑洞发电机研究进展

类星体被发现后,黑洞发电机研究成为阐述活动星系核（ＡＧＮ）的产能机制及其宇宙学演化的重要内容之一,主要涉及发电机物理和发电机输出特性两个方面,从广义相对论电动力学（ＧＲＥＤ）发电机理论和“膜模式”出发,在太阳发电机研究中引入引力磁项,形成和αωΩ黑洞发电机理论。发电机输出包括吸积盘辐射和喷流传播。辐射模型在于经典吸积模型的相对论扩展,已能说明ＳｇｒＡ＾＊的宽发射谱,燕可能根据谱特征确定发电机的构型     

A series expansion of the solid Earth tide effect on geopotential

We develop analytical series representing the main part of corrections to the geopotential coefficients caused by the solid Earth tides, where Love numbers are assumed to be frequency-independent. The series are compact, precise and valid over 1800 A.D.–2200 A.D. The maximum difference between the corrections given by the analytical series and their numerical values, obtained with use of the DE/LE-423 planetary/lunar ephemerides, does not exceed $0.7\times 10^{-12}$ . A new algorithm is proposed for calculating amplitudes of the additional variations of the geopotential coefficients for frequency dependence of Love numbers. It uses the representation of the Earth tide-generating potential in the standard HW95 format and takes into account the phase of tidal waves. Corrections of up to $2\times 10^{-12}$ to the published by the IERS Conventions (2010) amplitudes of the additional variations of the geopotential coefficients are suggested. Examples of use of the obtained series in analytical theories of motion of low-altitude STARLETTE and high-altitude ETALON-1 satellites are given.     

Determining parameters of Moon’s orbital and rotational motion from LLR observations using GRAIL and IERS-recommended models

The aim of this work is to combine the model of orbital and rotational motion of the Moon developed for DE430 with up-to-date astronomical, geodynamical, and geo- and selenophysical models. The parameters of the orbit and physical libration are determined in this work from lunar laser ranging (LLR) observations made at different observatories in 1970–2013. Parameters of other models are taken from solutions that were obtained independently from LLR. A new implementation of the DE430 lunar model, including the liquid core equations, was done within the EPM ephemeris. The postfit residuals of LLR observations make evident that the terrestrial models and solutions recommended by the IERS Conventions are compatible with the lunar theory. That includes: EGM2008 gravitational potential with conventional corrections and variations from solid and ocean tides; displacement of stations due to solid and ocean loading tides; and precession-nutation model. Usage of these models in the solution for LLR observations has allowed us to reduce the number of parameters to be fit. The fixed model of tidal variations of the geopotential has resulted in a lesser value of Moon’s extra eccentricity rate, as compared to the original DE430 model with two fit parameters. A mixed model of lunar gravitational potential was used, with some coefficients determined from LLR observations, and other taken from the GL660b solution obtained from the GRAIL spacecraft mission. Solutions obtain accurate positions for the ranging stations and the five retroreflectors. Station motion is derived for sites with long data spans. Dissipation is detected at the lunar fluid core-solid mantle boundary demonstrating that a fluid core is present. Tidal dissipation is strong at both Earth and Moon. Consequently, the lunar semimajor axis is expanding by 38.20 mm/yr, the tidal acceleration in mean longitude is \(-25.90 {{}^{\prime \prime }}/\mathrm{cy}^2\), and the eccentricity is increasing by \(1.48\times 10^{-11}\) each year.     

国际地球自转服务(IERS)评介

回顾了建立国际地球自转服务(IERS)的历史过程,简述了IERS机构各部门的设置及其作用。着重介绍了IERS重建的理由和过程以及新IERS机构的组成。最后给出IERS成果对现有科学研究的作用和意义。     

Tidal effects in differentiated viscoelastic bodies: a numerical approach

The majority of confirmed terrestrial exoplanets orbits close to their host stars and their evolution was likely altered by tidal interaction. Nevertheless, due to their viscoelastic properties on the tidal frequencies, their response cannot be described exactly by standardly employed constant-lag models. We therefore introduce a tidal model based on the numerical evaluation of a continuum mechanics problem describing the deformation of viscoelastic (Maxwell or Andrade) planetary mantles subjected to external force. We apply the method on a model Earth-size planet orbiting a low-mass star and study the effect of the orbital eccentricity, the mantle viscosity and the chosen rheology on the tidal dissipation, the complex Love numbers and the tidal torque. The number of stable spin states (i.e., zero tidal torque) grows with increasing mantle viscosity, similarly to the analytical model of Correia et al. (Astron Astrophys 571:A50, 2014) for homogeneous bodies. This behavior is only slightly influenced by the rheology used. Similarly, the Love numbers do not distinctly depend on the considered rheological model. The increase in viscosity affects the amplitude of their variations. The tidal heating described by the Maxwell rheology attains local minima associated with low spin-orbit resonances, with depth and shape depending on both the eccentricity and the viscosity. For the Andrade rheology, the minima at low resonances are very shallow and the tidal heating for all viscosities resembles a “fluid limit.” The tidal heating is the quantity influenced the most by the rheology, having thus possible impact on the internal thermal evolution.     

The interior structure of Mercury and its core sulfur content

Depth-dependent interior structure models of Mercury are calculated for several plausible chemical compositions of the core and of the mantle. For those models, we compute the associated libration amplitude, obliquity, tidal deformation, and tidal changes in the external potential. In particular we study the relation between the interior structure parameters for five different mantle mineralogies and two different temperature profiles together with two extreme crust density values. We investigate the influence of the core light element concentration, temperature, and melting law on core state and inner and outer core size. We show that a sulfur concentration above 10 wt% is unlikely if the temperature at the core-mantle boundary is above 1850 K and the silicate shell at least 240 km thick. The interior models can only have an inner core if the sulfur weight fraction is below 5 wt% for core-mantle boundary temperature in the 1850-2200 K range. Within our modeling hypotheses, we show that with the expected precision on the moment of inertia the core size can be estimated to a precision of about 50 km and the core sulfur concentration with an error of about 2 wt%. This uncertainty can only be reduced when more information on the mantle mineralogy of Mercury becomes available. However, we show that the uncertainty on the core size estimation can be greatly reduced, to about 25 km, if tidal surface displacements and tidal variations in the external potential are considered.     

Accurate analytical nutation series

In this letter we present the first accurate analytical nutation series, deduced from the Hamiltonian theory by the authors. They provide the highest accuracy ever obtained by any analytical nutation series, since the deviation in CEP (celestial ephemeris pole) offsets with respect to that of the IERS Conventions 1996 is kept below 1 mas in the time domain, in spite of still lacking ocanic corrections.     

Toward the construction of a new precession-nutation theory of nonrigid Earth

The accuracy of the rigid Earth solution SMART97 is 2?μ as over the time interval (1968, 2023), accuracy showed by the comparison with a numerical integration using the positions of the Moon, the Sun, and the planets given by DE403. To obtain a nonrigid Earth solution, we use the transfer function of Mathews et al. (2000) and , to keep the precision of our rigid Earth solution in the computation of the geophysical effects, we apply this transfer function to the Earth's angular velocity vector in order to avoid the inherent approximations of the classical methods. Moreover the perturbations of the third component of the angular velocity vector are taken into account. Lastly, we take into account, in an iterative process, the second order perturbations due to the geophysical effects. The results are compared with the Herring solution (1996) published in the IERS Conventions.     

A Hamiltonian theory for an elastic earth: First order analytical integration

In this article an approximate analytical integration is performed of the Hamiltonian corresponding to the rotational motion of an Earth whose elastic mantle is deformed by rotation and lunisolar attraction, using Deprit's perturbation method for the first order. Besides the usual terms, this Hamiltonian includes the perturbation of the kinetic energy and the elastic energy produced with the deformation, as well as their causes, the tidal and the centrifugal potential; these new terms have already been studied for the tidal deformation in previous articles (Getino and Ferrándiz, 1990a, 1990b). The effects of the deformation due to the centrifugal potential are studied in this article, following the same method as that used for the tidal deformation. Numerical tables of the periodic perturbations corresponding to the nutation in obliquity and longitude are obtained. As for the secular effects, a theoretical value of 457 days is obtained for Chandler's period.     

关于RSC(IERS)88 C01实现中的几个问题的讨论

本文讨论了IERS的RSC 88C01实现过程中在有关基本源选取和参考架定向确定等方面存在的问题,分析了它们对IERS天球参考架的影响,提出了建立和维持高精度天球参考架的几点建议。