电离层掩星反演的正则最大熵法

电离层掩星数据现已成为电离层观测数据的重要来源,对掩星数据的反演研究一直是掩星研究的热点.传统采用的改正TEC(1btal Electron Content)的Abel变换反演法为线性反演法,它会把测量误差带入反演结果中.为改善反演效果受测量误差的影响,引入两种非线性的反演方法一正则化和正则最大熵反演法.随后设计模拟试验,对3种方法进行验证、比较,得到正则最大熵反演法可以很好地减小测量误差的影响.     

非圆轨道GPS/LEO掩星反演地球大气参数的算法及讨论

在非圆轨道GPS和LEO卫星条件下，给出一种较为直接的GPS／LEO掩星反演地球大气参数技术中弯曲角序列的迭代算法，并在理论上对该迭代法的收敛性进行了严格的数学证明．利用GPS掩星反演模拟程序，定量估算了卫星圆轨道假设对GPS／LEO掩星反演地球大气参数的影响，并验证了在非圆轨道条件下各种迭代法的一致性．指出了文献中给出的级数展开迭代算法的不足之处．     

GPS/LEO掩星技术中超折射效应的修正

在GPS掩星探测地球大气技术中，Abel积分变换要求大气折射指数n是折射半径a的单值函数，当大气折射率的垂直梯度达到小于一个极限值dN／dr≈-0．16N-unit m^-1上述单值性不成立，称其为超折射，此时Abel积分变换不再适用,如果还是在形式上应用经典的Abel变换，在反演结果中就会产生负大气折射率偏差．描述了低对流层中超折射现象的物理特性和数学表示；在广义Abel积分变换的基础上，讨论了超折射层内和超折射层下的大气剖面反演算法；选择了一个简单的采样间隔内等大气折射率垂直梯度假设，对英国高分辨率无线电探空观测资料进行模拟计算，验证了负超折射与大气折射率偏差的关系，并提出的广义Abel变换的合理性.     

火星电离层无线电掩星探测仿真研究

对中俄联合火星星-星电离层掩星技术体制进行了分析和介绍,采用三维射线追踪方法对电离层掩星事件的电波观测值进行了模拟计算,并利用模拟的掩星观测数据进行了电子密度廓线反演,结果说明仿真算法可靠.利用仿真的方法,分别对掩星电波相位观测误差和卫星轨道误差等带来的反演误差进行了个例计算和分析,结果得到:5%周的相位测量误差对白天电离层掩星探测结果的影响可以忽略,而夜间电子密度测量的绝对误差小于4×108 m-3;卫星轨道误差对掩星的主要影响是导致电离层高度抬升或下降.结果表明,中俄联合火星电离层掩星探测技术体制先进,可望获得高精度的电子密度廓线;其技术体制也可以用于月球电离层环境的探测.     

卫星圆轨道假设对GPS无线电掩星反演地球大气参数的影响

给出GPS无线电掩星反演地球大气参数过程中计算大气折射角的解析表达式，以圆轨道假设下的大气折射角计算值为先验约束，采用迭代法对不引入圆轨道假设情况的大气折射角进行归算，在此基础，利用反演方法得到了引入和不引入圆轨道假定两种情况下大气参数（气压和温度）的差分序列，结果表明：卫星圆轨道假设对GPS无线电掩星反演大气参数的影响，在气压方面为1mbar左右，而在气温方面为1K左右，这一结果支持了目前无线电掩星定性误差估计研究中通常引入卫星圆轨道假设这个近似处理方法的合理性，同时也表明：若在高精度反演地球大气参数时，摒弃圆轨道是必要的。     

GPS掩星折射率剖面一维变分同化

近年来,GPS／LEO(全球定位系统／低地球轨道)卫星无线电掩星技术给出了地球大气探测的新途径．从LEO卫星观测到的掩星数据可以反演的地球大气的气压、水汽、温度等剖面;它们对气象和大气科学研究,是具有潜在价值的数据资源．掩星数据资料的同化技术可以有效地改进这些气象参数的剖面,从而改进目前的数值天气预报模式．在当前采用的一维变分同化反演技术中,可以用掩星观测资料的大气折射率或弯曲角剖面进行同化,来反演大气水汽和温度剖面以及海平面压强．作为独立自主开发的GPD／LEO掩星技术系统的一部分,以欧洲中尺度天气预报分析(ECMWF)资料为背景场,CHAMP 掩星观测得到的折射率剖面为观测值,采用Levenberg—Marquardt方法实行GPS掩星资料一维变分同化．在讨论中,用掩星观测点附近相应的探空气球资料来检验CHAMP掩星资料变分同化的结果．     

二维射线跟踪技术中的掩星面坐标系

在大气球对称假设下,可以从弯曲角观测资料通过Abel积分变换反演得到折射率廓线。在低层大气中,折射率水平不均匀性可能导致Abel反演在近地面的误差达到10%。避免该误差的一种有效方法是不通过Abel反演而采用非对称模式下二维射线跟踪算子直接将弯曲角同化到数值天气预报模式中。文章介绍了非对称模式下二维射线跟踪技术的背景、原理和数学模型,讨论了掩星面坐标系与地心惯性坐标系之间的转换关系。     

GPS/LEO掩星观测的变分同化技术

在简单介绍GPS/LEO掩星探测大气的发展历史和科学意义之后,详细阐述了反演的基本原理;分析了标准反演中存在的问题,并说明一维变分同化(1DVAR)在反演方法中的重要性;给出了一维变分同化中价值函数的求解,以及各种同化因子;简单介绍了对当前气象学中普遍使用的四维变分同化(4DVAR);重点讨论了各种同化方法,以及使用各种同化因子的优缺点。最后,通过CHAMP卫星的观测实例分析,验证了GPS数据在数值天气预报(NWP)中的作用,以及相对于标准反演法一维变分对气象要素的改进。     

GPS/LEO掩星技术中Abel积分变换的奇点问题

在GPS／LEO无线电掩星反演地球大气技术中，Abel积分是反演地球大气折射指数的最常用的方法，Abel积分存在积分奇点的问题，根据这一问题讨论了解决奇异积分的几种不同的积分方法，并提出了一种直接求Abel积分变换的解析解．通过数值积分模拟计算，对比了各种不同积分方法引入了的计算误差。     

掩星技术一维变分同化算法的改进

从LEO卫星观测到的掩星数据可以反演地球大气的气压、水汽、温度等廓线;它们对气象和大气科学研究是有价值的数据资源.掩星数据资料的同化技术可以有效地改进这些气象参数廓线,从而改进目前的数值天气预报精度.把掩星观测参数廓线用变分同化方法进入气象业务流程的最大困难是它的计算量太大,通过对一维变分同化价值函数进行改进和对迭代流程进行新的设计,避免了反复计算大维数矩阵的缺点,从而提高了变分同化的计算效率.在适用性讨论中,用背景场向量加上1个高斯白噪声作为真实值来检验CHAMP掩星资料变分同化的结果.     

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 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     

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,     

Research in Astronomy and Astrophysics CONTENTS

正~~     

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.     

Study on the Stressed Mirror Polishing with a Continuous Polishing Machine for Large Aperture Off-axis Aspheric Mirrors

A special stressed annular polishing technique is proposed to mill the off-axis aspheric sub-mirrors of a large segmented mirror with an annular polishing machine. Based on the basic principle of stressed annular polishing technique, a set of special stressing mechanisms are designed to convert milling the aspheric surfaces of sub-mirrors with different off-axis distances into milling the spherical surfaces with identical radii of curvature, so that they can be pol- ished simultaneously on a continuous polishing machine. It took about contin- uous 40 hours to polish a scaled-down mirror of the planning Chinese Future Giant Telescope (CFGT) using this technique. This mirror has the 330 mm di- ameter, 3.6 m off-axis distance, and the 21.6 m radius of curvature, and its max- imum asphericity is 16 micron. The experiment shows that this method has a high effciency, suits batch manufacturing, especially the batch manufacturing of aspheric sub-mirrors of the segmented primary mirror of an extremely large aperture telescope.     

Born–Infeld type phantom model in the ω–ω′ plane

In this paper, we investigate the dynamics of Born–Infeld (B–I) phantom model in the ω–ω′ plane, which is defined by the equation of state parameter for the dark energy and its derivative with respect to N (the logarithm of the scale factor a). We find the scalar field equation of motion in ω–ω′ plane, and show mathematically the property of attractor solutions which correspond to ω _φ ∼−1, Ω _φ =1, which avoid the “Big rip” problem and meets the current observations well.      

A Comparison of Several Image Locating Algorithms in Astronomical Instruments

In astronomical observations at optical wavelengths, a fast image tracking system can be adopted to reduce the effects of the atmospheric seeing and telescopic tracking error, and therefore improve the observing efficiency. Aiming at the need of astronomical observations, totally 5 kinds of algorithms in two categories were selected to make a comparative study on their accuracies and stabilities under different noise conditions by both numerical experiment and laboratory test. The results indicate that the normalized cross-correlation method and barycenter method have not only a higher accuracy but also a better reliability against interferences, they will be applied to the high-resolution spectrograph of the Xinglong 2.16 m telescope and the scienti?c instruments of the SONG (Stellar Observations Network Group) project, respectively.