非刚体地球的受迫章动奥波策项和倾斜模

讨论了非刚体地球受迫章动奥波策项与简正模表达式中倾斜模的关系。结果表明天球历书极章动中倾斜振项对应于角动量极的章动,在球历书极章动与角动量极的章动奥波策项之和。同时还给出了岁差速率与自转极的章动奥波策项间的数学关系。     

Forced nutations of a two-layer Earth model

In this paper we present a theory of the Earth rotation for a model composed of an inelastic mantle and a liquid core, including the dissipation in the core–mantle boundary (CMB). The main features of the theory are: (i) to be Hamiltonian, therefore the computation of some complex inner torques can be avoided; (ii) to be self-consistent and non-dependent on a previous rigid Earth theory, so there is no need to use transfer functions; (iii) to be analytical, the solution being derived by perturbation methods. Numerical nutation series deduced from the theory are compared with the IERS 96 empirical series, an accuracy better than 0.8 mas in providing celestial ephemeris pole (CEP) offsets .     

RDAN97: An Analytical Development of Rigid Earth Nutation Series Using the Torque Approach

New series of rigid Earth nutations for the angular momemtum axis, the rotation axis and the figure axis, named RDAN97, are computed using the torque approach. Besides the classical J2 terms coming from the Moon and the Sun, we also consider several additional effects: terms coming from J3 and J4 in the case of the Moon, direct and indirect planetary effects, lunar inequality, J2 tilt, planetary‐tilt, effects of the precession and nutations on the nutations, secular variations of the amplitudes, effects due to the triaxiality of the Earth, new additional out‐of‐phase terms coming from second order effect and relativistic effects. Finally, we obtain rigid Earth nutation series of 1529 terms in longitude and 984 terms in obliquity with a truncation level of 0.1 μ (microarcsecond) and 8 significant digits. The value of the dynamical flattening used in this theory is HD=(C-A)/C=0.0032737674 computed from the initial value pa=50′.2877/yr for the precession rate. These new rigid Earth nutation series are then compared with the most recent models (Hartmann et al., 1998; Souchay and Kinoshita, 1996, 1997; Bretagnon et al., 1997, 1998. We also compute a benchmark series (RDNN97) from the numerical ephemerides DE403/LE403 (Standish et al., 1995) in order to test our model. The comparison between our model (RDAN97) and the benchmark series (RDNN97) shows a maximum difference, in the time domain, of 69 μas in longitude and 29 μas in obliquity. In the frequency domain, the maximum differences are 6 μas in longitude and 4 μ as in obliquity which is below the level of precision of the most recent observations (0.2 mas in time domain (temporal resolution of 1 day) and 0.02 mas in frequency domain). This revised version was published online in July 2006 with corrections to the Cover Date.     

The morphologies of volcanic landforms at Central Elysium Planitia: Evidence for recent and fluid lavas on Mars

This paper focuses on physical parameters (flow rates and rheological properties) of lava flows observed in the Central Elysium Planitia (CEP) region of Mars. The flows are modeled as Newtonian fluids, using the Jeffrey's equation and the concept of Graetz number, or alternatively as Bingham fluids. In addition to these approaches, a theoretical model of the shape of shield volcanoes based on the solution for the porous flow of an unconfined aquifer is applied to 5 shields, providing independent quantifications of rheological variations between the shields. This analysis indicates that of the five volcanoes studied, two are partially buried by lava postdating their formation, a result which has been confirmed independently in one case by high resolution images. Our observations reveal that two types of lava flows may be found in the CEP region. The first group is composed of large lava flows with viscosities around ∼2.5×¹⁰⁵ Pa s or yield strengths ranging from 100 to 500 Pa. The second group includes small lava flows of the shield volcanoes and large leveed lava channels on the plains with viscosities below 10³ Pa s, or yield strengths less than 200 Pa. When compared with other volcanic regions on Mars investigated with similar approaches, these latter values are, at present, the lowest inferred for martian lava flows. Several hypotheses for the formation of these lavas are discussed in the context of CEP given that low viscosity can be the result of (1) high temperature, (2) low crystal content, (3) low Si abundance of the liquid phase, and/or (4) the presence of dissolved volatiles. Two scenarios are considered. In the first one, it is demonstrated that low viscosity lavas (of low silica content) can be produced in the context proposed by Schumacher and Breuer [Schumacher, S., Breuer, D., 2007. Geophys. Res. Lett. 34. L12202] for recent volcanism. However, geochemical maps derived from GRS measurements do not provide support for anomalously low silica concentrations in this region. In the second scenario, a water-rich magma is proposed, although arguments in favor of a water-rich mantle source below the CEP are not available at the present time.     

Long-Term Evolution of Orbits About A Precessing Oblate Planet: 1. The Case of Uniform Precession

It was believed until very recently that a near-equatorial satellite would always keep up with the planet’s equator (with oscillations in inclination, but without a secular drift). As explained in Efroimsky and Goldreich [Astronomy & Astrophysics (2004) Vol. 415, pp. 1187–1199], this misconception originated from a wrong interpretation of a (mathematically correct) result obtained in terms of non-osculating orbital elements. A similar analysis carried out in the language of osculating elements will endow the planetary equations with some extra terms caused by the planet’s obliquity change. Some of these terms will be non-trivial, in that they will not be amendments to the disturbing function. Due to the extra terms, the variations of a planet’s obliquity may cause a secular drift of its satellite orbit inclination. In this article we set out the analytical formalism for our study of this drift. We demonstrate that, in the case of uniform precession, the drift will be extremely slow, because the first-order terms responsible for the drift will be short-period and, thus, will have vanishing orbital averages (as anticipated 40 years ago by Peter Goldreich), while the secular terms will be of the second order only. However, it turns out that variations of the planetary precession make the first-order terms secular. For example, the planetary nutations will resonate with the satellite’s orbital frequency and, thereby, may instigate a secular drift. A detailed study of this process will be offered in a subsequent publication, while here we work out the required mathematical formalism and point out the key aspects of the dynamics. In this article, as well as in (Efroimsky 2004), we use the word ‘‘precession’’ in its most general sense which embraces the entire spectrum of changes of the spin-axis orientation -- from the long-term variations down to the Chandler Wobble down to nutations and to the polar wonder.     

On the precession constant: Values and constraints on the dynamical ellipticity; link with Oppolzer terms and tilt-over-mode

The luni-solar precession, derived by theoretical considerations from the precession of the equator, is one of the most important parameters for computing not only precession but also nutations, due to its relation to the dynamical flattening. In this paper, we review the numerical values of this parameter, from the geodynamical point of view as well as the astronomical point of view, from the observational point of view as well as from the theoretical point of view. In particular, we point out a difference of about 1 percent between the global Earth dynamical flattening derived from the astronomical observations and the values derived from the different geophysical computations. The nutation amplitudes depend on the Earth dynamical flattening and this dependence is amplified by a resonance at an important normal mode, the Tilt-Over-Mode (TOM). Since the astronomical point of view as well as the geophysical one are confronted, we also take the opportunity to make the link between the TOM and the expressions of the nutations of the different axes which, in turn, are related with one another by the Oppolzer terms. Both, the Oppolzer terms and the TOM originate from a reference frame tilt effect. In writing the link between the nutational motions of the different axes, and so, in writing the Oppolzer terms, we also make the link with the precessional motion.     

On the frequency and amplitude variations of the δ Scuti star CD−24 7599 (=XX Pyx)

We present 132 h of new time-series photometric observations of the δ Scuti star CD−24 7599 acquired during 86 nights from 1993 to 1996 to study its frequency and amplitude variations. By using all published observations we demonstrate that the three dominating pulsation modes of the star can change their photometric amplitudes within one month at certain times, while the amplitudes can remain constant within the measurement errors at other times. CD−24 7599 also exhibits frequency variations, which do not show any correspondence between the different modes.   The typical time-scale for the amplitude variations is found to be several hundred days, which is of the same order of magnitude as the inverse linear growth rates of a selected model. We find no evidence for periodic amplitude modulation of two of the investigated modes ( f ₂ and f ₃), but f ₁ may exhibit periodic modulation. The latter result could be spurious and requires confirmation. The observed frequency variations may either be continuous or reflect sudden frequency jumps. No evidence for cyclical period changes is obtained.   We exclude precession of the pulsation axis and oblique pulsation for the amplitude variations. Beating of closely spaced frequencies cannot explain the amplitude modulations of two of the modes, while it is possible for the third. Evolutionary effects, binarity, magnetic field changes or avoided crossings cannot be made responsible for the observed period changes. Only resonance between different modes may be able to explain the observations. However, at this stage a quantitative comparison is not possible. More observations, especially data leading to a definite mode identification and further measurements of the temporal behaviour of the amplitudes and frequencies of CD−24 7599, are required.     

Considerations concerning the non-rigid Earth nutation theory

This paper presents the reflections of the Working Group of which the tasks were to examine the non-rigid Earth nutation theory. To this aim, six different levels have been identified: Level 1 concerns the input model (giving profiles of the Earth's density and theological properties) for the calculation of the Earth's transfer function of Level 2; Level 2 concerns the integration inside the Earth in order to obtain the Earth's transfer function for the nutations at different frequencies; Level 3 concerns the rigid Earth nutations; Level 4 examines the convolution (products in the frequency domain) between the Earth's nutation transfer function obtained in Level 2, and the rigid Earth nutation (obtained in Level 3). This is for an Earth without ocean and atmosphere; Level 5 concerns the effects of the atmosphere and the oceans on the precession, obliquity rate, and nutations; Level 6 concerns the comparison with the VLBI observations, of the theoretical results obtained in Level 4, corrected for the effects obtained in Level 5.Each level is discussed at the state of the art of the developments.     

电离层延迟改正模型综述

电离层延迟改正模型通常可以分为广播星历用的预报模型、广域差分用的实时模型、后处理模型3类,不同应用要求需要选择不同的模型。主要比较分析了几种常用的电离层延迟改正模型: 用于广域差分中生成格网模型的三角级数模型、多项式模型、低阶球谐函数模型等都可以获得很好的改正效果,且这3种模型基本上是等价的;电离层延迟谐函数展开模型可以用来分析电离层长时间系列的变化特征;国际电离层参考模型IRI的改正精度一般可以达到60％的效果;而GPS 星历采用的Klobuchar模型的参数设置存在一些不足,对此提出了一些改进措施。     

The wavelet analysis of the period of solar activity

Using the wavelet technique, we analyzed the time series of solar radio fluxes at 2800 MHz as well as sunspot numbers and areas. The results are as follows: (1) These three data sets demonstrate that the most prominent period is 10.69 years and that all other periods are not obvious. (2) The wavelet power spectrum displays the changes of the power spectrum over the entire time-period range and shows the variations in the local power of a given period in a given time interval. Our analysis shows that periods shorter than one year are distinct only around solar activity maximum. (3) The time curves of the wavelet power at three periods (10.69 years, 5.11 years and 155.5 days) for the three times series are rather alike, with the same times of peaks and similar undulations. The shorter the period, the more frequent the fluctuations.     

Methods to correct the variable graduation error of the meridian circle

In the automated(computerized) meridian circle, the graduation error can be calibrated in a short time: a complete determination of the graduation error takes a few days, while a coarse measurement of the first ten dominant Fourier components of the graduation error takes only 15 min. Thus, we can monitor the annual and diurnal variations of the graduation. In our regular observations, the annual variation can be thus corrected for. This kind of correction seems to be necessary, judging from the observing accuracy of modern meridian observations. On the other hand, we could not detect a change of the graduation within one clear day as far as the dominant components of the diameter error are concerned. In our case we can therefore assume the graduation error to be constant within one day.     

Effective Sensitivity of A Non-uniform Phased Array of Short Dipoles

Short dipoles are a key element in new low frequency array antennas as proposed for LOFAR and other astronomical applications. Unfortunately standard texts on short dipole antennas are based on the effective area and do not lead to an astronomically useful sensitivity formulation in a straightforward manner. The concept of maximum effective area is applied to arrays of short dipoles and allows expressing the sensitivity as the ratio of this area over the effective sky brightness temperature as long as the output noise power is dominated by the antenna input radiation. For both quantities we only need to know the array directivity pattern that includes the mutual coupling effects when the actual loading conditions of the array elements are taken into account. Short dipole elements have a constant directivity pattern for frequencies below resonance, but they exhibit strong complex impedance variations that provide only narrow band performance when power matching is applied as required in transmit applications. However, in receive applications voltage or current sensing can be realized, for example with an active balun. Assisted by the steep increase of the sky brightness with wavelength for frequencies below 300 MHz, this can provide sky noise dominated performance over at least a three to one frequency range. Still the low frequency limit is determined by the amplifier noise contribution and the losses in the antenna and in the dielectric ground surrounding the elements. We show that for a sparse array with the elements non-uniformly distributed according to an exponential shell model, a constant sensitivity can be obtained over a frequency range of at least two octaves. In addition, such a configuration has a factor of six greater sensitivity than a rectangular array for a large part of the frequency band.     

Obliquity variations of terrestrial planets in habitable zones

We have investigated obliquity variations of possible terrestrial planets in habitable zones (HZs) perturbed by a giant planet(s) in extrasolar planetary systems. All the extrasolar planets so far discovered are inferred to be jovian-type gas giants. However, terrestrial planets could also exist in extrasolar planetary systems. In order for life, in particular for land-based life, to evolve and survive on a possible terrestrial planet in an HZ, small obliquity variations of the planet may be required in addition to its orbital stability, because large obliquity variations would cause significant climate change. It is known that large obliquity variations are caused by spin-orbit resonances where the precession frequency of the planet's spin nearly coincides with one of the precession frequencies of the ascending node of the planet's orbit. Using analytical expressions, we evaluated the obliquity variations of terrestrial planets with prograde spins in HZs. We found that the obliquity of terrestrial planets suffers large variations when the giant planet's orbit is separated by several Hill radii from an edge of the HZ, in which the orbits of the terrestrial planets in the HZ are marginally stable. Applying these results to the known extrasolar planetary systems, we found that about half of these systems can have terrestrial planets with small obliquity variations (smaller than 10°) over their entire HZs. However, the systems with both small obliquity variations and stable orbits in their HZs are only 1/5 of known systems. Most such systems are comprised of short-period giant planets. If additional planets are found in the known planetary systems, they generally tend to enhance the obliquity variations. On the other hand, if a large/close satellite exists, it significantly enhances the precession rate of the spin axis of a terrestrial planet and is likely to reduce the obliquity variations of the planet. Moreover, if a terrestrial planet is in a retrograde spin state, the spin-orbit resonance does not occur. Retrograde spin, or a large/close satellite might be essential for land-based life to survive on a terrestrial planet in an HZ.     

Diurnal and subdiurnal luni-solar nutations: comparisons and effects

We present a comparison of the diurnal and subdiurnal terms of the three last theories of rigid Earth's rotation: SMART97, RDAN97 and REN 2000. For a better interpretation of the observations, we characterize their contribution to the polar motion and we estimate the non-rigid effects, which are at the level of a few microarcseconds.     

对数周期偶极子天线的电特性分析

正在发展中的云南天文台０．５ＧＨｚ—１．５ＧＨｚ快速射电频谱分析仪的接收天线部分是计划在现有的十米口径抛物面天线的焦点上安装两付互相垂直的对数周期偶极子天线作为馈源。本文介绍了设计一个ＦＯＲＴＲＡＮ语言程序对对数周期偶极子天线电特性进行分析的过程，它可以在整个工作频段范围内计算出任意频率点上一副对数周期偶极子天线两个主平面内的方向图、波瓣宽度、增益和输入阻抗，并把结果以数据文件的形式记录下来，结合绘图软件，就可给出直观的结果，从而为对数周期偶极子天线的设计方案提供了可靠的理论判据。本文最后针对为云南天文台抛物面所设计的馈源作具体的分析。     

Observational manifestations and intrinsic properties of the RCR sources in terms of a unified model

We present a summary results of the study of radio sources showing significant variations of integral flux density using the data from the RATAN-600 surveys of 1980–1994 at a frequency of 7.6 cm. The majority of the detected variable sources have flat radio spectra, although there are also all other spectrum types found. Point and compact sources predominate, although all known morphological structures are found in the sample. Variability is detected both in quasars and galaxies. Using the catalog data, we found brightness variations in the optical and/or infrared ranges for a half of host objects of radio sources. We analyzed the properties of nonvariable and variable RCR sources. We compared the ratio of absolute magnitude to radio luminosity for sources with the active nucleus types determined from the optical data. It is found that this parameter is approximately the same for quasars with different radio luminosity. It isminimum for the strongest radio galaxies and grows up to the level characteristic of quasars with the decrease of radio luminosity. Considering that the ratio depends on obscuring properties of a dust torus, such behavior can be explained if we assume that the torus geometry and its optical depth depend on the source long. This parameter is slightly higher among variable sources than among nonvariable ones which counts in favor of the nucleus more open to an observer.     

A new exact method for line radiative transfer

We present a new method, the coupled escape probability (CEP), for exact calculation of line emission from multi-level systems, solving only algebraic equations for the level populations. The CEP formulation of the classical two-level problem is a set of linear equations , and we uncover an exact analytic expression for the emission from two-level optically thick sources that holds as long as they are in the 'effectively thin' regime. In a comparative study of a number of standard problems, the CEP method outperformed the leading line transfer methods by substantial margins.
The algebraic equations employed by our new method are already incorporated in numerous codes based on the escape probability approximation. All that is required for an exact solution with these existing codes is to augment the expression for the escape probability with simple zone-coupling terms. As an application, we find that standard escape probability calculations generally produce the correct cooling emission by the C  ii 158-μm line but not by the ³P lines of O  i .     

Correcting for the solar wind in pulsar timing observations: the role of simultaneous and low-frequency observations

The primary goal of pulsar timing array projects is to detect ultra-low-frequency gravitational waves. Pulsar data sets are affected by numerous noise processes including varying dispersive delays in the interstellar medium and from the solar wind. The solar wind can lead to rapidly changing variations that, with existing telescopes, can be hard to measure and then remove. In this paper we study the possibility of using a low frequency telescope to aid in such correction for the Parkes Pulsar Timing Array(PPTA) and also discuss whether the ultra-wide-bandwidth receiver for the FAST telescope is sufficient to model solar wind variations. Our key result is that a single wide-bandwidth receiver can be used to model and remove the effect of the solar wind. However, for pulsars that pass close to the Sun such as PSR J1022+1022, the solar wind is so variable that observations at two telescopes separated by a day are insufficient to correct the solar wind effect.     

关于天球参考报

章动序列计算和地球定向参数测定需要一个中间的天球参考极作参照，1984年，采用IAU1980章动理论，选取天球历书极作为参考极，利用改善岁差章动模型和由天文测地新技术确定地球定向参数实现的天球历书极，其精度可达0．1mas，随着理论和观测精度的提高，在微角秒量级下，章动和极移模型中周日和半周日成分分应被考虑，地球定向参数的高频成分已被测定，因此天球历书极的原先定义不再适用，需要更改，叙述了不同天球参考极的概念，天球历书极的定义，评述了天球历书极目前实现及其缺陷，介绍了新的天球参考极－天球中间极的定义及其实现。