Creation of model of quasi-trapped proton fluxes below Earth’s radiation belt

The object of investigation is the phenomenon of proton (from tens keV to several MeV) flux enhancement in near-equatorial region (L < 1.15) at altitude up to ∼1300 km (the storm-time equatorial belt). These fluxes are quite small but the problem of their origin is more interesting than the possible damage they can produce. The well known sources of these protons are radiation belt and ring current. The mechanism of transport is the charge-exchange on neutral hydrogen of exosphere and the charge-exchange on oxygen of upper atmosphere. Therefore this belt is something like the ring current projection to low altitudes. Using the large set of satellites data we obtain the average energy spectrum, the approximation of spectrum using kappa-function, the flux dependence on L, B geomagnetic parameters. On the basis of more than 30 years of experimental observations we made the empiric model that extends model of proton fluxes below 100 keV in the region of small L-values (L < 1.15). The model was realized as the package of programs integrated into COSRAD system available via Internet. The model can be used for revision of estimation of dose that low-orbital space devices obtain.     

FEM code simulation of the magnetospheric proton fluxes

We report a study of the numeric solution to the diffusive transport equation for energetic protons magnetically trapped in the Earth's equatorial magnetosphere. The analysis takes into account the pertinent physical processes in this region, including deceleration of protons by Coulomb collisional interactions with free and bound electrons, the charge exchange process, cosmic ray albedo neutron decay source, and electric and magnetic radial diffusion. These results were obtained using the Finite Element Method with magnetic moment and geomagnetic L-shell as free variables. Steady state boundary conditions were imposed at L=1 as zero distribution function and at L=7 with proton distribution function extracted from ATS 6 satellite observations. The FEM-code yields unidirectional proton flux in the energy range of 0.1–1000 MeV at the equatorial top of the geomagnetic lines, and the results are found to be in satisfactorily agreement with the empirical NASA AP-8 model proton flux within the energy range of 0.5–100 MeV. Below 500 keV, the empirical AP-8 model proton fluxes are several orders of magnitude greater than those computed with the FEM-code at L<3. This discrepancy is difficult to explain by uncertainties of boundary spectrum parameters or transport coefficients.     

A model for stochastic acceleration of electrons during geomagnetic storms

The theory of resonant diffusion is extended to fully relativistic plasmas, and we examine resonant interactions between electrons and electromagnetic R mode (whistler) and L-mode (EMIC) waves. Resonant diffusion curves are constructed for plasma parameters representative of the Earth's storm time magnetosphere, both inside and outside the plasmapause. EMIC waves can resonate with electrons > 1 MeV, but the energies remain nearly constant along the diffusion curves. Storm-time EMIC waves can induce rapid pitch—angle scattering, but the waves are ineffective for stochastic acceleration of elections. Substantial energy change can occur along the diffusion curves for interactions between resonant electrons and whistler—mode waves, especially in regions of low plasma density. Specifically, whistlers can accelerate electrons from energies near 100 keV to above 1 MeV outside the plasmapause. A model is proposed comprising energy diffusion by whistler-mode chorus and pitch-angle scattering by EMIC waves to account for the gradual acceleration of electrons over the region 4 ≤ L ≤ 6 during the recovery phase of a geomagnetic storm.     

A simple model for the interhemispheric coupling of the middle atmosphere circulation

The interhemispheric coupling of the middle atmosphere general circulation is characterized by a global anomaly pattern of the zonal-mean temperature. This pattern reflects an anomalous stratospheric and mesospheric residual circulation, in which a weaker (stronger) stratospheric winter circulation is linked to an upward (downward) shift of its upper mesospheric branch reaching from the summer to the winter pole. This phenomenon is robust in observational data and several middle atmosphere general circulation models. In the present study, the recently proposed mechanism of the interhemispheric coupling is unequivocally proven within the framework of a zonally symmetric model that excludes any additional effects due to resolved waves and non-zonally propagating gravity waves. Two simulations are conducted that differ in the strength of the polar vortex. A weaker polar vortex results in a downward shift of the winter mesospheric gravity wave drag. This leads to changes also in the summer upper mesosphere via a feedback solely between gravity wave breaking and the zonal-mean state. The accompanying temperature anomaly reproduces the pattern of the interhemispheric coupling.     

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

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

Modeling the (upper) solar atmosphere including the magnetic field

The atmosphere of the Sun is highly structured and dynamic in nature. From the photosphere and chromosphere into the transition region and the corona plasma-β changes from above to below one, i.e., while in the lower atmosphere the energy density of the plasma dominates, in the upper atmosphere the magnetic field plays the governing role – one might speak of a “magnetic transition”. Therefore the dynamics of the overshooting convection in the photosphere, the granulation, is shuffling the magnetic field around in the photosphere. This leads not only to a (re-)structuring of the magnetic field in the upper atmosphere, but induces also the dynamic reaction of the coronal plasma, e.g., due to reconnection events. Therefore the (complex) structure and the interaction of various magnetic patches is crucial to understand the structure, dynamics and heating of coronal plasma as well as its acceleration into the solar wind.

The present article will emphasize the need for three-dimensional modeling accounting for the complexity of the solar atmosphere to understand these processes. Some advances on 3D modeling of the upper solar atmosphere in magnetically closed as well as open regions will be presented together with diagnostic tools to compare these models to observations. This highlights the recent success of these models which in many respects closely match the observations.     

The hydroxyl radical as an indicator of SEP fluxes in the high-latitude terrestrial atmosphere

The low background values at nighttime of the mesospheric hydroxyl (OH) radical make it easier to single out the atmospheric response to the external solar forcing in Polar Regions. Because of the short lifetime of HO_x, it is possible to follow the trails of Solar Energetic Particle (SEP) events in the terrestrial atmosphere, as shown by Storini and Damiani (2008). The sensitivity of this indicator makes discernible not only extreme particle events with a flux peak of several thousand pfu [1 pfu = 1 particle/(cm² s sr)] at energies >10 MeV, but also those with lower flux up to about 300 pfu. Using data from the Microwave Limb Sounder (MLS) on board the EOS AURA satellite, we examined the correlation of OH abundance vs. solar proton flux for almost all the identified SEP events spanning from November 2004 to December 2006 (later on no more SEP events occurred during Solar Cycle no. 23). The channels at energies greater than 5 MeV and 10 MeV showed the best correlation values (r ∼ 0.90–0.95) at altitudes around 65–75 km whereas, as expected, the most energetic channels were most highly correlated at lower altitudes. Therefore, it is reasonably possible to estimate the solar proton flux from values of mesospheric OH (and viceversa) and it could be useful in studying periods with gaps in the records of solar particles.     

Modeling the time-intensity profile of solar flare generated particle fluxes in the inner heliosphere.

It is possible to model the time-intensity profile of solar particles expected in space after the occurrence of a significant solar flare on the sun. After the particles are accelerated in the flare process, if conditions are favorable, they may be released into the solar corona and then into space. The heliolongitudinal gradients observed in the inner heliosphere are extremely variable, reflecting the major magnetic structures in the solar corona which extend into space. These magnetic structures control the particle gradients in the inner heliosphere. The most extensive solar particle measurements are those observed by earth-orbiting spacecraft, and forecast and prediction procedures are best for the position of the earth. There is no consensus of how to extend the earth-based models to other locations in space. Local interplanetary conditions and structures exert considerable influence on the time-intensity profiles observed. The interplanetary shock may either reduce or enhance the particle intensity observed at a specific point in space and the observed effects are very dependent on energy.     

Monte Carlo study of secondary electrons and X-rays produced by different angular distributions of primary precipitating electrons interacting with the atmosphere

We present Monte Carlo simulations of the energy deposited by energetic electrons from the outer Van Allen belt precipitating into the atmosphere below 100 km. Consideration of the possible mechanism for the precipitation of outer belt electrons suggests that the electrons should be moving nearly horizontally when the enter the atmosphere; therefore we simulate and compare the precipitation of close to horizontally moving, vertically moving, and isotropically distributed electrons. We find that there should be a number of observable differences between the energy deposition for these different angular distributions of electrons.     

Some peculiarities of longitudinal distribution of proton fluxes at high latitudes

Dynamical features of proton fluxes at high and middle latitudes were studied based on data measured by Sun-synchronous low-altitude (800 km height) polar-orbiting three NOAA series satellites: POES 15, 16, and 17 during the geomagnetic storm on January, 21–22, 2005. Data from three satellites that passed the Northern hemisphere along different MLTs allow reconstructing the longitudinal distribution of the proton fluxes. Measurements of protons with energies of 30–80 keV and 80–240 keV (the ring current energy range) by 0- and 90-detectors were used to evaluate and compare the longitudinal asymmetry of proton flux distribution measured in the regions equatorward and poleward of the isotropic boundary. It was found that during all the phases of the geomagnetic storm distribution of the maximum flux of precipitating protons (0-detector data) is sufficiently asymmetric. The maximal flux position along MLT is moving from pre-midnight sector in quiet time to post-midnight one before and during SSC and moving back during recovery phase. The longitudinal distribution of precipitation maxima demonstrates the local increase in afternoon sector (approximately at 13:30 MLT) and decrease in the dusk one during SSC. These features are evident consequence of the magnetosphere compression. To identify the origin of the particles, the locations of maximum fluxes have been projected to the magnetosphere. It was determined that during geomagnetic storm main and recovery phases maximum fluxes were measured at latitudes poleward of the isotropic boundary. To evaluate the trapped particle flux asymmetry, the particle fluences (90-detector data) were calculated along the satellite orbit from L = 2 to the isotropic boundary. The total fluences of trapped particles calculated along the satellite orbit show regular asymmetry between dusk and dawn during main and recovery phases. The maximal intensity of proton fluxes of both investigated populations located poleward and equatorward of the isotropic boundary is achieved during SSC. The total flux measured during crossing the anisotropic region can be considered as a proxy for ring current injection rate.     

On the adjust of minor constituents density in middle atmosphere during solar proton events

Profiles of O₃ partial pressure and of other minor atmospheric constituents (NO, NO₂, HCL, HF and H₂O), observed in the middle atmosphere during Solar proton events (20.04.1998; 05.04.2000), were analysed. Conclusions were drawn that under SCR impact a short-term O₃ partial pressure increase and destruction of some freon constituents took place.     

Trapped proton fluxes measured on board LEO satellites in comparison with models

The proton fluxes from the low-Earth orbital satellites databases (NPOES-17 and CORONAS-F) were analyzed for the quiet geomagnetic period in April 2005. The satisfactory consent was found between the experimental and the AP8 model fluxes of the trapped protons with energy more than ∼10 MeV. At the same time, trapped proton fluxes with energy less than 10 MeV measured by LEO satellites were higher than the ones predicted by the AP8 model in the region of the SAA (drift shell L < 1.5).     

A middle atmosphere temperature reference model from satellite measurements

Temperature fields in the stratosphere and mesosphere have been derived from radiance measurements made by the Nimbus 5 SCR, the Nimbus 6 PMR, and the Nimbus 7 SAMS and LIMS radiometers. These instruments cover different latitude and height ranges and different times during the 1973–1983 period. The problems of combining different data sets are discussed, and examples from a proposed model atmosphere for the stratosphere and mesosphere are presented. The model is given in terms of zonal means and amplitude and phase of zonal waves 1 and 2 for temperature and geopotential height, as functions of latitude and pressure for each calendar month. Comparisons are made with the CIRA 1972 and the Koshelkov Southern Hemisphere models and with the SAMS results and in-situ rocket/radio sondes.     

A model for calculating the vertical distribution of the atmospheric electric potential in the exchange layer in a maritime clean atmosphere

A theoretical model is developed for calculating the vertical distribution of atmospheric electric potential in exchange layer of maritime clean atmosphere. The transport of space charge in electrode layer acts as a convective generator in this model and plays a major role in determining potential distribution in vertical. Eddy diffusion is the main mechanism responsible for the distribution of space charge in vertical. Our results show that potential at a particular level increases with increase in the strength of eddy diffusion under similar conditions. A method is suggested to estimate columnar resistance, the ionospheric potential and the vertical atmospheric electric potential distribution in exchange layer from measurements of total air-earth current density and surface electric field made over oceans. The results are validated and found to be in very good agreement with the previous aircraft measurements. Different parameters involved in the proposed methodology can be determined either theoretically, as in the present work, or experimentally using the near surface atmospheric electrical measurements or using some other surface-based measurement technique such as LIDAR. A graphical relationship between the atmospheric eddy diffusion coefficient and height of exchange layer obtained from atmospheric electrical approach, is reported.     

A numerical model characterizing internal gravity wave propagation into the upper atmosphere

A new two-dimensional, time-dependent and fully nonlinear model is developed to numerically simulate plane wave motions for internal gravity waves in a non-isothermal and windy atmosphere that accounts for the dissipation due to eddy and molecular processes. The atmosphere has been treated as a well mixed total gas with a constant mean molecular weight. The effects of Rayleigh friction and Newtonian cooling are applied near the upper boundary of the model to simulate the radiation conditions as well as to act as a sponge layer; lateral boundaries are periodic over a horizontal wavelength to simulate a horizontally infinite domain. The thermal excitation to initiate upward propagating waves is spatially localized in the troposphere and is a Gaussian function of time. A time-splitting technique is applied to the finite difference equations that are derived from the Navier–Stokes equations. The time integration for these highly coupled equations is performed using an explicit second order Lax–Wendroff scheme and an implicit Newton–Raphson scheme. The wave solutions derived from the model are found to be broadly agreeable with those derived from a Wentzel–Kramers–Brillouin theory.     

Proposal for a reference model of the middle atmosphere of the southern hemisphere

Available rocketsonde information has been used to compile tables of monthly mean temperature, pressure, density and zonal wind for the middle atmosphere of the southern hemisphere with the purpose of revising similar tables presented to COSPAR earlier. The altitude range is 25 to 80 km in steps of 5 km. The latitude range is 0° to 70°S with a 10° step. The compatability of different sets of temperature measurements is discussed. Mean values of temperature, pressure and zonal wind obtained for the southern hemisphere are compared with northern hemisphere model values. Large differences between the hemispheres (up to 20°C in temperature, 20–30% in pressure, 30–50 m/s in wind) imply that reference atmospheres such as CIRA should be complemented by southern hemisphere climatology.     

Effects of dust polarity and nonextensive electrons on the dust-ion acoustic solitons and double layers in earth atmosphere

Propagation of dustion acoustic solitary waves (DIASWs) and double layers is discussed in earth atmosphere, using the Sagdeev potential method. The best model for distribution function of electrons in earth atmosphere is found by fitting available data on different distribution functions. The nonextensive function with parameter $q=0.58$ provides the best fit on observations. Thus we analyze the propagation of localized waves in an unmagnetized plasma containing nonextensive electrons, inertial ions, and negatively/positively charged stationary dust. It is found that both compressive and rarefactive solitons as well as double layers exist depending on the sign (and the value) of dust polarity. Characters of propagated waves are described using the presented model.     

Comparing proton fluxes of central meridian SEP events with those predicted by SOLPENCO

We have developed an operational code, SOLPENCO, that can be used for space weather prediction schemes of solar energetic particle (SEP) events. SOLPENCO provides proton differential flux and cumulated fluence profiles from the onset of the event up to the arrival of the associated traveling interplanetary shock at the observer’s position (either 1.0 or 0.4 AU). SOLPENCO considers a variety of interplanetary scenarios where the SEP events develop. These scenarios include solar longitudes of the parent solar event ranging from E75 to W90, transit speeds of the associated shock ranging from 400 to 1700 km s⁻¹, proton energies ranging from 0.125 to 64 MeV, and interplanetary conditions for the energetic particle transport characterized by specific mean free paths. We compare the results of SOLPENCO with flux measurements of a set of SEP events observed at 1 AU that fulfill the following four conditions: (1) the association between the interplanetary shock observed at 1 AU and the parent solar event is well established; (2) the heliolongitude of the active region site is within 30° of the Sun–Earth line; (3) the event shows a significant proton flux increase at energies below 96 MeV; (4) the pre-event intensity background is low. The results are discussed in terms of the transit velocity of the shock and the proton energy. We draw conclusions about both the use of SOLPENCO as a prediction tool and the required improvements to make it useful for space weather purposes.     

Measurements of high-energy neutron and proton fluxes on-board “Mir-Spectr” orbital complex

The measurements of high-energy neutron (with energies 30–300 MeV) and proton (with energies 1–200 MeV) fluxes are being conducted on-board “Mir-Spectr” orbital complex. Neutrons are detected by the undirected (FOV 4π sr) scintillator spectrometer, consisting of 4 identical CsI(Tl) detector units (the effective area for neutrons 30 cm²). The gamma-quanta, which can be also detected by this instrument, are separated from neutrons by the analysis of the scintillator output pulse shape. To exclude registration of charged particles an anticoincidence plastic scintillator shield is realized in each detector unit. The proton fluxes are measured by the telescope based on 3 semiconductor detectors with small geometry factor (1 cm²×sr). As the first result of the experiment the upper limit of the integral flux of local and albedo neutrons in the equatorial region (L<1.1) was estimated. The results of this measurements can be useful for the radiation security. Also, the neutrons of solar flares can be detected in this experiment.     

A numerical model approximating extreme energetic electron events involved in the physical and chemical processes of the middle atmosphere

Relativistic electrons (with energies >150 keV) which originate in the outer radiation belt and detected by the Russian ‘Meteor’ series of satellites have been correlated with the atmospheric total ozone data compiled by almost 90 stations located around the world within the latitude zone 40°–70°N. In more than 60% of the stations examined we have detected a clear decrease of the ozone 3–5 days after the electron flux excess. A numerical model has been applied to approximate this effect based on relativistic electron initiated nitric oxides creation in the upper mesosphere with subsequent atmospheric transport (both vertical and horizontal) towards the upper stratosphere. A first attempt of local and temporal prediction of ozone depletion because of energetic electrons impact in the middle atmosphere has been illustrated.