Geoeffectiveness of interplanetary shocks during solar minimum (1995–1996) and solar maximum (2000)

Plasma and magnetic field parameter variations across fast forward interplanetary shocks are analyzed during the last solar cycle minimum (1995–1996, 15 shocks), and maximum year 2000 (50 shocks). It was observed that the solar wind velocity and magnetic field strength variation across the shocks were the parameters better correlated with Dst. Superposed epoch analysis centered on the shock showed that, during solar minimum, B _z profiles had a southward, long-duration variation superposed with fluctuations, whereas in solar maximum the B _z profile presented 2 peaks. The first peak occurred 4 hr after the shock, and seems to be associated with the magnetic field disturbed by the shock in the sheath region. The second peak occurred 19 hr after the shock, and seems to be associated with the ejecta fields. The difference in shape and peak in solar maximum (Dst peak =−50 nT, moderate activity) and minimum (Dst peak =−30 nT, weak activity) in average Dst profiles after shocks are, probably, a consequence of the energy injection in the magnetosphere being driven by different interplanetary southward magnetic structures. A statistical distribution of geomagnetic activity levels following interplanetary shocks was also obtained. It was observed that during solar maximum, 36% of interplanetary shocks were followed by intense (Dst≤−100 nT) and 28% by moderate (−50≤Dst <−100 nT) geomagnetic activity. During solar minimum, 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity, respectively. Thus, during solar maximum a higher relative number of interplanetary shocks might be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50–60% to be followed by intense/moderate geomagnetic activity.     

Statistical Comparison of Magnetic Clouds with Interplanetary Coronal Mass Ejections for Solar Cycle 23

We compare the number and characteristics of interplanetary coronal mass ejections (ICMEs) to those of magnetic clouds (MCs) by using in-situ solar wind plasma and magnetic field observations made at 1 AU during solar cycle 23. We found that ≈ 28% of ICMEs appear to contain MCs, since 103 magnetic clouds (MCs) occurred during 1995  – 2006, and 307 ICMEs occurred during 1996 – 2006. For the period between 1996 and 2006, 85 MCs are identified as part of ICMEs, and six MCs are not associated with ICMEs, which conflicts with the idea that MCs are usually a subset of ICMEs. It was also found that solar wind conditions inside MCs and ICMEs are usually similar, but the linear correlation between geomagnetic storm intensity (Dst _min) and relevant solar wind parameters is better for MCs than for ICMEs. The differences between average event duration (Δt) and average proton plasma β (〈β〉) are two of the major differences between MCs and ICMEs: i) the average duration of ICMEs (29.6 h) is 44% longer than for MCs (20.6 hours), and ii) the average of 〈β〉 is 0.01 for MCs and 0.24 for ICMEs. The difference between the definition of a MC and that for an ICME is one of the major reasons for these average characteristics being different (i.e., listed above as items i) and ii)), and it is the reason for the frequency of their occurrences being different.     

Coronal Holes and Solar Wind High-Speed Streams: II. Forecasting the Geomagnetic Effects

We present a simple method of forecasting the geomagnetic storms caused by high-speed streams (HSSs) in the solar wind. The method is based on the empirical correlation between the coronal hole area/position and the value of the Dst index, which is established in a period of low interplanetary coronal mass ejection (ICME) activity. On average, the highest geomagnetic activity, i.e., the minimum in Dst, occurs four days after a low-latitude coronal hole (CH) crosses the central meridian. The amplitude of the Dst dip is correlated with the CH area and depends on the magnetic polarity of the CH due to the Russell – McPherron effect. The Dst variation may be predicted by employing the expression Dst(t)=(−65±25×cos λ)[A(t ^*)]^0.5, where A(t ^*) is the fractional CH area measured in the central-meridian slice [−10°,10°] of the solar disc, λ is the ecliptic longitude of the Earth, ± stands for positive/negative CH polarity, and t−t ^*=4 days. In periods of low ICME activity, the proposed expression provides forecasting of the amplitude of the HSS-associated Dst dip to an accuracy of ≈30%. However, the time of occurrence of the Dst minimum cannot be predicted to better than ±2 days, and consequently, the overall mean relative difference between the observed and calculated daily values of Dst ranges around 50%.     

Periodic Appearance of Coronal Holes and the Related Variation of Solar Wind Parameters

We compared the variability of coronal hole (CH) areas (determined from daily GOES/SXI images) with solar wind (daily ACE data) and geomagnetic parameters for the time span 25 January 2005 until 11 September 2005 (late declining phase of solar cycle 23). Applying wavelet spectral analysis, a clear 9-day period is found in the CH time series. The GOES/SXI image sequence suggests that this periodic variation is caused by a mutual triangular distribution of CHs ∼120° apart in longitude. From solar wind parameters a 9-day periodicity was obtained as well, simultaneously with the 9-day period in the CH area time series. These findings provide strong evidence that the 9-day period in solar wind parameters, showing up as higher harmonic of the solar rotation frequency, is caused by the “periodic” longitudinal distribution of CHs on the Sun recurring for several solar rotations. The shape of the wavelet spectrum from the Dst index matches only weakly with that from the CH areas and is more similar to the wavelet spectrum of the solar wind magnetic field magnitude. The distinct 9-day period does not show up in sunspot group areas which gives further evidence that the solar wind modulation is strongly related to CH areas but not to active region complexes. The wavelet power spectra for the whole ACE data range (∼1998 – 2006) suggest that the 9-day period is not a singular phenomenon occurring only during a specific time range close to solar minimum but is occasionally also present during the maximum and decay phase of solar cycle 23. The main periods correspond to the solar rotation (27^d) as well as to the second (13.5^d) and third (9^d) harmonic. Electronic Supplementary Material The online version of this article () contains supplementary material, which is available to authorized users.     

The Dependence of the Geoeffectiveness of Interplanetary Flux Rope on Its Orientation,with Possible Application to Geomagnetic Storm Prediction

Interplanetary magnetic clouds (MCs) are one of the main sources of large non-recurrent geomagnetic storms. With the aid of a force-free flux rope model, the dependence of the intensity of geomagnetic activity (indicated by Dst index) on the axial orientation (denoted by θ and φ in GSE coordinates) of the magnetic cloud is analyzed theoretically. The distribution of the Dst values in the (θ, φ) plane is calculated by changing the axial orientation for various cases. It is concluded that (i) geomagnetic storms tend to occur in the region of θ<0°, especially in the region of θ≲−45°, where larger geomagnetic activity could be created; (ii) the intensity of geomagnetic activity varies more strongly with θ than with φ; (iii) when the parameters B ₀ (the magnetic field strength at the flux rope axis), R ₀ (the radius of the flux rope), or V (the bulk speed) increase, or |D| (the shortest distance between the flux rope axis and the x-axis in GSE coordinates) decreases, a flux rope not only can increase the intensity of geomagnetic activity, but also is more likely to create a storm, however the variation of n (the density) only has a little effect on the intensity; (iv) the most efficient orientation (MEO) in which a flux rope can cause the largest geomagnetic activity appears at φ∼0° or ∼ 180°, and some value of θ which depends mainly on D; (v) the minimum Dst value that could be caused by a flux rope is the most sensitive to changes in B ₀ and V of the flux rope, and for a stronger and/or faster MC, a wider range of orientations will be geoeffective. Further, through analyzing 20 MC-caused moderate to large geomagnetic storms during 1998 – 2003, a long-term prediction of MC-caused geomagnetic storms on the basis of the flux rope model is proposed and assessed. The comparison between the theoretical results and the observations shows that there is a close linear correlation between the estimated and observed minimum Dst values. This suggests that using the ideal flux rope to predict practical MC-caused geomagnetic storms is applicable. The possibility of the long-term prediction of MC-caused geomagnetic storms is discussed briefly.     

RING CURRENT DYNAMICS DURING THE 13–18 JULY 2000 STORM PERIOD

We study the development of the terrestrial ring current during the time interval of 13–18 July, 2000, which consisted of two small to moderate geomagnetic storms followed by a great storm with indices Dst=−300 nT and Kp=9. This period of intense geomagnetic activity was caused by three interplanetary coronal mass ejecta (ICME) each driving interplanetary shocks, the last shock being very strong and reaching Earth at ∼ 14 UT on 15 July. We note that (a) the sheath region behind the third shock was characterized by B _z fluctuations of ∼35 nT peak-to-peak amplitude, and (b) the ICME contained a negative to positive B _z variation extending for about 1 day, with a ∼ 6-hour long negative phase and a minimum B _z of about −55 nT. Both of these interplanetary sources caused considerable geomagnetic activity (Kp=8 to 9) despite their disparity as interplanetary triggers. We used our global ring current-atmosphere interaction model with initial and boundary conditions inferred from measurements from the hot plasma instruments on the Polar spacecraft and the geosynchronous Los Alamos satellites, and simulated the time evolution of H⁺, O⁺, and He⁺ ring current ion distributions. We found that the O⁺ content of the ring current increased after each shock and reached maximum values of ∼ 60% near minimum Dst of the great storm. We calculated the growth rate of electromagnetic ion cyclotron waves considering for the first time wave excitation at frequencies below O⁺ gyrofrequency. We found that the wave gain of O⁺ band waves is greater and is located at larger L shells than that of the He⁺ band waves during this storm interval. Isotropic pitch angle distributions indicating strong plasma wave scattering were observed by the imaging proton sensor (IPS) on Polar at the locations of maximum predicted wave gain, in good agreement with model simulations.     

Energy-Dependent Timing of Thermal Emission in Solar Flares

We report solar flare plasma to be multi-thermal in nature based on the theoretical model and study of the energy-dependent timing of thermal emission in ten M-class flares. We employ high-resolution X-ray spectra observed by the Si detector of the “Solar X-ray Spectrometer” (SOXS). The SOXS onboard the Indian GSAT-2 spacecraft was launched by the GSLV-D2 rocket on 8 May 2003. Firstly we model the spectral evolution of the X-ray line and continuum emission flux F(ε) from the flare by integrating a series of isothermal plasma flux. We find that the multi-temperature integrated flux F(ε) is a power-law function of ε with a spectral index (γ)≈−4.65. Next, based on spectral-temporal evolution of the flares we find that the emission in the energy range E=4 – 15 keV is dominated by temperatures of T=12 – 50 MK, while the multi-thermal power-law DEM index (δ) varies in the range of −4.4 and −5.7. The temporal evolution of the X-ray flux F(ε,t) assuming a multi-temperature plasma governed by thermal conduction cooling reveals that the temperature-dependent cooling time varies between 296 and 4640 s and the electron density (n _e) varies in the range of n _e=(1.77 – 29.3)×10¹⁰ cm⁻³. Employing temporal evolution technique in the current study as an alternative method for separating thermal from nonthermal components in the energy spectra, we measure the break-energy point, ranging between 14 and 21±1.0 keV.     

Evolution of Solar and Geomagnetic Activity Indices, and Their Relationship: 1960?�C?2001

We employ annually averaged solar and geomagnetic activity indices for the period 1960??C?2001 to analyze the relationship between different measures of solar activity as well as the relationship between solar activity and various aspects of geomagnetic activity. In particular, to quantify the solar activity we use the sunspot number R _s, group sunspot number R _g, cumulative sunspot area Cum, solar radio flux F10.7, and interplanetary magnetic field strength IMF. For the geomagnetic activity we employ global indices Ap, Dst and Dcx, as well as the regional geomagnetic index RES, specifically estimated for the European region. In the paper we present the relative evolution of these indices and quantify the correlations between them. Variations have been found in: i) time lag between the solar and geomagnetic indices; ii) relative amplitude of the geomagnetic and solar activity peaks; iii) dual-peak distribution in some of solar and geomagnetic indices. The behavior of geomagnetic indices is correlated the best with IMF variations. Interestingly, among geomagnetic indices, RES shows the highest degree of correlation with solar indices.     

Dynamic and Thermal Disappearance of Prominences and Their Geoeffectiveness

This paper is a qualitative study of 42 events of solar filament/prominence sudden disappearances (“disparitions brusques”; henceforth DBs) around two solar minima, 1985 – 1986 and 1994. The studied events were classified as 17 thermal and 25 dynamic disappearances. Associated events, i.e. coronal mass ejections (CMEs), type II bursts, evolution of nearby coronal holes, as well as solar wind speed, and geomagnetic disturbances are discussed. We have found that about 50% of the thermal DBs with adjacent (within 15° from the DB) coronal holes were associated with CMEs within a selected time window. All the studied thermal disappearances with adjacent coronal holes or accompanied by dynamic disappearances were associated with weak and medium geomagnetic storms. Also, nearly 64% of dynamic DBs were associated with CMEs. Ten (40%) dynamic disappearances were associated with intense geomagnetic storms, even when no CMEs was reported, six (24%) dynamic disappearances corresponded to extreme storms, and five (20%) corresponded to medium geomagnetic storms. The extreme geomagnetic storms appeared to be related to combined events, involving dynamic disappearances with adjacent coronal holes or including thermal disappearances. Furthermore, the geomagnetic activity (Dst index) increased if the source was close to the central meridian (±30°). The highest interplanetary magnetic field (B), longest duration, lowest southward direction B _z component, and lowest Dst were highly correlated for all studied events. The Sun – Earth transit time computed from the starting time of the sudden disappearance and the time its effect was measured at Earth was about 4.3 days and was mainly well correlated with the solar wind speed measured in situ (daily value).     

Variability in the Periodicity of 27 Days in Solar Indices

The short-term variation of solar indices, though typically near the solar rotation period of 27 days, can often deviate considerably from 27 days, in a wide range ∼ 19–33 days. The peak locations are within a day or two for all solar indices, indicating that the whole of the solar atmosphere is affected in a similar way. There are no systematic differences between the peaks of the chromosphere and the corona as such, but F10, X-rays, and coronal green line, which have uncertainties about their solar altitudes of origin, do show some differences (earlier peaks) as compared to other indices (chromospheric as well as coronal).     

Comparison of the Characteristics of Magnetic Clouds and Magnetic Cloud-Like Structures for the Events of 1995 – 2003

Using nine years of solar wind plasma and magnetic field data from the Wind mission, we investigated the characteristics of both magnetic clouds (MCs) and magnetic cloud-like structures (MCLs) during 1995 – 2003. A MCL structure is an event that is identified by an automatic scheme (Lepping, Wu, and Berdichevsky, Ann. Geophys. 23, 2687, 2005) with the same criteria as for a MC, but it is not usually identifiable as a flux rope by using the MC (Burlaga et al., J. Geophys. Res. 86, 6673, 1981) fitting model developed by Lepping, Jones, and Burlaga (Geophys. Res. Lett. 95(11), 957, 1990). The average occurrence rate is 9.5 for MCs and 13.6 for MCLs per year for the overall period of interest, and there were 82 MCs and 122 MCLs identified during this period. The characteristics of MCs and MCL structures are as follows: (1) The average duration, Δt, of MCs is 21.1 h, which is 40% longer than that for MCLs (Δt=15 h); (2) the average (minimum B _z found in MC/MCL measured in geocentric solar ecliptic coordinates) is −10.2 nT for MCs and −6 nT for MCLs; (3) the average Dst_min (minimum Dst caused by MCs/MCLs) is −82 nT for MCs and −37 nT for MCLs; (4) the average solar wind velocity is 453 km s⁻¹ for MCs and 413 km s⁻¹ for MCLs; (5) the average thermal speed is 24.6 km s⁻¹ for MCs and 27.7 km s⁻¹ for MCLs; (6) the average magnetic field intensity is 12.7 nT for MCs and 9.8 nT for MCLs; (7) the average solar wind density is 9.4 cm⁻³ for MCs and 6.3 cm⁻³ for MCLs; and (8) a MC is one of the most important interplanetary structures capable of causing severe geomagnetic storms. The longer duration, more intense magnetic field and higher solar wind speed of MCs, compared to those properties of the MCLs, are very likely the major reasons for MCs generally causing more severe geomagnetic storms than MCLs. But the fact that a MC is an important interplanetary structure with respect to geomagnetic storms is not new (e.g., Zhang and Burlaga, J. Geophys. Res. 93, 2511, 1988; Bothmer, ESA SP-535, 419, 2003).     

Relationship between the Parameters of Coronal Holes and High-Speed Solar Wind Streams over an Activity Cycle

The comparison of the brightness and area of coronal holes (CH) to the solar wind speed, which was started by Obridko et al. (Solar Phys. 260, 191, 2009a) has been continued. While the previous work was dealing with a relatively short time interval 2000 – 2006, here we have analyzed the data on coronal holes observed in the Sun throughout activity Cycle 23. A catalog of equatorial coronal holes has been compiled, and their brightness and area variations during the cycle have been analyzed. It is shown that CH is not merely an undisturbed zone between the active regions. The corona heating mechanism in CH seems to be essentially the same as in the regions of higher activity. The reduced brightness is the result of a specific structure with the magnetic field being quasi-radial at as low an altitude as 1.1R _⊙ or a bit higher. The plasma outflow decreases the measure of emission from CH. With an adequate choice of the photometric boundaries, the CH area and brightness indices display a fairly high correlation (0.6 – 0.8) with the solar wind velocity throughout the cycle, except for two years, which deviate dramatically – 2001 and 2007, i.e., the maximum and the minimum of the cycle. The mean brightness of the darkest part of CH, where the field lines are nearly radial at low altitudes, is of the order of 18 – 20% of the solar brightness, while the brightness of the other parts of the CH is 30 – 40%. The solar wind streams originate at the base of the coronal hole, which acts as an ejecting nozzle. The solar wind parameters in CH are determined at the level where the field lines are radial.     

Near-Earth Interplanetary Coronal Mass Ejections During Solar Cycle 23 (1996 – 2009): Catalog and Summary of Properties

In a previous study (Cane and Richardson, J. Geophys. Res. 108(A4), SSH6-1, 2003), we investigated the occurrence of interplanetary coronal mass ejections in the near-Earth solar wind during 1996 – 2002, corresponding to the increasing and maximum phases of solar cycle 23, and provided a “comprehensive” catalog of these events. In this paper, we present a revised and updated catalog of the ≈300 near-Earth ICMEs in 1996 – 2009, encompassing the complete cycle 23, and summarize their basic properties and geomagnetic effects. In particular, solar wind composition and charge state observations are now considered when identifying the ICMEs. In general, these additional data confirm the earlier identifications based predominantly on other solar wind plasma and magnetic field parameters. However, the boundaries of ICME-like plasma based on charge state/composition data may deviate significantly from those based on conventional plasma/magnetic field parameters. Furthermore, the much studied “magnetic clouds”, with flux-rope-like magnetic field configurations, may form just a substructure of the total ICME interval.     

Low-Energy Cutoffs in Electron Spectra of Solar Flares: Statistical Survey

The Reuven Ramaty High Energy Spectroscopic Imager (RHESSI) X-ray data base (February 2002 – May 2006) has been searched to find solar flares with weak thermal components and flat photon spectra. Using a regularized inversion technique, we determine the mean electron flux distribution from count spectra for a selection of events with flat photon spectra in the 15 – 20 keV energy range. Such spectral behavior is expected for photon spectra either affected by photospheric albedo or produced by electron spectra with an absence of electrons in a given energy range (e.g., a low-energy cutoff in the mean electron spectra of nonthemal particles). We have found 18 cases that exhibit a statistically significant local minimum (a dip) in the range of 13 – 19 keV. The positions and spectral indices of events with low-energy cutoff indicate that such features are likely to be the result of photospheric albedo. It is shown that if the isotropic albedo correction is applied, all low-energy cutoffs in the mean electron spectrum are removed, and hence the low-energy cutoffs in the mean electron spectrum of solar flares above ∼ 12 keV cannot be viewed as real features. If low-energy cutoffs exist in the mean electron spectra, their energies should be less than ∼ 12 keV.     

A Preliminary Estimate of the Size of the Coming Solar Cycle 24, based on Ohl’s Precursor Method

For many purposes (e.g., satellite drag, operation of power grids on Earth, and satellite communication systems), predictions of the strength of a solar cycle are needed. Predictions are made by using different methods, depending upon the characteristics of sunspot cycles. However, the method most successful seems to be the precursor method by Ohl and his group, in which the geomagnetic activity in the declining phase of a sunspot cycle is found to be well correlated with the sunspot maximum of the next cycle. In the present communication, the method is illustrated by plotting the 12-month running means aa(min ) of the geomagnetic disturbance index aa near sunspot minimum versus the 12-month running means of the sunspot number Rz near sunspot maximum [aa(min ) versus Rz(max )], using data for sunspot cycles 9 – 18 to predict the Rz(max ) of cycle 19, using data for cycles 9 – 19 to predict Rz(max ) of cycle 20, and so on, and finally using data for cycles 9 – 23 to predict Rz(max ) of cycle 24, which is expected to occur in 2011 – 2012. The correlations were good (∼+0.90) and our preliminary predicted Rz(max ) for cycle 24 is 142±24, though this can be regarded as an upper limit, since there are indications that solar minimum may occur as late as March 2008. (Some workers have reported that the aa values before 1957 would have an error of 3 nT; if true, the revised estimate would be 124±26.) This result of the precursor method is compared with several other predictions of cycle 24, which are in a very wide range (50 – 200), so that whatever may be the final observed value, some method or other will be discredited, as happened in the case of cycle 23.     

A New Forecasting Index for Solar Wind Velocity Based on EIT 284 Å Observations

Various solar wind forecasting methods have been developed during the past decade, such as the Wang?–?Sheeley model and the Hakamada?–?Akasofu?–?Fry Version 2 (HAFv2) model. Also, considerable correlation has been found between the solar wind speed v and the coronal hole (CH) area A _M on the visible side of the Sun, showing quantitative improvement of forecasting accuracy in low CME activity periods (e.g., Vr?nak, Temmer, and Veronig, Solar Phys. 240, 315, 2007a). Properties of lower layers of the solar atmosphere are good indications of the subsequent interplanetary and geomagnetic activities. We analyze the SOHO/EIT 284 Å images and construct a new forecasting factor (Pch) from the brightness of the solar EUV emission, and a good correlation is found between the Pch factor and the 3-day-lag solar wind velocity (v) probed by the ACE spacecraft. The main difference between the Pch and A _M factor is that Pch does not depend on the CH-boundary estimate and can reflect both the area and brightness of CH. A simple method of forecasting the solar wind speed near Earth in low CME activity periods is presented. Between Pch and v from 21 November until 26 December 2003, the linear correlation coefficient is R=0.89. For comparison we also analyze the data in the same period (DOY 25?–?125, 2005) as Vr?nak, Temmer, and Veronig (Solar Phys. 240, 315, 2007a), who used the CH areas A _M for predicting the solar wind parameters. In this period the correlation coefficient between Pch and v is R=0.70, whereas for A _M and v the correlation coefficient is R=0.62. The average relative difference between the calculated and the observed values is $\overline{|\delta|}\approx 12.15\%Various solar wind forecasting methods have been developed during the past decade, such as the Wang – Sheeley model and the Hakamada – Akasofu – Fry Version 2 (HAFv2) model. Also, considerable correlation has been found between the solar wind speed v and the coronal hole (CH) area A _M on the visible side of the Sun, showing quantitative improvement of forecasting accuracy in low CME activity periods (e.g., Vršnak, Temmer, and Veronig, Solar Phys. 240, 315, 2007a). Properties of lower layers of the solar atmosphere are good indications of the subsequent interplanetary and geomagnetic activities. We analyze the SOHO/EIT 284 ? images and construct a new forecasting factor (Pch) from the brightness of the solar EUV emission, and a good correlation is found between the Pch factor and the 3-day-lag solar wind velocity (v) probed by the ACE spacecraft. The main difference between the Pch and A _M factor is that Pch does not depend on the CH-boundary estimate and can reflect both the area and brightness of CH. A simple method of forecasting the solar wind speed near Earth in low CME activity periods is presented. Between Pch and v from 21 November until 26 December 2003, the linear correlation coefficient is R=0.89. For comparison we also analyze the data in the same period (DOY 25 – 125, 2005) as Vršnak, Temmer, and Veronig (Solar Phys. 240, 315, 2007a), who used the CH areas A _M for predicting the solar wind parameters. In this period the correlation coefficient between Pch and v is R=0.70, whereas for A _M and v the correlation coefficient is R=0.62. The average relative difference between the calculated and the observed values is . Furthermore, for the ten peaks during the analysis period, Pch and v show a correlation coefficient of R=0.78, and the average relative difference between the calculated and the observed peak values is . Moreover, the Pch factor can eliminate personal bias in the forecasting process, which existed in the method using CH area as input parameter, because CH area depends on the CH-boundary estimate but Pch does not. Until now the CH-boundary could not be easily determined since no quantitative criteria can be used to precisely locate CHs from observations, which led to differences in forecasting accuracy.     

The 156-day quasi-periodicity in solar indices

This paper investigates a series of daily solar indices: the sunspot number W (1900–2008), solar flux at 2800 MHz F _10.7 (1947–2008), and a number of X-ray flares N _x (1981–2008). The methods of Fourier and wavelet analysis are used to reveal the so-called 156-day Rieger-type periodicity (RTP). The W index is observed to have a statistically significant RTP amplitude in the neighborhood of the solar maxima in most of the solar cycles under study, except for cycles 14, 15, and 23. The 156-day peak is observed to have its largest power during the declining phase of cycle 16, at the maximum of cycle 21, and during the increasing phase of cycles 20 and 23. Statistically significant RTPs are also observed at the minima of cycles 17, 18 and 19. We conclude that there is no stable dependence between RTP and the solar cycle. The wavelet analysis shows that the pattern of the RTP time dependence for the F _10.7 index is almost identical to that of the W index. The correlation coefficient between the RTP curves is 0.95. The correlation coefficients for the pairs of indices W-N _x and F _10.7-N _x are 0.36 and 0.32, respectively. No time lags are found between the RTP starting points for different indices. Thus, the 156-day quasi-periodicity involves, almost simultaneously, events that occur in active regions of the solar atmosphere at different heights. This paper discusses the possible nature of RTP.     

Radial Evolution of Well-Observed Slow CMEs in the Distance Range 2 – 30 <Emphasis Type="Italic">R</Emphasis><Subscript>⊙</Subscript>

We performed a detailed analysis of 27 slow coronal mass ejections (CMEs) whose heights were measured in at least 30 coronagraphic images and were characterized by a high quality index (≥4). Our primary aim was to study the radial evolution of these CMEs and their properties in the range 2 – 30 solar radii. The instantaneous speeds of CMEs were calculated by using successive height – time data pairs. The obtained speed – distance profiles [v(R)] are fitted by a power law v = a(R−b) ^c . The power-law indices are found to be in the ranges a=30 – 386, b=1.95 – 3.92, and c=0.03 – 0.79. The power-law exponent c is found to be larger for slower and narrower CMEs. With the exception of two events that had approximately constant velocity, all events were accelerating. The majority of accelerating events shows a v(R) profile very similar to the solar-wind profile deduced by Sheeley et al. (Astrophys. J. 484, 472, 1997). This indicates that the dynamics of most slow CMEs are dominated by the solar wind drag.     

Relationships Among Magnetic Clouds, CMES, and Geomagnetic Storms

During solar cycle 23, 82 interplanetary magnetic clouds (MCs) were identified by the Magnetic Field Investigation (MFI) team using Wind (1995 – 2003) solar wind plasma and magnetic field data from solar minimum through the maximum of cycle 23. The average occurrence rate is 9.5 MCs per year for the overall period. It is found that some of the anomalies in the frequency of occurrence were during the early part of solar cycle 23: (i) only four MCs were observed in 1999, and (ii) an unusually large number of MCs (17 events) were observed in 1997, just after solar minimum. We also discuss the relationship between MCs, coronal mass ejections (CMEs), and geomagnetic storms. During the period 1996 – 2003, almost 8000 CMEs were observed by SOHO-LASCO. The occurrence frequency of MCs appears to be related neither to the occurrence of CMEs as observed by SOHO LASCO nor to the sunspot number. When we included “magnetic cloud-like structures” (MCLs, defined by Lepping, Wu, and Berdichevsky, 2005), we found that the occurrence of the joint set (MCs + MCLs) is correlated with both sunspot number and the occurrence rate of CMEs. The average duration of the MCL structures is ~40% shorter than that of the MCs. The MCs are typically more geoeffective than the MCLs, because the average southward field component is generally stronger and longer lasting in MCs than in MCLs. In addition, most severe storms caused by MCs/MCLs with Dst _min≤ −100 nT occurred in the active solar period.     

Stream Interactions and Interplanetary Coronal Mass Ejections at 5.3 AU near the Solar Ecliptic Plane

We have performed a survey of the characteristics of two types of large spatial-scale solar-wind structures, stream interaction regions (SIRs), and interplanetary coronal mass ejections (ICMEs), near 5.3 AU, using solar-wind observations from Ulysses. Our study is confined to the three aphelion passes of Ulysses, and also within ± 10° of the solar ecliptic plane, covering a part of 1992, 1997 – 1998, and 2003 – 2005, representing three slices of different phases of the solar activity cycle. Overall, there are 54 SIRs and 60 ICMEs in the survey. Many are merged in hybrid events, suggesting that they have undergone multiple interactions prior to reaching Jovian orbit. About 91% of SIRs occur with shocks, with 47% of such shocks being forward – reverse shock pairs. The solar-wind velocity sometimes stays constant or even decreases within the interaction region near 5.3 AU, in contrast with the gradual velocity increase during SIRs at 1 AU. Shocks are driven by 58% of ICMEs, with 94% of them being forward shocks. Some ICMEs seem to have multiple small flux ropes with different scales and properties. We quantitatively compare various properties of SIRs and ICMEs at 5.3 AU, and study their statistical distributions and variations with solar activity. The width, maximum dynamic pressure, and peak perpendicular pressure of SIRs all become larger than ICMEs. Dynamic pressure (P _dyn) is expected to be important for Jovian magnetospheric activity. We have examined the distributions of P _dyn of SIRs, ICMEs, and general solar wind, but these cannot explain the observed bimodal distribution of the location of the Jovian magnetopause. By comparing the properties of SIRs and ICMEs at 0.72, 1, and 5.3 AU, we find that the ICME expansion slows down significantly between 1 and 5.3 AU. Some transient and small streams in the inner heliosphere have merged into a single interaction region. Electronic Supplementary Material  The online version of this article () contains supplementary material, which is available to authorized users.