强地震活动事件期间热层大气密度异常变化的探测结果

地震活动一直是人类非常关注的自然灾害事件, 其对热层大气密度的影响还不是非常清楚. 2008年5月12日中国四川汶川发生震级8.0级强震事件, 随后, 在6月14日日本本州东部, 7月5日鄂霍次克海和7月19日日本本州东海岸发生震级7.0~7.6级强震事件, 利用期间中国星载大气密度探测器在630 km高度上就位探测的热层大气密度对探测结果进行综合分析, 结果表明, 强地震震源中心区域上空附近热层大气密度出现异常降变. 在时间上, 强地震发生前1~3天内就已开始出现大气密度降变, 强地震发生日附近降变达到谷值, 降变比达0.40~0.65. 这种降变的纬度区域范围震前位于强震中心所处纬度的±3º~±9º, 强震发生当日扩大到±8º~±20º.     

2007-2009谷年期间高热层大气密度的变化

在第23至第24太阳活动周的峰年之间,太阳活动谷年具有持续时间长,极低F_10.7太阳辐射通量（低至65）和超长期的零太阳黑子数记录等特点,因此是观测和研究在这种特殊背景下热层大气变化的极好机会.尤其是能充分理解和掌握在宁静环境下热层大气密度对弱太阳活动和小地磁扰动的响应特性.本文利用高度650 km以上星载大气密度探测器2007—2009年的连续探测数据进行分析,结果表明,在太阳辐射通量F_10.7极低值期间,较高热层大气密度对F_10.7的起伏具有更显著的响应变化.当F_10.7由70降至65时,日均大气密度会有4～5倍的显著降变,远大于通常大气模式中的降变值.同样在F_10.7极低值期间,较高热层大气密度对小地磁扰动也具有显著的响应增变,当日K_p指数之和由23增至30时,较高热层大气密度则会有80%～160%的强增变.     

基于轨道预报的热层模式评估

在太阳活动高低年的地磁平静/扰动环境下,利用不同热层大气模式J77,DTM78,MSIS00,JB2008和CHAMP加速度计反演密度,分析有无先验信息条件下的轨道预报误差.结果表明无先验信息的精密轨道预报中,热层模式的性能可能被弹道系数等参数偏差干扰,此时预报误差不能作为模式性能的评价标准.先验信息对轨道预报精度提升非常明显,尤其是地磁扰动期先进热层模式性能得以展现,轨道预报误差为无先验信息情况下的10%～25%.目前热层模式的主要缺陷存在于地磁扰动期.各模式之间的差异是：JB2008模式可以通过线性和单一频率周期项补偿,而J77及DTM78等模式还存在更多频率的误差.本文对不同情况下精密轨道预报的研究结果可为空间碎片碰撞预警等工程实践提供参考.     

一种基于温度参数的热层密度修正方法

热层大气的阻力效应是影响低轨航天器大量空间操作的重要因素, 尤其是经验密度模式, 其固有的至少15%的内符合误差已严重制约航天器轨道计算精度的提高. 针对广泛应用的经验密度模式, 选择物理背景简明、关联参数较少的JACCHIA71模式, 以地磁平静条件下的全球散逸层顶温度最小值T_c及125 km高度拐点温度T_x为对象, 建立密度相对于上述温度参数的条件方程, 推导密度相对于温度参数的解析偏导数, 并给出其最小二乘解. 同时, 利用CHAMP卫星数据对模式进行修正, 模式平均误差从40%降低至3%左右. 通过TG01飞行器的轨道预报比较, 修正前后轨道预报位置精度从2 km提升至1 km左右. 经过CHAMP卫星和TG01飞行器的实测数据检验, 验证了修正算法的正确性和有效性.     

基于电子密度同化的热层大气密度预报方法

采用热层电离层耦合模式TIEGCM和集合卡尔曼滤波同化方法,利用同化COSMIC电离层掩星电子密度数据优化热层电离层参量,并将模式预报的大气密度与CHAMP卫星大气密度数据进行对比,分别开展模拟和实测数据的同化预报实验.在模拟数据同化实验中,状态向量包含温度、风场和离子成分的实验结果表明,仅优化温度即可达到最优的热层大气密度预报效果.在实测数据同化实验中,将温度作为状态向量参数,优化结果表明,循环同化过程中模式预报的大气密度相对偏差的均方根误差在48h内从38%减小到27%,同化稳定时间至少需要30h.预报过程中大气密度预报效果的改善持续时间为34h.这表明电子密度同化能够改善热层大气密度的预报精度,设计的实验方案合理可行,可获得较长的预报时效.     

Use of two-line element data for thermosphere neutral density model calibration

Traditional empirical thermospheric density models are widely used in orbit determination and prediction of low-Earth satellites. Unfortunately, these models often exhibit large density errors of up to around 30% RMS. Density errors translate into orbit errors, adversely affecting applications such as re-entry operations, manoeuvre planning, collision avoidance and precise orbit determination for geodetic missions. The extensive database of two-line element (TLE) orbit data contains a wealth of information on satellite drag, at a sufficiently high spatial and temporal resolution to allow a calibration of existing neutral density models with a latency of one to two days. In our calibration software, new TLE data for selected objects is converted to satellite drag data on a daily basis. The resulting drag data is then used in a daily adjustment of density model calibration parameters, which modify the output of an existing empirical density model with the aim of increasing its accuracy. Two different calibration schemes have been tested using TLE data for about 50 objects during the year 2000. The schemes involve either height-dependent scale factors to the density or corrections to CIRA-72 model temperatures, which affect the density output based on a physical model. Both schemes have been applied with different spherical harmonic expansions of the parameters in latitude and local solar time. Five TLE objects, varying in perigee altitude between 280 and 530 km, were deliberately not used during calibration, in order to provide independent validation. Even with a single daily parameter, the RMS density model error along their tracks can already be reduced from the 30% to the 15% level. Adding additional parameters results in RMS errors lower than 12%.     

Thermosphere density and wind measurements in the equatorial region using a constellation of drag balance nanospacecraft

A mission for in situ thermosphere density and winds measurement is described, based on nanospacecraft equipped with a drag balance instrument (DBI) and a GPS receiver. The mission is based on nanosatellite clusters deployed in three orbital planes. In this study, clusters of 10 nanospacecraft are considered, leading to a mission based on a total of 30 nanospacecraft. The geometry analyzed is a symmetrical one, including an equatorial orbit and two orbits with the same inclination and opposing ascending nodes. The main idea is that, by combining the accurate information on the satellite inertial position and velocity provided by the GPS receiver and the drag acceleration intensity provided by the DBI, due to the orbits’ geometrical configuration, both atmospheric drag and wind can be resolved in a region close to the orbit nodes. Exploiting the Earth oblateness effect, a complete scan of the equatorial regions can be accomplished in the short mission lifetime typical of very low Earth orbit satellites, even in high solar activity peaks, when the expected nanospacecraft lifetime is about 40 days.     

一种NRLMSISE-00模型太阳辐射校正方法

根据空间天气的状态,调整大气模型的相关输入参数能够减小模型的计算误差.通过对比CHAMP卫星在轨大气密度探测数据与NRLMSISE-00模式的计算结果发现,通过调整F_10.7的输入,使轨道大气密度积分的模式计算结果与探测结果之间的误差达到最小,此时的F_10.7被称为理想F_10.7输入（F^*）.进一步的分析发现,F^*与太阳紫外辐射MgII指数存在很好的相关性,因此可以选择其他的太阳紫外辐射代理参数取代F_10.7,从而减小模型计算误差.本文采用神经网络技术,建立新的太阳紫外辐射代理参量F_euv与MgII,F_10.7等的对应模型,能够根据当日参数值计算F_euv.研究结果表明,新的代理参数能够有效减小NRLMSISE-00的计算误差.     

Limb Viewing Hyper Spectral Imager (LiVHySI) for airglow measurements onboard YOUTHSAT-1

The Limb Viewing Hyper Spectral Imager (LiVHySI) is one of the Indian payloads onboard YOUTHSAT (inclination 98.73°, apogee 817 km) launched in April, 2011. The Hyper-spectral imager has been operated in Earth’s limb viewing mode to measure airglow emissions in the spectral range 550–900 nm, from terrestrial upper atmosphere (i.e. 80 km altitude and above) with a line-of-sight range of about 3200 km. The altitude coverage is about 500 km with command selectable lowest altitude. This imaging spectrometer employs a Linearly Variable Filter (LVF) to generate the spectrum and an Active Pixel Sensor (APS) area array of 256 × 512 pixels, placed in close proximity of the LVF as detector. The spectral sampling is done at 1.06 nm interval. The optics used is an eight element f/2 telecentric lens system with 80 mm effective focal length. The detector is aligned with respect to the LVF such that its 512 pixel dimension covers the spectral range. The radiometric sensitivity of the imager is about 20 Rayleigh at noise floor through the signal integration for 10 s at wavelength 630 nm. The imager is being operated during the eclipsed portion of satellite orbits. The integration in the time/spatial domain could be chosen depending upon the season, solar and geomagnetic activity and/or specific target area. This paper primarily aims at describing LiVHySI, its in-orbit operations, quality, potential of the data and its first observations. The images reveal the thermospheric airglow at 630 nm to be the most prominent. These first LiVHySI observations carried out on the night of 21st April, 2011 are presented here, while the variability exhibited by the thermospheric nightglow at O(¹D) 630 nm has been described in detail.     

太阳极紫外辐射对热层大气密度的影响

对2001－2021年SOHO卫星的极紫外辐射测量数据,以及CHAMP,GRACE-A和SWARM-C卫星资料推导出的高分辨率大气密度数据进行统计分析,发现大气密度与极紫外测量值的相关系数大于密度与F_10.7指数的相关系数,证实极紫外辐射在不同地方时的影响程度存在显著差异,从而驱动大气密度的周日变化。利用三颗卫星的高度差异揭示极紫外辐射对大气密度的加热效应在350～500 km范围随着高度增加而减弱。统计得到极紫外辐射影响在地方时和纬度上的空间差异：对夏季半球的影响大于冬季半球;在白天,对中纬度地区的影响高于赤道和高纬度地区;在夜间,密度对辐射的斜率在夏季半球高纬度地区存在峰值,在冬季半球中纬度存在谷值,模型DTM2000和NRLMSISE00未能准确刻画。为了改进经验模型,提出基于球谐函数的拟合方法,优于主流模型周日效应采用的表达式,对周日效应建模和修正提供有益借鉴。利用昼夜间能量传输和热层大气经向环流机制探讨了统计结果的物理机制。     

EUV-TEC proxy to describe ionospheric variability using satellite-borne solar EUV measurements: First results

Primary photoionisation of major ionospheric constituents is calculated from satellite-borne solar EUV measurements. Number densities of the background atmosphere are taken from the NRLMSISE-00 climatology. From the calculated ionisation rates, a proxy termed EUV-TEC, which is based on the global total ionisation is calculated, and describes the ionospheric response to solar EUV and its variability. The proxy is compared against the global mean ionospheric total electron content (TEC) derived from GPS data. Results show that the EUV-TEC proxy provides a better overall representation of global TEC than conventional solar indices like F10.7 do. The EUV-TEC proxy may be used for scientific research, and to describe the ionospheric effects on radio communication and navigation systems.     

Variations of the O(D) and O(S) peak volume emission rates without direct solar effects

Airglow volume emission rates of the O(¹D) red line at 630.0 nm and the O(¹S) green line at 557.7 nm were measured by the Wind Imaging Interferometer (WINDII) on the Upper Atmospheric Research Satellite (UARS) during 1991–1997. Focus of this study is on the peak volume emission rates of the two airglows after removing the direct solar effect, which are referred to as the ‘dark’ peak emission rates. The main results are as follows. For the red line emission, at low and mid-latitudes the daytime variation does not have a clear pattern except an enhancement at dusk; during nighttime there is an enhancement in the equatorial region at 20–03 h, which has a semiannual variation with maxima at equinoxes; at solstices the daytime O(¹D) dark emission rate is stronger in winter than in summer. For both the green line E-and F-layers the distribution of the dark peak volume emission rate is symmetric about noon in all seasons, symmetric about the equator at equinoxes, and stronger in summer than in winter. The O(¹S) E-layer is profoundly affected by tides. For the first time the diurnal and semidiurnal amplitudes for the emission rates are derived using 24-h zero-sun data. The amplitude of the diurnal tide can be as large as 20% of the mean peak volume emission rate, and has maxima at the equator and about 40°N/S, and minima at about 20°N/S. The daily diurnal maximum is at noon at the equator but at midnight at 40°N/S. There is a clear semiannual variation of the diurnal amplitude in the equatorial region with maxima at equinoxes. The amplitude of the semidiurnal tide is mostly less than 10% of the mean peak volume emission rate with maximum amplitudes at noon and midnight. There is an annual variation of the semidiurnal amplitude at mid-latitudes peaking in summer. Aurorae appear in all three emission layers day and night. The green aurorae are brighter than the red aurorae, and the green E-layer aurorae are 2–3 times stronger than the F-layer aurorae. The green aurora has a clear midday gap in the F-layer and an afternoon gap in the E-layer. The red aurorae are particularly strong in the so-called cusp region at equinoxes.     

基于CHAMP卫星观测数据对热层大气密度的经验正交分析

对2003年(太阳活动较高年)至2007年(太阳活动低年) CHAMP卫星的热层大气密度观测数据进行了经验正交函数(EOF)分析, 得到了400 km高度上白天平均大气密度ρ的太阳活动周变化与年度变化等不同变化分量. 研究结果表明, ρ受太阳活动影响较大, 其太阳周变化分量与F10.7指数变化之间的相关系数可高达94.5 %; ρ的太阳周变化分量随纬度增加而减小, 且在中高纬地区, 南半球的值明显大于北半球的值, 在低纬地区则出现基本对称的双峰分布, 即赤道质量密度异常(EMA)结构. 在ρ的年变化中, 呈现出明显的季节变化, 即夏季低冬季高; 同时ρ的年变化幅度随太阳活动增加而增强, 随纬度增加而增强. 将本文结果与经验模式NRLMSISE00在观测条件下的输出数据进行对比, 发现两者的太阳周变化与年变化分量基本一致, 但本文观测数据的太阳周成分随纬度变化略小, 年变化幅度略大, 且NRLMSISE00模式不能再现EMA结构. 研究结果对揭示热层气候学变化特征具有重要意义.     

利用二维核回归估计的大气密度模式修正

传统经验大气密度模式预测大气密度存在的较大误差会引起低轨卫星轨道预报误差,对卫星的再入轨、控制计划、碰撞规避及精密定轨造成不利影响.利用天宫一号卫星探测数据,针对大气NRLMSISE-00模式计算的误差特点,在地磁相对平静(Ap ≤ 30)的时间段内,对相近地方时和纬度的模式误差分布进行分析发现,相近地方时和纬度的模式误差分布基本相同.利用二维核回归估计方法,对与预测点相近地方时和纬度的样本误差进行加权,估计预测点处的模式误差,进而按距离预测日期天数的长短,采用加权修正法对模式预测结果进行修正,修正后大气模式误差的均方差(RMS)由14.09%降至4.05%.研究结果表明,该修正方法可以显著提高大气密度预报精度.