Simulation of Dynamic Evolution of Ring Current Ion Flux by a Lunar Base Energetic Neutral Atom (ENA) Imaging

: The distribution of energetic ion ﬂux in the ring current region, such as a meteorological cumulonimbus cloud, stores up the particle energy for a geomagnetic substorm. It is helpful to study the geomagnetic substorm mechanism by using a lunar base ENA imaging simulation of the dynamic evolution of the ring current, and establishing the corresponding relationship between key node events of the substorm. Based on the previous observation experience and our simulation results of the dynamic evolution of the ring current, we propose a macroscopic model of substorms related to the dynamic evolution of ring currents and present the possibility of conﬁrming the causal sequence of some of those critical node events of substorms with the lunar base ENA imaging measurement. IBEX, operating in the ecliptic plane, may even give examples of the telemetry of ring current ion ﬂuxes through ENA measurements during substorms/quiets.


Introduction
The ring current ion flux spreads from 2 to 8 Earth radii (R E ), with the maximum at 4 Re, where ion energy ranges from tens of eV to tens of keV, accompanied by westward drift.The enhancement of the ring current is due to injection from the magnetotail plasma sheet during a geomagnetic storm.Its effect is to generate a magnetic field at the ground that is opposite to the geomagnetic field, causing a negative disturbance in the Dst index.The partial disturbance of the ring current produces substorms, which often affect the highlatitude auroral region.The plasma pressure and ion fluxes of the ring current increased dramatically during super substorms, while little change was observed for substorms with AE < 1000 nT [1].During substorms, the major terrestrial energetic neutral atom (ENA) emissions also originate from the ring current.
Ogasawara et al. [2] reported a continuous daylong observation of bright terrestrial ENA emissions with the Interstellar Boundary Explorer (IBEX) [3].The IBEX-Hi operates in a lunar resonance orbit with an apogee of ~48 Earth radius (R E ), and it can monitor ENAs from 0.3 to 6.0 keV for long periods of time.The correlation between the auroral electric jet (AE) index and ENA emissions was found through long-term monitoring.The results show that the temporal profile of ENA correlated with Earth's auroral activities.However, IBEX-Hi's single-pixel (6 • × 6 • ) scan observations make it difficult to distinguish the ENA emission source from the low-altitude auroral region or the ring current region, providing no further information about the trigger of the substorm.
There are many substorm models, and they all involve three key node phenomena, including the magnetotail lobe reconnection, cross-tail current disruption, and aurora brightening (see Figure 1).The widely accepted models are the near-Earth current disruption model [4][5][6][7] and the near-Earth neutral line model [8,9]; both include the three key phenomena mentioned above, but occur in different timing sequence.
There are many substorm models, and they all involve three key node phe including the magnetotail lobe reconnection, cross-tail current disruption, an brightening (see Figure 1).The widely accepted models are the near-Earth curren tion model [4][5][6][7] and the near-Earth neutral line model [8,9]; both include the phenomena mentioned above, but occur in different timing sequence.Among the three key node phenomena of the substorm, only auroral br (including the brightening of visible light, EUV, and ENA emission) is visible, w the spatiotemporal telemetry evidence of its occurrence.The other two items ha ibility and occur randomly at specific times and in limited airspace.There was n direct-measurement evidence to confirm their temporal and spatial existence.IBE served ENA scanning images of the magnetotail plasma sheet disconnection during the geomagnetic quiet period at the apogee of the lunar resonance orbit [1 brings a glimmer of light to the visual remote sensing of the ring current.
In Section 2, we discuss ENA imaging simulation of ring current dynamic e Section 3 describes a macro substorm model, which can be dynamically mon ENA imaging.In Section 4, the supporting evidence of previous observations fo posed model is given.Finally, the Section 5 gives a summary.

Lunar Orbit and Instrument Parameters for the Simulation
Average lunar orbit distance of 380,000 km, about 60 radii of the Earth (Re), respect to the Moon's orbit, has an angle of less than 1°.The space range of the spheres' ring current is between 2 and 8 RE.Therefore, the field of view of the in is designed to be 20° × 20°, and the angle resolution is better than 0.5°.It is p integrate 40 one-dimensional arrays, each consisting of 40 detectors, into a tw sional detector array.According to the lunar orbit positions (see Figure 2), th technical specifications of the simulation design of the lunar base ENA image FOV and angular resolution, are shown in Table 1.ENA generated by ring curr getic ions is the strongest omnidirectional emission source for particle imaging tion of the geo-magnetosphere in lunar orbit.The geometric factor of a single d the instrument is 0.0044 cm 2 sr, and the integration time is about 3 min, which the statistical requirements of ENA imaging [11].The ENA is suitable for the tra cle of ring current energetic ion distribution.The intensity and distribution of in different orbital positions are predicted by simulation, which is similar to the telemetry instruments, and the results are expected.Among the three key node phenomena of the substorm, only auroral brightening (including the brightening of visible light, EUV, and ENA emission) is visible, which has the spatiotemporal telemetry evidence of its occurrence.The other two items have no visibility and occur randomly at specific times and in limited airspace.There was never any direct-measurement evidence to confirm their temporal and spatial existence.IBEX-Hi observed ENA scanning images of the magnetotail plasma sheet disconnection (at 10 R E ) during the geomagnetic quiet period at the apogee of the lunar resonance orbit [10], which brings a glimmer of light to the visual remote sensing of the ring current.
In Section 2, we discuss ENA imaging simulation of ring current dynamic evolution.Section 3 describes a macro substorm model, which can be dynamically monitored by ENA imaging.In Section 4, the supporting evidence of previous observations for the proposed model is given.Finally, the Section 5 gives a summary.

Lunar Orbit and Instrument Parameters for the Simulation
Average lunar orbit distance of 380,000 km, about 60 radii of the Earth (R e ), 1 R e with respect to the Moon's orbit, has an angle of less than 1 • .The space range of the magnetospheres' ring current is between 2 and 8 R E .Therefore, the field of view of the instrument is designed to be 20 • × 20 • , and the angle resolution is better than 0.5 • .It is planned to integrate 40 one-dimensional arrays, each consisting of 40 detectors, into a twodimensional detector array.According to the lunar orbit positions (see Figure 2), the specific technical specifications of the simulation design of the lunar base ENA imager, include FOV and angular resolution, are shown in Table 1.ENA generated by ring current energetic ions is the strongest omnidirectional emission source for particle imaging exploration of the geo-magnetosphere in lunar orbit.The geometric factor of a single detector of the instrument is 0.0044 cm 2 sr, and the integration time is about 3 min, which can meet the statistical requirements of ENA imaging [11].The ENA is suitable for the tracer particle of ring current energetic ion distribution.The intensity and distribution of ENA flux in different orbital positions are predicted by simulation, which is similar to the previous telemetry instruments, and the results are expected.

Simulation Equations
The counts recorded in each pixel of an ENA image, C(δ, ε) with elevation δ and azimuth ε, are represented in the simulation equation, where ΔE is the energy range within, ΔT the integral time for the pixel, and ΔΩ the solid angle of the volume element pointing to the δ, ε pixel; jion represents ion differential flux at the integral volume element, A(δ, ε) the response function of a detector [12], and σ(E) the charge exchange cross sections [13]; n(r, φ, θ) is the exospheric neutral atomic density, where r represents the geocentric distance, φ the longitude, and θ the latitude; and dV is the volume element integral along the line of sight of the detector.The ion flux in the ring current region may be expressed in the form of [14] 1 Table 1.Instrument simulation design parameters [11].

Simulation Equations
The counts recorded in each pixel of an ENA image, C(δ, ε) with elevation δ and azimuth ε, are represented in the simulation equation, where ∆E is the energy range within, ∆T the integral time for the pixel, and ∆Ω the solid angle of the volume element pointing to the δ, ε pixel; j ion represents ion differential flux at the integral volume element, A(δ, ε) the response function of a detector [12], and σ(E) the charge exchange cross sections [13]; n(r, ϕ, θ) is the exospheric neutral atomic density, where r represents the geocentric distance, ϕ the longitude, and θ the latitude; and dV is the volume element integral along the line of sight of the detector.The ion flux in the ring current region may be expressed in the form of [14] where e = (1 + 1/κ) κ+1 ≈ 2.962 (k = 5.5), θ is the magnetic latitude, and where where ϕ s = 180 • , ξ = 0.73 for azimuth asymmetry of ring current ion flux.L is the magnetic shell index, where where Kp is the geomagnetic activity index, and the maximum ion fluxes and corresponding positions by the model above are shown in Table 2. Equatorial ion flux distributions under different Kp indices are shown in Figure 3. where where φs = 180°, ξ = 0.73 for azimuth asymmetry of ring current ion flux.L is the magnetic shell index, where ), ( ) ( , and where Kp is the geomagnetic activity index, and the maximum ion fluxes and corresponding positions by the model above are shown in Table 2. Equatorial ion flux distributions under different Kp indices are shown in Figure 3.The pitch angle function (fα) of the particles within the volume element is expressed in terms of αeq (equatorial pitch angle) such that The pitch angle function (f α ) of the particles within the volume element is expressed in terms of α eq (equatorial pitch angle) such that Here, the parameter K α (adjusting the pitch angle distribution width) corresponds to the Kp index from 0 to 7, which simulates the dynamic evolution of pitch angle from pancake to butterfly during a substorm [15][16][17][18], as shown in Figure 4. .
Here, the parameter Kα (adjusting the pitch angle distribution width) correspo the Kp index from 0 to 7, which simulates the dynamic evolution of pitch angle from cake to butterfly during a substorm [15][16][17][18], as shown in Figure 4.The minimum energy channel of the ENA imager starts at 4 keV, so it does rectly record the neutral gas (≤10 6 K) in the background.The ENA map recorded ENA imager is the energetic neutral atom produced after the charge exchange of th tral atom evaporated from the exosphere and energetic ions precipitated from th current region.
The neutral hydrogen densities provided by Chamberlain's model, adjust Rairden et.al. [19], are too low to be applicable in the ring current region.If conser of the radial flux of neutral atoms during a major geomagnetic storm is assumed, w however, add an additional term, (a0/r) 2 , to the Chamberlain model so that the dens be represented by , where r is the geocentric distance in RE, a0 = 1.78 RE represents the height of the exos and n0 = 1600 cm −3 represents the neutral density constant.This model is better fitt the yearly average hydrogen density in solar maximum conditions in the ring curr gion measured by Tinsley et al. [20], which should also be suitable in the substorm p

Simulation Results
We simulate the dynamic process of magnetospheric substorms with the incre the Kp index combined with ion pitch angle diffusion.The low-altitude polar EN hancement caused by substorms often appears before midnight and is accompan eastward drift.When the ENA-enhanced signal decays, it will be accompanied by ward movement [21,22].Thus, we mimic the substorm growth phase with the Kp ascending sequence of the simulation images at the full moon and the first quarte ures 5 and 6).The substorm recovery phase is simulated with the Kp index desce sequence of the simulation image in the last quarter (Figure 7).The minimum energy channel of the ENA imager starts at 4 keV, so it does not directly record the neutral gas (≤10 6 K) in the background.The ENA map recorded by the ENA imager is the energetic neutral atom produced after the charge exchange of the neutral atom evaporated from the exosphere and energetic ions precipitated from the ring current region.
The neutral hydrogen densities provided by Chamberlain's model, adjusted by Rairden et.al. [19], are too low to be applicable in the ring current region.If conservation of the radial flux of neutral atoms during a major geomagnetic storm is assumed, we may, however, add an additional term, (a 0 /r) 2 , to the Chamberlain model so that the density can be represented by where r is the geocentric distance in R E , a 0 = 1.78 R E represents the height of the exosphere, and n 0 = 1600 cm −3 represents the neutral density constant.This model is better fitted for the yearly average hydrogen density in solar maximum conditions in the ring current region measured by Tinsley et al. [20], which should also be suitable in the substorm period.

Simulation Results
We simulate the dynamic process of magnetospheric substorms with the increase in the Kp index combined with ion pitch angle diffusion.The low-altitude polar ENA enhancement caused by substorms often appears before midnight and is accompanied by eastward drift.When the ENA-enhanced signal decays, it will be accompanied by westward movement [21,22].Thus, we mimic the substorm growth phase with the Kp index ascending sequence of the simulation images at the full moon and the first quarter (Figures 5 and 6).The substorm recovery phase is simulated with the Kp index descending sequence of the simulation image in the last quarter (Figure 7).During geomagnetic quiet, Figures 5a, 6a and 7d, the pitch angle of ring current e ergetic ions is mostly 90° (pancake-shaped distribution), and the ENA emission source located in the equatorial ring current region and only propagates within the ecliptic plan with slow attenuation.At ~60 RE lunar orbit, the ENA imager can obtain about 200 coun in a 3 min integration time.Among them, the maximum ENA event count of a single pix is about 4-12, as shown in Table 3.The first quarter image of the ENAs in Figure 6a a distributed near L ≥ 7 and tend to expand along magnetic field lines.The equatorial rin current is far away from the Earth and the neutral gas density is low, resulting in a lo global energetic neutral index (GENI) by Ogasawara et al. [2].Almost no ENA emissio was measured in the auroral zone.With the increase in Kp index and the pitch angle dispersion, the ENA emissio source begins to shift along the magnetic field line to the low-altitude and high-latitud polar region, and the magnetic flux tubes converge in the auroral zone at low altitude.A the density of neutral gas is much higher than that in the equinoctial ring current, EN emission from the auroral zone is rapidly enhanced, and its GENI also increases corr spondingly, as shown in Figures 5b and 6b.After the dispersion of the pitch angle, th ring current energetic ions with continuous pitch angle distribution in the auroral zon During geomagnetic quiet, Figures 5a, 6a and 7d, the pitch angle of ring current energetic ions is mostly 90 • (pancake-shaped distribution), and the ENA emission source is located in the equatorial ring current region and only propagates within the ecliptic plane with slow attenuation.At ~60 R E lunar orbit, the ENA imager can obtain about 200 counts in a 3 min integration time.Among them, the maximum ENA event count of a single pixel is about 4-12, as shown in Table 3.The first quarter image of the ENAs in Figure 6a are distributed near L ≥ 7 and tend to expand along magnetic field lines.The equatorial ring current is far away from the Earth and the neutral gas density is low, resulting in a low global energetic neutral index (GENI) by Ogasawara et al. [2].Almost no ENA emission was measured in the auroral zone.With the increase in Kp index and the pitch angle dispersion, the ENA emission source begins to shift along the magnetic field line to the low-altitude and high-latitude polar region, and the magnetic flux tubes converge in the auroral zone at low altitude.As the density of neutral gas is much higher than that in the equinoctial ring current, ENA emission from the auroral zone is rapidly enhanced, and its GENI also increases correspondingly, as shown in Figures 5b and 6b.After the dispersion of the pitch angle, the ring current energetic ions with continuous pitch angle distribution in the auroral zone can produce almost omnidirectional ENA emissions.Here, we define the ENA emission brightening in the aurora zone as ENA aurora.
When the Kp index continues to increase and the pitch angle diffuses further, the energetic ion flux perpendicular to the magnetic field decreases, and the ENA emission source in the equatorial ring current begins to decrease, but the GENI further increases, as shown in Figures 5c and 6c.
When the pitch angle finally diffuses into a butterfly shape distribution (the energetic ion flux perpendicular to the magnetic field line decrease to a minimum), the ENA emission source basically transfers to the high-latitude aurora zone, and the GENI also reaches a maximum, as shown in Figures 5d and 6d.ENA simulation counts may increase rapidly with the increase in ion pitch angle diffusion, and both the ENA maximum and the sum counts (or GENI) increased by more than one order of magnitude, as shown in Table 3.
In Figure 7, the Kp index evolves in a reversed order, with the pitch distribution returning from butterfly to pancake.Lunar base ENA imaging measurements of the entire substorm process were not only characterized by GENI [2], but also visually monitored the evolution of the ring current ion flux through the distribution pattern of ENA emission sources.

Macroscopic Models of Magnetospheric Substorms
Based on possibility of the visual remote sensing of ring current ion flux by ENA imaging, we propose a macroscopic model of magnetospheric substorms dominated by solar wind dynamic pressure.We believe that the energetic ions that support geomagnetic activity are originally stored in the ring current region, as shown in Figures 5a and 6a.When the solar wind particle flow (bow shock) strikes, in the growth phase, the magnetopause dynamic pressure increases and the geomagnetic field is stretched tailward (bow shock attack marked in black, Figure 1).Due to the conservation of the first adiabatic invariant, the pitch angle of energetic ions diffuses within the ring current (current disruption marked in red, Figure 1), and then those ions precipitate along the field lines towards the low-altitude polar region (current wedge marked in green, Figure 1), thus completing the triggering of the aurora substorm (aurora marked in red, Figure 1), as shown from Figures 5b-d and 6b-d.At the same time, the geomagnetic field permeating the magnetotail solar wind particles stream constantly picks up energetic ions in the solar wind.During the recovery phase of the substorm, the solar wind particle flow passed away, the solar wind dynamic pressure decreases, and the magnetic field begins dipolarization (field depolarization marked in black, Figure 1).The magnetic field carries the energetic ions collected from the solar wind into the ring current (ion flow injection marked in green, Figure 1).This allows the energetic ions lost in the loss cone during the precipitation to be replenished, as shown from Figure 7a-d.Visual remote sensing of magnetic reconnection at the magnetotail lobe is not available at present.However, either the high-speed stream caused by magnetotail reconnection (reconnection marked in red, Figure 1) or Earthward convection caused by magnetic field dipolarization can contribute to ring current ion injection.
Except for the phenomenon of magnetotail lobe reconnection, the different ring current morphological changes corresponding to the whole substorm process can be monitored by ENA imaging remote sensing.In general, under solar wind pressure, the angular diffusion process of energy ion ejection is very fast in the growth stage, and the recovery of substorms driven by magnetic field dipole is slow.However, the high-speed flow of energetic ions generated by magnetotail reconnection will accelerate the process of the substorm recovery stage, and even trigger a continuous series of substorms.
According to the previous inversion logic of ENA imaging, all the mass spectra and energy spectra of ENA are usually converted into the corresponding flux distributions of energetic ions on the equatorial plane.By measuring the distribution of energetic ion flux (including mass spectrum and energy spectrum) in the ring current before the substorm, changes in energy and mass spectrum (such as changes in oxygen ion composition) during the substorm can be monitored so as to obtain targeted information about the source and acceleration of ions of different components.

Observation Basis of Macroscopic Model of Substorms
The previous mission for ENA observations, such as Astrid [23], IMAGE [24], TC-2 [25], TWINS [26], etc., have high orbital inclination.With the exception of IBEX operating on the ecliptic plane, there are almost no ENA telemetry signals directly from the ring current.The ENA brightening signal basically comes from the low-altitude polar regions, which is called ENA aurora here.However, the NUADU/TC-2 measurements and inversions of the ENA auroral evolutionary series, combined with the in situ measurements of other multi-satellite systems, reveal many characteristics about the evolution of the ring current [21,22].IBEX, operating within the ecliptic plane, even gave examples of ENA remote measurements of ring current ion fluxes during substorms/quiets [10].

Multi-Satellite Joint Observation of Causal Sequence of Substorm Events
High time resolution ENA imaging data recorded by the NUADU instrument aboard TC-2 are available throughout the major magnetic storm of 15 May 2005 [21].Ion flux distributions inversed from those ENA imaging measurements showed that the evolution of the ring current ion flux during a major storm was closely related to variations of two series of substorms.Comparisons of these retrieved distributions with in situ measurements of ion fluxes made aboard the spacecraft of the LANL constellation and with magnetic field measurements recorded aboard the spacecraft of the GOES constellation (both of which are located in the equatorial plane at ~6.6 R E in the environment of the outer radiation belt) showed covariations in the three data sets.The magnetic field response of the substorms occurred earlier at the geosynchronous altitude on the night side than was the case on the ground [21].It was found that the ENA aurora brightening occurred during geomagnetic field stretching in the growth phases of the substorms, rather than after Earthward-directed dipolarization events (Tables S1 and S2).This suggests that the energetic ions that cause the ENA aurora brightening were stored in the ring current region before.This work does not support the neutral line model of magnetotail reconnection triggering substorms.The ion flow injection/convection events that occur after this are simply the particle energy for the subsequent substorm.

Azimuth Evolution Direction of ENA Aurora
On 12 November 2004, TC-2/NUADU obtained the ENA image with the highest temporal resolution by using the close measurement at the perigee, and then reversed the ion flux distribution of the equatorial ring current [22].We found that in ENA images projected on the equatorial plane (i.e., inverted ring current ion flux), the brightening of the ENA auroras was usually accompanied by an eastward shift of the enhanced area of the ring current ion flux.On the contrary, the ENA aurora moved westward when it weakened [21,22].Due to the rotation of the Earth, the response sequence of the magnetospheric shell stretching tailward is from west to east, while the response sequence of magnetic field dipolarization is the opposite, stretching from east to west (Figure S1).However, the above variation in the geomagnetic field has regional characteristics, so it is hardly suitable to obtain in situ measurement to support specific geomagnetic field variations.Anyway, these observational results are indirect evidence that that magnetospheric magnetic disturbances trigger substorms.These IBEX-Hi ENA scanning images with gaps appearing at about −10 R E were originally thought to be remote sensing images of the plasma sheet disconnection events (Figures 2 and 3 of Mccomas et al. [10], Figures S2 and S4).However, the differential fluxes of ENA on both sides of the crack are similar, indicating that the emission mechanism of ENA on both sides is the same.During the geomagnetic quiet, the energetic ions in the ring current are confined by the magnetic field, and their pitch angle is about 90 • .The ENA produced by the exchange of charge between them and the surrounding neutral gas only propagates in the ecliptic plane and decays slowly.On the other hand, the energetic ions in the plasma sheet are almost in a state of thermal motion under the condition of high β value.The ENA produced by them has the characteristics of omnidirectional propagation and rapid attenuation.In Figure 4a [10].A substorm occurred on 28 October 2009 (Figure S3).At 3:45, the AE index (AE = 22 nT) began to grow, reaching a maximum at 11:15 (AE = 281 nT), and quieted down at 23:05 (AE = 44 nT) for about a whole day.During the substorm, the pitch angle of energetic ions in the ring current diffused, and the ENAs propagated in the ecliptic plane decreased, resulting in the observation illusion of plasma disconnection.Before midnight on the 28th, the geomagnetic activity returned to quiet, and the scanning FOV of IBEX-Hi could still cover the ring current at the far end of the magnetotail, so the enhanced signal of the ENAs was measured.The so-called plasma sheet disconnected can be understood as ENA signal weakening in the ring current region caused by the pitch angle diffusion of ring current ions during the substorm.Ogasawara et al. [2] found a positive correlation between the GENI and AE index.Unfortunately, due to the limited field of view, IBEX-Hi was not able to measure the large increase in ENA emission from the auroral zone during substorms.
IBEX orbit 52 was different.During ENA signal acquisition from 10:48 UT on 5 November 2009 to 02:23 UT on 7 November 2009, there was a small negative bay (about −6 nT) disturbance in the Dst index forecast (Figure S5) calculated by the ACE satellite real-time solar wind parameter measurement, and no substorm was generated.Although no substorms were formed, IBEX-Hi sensitively observed changes in the ion flux of the ring current and formed a transient crack in the ENA measurements, as shown in Figure 2 of Mccomas et al. [10].This indicates that the ENA emission response monitoring of the solar wind disturbance in the ring current region may be more sensitive than the AE index.

Conclusions
The terrestrial ring current is an aggregation area of energetic ions that are bound by the magnetic field of the magnetosphere with a pitch angle of about 90 • near the magnetotail equatorial plane.Similar to a cumulonimbus cloud in the meteorological sense, the spatial extent and the flux strength of the energetic ions are determinants of geomagnetic activity.The dynamic evolution of ENA emission in the ring current region and the global ENA emission flux are more sensitive to the solar wind disturbance response.During the solar maximum, the solar wind pressure surrounding the magnetopause is stronger, and the scale of the ring current is larger so that the geomagnetic activity is also stronger.The effective monitoring of the ring current will be the key data support for the space weather forecast.The Moon-based platform will be an ideal choice for ENA imaging monitoring of the terrestrial ring current.
Based on the motion characteristics of the ring current energetic ions, there is only one order of magnitude difference between the ENA counts generated in the equatorial plane and the auroral zone, which is within the effective range of the ENA imager measurements.The simulated ENA image can characterize the ring current ion flux distribution pattern of the magnetosphere.The lunar base two-dimensional ENA imager can be used to

Figure 1 .
Figure 1.Schematic diagram of the substorm model, including three key node phenomen in red), two connection processes (in green), and two driving factors (in black).

Figure 1 .
Figure 1.Schematic diagram of the substorm model, including three key node phenomena (marked in red), two connection processes (in green), and two driving factors (in black).

Figure 2 .
Figure 2. Schematic diagram of lunar orbit.Three simulation positions are labeled as first quarter, last quarter, and full moon.

Figure 2 .
Figure 2. Schematic diagram of lunar orbit.Three simulation positions are labeled as first quarter, last quarter, and full moon.

Figure 3 .
Figure 3. Simulated ion flux distribution of equatorial ring current before the ion pitch angle diffusion, where the Kp index is 0, 3, 5 and 7 from (A-D) respectively.

Figure 3 .
Figure 3. Simulated ion flux distribution of equatorial ring current before the ion pitch angle diffusion, where the Kp index is 0, 3, 5 and 7 from (A-D) respectively.

Figure 4 .
Figure 4. Pitch angle distributions of different Kp indices.

Figure 4 .
Figure 4. Pitch angle distributions of different Kp indices.

Figure 5 .
Figure 5. Simulation distribution of Luna based ENA imaging at full moon, where the order of K indices from (a-d) is the same as in Figure 3.

Figure 6 .
Figure 6.Simulation distribution of Luna based ENA imaging at the first quarter with the same K index order from (a-d) as Figure 5.

Figure 5 .Figure 5 .
Figure 5. Simulation distribution of Luna based ENA imaging at full moon, where the order of Kp indices from (a-d) is the same as in Figure 3.

Figure 6 .
Figure 6.Simulation distribution of Luna based ENA imaging at the first quarter with the same index order from (a-d) as Figure 5.

Figure 6 .
Figure 6.Simulation distribution of Luna based ENA imaging at the first quarter with the same Kp index order from (a-d) as Figure 5.

Figure 7 .
Figure 7. Simulation distribution of Luna based ENA imaging at the last quarter, the Kp index is 5, 3, and 0 from (a-d), in the reverse order of Figure 5.

Figure 7 .
Figure 7. Simulation distribution of Luna based ENA imaging at the last quarter, the Kp index is 7, 5, 3, and 0 from (a-d), in the reverse order of Figure 5.

4. 3 .
Examples of ENA Telemetry of Ring Current during Quiet/Substorms IBEX-Hi ENA scanning images from orbit 51 and 52 near apogee [10] (Figure S6) indicate that there are indeed ENA signals from the direction of magnetotail during geomagnetic quiet.These ENA signals with flux variations may serve as examples of the telemetry of ring current ion flux during geomagnetic substorms and quiet periods.
of Mccomas et al. [10] (Figure S7), IBEX-Hi's ENA count measurement, covering most of the magnetotail region from 26 October 2009 to 2 November 2009, shows that the ENA signal disappeared quickly at dusk on 29 October 2009, although the field of view of IBEX-Hi remained facing the plasma sheet until it entered the magnetotail on 2 November.Considering the position error caused by the width of the IBEX-Hi field of view, we believe that ENAs measured by IBEX-Hi at a distance of about 50 Re should be emitted from the ring current.The ENA flux sunk from 500 to 250 cm −2 sr −1 s −1 keV −1 in the direction of the magnetotail observed by IBEX-Hi at orbit 51 over the interval from 21:21 UT on 27 October 2009 to 13:40 UT on 29 October 2009, as shown in Figure 3 of Mccomas et al.

Table 3 .
Statistics of ENA counts.

Table 3 .
Statistics of ENA counts.