The radial distribution of gamma rays and cosmic rays in the outer Galaxy

The radial distribution of the high-energy (70 MeV-5 GeV) gamma-ray emissivity in the outer Milky Way is derived. The kinematics of HI are used to construct column-density maps in three galacto-centric distance ranges in the outer Galaxy. These maps are used in combination with COS-B gamma-ray data to determine gamma-ray emissivities in these distance ranges. A steep negative gradient of the emissivity for the 70 MeV-150 MeV energy range is found in the outer Galaxy. The emissivity for the 300 MeV-5 GeV range is found to be approximately constant (within 20%) and equal to the local value out to large (20 kpc) galacto-centric distances. These results imply a hardening of the gamma-ray spectrum with increrasing distance and for R > 16 kpc the spectrum is shown to be consistent with a π°-decay spectrum with the intensity expected from the local measurement of the cosmic-ray nuclei spectrum. The energy-dependent decrease is interpreted as a steep gradient in the cosmic-ray electron density and a near constancy of the nuclear component. The galactic origin of electrons with energies up to several hundreds of MeV is confirmed, while for cosmic-ray nuclei with energies of a few GeV either confinement in a large galactic halo or an extragalactic origin is suggested by the data.     

Study on solar sources and polar cap absorption events recorded in Antarctica

Particularly intense events occurred on the Sun in a period around minimum of solar activity during cycle 23. We investigated the characteristics of September 2005 and December 2006 events and the properties of the correlated observations of ionospheric absorption, obtained by a 30 MHz riometer installed at Mario Zucchelli Station (MZS-Antarctica), and of geomagnetic activity recorded at Scott Base (Antarctica). Solar events are studied using the characteristics of CMEs measured with SoHO/LASCO coronagraphs and the temporal evolution of solar energetic protons in different energy ranges measured by GOES 11 spacecraft.     

The radiation environment in the ISS-USLab measured by ALTEA: Spectra and relative nuclear abundances in the polar,equatorial and SAA regions

We present here the energy spectra relative to different geomagnetic regions as measured by the ALTEA (Anomalous Long Term Effects on Astronauts) detector in the International Space Station – USLab from August 2006 to July 2007.     

Launch and commissioning of the PAMELA experiment on board the Resurs-DK1 satellite

PAMELA is a satellite borne experiment designed to study with great accuracy cosmic rays of galactic, solar, and trapped nature in a wide energy range (protons: 80 MeV–700 GeV, electrons 50 MeV–400 GeV). Main objective is the study of the antimatter component: antiprotons (80 MeV–190 GeV), positrons (50 MeV–270 GeV) and search for antimatter (with a precision of the order of 10⁻⁸). The experiment, housed on board the Russian Resurs-DK1 satellite, was launched on June, 15th 2006 in a 350 × 600 km orbit with an inclination of 70°. The detector consists of a permanent magnet spectrometer core to provide rigidity and charge sign information, a Time-of-Flight system for velocity and charge information, a silicon–tungsten calorimeter and a neutron detector for lepton/hadron identification. An anticounter system is used off-line to reject false triggers coming from the satellite. In self-trigger mode the calorimeter, the neutron detector and a shower tail catcher are capable of an independent measure of the lepton (e⁺ + e⁻) component up to 2 TeV. In this work we focus on the first months of operations of the experiment during the commissioning phase.     

Evolution of magnetic fields and cosmic ray acceleration in supernova remnants

Observations show that the magnetic field in young supernova remnants (SNRs) is significantly stronger than can be expected from the compression of the circumstellar medium (CSM) by a factor of four expected for strong blast waves. Additionally, the polarization is mainly radial, which is also contrary to expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may help to explain these two observed features. They can increase the compression ratio to factors well over those of regular strong shocks by adding a relativistic plasma component to the pressure, and by draining the shock of energy when CRs escape from the region. The higher compression ratio will also allow for the contact discontinuity, which is subject to the Rayleigh–Taylor (R–T) instability, to reach much further out to the forward shock. This could create a preferred radial polarization of the magnetic field. With an Adaptive Mesh Refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three different configurations of the initial CSM magnetic field, and look at two different equations of state in order to look at the possible influence of a CR plasma component. The spectrum of CRs can be simulated using test particles, of which we also show some preliminary results that agree well with available analytical solutions.     

Mid 19th century minimum of galactic cosmic ray flux inferred from Ti in Allegan meteorite

Measurements of ⁴⁴Ti activity in meteorites show that the galactic cosmic ray (GCR) intensity has been declining in the interplanetary space during the past three centuries and has a component of cyclic variation, with periodicity of about 87 years [Taricco, C., Bhandari, N., Cane, D., et al. Galactic cosmic ray flux decline and periodicities in the interplanetary space during the last 3 centuries revealed by ⁴⁴Ti in meteorites. J. Geophys. Res. 111, A08102, 2006.]. In order to verify these results, we have measured ⁴⁴Ti activity in Allegan meteorite which fell in 1899 and in some other meteorites with better precision. The measurements confirm low cosmic ray flux and consequently high solar activity near the middle of 19th century.     

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.     

The charge-sign dependent effect in the solar modulation of cosmic rays

The centennial anniversary of the discovery of cosmic rays was in 2012. Since this discovery considerable progress has been made on several aspects related to galactic cosmic rays in the heliosphere. It is known that they encounter a turbulent solar wind with an imbedded heliospheric magnetic field when entering the Sun’s domain. This leads to significant global and temporal changes in their intensity inside the heliosphere, a process known as the solar modulation of cosmic rays. The prediction of a charge-sign dependent effect in solar modulation in the late 1970s and the confirmatory observational discoveries can also be considered as a milestone. A short review is given of these predictions based on theoretical and numerical modelling work, the observational confirmation and related issues.     

Long-term (50 years) measurements of cosmic ray fluxes in the atmosphere

Since the middle of 1957 till present time the group of researchers of P.N. Lebedev Physical Institute of the Russian Academy of Sciences has carried out the regular balloon borne measurements of charged particle fluxes in the atmosphere. The measurements are performed at polar (northern and southern) and middle latitudes and cover the interval of heights from the ground level up to 30–35 km. Standard detectors of particles (gas-discharged counters) have been used. More than 80,000 measurements of cosmic ray fluxes in the atmosphere have been performed to the present time. In the data analysis the geomagnetic field and the Earth’s atmosphere are used as cosmic ray spectrometers.     

Detection of the high energy component of Jovian electrons in Low Earth Orbit with the PAMELA experiment

The PAMELA experiment is devoted to the study of cosmic rays in Low Earth Orbit with an apparatus optimized to perform a precise determination of the galactic antimatter component of c.r. It is constituted by a number of detectors built around a permanent magnet spectrometer. PAMELA was launched in space on June 15th 2006 on board the Russian Resurs-DK1 satellite for a mission duration of 3 years. The characteristics of the detectors, the long lifetime and the orbit of the satellite, will allow to address several aspects of cosmic-ray physics. In this work we discuss the observational capabilities of PAMELA to detect the electron component above 50 MeV. The magnetic spectrometer allows a detailed measurement of the energy spectrum of electrons of galactic and Jovian origin. Long term measurements and correlations with Earth–Jupiter 13 months synodic period will allow to separate these two contributions and to measure the primary electron Jovian component, dominant in the 50–70 MeV energy range. With this technique it will also be possible to study the contribution to the electron spectrum of Jovian e⁻ reaccelerated up to 2 GeV at the Solar Wind Termination Shock.     

Applications and usage of the real-time Neutron Monitor Database

A high-time resolution Neutron Monitor Database (NMDB) has started to be realized in the frame of the Seventh Framework Programme of the European Commission. This database will include cosmic ray data from at least 18 neutron monitors distributed around the world and operated in real-time. The implementation of the NMDB will provide the opportunity for several research applications most of which will be realized in real-time mode. An important one will be the establishment of an Alert signal when dangerous solar cosmic ray particles are heading to the Earth, resulting into ground level enhancements effects registered by neutron monitors. Furthermore, on the basis of these events analysis, the mapping of all ground level enhancement features in near real-time mode will provide an overall picture of these phenomena and will be used as an input for the calculation of the ionization of the atmosphere. The latter will be useful together with other contributions to radiation dose calculations within the atmosphere at several altitudes and will reveal the absorbed doses during flights. Moreover, special algorithms for anisotropy and pitch angle distribution of solar cosmic rays, which have been developed over the years, will also be set online offering the advantage to give information about the conditions of the interplanetary space. All of the applications will serve the needs of the modern world which relies at space environment and will use the extensive network of neutron monitors as a multi-directional spectrographic detector. On top of which, the decreases of the cosmic ray intensity – known as Forbush decreases – will also be analyzed and a number of important parameters such as galactic cosmic ray anisotropy will be made available to the users of NMDB. A part of the NMDB project is also dedicated to the creation of a public outreach website with the scope to inform about cosmic rays and their possible effects on humans, technological systems and space-terrestrial environment. Therefore, NMDB will also stand as an informative gate on space research through neutron monitor’s data usage.     

On possible genesis of fractal dimensions in the turbulent pulsations of cosmic plasma – galactic cosmic rays – turbulent pulsation in planetary atmosphere system

The spectrum of turbulent pulsations induced in the atmosphere by the galactic cosmic rays is defined. A possible manifestation of genesis of fractal dimensions in the system of “spectrum of turbulent pulsations of cosmic plasma – galactic cosmic rays’ spectrum – spectrum of atmospheric turbulent pulsations” is analyzed.     

Modeling and experimental study of the 27-day variation of galactic cosmic-ray intensity for a solar wind velocity depending on heliolongitude

We develop a three-dimensional (3-D) model of the 27-day variation of galactic cosmic-ray (GCR) intensity with a spatial variation of the solar wind velocity. A consistent, divergence-free interplanetary magnetic field is derived by solving the corresponding Maxwell equations with a variable solar wind speed, which reproduces in situ observed experimental data for the time interval to be analyzed (24 August 2007–28 February 2008). We perform model calculations for the GCR intensity using the variable solar wind and the corresponding magnetic field. Results are compatible with experimental data; the correlation coefficient between our model predictions and observed 27-day GCR variation is 0.80 ± 0.05.     

Establishment of atmospheric weighted mean temperature model in the polar regions

Due to the special geographical location and extreme climate environment, the polar regions (Antarctic and Arctic) have an important impact on global climate change. Atmospheric weighted mean temperature (Tm) is a crucial parameter in the retrieval of precipitable water vapor (PWV) from the zenith wet delay (ZWD) of ground-based Global Navigation Satellite System (GNSS) signal propagation. In this paper, the correlation between weighted mean temperature and surface temperature (Ts) is studied firstly. It is shown that the correlation coefficients between Tm and Ts are 0.93 in the Antarctic and 0.94 in the Arctic. The linear regression Tm model and quadratic function Tm model of the Antarctic and the Arctic are established respectively using the radiosonde profiles of 12 stations in the Antarctic and 58 stations in the Arctic from 2008 to 2015. The accuracies of the linear regression Tm model, the quadratic function Tm model and GPT2w Tm model which is a state-of-the-art global Tm model are verified using the radiosonde profiles from 2016 to 2018 in the Antarctic and Arctic. Root Mean Square (RMS) errors of the linear regression Tm model, the quadratic function Tm model and GPT2w Tm model in the Antarctic are 3.07 K, 2.87 K and 4.32 K respectively, and those in the Arctic are 3.53 K, 3.38 K and 4.82 K, which indicates that the quadratic function Tm model has a higher accuracy compared to linear regression Tm model, and the accuracies of the two regional Tm models are better than that of GPT2w Tm model in the polar regions. In order to better evaluate the accuracy of Tm in the PWV retrieval, the PWV values of radiosondes are used for comparisons as the reference value. The RMS errors of PWV derived from the two Tm models are similar for 1.28 mm in the Antarctic and 1 mm in the Arctic respectively. In addition, the spatial and temporal variation characteristics of Tm are analyzed in the polar regions by spectral analysis of Tm data using fast Fourier transform. The results show that the Tm has obvious seasonality and annual periodicity in the polar regions, and the maximum difference between warm season and cold season is about 63 K. After comparing and analyzing the influences of latitude, longitude and elevation on the Tm in the polar regions, it is found that latitude and elevation have a greater influence on the Tm than the longitude. As the latitude and elevation increase, the Tm decreases, and vice versa in the polar regions.     

The scaler mode in the Pierre Auger Observatory to study heliospheric modulation of cosmic rays

The impact of the solar activity on the heliosphere has a strong influence on the modulation of the flux of low energy galactic cosmic rays arriving at Earth. Different instruments, such as neutron monitors or muon detectors, have been recording the variability of the cosmic ray flux at ground level for several decades. Although the Pierre Auger Observatory was designed to observe cosmic rays at the highest energies, it also records the count rates of low energy secondary particles (the scaler mode) for the self-calibration of its surface detector array. From observations using the scaler mode at the Pierre Auger Observatory, modulation of galactic cosmic rays due to solar transient activity has been observed (e.g., Forbush decreases). Due to the high total count rate coming from the combined area of its detectors, the Pierre Auger Observatory (its detectors have a total area greater than 16,000 m²) detects a flux of secondary particles of the order of ∼10⁸ counts per minute. Time variations of the cosmic ray flux related to the activity of the heliosphere can be determined with high accuracy. In this paper we briefly describe the scaler mode and analyze a Forbush decrease together with the interplanetary coronal mass ejection that originated it. The Auger scaler data are now publicly available.     

Fluxes and nuclear abundances of cosmic rays inside the magnetosphere using a transmission function approach

At 1 AU and outside the Earth’s magnetosphere, the relative abundances to protons for He (He/p), C (C/p) and Fe (Fe/p) nuclei were calculated using the observation data of AMS-01 (for p and He) and HEAO-3 (for C and Fe) above 0.8 GeV/nucleon. In addition, the transmission function (TF) for the GCR propagation inside the magnetosphere was evaluated using the IGRF and T96 (introduced by Tsyganenko and Stern) models to obtain permitted and forbidden trajectories inside the magnetosphere. The TF allowed one to derive the primary He-nuclei fluxes in the same geomagnetic regions of AMS-01 observations. These fluxes were found in good agreement with the observation data. Furthermore inside the magnetosphere in addition to the flux of helium, it allowed one to obtain those of the primary p, C, and Fe nuclei and the relative abundances of He, C and Fe nuclei to protons from the same observation data of AMS-01 and HEAO-3 above ≈0.8 GeV/nucleon. Up to a geomagnetic latitude of ≈45.84°, the relative isotopic abundances were found to depend on the mass number I_isot and, on average, range from a factor ≈2.31 up to ≈3.35 larger than those outside the magnetosphere at 1 AU. Thus, the magnetospheric isotopic/nuclear relative abundances differ from those inside the solar cavity and those in the interstellar space. The usage of the TF approach can allow one to determine the nuclear abundances in the magnetosphere at any geomagnetic latitude and, thus, any orbit, provided that the CR spectra are determined at 1 AU.     

Sawtooth substorms generated under conditions of the steadily high solar wind energy input into the magnetosphere: Relationship between PC, AL and ASYM indices

Periodicity in occurrence of magnetic disturbances in polar cap and auroral zone under conditions of steady and powerful solar wind influence on the magnetosphere is analyzed on the example of 9 storm events with distinctly expressed sawtooth substorms (N = 48). Relationships between the polar cap magnetic activity (PC-index), magnetic disturbances in the auroral zone (AL-index) and value of the ring current asymmetry (ASYM index) were examined within the intervals of the PC growth phase and the PC decline phase inherent to each substorm. It is shown that the substorm sudden onsets are always preceded by the PC growth and that the substorm development does not affect the PC growth rate. On achieving the disturbance maximum, the PC and AL indices are simultaneously fall down to the level preceding the substorm, so that the higher the substorm intensity, the larger is the AL and PC drop in the decline phase. The ASYM index increases and decreases in conformity with the PC and AL behavior, the correlation between ASYM and PC being better than between ASYM and AL. Level of the solar wind energy input into the magnetosphere determines periodicity and intensity of disturbances: the higher the coupling function E_KL, the higher is substorm intensity and shorter is substorm length. Taking into account the permanently high level of auroral activity and inconsistency of aurora behavior and magnetic onsets during sawtooth substorms, the conclusion is made that auroral ionosphere conductivity is typically high and ensures an extremely high intensity of field-aligned currents in R1 FAC system. The periodicity of sawtooth substorms is determined by recurrent depletions and restorations of R1 currents, which are responsible for coordinated variations of magnetic activity in the polar cap and auroral zone.     

A statistical analysis of low frequency geomagnetic field pulsations at two Antarctic geomagnetic observatories in the polar cap region

The aim of this study is to investigate the characteristics of low frequency (∼0.5–5 mHz) geomagnetic field fluctuations as recorded at two Antarctic stations within the polar cap: the Italian observatory Mario Zucchelli Station (TNB) and the French–Italian observatory Dome C (DMC) in order to investigate the spatial extension and propagation characteristics of the phenomena observed at very high latitude. The stations have approximately the same geographic latitude, but a very different corrected geomagnetic latitude, being DMC close to the geomagnetic pole and TNB closer to the auroral oval.     

Studies of polar mesosphere summer echoes with the EISCAT VHF and UHF radars: Information contained in the spectral shape

The nature of PMSE in the VHF and UHF frequency range is considered taking into account the shape of corresponding Doppler spectra. Assuming a turbulence-based model of PMSE it is argued that for cases where a VHF radar detects strong PMSE, the UHF radar could either detect enhanced coherent scattering caused by the same physical process as in the VHF (i.e., turbulence with large charged ice particles), there could be incoherent scattering modified by the charged ice particles, or there could be a mixture of both. In order to distinguish these cases a simple but robust method is introduced to characterize the shape of the Doppler spectra derived from observations at both frequencies. Spectral shapes are quantified with one simple fitting parameter of a generalized fit to the autocorrelation function (=Fourier transform of the Doppler spectrum). This parameter takes a value of 1 for a Lorentzian spectrum indicative of pure incoherent scatter from the D-region, a value of 2 for coherent scatter owing to turbulence, and a value of less than 1 for incoherent scatter modified by the presence of charged aerosol particles. This method is applicable to observations at altitudes between ∼70 and ∼90 km. Simultaneous observations with the EISCAT VHF and UHF radar are presented in which all three cases mentioned above are identified. For the case of incoherent scatter modified by the presence of charged aerosol particles we quantify the radius of the involved ice particles to exceed ∼5 nm. Most importantly, however, for the case where the UHF-signal exceeded the incoherent scatter signal significantly, the spectrum revealed a clear Gaussian shape indicative of a coherent scattering process with identical spectral width as for the VHF-observations. This finding gives strong support that both echoes are created by the same turbulence-based mechanism and not by different mechanisms as speculated by several previous authors.