On the application of the boundary element method in coronal magnetic field reconstruction

Solar magnetic field is believed to play a central role in solar activities and flares, filament eruptions as well as CMEs are due to the magnetic field re-organization and the interaction between the plasma and the field. At present the reliable magnetic field measurements are still confined to a few lower levels like in photosphere and chromosphere. Although IR technique may be applied to observe the coronal field but the technique is not well-established yet. Radio techniques may be applied to diagnose the coronal field but assumptions on radiation mechanisms and propagations are needed. Therefore extrapolation from photospheric data upwards is still the primary method to reconstruction coronal field. Potential field has minimum energy content and a force-free field can provide the required excess energy for energy release like flares, etc. Linear models have undesirable properties and it is expected to consider non-constant-alpha force-free field model. As the recent result indicates that the plasma beta is sandwich-ed distributed above the solar surface (Gary, 2001), care must be taken in modeling the coronal field correctly. As the reconstruction of solar coronal magnetic fields is an open boundary problem, it is desired to apply some technique that can incorporate this property. The boundary element method is a well-established numerical techniques that has been applied to many fields including open-space problems. It has also been applied to solar magnetic field problems for potential, linear force-free field and non-constant-alpha force-free field problems. It may also be extended to consider the non-force-free field problem. Here we introduce the procedure of the boundary element method and show its applications in reconstruction of solar magnetic field problems. This revised version was published online in August 2006 with corrections to the Cover Date.     

Magnetic Flux Emergence, Activity, Eruptions and Magnetic Clouds: Following Magnetic Field from the Sun to the Heliosphere

We present an overview of how the principal physical properties of magnetic flux which emerges from the toroidal fields in the tachocline through the turbulent convection zone to the solar surface are linked to solar activity events, emphasizing the effects of magnetic field evolution and interaction with other magnetic structures on the latter. We compare the results of different approaches using various magnetic observables to evaluate the probability of flare and coronal mass ejection (CME) activity and forecast eruptive activity on the short term (i.e. days). Then, after a brief overview of the observed properties of CMEs and their theoretical models, we discuss the ejecta properties and describe some typical magnetic and composition characteristics of magnetic clouds (MCs) and interplanetary CMEs (ICMEs). We review some individual examples to clarify the link between eruptions from the Sun and the properties of the resulting ejecta. The importance of a synthetic approach to solar and interplanetary magnetic fields and activity is emphasized.     

Coronal Observations of CMEs

CMEs have been observed for over 30 years with a wide variety of instruments. It is now possible to derive detailed and quantitative information on CME morphology, velocity, acceleration and mass. Flares associated with CMEs are observed in X-rays, and several different radio signatures are also seen. Optical and UV spectra of CMEs both on the disk and at the limb provide velocities along the line of sight and diagnostics for temperature, density and composition. From the vast quantity of data we attempt to synthesize the current state of knowledge of the properties of CMEs, along with some specific observed characteristics that illuminate the physical processes occurring during CME eruption. These include the common three-part structures of CMEs, which is generally attributed to compressed material at the leading edge, a low-density magnetic bubble and dense prominence gas. Signatures of shock waves are seen, but the location of these shocks relative to the other structures and the occurrence rate at the heights where Solar Energetic Particles are produced remains controversial. The relationships among CMEs, Moreton waves, EIT waves, and EUV dimming are also cloudy. The close connection between CMEs and flares suggests that magnetic reconnection plays an important role in CME eruption and evolution. We discuss the evidence for reconnection in current sheets from white-light, X-ray, radio and UV observations. Finally, we summarize the requirements for future instrumentation that might answer the outstanding questions and the opportunities that new space-based and ground-based observatories will provide in the future.     

The Pre-CME Sun

The coronal mass ejection (CME) phenomenon occurs in closed magnetic field regions on the Sun such as active regions, filament regions, transequatorial interconnection regions, and complexes involving a combination of these. This chapter describes the current knowledge on these closed field structures and how they lead to CMEs. After describing the specific magnetic structures observed in the CME source region, we compare the substructures of CMEs to what is observed before eruption. Evolution of the closed magnetic structures in response to various photospheric motions over different time scales (convection, differential rotation, meridional circulation) somehow leads to the eruption. We describe this pre-eruption evolution and attempt to link them to the observed features of CMEs. Small-scale energetic signatures in the form of electron acceleration (signified by nonthermal radio bursts at metric wavelengths) and plasma heating (observed as compact soft X-ray brightening) may be indicative of impending CMEs. We survey these pre-eruptive energy releases using observations taken before and during the eruption of several CMEs. Finally, we discuss how the observations can be converted into useful inputs to numerical models that can describe the CME initiation.     

CME Theory and Models

This chapter provides an overview of current efforts in the theory and modeling of CMEs. Five key areas are discussed: (1) CME initiation; (2) CME evolution and propagation; (3) the structure of interplanetary CMEs derived from flux rope modeling; (4) CME shock formation in the inner corona; and (5) particle acceleration and transport at CME driven shocks. In the section on CME initiation three contemporary models are highlighted. Two of these focus on how energy stored in the coronal magnetic field can be released violently to drive CMEs. The third model assumes that CMEs can be directly driven by currents from below the photosphere. CMEs evolve considerably as they expand from the magnetically dominated lower corona into the advectively dominated solar wind. The section on evolution and propagation presents two approaches to the problem. One is primarily analytical and focuses on the key physical processes involved. The other is primarily numerical and illustrates the complexity of possible interactions between the CME and the ambient medium. The section on flux rope fitting reviews the accuracy and reliability of various methods. The section on shock formation considers the effect of the rapid decrease in the magnetic field and plasma density with height. Finally, in the section on particle acceleration and transport, some recent developments in the theory of diffusive particle acceleration at CME shocks are discussed. These include efforts to combine self-consistently the process of particle acceleration in the vicinity of the shock with the subsequent escape and transport of particles to distant regions.     

The Localization of Particle Acceleration Sites in Solar Flares and CMES

We review the particular aspect of determining particle acceleration sites in solar flares and coronal mass ejections (CMEs). Depending on the magnetic field configuration at the particle acceleration site, distinctly different radiation signatures are produced: (1) If charged particles are accelerated along compact closed magnetic field lines, they precipitate to the solar chromosphere and produce hard X-rays, gamma rays, soft X-rays, and EUV emission; (2) if they are injected into large-scale closed magnetic field structures, they remain temporarily confined (or trapped) and produce gyrosynchrotron emission in radio and bremsstrahlung in soft X-rays; (3) if they are accelerated along open field lines they produce beam-driven plasma emission with a metric starting frequency; and (4) if they are accelerated in a propagating CME shock, they can escape into interplanetary space and produce beam-driven plasma emission with a decametric starting frequency. The latter two groups of accelerated particles can be geo-effective if suitably connected to the solar west side. Particle acceleration sites can often be localized by modeling the magnetic topology from images in different wavelengths and by measuring the particle velocity dispersion from time-of-flight delays.     

The Role of Magnetic Reconnection in CME Acceleration

Observations carried out from the coronagraphs on board space missions (LASCO/SOHO, Solar Maximum and Skylab) and ground-based facilities (HAO/Mauna Loa Observatory) show that coronal mass ejections (CMEs) can be classified into two classes based on their kinematics evolution. These two classes of CMEs are so-called fast and slow CMEs. The fast CME starts with a high initial speed that remains more or less constant; it is also called the constant-speed CME. On the other hand, the slow CME starts with a low initial speed, but shows a gradual acceleration; it is also called the accelerated and slow CME. Low and Zhang [Astrophys. J. 564, L53–L56, 2002] suggested that these two classes of CMEs could be a result of a difference in the initial topology of the magnetic fields associated with the underlying quiescent prominences. A normal prominence magnetic field topology will lead to a fast CME, while an inverse quiescent prominence results in a slow CME, because of the nature of the magnetic reconnection processes. In a recent study given by Wu et al. [Solar Phys. 225, 157–175, 2004], it was shown that an inverse quiescent prominence magnetic topology also could produce a fast CME. In this study, we perform a numerical MHD simulation for CMEs occurring in both normal and inverse quiescent prominence magnetic topology. This study demonstrates three major physical processes responsible for destabilization of these two types of prominence magnetic field topologies that can launch CMEs. These three initiation processes are identical to those used by Wu et al. [Solar Phys. 225, 157–175, 2004]. The simulations show that both fast and slow CMEs can be initiated from these two different types of magnetic topologies. However, the normal quiescent prominence magnetic topology does show the possibility for launching a reconnection island (or secondary O-line) that might be thought of as a “CME’’.     

What Are Special About Ground-Level Events?

Ground level events (GLEs) occupy the high-energy end of gradual solar energetic particle (SEP) events. They are associated with coronal mass ejections (CMEs) and solar flares, but we still do not clearly understand the special conditions that produce these rare events. During Solar Cycle 23, a total of 16 GLEs were registered, by ground-based neutron monitors. We first ask if these GLEs are clearly distinguishable from other SEP events observed from space. Setting aside possible difficulties in identifying all GLEs consistently, we then try to find observables which may unmistakably isolate these GLEs by studying the basic properties of the associated eruptions and the active regions (ARs) that produced them. It is found that neither the magnitudes of the CMEs and flares nor the complexities of the ARs give sufficient conditions for GLEs. It is possible to find CMEs, flares or ARs that are not associated with GLEs but that have more extreme properties than those associated with GLEs. We also try to evaluate the importance of magnetic field connection of the AR with Earth on the detection of GLEs and their onset times. Using the potential field source surface (PFSS) model, a half of the GLEs are found to be well-connected. However, the GLE onset time with respect to the onset of the associated flare and CME does not strongly depend on how well-connected the AR is. The GLE onset behavior may be largely determined by when and where the CME-driven shock develops. We could not relate the shocks responsible for the onsets of past GLEs with features in solar images, but the combined data from the Solar TErrestrial RElations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) have the potential to change this for GLEs that may occur in the rising phase of Solar Cycle 24.     

Multi-Wavelength Observations of CMEs and Associated Phenomena

This chapter reviews how our knowledge of CMEs and CME-associated phenomena has been improved, since the launch of the SOHO mission, thanks to multi-wavelength analysis. The combination of data obtained from space-based experiments and ground based instruments allows us to follow the space-time development of an event from the bottom of the corona to large distances in the interplanetary medium. Since CMEs originate in the low solar corona, understanding the physical processes that generate them is strongly dependant on coordinated multi-wavelength observations. CMEs display a large diversity in morphology and kinematic properties, but there is presently no statistical evidence that those properties may serve to group them into different classes. When a CME takes place, the coronal magnetic field undergoes restructuring. Much of the current research is focused on understanding how the corona sustains the stresses that allow the magnetic energy to build up and how, later on, this magnetic energy is released during eruptive flares and CMEs. Multi-wavelength observations have confirmed that reconnection plays a key role during the development of CMEs. Frequently, CMEs display a rather simple shape, exhibiting a well known three-part structure (bright leading edge, dark cavity and bright knot). These types of events have led to the proposal of the ‘`standard model’' of the development of a CME, a model which predicts the formation of current sheets. A few recent coronal observations provide some evidence for such sheets. Other more complex events correspond to multiple eruptions taking place on a time scale much shorter than the cadence of coronagraph instruments. They are often associated with large-scale dimming and coronal waves. The exact nature of these waves and the physical link between these different manifestations are not yet elucidated. We also discuss what kind of shocks are produced during a flare or a CME. Several questions remain unanswered. What is the nature of the shocks in the corona (blast-wave or piston-driven?) How they are related to Moreton waves seen in Hα? How they are related to interplanetary shocks? The last section discusses the origin of energetic electrons detected in the corona and in the interplanetary medium. “Complex type III-like events,”which are detected at hectometric wavelengths, high in the corona, and are associated with CMEs, appear to originate from electrons that have been accelerated lower in the corona and not at the bow shock of CMEs. Similarly, impulsive energetic electrons observed in the interplanetary medium are not the exclusive result of electron acceleration at the bow shocks of CMEs; rather they have a coronal origin.     

Modeling coronal streamers and their eruption

Numerical solutions of the time-dependent MHD equations are used to generate ambient coronal streamer structures in a corona characteristic of that near solar minimum. The streamers are then disrupted by slow photospheric shear motion at the base of magnetic field lines within the closed field region, which is currently believed to be responsible for producing at least some CMEs. In contrast to several other simulations of this phenomena, the polytropic index is maintained at a value of 5/3 through the addition of coronal heating. Observations are used as a guide in determining the thermodynamic structure and plasma beta in the ambient corona. For a shear speed of 2.5 km/sec, the streamer configuration evolves slowly for about 65 hours before erupting outward with the formation of a CME. The bright CME leading edge travels outward at a speed of about 240 km/sec, and the sheared field lines follow at a somewhat slower speed. A closed magnetic field region is ejected as the magnetic field lines that were opened by the CME reconnect and reform the streamer.     

S/WAVES: The Radio and Plasma Wave Investigation on the STEREO Mission

This paper introduces and describes the radio and plasma wave investigation on the STEREO Mission: STEREO/WAVES or S/WAVES. The S/WAVES instrument includes a suite of state-of-the-art experiments that provide comprehensive measurements of the three components of the fluctuating electric field from a fraction of a hertz up to 16 MHz, plus a single frequency channel near 30 MHz. The instrument has a direction finding or goniopolarimetry capability to perform 3D localization and tracking of radio emissions associated with streams of energetic electrons and shock waves associated with Coronal Mass Ejections (CMEs). The scientific objectives include: (i) remote observation and measurement of radio waves excited by energetic particles throughout the 3D heliosphere that are associated with the CMEs and with solar flare phenomena, and (ii) in-situ measurement of the properties of CMEs and interplanetary shocks, such as their electron density and temperature and the associated plasma waves near 1 Astronomical Unit (AU). Two companion papers provide details on specific aspects of the S/WAVES instrument, namely the electric antenna system (Bale et al., Space Sci. Rev., 2007) and the direction finding technique (Cecconi et al., Space Sci. Rev., 2007).     

The STEREO Mission: An Introduction

The twin STEREO spacecraft were launched on October 26, 2006, at 00:52 UT from Kennedy Space Center aboard a Delta 7925 launch vehicle. After a series of highly eccentric Earth orbits with apogees beyond the moon, each spacecraft used close flybys of the moon to escape into orbits about the Sun near 1 AU. Once in heliospheric orbit, one spacecraft trails Earth while the other leads. As viewed from the Sun, the two spacecraft separate at approximately 44 to 45 degrees per year. The purposes of the STEREO Mission are to understand the causes and mechanisms of coronal mass ejection (CME) initiation and to follow the propagation of CMEs through the inner heliosphere to Earth. Researchers will use STEREO measurements to study the mechanisms and sites of energetic particle acceleration and to develop three-dimensional (3-D) time-dependent models of the magnetic topology, temperature, density and velocity of the solar wind between the Sun and Earth. To accomplish these goals, each STEREO spacecraft is equipped with an almost identical set of optical, radio and in situ particles and fields instruments provided by U.S. and European investigators. The SECCHI suite of instruments includes two white light coronagraphs, an extreme ultraviolet imager and two heliospheric white light imagers which track CMEs out to 1 AU. The IMPACT suite of instruments measures in situ solar wind electrons, energetic electrons, protons and heavier ions. IMPACT also includes a magnetometer to measure the in situ magnetic field strength and direction. The PLASTIC instrument measures the composition of heavy ions in the ambient plasma as well as protons and alpha particles. The S/WAVES instrument uses radio waves to track the location of CME-driven shocks and the 3-D topology of open field lines along which flow particles produced by solar flares. Each of the four instrument packages produce a small real-time stream of selected data for purposes of predicting space weather events at Earth. NOAA forecasters at the Space Environment Center and others will use these data in their space weather forecasting and their resultant products will be widely used throughout the world. In addition to the four instrument teams, there is substantial participation by modeling and theory oriented teams. All STEREO data are freely available through individual Web sites at the four Principal Investigator institutions as well as at the STEREO Science Center located at NASA Goddard Space Flight Center.     

Heliospheric Physics: Linking the Sun to the Magnetosphere

Research into the heliospheric structure and its relation to the solar boundary is at an impasse. After successful predictions by Parker about the zeroth-order behavior of the heliospheric magnetic field and the solar wind, the heliospheric community struggles to make substantive progress toward a predictive model describing the connections between the Sun and its space environment, between the closed corona and the open corona extending to the planets. This is caused by our lack of understanding of the basic processes heating the corona and transporting open magnetic field. We detail the models used to describe this connectivity, from potential field source surface models to full MHD techniques. We discuss the current limitations of both approaches. Finally, we address a recent attempt to advance our understanding beyond these limitations. At this point in time the proposed theory remains controversial in the community, but it addresses important shortcomings of current approaches outlined above.     

Modelling of Solar Wind, CME Initiation and CME Propagation

Simulations of coronal mass ejections (CMEs) evolving in the interplanetary (IP) space from the Sun up to 1 AU are performed in the framework of ideal magnetohydrodynamics (MHD) by the means of a finite-volume, explicit solver. The aim is to quantify the effect of the background solar wind and of the CME initiation parameters, such as the initial magnetic polarity, on the evolution and on the geo-effectiveness of CMEs. First, three different solar wind models are reconstructed using the same numerical grid and the same numerical scheme. Then, different CME initiation models are considered: Magnetic foot point shearing and magnetic flux emergence. For the fast CME evolution studies, a very simple CME model is considered: A high-density and high-pressure magnetized plasma blob is superposed on a background steady state solar wind model with an initial velocity and launch direction. The simulations show that the initial magnetic polarity substantially affects the IP evolution of the CMEs influencing the propagation velocity, the shape, the trajectory (and thus, the geo-effectiveness).     

Magnetic Fields in the Inner Heliosphere

The structure of Heliospheric Magnetic Field (HMF) is a function of both the coronal conditions from which it originates and dynamic processes which take place in the solar wind. The division between the inner and outer regions of the heliosphere is the result of dynamic processes which form large scale structures with increasing heliocentric distance. The structure of the HMF is normally described in the reference frame based on Parker's geometric model, but is better understood as an extension of potential field models of the corona. The Heliospheric Current Sheet (HCS) separates the two dominant polarities in the heliosphere; its large scale geometry near solar minimum is well understood but its topology near solar maximum remains to be investigated by Ulysses. At solar minimum, Corotating Interaction Regions (CIRs) dominate the near-equatorial heliosphere and extend their influence to mid-latitudes; the polar regions of the heliosphere are dominated by uniform fast solar wind streams and large amplitude, long wavelength, mostly transverse magnetic fluctuations. Coronal Mass Ejections (CMEs) introduce transient variability into the large scale heliospheric structure and may dominate the inner heliosphere near solar maximum at all latitudes.     

Global Non-Potential Magnetic Models of the Solar Corona During the March 2015 Eclipse

Seven different models are applied to the same problem of simulating the Sun’s coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.     

Plasma Diagnostics of the Interstellar Medium with Radio Astronomy

We discuss the degree to which radio propagation measurements diagnose conditions in the ionized gas of the interstellar medium (ISM). The “signal generators” of the radio waves of interest are extragalactic radio sources (quasars and radio galaxies), as well as Galactic sources, primarily pulsars. The polarized synchrotron radiation of the Galactic non-thermal radiation also serves to probe the ISM, including space between the emitting regions and the solar system. Radio propagation measurements provide unique information on turbulence in the ISM as well as the mean plasma properties such as density and magnetic field strength. Radio propagation observations can provide input to the major contemporary questions on the nature of ISM turbulence, such as its dissipation mechanisms and the processes responsible for generating the turbulence on large spatial scales. Measurements of the large scale Galactic magnetic field via Faraday rotation provide unique observational input to theories of the generation of the Galactic field.     

Photospheric and Subphotospheric Dynamics of Emerging Magnetic Flux

Magnetic fields emerging from the Sun’s interior carry information about physical processes of magnetic field generation and transport in the convection zone. Soon after appearance on the solar surface the magnetic flux gets concentrated in sunspot regions and causes numerous active phenomena on the Sun. This paper discusses some properties of the emerging magnetic flux observed on the solar surface and in the interior. A statistical analysis of variations of the tilt angle of bipolar magnetic regions during the emergence shows that the systematic tilt with respect to the equator (the Joy’s law) is most likely established below the surface. However, no evidence of the dependence of the tilt angle on the amount of emerging magnetic flux, predicted by the rising magnetic flux rope theories, is found. Analysis of surface plasma flows in a large emerging active region reveals strong localized upflows and downflows at the initial phase of emergence but finds no evidence for large-scale flows indicating future appearance a large-scale magnetic structure. Local helioseismology provides important tools for mapping perturbations of the wave speed and mass flows below the surface. Initial results from SOHO/MDI and GONG reveal strong diverging flows during the flux emergence, and also localized converging flows around stable sunspots. The wave speed images obtained during the process of formation of a large active region, NOAA 10488, indicate that the magnetic flux gets concentrated in strong field structures just below the surface. Further studies of magnetic flux emergence require systematic helioseismic observations from the ground and space, and realistic MHD simulations of the subsurface dynamics.     

Mars Global Surveyor Observations of Solar Wind Magnetic Field Draping Around Mars

We examine the magnetic field in the martian magnetosheath due to solar wind draping. Mars Global Surveyor provided 3-D vector magnetic field measurements at a large range of altitudes, local times, and solar zenith angles as the spacecraft orbit evolved. We choose orbits with very clean signatures of draping to establish the nominal morphology of the magnetic field lines at local times of near-subsolar and near-terminator. Next, using a compilation of data from Mars Global Surveyor, we determine the average magnetic field morphology in the martian magnetosheath due to the solar wind interaction. The topology of the field is as expected from previous observations and predictions. The magnetic field magnitude peaks at low altitude and noon magnetic local time and decreases away from that point. The magnetic field has an inclination from the local horizontal of 5.6° on average in the dayside magnetosheath and 12.5° on the nightside. The inclination angle is closest to zero at noon magnetic local time and low altitude. It increases both upward and to later local times. The magnetic field in the induced magnetotail flares out from the Mars—Sun direction by 21°. Finally, we compare the observations to gasdynamic model predictions and find that the shocked solar wind flow in the martian magnetosheath can be treated as a gasdynamic flow with the magnetic pileup boundary as the inner boundary to the flow.