用武汉流星雷达观测象限仪座流星雨期间的流星速度

利用武汉流星雷达,首次成功地观测了象限仪座流星雨及流星雨期间的流星速度,讨论了利用单站全天空流星雷达观测流星雨的相关问题.从观测结果可以发现此次象限仪座流星雨发生在2004年1月4日的0000-0800LT,其中流星峰值出现在0400LT,而且通过流星雷达观测到的流星雨期间的流星回波平面推测得到的流星雨辐射点也与该流星雨的理论辐射点位置对应非常好.利用流星回波振幅的Fresnel振荡方法计算了此次流星雨期间观测到的流星的速度,分析了该流星速度的分布,这次流星雨期间观测到的流星速度主要集中在10-30 km/s,可以看出这种速度分布是由流星雨进入地球大气的初始速度和流星在大气中的减速过程共同决定的.最后研究了流星速度随高度的变化,并且由此讨论了地球大气对于流星体的减速作用.      

流星雷达系统相位差偏差的估计和校正

介绍了一种新的流星雷达系统的相位偏差估计和校正方法.利用流星回波的观测数据,用回波信号在各个接收通道之间的相位差,结合干涉式接收天线阵的几何关系,建立了各天线相位差测量值与偏差值之间的线性方程组,利用最小二乘法求解方程组,得到了流星雷达系统各个接收通道之间的相位差偏差估计值及校正后的流星回波到达角.与已有的流星雷达相位偏差估计和校正的方法相比,这种方法可以通过流星雷达的观测数据来计破算雷达系统各个接收天线通道之间的相位差偏差量,而不需要增加额外的硬件,实现了对观测数据的事后处理,可以方便地对已有数据进行校正.以2004年4-6月的武汉流星雷达观测数据为例,计算了流星雷达系统的偏差估计量,并用校正后的数据来计算流星回波的空间位置.结果表明,校正后流星回波数在各个方向上随高度的分布比校正前更符合统计分布.      

流星参数的数值分析和流星高度分布模型

流星余迹能够被后向散射雷达观测到, 利用观测结果, 可以分析和研究流星的空间分布和时间变化规律. 同时, 利用流星空间分布还可以进行空间碎片的研究. 基于标准理论, 对影响雷达回波功率的主要因素, 例如如双极扩散、余迹的初始半径、流星的有限速度, 以及雷达的脉冲重复频率在不同频率和速度下进行了数值分析和计算, 得到的流星衰减时间及双极扩散系数的观测结果与理论结果一致. 通过对昆明流星雷 达观测到的571632个流星进行统计分析, 得到了流星高度分布统计模型, 并利用该模型的分析结果与不同月份流星的观测数据进行对比, 结果比较一致.      

北极69°N和78°N空间碎片凝视探测的对比分析

利用北极69°N和78°N两套非相干散射雷达的首次空间碎片联合观测数据进行空间碎片参数（距离、速度、散射截面积、等效直径等）的对比分析,得出以下结论:两部雷达探测的碎片高度均主要分布在500～1100km和1400～1600km区间,但78°N雷达探测的碎片数量较多;空间碎片的径向速度均在-1.5～1.5km…-1区间,其中大部分为负值,说明在此次探测试验中碎片运动方向主要以远离雷达或地球为主;ESR雷达探测的空间碎片射截面积约为10-5～10-2m2,等效直径主要分布在4～10cm,而UHF雷达探测的空间碎片散射截面积约为10-6～10-2m2,等效直径主要分布在2～6cm,说明在同一高度上69°N雷达探测能力更强;经合理设置判据参数后得出重复检测次数,78°N雷达和69°N雷达分别有32次和14次重复检测,两部雷达共有4次重复检测.这些结果为空间碎片检测和建模提供了参考.      

全天空流星雷达相位差监测分析方法研究

介绍了一种直接利用流星观测数据, 根据流星的时空分布特性和天线的空间布阵关系建立数学模型, 对全天空流星雷达各天线通道的相位差进行监测分析和估算的新方法. 通过分析全天空流星雷达的流星观测数据, 获得了流星时空分布特性, 天线阵相位差变化对流星空间分布的影响, 特别是流星高度分布的标准差特性与各天线相位差的关系. 在此基础上, 模拟研究了利用流星高度分布的标准差来估算天线相位差的偏差, 并应用于中国三亚地区全天空流星雷达进行相位差监测分析和校正. 结果表明, 新方法无需任何附加硬件, 通过日常观测数据就能对某一通道的相位差变化或多个通道的相位差变化进行估算和分析,相位差监测精度优于2°. 对这些相位差变化进行校正, 可有效提高全天空流星雷达对流星的定位测量精度.      

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

利用中国廊坊站（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模式纬向风和经向风风速偏小.      

武汉上空偶发E层与流星雨的联系

利用全天空流星雷达和数字测高仪联合观测数据,初步分析了武汉(30°32′N,114°22′E)上空2002年3次大的流星雨(英仙座、狮子座、双子座)及2003年象限座流星雨爆发期间Es的出现率变化趋势及其与流星观测量的关系.结果表明,在流星雨爆发时,雷达观测到的流星数目会显著增加;流星雨峰值过后,Es出现率也会增大.流星数和Es出现率间的平均互相关函数有两个峰值,第一个峰值平均出现在流星雨爆发后的第2天前后,第二个峰值平均出现在流星雨后第6天前后.我们认为,第一个峰值对应于流星雨期间的大量流星雨电离产生的Es,第二个峰值对应于流星雨后沉积的金属离子在风剪切作用下形成的Es.     

长时间流星不均匀体母体特征个例研究

流星体坠入地球大气烧蚀电离产生流星等离子体尾迹,在等离子体不稳定性过程作用下产生流星不均匀体.利用光学视频和无线电雷达在低纬三亚开展流星体烧蚀和流星不均匀体综合探测结果,发展了一种获取流星不均匀体母体（流星体）特征参数的方法,并对2015年12月双子座流星雨期间观测的一次长持续时间流星不均匀体事例进行了分析,得到了其母体速度、质量和轨道等特征,结果显示产生这次流星不均匀体的流星体速度和轨道等具有双子座流星特点.该方法可应用于流星不均匀体及其母体特征研究.      

北纬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模式结果的气候变化特点符合很好,与其他纬度的平均纬向风气候变化特点基本类似.     

微流星消融对中高层区域金属层的影响

采用连续的微流星体消融模型,对小微粒的流星颗粒(半径＜100μm)消融过程进行分析,并估计对中高层区域(80～110 km)金属层的影响．无线电流星雷达与激光雷达的探测结果表明,半径在100μm以上的流星颗粒的消融物并不是金属层主要的源．本文从流星颗粒的动量方程和能量方程出发,计算不同半径下的流星微粒在不同的注入速度下的消融剖面,并分析了流星消融截止高度的变化趋势．结合长期的激光雷达观测,计算表明,10～4μ半径的流星微粒是普通金属层的主要影响因素．     

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.     

Antarctic meteor observations using the Davis MST and meteor radars

This paper presents the meteor observations obtained using two radars installed at Davis (68.6°S, 78.0°E), Antarctica. The Davis MST radar was installed primarily for observation of polar mesosphere summer echoes, with additional transmit and receive antennas installed to allow all-sky interferometric meteor radar observations. The Davis meteor radar performs dedicated all-sky interferometric meteor radar observations. The annual count rate variation for both radars peaks in mid-summer and minimizes in early Spring. The height distribution shows significant annual variation, with minimum (maximum) peak heights and maximum (minimum) height widths in early Spring (mid-summer). Although the meteor radar count rate and height distribution variations are consistent with a similar frequency meteor radar operating at Andenes (69.3°N), the peak heights show a much larger variation than at Andenes, while the count rate maximum-to-minimum ratios show a much smaller variation. Investigation of the effects of the temporal sampling parameters suggests that these differences are consistent with the different temporal sampling strategies used by the Davis and Andenes meteor radars. The new radiant mapping procedure of [Jones, J., Jones, W., Meteor radiant activity mapping using single-station radar observations, Mon. Not. R. Astron. Soc., 367(3), 1050–1056, doi: 10.1111/j.1365-2966.2006.10025.x, 2006] is investigated. The technique is used to detect the Southern delta-Aquarid meteor shower, and a previously unknown weak shower. Meteoroid speeds obtained using the Fresnel transform are presented. The diurnal, annual, and height variation of meteoroid speeds are presented, with the results found to be consistent with those obtained using specular meteor radars. Meteoroid speed estimates for echoes identified as Southern delta-Aquarid and Sextantid meteor candidates show good agreement with the theoretical pre-atmospheric speeds of these showers (41 km s⁻¹ and 32 km s⁻¹, respectively). The meteoroid speeds estimated for these showers show decreasing speed with decreasing height, consistent with the effects of meteoroid deceleration. Finally, we illustrate how the new radiant mapping and meteoroid speed techniques can be combined for unambiguous meteor shower detection, and use these techniques to detect a previously unknown weak shower.     

廊坊中频雷达风场FCA算法的研究及应用

中频雷达是获得临近空间60～100 km范围内大气风场的重要探测设备之一. 全相关分析方法是空分天线模式中频雷达的一种重要风场反演方法. 通过介绍全相关分析方法算法的一般理论, 根据廊坊中频雷达的天线配置, 分别研究了利用三副天线数据与四副天线数据的风场反演算法FCA-3A和FCA-4A, 讨论给出了全相关 分析方法的剔除标准. 用廊坊中频雷达实测的回波相位和振幅, 得到由该算法反演 的大气风场. 为验证算法的可靠性, 将其结果与ATRAD公司的反演结果进行比较. 结果显示大气纬向风和经向风随时间和高度的变化特性基本一致. FCA-3A与FCA-4A的风场反演结果基本一致, FCA-4A略好于FCA-3A.      

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

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

Plasma and electromagnetic wave simulations of meteors

Every day billions of meteoroids impact and disintegrate in the Earth’s atmosphere. Current estimates for this global meteor flux vary from 2000 to 200,000 tons per year, and estimates for the average velocity range between 10 km/s and 70 km/s. The basic properties of this global meteor flux, such as the average mass, velocity, and chemical composition remain poorly constrained. We believe much of the mystery surrounding the basic parameters of the interplanetary meteor flux exists for the following reason, the unknown sampling characteristics of different radar meteor observation techniques, which are used to derive or constrain most models. We believe this arises due to poorly understood radio scattering characteristics of the meteor plasma, especially in light of recent work showing that plasma turbulence and instability greatly influences meteor trail properties at every stage of evolution. We present our results on meteor plasmas simulations of head echoes using particle in cell (PIC) ions, which show that electric fields strongly influence early stage meteor plasma evolution, by accelerating ions away from the meteoroid body. We also present the results of finite difference time domain electromagnetic simulations (FDTD), which can calculate the radar cross section of the simulated meteor plasmas. These simulations have shown that the radar cross section depends in a complex manner on a number of parameters. These include the angle between radar and meteor entry, a large dependence on radar frequency, which shows that for a given meteor plasma size and density, the reflectivity as a function of probing radar frequency varies, but typically peaks below 100 MHz.     

The 16-day waves in the mesosphere and lower thermosphere over Wuhan (30.6°N, 114.5°E) and Adelaide (35°S, 138°E)

Winds from a meteor radar at Wuhan (30.6°N, 114.5°E) and a MF radar at Adelaide (35°S, 138°E) are used to study the 16-day waves in the mesosphere and lower thermosphere (MLT). The height range is 78–98 km at Wuhan and 70–98 km at Adelaide. By comparison, it is found that the zonal components at both sites are generally larger than the meridional ones, and eastward motion of the zonal background winds is favorable for the 16-day waves penetration to the MLT region. The zonal maximum amplitude appears in the autumn (September–October) around 86–98 km at Wuhan and in the winter months and early spring (July–October) around 72–82 km at Adelaide. Differences are found in wave amplitudes and time of appearance between the two years of 2002 and 2003. In 2003, the intensity of the wave amplitudes is relatively smaller than that for 2002 at both sites. The summer 16-day waves are comparatively weaker at Adelaide in both years, but stronger in 2002 at Wuhan near the mesopause and the lower thermosphere (86–98 km). The strong summer waves at Wuhan may come from the winter southern hemisphere.     

Comparison of mesopause region meteor radar winds,medium frequency radar winds and low frequency drifts over Germany

During 2004 and 2005 measurements of mesospheric/lower thermospheric (80–100 km) winds have been carried out in Germany using three different ground-based systems, namely a meteor radar (36.2 MHz) at the Collm Observatory (51.3°N, 13°E), a MF radar (3.18 MHz) at Juliusruh (54.6°N, 13.4°E) and the LF D1 measurements using a transmitter (177 kHz) at Zehlendorf near Berlin and receivers at Collm with the reflection point at 52.1°N, 13.2°E. This provides the possibility of comparing the results of different radar systems in nearly the same measuring volume. Meteor radar winds are generally stronger than the winds observed by MF and especially by LF radars. This difference is small near 80 km but increases with height. The difference between meteor radar and medium frequency radar winds is larger during winter than during summer, which might indicate an indirect influence of gravity waves on spaced antenna measurements.