Atmospheric pressure variations at extratropical latitudes associated with Forbush decreases of galactic cosmic rays

Changes of troposphere pressure associated with short-time variations of galactic cosmic rays (GCRs) taking place in the Northern hemisphere’s cold months (October–March) were analyzed for the period 1980–2006, NCEP/NCAR reanalysis data being used. Noticeable pressure variations during Forbush decreases of GCRs were revealed at extratropical latitudes of both hemispheres. The maxima of pressure increase were observed on the 3rd–4th days after the event onsets over Northern Europe and the European part of Russia in the Northern hemisphere, as well as on the 4th–5th days over the eastern part of the South Atlantic opposite Queen Maud Land and over the d’Urville Sea in the Southern Ocean. According to the weather chart analysis, the observed pressure growth, as a rule, results from the weakening of cyclones and intensification of anticyclone development in these areas. The presented results suggest that cosmic ray variations may influence the evolution of extratropical baric systems and play an important role in solar-terrestrial relationships.     

Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

Low-latitude geomagnetic response to the interplanetary conditions during very intense magnetic storms

The variations in the horizontal and declination components of the geomagnetic field in response to the interplanetary shocks driven by fast halo coronal mass ejections, fast solar wind streams from the coronal hole regions and the dynamic pressure pulses associated with these events are studied. Close association between the field-aligned current density (j_∥) and the fluctuations in the declination component (ΔD_ABG) at Alibag is found for intense storm conditions. Increase in the dawn-dusk interplanetary electric field (E_y) and ΔD_ABG are generally in phase. However, when the magnetospheric electric field is directed from dusk to dawn direction, a prominent scatter occurs between the two. It is suggested that low-latitude ground magnetic data may serve as a proxy for the interplanetary conditions in the solar wind.     

The solar wind control of the equatorial ionosphere dynamics during geomagnetic storms

The interplanetary magnetic field, geomagnetic variations, virtual ionosphere height h′F, and the critical frequency foF2 data during the geomagnetic storms are studied to demonstrate relationships between these phenomena. We study 5-min ionospheric variations using the first Western Pacific Ionosphere Campaign (1998–1999) observations, 5-min interplanetary magnetic field (IMF) and 5-min auroral electrojets data during a moderate geomagnetic storm. These data allowed us to demonstrate that the auroral and the equatorial ionospheric phenomena are developed practically simultaneously. Hourly average of the ionospheric foF2 and h′F variations at near equatorial stations during a similar storm show the same behavior. We suppose this is due to interaction between electric fields of the auroral and the equatorial ionosphere during geomagnetic storms. It is shown that the low-latitude ionosphere dynamics during these moderate storms was defined by the southward direction of the B_z-component of the interplanetary magnetic field. A southward IMF produces the Region I and Region II field-aligned currents (FAC) and polar electrojet current systems. We assume that the short-term ionospheric variations during geomagnetic storms can be explained mainly by the electric field of the FAC. The electric fields of the field-aligned currents can penetrate throughout the mid-latitude ionosphere to the equator and may serve as a coupling agent between the auroral and the equatorial ionosphere.     

Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere

The ionization profiles produced by galactic cosmic rays in the Earth atmosphere and ionosphere are obtained on the basis of Monte Carlo simulations. Cascade processes in the atmosphere are simulated using CORSIKA 6.52 code with FLUKA 2006 and QGSJET II hadronic interaction subroutines. Proton induced showers are considered using a realistic atmospheric model (US Standard Atmosphere). The energy deposit from different components is taken into account, namely electromagnetic, hadron and muon components. The curvature of the atmosphere is considered in the computer code. On the basis of the computed ionization yield function the ion pair production rate in the atmosphere is obtained for different conditions and locations. The model is applicable to the entire atmosphere, from ground level to upper atmosphere. Several applications of the obtained results are discussed. The Monte Carlo simulation model considers nuclear interactions below the altitude of 35 km. It is compared with analytical–numerical electron production rate model. The latter model which takes into account the electromagnetic interactions above altitudes of 35 km has two main regions of application: above 50 km (thin target model) and between 35 and 50 km (intermediate target model). A good agreement between the CORSIKA results and analytical–numerical model results is found above altitude of 35 km.     

Variations in cosmic ray cutoff rigidities during the great geomagnetic storm of November 2004

A very strong interplanetary and magnetospheric disturbance observed on 7–13 November 2004 can be regarded as one of the strongest events during the entire period of space observations. In this paper we report on the studies of cosmic ray cutoff rigidity variations during 7–13 November 2004 showing how storm conditions can affect the direct cosmic ray access to the inner magnetosphere. Effective cutoff rigidities have been calculated for selected points on the ground by tracing trajectories of cosmic ray particles through the magnetospheric magnetic field of the “storm-oriented” Tsyganenko 2003 model. Cutoff rigidity variations have also been determined by the spectrographic global survey method on the basis of experimental data of the neutron monitor network. Relations between the calculated and experimental cutoff rigidities and the geomagnetic Dst-index and interplanetary parameters have been investigated. Correlation coefficients between the cutoff rigidities obtained by the trajectory tracing method and the spectrographic global survey method have been found to be in the limits 0.76–0.89 for all stations except the low-latitude station Tokyo (0.35). The most pronounced correlation has been revealed between the cutoff rigidities that exhibited a very large variation of ∼1–1.5 GV during the magnetic storm and the Dst index.     

Depth dependency of neutron density produced by cosmic rays in the lunar subsurface

Depth dependency of neutrons produced by cosmic rays (CRs) in the lunar subsurface was estimated using the three-dimensional Monte Carlo particle and heavy ion transport simulation code, PHITS, incorporating the latest high energy nuclear data, JENDL/HE-2007. The PHITS simulations of equilibrium neutron density profiles in the lunar subsurface were compared with the measurement by Apollo 17 Lunar Neutron Probe Experiment (LNPE). Our calculations reproduced the LNPE data except for the 350–400 mg/cm² region under the improved condition using the CR spectra model based on the latest observations, well-tested nuclear interaction models with systematic cross section data, and JENDL/HE-2007.     

Strong geomagnetic activity forecast by neural networks under dominant southern orientation of the interplanetary magnetic field

The paper deals with the relation of the southern orientation of the north–south component _Bz$B_z$ of the interplanetary magnetic field to geomagnetic activity (GA) and subsequently a method is suggested of using the found facts to forecast potentially dangerous high GA. We have found that on a day with very high GA hourly averages of _Bz$B_z$ with a negative sign occur at least 16 times in typical cases. Since it is very difficult to estimate the orientation of _Bz$B_z$ in the immediate vicinity of the Earth one day or even a few days in advance, we have suggested using a neural-network model, which assumes the worse of the possibilities to forecast the danger of high GA – the dominant southern orientation of the interplanetary magnetic field. The input quantities of the proposed model were information about X-ray flares, type II and IV radio bursts as well as information about coronal mass ejections (CME). In comparing the GA forecasts with observations, we obtain values of the Hanssen–Kuiper skill score ranging from 0.463 to 0.727, which are usual values for similar forecasts of space weather. The proposed model provides forecasts of potentially dangerous high geomagnetic activity should the interplanetary CME (ICME), the originator of geomagnetic storms, hit the Earth under the most unfavorable configuration of cosmic magnetic fields. We cannot know in advance whether the unfavorable configuration is going to occur or not; we just know that it will occur with the probability of 31%.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

Nightside 135.6 nm emission enhancements of mid/low latitudes during geomagnetic storms as observed by the ionosphere PhotoMeter (IPM) on the Chinese meteorological satellite FY-3D

IPM has detected nightside 135.6 nm emission enhancements over a wide latitude range, from the sub-auroral latitudes to the equatorial regions during geomagnetic storms. Our work, presented in this paper, uses the data of IPM to understand these 135.6 nm emission enhancements during of geomagnetic storms and studies the variations of total electron content (TEC) and the F₂ layer peak electron density (NmF2) in the region of enhanced emissions. Middle and low latitude emission enhancements are presented during several medium storms in 2018. The variations of both the integrated electron content (IEC) derived from the nighttime OI 135.6 nm emission by IPM and TEC from the International GNSS Service (IGS) relative to the daily mean of magnetically quiet days of per each latitude bin (30°≦geographic latitude < 40°, 15°≦geographic latitude < 30°, 0°≦geographic latitude < 15°, ?15°≦geographic latitude < 0°, ?30°≦geographic latitude < -15°, ?40°≦geographic latitude < -30°) are investigated and show that on magnetically storm day, IEC by IPM always increases, while TEC from IGC may increase or decrease. Even if both increase, the increase of IEC is greater than that of TEC. From the comparison of IEC and TEC during magnetic storms, it can be seen that the enhancement of the nighttime 135.6 nm emissions is not entirely due to the ionospheric change. The time of IEC enhancements at each latitude bin is in good agreement, which mainly corresponds to the main phase time of the geomagnetic storm event and lasts until the recovery phase. The available ground-based ionosonde stations provide the values of NmF2 which match the 135.6 nm emissions measured by IPM in space and time. The variations of NmF2 squared can characterize the variations of the OI 135.6 nm emissions caused by O⁺ ions and electrons radiative recombination. The study results show that the OI 135.6 nm emission enhancements caused by O⁺ ions and electrons radiative recombination (where NmF2 squared increases) are obviously a contribution to the measured 135.6 nm emission enhancements by IPM. The contribution accounts for at least one of all contributions to the measured 135.6 nm emission enhancements by IPM. However, where the NmF2 squared provided by ionosonde decrease or change little (where the OI 135.6 nm emissions cause by O⁺ ions and electrons radiative recombination also decrease or change little), the emission enhancements measured by IPM at storm-time appear to come from the contributions of other mechanisms, such as energetic neutral atoms precipitation, or the mutual neutralization emission (O+ + O-→2O + h? (135.6 nm)) which also occupies a certain proportion in 135.6 nm airglow emission at night.     

Analytical approach to cosmic ray ionization by nuclei with charge Z in the middle atmosphere – Distribution of galactic CR effects

The effects of galactic and solar cosmic rays (CR) in the middle atmosphere are considered in this work. A new analytical approach for CR ionization by protons and nuclei with charge Z in the lower ionosphere and middle atmosphere is developed in this paper. For this purpose the ionization losses (dE/dh) according to the Bohr–Bethe–Bloch formula for the energetic charged particles are approximated in three different energy intervals. More accurate expressions for energy decrease E(h) and electron production rate profiles q(h) are derived. The obtained formulas allow comparatively easy computer programming. The integrand in q(h) gives the possibility for application of adequate numerical methods – such as Romberg method or Gauss quadrature, for the solution of the mathematical problem. On this way the process of interaction of cosmic ray particles with the upper, middle and lower atmosphere will be described much more realistically. Computations for cosmic ray ionization in the middle atmosphere are made. The full CR composition is taken into account: protons, Helium (α-particles), light L, medium M, heavy H and very heavy VH group of nuclei.