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).     

Stratospheric polar vortex as a possible reason for temporal variations of solar activity and galactic cosmic ray effects on the lower atmosphere circulation

Possible reasons for the temporal instability of long-term effects of solar activity (SA) and galactic cosmic ray (GCR) variations on the lower atmosphere circulation were studied. It was shown that the detected earlier ∼60-year oscillations of the amplitude and sign of SA/GCR effects on the troposphere pressure at high and middle latitudes (Veretenenko and Ogurtsov, Adv.Space Res., 2012) are closely related to the state of a cyclonic vortex forming in the polar stratosphere. The intensity of the vortex was found to reveal a roughly 60-year periodicity affecting the evolution of the large-scale atmospheric circulation and the character of SA/GCR effects. An intensification of both Arctic anticyclones and mid-latitudinal cyclones associated with an increase of GCR fluxes at minima of the 11-year solar cycles is observed in the epochs of a strong polar vortex. In the epochs of a weak polar vortex SA/GCR effects on the development of baric systems at middle and high latitudes were found to change the sign. The results obtained provide evidence that the mechanism of solar activity and cosmic ray influences on the lower atmosphere circulation involves changes in the evolution of the stratospheric polar vortex.     

Possible influence of cosmic rays on climate through thunderstorm clouds

With decreasing of cosmic ray (CR) intensity caused by increasing of solar activity (SA) or in some short periods of Forbush-decreases, the intensity of secondary CR relativistic electrons decreases and the probability of formation of thunderstorm clouds and discharges between clouds or between clouds and ground is also expected to decrease. This will influence on weather and climate. In this case is very important to have more detail information on the atmospheric electric field distribution in the atmosphere, additional to information what gave now electric field sensors (EFS) only in about one point near the ground. We show that CR not only influenced on atmospheric electric field phenomenon, but can give practically continuous information on the atmospheric electric field distribution in the atmosphere. We extend our theory of CR atmospheric electric field effect on electron–photon, muon and neutron component including different multiplicities. We take into account that about 0.07 of neutron monitor counting rate caused by negative soft muons captured by lead nucleons and formed mesoatoms with generation of several MeV energy neutrons from lead. In this case the neutron monitor or neutron super-monitor works as analyzer that detects muons of only one, negative sign. It is very important because the atmospheric electric field effect have opposite signs for positive and negative muons that main part of this effect in the muon telescope or in ionization chamber is compensated and we can observe only small part of total effect of one sign muons. On the basis of our general theory of CR atmospheric electric field effects with taking into account of negative soft muon acceleration and deceleration in the Earth atmosphere (in dependence of direction and intensity of electric field) we discuss the possibility of existing this effect in CR neutron monitor counting rate and in different multiplicities and calculate the expected effects in dependence of atmospheric electric field distribution in the atmosphere. We show that the comparison of observed effects with theoretically expected will give important information on the value of atmospheric electric field and its distribution in the atmosphere. We consider also the possible influence of secondary relativistic electrons of CR and relativistic electrons precipitated from the Earth’s radiation belts on thunderstorms and lightnings, and through this – on climate change.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.     

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.     

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.     

Regional and temporal variability of solar activity and galactic cosmic ray effects on the lower atmosphere circulation

In this work we studied the spatial and temporal structure of long-term effects of solar activity (SA) and galactic cosmic ray (GCR) variations on the lower atmosphere circulation as well as possible reasons for the peculiarities of this structure. The study revealed a strong latitudinal and regional dependence of SA/GCR effects on pressure variations in the lower troposphere which seems to be determined by specific features of baric systems formed in different regions. The temporal structure of SA/GCR effects on the troposphere circulation at high and middle latitudes is characterized by a roughly 60-year periodicity which is apparently due to the epochs of the large-scale atmospheric circulation. It is suggested that a possible mechanism of long-term effects of solar activity and cosmic ray variations on the troposphere circulation involves changes in the evolution of the polar vortex in the stratosphere of high latitudes, as well as planetary frontal zones.     

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.     

Solar rhythms in the characteristics of the Arctic frontal zone in the North Atlantic

Long-term changes of the Arctic frontal zone characteristics near the south-eastern coasts of Greenland were considered, the NCEP/NCAR reanalysis data being used. It was found that in the cold half of the year the temperature gradients in the layer 1000–500 hPa in the region under study reveal strong ∼10-yr and ∼22-yr periodicities that seem to be related to solar activity cycles. The results obtained suggest the influence of solar activity and cosmic ray variations on the structure of the temperature field of the troposphere resulting in the changes of the temperature contrasts in the Arctic frontal zone that, in turn, may affect the intensity of cyclogenesis at middle latitudes. The detected effects seem to indicate an important part of frontal zones in the mechanism of solar activity and cosmic ray variation influence on the development of extratropical baric systems. It is suggested that the variations of the temperature gradients revealed in the Arctic frontal zone are due to the radiative forcing of cloudiness changes which may be associated with geomagnetic activity and cosmic ray variations.     

27-day variations of the galactic cosmic ray intensity and anisotropy in the minimum of the 23rd solar activity cycle

We study the 27-day variations of the solar wind velocity, galactic cosmic ray (GCR) intensity and anisotropy in the last minimum epoch of solar activity (2007–2009, A < 0). The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy (A27A) in the current minimum epoch of solar activity (2007–2009, A < 0) is lesser than in previous positive polarity period as it is expected from the drift theory. So, polarity dependence rule for the 27-day variation of the GCR anisotropy is fully kept. It is a universal principle for the amplitudes of the 27-day variation of the GCR anisotropy. At the same time, the average amplitude of the 27-day variation of the GCR intensity (A27I) remains at the same level as for previous minimum epoch 1995–1997 (A > 0) showing by the same token an violation of its polarity dependence rule established earlier. We assume that this phenomenon could be generally related with the well established 27-day variation of the solar wind velocity being in anti-correlation with the similar changes of the 27-day variation of the GCR intensity. Generally, a character of the heliolongitudinal asymmetry of spatial large-scale structure of the solar wind velocity (SWV) established in the recent minimum epoch, preferentially pronounces in the behavior of the 27-day variation of the GCR intensity than anisotropy. The formation of the 27-day variation of the GCR anisotropy preferentially takes place in a restricted disk like local vicinity in the helioequatorial region, whilst the 27-day variation of the GCR intensity is formed in the global three dimensional vicinity of the heliosphere.     

Rigidity spectrum of the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity

We study the temporal evolution of the power rigidity spectrum of the first (27 days) and the second (14 days) harmonics of the 27-day variation of the galactic cosmic ray intensity measured by neutron monitors in the period of 1965–2002. The rigidity spectrum of these variations can be approximated by a power law. We show the rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity have similar time profiles. These spectra are hard (γ ≈ 0.5 ± 0.1) and soft (γ ≈ 1.1 ± 0.2) during solar maximum and minimum activity, respectively. We ascribe this to the alternation of the sizes of the modulation regions responsible for the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity. Especially, the average radial sizes of the modulation regions of the 27-day variation and the heliolatitudinal extension of the heliolongitudinal asymmetry are smaller during solar minimum than during solar maximum. We show also, that the temporal changes of the power rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity are in a negative correlation with the changes of the rigidity spectrum of the corresponding 11-year variation.     

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

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

North–South asymmetry in cosmic ray fluxes as measured in the stratosphere and in selected solar wind parameters in the near-Earth space

North–South asymmetry in the cosmic ray fluxes as resulted from the long-term balloon measurements in the northern and southern polar stratosphere does not agree with that found from the neutron monitor data. In order to reveal possible sources of the observed asymmetry, selected interplanetary parameters were examined. North–South asymmetry relative to the heliospheric neutral sheet was considered for solar wind velocity, plasma density and some other solar plasma parameters. It is shown that North–South asymmetry of the solar wind velocity and plasma density depends on the Earth’s heliolatitude and the phase of the 11-year solar activity cycle. This may be relevant to the results of cosmic ray measurements in the stratosphere.     

Cosmic ray variations of solar origin in relation to human physiological state during the December 2006 solar extreme events

There is an increasing amount of evidence linking biological effects to solar and geomagnetic disturbances. A series of studies is published referring to the changes in human physiological responses at different levels of geomagnetic activity. In this study, the possible relation between the daily variations of cosmic ray intensity, measured by the Neutron Monitor at the Cosmic Ray Station of the University of Athens (http://cosray.phys.uoa.gr) and the average daily and hourly heart rate variations of persons, with no symptoms or hospital admission, monitored by Holter electrocardiogram, is considered. This work refers to a group of persons admitted to the cardiological clinic of the KAT Hospital in Athens during the time period from 4th to 24th December 2006 that is characterized by extreme solar and geomagnetic activity. A series of Forbush decreases started on 6th December and lasted until the end of the month and a great solar proton event causing a Ground Level Enhancement (GLE) of the cosmic ray intensity on 13th December occurred. A sudden decrease of the cosmic ray intensity on 15th December, when a geomagnetic storm was registered, was also recorded in Athens Neutron Monitor station (cut-off rigidity 8.53 GV) with amplitude of 4%. It is noticed that during geomagnetically quiet days the heart rate and the cosmic ray intensity variations are positively correlated. When intense cosmic ray variations, like Forbush decreases and relativistic proton events produced by strong solar phenomena occur, cosmic ray intensity and heart rate get minimum values and their variations, also, coincide. During these events the correlation coefficient of these two parameters changes and follows the behavior of the cosmic ray intensity variations. This is only a small part of an extended investigation, which has begun using data from the year 2002 and is still in progress.     

Lunar cosmic ray radiation environments during Luna and Lunar Reconnaissance Orbiter missions

The RV-2N-series instruments onboard Luna missions and the Cosmic Ray Telescope for the Effects of Radiation (CRaTER) instrument onboard Lunar Reconnaissance Orbiter (LRO) were designed to characterize the global lunar radiation environment and its biological impacts by measuring cosmic ray (CR) intensity. In this study, we have shown that the RV-2N-series instruments onboard of Russian Luna missions and the CRaTER reliably detect both background CRs and solar proton events (SPEs) in the lunar radiation environment using the proton intensity measured by the RV-2N-series onboard Luna missions out of the Russian Luna program for the exploration of the Moon (November 1970–August 1975) and the CR intensity on the Moon observed by the CRaTER (June 2009–March 2011). Those were compared with the CR intensities observed by neutron monitors (McMurdo, Thule, Oulu) on the Earth. The sunspot number is used as the index of solar activity (NOAA National Geophysical Data Center). As a result, the background CR intensities on the Moon turned out to have a good anti-correlation with the solar activity. We have also identified the proton intensity increasing events on the Moon which have the similar profiles to those observed by neutron monitors on the Earth. Most of these events show the significant increase of proton intensities in the lunar radiation environment when the SPEs associated with solar eruptions are verified. Therefore, most of the proton intensity increasing events are associated with the energetic solar particles in the lunar environment.     

The dynamic heliosphere,solar activity,and cosmic rays

This brief review addresses the relation between solar activity, cosmic ray variations and the dynamics of the heliosphere. The global features of the heliosphere influence what happens inside its boundaries on a variety of time-scales. Galactic and anomalous cosmic rays are the messengers that convey vital information on global heliospheric changes in the manner that they respond to these changes. By observing cosmic rays over a large range of energies at Earth, and with various space detectors, a better understanding is gained about space weather and climate. The causes of the cosmic ray variability are reviewed, with emphasis on the 11-year and 22-year cycles, step modulation, charge-sign dependent modulation and particle drifts. Advances in this field are selectively discussed in the context of what still are some of the important uncertainties and outstanding issues.     

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.     

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.     

Propagation modes of energetic charged particles in the heliosphere

Several years ago, the anisotropic diffusion and convective transport accompanied by adiabatic deceleration were considered as the principal means for cosmic ray propagation. Particles of relatively small energies (～ 1 MeV) can propagate along the force lines of the magnetic field without scattering at distances of several astronomical units in the quiet heliosphere. The theory describing the 11-year variation of galactic cosmic ray intensity and the propagation of solar cosmic rays was founded on this basis. However, the anomalies of the 11-year variation of galactic cosmic ray intensity in 1969–1971 revealed the necessity to take into account the influence of the general electromagnetic field of the heliosphere giving rise to a rapid magnetic drift of particles. The particles drift either from the magnetic axis to the ecliptic plane (in the cycle of 1969–1980) or in the opposite direction depending on the sign of the general magnetic field of the sun. The neutral layers along which the drift velocity is comparable to the particle velocity is of great significance. However, in the presence of sector structure, the time of particle propagation along the neutral layer from the boundary of the modulation region to the earth orbit is substantially increased. Thus a marked adiabatic deceleration is here possible. The time delay observed in the recovery of proton intensities at various energies can be explained in terms of a transient phase of the interplanetary field following the polarity reversal.