Measuring Coronal Magnetic Fields with Remote Sensing Observations of Shock Waves

Bemporad, Alessandro; Susino, Roberto; Frassati, Federica; Fineschi, Silvano

doi:10.3389/fspas.2016.00017

OPINION article

Front. Astron. Space Sci., 27 May 2016

Sec. Stellar and Solar Physics

Volume 3 - 2016 | https://doi.org/10.3389/fspas.2016.00017

This article is part of the Research Topic Coronal Magnetometry View all 15 articles

Measuring Coronal Magnetic Fields with Remote Sensing Observations of Shock Waves

$\r\nAlessandro Bemporad*$ Alessandro Bemporad^*

Roberto Susino

Federica Frassati

Silvano Fineschi

INAF, Turin Astrophysical Observatory, Pino Torinese, Italy

Introduction

State of the Art

Our limited knowledge of the magnetic fields structuring in the solar corona represents today the main hurdle in our understanding of its structure and dynamic. Over the last decades significant efforts have been dedicated to measure these fields, by approaching the problem on many different sides and in particular: (i) by improving our theoretical understanding of the modification (via Zeeman and Hanle effects) induced by these fields on the polarization of coronal emission lines, (ii) by developing new instrumentation to measure directly with spectro-polarimeters these modifications, (iii) by improving the reliability of the extrapolated coronal fields starting from photospheric measurements, (iv) by developing new techniques to analyse existing remote sensing data and infer properties of these fields, or by combining all these different approaches (e.g., Chifu et al., 2015).

In this paper we focus on the fourth method, discussing the advantages and disadvantages of techniques recently developed. Over the last few years it has been shown that coronagraphic white light (WL) observations of major Coronal Mass Ejections (CMEs) show the presence of hemispherical regions expanding ahead of the CME fronts (Figure 1) where faint increases in the WL intensity are detected (e.g., Vourlidas et al., 2003; Ontiveros and Vourlidas, 2009). More recently similar features have also been detected in the early stages of CME developments observed with EUV imagers (e.g., Kozarev et al., 2011). Usually these slightly brighter regions are interpreted as the coronal plasma compressed by the transit of the shock wave driven by the super-alfvénic expansion of the CME, a region called “shock sheath.” This interpretation has also been supported by MHD simulations (e.g., Manchester et al., 2008) and forward modeling of synthetic observations (Vourlidas et al., 2013) showing the formation of an arch shaped feature very similar to what is observed in WL coronagraphs. In this scenario, it has been recently demonstrated that measurements of coronal fields can be successfully derived from remote sensing observations of interplanetary shocks with at least two techniques, briefly described in the next Sections.

FIGURE 1

Figure 1. The expansion of a CME-driven shock wave as observed in the LASCO/C2 WL coronagraph (Bemporad and Mancuso, 2010).

The Shock Standoff Distance Technique

Gopalswamy and Yashiro (2011) first applied to remote sensing observations of CMEs this technique based on the measurement of the so called shock “standoff distance” Δ, i.e., the distance between the shock driver (i.e., the expanding flux-rope) and the shock wave. Given Δ and the radius of curvature R of the flux-rope, a semi-empirical formula relates the Δ/R ratio with the value of the shock sonic Mach number M_s. The M_s value derived from the observed Δ/R ratio can be combined with the measured shock speed v_sh to provide a measurement of the upstream Alfvén speed v_A, by assuming that the Alfvéninc Mach number M_A ~ M_s, and by assuming a value for the pre-shock solar wind speed v_out. Given also the measured pre-shock coronal electron density n_e, the magnetic field strength is provided by B = v_A (μ m_p n_e)^1∕2. Various authors applied this method to shocks observed with WL coronagraphs (Kim et al., 2012), EUV full disk imagers (Gopalswamy et al., 2012) and heliospheric imagers (Poomvises et al., 2012), providing measurements of the coronal field up to heliocentric distances larger than 120 solar radii.

Nevertheless, this technique has various limitations. As pointed out by Savani et al. (2012) there are some uncertainties in the expression of the semi-empirical relationship between the Δ/R ratio with M_s which has two expressions:

\begin{array}{l} \frac{Δ}{R} = k \frac{(γ - 1) M_{s}^{2} + 2}{(γ + 1) M_{s}^{2}} and \frac{Δ}{R} = k \frac{(γ - 1) M_{s}^{2} + 2}{(γ + 1) (M_{s}^{2} - 1)}, \end{array}

where, the first expression was originally derived by Priest (1984) for the hydrodynamic (HD) shock case, while the second one is an adjustment made on intuitive basis by Farris and Russell (1994) to avoid the fact that the first expression is not valid under low Mach number regime (M_s < 3). A second problem is that (as discussed in Savani et al., 2012) the value of the constant k can vary depending on the oblateness of the shock driver (hence of the CME flux rope), which is usually a poorly constrained parameter (because of projection effects and unknown 3D shape of the CME) and is also expected to change as function of time/distance (Savani et al., 2011). The third problem is that (even though more appropriate MHD formulas were provided by Priest and Forbes, 2007) the above expression holds for the HD case, and can be applied in the magneto-hydrodynamic (MHD) case if one assumes that M_A ~ M_s, hence v_A ~ v_s, while usually in the solar corona v_A >> v_s. The HD solution is a good approximation for the MHD case when the magnetic pressure is much smaller than the plasma pressure (β >> 1), but this is not the case for the lower corona and is true only for the interplanetary medium, hence for the analysis of heliospheric images or in situ measurements (Fairfield et al., 2001). Moreover, the application of the above method to CMEs observed with WL coronagraphs requires the assumption of another free parameter, the solar wind speed, to convert the measured shock speed into an upstream speed. Also, this technique is able to provide only a 1D radial profile of the magnetic field strength in the direction of propagation of the shock.

The WL-UV Technique

To overcome these limitations, an alternative technique to derive coronal magnetic fields from shocks has been developed by the Group in Turin. In particular, Bemporad and Mancuso (2010) first demonstrated that WL and UV observations of a coronal shock (Figure 1) can be combined to infer the magnetic field strength and deflection across the shock front. In the technique pre- and post-shock densities are derived as usual from WL, while pre-shock solar wind velocities and plasma temperatures are derived from UV with standard techniques. The WL images are also used to estimate the shock speed and the inclination angle ϑ_Bn between the normal to the shock surface (projected on the plane of the sky) and the pre-shock magnetic field, assuming that above two solar radii it's radial. Given the shock compression ratio, shock inclination and velocity derived from WL, and the pre-shock plasma temperature and outflow velocity derived from UV, the MHD Rankine-Hugoniot equations written for the general case of an oblique shock provide not only the post-shock plasma temperatures, but also the post-shock outflow velocities, pre- and post-shock magnetic field strengths, as well as the velocity and magnetic field vector deflections across the shock. This technique, being a combination of WL and UV data, is very promising for application on future observations by the Metis coronagraph (Antonucci et al., 2012; Fineschi et al., 2012) on-board the forthcoming ESA-Solar Orbiter mission.

In a following work (Bemporad and Mancuso, 2013) we introduced an empirical relationship to derive the Alfvénic Mach number M_A∠ for the general case of an oblique shock directly from WL observations of shocks as

\begin{array}{l} M_{A ∠} = \sqrt{{(M_{A ⊥} sin ϑ_{B n})}^{2} + {(M_{A ∕ ∕} cos ϑ_{B n})}^{2}} \end{array}

where M_A_⊥ and M_A∕∕ are the Alfvénic Mach numbers for the special cases of perpendicular and parallel shocks, respectively, that are easily estimated from the shock compression ratio by assuming low plasma β condition. This empirical relationship, tested observationally (Bemporad et al., 2014) and numerically (Bacchini et al., 2015), allowed the authors to derive for the first time the magnetic field strength in the corona crossed by a shock over a huge region covering 10 solar radii in altitudes and 110° in latitudes (Susino et al., 2015).

This technique has also some uncertainties. First, it assumes that the pre-shock magnetic field is radial, an hypothesis which is more reliable above coronal streamer regions and usually not acceptable above polar regions, where super-radial expansion of the solar wind occurs. Second, it assumes that all the plasma parameters are derived on the plane of the sky, and (as it happens usually for remote sensing observations) it is hard to estimate the effect of the line of sight (LOS) integration in the determination of the shock compression ratio and all the other parameters derived from this quantity. Issues related with LOS integration are present in all kinds of remote sensing observations of optically thin plasmas. Nevertheless, we point out that in the above determinations of the shock compression ratios from WL the authors showed how to take into account different LOS extensions of the shocked plasmas at different altitudes (see e.g., Bemporad and Mancuso, 2010), thus improving the simplistic usual assumption by previous authors of a constant LOS thickness (e.g., Ontiveros and Vourlidas, 2009).

Present Debates on Magnetic Field Measurements with Shock

After a first excitement by part of the community for the possibility to directly observe the propagation of shocks in the corona, the validity of these results are becoming more debated. It has been pointed out that the observed hemispherical WL features expanding ahead of CMEs could not be unambiguously differentiated from the front of compression waves, due to the pileup of coronal plasma lying above the expanding CME. The hypothesis that these features are really shocks is better supported when at the same time at least one of these observational features is also reported:

1. remote sensing observation in radio data of a type-II burst;

2. remote sensing observation on EUV-UV spectroscopic data of non-thermal line broadening ahead of the CME front;

3. in situ observation of a gradual SEP event.

Another observational property supporting the real formation of a CME-driven shock is provided directly by the speed of the CME front as observed in WL images. In fact, the formation of a shock is expected only when the driver (i.e., the expanding CME) is moving in the corona faster than the local Alfvén, sound and/or magneto-sonic speeds. In agreement with this, faster CMEs are also statistically more associated with type-II radio bursts (Gopalswamy et al., 2010) and a good correlation exists between CME speeds and SEP fluxes (Kahler and Vourlidas, 2013). In any case, the detection of one or more of these features is not telling us where the shock is formed in the corona, and the possibility that the observed WL feature is not the real shock cannot be completely ruled out.

A criticism often also pointed out is that, if these hemispherical features are really the shocks, then the intensity variation plotted perpendicular to the shock should contain a steep discontinuity, similar to those observed in density measurements acquired by in situ spacecraft at shock transit. On the other hand, radial plots of the WL intensity across these hemispherical features can show a gradual increase of the WL, hence a gradual increase in the electron column density of the plasma. Nevertheless, this is only an apparent inconsistency due to a combination of the 3D geometry of the shock surface and the integration along the LOS through the optically thin coronal plasma. In fact, the shock surface will have an almost hemispherical shape (as demonstrated by 3D numerical simulations), and the location of the shock boundary in the 2D projected view will correspond to the pixels where the lines of sight graze this surface. Then, moving radially into the shock sheath, the fraction of the LOS intercepting the sheath region will progressively increase, leading to the observed progressive increase in the WL intensity, as nicely shown by Manchester et al. (2008).

More in general, the above techniques require many assumptions to derive the magnetic fields from the observed WL images, making it difficult to estimate the real accuracy of these measurements. We suggest that these uncertainties will be better constrained by combining data analysis with forward modeling of synthetic observations (e.g., Gibson, 2015).

Conclusions and Future Perspectives

Despite the uncertainties discussed here, both the “standoff” and “WL-UV” techniques (recently applied for inter-comparison to the same event by Susino et al., 2015) are now providing reasonable measurements of the coronal magnetic fields over broad intervals in altitudes and latitudes never reached before. Comparison with other observational methods and the analysis of synthetic data will likely help us to improve and optimize these techniques in the future. A limit of these techniques is that field measurements are provided only after the eruption responsible for the shock wave: hence these methods give an “a posteriori” knowledge of the pre-CME coronal fields, and likely will be applicable for forecasting purposes only when statistical analyses will be carried out.

Author Contributions

All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

RS acknowledges support from ASI contract I/013/12/0-1; FF acknowledge support from INAF Ph.D. Grant.

References

Antonucci, E., Fineschi, S., Naletto, G., Romoli, M., Spadaro, D., Nicolini, G., et al. (2012). Multi element telescope for imaging and spectroscopy (metis) coronagraph for the solar orbiter mission, in Proceedings of SPIE, Vol. 8443 (Amsterdam).

Google Scholar

Bacchini, F., Susino, R., Bemporad, A., and Lapenta, G. (2015). Plasma physical parameters along CME-driven Shocks. II. Observation-simulation comparison. Astrophys. J. 809, 58–70. doi: 10.1088/0004-637X/809/1/58