Organization of Energetic Particles by the Solar Wind Structure During the Declining to Minimum Phase of Solar Cycle 23

We investigate the organization of the low energy energetic particles (≤1 MeV) by solar wind structures, in particular corotating interaction regions (CIRs) and shocks driven by interplanetary coronal mass ejections, during the declining-to-minimum phase of Solar Cycle 23 from Carrington rotation 1999 to 2088 (January 2003 to October 2009). Because CIR-associated particles are very prominent during the solar minimum, the unusually long solar minimum period of this current cycle provides an opportunity to examine the overall organization of CIR energetic particles for a much longer period than during any other minimum since the dawn of the Space Age. We find that the particle enhancements associated with CIRs this minimum period recurred for many solar rotations, up to 30 at times, due to several high-speed solar wind streams that persisted. However, very few significant CIR-related energetic particle enhancements were observed towards the end of our study period, reflecting the overall weak high-speed streams that occurred at this time. We also contrast the solar minimum observations with the declining phase when a number of solar energetic particle events occurred, producing a mixed particle population. In addition, we compare the observations from this minimum period with those from the previous solar cycle. One of the main differences we find is the shorter recurrence rate of the high-speed solar wind streams (～10 solar rotations) and the related CIR energetic particle enhancements for the Solar Cycle 22 minimum period. Overall our study provides insight into the coexistence of different populations of energetic particles, as well as an overview of the large-scale organization of the energetic particle populations approaching the beginning of Solar Cycle 24.     

The Rotation Rates of Solar Magnetic Fields During Solar Cycles 21 – 23

Employing the synoptic maps of the photospheric magnetic fields from the beginning of solar cycle 21 to the end of 23, we first build up a time – longitude stackplot at each latitude between ±35°. On each stackplot there are many tilted magnetic structures clearly reflecting the rotation rates, and we adopt a cross-correlation technique to explore the rotation rates from these tilted structures. Our new method avoids artificially choosing magnetic tracers, and it is convenient for investigating the rotation rates of the positive and negative fields by omitting one kind of field on the stackplots. We have obtained the following results. i) The rotation rates of the positive and negative fields (or the leader and follower polarities, depending on the hemispheres and solar cycles) between latitudes ±35° during solar cycles 21–23 are derived. The reversal times of the leader and follower polarities are usually not consistent with the years of the solar minimum, nevertheless, at latitudes ±16°, the reversal times are almost simultaneous with them. ii) The rotation rates of the three solar cycles averaged over each cycle are calculated separately for the positive, negative and total fields. The latitude profiles of rotation of the positive and negative fields exhibit equatorial symmetries with each other, and those of the total fields lie between them. iii) The differences in rotation rates between the leader and follower polarities are obtained. They are very small near the equator, and increase as latitude increases. In the latitude range of 5° – 20°, these differences reach 0.05 deg day⁻¹, and the mean difference for solar cycle 22 is somewhat smaller than cycles 21 and 23 in these latitude regions. Then, the differences reduce again at latitudes higher than 20°.     

CRRES Measurements of Energetic Helium During the 1990–1991 Solar Maximum

During the 1990–1991 solar maximum, the CRRES satellite measured helium from 38 to 110 MeV n^–1, with isotopic resolution, during both solar quiet periods and a number of large solar flares, the largest of which were seen during March and June 1991. Helium differential energy spectra and isotopic ratios are analyzed and indicate that (1) the series of large solar energetic particle (SEP) events of 2–22 June display characteristics consistent with CME-driven interplanetary shock acceleration; (2) the SEP events of 23–28 March exhibit signatures of both CME-driven shock acceleration and impulsive SEP acceleration; (3) below about 60 MeV n^–1, the helium flux measured by CRRES is dominated by solar helium even during periods of least solar activity; (4) the solar helium below 60 MeV n^–1 is enriched in ³He, with a mean ³He/⁴He ratio of about 0.18 throughout most of the CRRES mission `quiet' periods; and (5) an association of this solar component with small CMEs occurring during the periods selected as solar `quiet' times.     

Differential Rotation of Strong Magnetic Flux During Solar Cycles 21 – 23

In the present investigation we measure the differential rotation of strong magnetic flux during solar cycles 21 – 23 with the method of wavelet transforms. We find that the cycle-averaged synodic rotation rate of strong magnetic flux can be written as ω=13.47−2.58sin ² θ or ω=13.45−2.06sin ² θ−1.37sin ⁴ θ, where θ is the latitude. They agree well with the results derived from sunspots. A north–south asymmetry of the rotation rate is found at high latitudes (28°<θ<40°). The strong flux in the southern hemisphere rotates faster than that in the northern hemisphere by 0.2 deg day⁻¹. The asymmetry continued for cycles 21 – 23 and may be a secular property.     

The “FIP Effect” and the Origins of Solar Energetic Particles and of the Solar Wind

We find that the element abundances in solar energetic particles (SEPs) and in the slow solar wind (SSW), relative to those in the photosphere, show different patterns as a function of the first ionization potential (FIP) of the elements. Generally, the SEP and SSW abundances reflect abundance samples of the solar corona, where low-FIP elements, ionized in the chromosphere, are more efficiently conveyed upward to the corona than high-FIP elements that are initially neutral atoms. Abundances of the elements, especially C, P, and S, show a crossover from low to high FIP at \({\approx}\,10~\mbox{eV}\) in the SEPs but \({\approx}\,14~\mbox{eV}\) for the solar wind. Naively, this seems to suggest cooler plasma from sunspots beneath active regions. More likely, if the ponderomotive force of Alfvén waves preferentially conveys low-FIP ions into the corona, the source plasma that eventually will be shock-accelerated as SEPs originates in magnetic structures where Alfvén waves resonate with the loop length on closed magnetic field lines. This concentrates FIP fractionation near the top of the chromosphere. Meanwhile, the source of the SSW may lie near the base of diverging open-field lines surrounding, but outside of, active regions, where such resonance does not exist, allowing fractionation throughout the chromosphere. We also find that energetic particles accelerated from the solar wind itself by shock waves at corotating interaction regions, generally beyond 1 AU, confirm the FIP pattern of the solar wind.     

An Extreme Solar Event of 20 January 2005: Properties of the Flare and the Origin of Energetic Particles

The famous extreme solar and particle event of 20 January 2005 is analyzed from two perspectives. Firstly, using multi-spectral data, we study temporal, spectral, and spatial features of the main phase of the flare, when the strongest emissions from microwaves up to 200 MeV gamma-rays were observed. Secondly, we relate our results to a long-standing controversy on the origin of solar energetic particles (SEP) arriving at Earth, i.e., acceleration in flares, or shocks ahead of coronal mass ejections (CMEs). Our analysis shows that all electromagnetic emissions from microwaves up to 2.22 MeV line gamma-rays during the main flare phase originated within a compact structure located just above sunspot umbrae. In particular, a huge (≈ 10⁵ sfu) radio burst with a high frequency maximum at 30 GHz was observed, indicating the presence of a large number of energetic electrons in very strong magnetic fields. Thus, protons and electrons responsible for various flare emissions during its main phase were accelerated within the magnetic field of the active region. The leading, impulsive parts of the ground-level enhancement (GLE), and highest-energy gamma-rays identified with π ⁰-decay emission, are similar and closely correspond in time. The origin of the π ⁰-decay gamma-rays is argued to be the same as that of lower-energy emissions, although this is not proven. On the other hand, we estimate the sky-plane speed of the CME to be 2 000 – 2 600 km s⁻¹, i.e., high, but of the same order as preceding non-GLE-related CMEs from the same active region. Hence, the flare itself rather than the CME appears to determine the extreme nature of this event. We therefore conclude that the acceleration, at least, to sub-relativistic energies, of electrons and protons, responsible for both the major flare emissions and the leading spike of SEP/GLE by 07 UT, are likely to have occurred nearly simultaneously within the flare region. However, our analysis does not rule out a probable contribution from particles accelerated in the CME-driven shock for the leading GLE spike, which seemed to dominate at later stages of the SEP event. S.N. Kuznetsov deceased 17 May 2007.     

3D Coronal Structures and Magnetic Field During the Total Solar Eclipse of 29 March 2006

The good quality of the observing sequence of about 60 photographs of the white-light corona taken during the total solar eclipse observations on 29 March 2006, in Al Sallum, Egypt, enable us to use a new method of image processing for enhancement of the fine structure of coronal phenomena. We present selected magnetic-field lines derived for different parameters of the extrapolation model. The coincidence of the observed coronal white-light fine structures and the computed field-line positions provides a 3D causal relationship between coronal structures and the coronal magnetic field.     

Observations of Double-Periodic X-Ray Emission in Interacting Systems of Solar Flare Loops

We investigate the M1.8 solar flare of 20 October 2002. The flare was accompanied by quasi-periodic pulsations (QPP) of both thermal and nonthermal hard X-ray emissions (HXR) observed by RHESSI in the 3?–?50 keV energy range. Analysis of the HXR time profiles in different energy channels made with the Lomb periodogram has indicated two statistically significant time periods of about 16 and 36 s. The 36 s QPP were observed only in the nonthermal HXR emission in the impulsive phase of the flare. The 16 s QPP were found in thermal and nonthermal HXR emission both in the impulsive and in the decay phases of the flare. Imaging analysis of the flare region, the determined time periods of the QPP, and the estimated physical parameters of the flare loops allowed us to interpret the observed QPP in terms of MHD oscillations excited in two spatially separated, but interacting systems of flaring loops.     

Solar Cycle Variation of Magnetic Flux Ropes in a Quasi-Static Coronal Evolution Model

The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free “potential-field” extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helicity due to flux emergence and large-scale photospheric motions. This helicity concentrates into twisted magnetic flux ropes, which may lose equilibrium and be ejected. Here, we consider how the magnetic structure predicted by this model – in particular the flux ropes – varies over the solar activity cycle, based on photospheric input data from six periods of cycle 23. The number of flux ropes doubles from minimum to maximum, following the total length of photospheric polarity inversion lines. However, the number of flux rope ejections increases by a factor of eight, following the emergence rate of active regions. This is broadly consistent with the observed cycle modulation of coronal mass ejections, although the actual rate of ejections in the simulation is about a fifth of the rate of observed events. The model predicts that, even at minimum, differential rotation will produce sheared, non-potential, magnetic structure at all latitudes.     

Spectral Flattening During Solar Radio Bursts At Cm–mm Wavelengths and the Dynamics of Energetic Electrons in a Flare Loop

Until recently, most of the information on particle acceleration processes in solar flares has been obtained from hard X-ray and cm-microwave observations. As a rule they provide information on electrons with energies below 300 keV. During recent years it became possible to measure the gamma-ray and millimeter radio emission with improved sensitivities. These spectral ranges carry information on much higher energy electrons. We studied the temporal and spectral behaviour of the radio burst emission at centimeter-millimeter wavelengths (8–50 GHz) by using the data from the patrol instruments of IAP (Bern University). We have analyzed more than 20 impulsive and long duration radio bursts (of 10 s to several 100 s duration).The main finding of the data analysis is the presence of spectral flattening throughout the bursts, which occurs always during the decay phase of flux peaks, at frequencies well above the spectral peak frequency and independently of burst duration. Furthermore, for some of the bursts, the flux maxima at higher frequencies are delayed. These findings can serve as evidence of the hardening of the electron spectrum at energies above some hundreds of keV during the decay phase of cm–mm flux peaks. As a most likely reason for such a hardening we consider Coulomb collisions of energetic electrons continuously injected and trapped in a flaring loop.     

Investigation of the Solar Differential Rotation of Compact Magnetic Elements for 1966 – 1986

The differential rotation of compact magnetic elements during activity cycles 20 and 21 (1966 – 1986) is studied by using solar synoptic charts. For each hemisphere the compact magnetic elements with the polarity of the circumpolar magnetic field have larger rotation rates than the elements with the opposite polarity. This difference in rotation rates is present during the whole cycle except during the polarity reversal of the circumpolar field.     

Corrections to Solar Thermal Structure when a Turbulent Magnetic Field is Included

Correction of non-ideal effect due to a magnetic fluctuating tensor is derived from the ideal MHD equations. The inclusion of a magnetic turbulent field leads to modifications of the hydrostatic equilibrium equation and thermodynamical variables such as the temperature T, the adiabatic exponent γ, the adiabatic temperature gradient △↓ad and the temperature gradient △↓. In particular, the modifications in the adiabatic and radiative temperature gradients will result in a change in the Schwarzchild criterion, hence in the location of the base of the convective zone. Incorporating the modifications, we construct a modified thermodynamical equilibrium structure of the Sun.     

Interpreting Helioseismic Structure Inversion Results of Solar Active Regions

Helioseismic techniques such as ring-diagram analysis have often been used to determine the subsurface structural differences between solar active and quiet regions. Results obtained by inverting the frequency differences between the regions are usually interpreted as the sound-speed differences between them. These in turn are used as a measure of temperature and magnetic-field strength differences between the two regions. In this paper we first show that the “sound-speed” difference obtained from inversions is actually a combination of sound-speed difference and a magnetic component. Hence, the inversion result is not directly related to the thermal structure. Next, using solar models that include magnetic fields, we develop a formulation to use the inversion results to infer the differences in the magnetic and thermal structures between active and quiet regions. We then apply our technique to existing structure inversion results for different pairs of active and quiet regions. We find that the effect of magnetic fields is strongest in a shallow region above 0.985R _⊙ and that the strengths of magnetic-field effects at the surface and in the deeper (r<0.98R _⊙) layers are inversely related (i.e., the stronger the surface magnetic field the smaller the magnetic effects in the deeper layers, and vice versa). We also find that the magnetic effects in the deeper layers are the strongest in the quiet regions, consistent with the fact that these are basically regions with weakest magnetic fields at the surface. Because the quiet regions were selected to precede or follow their companion active regions, the results could have implications about the evolution of magnetic fields under active regions.     

Magnetic Discontinuities in the Near-Earth Solar Wind: Evidence of In-Transit Turbulence or Remnants of Coronal Structure?

Fluctuations in the solar wind plasma and magnetic field are well described by the sum of two power law distributions. It has been postulated that these distributions are the result of two independent processes: turbulence, which contributes mainly to the smaller fluctuations, and crossing the boundaries of flux tubes of coronal origin, which dominates the larger variations. In this study we explore the correspondence between changes in the magnetic field with changes in other solar wind properties. Changes in density and temperature may result from either turbulence or coronal structures, whereas changes in composition, such as the alpha-to-proton ratio are unlikely to arise from in-transit effects. Observations spanning the entire ACE dataset are compared with a null hypothesis of no correlation between magnetic field discontinuities and changes in other solar wind parameters. Evidence for coronal structuring is weaker than for in-transit turbulence, with only ∼ 25% of large magnetic field discontinuities associated with a significant change in the alpha-to-proton ratio, compared to ∼ 40% for significant density and temperature changes. However, note that a lack of detectable alpha-to-proton signature is not sufficient to discount a structure as having a solar origin.     

Modeling of Local Magnetic Field Enhancements within Solar Flux Ropes

To model and study local magnetic-field enhancements in a solar flux rope we consider the magnetic field in its interior as a superposition of two linear (constant α) force-free magnetic-field distributions, viz. a global one, which is locally similar to a part of the cylinder, and a local torus-shaped magnetic distribution. The newly derived solution for a toroid with an aspect ratio close to unity is applied. The symmetry axis of the toroid and that of the cylinder may or may not coincide. Both the large and small radii of the toroid are set equal to the cylinder’s radius. The total magnetic field distribution yields a flux tube which has a variable diameter with local minima and maxima. In principle, this approach can be used for the interpretation and analysis of solar-limb observations of coronal loops.     

Average Thickness of Magnetosheath Upstream of Magnetic Clouds at 1 AU versus Solar Longitude of Source

Starting with a large number (N=100) of Wind magnetic clouds (MCs) and applying necessary restrictions, we find a proper set of N=29 to investigate the average ecliptic plane projection of the upstream magnetosheath thickness as a function of the longitude of the solar source of the MCs, for those cases of MCs having upstream shock waves. A few of the obvious restrictions on the full set of MCs are the need for there to exist a driven upstream shock wave, knowledge of the MC’s solar source, and restriction to only MCs of low axial latitudes. The analysis required splitting this set into two subsets according to average magnetosheath speed: slow/average (300 – 500 km s⁻¹) and fast (500 – 1100 km s⁻¹) speeds. Only the fast set gives plausible results, where the estimated magnetosheath thickness (ΔS) goes from 0.042 to 0.079 AU (at 1 AU) over the longitude sector of 0° (adjusted source-center longitude of the average magnetic cloud) to 40° off center (East or West), based on N=11 appropriate cases. These estimates are well determined with a sigma (σ) for the fit of 0.0055 AU, where σ is effectively the same as (chi-squared) for the appropriate quadratic fit. The associated linear correlation coefficient for ΔS versus |Longitude| was very good (c.c.=0.93) for the fast range, and ΔS at 60° longitude is extrapolated to be 2.7 times the value at 0°. For the slower speeds we obtain the surprising result that ΔS is typically more-or-less constant at 0.040±0.013 AU at all longitudes, indicating that the MC as a driver, when moving close to the normal solar wind speed, has little influence on magnetosheath thickness. In some cases, the correct choice between two candidate solar-source longitudes for a fast MC might be made by noting the value of the observed ΔS just upstream of the MC. Also, we point out that, for the 29 events, the average sheath speed was well correlated with the quantity ΔV[=(〈V _MC〉−〈V _UPSTREAM〉)], and also with both 〈V _MC〉 and 〈V _MC,T〉, where 〈V _MC〉 is the first one-hour average of the MC speed, 〈V _MC,T〉 is the average MC speed across the full MC, and 〈V _UPSTREAM〉 is a five-hour average of the solar wind speed just upstream of the shock.     

Deflection of Coronal Rays by Remote CMEs: Shock Wave or Magnetic Pressure?

We analyze five events of the interaction of coronal mass ejections (CMEs) with the remote coronal rays located up to 90° away from the CME as observed by the SOHO/LASCO C2 coronagraph. Using sequences of SOHO/LASCO C2 images, we estimate the kink propagation in the coronal rays during their interaction with the corresponding CMEs ranging from 180 to 920 km s⁻¹ within the interval of radial distances from 3 R _⊙ to 6 R _⊙. We conclude that all studied events do not correspond to the expected pattern of shock wave propagation in the corona. Coronal ray deflection can be interpreted as the influence of the magnetic field of a moving flux rope within the CME. The motion of a large-scale flux rope away from the Sun creates changes in the structure of surrounding field lines, which are similar to the kink propagation along coronal rays. The retardation of the potential should be taken into account since the flux rope moves at a high speed, comparable with the Alfvén speed.     

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.     

Phase Space Analysis: The Equilibrium of the Solar Magnetic Cycle

We present the results of a statistical study of the solar cycle based on the analysis of the superficial toroidal magnetic field component phase space. The magnetic field component used to create the embedded phase space was constructed from monthly sunspot number observations since 1750. The phase space was split into 32 sections (or time instants) and the average values of the orbits on this phase space were calculated (giving the most probable cycle). In this phase space it is shown that the magnetic field on the Sun’s surface evolves through a set of orbits that go around a mean orbit (i.e., the most probable magnetic cycle that we interpret as the equilibrium solution). It follows that the most probable cycle is well represented by a van der Pol oscillator limit curve (equilibrium solution), as can be derived from mean-field dynamo theory. This analysis also retrieves the empirical Gnevyshev – Ohl’s rule between the first and second parts of the solar magnetic cycle. The sunspot number evolution corresponding to the most probable cycle (in phase space) is presented.     

Short-Term Periodicities in Sunspot Activities During the Descending Phase of Solar Cycle 23

In the present study, the short-term periodicities in the daily data of the sunspot numbers and areas are investigated separately for the full disk, northern, and southern hemispheres during Solar Cycle 23 for a time interval from 1 January 2003 to 30 November 2007 corresponding to the descending and minimum phase of the cycle. The wavelet power spectrum technique exhibited a number of quasi-periodic oscillations in all the datasets. In the high frequency range, we find a prominent period of 22 – 35 days in both sunspot indicators. Other quasi-periods in the range of 40 – 60, 70 – 90, 110 – 130, 140 – 160, and 220 – 240 days are detected in the sunspot number time series in different hemispheres at different time intervals. In the sunspot area data, quasi-periods in the range of 50 – 80, 90 – 110, 115 – 130, 140 – 155, 160 – 190, and about 230 days were noted in different hemispheres within the time period of analysis. The present investigation shows that the well-known “Rieger periodicity” of 150 – 160 days reappears during the descending phase of Solar Cycle 23, but this is prominent mainly in the southern part of the Sun. Possible explanations of these observed periodicities are delivered on the basis of earlier results detected in photospheric magnetic field time series (Knaack, Stenflo, and Berdyugina in Astron. Astrophys. 438, 1067, 2005) and solar r-mode oscillations.