武汉中层、低热层大气角谱中频雷达观测

武汉中频雷达是利用中层、低热层中电子密度不均匀体的散射来测量大气的水平风场和电子密度剖面。雷达在计算风场的过程中可得到一些该层大气中电子密度不均匀体的寿命和空间尺度等参量值。结合这些参量和大气风场值可计算得到大气的角谱。其计算方法包括全相关分析技术的谱宽法和空间相关法。前者计算的值被认为是大气角谱的上限值。应用这两种方法,利用2001年2月9日武汉中频雷达的观测数据,对武昌上空中层、低热层大气的角谱进行了计算。结果得到大气角谱随高度增加略微增加,如在68km为6°,90km达最大为11°,其平均值为9.2°。利用2001年2月4—10日86km高度上的数据,得到一个7天平均的日变化曲线,发现大气角谱值在本地夜晚最小。     

昆明MF雷达及初步探测结果

MF雷达是中层大气风场及低电离层电子密度观测的重要手段. 简要介绍了新建昆明MF雷达设备的工作原理、结构及工作模式, 并对观测结果进行初步分析. 对2009年1月观测数据的分析表明, 该月纬向风场最大可达80 m/s, 经向风场则较小, 一般不超过40 m/s, 且二者均呈现一定波动性. 相应的LS谱及谐波拟合分析表明,周日潮汐是80~100 km高度大气风场的主要扰动成分, 其振幅随高度改变, 相位向下传播, 且周日潮汐经向分量相位超前于纬向分量相位. 此外, 分析了MF雷达白天的电子密度观测结果, 并与IRI2000进行比较, 发现两者在变化趋势上有非常好的一致性, 但雷达观测结果小于IRI2000给出的参考值.      

中高层大气风场对地磁活动的响应特征

统计研究漠河、北京、武汉流星雷达观测到的2012-2018年80～100 km高度的风场数据，比较在地磁平静期(Kp≤2)和地磁扰动期(Kp≥4)的日平均风场数据，得到在地磁活动期风场的变化特征。研究结果表明，在地磁扰动时风场变化具有季节差异和纬度差异。地磁扰动期间，纬向风在较高纬度地区倾向于中间层西风增强，低热层东风增强，纬度较低地区倾向于东风增强。春季，地磁活动对纬向风的影响没有纬度差异，在夏冬季随着纬度的降低中间层东风增强明显。地磁活动对经向风的影响具有季节差异，对春冬季节的影响强于夏秋季节。研究表明，地磁活动对纬向风的影响可达9 m·s^–1左右，对经向风的影响可达5 m·s^–1左右。地磁活动对中性大气风场的影响可达80 km。     

中层小尺度风切变的观测

十个中层大气风剖面于1988年夏季在北欧的Andφya(69°N)上空用装有铝箔的火箭测得。使用了一种新型铝箔,使其测量的高度范围覆盖103-85km,测量的高度分辨率为25m.观测的风切变剖面显示了一多层结构,一般由3-9个切变层组成,切变层厚度通常小于200m.观测到高达40-90m/s/km的强风切变,且如此大的切变存在于所测区间内不同高度上。在一个连续测量的五个风切变剖面中,有寿命至少2.5小时,相位向下运动速度为0.4m/s,以及最大切变振幅为180m/s/km的强风切变,显示了稳定和持续特性。      

大气重力波引起的突发钠层的模拟研究

应用钠原子和中性大气分子的质量连续性方程来模拟突发纳层（SSL），垂直风场采用接近实际大气重力波的正弦行波模式，结果较好地反映了SSL的形成过程。SSL的形成时刻在5－15min之间并可持续到30min之后，形成高度大约在90－100km之间，峰宽为0.5-2km之间，这些都与实际观测SSL的特点相符，同时还进一步地研究了当重力波参数（主要指垂直波长和周期）、风速以及常态钠层半宽度发生变化时SSL的变化趋势。     

重力波对中间层和低热层大气环流的影响

利用β通道准地转近似大气平均运动方程组,采用重力波线性饱和参量化方法,定性地研究大气重力 波对中间层和低热层大气环流的作用.模拟计算得到,大气重力波对平均东西风速可产生100m·s-1／d左右的作 用力和产生120 m2/s的湍流扩散,这些作用平衡了Coriolis扭力,导致大气的平均东西风速大大偏离辐射平衡 风场,中层顶附近的平均东西风速在冬季(夏季)甚至反转为东风(西风).平均东西风速计算结果与冬季和夏季 中频雷达东西风速观测值和大气模式剖面等大致一致.     

中国廊坊中间层和低热层大气平均风观测模拟

利用中国廊坊站（39.4°N,116.7°E）流星雷达在2012年4月1日至2013年3月31日的水平风场观测数据,分析廊坊上空80~100km的中间层与低热层（Mesosphere and Lower Thermosphere,MLT）大气平均纬向风和经向风的季节变化特征.结果表明平均纬向风和经向风都表现出明显的季节变化特征.平均纬向风在冬季MLT盛行西风,极大值位于中间层顶,随高度增加西风减弱;在夏季中间层为东风,低热层为强西风,风向转换高度约为82km.平均经向风在冬季以南风为主,在夏季盛行北风.纬向风和经向风在春秋两季主要表现为过渡阶段.流星雷达观测结果与WACCM4模式和HWM93模式模拟的气候变化特点基本一致,但WACCM4模式纬向风和经向风风速偏大,而HWM93模式纬向风和经向风风速偏小.      

利用Aura卫星资料计算全球中层大气背景风场

利用2005年第二代地球观测卫星系统(EOS)的Aura卫星上MLS观测的压强、温度、密度等数据,推算出全球中层大气的平均背景风场,分析了中层大气风场随时间和高度变化的特点.与武汉流星雷达及澳大利亚Adelaide台站观测的比较结果显示,用Aura数据推算出来的风场与实际观测比较符合.与HWM-93模式的比较显示,Aura风场随时间和空间变化的总体趋势与HWM-93基本吻合.特别是在80km以下的高度范围内Aura数据与HWM-93数据符合得比较好;在80 km以上的高度,Aura所算得的风值与HWM-93风值的差别逐渐增大,Aura风值普遍比HWM-93的要大.      

武汉上空中层和低热层大气潮汐的流星雷达观测

武汉流星雷达是2002年元月建成的我国第一部全天空流星雷达，本文对2002年2月19日到7月31日流星雷达观测的潮汐的讨论表明，武汉中层顶以周日潮汐为潮汐运动的主要分量，它的强度远大于半日潮汐，周日潮汐和半日潮汐的波源都在80km以下．周日潮汐分量在3、4月份最强，并且经向分量略强于纬向分量．两个分量的峰值在约95km处出现，分别达到44m／s和60m／s．半日潮的最大值24m／s出现在4月初约93km处．周日潮汐和半日潮汐的振幅和相位随时间呈现出拟周期变化的特征，这可能是潮汐与行星波非线形相互作用的结果．观测结果与GSWM模型的比较表明，GSWM模型在相位随高度变化趋势上与观测结果一致，但模型的周日潮相位比观测约超前1—2h，半日潮相位约滞后1—4h．在周日潮汐较强的月份，模型与观测有较大的差异，观测的幅度通常在95km附近有极大值，而模型并没有极大值．GSWM模型对半日潮的幅度的估计通常过小，观测的半日潮汐幅度有时甚至超过模型值的一倍以上．     

北纬30°N中间层和低热层大气平均风中频雷达观测

利用武汉(30.5°N,114.4°W)中频雷达在2001年1月1日至3月18日、2002年2月5日至3月18日、2002年10月18日至12月31日期间,日本Yamagawa(31.2°N,130.6°W)中频雷达在1997年1月1日至10月17日期间的数据,分析北纬30°N地区上空60-98km高度的中间层、低热层大气平均风的变化规律,结果表明平均纬向风和经向风都具有明显的季节变化.平均纬向风在冬季基本为西风,随高度增加,西风减弱,甚至在上部会出现微弱的东风;夏季中间层表现为强烈的东风,低热层则为西风,风向转换高度在80km附近;春季和秋季为转换季节,在春季出现舌状东风结构.80km附近的平均经向风场在冬季以南风为主,在夏季则以北风为主.不同年份的平均风场存在年际变化性,但其气候变化特点非常相似.中频雷达观测结果与HWM93模式结果的气候变化特点符合很好,与其他纬度的平均纬向风气候变化特点基本类似.     

Mesosphere/lower thermosphere wind regime parameters using a newly installed SKiYMET meteor radar at Kazan (56°N, 49°E)

New meteor radar (MR) horizontal wind data obtained during 2015–2018 at Kazan (56°N, 49°E) are presented. The measurements were carried out with a state-of-the-art SKiYMET meteor radar. Monthly mean vertical profiles of zonal and meridional components of the prevailing wind speeds, also amplitudes and phases of the components of diurnal (DT) and semidiurnal tide (SDT) winds are displayed as contour plots for a mean calendar year over the four recent years and compared with distributions of these parameters provided by the previous multiyear (1986–2002) meteor radar (MR) measurements at Kazan and by the recent HWM07 empirical model. The analysis shows that the SKiYMET zonal and meridional prevailing wind speeds are generally in good agreement, sharing the same seasonal features, with the earlier MR seasonal winds. Comparisons with the HWM07 model are not favourable: eastward solstitial cells as modelled are significantly larger, >30?m/s compared to 15–20?m/s. Also, reversal lines are too variable with height, and the positions of modelled cells (positive and negative) are unlike those of either MRs at Kazan or other MLT radars. Both MR systems provide the large SDT amplitudes, approximately 30?m/s and vertical wavelengths, approximately 55?km, for both components at middle latitudes in winter. They also show the well known strong SDT September feature (heights 85–100?km, the vertical wavelength ～55–60?km), and the weak summer SDT for 80–91?km. HWM07 shows unrealistic amplitudes and phases above 90?km by height and month: minimal amplitudes in equinoxes and no September feature.The weak DT of middle to high latitudes provide similar amplitude and phase structures from both MRs, 1986–2002 and 2015–2017: largest amplitudes (10–12 or 8–10?m/s) for the evanescent meridional tide in summer, peaking in late July; weakest (0–2, 2–4?m/s) at 80 to 92–96?km, when the tide is vertically propagating (January, February, November, December) with a vertical wavelength near 40?km. Again, HWM07 differs in amplitude and phase structures: showing peak amplitudes in equinoxes: April, 15?m/s at 88?km; October, 21?m/s at 89?km.Coupling of the MR wind parameters with the ERA5 wind parameters is studied for a case in 2016. It is shown that the prevailing winds and DT amplitudes and phases of both datasets can be simply linked together, but that the ERA5 SDT amplitudes are significantly underestimated at the top model levels of the ERA5 reanalysis project.     

武汉突发E层的季节变化及其与中层风的关系

通过对2018年武汉站（114.61°E,30.53°N）数字测高仪记录的数据统计分析,研究武汉地区突发E层的特征。研究发现:2018年武汉地区电离层突发E层的临界频率在夏季最高,在冬季有一个次增强现象,在春秋季较低;在正午前最大,在日落后出现小幅度提升,在日出前最低。利用指定动态全球大气气候扩展模型SD-WACCM-X模拟出2018年武汉上空90～140 km高度的平均风场,探讨突发E层与背景风之间的关系,揭示突发E层的形成机制。结果表明,半日潮汐分量可能诱导突发E层临界频率的半日变化,周日潮汐分量诱导突发E层临界频率的周日变化;突发E层的强度可能与纬向风场120 km高度处的风剪切有关。      

2009年7月22日日全食电离层效应观测

首次利用廊坊中频雷达和武汉、嘉兴、廊坊等三站GPS对2009年7月22日日全食电离层效应进行了观测．日食期间,中频雷达D层78km高度上电子密度减小了约67％,电子密度为200cm^-3的高度上升了近10km,GPS／TEC减小了1TECU左右,其变化的最大相位与日食最大相位几乎同步;日食后,观测到周期为2个多小时的电离层扰动现象．     

北纬30°中层顶区域钠与铁原子层的结构和年际变化

基于武汉大学Na和Fe激光雷达在2004年1月至2011年12月期间的观测数据，得到武汉上空中层顶区域Na和Fe原子层的平均特性、夜间变化和季节变化特征.Na层平均质心高度为91.36km，平均RMS（均方根）宽度为4.64km.Fe层平均质心高度为88.99km，平均RMS宽度为4.57km.在充分考虑金属层夜间变化和季节变化对数据样本影响的基础上，获取了Na层和Fe层结构在此期间的年际变化特征.对Na层和Fe层质心高度及RMS宽度的年际变化进行线性拟合，发现Na层和Fe层在此期间均相对稳定，Na层质心高度在近8年间仅有约58m的下降，变化率为-7.91m·a-1，Na层RMS宽度减小约151m，变化率为-20.60m·a-1.同期，Fe层的质心高度下降了约230m，变化率为-31.36m·a-1，Fe层RMS宽度则有所增大，变化率为21.01m·a-1.      

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

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

Study of ionospheric D region changes during solar flares using MF radar measurements

During solar flares, the X-ray radiation suddenly increases, resulting in an increase in the electron density of the atmospheric D region and a strong absorption of short-wave radio waves. Based on Langfang medium frequency (MF) radar, this paper analyzed the variation characteristics of D region in the lower ionosphere from 62 km to 82 km. The analysis focused on multiple C-level and M-level solar flare events before and after the large-scale flare event at 11:53 (UT) on September 6, 2017. The results show that it is difficult to detect the electron density over 70 km in Langfang during solar flares, but the electron density value can be obtained as low as 62 km, and the stronger the flare intensity, the lower the detectable electron density height. Besides, the equal electron density height, the received power of X and O waves will also be significantly reduced during the flares, and the reduction of equal electron density height has a weak linear relationship with flare intensity.