皖小行星轨道演化的研究现状

近地小行星由于地球的特殊关系为世人所瞩目,对其研究始于５０年代,近年来已成为太阳系动力学中的一个热门课题。全面地介绍了近年来关于近地小行星的研究状况,包括对近地小行星轨道演化的研究的方法和主要结果。     

共振在近地小行星轨道演化中的作用

本文采用了合理的力学模型和可靠的计算方法,研究了近地小行星轨道运动,揭示了共振（轨道共振和长期共振）在其轨道演化中的作用,并初步探讨了近地小行星的起源问题。     

近地小行星轨道演化的研究

主要阐述近年来在近地小行星轨道演化研究工作中所获得的一些基本结果，即合理的力学模型和相应的有效算法，并以实际预报算例（近地小行星与地球的交会状态）与有关权威性的结果作了比较，证实这些研究结果确实是可信的。在给出的力学模型中，考虑了所有可能影响近地小行星运动的力学因素，包括各大天体和较大的主带小行星的引力作用、有关天体的扁率影响以及源于太阳引力的后牛顿效应。而在计算方法中，合理地处理了变步长问题和月球位置量这种相对而言的快变化问题，使得数值求解一个高维方程组时，对各天体而言，可采用同一步长进行 积分，避免了求解过程中的复杂性。     

短期威胁小行星

为更明确近地小行星撞击地球威胁的监测预警需求, 提出了"短期威胁小行星"的概念, 即未来100yr内可能对地球造成撞击威胁且等效直径大于10m的近地小行星. 以目前已发现的756颗短期威胁小行星为基础, 分析短期威胁小行星的轨道分布特点, 研究显示其与一般近地小行星的轨道分布存在差异, 短期威胁小行星的轨道半长轴更集中于1au, 轨道面更集中于黄道面. 基于近地小行星的数量模型, 初步建立了短期威胁小行星的数量估计模型, 并预估了未来100yr内存在撞击可能的短期威胁小行星的总体数量. 短期威胁小行星的特定研究对制定近地小行星搜巡监测策略有重要意义.     

小行星网络数据库发展评述

评述了近年来国际上与小行星相关的网络数据库的发展概况和意义。概要介绍了小行星轨道数据库、小行星测光数据库、小行星红外数据库、近地小行星数据库和小行星综合数据库的开发成果和应用价值，以及基于小行星数据库的历表服务；展望了小行星网络数据库的发展趋势。     

近地小行星(10302) 1989 ML和(4660) Nereus的轨道特征及演化分析

近地小行星(10302) 1989 ML和(4660) Nereus作为下一代深空探测的候选目标一直备受关注. 在考虑太阳系主要天体的动力学背景下, 通过计算最大Lyapunov指数(MLE)及MEGNO (Mean Exponential Growth factor of Nearby Orbits)指数讨论它们的稳定性. 同时, 对每个小行星, 在其观测误差范围内按多元正态分布各选取1000个克隆粒子, 通过统计分析显示这两个小行星在10万年内可能的运动范围, 给出半长径-偏心率空间中的出现次数分布图, 并统计小行星与地球或其他大行星之间的密近交汇及碰撞的概率. 此外还对这两个小行星的标称轨道进行长期共振、Kozai共振及平运动共振的动力学分析. 综上得出结论, 1989 ML处在平运动共振主导的区域, 发生密近交汇的概率较小, 从而其轨道相对较稳定; 而Nereus处在地球的密近交汇区域, 轨道极不稳定.     

近地小行星发现数目与发现场景的统计分析

近地小行星是一类可能对地球安全造成潜在威胁的太阳系小天体, 目前绝大部分的近地小行星是由地基望远镜发现的, 且数目仍在不断增加. 为了对我国未来开展近地小行星发现监测提供参考和借鉴, 利用国际小行星中心公开的数据库对所有近地小行星首次发现时刻的观测资料开展了多维度统计分析. 发现望远镜探测能力的限制会对近地小行星的发现造成选择效应, 导致不同轨道类型近地小行星发现的相对比例逐年变化且与直径有关. 另外, 结合数值模拟获得的轨道数据, 对近地小行星首次发现时的观测场景进行了还原, 获得了发现时刻近地小行星位置在不同天球坐标系的分布, 分析了其分布特征与季节、测站纬度和小行星直径的依赖关系. 最后, 通过分析数据定量考察了太阳、月球和银道面对近地小行星发现的影响, 发现地基望远镜一般难以发现来自太阳方向90$^\circ$范围内直径140m以下的近地小行星, 并且随着小行星直径的减小该限制范围也将变大; 月光污染对近地小行星发现的影响也非常显著, 望月前后几天的观测限制可导致约29%的目标无法被发现, 而且分析表明农历上半月发现的目标一般比下半月发现的更难以被跟踪观测; 银道面特别是银心方向会对近地小行星发现产生影响, 使得黄道面附近存在与季节相关的观测``盲区''.     

“NEAR-苏梅克”造访“爱神”

美国宇航局于1996年2月17日发射了一艘太空船,主要目标是与国际编号第433号小行星“爱神”交会,对其展开为期近一年的近距离考察,故被命名为“近地小行星交会探测器”,英文缩写就是NEAR,后又更名为“NEAR-苏梅克”。近地小行星有可能穿越地球轨道,对地球有一定的威胁,因此,它们是天文学家严密监视的对象。第433号小行     

宇宙信息

WISE首次发现近地小行星 美国宇航局的大视场红外巡天探测器（WISE）新发现了它的首颗近地小行星,而预计在整个工作寿命中它总共会发现大约数百颗。：这颗直径1千米的小行星编号为2010 AB78,发现于2010年1月12日。发现时它距离地球大约1．58亿千米,不会对地球构成威胁,在未来的几个世纪中它也不会接近地球。这颗小行星在一条椭圆轨道上绕太阳运动,而它的轨道面也并不和地球公转轨道面重合。     

天体力学的几个热点问题

阐述当今天体力学前沿课题中的几个热点问题，近地小行星与地球的交会及其动力演化，航天器定轨新手段中的得一星跟踪自主定轨方法以及星际探测中的轨道力学问题。     

On the theory of the oblateness of the Sun

In this paper we first emphasize why it is important to know the successive zonal harmonics of the Sun's figure with high accuracy: mainly fundamental astrometry, helioseismology, planetary motions and relativistic effects. Then we briefly comment why the Sun appears oblate, going back to primitive definitions in order to underline some discrepancies in theories and to emphasize again the relevant hypotheses. We propose a new theoretical approach entirely based on an expansion in terms of Legendre's functions, including the differential rotation of the Sun at the surface. This permits linking the two first spherical harmonic coefficients (J ₂ and J ₄) with the geometric parameters that can be measured on the Sun (equatorial and polar radii). We emphasize the difficulties in inferring gravitational oblateness from visual measurements of the geometric oblateness, and more generally a dynamical flattening. Results are given for different observed rotational laws. It is shown that the surface oblateness is surely upper bounded by 11 milliarcsecond. As a consequence of the observed surface and sub-surface differential rotation laws, we deduce a measure of the two first gravitational harmonics, the quadrupole and the octopole moment of the Sun: J ₂=−(6.13±2.52)×10⁻⁷ if all observed data are taken into account, and respectively, J ₂=−(6.84±3.75)×10⁻⁷ if only sunspot data are considered, and J ₂=−(3.49±1.86)×10⁻⁷ in the case of helioseismic data alone. The value deduced from all available data for the octopole is: J ₄=(2.8±2.1)×10⁻¹². These values are compared to some others found in the literature. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005238718479     

Accounting for the viscosity of the fluid core in the problem of the physical libration of the moon

A two-component theoretical model of the physical libration of the Moon in longitude is constructed with account taken of the viscosity of the core. In the new version, a hydrodynamic problem of motion of a fluid filling a solid rotating shell is solved. It is found that surfaces of equal angular velocity are spherical, and a velocity field of the fluid core of the Moon is described by elementary functions. A distribution of the internal pressure in the core is found. An angular momentum exchange between the fluid core and solid mantle is described by a third-order differential equation with a right-hand side. The roots of a characteristic equation are studied and the stability of rotation is proved. A libration angle as a function of time is found using the derived solution of the differential equation. Limiting cases of infinitely large and infinitely small viscosity are considered and an effect of lag of a libration phase from a phase of action of an external moment of forces is ascertained. This makes it possible to estimate the viscosity and sizes of the lunar fluid core from data of observations.     

Calculation of the global system of the field-aligned currents on the dayside of the magnetosphere

Using the well-known equation for the normal component of the current which exist near the tangential discontinuity in the plasma in the case of the frozen-in magnetic field, and supposing that the current closes in the ionosphere in the auroral oval in the region 1, one calculates and compares with the data of observations the dependence of the density of the field-aligned current at the level of the ionosphere on the local time.     

Rotation of the elastic Earth: the role of the angular-velocity-dependence of the elasticity-caused perturbation

We calculate the so-called convective term, which shows up in the expression for the angular velocity of the elastic Earth, within the Andoyer formalism. The term emerges due to the fact that the elasticity-caused perturbation depends not only on the instantaneous orientation of the Earth but also on its instantaneous angular velocity. We demonstrate that this term makes a considerable contribution into the overall angular velocity. At the same time the convective term turns out to be automatically included into the correction to the nutation series due to the elasticity, if the series is defined by the perturbation of the figure axis (and not of the rotational axis) in accordance with the current IAU resolution. Hence it is not necessary to take the effect of the convective term into consideration in the perturbation of the elastic Earth as far as the nutation is related to the motion of the figure axis.     

Impact of the Quadrupole Moment of the Sun on the Dynamics of the Earth-Moon System

Range of values of the Sun's mass quadrupole moment of coefficient J₂ arising both from experimental and theoretical determinations enlarge across literature on two orders of magnitude, from around 10^-7 until to 10^-5. The accurate knowledge of the Moon's physical librations, for which the Lunar Laser Ranging data reach an outstanding precision level, prove to be appropriate to reduce the interval of J₂ values by giving an upper bound of J₂. A solar quadrupole moment as high as 1.1 10^-5 given either from the upper bounds of the error bars of the observations, or from the Roche's theory, is not compatible with the knowledge of the lunar librations accurately modeled and observed with the LLR experiment. The suitable values of J₂ have to be smaller than 3.0 10^-6. As a consequence, this upper bound of 3.0 10^-6 is accepted to study the impact of the Sun's quadrupole moment of mass on the dynamics of the Earth-Moon system. Such as effect (with J₂ = 5.5±1.3 × 10^-6) has been already tested in 1983 by Campbell & Moffat using analytical approximate equations, and thus for the orbits of Mercury, Venus, the Earth and Icarus. The approximate equations are no longer sufficient compared with present observational data and exact equations are required. As if to compute the effect on the lunar librations, we have used our BJV relativistic model of solar system integration including the spin-orbit coupled motion of the Moon. The model is solved by numerical integration. The BJV model stems from general relativity by using the DSX formalism for purposes of celestial mechanics when it is about to deal with a system of n extended, weakly self-gravitating, rotating and deformable bodies in mutual interactions. The resulting effects on the orbital elements of the Earth have been computed and plotted over 160 and 1600 years. The impact of the quadrupole moment of the Sun on the Earth's orbital motion is mainly characterized by variations of , , and . As a consequence, the Sun's quadrupole moment of mass could play a sensible role over long time periods of integration of solar system models. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the role of the motion of the photospheric footpoints of the magnetic field lines

By means of a simple relation between the velocity v of the fluid particle and the velocity v_f of the photospheric footpoint of the magnetic field line v_z and B_z being respectively the components of v and the magnetic field B normal to the photospheric surface, it is shown formally that through the phtospheric surface the transport of all the quantities attributed to the magnetic field, such as the magnetic flux, the magnetic energy and the helicity, is independent of v_z, and v_f is the only kinematical quantity on which the transport depends. In addition, in the neighborhood of the neutral line the velocity v_l of the moving curve of constant B_z is found to be equal approximately to the component of v or v_f in the direction of v_l. Since v_l can be measured or extimated, so can the components of v and v_f near the neutral line.     

Determination of the temperature of the solar corona from the spectrum of the electron-scattering continuum

When the K-corona is formed by the scattering of photospheric radiation from free electrons, the Fraunhofer lines are greatly broadened by the thermal motions of the hot electrons. This paper discusses the possibility of measuring the coronal electron temperature from the residual depressions in the K-coronal spectrum. If the ratio of the intensities at 4100 Å and 3900 Å can be measured to an accuracy of ±1%, the coronal temperature can be inferred to an accuracy of ±0.2 MK. The temperature of a coronal inhomogeneity may also be measured by this method, provided the position angle is known.Now at Fraunhofer Institute, Freiburg, Germany.     

On the systematic errors of the determinations of the absolute orientation of the meridian marks

The classical method of determination of the absolute azimuth (or Bessel's parameter n) can secure sufficiently precision for RA from observations of stars at high geographical latitudes during polar night only.