Source investigation of impulsive He-rich particle events

We have investigated the source characteristic and coronal magnetic field structure of six impulsive solar energetic particle (SEP) events selected from Wang et al. [Wang, Y.-M., Pick, M., Mason, G.M. Coronal holes, jets, and the origin of ³He-rich particle events. ApJ 639, 495, 2006] and Pick et al. [Pick, M., Mason, G.M., Wang, Y.-M., Tan, C., Wang, L. Solar source regions for ³He-rich solar energetic particle events identified using imaging radio, optical, and energetic particle observations. ApJ 648, 1247, 2006]. Some results are obtained: first, 2 events are associated with wide (≈100°) CMEs (hereafter wide CME events), another 4 events are associated with narrow (?40°) CMEs (hereafter narrow CME events); second, the coronal magnetic field configuration of narrow CME events appear more simple than that of the wide CME events; third, the photospheric magnetic field evolutions of all these events show new emergence of fluxes, while one case also shows magnetic flux cancellation; fourth, the EUV jets usually occurred very close to the footpoint of the magnetic field loop, while meter type III bursts occurred near or at the top of the loop and higher than EUV jets. Furthermore, the heights of type III bursts are estimated from the result of the coronal magnetic field extrapolations.     

Particle acceleration and gamma rays in solar flares: Recent observations and new modeling

Experiments on SMM, GAMMA, Yohkoh, GRANAT, Compton GRO, INTEGRAL, RHESSI and CORONAS-F satellites over the past three decades have provided copious data for fundamental research relating to particle acceleration, transport and energetics of flares and to the ambient abundance of the solar corona, chromosphere and photosphere. We summarize main results of solar gamma-astronomy (including some results of several joint Russian–Chinese projects) and try to appraise critically a real contribution of those results into modern understanding of solar flares, particle acceleration at the Sun and some properties of the solar atmosphere. Recent findings based on the RHESSI, INTEGRAL and CORONAS-F measurements (source locations, spectrum peculiarities, ³He abundance etc.) are especially discussed. Some unusual features of extreme solar events (e.g., 28 October 2003 and 20 January 2005) have been found in gamma-ray production and generation of relativistic particles (solar cosmic rays, or SCR). A number of different plausible assumptions are considered concerning the details of underlying physical processes during large flares: (1) existence of a steeper distribution of surrounding medium density as compared to a standard astrophysical model (HSRA) for the solar atmosphere; (2) enhanced content of the ³He isotope; (3) formation of magnetic trap with specific properties; (4) prevailing non-uniform (e.g., fan-like) velocity (angular) distributions of secondary neutrons, etc. It is emphasized that real progress in this field may be achieved only by combination of gamma-ray data in different energy ranges with multi-wave and energetic particle observations during the same event. We especially note several promising lines for the further studies: (1) resonant acceleration of the ³He ions in the corona; (2) timing of the flare evolution by gamma-ray fluxes in energy range above 90 MeV; (3) separation of gamma-ray fluxes from different sources at/near the Sun (e.g., different acceleration sources/episodes during the same flare, contribution of energetic particles accelerated by the CME-driven shocks etc.); (4) asymmetric magnetic geometry and new magnetic topology models of the near-limb flares; (5) modeling of self-consistent time scenario of the event.     

Propagation of energetic particles in the solar wind

The characteristics of solar energetic particle events as observed in interplanetary space depend on many physical processes acting at the source and during the transport from the source to the observer. These processes, such as acceleration at the Sun and the propagation near the Sun and in interplanetary space depend, in general, on both the particle velocity and rigidity. Thus, the evaluation of both the nuclear charge and/or atomic mass and the ionic charge of heavy ions turns out to be essential for the interpretation of the physical parameters observed, such as the energy spectra and the compositional variations during individual solar energetic particle events. In this paper recent results on the direct determination of the charge states of He, C, O, and Fe will be summarized. Using these results the compositional variations during individual solar particle events will be discussed. It will be shown that ratio changes by a factor of ～ 10 during the onset phase of solar particle events, as frequently observed, could be explained not only by rigidity dependent interplanetary propagation, but also by rigidity dependent diffusive propagation in the corona. However, there is now increasing experimental evidence that also other processes, such as compositional variations at the source and discontinuities of the interplanetary magnetic field, separating two different particle populations, may be important. Thus the picture emerges that these variations do not have a unique explanation but rather that each event has to be investigated individually.     

Temporal profiles of solar energetic particle events from SOHO/EPHIN data

Temporal profiles of energetic ions and electrons observed at 1 AU during solar energetic particle events are mainly determined by particle injection features, the observer location relative to the source region at the Sun, and the interplanetary space plasma and field conditions during particle transport. In this work, temporal profiles of 18 solar energetic particle events have been analyzed by fitting a pulse function to them in order to find a set of parameters which can be used to characterize the events.     

Asymmetric precipitation in a coronal loop as explanation of a singular observed spectrum

Almost 10 years of solar submillimeter observations have shown new aspects of solar activity, such as the presence of rapid solar spikes associated with the launch of coronal mass ejections and an increasing submillimeter spectral component in flares. We analyse the singular microwave–submillimeter spectrum of an M class solar flare on 20 December, 2002. Flux density observations measured by Sun patrol telescopes and the Solar Submillimeter Telescope are used to build the radio spectrum, which is fitted using Ramaty’s code. At submillimeter frequencies the spectrum shows a component different from the microwave classical burst. The fitting is achieved proposing two homogeneous sources of emission. This theoretical fitting is in agreement with differential precipitation through a magnetically asymmetric loop or set of loops. From a coronal magnetic field model we infer an asymmetric magnetic structure at the flare location. The model proposed to quantify the differential precipitation rates due to the asymmetry results in a total precipitation ratio Q₂/Q₁≈10⁴–10⁵, where Q₁(Q₂) represents the total precipitation in the loop foot with the high (low) magnetic field intensity. This ratio agrees with the electron total number ratio of the two sources proposed to fit the radio spectrum.     

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.     

Probing heliospheric diffusion coefficients with solar energetic particles

The dynamics of solar particle events provide a direct link to the understanding of properties of wave–particle interactions, and to the nature of the solar wind fluctuations. Depending on their energy, the often simultaneously observed electrons, protons and ions interact with different wavenumber ranges of the fluctuations, and are sensitive to various aspects of the dynamical nature of the solar wind turbulence. In general, the evolution of particle events is also sensitive to the spatial variation of the transport parameters between the Sun and a few AU. Together with in situ plasma and magnetic field observations this information can be used to extrapolate the properties of transport parameters into the more distant Heliosphere. Recent developments in the theory of parallel transport of energetic particles, and examples for the modelling of solar particle events and the derivation of transport parameters are considered. A dynamical quasi-linear theory is presented which gives special emphasis to the geometry and dynamic nature of the fluctuations, and which is able to provide particle mean free paths solely from observed plasma parameters, in good agreement with those derived by the modelling. Possibilities to apply the above results to the study of other energetic particle processes in the Heliosphere are discussed.     

SEP events and multi-spacecraft observations: Constraints on theory

Fifteen solar energetic particle (SEP) events have been analyzed using proton flux data recorded by the Helios 1, Helios 2, and IMP 8 spacecraft in the energy range ∼4–40 MeV during 1974–1982. For each of the events at least two of the spacecraft have their nominal magnetic footpoint within 20° in heliocentric longitude from each other. The SEP events are sub-grouped as a function of their heliocentric longitudinal separation and heliocentric radial distance from the SEP associated solar flare and several case studies are presented in this paper. Main results concerning their usage in estimating the SEP radial dependence are given. Moreover, we investigate the behavior of the third not connected spacecraft in order to study the dependence of the proton flux as a function of flare location. It is found that the contribution of the longitudinal gradient in determining variations in the SEP proton flux is particularly relevant for spacecraft having their magnetic connection footpoint separated from the flare between 30° and 50°.     

Relation of decline characteristics of 2–4.6 MeV protons in SEP events to solar wind parameters

The shape of the particle flux decline in solar energetic particle (SEP) events is of particular importance in understanding the propagation of energetic particles in the interplanetary medium. Power-law time profiles indicate the dominance of diffusive propagation, whereas exponential-law decline emphasizes convection transport and adiabatic deceleration. Values obtained theoretically for the decay time in the latter case are reasonably close to the fitted slopes in nearly half of all events when the solar wind speed stays nearly constant. Dependencies of characteristic decay time τ and spectral index γ on environmental plasma parameters are considered. Parts of exponential-law declines when solar wind speed: (a) decreases with time, (b) is constant, and (c) increases with time through the interval are analyzed separately. Both average values and dispersions of size distributions of τ for these three groups markedly differ in accordance with theoretical expectations.     

A question raised from the observation of dynamic cusp formation: When and where does particle acceleration occur?

We present observations of a C9.4 flare on 2002 June 2 in EUV (TRACE) and X-rays (RHESSI). The multiwavelength data reveal: (1) the involvement of a quadrupole magnetic configuration; (2) loop expansion and ribbon motion in the pre-impulsive phase; (3) gradual formation of a new compact loop with a long cusp at the top during the impulsive phase of the flare; (4) appearance of a large, twisted loop above the cusp expanding outward immediately after the hard X-ray peak; and (5) X-ray emission observed only from the new compact loop and the cusp. In particular, the gradual formation of an EUV cusp feature is very clear. The observations also reveal the timing of the cusp formation and particle acceleration: most of the impulsive hard X-rays (>25 keV) were emitted before the cusp was seen. This suggests that fast reconnection occurred during the restructuring of the magnetic configuration, resulting in more efficient particle acceleration, while the reconnection slowed after the cusp was completely formed and the magnetic geometry was stabilized. This observation is consistent with the observations obtained with Yohkoh/Soft X-ray Telescope (SXT) that soft X-ray cusp structures only appear after the major impulsive energy release in solar flares. These observations have important implications for the modeling of magnetic reconnection and particle acceleration.     

Power-law spectra of energetic electrons in solar flares from the maximum entropy and dimensional considerations

The maximum entropy formalism and dimensional analysis are used to derive a power-law spectrum of accelerated electrons in impulsive solar flares, where the particles can contain a significant fraction of the total flare energy. Entropy considerations are used to derive a power-law spectrum for a particle distribution characterised by its order of magnitude of energy. The derivation extends an earlier one-dimensional argument to the case of an isotropic three-dimensional particle distribution. Dimensional arguments employ the idea that the spectrum should reflect a balance between the processes of energy input into the corona and energy dissipation in solar flares. The governing parameters are suggested on theoretical grounds and shown to be consistent with solar flare observations. The flare electron flux, differential in the non-relativistic electron kinetic energy E, is predicted to scale as $E^{-3}$. This scaling is in agreement with RHESSI measurements of the hard X-ray flux that is generated by deka-keV electrons, accelerated in intense solar flares.     

太阳宇宙线源物质的高色球层模式和重离子丰度过量机制的探讨

本文利用太阳能量粒子事件中重离子平均丰度过量的资料,计算得到太阳能量粒子源物质的温度,提出了描述太阳宇宙线能量粒子源物质的新模式——高色球层模式;太阳耀斑观测确定,太阳宇宙线耀斑的加速区一般最可能出现在低日冕甚至高达几万公里的高度,从而,太阳宇宙线的源和加速区通常不位于同一区域;进而提出了描述太阳能量粒子事件中重离子丰度过量的可能机制——其源物质是通过太阳黑子的冻结型无力场从高色球层输送到活动区,形成耀斑前加速区内重离子丰度大和耀斑后宇宙线中重元素丰度的过量.      

Model of a fluxtube with a twisted magnetic field in the stratified solar atmosphere

We build a single vertical straight magnetic fluxtube spanning the solar photosphere and the transition region which does not expand with height. We assume that the fluxtube containing twisted magnetic fields is in magnetohydrostatic equilibrium within a realistic stratified atmosphere subject to solar gravity. Incorporating specific forms of current density and gas pressure in the Grad–Shafranov equation, we solve the magnetic flux function, and find it to be separable with a Coulomb wave function in radial direction while the vertical part of the solution decreases exponentially. We employ improved fluxtube boundary conditions and take a realistic ambient external pressure for the photosphere to transition region, to derive a family of solutions for reasonable values of the fluxtube radius and magnetic field strength at the base of the axis that are the free parameters in our model. We find that our model estimates are consistent with the magnetic field strength and the radii of Magnetic bright points (MBPs) as estimated from observations. We also derive thermodynamic quantities inside the fluxtube.     

How far are we from a ‘Standard Model’ of the solar dynamo?

Over the last few years, dynamo theorists seem to be converging on a basic scenario as to how the solar dynamo operates. The strong toroidal component of the magnetic field is produced in the tachocline, from where it rises due to magnetic buoyancy to produce active regions at the solar surface. The decay of tilted bipolar active regions at the surface gives rise to the poloidal component, which is first advected poleward by the meridional circulation and then taken below the surface to the tachocline where it can be stretched to produce the toroidal component. The mathematical formulation of this basic model, however, involves the specification of some parameters which are still uncertain. We review these remaining uncertainties which have resulted in disagreements amongst various research groups and have made it impossible to still arrive at something that can be called a standard model of the solar dynamo.     

Statistical properties of the most powerful solar and heliospheric disturbances

We present and discuss here the first version of a data base of extreme solar and heliospheric events. The data base contains now 87 extreme events mostly since 1940. An event is classified as extreme if one of the three critical parameters passed a lower limit. The critical parameters were the X-ray flux (parameter R), solar proton flux (parameter S) and geomagnetic disturbance level (parameter G). We find that the five strongest extreme events based on four variables (X-rays SEP, Dst, Ap) are completely separate except for the October 2003 event which is one the five most extreme events according to SEP, Dst and Ap. This underlines the special character of the October 2003 event, making it unique within 35 years. We also find that the events based on R and G are rather separate, indicating that the location of even extreme flares on the solar disk is important for geomagnetic effects. We also find that S = 3 events are not extreme in the same sense as R > 3 and G > 3 events, while S = 5 events are missing so far. This suggests that it might be useful to rescale the classification of SEP fluxes.     

PPS-87: a new event oriented solar proton prediction model.

A new event-oriented solar proton prediction model has been developed and implemented at the USAF Space Environment forecast facility. This new model generates predicted solar proton time-intensity profiles for a number of user adjustable energy ranges and is also capable of making predictions for the heavy ion flux. The computer program is designed so a forecaster can select inputs based on the data available in near real-time at the forecast center as the solar flare is occurring. The predicted event amplitude is based on the electromagnetic emission parameters of the solar flare (either microwave or soft X-ray emission) and the solar flare position on the sun. The model also has an update capability where the forecaster can normalize the prediction to actual spacecraft observations of spectral slope and particle flux as the event is occurring in order to more accurately predict the future time-intensity profile of the solar particle flux. Besides containing improvements in the accuracy of the predicted energetic particle event onset time and magnitude, the new model converts the predicted solar particle flux into an expected radiation dose that might be experienced by an astronaut during EVA activities or inside the space shuttle.     

Mode selection mechanism and magnetic diffusivity in a non-axisymmetric solar dynamo model

The generation of solar non-axisymmetric magnetic fields is studied based on a linear α²–Ω dynamo model in a rotating spherical frame. The model consists of a solar-like differential rotation, a magnetic diffusivity varied with depth, and three types of α-effects with different locations, i.e. the tachocline, the whole convective zone and the sub-surface. Some comparisons of the critical α-values of axisymmetric (m = 0) and longitude-dependent modes (m = 1,2,3) are presented to show the roles of the magnetic diffusivity in the problem of modes selection. With the changing of diffusivity intensity for the given solar differential rotation system, the dominant mode possibly changes likewise and the stronger the diffusivity is, the easier the non-axisymmetric modes are excited. The influence of the diffusivity and differential rotation on the configurations of the dominant modes are also presented.     

Modeling the effects of pitch-angle scattering processes on the transport of solar energetic particles along the interplanetary magnetic field

We present a Monte-Carlo technique to study the time-dependent transport of energetic particles in the interplanetary medium. We use the guiding center approximation between discrete finite pitch-angle scatterings to quantify the competing effects of focusing and pitch-angle scattering on energetic particles propagating along a Parker spiral magnetic field. We consider that the pitch-angle scattering process is produced by small-scale magnetic field irregularities frozen in the expanding solar wind. We also include the effects of both solar wind convection and adiabatic deceleration. We use a joint probability distribution P(h, μ′) = p(h; μ′)q(μ′; μ) to describe both the distance traveled by the particle between two scattering processes and the change in the particle pitch-angle after a scattering process. Here, p(h; μ′) is the conditional probability that the particle travels a distance h along the field line before the next scattering if it had a pitch-angle cosine μ′ after the previous scattering, and q(μ′; μ) is the conditional probability for the pitch-angle cosine μ′ if the pitch-angle cosine was μ before the scattering. We consider several functional forms to describe the processes of pitch-angle scattering, such as an isotropic scattering without any memory of the initial particle’s pitch-angle or processes in which the scattering result depends upon the initial particle’s pitch-angle. The results of our simulations are pitch-angle distributions and time-intensity profiles that can be directly compared to spacecraft observations. Comparison of our simulations with near-relativistic (45–290 keV) electron events observed by the Electron, Proton and Alpha Monitor on board the Advanced Composition Explorer allows us to estimate both the time dependence of the injection of near-relativistic electrons into the interplanetary medium and the conditions for electron propagation along the interplanetary magnetic field.     

Modeling of “gradual” solar energeticparticle events using a stochastic differential equation method

We have modeled “gradual” solar energetic particle events through numerical simulations using a StochasticDifferential Equation (SDE) method. We consider that energetic particle events are roughly divided into two groups: (1) where the shock was driven by coronal mass ejections (CMEs) associated with large solar flares, and (2) where they have no related solar events apart from the CMEs. (The detailed classification of energetic particle events was discussed in our previous paper.) What we call “gradual” solar energetic particle events belong to the former group. Particles with energies greater than 10 MeV are observed within several hours after the occurrence of flares and CMEs in many gradual events. By applying the SDE method coupled with particle splitting to diffusive acceleration, we found that an injection of high energy particles is necessary for early enhancement of such a high-energy proton flux and that it should not be presumed that the solar wind particles act as the seed population.     

Variability of spectra in large solar energeticparticle events

Using data from the Solar Isotope Spectrometer on the Advanced Composition Explorer obtained during 36large solar energetic particle events (SEPs) that occurred during 1997–2002 we have examined the spectral characteristics of oxygen and iron. Based on the shape of the oxygen spectrum during the decay phase following the peak in particle intensity, each SEP event was categorized as either exponential (7 events) or power law (29 events). We find that the exponential events were typically the larger events (in terms of peak oxygen intensity) and had Fe/0 ratios that strongly decreased with increasing energy.Event-averaged Fe/0 ratios (integrated over 12 to 60 MeV/nucleon) were at or below coronal abundances for nearly all these events, while the ratios obtained in the power law events were typically enhanced over coronal values. The majority of the power law events had the same spectral index for both 0 and Fe resulting in an Fe/0 ratio independent of energy. However 6 of the 29 power law events had Fe/0 ratios that increased with increasing energy due to an Fe spectral index less negative than that of 0. We consider simple diffusion theory in an effort to understand the nature of these events and obtain preliminary but promising results.