Trajectory and orbit of the Tunguska meteorite revisited

Abstract— A critical survey is presented of all determinations of the azimuth and inclination of the Tunguska meteorite's trajectory based either on eyewitness testimonies or on the mathematical treatment of the forest-leveling field in the area of the catastrophe. The eyewitness testimonies collected in the neighborhood of the Nizhnyaya Tunguska River indicate the most probable azimuth of the trajectory projection to be 104° from the north to the east, which is close to the most recent azimuth estimate from the forest-leveling field, 99°. For the most part of the trajectory, its inclination could not exceed 15°. However, it is seen from aerodynamic calculations that the combined action of the gravity field and a nonzero aerodynamic lift could increase the inclination to 40° as the end of the trajectory was approached. Meteoroid orbits are calculated for a broad family of trajectories with azimuths ranging from 99° (Fast et al, 1976) to 137° (Krinov, 1949) and geocentric velocities ranging from 25 to 40 km/s. Orbits with large azimuth values (120° and larger) are shown to belong to the asteroidal type. They are succeeded by the orbits of short-period and long-period comets, whereas very small azimuth values and large geocentric velocities correspond to the region of hyperbolic orbits. Certain restrictions on the possible trajectory azimuths and geocentric velocities of the Tunguska body are imposed by this study.     

THE TUNGUSKA EVENT AND THE METEORS WITH TERMINAL FLARES

The comparison of the Tunguska body explosion with the effect of terminal flares of meteors and fireballs leads us to the conclusion that these events are of a similar nature but differ only by their scale. We consider that the dynamics of progressive breaking and evaporation of meteoric bodies during their entry into the terrestrial atmosphere could explain the terminal burst. An extremely porous body model for the Tunguska meteorite was analysed and rejected as unsatisfactory. The realistic values of the initial velocity (～30 km/sec) and of the inclination angle for the Tunguska's trajectory (5–15°) give orbital elements not in contradiction with the cometary origin of the Tunguska body.     

Comment on “Evidence for a very high carbon/iridium ratio in the Tunguska impactor” by K. L. Rasmussen,H. J. F. Olsen,R. Gwozdz and E. M. Kolesnikov

Abstract— We discuss possible evidence for a dilution of ¹⁴C caused by the Tunguska impact event, proposed by Rasmussen et al. (1999). The results presented in that paper and other available information do not support this hypothesis.     

MASS FLUX OF ASTEROIDAL ORIGIN METEOROIDS ON PERIODIC COMET NUCLEI

We examine the potential contamination of cometary nuclei through impacts from asteroidal origin meteoroids. The paper uses a simple model and has the goal of determining whether asteroidal contamination is potentially significant. We assume a meteoroid power law mass distribution with index values in the range from s=1.83 to s=2.09. We used maximum and minimum models which we believe will bracket the true meteoroid mass distribution. We identify those comets which are expected to be most significantly contaminated, and find values of up to 3.6 kg of asteroidal meteoroid impact per square meter of the cometary surface per orbital revolution. This is less than the expected mass loss per perihelion passage for most comets. Therefore any remnant effects of the contamination will depend on the penetration depth of the meteoroids in the cometary nucleus, and possibly on the distribution of active and inactive areas on cometary nuclei. We present a simple model which suggests that even small meteoroids will embed relatively deeply into a cometary nucleus.     

Determination of the height of the “meteoric explosion”

When cosmic bodies of asteroidal and cometary origin, with a size from 20 to approximately 100 m, enter dense atmospheric layers, they are destroyed with a large probability under the action of aerodynamic forces and decelerated with the transfer of their energy to the air at heights from 20–30 to several kilometers. The forming shock wave reaches the Earth’s surface and can cause considerable damage at great distances from the entry path similar to the action of a high-altitude explosion. We have performed a numerical simulation of the disruption (with allowance for evaporation of fragments) and deceleration of meteoroids having the aforesaid dimensions and entering the Earth’s atmosphere at different angles and determined the height of the equivalent explosion point generating the same shock wave as the fall of a cosmic body with the given parameters. It turns out that this height does not depend on the velocity of the body and is approximately equal to the height at which this velocity is reduced by half. The obtained results were successfully approximated by a simple analytical formula allowing one to easily determine the height of an equivalent explosion depending on the dimensions of the body, its density, and angle of entry into the atmosphere. A comparison of the obtained results with well-known approximate analytical (pancake) models is presented and an application of the obtained formula to specific events, in particular, to the fall of the Chelyabinsk meteorite on February 15, 2013, and Tunguska event of 1908, is discussed.     

Authors' Reply

Abstract— Jull et al. propose an alternative interpretation of our depth vs. ¹⁴C data measured on a peat core from the central Tunguska impact site (Rasmussen et al., 1999). We find that the proposed alternative is untenable.     

Nitrate in the Greenland ice sheet in the years following the 1908 tunguska event

The Tunguska event on 30 June 1908 has been subjected to much speculation within different fields of research. Publication of the results of the 1961 expedition to the Tunguska area (Florensky, 1963) supports that a cometary impact caused the event. Based on this interpretation, calculations of the impactor energy release and explosion height have been reported by Ben-Menahem (1975), and velocity, mass, and density of the impactor by Petrov and Stulov (1975). Park (1978) and Turco et al., 1981, Turco et al., 1982, used these numbers to calculate a production of ca. 30 × 10⁶ tons of NO during atmospheric transit. This paper presents a high-resolution study of nitrate concentration in the Greenland ice sheet in ca. 10 years covering the Tunguska event. No signs of excess nitrate are found in three ice cores from two different sites in Greenland in the years following the Tunguska event. By comparing these results with results for other aerosols generally found in the ice, the lack of excess NO₃^? following the Tunguska event can be interpreted as indicating that the impactor nitrate production calculated by Park (1978) and Turco et al., 1981, Turco et al., 1982 are 1–2 orders of magnitude too high. To explain this it is suggested, from other lines of reasoning, that the impactor density determined by Petrov and Stulov (1975) probably is too low.     

Migration of Comets

Based on one of the particular cases of twice averaged model Hill problem with the allowance for the oblateness of the central body a quadrature is derived for the determination of the migration time of cometary nuclei from various cometary reservoirs and a (14-th order) algebraic equation for the determination of the initial conditions that allow the escape of the cometary nucleus (which at the initial instant of time moves in an orbit with arbitrary eccentricity (0 < e < 1) and inclination (0° < i < 180°) deeply inside the sphere o f action of the central body) from the sphere of action of the central body or its impact onto the central body. We analyze the shape of the boundaries of the hypothetical cometary reservoirs and the method of searching for regions of high concentration of interstellar particles in the Solar System.     

Evidence for a very high carbon/iridium ratio in the Tunguska impactor

Abstract— We have measured excess Ir and depletion of ¹⁴C, two independent indicators of cosmic material, in peat cores from the central Tunguska impact site. Both Ir and ¹⁴C show pronounced anomalies in the same stratigraphical depth interval. We have estimated an integral deposition of nonradioactive cosmogenic C of 6.8 ± 1.0 mg C cm^?2, and an integrated Ir deposition of 5.9 ± 1.2 pg Ir cm^?2. The very high C/Ir ratio and a deduced δ¹³C value of +55 ± 10% relative to V Pee Dee Belemnite (VPDB) of the impactor material found in this study points towards a cometary type impactor, rather than a chondritic or achondritic asteroidal type impactor.     

The thermal history of interplanetary dust particles collected in the Earth's stratosphere

Abstract— Fragments of 24 individual interplanetary dust particles (IDPs) collected in the Earth's stratosphere were obtained from NASA's Johnson Space Center collection and subjected to pulse-heating sequences to extract He and Ne and to learn about the thermal history of the particles. A motivation for the investigation was to see if the procedure would help distinguish between IDPs of asteroidal and cometary origin. The use of a sequence of short-duration heat pulses to perform the extractions is an improvement over the employment of a step-heating sequence, as was used in a previous investigation. The particles studied were fragments of larger parent IDPs, other fragments of which, in coordinated experiments, are undergoing studies of elemental and mineralogical composition in other laboratories. While the present investigation will provide useful temperature history data for the particles, the relatively large size of the parent IDPs (～40 μm in diameter) resulted in high entry deceleration temperatures. This limited the usefulness of the study for distinguishing between particles of asteroidal and cometary origin.     

Evidence from IRAS for a very young, partially formed dust band

The relative proportions of asteroidal and cometary materials in the zodiacal cloud is an ongoing debate. The determination of the asteroidal component is constrained through the study of the Solar System dust bands (the fine-structure component superimposed on the broad background cloud), since they have been confidently linked to specific, young, asteroid families in the main belt. The disruptions that produce these families also result in the injection of dust into the cloud and thus hold the key to determining at least a minimum value for the asteroidal contribution to the zodiacal cloud. There are currently known to be at least three dust band pairs, one at approximately 9.35° associated with the Veritas family and two central band pairs near the ecliptic, one of which is associated with the Karin subcluster of the Koronis family. Through careful co-adding of almost all the pole-to-pole intensity scans in the mid-infrared wavebands of the Infrared Astronomical Satellite (IRAS) data set, we find strong evidence for a partial Solar System dust band, that is, a very young dust band in the process of formation, at approximately 17° latitude. We think this is a confirmation of the M/N partial band pair first suggested by Sykes [1988. IRAS observations of extended zodiacal structures. Astrophys. J. 334, L55-L58]. The new dust band appears at some but not all ecliptic longitudes, as expected for a young, partially formed dust band. We present preliminary modeling of the new, partial dust band which allows us to put constraints on the age of the disruption event, the inclination and node of the parent body at the time of disruption, and the quantity of dust injected into the zodiacal cloud.     

Explosion of a comet nucleus fragment in the earth’s atmosphere

This paper analyzes data on thermal explosions of large meteoroids in the earth’s atmosphere. The cumulative function of flux of space bodies is corrected with regard to the explosion height, which is determined, according to our approach, by maximum braking. As a result, the integral function of flux in the work [Brown, P., Spalding, R.E., ReVelle, D.O., et al., The Flux of Small Near-Earth Objects Colliding with the Earth, Nature, 2002, vol. 420, pp. 314–316] is consistent with the one we derived earlier. It is found that at least one phenomenon of those discussed in the paper by Brown et al. is a result of explosion of a comet nucleus fragment. It is shown that the Tunguska phenomenon cannot be explained within a monolithic body model.     

On the magnetic field in the tail of Comet Halley

The method of evaluation of the cometary magnetic field proposed by Podgornyet al. (1980) is shown not to be self-consistent. An alternative method is discussed.     

A clast of Bali‐like oxidized CV material in the reduced CV chondrite breccia Vigarano

Abstract— The CV (Vigarano‐type) chondrites are a petrologically diverse group of meteorites that are divided into the reduced and the Bali‐like and Allende‐like oxidized subgroups largely based on secondary mineralogy (Weisberg et al., 1997; Krot et al., 1998b). Some chondrules and calcium‐aluminum‐rich inclusions (CAIs) in the reduced CV chondrite Vigarano show alteration features similar to those in Allende: metal is oxidized to magnetite; low‐Ca pyroxene, forsterite, and magnetite are rimmed and veined by ferrous olivine (Fs_40–50); and plagioclase mesostases and melilite are replaced by nepheline and sodalite (Sylvester et al., 1993; Kimura and Ikeda, 1996, 1997, 1998). Our petrographic observations indicate that Vigarano also contains individual chondrules, chondrule fragments, and lithic clasts of the Bali‐like oxidized CV materials. The largest lithic clast (about 1 times 2 cm in size) is composed of opaque matrix, type‐I chondrules (400–2000 μm in apparent diameter) surrounded by coarse‐grained and fine‐grained rims, and rare CAIs. The matrix‐chondrule ratio is about 1.1. Opaque nodules in chondrules in the clast consist of Cr‐poor and Cr‐rich magnetite, Ni‐ and Co‐rich metal, Ni‐poor and Ni‐rich sulfide; low‐Ni metal nodules occur only inside chondrule phenocrysts. Chromium‐poor magnetite is preferentially replaced by fayalite. Chondrule mesostases are replaced by phyllosilicates; low‐Ca pyroxene and olivine phenocrysts appear to be unaltered. Matrix in the clast consists of very fine‐grained (<1 μm) ferrous olivine, anhedral fayalite grains (Fa_80–100), rounded objects of porous Ca‐Fe‐rich pyroxenes (Fs_10–50Wo₅₀), Ni‐poor sulfide, Ni‐ and Co‐rich metal, and phyllosilicates; magnetite is rare. On the basis of the presence of the Bali‐like lithified chondritic clast—in addition to individual chondrules and CAIs of both Bali‐like and Allende‐like materials—in the reduced CV chondrite Vigarano, we infer that (1) all three types of materials were mixed during regolith gardening on the CV asteroidal body, and (2) the reduced and oxidized CV materials may have originated from a single, heterogeneously altered asteroid.     

On the formation of cometary nuclei in dense globules

Evolution of cometary orbits by planetary perturbations, weakly hyperbolic original orbits of comets calculated by Marsdenet al. (1978) are taken to indicate the interstellar origin of comets, and the possible formation of cometary nuclei in interstellar globules is discussed. The process is sedimentation of dust grains. It is shown that if a globule is at 40 K, its lifetime is sufficiently long to allow the sedimentation.     

Gravitational effects after the impact disruption of a minor planet: geometrical properties and criteria for the reaccumulation

After the catastrophic disruption of a planetary body the fragments move according to their mutual gravitational attraction, finally resulting into a more or less massive reaccumulation, as well as into the formation of binary or multiple systems. In this paper we analyze this process by means of the outcomes of the semi-empirical model of catastrophic impacts described by Paolicchi et al., 1996and applied to impacts at planetary sizes by Paolicchi et al., 1993and Doressoundiram et al., 1997. It is possible to identify the location in the parent body of fragments which are going to be reaccumulated or ejected, or to form binaries.Moreover, we compare the results of numerical integrations with three analytical predicting criteria existing in the literature, as well as with a new one, based on the definition of iso-velocity surfaces. We show that: (a) two of the criteria presented and used in the literature may lead to severe over- or underestimates of the amount of reaccumulated mass; (b) the new criterion introduced here and the previous one described by Petit and Farinella, 1993are capable of giving an accurate estimate in many cases, but are less effective when the reaccumulation is limited and not strongly concentrated onto a single big attractor; (c) the region in the target where the analytic criteria fail is generally the transition region between escaping and reaccumulating. This is the same region where other interesting phenomena, such as the formation of binaries, take place.     

Possible effects of secular resonances in Phobos and Triton

Due to the tides, the orbits of Phobos and Triton are contracting. While their semi major axes are decreasing, several possibilities of secular resonances involving node, argument of the pericenter and mean motion of the Sun will take place. In the case of Mars, if the obliquity (ε), during the passage through some resonances, is not so small, very significant variations of the inclination will appear. In one case, capture is almost certain provided that ε?20°. For Triton there are also similar situations, but capture seems to be not possible, mainly because in S₁ state, Triton's orbit is sufficiently inclined (far) with respect to the Neptune's equator. Following Chyba et al. (Astron. Astrophys. 219 (1989) 123), a simplified equation that gives the evolution of the inclination versus the semi major axis, is derived. The time needed for Triton crash onto Neptune is longer than that one obtained by these authors, but the main difference is due to the new data used here. In general, even in the case of non-capture passages, some significant jumps in inclination and in eccentricities are possible.     

Quantitative Analysis of H2OComa Images Using a Multiscale MHD Model with Detailed Ion Chemistry

The observation of ions created by ionization of cometary gas, either by ground-based observations or byin situmeasurements can give us useful information about the gas production and composition of comets. However, due to the interaction of ions with the magnetized solar wind and their high chemical reactivity, it is not possible to relate measured ion densities (or column densities) directly to the parent gas densities. In order to quantitatively analyze measured ion abundances in cometary comae it is necessary to understand their dynamics and chemistry. We have developed a detailed ion–chemical network of cometary atmospheres. We include production of ions by photo- and electron impact-ionization of a background neutral atmosphere, charge exchange of solar wind ions with cometary atoms/molecules, reactions between ions and molecules, and dissociative recombination of molecular ions with thermal electrons. By combining the ion–chemical network with the three-dimensional plasma flow as computed by a new fully three-dimensional MHD model of cometary plasma environments (Gombosiet al.1996) we are able to compute the density of the major cometary ions everywhere in the coma. The input parameters for our model are the solar wind conditions (density, speed, temperature, magnetic field) and the composition and production rate of the gas. We applied our model to Comet P/Halley in early March 1986, for which the input parameters are reasonably well known. We compare the resulting column density of H₂O⁺with ground-based observations of H₂O⁺from DiSantiet al.(1990). The results of our model are in good agreement with both the spatial distribution and the absolute abundance of H₂O⁺and with their variations with the changing overall water production rate between two days. The results are encouraging that it will be possible to obtain production rates of neutral cometary constituents from observations of their ion products.     

Impact origin of the Vesta family

Abstract— The compelling petrographic link (Consolmagno and Drake, 1977; Gaffey, 1983) between basaltic achondrite meteorites and the ～530 km diameter asteroid 4 Vesta has been tempered by a perceived difficulty in launching rocks from this asteroid's surface at speeds sufficient to bring them to Earth (Wasson and Wetherill, 1979) without obliterating Vesta's signature crust. I address this impasse in response to recent imaging (Zellner et al, 1996; Dumas and Hainaut, 1996) of a ～450 km impact basin across Vesta's southern hemisphere (Thomas et al., 1997) and model the basin-forming collision using a detailed two-dimensional hydrocode with brittle fracture including self-gravitational compression (cf., Asphaug and Melosh, 1993). A ～42 km diameter asteroid striking Vesta's basaltic crust (atop a denser mantle and iron core) at 5.4 km/s launches multikilometer fragments up to ～600 m/s without inverting distal stratigraphy, according to the code. This modeling, together with collisional, dynamical, rheological and exposure-age timescales (Marzari et al., 1996; Welten et al., 1996), and observations of V-type asteroids (Binzel and Xu, 1993) suggests a recent (<～1 Ga) impact origin for the Vesta family and a possible Vesta origin for Earth-approaching V-type asteroids (Cruik-shank et al., 1991).     

Stability of protoplanets around RU Lupi

It is shown that the protoplanetary dust condensations around the T-Tauri star RU Lup, proposed by Gahmet al. (1974), would be unstable to tidal disruption if they are bound only by gravitational forces. There are extreme conditions under which the condensations would be stable, but such conditions would easily be verifiable observationally.