Spacecraft potential variations of the Swarm satellites at low Earth orbital altitudes

In this study, we provide the first detailed analysis of variations in the spacecraft potential (Vs) of the three Swarm satellites, which are flying at about 400-500 km. Unlike previous studies that have investigated extreme charging events, usually with spacecraft potentials as negative as −100 V, this study is focused on variations of Swarm Vs readings, which fall within a few negative volts. The Swarm observations show that spacecraft at low Earth orbital (LEO) altitudes are charged only slightly negatively, varying between −7 V and 0 V, with the majority of recorded potentials at these altitudes clustering close to −2 V. However, a second peak of Vs data is found at −5.5 V, though the event numbers for these more-negative observations are less, by an order of magnitude, than for incidents near the −2 V peak. These two distinct Vs peaks suggest two different causes. We have thus divided the Swarm spacecraft Vs data into two categories: less-negatively charged (−5 < Vs < 0 V) and more-negatively-charged (−6.5 < Vs < −5 V). These two Vs categories exhibit different spatial and temporal distributions. The Vs observations in the first category remain relatively closer to 0 V above the magnetic equator, but become much more negative at low and middle latitudes on the day side;at high latitudes, these first-category Vs readings are relatively more-negative during local summer. Second-category Vs events cluster into two bands at the middle latitudes (between ±20°-50° magnetic latitude), but with slightly more negative readings at the South Atlantic Anomaly (SAA) region;at high latitudes, these rarer but more-negative second-category Vs events exhibit relatively more-negative values during local winter, which is opposite to the seasonal pattern seen in the first category. By comparing Vs data to the distributions of background plasma density at Swarm altitudes, we find for the first category that more-negative Vs readings are recorded at regions with higher background plasma density, while for the second category the more-negative Vs data are observed at regions with lower background plasma density. This can be explained as follows: the electron and ion fluxes incident on Swarm surface, whose differences determine the potential of Swarm, are dominated by the background “cold” plasma (due to ionization) and “hot” plasma (due to precipitated particles from magnetosphere) for the two Vs categories, respectively.

A statistical analysis of the Kappa-type energy spectrum distribution of radiation belt electrons observed by Van Allen Probes

The energy spectrum of energetic electrons is a key factor representing the dynamic variations of Earth’s Van Allen radiation belts.Increased measurements have indicated that the commonly used Maxwellian and Kappa distributions are inadequate for capturing the realistic spectral distributions of radiation belt electrons.Here we adopt the Kappa-type(KT)distribution as the fitting function and perform a statistical analysis to investigate the radiation belt electron flux spectra observed by the Van Allen Probes.By calculating the optimal values of the key KT distribution parameters(i.e.,κandθ2)from the observed spectral shapes,we fit the radiation belt electron fluxes at different L-shells under different geomagnetic conditions.In this manner,we obtain typical values of the KT distribution parameters,which are statistically feasible for modeling the radiation belt electron flux profiles during either geomagnetically quiet or active periods.A comparison of the KT distribution model results with those using the Maxwellian or Kappa distribution reveals the advantage of the KT distribution for studying the overall properties of the radiation belt electron spectral distribution,which has important implications for deepening the current understanding of the radiation belt electron dynamics under evolving geomagnetic conditions.

Effect of interplanetary magnetic field B_(x)on the polar electrojets as observed by CHAMP and Swarm satellites

Based on 16 years of magnetic field observations from CHAMP and Swarm satellites,this study investigates the influence of the Interplanetary Magnetic Field(IMF)Bx component on the location and peak current density of the polar electrojets(PEJs).We find that the IMF Bx displays obvious local time,seasonal,and hemispherical effects on the PEJs,as follows:(1)Compared to other local times,its influence is weakest at dawn and dusk.(2)In the midnight sectors of both hemispheres,the IMF Bx tends to amplify the westward PEJ when it is<0 in the Northern Hemisphere and when it is>0 in the Southern Hemisphere;this effect is relatively stronger in the local winter hemisphere.(3)At noontime,the IMF Bx intensifies the eastward current when it is<0 in the Northern Hemisphere;in the Southern Hemisphere when it is>0,it reduces the westward current;this effect is notably more prominent in the local summer hemisphere.(4)Moreover,the noontime eastward current shifts towards higher latitudes,while the midnight westward current migrates towards lower latitudes when IMF Bx is<0 in the Northern Hemisphere and when it is>0 in the Southern Hemisphere.

A data assimilation-based forecast model of outer radiation belt electron fluxes

Because radiation belt electrons can pose a potential threat to the safety of satellites orbiting in space,it is of great importance to develop a reliable model that can predict the highly dynamic variations in outer radiation belt electron fluxes.In the present study,we develop a forecast model of radiation belt electron fluxes based on the data assimilation method,in terms of Van Allen Probe measurements combined with three-dimensional radiation belt numerical simulations.Our forecast model can cover the entire outer radiation belt with a high temporal resolution(1 hour)and a spatial resolution of 0.25 L over a wide range of both electron energy(0.1-5.0 MeV)and pitch angle(5°-90°).On the basis of this model,we forecast hourly electron fluxes for the next 1,2,and 3 days during an intense geomagnetic storm and evaluate the corresponding prediction performance.Our model can reasonably predict the stormtime evolution of radiation belt electrons with high prediction efficiency(up to~0.8-1).The best prediction performance is found for~0.3-3 MeV electrons at L=~3.25-4.5,which extends to higher L and lower energies with increasing pitch angle.Our results demonstrate that the forecast model developed can be a powerful tool to predict the spatiotemporal changes in outer radiation belt electron fluxes,and the model has both scientific significance and practical implications.

Proton pitch angle distributions in the Martian induced magnetosphere: A survey of Tianwen-1 Mars Ion and Neutral Particle Analyzer observations

The pitch angle distributions of ions and electrons can be affected by various processes;thus,they can serve as an important indicator of the physical mechanisms driving the dynamics of space plasmas.From observations from the Mars Ion and Neutral Particle Analyzer onboard the Tianwen-1 orbiter,we calculated the pitch angle distributions of protons in the Martian induced magnetosphere by using information from the magnetohydrodynamically simulated magnetic field,and we statistically analyzed the spatial occurrence pattern of different types of pitch angle distributions.Even though no symmetrical features were seen in the dataset,we found the dominance of the field-aligned distribution type over the energy range from 188 to 6232 eV.Maps of the occurrence rate showed the preferential presence of a trapped-like distribution at the lower altitudes of the surveyed nightside region.Although our results are more or less restricted by the adopted magnetic field,they indicate the complexity of the near-Mars proton pitch angle distributions and infer the possibility of wave–particle interactions in the Martian induced magnetosphere.

Cross-satellite calibration of high-energy electron fluxes measured by FengYun-4A based on Arase observations

We use the High-energy Electron Experiments(HEP)instrument onboard Arase(ERG)to conduct an energy-dependent cross-satellite calibration of electron fluxes measured by the High Energy Particle Detector(HEPD)onboard FengYun-4A(FY-4A)spanning from April 1,2017,to September 30,2019.By tracing the two-dimensional magnetic positions(L,magnetic local time[MLT])of FY-4A at each time,we compare the datasets of the conjugate electron fluxes over the range of 245–894 keV in 6 energy channels for the satellite pair within different sets of L×MLT.The variations in the electron fluxes observed by FY-4A generally agree with the Arase measurements,and the percentages of the ratios of electron flux conjunctions within a factor of 2 are larger than 50%.Compared with Arase,FY-4A systematically overestimates electron fluxes at all 6 energy channels,with the corresponding calibration factors ranging from 0.67 to 0.81.After the cross-satellite calibration,the electron flux conjunctions between FY-4A and Arase show better agreement,with much smaller normalized root mean square errors.Our results provide a valuable reference for the application of FY-4A high-energy electron datasets to in-depth investigations of the Earth’s radiation belt electron dynamics.

Overshoot phenomena:Observation and simulation

High-power O-mode radio waves can excite artificial instabilities in the F region,according to experiments conducted at the European Incoherent Scatter Science Association(EISCAT)heating facility.The main instabilities include the parametric decay instability(PDI),oscillating two-stream instability(OTSI),and thermal parametric instability(TPI).The PDI and OTSI not only compete with each other,but also compete with the TPI,leading to a two-stage overshoot phenomenon:a miniovershoot occurs on a millisecond time scale after pump-on,followed by the main overshoot.We gain insight into the miniovershoot via a generalized Zakharov model,whereas the main overshoot can be observed as an enhanced plasma line overshoot phenomenon in incoherent scatter radar spectra.We can also observe that the zero-frequency ion line exists only in the initial heating period after a cold start and that the upshifted and downshifted ion lines behave irregularly in the spectra.The simulation results show that competition between the PDI and OTSI leads to an initial peak,which we named the pre-miniovershoot.The following processes,namely ion density caviton generation,and collapse and cascade in the development of the PDI,contribute to the miniovershoot phenomenon.

On the loss mechanisms of radiation belt electron dropouts during the 12 September 2014 geomagnetic storm

Radiation belt electron dropouts indicate electron flux decay to the background level during geomagnetic storms,which is commonly attributed to the effects of wave-induced pitch angle scattering and magnetopause shadowing.To investigate the loss mechanisms of radiation belt electron dropouts triggered by a solar wind dynamic pressure pulse event on 12 September 2014,we comprehensively analyzed the particle and wave measurements from Van Allen Probes.The dropout event was divided into three periods:before the storm,the initial phase of the storm,and the main phase of the storm.The electron pitch angle distributions(PADs)and electron flux dropouts during the initial and main phases of this storm were investigated,and the evolution of the radial profile of electron phase space density(PSD)and the(μ,K)dependence of electron PSD dropouts(whereμ,K,and L^*are the three adiabatic invariants)were analyzed.The energy-independent decay of electrons at L>4.5 was accompanied by butterfly PADs,suggesting that the magnetopause shadowing process may be the major loss mechanism during the initial phase of the storm at L>4.5.The features of electron dropouts and 90°-peaked PADs were observed only for>1 MeV electrons at L<4,indicating that the wave-induced scattering effect may dominate the electron loss processes at the lower L-shell during the main phase of the storm.Evaluations of the(μ,K)dependence of electron PSD drops and calculations of the minimum electron resonant energies of H+-band electromagnetic ion cyclotron(EMIC)waves support the scenario that the observed PSD drop peaks around L^*=3.9 may be caused mainly by the scattering of EMIC waves,whereas the drop peaks around L^*=4.6 may result from a combination of EMIC wave scattering and outward radial diffusion.

Review of ionospheric irregularities and ionospheric electrodynamic coupling in the middle latitude region

This paper briefly reviews ionospheric irregularities that occur in the E and F regions at mid-latitudes. Sporadic E(ES) is a common ionospheric irregularity phenomenon that is first noticed in the E layer. ES mainly appears during daytime in summer hemispheres, and is formed primarily from neutral wind shear in the mesosphere and lower thermosphere(MLT) region. Field-aligned irregularity(FAI) in the E region is also observed by Very High Frequency(VHF) radar in mid-latitude regions. FAI frequently occurs after sunset in summer hemispheres, and spectrum features of E region FAI echoes suggest that type-2 irregularity is dominant in the nighttime ionosphere. A close relationship between ES and E region FAI implies that ES may be a possible source of E region FAI in the nighttime ionosphere. Strong neutral wind shear, steep ES plasma density gradient, and a polarized electric field are the significant factors affecting the formation of E region FAI. At mid-latitudes, joint observational experiments including ionosonde, VHF radar, Global Positioning System(GPS) stations, and all-sky optical images have revealed strong connections across different scales of ionospheric irregularities in the nighttime F region, such as spread F(SF), medium-scale traveling ionospheric disturbances(MSTID), and F region FAI.Observations suggest that different scales of ionospheric irregularities are generally attributed to the Perkins instability and subsequently excited gradient drift instability. Nighttime MSTID can further evolve into small-scale structures through a nonlinear cascade process when a steep plasma density gradient exists at the bottom of the F region. In addition, the effect of ionospheric electrodynamic coupling processes, including ionospheric E-F coupling and inter-hemispheric coupling on the generation of ionospheric irregularities, becomes more prominent due to the significant dip angle and equipotentiality of magnetic field lines in the mid-latitude ionosphere. Polarized electric fields can map to different ionospheric regions and excite plasma instabilities which form ionospheric irregularities. Nevertheless,the mapping efficiency of a polarized electric field depends on the ionospheric background and spatial scale of the field.

A detailed investigation of low latitude tweek atmospherics observed by the WHU ELF/VLF receiver:2.Occurrence features and associated ionospheric parameters

As a companion paper to Zhou RX et al.(2020),this study describes application of the automatic detection and analysis module to identify all the tweek atmospherics detectible in the WHU ELF/VLF receiver data collected at Suizhou station during the period of 3 February through 29 February 2016.Detailed analysis of the identified low-latitude tweek events reveals that the occurrence rate varies considerably—from 800 to 6000 tweeks per day,and exhibits a strong diurnal and local time dependence,the peak occurring before local midnight.The diurnal variation of identified tweeks was similar to that of the lightning data obtained by the World-Wide Lightning Location Network(WWLLN)..Estimates of the propagation distance and ionospheric reflection height of tweek atmospherics suggest that the majority(~92%)of the low latitude tweeks originate from the lightning activity within a radius of 4000 km and that they are very likely to reflect from the lower ionospheric D-region at the height range of 75–85 km.At these lower ionospheric reflection altitudes,~74%of the corresponding electron densities from the tweek spectral measurements are within 24.5–27.5 cm^-3.The daily variation of estimated D-region electron densities in the considered period(February 2016)also exhibits a small overall increasing trend from early to later in the month.

Evidence of X-mode heating suppressing O-mode heating

In this study,we present three experiments carried out at the EISCAT(European Incoherent Scatter Scientific Association)heating facility on October 29 and 30,2015.The results from the first experiment showed overshoot during the O-mode heating period.The second experiment,which used cold-start X-mode heating,showed the generation of parametric decay instability,whereas overshoot was not observed.The third experiment used power-stepped X-mode heating with noticeable O-mode wave leakage.Parametric decay instability and oscillating two-stream instability were generated at the O-mode reflection height without the overshoot effect,which implies suppression of the thermal parametric instability with X-mode heating.We propose that the electron temperature increased because X-mode heating below the upper hybrid height decreased the growth rate of the thermal parametric instability.

A detailed investigation of low latitude tweek atmospherics observed by the WHU ELF/VLF receiver:Ⅰ. Automatic detection and analysis method

As a dispersive wave mode produced by lightning strokes, tweek atmospherics provide important hints of lower ionospheric(i.e., D-region) electron density. Based on data accumulation from the WHU ELF/VLF receiver system, we develop an automatic detection module in terms of the maximum-entropy-spectral-estimation(MESE) method to identify unambiguous instances of low latitude tweeks.We justify the feasibility of our procedure through a detailed analysis of the data observed at the Suizhou Station(31.57°N, 113.32°E) on17 February 2016. A total of 3961 tweeks were registered by visual inspection;the automatic detection method captured 4342 tweeks, of which 3361 were correct ones, producing a correctness percentage of 77.4%(= 3361/4342) and a false alarm rate of 22.6%(= 981/4342).A Short-Time Fourier Transformation(STFT) was also applied to trace the power spectral profiles of identified tweeks and to evaluate the tweek propagation distance. It is found that the fitting accuracy of the frequency–time curve and the relative difference of propagation distance between the two methods through the slope and through the intercept can be used to further improve the accuracy of automatic tweek identification. We suggest that our automatic tweek detection and analysis method therefore supplies a valuable means to investigate features of low latitude tweek atmospherics and associated ionospheric parameters comprehensively.

Non-storm erosion of MeV electron outer radiation belt down to L^(*)<4.0 associated with successive enhancements of solar wind density

We report an unusual non-storm erosion event of outer zone MeV electron distribution during three successive solar wind number density enhancements(SWDEs)on November 27-30,2015.Loss of MeV electrons and energy-dependent narrowing of electron pitch angle distributions(PAD)first developed at L^(*)=5.5 and then moved down to L^(*)<4.According to the evolution of the electron phase space density(PSD)profile,losses of electrons with small pitch angles at L^(*)>4 during SWDE1 are mainly due to outward radial diffusion.However during SWDE2&3,scattering loss due to EMIC waves is dominant at 4<L^(*)<5.As for electrons with large pitch angles,outward radial diffusion is the primary loss mechanism throughout all SWDEs which is consistent with the incursion of the Last Closed Drift Shell(LCDS).The inner edge of EMIC wave activity moved from L^(*)~5 to L^(*)~4 and from L~6.4 to L~4.2 from SWDE1 to SWDE2&3,respectively,observed by Van Allen Probes and by ground stations.This is consistent with the inward penetration of anisotropic energetic protons from L^(*)=4.5 to L^(*)=3.5,suggesting that the inward extension of EMIC waves may be driven by the inward injection of anisotropic energetic protons from the dense plasma sheet.

An empirical model of the global distribution of plasmaspheric hiss based on Van Allen Probes EMFISIS measurements

Using wave measurements from the EMFISIS instrument onboard Van Allen Probes,we investigate statistically the spatial distributions of the intensity of plasmaspheric hiss waves.To reproduce these empirical results,we establish a fitting model that is a thirdorder polynomial function of L-shell,magnetic local time(MLT),magnetic latitude(MLAT),and AE*.Quantitative comparisons indicate that the model’s fitting functions can reflect favorably the major empirical features of the global distribution of hiss wave intensity,including substorm dependence and the MLT asymmetry.Our results therefore provide a useful analytic model that can be readily employed in future simulations of global radiation belt electron dynamics under the impact of plasmaspheric hiss waves in geospace.

Hemispheric Distribution of Lower-band Chorus Waves Observed by Van Allen Probes

Whistler mode chorus waves are important electromagnetic emissions due to their dual roles in acceleration and loss processes of Earth's radiation belt electrons.A detailed global survey of lower-band chorus is performed using EMFISIS data from Van Allen Probes in near-equatorial orbits.In addition to the confirmation of the positive correlation of chorus wave intensities to geomagnetic activity and dayside-nightside distribution asymmetry of wave amplitude and occurrence probability,the analysis results find that in statistics lower-band chorus emissions exhibit higher wave occurrence rates and larger normalized peak wave frequencies in the magnetically northern hemisphere but somehow stronger peak wave intensities in the magnetically southern hemisphere.While overall the differences between the two magnetically hemispheric distributions tend to be not significant,it is important to establish the magnetically hemispheric distribution profiles of lowerband chorus with respect to L-shell,magnetic local time,and geomagnetic latitude for improved understanding of chorus-induced dynamics of radiation belt electrons.

Triangulation of red sprites observed above a mesoscale convective system in North China

The triangulation of red sprites was obtained, based on concurrent observations over a mesoscale convective system(MCS) in North China from two stations separated by about 450 km. In addition, broadband sferics from the sprite-producing lightning were measured at five ground stations, making it possible to locate and identify the individual causative lightning discharges for different elements in this dancing sprite event. The results of our analyses indicate that the sprites were produced above the trailing stratiform region of the MCS, and their parent strokes were located mainly in the peripheral area of the stratiform. The lateral offset between sprites and causative strokes ranges from a few km to more than 50 km. In a particularly bright sprite, with a distinct halo feature and streamers descending down to an altitude of approximately 48 km, the sprite current signal identified in the electric sferic, measured at a range of about 1,110 km, peaked at approximately 1 ms after the return stroke.

Response of photoelectron peaks in the Martian ionosphere to solar EUV/X-ray irradiance

An important population of the dayside Martian ionosphere are photoelectrons that are produced by solar Extreme Ultraviolet and X-ray ionization of atmospheric neutrals.A typical photoelectron energy spectrum is characterized by a distinctive peak near 27 eV related to the strong solar HeⅡ emission line at 30.4 nm,and an additional peak near 500 eV related to O Auger ionization.In this study,the extensive measurements made by the Solar Wind Electron Analyzer on board the recent Mars Atmosphere and Volatile Evolution spacecraft are analyzed and found to verify the scenario that Martian ionosphere photoelectrons are driven by solar radiation.We report that the photoelectron intensities at the centers of both peaks increase steadily with increasing solar ionizing flux below 90 nm and that the observed solar cycle variation is substantially more prominent near the O Auger peak than near the HeⅡ peak.The latter observation is clearly driven by a larger variability in solar irradiance at shorter wavelengths.When the solar ionizing flux increases from 1 mW·m^-2 to 2.5 mW·m^-2,the photoelectron intensity increases by a factor of 3.2 at the HeⅡ peak and by a much larger factor of 10.5 at the O Auger peak,both within the optically thin regions of the Martian atmosphere.

Test particle simulations of resonant interactions between energetic electrons and discrete, multi-frequency artificial whistler waves in the plasmasphere

Modulated high frequency （HF） heating of the ionosphere provides a feasible means of artificially generating ex- tremely low frequency （ELF）/very low frequency （VLF） whistler waves, which can leak into the inner magnetosphere and contribute to resonant interactions with high energy electrons. Combining the ray tracing method and test particle simulations, we evaluate the effects of energetic electron resonant scattering driven by the discrete, multi-frequency arti- ficially generated ELF/VLF waves. The simulation results indicate a stochastic behavior of electrons and a linear profile of pitch angle and kinetic energy variations averaged over all test electrons. These features are similar to those associated with single-frequency waves. The computed local diffusion coefficients show that, although the momentum diffusion of relativistic electrons due to artificial ELF/VLF whistlers with a nominal amplitude of ~ 1 pT is minor, the pitch angle scattering can be notably efficient at low pitch angles near the loss cone, which supports the feasibility of artificial triggering of multi-frequency ELF/VLF whistler waves for the removal of high energy electrons from the magnetosphere. We also investigate the dependences of diffusion coefficients on the frequency interval （△f） of the discrete, multi-frequency waves. We find that there is a threshold value of Af for which the net diffusion coefficient of multi-frequency whistlers is inversely proportional to △f （proportional to the frequency components Nw） when △f is below the threshold value but it remains unchanged with increasing Af when △f is larger than the threshold value. This is explained as being due to the fact that the resonant scattering effect of broadband waves is the sum of the effects of each frequency in the ＇effective frequency band＇. Our results suggest that the modulation frequency of HF heating of the ionosphere can be appropriately selected with reasonable frequency intervals so that better performance of controlled precipitation of high energy electrons in the plasmasphere by artificial ELF/VLF whistler waves can be achieved.

Numerical simulation of chorus-driving acceleration of relativistic electrons at extremely low L-shell during geomagnetic storms

During 2018 major geomagnetic storm,relativistic electron enhancements in extremely low L-shell regions(reaching L∼3)have been reported based on observations of ZH-1 and Van Allen probes satellites,and the storm is highly likely to be accelerated by strong whistler-mode waves occurring near very low L-shell regions where the plasmapause was suppressed.It is very interesting to observe the intense chorus-accelerated electrons locating in such low L-shells and filling into the slot region.In this paper,we further perform numerical simulation by solving the two-dimensional Fokker-Planck equation based on the bounce-averaged diffusion rates.Numerical results demonstrate the evolution processes of the chorus-driven electron flux and confirm the flux enhancement in low pitch angle ranges(20◦-50◦)after the wave-particle interaction for tens of hours.The simulation result is consistent with the observation of potential butterfly pitch angle distributions of relativistic electrons from both ZH-1 and Van Allen probes.

Simultaneous occurrence of four magnetospheric wave modes and the resultant combined scattering effect on radiation belt electrons

We report a representative concurrent event of four wave modes at L≈5.0,including electrostatic electron cyclotron harmonic(ECH)waves,exohiss,magnetosonic(MS)waves,and electromagnetic ion cyclotron(EMIC)waves,based on the observations from Van Allen Probe A on October 15,2015.The diffusion coefficients induced by these waves are calculated by using both the Full Diffusion Code and test particle simulations.Moreover,the scattering effects of these waves on energetic electrons are simulated by using a two-dimensional Fokker-Planck diffusion model.The results show that ECH waves mainly scatter low-pitch-angle(<20°)electrons at 0.1-10 keV;exohiss can significantly scatter hundreds of kiloelectron volt electrons to form a reversed energy spectrum;MS waves mainly affect high-pitch-angle electrons(>60°);and EMIC waves scatter only>5 MeV electrons.The combined scattering effects of exohiss and MS waves are stronger than those of exohiss alone.The top-hat pitch angle distributions produced by exohiss are relaxed after adding the effect of MS waves.Because the energies of electrons scattered by ECH waves and EMIC waves are much lower and higher than those scattered by exohiss and MS waves,respectively,the combined scattering effects with the addition of ECH and EMIC waves show little difference from the results for the combination of MS waves and exohiss.These results suggest that distinct wave modes can occur simultaneously and scatter electrons in combination or individually,which requires careful consideration in future global simulations of the complex dynamics of radiation belt energetic electrons.