基于STK的GPS/LEO无线电掩星技术仿真研究

GPS/LEO无线电掩星技术反演地球大气参数剖面已经具有较高的精度. 国外开展了多个GPS/LEO掩星项目, 但中国还尚未深入进行相关的实验, 这制约了中国掩星技术的发展. 本文提出基于STK进行GPS/LEO掩星技术研究的方法; 根据GPS/LEO掩星的原理, 推导出掩星事件发生的条件和掩星切点的计算公式; 利用STK对掩星过程进行模拟, 得到掩星数据. 在大气球对称假设和大气模型已知的条件下, 反演得到中性大气折射指数. 通过比较模型和反演数据, 表明反演数据精度较高, 验证了利用STK模拟GPS/LEO掩星实验方法的可行性.      

GRO和LRO掩星事件模拟研究

GRO（Global Navigation Satellite System Radio Occultation）和LRO（Low Earth Orbit Radio Occultation）联合组网探测地球大气是无线电掩星探测技术的主要发展方向.本文根据掩星事件的数学判据,仿真分析了LEO卫星主要轨道参数对GRO和LRO掩星事件数量和全球分布情况的影响.研究表明:卫星轨道越低GRO掩星事件越多;轨道倾角在30°和75°之间时,GRO掩星事件较多,全球覆盖率也较大;利用极轨卫星进行LRO掩星探测时,LRO掩星事件较均匀地分布在各纬度带.研究成果对GRO和LRO联合星座设计具有参考价值.      

LEO-LEO掩星事件持续时间及卫星间相对角速度仿真分析

通过对LEO-LEO掩星事件仿真,分析卫星轨道参数对掩星事件持续时间和卫星间相对角速度的影响.结果表明:对掩星事件持续时间影响较大的LEO-LEO卫星轨道参数依次为轨道倾角、升交点赤经、轨道高度,近地点角距对掩星事件持续时间影响不大;LEO-LEO卫星轨道倾角互补时,掩星事件持续时间最短,为98s;两颗LEO卫星轨道倾角之和或者升交点赤经之差在120°-240°范围内,掩星持续时间约100~150s;两颗LEO卫星轨道越高,掩星持续时间越短;掩星发生时,LEO卫星间水平方向相对角速度最大值约为0.14（°）·-1,俯仰方向相对角速度最大值约为0.078（°）·s-1.      

掩星大气探测系统误差源分析

根据掩星大气探测反演原理,研究了影响系统产品精度的误差源,并对多普勒观测误差、天线相位中心变化误差、多路径误差、卫星质心变化误差、导航卫星钟飘、卫星姿态变化误差和精密定轨速度误差进行分析,根据工程实现可行性提出了系统误差分配建议,为掩星大气探测系统设计提供参考.      

ROPP反演软件算法及其精度分析

介绍了ROPP反演软件中无线电掩星反演的算法与精度分析. 采用COSMIC卫星2008年1 月1日全天的附加相位数据, 反演得到折射率、温度、压强与湿度等参数, 并与CDAAC 相应结果进行对比. 实验结果表明, 在30km高度以下, 折射率、压强和湿度的相对 误差在2%以内, 温度误差不超过2K.      

FY-3D卫星的北斗掩星分布特征与误差特性

利用2018年1－3月FY-3D卫星的掩星折射率数据,研究了北斗导航卫星系统的掩星分布特点、数据精度以及误差统计特征。北斗导航卫星系统同步静止轨道掩星沿卫星轨道呈弧状分布在南北两极地区,倾斜轨道掩星在东西半球低纬度地区分别形成一小一大两个空洞,中地球轨道掩星则全球均匀分布。北斗掩星折射率数据精度在探测核心区域,即12~32 km范围内,与ERA5再分析资料计算的折射率相比,平均偏差的标准差约为1.5%,在核心区外,标准差从1.5%逐渐增大到6%。静止轨道掩星的平均偏差在高层略大于倾斜轨道和中地球轨道掩星。下降掩星在20 km以上区域的标准差大于上升掩星,20 km以下区域小于上升掩星。高纬地区北斗掩星标准差最小,低纬地区最大,对流层中下层尤其明显。分析结果表明,北斗掩星的数据精度和误差特征与GPS掩星数据相似。      

地基行星大气掩星观测资料流程和数据整理

在中国火星探测萤火一号(YH-1)计划中, 包括了地基掩星观测反演火星大气的科研任务. 观测资料整理是反演流程的第一步. 本文描述了地基火星大气掩星观测处理软件系统的观测数据流程和观测资料整理模块,并详细介绍了观测资料整理模块的结构和功能, 其中包括时间系统转换、历表插值、坐标系变换、信号时延改正以及掩星平面建立. 利用行星数据系统公布的火星快车无线电科学数据和由SPICE得到的地球、火星历表以及火星快车的轨道数据, 结合本文的算法, 得到了一些实验结果.      

基于X射线掩星探测反演行星大气密度

X射线掩星是一种常见的天文现象，基于X射线掩星探测的大气密度反演是一种涉及学科交叉的新方法，其通过处理高能X射线天体辐射源的掩星观测数据实现大气密度的反演，基本原理为X射线在大气中传播时，X射线光子被大气中的原子（包括分子中的原子）吸收和散射，从而导致X射线强度发生衰减，根据衰减后X射线信号的强度反演对应的密度廓线。本文根据X射线掩星探测的应用需求，论证了基于X射线掩星实现大气密度反演的新方法，重点介绍了光变曲线拟合和能谱拟合两种地球中高层大气密度反演算法，分析了X射线掩星探测反演大气密度的研究进展和研究方法，对基于X射线掩星反演大气密度的优点进行分析和讨论，进而对X射线掩星探测的应用场景进行展望。结果表明，作为一种新型中高层大气密度测量手段，X射线掩星探测可对中高层大气密度实现有效探测，弥补了目前中高层大气密度实测数据的不足。     

利用卫星两行轨道根数反演热层密度

两行轨道根数（TLEs）是基于一般摄动理论产生的用于预报地球轨道飞行器位置和速度的一组轨道参数,通过求解大气阻力微分方程,可反演出热层大气密度. 本文选取近圆轨道CHAMP卫星和椭圆轨道Explorer8卫星,以两行轨道根数数据为基础,计算反弹道系数,并根据不同轨道特征采用两种不同反演方法对热层大气密度进行研究. 结果表明,这两种方法反演得到的大气密度与实测值均符合较好,其中CHAMP卫星的反演结果和经验模式值相对于实测值的误差分别为7.94%和13.94%,Explorer8卫星的误差分别为9.04%和14.32%. 相比模式值,利用两行轨道根数数据反演的热层大气密度更接近于实测值,说明该方法可以作为获取大量可靠大气密度数据的一种有效途径.      

低轨航天器弹道系数估算及热层大气模型误差分析

利用低轨（LEO）航天器在轨期间两行轨道根数（TLEs）数据,结合经验大气密度模型NRLMSISE00,反演计算得到其在轨期间的弹道系数B’,以31年B’的平均值代替弹道系数真值,分别通过标准球形目标卫星对比以及物理参数基本相同的非球形目标卫星对比,对弹道系数真值进行了检验;利用不同外形目标卫星弹道系数在不同太阳活动周内的变化规律,结合太阳和地磁活动变化,估计经验大气密度模型的误差分布. 结果表明,利用反演弹道系数31年的平均值来代替真值,其在理论值的正常误差范围内;大气密度模型误差在210～526km高度范围内存在相同的变化趋势,且模型误差随高度增加而增大;在短周期内B’变化与太阳活动指数F_10.7存在反相关性;密度模型不能有效模拟2008年出现的大气密度异常低. 以上结果表明,经验大气密度模型结果需要修正,尤其是在太阳活动峰年和谷年,此外,磁暴期间模型误差的修正对卫星定轨和轨道预报等也具有重要意义.      

Effects of solar activity variations on the low latitude topside nighttime ionosphere

We have studied the topside nighttime ionosphere of the low latitude region using data obtained from DMSP F15, ROCSAT-1, KOMPSAT-1, and GUVI on the TIMED satellite for the period of 2000–2004, during which solar activity decreased from its maximum. As these satellites operated at different altitudes, we were able to discriminate altitude dependence of several key ionospheric parameters on the level of solar activity. For example, with intensifying solar activity, electron density was seen to increase more rapidly at higher altitudes than at lower altitudes, implying that the corresponding scale height also increased. The density increased without saturation at all observed altitudes when plotted against solar EUV flux instead of F10.7. The results of the present study, as compared with those of previous studies for lower altitudes, indicate that topside vertical scale height increases with altitude and that, when solar activity increases, topside vertical scale height increases more rapidly at higher altitudes than at lower altitudes. Temperature also increased more rapidly at higher altitudes than at lower altitudes as solar activity increased. In addition, the height of the F2 peak was seen to increase with increasing solar activity, along with the oxygen ion fraction measured above the F2 peak. These results confirm that the topside ionosphere rises and expands with increasing solar activity.     

Mapping the earth's magnetic and gravity fields from space: Current status and future prospects

Orbital potential field measurements are sensitive to regional variations in earth density and magnetization that occur over scales of a few hundred kilometers or greater. Global field models currently available are able to distinguish gravity variations of ±5 milligal over distances of ～1,000 km and magnetic variations of ±6 gamma over distances of ～300 km at the earth's surface. Regional variations in field strength have been detected in orbital measurements that are not apparent in higher resolution, low altitude surveys. NASA is presently studying a spacecraft mission known as GRAVSAT/MAGSAT, which would be the first satellite mission to perform a simultaneous survey of the earth's gravity and magnetic fields at low orbital altitudes. GRAVSAT/MAGSAT has been proposed for launch during the latter nineteen-eighties, and it would measure gravity field strength to an accuracy of 1 milligal and magnetic field strength to an accuracy of 2 gamma (scalar)/5 gamma (vector components) over a distance of roughly 100 km. Even greater improvements in the accuracy and spatial resolution of orbital surveys are anticipated during the nineteen-nineties with the development of potential field gradiometers and a tethered satellite system that can be deployed from the Space Shuttle to altitudes of 120 km above the earth's surface.     

A Method of Inversing the Peak Density of Atomic Oxygen Vertical Distribution in the MLT Region From the OI (557.7nm) Night Airglow Intensity

In this paper, using the MSISE-90 model as the reference atmosphere, we discuss the feasibility and method of deducing the peak densities of the undisturbed atomic oxygen profiles in the MLT region (the mesosphere and lower thermosphere region) from OI (557.7 nm) night airglow intersities. The peak densities for different seasons, latitudes and longitudes are deduced from OI (557.7nm) airglow intensities through this expression. We analyze the features of inversion relative errors and discuss the influence of the variations in temperature on inversion errors. The results indicate that all inversion errors are less than 5% except for those at high altitudes in the summer hemisphere. And the impact of the variations in temperature on errors is not significant.      

The ratio of hazardous meteoroids to orbital debris in near-Earth space

Orbital debris is known to pose a substantial threat to Earth-orbiting spacecraft at certain altitudes. For instance, the orbital debris flux near Sun-synchronous altitudes of 600–800 km is particularly high due in part to the 2007 Fengyun-1C anti-satellite test and the 2009 Iridium-Kosmos collision. At other altitudes, however, the orbital debris population is minimal and the primary impactor population is not man-made debris particles but naturally occurring meteoroids. While the spacecraft community has some awareness of the risk posed by debris, there is a common misconception that orbital debris impacts dominate the risk at all locations. In this paper, we present a damage-limited comparison between meteoroids and orbital debris near the Earth for a range of orbital altitude and inclination, using NASA’s latest models for each environment. Overall, orbital debris dominates the impact risk between altitudes of 600 and 1300 km, while meteoroids dominate below 270 km and above 4800 km.     

Enhancing the kinematic precise orbit determination of low earth orbiters using GPS receiver clock modelling

Clock error estimation has been the focus of a great deal of research because of the extensive usage of clocks in GPS positioning applications. The receiver clock error in the spacecraft orbit determination is commonly estimated on an epoch-by-epoch basis, along with the spacecraft’s position. However, due to the high correlation between the spacecraft orbit altitude and the receiver clock parameters, estimates of the radial component are degraded in the kinematic approach. Using clocks with high stability, the predictable behaviour of the receiver oscillator can be exploited to improve the positioning accuracy, especially for the radial component. This paper introduces two GPS receiver clock models to describe the deterministic and stochastic property of the receiver clock, both of which can improve the accuracy of kinematic orbit determination for spacecraft in low earth orbit. In particular, the clock parameters are estimated as time offset and frequency offset in the two-state model. The frequency drift is also estimated as an unknown parameter in the three-state model. Additionally, residual non-deterministic random errors such as frequency white noise, frequency random walk noise and frequency random run noise are modelled. Test results indicate that the positioning accuracy could be improved significantly using one day of GRACE flight data. In particular, the error of the radial component was reduced by over 40.0% in the real-time scenario.     

近圆轨道地球观察卫星轨道调整动态分析

本文讨论了使用小推力在有限弧段上调整地球观察卫星轨道偏差的动力学过程。首先,建立调整模型及与之相应的控制方程;其次,进行了轨道变换过程的分析;最后,讨论了推力偏斜的影响。所得结果可供同类卫星轨道调整技术方案设计时参考。      

由电离层电位所建立的大气电状态

本文在已知电离层电位分布下, 解析地计算了大气电位, 电场和电流强度的全球分布.结果表明, 在大气导电率随高度呈指数增加的情况下, 100km高度上的电离层电位, 几乎无衰减地扫到25km以下.大气电场较强的区域主要在20km以下的低层大气区, 其垂直分量比水平分量大4个数量级.而中高层大气电场较弱, 且两分量量级相当.本文还提出了一种考虑地面形状对大气电场影响的解析方法.      

Summary of current radiation dosimetry results on manned spacecraft.

Measurements of radiation exposures aboard manned space flights of various altitudes, orbital inclinations and durations were performed by means of passive radiation detectors, thermoluminescent detectors (TLD's), and in some cases by active electronic counters. The TLD's and electronic counters covered the lower portion of the LET (linear energy transfer) spectra, while the nuclear track detectors measured high-LET produced by HZE particles. In Spacelab (SL-1), TLD's recorded a range of 102 to 190-millirad, yielding an average low-LET dose rate of 11.2 mrad per day inside the module, about twice the dose rate measured on previous space shuttle flights. Because of a higher inclination of the SL-1 orbit (57 degrees versus 28.5 degrees for previous shuttle flights), substantial fluxes of highly ionizing HZE particles were also observed, yielding an overall average mission dose-equivalent of about 135 millirem, about three times higher than measured an previous shuttle missions. A dose rate more than an order of magnitude higher than for any other space shuttle light was obtained for mission STS-41C, reflecting the highest orbital altitude to date of 519 km.     

不同上边界条件下的极区电离层数值模拟

利用一维自洽的极区电离层模型,研究了沿磁力线方向不同电离层-磁层耦合条件下极区电离层的响应.此模型在110-610km的电离层空间区域内,综合求解描述极区电离层的连续性方程、动量方程和能量方程,以得到电离层数值解.研究发现,上边界条件在200 km以上的高度能显著地影响电离层参量的形态.较高的O+上行速度对应较低的F层峰值和较高的电子温度.不同边界O+上行速度对应的温度高度剖面完全不同.200km以上电子温度高度剖面不但由来自磁层的热流通量所控制,同时还受到场向O+速度的影响.对利用电离层模型研究电离层内部物理过程提出了建议.      

Seasonal variations in vertical distribution of meteor decay time as observed from meteor radars at 8.5°N and 80°N

The decay times of meteor radar echoes have been used for decades to investigate characteristics of the mesosphere and lower thermosphere (MLT) region. As the meteor echo decay time depends on background atmospheric parameters, in the present communication, we examine the seasonal variation of the vertical distributions of underdense meteor echo decay times with respect to echo strength. Observations from two similar radars located at two distinct geographical locations, Thumba (8.5°N, 77°E) and Eureka (80°N, 85.8°W) were used for the present study. Here, the radar received signal power is categorized into strong and weak echoes and vertical profiles of their decay times are constructed. It has been noticed that the monthly mean decay time vertical profile turning altitude (i.e., inflection point) varies in the range of 80–87?km of altitude depending on latitude. The turning altitude is observed at relatively lower heights in the winter than in summer at both the latitudes. The present analysis shows that the meteor decay time below the mean turning altitude follows a decreasing trend with decreasing altitude, which is quite distinct to the behaviour of ambipolar diffusion. It is also observed that there is a difference in mean decay time of strong and weak echoes below 90?km of altitude, which is very prominently seen at lower altitudes. This difference shows a seasonal pattern at high latitude, but does not show any seasonal variation at low latitude. The present results are discussed in light of current understanding of the meteor decay time.