Harmonic analysis of the ionospheric electron densities retrieved from FORMOSAT-3/COSMIC radio occultation measurements

In order to investigate the regular variations of the ionosphere, the least-squares harmonic estimation is applied to the time series of ionospheric electron densities derived from about five years of Global Positioning System radio occultation observations by FORMOSAT-3/COSMIC satellites. The analysis is done for different latitudes and altitudes in the region of Iran. The least-squares harmonic estimation is found to be a powerful tool for the frequency analysis of the completely unevenly spaced time series of radio occultation measurements. Although the obtained results are slightly different from the exact expected cycles (i.e. annual and diurnal components with their Fourier decompositions, and the 27-day period) due to the low horizontal resolution of radio occultation measurements, high vertical resolution of the observations enables us to detect not only the total electron content variations but also periodic patterns of electron densities at different altitudes of the ionosphere. The dominant diurnal and annual signals together with their Fourier series decompositions are obtained, which are consistent with the previous analyses on the total electron content. In the equatorial anomaly band, the annual component is weaker than its Fourier decomposition periods. In particular, the semiannual period dominates the annual component, indicating the relationship between the semiannual variation of the electron densities and the ionospheric equatorial anomaly. From detection of the phases of the components, it is revealed that the annual signal generally has its maximum value in summers at high altitudes, and in the winters at low altitudes. This is probably due to the higher [O/N₂] ratios in winter than in the summer in the lower ionosphere. Furthermore, the semiannual component mostly peaks around solstices or about a month before/after them.     

基于不同天文标准计算地球引力对卫星轨道的影响

建立了一个考虑地球非球形引力摄动的卫星轨道摄动模型并用C＋＋语言编写程序进行轨道计算,其积分器采用Runge-Kutta-Fehlberg7（8）。建模过程中,对岁差、章动等量的计算分别考虑了美国海军天文台1981年和2005年的两种天文标准。通过与STK对比,发现无论采用哪一种标准,轨道传播两天后的位置误差都不会超过分米量级,速度误差不会超过毫米每秒的量级。说明了建模和编程计算的合理性,以及采用新旧两种天文标准计算卫星轨道的差别甚小。     

VLBI for better gravimetry in SELENE

The Japanese lunar explorer SELENE (SElenological and Engineering Explorer), to be launched in 2007, will for the first time utilize VLBI observations in lunar gravimetry investigations. This will particularly improve the accuracy to which the low degree gravitational harmonics and the gravity field near the limb can be measured, and when combined with Doppler measurements will enable three-dimensional information to be extracted. Differential VLBI Radio sources called VRAD experiment involves two on-board sub-satellites, Rstar and Vstar. These will be observed using differential VLBI to measure the trajectories of the satellites with the Japanese network named VERA (VLBI Exploration of Radio Astrometry) and an international VLBI network.     

Simulation study for the determination of the lunar gravity field from PRARE-L tracking onboard the German LEO mission

A simulation study has been performed at GFZ Potsdam, which shows the anticipated improvement of the lunar gravity field model with respect to current (LP150Q model) or near-future (SELENE) knowledge in the framework of the planned German Lunar Explorations Orbiter (LEO) mission, based on PRARE-L (Precise Range And Range-rate Equipment – Lunar version) Satellite-to-Satellite (SST) and Satellite-Earth-Satellite (SEST) tracking observations. It is shown that the global mean error of the lunar gravity field can be reduced to less than 0.1 mGal at a spatial resolution of 50 km. In the spectral domain, this means a factor of 10 (long wavelengths) and some 100 (mid to short wavelengths) improvement as compared to predictions for SELENE or a factor of 1000 with respect to LP150Q. Furthermore, a higher spatial resolution of up to 28 km seems feasible and would correspond to a factor of 2–3 improvement of SELENE results. Moreover, PRARE-L is expected to derive the low-degree coefficients of the lunar gravity field with unprecedented accuracy. Considering long mission duration (at least 1 year is planned) this would allow for the first time a precise direct determination of the low-degree tidal Love numbers of the Moon and, in combination with high precision SEST, would provide an experimental basis to study relativistic effects such as the periselenium advance in the Earth–Moon system.     

“火卫1”轨道预报与动力学分析

使用数值积分的方法对"火卫1"(Phobos)轨道进行预报计算,介绍了计算"火卫1"轨道时考虑的动力模型与数值积分方法,分析了不同动力模型及参量对轨道计算结果的影响,并确定最终选择的模型。最后将计算结果与不同机构发布的最新历表进行比较,计算位置和轨道根数的差异。得到的结果与各星历之间差异的量级相同,进而验证了相关模型和算法的精度和可靠性。该轨道预报算法模型与分析结果将对"火卫1"探测任务提供参考。     

一种空间目标RCS序列的姿态异常检测方法

准确、快速判断空间目标姿态运动模式异常,对于空间目标监测具有重要意义。针对空间目标雷达散射截面 (Radar Cross Section,RCS)序列,提出一种基于小波包分解（Wavelet Packet Decomposition,WPD）能量谱特征的无监督机器学习异常检测方法,并采用单类支持向量机(One Class Support Vector Machine,OCSVM)验证异常检测效果。设置了几种典型异常场景进行仿真分析,试验结果表明,该方法能有效检测出三轴稳定类空间目标发生失稳旋转的姿态异常。相比于传统统计参数特征、小波变换统计参数特征及能量特征的姿态判别方法,具有检测概率高、鲁棒性好的特点。     

动载作用下管内汽水两相流传热特性

设计搭建了旋转平台,采用电加热的方法对过载状态下径向水平管内汽水两相流流动及传热特性进行了实验研究.实验测试段为内径20 mm、长400 mm的有机玻璃圆管.通过调节转台转速来实现不同过载量,进而研究过载状态下汽水两相流传热机理.结果表明:随着过载的增大,管内两相流压力增大,流阻急剧增大,流量减小,空隙率减小,换热量减小,管外传热系数增大.      

新疆于田7.2 级地震前的电离层电磁扰动

利用法国DEMETER 卫星观测的电磁、等离子体等参量的观测数据分析了2008 年3 月21日新疆于田7.2 级地震前离子温度Ti、离子密度Ni、甚低频(VLF) 磁场等出现的异常变化. Ti分析显示震中区上空出现的突变信号在更大空间多次重复出现, 作为地震异常的信度不高. 通过多次重访轨道对比发现H+ 离子密度整体变化形态一致, 但夏秋季节的峰值变化幅度是冬春季节的1 倍. 通过把2008 年的数据与2007 年和2006 年同期的重访轨道数据进行对比, 发现2008 年在于田地震前1 个月左右Ni 在峰值区呈现了显著增强的变化趋势, 并一直持续到地震前, 2008 年2 月以后的峰值幅度是其他两年同期的~1 倍多. 对于磁场频谱数据, 抽取了震中区域上空2000 km 范围内多个轨道的单频(400 Hz)信息, 并在研究区域统一纬度空间(0~60°N) 进行了时间序列对比, 结果发现地震前2 日内震中区域上空400 Hz 谱密度增加了近一个数量级, 反映空间存在明显的VLF 电磁辐射现象. 对于3 月20 日极低频段(ELF) 三分量电磁场的矢量分析显示观测时段内存在一些左旋极化电磁信号, 反映了空间电离层离子密度变化伴随ELF 电磁辐射信息. 综合分析认为, 这次地震前各参量的变化过程与目前地震孕育及电离层耦合机理相吻合.      

Variability study of the crest-to-trough TEC ratio of the equatorial ionization anomaly around 120°E longitude

In this paper, latitudinal profiles of the vertical total electron content (TEC) deduced from the dual-frequency GPS measurements obtained at ground stations around 120°E longitude were used to study the variability of the equatorial ionization anomaly (EIA). The present study mainly focuses on the analysis of the crest-to-trough TEC ratio (TEC-CTR) which is an important parameter representing the strength of EIA. Data used for the present study covered the time period from 01 January, 1998 to 31 December, 2004. An empirical orthogonal function analysis method is used to obtain the main features of the TEC-CTR’s diurnal and seasonal variations as well as its solar activity level dependency. Our results showed that: (1) The diurnal variation pattern of the TEC-CTR at 120°E longitude is characterized by two remarkable peaks, one occurring in the post-noon hours around 13–14 LT, and the other occurring in the post-sunset hours around 20–21 LT, and the post-sunset peak has a much higher value than the post-noon one. (2) Both for the north and south crests, the TEC-CTR at 120°E longitude showed a semi-annual variation with maximum peak values occurring in the equinoctial months. (3) TEC-CTR for the north crest has lower values in summer than in winter, whereas TEC-CTR for the south crest does not show this ‘winter anomaly’ effect. In other words, TEC-CTR for both the north and south crests has higher values in the northern hemispheric winter than in the northern hemispheric summer. (4) TEC-CTR in the daytime post-noon hours (12–14 LT) does not vary much with the solar activity, however, TEC-CTR in the post-sunset hours (19–21 LT) shows a clear dependence on the solar activity, its values increasing with solar activity.     

Second post-Newtonian approximation of scalar-tensor theory of gravity

Deep space laser ranging missions like ASTROD I (Single-Spacecraft Astrodynamical Space Test of Relativity using Optical Devices) and ASTROD, together with astrometry missions like GAIA and LATOR will be able to test relativistic gravity to an unprecedented level of accuracy. More precisely, these missions will enable us to test relativistic gravity to 10^-7–10^-9${10}^{-7}''–''{10}^{-9}$ of the size of relativistic (post-Newtonian) effects, and will require second post-Newtonian approximation of relevant theories of gravity. The first post-Newtonian approximation is valid to 10^-6${10}^{-6}$ and the second post-Newtonian approximation is valid to 10^-12${10}^{-12}$ in terms of post-Newtonian effects in the solar system. The scalar-tensor theory is widely discussed and used in tests of relativistic gravity, especially after the interests in inflation models and in dark energy models. In the Lagrangian, intermediate-range gravity term has a similar form as cosmological term. Here we present the full second post-Newtonian approximation of the scalar-tensor theory including viable examples of intermediate-range gravity. We use Chandrasekhar’s approach to derive the metric coefficients and the equation of the hydrodynamics governing a perfect fluid in the second post-Newtonian approximation in scalar-tensor theory; all terms inclusive of O(c^-4)$\mathrm{O}(c^{-4})$ are retained consistently in the equations of motion.