The Velocity Distribution of Meteoroids at the Earth as Measured by the Canadian Meteor Orbit Radar (CMOR)

The velocity distribution of meteoroids at the Earth is measured using a time-of-flight measurement technique applied to data collected by the CMOR radar (29.85 MHz). Comparison to earlier velocity measurements from the Harvard Radio Meteor Project suggests that HRMP suffered from biases which underestimated the number of fragmenting meteoroids. This bias results in a systematic underestimation of the numbers of higher velocity meteoroids. Other works (cf. Taylor and Elford, 1998) have also found additional biases in the HRMP which suggest the original HRMP meteoroid velocity analysis may have underestimated the fraction of high velocity meteors by factors up to 10⁴.     

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 Maribo CM2 meteorite fall—Survival of weak material at high entry speed

High entry speed (>25 km s^?1) and low density (<2500 kg m^?3) are the two factors that lower the chance of a meteoroid to drop meteorites. The 26 g carbonaceous (CM2) meteorite Maribo recovered in Denmark in 2009 was delivered by a bright bolide observed by several instruments across northern and central Europe. By reanalyzing the available data, we confirmed the previously reported high entry speed of (28.3 ± 0.3) km s^?1 and trajectory with slope of 31° to the horizontal. In order to understand how such a fragile material survived, we applied three different models of meteoroid atmospheric fragmentation to the detailed bolide light curve obtained by radiometers located in Czech Republic. The Maribo meteoroid was found to be quite inhomogeneous with different parts fragmenting at different dynamic pressures. While 30–40% of the (2000 ± 1000) kg entry mass was destroyed already at 0.02 MPa, another 25–40%, according to different models, survived without fragmentation up to the relatively large dynamic pressures of 3–5 MPa. These pressures are only slightly lower than the measured tensile strengths of hydrated carbonaceous chondrite (CC) meteorites and are comparable with usual atmospheric fragmentation pressures of ordinary chondritic (OC) meteoroids. While internal cracks weaken OC meteoroids in comparison with meteorites, this effect seems to be absent in CC, enabling meteorite delivery even at high speeds, though in the form of only small fragments.     

Watching meteors on Triton

The thin atmosphere of Neptune's moon Triton is dense enough to ablate micrometeoroids as they pass through. A combination of Triton's orbital velocity around Neptune and its orbital velocity around the Sun gives a maximum meteoroid impact velocity of approximately 19 km s⁻¹, sufficient to heat the micrometeoroids to visibility as they enter. The ablation profiles of icy and stony micrometeoroids were calculated, along with the estimated brightness of the meteors. In contrast to the terrestrial case, visible meteors would extend very close to the surface of Triton. In addition, the variation in the meteoroid impact velocity as Triton orbits Neptune produces a large variation in the brightness of meteors with orbital phase, a unique Solar System phenomenon.     

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 structure from radar head echoes

The Advanced Research Project Agency (ARPA) Long Range Tracking and Instrumentation Radar has recorded thousands of head echoes from small meteoroids, which include detailed trajectory information as well as ionization measurements. In total, 25 complete ionization curves have been matched using a detailed model of meteoroid ablation, though the solutions are not necessarily unique. While measurements of the spread along the trajectory of the echoes indicate that most meteors in this size range do not have large separations among fragments, the ionization curves are consistent with fragmenting bodies in the most cases. Very precise radar measurements of meteors can be a valuable source of data on the chemical and physical properties of small meteoroids.     

The hyperthermal ionization and high absolutemeteor velocities observed with HPLA radars

High Power Large Aperture (HPLA) radars generally observe very high meteor velocities averaging over 50 km s⁻¹. There are only a few events recorded around 30 km s⁻¹, while meteors at 20 km s⁻¹ or slower are very rare. This is a clear and debated contradiction to specular meteor radar results. A high plasma density condition contributes, but the dominating phenomenon is the hyperthermal ionization mechanism due to chemical dynamics of the ionization process. The observed high velocities can be explained in terms of high hyperthermal ionization cross-sections for collisions between ablated meteoroid metal atoms such as Na and/or Fe and atmospheric species.     

The first confirmed Perseid lunar impact flash

The first confirmed lunar impact flash due to a non-Leonid meteoroid is reported. The observed Perseid meteoroid impact occurred at 18^h28^m27^s on August 11, 2004 (UT). The selenographic coordinates of the lunar impact flash are 48±1° N and 72±2° E, and the flash had a visual magnitude of ca. 9.5 with duration of about 1/30 s. The mass of the impactor is estimated to have been 12 g based on a nominal model with conversion efficiency from kinetic to optical energy of 2×¹⁰⁻³. Extrapolation of a power law size-frequency distribution fitting the sub-centimeter Perseid meteoric particles to large meteoroids suggests that several flashes should have been observed at this optical efficiency. The detection of only one flash may indicate that the optical efficiency for Perseid lunar impact is much lower, or that the slope of the size distribution differs between large meteoroids and typical sub-centimeter meteoric particles.     

Sputtering and high altitude meteors

Conventional ablation theory assumes that a meteoroid undergoes intensive heating during atmospheric flight and surface atoms are liberated through thermal processes. Our research has indicated that physical sputtering could play a significant role in meteoroid mass loss. Using a 4th order Runge-Kutta numerical integration technique, we tabulated the mass loss due to the two ablation mechanisms and computed the fraction of total mass lost due to sputtering. We modeled cometary structure meteoroids with masses ranging from 10⁻¹³ to 10⁻³ kg and velocities ranging from 11.2 to 71 km s⁻¹. Our results indicate that a significant fraction of the mass loss for small, fast meteors is due to sputtering, particularly in the early portion of the light curve. In the past 6 years evidence has emerged for meteor luminosity at heights greater than can be explained by conventional ablation theory. We have applied our sputtering model and find excellent agreement with these observations, and therefore suggest that sputtered material accounts for the new type of radiation found at great heights.     

Decreased values of cosmic dust number density estimates in the Solar System

Experiments to investigate the effect of impacts on side-walls of dust detectors such as the present NASA/ESA Galileo/Ulysses instrument are reported. Side walls constitute 27% of the internal area of these instruments, and increase field of view from 140° to 180°. Impact of cosmic dust particles onto Galileo/Ulysses Al side walls was simulated by firing Fe particles, 0.5-5 μm diameter, 2-50 km s⁻¹, onto an Al plate, simulating the targets of Galileo and Ulysses dust instruments. Since side wall impacts affect the rise time of the target ionization signal, the degree to which particle fluxes are overestimated varies with velocity. Side-wall impacts at particle velocities of 2-20 km s⁻¹ yield rise times 10-30% longer than for direct impacts, so that derived impact velocity is reduced by a factor of ∼2. Impacts on side wall at 20-50 km s⁻¹ reduced rise times by a factor of ∼10 relative to direct impact data. This would result in serious overestimates of flux of particles intersecting the dust instrument at velocities of 20-50 km s⁻¹. Taking into account differences in laboratory calibration geometry we obtain the following percentages for previous overestimates of incident particle number density values from the Galileo instrument [Grün et al., 1992. The Galileo dust detector. Space Sci. Rev. 60, 317-340]: 55% for 2 km s⁻¹ impacts, 27% at 10 km s⁻¹ and 400% at 70 km s⁻¹. We predict that individual particle masses are overestimated by ∼10-90% when side-wall impacts occur at 2-20 km s⁻¹, and underestimated by ∼10-¹⁰² at 20-50 km s⁻¹. We predict that wall impacts at 20-50 km s⁻¹ can be identified in Galileo instrument data on account of their unusually short target rise times. The side-wall calibration is used to obtain new revised values [Krüger et al., 2000. A dust cloud of Ganymede maintained by hypervelocity impacts of interplanetary micrometeoroids. Planet. Space Sci. 48, 1457-1471; 2003. Impact-generated dust clouds surrounding the Galilean moons. Icarus 164, 170-187] of the Galilean satellite dust number densities of 9.4×¹⁰⁻⁵, 9.9×¹⁰⁻⁵, 4.1×¹⁰⁻⁵, and 6.8×¹⁰⁻⁵ m⁻³ at 1 satellite radius from Io, Europa, Ganymede, and Callisto, respectively. Additionally, interplanetary particle number densities detected by the Galileo mission are found to be 1.6×¹⁰⁻⁴, 7.9×¹⁰⁻⁴, 3.2×¹⁰⁻⁵, 3.2×¹⁰⁻⁵, and 7.9×¹⁰⁻⁴ m⁻³ at heliocentric distances of 0.7, 1, 2, 3, and 5 AU, respectively. Work by Burchell et al. [1999b. Acceleration of conducting polymer-coated latex particles as projectiles in hypervelocity impact experiments. J. Phys. D: Appl. Phys. 32, 1719-1728] suggests that low-density “fluffy” particles encountered by Ulysses will not significantly affect our results—further calibration would be useful to confirm this.     

A survey of meteor spectra and orbits: evidence for three populations of Na-free meteoroids

We present a survey of 97 spectra of mainly sporadic meteors in the magnitude range +3 to −1, corresponding to meteoroid sizes 1-10 mm. For the majority of the meteors, heliocentric orbits are known as well. We classified the spectra according to relative intensities of the lines of Mg, Na, and Fe. Theoretical intensities of these lines for a chondritic composition of the meteoroid and a wide range of excitation and ionization conditions were computed. We found that only a minority of the meteoroids show chondritic composition. Three distinct populations of Na-free meteoroids, each comprising ∼10% of sporadic meteoroids in the studied size range, were identified. The first population are meteoroids on asteroidal orbits containing only Fe lines in their spectra and possibly related to iron-nickel meteorites. The second population are meteoroids on orbits with small perihelia (q?0.2 AU), where Na was lost by thermal desorption. The third population of Na-free meteoroids resides on Halley type cometary orbits. This material was possibly formed by irradiation of cometary surfaces by cosmic rays in the Oort cloud. The composition of meteoroids on Halley type orbits is diverse, probably reflecting internal inhomogeneity of comets. On average, cometary dust has lower than chondritic Fe/Mg ratio. Surprisingly, iron meteoroids prevail among millimeter-sized meteoroids on typical Apollo-asteroid orbits. We have also found varying content of Na in the members of the Geminid meteoroid stream, suggesting that Geminid meteoroids were not released from their parent body at the same time.     

A technique for calculating meteor plasma density and meteoroid mass from radar head echo scattering

Large-aperture radars detect the high-density plasma that forms in the vicinity of a meteoroid and moves approximately at its velocity; reflections from these plasmas are called head echoes. To determine the head plasma density and configuration, we model the interaction of a radar wave with the plasma without using assumptions about plasma density. This paper presents a scattering method that enables us to convert measurements of radar cross-section (RCS) from a head echo into plasma density by applying a spherical scattering model. We use three methods to validate our model. First, we compare the maximum plasma densities determined from the spherical solution using 30 head echoes detected simultaneously at VHF and UHF. Second, we use a head echo detected simultaneously at VHF, UHF and L-band to compare plasma densities at all frequencies. Finally, we apply our spherical solution to 723 VHF head echoes and calculate plasma density, line density and meteoroid mass in order to compare these values with those obtained from a meteoroid ablation and ionization model. In all three comparisons, our results show that the spherical solution produces consistent results across a wide frequency range and agrees well with the single-body ablation model.     

Deconvoluting measurement uncertainty from the meteor speed distribution

Debiasing the velocity distribution of meteors observed by the Canadian Meteor Orbit Radar (CMOR) yields a distribution with large numbers of slow meteors. The distribution also contains significant numbers of hyperbolic meteors, in conflict with the expectation that interstellar meteors should be rare. In Moorhead et al. (2017a), we noted that measurement uncertainties were possibly smoothing the speed distribution and redistributing meteors to the extreme ends of the speed distribution. In this report, we use techniques analogous to image sharpening to remove the blurring caused by measurement uncertainties. The deconvolved speed distribution appears to have no meteors slower than 14 km s⁻¹ and none faster than 74 km s⁻¹. The result is to substantially raise the characteristic velocity of incoming meteoroids from 12.9 to 20.0 km s⁻¹.     

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.     

OPTICAL PREDICTIONS FOR HIGH GEOCENTRIC VELOCITY METEORS

In this study we numerically modelled the atmospheric ablation and luminosity of cometary structure meteoroids with geocentric velocities from 71 to 200 km/s. We considered meteoroid masses ranging from 10⁻¹³ to 10⁻⁶ kg. Expected heights of ablation and maximum luminosity absolute magnitudes are determined. Height and trail length values are used to calculate the angle traversed in a single video frame. It is found that for pre-atmospheric meteoroid masses of greater than 10⁻⁸ kg, high geocentric velocity meteors should be detectable with current electro-optical technology if properly optimised.     

Release of neutral sodium atoms from the surface of Mercury induced by meteoroid impacts

Meteoroid impact has been shown to be a source of sodium, and most likely of other elements, on the Moon. The same process could be also relevant for Mercury. In this work we calculate the vapor and neutral Na production rates on Mercury due to the impacts of meteoroids in the radius range of 10⁻⁸-10⁻¹ m. We limit our calculations to this size range, because meteoroids with radius larger than 10⁻¹ m have not to be found important for the daily production of the exosphere. This work is based on a new dynamical model of the meteoroid flux at the heliocentric distance of Mercury, regarding objects in the size range 10⁻²-10⁻¹ m. This size range, never investigated before, is not affected by nongravitational forces, such as the Poynting-Robertson effect, which is dominant for particles smaller than 10⁻² m. In order to evaluate the release of neutral sodium atoms also for smaller meteoroids we have used the distribution reported by M.J. Cintala [1992. Impact-induced thermal effects in the lunar and mercurian regoliths. J. Geophys. Res. 97, 947-973] calculated for particle size range 10⁻⁸-10⁻³ m. We have extrapolated this distribution up to 10⁻² m and we have based the impact calculations on a new surface composition assuming 90% plagioclase and 10% pyroxene. The results of our model are that (i) the total mass of vapor produced by the impact of meteoroids in the size range 10⁻⁸-10⁻¹ m is 4.752×¹⁰⁸ g per year, and (ii) the production rate of neutral sodium atoms is 1.5×¹⁰²² s⁻¹.     

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.     

A search for large meteoroids in the Perseid stream

Abstract— We report on two surveys conducted during the times of Perseid shower maximum in 1997 and 1998. The first survey entailed the video monitoring of the Moon's disk with the intent of recording the optical flashes that should result when large meteoroids strike the lunar surface. The second survey consisted of a combination video camera and very low frequency (VLF) radiowave receiver system capable of detecting electrophonic meteors during their ablation in the Earth's atmosphere. Using standard ablation theory, we find that for a Perseid meteoroid to be capable of generating electrophonic sounds, it must have an initial mass in excess of 495 kg. We also find, as a result of the surveys, an upper limit of 2 × 10^?17 m^?2 s^?1 to the flux of electrophonic Perseid meteors entering the Earth's atmosphere. Although our study indicates that large, meter-sized meteoroids must, at best, be sparsely distributed within the Perseid stream, we briefly discuss some tantalizing lines of evidence, found from within the astronomical literature, that hint at their true existence.     

Determining meteoroid bulk densities using a plasma scattering model with high-power large-aperture radar data

We present an improved technique for calculating bulk densities of low-mass (<1 g) meteoroids using a scattering model applied to the high-density plasma formed around the meteoroid as it enters Earth’s atmosphere. These plasmas, referred to as head echoes, travel at or near the speed of the meteoroid, thereby allowing the determination of the ballistic coefficient (mass divided by physical cross-section), which depends upon speed and deceleration. Concurrently, we apply a scattering model to the returned signal strength of the head echo in order to correlate radar-cross-section (RCS) to plasma density and meteoroid mass. In this way, we can uniquely solve for the meteoroid mass, radius and bulk density independently. We have applied this new technique to head echo data collected in 2007 and 2008 simultaneously at VHF (160 MHz) and UHF (422 MHz) at ALTAIR, which is a high-power large-aperture radar located on the Kwajalein Atoll. These data include approximately 20,000 detections with dual-frequency, dual-polarization, and monopulse (i.e. angle) returns. From 2000 detections with the smallest monopulse errors, we find a mean meteoroid bulk density of 0.9 g/cm³ with observations spanning almost three orders of magnitude from 0.01 g/cm³ to 8 g/cm³. Our results show a clear dependence between meteoroid bulk density and altitude of head echo formation, as well as dependence between meteoroid bulk density and 3D speed. The highest bulk densities are detected at the lowest altitudes and lowest speeds. Additionally, we stipulate that the approximations used to derive the ballistic parameter, in addition to neglecting fragmentation, suggest that the traditional ballistic parameter must be used with caution when determining meteoroid parameters.     

On the shape of meteor light curves, and a fully analytic solution to the equations of meteoroid ablation

The shape and characteristics (beginning and end heights, and height of maximum brightness) of meteor light curves are investigated under the constraint that the surface area S that a meteoroid presents to the oncoming air flow varies as a power law in the meteoroid mass m such that   S ∼ m ^α  . We investigate the meteoroid ablation for a range of values of α, and find that the  α= 1  condition allows for a fully analytic solution to the coupled differential equations of meteoroid ablation when the density profile is that of an isothermal atmosphere. The possible geometrical properties of Geminid meteoroids are discussed in terms of the  α= 1  ablation model and it is shown that they are consistent with being derived from an asteroidal, rather than cometary, parent body.