On the possibility of a secondary explosion in the antipode of a large meteorite impact point

At large distances, due to atmospheric absorption and the dispersion of high-frequency components, the airwaves from the fall of large meteorites or heavy-yield explosions are transformed into an infrasonic wave train propagating over large distances via atmospheric sound channels. In approaching the antipode, the amplitude of infrasonic oscillations increases significantly and the nonlinear effects may trigger the formation of a blast wave, that is, another explosion. The condition which allows such a phenomenon to happen was obtained in this study. Infrasonic waves from the Tunguska fall event and waves generated by the largest nuclear explosions were considered in this study.     

Propagation of auroral infrasonic waves studied by ray-tracings

On the basis of a ray tracing method the propagation and the attenuation of an auroral infrasonic wave are studied. Relations between the direct and reflected waves recorded at the Syowa Station, Antarctica, are clarified with regard to; (1) the delay time, (2) the intensity ratio, and (3) trace velocities. The time required for a wave to travel from the source to the ground is calculated as a function of a source altitude. The retardation time of the wave arrival behind the zenith crossing of the source current is deduced. A method is proposed for estimating the altitude of a source current from the retardation and a trace velocity of the wave. It is concluded that the existence of a supersonic equatorward motion of an electrojet which continues for a certain distance is necessary for the observation of auroral infrasonic waves. This distance must exceed at least 60 km equatorwards from the zenith to enable the direct wave to be observed and with total length of 930 km to enable the reflected wave to be observed. From these conditions it is also concluded that the infrasonic wave is not seen in mid latitudes and the reflected wave is a rare phenomenon.     

Parameters of Acoustic Signals Generated by the Atmospheric Meteoroid Explosion over Romania on January 7, 2015

System spectral analysis of temporal variations in the level of acoustic signals recorded at a number of European infrasound stations is carried out. The prevailing periods that varied within 3–5 s are found. Initial kinetic and acoustic energies of the Romanian meteoroid and infrasound parameters are calculated: celerity (280 m/s), acoustic efficiency, and stratospheric wind velocity (about 20 m/s). Parameters of cylindrical and explosive shock waves are found: duration, characteristic size, and amplitude. The dependence of the attenuation of infrasonic waves on distance is estimated. The results of the estimates are in good agreement with the results of the observations.     

An auroral infrasonic substorm investigation using Chatanika radar and other geophysical sensors

Characteristics of the supersonic auroral arcs within the 0905 UT 2 April 1973 substorm were determined using data from (1) all-sky cameras; (2), surface magnetometers, (3) multispectral scanning photometers, (4) 30MHz riometers, (5) Chatanika incoherent-scatter radar, (6) Homer auroral radar, and (7) infrasonic microphone arrays at College and Stevens Village in Alaska. These data were analyzed to determine the properties of an auroral electrojet arc that generates auroral infrasonic waves (AIW).

An arc that was show to be the source of an AIW was found to have the following characteristics: (1) a velocity of 500 m/sec traveling from an azimuth of 350°; (2) an intensity in 4278 A of 26 Kr, (3) a maximum electron density of 2.8 × 10⁶ el/cm⁶ at 100km height, (4) an equivalent westward line current of 2.8 × 10⁶ A, (5) orientation of ΔH parallel to the AIW direction of travel and perpendicular to the arc's long axis, (6) a characteristic energy of the primary auroral electron spectrum of 3.0keV, and (7) an energy deposition rate for the auroral pdarticles of 100 erg/cm² sec.     

The Morávka meteorite fall: 2. Interpretation of infrasonic and seismic data

Abstract— The sound production from the Morávka fireball has been examined in detail making use of infrasound and seismic data. A detailed analysis of the production and propagation of sonic waves during the atmospheric entry of the Morávka meteoroid demonstrates that the acoustic energy was produced both by the hypersonic flight of the meteoroid (producing a cylindrical blast wave) and by individual fragmentation events of the meteoroid, which acted as small explosions (producing quasispherical shock waves). The deviation of the ray normals for the fragmentation events was found to be as much as 30° beyond that expected from a purely cylindrical line source blast. The main fragmentation of the bolide was confined to heights above 30 km with a possible maximum in acoustic energy production near 38 km. Seismic stations recorded both the direct arrival of the airwaves (the strongest signal) as well as air‐coupled P‐waves and Rayleigh waves (earlier signals). In addition, deep underground stations detected the seismic signature of the fireball. The seismic data alone permit reconstruction of the fireball trajectory to a precision on the order of a few degrees. The velocity of the meteoroid is much less well‐determined by these seismic data. The more distant infrasonic station detected 3 distinct signals from the fireball, identified as a thermospheric return, a stratospheric return, and an unusual mode propagating through the stratosphere horizontally and then leaking to the receiver.     

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.     

用粒子模拟方法研究离子声波

<正>向传播朗缪尔波被离子声波散射是太阳射电Ⅲ型暴基波和谐波激发的重要过程.使用粒子模拟方法对电子束流激发朗缪尔波的过程进行了模拟,同时对产生的反向朗缪尔波、朗缪尔波2次谐波和朗缪尔波通过非线性过程产生的离子声波的性质进行了分析研究.为了更好地研究离子声波,模拟时单独计算了由离子扰动引起的电场.模拟计算了不同初始参数下产生的离子声波强度,发现离子的温度和质量对离子声波的产生有重要作用,验证了反向朗缪尔波与离子声波的相关性.同时在模拟中验证了朗缪尔波的衰变过程,确认了离子声波对反向朗缪尔波的放大作用.     

A Simulation Study of the Plasma Wave Enhancements in the Earth’s Equatorial Plasmasphere

In the equatorial plasmasphere, plasma waves are frequently observed. To improve our understanding of the mechanism generating plasma waves from instabilities, a comparison of observations, linear growth-rate calculations, and simulation results is presented. To start the numerical experiments from realistic initial plasma conditions, we use the initial parameters inferred from observational data obtained around the plasma-wave generation region by the Akebono satellite. The linear growth rates of waves of different modes are calculated under resonance conditions, and compared with simulation results and observations. By employing numerical experiments by a particle code, we first show that upper hybrid-, Z-, and whistler-mode waves are excited through instabilities driven by a ring-type velocity distribution. The simulation results suggest a possibility that energetic electrons with energies of some tens of keV confined around the geomagnetic equator are responsible for the observed enhancements of Z- and whistler-mode waves. While the comparison between linear growth-rate calculations and observations shows the different tendency of wave amplitude of Z-mode and whistler-mode waves, the wave amplitude of these wave modes in the simulation results is consistent with the observation.     

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.     

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.     

Plasma waves caused by transient heat conduction in a coronal loop and anomalous resistivity

Numerical simulation is made of the transient heat conduction during local heating in a model coronal magnetic loop with an axial electric current. It is assumed that a segment near the top of the normal coronal loop is heated to above 10⁷ K by a sufficiently small heat input as compared with the total flare energy. A hump appears in the velocity distribution of electrons moving down the temperature gradient with speeds slightly below the thermal one. Consequently, electron plasma waves are excited. The high intensity of the waves persists in the upper region of the loop for more than a second until the termination of the simulation. The energy density of the plasma waves normalized with respect to thermal density is 10^–3.5 at maximum. A theoretical estimate gives an anomalous resistivity 5 orders of magnitude greater than an initial value. Based on the above result, we propose a model for impulsive loop flares.     

Formation of fast magnetosonic shock waves during a two-current-loop collision in solar flares

We present a numerical simulation of the fast magnetosonic shock wave formation during a two-current-loop collision by using a magnetohydrodynamical model. It is shown that the rarefaction waves are generated in the initial stage when the two current loops start to collide. After the rarefaction waves propagate away from the excited region, the fast magnetosonic waves with density enhancement can be produced for the simulation when the current strength of the loop is weak. As the current becomes strong enough, the magnetosonic shock waves can be generated in the direction perpendicular to that of the two-loop collision.     

Cross-scale coupling at a perpendicular collisionless shock

A full particle simulation study is carried out on a perpendicular collisionless shock with a relatively low Alfven Mach number (M_A = 5). Recent self-consistent hybrid and full particle simulations have demonstrated ion kinetics are essential for the non-stationarity of perpendicular collisionless shocks, which means that physical processes due to ion kinetics modify the shock jump condition for fluid plasmas. This is a cross-scale coupling between fluid dynamics and ion kinetics. On the other hand, it is not easy to study cross-scale coupling of electron kinetics with ion kinetics or fluid dynamics, because it is a heavy task to conduct large-scale full particle simulations of collisionless shocks. In the present study, we have performed a two-dimensional (2D) electromagnetic full particle simulation with a “shock-rest-frame model”. The simulation domain is taken to be larger than the ion inertial length in order to include full kinetics of both electrons and ions. The present simulation result has confirmed the transition of shock structures from the cyclic self-reformation to the quasi-stationary shock front. During the transition, electrons and ions are thermalized in the direction parallel to the shock magnetic field. Ions are thermalized by low-frequency electromagnetic waves (or rippled structures) excited by strong ion temperature anisotropy at the shock foot, while electrons are thermalized by high-frequency electromagnetic waves (or whistler mode waves) excited by electron temperature anisotropy at the shock overshoot. Ion acoustic waves are also excited at the shock overshoot where the electron parallel temperature becomes higher than the ion parallel temperature. We expect that ion acoustic waves are responsible for parallel diffusion of both electrons and ions, and that a cross-scale coupling between an ion-scale mesoscopic instability and an electron-scale microscopic instability is important for structures and dynamics of a collisionless perpendicular shock.     

Nonlinear electrostatic coherent structures: solitary and shock waves in a dissipative, nonplanar multi-component quantum plasma

The nonlinear propagation of ion-acoustic solitary and shock waves in a dissipative, nonplanar quantum plasma comprised of electrons, positrons, and ions are studied. A modified Korteweg-de Vries Burgers equation is derived in the limit of low frequency and long wavelength by taking into account the kinematic viscosity among the plasma constituents. It is shown that this plasma system supports the propagation of both compressive and rarefactive nonlinear waves. The effects of variation of various plasma parameters on the time evolution of nonplanar solitary waves, the profile of shock waves, and the nonlinear structure induced by the collision of solitary waves are discussed. It is found that these parameters have significant effects on the properties of nonlinear waves in cylindrical and spherical geometries, and these effects for compressive and rarefactive nonlinear waves are obviously different.     

The emission and absorption of waves by charged particles in magnetized plasmas

A semiclassical theory describing the emission and absorption of waves is applied to the interaction of charged particles with waves in magnetized plasmas. Spontaneous emission of all cold plasma wave modes is calculated in detail. The method gives the absorption coefficient for the waves and a diffusion equation in momentum space for the particles describing the effects of the induced processes.Coefficients describing the systematic change of particle parameters follow from the diffusion equation. Applications of astrophysical interest are outlined.     

A new numerical method for simulating the solar wind Alfvén waves: Development of the Vlasov-MHD model

A novel scheme of plasma simulation particularly suited for computing the one-dimensional nonlinear evolution of parallel propagating solar wind Alfvén waves is presented. The scheme is based on the Vlasov and the MHD models, for solving the longitudinal and the transverse components, respectively. As long as the nonlinearity is not very large (so that the longitudinal and transverse components are well separated), our Vlasov-MHD model can correctly describe evolution of finite amplitude parallel Alfvén waves, which are typical in the solar wind, both in the linear and nonlinear stages. The present model can be applied to discussions of phenomena where the parallel Alfvén waves play major roles, for example, the solar coronal heating and solar wind acceleration by the Alfvén waves propagating from the photosphere.     

Compound Effect of Alfvén Waves and Ion-Cyclotron Waves on Heating/Acceleration of Minor Ions via the Pickup Process

A scenario is proposed to explain the preferential heating of minor ions and differential-streaming velocity between minor ions and protons observed in the solar corona and in the solar wind. It is demonstrated by test-particle simulations that minor ions can be nearly fully picked up by intrinsic Alfvén-cyclotron waves observed in the solar wind based on the observed wave energy density. Both high-frequency ion-cyclotron waves and low-frequency Alfvén waves play crucial roles in the pickup process. A minor ion can first gain a high magnetic moment through the resonant wave–particle interaction with ion-cyclotron waves, and then this ion with a large magnetic moment can be trapped by magnetic mirror-like field structures in the presence of the low-frequency Alfvén waves. As a result, the ion is picked up by these Alfvén-cyclotron waves. However, minor ions can only be partially picked up in the corona because of the low wave energy density and low plasma β. During the pickup process, minor ions are stochastically heated and accelerated by Alfvén-cyclotron waves so that they are hotter and flow faster than protons. The compound effect of Alfvén waves and ion-cyclotron waves is important in the heating and acceleration of minor ions. The kinetic properties of minor ions from simulation results are generally consistent with in-situ and remote features observed in the solar wind and solar corona.     

Kadomtsev–Petviashvili equation for solitary waves in warm dense astrophysical electron-positron-ion plasmas

Theoretical investigation is carried out to understand the dynamics and stability of three dimensional ion solitary waves propagating in dense plasma comprising of ultra-relativistic degenerate electrons and positrons and warm ions. A linear dispersion relation is derived which shows a strong dependence of wave on positron concentration (through the change of density balance) and ion-to-degenerate electron temperature ratio. A nonlinear Kadomtsev-Petviashvili equation is derived by employing the reductive perturbation technique and solved analytically and the conditions for existence of stable solitary waves are found. The analysis reveals that only compressive solitary waves exist in the system. Effects of the change of density balance and Fermi temperature ratios are studied in detail, both analytically and numerically. Furthermore, the conditions for stable solitary waves are discussed by using energy consideration method. The numerical results are also presented by using the parameters consistent with the degenerate and ultrarelativistic astrophysical plasmas.     

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.     

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