Short-duration Electron Precipitation Studied by Test Particle Simulation
Short-duration Electron Precipitation Studied by Test Particle Simulation
Journal of Astronomy and Space Sciences. 2015. Dec, 32(4): 317-325
Copyright © 2015, The Korean Space Science Society
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : November 18, 2015
  • Accepted : December 04, 2015
  • Published : December 15, 2015
Export by style
Cited by
About the Authors
Jaejin, Lee
University of Science and Technology, Daejeon 34113, Korea
Kyung-Chan, Kim
Korea Astronomy and Space Science Institute, Daejeon 34055, Korea
Jong-Gil, Lee
University of Science and Technology, Daejeon 34113, Korea

Energy spectra of electron microbursts from 170 keV to 340 keV have been measured by the solid-state detectors aboard the low-altitude (680 km) polar-orbiting Korean STSAT-1 (Science and Technology SATellite). These measurements have revealed two important characteristics unique to the microbursts: (1) They are produced by a fast-loss cone-filling process in which the interaction time for pitch-angle scattering is less than 50 ms and (2) The e-folding energy of the perpendicular component is larger than that of the parallel component, and the loss cone is not completely filled by electrons. To understand how wave-particle interactions could generate microbursts, we performed a test particle simulation and investigated how the waves scattered electron pitch angles within the timescale required for microburst precipitation. The application of rising-frequency whistler-mode waves to electrons of different energies moving in a dipole magnetic field showed that chorus magnetic wave fields, rather than electric fields, were the main cause of microburst events, which implied that microbursts could be produced by a quasi-adiabatic process. In addition, the simulation results showed that high-energy electrons could resonate with chorus waves at high magnetic latitudes where the loss cone was larger, which might explain the decreased e-folding energy of precipitated microbursts compared to that of trapped electrons.
Electron microbursts represent the precipitation of electrons during periods of less than ~1 sec at L = 4-8 and 06-18 local times (Parks et al. 1965 ; Parks 1978) . They were discovered by X-ray balloon-borne experiments in the early sixties (Anderson & Milton 1964) . Microbursts in the 20-100 keV energy range can be characterized by exponential energy spectra with e-folding energies E 0 ~5-20 keV (Anderson et al. 1966 ; Rosenberg et al. 1990 ; Reinard et al. 1997) .
Cyclotron resonance between electrons and magnetic chorus waves has been suggested as a possible mechanism for microburst generation (Rosenberg et al. 1990 ; Skoug et al. 1996 ; Lakhina et al. 2010 ; Lee et al. 2014) . Here, the term chorus refers to bursty whistler-mode waves with elements of ~0.1-1 sec. The elements comprise coherent sub-elements of ~5-10 msec durations with amplitudes of 0.2-0.3 nT (Tsurutani et al. 2009) . The timescales of chorus waves are similar to those of electron microbursts, and the coincident occurrence of microbursts and VLF chorus waves supports the idea that microbursts are produced by resonant interactions between electrons and cyclotron waves (Rosenberg et al. 1981) .
Imhof et al. (1992) and Nakamura et al. (2000) reported observations of impulsive electron precipitations with energies exceeding 1 MeV by spacecraft; these precipitations were named relativistic microbursts. However, relativistic microbursts have not been observed by balloon-borne experiments, and at this time it is unknown whether these microburst phenomena are related or originate from the same source.
Quasi-linear theory has been widely used to study pitch angle scattering by wave-particle interactions. In quasilinear theory, the diffusion equation is solved by assuming stochastic particle scattering caused by a succession of small-amplitude and random-phase waves (Lyons 1974 ; Summers 2005 ; Albert 2005 ; Glauert & Horne 2005) . Quasi-linear theory provides a good description of the average properties of the diffusion processes, but omits particle trapping and highly nonlinear effects (Horne et al. 2003) . Test particle simulation methods, in which the trajectories and pitch angles are calculated stepwise for each particle, include these non-linear effects, but the required calculation time to simulate a complete particle distribution is impractically large (Chang & Inan 1985 ; Inan 1987 ; Rosenberg et al. 1990) .
Bortnik & Thorne (2007) simulated electron microbursts generated by whistler-mode chorus wave interactions. They calculated the time-averaged pitch-angle changes of resonant particles from an equation proposed by Bell (1984) and estimated the precipitating electron flux by a convolution method. The simulation showed a short-timescale (~0.6 sec) electron precipitation with a peak energy of ~100–300 keV. Hikishima et al. (2010) used a self-consistent full-particle simulation to show that microburst-like precipitations of electrons with energies of 10–100 keV could accompany the generation of discrete chorus-wave emissions. These studies successfully reproduced microburst-like structures by resonant interactions between electrons and chorus waves, supporting the proposed generation mechanism of resonant interactions for microbursts. However, these studies did not address the detailed physical nature of the operations of the interactions in the production of microbursts.
In this study, we used a test particle simulation to investigate the details of resonant wave-particle interactions for electrons in the energy range of ~170–360 keV, including the differences between the fluxes in the parallel and perpendicular directions. We then compared the test particle simulation results with observations made by the Korean STSAT-1 (Science and Technology SATellite), showing that the test particle simulation could reproduce several observed features. These results further support the idea that wave-particle interactions involving cyclotron resonance can generate fast pitch-angle diffusion and thus produce microburst-like precipitations.
Two Solid-State Telescopes (SSTs) on STSAT-1 were launched into a sun-synchronous low-altitude (680 km) orbit on September 27, 2003. The first and second SSTs were aimed perpendicular and upward or parallel, respectively to the direction of the Earth’s magnetic field. The SSTs consisted of 300-μm-thick Si detectors designed to measure electrons in the energy spectra of ~190–360 keV (perpendicular) and ~170–330 keV (parallel). Each SST had a 33.9° Field of View (FOV) with a geometric factor of 0.045cm2 sr. To block ionic and UV radiation from the sun, an Al-coated Lexan foil was placed in front of each detector. Each detector had 30 energy channels with time resolutions of 50 msec, sufficient to resolve microburst structures (Lee et al. 2005 ; Lee et al. 2012) .
STSAT-1 operated successfully over approximately 300 passes through the auroral zone until March 31, 2005. During this period, electron microbursts were detected on six different occasions. An important feature of these observations was the absence of any time delay between the appearance of trapped (perpendicular) and precipitated (parallel) microbursts, which indicated that microbursts were generated by prompt loss-cone filling within a timescale of ~50 msec (the time resolution of the SSTs) or less. If we define the diffusion coefficient ( D ) in terms of time Δt for the electron pitch angles ( α ) to diffuse in a Gaussian distribution, with a standard deviation of Δα = ( D Δt ) 1/2 , and furthermore assume that microbursts are generated on the equator where the B-field is ~100 nT, we show that electrons must diffuse approximately 2.4° in 50 msec, with diffusion coefficients exceeding 3.5 × 10 -2 rad 2 /sec (Lee et al. 2005) .
Electron microbursts have been suggested to generate from cyclotron resonance with chorus waves (Rosenberg et al. 1981) . Whistler-mode chorus waves possess circularly polarized magnetic components and elliptically polarized electric components (Verkhoglyadova et al. 2009 ; Verkhoglyadova et al. 2010) . For simplicity, we will assume that both electric and magnetic wave components are circularly polarized. Fig. 1 depicts schematically the interactions of whistler-mode waves with electrons. Both electric ( E w ) and magnetic ( B w ) wave fields can scatter the pitch angles of electrons. However, the specific ways in which they interact with electrons differ significantly.
PPT Slide
Lager Image
Schematics showing interactions between electrons and cyclotron waves. (a) The sequence of particle interactions with wave B-field (from left to right). (b) The sequence of particle interactions with wave E-field.
Let us first consider a magnetic-component wave interaction, shown in Fig. 1 (a). B w is the magnetic field of the wave, B 0 is a static background magnetic field, and V e indicates the electron velocity vector. At the initial state, assuming the electron moves parallel to the background magnetic field where the pitch angle with respect to the background magnetic field is zero, however the actual pitch-angle of the electron is θ , because the total magnetic field is the vector sum of B w and B 0 . While an electron rotates through 180° with respect to the total B-field by gyro motion, B w also rotates relative to the background magnetic field. At this point, the electron has a pitch angle of 2 θ with respect to the background magnetic field. After rotating through another 180°, the pitch angle of the electron increases to 4 θ . While a conceptual understanding of the interactions of electrons with magnetic wave fields is simple, the quantitative effects of these interactions are more difficult to describe, because we cannot solve the equation of motion analytically. However, in this simple diagram, we see that the pitch angle can be changed more quickly when the background magnetic field is smaller, or when the wave magnetic field is stronger. Thus, the wave magnetic field scattering process is expected to be efficient in the equatorial region or in minimum B pockets where the magnetic field is minimum (Tsurutani et al. 2009) . For a wave magnetic field B w (~100 pT) much smaller than B 0 (~100 nT at the equator), the electron has a gyro-frequency of the order ~10 kHz. One can see that the pitch angle of the electron can change by several degrees in a few msec.
Meanwhile, the electric wave component E w shown in Fig. 1 (b) can accelerate or decelerate electrons. The perpendicular velocity increases or decreases continuously, changing the pitch angle as the resonance condition is satisfied. The electron pitch angle at time t can be expressed by Eq. (1), where α 0 is the initial pitch angle, e is the electron charge, me is the electron mass, and ve is the parallel component of the electron velocity vector. This equation shows that the pitch angle changes more quickly when the parallel electron velocity is smaller, which indicates that electric field waves can more efficiently scatter the electron pitch angles of lower-energy electrons.
PPT Slide
Lager Image
We now show the results and discuss electron pitch-angle variations as calculated by the test particle simulation in a uniform background magnetic field. The assumption of a uniform magnetic field of course does not reflect the real situation. However, this simple first-step approach provides an idea of the timescales of pitch-angle scattering. In Section 6, we extend this simulation model to the dipole field for a more realistic microburst generation.
PPT Slide
Lager Image
We solve numerically the equation of motion expressed in Eq. (2),where me is the electron mass, γ is the relativistic factor, q is the charge,and v is the electron velocity vector. The circularly polarized magnetic and electric fields. The wave particle interaction requires a resonance condition;Eq. (3) must be satisfied:
PPT Slide
Lager Image
where k|| and v|| are the parallel components of the wave vector and particle velocity vector, respectively, relative to the background magnetic field B 0 , Ω c is the electron cyclotron frequency, γ is the relativistic factor, and n is an integer denoting the cyclotron harmonic. For electrons at 300 keV interacting with coherent whistler-mode waves at 1 kHz propagating parallel to B 0 , the cyclotron frequency is ~15.1 kHz when B 0 = 542 nT (the plasma frequency is 260 kHz). Note that both electrons and waves should move toward each other to satisfy the resonance equation of Eq. (3).
Fig. 2 shows the simulation results of electrons interacting with whistler-mode magnetic and electric wave fields. The red line corresponds to the pitch-angle variation when only magnetic waves of 0.1 nT amplitude are applied to the 300 keV electrons having initial pitch angles of 0°. As time increases, the pitch angle oscillates between a maximum of 7.5° and minimum of 0°. In Fig. 1 (a), the electron total velocity was unchanged during the wave interaction, because the process was adiabatic. If the pitch angle were increased, then the decrease in parallel velocity would break the resonance condition of Eq. (3). Therefore, the pitch angles are confined to a limited range, within which they oscillate as shown by the red line in Fig. 2 .
PPT Slide
Lager Image
Electron (300 keV) pitch-angle change by circularly polarized magnetic (0.1 nT) and electric (3 mV/m) wave fields. While magnetic wave interaction produces pitch-angle oscillation, the electric wave field linearly increases the pitch-angle.
When 3 mV/m electric waves are applied, the pitch angle increases linearly, as shown by the blue line in Fig. 2 . In the schematic example of electric wave field interaction shown in Fig. 1 (b), the parallel velocity component is unchanged and the resonance condition is satisfied during the interaction. In this case, only the perpendicular velocity component is changed, and the pitch angle increases linearly in the non-relativistic case. In a relativistic situation, the increase of velocity changes the relativistic factor γ in Eq. (3), and also breaks the resonance condition.
Between magnetic and electric wave fields, one may ask which wave component contributes more dominantly to generate microbursts. This depends on the amplitudes and interaction times of the waves. In the cold-plasma model, the electric and magnetic waves are related by Eq. (4):
PPT Slide
Lager Image
where ω is the wave frequency and k is the wave vector derived from the dispersion relation of the whistler-mode wave. In this equation, a magnetic wave field strength of 0.1 nT corresponds to an electric field of ~ 3mV/m, similar to the condition described above. If the interaction time is less than 13 msec ( Fig. 2 ), the magnetic wave fields dominate pitch-angle diffusion, while electric fields dominate only for long-term interactions. Therefore, the interaction time is important in determining which waves dominate the generation of microbursts.
Fig. 2 shows the results of pitch-angle variation for electrons with initial pitch angles of 0°. If the initial pitch angle is not 0°, the motion of the electron depends on the gyro-phase. Fig. 3 shows the results of a test particle simulation in which αINITIAL is the initial pitch angle of the electrons and αFINAL is the final pitch angle after interacting with waves for 3 msec (see Eq. 2). Black dots indicate electrons of different gyro-phases. In the absence of wave interaction, αFINAL would be equal to αINITIAL , as indicated by the black diagonal line in Fig. 3 . With wave interaction, αFINAL depends on the gyro-phases and spread as the interaction time increases. Note the 0° initial pitch angle has no dependence on the gyro-phases and moves to a final pitch angle of 2.1° ( Fig. 3 ). An electron having an initial pitch angle of 2.1° can be changed to one having a final pitch angle of 0°. Hence, if αFINAL is less than the loss-cone angle, the electrons are lost into the Earth’s atmosphere.
PPT Slide
Lager Image
Pitch-angle diffusion simulated by test particle model. αINITIAL is the initial pitch angle of electrons and αFINAL is the final pitch angle after interacting with waves for 3 ms. The pitch-angle scattering by whistler-mode waves depends on the gyrophase of electrons.
Before we show the test particle simulation results in a dipole magnetic field, we discuss first how the resonance condition is satisfied in a dipole field. In Eq. (3), the gyrofrequency Ω c is a function of the magnetic field. The wave vector k can be expressed as a function of w from the dispersion relation Eq. (5), where ωp is the plasma frequency depending on plasma density. If we assume that the plasma content of a magnetic flux tube is constant along the length of the flux tube, then the plasma density depends linearly on the magnetic field. In this case, ωp is also a function of the magnetic field, as shown in Eq. (6), where np is the thermal plasma density. For the whistler-mode waves of 1 kHz, the resonance condition can be expressed as a function of the magnetic field, and we can solve Eq. (3) for magnetic latitude and for fixed-energy electrons in a dipole field. The assumption that the plasma density depends on the background geomagnetic field may be false. However, whistler-mode waves do propagate through magnetic flux tubes, and the plasma density can be expressed by magnetic latitude.
PPT Slide
Lager Image
PPT Slide
Lager Image
Fig. 4 shows a plot of
PPT Slide
Lager Image
as a function of the magnetic latitude when 1 kHz whistler-mode waves are applied to electrons with energies of 10 keV, 100 keV, 1 Mev, and 10 MeV. Here, we assume a plasma density of 2 cm – 3 at the equatorial region, L of 5.5, and anti-parallel electron movement relative to k. Notably, when
PPT Slide
Lager Image
the resonance condition is satisfied, and wave-particle interaction occurs. As shown in Fig. 4 , electrons at 10 keV do not satisfy the resonance condition, so their pitch angles are unchanged by wave interactions. However, electrons at 100 keV interact with waves in the equatorial region where the loss cone has a small angle, and electrons can be scattered into the loss cone. This may explain why microbursts are not detected at 10 keV, while 100 keV electron microbursts have been observed (Datta et al. 1997 ; Lee et al. 2005) .
PPT Slide
Lager Image
Resonance condition as a function of the magnetic latitude for waves interacting with electrons at 10 keV, 100 keV, 1 MeV, and 10 MeV.
Electrons with energies in the MeV range can interact with waves at high latitudes, where the loss cone angle is larger. In this case, the electrons would have longer angular paths to reach the point at which the pitch angle is zero. Therefore, it is difficult for high-energy electrons to be scattered into the loss cone; higher fluxes would be observed for microbursts with larger pitch angles outside the loss cone than that observed for microburst electrons traveling parallel to the magnetic field direction.
In our test particle simulation, we assume that waves are generated at the equatorial region of L = 5.5 and propagate pole-ward along the direction of the magnetic field. The waves belong to the whistler mode, with rising frequencies similar to those of chorus waves:
PPT Slide
Lager Image
where Bm is the wave magnetic field,
PPT Slide
Lager Image
is the wave frequency, t is time, k is the wave vector calculated from Eq.(6), z is the distance the waves travel, and φ is the wave phase. The wave frequency varies from 800 Hz to 1,200Hz during a time interval of 250 msec. Notably, lower-frequency waves can interact with higher-energy electrons at the same magnetic latitude (see Eq.(3)). This implies that low-frequency waves interact with high-energy electrons at low latitudes and can induce pitch-angle diffusion in the electrons. In our simulation, we assume the application of rising-tone waves to simulate interactions with chorus waves. However, any whistler-mode wave, including falling-tone or broadband, can similarly diffuse electrons at similar amplitudes.
Chorus waves can have amplitudes of 200-300 pT for timescales of~5-10 msec (Tsurutani et al. 2009) . We assume constant wave amplitudes of 150 pT for all wave frequencies: the waves move at a group velocity expressed by Eq.(8) derived from the cold plasma model. We assume that the chorus propagates parallel to the magnetic field. It has been shown by ray-tracing methods that off-axis chorus waves will propagate obliquely relative to the magnetic field (Bortnik et al. 2007) . We will follow this study with further research to examine more complex cases in detail.
PPT Slide
Lager Image
Low plasma density produces a small wave vector k in Eq. (5), and provides an environment in which high-energy electrons can interact with waves at low latitudes. While plasma density data is not available during our microburst events, it is known that the plasma density is ~1 cm -3 at L = 6 (Carpenter & Anderson 1992) . We assume a plasma density of 2 cm -3 at the equator; the density depends linearly on the magnetic field. The electric wave field is derived from Eq. (4).
Fig. 5 shows an example of electron pitch-angle changes due to wave-particle interactions in a dipole field. When the wave packet described by Eq. 7 arrives at the magnetic latitude of 25° from the equator, electrons begin interacting with the wave front at 800 Hz. We assume the electrons have initial pitch angles of 5°. At the beginning of the interaction (~25°), the electron pitch angles oscillate at small amplitude because the resonance condition is not fully satisfied. When the electrons move to the lower magnetic latitude of 21°, the wave frequency is increased slightly, thus satisfying the resonance condition. In this region, the pitch angles change dramatically over time scales of 7 msec; the amount of pitch angle variation depends on the gyro-phases of the electrons. Fig. 5 shows the trajectories of four electrons having different gyro-phases.
PPT Slide
Lager Image
Changes of electron pitch angle by interactions with waves in a dipole field. When a wave packet arrives at the magnetic latitude of 25°, electrons having an initial pitch angle of 5° begin interacting with the waves.
When the electrons reach the magnetic latitude of 18°, the resonance condition is not satisfied; the pitch angles show small-amplitude oscillations. After the electrons pass through the latitude of 6°, where the waves are terminated at the frequency of 1,200 Hz, they move adiabatically along the magnetic field. Some electrons of various pitch angles and phases can become aligned with the background magnetic field and precipitate into the Earth’s atmosphere. We suggest this mechanism for electron microburst generation.
In order to reproduce the observed electron microburst structures, our simulation used a large number of test particles in a dipole field. As the waves traveled along the magnetic field, we calculated the electron motion for particles of six different energies of 170, 200, 230, 260, 290, and 320 keV. We considered electrons having pitch angles only approaching that of the loss cone, from pitch angles of 3.3° to 10.3° at the magnetic equator. Each electron had a different initial gyro-phase. Each set of test particle electrons began to interact with the waves at different magnetic latitudes, from 5° to 45° at intervals of 1°. The total number of test particles used in this simulation was the six energies × 720 gyro-phases × 280 pitch angles × 40 latitudes. We calculated the motion of each electron with an adaptive step-size control using the Runge-Kutta method. A parallel computing machine having 32 processing cores was used for the calculations.
Information on the pitch-angle distribution when microbursts occur is not available, but the distribution could depend on the location and magnetic activity (West et al. 1973). We assumed the electron pitch-angle distribution at the equator as shown by the black solid line in Fig. 6 , where the electron loss-cone angle is 3.3° (corresponding to a pitch angle of 90° at an altitude of 100 km). The electron flux increased linearly up to the pitch angle of 8°. Presumably, the equatorial pitch angle corresponding to 90° at the satellite altitude of 680 km is located close to the boundary of the loss cone. When electrons are diffused into the loss cone by wave-particle interactions, as indicated by the dotted line, the electron flux increases would be observed by both the parallel and perpendicular detectors on STSAT-1. We have indicated the parallel and perpendicular detectors as blue and red rectangles, respectively.
PPT Slide
Lager Image
Electron pitch-angle distribution on the equator. The solid line shows the initial pitch-angle distribution; the dotted line shows the distribution after pitch-angle scattering. Blue and red rectangles represent the parallel and perpendicular detector windows, respectively, at an altitude of 680 km.
After the electrons interact with the waves, they move adiabatically along the magnetic field, ultimately reaching the STSAT-1 detectors, whose FOV can detect electrons with pitch angles of ~0–17° and ~73–90° for the parallel and perpendicular detectors, respectively, at 680 km altitude. Fig. 7 shows the simulation results. The top two panels show the perpendicular and parallel electron spectrograms, and the subsequent panels show the e-folding energies and total electron fluxes. The simulation results agree well with the STSAT-1 observations. The time shown (at the bottom) is the sum of three time intervals: the time for wave propagation from the equator to the magnetic latitude at which electrons begin interacting with waves ( Tw ), the wave-particle interaction time ( Ti ), and the time for electrons to move from the latitude at which the wave and electron interactions are terminated to the altitude of 680 km ( Tt ). The electron flux is accumulated for a time interval of 50 msec, equal to the time resolution of the STSAT-1 detectors. As shown in Fig. 7 , the electron fluxes are enhanced from 0.3 to 0.8 sec, a signature similar to the electron microbursts reproduced by the simulation. The STSAT-1 observation showing the simultaneous increase of microburst electron fluxes in electrons traveling both perpendicular and parallel to the magnetic field is interpreted as a result of the fast loss-cone filling process; we find similar electron flux increases in the simulation results. Similarly, the simulations also reproduce the result showing the higher e-folding energy of perpendicular microbursts compared to that of microbursts in the parallel direction. In addition, we also see an anisotropic pitch-angle distribution.
PPT Slide
Lager Image
An electron microburst-like precipitation obtained by the test particle simulation. From top to bottom: perpendicular and parallel electron spectrograms, e-folding energy, and total flux summed from differential fluxes shown in the top and second panels.
The test particle simulation can reproduce the microburstlike characteristics observed by STSAT-1. The important findings of our simulation study are the effects of duration and frequency bandwidth of the chorus waves on the microburst precipitation timescale, and the dominance of magnetic wave fields in the pitch-angle scattering, which indicates that the process is quasi-adiabatic.
The simulated electron microburst has a timescale of ~0.5 sec, consistent with those observed by STSAT-1. The main factor determining microburst timescale is the wave occurrence duration (250 msec in this simulation). However, the frequency bandwidth also significantly affects the microburst timescale. If we were to use monochromatic waves instead of multi-spectral waves, the resonance condition would be satisfied only for electrons with a specific energy at specific magnetic latitudes. In this monochromatic case, the electron microburst duration at fixed electron energy would be equal to the duration of the waves. However, if the waves have a frequency bandwidth, mono-energetic electrons can interact with the waves at multiple latitude points. In this case, the microburst duration is increased by the time required for the wave to travel through the interacting region. Thus, the observed electron microburst duration of ~0.25–1 sec appears to result from a combination of wave duration and frequency bandwidth.
Both electric and magnetic wave fields can cause pitch-angle scattering. In our simulation, we found the contribution of electric wave fields for pitch angle changes is ~20%, because the wave-particle interaction time of ~7 msec is very short. This short interaction time results from the non-uniform background magnetic field (a dipole field in our simulation). The dominant contribution of the magnetic wave field suggests that electron microbursts may be generated by quasi-adiabatic processes.
It has been suggested that chorus waves might be crucial to the acceleration of radiation-belt electrons (Boltnik & Thorne 2007; Shprits et al. 2006). While the proposed short wave-particle interaction for generating microbursts does not include significant electron acceleration, it is possible that multi-interactions with chorus waves could accelerate the electrons.
The simulated microburst ( Fig. 7 ) has a fast-rising and slow-falling shape, while the STSAT-1 observations generally show slow-rising and fast-falling microburst structures. Various types and shapes of microbursts have been observed; Datta et al. (1997) showed that the different shapes of these microbursts could correlate to variations of the diffusion coefficient. The diffusion rate depends on the wave amplitude in the wave-particle interaction process, implying that microburst shapes are determined by whistler-mode wave structures. In our simulation, we applied waves of constant amplitude, yet chorus waves have substructures with amplitudes that might be changed during propagation (Tsurutani et al. 2009) .
The simulation results indicate that high-energy electrons can interact with waves at high latitudes, where the loss cone is large and the background magnetic field is strong ( Fig. 4 ). In addition, the interaction time is shorter because of the higher parallel speeds. At high latitudes, the strong background magnetic field reduces the magnetic wave field effect, as shown in Fig. 1 (a). The electric wave field interaction is independent of the background magnetic field, yet only electrons having large parallel velocities can satisfy the resonance condition at high latitudes. Thus, electric and magnetic wave field interactions are not effective for scattering the pitch angles of electrons at high latitudes.
This may explain why only low-energy (<100 keV) microbursts have been observed by balloon-borne experiments that detect the X-rays produced by low-energy field-aligned electrons. Notably, relativistic-energy microbursts have only been observed thus far in situ on spacecraft-borne detectors. These spacecraft did not utilize altitude-control systems, and the detectors measured both high-energy electrons in the perpendicular direction and lower-energy electrons in the parallel direction. Lee et al. (2005) showed that the e-folding energy of microbursts increased following magnetic storms, allowing in-situ measurements to detect relativistic microbursts. It is thus possible that low-energy (< 100 keV) and relativistic electron microbursts arise from the same source.
Our test particle simulation model is not self-consistent. Moreover, we assumed constant wave amplitudes with a simple cold plasma model. Nevertheless, our model successfully simulated microburst-like precipitation and electron energy spectra similar to those measured by STSAT-1. This model will be improved to describe the detailed temporal and spatial structures of electron microbursts in future studies. However, at the present time, little information is available on the spatial scales of relativistic microbursts. Future measurements of spatial scales with a constellation of small satellites (for example, CubeSats) can provide important information on the exact spatial and temporal scales of relativistic microbursts.
This work was supported by NMSC (National Meteorological Satellite Centre) of KMA (Korea Meteorological Administration) with the research project of “Geostationary Meteorological Satellite Ground Segment Development”.
Albert JM (2005) Evaluation of quasi-linear diffusion coefficients for whistler mode waves in a plasma with arbitrary density ratio J. Geophys. Res. 110 A03218 -
Anderson KA , Milton DW (1964) Balloon observations of x-rays in the auroral zone J. Geophys. Res. 69 4457 - 4479    DOI : 10.1029/JZ069i021p04457
Anderson KA , Chase LM , Hudson HS , Lampton M , Milton DW (1966) Balloon and rocket observations of auroral-zone microburst J. Geophys. Res. 71 4617 - 4589    DOI : 10.1029/JZ071i019p04617
Bell TF (1984) The Nonlinear Gyroresonance Interaction Between Energetic Electrons and Coherent VLF Waves Propagating at an Arbitrary Angle with Respect to the Earth’s Magnetic Field J. Geophys. Res. 89 905 - 918    DOI : 10.1029/JA089iA02p00905
Bortnik J , Thorne RM (2007) The dual role of ELF/VLF chorus waves in the acceleration and precipitation of radiation belt electrons J. Atmos. Sol. Terr. Phys. 69 378 - 386    DOI : 10.1016/j.jastp.2006.05.030
Bortnik J , Thorne RM , Meredith NP (2007) Modeling the propagation characteristics of chorus using CRRES suprathermal electron fluxes J. Geophys. Res. 112 A08204 -
Chang HC , Inan US (1985) Test Particle Modeling of Wave-Induced Energetic Electron Precipitation J. Geophys. Res. 90 6409 - 6418    DOI : 10.1029/JA090iA07p06409
Carpenter DL , Anderson RR (1992) An ISEE/Whistler Model of Equatorial Electron Density in the Magnetosphere J. Geophys. Res. 97 1097 - 1108    DOI : 10.1029/91JA01548
Datta S , Skoug RM , McCarthy MP , Parks GK (1997) Modeling of microburst electron precipitation using pitch angle diffusion theory J. Geophys. Res. 102 17325 - 17224    DOI : 10.1029/97JA00942
Glauert SA , Horne RB (2005) Calculation of pitch angle and energy diffusion coefficients with the PADIE code J. Geophys. Res. 110 A04206 -
Hikishima M , Omura Y , Summers D (2010) Microburst precipitation of energetic electrons associated with chorus wave generation Geophys. Res. Lett. 37 L07103 -    DOI : 10.1029/2010GL042678
Horne RB , Glauert SA , Thorne RM (2003) Resonant diffusion of radiation belt electrons by whistler-mode chorus Geophys. Res. Lett. 30 1493 -    DOI : 10.1029/2003GL016963
Imhof WL , Voss HD , Mobilia J , Datlowe DW , Gaines EE (1992) Relativistic electron microbursts J. Geophys. Res. 97 13829 - 13837    DOI : 10.1029/92JA01138
Inan US (1987) Gyroresonant pitch angle scattering by coherent and incoherent whistler mode waves in themagnetosphere J. Geophys. Res. 92 127 - 142    DOI : 10.1029/JA092iA01p00127
Lakhina GS , Tsurutani BT , Verkhoglyadova OP , Pickett JS (2010) Pitch angle transport of electrons due to cyclotron interactions with the coherent chorus subelements J. Geophys. Res. 115 A00F15 -
Lee DH , Lee DY , Shin DK , Kim JH , Cho JH (2014) A Statistical Test of the Relationship Between Chorus Wave Activation and Anisotropy of Electron Phase Space Density J. Astron. Space Sci. 31 295 - 301    DOI : 10.5140/JASS.2014.31.4.295
Lee JJ , Parks GK , Min KW , Kim HJ , Park J (2005) Energy spectra of ~170-360 keV electron microbursts measured by the Korean STSAT-1 Geophys. Res. Lett. 32 L13106 -    DOI : 10.1029/2005GL022996
Lee JJ , Parks GK , Lee E , Tsurutani BT , Hwang J (2012) Anisotropic pitch angle distribution of ~100 keV microburst electrons in the loss cone: measurements from STSAT-1 Ann. Geophys. 30 1567 - 1573    DOI : 10.5194/angeo-30-1567-2012
Lyons LR (1974) Pitch angle and energy diffusion coefficients from resonant interactions with ion-cyclotron and whistler waves J. Plasma Phys. 12 417 -    DOI : 10.1017/S002237780002537X
Nakamura R , Isowa M , Kamide Y , Baker DN , Blake JB (2000) SAMPEX observations of precipitation bursts in the outer radiation belt J. Geophys. Res. 105 15875 - 15885    DOI : 10.1029/2000JA900018
Parks GK (1978) Microburst precipitation phenomena J. Geomag. Geoelectr. 30 327 - 341    DOI : 10.5636/jgg.30.327
Parks GK , Hudson HS , Milton DW , Anderson KA (1965) Spatial Asymmetry and Periodic Time Variations of X-Ray Microbursts in the Auroral Zone J. Geophys. Res. 70 4976 - 4978    DOI : 10.1029/JZ070i019p04976
Reinard AA , Skoug RM , Datta S , Parks GK (1997) Energy spectral characteristics of auroral electron microburst precipitation Geophys. Res. Lett. 24 611 - 614    DOI : 10.1029/97GL00377
Rosenberg TJ , Siren JC , Matthews DL , Marthinsen K , Holtet JA (1981) Conjugacy of electron microbursts and VLF chorus J. Geophys. Res. 86 5819 - 5832    DOI : 10.1029/JA086iA07p05819
Rosenberg TJ , Wei R , Detrick DL (1990) Observations and modeling of wave-induced microburst electron precipitation J. Geophys. Res. 95 6467 - 6475    DOI : 10.1029/JA095iA05p06467
Skoug RM , Datta S , McCarthy MP , Parks GK (1996) A cyclotron resonance model of VLF chorus emissions detected during electron microburst precipitation J. Geophys. Res. 101 21481 - 21492    DOI : 10.1029/96JA02007
Summers D , Mace RL , Hellberg MA (2005) Pitch-angle scattering rates in planetary magnetospheres J. Plasma Phys. 71 237 - 250    DOI : 10.1017/S0022377804003186
Tsurutani BT , Verkhoglyadova OP , Lakhina GS , Yagitani S (2009) Properties of dayside outer zone chorus during HILDCAA events: Loss of energetic electrons J. Geophys. Res. 114 A03207 -
Verkhoglyadova OP , Tsurutani BT , Omura Y , Yagitani S (2009) Properties of dayside nonlinear rising tone chorus emissions at large L observed by GEOTAIL Earth Planets Space 61 625 - 628    DOI : 10.1186/BF03352937
Verkhoglyadova OP , Tsurutani BT , Lakhina GS (2010) Properties of obliquely propagating chorus J. Geophys. Res. 115 A00F19 -
West HI , Buck RM , Walton JR (1973) Electron Pitch Angle Distributions throughout the Magnetosphere as Observed on Ogo 5 J. Geophys. Res. 78 1064 - 1082    DOI : 10.1029/JA078i007p01064