Analyzing Radar Meteor Trail Echoes using the Fresnel Transform Technique: A Signal Processing Viewpoint

Several techniques have been proposed for measuring speeds of meteoroids observed using radars. A recent technique involves the use of Fresnel transforms to accurately determine the speed of a meteoroid producing the trail. We follow a numerical modeling approach to analyze this technique in detail. Our studies indicate that high sensitivity to background noise levels might be a possible shortcoming of the Fresnel transform method. A matched filtering approach is presented as an alternative to alleviate this sensitivity to the noise problem. Performance of the two techniques is compared using numerical modeling and data from a 30 MHz radar.     

Plasma and Electromagnetic Simulations of Meteor Head Echo Radar Reflections

Recently, meteor head echo detections from high powered large aperture radars (HPLA) have brought new measurements to bear on the study of sporadic interplanetary meteors. These same observations have demonstrated an ability to observe smaller meteoroids without some of the geometrical restrictions of specular radar techniques. Yet incorporating data from various radar reflection types and from different radars into a single consistent model has proven challenging. 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. In order to overcome some of the unknown relationships between meteoroid characteristics (such as mass and velocity) and the resulting head echo radar cross-sections (RCS), we present our results on meteor plasma simulations of head echo plasmas 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 at speeds as large as several kilometers per second. We also present the results of finite difference time domain electromagnetic simulations (FDTD), which can calculate the radar cross-section of the simulated meteor plasma electron distributions. These simulations have shown that the radar cross-section depends in a complex manner on a number of parameters. In this paper we demonstrate that for a given head echo plasma the RCS as a function of radar frequency peaks at sqrt (2*peak plasma frequency) and then decays linearly on a dB scale with increasing radar frequency. We also demonstrate that for a fixed radar frequency, the RCS increases linearly on a dB scale with increasing head echo plasma frequency. These simulations and resulting characterization of the head echo radar cross-section will both help relate HPLA radar observations to meteoroid properties and aid in determining a particular radar facility’s ability to observe various meteoroid populations.     

Meteor velocities: A new look at an old problem

Meteoroids that orbit the Sun encounter the Earth with speeds between 11 and 74 km/sec. However, the distribution of the velocities of meteoroids between these limits is not well known. The uncertainty is caused by the difficulty in measuring the true flux of meteors at the extrema of the velocity distribution. Whilst the most comprehensive measurements of meteor flux are those obtained using radio techniques, meteors with speeds > 50 km/sec occur at heights where the effects of initial radius of the trail and diffusion significantly reduce the radio reflection from the trails; on the other hand the high dependence of the collisional ionization probability on velocity (to the power 3.5) significantly inhibits the detection of meteors with speeds < 20 km/sec. Recent developments in meteor radar systems are now making it possible to measure the velocity of meteors at the extrema of the distribution. For meteoroids ablating at heights between 100 and 120 km the speed of entry can be measured at 2 and 6 MHz using a radar with a 1 km diameter array located near Adelaide; these observations will commence early in 1995. In the meantime a 54 MHz MST radar is being operated at a pulse repetition frequency of 1024 Hz to search for the presence of interstellar (speed > 74 km/sec) meteors. Both these radars exploit the phase information available prior to the closest-approach (t_o) point.     

The initial train radius of sporadic meteors

The height distributions, velocity distributions and flux measurements of underdense echoes determined from meteor radar observations are significantly affected by the attenuation associated with the initial radius of meteor trains. Dual-frequency radar observations of a very large set of sporadic radar meteors at 29 and 38 MHz yield estimates of the initial train radius and its dependence on height and meteoroid speed as determined by the time-delay method. We provide empirical formulae that can be used to correct meteoroid fluxes for the effect of initial train radius at other radio frequencies.     

Radar Campaign to Determine the Dependence of Initial Radii of Meteor Plasma Trains on Trajectory and Orbit

Experimental and theoretical work on the transverse dimensions of meteoric plasma trains have not converged to provide generally accepted values especially uncertain is the dependence of the train radii on meteor speeds. The roles of the meteoroid structure, fragmentation and plasma processes such as ion–electron instabilities need establishing. Knowledge of the quantitative spatial distribution of plasma in meteor trains is essential for a correct interpretation of fluxes and orbital characteristics. A current project is described which employs the AMOR 26 MHz radar facility in conjunction with a frequency managed radar operating at longer wavelengths designed to measure the ionization train radii, heights, atmospheric speeds and orbits of individual meteors.     

Meteor Outburst Profiles and Cometary Ejection Models

The spatial structure of meteor streams, and the activity profiles of their corresponding meteor showers, depend firstly on the distribution of meteoroid orbits soon after ejection from the parent comet nucleus, and secondly on the subsequent dynamical evolution. The latter increases in importance as more time elapses. For younger structures within streams, notably the dust trails that cause sharp meteor outbursts, it is the cometary ejection model (meteoroid production rate as a function of time through the several months of the comet’s perihelion return, and velocity distribution of the meteoroids released) that primarily determines the shape and width of the trail structure. This paper describes how a trail cross section can be calculated once an ejection model has been assumed. Such calculations, if made for a range of ejection model parameters and compared with observed parameters of storms and outbursts, can be used to constrain quantitatively the process of meteoroid ejection from the nucleus, including the mass distribution of ejected meteoroids.     

Very small bodies in the solar system: The Earth's atmosphere as detector

Particles of mass less than about 1 gm are a minor fraction of the total matter impinging on the Earth averaged over millennia time scales. However, these particles dominate during a single particular year and produce the most obvious evidence of incoming extra-terrestrial matter in the form of ablation trails in the atmosphere which are visible at night as meteors.Observations of meteors give astronomical information on the composition, structure, and cometary associations of the particles. The composition is deduced from optical spectra of meteors, whilst telescopic studies of the trails during formation give information on the physical structure of the particles. Any cometary associations are deduced from measurement of meteor orbits determined photographically, using television, or by radar.Meteors occur in the atmosphere at heights from about 70 to 120 km. Optical observations are restricted to night-time and usually under conditions of low moonlight. A typical television based detector can record +8^M meteors with a sporadic rate of 15–20 per hour and velocities accurate to about 3%. The luminosity of the trail is strongly dependent on the velocity of the meteoroid (to about the third power).Radar observations of meteors are unrestricted by weather or time of day, and can readily detect meteors at least two orders of magnitude smaller in mass than those detectable optically. Again the observations are heavily biased toward the higher velocities as the electron line density varies approximately asV ^3.5. However, the higher the velocity of the meteoroid the greater the height of the meteor trail, and the reduced probability of radar detection due to rapid diffusion of the trail. Thus radar observations tend to select meteors in the intermediate velocity range 30–40 km s^–1.     

MARSIS surface reflectivity of the south residual cap of Mars

The south residual cap of Mars is commonly described as a thin and bright layer of CO₂-ice. The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a low-frequency radar on board Mars Express operating at the wavelength between 55 and 230 m in vacuum. The reflection of the radar wave on a stratified medium like the residual cap can generate interferences, causing weaker surface reflections compared to reflections from a pure water ice surface. In order to understand this anomalous low reflectivity, we propose a stratified medium model, which allows us to estimate both the thickness and the dielectric constant of the optically thin slab. First, we consider the residual cap as single unit and show that the decrease in the reflected echo strength is well explained by a mean thickness of 11 m and a mean dielectric constant of 2.2. This value of dielectric constant is close to the experimental value 2.12 for pure CO₂-ice. Second, we study the spatial variability of the radar surface reflectivity. We observe that the reflectivity is not homogeneous over the residual cap. This heterogeneity can be modeled either by variable thickness or variable dielectric constant. The surface reflectivity shows that two different units comprise the residual cap, one central unit with high reflectivity and surrounding, less reflective units.     

Radio and Meteor Science Outcomes From Comparisons of Meteor Radar Observations at AMISR Poker Flat, Sondrestrom, and Arecibo

Radio science and meteor physics issues regarding meteor “head-echo” observations with high power, large aperture (HPLA) radars, include the frequency and latitude dependency of the observed meteor altitude, speed, and deceleration distributions. We address these issues via the first ever use and analysis of meteor observations from the Poker Flat AMISR (PFISR: 449.3 MHz), Sondrestrom (SRF: 1,290 MHz), and Arecibo (AO: 430 MHz) radars. The PFISR and SRF radars are located near the Arctic Circle while AO is in the tropics. The meteors observed at each radar were detected and analyzed using the same automated FFT periodic micrometeor searching algorithm. Meteor parameters (event altitude, velocity, and deceleration distributions) from all three facilities are compared revealing a clearly defined altitude “ceiling effect” in the 1,290 MHz results relative to the 430/449.3 MHz results. This effect is even more striking in that the Arecibo and PFISR distributions are similar even though the two radars are over 2,000 times different in sensitivity and at very different latitudes, thus providing the first statistical evidence that HPLA meteor radar observations are dominated by the incident wavelength, regardless of the other radar parameters. We also offer insights into the meteoroid fragmentation and “terminal” process.     

A meteoroid stream survey using the Canadian Meteor Orbit Radar: I. Methodology and radiant catalogue

Using a meteor orbit radar, a total of more than 2.5 million meteoroids with masses ∼10⁻⁷ kg have had orbits measured in the interval 2002-2006. From these data, a total of 45 meteoroid streams have been identified using a wavelet transform approach to isolate enhancements in radiant density in geocentric coordinates. Of the recorded streams, 12 are previously unreported or unrecognized. The survey finds >90% of all meteoroids at this size range are part of the sporadic meteoroid background. A large fraction of the radar detected streams have q<0.15 AU suggestive of a strong contribution from sungrazing comets to the meteoroid stream population currently intersecting the Earth. We find a remarkably long period of activity for the Taurid shower (almost half the year as a clearly definable radiant) and several streams notable for a high proportion of small meteoroids only, among these a strong new shower in January at the time of the Quadrantids (January Leonids). A new shower (Epsilon Perseids) has also been identified with orbital elements almost identical to Comet 96P/Machholz.     

Meteor head echo radar data: Mass-velocity selection effects

High-power, large-aperture (HPLA) radars detect the plasma that forms in the vicinity of a meteoroid and moves approximately at its velocity; reflections from these plasmas are called head echoes. For over a decade, HPLA radars have been detecting head echoes with peak velocity distributions >50 km/s. These results have created some controversy within the field of meteor physics because previous data, including spacecraft impact cratering studies, optical and specular meteor data, indicate that the peak of the velocity distribution to a set limiting mass should be <20 km/s [Love, S.G., Brownlee, D.E., 1993. Science 262, 550-553]. Thus the question of whether HPLA radars are preferentially detecting high-velocity meteors arises. In this paper we attempt to address this question by examining both modeled and measured head echo data using the ALTAIR radar, collected during the Leonid 1998 and 1999 showers. These data comprise meteors originating primarily from the North Apex sporadic meteor source. First, we use our scattering theory to convert measured radar-cross-section (RCS) to electron line density and mass, as well as to convert modeled electron line density and mass to RCS. We subsequently compare the dependence between mass, velocity, mean-free-path, RCS and line density using both the measured and modeled data by performing a multiple, linear regression fit. We find a strong correlation between derived mass and velocity and show that line density is approximately proportional to mass times velocity^3.1. Next, we determine the cumulative mass index using subsets of our data and use this mass index, along with the results of our regression fit, to weight the velocity distribution. Our results show that while there does indeed exist a bias in the measured head echo velocity distribution, it is smaller than those calculated using traditional specular trail data due to the different scattering mechanism, and also includes a bias against the low-mass, very high-velocity meteoroids.     

Meteoroid Bulk Density Determination Using Radar Head Echo Observations

We present the results of a study of meteoroid bulk densities determined from meteor head echoes observed by radar. Meteor observations were made using the Advanced Research Projects Agency Long-Range Tracking And Instrumentation Radar (ALTAIR). ALTAIR is particularly well suited to the detection of meteor head echoes, being capable of detecting upwards of 1000 meteor head echoes per hour. Data were collected for 19 beam pointings and are comprised of approximately 70 min. of VHF observations. During these observations the ALTAIR beam was directed largely at the north apex sporadic source. Densities are calculated using the classical physical theory of meteors. Meteoroid masses are determined by applying a full wave scattering theory to the observed radar cross-section. Observed meteoroids are predominantly in the 10⁻¹⁰ to 10⁻⁶ kg mass range. We find that the vast majority of meteoroid densities are consistent with low density, highly porous objects as would be expected from cometary sources. The median calculated bulk density was found to be 900 kg/m³. The orbital distribution of this population of meteoroids was found to be highly inclined.     

Solar flare activity: Evidence for large-scale changes in the past

An analysis of radar and photographic meteor data and of spacecraft meteoroid penetration data indicates that there probably has not been a large increase in meteoroid impact rates in the last 10⁴ yr. The solar flare tracks observed in the glass linings of meteoroid impact pits on lunar rock 15205 are therefore reanalyzed assuming a meteoroid flux that is constant in time. Based on this assumption, the data suggest that the production rate of Fe-group solar flare tracks may have varied by as much as a factor of 50 on a time scale of about 10⁴ yr. No independently obtained data are known to require conflict with this interpretation. Confidence in this conclusion is somewhat qualified by the experimental and analytical uncertainties involved, but the conclusion nevertheless remains the present “best” explanation for the observed data trends.     

THE METEOR FLUX: IT DEPENDS HOW YOU LOOK

In this paper, we use radar observations from a 50 MHz radar stationed near Salinas, Puerto Rico, to study the variability of specular as well as non-specular meteor trails in the E-region ionosphere. The observations were made from 18:00 to 08:00 h AST over various days in 1998 and 1999 during the Coqui II Campaign [Urbina et al., 2000, Geophys. Rev. Lett. 27, 2853–2856]. The radar system had two sub-arrays, both produced beams pointed to the north in the magnetic meridian plane, perpendicular to the magnetic field, at an elevation angle of approximately 41 degrees. The Coqui II radar is sensitive to at least two types of echoes from meteor trails: (1) Specular reflections from trails oriented perpendicular to the radar beam, and (2) scattering, or, non-specular reflections, from trails deposited with arbitrary orientations. We examine and compare the diurnal and seasonal variability of echoes from specular and non-specular returns observed with the Coqui II radar. We also compare these results with meteor head echo observations made with the Arecibo 430 MHz radar. We use common region observations of these three types of meteor echoes to show that the diurnal and seasonal variability of specular trails, non-specular trails, and head echoes are not equivalent. The implications of these results on global meteor mass flux estimates obtained from specular meteor observations remains to be examined.     

The meteor radar response function: Theory and application to narrow beam MST radar

The meteor radar response function is an important tool for analyzing meteor backscatter observed by radar systems. We extend previous work on the development of the response function to include a non-uniform meteor ionization profile, provided by meteor ablation theory, in contrast to what has been assumed in the past. This has the advantage that the height distribution of meteors expected to be observed by a radar meteor system may be accurately modeled. Such modeling leads to meteor height distributions that have implications for the composition of those meteoroids ablating at high altitudes which may be observed by “non-traditional” meteor radars operating at MF/HF. The response function is then employed to investigate meteor backscatter observed by narrow beam MST radars which in recent years have been used increasingly to observe meteors.     

Real-Time Flux Density Measurements of the 2011 Draconid Meteor Outburst

During the 2011 outburst of the Draconid meteor shower, members of the Video Meteor Network of the International Meteor Organization provided, for the first time, fully automated flux density measurements in the optical domain. The data set revealed a primary maximum at 20:09 UT ± 5 min on 8 October 2011 (195.036° solar longitude) with an equivalent meteoroid flux density of (118 ± 10) × 10^?3/km²/h at a meteor limiting magnitude of +6.5, which is thought to be caused by the 1900 dust trail. We also find that the outburst had a full width at half maximum of 80 min, a mean radiant position of α = 262.2°, δ = +56.2° (±1.3°) and geocentric velocity of v_geo = 17.4 km/s (±0.5 km/s). Finally, our data set appears to be consistent with a small sub-maximum at 19:34 UT ±7 min (195.036° solar longitude) which has earlier been reported by radio observations and may be attributed to the 1907 dust trail. We plan to implement automated real-time flux density measurements for all known meteor showers on a regular basis soon.     

Meteor Showers Associated with the Taurid Complex Asteroids

Recent theoretical and observational work has shown that the asteroids belonging to the Taurid meteoroid complex have a cometary nature. If so, then they might possess related meteoroid streams producing meteor showers in the Earth atmosphere. We studied the orbital evolution of ten numbered Taurid complex asteroids by the Halphen-Goryachev method. It turned out that all of these asteroids are quadruple crossers relative to the Earth's orbit. Therefore their proposed meteoroid streams may in theory each produce four meteor showers. The theoretical orbital elements and geocentric radiants of these showers are determined and compared with the available observational data. The existence of the predicted forty meteor showers of the ten Taurid complex asteroids is confirmed by a search of the published catalogues of observed meteor shower radiants and orbits, and of the archives of the IAU Meteor Data Center (Lund). The existence of meteor showers associated with the Taurid Complex Asteroids confirms that, most likely, these asteroids are extinct comets. This revised version was published online in July 2006 with corrections to the Cover Date.     

Techniques For Measuring Radar Meteor Speeds

We compare the results from the application of four different methods to determine the speed of meteoroids from single station radar data. The methods used are the pre-t ₀ amplitude, post-t ₀ amplitude, pre-t ₀ phase and the Fresnel transform (FT) methods. Speeds from the first three methods are compared to the FT method since, requiring the use of the entire records of both the amplitude and phase data, this method is the most accurate of the four.     

Geminid Meteor shower activity 2003–2005 as observed by Gadanki radar

The distribution of meteor signals reflected from a backscatter radar is considered according to their duration. This duration time (T) is used to classify the meteor echoes and to calculate the mass index (S) of different meteoroids of shower plus sporadic background. Observational data on particle size distribution of the Geminid meteor shower are very scarce, particularly at low latitudes. In this paper the observational data from Gadanki radar (13.46°N, 79.18°E) have been used to determine the particle size distribution and the number density of meteoroids inside the stream of the Geminid meteor shower. The mean variation of meteor number density across the stream has been determined for three echo duration classes, T<0.4, T=0.4–1 and T>1 s. We are more interested in the appearance of echoes of various durations and therefore meteors of various masses in order to understand more on the filamentary structure of the stream. It is observed that the faint particle flux peaks earlier than the larger particles. We found a decreasing trend in the mass index values from the day of peak activity to the next observation days. The mass index profile was found to be U-shaped with a minimum value near the time of peak activity. The observed minimum s values are 1.64±0.05 and 1.65±0.04 in the years 2003 and 2005, respectively. The activity of the shower indicates the mass segregation of meteoroids inside the stream. Our results are best comparable with the “scissors” structure model of the meteoroid stream formation of Ryabova [2007. Mathematical modeling of the Geminid meteoroid stream. Mon. Not. R. Astron. Soc. 375, 1371–1380] by considering the asteroid 3200 Phaethon as an extinct comet.     

CAMS: Cameras for Allsky Meteor Surveillance to establish minor meteor showers

First results are presented from a newly developed meteoroid orbit survey, called CAMS – Cameras for Allsky Meteor Surveillance, which combines meteor detection algorithms for low-light video observations with traditional video surveillance tools. Sixty video cameras at three stations monitor the sky above 31° elevation. Goal of CAMS is to verify meteor showers in search of their parent comets among newly discovered near-Earth objects.This paper outlines the concept of operations, the hardware, and software methods used during operation and in the data reduction pipeline, and accompanies the data release of the first batch of meteoroid orbits. During the month of November 2010, 2169 precisely reduced meteoroid trajectories from 17 nights have an error in the apparent radiant of the trajectory <2° and error in speed <10%. Median values of the error are 0.31° and 0.53 km/s, respectively, sufficient to resolve the intrinsic dispersion of annual meteor showers and resolve minor showers from the sporadic background. The limiting visual magnitude of the cameras is +5.4, recording meteors of +4 magnitude and brighter, bright enough to stand out from the mostly fainter sporadic meteors detected as under dense radar echoes.CAMS readily detected all established showers (6) active during the clear nights in November. Of the showers that needed confirmation, we confirm the theta Aurigids (THA, IAU#390), the chi Taurids (CTA, IAU#388), and the omicron Eridanids (OER, IAU#338). We conclude that the iota November Aurigids (IAR, IAU#248) are in fact the combined activity of the theta Aurigids and chi Taurids, and this shower should be dismissed from the list. Finally, there is also a clustering consistent with the zeta Cancrids (ZCN, IAU#243), but we cannot exclude that this is lower perihelion dust belonging to the Orionid shower.Data are submitted to the IAU Meteor Data Center on a semi-regular basis, and can be accessed also at http://cams.seti.org.