This site uses cookies. By continuing to use this site you agree to our use of cookies. To find out more, see our Privacy and Cookies policy.

Prediction of the Sun's Coronal Magnetic Field and Forward-modeled Polarization Characteristics for the 2019 July 2 Total Solar Eclipse

, , , , and

Published 2020 February 11 © 2020. The American Astronomical Society. All rights reserved.
, , Citation Soumyaranjan Dash et al 2020 ApJ 890 37 DOI 10.3847/1538-4357/ab6a91

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/890/1/37

Abstract

On 2019 July 2 a total solar eclipse—visible across parts of the Southern Pacific Ocean, Chile, and Argentina—enabled observations of the Sun's corona. The structure and emission characteristics of the corona are determined by underlying magnetic fields, which also govern coronal heating and solar eruptive events. However, coronal magnetic field measurements remain an outstanding challenge. Coronal magnetic field models serve an important purpose in this context. Earlier work has demonstrated that the large-scale coronal structure is governed by surface flux evolution and memory buildup, which allows for its prediction on solar rotational timescales. Utilizing this idea and based upon a 51 day forward run of a predictive solar surface flux transport model and a potential field source surface model, we predict the coronal structure of the 2019 July 2 solar eclipse. We also forward model the polarization characteristics of the coronal emission. Our prediction of two large-scale streamer structures and their locations on the east and west limbs of the Sun match eclipse observations reasonably well. We demonstrate that the Sun's polar fields strongly influence the modeled corona, concluding that accurate polar field observations are critical. This study is relevant for coronal magnetometry initiatives envisaged with the Daniel K. Inouye Solar Telescope, Coronal Multichannel Polarimeter and upcoming space-based instruments such as Solar Orbiter, Solar Ultraviolet Imaging Telescope and the Variable Emission Line Coronagraph on board the Indian Space Research Organisation's Aditya-L1 space mission.

Export citation and abstract BibTeX RIS

1. Introduction

The Sun's surface and its outer atmosphere—the corona—forms the time-dependent inner boundary of the heliosphere, and therefore, forces it. Magnetic field dynamics in the corona including magnetic reconnection and heating is induced by photospheric flux emergence and evolution. In turn this coronal dynamics spawn the solar wind, solar flares, coronal mass ejections (CMEs) and energetic particle flux and high energy radiation which collectively create hazardous space environmental conditions. Understanding the origin of space weather at the Sun and building predictive capabilities is a high priority goal in the space sciences (Schrijver et al. 2015). Although the origin of space weather can be traced back to the dynamics of magnetic fields in the Sun's corona, direct observations of coronal magnetic fields is challenging because of the low photon flux associated with the tenuous coronal plasma. Thus computational approaches based on our theoretical understanding of how the Sun's magnetic fields emerge through the surface, evolve, and permeate the corona are essential to model and understand coronal dynamics.

The Alfvén speed in the corona (the speed at which magnetic disturbances propagate) is much faster (of the order of 1000 km s−1) than the large-scale surface motions that drive coronal evolution (1–3 km s−1). Consequently, a reasonable approach in theoretical models is to assume that the coronal magnetic field distribution evolves quasi-statically in response to photospheric forcing (Cook et al. 2009). These models can be categorized into four broad classes: potential field source surface extrapolation models, force-free models, magnetohydrostatic models, and full magnetohydrodynamic (MHD) models. While the first three primarily provide the three-dimensional close-to-equilibrium magnetic field structure in the solar corona, the last approach can self-consistently provide the magnetic field as well as thermal characteristics. All these modeling approaches have their own advantages and limitations and often a particular model is chosen based on the balance of convenience and sophistication necessary for a particular problem. For a detailed account on coronal magnetic field models see Mackay & Yeates (2012) and Yeates et al. (2015).

Testing these theoretical modeling approaches through "true" predictions and constraining them through observations is best achieved at the time of solar eclipses when lunar occultation masks the bright solar disk revealing the Sun's faint coronal structure. In an earlier work, Nandy et al. (2018) postulated a novel methodology for long-term coronal field prediction. Posteclipse assessment demonstrated that they successfully predicted all the major large-scale structures of the great American solar eclipse (2017 August 21) including a pseudostreamer that was not captured in an alternative prediction methodology (Mikić et al. 2018).

Nandy et al. (2018) used a data-driven Predictive Solar Surface Flux Transport (PSSFT) model developed at the Center of Excellence in Space Sciences India (CESSI; Bhowmik & Nandy 2018) to first predict the surface magnetic field distribution and then feed this into a potential field extrapolation model to obtain the prediction of the coronal magnetic field structure. The PSSFT model simulates the evolution of the surface magnetic field driven by the emergence of tilted active regions and redistribution of associated magnetic flux mediated via large-scale plasma flows and supergranular diffusion. This process is known as the Babcock–Leighton (B-L) mechanism (Babcock 1961; Leighton 1969; Wang et al. 1989; van Ballegooijen et al. 1998; Schrijver 2001; Sheeley 2005; Cameron et al. 2010; Mackay & Yeates 2012) in the context of dynamo theory and bridges solar internal dynamics with coronal field evolution.

The CESSI PSSFT model used in Nandy et al. (2018) is driven by observed sunspot data over century timescale and calibrated with polar flux observations spanning multiple solar cycles (Muñoz-Jaramillo et al. 2012; Bhowmik & Nandy 2018). The slow evolution of the photospheric magnetic field aids in building up a long-term memory (due to mostly deterministic surface flux transport processes) in the PSSFT model. This enables the model to make long-term predictions of the large-scale surface field distribution including the high latitude polar fields (Bhowmik & Nandy 2018)—which can be used as inputs for making coronal field predictions. An accurately predicted surface magnetic field map ensures evaluation of large-scale structures of the coronal magnetic field (Schrijver & De Rosa 2003) with better precision. We note that given the large errors in polar field measurements due to projection effects from plane-of-ecliptic observations, such predictive surface field evolution models are critical to accurate coronal field predictions.

Under low plasma-β (i.e., magnetic pressure dominating over gas pressure) and current-free approximations the potential field source surface (PFSS) extrapolation technique (Altschuler & Newkirk 1969; Schatten et al. 1969) may be employed to simulate coronal magnetic fields. We note that dynamic coronal simulations based on magnetofrictional approaches transmit this surface memory to the corona (Yang et al. 1986; van Ballegooijen et al. 2000; Yeates 2014), while PFSS-like extrapolation relies on the memory of the surface field itself. Dynamic MHD models or magnetofrictional models are certainly more realistic approaches to simulating the coronal field, especially during epochs when a large number of active region associated highly non-potential structures are present. However, it is absolutely critical that the major driver of coronal models, i.e., the Sun's surface magnetic field map is accurate. Our methodology lays more emphasis on this latter aspect. While PFSS model-based extrapolations assume a current-free corona and therefore has some limitations, it does a reasonable job at extracting the large-scale coronal field. Additionally, forward-modeled polarization characteristics inferred with such a model—which are computationally very efficient—has the potential to provide near real time validation of coronal field observations utilizing magnetometry technique.

Here, we utilize the Nandy et al. (2018) prediction methodology—a combination of the PSSFT and PFSS models—to predict the large-scale coronal magnetic structure of the 2019 July 2 total solar eclipse. A short research note with the predicted coronal magnetic structure is communicated in Dash et al. (2019). This work contains description of the methodology with detailed in-depth analysis of the coronal magnetic field structure and topology including the possibility of a pseudostreamer appearing in the future. Additionally, we present forward-modeled polarization maps for the coronal emission expected during the eclipse which may be compared with coronal magnetometry observations. It is thought that compromised polar field observations lead to a discrepancy in the orientation of the streamer tips which is handled by adding an apparent missing flux near the poles in some predictive approaches (Riley et al. 2019). To test this idea we further perform a parameter space study to explore how varying flux content in the Sun's polar regions impact the orientation of coronal streamers. Finally, we compare our predictions with the white-light corona, extended coronal streamers and plumes observed during the eclipse. The rest of the paper is structured as follows: a brief description of the computational models is provided in Section 2; results are presented in Section 3; a comparison between the simulation-generated and observed corona on the day of eclipse is given in Section 4; concluding discussions follow in Section 5.

2. Numerical Models and Sunspot Data Input

Evolution of magnetic field on the solar surface is governed by the magnetic induction equation. As photospheric magnetic field is predominantly in the radial direction (Solanki 1993), in our PSSFT model we solve only the radial component (Br) of the induction equation which in spherical polar coordinates is given by

Equation (1)

The symbols, θ and ϕ represent colatitude and longitude, R is the solar radius. The large-scale velocity fields, the differential rotation and meridional circulation on the solar surface are denoted by ω(θ) and v(θ), respectively. These plasma flows are modeled using empirical functions (Bhowmik & Nandy 2018), which are observationally verified (Snodgrass 1983; van Ballegooijen et al. 1998). The parameter ηh is the effective diffusion coefficient associated with the turbulent motion of supergranules and S(θ, ϕ, t) is the source term describing the emergence of new sunspots. We note that any new spot is included in the PSSFT model when it has the maximum area coverage on the solar surface during its lifetime (or period of visibility on the solar disk).

The century-scale PSSFT simulation utilizes sunspot data recorded by Royal Greenwich Observatory (RGO) and United States Air Force (USAF)/National Oceanic and Atmospheric Administration (NOAA) for the period 1913 August until 2016 September (Bhowmik & Nandy 2018). The database provides essential details associated with active regions such as position on the solar surface, area coverage, and the corresponding time, etc. From 2016 October onward, the observed sunspot data used in the PSSFT model is acquired from the Helioseismic and Magnetic Imager (HMI) on board NASA's Solar Dynamic Observatory (SDO). The sunspot area given by HMI is scaled down by a constant factor for consistency with the RGO-NOAA/USAF recorded area; this cross-instrument calibration is necessary as the original calibration of the century-scale PSSFT simulation is based on the later database. The requirement of such cross-instrumental calibration arises due to the different techniques employed for recording sunspot data (Muñoz-Jaramillo et al. 2015). The scaling factor is evaluated by comparing the area associated with 1174 active regions recorded in both the NOAA/USAF and HMI database for an overlapping period of six years (2010 May–2016 September). We measure the ratio between the area recorded by the two databases and use a Gaussian fit to the distribution for calculating the mean. The mean ratio (0.244) is utilized to scale down the HMI-recorded area. The magnetic flux of the active region is evaluated based on a linear empirical relation (Sheeley 1966; Dikpati et al. 2006): Φ(A) = 7.0 × 1019A maxwells, where A is the area in unit of microhemispheres. We assume all active regions appearing on the solar surface are ideal β-spots where the magnetic flux is equally distributed between the leading and following polarity spots. The associated tilt angle is decided by Joy's law with a cycle-dependent refinement (see Equation (5) of Bhowmik & Nandy 2018).

The last observed active region included in the PSSFT model is AR12741—which attained its maximum area coverage on 2019 May 12. The PSSFT model is then forward run until the day of the eclipse on 2019 July 2 (assuming no new sunspot emergence until then) to generate the predicted surface magnetic field distribution.

To model the corresponding coronal magnetic field structure we use an PFSS extrapolation technique with the predicted surface magnetic field used as the bottom boundary of the computational domain. The extrapolation is extended up to the source surface (2.5 R) beyond which we assume the magnetic field to become radial (Davis 1965). We utilize the PFSS extrapolation model developed by Yeates (2018).4

3. Results

3.1. Prediction of the Coronal Magnetic Field Structure

The PSSFT predicted surface magnetic field distribution corresponding to 2019 July 2 is depicted in Figure 1(a). The radial component of the magnetic field (Br) obtained from forward running the PSSFT simulation is plotted as a function of latitude and longitude. The region between the east and west limbs represents the surface magnetic field map on the solar disk during the eclipse. As the Sun has entered the minimum phase of solar cycle 24, we witness strong concentrations of unipolar magnetic flux of opposite polarities at both polar caps. This is indicative of a large-scale magnetic configuration with dominant dipolar characteristics. Note that the locally confined clusters of magnetic field near the equatorial region (in this and the following figures) correspond to the residual flux of the emerged active regions and should not be mistaken as new spot emergence. Figure 1(b) represents the polarity distribution of the surface magnetic field which is evaluated by calculating the quantity, $P(\theta ,\phi )={B}_{r}(\theta ,\phi )/| {B}_{r}(\theta ,\phi )| $. This image represents how the magnetic polarity is distributed in the large-scale structures on the surface.

Figure 1.

Figure 1. (a) Predicted solar surface magnetic field (Br) for 2019 July 2 obtained from the PSSFT model. The colorbar represents the field strength saturated to ±10 G. (b) Saturated image to indicate the polarity (and inversion lines) of the surface field distribution. The solid black line represents the central meridian, the dashed blue and maroon lines correspond to the east and west limb of the visible solar disk on 2019 July 2. The rectangular box on the east limb represents an area of interest which has the possibility of hosting a pseudostreamer.

Standard image High-resolution image

The predicted surface magnetic field is then utilized in the PFSS extrapolation model as its lower boundary condition to obtain the global coronal magnetic field. The surface field is mapped on the circular disk of Figure 2(a), where the white curves correspond to the magnetic polarity inversion lines. The PFSS model generated coronal field lines near the limb are extracted to represent the plane-of-sky corona expected on the day of the eclipse. The locations marked as Regions 1 and 3 correspond to the north and the south pole. Based on the polarity of the magnetic field at the foot-points the open field lines are color-coded in light red (radially outward) and cyan (radially inward) while all closed lines are colored in black. The predicted coronal field has two prominent streamer structures, one on the west limb (whose tip or cusp is denoted as Region 2) and the other on the east limb (Region 4). The set of closed black curves in Figure 2(a) corresponding to large-scale closed loop structures separate open flux (coronal holes) of opposite polarities—the defining characteristic of helmet streamers (Wang et al. 2007). Although the foot-points of the closed magnetic field lines associated with streamers are extended across the solar equator and major portion of the two hemispheres (north and south) they still have a directional sense. The cusp of the streamer on the east limb (Region 4) is centered below the solar equatorial plane. The cusp of the streamer on the west limb (Region 2) is centered above the equatorial plane. These kind of streamer structures are known to form extended coronal plumes visible in large-angle coronagraphs such as the Large Angle and Spectrometric COronagraph (LASCO) instrument on board the Solar and Heliospheric Observatory spacecraft (SOHO). The presence of closed field line structures on the east limb (somewhat inclined to the plane of sky) on the northern edge of the streamer (Region 5) makes the appearance of this streamer somewhat diffused with the northern edge extended to higher latitudes, with the likelihood of a corresponding coronal plume overlying Region 5.

Figure 2.

Figure 2. (a) Predicted large-scale coronal structure of the 2019 July 2 solar eclipse with open field lines denoted in light red (radially outward) and cyan (radially inward). Closed field lines are denoted in black. The white curves on the solar disk correspond to polarity inversion lines. Regions marked 1 and 3 denote open field lines (coronal holes) at the polar regions. Regions marked 2 and 4 are the cusps of streamers which are predicted to occur on the west (right) and east (left) limbs, respectively, indicating the inclination at which extended coronal plumes may be expected in large-angle coronagraph observations. Region 5 corresponds to the additional closed field lines adjacent to the helmet streamer on the east limb. Region 6 overlies a surface field distribution which is symptomatic of a pseudostreamer host but is unlikely to generate one by 2019 July 2 (see the text). (b) A representative white-light corona synthetically constructed from the magnetic field structure to indicate the plausible white light appearance of the corona during the total solar eclipse.

Standard image High-resolution image

The large-scale plasma motion or supergranular convection on the photosphere can generate foot-point motion resulting in the rise of the closed loops within the helmet streamers. This can trigger a three-dimensional reconnection with the overlying open field lines. Consequently, plasma materials get energized resulting in heating and emission. In addition, Thomson scattering of photons from regions with enhanced charged particle density (such as in magnetic loops) contribute to the overall appearance of the white light corona (structured by magnetic fields). These are processes we cannot capture with our simplistic coronal field model, however, our coronal magnetic field may be utilized to generate a "synthetic" white light corona. Figure 2(b) is a representation of the white-light corona based on the simulated magnetic field distribution. Here the magnetic field lines are plotted using a single color (white) wherein closed field lines are assigned more weight as compared to the open ones. Additionally, an inverse r2 filter is used on the resulting image to generate the "synthetic" corona. This image may be taken as a qualitative guide to the white light corona (within 2.5 R) observed on 2019 July 2.

3.2. Coronal Magnetometry: Forward-modeled Coronal Polarization Characteristics

The Stokes I, Q, U, V polarization vectors have shown considerable potential to be used as a direct diagnostic of coronal magnetism (Lin et al. 2004). The line-of-sight magnetic field strength provided by the longitudinal Zeeman effect is manifested in circular polarization, V. The linear polarization represented by Stokes vectors Q and U originate due to resonance scattering of photons by the electrons in the corona. These Stokes vectors Q and U contain information about the direction of the magnetic field projected onto the plane of sky of the corona. To generate the synthetic Stokes data (Judge et al. 2006) we utilize the FORWARD tool set (Gibson et al. 2016) using a spherically symmetric hydrostatic temperature and density profile weighted differently to open and closed field lines for the background corona (Gibson et al. 1999) which is incorporated using the TOPOLOGY module. It is important to note that the normalized polarization vectors are independent of the density of the medium. FORWARD uses the Coronal Line Emission polarimetry code developed by Judge & Casini (2001) to synthesize Stokes (I, Q, U, V) line profiles. When the magnetic field is oriented at the Van Vleck angle = 54fdg74 relative to the radial direction (in the solar coordinate system) the linear polarization ($L=\sqrt{{Q}^{2}+{U}^{2}}$) becomes zero (Van Vleck 1925). We refer to the regions where this occur as "Van Vleck nulls." The Stokes V/I profile gives the line-of-sight intensity-weighted average magnetic field strength (BLOS).

In Figure 3 we present the forward-modeled polarimetric maps derived from the predicted coronal magnetic field structure for 2019 July 2. Figure 3(a) shows the linear polarization vectors (blue lines) corresponding to emission from the Fe xiii transition at 1074.7 nm. The magnetic field vectors in the plane of sky are denoted as red arrows. As expected, the linear polarization vectors correctly identify the direction of the plane-of-sky magnetic field, including the diverging and converging radial fields from the north and south solar poles and the curvature underlying the closed streamer belts. The degree of linear polarization L/I projected onto the plane of sky is shown in Figure 3(b), wherein, the dark regions corresponding to Van Vleck nulls denote curved field lines of closed streamer loops (oriented at the Van Vleck angle to the local radial direction); note the correspondence with Figure 2(a). Intriguingly, the low-lying double-loop structure evident at high latitudes in the northeast limb hints at the existence of a pseudostreamer (Rachmeler et al. 2014), to which we come back later. Figure 3(c) depicts the circular polarization given by Stokes V/I which is proportional to the line-of-sight magnetic field strength. Note that Stokes $V\propto {B}_{\mathrm{LOS}}\cos \theta $, where θ is the polar angle of the magnetic field relative to the line of sight (Casini & Judge 1999). Therefore, its sign indicates the direction of the field toward or away from the line of sight. Blue shading indicates line-of-sight integrated magnetic fields directed away from the observer (while red denotes field directed toward the observer). The plot indicates the presence of clock-wise closed loops (as viewed from solar north) connecting the positive polarity patch near the east limb to a negative polarity patch behind the limb which we confirm from the full 3D coronal magnetic field structure (not shown here).

Figure 3.

Figure 3. Forward-modeled coronal polarization maps. (a) Plane-of-sky magnetic field vectors (red arrows) and linear polarization vectors (blue lines) with polarization brightness in the background (b) Degree of linear polarization, L/I. (c) The circular polarization given by Stokes V/I.

Standard image High-resolution image

Such predicted polarization characteristics can, on the one hand, aid in the interpretation of coronal magnetometry observations made during solar eclipses or with coronagraph instruments. On the other hand, the observations themselves can help constrain coronal field magnetic field models which underlie the forward-modeled polarization characteristics.

3.3. Evaluating the Possibility of a Pseudostreamer

The coronal magnetic field generated from the PSSFT−PFSS coupled model has a narrow collimated structure (marked as Region 6 on the northeast limb of Figure 2(a)) with a localized void and very low-lying field loops closing near the surface. Such a magnetic configuration is indicative of a pseudostreamer which has not quite matured to be visible. Visible pseudostreamers, in general, materialize under certain conditions in the coronal regions which overlie surface magnetic field distribution where a narrow region of one polarity separates two surrounding opposite polarity open flux patches (Wang et al. 2007; Rachmeler et al. 2014; Abbo et al. 2015). The surface polarity configuration of this region of interest is highlighted within the rectangular box of the surface magnetic field in Figure 1(b). We observe a narrow region of negative polarity separating two positive polarity patches on the solar surface and two associated polarity inversion lines. This region is centered around 42 N on the east limb. Could this region with pseudostreamer favorable surface magnetic field distribution mature into a visible pseudostreamer?

A basic pseudostreamer configuration is characterized by two polarity inversion lines under the cusp of the streamer and a pair of loop arcades (consisting of closed field regions) within the streamer. When conditions are favorable (e.g., induced foot-point motions due to shearing or differential rotation) the closed loops can rise and undergo two types of interchange reconnection with the surrounding open field lines. One possibility is that closed field lines at the outer edge of any arcade can undergo three-dimensional reconnection with adjacent open field lines. Another possibility is that sheared closed field lines can rise and expand into the corona where they reconnect with opposite polarity open field lines across the null point (near the source surface). The former necessitates loop arcades with field strength similar to the adjacent open field and the latter necessitates loop arcades which are high enough to form an x-point with the adjacent open field lines.

We investigate the coronal magnetic field configuration and the underlying surface magnetic field distribution in this region of interest in Figure 4. We find that the weak field closed loops (black) are extremely low-lying and have not yet matured into arcade-like structures (see Figure 4(a)). Given the weak underlying fields at this time they are also unlikely to reach higher altitudes to form an x-point with the open field lines (light red). The possibility of interchange reconnection that could trigger a visibly stable pseudostreamer by 2019 July 2, is therefore, low. Intriguingly, however, a longer forward run of the PSSFT model for a further solar rotational timescale indicates that transport of positive radial flux toward the north pole is still ongoing which is expected to increase the field strength of the open field patch just south of the region of interest; note the increasing height of the local peak indicated by the arrow in Figure 4(b). This flux pile-up could give rise to a visible pseudostreamer under favorable conditions at a date later than the day of the eclipse.

Figure 4.

Figure 4. (a) The predicted coronal magnetic field line distribution around "Region 6" (with reference to Figure 2). The weak field closed loops denoted in black are extremely low-lying and enclosed by open field lines denoted in light red. (b) The black curve represents the surface magnetic field strength along the east (left) limb on the day of the eclipse. The magenta curve denotes the PSSFT model predicted field strength along the same limb one solar rotation later indicating the possibility of increased positive radial field transport to the edge of the region of interest.

Standard image High-resolution image

3.4. Dependence of Streamer Tilt Angles on Surface Magnetic Field Distribution

It is evident that in the SFT–PFSS framework the large-scale coronal structures such as helmet streamers (HS) are primarily influenced by the surface magnetic field distribution. We explore this interconnection further to explore how variations in the polar flux impact the tilt angle of the helmet streamers. We consider two additional cases with the polar cap magnetic flux being (artificially) halved and doubled (within 30 polar caps) relative to the primary PSSFT simulation used for our prediction. We calculate the tilt angle of the HS by measuring the slope of the line joining the tip the structure and the center of the solar disk. Table 1 presents the dependence of the tilt angle of the streamers on the the magnetic dipole moment for these three cases. This dependence is plotted in Figure 5 which indicates that with increasing polar field strength the tips of the streamers are more likely to bend toward the equatorial plane of the Sun. It is important to note that the buildup of the surface magnetic field is primarily driven by the B-L mechanism (Bhowmik & Nandy 2018). Observed large-scale magnetic flux patches are inserted as ideal β-spots in our simulation to drive the flux transport process. Using more realistic observed active region configurations as well as newly emerged small scale (non-active region) flux patches may further improve the accuracy of the simulated flux transport dynamics.

Figure 5.

Figure 5. Plot showing the variation of the tilt angle of the streamer with the dipole moment. Filled cross and star denotes the properties of streamers on the northwest limb and southeast limb, respectively. The marker groups colored orange, green, and red correspond to half polar field strength, original polar field strength and double polar field strength of the SFT simulations, respectively.

Standard image High-resolution image

Table 1.  Variation of Simulated Streamer Tilt Angle with the Sun's Magnetic Dipole Moment

  HS Tilt Angle HS Tilt Angle Dipole Moment
  at Southeast Limb at Northwest Limb (Mx)
SFT polar field halved −23fdg68 14fdg19 3.4418e+22
SFT original −24fdg77 11fdg47 5.1204e+22
SFT polar field doubled −4fdg63 9fdg72 8.4901e+22

Download table as:  ASCIITypeset image

4. Comparison of the Model Prediction with Eclipse Observations

A comparative study between the observed and predicted corona completes our analyses. We utilize coronal images provided by K-coronagraph (K-Cor) from COronal Solar Magnetism Observatory (COSMO) facility suite at Mauna Loa Solar Observatory (MLSO) and the LASCO C2 instrument for this purpose. During a total solar eclipse, large-scale magnetic structures in the lower corona which are otherwise not visible, becomes observable. In Figure 6 we present a composite image which depicts the PSSFT model predicted coronal magnetic field lines overlaid on COSMO K-Cor white light observation for 2019 July 2 until 3R. A LASCO C2 image is superimposed on this image until 6R to track the helmet streamers and polar plumes in the extended corona. The bright regions in the MLSO and LASCO observations indicate the presence of streamer like structures. The closed magnetic loops, as predicted by PSSFT, envelope the observed helmet streamers (marked as Region 2 and 4 in Figure 2). Extended coronal plumes observed by LASCO C2 are observed to emanate from the cusp regions of these streamers.

Figure 6.

Figure 6. Our predicted large-scale coronal magnetic field overlaid on the white-light corona observed by the MLSO/K-Cor instrument (up to 3.0 R). A LASCO C2 coronagraph image is plotted until 6.0 R for the 2019 July 2 eclipse. The PFSS generated coronal magnetic field lines plotted in red and black correspond to open and closed magnetic fields, respectively.

Standard image High-resolution image

5. Concluding Discussions

In summary, here we predict the Sun's coronal magnetic field structure expected to be observed during the total solar eclipse of 2019 July 2. We also present forward-modeled polarization characteristics of the coronal magnetic field that should inform and aid in the interpretation of observations during the eclipse.

It is our belief that the usage of the PSSFT model allows for better prediction of the surface magnetic field, especially at high latitudes enabled via surface plasma flux transport processes. Such models assimilating surface magnetic field data provide ideal boundary conditions for simulating and predicting the large-scale coronal structure. We note that our scheme is not perfect and the PFSS model for coronal field extrapolation has certain limitations. We have ignored any non-linear effects on the plasma transport in the surface flux evolution (whose impact is expected to be minimal during declining phases of the cycle). Also PFSS extrapolations for simulating the corona cannot account for current carrying large-scale sheared structures. Neither can we account for the impact of the heliospheric current sheet near the source surface and beyond, interactions with which can deflect the cusps of streamer belts and coronal plumes at larger distances from the Sun. Nonetheless, as outlined earlier, compared with more complex approaches this methodology is a more efficient means for predicting the global coronal structure fast and validating it in near real time with coronal magnetometry observations. Comparison of our PSSFT–PFSS models based prediction with posteclipse observations yield reasonable agreement with the observed large-scale corona during the 2019 July 2 solar eclipse, further reiterating the potential of this approach.

For our study, we have used a spherically symmetric density profile weighted by the nature of the magnetic field lines (open/closed) to compute the polarization vectors. This can be further improved by using observational plasma parameters e.g., density, pressure, temperature distribution in the corona which are currently not available. However, utilizing them wont affect the calculation of normalized Stokes vectors as it is dependent only on the line-of-sight magnetic field (Casini & Judge 1999; Lin & Casini 2000). Hence the Stokes vectors plotted in Figure 3(a) (marked as blue lines) and the degree of linear polarization plotted in Figure 3(b) are expected to remain unchanged.

Our analysis of the dependence of the inferred location of helmet streamers on the Sun's polar field strength (see also Riley et al. 2019) indicates that (projection) errors in polar observations can compromise the inference of coronal structures. In this context, either accurate surface flux transport simulations, or out-of-ecliptic plane observations of the Sun's polar field (such as that envisaged by the Solar Orbiter spacecraft) will be critical for coronal field simulations.

The coronal magnetic fields, which evolve in response to driving from the solar surface, govern spatial and temporal variations of the slow and fast components of the solar wind and heliospehric open flux. They also spawn solar flares and CMEs which have severe space weather impacts. Several future facilities are focusing on coronal magnetometry and gearing up to return coronal magnetic field measurements. These include the ground-based facilities DKIST and CoMP (Tomczyk et al. 2008). ISRO's Aditya-L1 space mission, currently under development, will fly the Solar Ultraviolet Imaging Telescope (SUIT; Tripathi et al. 2017) and the Variable Emission Line Coronagraph (VELC; Prasad et al. 2017) instruments which would simultaneously observe filament-prominence-arcade systems and coronal magnetic fields. We expect these multi-viewpoint, multi-wavelength observations to revolutionize the field of coronal magnetometry. Collective endeavors combining theoretical modeling with these coronal observations are expected to lead to refined data-driven operational forecasting models for solar activity induced space weather.

We are thankful to Sarah Gibson, Anthony Yeates, and the anonymous referee for helpful comments. CESSI is funded by the Ministry of Human Resource Development, Government of India, under the Frontier Areas of Science and Technology (FAST) scheme. S.D. acknowledges funding from the INSPIRE program of the Department of Science and Technology, Government of India. We acknowledge utilization of data from the NASA/SDO HMI instrument maintained by the HMI team and the Royal Greenwich Observatory/USAF-NOAA active region database compiled by David H. Hathaway. The authors acknowledge the use of data from the SOHO/LASCO and MLSO/COSMO K-Cor instruments. Additional graphics and text related to this prediction are available online at http://www.cessi.in/solareclipse2019.

Footnotes

Please wait… references are loading.
10.3847/1538-4357/ab6a91