Solar Spicules: A Review of Recent Models and Targets for Future Observations – (Invited Review)

Since their discovery over 100 years ago, there have been many suggestions for the origin and development of solar spicules. Because the velocities of spicules are comparable to the sound and Alfvén speeds of the low chromosphere, linear theory cannot fully describe them. Consequently, detailed tests of theoretical ideas had to await the development of computing power that only became available during the 1970s. This work reviews theories for spicules and spicule-like features over approximately the past 25 years, with an emphasis on the models based on nonlinear numerical simulations. These models have given us physical insight into wave propagation in the solar atmosphere, and have helped elucidate how such waves, and associated shock waves, may be capable of creating motions and structures on magnetic flux tubes in the lower solar atmosphere. So far, however, it has been difficult to reproduce the most-commonly-quoted parameters for spicules with these models, using what appears to be the most suitable input parameters. A key impediment to developing satisfactory models has been the lack of reliable observational information, which is a consequence of the small angular size and transient lifetime of spicules. I close with a list of key observational questions to be addressed with space-based satellites, such as the currently operating TRACE satellite, and especially the upcoming Solar-B mission. Answers to these questions will help determine which, if any, of the current models correctly explains spicules.     

The coupling of solar convection and rotation – (Invited Review)

In recent years, helioseismology has provided an unprecedented look at the dynamics of the solar interior. These new insights have been accompanied by tremendous advances in high-performance computing technology, prompting increasingly sophisticated and realistic numerical models of solar convection. Among the most important helioseismic constraints on global-scale convection models is the mean differential rotation profile of the solar envelope, which is established by convection under the influence of rotation. The highly turbulent nature of solar convection makes this rotational influence difficult to determine and model. I will begin this review by discussing the solar rotation profile inferred from helioseismic measurements and various theoretical and numerical approaches to account for it. Computational constraints limited early numerical models to relatively laminar flow regimes but more recent investigations have begun to explore the distinct nature of turbulent convection. After a brief overview of empirical and numerical results on the related Rayleigh-Bernard system, I will outline the current state of numerical modeling of turbulent convection in rotating, stratified fluids, first in Cartesian and then in spherical geometries. The emphasis throughout will be on how rotation influences the structure, evolution, and transport processes of turbulent convection and what type of differential rotation can result.     

Time-Distance Inversion Methods and Results – (Invited Review)

The current interpretations of the travel-time measurements in quiet and active regions on the Sun are discussed. These interpretations are based on various approximations to the 3-D wave equation such as the Fermat principle for acoustic rays and the Born approximation. The ray approximation and its modifications have provided the first view of the 3-D structures and flows in the solar interior. However, more accurate and computationally efficient approximations describing the relation between the wave travel times and the internal properties are required to study the structures and flows in detail. Inversion of the large three-dimensional datasets is efficiently carried out by regularized iterative methods. Some results of time-distance inversions for emerging active regions, sunspots, meridional flows and supergranulation are presented. An active region which emerged on the solar disk in January 1998, was studied from SOHO/MDI for eight days, both before and after its emergence at the surface. The results show a complicated structure of the emerging region in the interior, and suggest that the emerging flux ropes travel very quickly through the depth range of our observations. The estimated speed of emergence is about 1.3 km s^–1. Tomographic images of a large sunspot reveal sunspot `fingers' - long narrow structures at a depth of about 4 Mm, which connect the sunspot with surrounding pores of the same polarity.     

Conditions for the Formation and Maintenance of Filaments – (Invited Review)

Observational conditions for the formation and maintenance of filaments are reviewed since 1989 in the light of recent findings on their structure, chirality, inferred magnetic topology, and mass flows. Recent observations confirm the necessary conditions previously cited: (1) their location at a boundary between opposite-polarity magnetic fields (2) a system of overlying coronal loops, (3) a magnetically-defined channel beneath, (4) the convergence of the opposite-polarity network magnetic fields towards their common boundary within the channel and (5) cancellation of magnetic flux at the common polarity boundary. Evidence is put forth for three additional conditions associated with fully developed filaments: (A) field-aligned mass flows parallel with their fine structure (B) a multi-polar background source of small-scale magnetic fields necessary for the formation of the filament barbs and (C) a handedness property known as chirality which requires them to be either of two types, dextral or sinistral. One-to-one relationships have been established between the chirality of filaments and the chirality of their filament channels and overlying coronal arcades. These findings reinforce earlier evidence that every filament magnetic field is separate from the magnetic field of the overlying arcade but both are parts of a larger magnetic field system. The larger system has at least quadrupolar footprints in the photosphere and includes the filament channel and subphotospheric magnetic fields, This ‘systems’ view of filaments and their environment enables new perspectives on why arcades and channels are invariable conditions for their existence. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005026814076     

Helioseismic Holography of Active-Region Subphotospheres – (Invited Review)

The development of solar acoustic holography has opened a major new diagnostic avenue in local helioseismology. It has revealed `acoustic moats' surrounding sunspots, `acoustic glories' surrounding complex active regions, and `acoustic condensations' suggesting the existence of significant seismic anomalies up to 20 Mm beneath active-region photospheres. Phase-sensitive seismic holography is now yielding high-resolution maps of sound travel-time anomalies caused by magnetic forces in the immediate subphotosphere, apparent thermal enhancements in acoustic moats, and Doppler signatures of subsurface flows. It has given us the first seismic images of a solar flare, and has uncovered a remarkable anomaly in the statistical distribution of seismic emission from acoustic glories. Seismic holography will probably give us the means for early detection of large active regions on the far-surface of the Sun, and possibly of deep subsurface activity as well. This powerful diagnostic now promises a new insight into the hydromechanical and thermal environments of the solar interior in the local perspective.     

Solar Magnetoconvection – (Invited Review)

In recent years the study of how magnetic fields interact with thermal convection in the Sun has made significant advances. These are largely due to the rapidly increasing computer power and its application to more physically relevant parameters regimes and to more realistic physics and geometry in numerical models. Here we present a survey of recent results following one line of investigations and discuss and compare the results of these with observed phenomena.     

Aspects of Three-Dimensional Magnetic Reconnection – (Invited Review)

In this review paper we discuss several aspects of magnetic reconnection theory, focusing on the field-line motions that are associated with reconnection. A new exact solution of the nonlinear MHD equations for reconnective annihilation is presented which represents a two-fold generalization of the previous solutions. Magnetic reconnection at null points by several mechanisms is summarized, including spine reconnection, fan reconnection and separator reconnection, where it is pointed out that two common features of separator reconnection are the rapid flipping of magnetic field lines and the collapse of the separator to a current sheet. In addition, a formula for the rate of reconnection between two flux tubes is derived. The magnetic field of the corona is highly complex, since the magnetic carpet consists of a multitude of sources in the photosphere. Progress in understanding this complexity may, however, be made by constructing the skeleton of the field and developing a theory for the local and global bifurcations between the different topologies. The eruption of flux from the Sun may even sometimes be due to a change of topology caused by emerging flux break-out. A CD-ROM attached to this paper presents the results of a toy model of vacuum reconnection, which suggests that rapid flipping of field lines in fan and separator reconnection is an essential ingredient also in real non-vacuum conditions. In addition, it gives an example of binary reconnection between a pair of unbalanced sources as they move around, which may contribute significantly to coronal heating. Finally, we present examples in TRACE movies of geometrical changes of the coronal magnetic field that are a likely result of large-scale magnetic reconnection. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005248007615     

Waves and Oscillations in the Corona – (Invited Review)

It has long been suggested on theoretical grounds that MHD waves must occur in the solar corona, and have important implications for coronal physics. An unequivocal identification of such waves has however proved elusive, though a number of events were consistent with an interpretation in terms of MHD waves. Recent detailed observations of waves in events observed by SOHO and TRACE removes that uncertainty, and raises the importance of MHD waves in the corona to a higher level. Here we review theoretical aspects of how MHD waves and oscillations may occur in a coronal medium. Detailed observations of waves and oscillations in coronal loops, plumes and prominences make feasible the development of coronal seismology, whereby parameters of the coronal plasma (notably the Alfvén speed and through this the magnetic field strength) may be determined from properties of the oscillations. MHD fast waves are refracted by regions of low Alfvén speed and slow waves are closely field-guided, making regions of dense coronal plasma (such as coronal loops and plumes) natural wave guides for MHD waves. There are analogies with sound waves in ocean layers and with elastic waves in the Earth's crust. Recent observations also indicate that coronal oscillations are damped. We consider the various ways this may be brought about, and its implications for coronal heating.     

Towards Understanding Solar Convection and Activity – (Invited Review)

The dynamics of the large-scale eddies which advect angular momentum through the convection zone is controlled in a significant way by the boundary conditions, which, if they are not modelled adequately, do not lead to a distribution of angular velocity that is consistent with observation. The transition boundary layer separating the convection zone from the radiative interior is thought to play a critical role in controlling the magnetic field in the convection zone, and is probably not wholly irrelevant to understanding the cycle of solar activity.     

Sunspot Oscillations: A Review – (Invited Review)

The current state of our knowledge, and ignorance, of the nature of oscillations in sunspots is surveyed. An effort is made to summarize the robust aspects of both the observational and theoretical components of the subject in a coherent, and common, conceptual framework. Detailed discussions of the various controversial issues are avoided except in instances where new viewpoints are advanced. Instead, extensive references are made to the growing literature on the subject, and generous explanatory remarks are made to guide the reader who wishes to delve more deeply into the underpinnings of the subject matter.     

Basic Principles of Solar Acoustic Holography – (Invited Review)

We summarize the basic principles of holographic seismic imaging of the solar interior, drawing on familiar principles in optics and parallels with standard optical holography. Computational seismic holography is accomplished by the phase-coherent wave-mechanical reconstruction of the p-mode acoustic field into the solar interior based on helioseismic observations at the solar surface. It treats the acoustic field at the solar surface in a way broadly analogous to how the eye treats electromagnetic radiation at the surface of the cornea, wave-mechanically refocusing radiation from submerged sources to render stigmatic images that can be sampled over focal surfaces at any desired depth. Holographic diagnostics offer a straight-forward assessment of the informational content of the observed p-mode spectrum independent of prospective physical models of the local interior anomalies that it represents. Computational holography was proposed as the optimum approach whereby to address the severe diffraction effects that confront standard tomography in the solar p-mode environment. It has given us a number of remarkable discoveries in the last two years and now promises a new insight into solar interior structure and dynamics in the local perspective. We compare the diagnostic roles of simple acoustic-power holography and phase-sensitive holography, and anticipate approaches to solar interior modeling based on holographic signatures. We identify simple computational principles that, applied to high-quality helioseismic observations, make it easy for prospective analysts to produce high-quality holographic images for practical applications in local helioseismology.     

A comparison of differential rotation measurements – (Invited Review)

Observers have long measured solar rotation with different techniques and obtained different results. This paper compares differential rotation measurements from four techniques: Doppler shift, Doppler feature tracking, magnetic feature tracking, and p-mode splittings. The different rotation rates measured by the first three techniques are interpreted as rotation rates of solar phenomena which depend on the properties and depth of that which is measured. This interpretation is supported by comparison with rotation measurements obtained from p-mode splittings except for Doppler features. The rotation rate of the plasma corresponds to the surface rate obtained by inversions; the rates of magnetic features correspond to the rotation rate at various depths within the convection zone. Supergranulation rotates at a rate greater than the maximum rotation rate within the convection zone, suggesting that supergranules are not simple convection cells anchored at a particular depth.     

Acoustic Imaging of Solar Active Regions – (Invited Review)

Acoustic imaging is a new method to construct the acoustic signal at a point on the solar surface or in the solar interior with the signals measured at the solar surface. The constructed signals contain both intensity information and phase information. The intensity is computed by summing the squared amplitude of the constructed signal over time. The phase of constructed signals can be studied by the cross-correlation function between the time series constructed with ingoing waves and outgoing waves. The location of the envelope peak of the cross-correlation function and the phase of the cross-correlation function contain different information on the physical conditions of the plasma along the wave path. From the constructed signals, one can form the two-dimensional outgoing intensity map, absorption map, phase-shift map, and envelope-shift map of a target region at different focal depths. The perturbed physical conditions caused by the magnetic fields of active regions manifest in these maps. The outgoing intensity is lower in magnetic regions than the quiet Sun. The group travel time and phase travel time are smaller in magnetic regions than in the quiet Sun. In this paper, we review the studies of active regions, including emerging flux regions, with acoustic imaging.     

Fluid Dynamics and MHD of the Solar Convection Zone and Tachocline: Current Understanding and Unsolved Problems – (Invited Review)

We review recent progress and define unanswered scientific questions in five related topics: granulation- to supergranulation-scale convection and magnetic structures; global convection and circulation; the rise of magnetic flux tubes to the photosphere, and their injection into the base of the convection zone; tachocline fluid dynamics and MHD; and the solar dynamo. We close with a set of observational `targets' for helioseismologists to aim for.     

Statistical Study of ICMEs and Their Sheaths During Solar Cycle 23 (1996 – 2008)

We have created a new catalog of 325 interplanetary coronal mass ejections (ICMEs) using their in-situ plasma signatures from 1996 to 2008; this time period includes Solar Cycle 23. The data set came from the OMNI near-Earth database. The one-minute resolution data that we used include magnetic-field strength, solar-wind speed, proton density, proton temperature, and plasma β. We compared this new catalog with other published catalogs. For every event, we indicated the presence of an ICME-driven shock. We identified the boundaries of ICMEs and their sheaths, and examined the statistical properties of characteristic parameters. We derived the duration and radial width of ICMEs and sheaths in the region near Earth. The statistical analysis of all events shows that, on average, sheaths travel faster than ICMEs, which indicates the expansion of CMEs in the interplanetary medium. They have higher mean magnetic-field strength values than ICMEs, and they are denser. They have higher mean proton temperature and plasma β than ICMEs, but they are smaller than ICMEs and last for a shorter time. The events were divided into different categories according to whether they included a shock and according to the phase of Solar Cycle 23 in which they are observed, i.e. ascending, maximum, or descending phase. We compared the different categories. We present a catalog of events available to the scientific community that studies ICMEs, and show the distribution and statistical properties of various parameters during these phenomena that govern the solar wind, the interplanetary medium, and space weather.     

Latitude Dependence of the Variations of Sunspot Group Numbers (SGN) and Coronal Mass Ejections (CMEs) in Cycle 23

The 12-month running means of the conventional sunspot number Rz, the sunspot group numbers (SGN) and the frequency of occurrence of Coronal Mass Ejections (CMEs) were examined for cycle 23 (1996 – 2006). For the whole disc, the SGN and Rz plots were almost identical. Hence, SGN could be used as a proxy for Rz, for which latitude data are not available. SGN values were used for 5° latitude belts 0° – 5°, 5° – 10°, 10° – 15°, 15° – 20°, 20° – 25°, 25° – 30° and > 30°, separately in each hemisphere north and south. Roughly, from latitudes 25° – 30° N to 20° – 25° N, the peaks seem to have occurred later for lower latitudes, from latitudes 20° – 25° N to 15° – 20° N, the peaks are stagnant or occur slightly earlier, and then from latitudes 15° – 20° N to 0° – 5° N, the peaks seem to have occurred again later for lower latitudes. Thus, some latitudinal migration is suggested, clearly in the northern hemisphere, not very clearly in the southern hemisphere, first to the equator in 1998, stagnant or slightly poleward in 1999, and then to the equator again from 2000 onwards, the latter reminiscent of the Maunder butterfly diagrams. Similar plots for CME occurrence frequency also showed multiple peaks (two or three) in almost all latitude belts, but the peaks were almost simultaneous at all latitudes, indicating no latitudinal migration. For similar latitude belts, SGN and CME plots were dissimilar in almost all latitude belts except 10° – 20° S. The CME plots had in general more peaks than the SGN plots, and the peaks of SGN often did not match with those of CME. In the CME data, it was noticed that whereas the values declined from 2002 to 2003, there was no further decline during 2003 – 2006 as one would have expected to occur during the declining phase of sunspots, where 2007 is almost a year of sunspot minimum. An inquiry at GSFC-NASA revealed that the person who creates the preliminary list was changed in 2004 and the new person picks out more weak CMEs. Thus a subjectivity (overestimates after 2002) seems to be involved and hence, values obtained before and during 2002 are not directly comparable to values recorded after 2002, except for CMEs with widths exceeding 60°.     

Spectral Analysis and Atmospheric Models of Microflares

By use of the high-resolution spectral data obtained with THEMIS on 2002 September 5, the spectra and characteristics of five well-observed microflares have been analyzed. Our results indicate that some of them are located near the longitudinal magnetic polarity inversion lines. All the microflares are accompanied by mass motions. The most obvious characteristic of the Ha microflare spectra is the emission at the center of both Ha and Can 8542 A lines. For the first time both thermal and non-thermal semi-empirical atmospheric models for the conspicuous and faint microflares are computed. In computing the non-thermal models, we assume that the electron beam resulting from magnetic reconnection is produced in the chromosphere, because it requires lower energies for the injected particles. It is found there is obvious heating in the low chromosphere. The temperature enhancement is about 1000-2200 K in the thermal models. If the non-thermal effects are included, then the required temperature increase can be reduced by 100-150 K. These imply that the Ha microflares can probably be produced by magnetic reconnection in the solar lower atmosphere. The radiative and kinetic energies of the Ha microflares are estimated and the total energy is found to be 1027 - 4×1028 erg.     

Interplanetary and Geomagnetic Consequences of Interacting CMEs of 13 – 14 June 2012

We report on the kinematics of two interacting CMEs observed on 13 and 14 June 2012. The two CMEs originated from the same active region NOAA 11504. After their launches which were separated by several hours, they were observed to interact at a distance of \(100~R_{\odot}\) from the Sun. The interaction led to a moderate geomagnetic storm at the Earth with minimum \(\mathrm{D}_{\mathrm{st}}\) index of approximately ?86 nT. The kinematics of the two CMEs is estimated using data from the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) instrument onboard the Solar Terrestrial Relations Observatory (STEREO). Assuming a head-on collision scenario, we find that the collision is inelastic in nature. Further, the signatures of their interaction are examined using the in situ observations obtained by Wind and the Advance Composition Explorer (ACE) spacecraft. It is also found that this interaction event led to the strongest sudden storm commencement (SSC) (\({\approx\,}150\) nT) of the present Solar Cycle 24. The SSC was of long duration, approximately 20 hours. The role of interacting CMEs in enhancing the geoeffectiveness is examined.