Infrasonic Observations of Meteoroids: Preliminary Results from a Coordinated Optical-radar-infrasound Observing Campaign

Recent observations using the newly installed Elginfield infrasound array in coordination with the Southern Ontario all-sky meteor camera network and Canadian Meteor Orbit Radar (CMOR) has shown that the number of meteors producing infrasound at the Earth’s surface is more frequent than previously thought. These data show the flux of meteoroids capable of producing infrasound at the ground is at least 1/month and is limited to meteors with peak visual brightness above −2. Comparisons to current meteor infrasound theory show excellent agreement with amplitude and period predictions for weakly non-linear shock waves using a realistic vertically inhomogeneous atmosphere. Similar predictions show isothermal assumptions underestimate the amplitude by orders of magnitude.     

Heating of Jupiter’s thermosphere by the dissipation of upward propagating acoustic waves

Thunderstorms in Jupiter’s atmosphere are likely to be prodigious generators of acoustic waves, as are thunderstorms in Earth’s atmosphere. Accordingly, we have used a numerical model to study the dissipation in Jupiter’s thermosphere of upward propagating acoustic waves. Model simulations are performed for a range of wave periods and horizontal wavelengths believed to characterize these acoustic waves. The possibility that the thermospheric waves observed by the Galileo Probe might be acoustic waves is also investigated. Whereas dissipating gravity waves can cool the upper thermosphere through the effects of sensible heat flux divergence, it is found that acoustic waves mainly heat the Jovian thermosphere through effects of molecular dissipation, sensible heat flux divergence, and Eulerian drift work. Only wave-induced pressure gradient work cools the atmosphere, an effect that operates at all altitudes. The sum of all effects is acoustic wave heating at all heights. Acoustic waves and gravity waves heat and cool the atmosphere in fundamentally different ways. Though the amplitudes and mechanical energy fluxes of acoustic waves are poorly constrained in Jupiter’s atmosphere, the calculations suggest that dissipating acoustic waves can locally heat the thermosphere at a significant rate, tens to a hundred Kelvins per day, and thereby account for the high temperatures of Jupiter’s upper atmosphere. It is unlikely that the waves detected by the Galileo Probe were acoustic waves; if they were, they would have heated Jupiter’s thermosphere at enormous rates.     

Parameters of infrasonic waves generated by the Chelyabinsk meteoroid on February 15, 2013

Data from infrasound stations are used to determine the basic parameters of infrasonic waves generated during the passage and airburst of the Chelyabinsk space body: time lag, duration, spectral structure, dispersion dependence, and celerity. A simulation is performed for the parameters of the infrasonic waves. A comparison between the simulation and observation results shows a good agreement.     

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.     

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.     

An entry model for the Tagish Lake fireball using seismic,satellite and infrasound records

Abstract— We present instrumental observations of the Tagish Lake fireball and interpret the observed characteristics in the context of two different models of ablation. From these models we estimate the pre‐atmospheric mass of the Tagish Lake meteoroid to be ?56 tonnes and its porosity to be between 37 and 58%, with the lowest part of this range most probable. These models further suggest that some 1300 kg of gram‐sized or larger Tagish Lake material survived ablation to reach the Earth's surface, representing an ablation loss of 97% for the fireball. Satellite recordings of the Tagish Lake fireball indicate that 1.1 times 10¹² J of optical energy were emitted by the fireball during the last 4 s of its flight. The fraction of the total kinetic energy converted to light in the satellite pass band is found to be 16%. Infrasonic observations of the airwave associated with the fireball establish a total energy for the event of 1.66 ± 0.70 kT TNT equivalent energy. The fraction of this total energy converted to acoustic signal energy is found to be between 0.10 and 0.23%. Examination of the seismic recordings of the airwave from Tagish Lake have established that the acoustic energy near the sub‐terminal point is converted to seismic body waves in the upper‐most portion of the Earth's crust. The acoustic energy to seismic energy coupling efficiency is found to be near 10^?6 for the Tagish Lake fireball. The resulting energy estimate is near 1.7 kT, corresponding to a meteoroid 4 m in diameter. The seismic record indicates extensive, nearly continuous fragmentation of the body over the height intervals from 50 to 32 km. Seismic and infrasound energy estimates are in close agreement with the pre‐atmospheric mass of 56 tonnes established from the modeling. The observed flight characteristics of the Tagish Lake fireball indicate that the bulk compressive strength of the pre‐atmospheric Tagish Lake meteoroid was near 0.25 MPa, while the material compressive strength (most appropriate to the recovered meteorites) was closer to 0.7 MPa. These are much lower than values found for fireballs of ordinary chondritic composition. The behavior of the Tagish Lake fireball suggests that it represents the lowest end of the strength spectrum of carbonaceous chondrites or the high end of cometary meteoroids. The bulk density and porosity results for the Tagish Lake meteoroid suggest that the low bulk densities measured for some small primitive bodies in the solar system may reflect physical structure dominated by microporosity rather than macroporosity and rubble‐pile assemblages.     

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.     

Arbitrary Amplitude Ion Acoustic Solitary Waves in An Unmagnetized Two Electron Population Ultra-Relativistic Dense Plasmas

The nonlinear wave structure of arbitrary amplitude ion acoustic solitary waves (IASWs) are studied in the Sagdeev’s pseudopotential framework for an ultra-relativistic degenerate dense plasma comprising cold and hot electrons and inertial ultra-cold ions. By employing standard normal-mode analysis the dispersion relation for linear waves is studied. The numerical results are presented to understand the features of ion acoustic solitary wave structures. It is shown that the present plasma model supports IASWs having positive potential well. Also, it is found that the small amplitude rarefactive double layer solution can exist in such a plasma system in some parametric region. It is shown that solitary structures and double layers are affected by relevant plasma parameters.     

Seismology of the solar atmosphere

We describe a new instrument for seismically probing the properties of the Sun's lower atmosphere, and present some first results from an observational campaign carried out at the geographic South Pole during the austral summer of 2002/2003. A preliminary analysis of the data (simultaneous, high-cadence observations of the velocity signals from the photosphere and low chromosphere) shows that the well-known suppression of acoustic power in regions of strong magnetic field, and enhancement of high-frequency power around active regions (acoustic halos), are both consistent with a spreading out of the magnetic field lines with increasing height in the atmosphere. The data have also revealed some unexpected wave behavior. First, evanescent-like waves are found at frequencies substantially above the acoustic cut-off frequency in regions of intermediate magnetic field. Second, upward- and downward-propagating waves are detected in areas of strong magnetic field such as sunspots and plage: even at frequencies below the acoustic cut-off frequency. Third, the wave behavior in regions of strong magnetic field can change over periods of a few hours from propagating to evanescent. While we have no concrete explanation for the first two results, the latter result opens up the question of whether sound waves are involved in short-term events such as flares or CME's.     

The Villalbeto de la Peña meteorite fall: I. Fireball energy,meteorite recovery,strewn field,and petrography

Abstract— An impressive daylight fireball was observed from Spain, Portugal, and the south of France at 16h46m45s UTC on January 4, 2004. The meteoroid penetrated into the atmosphere, generating shock waves that reached the ground and produced audible booms. The associated airwave was recorded at a seismic station located 90 km north of the fireball trajectory in Spain, and at an infrasound station in France located 750 km north‐east of the fireball. The absolute magnitude of the bolide has been determined to be ?18 ± 1 from a casual video record. The energy released in the atmosphere determined from photometric, seismic, and infrasound data was about 0.02 kilotons (kt). A massive fragmentation occurred at a height of 28 ± 0.2 km, resulting in a meteorite strewn field of 20 × 6 km. The first meteorite specimen was found on January 11, 2004, near the village of Villalbeto de la Peña, in northern Palencia (Spain). To date, about 4.6 kg of meteorite mass have been recovered during several recovery campaigns. The meteorite is a moderately shocked (S4) L6 ordinary chondrite with a cosmic‐ray‐exposure age of 48 ± 5 Ma. Radioisotope analysis shows that the original body had a mass of 760 ± 150 kg, which is in agreement with the estimated mass obtained from photometric and seismic measurements.     

Stardust—An artificial,low‐velocity “meteor” fall and recovery: 15 January 2006

Abstract— On January 15, 2006, Stardust, a man‐made space capsule, plummeted to Earth for a soft landing after spending seven years in space. Since the expected initial speed of the body was about 12.9 km/s, a four‐element ground‐based infrasound array was deployed to Wendover, Nevada, USA, to measure the hypersonic booms from the re‐entry. At a distance of ～33 km from the nominal trajectory, we easily recorded the weak acoustic arrivals and their continued rumbling after the main hypersonic boom arrival. In this paper, we report on subsequent analyses of these data, including an assessment of the expected entry characteristics (dynamics, energetics, ablation and panchromatic luminosity, etc.) on the basis of a bolide/meteor/fireball entry model that was specifically adapted for modeling a re‐entering man‐made object. Throughout the infrasonic data analyses, we compared our results for Stardust to those previously obtained for Genesis. From the associated entry parameters, we were also able to compute the kinetic energy density conservation properties for the propagating line source blast wave and compared the inviscid theoretical predictions against observed ground‐based infrasound amplitude and wave period data as a function of range. Finally, we made a top‐down bottom‐up assessment of the line source wave normals propagating downward into the complex temperature/sound speed and horizontal wind speed environment during January 15, 2006. This assessment proved to be generally consistent with the signal processing analysis and with the observed time delay between the known Stardust entry and the time of observations of infrasound signals, and so forth.     

Atmospheric acoustics of Titan, Mars, Venus, and Earth

Planetary atmospheres are complex dynamical systems whose structure, composition, and dynamics intimately affect the propagation of sound. Thus, acoustic waves, being coupled directly to the medium, can effectively probe planetary environments. Here we show how the acoustic absorption and speed of sound in the atmospheres of Venus, Mars, Titan, and Earth (as predicted by a recent molecular acoustics model) mirror the different environments. Starting at the surface, where the sound speed ranges from ∼200 m/s for Titan to ∼410 m/s for Venus, the vertical sound speed profiles reveal differences in the atmospheres' thermal layering and composition. The absorption profiles are relatively smooth for Mars, Titan, and Earth while Venus stands out with a noticeable attenuation dip occurring between 40 and 100 km. We also simulate a descent module sampling the sound field produced by a low-frequency “event” near the surface noting the occurrence of acoustic quiet zones.     

彗木相撞中的波动

1994年7月18日至24日期间，彗星苏梅克-列维9(SL-9)的超过20块碎片与木星发生了相撞．哈勃空间望远镜(HST)拍摄到的图像揭示了木星大气对撞击的动力学响应．具有重要意义的是观测到5个撞击点周围的圆环，它们以450米／秒的常速度向外运动．环的圆形性表明它们是波．因为对于不同大小的撞击，波速是常量，可以推断出传播速度与爆炸能量无关．这意味着这些波动是线性波．评述现行理论所使用的3类候选波，亦即惯性引力波、声波和地震波，介绍的重点是前面两种．     

On the speed and shape of electron acoustic solitary waves

Electron acoustic solitary waves in a collisionless plasma consisting of a cold electron fluid and non-thermal hot electrons are investigated by a direct analysis of the field equations. The Sagdeev potential is obtained in terms of electron acoustic speed by simply solving an algebraic equation. It is found that the amplitude and width of the electron acoustic solitary waves as well as the parametric regime where the solitons can exist are very sensitive to the population of energetic non-thermal hot electrons. The soliton and double layer solutions are obtained as a small-amplitude approximation. The effect of non-thermal hot electrons is found to significantly change the properties of the electron acoustic solitary waves (EAWs). A comparison with the Viking Satellite observations in the day side auroral zone is also discussed.     

Solar atmospheric heating

The solar atmosphere is heated by a flux of mechanical waves propagating in one or more of the modes: acoustic, Alfvén and gravitational.The acoustic theory is compared with observational data and found inadequate. First, the theory is based quantitatively on the Böhm-Vitense convection zone model, and large-scale convective motions (supergranulation) and magnetic fields (unipolar regions) show that convection has another form. On the other hand, when granular motions are invoked the energy flux is too small. Second, atmospheric heating is localized in faculae, and enhanced acoustic flux beneath these regions is no longer explicable. Finally, the short periods of 10–30 s invoked recently appear inexplicable. Objections to the gravitational wave heating process are given briefly.Previous objections to Alfvén waves as an energy source followed from the belief that fields were generally uniform and of strength 50 G, now known to be incorrect. Models of Alfvén wave generation are based on (i) granule eddy motions, (ii) overstable oscillations in subsurface flux tubes and sunspot flux ropes, and (iii) supergranule motions, both horizontal and vertical.The first provides waves which propagate along thin flux tubes oscillating as taut wires in a compressible fluid; they may explain mottles, fibrils and other small emission features. The second may explain the enormous dissipation in spot groups, including flares. The third has been invoked earlier to explain spicules, and may have effects in the solar wind.     

Effect of electron temperature on small-amplitude electron acoustic solitary waves in non-planar geometry

Effects of electron temperature on the propagation of electron acoustic solitary waves in plasma with stationary ions, cold and superthermal hot electrons is investigated in non-planar geometry employing reductive perturbation method. Modified Korteweg–de Vries equation is derived in the small amplitude approximation limit. The analytical and numerical calculations of the KdV equation reveal that the phase velocity of the electron acoustic waves increases as one goes from planar to non planar geometry. It is shown that the electron temperature ratio changes the width and amplitude of the solitary waves and when electron temperature is not taken into account,our results completely agree with the results of Javidan & Pakzad (2012). It is found that at small values of \(\tau \), solitary wave structures behave differently in cylindrical (\(\text {m} = 1\)), spherical (\(\text {m} = 2\)) and planar geometry (\(\text {m} = 0\)) but looks similar at large values of \(\tau \). These results may be useful to understand the solitary wave characteristics in laboratory and space environments where the plasma have multiple temperature electrons.     

Dynamical implications of Si iv line profiles from OSO-8 observations

The widths and rms variations of line centers for 1393 profiles from the Si iv ion in the solar transition zone, as observed by OSO-8, are analyzed to give the amplitudes and periods of three postulated types of disturbance: sinusoidal acoustic waves; sinusoidal acoustic shocks; and normally propagated MHD waves. All three assumptions lead to mean intervals between disturbances of 40–50 s, but acoustic waves and shocks are rules out because they predict intensifications of brightness far in excess of anything observed. MHD disturbances in magnetic fields 2 G are consistent with observations, but normally propagated disturbances also give too great an intensification in brightness. Disturbances with wave vectors at an angle to the magnetic field are suggested.     

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.     

Nonthermal Velocities in the Solar Transition Zone and Corona

Nonthermal velocities are presented for spectral lines covering the temperature range 10 4–10 6 K, measured from high-spectral-resolution data for several solar features observed at the limb by the high resolution telescope and spectrograph (HRTS), including a coronal hole, quiescent regions and several small-scale active regions. These results are compared with predictions based on acoustic waves and heating via Alfvén waves. It is likely that more than one mechanism is operating simultaneously, in particular, resonant Alfvén wave heating, which is very sensitive to background plasma motions.     

彗木相撞中的惯性引力波

苏梅克-列维9号彗星(SL9)与木星相撞后，在木星上观测到的以常速度(～450m／s)向外扩展的圆环意味着这是碰撞在木星大气中引起的线性波动．我们选取：非旋转、无粘性、密度分层、不可压缩的木星大气模型，而且木星大气以水平速度U=b az运动；给出初始扰动压力P(r；0)作为碰撞的初始条件，用流体力学方程组求解了彗木相撞中的惯性引力波．结果表明：当木星大气以速度U=U0(～170m／s)运动时，彗星碎片的大部分能量都用来产生内波，同时还得到彗星碎片的撞击深度H与水平相速Vp的关系式．当木星大气以速度为U=b az运动时，木星大气的扰动能量不再是在动能和势能间均分。