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.     

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

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

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.     

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.     

The physical mechanism of the solar variability influence on electrical and climatic characteristics of the troposphere

Possible mechanisms of solar–climatic connections, which may be of importance over short and long time intervals, are discussed. The variations of energetic balance of Earth’s climatic system for the last 50 years are estimated. It is ascertained that the imbalance between the flux of solar energy that comes to the Earth and radiates to space is of 0.1% for the last ten years. The significance is analyzed for the possible influence of variations of solar constant upon the energetic balance of the atmosphere.The physical mechanism of the influence of solar activity on climatic characteristics and the atmospheric circulation is suggested and theoretically substantiated. The mechanism is based on the redistribution in lower-troposphere of condensation nuclei by the vertical electric field. This electric field is determined by the ionosphere–Earth electric potential, which in the Polar Regions is controlled not only by tropical thunderstorms and by the galactic cosmic-ray intensity but also by solar cosmic-ray fluxes. The height redistribution in the atmosphere of condensation nuclei with a change of the electric field of the atmosphere is accompanied by a change in total latent heat (phase transition of water vapor), by changes in radiation balance, and by subsequent changes of the thermobaric field of troposphere. The results of analysis of thermobaric field variations for the periods of invasion of abnormally powerful solar cosmic ray fluxes and magnetic storms confirm the reality of manifestation of heliogeophysical disturbances.     

The influence of interplanetary shocks on solar protons measured in the stratosphere.

Since the beginning of the 22nd solar cycle twenty solar proton events were observed by the regular balloon measurements of cosmic rays. Temporal changes of intensities and energy spectra of solar protons with energy 100-500 MeV were obtained. The strong influence of interplanetary shock waves on the proton flux characteristics near the Earth was observed. Possible effects of solar proton transport in the vicinity of shock fronts are discussed to explain the observational data.     

Dynamics of solar cosmic ray events: Processes at large heliocentric distances (⪢1 AU)

Observations of solar cosmic ray events far from the sun (?1 AU) became possible after the launch of Pioneer 10 in 1972. Four spacecraft have now travelled beyond the orbit of Jupiter - Pioneer 10/11 and Voyager 1/2 — and are producing a growing body of distant observations of solar cosmic ray events. Initial studies using Pioneer 10/11 data out to ～6 AU interpreted flare particle observations in terms of a diffusion model, including the effects of convection and adiabatic energy loss. This model enjoyed general success in explaining the time-intensity profiles in cases where the spacecraft connection longitude at the sun did not change significantly with time. The results implied that the radial diffusion coefficient (K_r) increased slowly with distance over that radial range. More recent results at larger distances imply that K_r may begin to decrease beyond ～5 AU. It is not yet clear whether the standard diffusion model will be adequate to explain solar events well beyond 5 AU. The fact that large events at very large distances can last up to two solar rotations implies that solar wind stream structure will also play a role in the event dynamics. In general, however, observations at large distances offer perhaps the best hope of separating interplanetary propagation effects from coronal storage and propagation effects which frequently dominate observed event profiles at 1 AU.     

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.     

Effects of solar modulation on the cosmic ray positron fraction

We implemented a 2D Monte Carlo model to simulate the solar modulation of galactic cosmic rays. The model is based on the Parker’s transport equation which contains diffusion, convection, particle drift and energy loss. Following the evolution in time of the solar activity, we are able to modulate a local interstellar spectrum (LIS), that we assumed isotropic beyond the termination shock, down to the Earth position inside the heliosphere. In this work we focused our attention to the cosmic ray positron fraction at energy below ∼10 GeV, showing how the particle drift processes could explain different results for AMS-01 and PAMELA. We compare our modulated spectra with observations at Earth, and then make a prediction of the cosmic ray positron fraction for the AMS-02 experiment.     

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.     

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

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.     

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 atmospheric cosmic- and solar energetic particle radiation environment at aircraft altitudes

Galactic cosmic rays interact with the solar wind, the earth's magnetic field and its atmosphere to produce hadron, lepton and photon fields at aircraft altitudes. In addition to cosmic rays, energetic particles generated by solar activity bombard the earth from time to time. These particles, while less energetic than cosmic rays, also produce radiation fields at aircraft altitudes which have qualitatively the same properties as atmospheric cosmic rays. We have used a code based on transport theory to calculate atmospheric cosmic-ray quantities and compared them with experimental data. Agreement with these data is seen to be good. We have then used this code to calculate equivalent doses to aircraft crews. We have also used the code to calculate radiation doses from several large solar energetic particle events which took place in 1989, including the very large event that occurred on September 29^th and 30^th of that year. The spectra incident on the atmosphere were determined assuming diffusive shock theory.     

50–200 GV cosmic rays at various heliolatitudes

We consider regular motion of 50 – 200 GV particles in a large-scale interplanetary magnetic field model which contains a wavy neutral sheet responsible for the sector-structure. Numerical calculations based upon energy losses along various trajectories are carried out to obtain the predicted omni-directional density and anisotropy of cosmic rays at various solar latitudes. A marked difference is found between odd and even solar cycles. The post-1969 field configuration gives small radial and large latitudinal gradient: cosmic ray density increases toward the poles. The latitudinal gradient turns out smaller and of opposite sense for the pre-1969 epoch. Anisotropy changes dramatically as we move off the solar equator: corotation appears to be restricted to low latitudes.     

Modeling the solar cosmic ray event of 13 December 2006 using ground level neutron monitor data

In order to understand the physics under extreme solar conditions such as those producing ground level enhancements of solar cosmic rays, it is important to use accurate and reliable models. The NM-BANGLE Model is a new cosmic ray model which couples primary solar cosmic rays at the top of the Earth’s atmosphere with the secondary ones detected at ground level by neutron monitors during GLEs. This model calculates the evolution of several GLE parameters such as the solar cosmic ray spectrum, anisotropy and particle flux distribution, revealing crucial information on the energetic particle propagation and distribution. The total output of the NM-BANGLE Model is a multi-dimensional GLE picture that gives an important contribution to revealing the characteristics of solar energetic particle events recorded at ground level. In this work, the results of the NM-BANGLE Model application to the recent GLE of 13 December 2006 are presented and discussed. Moreover, a comparison with the extreme event of 20 January 2005 (GLE69) has been realized.     

The atmospheric radiation response to solar-particle-events.

High-energy solar particles, produced in association with solar flares and coronal mass ejections, occasionally bombard the earth's atmosphere. resulting in radiation intensities additional to the background cosmic radiation. Access of these particles to the earth's vicinity during times of geomagnetic disturbances are not adequately described by using static geomagnetic field models. These solar fluxes are also often distributed non uniformly in space, so that fluxes measured by satellites obtained at great distances from the earth and which sample large volumes of space around the earth cannot be used to predict fluxes locally at the earth's surface. We present here a method which uses the ground-level neutron monitor counting rates as adjoint sources of the flux in the atmosphere immediately above them to obtain solar-particle effective dose rates as a function of position over the earth's surface. We have applied this approach to the large September 29-30, 1989 ground-level event (designated GLE 42) to obtain the magnitude and distribution of the solar-particle effective dose rate from an atypically large event. The results of these calculations clearly show the effect of the softer particle spectra associated with solar particle events, as compared with galactic cosmic rays, results in a greater sensitivity to the geomagnetic field, and, unlike cosmic rays, the near-absence of a "knee" near 60 degrees geomagnetic latitude.     

Peculiarities of galactic cosmic ray Forbush effects during October–November 2003

Features of two successive Forbush effects of the galactic cosmic ray intensity in October–November 2003 have been studied based on the neutron monitors data. The rigidity spectrum of the galactic cosmic ray intensity in the course of the first Forbush effect (22–27 October) is gradually hardening, while the rigidity spectrum of the second Forbush effect (28 October–10 November) from the starting moment is very hard. As far, the energy range of the turbulence of the interplanetary magnetic field is in general responsible for the diffusion of galactic cosmic ray particles of the energy 5–50 GeV (to which neutron monitors are sensitive), we postulate that the gradually hardening (from day to day) of the rigidity spectrum of the first Forbush effect is associated with the enhancement of the power spectral density in the energy range of the interplanetary magnetic field turbulence caused by the large scale irregularities generated due to the interaction of the extending high speed disturbances with the background solar wind. The very hard rigidity spectrum (from the starting moment) of the second Forbush effect is generally associated with the well established new structure of the energy range of the interplanetary magnetic field turbulence enriched by the already created large scale irregularities. The gradually softening of the rigidity spectrum during the recovery phase of the second Forbush effect confirms that the disturbed interplanetary magnetic field turbulence step by step returns to the initial state.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.