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

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

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.     

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

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

A new narrow beam Doppler radar at 3 MHz for studies of the high-latitude middle atmosphere

A Doppler radar at 3.17 MHz has been installed at Saura close to the Andøya Rocket Range as part of the ALOMAR observatory at Andenes, Norway in summer 2002 to improve the ground based capabilities for measurements of small scale features and turbulence in the mesosphere. The main feature of the new Saura MF radar is the transmitting/receiving antenna which is arranged as a Mills Cross of 29 crossed half-wave dipoles with a minimum beam width of about 7°. Each dipole is fed by its own transceiver, and the individual phase control of the 58 transceiver modules on transmission and reception provides high flexibility in beam forming and pointing as well as transmission switching between ordinary and extraordinary mode circular polarisation. In addition, beams with different widths at the same pointing angle can be formed. For multiple receiver applications (spaced antenna wind measurements, all-sky meteor detections) four independent receiving channels are available.     

Creating an artificial Geminid meteor shower: Correlation between ejecta velocity and observability

One of the interesting arguments for a space impact mission to asteroid 3200 Phaethon is to create an artificial Geminid meteor shower. In this work we investigate the artificial shower’s dates of observability and dependence on ejecta velocity using dust trail theory. We find that when the dust ejecta velocities are 200 m/s the artificial meteor showers start to be visible in 2204 and continue for about 30 years. If the dust ejecta velocity is 20 m/s they only last 10 years from 2215 to 2225. Thus, the onset of artificial shower activity begins sooner and lasts longer with higher ejecta velocities. To produce an artificial meteor shower with 3200 Phaethon as the parent will require higher impact energy than the Deep Impact spacecraft delivered to 9P/Tempel 1.     

廊坊中频雷达流星观测及初步结果

中频雷达用来开展夜间100km高度以上的流星观测,获得流星随时间、高度、方位的分布情况及流星体速度、流星辐射点、流星余迹径向速度等参数,其探测数据可用于流星天文学、中层大气动力学等领域的研究.利用2017年11月16日12:00UT-22:00UT期间廊坊观测站（39.4°N,116.7°E）的中频雷达数据,首次开展了中国中纬度地区夜间流星观测实验,共检测到94个流星回波信号,集中分布在97～115km高度范围内,平均高度为106.5km,计算得到了流星回波的双极扩散系数、方位分布等相关参数,并与国外中频雷达流星探测结果进行了初步比较.      

Meteoroid bombardment as a source of the lunar exosphere

The column densities of impact-produced metal atoms in the exosphere during the peaks of activity of the main meteor showers – Geminids, Quadrantids and Perseids – and during quiet periods are estimated. The Na supply rate is estimated to be 2 × 10⁴, 3 × 10³, 10⁴, and 2 × 10⁴ atoms cm⁻² s⁻¹ for sporadic meteoroids, Perseid, Geminid, and Quadrantid meteor showers, respectively. A low upper limit on Ca in the lunar exosphere is explained by the condensation of Ca into dust grains during expansion of the cooling impact-produced vapor cloud. The chemical composition of gas-phase species released to the lunar exosphere during meteoroid impacts has been estimated. Most impact-produced molecules that contain metals are destroyed by solar photons while on ballistic trajectories. Energies of Na, K, Ca, and Mg atoms produced via photolysis of the respective monoxides are estimated to be 0.4, 0.35, 0.6, and 0.45 eV, respectively. The relative content of impact-produced Na and K atoms is maximal at altitudes of about 1000–2000 km and during the main meteor showers, lunar eclipses, and passages of the Moon through the Earth’s magnetosphere.     

Asteroid and comet hazard: Identification problem of observed space objects with the parental bodies

This article focuses on the genetic identification of observed small cosmic bodies with alleged parental bodies; namely, comets, asteroids and meteoroid swarms. There is a problem of the upper D-value limit as a measure of proximity between the orbits of the bodies in the five-dimensional phase space (Southworth and Hawkins, 1963). In the study of genetic relationships of the comet and meteor complexes, the D value is usually taken as equal to 0.2 for all meteor showers. However, the upper D limit should be investigated for each meteoroid complex. For example, such investigation was performed for the Taurid meteor complex (Porub?an et al., 2006). In this paper, the upper D-criterion limit value was investigated for the Perseid meteor shower. The 1862 III Swift–Tuttle comet is its parental comet.     

An intercomparison of meteor radar measurements at the South Pole using two different processing systems

A new meteor radar system was installed at the Amundsen–Scott station South Pole in 2001 to further the understanding of the dynamics of the Antarctic region. The antenna array consists of four yagis pointed along the 0°, 90°, 180°, and 270° meridians and five folded crossed dipoles arranged in a cross configuration and operating as an interferometer to provide position measurements for the detected radio meteors. The four yagis are time division multiplexed and used for both transmitting and receiving while the five folded crossed dipoles are only used for reception. The current arrangement of data acquisition (DAQ) systems at the South Pole allows the collection of meteors in a configuration similar to the previous meteor radar system that operated at the South Pole in the mid 1990s while also using an interferometer to accurately determine the meteor positions in the sky, which enables the determination of the vertical structure of the observed waves. This has been accomplished through the use of two DAQ and post-processing systems: COBRA (Colorado Obninsk radar) connected to the yagis and MEDAC (meteor echo detection and collection) connected to the folded crossed dipoles. With two separate DAQ systems operating in parallel we have the ability to directly compare the results and understand the inherent variability in the derived scientific results based on different system architectures and processing assumptions. The impact of operating a system without an interferometer on the amplitudes and phases of the observed wave components is considered. We find that the lack of altitude resolution of the COBRA DAQ system leads to an underestimation of the amplitude of the s = 1 component of the semidiurnal tide of ∼20% during the summer months.     

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.     

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

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

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.     

Redefinition of the meteor storm from the point of view of spaceflight security

ZHR¹ is defined as number of meteoroids passing 1000 km² zenith area per hour which can produce craters no less than 1 cm in diameter on an aluminum material. The relationship between ZHR and ZHR¹ is deduced. We evaluate the strong meteor showers since 1990s by ZHR¹ instead of ZHR and find the Giacobinid (Draconid) shower in 1998 is much stronger than the Leonid shower in 1999, 2001 and 2002.     

Japanese research project on Arctic and Antarctic observations of the middle atmosphere

An all-sky optical imager is in routine observation at the South Pole. Monochromatic images of aurora and air glow at N₂⁺ 427.8nm, OI 557.7nm, OI 630nm and OH 730nm are supplying significant information on the magnetospheric process in the polar cap and cusp/cleft region along with atmospheric wave signature at this particular point. SuperDARN radars in Antarctica make observations over the South Pole.

At Syowa Station, Antarctica, a multi-instrumental observation project is now being implemented for the study of the polar upper atmosphere from the mesosphere to the thermosphere, where complex physical and chemical processes take place making the region very attractive for scientific research. Two HF radars, which are part of SuperDARN radars, have been already installed and started observations. By the end of 1999, all-sky imagers, photo meters, a Na temperature Lidar, an MF radar and a Fabry-Perot interferometer will be introduced and start collecting various physical parameters on a routine basis.

In the Arctic region, we are planning to deploy coordinated ground-based observations with optical, radio and radar sensing of the polar middle and upper atmosphere in conjunction with EISCAT radars. Scientific goals are versatile to shed light on the tangled coupling processes in response to magnetospheric disturbances from above and bi-lateral interactions with high-density lower atmospheric layers. These are outlined in this paper.     

Gravity waves and wind-shear in the MLT at 23°S

We have used the technique suggested by Hocking [Hocking, W. A new approach to momentum flux determinations using SKiYMET meteor radars. Ann. Geophys. 23, 2005.] to derive short period wind variances in the 80–100 km region from meteor radar data. We find that these fluctuating winds, assumed to correspond to gravity waves and turbulence, are closely correlated with the vertical shear of the horizontal tidal winds. This close correlation suggests that in situ wind shear may be a major source of gravity waves and turbulence in the MLT. If this is the case, gravity waves generated in the troposphere and propagating up to the MLT region, generally assumed to constitute an important influence on the climatology of the region, may be a less important source of energy and momentum in the 80–100 km region than has been hitherto believed.     

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.     

Measurement of turbulent kinetic energy dissipation rates in the mesosphere by a 3 MHz Doppler radar

A new narrow beam Doppler radar operating at 3.17 MHz has been installed close to the Andøya Rocket Range in Andenes, Norway in summer 2002 in order to improve the ground based capabilities for measurements of turbulence in the mesosphere. The main feature of the radar is a Mills Cross transmitting/receiving antenna consisting of 29 crossed half-wave dipoles. In combination with the modular transceiver system this provides high flexibility in beam forming and pointing. In general, vertical and oblique beams with a minimum one way half-power full-beam width (HPFW) of 6.6° are used. The observations are usually performed with a height resolution of 1 km and with off-zenith beams at 7.3° directed towards NW, NE, SE, and SW. Turbulence intensities have been estimated from the width of the observed signal spectra using an computationally intensive correction method which requires precise knowledge of the antenna radiation pattern. The program uses real-time measurements of the wind field in all determinations. Turbulent kinetic energy dissipation rates based on radar observations are presented and compared with corresponding climatological summer and winter profiles from rocket measurements, as well as with single profiles from model runs for selected periods from September 2003 to Summer 2004. The mean turbulent kinetic energy dissipation rates based on these radar measurements are about 5 mW/kg at 60 km altitude and about 20 mW/kg at 80 km, in reasonable agreement with mean turbulence intensities obtained from previous rocket soundings at Andenes.     

Debris flux comparisons from the Goldstone Radar,Haystack Radar,and Hax Radar prior,during, and after the last solar maximum

The continual monitoring of the low Earth orbit (LEO) debris environment using highly sensitive radars is essential for an accurate characterization of these dynamic populations. Debris populations are continually evolving since there are new debris sources, previously unrecognized debris sources, and debris loss mechanisms that are dependent on the dynamic space environment. Such radar data are used to supplement, update, and validate existing orbital debris models. NASA has been utilizing radar observations of the debris environment for over a decade from three complementary radars: the NASA JPL Goldstone radar, the MIT Lincoln Laboratory (MIT/LL) Long Range Imaging Radar (known as the Haystack radar), and the MIT/LL Haystack Auxiliary radar (HAX). All of these systems are highly sensitive radars that operate in a fixed staring mode to statistically sample orbital debris in the LEO environment. Each of these radars is ideally suited to measure debris within a specific size region. The Goldstone radar generally observes objects with sizes from 2 mm to 1 cm. The Haystack radar generally measures from 5 mm to several meters. The HAX radar generally measures from 2 cm to several meters. These overlapping size regions allow a continuous measurement of cumulative debris flux versus diameter from 2 mm to several meters for a given altitude window. This is demonstrated for all three radars by comparing the debris flux versus diameter over 200 km altitude windows for 3 nonconsecutive years from 1998 to 2003. These years correspond to periods before, during, and after the peak of the last solar cycle. Comparing the year to year flux from Haystack for each of these altitude regions indicate statistically significant changes in subsets of the debris populations. Potential causes of these changes are discussed. These analysis results include error bars that represent statistical sampling errors.     

Airglow OH emission height inferred from the OH temperature and meteor trail diffusion coefficient

Simultaneous observations of the airglow OH(6,2) band rotational temperature, TOH, and meteor trail ambipolar diffusion coefficient, D, were carried out at Shigaraki (35°N, 136°E), during PSMOS 2003 Campaign, January 28 to February 8, 2003. The OH emission height was estimated by cross correlation analysis of the TOH and D nocturnal variations. A good correlation between TOH and D was obtained at 85 km of altitude. From the nocturnal variations of TOH and D, it is found that the OH emission peak height varied from 88 km before the midnight to 84 km in the early morning. The height variation could be caused by an atmospheric tidal effect in the emission height.