ORIGINAL RESEARCH article
Sec. Space Physics
Exploiting White-Light Observations to Improve Estimates of Magnetic Connectivity
- 1IRAP, CNRS, CNES, Université Toulouse III–Paul Sabatier, Toulouse, France
- 2LDE3, DAp/AIM, CEA/IRFU, CNRS/INSU, Gif-sur-Yvette, France
The Solar Orbiter (SolO) and Parker Solar Probe missions have opened up new challenges for the heliospheric scientific community. Their proximity to the Sun and their high quality measurements allow us to investigate, for the first time, potential sources for the solar wind plasma measured in situ. More accurate estimates of magnetic connectivities from spacecraft to the Sun are required to support science and operations for these missions. We present a methodology to systematically compare coronal and heliospheric models against white-light (WL) observations. WL images from the SOlar and Heliospheric Observatory (SoHO) are processed to unveil the faint structures of the K-corona. Images are then concatenated over time and are projected into a Carrington synoptic map. Features of interest such as the Streamer Belt (SB) are reduced to simplified geometric objects. Finally, a metric is defined to rank models according to their performance against WL observations. The method has been exploited to reproduce magnetic sectors from WL observations. We tested our results against one year of in situ magnetic polarity measurements taken at near one AU from the Advanced Composition Explorer (ACE) and the Solar TErrestrial RElations Observatory (STEREO-A). We obtained a good correlation that emphasizes the relevance of using WL observations to infer the shape of the sector structure. We show that WL observations provide additional constraints to better select model parameters such as the input photospheric magnetic map. We highlight the capability of this technique to systematically optimize coronal and heliospheric models using continuous and near-real-time WL observations. Several relevant practical applications are discussed, which should allow us to improve connectivity estimates.
A new era in heliophysics has begun with the recent launches of the Solar Orbiter (SolO) and Parker Solar Probe (PSP) missions. The comprehensive instrumentation onboard these spacecraft and their proximity to the Sun will allow for detailed investigations of the relation between the solar wind properties measured in situ (IS) and their solar source regions observed by remote-sensing (RS) instruments. This IS-RS connection is typically established by considering magnetic field lines that connect to the same solar wind source region at the Sun. They can also act to channel energetic particles. Whichever the chosen model, Magneto-Hydro-Dynamic (MHD) or Potential Field Source Surface (PFSS) models, a common critical parameter is the input photospheric magnetic field map. Many magnetograms are now made available to the scientific community and can serve as the initial magnetic boundary condition for solar wind simulations. Some of them not only offer a single set but an ensemble of possibilities such as those computed by the Air Force Data Assimilative Photospheric Flux Transport (ADAPT) model, (see e.g. Arge et al., 2010; Arge et al., 2013). The ADAPT maps are global magnetograms of the photospheric magnetic flux and are produced using data assimilation techniques along with a magnetic flux transport model (see Worden and Harvey, 2000). Each ADAPT map has different realizations of the status of the photospheric magnetic field at a certain time. Selecting among so many magnetograms is a challenging task that is often done arbitrarily. In this work, we present a novel technique that it is designed to automatically test the reliability of the different realizations of the ADAPT maps. The developed method is based on a comparison of the global models of the solar atmosphere with the location of streamers in synoptic maps from white-light (WL) observations. The method is versatile by allowing us to perform such comparisons using any magnetograms and any type of global models.
The convenience of Carrington maps derived from WL observations has been demonstrated in many studies. Some of the earliest WL synoptic maps were constructed with K-corona images from the CORONASCOPE II and SOLWIND missions (Bohlin, 1970; Wang and Sheeley, 1992), as well as from the Mauna Loa Solar Observatory in Hawaï (Hansen et al., 1976). Later on, synoptic observations from the (WL) coronagraph onboard Skylab have been exploited to model the large scale coronal density structures by assuming that the peaks in brightness mark the location of the heliospheric current sheet (HCS) (Guhathakurta et al., 1996). The continuous monitoring of the solar corona by the Solar and Heliospheric Observatory (SoHO) has then enabled a more systematic comparison between the location of streamers with the magnetic topology of the solar corona derived from PFSS calculations (Wang et al., 1998; Wang et al., 2000; Wang et al., 2007) and global coronal models (Gibson et al., 2003; Thernisien and Howard, 2006; de Patoul et al., 2015; Pinto and Rouillard, 2017). Rotational tomography techniques have been developed recently to correct for line of sight effects and convert WL observations into synoptic density maps (Morgan and Cook, 2020). Other techniques have involved the combination of coronagraph images obtained from multiple vantage points, (e.g. SoHO and STEREO) simultaneously to derive synchronic circumsolar maps of helmet streamers (Sasso et al., 2019). WL Carrington maps have also been exploited to track the origin of the solar wind measured in situ at PSP (see e.g., Rouillard et al., 2020a; Rouillard et al., 2020c; Griton et al., 2021) or to study the fine structure of coronal rays observed remotely at PSP (see Poirier et al., 2020). The deflection of streamers during solar storms can also be investigated using WL Carrington maps (see e.g., Kouloumvakos et al., 2020). The assumption in these comparisons is that the bright helmet streamers mark the location of the (HCS). We will make a similar assumption throughout this work and in particular we will test its reliability in Section 3.
In Section 2, we show in detail the methodology followed to perform the systematic assessment of coronal models against WL observations. We exploit our automated detection of the streamer belt (SB) to infer the distribution of magnetic sectors by assuming that the core of the SB hosts the HCS. In Section 3, we then directly compare the sector structure inferred at the Sun with the one measured in situ by two space-based observatories orbiting one AU. Finally, in Sections 4, 5 we highlight the successes and also discuss the limitations of our automated technique for several practical applications for space weather related topics. We also discuss a connection with PSP and SolO observations.
2 Data and Methods
In this section we present the technique developed to compare coronal models against WL observations using an automated and systematic approach. We start with a description of the methods used to pre-process the coronagraph images (sub Section 2.1) and to construct the WL synoptic maps (sub Section 2.2). In sub Section 2.3 we present a technique used to identify and characterize streamer signatures in WL synoptic maps. Finally, in sub Section 2.4 we define a metric that can be used to score coronal models against WL observations and we compare a coronal model with observations using a test case.
2.1 Pre-Processing the White-Light Images
In this study we exploit WL observations from the Large Angle and Spectrometric Coronagraph (LASCO) C2 (LC2) coronagraph on-board SoHO. SoHO is a collaborative ESA-NASA mission launched on December 2, 1995, and which is still operating at the Lagrange point L1. The SoHO spacecraft is equipped with an imaging suite composed of three coronagraphs (see Brueckner et al., 1995), two of which (namely the LASCO, C2, and C3) are still operating and observe continuously the corona in visible light. The LC2 coronagraph has a field-of-view (FOV) extending from
The WL emission observed by the coronagraphs comes from several distinct contributions. The dominant one comes from the photospheric light that is scattered by dust particles; this forms the F-corona. Dust particles are ubiquitous around the Sun but a number of theoretical studies predict the presence of a “dust-free” zone (see e.g., Lamy, 1974; Mukai et al., 1974), which might have been detected for the first time by the Wide-field Imager for Solar PRobe (WISPR) onboard PSP (see Howard et al., 2019). Another contribution is the photospheric light scattered by coronal electrons; this forms the so-called K-corona that is of interest for this study. WL images of the K-corona can help to locate streamers and streamer rays, where a subset of the slow solar wind forms. The brightest signatures in WL images come from regions of dense plasma such as the streamers.
Coronagraphic images must be processed beforehand to remove the F-corona. An estimate of the F-corona brightness distribution is derived by considering LC2 images over several days to a month, this image is then subtracted from the raw image to reveal the fine details of the K-corona. This technique has been used for several decades and is adequate for spacecraft moving slowly along their orbit such as SoHO located at one AU. It is however no longer suitable for fast-moving observatories such as WISPR onboard PSP, requiring more sophisticated techniques (see Stenborg et al., 2018; Hess et al., 2020). For SoHO LASCO we exploit the “monthly min” background images that are provided through https://sohowww.nascom.nasa.gov/.
In Figure 1 we show an example of a processed LASCO C2 image on the 2020/10/03 at 12:00. This image reveals the K-corona composed on that day of two large streamer rays one off the East (left) and West (right) limbs. These streamer rays emanate from the tip of helmet streamers and extend to the high corona.
FIGURE 1. A processed LASCO C2 image on the 2020-10-03 at 12:00. The green circle depicts the band of pixels that are extracted at
2.2 Construction of the White-Light Synoptic Maps
To visualize the complete SB we construct a synoptic map given in a heliographic latitude vs. Carrington longitude format by considering LC2 images obtained over a complete solar rotation. The SB is the name given to the 3-D shape formed by streamer rays over the full solar rotation. Each LC2 image is first projected onto the helioprojective sphere, which is defined as the sphere of radius the Sun—observer distance and with the observer at its center. Once a coordinate is associated to each image pixel, a band of pixels of circular shape is extracted from each image at a specific heliocentric radial distance (green line in Figure 1), which will be called hereafter the height of the Carrington map. After a coordinate transformation into the Carrington frame, the extracted pixels are mapped onto a latitude vs. longitude map (green line in Figure 2). As images are concatenated over time, pixels fill up the synoptic map until the latter covers a full solar rotation. The full synoptic map can then be continually updated as new images become available, thus providing a dynamic synoptic map.
FIGURE 2. White-light synoptic (or Carrington) map at a height of
The height at which pixels are extracted from the images plays a significant role into building a synoptic map. A height too close to the inner edge of the coronagraph would give a signal that is too saturated, whereas WL features are not so well resolved near the outer edge. Furthermore, streamer rays get narrower with increasing heliocentric distance, resulting in a change of the SB thickness in synoptic maps. A synoptic map with a SB that does not evolve anymore with increasing height is therefore preferred. This condition is usually satisfied at a height greater than
Although WL synoptic maps such as Carrington maps have been widely used in the past, their interpretation can be tricky. Relating bright WL signatures to dense coronal structures is not trivial. The major difficulty arises from the fact that a WL detector integrates light emitted from a range of locations extending along the line-of-sight. This induces some uncertainty in associating specific structures modeled in 3-D, such as streamers, with their WL signature. The combined effects of the inverse square falloff with heliocentric radial distance, of both the solar wind density and of the photospheric light intensity, mean that most of the WL contribution is expected to come from electrons located near the “Thomson sphere” (e.g. Vourlidas and Howard, 2006; Howard and Tappin, 2009). This sphere is geometrically defined as a sphere of diameter, the Sun—observer distance, and whose center is located at half this distance. A later study from Howard and DeForest (2012) showed that contributions to the total brightness are actually distributed over a broader region called the “Thomson plateau”, with the “Thomson sphere” standing in the middle. A recent study by Morgan and Cook (2020) tackled this line-of-sight integration problem by using rotational tomography which exploits the Sun’s rotation to resolve spatial density structures. This is a highly sophisticated method that requires additional modeling assumptions. We do not exploit such a method in order to keep our automated technique purely based on WL observations without any additional modeling aspects. We are aware of these line-of-sight related uncertainties, which are expected to occur mostly at places of large folds in the SB and mainly in the longitudinal direction.
The WL synoptic map shown in Figure 2 depicts a coronal configuration typical of a low solar activity, with a SB that remains flat most of the time. Only a small folding of the SB can be seen at 70 deg Carrington longitude along with a drift away of the SB toward higher latitudes. PFSS reconstructions of the coronal magnetic field suggest that, in this case, this bright excursion is part of an arch from a pseudo-streamer. But other inverted arch features, called “bananas” in Gibson et al. (2003), can also be observed in WL synoptic maps from the projection of non-equatorial streamers due to coronal rotation (see also Wang and Sheeley, 1992; Wang et al., 2000; Wang et al., 2007). These projection effects are particularly accentuated for observatories that are located much closer to the Sun such as WISPR onboard PSP (see e.g., Liewer et al., 2019; Poirier et al., 2020).
2.3 Identifying the Features of Interest in Synoptic Maps
Our aim is to exploit WL observations to determine if the HCS derived by a particular magnetic model of the solar atmosphere reproduces the position of bipolar streamers well. Therefore, our algorithm relies on the clear identification of coronal rays that originate from the tip of bipolar streamers. However, coronal rays observed in coronagraph images can have other origins. For instance they can emanate from other types of streamers called pseudo-streamers, which are particularly common during high solar activity. In contrast to bipolar streamers, pseudo-streamers form in open field regions of the same polarity and therefore cannot be related to the origin of the HCS of interest here. Fortunately, streamers and pseudo-streamers usually have different radial extensions in the corona, and there can be used to discriminate between the two. Pseudo-streamers are expected to be smaller structures compared to streamers, with a cusp located lower in the corona than bipolar streamers (Wang et al., 2007). We have found that at a heliocentric distance of
We also apply a 2-D median filter to reduce the effects of “hot pixels” and small artifacts related, (e.g. remaining cosmic rays), to ease the automatic identification and tracing of the SB. This tracing is performed by using a specific brightness threshold value (dark solid line in Figure 2). During the building process of the WL map, a dynamic scaling of the brightness is performed. Brightness is normalized by the maximal value of the current band of pixels that is extracted from the image. A single threshold value can therefore be set up to trace the SB contour in a consistent manner over the full solar rotation. Another advantage is that it allows a coherent scientific interpretation since the maximal brightness is likely to take place around the HCS.
The SB identification can be compromised in certain regions where WL emissions from the SB become faint as seen near longitudes 60° and 200° in Figure 2. This dimming happens when solar rotation brings the local tangent plane to the belt at high angles to the line of sight of the observer. Indeed, when the angle between the plane’s normal vector and the line of sight decreases, the section of the belt scattering photons back to the observer becomes much smaller. This effect becomes very common when the Sun gets more active with a dominant quadrupolar structure and the belt undergoes significant latitudinal excursions. In such configurations, the SB appears no longer continuous but rather fragmented, visible as separate SB bands in Carrington maps. These regions of large folding are difficult to address in an automated manner and they are consequently removed from the comparison with the models.
As discussed earlier, coronal rays from streamers are expected to be a good proxy to visualize where solar magnetic field reversals take place in the outer corona and inner heliosphere. The SB identified above already outlines the region where we expect those reversals. But we can go further by assuming that magnetic polarity changes are likely to occur where the plasma is densest within the core of the SB, and hence where the streamer total brightness reaches a maximum. Indeed, the HCS is usually measured in situ during the passage of the thin Heliospheric Plasma Sheet (HPS), a region of the solar wind that is usually associated with an order of magnitude increase in density (Winterhalter et al., 1994). The thickness of the HPS measured in situ tends to be thinner than streamers; the first WISPR images taken by PSP suggest that the HPS is a smaller but very bright region of the streamer located right at its core (Poirier et al., 2020; Rouillard et al., 2020c). An algorithm was therefore set up to trace the line passing through all the peaks of total brightness. This line covers a full solar rotation and is hereafter called the Streamer-Maximum-Brightness (SMB) line (dark dashed line in Figure 2). In Section 3, we will show a very good agreement between the estimated times when spacecraft are expected to cross the SMB and the crossings of HCS detected in situ.
2.4 Defining a Metric to Compare the Coronal Models With Observations
The technique presented in this paper aims to evaluate global coronal models against WL observations in an automated and systematic manner. For this purpose one needs to define a metric, a single quantity that evaluates how well the models agree with the WL observations. We define a metric that is comprehensive and which aims to compare distinct models to each other. We will illustrate the final metric by applying our automated technique to several extrapolations produced by a PFSS model. We choose a PFSS model rather than MHD models for simplicity, but the automated technique is applicable to any other types of models such as full 3-D MHD as long as they provide the location of the HCS (see e.g. Badman et al., 2021).
The position of the neutral line (NL) is therefore a necessary input from the model to be evaluated. For a PFSS model, the NL location is obtained at the outer boundary of the extrapolation domain, which is the source surface where magnetic field lines are assumed to be radial and open into the interplanetary medium. The height of this source surface is generally
Angular distances between the model NL and the SMB line are calculated for a maximum of 180 points along the streamer belt. The SB contour is then exploited in order to normalize these angular distances by the observed thickness of the streamer. In regions where the outline of the SB is not clearly determined, usually due to large folds (see sub Section 2.2), the normalized angular distances are not computed. This ensures that uncertain regions where the SMB is poorly observed do not affect the final result.
The next step is to combine the modeled NL—SB angular distances derived over 360 degrees of longitude into a global score that tells us how well the model NL matches the SB location. A mean value
In addition to the mean value
A final parameter called the confidence score, C, is computed by combining both
The gain function has been designed so that a mean value
In Figure 3 we compare several PFSS extrapolations against the SMB line and SB contour identified earlier in Figure 2. Results from the metric defined above are given in the legend including the mean score
FIGURE 3. Comparison of several NLs against the SMB line and SB contour identified in Figure 2. NLs computed from the PFSS model using different
FIGURE 4. (A): 2-D colored map of confidence scores C as a function of both the gain and penalty scores. (B) and (C): gain and penalty scores as a function of the mean value
We illustrate our scoring technique in Figure 5. Confidence scores of multiple PFSS extrapolations can easily be compared with each other, allowing a quick selection of the PFSS model parameters that best fit the SB location on a specific date. In this example, the GONG/mrzqs (“ZGONG”) magnetogram with a RSS at
FIGURE 5. Summary plot of PFSS models performances according to WL, using GONG/mrzqs (labeled “ZGONG”) and GONG/ADAPT (labeled “AGONG_XX”) magnetograms and three distinct
3 Correlation With Magnetic Sectors Measured in situ at One AU
As discussed above, the brightest WL emissions are likely associated with the very dense plasma located in the HPS that encloses the HCS. This assumption can be tested by connecting WL observations with direct in situ measurements of the magnetic sector structure. In the early 70 s, a strong correlation was already found between the brightest regions of the K-corona and the sector boundaries measured in situ at one AU (Hansen et al., 1974; Howard and Koomen, 1974). We check again this correlation by developing a method to assess in a statistical manner, the magnetic field reversals measured in situ vs. the SMB lines derived from the WL synoptic maps. We exploited one year of magnetic field measurements from the Advanced Composition Explorer (ACE) and the Solar TErrestrial RElations Observatory (STEREO-A).
A sector structure can be derived from the WL synoptic maps by assuming that the SMB line is the location of the neutral line at
WL-derived magnetic sector maps can then be compared to sectors measured directly in situ. In order to compare one year of data in a convenient way we bin the data in Bartels map format. Such maps are shown in Figure 6 for STEREO-A [panels (A) and (C)] and ACE [panels (B) and (D)]. They have been constructed using 16-h averaged in situ measurements [panels (C) and (D)] and WL observations [panels (A) and (B)]. Bartels maps display successive Carrington rotations along the ordinate, and time during a specific rotation along the x axis (here the Carrington longitude). With a Solar rotation period of 27.27 days as viewed from Earth (synodic period), the map is progressively filled up from right to left and from top to bottom. Bartels maps are useful to visualize at a glance magnetic sectors and their variability over consecutive Carrington rotations.
FIGURE 6. Bartel’s diagrams constructed with the SMB lines from the WL observations, for STEREO-A (A) and ACE(B). Bartel’s diagrams using direct in situ polarity measurements taken over the full year 2019 at STEREO-A (C) and ACE(D). Blue and red colors indicate positive and negative polarities respectively (dark colors correspond to data gaps).
Data are color-coded in blue (and red) for positive (and negative) polarities. Figure 6 brings out a good correspondence between in situ and WL derived Bartel’s diagrams with very similar sector structures, with
The presented metric can be of interest for many direct applications. In this section we present preliminary works that illustrate several of these possibilities.
4.1 Understanding in Situ Measurements at PSP
Such a comparison can be repeated in the context of the SolO and PSP missions. Their proximity to the Sun offer the opportunity to probe in greater detail the relation between HCS crossings measured in situ and their origin in WL. Since those measurements have been taken in the heliosphere, a projection of these data back to the solar corona is needed. We use the same ballistic mapping as employed before, with the nominal Parker spiral extending from PSP down to the height of the WL synoptic map at
FIGURE 7. History of magnetic polarities measured in situ at PSP from 21-10-2018 to 04-12-2018. In the background a WL synoptic map from LASCO C2 at
4.2 Constraining PFSS Coronal Models in a Systematic Manner
PFSS coronal models usually set a source surface height at
4.3 Constraining MHD Simulations
One of the most critical parameters in any coronal or heliospheric model is the magnetogram used for the inner boundary of the simulated domain. Whatever the accuracy of the model, the magnetogram will have a major impact on the final outcome. The technique presented in this work can be used to select which combination of magnetograms and models provide the best modeling of the streamer belt. MHD models are more challenging to operate because more constraints are required compared to PFSS models. In addition to the magnetic field, one needs to set the plasma bulk properties at the inner boundary of the simulation domain. The well-known Wang-Sheeley-Arge (WSA) model established an empirical law to estimate the solar wind velocity based on the expansion factor of open magnetic field lines between the photosphere and the source surface (see Wang and Sheeley, 1990; Arge and Pizzo, 2000). We have recently exploited our SMB identification technique to derive an alternative empirical law relating the angular distance of a point situated anywhere on the Carrington map to the SMB, to the plasma density and velocity that we should expect. In situ measurements from the Ulysses, ACE and WIND missions have been exploited to build empirical laws for the solar wind density and velocity, function of the angular distance to the SMB, this will be presented in a future study. In this relation, regions situated close to the SMB will produce a slower and denser solar wind than regions located far from the SMB. Examples of such derived maps are given in Figure 8 and have been injected as solar wind density and velocity 2-D maps in MHD heliospheric simulations extending typically from 21.5 to 215 solar radii such as the ENLIL model (see Odstrcil, 2003). The inner boundary of the MHD simulation can then be updated on a, (e.g., daily) basis as new WL observations are available and the algorithm presented here updates the position of the SMB. Results on this topic are still preliminary and will be presented in a dedicated study. A major benefit of this method is that it allows to dynamically constrain heliospheric models, therefore improving forecast capabilities for space weather purposes.
FIGURE 8. 2-D maps of solar wind density (A) and radial velocity (B) determined from the SMB line extracted on the 2019-06-15 at 12:00 UTC from LASCO-C2 WL observations at
4.4 Improving Sun to S/C Connectivity Estimates
Coronal and heliospheric models are essential in the context of SolO and PSP science because they allow us to estimate the possible origins of the solar wind plasma measured in situ at those spacecraft or connect energetic particles to potential accelerators at the Sun, (e.g. flares, shocks). A web-based service tool has been put in place by the Solar Orbiter Data Analysis Working Group (MADAWG) to provide support to SolO operations (Rouillard et al., 2020b). This tool, called the Magnetic Connectivity Tool, is directly accessible to the public at this website1. Forecast and past magnetic connectivities among many observatories such as PSP, SolO, BepiColombo, STEREO, and ACE are available via this web interface. Different magnetograms, coronal and heliospheric models can also be specified by the user. Recently, we implemented the technique presented in this paper to automatically select the combination of PFSS reconstructions—magnetograms that best fit WL observations to improve connectivity estimates. The selection is performed among several distinct sources of magnetograms such as those produced by the National Solar Observatory (NSO) and the Wilcox Solar Observatory (WSO), including those computed with the ADAPT model.
The most critical step in the methodology presented in this work is having a correct identification of the SMB line. We have noticed several difficulties that can make this identification challenging and sometimes even produce erroneous estimates.
We especially pointed out cases where there are large folds in the SB. Since our method does not resolve the line-of-sight (as in tomographical approaches), our metric can not give an accurate estimate of the HCS location where the SB is significantly warped. Nearby SB folds, inverted arches can also appear which result from the projection of non-equatorial streamer rays along the line-of-sight. Regions of large warping of the SB are in general omitted in our metric because their WL signature is too faint to be detected by the algorithm. However, we plan to incorporate an actual detection of SB folds so that they can be systematically removed from our metric.
Furthermore, highly structured coronal configurations that are typically observed during Solar maximum may require further adjustments of the identification technique, which has not yet been fully tested for those particular cases. Such extreme cases where a significant fraction of the SB is highly warped may affect the reliability of our metric.
In addition, the presence of pseudo-streamers can make the identification of the SMB ambiguous because they can be on occasion brighter than the bipolar streamers, especially in regions where the SB is warped. Distinguishing pseudo-streamers from bipolar streamers is therefore tricky. We found that lifting the height of the WL synoptic maps to
Another approach to separate pseudo-streamers from bipolar streamers could be designed by exploiting in situ measurements, as was done in Section 3 and in sub Section 4.1. A combination of our WL-based technique with a back-mapping of the sector structure measured in situ could be devised in an automated manner. It was attempted on a case by case basis by Badman et al. (2021) and showed some great potential. It can also provide a local assessment of the models and so diminish the impact of the global score provided by the WL-based approach. This is a possibility that we will attempt in the future.
Finally, another important point relates to the cadence with which WL Carrington maps are updated. We have found that discrepancies between the SMB line and the sector structure measured in situ from ACE, STEREO (see Section 3) or PSP (see sub Section 4.1) are often caused by the updating of the WL synoptic maps from off-limb observations. The rate of update in the map at a particular Carrington longitude and from a single viewpoint is set by the time a particular longitude rotates from the East to the West limb. From one AU this imposes an update every 14 days. This cadence can be greatly increased by combining coronagraphic observations from other vantage points than LASCO. This was the motivation for the development of merged synoptic maps from SoHO LASCO and STEREO SECCHI as performed by Sasso et al. (2019). Our next step is to incorporate these synchronic maps built from several vantage points.
We introduced a novel automated technique that provides an important constraint to coronal and heliospheric models based on the exploitation of WL observations. WL images from the LASCO C2 coronagraph are processed and converted into synoptic maps called Carrington maps, which are dynamically updated as new images are available. Features of interest such as the SB and the SMB are automatically identified by an algorithm. Future applications of this technique will include other vantage points including data from STEREO, SolO, PSP and the Polarimeter to Unify the Corona and Heliosphere (PUNCH). Reducing the initial synoptic maps to these two simple geometric elements allows us to perform direct comparisons with coronal models. Therefore a metric has been defined to quantify these comparisons with a score that allows us to assess automatically the quality of the models with respect to WL observations.
For illustrative purposes we applied this automated technique to PFSS models only and on a subset of magnetograms. But the technique can easily be applied to other coronal models including MHD models, and with other sources of photospheric magnetic maps.
Finally, we emphasized the potential of this technique for several practical applications. By providing additional and systematic constraints to coronal models, connectivity estimates can be improved that will support SolO and PSP related science topics. By providing additional constraints to MHD heliospheric models, the WL based optimization technique can also help to improve forecast capabilities for space weather applications.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
NP, AR, and AK developed the pipeline from the SB identification in WL synoptic maps to the final models scoring. AK developed the pipeline that processes raw white-light images and produces Carrington synoptic maps. AP, MI, AR, and NP performed the statistical comparison between in situ and white-light derived magnetic sectors. NF produced Figure 7 (and the corresponding animation) of PSP in situ magnetic polarities vs. WL synoptic maps. RP, VR, AR, and NP carried on the preliminary tests using the SMB line to constrain the plasma bulk properties at the inner boundary of MHD heliospheric simulations. RP and MI developed the PFSS model based on Wang and Sheeley (1992) using a spherical harmonics method. MI and MA developed the code and web interface of the Magnetic Connectivity Tool. NP wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
AK acknowledge financial support from the ANR project SLOW_SOURCE (ANR-18-ERC1-0006-01), COROSHOCK (ANR-17-CE31-0006-01), and FP7 HELCATS project https://www.helcats-fp7.eu/ under the FP7 EU contract number 606692. The work of AR, NP, VR and MA was funded by the ERC SLOW_SOURCE project (SLOW_SOURCE - DLV-819189). The IRAP team acknowledges support from the French space agency (Centre National des Etudes Spatiales; CNES; https://cnes.fr/fr) that funds the plasma physics data center (Centre de Données de la Physique des Plasmas; CDPP; http://cdpp.eu/), the Multi Experiment Data & Operation Center (MEDOC; https://idoc.ias.u-psud.fr/MEDOC) and the space weather team in Toulouse (Solar-Terrestrial Observations and Modelling Service; STORMS; https://stormsweb.irap.omp.eu/). This includes funding for the data mining tools AMDA (http://amda.cdpp.eu/), CLWEB (http://clweb.cesr.fr/) and the propagation tool (http://propagationtool.cdpp.eu).
Conflict of Interest
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.
We acknowledge the scientific and engineering teams who designed, built and share the data of the SoHO LASCO C2, ACE MAG, STEREO-A MAG and PSP FIELDS instruments. The authors are grateful to the Referees for constructive comments and recommendations which helped to improve the readability and quality of the paper.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2021.684734/full#supplementary-material
Arge, C. N., Henney, C. J., Hernandez, I. G., Toussaint, W. A., Koller, J., and Godinez, H. C. (2013). “Modeling the corona and Solar Wind Using ADAPT Maps that Include Far-Side Observations,” in Solar Wind 13. Editors G. P. Zank, J. Borovsky, R. Bruno, J. Cirtain, S. Cranmer, H. Elliottet al. (vol. 1539 of American Institute of Physics Conference Series), 11–14. doi:10.1063/1.4810977
Arge, C. N., Henney, C. J., Koller, J., Compeau, C. R., Young, S., MacKenzie, D., et al. (2010). “Air Force Data Assimilative Photospheric Flux Transport (ADAPT) Model,” in Twelfth International Solar Wind Conference. Editors M. Maksimovic, K. Issautier, N. Meyer-Vernet, M. Moncuquet, and F. Pantellini (vol. 1216 of American Institute of Physics Conference Series), 343–346. doi:10.1063/1.3395870
Badman, S. T., Bale, S. D., Martínez Oliveros, J. C., Panasenco, O., Velli, M., Stansby, D., et al. (2020). Magnetic Connectivity of the Ecliptic Plane within 0.5 au: Potential Field Source Surface Modeling of the First Parker. Solar Probe Encounter 246, 23. doi:10.3847/1538-4365/ab4da7
de Patoul, J., Foullon, C., and Riley, P. (2015). 3D Electron Density Distributions in the Solar Corona during Solar Minima: Assessment for More Realistic Solar Wind Modeling. 814, 68. doi:10.1088/0004-637X/814/1/68
Gibson, S. E., Foster, D. J., Guhathakurta, M., Holzer, T., and Cyr, St. O. C. (2003). Three-dimensional Coronal Density Structure: 1. Model. J. Geophys. Res. (Space Physics) 108, 1444. doi:10.1029/2003JA009994
Howard, R. A., Vourlidas, A., Bothmer, V., Colaninno, R. C., DeForest, C. E., Gallagher, B., et al. (2019). Near-Sun Observations of an F-corona Decrease and K-corona fine Structure. 576, 232–236. doi:10.1038/s41586-019-1807-x
Kouloumvakos, A., Rouillard, A. P., Share, G. H., Plotnikov, I., Murphy, R., Papaioannou, A., et al. (2020). Evidence for a Coronal Shock Wave Origin for Relativistic Protons Producing Solar Gamma-Rays and Observed by Neutron Monitors at Earth. 893, 76. doi:10.3847/1538-4357/ab8227
Liewer, P., Vourlidas, A., Thernisien, A., Qiu, J., Penteado, P., Nisticò, G., et al. (2019). Simulating White Light Images of Coronal Structures for WISPR/Parker Solar Probe: Effects of the Near-Sun Elliptical Orbit. 294, 93. doi:10.1007/s11207-019-1489-4
Panasenco, O., Velli, M., D’Amicis, R., Shi, C., Réville, V., Bale, S. D., et al. (2020). Exploring Solar Wind Origins and Connecting Plasma Flows from the Parker Solar Probe to 1 au: Nonspherical Source Surface and Alfvénic Fluctuations. 246, 54. doi:10.3847/1538-4365/ab61f4
Poirier, N., Kouloumvakos, A., Rouillard, A. P., Pinto, R. F., Vourlidas, A., Stenborg, G., et al. (2020). Detailed Imaging of Coronal Rays with the Parker Solar Probe. 246, 60. doi:10.3847/1538-4365/ab6324
Rouillard, A. P., Kouloumvakos, A., Vourlidas, A., Kasper, J., Bale, S., Raouafi, N.-E., et al. (2020a). Relating Streamer Flows to Density and Magnetic Structures at the Parker Solar Probe. 246, 37. doi:10.3847/1538-4365/ab579a
Rouillard, A. P., Pinto, R. F., Vourlidas, A., De Groof, A., Thompson, W. T., Bemporad, A., et al. (2020b). Models and Data Analysis Tools for the Solar Orbiter mission. 642, A2. doi:10.1051/0004-6361/201935305
Rouillard, A. P., Poirier, N., Lavarra, M., Bourdelle, A., Dalmasse, K., Kouloumvakos, A., et al. (2020c). Modeling the Early Evolution of a Slow Coronal Mass Ejection Imaged by the Parker Solar Probe. 246, 72. doi:10.3847/1538-4365/ab6610
Sasso, C., Pinto, R. F., Andretta, V., Howard, R. A., Vourlidas, A., Bemporad, A., et al. (2019). Comparing Extrapolations of the Coronal Magnetic Field Structure at 2.5 R with Multi-Viewpoint Coronagraphic Observations. 627, A9. doi:10.1051/0004-6361/201834125
Keywords: white-light imagery, modeling, space weather, sun: slow solar wind, sun: magnetic fields, sun: coronal streamers
Citation: Poirier N, Rouillard AP, Kouloumvakos A, Przybylak A, Fargette N, Pobeda R, Réville V, Pinto RF, Indurain M and Alexandre M (2021) Exploiting White-Light Observations to Improve Estimates of Magnetic Connectivity. Front. Astron. Space Sci. 8:684734. doi: 10.3389/fspas.2021.684734
Received: 23 March 2021; Accepted: 05 May 2021;
Published: 25 May 2021.
Edited by:Zoltan Voros, Austrian Academy of Sciences, Austria
Reviewed by:Pete Riley, Predictive Science, United States
Gordon James Duncan Petrie, National Solar Observatory, United States
Huw Morgan, Aberystwyth University, United Kingdom
Copyright © 2021 Poirier, Rouillard, Kouloumvakos, Przybylak, Fargette, Pobeda, Réville, Pinto, Indurain and Alexandre. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nicolas Poirier, firstname.lastname@example.org