Whistler studies of the plasmasphere shape and dynamics

Whistler studies of the plasmapause/plasmasphere are traced from their beginnings during the IGY through the early 1960's, when extensive data from Antarctica became available. Highlights of this period include discovery of the ‘knee’ in the equatorial electron density profile, initial comparisons with results from the Lunik probes, identification of magnetic storm effects, and discovery of the duskside bulge, or region of larger plasmasphere radius, as well as smaller-scale (Δφ ≈ 20°) variations in radius with longitude. In the mid-1960's, whistlers provided the first evidence of cross-L plasma drift patterns in the outer plasmasphere. From a present day perspective, the plasmasphere is seen as a region penetrated, perhaps most efficiently in the dusk sector, by the unsteady component of high latitude electric fields. In the pre-dawn sector, post substorm outward drifts may be an aftereffect of the shielding of the plasmasphere against the steadier components of the substorm electric fields. The available indirect whistler evidence of plasmasphere erosion during large disturbances suggests that erosion occurs primarily in the dusk-premidnight sector.     

Ion cyclotron waves observed near the plasmapause

Pc2 electromagnetic ion cyclotron waves at 0.1 Hz, near the oxygen cyclotron frequency, have been observed by ISEE-1 and -2 between L = 7.6 − 5.8 on an inbound near equatorial pass in the dusk sector. The waves occurred in a thick plasmapause of width ⋍ 1 R_e and penetrated ⋍ 1 R_e into the plasmasphere. Wave onset was accompanied by significant increases in the thermal (0–100 eV) He⁺ and the warm (0.1–16 keV/e) O⁺ and He⁺ heavy ion populations. Wave polarization is predominantly left-handed with propagation almost parallel to the ambient magnetic field, and the spectral slot and polarization reversal predicted by multicomponent cold plasma propagation theory are identified in the wave data. The results are considered an example of wave-particle interactions occurring during the outer plasmasphere refilling process at the time of the substorm recovery phase.     

Towards multi-dimensional space weather modeling for energetic oxygen ions in the earth''s inner magnetosphere: Equilibrium configuration

We report the first 3+1 dimensional model development for energetic atomic oxygen ions in the Earth's radiation belts. Energetic Oxygen ions cans be supplied to the Earth's Inner magnetosphere from the sun (as a component of solar wind and solar energetic particles), from anomalous cosmic rays, and from acceleration processes acting on ionospheric atomic oxygen ions. We have built a multi-dimensional oxygen ion model in the following free parameters: geomagnetic L-shell, the magnetic moment, the second adiabatic invariant, and the discrete charge state number. Quiet time, steady state oxygen ion distributions have been obtained numerically from an assumed outer radiation zone boundary condition at L=7, average values of the radial diffusion coefficients, and standard values for the exospheric neutral densities due to the MSIS-86 upper atmosphere and exosphere neutral thermal particle density model. Average distributions of free electrons in the plasmasphere were also assumed with a mean plasmapause location just beyond L=4. We included the six lowest ionic charge states of atomic oxygen (¹⁶O) based on an existing charge exchange cross section compilation by Spjeldvik and Fritz (1978). Computed oxygen ion distributions include the resulting equilibrium structure of energy oxygen ions between 10 KeV and 100 MeV.     

连续磁活动对等离子体层演化的影响

利用嫦娥三号极紫外相机观测的2014年2月21日等离子体层极紫外对数图像,分析了一系列磁活动状态下等离子体层晨侧视角的演化.由等离子体层质子的相空间分布,模拟了 2014年2月18-22日发生一系列磁暴事件时等离子体层在磁赤道面的演化.通过观测与模拟发现,等离子体层实际的填充速度大于模拟时等离子体层的填充速度.推测昏侧...     

Plasmasphere thermal structure as measured by ISEE-1 and DE-1

Characteristics of plasmaspheric ion thermal structure are presented from a statistical survey of low-energy of (< 50 eV) ion measurements made by the Retarding Ion Mass Spectrometer (RIMS) on the DE-1 satellite. Morning and evening results are compared to illustrate diurnal trends. Typical day side temperatures range from about 4000 K in the inner plasmasphere to over 10000 K in the outer plasmasphere, while corresponding evening side temperatures range from near 2000 K to over 10000 K. Magnetic activity is found to affect the morning and evening sides somewhat differently. Temperatures are found to remain constant or increase with altitude along magnetic field lines, depending on local time and L shell. Thermal equilibrium between H+ and He+ prevails to a high degree throughout the plasmasphere. Ion temperatures from the Plasma Composition Experiment (PCE) on ISEE-1 are generally consistent with those from DE-1/RIMS, but are lower and tend to indicate more large scale structure on the day side.     

Detection of suprathermal ionospheric O+ ions inside the plasmasphere

The plasma diagnostic experiments on the AUREOL-3 satellite have revealed flows of low energy 0⁺ ions deep inside the night plasmasphere during a large substorm. Flux gradients of the 0⁺ ions were accompanied by enhancements of ELF electric field noise. The appearance of suprathermal ions at $\text{L ? 2.5 – 3}$ is interpreted within the framework of electrostatic ion-cyclotron acceleration of ionospheric ions in the diffuse auroral zone /12/ followed by a radial displacement of these ions inside the plasmasphere driven by azimuthal electric fields during substorm activity. Electrostatic oscillations observed inside the plasmasphere are apparently associated with gradient instability at the sharp boundaries of suprathermal ion flows.     

Structure and properties of the Earth's plasmasphere

The widely used concept of the plasmapause as the last closed electric field equipotential in the equatorial plane of the magnetosphere is oversimplified. The field aligned plasma motions are of substantial importance in the plasmapause formation and should be taken into account. Distributions of the main plasma parameters measured from the Prognoz-5 satellite are presented. The diurnal variations of the plasmapause height and the plasmasphere thermal properties are discussed.     

Sources and losses of ring current ions: An update

The cleft ion fountain has been identified as a prodigious source of upflowing suprathermal ionospheric plasma. Modeling efforts have traced the path of these ions from the polar ionosphere along trajectories where the ions are energized to keV energies and deposited in the near earth plasma sheet. Mass and energy dispersion of these ions accounts in a natural way for the observed variation in heavy ion content of the plasma sheet. Observations of ion composition in the plasma sheet by the AMPTE and ISEE spacecraft establish that ionospheric ions dominate in the near earth plasma sheet but solar wind ions become significant tailward. The heavy ion content of the plasma sheet increases with both solar cycle and magnetic activity. Direct injection of ionospheric ions into the ring current has been observed in the outer plasmasphere. Several mechanisms for the direct injection of ions from the plasmasphere and ionosphere into the ring current have appeared. Estimation of ionospheric source strengths and residence times have led to an estimate of the magnetospheric densities that would result solely from an ionospheric outflow populating the magnetosphere. Estimated densities were quite reasonable even without inclusion of a solar wind source of ions. Ring current ions decay primarily via charge exchange with the hydrogen geocorona, however, the roles of pitch angle diffusion and Coulomb collisions in this decay process are being clarified.

Modeling and observations of ENA by the 1SEE1 spacecraft has led to a re-affirmation of the dominant role of charge exchange in ring current decay. Ion cyclotron waves contribute to ring current decay in the dusk bulge region. The role of low frequency. (< 1 Hz) ion cyclotron waves in the plasmasphere is still unclear. Other wave modes may be responsible for the pitch angle diffusion and subsequent loss of ring current ions. Coulomb collisional energy losses from ring current O⁺ to thermal electrons are sufficient to power SAR arcs and represent an energy sink for ring current O⁺ within the plasmasphere. Coulomb collisions may be important for decay of low energy (< 10 KeV) ring current ions in the plasmasphere.     

A study of the plasmasphere density distribution using computed tomography methods from the EUV radiation data

In this paper, we introduce a new topic for imaging plasmasphere density with CT techniques. Both parallel-beam and fan-beam reconstruction formulas are derived with a special weighting function. Different from the traditional X-ray CT, the extreme ultraviolet (EUV) radiation may be totally absorbed by the Earth. Only the portion between the Earth and the detector enters the EUV sensor and be measured. The key step in our derivation is to remove the complete absorption effect with a specific weighting function. Computer-simulation studies are present to verify our formulas for the plasmasphere density.     

Plasma transport in the plasmasphere

The plasmasphere is filled with very low energy plasma upwelling from the topside ionosphere. The field-aligned distribution of this thermal ionospheric plasma is controlled by the gravitational and centrifugal potential distribution. There are two extreme types of hydrostatic plasma distribution in this field-aligned potential : the Diffusive Equilibrium distribution and the Exospheric Equilibrium distribution corresponding respectively to a saturated and to an almost empty magnetic flux tube. As a result of pitch angle scattering by Coulomb collisions an increasing number of ions escaping from the ionosphere are stored on trapped orbits with mirror points at high altitudes in the low density region. As a result of collisions the field-aligned density distribution gradually changes from exospheric equilibrium with a highly anisotropic pitch angle (cigar like) distribution to a diffusive equilibrium with a nearly isotropic pitch angle distribution. It is shown that the suprathermal ions become anisotropic much more slowly than ions of energies smaller than 1 eV. The Coulomb collision times have been estimated for flux tubes at different L values. A numerical simulation of the flux tube refilling process has been presented. The diurnal variation of the equatorial plasma density has been illustrated for plasma elements convected along drift paths which have a large dawn- dusk asymetry. The formation of a Light Ion Trough is discussed. Finally, evidence has also been given for the existence of a ‘plasmaspheric wind’ corresponding to a slow subsonic and continuous radial expansion of the plasma stored in the plasmasphere.     

Specification of the Earth’s plasmasphere with data assimilation

In this paper we report on initial work toward data assimilative modeling of the Earth’s plasmasphere. As the medium of propagation for waves which are responsible for acceleration and decay of the radiation belts, an accurate assimilative model of the plasmasphere is crucial for optimizing the accurate prediction of the radiation environments encountered by satellites. On longer time-scales the plasmasphere exhibits significant dynamics. Although these dynamics are modeled well by existing models, they require detailed global knowledge of magnetospheric configuration which is not always readily available. For that reason data assimilation can be expected to be an effective tool in improving the modeling accuracy of the plasmasphere. In this paper we demonstrate that a relatively modest number of measurements, combined with a simple data assimilation scheme, inspired by the ensemble Kalman filtering data assimilation technique does a good job of reproducing the overall structure of the plasmasphere including plume development. This raises hopes that data assimilation will be an effective method for accurately representing the configuration of the plasmasphere for space weather applications.     

The location and intrinsic luminosity of the hidden nucleus of NGC 1068

Two key problems on the ‘hidden’ nucleus of NGC 1068 are discussed, based on the archival Hubble Space Telescope images. First, we discuss the accurate location of the nucleus. We have found that the most probable location is only ∼ 0.″1 (∼ 7 pc) south of the UV brightest cloud. Second, we consider the intrinsic luminosity of the hidden nucleus. We show that its lower limit is as large as ∼ 2 × 10⁴⁵ erg/sec, suggesting that the luminosity is almost at the Eddington limit.     

多离子空间等离子体中快磁声波与粒子的回旋共振特性

快磁声波是空间等离子体中一种接近垂直传播的右旋极化电磁波,能够在等离子体层内外传播.快磁声波与带电粒子的回旋共振相互作用能够导致高能电子随机加速和投掷角扩散、能量质子投掷角扩散等,从而影响辐射带高能带电粒子的动态过程.分别基于完整的色散关系和高密度近似的色散关系,在不同空间等离子体条件下研究多离子空间等离子体中不同传播角的快磁声波色散曲线,并计算了快磁声波与H+,He+和O+离子的最小共振能量.结果表明,当传播角较小时,采用高密度近似与采用完整色散关系计算的离子最小共振能量没有太大差别.在中低密度中强磁场空间等离子体中,传播角≥ 88°时高密度近似色散关系会带来很大的误差,因此应利用完整色散关系计算最小共振能量.      

Wave emissions at half electron gyroharmonics in the equatorial plasmasphere region: CLUSTER observations and statistics

Intense (n + 1/2) f_ce emissions are a common phenomenon observed in the terrestrial inner magnetosphere. One of their interests is their possible effect in the pitch angle scattering of plasmasheet keV-electron, leading to diffuse auroras. In this paper, we present CLUSTER’s point of view about this topic, in the equatorial region of the plasmasphere, via a statistical study using 3 years of data. Spectral characteristics of these waves, which represent an important clue concerning their generation mechanism, are obtained using WHISPER data near perigee. Details on the wave spectral signature are shown in an event study, in particular their splitting in fine frequency bands. The orbit configuration of the four spacecraft offers a complete sampling on all MLT sectors. A higher occurrence rate of the emissions in the dawn sector and their confinement to the geomagnetic equator, pointed out in previous studies, are confirmed and described with additional details. The proximity of emission sites, both to the plasmapause layer and to the geomagnetic equator surface, seems to be of great importance in the behaviour of the (n + 1/2) f_ce wave characteristics. Our study indicates for the first time, that both the intensity of (n + 1/2) f_ce emissions, and the number of harmonic bands they cover, are increasing as the observation point is located further away outside from the plasmapause layer. Moreover, a study of the wave intensity in the first harmonic band (near 3/2 f_ce) shows higher amplitude for these emissions than previous published values, these emissions can play a role in the scattering of hot electrons. Finally, geomagnetic activity influence, studied via time series of the D_st index preceding observations, indicates that (n + 1/2) f_ce emission events are observed at CLUSTER position under moderate geomagnetic activity conditions, no specific D_st time variation being required.     

Initial nighttime ionospheric observations with advanced ionospheric probe onboard FORMOSAT-5

FORMOSAT-5 satellite was launched into a sun-synchronous orbit at 720 km altitude with 98.28° inclination on 25 August 2017. The onboard scientific payload, Advanced Ionospheric Probe (AIP) is capable of measuring topside ionospheric ion density, cross-track flow velocities, ion composition and temperature, and electron temperature. Initial observations of nighttime midlatitude ionospheric density and vertical flow velocity variations at 2230 LT sector during a few quiet magnetic days in December 2017 are studied here. Longitudinal density variations in the equatorward edge of midlatitude ionospheric trough (MIT) region are noticed. Accompanied with this density variation, the vertical flow velocities also behave differently. Although the density difference has been stated due to zonal wind effect related to the declination of the geomagnetic field lines, the vertical flow velocity variation seems to play the opposite role. All these density and vertical flow observations in the northern winter hemisphere can only be explained by the longitudinal differences in the diffusion velocity coming down from the protonsphere (plasmasphere). In addition, the hemispheric asymmetry in the vertical flow velocity can also be explained by the interaction between the topside ionosphere and the protonsphere. The observed vertical flow variations near MIT at different longitudes should present a new potential tool for the study of MIT formation.     

A new inversion method of plasma density distribution of plasmasphere in the geomagnetic equatorial plane from IMAGE data

The plasma density distribution of plasmasphere in the geomagnetic equatorial plane can help us study the magnetosphere like plasmasphere, ionosphere and their kinetics. In this paper, we introduce a new inversion method, GE-ART, to calculate the plasma density distribution in the geomagnetic equatorial plane from the Extreme Ultraviolet (EUV) data of IMAGE satellite under the assumption that the plasma density is constant along each geomagnetic field line. The new GE-ART algorithm was derived from the traditional Algebraic Reconstruction Techniques (ART) in Computed Tomography (CT) which was different from the several existing methods. In this new method, each value of the EUV image data was back-projected evenly to the geomagnetic field lines intersected by this EUV sight. A 3-D inversion matrix was produced by the contributions of all the voxels contained in the plasmasphere covered by the EUV sensor. That is, we considered that each value of the EUV image data was relative to the plasma densities of all the voxels passed through by the corresponding EUV radiation, which is the biggest difference to all the existing inversion methods. Finally, the GE-ART algorithm was evaluated by the real EUV data from the IMAGE satellite.     

Conical ion distributions at one earth radius — observations from the acceleration region?

Thermal ion measurements from the Retarding Ion Mass Spectrometer (RIMS) on Dynamics Explorer 1 (DE 1) in the night side auroral region were surveyed for evidence of ion acceleration. The RIMS measurements showed evidence for ion acceleration in the 2–10,000 km altitude range, with ion distributions peaked near 90°, and with temperatures of 1 to 10 eV. Two illustrations of the RIMS data for such observations are given here. The conical distributions are found at the low latitude edge of the auroral region, just outside the plasmapause. In the first example, the three major ion species (H+, He+, and O+) show evidence of acceleration. The angular distributions are peaked at different pitch angles, indicating that the different species have been accelerated at different altitudes. The H+ flux is higher than the O+ flux in this first example, in the RIMS energy range (0–50 eV). This is apparently typical of the RIMS observations on the night side. In the second example, only O+ is transversely accelerated.     

Relationship of the diffuse aurora and SAR arc dynamics to substorms and storms

Investigation results of a diffuse aurora (DA) and stable auroral red (SAR) arc dynamics based on spectrophotometric observations at the Yakutsk meridian (199°E geomagnetic longitude) are presented. The relationship of an equatorward extension of DA in the 557.7 nm emission to a substorm growth phase during the magnetospheric convection intensification after the turn of IMF B_Z to the south is shown. The formation of SAR arc during the substorm expansion phase is investigated. The association of SAR arc dynamics with the development of asymmetric ring current (substorm injection) during the main phase of a storm is analyzed. It is shown how the pulsating precipitations of energetic ring current particles develop in the outer plasmasphere based on photometric observations.     

A simplified D-region ion chemistry scheme and its possible use for IRI lower ionosphere modelling

Ion composition of the D region is principally characterized by the existence of two distinct regions of predominant molecular ions and predominant cluster ions, separated from each other by a rather sharp ‘transition height’, which is proposed to be included in the IRI as an additional parameter, supplementing the electron density models. It is possible to predict the position of this ‘transition height’ at a given place and time with the aid of a simplified ion chemistry scheme which is shown to be satisfactorily compatible with experimental ion composition data available in the literature. Our suggested method of this prediction makes use of the (IRI or experimental) electron density profile at the location and season in question, together with an effective clustering rate coeeficient calculated from corresponding temperature and density profiles taken from a suitable reference model of the neutral atmosphere.