Diagnosing the storm-time ring current and other magnetospheric plasma domains using relativistic electron data

Ring current ions and relativistic electrons simultaneously measured on board MOLNIYA-1 are analyzed in comparison with the ground-based magnetometer data for the period of a strong magnetic storm (|D_st|max≈230 nT). Injection of >500 keV electrons into the slot region (L≈3) near equatorial plane is occurred on time scale ≈1 hour, when, during the magnetic storm maximum, the extreme low-latitude position of auroral electrojets is reached and ring current becomes more symmetrical. Positions of both the ring current maximum and electron intensity maximum (L_max) are consistent to our previous result: |D_st|_max = 2.75 • 10⁴/L⁴_max. An extreme storm-time low-latitude position of the west electrojet center (for amplitudes of |D_st|_max up to 600 nT) is shown to be in a good consistence with this empirical dependence. It is supposed the trapped radiation boundary collapses down to L≈L_max in the course of the storm main phase.     

An investigation of the spatial variations of the energetic trapped electrons at low altitudes

We have analyzed the trapped electron data (0.19–3.2 MeV) taken by the Japanese OHZORA satellite operated at 350–850 km altitude in polar orbit during 1984–1987 near solar minimum. The electron observations reveal all the global attributes of the quiet-time electron radiation belts, such as the South Atlantic Anomaly, the electron “slot”, and the outer radiation belt regions. The electron data are in general agreement with the NASA AE-8 electron model, but there are differences, particularly with respect to distinctive local-time variations in the slot region. In this paper, we present results from analyses of variations of the electron pitch angle distributions with local time, L-shell and altitude.     

Dynamics of the low altitude secondary proton radiation belt

At the interface between the upper atmosphere and the radiation belt region, there exists a secondary radiation belt consisting mainly of energetic ions that have become neutralized in the ring current and the main radiation belt and then re-ionized by collisions in the inner exosphere. The time history of the proton fluxes in the 0.64 – 35 MeV energy range was traced in the equatorial region beneath the main radiation belts during the three year period from 21 February 1984 to 26 March 1987 using data obtained with the HEP experiment on board the Japanese OHZORA satellite. During most of this period a fairly small proton flux of −1.2 cm⁻² s⁻¹ sr⁻¹ was detected on geomagnetic field lines in the range 1.05 < L < 1.15. We report a few surprisingly deep and rapid flux decreases (flux reduction by typically two orders of magnitude). These flux decreases were also long in duration (lasting up to three months). We also registered abrupt flux increases where the magnitude of the proton flux enhancements could reach three orders of magnitude with an enhancement duration of 1–3 days. Possible reasons for these unexpected phenomena are discussed.     

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

Experimental investigations of quasistable radiation belts formed after solar proton events in September–October 1989 and March 1991 based on measurements made by “Liulin” dosimeter-radiometer on board the “MIR” space station

Since 1988 high sensitivity dosimeter-radiometer “Liulin” has been installed on board the MIR space station. Device measured absorbed dose rate and flux of penetrating particles. Results of measurements showed that after powerful solar proton events (SPE) September–October, 1989 and March, 1991 additional quasistable radiation belts were formed in the near earth space within the interval L=1.8−3.0. These “new” belts were observed as an additional maximums in flux (and sometimes dose) channels when crossing the SAA region. “New” belts were quasi stable and existed at least several months, decaying slightly after SPE. Dose to flux ratio analysis showed that major components of these belts were energetic electrons and protons arising in connection with preceding SPEs.     

Energetic neutron and gamma-ray spectra under the earth radiation belts according to “SALUTE-7”-“KOSMOS-1686” orbital complex and “CORONAS-i” satellite data

The spectra of neutrons >10 MeV and gamma-rays 1.5–100 MeV under the Earth Radiation Belts, restored from the data, obtained onboard orbital complex “SALUTE-7”-“KOSMOS-1686”, are presented. The spectra shapes are similar to those for albedo neutrons and gamma-rays, but absolute values of their fluxes (0.2 cm⁻² s⁻¹ for neutrons, 0.8 cm⁻² s⁻¹ for gamma-rays at the equator and 1.2 cm⁻² s⁻¹, 1.9 cm⁻² s⁻¹, accordingly, at L=1.9) are several times as large. It is possibly explained by the fact that most of the detected particles were produced by the cosmic ray interactions with the orbital complex matter. Neutron and gamma-ray fluxes obtained from “CORONAS-I” data are near those for albedo particles.     

Classification of solar wind structures and intense geomagnetic storm alarms with self-organizing maps

Intense geomagnetic storms (D_st < −100 nT) usually occur when a large interplanetary duskward-electric field (with Ey > 5 mV m⁻¹) lasts for more than 3 h. In this article, a self-organizing map (SOM) neural network is used to recognize different patterns in the temporal variation of hourly averaged Ey data and to predict intense storms. The input parameters of SOM are the hourly averaged Ey data over 3 h. The output layer of the SOM has a total of 400 neurons. The hourly Ey data are calculated from solar wind data, which are provided by NSSDC OMNIWeb and ACE spacecraft and contain information on 143 intense storms and a fair number of moderate storms, weak storms and quiet periods between September 3, 1966 and June 30, 2002. Our results show that SOM is able to classify solar wind structures and therefore to give timely intense storm alarms. In our SOM, 21 neurons out of 400 are identified to be closely associated with the intense storms and they successfully predict 134 intense storms out of the 143 ones selected. In particular, there are 14 neurons for which, if one or more of them are present, the occurrence probability of intense storms is about 90%. In addition, several of these 14 neurons can predict big magnetic storm (D_st  −180 nT). In summary, our method achieves high accuracy in predicting intense geomagnetic storms and could be applied in space environment prediction.     

Geminga pulsar observations with gamma-telescope Gamma-1

The Geminga light curve obtained with the “Gamma-1” telescope features two peaks separated by 0.5 ± 0.03 period. The light curve is pronounced for γ-quanta energies higher than 400 MeV. The pulsed flux upper limit (1σ) in the energy interval 50 – 300 MeV is 6·10⁻⁷ cm⁻²sec⁻¹. For energies >300 MeV the pulsed component power law spectrum has an exponent 1.1 _−0.3^+1.1 and an integral flux (1.1±0.3)·10⁻⁶ cm⁻²sec⁻¹.     

FEM code simulation of the magnetospheric proton fluxes

We report a study of the numeric solution to the diffusive transport equation for energetic protons magnetically trapped in the Earth's equatorial magnetosphere. The analysis takes into account the pertinent physical processes in this region, including deceleration of protons by Coulomb collisional interactions with free and bound electrons, the charge exchange process, cosmic ray albedo neutron decay source, and electric and magnetic radial diffusion. These results were obtained using the Finite Element Method with magnetic moment and geomagnetic L-shell as free variables. Steady state boundary conditions were imposed at L=1 as zero distribution function and at L=7 with proton distribution function extracted from ATS 6 satellite observations. The FEM-code yields unidirectional proton flux in the energy range of 0.1–1000 MeV at the equatorial top of the geomagnetic lines, and the results are found to be in satisfactorily agreement with the empirical NASA AP-8 model proton flux within the energy range of 0.5–100 MeV. Below 500 keV, the empirical AP-8 model proton fluxes are several orders of magnitude greater than those computed with the FEM-code at L<3. This discrepancy is difficult to explain by uncertainties of boundary spectrum parameters or transport coefficients.     

Interplanetary shocks and sudden impulses during solar maximum (2000) and solar minimum (1995–1996)

In this work a study is performed on the correlation between fast forward interplanetary shock parameters at 1 Astronomical Unit and sudden impulse (SI) amplitudes in the H-component of the geomagnetic field, for periods of solar activity maximum (year 2000) and minimum (year 1995–1996). Solar wind temperature, density and speed, and total magnetic field, were taken to calculate the static pressures (thermal and magnetic) both in the upstream and downstream sides of the shocks. The variations of the solar wind parameters and pressures were then correlated with SI amplitudes. The solar wind speed variations presented good correlations with sudden impulses, with correlation coefficients larger than 0.70 both in solar maximum and solar minimum, whereas the solar wind density presented very low correlation. The parameter better correlated with SI was the square root dynamic pressure variation, showing a larger correlation during solar maximum (r = 0.82) than during solar minimum (r = 0.77). The correlations of SI with square root thermal and magnetic pressure were smaller than with the dynamic pressure, but they also present a good correlation, with r > 0.70 during both solar maximum and minimum. Multiple linear correlation analysis of SI in terms of the three pressure terms have shown that 78% and 85% of the variance in SI during solar maximum and minimum, respectively, are explained by the three pressure variations. Average sudden impulse amplitude was 25 nT during solar maximum and 21 nT during solar minimum, while average square root dynamic pressure variation is 1.20 and 0.86 nPa^1/2 during solar maximum and minimum, respectively. Thus on average, fast forward interplanetary shocks are 33% stronger during solar maximum than during solar minimum, and the magnetospheric SI response has amplitude 20% higher during solar maximum than during solar minimum. A comparison with theoretical predictions (Tsyganenko’s model corrected by Earth’s induced currents) of the coefficient of sudden impulse change with solar wind dynamic pressure variation showed excellent agreement, with values around 17 nT/nPa^1/2.     

Near infrared observations of comet Halley from Vega 2

Spatial distribution of the continuum radiation in the range of 0.95–1.9 μm presumes total dust production rate of the comet of 10ρ tonne s⁻¹ (ρ is the dust material density) and its angular distribution proportional cos . Observations of the water vapor band at 1.38 μ m reveal strong jets, their time shift from the dust jet measured in situ is consistent with gas velocity of 0.82±0.1 km s⁻¹ and dust velocity of 0.55±0.08 km s⁻¹. The OH vibrational-rotational bands observed are excided directly via photolysis of water vapor. Water vapor production rate deduced from the H₂O band and OH band intensities is 8×10²⁹ s⁻¹. Intensity of the CN(0,0) band result in the CN column density of 9×10¹² cm⁻², i.e. larger by a factor of 3 than given by the violet band.     

INTERBALL statistical study of ion flow fluctuations in the plasma sheet

We use ion distribution measurements with CORALL instrument on-board the INTERBALL/Tail spacecraft to study plasma flows in the mid-tail (−9> X> −27 R_E) plasma sheet. Three velocity components computed every 2 minutes exhibit two types of velocity variations: Earthward bursty bulk flows (BBFs) and random flow fluctuations. Their properties are in a good agreement with the observations of the ISEE-2 spacecraft (Borovsky et al., 1997). The INTERBALL/Tail spacecraft configuration favors measurements of V_z component, in contrast to previous experiments in which only V_x and V_y were measured reliably. In the outer part of the plasma sheet V_y and V_z fluctuations were close to each other (variances σ(V_y) and σ(V_z) were about 160 and 110 km/s, respectively), but in the inner part at the dusk flank amplitude of V_y fluctuations increased and was 2 times higher than that of V_z component. This asymmetry of fluctuations should be taken into account during modern theoretical analysis and simulations.     

Response of the low latitude geomagnetic field to the major proton event of November 2001

A major solar flare eruption occurred at 16:20 UT on 4 November 2001, followed by strong solar radiation storm and proton event recorded by the SOHO and other interplanetary satellites. Coronal mass ejection associated with the flare event triggered an interplanetary shock, which impacted the geomagnetic field after about 33 h. The shock impact was quite intense to produce a SSC magnitude of 80 nT in the low latitude ground magnetic records followed by sharp and deep main phase (Dst −300 nT) in the first stage, following the density (N_p) enhancement. High time resolution digital magnetic field data from the equatorial and low latitude stations in India are analyzed to study the influence of various IP parameters on the intensity and duration of the magnetic storm. A double step storm was found to be in progress caused by the multiple injections. During the period of recovery, after a period of 8 h, a third stage of depression in the ground magnetic field was set in, which corresponded to the southward directed Bz. The energy transfer processes associated with the event is presented.     

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

CMB anisotropy testing from small satellites

We review the advantages and possibilities of small satellites. New results of data reduction of the satellite-borne experiment RELICT-1 are presented. For the inflation spectrum of primordial perturbations we obtained the estimate for quadrupole component 6·10⁻⁶ <ΔT₂/T<3.3·10⁻⁶. The RELICT-2 mission will provide a possibility of measurement of CBR anisotropy down to the level less than ΔT₂/T = 10⁻⁸. We present the results of engineering testing of RELICT-2 measurement system and discuss ways of improving of the radiometers sensitivities.     

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.     

Spatial coherence and electromagnetic wave generation during electron beam experiments in space

A two-dimensional electromagnetic simulation model is used to investigate the production of whistler waves in connection with electron beam experiments in space. The spectrum is observed to peak near 0.7 ω_e, and the conversion efficiency of beam energy to whistler waves is about 5 × 10⁻⁵. The whistlers can be trapped in a density trough extending out from the spacecraft and experience ducted propagation.     

Solar cycle variation of the low-altitude trapped proton flux

Under NASA's Space Environment Effects (SEE) program, we are developing new models for the low-altitude (250–1000 km, L < 1.5) trapped radiation environment based on data from the TIROS/NOAA polar orbiting spacecraft. The unique features of this data base and model include the long time series (more than one complete solar cycle) obtained from the TIROS/NOAA data and the use of a coordinate system more applicable to the low-altitude environment. The data show a strong variation (as much as a factor of 10) over the solar cycle and a hysteresis effect between the rising and falling portions of the solar cycle. Both the solar cycle variation and the hysteresis are functions of L. In addition to the hysteresis effect, the flux during a given cycle appears to be a function of the previous cycle. Superimposed on the gradual variation over the solar cycle, transient effects, correlated with solar particle events (SPEs), can be clearly seen. Comparison with the AP8 models shows that the measured flux is a factor of 2–3 higher than the model. These data have important implications for the development and use of trapped radiation models, and will also contribute to our knowledge of the source and loss mechanisms at work in the inner zone.