Issue 
A&A
Volume 659, March 2022



Article Number  A120  
Number of page(s)  9  
Section  Celestial mechanics and astrometry  
DOI  https://doi.org/10.1051/00046361/202140810  
Published online  29 March 2022 
Distribution of dust ejected from the lunar surface into the EarthMoon system
^{1}
Department of Engineering Mechanics, Shandong University,
250061
Jinan,
PR China
^{2}
Institute of Geological Sciences, Freie Universität Berlin,
Berlin, Germany
^{3}
Space Physics and Astronomy Research Unit, University of Oulu,
90014
Oulu, Finland
^{4}
School of Aeronautics and Astronautics, Sun Yatsen University,
Shenzhen Campus,
518107
Shenzhen, PR China
email: liuxd36@mail.sysu.edu.cn
Received:
15
March
2021
Accepted:
30
January
2022
Aims. An asymmetric dust cloud was detected around the Moon by the Lunar Dust Experiment on board the Lunar Atmosphere and Dust Environment Explorer mission. We investigate the dynamics of the grains that escape the Moon and their configuration in the EarthMoon system.
Methods. We use a plausible initial ejecta distribution and mass production rate for the ejected dust. Various forces, including the solar radiation pressure and the gravity of the Moon, Earth, and Sun, are considered in the dynamical model, and direct numerical integrations of trajectories of dust particles are performed. The final states, the average life spans, and the fraction of retrograde grains as functions of particle size are computed. The number density distribution in the EarthMoon system is obtained through longterm simulations.
Results. The average life spans depend on the size of dust particles and show a rapid increase in the size range between 1 and 10 μm. About 3.6 × 10^{−3} kg s^{−1} (~2%) particles ejected from the lunar surface escape the gravity of the Moon, and they form an asymmetric torus between the Earth and the Moon in the range [10 R_{E}, 50 R_{E}], which is offset toward the direction of the Sun. A considerable number of retrograde particles occur in the EarthMoon system.
Key words: planets and satellites: rings / zodiacal dust / Moon / Earth / celestial mechanics / methods: numerical
© ESO 2022
1 Introduction
In the past, studies on the configuration of circumplanetary dust mainly focused on the grain particles around the giant planets and Mars (e.g., Tiscareno & Murray 2018; Spahn et al. 2019). The dynamics of dust particles in the EarthMoon system has been less explored. Although the Munich Dust Counter did not detect the lunar ejecta dust cloud (Iglseder et al. 1996), a permanent and asymmetric dust cloud engulfing the Moon was indeed detected by the Lunar Dust Experiment (LDEX; Horányi et al. 2015). From the LDEX measurements, the size distribution of the dust grains was derived, and the density of particles was found to drop with altitude and to vary azimuthally around the Moon, reaching a peak at 5–7 h lunar local time. Six previously identified meteoroid populations, helion, antihelion, apex, antiapex, northern toroidal, and southern toroidal, were found to be plausible interplanetary projectiles that generate the lunar dust cloud via highspeed impacts (Szalay et al. 2019). The asymmetric nature of this dust cloud, the result of the asymmetric impactor flux, was analyzed, and a consistent mass production rate per unit surface as a function of impact direction was derived. Szalay & Horányi (2016) also analyzed particles that fall back to the lunar surface, which follow a cumulative size distribution with exponent 2.7. The dynamics of dust particles ejected from the Moon was previously studied by Colombo et al. (1966), and dust orbiting the nearEarth environment was studied by Peale (1966).
Investigating the lunar dust distributions at different times, heights, and positions was ever one of optional tasks of Chinese Lunar Exploration Mission, Chang’e 4 (Wang & Liu 2016). At the end of 2020, Chang’e 5 launched successfully. The dust distribution can be obtained from dust dynamics modeling; as such, the total dust flux received by the spacecraft can be calculated, and the dust hazard for a specific mission can be evaluated. Space activities in the EarthMoon system are much more frequent compared with other planets, and thus the study of dust particles is of great significance for assessing the space environment and ensuring the security of explorations.
In this work, we focus on the dynamics and distribution of particles ejected from the surface of the Moon that escape the Moon’s gravity. The paper is structured as follows. In Sect. 2 the dynamical model for the motion of particles ejected from the lunar surface is presented, taking the solar radiation pressure and gravity of the Sun, the Earth, and the Moon into consideration. In Sect. 3 the mass production rate and the initial ejecta distribution on the surface of the Moon are deduced. In Sect. 4 the detailed simulation results are presented, with the distribution of sinks and the average life spans of particles as functions of size. The dust number density distribution in the Earthcentered inertial frame is also presented. Finally, the distributions and the evolution of osculating orbital elements of dust particles are given.
2 Dynamical model
Dust particles in interplanetary and circumplanetary space can be affected by various forces, including solar gravity, solar radiation pressure, PoyntingRobertson drag, the Lorentz force, and the gravity of the planets and other bodies in the Solar System. For particles ejected from the surface of the Moon, solar radiation pressure and the gravity of the Moon, Earth, and Sun are considered in this paper. We estimated the strength of the perturbation induced by the Lorentz force for typical strengths of the magnetic field and grain charges in the interplanetary magnetic field, the Earth’s plasma sphere, and the geomagnetic tail of the Earth. We find that outside 10 R_{E} the strength of the Lorentz force is small compared to the perturbations induced by solar radiation pressure, by at least two orders of magnitude for grains larger than 10 μm; even for grains of 1 μm, solar radiation pressure dominates by more than an order of magnitude. An Earthcentered J 2000 inertial frame Oxyz is used in our simulation. Here, the x axis points in the direction of vernal equinox at the J2000 epoch, the z axis is the normal of the Earth’s equatorial plane (north), and the y axis is determined by the righthanded rule.
Physical parameters for the Sun, Earth, and Moon used in this paper are shown in Table 1. All the ephemerides of the Sun, Earth, and Moon in our model are taken from NAIF SPICE toolkit^{1}.
The equation of motion of one particle ejected from the surface of the Moon reads (1)
where r is the Earthcentric radius vector of thedust particle and is the acceleration caused by the gravity of the Earth, (2)
Here, G is the gravitational constant, M_{E} is the mass of the Earth, R_{E} is the reference radius of Earth, θ is the colatitude in an Earthcentered bodyfixed frame, and P_{2} is the Legendre function of degree 2. In our simulation only the seconddegree zonal harmonic, J_{2} ≈ 1.082 × 10^{−3} (Pavlis et al. 2012), is considered. The variable denotes the acceleration due to the solar radiation pressure (Burns et al. 1979), (3)
where Q_{c} = 1.36 × 10^{3} W∕m^{2} is the solar radiation energy flux at one astronomical unit (AU) distance and Q_{pr} is the solar radiation pressure efficiency factor, which depends on the size and material of the dust particle. The value of Q_{pr} was calculated from Mie theory for spherical particles using optical constants for silicates taken from Mukai (1989) (see our Fig. 1). The symbol ρ_{g} = 3500 kg m^{−3} denotes volumetric mass density for lunar dust (Solomon 1974), r_{g} is the radius of the particle, c is the speed of light, r_{S} is the vector from the Earth to the Sun, and is the unit vector from the Sun to the particle.
The acceleration by the gravity of Sun and Moon reads (e.g., Murray & Dermott 1999; Liu & Schmidt 2019) (4)
where M_{S} and M_{M} are the masses of the Sun and the Moon. The symbol r_{pS} denotes the vector from the dust particleto the Sun, r_{pM} is the vector from the dust particle to the Moon, and r_{M} is the vectorfrom the Earth to the Moon.
In Sect. 4 we analyze the contribution of different perturbation forces to the evolution of eccentricity and inclination, d e∕d t and d i∕d t, which can be calculated from the perturbation equations as (Murray & Dermott 1999) (5) (6)
Here, , , and are the magnitudes of the radial, tangential, and normal components of the acceleration, respectively, which we determined from the equations of motion (Eq. (1)). The symbols a, e, i, ω, f, and E denote the osculating semimajor axis, eccentricity, inclination, argument of periapsis, true anomaly, and eccentric anomaly, respectively. We used the standard expressions from the twobody problem for radial distance, (7)
and angular momentum per unit mass, (8)
The integration of a trajectory is terminated in our simulation if the grain hits the surface of the Earth or the Moon or if it escapes from the EarthMoon system. We did not consider the effect of the Earth’s atmosphere. Collisions with the Earth or Moon could be easily missed due to the discrete time steps of the integrator. Therefore, cubic Hermite interpolation is used in our model to ensure the detection of collisions with the Earth or Moon if they occur between two consecutive time steps of the integrator (Liu et al. 2016).
Physical and orbital properties of celestial bodies.
Fig. 1 Solar radiation pressure efficiency, Q_{pr}, for spherical silicate particles calculated from Mie theory. 
3 Mass production rate and initial ejecta mass distribution
For the starting velocity we adopted a simple model in which the ejection direction is perpendicular to the lunar surface and all the particles have the same velocity distribution regardless of size. We assumed that the starting velocity follows a power law with exponent q. By normalizing , we obtain the form of the initial velocity distribution, (9)
where q = 2.2 was used, suitable for a regolithcovered surface (Krivov et al. 2003), v_{0} = 130 m s^{−1} is the minimum speed for lunar ejecta (Horányi et al. 2015), and v_{max} = 2v_{esc} ≈ 4.68 km s^{−1} is adopted for the maximum ejection speed.
Via normalization to the total mass production rate, , the distribution of mass can be expressed as (10)
where α = 0.9 was inferred from the LDEX data (Horányi et al. 2015), m_{min} is the minimum ejecta mass, and m_{max} is the maximum ejecta mass. The ejecta sizedistribution is obtained from the transformation (e.g., Liu & Schmidt 2018), (11)
where .
To apply these distributions of mass or radius, we must specify the total mass production rate, . For the first step, the mass production rate per unit surface is (12)
where F_{p} is the number flux of impactors with characteristic mass m_{p}, treated here as a parallel beam that hits the surface with an angle ϕ from the surface normal, and cosϕ is the projection area factor. The yield, Y_{p}, is defined as the ratio of the ejecta particles’ mass to the impactors’ mass, which is also a function of the material of the target surface as well as the mass and velocity of the impactors (Koschny & Grün 2001a).
Considering the effects of oblique impacts from experiments by Gault (1973) and the dependence on projectile mass and velocity, we have (13)
The mass production on the surface of the Moon varies with local time because the angle ϕ varies for projectile populations that approach the EarthMoon system from different directions. Six meteoroid populations, helion (HE), antihelion (AH), apex (AP), antiapex (AA), northern toroidal (NT), and southern toroidal (ST), are found to be plausible sources of highspeed impacts (Szalay et al. 2019). Following Szalay & Horányi (2015), and taking Eqs. (13) and (14) into account, the mass production rate per unit surface reads (14)
Here, C is a normalization constant, s labels the different projectile populations, F_{s}, m_{s}, and v_{s} are the number flux, characteristic mass, and impact velocity for each source, and Θ_{H} denotes the Heaviside function, which accounts for the fact that a given projectile population can only reach one hemisphere of the Moon. We followed Szalay et al. (2019) and used relative contributions w_{s} = (0.198, 0.198, 0.303, 0.025, 0.138, 0.138) for the HE, AH, AP, AA, NT, and ST sources of meteoroids. The symbol ϕ_{s} is the impact angle for each source from the surface normal, (15)
Here, θ is the colatitude of the surface element, θ_{s} is the colatitude for sources, θ_{NT} = 30°, θ_{ST} = 150° with φ_{NT} = φ_{ST} = 0° (Szalay et al. 2019), φ is defined as the longitude from the apex of the motion of the EarthMoon system for the surface element, φ_{s} is the sourcespecific longitude of approach in the ecliptic plane (HE, AH, AP, and AA), φ_{HE} = 65°, φ_{AH} = 295°, φ_{AP} = 0°, and φ_{AA} = 180° (Szalay & Horányi 2016). Therefore, Eq. (14) can be written as a function of the longitude, φ, and colatitude, θ, on the Moon. The total mass production rate, , for the Moon was estimated to be approximately 0.2 kg s^{−1} (Szalay et al. 2019; Pokornỳ et al. 2019). By normalization, , the normalization constant, C, in Eq. (14) can be fixed.
In our simulation, 200 starting positions on the surface of the Moon, including 20 longitudes and 10 colatitudes, are selected, which divides the whole surface of the Moon into 20 × 10 elements. Longitudes, φ, are equidistant from 0 to 2π, and colatitudes, θ, vary from 0 to π. Accordingly, the mass production rate for a given surface element (i, j) reads (16)
where i is the longitude index from 1 to 20, j is the colatitude index from 1 to 10, R_{M} is the radius of the Moon, and are boundary longitudes of surface element (i, j), and and are boundary colatitudes of surface element (i, j). The distribution of mass production rates is shown in Fig. 2.
Correspondingly, the size distribution for surface element (i, j) reads (17)
4 Simulation results
To cover the size distribution of ejecta, we selected nine dust particle radii: 0.2, 0.5, 1, 2, 5, 10, 20, 50, and 100 μm. Here, r_{min} = 0.2 μm and r_{max} = 100 μm, so and. For the initial velocity, ten velocities in the range [0.95 v_{esc}, 2 v_{esc}] were used, where v_{esc} is the escape velocity of the Moon. We note that the velocity distribution is normalized to unity in the range [v_{0}, v_{max} = 2v_{esc}], where v_{0} = 130 m s^{−1} is the minimal starting velocity of ejecta. Integrating over the velocity distribution in the interval [0.95v_{esc}, 2v_{esc}] then removes in our normalization the grains in [v_{0}, 0.95v_{esc}], which all fall back to the surface. A small fraction of grains in the range [0.95v_{esc}, v_{esc}] escape the Moon due to threebody effects. They are considered in our evaluation of the simulations.
To average over seasonal effects, we assumed that particles start at 12 different times, covering equidistantly one period of the ascending node precession of the Moon. We verified that this choice also covers the different phases of the Moon fairly uniformly. The number of grains necessary to cover the ranges of radius, velocity, time, colatitude, and longitude are 9, 10, 12, 10, and 20, respectively. Thus, the total number of particles is 216 000, which requires a huge amount of CPU time. The longterm simulations were carried out on the supercomputer located at the Finnish CSCIT Center for Science.
Fig. 2 Mass production rate per surface element on the lunar surface. Zero longitude is defined by the direction of the apex of the motion of the EarthMoon system. 
Fate of particles ejected from the lunar surface after 100 yr (fraction of 24 000 particles for each size).
4.1 Particle sinks and particle lifetimes
The final state of particles ejected from the surface of the Moon after 100 yr is shown in Table 2. Most dust particles (>88%) ultimately leave the EarthMoon system after 100 Earth years,regardless of size.
Only a very small fraction of grains (>2 μm) remain in orbit after 100 yr. We do not integrate these trajectories further, to save CPU time. The fraction of grains that rapidly reimpact the Moon within one orbital period is approximately 22%, which only weakly depends on grain size. We note that this fraction depends on our choice of v_{min} = 0.95v_{esc}, so the fraction of grains falling back to the Moon within one orbital period is not listed in Table 2. The grains that reimpact the Moon after a longer time are listed in Table 2 in the column “Hit Moon”. Small particles are more sensitive to solar radiation pressure. The fraction of grains that hit the Moon increases with grain size. The fraction that hit the Earth grows from 1.2 to 10% at a grain size of 2 μm and then begins to drop to a value of 3% for 100 μm grains.
Figure 3a shows the average life spans as a function of the size of the dust particles. Small particles, including those sized 0.2, 0.5, and 1 μm, have short lifetimes, less than 0.1 yr. For particles between 1 and 10 μm, the life span increases rapidly with grain size. Large particles with sizes bigger than 10 μm have longer lifetimes, approximately 0.7 Earth years. The parameter β denotes the ratio of solar radiation pressure to gravity (Burns et al. 1979): (18)
As shown in Fig. 3b, the effect of solar radiation pressure is large for small particles, which induces high eccentricities, such that small grains are rapidly expelled from the EarthMoon system and the life spans of those particles are short. With increasing grain size, β decreases rapidly and becomes negligible for grains of 10 μm or larger. Thus, the lifetimes of these particles are much longer.
Fig. 3 Average lifetimes and β for grains of different sizes. Panel a: average lifetimes. Panel b: β. 
4.2 Steadystate configuration
A cylindrical grid is utilized in the simulation process, where , ϕ_{c} = atan2(y, x) and z_{c} = z. Indexes i_{cell}, j_{cell}, and k_{cell} are used to label grid cells. Hence, the number density of a given grid cell (i_{cell}, j_{cell}, k_{cell}) reads (19)
For all trajectories, the grain positions are stored at equal time intervals, Δt. The symbol n(i_{cell}, j_{cell}, k_{cell}, r_{g}, v, i, j) is the plain number of dust particles with radius r_{g}, starting velocity v, and starting surface element (i, j) recorded in the cell (i_{cell}, j_{cell}, k_{cell}), V (i_{cell}, j_{cell}, k_{cell}) is the volume of the cell (i_{cell}, j_{cell}, k_{cell}), and N_{start}(r_{g}, v, i, j) is the number of grains with a given starting radius r_{g}, velocity v, and starting element (i, j).
In a similar manner we can obtain the geometric optical depth covered by the grains in a cell, by integrating over in Eq. (21) and dividing by, instead of volume, the surface of the cell (along a given line of sight). The normal optical depth and number density distribution are illustrated in Fig. 4. Figure 4a shows the normal optical depth of the dust in a Sunoriented frame, where the positive x axis always points in the direction of the Sun. Particles occupy a torus between the Earth and the Moon with an outer edge that is separated by roughly one lunar Hill radius from the orbit of the Moon. The torus is asymmetric and slightly offset toward the Sun. From Fig. 4b, we can see that dust particles are distributed nearly symmetrically about the equator of the Earth, covering a wide vertical range of up to tens of Earth radii. We note, however, that the configuration shown in Fig. 4b was evaluated in an inertial frame, averaging over grains ejected from the Moon at 12 times equidistantly over one period of the ascending node precession of the Moon, and over an evolution time of more than 100 yr. Thus, the effects of the inclination of the Sun and the Moon are averaged out. The instantaneous dust configuration is expected to differ somewhat from the one shown in Fig. 4b because the instantaneous Laplace plane of the dynamical problem does not coincidewith the equatorial plane of the Earth.
The differential size distribution in the torus between the Earth and the Moon is shown in Fig. 5. The size distribution in the torus after longterm simulations remains close to the initial ejecta distribution but has flattened in the size range of 1–10 μm. This can be attributed to the rapidly changing life span for grains in this size range (see Fig. 3a).
Fig. 4 Optical depth and number density of particles ejected from the surface of the Moon. Panel a: normal geometric optical depth. Panel b: azimuthaveraged number density in the ρ_{c} − z_{c} plane. The red circle is the Earth, and the red line denotes the orbit of the Moon. 
Fig. 5 Steadystate differential size distribution for the torus from simulations. The red line denotes the initial ejecta size distribution from normalized to unity, and the blue line denotes the steadystate size distribution in the torus from Fig. 4. 
4.3 Analysis of orbital elements
Figure 6 shows the distribution of eccentricity for 0.2 μm particles. The eccentricities are dispersed from 0 to 15, that is, a large number (~80%) of particles are hyperbolic, which explains their short life spans. These grains leave the EarthMoon system in a very shorttime and contribute little to the number density. The evolution of eccentricity for a typical 0.2 μm particle with a lifetime of 0.12 yr is shown in Fig. 7a. The initial eccentricity is near 0.5, that is, this particle starts from an elliptical orbit. Initially, the eccentricity fluctuates slightly within a few days and then grows continually, representing a typical evolution from an ellipse to a hyperbola, which is a consequenceof the strong effect of solar radiation pressure. Figure 7b depicts the contribution of the perturbations induced by lunar gravity, solar gravity, and radiation pressure on the growth rate d e∕dt, calculated from Eq. (5). Solar gravity has a negligible influence on de∕dt. Initially, thede∕dt due to solar radiation pressure is a small negative or positive value. After 0.08 yr, solar radiation pressure dominates and d e∕dt remains positive. As a result, the eccentricity increases monotonically, and finally, the orbit becomes hyperbolic.
Large particles remain bound to the EarthMoon system for a much longer time, and they make a significant contributionto the dust population. The distributions of the semimajor axis, eccentricity, and inclination of all 100 μm particles are shown in Fig. 8. The semimajor axis peaks around 20 R_{E}, mainly ranging from 10 R_{E} to 50 R_{E}. The eccentricity has a centered distribution below unity, peaking near 0.7. The inclination exhibits two peaks, at values of 20° and 140°, that is, a significant fraction of grains are on retrograde orbits.
We also calculated the distributions of the solar angle for 1 and 100 μm grains (Fig. 9). The solar angle is defined as the angle between the grain’s orbital pericenter and the Sun as seen from Earth (Hamilton 1993): (20)
Here, Ω is the longitude of the ascending node, ω is the argument of periapsis, and ψ_{⊙} is the solar longitude in the inertial frame.
From Fig. 9, the solar angle distribution is approximately symmetric about ϕ_{⊙} = 180°. The solar angle of 1 μm particles peaks at 180°, that is, the population of these grains is offset toward the Sun. For 100 μm particles, the curve is smoother but shows two peaks, at 0° and 180°, which indicates that the offset of the torus formed by 100 μm particles is much smaller compared to that of the 1 μm grains. For particles of different sizes, the distribution exhibits different kinds of asymmetry, with offsets toward or away from the Sun. However, averaging over grain size, the torus is slightly offset toward the Sun (see Fig. 4a).
Figures 10 and 11 show the evolution of radial distance and inclination for a single 100 μm particle. This is one of the rare particles in our simulations that still remained in orbit after 100 yr. The grain covers a range between 10 R_{E} and 50 R_{E}, which is consistent with the distribution in Fig. 4a. The inclination oscillates from 0° to 60°. This particle is always moving on a prograde orbit within its life span. Such particles make a large contribution to the distribution of inclination (see Fig. 8c).
A considerable part of the particles in our simulations are on a retrograde orbit. The fractions of retrograde particles for different sizes are shown in Fig. 12. Some grains are ejected from the Moon directly into retrograde orbits (Colombo et al. 1966), but a larger number of grains evolve from a prograde to a retrograde orbit in the longterm simulation. To show this, a 1 μm particle is used here as an example. The evolution of inclination of this particle and contributions of specific forces to d i∕dt calculated from Eq. (6) are shown in Figs. 13 and 14. From Fig. 13, we see that the initial inclination of the grain is near 60°, and the final one is 170°. This particle evolves from a prograde orbit to a retrograde orbit multiple times. As can be seen in Fig. 14, the effects of gravity of the Moon and Sun on the evolution of d i∕dt are much smaller than the one induced by the solar radiation pressure, although there are some peaks from encounters with the Moon. The approximate change in inclination can be evaluated via the integration of the contribution to d i∕dt due to solar radiation pressure, which is about 100°.
Fig. 6 Distribution of eccentricity for 0.2 μm particles. 
Fig. 7 Evolution of eccentricity and de∕dt for one 0.2 μm particle. Panel a: evolution of eccentricity. Panel b: evolution of de∕dt. 
Fig. 8 Distribution of orbital elements for 100 μm particles ejected from the surface of the Moon. Panel a: distribution of the semimajor axis. Panel b: distribution of eccentricity. Panel c: distribution of inclination. 
Fig. 9 Distribution of solar angle for 1 and 100 μm particles ejected from the surface of the Moon. 
Fig. 10 Evolution of r for a 100 μm particle. 
Fig. 11 Evolution of inclination for a 100 μm particle. 
Fig. 12 Fraction of retrograde particles. 
Fig. 13 Evolution of inclination for a 1 μm particle. 
5 Conclusions
In this work we have derived the steadystate configuration of particles in the EarthMoon system that have been ejected from the lunar surface in hypervelocity impacts of micrometeoroids (Horányi et al. 2015). A variety of forces, including the solar radiation pressure and the gravity of the Sun, the Earth, and the Moon, are considered in a numerical exploration of the trajectories of the ejected particles in the system. Expanding on previous work in the literature (Szalay & Horányi 2015, 2016; Szalay et al. 2019), we calculated the initial ejecta distribution and mass production rate on the lunar surface, which we used as initial conditions for the integrations. The final states (sinks), the average life spans, and the fraction of retrograde grains have been derived as functions of grain size. We also show examples for the evolution and distribution of orbital elements of particles. Small particles tend to rapidly evolve into hyperbolic orbits (on timescales of weeks), while most large particles remain in elliptical orbits for a longer time (up to a year). Manyof the grains develop high inclinations, and we find that a substantial fraction evolve into retrograde orbits.
From our longterm integrations we find that about 3.6 × 10^{−3} kg s^{−1} (1.8%) particles escape from the lunar gravity. These particles form a tenuous broad torus, the densest part of which roughly spans a distance of 40 R_{E} from the Earth. We also estimated the mass rate of particles coming into the Earth’s atmosphere, , which is a nonnegligible value. As the space activities in the EarthMoon system are much more frequent compared with other planets, dust around the Earth seriously threatens not only the mechanical structure of spacecraft but potentially also the health of astronauts. In Apollo missions, astronauts were exposed to the lunar dust environment, and respiratory, dermal, and ocular irritations from lunar dust were reported (Turci et al. 2015). NASA established the Lunar Airborne Dust Toxicity Advisory Group (LADTAG) to assess the risk degree for spacecraft and astronauts in lunar missions (KhanMayberry 2008). Thus, analysis of the characters for particles in the torus around the Earth is of great significance for assessing the space environment and ensuring thesecurity of explorations.
We characterize the torus in terms of the normal geometric optical depth and number density. Owing to the effect of solar radiation pressure, the torus is mildly offset toward the solar direction. The peak normal optical depth we obtain is on the order of 10^{−12}, which is small compared to the optical depth of the Thebe extension of Jupiter’s gossamer rings, ~ 10^{−9} (De Pater et al. 2018), and the limits inferred for the yet undetected Phobos ring of Mars, ~ 10^{−8} (Showalter et al. 2006). The peak number density we find for the lunar torus is on the order of 10^{−9} m^{−3}, which is about four orders of magnitude smaller than the density of the Deimos ring inferred for radii of 15 μm by Juhász & Horányi (1995) and four and five orders of magnitude smaller than the number for the Deimos and Phobos rings (≥ 0.5 μm) determined by Liu & Schmidt (2021), respectively. For an in situ detector on a circular orbit at a distance of 10 R_{E} pointing in the prograde orbit direction, this translates into a flux of about 50 grains per m^{2} and per year (taking into account that a fraction of the grains are on retrograde orbits). From the Interplanetary Meteoroid Engineering Model (IMEM; Dikarev et al. 2005). we find that this is a factor of several smaller than the flux of micronsized interplanetary micrometeoroids at 1 AU. The directional distributions of the lunar grains and the interplanetary particles, however, would be different.
Our results on the grain lifetimes, their sinks, and their typical orbital evolution are robust because the relevant perturbation forces are taken into account and the respective parameters (masses of the gravitational perturbers, the solar constant, and the solar radiation pressure efficiency factor) are known up to small uncertainties. Large uncertainties, however, remain in our calibration of number density and optical depth. On one hand, our simulation results are based on the mass production rate of reported by Szalay et al. (2019). If we follow the method used in the literature for an analysis of ejecta from the Galilean moons (Krivov et al. 2003; Sremčević et al. 2003), employ the mass flux at 1AU from the IMEM model (Dikarev et al. 2005), and use the empirical expression for the yield suggested by Koschny & Grün (2001a), we obtain a mass production rate that is about one order of magnitude higher. Although in our analysis we prefer to use the value from the direct measurement, this suggests that the uncertainty in the calibration of our model is at least one order of magnitude. On the other hand, we assume a power law distribution for the ejection velocity, whose high velocity tail determines the fraction of ejecta that escape the Moon in our model. In the literature, a different functional form was fitted successfully to the bound particles detected by LDEX (Szalay & Horányi 2016). If the tail of the velocity distribution of the lunar ejecta turns out to be different from the one assumed in our paper, then the dust densities derived from our model will also change.
Fig. 14 Variation in di∕dt due to various forces for a 1 μm particle. Panel a: variation in di∕dt due to lunar gravity. Panel b: variation in di∕dt due to solar gravity. Panel c: variation in di∕dt due to solar radiation pressure. 
Acknowledgements
This work was supported by National Natural Science Foundation of China (No. 12002397) and by the China Scholarship Council (CSC, 201906220134). We acknowledge CSC – IT Center for Science for the allocation of computational resources.
References
 Burns, J. A., Lamy, P. L., & Soter, S. 1979, Icarus, 40, 1 [Google Scholar]
 Colombo, G., Shapiro, I. I., & Lautman, D. A. 1966, J. Geophys. Res., 71, 5719 [NASA ADS] [CrossRef] [Google Scholar]
 De Pater, I., Hamilton, D. P., Showalter, M. R., Throop, H. B., & Burns, J. A. 2018, Planetary Ring Systems, The Rings of Jupiter, eds. M. S. Tiscareno, & C. D. Murray (Cambridge: Cambridge University Press), 125 [CrossRef] [Google Scholar]
 Dikarev, V., Grün, E., Baggaley, J., et al. 2005, Adv. Space Res., 35, 1282 [NASA ADS] [CrossRef] [Google Scholar]
 Gault, D. E. 1973, The Moon, 6, 32 [NASA ADS] [CrossRef] [Google Scholar]
 Hamilton, D. P. 1993, Icarus, 101, 244 [NASA ADS] [CrossRef] [Google Scholar]
 Horányi, M., Szalay, J., Kempf, S., et al. 2015, Nature, 522, 324 [CrossRef] [Google Scholar]
 Iglseder, H., Uesugi, K., & Svedhem, H. 1996, Adv. Space Res., 17, 177 [NASA ADS] [CrossRef] [Google Scholar]
 Juhász, A., & Horányi, M. 1995, J. Geophys. Res. Planets, 100, 3277 [CrossRef] [Google Scholar]
 KhanMayberry, N. 2008, Acta Astron., 63, 1006 [NASA ADS] [CrossRef] [Google Scholar]
 Koschny, D., & Grün, E. 2001a, Icarus, 154, 391 [NASA ADS] [CrossRef] [Google Scholar]
 Krivov, A. V., Sremčević, M., Spahn, F., Dikarev, V. V., & Kholshevnikov, K. V. 2003, Planet. Space Sci., 51, 251 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, X., & Schmidt, J. 2018, A&A, 609, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Liu, X., & Schmidt, J. 2019, Astrodynamics, 3, 17 [CrossRef] [Google Scholar]
 Liu, X., & Schmidt, J. 2021, MNRAS, 500, 2979 [Google Scholar]
 Liu, X., Sachse, M., Spahn, F., & Schmidt, J. 2016, J. Geophys. Res. Planets, 121, 1141 [NASA ADS] [CrossRef] [Google Scholar]
 Mukai, T. 1989, Evolution of Interstellar Dust and Related Topics, Proceedings of the International School of Physics (North Holland: Elsevier Science) [Google Scholar]
 Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge: Cambridge University Press) [Google Scholar]
 Pavlis, N. K., Holmes, S. A., Kenyon, S. C., & Factor, J. K. 2012, J. Geophys. Res. Solid Earth, 117, B04406 [Google Scholar]
 Peale, S. J. 1966, J. Geophys. Res., 71, 911 [NASA ADS] [CrossRef] [Google Scholar]
 Pokornỳ, P., Janches, D., Sarantos, M., et al. 2019, J. Geophys. Res. Planets, 124, 752 [CrossRef] [Google Scholar]
 Showalter, M. R., Hamilton, D. P., & Nicholson, P. D. 2006, Planet. Space Sci., 54, 844 [NASA ADS] [CrossRef] [Google Scholar]
 Solomon, S. C. 1974, The Moon, 9, 147 [NASA ADS] [CrossRef] [Google Scholar]
 Spahn, F., Sachse, M., Seiß, M., et al. 2019, Space Sci. Rev., 215, 1 [CrossRef] [Google Scholar]
 Sremčević, M., Krivov, A. V., & Spahn, F. 2003, Planet. Space Sci., 51, 455 [CrossRef] [Google Scholar]
 Szalay, J. R., & Horányi, M. 2015, Geophys. Res. Lett., 42, 5141 [NASA ADS] [CrossRef] [Google Scholar]
 Szalay, J. R., & Horányi, M. 2016, Geophys. Res. Lett., 43, 4893 [NASA ADS] [CrossRef] [Google Scholar]
 Szalay, J. R., Pokornỳ, P., Sternovsky, Z., et al. 2019, J. Geophys. Res. Planets, 124, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Tiscareno, M. S., & Murray, C. D. 2018, Planetary Ring Systems: Properties, Structure, and Evolution (Cambridge: Cambridge University Press) [CrossRef] [Google Scholar]
 Turci, F., Corazzari, I., Alberto, G., Martra, G., & Fubini, B. 2015, Astrobiology, 15, 371 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Q., & Liu, J. 2016, Acta Astron., 127, 678 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Fate of particles ejected from the lunar surface after 100 yr (fraction of 24 000 particles for each size).
All Figures
Fig. 1 Solar radiation pressure efficiency, Q_{pr}, for spherical silicate particles calculated from Mie theory. 

In the text 
Fig. 2 Mass production rate per surface element on the lunar surface. Zero longitude is defined by the direction of the apex of the motion of the EarthMoon system. 

In the text 
Fig. 3 Average lifetimes and β for grains of different sizes. Panel a: average lifetimes. Panel b: β. 

In the text 
Fig. 4 Optical depth and number density of particles ejected from the surface of the Moon. Panel a: normal geometric optical depth. Panel b: azimuthaveraged number density in the ρ_{c} − z_{c} plane. The red circle is the Earth, and the red line denotes the orbit of the Moon. 

In the text 
Fig. 5 Steadystate differential size distribution for the torus from simulations. The red line denotes the initial ejecta size distribution from normalized to unity, and the blue line denotes the steadystate size distribution in the torus from Fig. 4. 

In the text 
Fig. 6 Distribution of eccentricity for 0.2 μm particles. 

In the text 
Fig. 7 Evolution of eccentricity and de∕dt for one 0.2 μm particle. Panel a: evolution of eccentricity. Panel b: evolution of de∕dt. 

In the text 
Fig. 8 Distribution of orbital elements for 100 μm particles ejected from the surface of the Moon. Panel a: distribution of the semimajor axis. Panel b: distribution of eccentricity. Panel c: distribution of inclination. 

In the text 
Fig. 9 Distribution of solar angle for 1 and 100 μm particles ejected from the surface of the Moon. 

In the text 
Fig. 10 Evolution of r for a 100 μm particle. 

In the text 
Fig. 11 Evolution of inclination for a 100 μm particle. 

In the text 
Fig. 12 Fraction of retrograde particles. 

In the text 
Fig. 13 Evolution of inclination for a 1 μm particle. 

In the text 
Fig. 14 Variation in di∕dt due to various forces for a 1 μm particle. Panel a: variation in di∕dt due to lunar gravity. Panel b: variation in di∕dt due to solar gravity. Panel c: variation in di∕dt due to solar radiation pressure. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.