从地面发射月球探测器的窗口选择

典型的月球探测器飞行轨道包括地球停泊轨道段、地月转移轨道段、月球卫星轨道段和着月轨道段。首先介绍了设计从地面发射月球探测器轨道典型的约束条件;然后,借助于二体假设,建立解析表达式,分析各种约束对窗口选择的影响,给出了各轨道段概略的飞行时间和粗窗口;最后,利用精确的探测器轨道动力学模型,计算精窗口,并给出了一则算例,所得结论可为月球探测器轨道发射、轨道设计提供依据。     

经过月球旁近的低能地月转移轨道

基于日、地、月构成的双圆问题(BCP,Bicircular Problem)研究了经过月球旁近的低能地月转移轨道,总结了这些轨道在相空间的分布特点.首先基于BCP模型,利用BCP系统的不变流形,搜索出经过月球旁近的低能地月转移轨道.然后把时间作为非自治系统相空间的增广维度,给出了能够反映出转移轨道在增广相空间分布情况的状态空间图,研究表明转移轨道以族的形式分布于相空间中,并且任意时刻都可以作为此类轨道的出发时刻.最后分析了不同转移轨道族各自速度增量、飞行时间以及系统能量的变化规律,分别得到了速度增量最优轨道族和飞行时间最优轨道族.     

木星探测轨道分析与设计

研究了与木星探测相关的轨道设计问题.重点关注木星探测轨道与火星、金星等类地行星探测轨道的不同及由此带来的轨道设计难点.首先分析了绕木星探测任务轨道的选择.建立近似模型讨论了向木星飞行需要借助多颗行星的多次引力辅助,对地木转移的多种行星引力辅助序列,使用粒子群算法搜索了2020年至2025年之间的燃料最省飞行方案并对比得到了向木星飞行较好的引力辅助方式为金星-地球-地球引力辅助.结合多任务探测,研究了航天器在飞向木星途中穿越主小行星带飞越探测小行星的轨道设计.最后,给出2023年发射完整的结合引力辅助与小行星多次飞越的木星探测轨道设计算例.     

赋形旋转抛物面天线母线自适应分段拟合

赋形旋转抛物面天线母线拟合常用的方法包括多项式整体拟合和等间隔分段拟合,多项式整体拟合结果阶次高、计算量大,且高阶多项式在边缘处拟合的结果容易振荡;而等间隔分段拟合则是对母线数据按等间隔分段,分段方式具有一定的盲目性,容易造成拟合参数多、光滑性差等问题.针对上述问题,提出了一种基于母线拟合残差分布的自适应分段拟合方法,该方法包括初始整体拟合和分段拟合两步,初始整体拟合用以确定各离散点的拟合残差分布,分段拟合先参照残差分布情况对离散点数据分段,然后采用低阶多项式对各段数据进行拟合.经过实例拟合对比,该方法可避免高阶拟合的不稳定性,减少了分段数,更适用于赋形旋转抛物面天线母线的拟合.     

天体测量法探测系外行星

在目前已发现的系外行星中,绝大多数是由视向速度法和凌星法探测得到的,天体测量法仅发现了1颗.gaia卫星数据即将发布,天体测量法将逐步成为系外行星探测的重要方法之一.基于天体测量法给出的恒星位置参数序列,讨论了在求解行星质量和轨道参数时涉及的动力学条件方程计算问题,给出了具体微分改正公式,同时也进行了必要的仿真模拟计算.建立的方法可以较容易地推广到多行星系统.     

基于二叉树的太阳系小天体数值历表建立方法

在小天体深空探测任务中,经常需要调用力模型较完备的精确小天体数值历表,然而由于部分小天体可能有较大的轨道偏心率以及存在近距离飞越大行星的可能,导致传统的等区间分段切比雪夫多项式插值方法不适用于构建这些小天体的数值历表.在传统方法基础上采用结合误差自适应的二分法对原始轨道进行分段,并采用二叉树结构来组织所有的分段区间,能解决历表读取中的系数反查问题.以大偏心率近地小行星4179 Toutatis和3200 Phaethon、大偏心率彗星2P/Encke以及将近距离飞越地球的99942 Apophis小行星为例,通过数值仿真证实了该算法的有效性.由于该方法没有对小天体的轨道性质做任何要求,因此能够作为构建太阳系小天体数值历表的通用方法,可应用于大量目标小天体的精确轨道设计中.     

SGP4/SDP4模型精度分析

本文基于最新发布的SGP4/SDP4(Simplified General Perturbation Version 4/Simplified Deep-space Perturbation Version 4)模型设计了一套定轨方案,从空间目标库中挑选出不同类型和轨道参数的1120个目标进行计算,定量给出了SGP4/SDP4模型处理不同类型空间目标的定轨预报精度.结果表明:近地目标定轨精度为百米量级;半同步和同步轨道定轨精度平均为0.7和1.9km.椭圆轨道目标的定轨精度与偏心率有关,除少数e>0.8的椭圆轨道目标,绝大多数椭圆轨道目标定轨误差均小于10km.用SGP4/SDP4模型对近地目标预报3天,半同步轨道预报30天,同步轨道预报15天,椭圆轨道预报1天,预报误差一般不超过40km.     

基于无量纲的轨道参数开展卫星性质统计研究

目前已发现了285颗围绕太阳系八大行星公转的卫星, 它们的轨道和物理性质呈现了丰富多样性. 目前为止, 几乎所有的卫星研究工作都基于单个卫星系统或者卫星群, 似乎缺少统一的研究. 提出了一个新的与行星性质无关、只与恒星半径有关的轨道参数n, 定义为以太阳半径为单位的轨道半长轴的自然对数. 不同行星的卫星的n值都存在双极分布, 绝大部分卫星在$n\gtrsim2$区间, 其次在$n\lesssim-1$区间, 位于中间区域的行星则很少. 从卫星物理参数和轨道参数与n的关系中发现, 分属六大行星的卫星有明显的共同特征. 首先, 轨道偏心率和轨道倾角偏大的卫星的n值都在3.5左右, 它们都是巨行星的不规则卫星. 其次, n值在-1和1之间的卫星绝大部分体积大、质量大、反照率高、自转速度慢. 从文献中找到11颗系外卫星候选体, 获得了它们轨道n值和卫星质量, 发现后者也是在-1< n< 1区间最大,其他区间偏小.这些统一的 规律暗示,太阳系内不同行星的卫星形成机制以及太阳系外卫星的形成机制可能一样或类似.     

卫星跟踪卫星模式中轨道参数需求分析

首次基于半解析法利用GRACE(Gravity Recovery and Climate Experiment)双星K波段星间速度误差、GPS接收机轨道误差和加速度计非保守力误差影响累计大地水准面精度的联合模型开展了卫星跟踪卫星模式中轨道参数的需求分析.建议我国将来首颗重力卫星的平均轨道高度设计为400 km和平均星间距离设计为220 km较优.此研究不仅为我国将来卫星重力测量计划中轨道参数的优化选取以及全球重力场精度的有效和快速估计提供了理论基础和计算保证,同时对将来国际GRACE Follow-On地球重力测量计划和GRAIL(Gravity Recovery and Interior Laboratory)月球重力探测计划的发展方向具有一定的指导意义.     

月球探测中转移轨道误差分析和中途修正计算

针对返回型月球探测器,基于月球探测中转移轨道的动力学特征和误差传递矩阵的性质,分别对地月转移轨道和月地转移轨道的误差传递特点进行了研究,并根据误差传递矩阵给出估算第1次中途轨道修正速度增量的线性公式.通过具体算例,给出在实际力学模型下月球探测中转移轨道误差传递性质,讨论了目标点和目标轨道两种不同的轨道修正方法的特点和适用情形,并结合再入约束条件对月地转移轨道第2次中途修正进行了分析和计算.     

Instructions for authors

正Research in Astronomy and Astrophysics publishes original research papers and reviews 011 all liranchcs of astronoiriy and astrophysics.Reviews arc by invitation only.Important new results that require rapid publication can be submitted as a Letter(Letters must be restricted in length     

A thermodynamic model of high temperature lava vaporization on Io

We modified the MAGMA chemical equilibrium code developed by Fegley and Cameron (1987, Earth Planet. Sci. Lett. 82, 207-222) and used it to model vaporization of high temperature silicate lavas on Io. The MAGMA code computes chemical equilibria in a melt, between melt and its equilibrium vapor, and in the gas phase. The good agreement of MAGMA code results with experimental data and with other computer codes is demonstrated. The temperature-dependent pressure and composition of vapor in equilibrium with lava is calculated from 1700 to 2400 K for 109 different silicate lavas in the ONaKFeSiMgCaAlTi system. Results for five lavas (tholeiitic basalt, alkali basalt, Barberton komatiite, dunite, and a molten type B1 Ca, Al-rich inclusion) are discussed in detail. The effects of continuous fractional vaporization on chemistry of these lavas and their equilibrium vapor are presented. The predicted abundances (relative to Na) of K, Fe, Si, Al, Ca, and Ti in the vapor equilibrated with lavas at 1900 K are lower than published upper limits for Io's atmosphere (which do not include Mg). We predict evaporative loss of alkalis, Fe, and Si during volcanic eruptions. Sodium is more volatile than K, and the Na/K ratio in the gas is decreased by fractional vaporization. This process can match Io's atmospheric Na/K ratio of 10±3 reported by Brown (2001, Icarus 151, 190-195). Silicon monoxide is an abundant species in the vapor above lavas. Spectroscopic searches are recommended for SiO at IR and mm wavelengths. Reactions of metallic vapors with S- and Cl-bearing volcanic gases may form other unusual gases including MgCl₂, MgS, MgCl, FeCl₂, FeS, FeCl, and SiS.     

Instructions for authors

正Research in Astronomy and Astrophysics publishes original research papers and reviews on all branches of astronomy and astrophysics.Reviews are by invitation only.Important new results that require rapid publication can be submitted as a Letter(Letters must be restricted in length to 6 printed pages).Authors who submit a paper are expected to be able to certify that the paper is original work,     

Alkali and halogen chemistry in volcanic gases on Io

We use chemical equilibrium calculations to model the speciation of alkalis and halogens in volcanic gases emitted on Io. The calculations cover wide temperature (500-2000 K) and pressure (10⁻⁶ to 10⁺¹ bars) ranges, which overlap the nominal conditions at Pele (T=1760 K, P=0.01 bars). About 230 compounds of 11 elements (O, S, Li, Na, K, Rb, Cs, F, Cl, Br, I) are considered. The elemental abundances for O, S, Na, K, and Cl are based upon observations. CI chondritic elemental abundances relative to sulfur are used for the other alkalis and halogens (as yet unobserved on Io). We predict the major alkali species in Pele-like volcanic gases and the percentage distribution of each alkali are LiCl (73%), LiF (27%); NaCl (81%), Na (16%), NaF (3%); KCl (91%), K (5%), KF (4%); RbCl (93%), Rb (4%), RbF (3%); CsCl (92%), CsF (6%), Cs (2%). Likewise the major halogen species and the percentage distribution of each halogen are NaF (88%), KF (10%), LiF (2%); NaCl (89%), KCl (11%); NaBr (89%), KBr (10%), Br (1%); NaI (61%), I (30%), KI (9%). We predict the major halogen condensates and their condensation temperatures at P=0.01 bar are NaF (1115 K), LiF (970 K); NaCl (1050 K), KCl (950 K); KBr (750 K), RbBr (730 K), CsBr (645 K); and solid I₂ (200 K). We also model disequilibrium chemistry of the alkalis and halogens in the volcanic plume. Based on this work and our prior modeling for Na, K, and Cl in a volcanic plume, we predict the major loss processes for the alkali halide gases are photolysis and/or condensation onto grains. Their estimated photochemical lifetimes range from a few minutes for alkali iodides to a few hours for alkali fluorides. Condensation is apparently the only loss process for elemental iodine. On the basis of elemental abundances and photochemical lifetimes, we recommend searching for gaseous KCl, NaF, LiF, LiCl, RbF, RbCl, CsF, and CsCl around volcanic vents during eruptions. Based on abundance considerations and observations of brown dwarfs we also recommend a search of Io's extended atmosphere and the Io plasma torus for neutral and ionized Li, Cs, Rb, and F.     

The climatic and physiographic controls of the eastern Mediterranean over the late Pleistocene climates in the southern Levant and its neighboring deserts

Modern-day synoptic-scale eastern Mediterranean climatology provides a useful context to synthesize the diverse late Pleistocene (60–12 ka) paleohydrologic and paleoenvironmental indicators of past climatic conditions in the Levant and the deserts to its south and east. We first critically evaluate, extract, and summarize paleoenvironmental and paleohydrologic records. Then, we propose a framework of eastern Mediterranean atmospheric circulation features interacting with the morphology and location of the southeast Mediterranean coast. Together they strongly control the spatial distribution of rainfall and wind pattern. This cyclone–physiography interaction enforces the observed rainfall patterns by hampering rainfall generation south and southeast of the latitude of the north Sinai coast, currently at 31°15′.The proposed framework explains the much-increased rains in Lebanon and northern Israel and Jordan as deduced from pollen, rise and maintenance of Lake Lisan, and speleothem formation in areas currently arid and semiarid. The proposed framework also accounts for the southward and eastward transition into semiarid, arid, and hyperarid deserts as expressed in thick loess accumulation at the deserts' margins, dune migration from west to east in the Sinai and the western Negev, and the formation of hyperarid (< 80 mm yr^− 1) gypsic–salic soils in the southern Negev and Sinai. Our climatic synthesis explains the hyperarid condition in the southern Negev, located only 200–250 km south of the much-increased rains in the north, probably reflecting a steeper rainfall gradient than the present-day gradient from the wetter Levant into its bordering southern and eastern deserts.At present, the rainiest winter seasons in Lebanon and northern and central Israel are associated with more frequent (+ 20%), deeper Cyprus Lows traversing the eastern Mediterranean at approximately the latitude of southern Turkey. Even these wettest years in northern Israel do not yield above average annual rainfall amounts in the hyperarid southern Negev. This region is mainly influenced by the Active Red Sea Troughs that produce only localized rains. The eastern Mediterranean Cyprus Lows also produce more dust storms and transport higher amounts of suspended dust to the loess area than any other atmospheric pattern. Concurrent rainfall and dust are essential to the late Pleistocene formation of the elongated thick loess zone along the desert northern margin. Even with existing dust storms, the lack of rain and very sparse vegetation account for the absence of late Pleistocene loess sequences from the southern Negev and the formation of hyperarid soils.When the north Sinai coast shifted 30–70 km northwest due to last glacial global sea level lowering, the newly exposed coastal areas supplied the sand and dust to these active eastern Mediterranean cyclones. This enforced the latitude of the northern boundary of the loess zone to be directly due east of the LGM shoreline. This shift of coast to the northwest inhibited rainfall in the southern Levant deserts and maintained their hyperaridity. Concurrently, frequent deep eastern Mediterranean Cyprus Lows were funneled along the northern Mediterranean increasing (probably doubling) the rains in central and northern Israel, Lebanon, southwestern Syria and northern Jordan. These storms and rains formed lakes, forests, and speleothems only a short distance north of the deserts in the southern Levant.     

Active hot star shells

Interferometry in the visible provides milliarcsecond spatial resolution and thus has been used for studying the circumstellar environment of active hot stars. In this paper I will illustrate how the visibility modulus and phase can be used to better constrain the physics of Be disks through results from the VLA, the MkII and the GI2T interferometers. I will insist on the importance and the potential of coupling high angular resolution with high spectral resolution to the study of Be shells. Finally I will present a possible study of the circumstellar disk of Be stars using the VLTI. This revised version was published online in July 2006 with corrections to the Cover Date.     

Fitting X-ray Afterglow Light curves of Gamma-ray Bursts by Using the Magnetar Energy Injection Model

The central compact object for some gamma-ray bursts (GRBs) may be a strongly magnetized millisecond pulsar. It can inject energy to the outer shock of the GRB by through the magnetic dipole radiation, and therefore causes the shallow decay of the early afterglow. Recently, from a large number of GRB X-ray afterglows observed by Swift/XRT(X-ray telescope), it is revealed that many of them exhibit the shallow decay about 10²∼10⁴ s after the burst prompt emission. We have fitted the X-ray afterglow light curves of 11 GRBs by using the energy injection model of a magnetar with the rotation period in the millisecond order of magnitude. The obtained result shows the validity and universality of the magnetar energy injection model in explaining the shallow decay of afterglows, and simultaneously provides some constraints on the magnetic field strength and rotation period of the central magnetar.     

Natural superrotation of a rarefied planetary atmosphere

We consider the kinetics of a rarefied rotating planetary atmosphere. The spatial distributions of the atmospheric-gas density and mean angular velocity were determined by analyzing the exact solution of the two-dimensional kinetic equation. We show that the angular velocity of the gas at some distance from the planet could be higher than that in the initial layer starting from which the atmosphere is rarefied. Our model calculations elucidate the superrotation mechanism under consideration.