skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Atomic‐Scale Simulations of Meteor Ablation
Abstract Meteoroids smaller than a microgram constantly bombard the Earth, depositing material in the mesosphere and lower thermosphere. Meteoroid ablation, the explosive evaporation of meteoroids due to erosive impacts of atmospheric particles, consists of sputtering and thermal ablation. This paper presents the first atomic‐scale modeling of sputtering, the initial stage of ablation where hypersonic collisions between the meteoroid and atmospheric particles cause the direct ejection of atoms from the meteoroid surface. Because meteoroids gain thermal energy from these particle impacts, these interactions are important for thermal ablation as well. In this study, a molecular dynamics simulator calculates the energy distribution of the sputtered particles as a function of the species, velocity, and angle of the incoming atmospheric particles. The sputtering yield generally agrees with semi‐empirical equations at normal incidence but disagrees with the generally accepted angular dependence.Λ, the fraction of energy from a single atmospheric particle impact incorporated into the meteoroid, was found to be less than 1 and dependent on the velocity, angle, atmospheric species, and meteoroid material. Applying this newΛto an ablation model results in a slower meteoroid temperature increase and mass loss rate as a function of altitude. This alteration results in changes in the expected electron line densities and visual magnitudes of meteoroids. Notably, this analysis leads to the prediction that meteoroids will generally ablate 1–4 km lower than previously predicted. This affects analysis of radar and visual measurements, as well as determination of meteoroid mass.  more » « less
Award ID(s):
1755020 1755350
PAR ID:
10449781
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Space Physics
Volume:
125
Issue:
9
ISSN:
2169-9380
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract High‐power large‐aperture radar instruments observe numerous meteor head echoes per minute. Head echoes result from reflections of radio waves from plasma surrounding meteoroids as they enter Earth's atmosphere. Knowledge of the spatial distribution of electrons in this plasma is essential to determining the mass loss rate of the meteor as a function of its measured radar cross‐section. Prior work applies theoretical and computational methods to determine the electron density distribution, but assumes the meteoroid emits neutral particles uniformly across its surface. In this paper, a numerical surface ablation model demonstrates that meteoroid mass loss may occur preferentially in the direction facing the oncoming atmosphere. Specifically, meteoroid mass loss becomes proportional to the frontal surface area facing the freestream atmosphere in the limit of high Biot number, but remains isotropic in the limit of low Biot number. Meteoroid rotation has a small effect on the direction of ejected mass, but the effect is insignificant compared to variation in meteoroid properties that affect the Biot number. This result informs our computational meteor plasma model, in which we compare the effect of meteoroid vaporization on the plasma distribution in the limits of low versus high Biot number. The resulting electron density profiles demonstrate order‐of‐magnitude agreement between each other, with peak difference of 70% immediately upstream of the meteoroid. This implies that the directional distribution of vaporizing neutrals likely does not significantly influence head echo observations, lending credence to existing work that assumes isotropic ablation. 
    more » « less
  2. Abstract Angular response functions are derived for four electron channels and six proton channels of the SEM‐2 MEPED particle telescopes on the POES and MetOp satellites from Geant4 simulations previously used to derive the energy response. They are combined with model electron distributions in energy and pitch angle to show that the vertical 0° telescope, intended to measure precipitating electrons, instead usually measures trapped or quasi‐trapped electrons, except during times of enhanced pitch angle diffusion. A simplified dynamical model of the radiation belt electron distribution near the loss cone, as a function of longitude, energy, and pitch angle, that accounts for pitch angle diffusion, azimuthal drift, and atmospheric backscatter is fit to sample MEPED electron data atL = 4during times of differing diffusion rates. It is then used to compute precipitating electron flux, as function of energy and longitude, that is lower than would be estimated by assuming that the 0° telescope always measures precipitating electrons. 
    more » « less
  3. Abstract To date, there is no consensus on the probability distribution of particle velocities during bedload transport, with some studies suggesting an exponential‐like distribution while others a Gaussian‐like distribution. Yet, the form of this distribution is key for the determination of sediment flux and the dispersion characteristics of tracers in rivers. Combining theoretical analysis of the Fokker‐Planck equation for particle motions, numerical simulations of the corresponding Langevin equation, and measurements of motion in high‐speed imagery from particle‐tracking experiments, we examine the statistics of bedload particle trajectories, revealing a two‐regime distance‐time (L‐Tp) scaling for the particle hops (measured from start to stop). We show that particles of short hop distances scale asL~giving rise to the Weibull‐like front of the hop distance distribution, while particles of long hop distances transition to a different scaling regime ofL~Tpleading to the exponential‐like tail of the hop distance distribution. By demonstrating that the predominance of mostly long hop particles results in a Gaussian‐like velocity distribution, while a mixture of both short and long hop distance particles leads to an exponential‐like velocity distribution, we argue that the form of the probability distribution of particle velocities can depend on the physical environment within which particle transport occurs, explaining and unifying disparate views on particle velocity statistics reported in the literature. 
    more » « less
  4. ABSTRACT This work presents the result of sporadic meteor radiant density distribution using the Arecibo 430 MHz incoherent scatter radar (ISR) located in Puerto Rico for the first time. Although numerous meteor studies have been carried out using the Arecibo ISR, meteoroid radiant density distribution has remained a mystery as the Arecibo radar cannot measure vector velocity. A numerical orbital simulation algorithm using dynamic programming and stochastic gradient descent is designed to solve the sporadic meteoroid radiant density and the corresponding speed distributions of the meteors observed at Arecibo. The data set for the algorithm comprises over 250 000 meteors from Arecibo observations between 2009 and 2017. Five of the six recognized sporadic meteor sources can be identified from our result. There is no clearly identifiable South Apex source. Instead, there is a broad distribution between +/−30° ecliptic latitude, with the peak density located in the North Apex direction. Our results also indicate that the Arecibo radar is not sensitive to meteors travelling straight into or perpendicular to the antenna beam but is most sensitive to meteors with an arrival angle between 30° and 60°. Our analysis indicates that about 75 per cent of meteoroids observed by the Arecibo radar travel in prograde orbits when the impact probability is considered. Most of the retrograde meteoroids travel in inclined low-eccentricity orbits. 
    more » « less
  5. Abstract Accurate specification of ionization production by energetic electron precipitation is critical for atmospheric chemistry models to assess the resultant atmospheric effects. Recent model‐observation comparison studies have increasingly highlighted the importance of considering precipitation fluxes in the full range of electron energy and pitch angle. However, previous parameterization methods were mostly proposed for isotropically precipitation electrons with energies up to 1 MeV, and the pitch angle dependence has not yet been parameterized. In this paper, we first characterize and tabulate the atmospheric ionization response to monoenergetic electrons with different pitch angles and energies between∼3 keV and∼33 MeV. A generalized method that fully accounts for the dependence of ionization production on background atmospheric conditions, electron energy, and pitch angle has been developed based on the parameterization method of Fang et al. (2010,https://doi.org/10.1029/2010GL045406). Moreover, we validate this method using 100 random atmospheric profiles and precipitation fluxes with monoenergetic and exponential energy distributions, and isotropic and sine pitch angle distributions. In a suite of 6,100 validation tests, the error in peak ionization altitude is found to be within 1 km in 91% of all the tests with a mean error of 2.7% in peak ionization rate and 1.9% in total ionization. This method therefore provides a reliable means to convert space‐measured precipitation energy and pitch angle distributions into ionization inputs for atmospheric chemistry models. 
    more » « less