Fragmentation model analysis of the observed atmospheric trajectory of the Tagish Lake fireball

Abstract— A recently published meteoroid fragmentation model (FM) was applied to observational data on the Tagish Lake meteoric fireball. An initial mass of 56,000 kg, derived from seismic and infrasound data by Brown et al. (2002), proved to be consistent with a very low value of intrinsic ablation coefficient of 0.0009 s² km^?2. The average residual of the best fit to the observed light curve was ±0.10 stellar magnitude. The apparent ablation coefficient varied from 0.0009 to 1.52 s² km^?2 with an average value of 0.054 s² km^?2 (determined by the gross fragmentation [GF] model). The FM found 33 individual fragmentation events during the penetration of the 56,000 kg initial mass of the Tagish Lake meteoroid through the atmosphere, with five of the events fragmenting more than 10% of the instantaneous mass of the main body. The largest event fragmented 88% of the mass of the main body at a height of 34.4 km. The velocity of the main body mass of 2660 kg at a height of 29.2 km (the last observed light) was 13.1 km/s. Strong fragmentation at heights lower than 29.2 km is very probable. The extreme fragmentation process of the Tagish Lake meteoroid puts its classification well outside the IIIB type in the direction of less cohesive bodies. The light curve could not be explained at all by making use of only the apparent ablation coefficient and apparent luminous efficiency.     

Recent Advances in Bolide Entry Modeling:A Bolide Potpourri*

In this paper, we will review recent research on numerous aspects of bolide entry into a planetary atmosphere, including such topics as the entry dynamics, energetics, ablation, deceleration, fragmentation, luminosity, mechanical wave generation processes, a total (panchromatic) power budget including differential and integral efficiencies vs. time, etc. Fragmentation, triggered by stagnation pressures exceeding the bolide breaking strength, has been included with subsequent wake behavior in either a collective or non-collective behavior limit. We have also utilized the differential panchromatic luminous efficiency of ReVelle and Ceplecha (2002c, Proceedings of Asteroids, Comets, Meteors ACM 2002, 29 July–2August, 285–288) to compute bolide luminosity. In addition we also introduce the concept of the differential and integral acoustic/infrasonic efficiency and generalized it to the case of mechanical wave efficiency including internal atmospheric gravity waves generated during entry. Unlike the other efficiencies which are assumed to be a constant multiple of the luminous efficiency, the acoustic efficiency is calculated independently using a 'first principles' approach. All of these topics have been pursued using either a homogeneous or a porous meteoroid model with great success. As a direct result, porosity seems to be a rather good possibility for explaining anomalous meteoroid behavior in the atmosphere.     

The beginning heights and light curves of high‐altitude meteors

Abstract— In this paper, we provide an overview of meteors with high beginning height. During the recent Leonid meteor storms, as well as within the regular double station video observations of other meteor showers, we recorded 164 meteors with a beginning height above 130 km. We found that beginning heights between 130 and 150 km are quite usual, especially for the Leonid meteor shower. Conversely, meteors with beginning heights above 160 km are very rare even among Leonids. From the meteor light curves, we are able to distinguish two different processes that govern radiation of the meteors at different altitudes. Light curves vary greatly above 130 km and exhibit sudden changes in meteor brightness. Sputtering from the meteoroid surface is the dominating process during this phase of the meteor luminous trajectory. Around 130 km, the process switches to ablation and the light curves become similar to the light curves of standard meteors. The sputtering model was successfully applied to explain the difference in the beginning heights of high‐altitude Leonid and Perseid meteors. We show also that this process in connection with high altitude fragmentation could explain the anomalously high beginning heights of several relatively faint meteors.     

Recent Advances in Bolide Entry Modeling: A Bolide Potpourri

In this paper, we will review recent research on numerous aspects of bolide entry into a planetary atmosphere, including such topics as the entry dynamics, energetics, ablation, deceleration, fragmentation, luminosity, mechanical wave generation processes, a total (panchromatic) power budget including differential and integral efficiencies versus time, etc. Fragmentation, triggered by stagnation pressures exceeding the bolide breaking strength, has been subsequently included in either a collective or non-collective wake behavior limit. We have also utilized the differential panchromatic luminous efficiency of ReVelle and Ceplecha (2002) to compute bolide luminosity. In addition we also introduce the concept of the differential and integral acoustic/infrasonic efficiency and generalized it to the case of mechanical wave efficiency including internal atmospheric gravity waves generated during entry. Unlike the other efficiencies which are assumed to be a constant multiple of the luminous efficiency, the acoustic efficiency is calculated independently using a “first principles” approach. All of these topics have been pursued using either a homogeneous or a porous meteoroid model with great success. As a direct result, porosity seems to be a rather good possibility for explaining anomalous meteoroid behavior in the atmosphere. Invited Paper Presented at Meteoroids 2004; Presented at University of Western Ontario, London, Ontario, Canada, August 16–20, 2004     

Approximation of the observed motion of bolides by the analytical solution of the equations of meteor physics

A great volume of data has been accumulated thus far related to the photoregistration of the paths of meteor bodies in the terrestrial atmosphere. Most images have been obtained by four bolide networks, which operate in the USA, Canada, Europe, and Spain in different time periods. The approximation of the actual data using theoretical models makes it possible to achieve additional estimates, which do not directly follow from the observations. In the present study, we suggest an algorithm to find such parameters of the theoretical relationship between the height and the velocity of the bolide motion that help to fit observations along the luminous part of the trajectories in the best way. The main difference from previous studies is that the given observations are approximated using the analytical solution of the equations of meteor physics. The model presented in this study was applied to a number of bright meteors observed by the Canadian camera network and by the US Prairie network and to the Benésov bolide, which is one of the largest fireballs registered by the European network. The correct mathematical modeling of meteor events in the atmosphere is necessary for further estimates of the key parameters, including the extra-atmospheric mass, the ablation coefficient, and the effective enthalpy of evaporation of entering bodies. In turn, this information is needed by some applications, namely, those aimed at studying the problems of asteroid and comet security, to develop measures of planetary defense, and to determine the bodies that can reach Earth’s surface.     

Detailed data for 259 fireballs from the Canadian camera network and inferences concerning the influx of large meteoroids

Abstract— We present data for 259 meteoric fireballs observed with the Canadian camera network, including velocities, heights, orbits, luminosities along each trail, estimates of preatmospheric masses and surviving meteorites (if any) as well as membership in meteor showers. Some 213 of the events comprise an unbiased sample of the 754 fireballs observed in a total of 1.51 × 10¹⁰ km² h of clear-sky observations. The number of fireballs and the amount of clear sky in which they were recorded are given for each day of the year. We find at least 37% of the unbiased sample are members of some 15 recognized meteor showers. Preatmospheric masses, based on an assumed luminous efficiency of 0.04 for velocities >10 km s^?1, range from 1 g for some very fast fireballs up to hundreds of kilograms for the largest events. We present plots and equations for the flux, as a function of initial mass, for the entire group of fireballs and for some subgroups: meteorite-dropping objects; meteor shower members; groups that appear to be mainly of asteroidal or cometary origin; and for very fast objects. For masses of a few kilograms, asteroidal objects outnumber cometary ones. Cometary objects attain greater peak brightness than asteroidal ones of equal mass largely due to higher velocity, but also because they fragment more severely. For 66 fireballs, we estimate the meteoroid density using photometric and dynamic masses. Presumed cometary objects have typical densities near 1.0, while asteroidal values show two groups that suggest meteoroids similar to carbonaceous and ordinary chondrites. Our basic data may be used by others for further studies or to reexamine our results using assumptions different from those employed in this paper.     

Airborne intensified charge-coupled device observations of the 1998 Leonid shower

Abstract— We have used dual coaxial microchannel plate image-intensified monochrome charge-coupled device (CCD) detectors run at standard NTSC frame rates (30 frames per second, fps) to study the Leonid meteor shower on 1998 November 17 from an airborne platform at an altitude of ～13 km. These observations were part of NASA's 1998 Leonid multi-instrument aircraft campaign (MAC). The observing systems had fields of view (width) of 16.3° and 9.5°, and limiting stellar sensitivities of +8.3^m and +8.9^m. During 12 h of recording, 230 meteors were detected, of which 65 were Leonid meteors. Light curves are presented for 53 of these meteors. The magnitudes at peak brightness of the meteors investigated were generally in the range from +4.0^m to +6.0^m. The mass distribution indices for the two samples are 1.67 and 1.44, with the former being based on the wider field of view dataset. The light curves were skewed with the brightest point towards the beginning of the meteor trail. The F parameter for points one magnitude below maximum luminosity had a mean value of 0.47 for the wider field system and 0.37 for the more sensitive narrower field system. We provide leading and trailing edge light curve slopes for each meteor as another indication of light curve shape. There were few obvious flares on the light curves, indicating that in-flight fragmentation into a large number of grains is not common. There is variability in light curve shape from meteor to meteor. The light curves are inconsistent with single, compact body meteor theory, and we interpret the data as indicative of a two-component dustball model with metal or silicate grains bonded by a lower boiling point, possibly organic, substance. The variation in light curve shape may be indicative of differences in mass distribution of the constituent grains. We provide trail length vs. magnitude data. There is only a slight hint of a bend at +5^m in the data, representing the difference between meteors that have broken into a cluster of grains prior to grain ablation, and those that continue to fragment during the grain ablation phase. Two specific meteors show interesting light curve features. One meteor is nebulous in appearance, with significant transverse width. The apparent light production region extends for 450 m from the center of the meteor path. Another meteor has several main fragments, and evidence of significant separated fragments. We offer several suggestions for improvements for the 1999 Leonid MAC light curve experiment.     

Reanalysis of the Historic AFTAC Bolide Infrasound Database

We have recently digitized and partially reanalyzed the historic bolide infrasonic database. These 10 events were originally detected by the U.S. Air Force Technical Applications Center (AFTAC) from ∼1960 to 1974. In this paper we present the first preliminary reanalysis results for two of the 10 bolide events, namely the Revelstoke bolide of 3/31/1965 as well as the Prince Edward Islands (P.E.I). S. African bolide of 8/03/1963, which were among the largest bolides detected during this time period. These bolides have been investigated initially since they are most likely to have had a significant effect on the computed global influx rate of ReVelle (Global Infrasonic Monitoring of Large Bolides, pp 483–490, 2001) as indicated in Brown et al. (Nature, 420:314–316, 2002). We are in the process of recomputing all relevant infrasonic propagation quantities such as plane wave back azimuth, signal velocities, power spectra, spectrograms, as well as energy estimates using multiple techniques. In a future paper we will present a complete digital reanalysis of the AFTAC bolide infrasonic data and its final resulting global bolide influx implications.     

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.     

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.     

Peculiarities of the Fragmentation of Benešov's Bolide

One of the previously published models for the fragmentation of meteoroids (Baldwin and Sheaffer, 1971) is used to analyze the observed segment of the trajectory of Beneov's bolide (Beneov, EN 070591) of the European Fireball Network. An important element of these observations is the separation of luminous fragments from the main body at the final radiation stage of the bolide as a whole. The application of the least-squares method with the use of various models as trial functions shows that the model of a single body with ablation provides the best fit to the observed trajectory.     

Acoustic-Gravity Waves from Bolide Sources

We have developed a new approach to modeling the acoustic-gravity wave (AGW) radiation from bolide sources. This first effort involves entry modeling of bolide sources that have available satellite data through procedures developed in ReVelle (Earth Moon Planets 95, 441–476, 2004a; in: A. Milani, G. Valsecchi, D. Vokrouhlicky (eds) NEO Fireball Diversity: Energetics-based Entry Modeling and Analysis Techniques, Near-earth Objects: Our Celestial Neighbors (IAU S236), 2007b). Results from the entry modeling are directly coupled to AGW production through line source blast wave theory for the initial wave amplitude and period at (at 10 blast wave radii and perpendicular to the trajectory). The second effort involves the prediction of the formation and or dominance of the propagation of the atmospheric Lamb, edge-wave composite mode in a viscous fluid (Pierce, J. Acoust. Soc. Amer. 35, 1798–1807, 1963) as a function of the source energy, horizontal range and source altitude using the Lamb wave frequency that was deduced directly during the entry modeling and that is used as a surrogate for the source energy. We have also determined that Lamb wave production by bolides at close range decreases dramatically as either the source energy decreases or the source altitude increases. Finally using procedures in Gill (Atmospheric-Ocean Dynamics, 1982) and in Tolstoy (Wave Propagation, 1973), we have analyzed two simple dispersion relationships and have calculated the expected dispersion for the Lamb edge-wave mode and for the excited, propagating internal acoustic waves. Finally, we have used the above formalism to fully evaluate these techniques for four large bolides, namely: the Tunguska bolide of June 30, 1908; the Revelstoke bolide of March 31, 1965; the Crete bolide of June 6, 2002 and the Antarctic bolide of September 3, 2004. Due to page limitations, we will only present results in detail for the Revelstoke bolide.     

Perseid meteor stream 1991 – 1993 from TV observations in Kiev

Some results of the double station television meteor observations provided in Kiev during the Perseid shower period in years 1991 – 1993 are presented. The dependence of the beginning hights of meteors on their initial velocities based on 57 best reduced double station meteors is constructed. Exceptionaly great beginning heights of meteors are discussed, as well as the exceptional case of a Perseid bolide.     

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.     

The Pribram,Lost City,Innisfree, and Neuschwanstein falls: An analysis of the atmospheric trajectories

To date, several meteorites have been found for which their flight in the atmosphere was recorded by special fireball camera networks. Because of this, a thorough analysis of the instrumentally registered falls is of current importance. For such fireballs, not only the high-quality photo images of the motion in the atmosphere exist, but also the density and the shape of the meteor body fragments reached the Earth’s surface are known for sure. In the present study, for the Innisfree, Lost City, and Pribram fireballs, new models of the entry to the atmosphere have been built. The values of the ballistic coefficient and the mass-loss parameter providing the best approximation for the observations of the luminous trajectory segment with the analytical solution of the meteor physics equations have been obtained. From recent results of the numerical experiments on the supersonic airflow of bodies of various shapes, the preatmospheric masses of the fireballs, as well as the dynamic estimates of the mass at the other trajectory points, were obtained. In particular, the terminal mass of the fireballs in the lower segment of the analyzed trajectories is in good agreement with the total mass of the meteorite material recovered in all of the cases considered. Moreover, to calculate the acceleration of the meteor bodies, a new analytical formula has been suggested, which allows the obtained theoretical time dependencies of the velocity and altitude to be compared with the observational data.     

The extra-atmospheric mass of small meteoroids of the Prairie and Canada bolide camera networks

The existing methods for determining the extra-atmospheric mass of meteor bodies from observations of their movement in the atmosphere allow a certain arbitrariness. Active attempts to overcome the discrepancy between the results of calculations based on different approaches often lead to physically incorrect conclusions. A way out is to laboriously accumulate the estimates and computation results and to consistently remove ambiguities. To correctly interpret the observed brightness of a meteor, one should use contemporary methods and the results of physical studies of the emitting gas. In the present work, the extra-atmospheric masses of small meteoroids of the Prairie and Canada bolide camera networks were calculated from the observed braking. It turned out that, in many cases, the conditions of movement of meteor bodies in the atmosphere corresponded to a free molecular airflow about a body. The so-called dynamic mass of the bodies was estimated from the real densities of the meteoroid material, which corresponded to monolithic water ice and stone, and for the proper values of the product of the drag coefficient and shape factor. When producing the trial function for the body trajectories in the “velocity-altitude” variables, we did not allow for fragmentation explicitly, since it is less probable for small meteoroids than for large ones. As before, our estimates differ substantially from the photometric masses published in the corresponding tables.     

Entry dynamics and acoustics/infrasonic/seismic analysis for the Neuschwanstein meteorite fall

Abstract— We have analyzed several types of data associated with the well‐documented fall of the Neuschwanstein meteorites on April 6, 2002 (a total of three meteorites have been recovered). This includes ground‐based photographic and radiometer data as well as infrasound and seismic data from this very significant bolide event (Spurný et al. 2002, 2003). We have also used these data to model the entry of Neuschwanstein, including the expected dynamics, energetics, panchromatic luminosity, and associated fragmentation effects. In addition, we have calculated the differential efficiency of acoustical waves for Neuschwanstein and used these values to compare against the efficiency calculated using available ground‐based infrasound data. This new numerical technique has allowed the source height to be determined independent of ray tracing solutions. We have also carried out theoretical ray tracing for a moving point source (not strictly a cylindrical line emission) and for an infinite speed line source. In addition, we have determined the ray turning heights as a function of the source height for both initially upward and downward propagating rays, independent of the explicit ray tracing (detailed propagation path) programs. These results all agree on the origins of the acoustic emission and explicit source heights for Neuschwanstein for the strongest infrasonic signals. Calculated source energies using more than four different independent approaches agree that Neuschwanstein was certainly <500 kg in initial mass, given the initial velocity of 20.95 km/s, resulting in an initial source energy ≤0.0157‐0.0276 kt TNT equivalent (4.185 times 10¹² J). Local source energies at the calculated infrasonic/seismic source altitudes are up to two orders of magnitude smaller than this initial source energy.     

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.     

Multi-Instrument Observations of Bright Meteors in the Czech Republic

We present detailed data on 8 bright meteors recorded simultaneously by different observational techniques. All meteors were recorded by all-sky cameras at the Czech stations of the European Fireball Network and by image intensified TV cameras placed at Ondrejov and Kunzak observatories. As well as direct photographic and LLLTV recordings, most of meteors were recorded also by the spectral TV camera and some also by photographic spectral cameras. For 6 cases, lightcurves from radiometers with very high time resolution (1200 s⁻¹) are also available. From all these detections we found a significant difference between TV and photographic beginning heights. TV beginnings are in average about 40 km higher than the photographic ones. We found that meteor brightness is up to 2 magnitudes higher in the photographic system than in the TV system. This difference for high velocity meteors is mainly caused by the presence of strong Ca⁺ lines in the blue part of the spectrum, where the image intensifier is only marginally sensitive. At heights above 110 km, the Na line is usually brighter than the Mg line, while at lower heights both lines have comparable brightness. In one of two captured spectra of short duration luminous trains, a small initial brightening of the Mg and Na lines caused by recombination processes was observed.     

Genesis—An artificial,low velocity “meteor” fall and recovery: September 8, 2004

Abstract— On September 8, 2004, Genesis, a manmade space capsule, plummeted to Earth after almost three years in space. A ground‐based infrasound array was deployed to Wendover, Nevada, to measure the “hypersonic boom” from the reentry, since the expected initial reentry speed of the body was about 11 km/sec. Due to the complete failure of its dual parachute system, we had a unique opportunity to assess the degree of reliability of our previously developed relations for natural meteors and bolides to analyze this well‐characterized manmade body. At ?20–50 km from the nominal trajectory, we succeeded in recording over two minutes of infrasonic signals from Genesis. Here we report on subsequent analyses of these infrasonic data, including an assessment of the expected entry characteristics on the basis of a bolide/meteor/fireball entry model specifically adapted to modeling reentering manmade objects. From these simulations, we were able to evaluate the line source blast wave relaxation radius, the differential acoustic efficiency, etc., to compute an approximate total power balance during entry. Next, we analyzed the detailed signals arriving from Genesis using a numerical, signal detection and wave processing software package (Matseis/Infra_Tool). We established the initial and subsequent arrivals and evaluated its plane wave back azimuths and elevation arrival angles and the degree of maximum, pair‐wise cross‐correlation, its power spectrum, spectrogram analysis, standard seismic f‐k analysis, etc. From the associated entry parameters, we computed the kinetic energy density conservation properties for the propagating line source blast waves and compared these predictions against observed ground‐based infrasound amplitude and wave period data as a function of range. We discovered that previously computed differential acoustic efficiencies were unreliable at Mach numbers below about 10. This is because we had assumed that a line source explosion was applicable, whereas at very low Mach numbers, typical of recovered meteorites, the detailed source characteristics are closer to those of supersonic objects. When corrections for these unphysical, very high efficiencies were made, agreement between theory and observations improved. We also made an assessment for the energy of the blast wave source from the ground‐based infrasound data using several other techniques that were also adapted from previous bolide studies. Finally, we made a top‐down‐bottom‐up assessment of the line source wave normals propagating via refraction downward into the complex middle atmospheric environment. This assessment proved to be generally consistent with the digital signal processing analysis and with the observed time delay between the known Genesis reentry and the infrasonic observations.