Atmospheric shock waves after impacts of cosmic bodies up to 1000 m in diameter

The results of numerical modeling of meteoroids' interaction with Earth's atmosphere are presented. We model the entry in two dimensions and then interpolate the results into a 3‐D model to calculate interaction of shock waves with the surface. Maximum shock pressures, wind speeds, and areas subjected to substantial overpressure are calculated for oblique impacts of asteroids and comets. We show that vertical impacts produce a smaller damage zone on the surface than oblique ones. Damage caused by shock waves covers an order of magnitude larger area than any other hazardous effects. The function of energy release in the atmosphere, which is traditionally used in meteoritics, has a limited application if cosmic bodies are larger than tens of meters in diameter: at each time moment energy is smoothed along a substantial length of the trajectory; both emitted radiation (routinely used for calibration of semi‐analytical models) and shock wave amplitude are complex functions of temperature–density distributions in atmosphere.     

SUPRACENTER: Locating fireball terminal bursts in the atmosphere using seismic arrivals

Abstract— Terminal bursts and fragmentations of meteoritic fireballs in the atmosphere may now be accurately located in four dimensions (three spatial + temporal) using seismic arrival times of their acoustic waves recorded by seismometer, camera, microphone, and/or infrasound stations on the ground. A computer program, SUPRACENTER, calculates travel times by ray tracing through realistic atmospheres (that include winds) and locates source positions by minimization of travel time residuals. This is analogous to earthquake hypocenter location in the solid Earth but is done through a variably moving medium. Inclusion of realistic atmospheric ray tracing has removed the need for the simplifying assumption of an isotropic atmosphere or an approximation to account for “wind drift.” This “drift” is on the order of several km when strong, unidirectional winds are present in the atmosphere at the time of a fireball's occurrence. SUPRACENTER‐derived locations of three seismically recorded fireballs: 1) the October 9, 1997 El Paso superbolide; 2) the January 25, 1989 Mt. Adams fireball; and 3) the May 6, 2000 Morávka fireball (with its associated meteorite fall), are consistent with (and, probably, an improvement upon) the locations derived from eyewitness, photographic, and video observations from the respective individual events. If direct acoustic seismic arrivals can be quickly identified for a fireball event, terminal burst locations (and, potentially, trajectory geometry and velocity information) can be quickly derived, aiding any meteorite recovery efforts during the early days after the fall. Potentially, seismic records may yield enough trajectory information to assist in the derivation of orbits for entering projectiles.     

A new method of meteor trajectory determination applied to multiple unsynchronized video cameras

Abstract– A new approach has been formulated to solve for the straight line trajectory of a meteor through the atmosphere when given multiple camera views of the meteor’s luminous track. Using a motion propagation model in 3‐D space plus time, and iteratively solving for all free model parameters simultaneously, one can obtain a fully coupled solution to the apparent radiant direction, 3‐D begin position, atmospheric entry speed, deceleration terms, and timing offsets when using data from unsynchronized video cameras. A Monte Carlo component adds empirical error estimation for each of the key model parameters computed. This multiparameter fitting method extends the allowable collection geometries for meteor trajectory estimation to lower convergence angles between camera‐meteor‐camera lines of sight and smaller site separation distances.     

“Comet‐tail” ejecta streaks: A predicted cratering landform unique to Titan

Abstract— A model for an impact ejecta landform peculiar to Saturn's moon Titan is presented. Expansion of the ejecta plume from moderate‐sized craters is constrained by Titan's thick atmosphere. Much of the plume is collimated along the incoming bolide's trajectory, as was observed for plumes from impacts on Jupiter of P/Shoemaker‐Levy‐9, but is retained as a linear, diagonal ejecta cloud, unlike on Venus where the plume “blows out.” On Titan, the blowout is suppressed because the vertically‐extended atmosphere requires a long wake to reach the vacuum of space, and the modest impact velocities mean plume expansion along the wake is slow enough to allow the wake to close off. Beyond the immediate ejecta blanket around the crater, distal ejecta is released into the atmosphere from an oblique line source: this material is winnowed by the zonal wind field to form streaks, with coarse radar‐bright particles transported less far than fine radar‐dark material. Thus, the ejecta form two distinct streaks faintly reminiscent of dual comet tails, a sharply W‐E radar‐dark one, and a less swept and sometimes comma‐shaped radar‐bright one.     

Orbit and origin of the LL7 chondrite Dishchii'bikoh (Arizona)

The trajectory and orbit of the LL7 ordinary chondrite Dishchii'bikoh are derived from low‐light video observations of a fireball first detected at 10:56:26 UTC on June 2, 2016. Results show a relatively steep ~21° inclined orbit and a short 1.13 AU semimajor axis. Following entry in Earth's atmosphere, the meteor luminosity oscillated corresponding to a meteoroid spin rate of 2.28 ± 0.02 rotations per second. A large fragment broke off at 44 km altitude. Further down, mass was lost to dust during flares at altitudes of 34, 29, and 25 km. Surviving meteorites were detected by Doppler weather radar and several small 0.9–29 g meteorites were recovered under the radar reflection footprint. Based on cosmogenic radionuclides and ground‐based radiometric observations, the Dishchii'bikoh meteoroid was 80 ± 20 cm in diameter assuming the density was 3.5 g/cm³. The meteoroid's collisional history confirms that the unusual petrologic class of LL7 does not require a different parent body than three previously observed LL chondrite falls. Dishchii'bikoh was ejected 11 Ma ago from parent body material that has a 4471 ± 6 Ma U‐Pb age, the same as that of Chelyabinsk (4452 ± 21 Ma). The distribution of the four known pre‐impact LL chondrite orbits is best matched by dynamical modeling if the source of LL chondrites is in the inner asteroid belt in a low inclined orbit, with the highly inclined Dishchii'bikoh being the result of interactions with Earth before impacting.     

A model for Type-IV emission in the solar burst of June 9, 1959 at decametric wavelengths

A major radio burst at decametric frequencies at 1638 UT on June 9, 1959 is apparently a Type-IV continuum burst of the kind that drifts from high to low frequencies. We present observations of flux variations and East-West positions of the emission at both 18 and 38 MHz. The burst moves outward at a speed of about 4700 km·sec^?1 at each frequency to a height of about 3 R _⊙ from the sun's center and then returns to the sun. This behavior is not simultaneous at 18 MHz and 38 MHz; the outward moving phase of 18 MHz emission occurs during the return phase of 38 MHz. We suspect that a solitary Alfvén wave or shock traverses the outer corona at the time of this burst. Relativistic electrons created low in the solar atmosphere travel freely along radial lines of force up to the coronal-streamer heights of the Type-IV burst. Upon encountering the shock, the electrons emit locally intense synchrotron emission, and pass through the shock on out into interplanetary space. This model appears to be consistent with other shock front phenomena in interplanetary space and the corona. Finally, the Razin effect (Boischot and Clavelier, 1967) suggests that low-frequency cut-offs in Type-IV bursts ought to be quite constant in frequency, and not higher than between 0.4 to 4 MHz.     

On penetration depth of the Shoemaker-Levy 9 fragments into the Jovian atmosphere

The actual penetration depth of the Shoemaker-Levy 9 fragments into the Jovian atmosphere is still an open question. From fundamental equations of meteoric physics with variable cross-section, a new analytic model of energy release of the fragments is presented. In use of reasonable parameters, a series of results are calculated for different initial mass of the fragments. The results show that the largest fragment explodes above pressure levels of 3 bars and does not penetrate into the H₂O cloud layer of the Jovian atmosphere, and that airburst of smaller fragments occur even above the upper cloud layer.     

An atmospheric blast/thermal model for the formation of high‐latitude pedestal craters

Abstract— Although tenuous, the atmosphere of Mars affects the evolution of impact‐generated vapor. Early‐time vapor from a vertical impact expands symmetrically, directly transferring a small percentage of the initial kinetic energy of impact to the atmosphere. This energy, in turn, induces a hemispherical shock wave that propagates outward as an intense airblast (due to high‐speed expansion of vapor) followed by a thermal pulse of extreme atmospheric temperatures (from thermal energy of expansion). This study models the atmospheric response to such early‐time energy coupling using the CTH hydrocode written at Sandia National Laboratories. Results show that the surface surrounding a 10 km diameter crater (6 km “apparent” diameter) on Mars will be subjected to intense winds (?200 m/s) and extreme atmospheric temperatures. These elevated temperatures are sufficient to melt subsurface volatiles at a depth of several centimeters for an ice‐rich substrate. Ensuing surface signatures extend to distal locations (?4 apparent crater diameters for a case of 0.1% energy coupling) and include striations, thermally armored surfaces, and/or ejecta pedestals—all of which are exhibited surrounding the freshest high‐latitude craters on Mars. The combined effects of the atmospheric blast and thermal pulse, resulting in the generation of a crater‐centered erosion‐resistant armored surface, thus provide a new, very plausible formation model for high‐latitude Martian pedestal craters.     

Theoretical study of partial frequency redistribution function in irradiated, moving atmospheres of close binary components

We have studied the effects of partial frequency redistribution function with angle-averaged   R _II-A  in irradiated and moving atmospheres of close binary components. We have considered the atmospheric extension of the primary component to be twice the radius of the primary component in a close binary system. We have considered two cases: (i) when the atmosphere is at rest and (ii) when the atmosphere is moving. In both the cases, we have computed the line profiles along the line of sight for a given optical depth. The irradiation from the secondary component is assumed to be one, five and 10 times the self-radiation. The line fluxes in the line of sight are calculated by using the total source functions due to self-radiation of the primary component and due to the irradiation from the secondary component. We have noted double-peaked emission lines in the case of a static medium and a reduction of emission peaks in the case of velocity field.     

一个微波爆发的源区磁场和高能电子的能流估计

本文分析了１９９３年１０月２日０７：３９：４０－０７：４１：００ＵＴ时段太阳产生的一个多脉冲微波暴的观测，认为它是由多个脉冲爆发叠加在一个慢变爆发背景上组成的．根据谱分析和利用我们的日冕磁场诊断公式［１］，第一次获得了一个爆发源区的磁场强度和高能电子的信息，其主要结果是：（１）脉冲爆发分量在光薄部分的射电谱指数的平均值比慢变爆发背景的值小１，即前者的谱比后者的硬．在１９．６ＧＨｚ上的亮度温度前者比后者高６倍．（２）从脉冲爆发分量和慢变爆发背景分量推断的源区磁场平均值分别为１５８和５３１高斯，且发现在爆发期间，慢变暴源区磁场强度随时间圣马鞍形变化，在极大相的值比脉冲相和下降相低约５０％（３）产生脉冲暴分量的高能电子的柱密度ＮＬ和数密度Ｎ（＞Ｅ０）分别为慢变暴分量的４％和８％，但它们所携带的能流和发射系数要比慢变爆发分量的值高１倍和８倍！表明这两种爆发成份可能分别来自能谱不同的两群电子在不同爆发源区的辐射．     

The rate of small impacts on Earth

Abstract— Asteroids tens to hundreds of meters in diameter constitute the most immediate impact hazard to human populations, yet the rate at which they arrive at Earth's surface is poorly known. Astronomic observations are still incomplete in this size range; impactors are subjected to disruption in Earth's atmosphere, and unlike the Moon, small craters on Earth are rapidly eroded. In this paper, we first model the atmospheric behavior of iron and stony bodies over the mass range 1–10¹² kg (size range 6 cm‐1 km) taking into account deceleration, ablation, and fragmentation. Previous models in meteoritics deal with rather small masses (<10⁵–10⁶ kg) with the aim of interpreting registered fireballs in atmosphere, or with substantially larger objects without taking into account asteroid disruption to model cratering processes. A few earlier attempts to model terrestrial crater strewn fields did not take into account possible cascade fragmentation. We have performed large numbers of simulations in a wide mass range, using both the earlier “pancake” models and also the separated fragments model to develop a statistical picture of atmosphere‐bolide interaction for both iron and stony impactors with initial diameters up to ?1 km. Second, using a compilation of data for the flux at the upper atmosphere, we have derived a cumulative size‐frequency distribution (SFD) for upper atmosphere impactors. This curve is a close fit to virtually all of the upper atmosphere data over 16 orders of magnitude. Third, we have applied our model results to scale the upper atmosphere curve to a flux at the Earth's surface, elucidating the impact rate of objects <1 km diameter on Earth. We find that iron meteorites >5 times 10⁴ kg (2.5 m) arrive at the Earth's surface approximately once every 50 years. Iron bodies a few meters in diameter (10⁵–10⁶ kg), which form craters ?100 m in diameter, will strike the Earth's land area every 500 years. Larger bodies will form craters 0.5 km in diameter every 20,000 years, and craters 1 km in diameter will be formed on the Earth's land area every 50,000 years. Tunguska events (low‐level atmospheric disruption of stony bolides >10⁸ kg) may occur every 500 years. Bodies capable of producing hazardous tsunami (?200 m diameter projectiles) should strike the Earth's surface every ?100,000 years. This data also allows us to assess the completeness of the terrestrial crater record for a given area over a given time interval.     

Thermal processing of chondrule precursors in planetesimal bow shocks

Abstract— We test the hypothesis that chondrules (and Type B and C calcium-aluminum-rich inclusions, CAIs) originated during passage of precursors through bow shocks upstream of planetesimals moving supersonically relative to nebula gas. A two-dimensional piecewise parabolic method (PPM) hydrocode, supplemented by a one-dimensional adiabatic shock model, is employed to simulate the postshock gas density, temperature, and velocity fields for given planetesimal sizes, velocities, and ambient nebular densities and temperatures. Thermal histories of incident silicate particles are calculated in the free molecular flow approximation by integration of the one-dimensional equations of gas-grain energy and momentum transfer. For gas number densities >10¹⁴ cm^?3, Mach numbers in the range of 4 to 5 are sufficient to melt isolated spherical particles with radii in the range 0.05 to 0.5 mm during passage of shocked gas thicknesses of 25–35 km. Minimum gas-planetesimal relative velocities are in the range 5.5–7 km/s, implying orbital eccentricities >0.2 and/or inclinations >15°. Melting of centimeter-sized CAI precursors requires either higher Mach numbers (6–7) or ambient gas densities >10¹⁵ cm^?3. For a constant radial distribution of planetesimal orbital eccentricities and inclinations, the model predicts more efficient melting of precursor particles at decreasing radial distances from the Sun where planetesimal velocities are largest. In order to process a significant fraction of solids in the nebula, planetesimals near ～2.5 AU during the chondrule formation epoch must have had a range of eccentricities and inclinations comparable to those presently observed in the residual asteroid belt. The most likely energy source for maintaining the necessary gas-planetesimal relative velocities is external gravitational perturbations associated with the forming outer planets, primarily Jupiter.     

Estimate of plasma parameters in a coronal loop by means of a fiber burst

Analyzing a fiber burst occurring in the event on 19 August, 1981, plasma parameters of the burst source volume are estimated by means of Kuijpers' (1975) fiber burst model. The derived height dependence of the density in a coronal loop agrees with an isothermal barometric loop atmosphere.     

On the Number Density of Sunyaev-Zeldovich Clusters of Galaxies

We study the conversion of a neutron star to a strange star as a possible energy source for gamma-ray bursts. We use different recent models for the equation of state of neutron star matter and strange quark matter. We show that the total amount of energy liberated in the conversion is in the range of &parl0;1-4&parr0;x1053 ergs (1 order of magnitude larger than previous estimates) and is in agreement with the energy required to power gamma-ray burst sources at cosmological distances.     

The orbit,atmospheric dynamics,and initial mass of the Park Forest meteorite

Abstract— The fireball accompanying the Park Forest meteorite fall (L5) was recorded by ground‐based videographers, satellite systems, infrasound, seismic, and acoustic instruments. This meteorite shower produced at least 18 kg of recovered fragments on the ground (Simon et al. 2004). By combining the satellite trajectory solution with precise ground‐based video recording from a single site, we have measured the original entry velocity for the meteoroid to be 19.5 ± 0.3 km/s. The earliest video recording of the fireball was made near the altitude of 82 km. The slope of the trajectory was 29° from the vertical, with a radiant azimuth (astronomical) of 21° and a terminal height measured by infrared satellite systems of 18 km. The meteoroid's orbit has a relatively large semi‐major axis of 2.53 ± 0.19 AU, large aphelion of 4.26 ± 0.38 AU, and low inclination. The fireball reached a peak absolute visual magnitude of ?22, with three major framentation episodes at the altitudes of 37, 29, and 22 km. Acoustic recordings of the fireball airwave suggest that fragmentation was a dominant process in production of sound and that some major fragments from the fireball remained supersonic to heights as low as ?10 km. Seismic and acoustic recordings show evidence of fragmentation at 42, 36, 29, and 17 km. Examination of implied energies/initial masses from all techniques (satellite optical, infrasound, seismic, modeling) leads us to conclude that the most probable initial mass was (11 ± 3) × 10³ kg, corresponding to an original energy of ?0.5 kt TNT (2.1 times 10¹² J) and a diameter of 1.8 m. These values correspond to an integral bolometric efficiency of 7 ± 2%. Early fragmentation ram pressures of <1 MPa and major fragmentations occurring with ram pressures of 2–5 MPa suggest that meter‐class stony near‐Earth asteroids (NEAs) have tensile strengths more than an order of magnitude lower than have been measured for ordinary chondrites. One implication of this observation is that the rotation period for small, fast‐rotating NEAs is likely to be >30 seconds.     

A slowly moving plasmoid associated with a filament eruption

We report the imaging observations of a slowly moving type IV burst associated with a filament eruption. This event was preceded by weak type III burst activity and was accompanied by a quasi-stationary continuum that persisted for several hours. The starting times and speeds of moving type IV burst and the erupting filament are nearly the same, implying a close physical relation between the two. The moving type IV burst is interpreted as gyrosynchrotron emission from a plasmoid containing a magnetic field of 1–2 G and nonthermal electrons of density 10⁵–10⁶ cm^–3 with a relatively low average energy of 50 keV.     

A novel approach to fireball modeling: The observable and the calculated

Estimating the mass of a meteoroid passing through the Earth's atmosphere is essential to determining potential meteorite fall positions. High‐resolution fireball images from dedicated camera networks provide the position and timing for fireball bright flight trajectories. There are two established mass determination methods: the photometric and the dynamic. A new approach is proposed, based on the dynamic method. A dynamic optimization initially constrains unknown meteoroid characteristics which are then used in a parametric model for an extended Kalman filter. The extended Kalman filter estimates the position, velocity, and mass of the meteoroid body throughout its flight, and quantitatively models uncertainties. Uncertainties have not previously been modeled so explicitly and are essential for determining fall distributions for potential meteorites. This two‐step method aims to automate the process of mass determination for application to any trajectory data set and has been applied to observations of the Bunburra Rockhole fireball. The new method naturally handles noisy raw data. Initial and terminal bright flight mass results are consistent with other works based on the established photometric method and cosmic ray analysis. A full analysis of fragmentation and the variability in the heat‐transfer coefficient will be explored in future versions of the model.     

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.     

A solar super‐flare as cause for the 14C variation in AD 774/5?

We present further considerations regarding the strong ¹⁴C variation in AD 774/5. For its cause, either a solar super‐flare or a short gamma‐ray burst were suggested. We show that all kinds of stellar or neutron star flares would be too weak for the observed energy input at Earth in AD 774/5. Even though Maehara et al. (2012) present two super‐flares with ∼10³⁵ erg of presumably solar‐type stars, we would like to caution: These two stars are poorly studied and may well be close binaries, and/or having a M‐type dwarf companion, and/or may be much younger and/or much more magnetic than the Sun – in any such case, they might not be true solar analog stars. From the frequency of large stellar flares averaged over all stellar activity phases (maybe obtained only during grand activity maxima), one can derive (a limit of) the probability for a large solar flare at a random time of normal activity: We find the probability for one flare within 3000 years to be possibly as low as 0.3 to 0.008 considering the full 1σ error range. Given the energy estimate in Miyake et al. (2012) for the AD 774/5 event, it would need to be ∼2000 stronger than the Carrington event as solar super‐flare. If the AD 774/5 event as solar flare would be beamed (to an angle of only ∼24°), 100 times lower energy would be needed. A new AD 774/5 energy estimate by Usoskin et al. (2013) with a different carbon cycle model, yielding 4 ot 6 time lower ¹⁴C production, predicts 4–6 times less energy. If both reductions are applied, the AD 774/5 event would need to be only ∼4 times stronger than the Carrington event in 1859 (if both had similar spectra). However, neither ¹⁴C nor ¹⁰Be peaks were found around AD 1859. Hence, the AD 774/5 event (as solar flare) either was not beamed that strongly, and/or it would have been much more than 4‐6 times stronger than Carrington, and/or the lower energy estimate (Usoskin et al. 2013) is not correct, and/or such solar flares cannot form (enough) ¹⁴C and ¹⁰Be. The 1956 solar energetic particle event was followed by a small decrease in directly observed cosmic rays. We conclude that large solar super‐flares remain very unlikely as the cause for the ¹⁴C increase in AD 774/5. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Interplanetary baseline observations of type III solar radio bursts

Simultaneous observations of type III radio bursts from spacecraft separated by 0.43 AU have been made using the solar orbiters HELIOS-A and HELIOS-B. The burst beginning at 19:22 UT on March 28, 1976 has been located from the intersection of the source directions measured at each spacecraft, and from burst arrival time differences. The source positions range from 0.03 AU from the Sun at 3000 kHz to 0.08 AU at 585 kHz. The electron density along the burst trajectory, and the exciter velocity (=0.13c) were determined directly, without the need to assume a density model as has been done with single-spacecraft observations. The separation of HELIOS-A and -B has also provided the first measurements of burst directivity at low frequencies. For the March 28 burst the intensity observed from near the source longitude (HELIOS-B) was 3–10 dB greater than that from 60° west of the source (HELIOS-A).