Aura/MLS与TIMED/SABER观测全球重力波特性

大气重力波是临近空间环境主要大气波动之一,对全球环流具有重要影响。卫星上搭载的临边探测器能够探测临近空间大气温度,可用于临近空间大气重力波研究。利用2012－2014年Aura的微波临边探测器（MLS）和TIMED的红外临边探测器（SABER）的探测数据,对20~50 km高度的大气重力波扰动分布特征开展了分析研究,两种观测重力波活动基本一致,重力波随季节、纬度及高度的变化显著。冬季半球高纬度重力波扰动较强,赤道和夏季半球近赤道地区上空也存在明显重力波活动区域,夏季半球高纬度重力波扰动最弱。重力波扰动强度随高度增加。TIMED/SABER重力波扰动强度数值比 Aura/MLS略强。      

一种全球临近空间大气密度建模方法及应用

基于TIMED/SABER卫星2002—2018年观测的20~100 km大气密度数据，统计获得多年月平均值和标准偏差的全球网格数据。利用网格数据，分析了大气密度的变化特征。以网格数据为基准，计算了USSA76的相对偏差，分析了USSA76相对偏差的分布特征。以网格数据为驱动，将大气密度表征为平均值与大尺度扰动量和小尺度扰动量的加和，大尺度扰动和小尺度扰动分别采用余弦函数和一阶自回归模型表征，初步建立了全球临近空间大气密度模型。通过对比模型仿真值与激光雷达观测值，表明模型仿真值与观测值具有较好的吻合度，验证了建模方法的可行性。利用蒙特卡罗方法可再现给定轨迹上所有可能的大气状态。      

神舟4号大气成分探测的新结果

神舟4号(SZ-4)大气成分探测器搭载在SZ-4留轨舱上于2002年12月30日发射入轨，在在轨运行的前3个多月中，正值地球南半球处于夏半球季节，并发生了多次中低强度的地磁扰动事件，SZ-4大气成分探测器测得了轨道舱运行高度上(330-362km附近)大气成分的响应变化和异常现象新结果．探测数据表明，中低强度的地磁扰动事件也能引起热层大气中主要成分O和N2的数密度值增高的响应变化．同样在进入地磁扰动峰期后较高纬度处出现了N2的异常增变和O的异常降变以及大气密度异常扰动的现象，但此期间所出现异常现象的地域与SZ-3和SZ-2大气成分探测器探测结果相反，它仅位于南半球较高纬度地区．     

基于智能图搜索的滑翔式高超声速飞行器路径规划方法

研究了高超声速飞行器在线路径规划问题。首先建立滑翔式高超声速再入飞行器动力学模型,在此基础上综合考虑高超声速飞行器在飞行过程中的动压约束、热流约束、过载约束等,对动力学模型进行解析计算,分析再入走廊,并进一步分析形成滑翔式高超声速飞行器的覆盖区域。基于飞行器覆盖区域并结合工程应用需求,提出速度分段的高超声速飞行器可达区域计算方法,确定在不同终端速度约束下的高超声速飞行器可抵达区域。并以此为基础,以各个速度分段的终端速度为节点,提出基于启发式路径搜索算法的分段路径规划方法,通过变步长图搜索算法可在飞行过程中在线确定各个速度分段区间的终点,完成分段路径搜索规划。最后通过数值仿真验证了所提出路径规划方法的有效性,结果表明,所提出的方法计算效率高,能够有效规划高超声速飞行器的机动飞行轨迹,具有一定的工程应用前景。     

高超声速飞行器热流密度/分层温度/碳化层研究

受试验设备能力限制,地面风洞无法完全模拟高超声速飞行器临近空间热环境。文章采用在飞行器表面开孔安装长时耐高温热流传感器直接测量热流密度的方法,国内首次获得Ma12以上高超声速飞行器表面热流密度时变数据和边界层转捩特征。实测热流值与理论预示值规律相同,两者偏差小于20%。针对树脂基材料导热微分方程中虽考虑了热解吸热项,但未考虑导热系数随温度变化情况,采用在树脂基材料导热微分方程中加入物性参数随温度变化项的方法,计算了飞行器热防护结构内部分层温度和碳化层厚度,并与实测结果进行了比较,不考虑树脂热解特性和材料物性参数随温度变化,理论值高于实测值,最大偏差275~320℃;考虑热解特性和物性参数随温度变化情况,计算值与实测值最大偏差小于70℃。     

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

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

考虑射程损失的高超声速飞行器突防弹道分析与设计方法

研究了高超声速飞行器突防机动造成速度损失进而影响射程的问题。高超声速飞行器在面对拦截威胁时需要靠自身机动进行躲避,而机动躲避将会使得飞行状态产生偏离,进而影响高超声速飞行器射程。针对此问题,首先分别建立弹道式和滑翔式高超声速飞行器运动学和动力学模型,然后考虑突防机动过程中造成的速度、弹道倾角等弹道参数变化,分别对弹道式和滑翔式高超声速飞行器的射程与弹道参数关系模型进行构建,得到突防机动-射程变化关系。之后, 通过数值仿真对突防机动-射程变化二者之间的关系进行验证,结果表明弹道式和滑翔式高超声速飞行器的机动均会导致射程变化,且变化规律与理论分析基本一致,验证了所提出突防机动-射程变化模型的正确性和有效性。最后,基于突防机动-射程变化模型,针对两种高超声速飞行器分别给出对射程变化影响较小的突防机动策略,为提升飞行器飞行性能提供理论和方法基础。     

高超声速翼面气动热与静气动弹性综合分析

研究了考虑热效应影响的高超声速飞行器的静气动弹性问题.建立了考虑热效应影响的高超声速飞行器从气动热计算到静气动弹性分析的快速分析方法,该方法采用分层求解思路对热气动弹性问题进行解耦,利用边界层外的无粘数值求解和边界层内的工程算法相结合的方法,计算高超声速飞行器的表面热流,获得模型的温度场,并基于该温度场计算结构热刚度矩阵.以高超声速飞行器小展弦比翼面为研究对象,计算了该模型在3种飞行状态下的温度场、热应力、热刚度和静气动弹性特性,从而验证了热静气动弹性快速分析方法的可行性和适用性.     

TIMED/SABER与AURA/MLS临近空间探测温度数据比较

利用AURA/MLS数据（V4.2）和TIMED/SABER数据（V2.0）对20~92km高度的大气温度进行比较分析,计算AURA/MLS数据与TIMED/SABER数据的温度绝对偏差,并分析平均温度偏差在不同季节中随经度、纬度和高度的变化特征.结果表明:20~80km高度的平均温度偏差在±6K以内,相对偏差在3%以内;80~90km高度平均温度偏差减小至-10K以下,相对偏差在9%以内.中低纬度地区平均温度偏差廓线的变化趋势一致,从20km高度的-3K左右的负偏差逐渐增加,在45~50km高度的平流层顶处有较明显的3K左右的正偏差峰值.平均温度偏差随纬度变化明显,随经度变化很小.研究结果可为卫星数据的应用提供参考依据.      

基于Cluster观测的磁尾场向电流在南北半球投影位置的分布

利用Cluster四颗卫星的磁场探测数据计算磁尾场向电流并投影到极区电离层,研究其投影位置在南北半球的分布规律,统计过程中去除了强磁暴（磁暴主相Dst<–100 nT）期间的场向电流事件。结果显示:磁尾场向电流事件在极区投影位置的纬度分布具有明显的南北半球不对称性,北半球为单峰结构,南半球为双峰结构。在北半球投影到较低纬度（<64°）的场向电流事件数目明显多于南半球,并且所能达到的最低纬度更低;在南半球投影到较高纬度（>74°）的场向电流事件数目明显多于北半球,并且所能达到的最高纬度更高。地磁平静条件下（|AL|<100 nT）,磁尾场向电流密度随磁地方时（MLT）呈递增趋势,这一结果与低高度卫星在极区对I区场向电流的探测结果符合很好。研究结果表明,磁尾场向电流投影位置的纬度分布呈现出明显的南北不对称性,这与南北半球磁尾场向电流的空间分布以及磁层中磁场结构具有密切关系。      

Hemispheric,seasonal and latitudinal dependence of storm-time ionosphere during low solar activity period

Analysis of the seasonal, hemispheric and latitudinal variation of the ionospheric F2 peak during periods of disturbed geomagnetic conditions in 2011, a year of low solar activity, had been studied using hourly data obtained from low- and mid-latitude ionosonde stations. Our results showed an enhancement in F2-layer maximum electron density (NmF2) at daytime over low latitudes. For the mid-latitude stations, NmF2 depletion pre-dominates the daytime and overturned at nighttime. In general, the variation in terms of magnitude is higher in the low-latitude than at mid-latitude. The nighttime decrease in NmF2 is accompanied by a corresponding F2 peak height (hmF2) increase and overturned at daytime. The hmF2 response during the equinoctial months is lower than the solstices. NmF2 shows distinct seasonal, hemispheric and latitudinal dependence in its response. Appearance of a significant ionospheric effect in southern hemisphere is higher than in the northern hemisphere, and is more pronounced in the equinoxes at low latitudes. At mid-latitudes, the ionospheric effect is insignificant at both hemispheres. A negative ionospheric response dominates the whole seasons at the mid-latitude except for March equinox. The reverse is the case for the hmF2 observation. The amplitudes of both the NmF2 and hmF2 increase with increasing latitude and maximize in the southern hemisphere in terms of longitude.     

MSIS-86和IRI-86模式的扩展应用

本文简要地介绍了MSIS-86和IRI-86两个模式,并用此两模式所给的参数推算出热层中性大气风系、碰撞频率和电离层电导率的时空分布。      

Proposal for a reference model of the middle atmosphere of the southern hemisphere

Available rocketsonde information has been used to compile tables of monthly mean temperature, pressure, density and zonal wind for the middle atmosphere of the southern hemisphere with the purpose of revising similar tables presented to COSPAR earlier. The altitude range is 25 to 80 km in steps of 5 km. The latitude range is 0° to 70°S with a 10° step. The compatability of different sets of temperature measurements is discussed. Mean values of temperature, pressure and zonal wind obtained for the southern hemisphere are compared with northern hemisphere model values. Large differences between the hemispheres (up to 20°C in temperature, 20–30% in pressure, 30–50 m/s in wind) imply that reference atmospheres such as CIRA should be complemented by southern hemisphere climatology.     

南大西洋异常区高能质子通量南北向差异

利用2012年8月NOAA N15,N16,N17,N18卫星中能质子和电子探测器（MEPED）>6.9MeV质子能档的数据,研究了800km高度处南大西洋异常区（SAA）质子分布的南北方向差异.数据显示,质子计数率沿经线随纬度的变化近似满足高斯分布,且向南运动的质子多于向北运动的质子.MEPED具有0°和90°两个探测方向,可对质子投掷角分布进行估算,从而消除在探测南向和北向质子时探测器方向角偏差所造成的影响.结果表明,在经度40°W、纬度13°S至23°S位置上,仍然存在南向质子比北向质子多20%~30%的现象.其原因可能与地磁异常使南向质子镜像点更低,更多质子散射进入大气层沉降有关.分析还发现,南大西洋异常区质子通量存在地方时变化.日侧比夜侧强约20%,这可能是由于地磁场日变化引起的.      

顶部电离层总离子密度经度的变化

利用DMSP F13卫星1996-2005年共10年的观测数据,研究地磁中低纬地区黄昏时段(18:00 LT)顶部电离层总离子密度经度变化的季节、地磁纬度和太阳活动变化特征.结果表明总的经度变化在低纬地区与中纬地区具有明显不同特征.不同经度结构的季节变化均以年变化为主,但纬度分布具有明显差异.一波结构主要集中在中纬地区,且南半球明显强于北半球;二波结构南北半球不对称性非常明显;三波结构和四波结构均为低纬地区明显强于中纬地区.通过分析不同波结构对总经度变化的贡献发现,一波结构在南半球中纬地区贡献最大,二波结构在12月前后的15°N附近贡献较大,三波结构和四波结构仅在低纬地区有较强贡献.在不同太阳活动条件下,不同波结构的贡献率有明显变化.      

基于COSMIC掩星精密定轨数据的等离子体层电子含量研究

等离子体层是日地环境重要的组成部分.本文利用COSMIC掩星精密定轨数据经处理后得到的podTec文件获取等离子体层电子含量（PEC）对等离子体层进行研究.将podTec数据进行处理后获得的PEC（pod-PEC）和IRI-Plas经验模型提供的PEC （IRI-PEC）进行对比，发现pod-PEC与IRI-PEC符合得较好.在低（0°—20°）、中（20°—50°）、高（50°—90°）修正地磁纬度带下，分析了COSMIC在太阳活动极大年（2014年）3，6，9和12月的pod-PEC，得到如下结论:PEC随着纬度升高而逐渐减少，且3，9月PEC在中低纬关于磁赤道的南北对称性较好，6月北半球各纬度带的PEC明显高于南半球同一纬度带的值，而12月情况则完全相反，南半球中纬的PEC甚至会等于北半球低纬的PEC值；PEC在白天高而晚上低，高纬地区的PEC昼夜变化不明显；PEC具有明显的季节性.对于北半球，一年中PEC最大值出现在春季，冬秋季次之，夏季最低，具有明显的年度异常现象.      

A statistical study of the mid-latitude NmF2 winter anomaly

Hourly values of the F2-layer peak density, NmF2, measured by 62 ionosonde stations from 1957 to 2008 at middle geomagnetic latitudes of the northern and southern hemispheres are used in a statistical study of the F2-region winter anomaly. We analyze a maximum value, R, of the winter/summer geomagnetically quiet daytime NmF2 ratio over each ionosonde for approximately the same winter and summer solar activity conditions. The conditional probability of an occurrence of R in an interval of R, the most frequent value of R, the mean expected value of R, and the conditional probability to observe the F2-region winter anomaly are calculated and studied for low, moderate, and high solar activity conditions. It is found that northern hemisphere stations show significant cross-correlations of winter anomaly statistics with geomagnetic latitude for moderate and high solar activities.