Influence of the Precipitating Energetic Particles on Atmospheric Chemistry and Climate

We evaluate the influence of the galactic cosmic rays (GCR), solar proton events (SPE), and energetic electron precipitation (EEP) on chemical composition of the atmosphere, dynamics, and climate using the chemistry-climate model SOCOL. We have carried out two 46-year long runs. The reference run is driven by a widely employed forcing set and, for the experiment run, we have included additional sources of NO _x and HO _x caused by all considered energetic particles. The results show that the effects of the GCR, SPE, and EEP fluxes on the chemical composition are most pronounced in the polar mesosphere and upper stratosphere; however, they are also detectable and statistically significant in the lower atmosphere consisting of an ozone increase up to 3?% in the troposphere and ozone depletion up to 8?% in the middle stratosphere. The thermal effect of the ozone depletion in the stratosphere propagates down, leading to a warming by up to 1?K averaged over 46?years over Europe during the winter season. Our results suggest that the energetic particles are able to affect atmospheric chemical composition, dynamics, and climate.     

The effect of model assumptions on computations of cosmic ray induced ionization in the atmosphere

The longitudinal profile of atmospheric cascades is sensitive to the energy, mass of the primary particle and to atmospheric state. In this work are compared ionization yield functions Y for winter, summer and US standard profiles of Earth's atmosphere. The various profiles are obtained on the basis of CORSIKA 6.52 code simulations using FLUKA 2006 and QGSJET II hadronic interaction subroutines. The energy deposit of proton induced cascade processes in the atmosphere is calculated for different types of atmospheres. The ion pair production in the atmosphere and the contribution of the different shower components, precisely the electromagnetic, muon and hadronic is estimated according applied atmospheric types. In addition simulations with different hadronic interaction models GHEISHA 2002, FLUKA 2006 and QGSJET II are carried out. The ion pair production in the atmosphere and the contribution of different shower components is estimated according the assumed hadronic interaction models. The yield function Y for total ionization, respectively, for the different components is compared. The observed differences are widely discussed. General conclusion concerning the application of various atmospheric profiles and hadron interaction models is carried out.     

Variations of solar radiation input to the lower atmosphere associated with different helio/geophysical factors

The influence of helio/geophysical factors on the solar energy input to the lower atmosphere has been studied at the network of actinometric stations of Russia in different latitudinal belts. It was found that there are appreciable changes in the half-yearly values of total radiation associated with galactic cosmic ray (GCR) variations in the 11-yr solar cycle, the increase of GCR flux being accompanied by a decrease of the total radiation at higher latitudes and by its increase at lower latitudes. Auroral phenomena and solar flare activity are likely to affect the solar radiation input to the high-latitudinal belt together with GCR variations, the increase of both these factors resulting in the decrease of total radiation. The changes found in the total radiation fluxes in the lower atmosphere seem to be related to the cloud cover variations associated with the solar and geophysical phenomena under study. The variations of the solar radiation input in the 11-yr-cycle amounting to ±4–6% may be an important factor affecting tropospheric dynamics.     

Cosmogenic21Ne and22Ne depth profiles in chondrites

The production rate profiles of²¹Ne and²²Ne as a function of depth in meteoroids due to spallation by solar flare cosmic rays (SCR) and galactic cosmic rays (GCR) are calculated and their dependence on size and composition of meteoroids has been evaluated. The GCR production rate at a given depth increases with size for radii<25cm and then decreases whereas the²²Ne²¹Ne ratio (NeR) generally decreases with size and depth. The calculated GCR production rates and NeR are consistent with the measurements in several Chondrites. A plot of track production rate vs. NeR shows that some chondrites have NeR values smaller than those expected for their sizes. Thes obeervation suggestsat least a two-stage irradiation for such meteorites; the meteoroid exposure as a small body in the interplanetary space must have been preceded by exposure under deep shielding, possibly in its parent body.     

Study of spatial and temporal structure of long-term effects of solar activity and cosmic ray variations on the lower atmosphere circulation

The spatial and temporal structure of the effects of solar activity (SA) and galactic cosmic ray (GCR) flux variations on the lower atmosphere circulation has been studied based on NCEP/NCAR reanalysis archive for 1948–2006 and MSLP (Climatic Research Unit, UK) data for 1873–2000. It has been shown that the GCR effects on pressure variations are characterized by a strong latitudinal and regional dependence, which is determined by specific features of the tropospheric circulation in the studied regions. The distribution of the correlation coefficients for mean yearly values of atmospheric pressure with the GCR flux intensity is closely related to the position of the main climatological fronts. The periodic (∼60 years) changes in the correlation sign of the pressure at high and middle latitudes with Wolf numbers have been revealed. It has been suggested that the changes of the sign of SA/GCR effects on atmospheric pressure are caused by the changes of the macrocirculation epochs, which, in turn, may be related to large-scale processes on the Sun.     

Galactic Cosmic Ray Intensity in the Upcoming Minimum of the Solar Activity Cycle

During the prolonged and deep minimum of solar activity between cycles 23 and 24, an unusual behavior of the heliospheric characteristics and increased intensity of galactic cosmic rays (GCRs) near the Earth’s orbit were observed. The maximum of the current solar cycle 24 is lower than the previous one, and the decline in solar and, therefore, heliospheric activity is expected to continue in the next cycle. In these conditions, it is important for an understanding of the process of GCR modulation in the heliosphere, as well as for applied purposes (evaluation of the radiation safety of planned space flights, etc.), to estimate quantitatively the possible GCR characteristics near the Earth in the upcoming solar minimum (~2019–2020). Our estimation is based on the prediction of the heliospheric characteristics that are important for cosmic ray modulation, as well as on numeric calculations of GCR intensity. Additionally, we consider the distribution of the intensity and other GCR characteristics in the heliosphere and discuss the intercycle variations in the GCR characteristics that are integral for the whole heliosphere (total energy, mean energy, and charge).     

Cosmic ray-induced ionization in the atmosphere: spatial and temporal changes

Detailed calculations of the time-variable spatial distribution of cosmic ray-induced ionization of the lower atmosphere are presented using a physical model. Using the differential energy spectrum of cosmic rays obtained from the worldwide neutron monitor network since 1951 and taking into account also the slow changes in the geomagnetic dipole, we have calculated the corresponding 3D (geographical coordinates and altitude) equilibrium ion concentration in the lower atmosphere as a function of time for the period 1951–2000. A comparison to the results of measurements validates the calculation method, as the calculated cosmic ray-induced ionization reproduces in general the observed altitudinal and latitudinal profiles of the ion concentration. The results of the present work provide a basis for a quantitative study of the solar–terrestrial relationships on long time scales.     

Oscillations of galactic cosmic rays and solar indices before the arrival of relativistic solar protons

Using modern wavelet analysis techniques, we have made an attempt to search for oscillations of intensity of galactic cosmic rays (GCR), sunspot numbers (SS) and magnitudes of coronal index (CI) implying that the time evolution of those oscillations may serve as a precursor of Ground Level Enhancements (GLEs) of solar cosmic rays (SCR). From total number of 70 GLEs registered in 1942–2006, the four large events — 23 February 1956, 14 July 2000, 28 October 2003, and 20 January 2005 — have been chosen for our study. By the results of our analysis, it was shown that a frequency of oscillations of GCR decreases as time approaches to the event day. We have also studied a behaviour of common periodicities of GCR and SCR within the time interval of individual GLE. The oscillations of GLE occurrence rate (OR) at different stages of the solar activity (SA) cycle is of special interest. We have found some common periodicities of SS and CI in the range of short (2.8, 5.2, 27 and 60 days), medium (0.3, 0.5, 0.7, 1.3, 1.8 and 3.2 years) and long (4.6 and 11.0 years) periods. Short and medium periodicities, in general, are rather concentrated around the maxima of solar cycles and display the complex phase relations. When comparing these results with the behaviour of OR oscillations we found that the period of 11 years is dominating (controlling); it is continuous over the entire time interval of 1942–2006, and during all this time it displays high synchronization and clear linear ratios between the phases of oscillations of η, SS and CI. It implies that SCR generation is not isolated stochastic phenomena characteristic exclusively for chromospheric and/or coronal structures. In fact, this process may have global features and involve large regions in the Sun’s atmosphere.     

Influence of energetic particles on atmospheric ionization

Atmospheric ionization rate for the incident electrons and heavy nuclei (protons and α-particles) of various origin is simulated using GEANT4 code. The photon-nuclei cascade generated by a primary particle is included in this code. “Secondary” bremsstrahlung effect from electrons causes peak ionization at low altitudes along with the more powerful peak at high altitudes. The contribution of helium to the total ionization rate, which turns out to be ∼10% of that of protons, is discussed.     

Global magnetic field of the Sun and long-term variations of galactic cosmic rays

The paper deals with the relation of long-term variations of 10 GV galactic cosmic rays (GCR) to the global solar magnetic field and solar wind parameters. This study continues previous works, where the tilt of the heliospheric current sheet (HCS) and other solar-heliospheric parameters are successfully used to describe long-term variations of cosmic rays in the past two solar cycles. The novelty of the present work is the use of the HCS tilt and other parameters reconstructed from Hα observations of filaments for the period when direct global solar magnetic field observations were unavailable. Thus, we could extend the GCR simulation interval back to 1953. The analysis of data for 1953–1999 revealed a good correlation (the correlation coefficient >0.88) between the solar-heliospheric parameters and GCR in different cycles of solar activity. Moreover, the approach applied makes it possible to describe the behavior of cosmic rays in the epochs of solar maxima, which could not be done before. This indicates both the adequacy of the model and the reliability of the reconstructed global solar magnetic field parameters.     

The rigidity spectrum of the harmonics of the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity: 1965–2002

We study temporal changes of the rigidity (R) spectrum of the harmonics of the 27-day variation of the galactic cosmic ray (GCR) intensity using neutron monitors (NM) data for the period 1965–2002. We show that the rigidity spectrum of the third harmonic (9 days) of the 27-day variation of the GCR intensity changes in a similar way as the spectra of the first and second harmonics, being hard in the maximum epochs and soft in the minimum epochs of solar activity. We ascribe this finding to the alternation of the sizes of the modulation regions of the 27-day variation of the GCR intensity in different epochs of solar activity. The average size of the vicinity of the corotating interaction regions, causing the 27-day variation of the GCR intensity, is less in the minimum epochs than in the maximum epochs of solar activity. A vicinity of the corotating interaction regions of larger size involves in modulation higher rigidity particles of GCR than the vicinity of smaller size; thus, this statement can be considered as one of the reasons leading to the hardening of the rigidity spectrum of the harmonics of the 27-day variation of the GCR intensity in maximum epochs compared with minimum epochs of solar activity.We also show that the temporal changes of the power rigidity spectrum of the third harmonic of the 27-day variation of the GCR intensity are negatively correlated with the rigidity spectrum of the 11-year variation of the galactic cosmic ray intensity.We found a recurrence in the temporal changes of the amplitudes of the first harmonic of the 27-day variation of the GCR intensity and in some parameters of solar activity and solar wind.     

The proposed connection between clouds and cosmic rays: cloud behaviour during the past 50–120 years

Several authors have suggested that a link exists between the flux of galactic cosmic rays (GCR) and cloudiness. Here we review the evidence for such a connection from studies of cloud factors using both satellite and ground-based data. In particular, we search for evidence for the low cloud decrease predicted by the rising levels of solar activity and the low cloud-cosmic ray flux correlation indicated by satellite data. Sunshine and synoptic cloud records both indicate that the global total cloud cover has increased during the past century. This increase in total cloud cover argues against a dominating role by solar activity (via GCR) over cloud formation on centennial time scales. Either the predicted low cloud decrease has not occurred or the medium-high level cloud has increased to a greater extent than low cloud has decreased.As there is no accurate long term data available on low cloud behaviour during the last century, we are not able to totally dismiss the link between GCR and cloudiness, but we list a number of arguments for and against the proposed cosmic ray-cloud connection.     

Changes in the chemical composition of the atmosphere in the polar regions of the Earth after solar proton flares (3d modeling)

The paper presents the results of numerical photochemical simulations of the impact of the most powerful solar proton flares during the 23rd solar cycle on the ozonosphere in the polar regions of the Earth. A global 3D photochemical model, CHARM, developed at Central Aerological Observatory (CAO) was used in the simulations. The model introduces an additional source of nitrogen atoms and OH radicals. These components are formed due to the ionization effect of solar protons in the Earth’s atmosphere. The ionization rate was determined from data on proton fluxes measured by GOES satellites. The production rate of additional NO _x and HО _x molecules per ion pair was based on published theoretical studies. It is shown that the most intense flares in the 23rd solar cycle (2000, 2001, and 2003) destroyed ozone in the mesosphere to a great extent (sometimes completely, for example, during the July 14, 2000, event). It is found that the response of ozone to solar proton events follows a seasonal pattern. For the first time, the long-term effect of solar proton events is identified; it is approximately one year.     

Formation of large-scale disturbances in the upper atmosphere caused by acoustic gravity wave sources on the Earth’s surface

The results of a model study of the acoustic gravity wave (AGW) propagation from the Earth’s surface to the upper atmospheric altitudes have been considered. Numerical calculations have been performed using a nonhydrostatic model of the atmosphere, which takes into account nonlinear and dissipative processes originating when waves propagate upward. The model source of atmospheric disturbances has been specified in an area localized on the Earth’s surface. The disturbance source frequency spectrum includes harmonics at frequencies of 0.5ωg-1.5ωg (ωg is the Brunt-Väisälä frequency near the Earth’s surface). The calculations indicated that AGW propagation and dissipation over the source result in the fact that the region of large-scale spatial disturbances of the upper atmosphere mean state is formed at ～200 km altitudes. This region substantially affects AGW propagation and results in waveguide propagation of AGWs with periods shorter than the Väisälä-Brunt period at the altitude of a disturbed atmosphere. The dissipation of AGWs propagating in such a waveguide results in a waveguide horizontal expansion. The extension of the disturbed region of the mean state of the upper atmosphere and, consequently, the waveguide length can reach ～1000 km, if the AGW ground source operates for ～1 h. The physical mechanism by which large-scale disturbances are formed in the upper atmosphere, based on the propagation and dissipation of AGWs with periods shorter than the Väisälä-Brunt period in the upper atmosphere, explains why these disturbances are rapidly generated and localized above AGW sources located on the Earth’s surface or in the lower atmosphere.     

Impact of different energies of precipitating particles on NOx generation in the middle and upper atmosphere during geomagnetic storms

Energetic particle precipitation couples the solar wind to the Earth's atmosphere and indirectly to Earth's climate. Ionisation and dissociation increases, due to particle precipitation, create odd nitrogen (NO_x) and odd hydrogen (HO_X) in the upper atmosphere, which can affect ozone chemistry. The long-lived NO_x can be transported downwards into the stratosphere, particularly during the polar winter. Thus, the impact of NO_x is determined by both the initial ionisation production, which is a function of the particle flux and energy spectrum, as well as transport rates. In this paper, we use the Sodankylä Ion and Neurtal Chemistry (SIC) model to simulate the production of NO_x from examples of the most representative particle flux and energy spectra available today of solar proton events (SPE), auroral energy electrons, and relativistic electron precipitation (REP). Large SPEs are found to produce higher initial NO_x concentrations than long-lived REP events, which themselves produce higher initial NO_x levels than auroral electron precipitation. Only REP microburst events were found to be insignificant in terms of generating NO_x. We show that the Global Ozone Monitoring by Occultation of Stars (GOMOS) observations from the Arctic winter 2003–2004 are consistent with NO_x generation by a combination of SPE, auroral altitude precipitation, and long-lived REP events.     

Solar proton activity during cycle 23 and changes in the ozonosphere: Numerical simulation and analysis of observational data

Using the data on solar proton fluxes measured on board the GOES satellites, the most powerful solar proton events (SPEs) of solar cycle 23 are selected, and ionization rates in the atmosphere in these periods at high latitudes of the Northern Hemisphere are calculated. Assuming that each ion pair formed at the retardation of solar protons in the atmosphere leads to the formation of 1.25 molecules of nitric oxide, 2.0 molecules of the OH radical, and one oxygen atom, changes in the content of ozone, nitrogen and other compounds were calculated using a photochemical model. The calculations showed that the strongest ionization and destruction of ozone was caused by SPEs that occurred on July 14, 2000; November 8, 2000; November 4, 2001; and October 28, 2003. The results can form the basis for compiling the catalog of changes in ionization and ozone in the atmosphere caused by solar proton activity.     

Equatorial temperature and wind anomaly (ETWA)—a review

We review the observational evidence for describing the characteristics of the equatorial temperature and zonal wind anomaly (ETWA) in the low-latitude thermosphere in solar maximum and minimum periods. In spite of some new results on ETWA in the last decade, including its discovery, there is no satisfactory explanation in our understanding of the phenomenon. The two suggested mechanisms for heating at the crests of the equatorial ionization anomaly (EIA) to form the equatorial temperature anomaly (ETA) are due to (1) the ion-drag on the zonal winds resulting in the transfer of kinetic energy into heat energy and (2) the exothermic chemical reactions involved in the dominant O⁺ion re-combinations. To verify which of the two suggested mechanisms is the most effective in causing ETA, it is necessary to measure simultaneously a few parameters in situ by the satellite-borne instruments. They are (1) the electron density (N_e) and temperature (T_e), (2) the molecular and atomic ion densities (N_ij) and ion temperatures (T_i), (3) the gas temperatures (T) and densities of the gas constituents, (4) the vector winds or at least the zonal (Z) and vertical (V) wind components and (5) the drift velocities of the ionization. These together with the simultaneous ground-based measurements, will resolve identifying not only the dominant mechanism(s) for ETWA, but also the processes responsible for the enigmatic phenomena, such as the equatorial spread-F (ESF), the midnight temperature maximum (MTM) and the possible role of the EIA in their occurrences.     

Helium isotopes in the earth's interior and in the atmosphere: A degassing model of the earth

A first-order degassing model was applied to describe the evolution of helium content and isotope composition in the earth and in the atmosphere. The main events described by the model are: (1) the earth-trapped primordial rare gases at the moment of its accretion; (2) later, the solid earth lost primordial and radiogenic rare gases, and (3) they were accumulated in the atmosphere; (4) in addition,³He was formed in the atmosphere due to cosmic irradiation, accretion from solar wind, etc.; (5)³He and⁴He dissipated into space at different loss rates.Study of this model confirms the concept that some of primordial helium is retained in the interior of the earth; terrestrial helium (³He/⁴He～ 2 × 10^?5) was most probably formed as a mixture of primordial (³He/⁴He= 3 × 10^?4) and radiogenic (³He/⁴He～ 3 × 10^?8) helium. For achondritic concentrations of heavy radioactive elements (U= 2.25 × 10^?8g/g) the calculated⁴He flux from the earth is equal to 5.7 × 10⁶ at cm^?2 sec^?1. The corresponding³He flux is about 114 at cm^?2 sec^?1. In discussing the aeronomic problem of helium it is necessary to take into account that the earth is the main source of the light helium isotope.     

Acoustic energy transfer to the upper atmosphere from sinusoidal sources and a role of nonlinear processes

In the framework of solving the problem of acoustic energy transfer from near-surface sources through the upper atmosphere, the propagation of sinusoidal signals of different origin is studied. All calculations are made by means of a model that takes into account the inhomogeneity of the atmosphere, nonlinear effects, absorption, divergence of wave front due to long-range acoustic wave propagation, etc., but does not include the effect of gravity. Infrasonic waves of various periods and their absorption at various heights of the atmosphere are investigated. The calculations show that a sinusoidal signal is destroyed by nonlinear processes during its upward propagation; it transforms into two, initial and final, impulses. The location of the “transformation zone” depends on frequency; its height increases with decreasing frequency. The acoustic waves can heat the upper atmosphere, for example, waves with a period of 3 min generated by thunderstorms can heat the atmosphere by up to ΔT_a=13.08 K/day in the region of 323–431 km. The efficiency of a point artificial emitter is too weak to heat the atmosphere significantly.     

Ball lightning as a radiation bubble

Summary A model for ball lightning is presented and its properties discussed. The model is that of a microwave radiation field contained within a plasma dielectric sphere, resonant at a frequency much greater than the electron-molecule collision frequency. Calculations are made of the energy stored in the microwave field, the electron temperature, the rate of energy loss due to ionization, and the effects of recombination. It is concluded (a) that a self-consistent set of conditions can be obtained only when the neutral density within the sphere is much lower than atmospheric, and (b) that the microwave field strength required is of the order of 10⁹V/m. Under these conditions the radiation pressure becomes comparable to that of the atmosphere; the second case investigated is when these two become equal. It is then found that this radiation bubble appears to satisfy the requirements imposed on energy storage and dissipation, recombination and resonant frequency. The electron density in such a bubble is found to be typically 10¹¹ cm^–3 and the stored energy typically 10³ Joules. A discussion is given of some of the problems associated with the model-formation, stability, neutral number density, and the problem of hydrostatic equilibrium.