Effect of random nature of cosmic ray sources – Supernova remnants – on cosmic ray intensity fluctuations,anisotropy, and electron energy spectrum

The random nature of sources (the supernova remnants) leads to the fluctuations of cosmic ray intensity in space and time. We calculate the expected fluctuations in a flat-halo diffusion model for particles with energies from 0.1 to 10³ TeV. The data on energy spectra and anisotropy of very high energy protons, nuclei and electrons, and the astronomical data on supernova remnants, the potential sources of cosmic rays, are used to constrain the value of the cosmic-ray diffusion coefficient and its dependence on energy.     

Observations of non-migrating tides and ionospheric perturbations of O(D) airglow by ISUAL instrument

In this paper, we present the spatial variations of O(¹D) airglow observed by the ISUAL (Imager of Sprites and Upper Atmospheric Lightning) instrument on board the FORMOSAT-2 satellite. With a CCD camera and a 630 nm filter, ISUAL can measure global atmospheric emissions lying between the heights of 80 and 300 km. In days of 3–6 September 2008 and 25–27 February 2009, ISUAL has measured the emissions of O(¹D) airglow with results showing strong longitudinal peak-3 and peak-4 structures. The Lomb-Scargle analyses for these two cases show periods of longitudes of 120° and 90° supporting the DE2 and DE3 non-migrating tides. The 630 nm emissions are enhanced in equatorial regions and are lying along the equator. Over Africa its intensity can sometimes increase up to 80% relative to other longitudes. The perturbation is so strong that non-migrating tides are erased. A case of bimodal distribution with strong emissions at latitudes in equator and mid-latitude in geographic coordinates was observed.     

Study of the 27-day variations of the galactic cosmic ray intensity and anisotropy

We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.     

Prediction of expected global climate change by forecasting of galactic cosmic ray intensity time variation in near future based on solar magnetic field data

A method of prediction of expected part of global climate change caused by cosmic ray (CR) by forecasting of galactic cosmic ray intensity time variation in near future based on solar activity data prediction and determined parameters of convection-diffusion and drift mechanisms is presented. This gave possibility to make prediction of expected part of global climate change, caused by long-term cosmic ray intensity variation. In this paper, we use the model of cosmic ray modulation in the Heliosphere, which considers a relation between long-term cosmic ray variations with parameters of the solar magnetic field. The later now can be predicted with good accuracy. By using this prediction, the expected cosmic ray variations in the near Earth space also can be estimated with a good accuracy. It is shown that there are two possibilities: (1) to predict cosmic ray intensity for 1–6 months by using a delay of long-term cosmic ray variations relatively to effects of the solar activity and (2) to predict cosmic ray intensity for the next solar cycle. For the second case, the prediction of the global solar magnetic field characteristics is crucial. For both cases, reliable long-term cosmic ray and solar activity data as well as solar magnetic field are necessary. For solar magnetic field, we used results of two magnetographs (from Stanford and Kitt Peak Observatories). The obtained forecasting of long-term cosmic ray intensity variation we use for estimation of the part of global climate change caused by cosmic ray intensity changing (influenced on global cloudiness covering).     

Long-term variations of the surface pressure in the North Atlantic and possible association with solar activity and galactic cosmic rays

Long-term variations of the surface pressure in the North Atlantic for the period 1874–1995 (Mean Sea Level Pressure archive, Climatic Research Unit, UK) were compared with indices of solar and geomagnetic activity and the galactic cosmic ray (GCR) variations characterized by the concentration of the cosmogenic isotope ¹⁰Be. A periodicity of ∼80 yrs close to the Gleissberg cycle in the intensity of the 11-yr solar cycles was found in the pressure variations at middle latitudes (45–65°N) in the cold half of the year, which is the period of intensive cyclogenesis. It was shown that a long-term increase of pressure in this region coincided with a secular rise of solar/geomagnetic activity which was accompanied by a decrease in GCR intensity. Long-term decreases of pressure were observed during the periods of low (or decreasing) intensities of sunspot cycles. Similar features were also found in the spectral characteristics of geomagnetic activity indices, GCR intensity and pressure at middle latitudes on the quasi-decadal time scale. Effects of solar activity/GCR variations on the surface pressure seem to be more pronounced in the North Atlantic zone of intensive cyclogenesis (near the eastern coasts of North America). The results obtained suggest possible links between long-term variations in cyclonic activity at extratropical latitudes of the North Atlantic and solar activity/GCR variations on the time scales from ∼10 to ∼100 yrs.     

Application of HF Doppler measurements for the investigation of internal atmospheric waves in the ionosphere

We present results of the spectral analysis of data series of Doppler frequency shifted signals reflected from the ionosphere, using experimental data received at Kazan University, Russia. Spectra of variations with periods from 1 min to 60 days have been calculated and analyzed for different scales of periods. The power spectral density for spring and winter differs by a factor of 3–4. Local maxima of variation amplitude are detected, which are statistically significant. The periods of these amplitude increases range from 6 to 12 min for winter, and from 24 to 48 min for autumn. Properties of spectra for variations with the periods of 1–72 h have been analyzed. The maximum of variation intensity for all seasons and frequencies corresponds to the period of 24 h. Spectra of variations with periods from 3 to 60 days have been calculated. The maxima periods of power spectral density have been detected by the MUSIC method for the high spectral resolution. The detected periods correspond to planetary wave periods. Analysis of spectra for days with different level of geomagnetic activity shows that the intensity of variations for days with a high level of geomagnetic activity is higher.     

Delineation of possible influence of solar variability and galactic cosmic rays on terrestrial climate parameters

The present paper has investigated the associations of solar activity (SA), represented by total solar irradiance (TSI), galactic cosmic rays (GCR) and terrestrial climate parameters in particular the global cloudiness and global surface temperature. To that end, we have analysed thirty five years (1983–2018) data of these parameters and have applied the Granger-causality test in order to assess whether there is any potential predictability power of one indicator to the other. The correlations among the involved parameters are tested using Vector Auto Regression (VAR) model and variance decomposition method. As a result of the above analysis, we have found that the TSI is an important factor and has contributed about 8.77 ± 0.42% in the cosmic ray intensity variations. In case of cloud cover variations, the other three parameters (TSI, cosmic ray and global surface temperature) have played a significant role. Further, the TSI changes have contributed 1.68 ± 0.03% fluctuations in the variance of the cloud cover while the cosmic ray intensity and global surface temperature have contributed about 4.89 ± 0.08% and 10.87 ± 1.41% respectively. In case of the global surface temperature anomaly both TSI and cloud covers have contributed about 5.07 ± 0.47% and 14.42 ± 2.13% fluctuations respectively. Additionally, we have also assessed the impact of internal climate oscillations like multivariate ENSO index (MEI), north Atlantic oscillations (NAO) and quasi biennial oscillations (QBO) on cloud cover variations. The contribution of these internal oscillations e.g. ENSO, NAO and QBO in cloud cover variation were reported as 7.48 ± 1.02%, 5.51 ± 0.16% and 1.36 ± 0.43% respectively.     

Relativistic solar protons in the ground level event of 23 February 1956: New study

In this study we applied again to the outstanding solar particle event of 23 February 1956, the largest one in the entire history of observations of solar cosmic rays. Due to significant improvement of the analysis/modeling techniques and new understanding of physical processes in the solar atmosphere and interplanetary space, a possibility arises to interpret the old data in the light of modern concept of multiple particle acceleration at/near the Sun. In our new analysis the data of available then neutron monitors and muon telescopes are used. The technique of the analysis includes: (a) calculation of asymptotic cones of ground-based detectors; (b) modeling of cosmic ray detector responses at variable parameters of the flux of solar relativistic protons; (c) determination of primary solar proton parameters outside magnetosphere by comparison of computed responses with observations. Certain evidence was obtained that the flux of relativistic solar protons consisted of two distinct components: prompt and delayed ones. The prompt component with exponential energy spectrum caused a giant impulse-like increase at a number of European cosmic ray stations. The delayed component had a power-law spectrum and was a cause of gradual increase at cosmic ray stations in the North American region. A numerical simulation of the proton acceleration in the vicinity of the magnetic reconnection region brings to the proton spectrum with exponential dependence on energy. This agrees with observational data for the prompt component. It is also shown that the huge increase in ∼5000% on neutron monitors was due to the prompt component only with the exponential proton spectrum. The power-law spectrum of comparable intensity gave considerably smaller effect.     

The structure of the solar cycle maximum phase in the galactic cosmic ray 11-year variation

Two phenomena connected with the maximum phase of the 11-year solar cycle in the galactic cosmic ray intensity – the change in the energy dependence of the intensity variations and the double-peak structure in the intensity modulation time profile – are considered for the last five solar cycles (Nos. 19–23). The distinct 22-year cycle in the magnitude of the so called energy hysteresis is observed.The periods of the solar cycle maximum phase in the galactic cosmic ray intensity, characterized by the specific energy dependence of the intensity, are estimated. It is found that the double-peak structures belonging to the solar cycle maximum phase and those around it are very similar both in the amplitude and in its energy dependence.     

基于小波与交叉小波分析的太阳黑子与宇宙线相关性研究

利用小波分析和交叉小波分析方法, 根据太阳黑子数以及Huancayo和Climax两个测站的月均宇宙线数据, 分析了两个测站的月均宇宙线周期变化, 同时利用太阳黑子数R₁₂对Climax站宇宙线流量进行预测研究. 小波分析结果表明, 太阳黑子与宇宙线除存在显著的11年周期外, 太阳活动高年期间还存在1～6个月尺度的周期特性, 在第22太阳周活动高年时还出现了6～8和1～22个月的变化周期; 交叉小波分析结果表明, 在130个月左右的周期上宇宙线与太阳黑子具有显著的负相关性, 并且宇宙线的变化滞后太阳黑子约8个月; 分别采用预测时刻和8个月前的太阳黑子数, 预测相对误差为3.8912%和3.2386%. 本文方法同样适用于估算其他空间天气参量之间的周期和相关性, 提高各种空间天气参量的预测或预报精度.      

Doppler observations of infrasonic waves of meteorological origin at ionospheric heights

Ionospheric effects of meteorological origin observed by the continuous HF Doppler sounder over the Czech Republic are reported in this paper. We focused on detection of waves of periods 1–10 min. We discuss the influence of dynamics and intensity of active weather systems on the occurrence of short period waves and dependence of the observed ionospheric effects on the height of reflection of the sounding radio wave. We observed 3–5 min waves during a severe weather event in summer and 2.5–4 min waves during a severe weather event in winter. We excluded possible geomagnetic origin of these oscillations by the analysis of fluctuations of the local geomagnetic field. In eight cases of 10, wave activity in the analysed period range was not significantly increased comparing to quiet days. The intensity of weather systems as well as the location of potential sources of waves towards the points of HF Doppler shift observation influence significantly the occurrence of infrasonic waves in the ionosphere. The results in Central Europe differ considerably from those previously obtained in North America. As a possible reason, we discuss different intensity and dynamics of weather systems in both regions.     

Estimation of long-term cosmic ray intensity variation in near future and prediction of their contribution in expected global climate change

On the basis of results obtained in our paper [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], we determine: the dimension of the Heliosphere (modulation region), radial diffusion coefficient and other parameters of convection–diffusion; drift mechanisms of long-term variations of cosmic ray (CR) dependence on particle energy; level of solar activity (SA); and generally, the solar magnetic field. We obtain this important information on the basis of CR and SA data in the past, taking into account the theory of convection–diffusion and global drift modulation of galactic CR in the Heliosphere. By using these results and other regularly published predictions of expected SA variation in the near future, as well as predictions of the next SA cycle, we may make predictions of long-term cosmic ray intensity variation expected in the near future (up to 10–12 years). In [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], properties of connections between long-term variation in CR intensity and some part of a global climate change were estimated, controlled by solar activity through CR. We show that in this way we may make predictions of some part of a global climate change expected in the near future (up to 10–12 years and maybe more, depending upon the period during which definite predictions of SA can be made), controlled by solar activity through CR. In this case, estimations of expected long-term changes in the planetary distribution of cutoff rigidities, which also influence CR intensity, as well as CR-influenced effects on global climate variation, become important.     

Research by the Tasmanian cosmic ray group during the International Geophysical Year

Systematic recording of the cosmic radiation commenced in Hobart in 1946 and at Mawson in Antarctica in 1955, making these two of the longest running cosmic ray observatories in the world. For the IGY, observations were also made at a sub-Antarctic island and near the equator, and an airborne survey of the nucleonic component was made from Geomagnetic Latitude −60°, south of Australia, to Japan and back. At Hobart there were neutron monitors, vertical and inclined muon telescopes, an ionization chamber, and two muon telescopes at ∼40 m of water equivalent underground. The research based on these and other observations determined the energy dependence of the Forbush and 11-year variations and concentrated, in particular, on understanding the anisotropic nature of galactic cosmic rays up to 150 GeV; the anisotropies in the onset phase of Forbush decreases; and the anisotropies in solar cosmic ray events. An investigation was initiated to calculate the trajectories and cutoff rigidities of cosmic rays in a high order simulation of the geomagnetic field. This was completed in 1959–60.     

23rd solar cycle: Solar activity parameters and associated Forbush decreases

We report Forbush decreases (FD) in cosmic ray intensity from January 1996 to December 2008, the whole Solar Cycle 23rd. Statistical analysis is done for only 152 events for which associated solar flare position, flare classes, and Coronal Mass Ejections (CME) speed are given. We applied FD parameters taken from the Forbush Effects and Interplanetary Disturbances databases maintained by the Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radiowave Propagation (IZMIRAN), obtained by processing the data of the worldwide neutron monitor network using the global survey method (GSM) (A. Belov et al., 2018). For the said number of events, we examine their effect on interplanetary space and the decrease of the galactic cosmic rays (GCR) near Earth. We found that the 11–20° latitudinal belt shows more FD- associated flare events than the other latitudinal belts, and on this belt, the Southern hemisphere is more active. The results reveal that FDs and solar flares are well correlated. Statistical analysis is carried out for the magnitude of the CR decrease with solar and geomagnetic parameters.     

低平流层温度对银河系宇宙线强度变化的响应

选用芬兰中子堆监测站从lop年7月至1989年3月的宇宙线强度月平均曲线数据,和在相对应期间我国境内海拉尔等4个高空气象观测站的30hPa和100hPa等位势高度上常规气温观测月平均资料,经统计分析,取得初步结果:低平流层气温与银河系宇宙线辐射强度之间呈现负相关的特性;与银河系宇宙线辐射强度峰区和谷区相对应会在低平流层中引起幅值为2-3℃的温度变化;不同纬度观测站所得观测数据表明,上述增温幅值随纬度降低而减少。      

Measurements of the radiation quality factor Q at aviation altitudes during solar minimum (2006–2008)

In radiation protection, the Q-factor has been defined to describe the biological effectiveness of the energy deposition or absorbed dose to humans in the mixed radiation fields at aviation altitudes. This particular radiation field is generated by the interactions of primary cosmic particles with the atoms of the constituents of the Earth’s atmosphere. Thus the intensity, characterized by the ambient dose equivalent rate H∗(10), depends on the flight altitude and the energy spectra of the particles, mainly protons and alpha particles, impinging on the atmosphere. These charged cosmic projectiles are deflected both by the interplanetary and the Earth’s magnetic field such that the corresponding energy spectra are modulated by these fields. The solar minimum is a time period of particular interest since the interplanetary magnetic field is weakest within the 11-year solar cycle and the dose rates at aviation altitudes reach their maximum due to the reduced shielding of galactic cosmic radiation. For this reason, the German Aerospace Center (DLR) performed repeated dosimetric on-board measurements in cooperation with several German airlines during the past solar minimum from March 2006 to August 2008. The Q-factors measured with a TEPC range from 1.98 at the equator to 2.60 in the polar region.     

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.     

Observed solar cycle modulation of galactic cosmic rays over a range of rigidities

We have studied the long-term, steady-state, solar cycle modulation of galactic cosmic ray intensity for seven cycles (17–23). Our analysis is based on the data obtained with a variety of detectors on earth (neutron monitors of the global network and muon detectors) as well as telescopes flown on high altitude balloons and on-board near-earth satellites. The median rigidity of response for these detectors to galactic cosmic ray spectrum lies in the range 1–70 GV. We correlate cosmic ray data to sunspot numbers, Ap, solar wind bulk speed (V), magnetic field (B), as well as to the cycle maximum (M), minimum (m), and the epochs of the solar polar field reversals. This enables us to derive the rigidity dependence of observations, and helps us to define the characteristics of the modulation function in the heliosphere.     

The IGY and beyond: A brief history of ground-based cosmic-ray detectors

The state of art of ground-based cosmic-ray research from its discovery to present is reviewed. After discovery of cosmic rays by Hess in 1912, the nature of the primary and secondary radiation was established from recordings by a variety of instruments, sensitive to various components of cosmic rays and operated at different latitudes, longitudes and altitudes, including instruments carried by balloons. The IGY formalized international co-operation and coordinated study of cosmic rays, which is vital for meaningful interpretation of cosmic-ray data. Data collected at different geographic locations require an effective cutoff rigidity as a data ordering parameter. This parameter is obtained from tracing trajectories of primary cosmic rays in the Earth’s magnetic field. After 50 years the world’s neutron monitor network remains still the backbone for studying intensity variations of primary cosmic rays in the rigidity ranges between 1 and 15 GV, associated with transport and with transient events. Also the penetrating muon and neutrino components of secondary cosmic rays have a long history of recording and fundamental problem investigations. Valuable data about composition and spectrum of primary cosmic rays in ever increasing high-energy regions have been obtained during the years of investigations with various configurations and types of extensive air shower detectors. The culture of personal involvement of the physicist in carrying out experiments and data acquisition characterized the continued vitality of cosmic-ray investigations ranging from its atmospheric, geomagnetic and heliospheric transport through to its solar and astrophysical origins.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.