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 streams associated with Jupiter family comets

Meteor showers have been observed for a considerable time, and the cause, meteoroids from a meteoroid stream ablating in the Earth's atmosphere, has also been understood for centuries. The connection between meteoroid streams and comets was also established 150 years ago. Since that time our ability both to understand the physics and to numerically model the situation has steadily increased. We will review the current state of knowledge. However, just as there are differences between the behaviour of long period comets, Halley family comets and Jupiter family comets, so also differences exist between the associated meteoroid streams. Streams associated with Jupiter family comets show much more variety in their behaviour, driven by the gravitational perturbations from Jupiter. The more interesting showers associated with Jupiter family comets will be discussed individually.     

Meteor outbursts from long-period comet dust trails

Long-period comets have narrow one-revolution old dust trails that can cause meteor outbursts when encountered by Earth. To facilitate observing campaigns that will characterize and perhaps help find Earth-threatening, long-period comets from their trace of meteoric debris, we use past accounts of outbursts from 14 different showers to calculate the future dust trail positions near Earth’s orbit. We also examine known near-Earth, long-period comets and identify five potential new showers, which can be utilized to learn more about these objects. We demonstrate that it is the one-revolution trail that is responsible for meteor outbursts. A method that calculates in what year these showers are likely to return and at what hour is presented. The calculations improve on earlier approximate methods that used the Sun’s reflex motion to gauge the trail motion relative to Earth’s orbit.     

Dust Trails of 8P/Tuttle and the Unusual Outbursts of the Ursid Shower

We calculate the position of dust trails from comet 8P/Tuttle, in an effort to explain unusual Ursid meteor shower outbursts that were seen when the comet was near aphelion. Comet 8P/Tuttle is a Halley-type comet in a 13.6-year orbit, passing just outside of Earth's orbit. We find that the meteoroids tend to be trapped in the 12:14 mean motion resonance with Jupiter, while the comet librates in a slightly shorter period orbit around the 13:15 resonance. It takes 6 centuries to decrease the perihelion of the meteoroid orbits enough to intersect Earth's orbit, during which time the meteoroids and comet separate in mean anomaly by 6 years, thus explaining the 6-year lag between the comet's return and Ursid outbursts. The resonances also prevent dispersion along the comet orbit and limit viewing to only one year in each return. We identified past dust trail encounters with dust trails from 1392 (Dec. 1945) and 1378 (Dec. 1986) and predicted another outburst on 2000 December 22 at around 7:29 and 8:35 UT, respectively, from dust trails dating to the 1405 and 1392 returns. This event was observed from California using video and photographic techniques. At the same time, five Global-MS-Net stations in Finland, Japan, and Belgium counted meteors using forward meteor scatter. The outburst peaked at 8:06±07 UT, December 22, at zenith hourly rate ∼90 per hour, and the Ursid rates were above half peak intensity during 4.2 h. We find that most Ursid orbits do scatter around the anticipated positions, confirming the link with comet 8P/Tuttle and the epoch of ejection. The 1405 and 1392 dust trails appear to have contributed similar amounts to the activity profile. Some orbits provide a hint of much older debris being present as well. This work is the strongest evidence yet for the relevance of mean motion resonances in Halley-type comet dust trail evolution.     

Meteoroid streams and the sporadic background

Summary There is a general agreement that meteoroid streams form through the ejection of dust grains, or meteoroids, up to a few centimeters in size from comets and possibly asteroids. After ejection these meteoroids are subject to forces arising from Solar radiation and the gravitational fields of the planets. Meteoroids may also break up into smaller ones through collisions and other effects. In many cases meteor showers have been observed for millennia, with material being fed into the stream throughout this period from the parent and material lost through the external effects mentioned. Much of the lost material forms the general sporadic background. This paper will review our state of knowledge of the processes involved above and will also aim to give some insight into the structure of the sporadic background     

Venus-intercepting meteoroid streams

This study is motivated by the possibility of determining the large-body meteoroid flux at the orbit of Venus. Towards this end, we attempt to estimate the times at which enhanced meteoric activity might be observed in the planet's atmosphere. While a number of meteoroid streams are identified as satisfying common Earth and Venus intercept conditions, it is not clear from the Earth-observed data if these streams contain large-body meteoroids. A subset of the Taurid Complex objects may produce fireball-rich meteor showers on Venus. A total of 11 short-period, periodic comets and 46 near-Earth asteroids approach the orbit of Venus to within 0.1 au, and these objects may have associated meteoroid streams. Comets 27P/Crommelin and 7P/Pons–Winnecke are identified as candidate parents to possible periodic meteor showers at the orbit of Venus.     

A micrometeor component of the 1998 Leonid shower

Most astronomers expected a significant meteor shower associated with the Leonid meteoroid stream to appear in 1998 and 1999. An enhanced shower was widely observed in both years, and details can be found in many published articles. In 1998, one remarkable feature was the appearance of a strong component, rich in bright meteors, which appeared about 16 h before the expected maximum of the main shower, but another observed feature was an abnormal peak in the ionosphere characteristic value f _b E _s which was detected about 18 h after the main shower. A very high value of f _b E _s persisted for over an hour. The likely explanation is that the ionosphere was bombarded by an additional swarm of meteoroids, much smaller than those that produce a visible trail or an ionization trail that can be picked up by radio detectors. The different dynamical behaviours between small and large meteoroids are investigated and, in consequence, an explanation for the observed phenomena is offered and 1933 is suggested as being the likely ejection time.     

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.     

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.     

Searching for Light Curve Evidence of Meteoroid Structure and Fragmentation

We used light curve analysis to search for evidence of the dustball meteoroid model. Leonid, Taurid, Alpha Monocerotid and sporadic meteors from November 2003 were observed and analyzed using uniform methodology. Meteors from these four sources were examined for evidence of fragmentation by examining light curve shape and searching for light curve irregularities. Differences in meteoroid structure should be reflected by differences in meteor light curves. The resulting meteor light curve F-parameter values showed no statistically significant differences between the meteors from the various cometary showers or the sporadic meteors. The F-parameter values also suggest that the meteoroids from these sources do not follow a single body ablation model, which suggests that all four sources produce meteoroids with a dustball structure.     

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.     

Quantitative model of the release of sodium from meteoroids in the vicinity of the Sun: Application to Geminids

A considerable depletion of sodium was observed in Geminid meteoroids. To explain this phenomenon, we developed a quantitative model of sodium loss from meteoroids due to solar heating. We found that sodium can be lost completely from Geminid meteoroids after several thousands of years when they are composed of grains with sizes up to ∼100 μm. The observed variations of sodium abundances in Geminid meteor spectra can be explained by differences in the grain sizes among these meteoroids. Sodium depletions are also to be expected for other meteoroid streams with perihelion distances smaller than ∼0.2 AU. In our model, the meteoroids were represented by spherical dust-balls of spherical grains with an interconnected pore space system. The grains have no porosity and contain usual minerals known from meteorites and IDP's, including small amount of Na-bearing minerals. We modeled the sequence of three consecutive processes for sodium loss in Geminid meteoroids: (i) solid-state diffusion of Na atoms from Na-bearing minerals to the surface of grains, (ii) thermal desorption from grain surfaces and (iii) diffusion through the pore system to the space. The unknown material parameters were approximated by terrestrial analogs; the solid-state diffusion of Na in the grains was approximated by the diffusion rates for albite and orthoclase.     

The origin of the June Bootid outburst in 1998 and determination of cometary ejection velocities

Numerical integrations are used to show that the main contribution to the outburst observed in the June Bootid meteor shower in 1998 was a subset of meteoroids released from the parent comet, 7P/Pons–Winnecke, at its 1825 return. A substantial part of the June Bootid stream is in 2:1 resonance with Jupiter. This inhibits chaotic motion, allowing structures in the stream to remain compact enough over centuries that meteor outbursts can still be produced. Circumstances of ejection in 1825 are calculated that exactly result in orbits capable of producing meteors at the observed time in 1998. Required ejection velocities are  10–20 m s^-1  .     

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.     

Accuracy of meteoroid speeds determined using a Fresnel transform procedure

New methods of determining meteor speeds using radar are giving results with an accuracy of better that 1%. It is anticipated that this degree of precision will allow determinations of pre-atmospheric speeds of shower meteors as well as estimates of the density of the meteoroids. The next step is to determine under what conditions these new measurements are reliable.Errors in meteoroid speeds determined using a Fresnel transform procedure applied to radar meteor data are investigated. The procedure determines the reflectivity of a meteor trail as a function of position, by application of the Fresnel transform to the time series of a radar reflection from the trail observed at a single detection station. It has previously been shown that this procedure can be used to determine the speed of the meteoroid, by finding the assumed speed that gives a reflectivity image that best meets physical expectations. It has also been shown that speeds determined by this method agree with those from the well established “pre-t_o phase” method when applied to reflections with a high signal to noise ratio. However, there is a discrepancy between the two methods for weaker reflections. A method to investigate the discrepancy is described and applied, with the finding that the speed determined by using the Fresnel transform procedure is more accurate for weaker reflections than that given by the “pre-t_o phase” method.     

A dynamical model of the sporadic meteoroid complex

Sporadic meteoroids are the most abundant yet least understood component of the Earth's meteoroid complex. This paper aims to build a physics-based model of this complex calibrated with five years of radar observations. The model of the sporadic meteoroid complex presented here includes the effects of the Sun and all eight planets, radiation forces and collisions. The model uses the observed meteor patrol radar strengths of the sporadic meteors to solve for the dust production rates of the populations of comets modeled, as well as the mass index. The model can explain some of the differences between the meteor velocity distributions seen by transverse versus radial scatter radars. The different ionization limits of the two techniques result in their looking at different populations with different velocity distributions. Radial scatter radars see primarily meteors from 55P/Tempel-Tuttle (or an orbitally similar lost comet), while transverse scatter radars are dominated by larger meteoroids from the Jupiter-family comets. In fact, our results suggest that the sporadic complex is better understood as originating from a small number of comets which transfer material to near-Earth space quite efficiently, rather than as a product of the cometary population as a whole. The model also sheds light on variations in the mass index reported by different radars, revealing it to be a result of their sampling different portions of the meteoroid population. In addition, we find that a mass index of s=2.34 as observed at Earth requires a shallower index (s=2.2) at the time of meteoroid production because of size-dependent processes in the evolution of meteoroids. The model also reveals the origin of the 55° radius ring seen centered on the Earth's apex (a result of high-inclination meteoroids undergoing Kozai oscillation) and the central condensations seen in the apex sources, as well as providing insight into the strength asymmetry of the helion and anti-helion sources.     

The ejection velocity of meteoroids from cometary nuclei deduced from observations of meteor shower outbursts

The ejection velocities of meteoroids belonging to the Leonid and Perseid meteoroid streams are deduced from the observed differences between the longitude of the ascending node of the outburst meteoroids and that of the parent comet. The difference is very sensitive to the true anomaly of the ejection point, as well as the ejection velocity, and probable values for both are discussed.     

The P/Halley Stream: Meteor Showers on Earth, Venus and Mars

We have simulated the formation and evolution of comet 1P/Halley’s meteoroid stream by ejecting particles from the nucleus 5000 years ago and propagating them forward to the present. Our aim is to determine the existence and characteristics of associated meteor showers at Mars and Venus and compare them with 1P/Halley’s two known showers at the Earth. We find that one shower should be present at Venus and two at Mars. The number of meteors in those atmospheres would, in general, be less than that at the Earth. The descending node branch of the Halley stream at Mars exhibits a clumpy structure. We identified at least one of these clumps as particles trapped in the 7:1 mean motion resonance with Jupiter, potentially capable of producing meteor ourbursts of ZHR∼1000 roughly once per century.     

NEOCAM: The Near Earth Object Chemical Analysis Mission

The prime measurement objective of the Near Earth Object Chemical Analysis Mission (NEOCAM) is to obtain the ultraviolet spectra of meteors entering the terrestrial atmosphere from ∼125 to 300 nm in meteor showers. All of the spectra will be collected using a slitless ultraviolet spectrometer in Earth orbit. Analysis of these spectra will reveal the degree of chemical diversity in the meteors, as observed in a single meteor shower. Such meteors are traceable to a specific parent body and we know exactly when the meteoroids in a particular shower were released from that parent body (Asher, in: Arlt (ed.) Proc. International Meteor Conference, 2000; Lyytinen and van Flandern, Earth Moon Planets 82–83:149–166, 2000). By observing multiple apparitions of meteor showers we can therefore obtain quasi-stratigraphic information on an individual comet or asteroid. We might also be able to measure systematic effects of chemical weathering in meteoroids from specific parent bodies by looking for correlations in the depletions of the more volatile elements as a function of space exposure (Borovička et al., Icarus 174:15–30, 2005). By observing the relation between meteor entry characteristics (such as the rate of deceleration or breakup) and chemistry we can determine if our meteorite collection is deficient in the most volatile-rich samples. Finally, we can obtain a direct measurement of metal deposition into the terrestrial stratosphere that may act to catalyze atmospheric chemical reactions.     

On the mechanisms leading to orphan meteoroid streams

We analyse several mechanisms capable of creating orphan meteoroid streams (OMSs) for which a parent has not been identified. OMSs have been observed as meteor showers since the XIXth century and by the IRAS satellite in the 1980s. We find that the process of close encounters with giant planets (particularly Jupiter) is the most efficient mechanism to create them: only a limited section of the stream is perturbed and follows the parent body on its new orbit, while the majority of the meteoroids remain in their pre-encounter orbit or in an intermediate state, breaking the link with their parent body. Cometary non-gravitational forces can also contribute to the process since they cause the comet to drift away from its stream. However, they are not sufficient by themselves to produce an OMS. Resonances can either split or confine a stream over a long time (>1000 yr). Some meteoroid streams may look like OMSs since their parent comet is dormant or not observable (e.g. long period). Even if new techniques succeed in linking minor objects to meteoroid streams, OMSs will still exist simply because cometary nuclei are subject to complete disruption leading to their disappearance.