太阳高能粒子事件强度与相关双CME事件关系的两个事例

太阳高能粒子事件常伴随太阳耀斑和日冕物质抛射事件(Coronal Mass Ejections,CME)出现,由于太阳高能粒子事件的关键因素是双CME的相互作用,利用SOHO卫星观测的高能粒子强度、耀斑强度以及CME的相对高度与时间,通过高度与时间拟合得到的速度,分析了2001年4月15日和2005年1月20日的太阳高能粒子事件强度与相关双CME事件的关系,发现这两个太阳高能粒子事件中E ≥ 10MeV质子的强度与双CME事件无关.因此在这两次太阳高能粒子事件早期,E ≥ 10MeV质子的强度只与相关太阳耀斑和CME有关.      

Solar proton events recorded in the stratosphere during cosmic ray balloon observations in 1957–2008

Long-term balloon observations have been performed by the Lebedev Physical Institute since 1957 up to the present time. The observations are taken several times a week at the polar and mid latitudes and allow us to study dynamics of galactic and solar cosmic ray as well as secondary particle fluxes in the atmosphere and in the near-Earth space. Solar energetic particles (120) – mostly protons – (SEP) events with >100 MeV proton intensity above 1 cm⁻² s⁻¹ s⁻¹ were recorded during 1958–2006. Before the advent of the SEP monitoring on spacecraft these results constituted the only homogeneous series of >100 MeV SEP events. The SEP intensities and energy spectra inferred from the Lebedev Physical Institute observations are consistent with the results taken in the adjacent energy intervals by the spacecraft and neutron monitors. Joint consideration of the SEP events series recorded by balloons and by neutron monitors during solar cycles 20–23 makes it possible to restore the probable number of events in solar cycle 19, which was not properly covered by observations. Some correlation was found between duration of SEP event production in a solar cycle and sunspot cycle characteristics.     

1996-2002年太阳耀斑的统计分析

分析了1996-2002年南北半球的太阳黑子相对数和南北半球太阳X射线耀斑级别(简称Imp)≥M1．0的太阳X射线耀斑的特征和不对称性．分析结果表明，南北半球的太阳耀斑活动的程度交替上升，在2001年7月以前北半球的太阳耀斑活动强于南半球，2001年7月开始耀斑活动逐渐以南半球为主．本文还逐月分析了1996—2001年南北半球的耀斑指数．2000年7月为第23周太阳指数最大的一个月，与第23周太阳黑子相对数最大月均值吻合．     

太阳黑子数及Ap指数周期变化特征的小波分析

应用Morlet小波变换方法从多个变化尺度上对1932—2000年的太阳黑子数及Ap指数的变化特征进行分析．(1)太阳黑子数存在准11年、准32年的周期变化特征及Ap指数存在准32年、准11年、准6个月、准27天和准13．9天的周期变化特征；(2)太阳黑子数及Ap指数有着相似的准11年周期变化，但Ap指数极值的出现要比太阳黑子数极值出现滞后1—2年；(3)Ap指数准27天的周期变化在太阳黑子活动高、低年不同，在太阳活动低年，Ap指数有着较稳定的准27天周期变化，但在太阳活动高年，27天的周期变化几乎消失，这种周期变化的消失和出现时间可在Morlet小波变换图中体现出来。     

Relative ion Fe,C and O abundances in quiet time particle fluxes in the 23 SC

Using ACE and SOHO data the origin of quiet-time low-energy particle fluxes at 1 AU is studied in the 23rd solar cycle. One of the selection criteria of quiet-time periods is to demand that H/He < 10 provided that periods with noticeable contribution of remnants of gradual events have been excluded from consideration. Our results suggest different origin of 0.03–3 MeV/nucleon particles – different seed populations accelerated and different acceleration processes. During the ascending, maximum and descending phases of solar activity quiet-time ions consist of coronal particles accelerated to suprathermal energies in about a half of the quiet periods, the rest of quiet-time fluxes originates from particle acceleration in processes similar to those in small impulsive solar flares rich in Fe. At solar minimum the bulk solar wind particles serve as seed population.     

Solar particle precipitation into the polar atmosphere and their dependence on hemisphere and local time

The precipitation of solar energetic particles, protons as well as electrons, at high latitudes is commonly assumed to be homogeneous across both polar caps. Using Low-Earth Orbit POES (Polar Orbiting Environmental Satellites) we determine particle penetration ratios into the polar atmosphere for protons ranging from about 0.1 MeV to 500 MeV and for electrons spanning about one order of magnitude in energy with a maximum of 0.3 MeV. Based on power law fits for the POES spectrum we show, that for energies interesting for middle and lower atmosphere chemistry, particle flux over the poles is comparable in magnitude to flux at the geostationary orbit or at L1 in interplanetary space. The time period under study are the solar energetic particle (SEP) event series of October/November 2003 and January 2005.     

Sunspot activity and cosmic ray modulation at 1 a.u. for 1900–2013

The descent of sunspot cycle 23 to an unprecedented minimum of long duration in 2006–2009 led to a prolonged galactic cosmic ray (GCR) recovery to the highest level observed in the instrumental era for a variety of energetic charged particle species on Earth, over a wide range of rigidities. The remarkable GCR increase measured by several ground-based, balloon-borne, and detectors on a satellite is described and discussed. It is accompanied by a decrease in solar wind velocity and interplanetary magnetic field at 1 a.u., reaching the lowest values since measurements of the solar wind began in October 1963; the solar polar field strength (μT) measured at the Wilcox Solar Observatory (WSO) is also significantly reduced compared to prior cycles since the start of the program in 1976, the polar field in the northern hemisphere reversed in June 2012 and again in February 2014, that in the southern hemisphere reversed in July 2013. If updates of WSO data confirm the second reversal in northern solar hemisphere, it would pose a serious challenge to the Dynamo Theory. The long-term change in solar behavior may have begun in 1992, perhaps earlier. The physical underpinnings of these solar changes need to be understood and their effect on GCR modulation processes clarified. The study discusses the recent phenomena in the context of GCR modulation since 1900. These happenings affected our empirical predictions for the key parameters for the next two sunspot cycles (they may be progressively less active than sunspot cycle 24) but it enhanced support for our prediction that solar activity is descending into a Dalton-like grand minimum in the middle of the twentyfirst century, reducing the frequency of the coronal mass ejections; they determine the space weather affecting the quality of life on Earth, radiation dose for hardware and human activities in space as well as the frequency of large Forbush decreases at 1 a.u.     

Galactic cosmic ray flux simulation and prediction.

A dynamic galactic cosmic ray model is proposed to quantitatively describe the z=1-28 ions and electrons of E=10-10(5) MeV/nucleon and their particle flux variations around the Earth's orbit and beyond the Earth's magnetosphere due to diverse large-scale variations of solar activity factors. The variations of large-scale heliospheric magnetic fields and the galactic cosmic ray flux variation time delays relative to solar activity variations are simulated. The lag characteristics and sunspot number predictions having been determined in detail, the model can be used to predict galactic cosmic ray flux levels.     

Morphological studies on ionospheric VHF scintillations over an Indian low latitude station during a solar cycle period (2001–2010)

The amplitude scintillations data recorded at 244 MHz from the geostationary satellite, FLEETSAT (73°E) at a low latitude station, Waltair (17.7°N, 83.3°E) during the ten year period of high to low solar activity from 2001 to 2010 is considered to study the occurrence characteristics of the VHF scintillations. A close association between the intense scintillations on VHF signals during pre-midnight hours, associated with range type of spread-F on ionograms and a relatively weak and slow fading scintillations during post-midnight hours associated with frequency type of spread-F is observed during the relatively high sunspot years from 2001 to 2004, whereas during the low sunspot years from 2005 to 2010 the scintillation activity as well as spread-F activity are found to be minimum. During both the high and low sunspot years, it is observed that the maximum scintillation activity occurs during equinoctial months followed by winter with the minimum occurrence during summer months. The annual mean percentage occurrence of scintillations is found to be clearly associated with the variations in the annual mean sunspot number. The nocturnal variations in the occurrence of scintillations show the onset of scintillation activity starts from 19:00 h LT with maximum of occurrence around 21:00 h LT. A clear semiannual variation in the occurrence of scintillations is observed during pre-midnight hours with two peaks in equinoctial months of March/April and September/October. The number of scintillation patches observed is found to be more during pre-midnight hours compared to those during post-midnight hours. The most probable scintillation patch duration lies around 30 min. Further, it is also found that the number of scintillation patches with durations of 60 min and more decreases with the increase in the patch duration. It is also observed in general that the scintillation activity is inhibited during geomagnetic disturbed days.     

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.     

Synoptic magnetic field in cycle 23 and in the beginning of the cycle 24

The SOHO/MDI data provide the uniform time series of the synoptic magnetic maps which cover the period of the cycle 23 and the beginning of the cycle 24. It is very interesting period because of the long and deep solar minimum between the cycles 23 and 24. Synoptic structure of the solar magnetic field shows variability during solar cycles. It is known that the magnetic activity contributes to the solar irradiance. The axisymmetrical distribution of the magnetic flux (Fig. 3c) is closely associated with the ‘butterfly’ diagram in the EUV emission (Benevolenskaya et al., 2001). And, also, the magnetic field (B_∥) shows the non-uniform distributions of the solar activity with longitude, so-called ‘active zones’, and ‘coronal holes’ in the mid-latitude. Polar coronal holes are forming after the solar maxima and they persist during the solar minima. SOHO/EIT data in the emission of Fe XII (195 Å) could be a proxy for the coronal holes tracking. The active longitudinal zones or active longitude exist due to the reappearance of the activity and it is clearly seen in the synoptic structure of the solar cycle. On the descending branch of the solar cycle 23 active zones are less pronounced comparing with previous cycles 20, 21 and 22. Moreover, the weak polar magnetic field precedes the long and deep solar minimum. In this paper we have discussed the development of solar cycles 23 and 24 in details.     

中等磁暴太阳周分布的统计研究

依据实际观测的中等磁暴数据,统计分析了中等磁暴的太阳周分布.分析结果表明,在一个太阳活动周内,每年中等磁暴随时间的变化出现多个峰值,其中,最大峰值均出现在太阳活动周的下降段,即中等磁暴的峰值比太阳黑子数平滑年均值的峰值要滞后,滞后的时间为2～3年.超过70%的中等磁暴出现在太阳活动周的下降段,这表明绝大多数中等磁暴出现在太阳活动周的下降段.通过对中等磁暴平滑月均值与太阳黑子数平滑月均值相位差的计算分析发现,中等磁暴峰值出现的时间比太阳黑子数峰值出现的时间要滞后,不同太阳活动周中等磁暴峰值出现的时间与太阳黑子数峰值时间滞后的程度不同.      

Validation of IRI-2012 TEC model over Ethiopia during solar minimum (2009) and solar maximum (2013) phases

This paper investigates the capacity of the latest version of the International Reference Ionosphere (IRI-2012) model in predicting the vertical Total Electron Content (vTEC) over Ethiopian regions during solar minimum (2009) and solar maximum (2013) phases. This has been carried out by comparing the IRI-2012 modeled and experimental vTEC inferred from eight ground based dual frequency GPS (Global Positioning System) receivers installed recently at different regions of the country. In this work, the diurnal, monthly and seasonal variation in the measured vTEC have been analyzed and compared with the IRI-2012 modeled vTEC. During the solar minimum phase, the lowest and highest diurnal peak of the experimental vTEC are observed in July and October, respectively. In general, the diurnal variability of vTEC has shown minimum values around 0300 UT (0600 LT) and maximum values between around 1000 and 1300 UT (1300 and 1600 LT) during both solar activity phases. Moreover, the maximum and minimum monthly and seasonal mean hourly vTEC values are observed in October and July and in the March equinox and June solstice, respectively. It is also shown that the IRI-2012-model better predicts the diurnal vTEC in the time interval of about 0000–0300 UT (0300–0600 LT) during the solar minimum phase. However, the model generally overestimates the diurnal vTEC except in the time interval of about 0900–1500 UT (1200–1800 LT) during the solar maximum phase. The overall result of this work shows that the diurnal vTEC prediction performance of the model is generally better during the solar minimum phase than during solar maximum phase. Regarding the monthly and seasonal prediction capacity of the model, there is a good agreement between the modeled and measured monthly and seasonal mean hourly vTEC values in January and December solstice, respectively. Another result of the work depicts that unlike the GPS–TEC the IRI-2012 TEC does not respond to the effect resulted from geomagnetic storms.     

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.     

Solar proton enhancements in different energy channels and coronal mass ejections during the last solar cycle

The main properties of 11622 coronal mass ejections (CMEs) observed by the Solar and Heliospheric Observatory (SOHO) mission’s Large Angle and Spectrometric Coronagraph (LASCO-C2) from January 1996 through December 2006 are considered. Moreover, the extended database of solar proton enhancements (SPEs) with proton flux >0.1 pfu at energy >10 MeV measured at the Earth’s orbit is also studied. A comparison of these databases gives new results concerning the sources and acceleration mechanisms of solar energetic particles. Specifically, coronal mass ejections with width >180° (wide) and linear speed >800 km/s (fast) seem they have the best correlation with solar proton enhancements. The study of some specific solar parameters, such as soft X-ray flares, sunspot numbers, solar flare index etc. has showed that the soft X-ray flares with importance >M5 may provide a reasonable proxy index for the SPE production rate. From this work, it is outlined that the good relation of the fast and wide coronal mass ejections to proton enhancements seems to lead to a similar conclusion. In spite of the fact that in the case of CMEs the statistics cover only the last solar cycle, while the measurements of SXR flares are extended over three solar cycles, it is obvious for the studied period that the coronal mass ejections can also provide a good index for the solar proton production.     

Experimental data and analysis of the October 2003 Forbush decrease

On October 28, 2003 an Earthward-directed coronal mass ejection (CME) was observed from SOHO/LASCO imagery in conjunction with an X17 solar flare. The CME, traveling at nearly 2000 km/s, impacted the Earth on October 29, 2003 causing ground-based particle detectors to register a counting rate drop known as a Forbush decrease. The CME was not only responsible for affecting the rate of cosmic rays, but also caused anisotropies in their direction of incidence. Data from Project GRAND, an array of proportional wire chambers which detects secondary muons, are presented during the time of this Forbush decrease.     

Ground-level enhancements during solar cycle 23: results from SVIRCO,LOMNICKY STIT and LARC neutron monitors

The list of Ground-Level Enhancements (GLEs) occurring during the on-going solar cycle (no. 23) is used to show that they are absent during a time interval of the maximum phase in which a relative reduction in the sunspot area extent occurs. This period refers to October 2000 to March 2001 and it is associated with the concept of the Gnevyshev Gap, introduced 10 years ago by the Rome cosmic ray group. Also, effects of the Gnevyshev Gap in the solar cycle modulation of the atmospheric attenuation coefficient for Rome neutron monitor during the past solar cycle (no. 22) are illustrated. Moreover, it is suggested that GLE data files should be prepared by using the appropriate attenuation coefficient for each level of solar activity.     

Study of the 28 October 2003 and 20 January 2005 solar flares by means of 2.223 MeV gamma-emission line

By the data on intensity-time profiles of the neutron capture line of 2.223 MeV we have studied some characteristics of two solar flares, 28 October 2003 and 20 January 2005 (INTEGRAL and CORONAS-F observations, respectively). The SINP code was applied making allowance for the main processes of neutron interactions and deceleration in the solar plasma, character of neutron source, losses of neutrons and density model of the solar atmosphere. Comparison of the computed time profiles of 2.223 MeV line with observed ones for the flare of 28 October 2003 confirms the results obtained earlier for three other flares. Namely, the effect of density enhancement (EDE) in the sub-flare region, as well as the variations (hardening) of accelerated particle spectrum in the course of the event have been confirmed. The usual modeling procedure by the SINP code, however, seems to be inapplicable to the event of 20 January 2005. Possible causes of density enhancements during some flares and peculiarities of the 20 January 2005 flare are discussed.     

Energetic particle signatures of geoeffective coronal mass ejections

We have studied statistically associations of moderate and intense geomagnetic storms with coronal mass ejections (CMEs) and energetic particle events. The goal was to identify specific energetic particle signatures, which could be used to improve the predictions of the geoeffectiveness of full and partial halo CMEs. Protons in the range 1–110 MeV from the ERNE experiment onboard SOHO are used in the analysis. The study covers the time period from August 1996 to July 2000. We demonstrate the feasibility of energetic particle observations as an additional source of information in evaluating the geoeffectiveness of full and partial halo CMEs. Based on the observed onset times of solar energetic particle (SEP) events and energetic storm particle (ESP) events, we derive a proxy for the transit times of shocks driven by the interplanetary counterparts of coronal mass ejections from the Sun to the Earth. For a limited number of geomagnetic storms which can be associated to both SEP and ESP signatures, we found that this transit time correlates with the strength of geomagnetic storms.     

Low latitude nighttime ionospheric vertical E × B drifts at African region

The nighttime vertical E × B drifts velocities of the F2-region were inferred from the hourly hmF2 values obtained from ionosonde data over an African equatorial station, Ilorin (8.50^oN, 4.68^oE; dip lat. 2.95^o) during period of low solar activity. For each season, the plasma drift Vz is characterized by an evening upward enhancement, then by a downward reversal at 1900 LT till around 0000 LT, except for June solstice. This was explained using the Rayleigh–Taylor (R-T) instability mechanism. The occasional drift differences in Vz obtained by inferred and direct measurement over Ilorin and Jicamarca, respectively are reflective of the importance of chemistry and divergent transport system due to both the E region electric and magnetic fields instead of simple motions. The pre-reversal enhancement (PRE) magnitude is higher during the equinoctial months than the solsticial months over Jicamarca, highest during December solstice and the equinoctial months over Ilorin, suggesting the dominance of higher E × B fountain during equinoxes at both stations. The lowest PRE magnitude was in June solstice. The appearance of post-noon peak in NmF2 around 1700 LT is highest during the equinoctial months and lowest during the solsticial period. A general sharp drop in NmF2 around 1800 LT is distinct immediately after sunset, lowest during June solstice and highest in March equinox. Our result suggests that between 0930 and 2100 LT, the general theory that vertical drifts obtained by digisonde measurements only match the E × B drift if the F layer is higher than 300 km is reliable, but does not hold for the nighttime period of 2200–0600 LT under condition of solar minima. Hence, the condition may not be sufficient for the representation of vertical plasma drift at nighttime during solar minima. This assertion may still be tentative, as more equatorial stations needed to be studied for better confirmation.