Duncan Hendry Mackay

The Sun’s Open Flux plays a vital role in heliospheric physics and for the origin of the solar wind. A key problem faced by global magnetic field models of the solar corona is that they consistently underestimate the amount of Open Flux compared to in situ measurements. The purpose of the present paper is to consider if statistically averaged helicity condensation (SAHC) can reduce this discrepancy by increasing the simulated Open Flux found in potential field source surface (PFSS) models. We consider each Carrington rotation from 1646–2149 (1976 September to 2014 April). For each rotation, we hold the photosphere fixed, construct a PFSS model, and inject nonpotentiality into the coronal field by applying SAHC at different angular rotation rates, ω. When ω exceeds 2.5 × 10⁻⁶ s⁻¹, both ambient and sporadic enhancements of Open Flux are found within the Carrington rotation timescale. As ω increases, the amount of ambient Open Flux increases, along with the frequency of sporadic events. These enhancements are accompanied by increased footpoint area of open field lines, where the open field locations at each pole expand by up to 3^∘–6^∘ of latitude. A solar cycle dependence is found with higher increases in ambient Open Flux at solar maximum compared to solar minimum. Comparing our results with 27 day averaged OMNI spacecraft measurements, our highest value of ω = 10 × 10⁻⁶ s⁻¹ can resolve just under 50% of the discrepancy between in situ Open Flux measurements and PFSS models. Although helicity condensation cannot fully resolve the discrepancy, it can significantly contribute toward reducing it.

Whitepaper #287 in the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. Main topics: space weather applications; basic research. Additional topics: space weather research/applications/operations pipeline; solar physics; space-based missions/projects; […]

Whitepaper #333 in the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. Main topics: space weather applications; basic research. Additional topics: space weather research/operations/research loop; space weather research/applications/operations pipeline; […]

The solar corona exhibits a wide range of phenomena, from highly non-potential objects such as solar filaments to near-potential structures such as the open magnetic flux. For any global model to be useful in space weather applications, the model must on a single day reproduce all of these phenomena in the same simulation, using the same set of coronal physics and parameters. The purpose of the present paper is to evaluate whether the evolving magnetofrictional model can achieve this goal. Twenty-eight separate simulations are analyzed, where each tries to reproduce both highly non-potential and near-potential phenomena observed in the solar corona on the same day. This day is chosen to be 2015 March 20, the date of the solar eclipse. The study evaluates how the cadence of bipole data, ideal or nonideal coronal physics, and the variety of helicity injection mechanisms affect the accuracy of the simulations. Results show that significantly better agreement arises when using Advective Flux Transport (AFT) synoptic maps to drive the simulations, as compared to 27-day Carrington rotation synoptic maps. Using the nonideal effect of hyperdiffusion leads to the worst agreement with all coronal phenomena. Alternatively, when running either ideal or ohmic diffusion simulations with helicity condensation or bipoles with a self-helicity, a good agreement with both on-disk and limb structures can be found. This suggests that future studies aiming to simulate the corona and reproduce multiple phenomena on a given day should use data products such as AFT and avoid using the nonideal physics of hyperdiffusion.

Context. Our study investigates a mini-filament eruption associated with cancelling magnetic fluxes. The eruption originates from a small-scale loop complex commonly known as a coronal bright point (CBP). The event is uniquely recorded in both the imaging and spectroscopic data taken with the Interface Region Imaging Spectrograph (IRIS). Aims. The investigation aims to gain a better understanding of the physical processes driving these ubiquitous small-scale eruptions. Methods. We analysed IRIS spectroscopic and slit-jaw imaging observations as well as images taken in the extreme-ultraviolet channels of the Atmospheric Imaging Assembly (AIA) and line-of-sight magnetic-field data from the Helioseismic Magnetic Imager (HMI) on board the Solar Dynamics Observatory. As the observations can only indicate the possible physical processes at play, we also employed a non-linear force-free field (NLFFF) relaxation approach based on the HMI magnetogram time series. This allowed us to further investigate the evolution of the magnetic-field structures involved in the eruption process. Results. We identified a strong small-scale brightening as a micro-flare in a CBP, recorded in emission from chromospheric to flaring plasmas. The mini-eruption is manifested via the ejection of hot (CBP loops) and cool (mini-filament) plasma recorded in both the imaging and spectroscopic data. The micro-flare is preceded by the appearance of an elongated bright feature in the IRIS slit-jaw 1400 Å images, located above the polarity inversion line. The micro-flare starts with an IRIS pixel size brightening and propagates bi-directionally along the elongated feature. We detected, in both the spectral and imaging IRIS data and AIA data, strong flows along and at the edges of the elongated feature; we believe that these represent reconnection outflows. Both edges of the elongated feature that wrap around the edges of the erupting MF evolve into a J-type shape, creating a sigmoid appearance. A quasi-separatrix layer (QSL) is identified in the vicinity of the polarity inversion line by computing the squashing factor, Q, in different horizontal planes of the NLFFF model. Conclusions. This CBP spectro-imaging study provides further evidence that CBPs represent downscaled active regions and, as such, they may make a significant contribution to the mass and energy balance of the solar atmosphere. They are the sources of all range of typical active-region features, including magnetic reconnection along QSLs, (mini-)filament eruptions, (micro-)flaring, reconnection outflows, etc. The QSL reconnection site has the same spectral appearance as the so-called explosive events identified by strong blue- and red-shifted emission, thus providing an answer to an outstanding question regarding the true nature of this spectral phenomenon.

An important element of 3D data-driven simulations of solar magnetic fields is the determination of the horizontal electric field at the solar photosphere. This electric field is used to drive the 3D simulations and inject energy and helicity into the solar corona. One outstanding problem is the localisation of the horizontal electric field such that it is consistent with Ohm’s law. Yeates (Astrophys. J.836(1), 131, 2017) put forward a new “sparse” technique for computing the horizontal electric field from normal-component magnetograms that minimises the number of non-zero values. This aims to produce a better representation of Ohm’s law compared to previously used “non-sparse” techniques. To test this new approach we apply it to active region (AR) 10977, along with the previously developed non-sparse technique of Mackay, Green, and van Ballegooijen (Astrophys. J.729(2), 97, 2011). A detailed comparison of the two techniques with coronal observations is used to determine which is the most successful. Results show that the non-sparse technique of Mackay, Green, and van Ballegooijen (2011) produces the best representation for the formation and structure of the sigmoid above AR 10977. In contrast, the Yeates (2017) approach injects strong horizontal fields between spatially separated, evolving magnetic polarities. This injection produces highly twisted unphysical field lines with significantly higher magnetic energy and helicity. It is also demonstrated that the Yeates (2017) approach produces significantly different results that can be inconsistent with the observations depending on whether the horizontal electric field is solved directly or indirectly through the magnetic vector potential. In contrast, the Mackay, Green, and van Ballegooijen (2011) method produces consistent results using either approach. The sparse technique of Yeates (2017) has significant pitfalls when applied to spatially resolved solar data, where future studies need to investigate why these problems arise.

Context. Coronal mass ejections (CMEs) that exhibit weak or no eruption signatures in the low corona, known as stealth CMEs, are problematic as upon arrival at Earth they can lead to geomagnetic disturbances that were not predicted by space weather forecasters. Aims. We investigate the origin and eruption of a stealth event that occurred on 2015 January 3 that was responsible for a strong geomagnetic storm upon its arrival at Earth. Methods. To simulate the coronal magnetic field and plasma parameters of the eruption we use a coupled approach. This approach combines an evolutionary nonlinear force-free field model of the global corona with a MHD simulation. Results. The combined simulation approach accurately reproduces the stealth event and suggests that sympathetic eruptions occur. In the combined simulation we found that three flux ropes form and then erupt. The first two flux ropes, which are connected to a large AR complex behind the east limb, erupt first producing two near-simultaneous CMEs. These CMEs are closely followed by a third, weaker flux rope eruption in the simulation that originated between the periphery of AR 12252 and the southern polar coronal hole. The third eruption coincides with a faint coronal dimming, which appears in the SDO/AIA 211 Å observations, that is attributed as the source responsible for the stealth event and later the geomagnetic disturbance at 1 AU. The incorrect interpretation of the stealth event being linked to the occurrence of a single partial halo CME observed by LASCO/C2 is mainly due to the lack of STEREO observations being available at the time of the CMEs. The simulation also shows that the LASCO CME is not a single event but rather two near-simultaneous CMEs. Conclusions. These results show the significance of the coupled data-driven simulation approach in interpreting the eruption and that an operational L5 mission is crucial for space weather forecasting.

The coronal magnetic field evolution of 20 bipolar active regions (ARs) is simulated from their emergence to decay using the time-dependent nonlinear force-free field method of Mackay, Green, and van Ballegooijen (Astrophys. J. 729, 97, 2011). A time sequence of cleaned photospheric line-of-sight magnetograms, which covers the entire evolution of each AR, is used to drive the simulation. A comparison of the simulated coronal magnetic field with the 171 and 193 Å observations obtained by the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA), is made for each AR by manual inspection. The results show that it is possible to reproduce the evolution of the main coronal features such as small- and large-scale coronal loops, filaments and sheared structures for 80% of the ARs. Varying the boundary and initial conditions, along with the addition of physical effects such as Ohmic diffusion, hyperdiffusion and a horizontal magnetic field injection at the photosphere, improves the match between the observations and simulated coronal evolution by 20%. The simulations were able to reproduce the build-up to eruption for 50% of the observed eruptions associated with the ARs. The mean unsigned time difference between the eruptions occurring in the observations compared to the time of eruption onset in the simulations was found to be ≈5 hrs. The simulations were particularly successful in capturing the build-up to eruption for all four eruptions that originated from the internal polarity inversion line of the ARs. The technique was less successful in reproducing the onset of eruptions that originated from the periphery of ARs and large-scale coronal structures. For these cases global, rather than local, nonlinear force-free field models must be used. While the technique has shown some success, eruptions that occur in quick succession are difficult to reproduce by this method and future iterations of the model need to address this.

Context. We report on the third part of a series of studies on eruptions associated with small-scale loop complexes named coronal bright points (CBPs). Aims. A single case study of a CBP in an equatorial coronal hole with an exceptionally large size is investigated to expand on our understanding of the formation of mini-filaments, their destabilisation, and the origin of the eruption triggering the formation of jet-like features recorded in extreme ultraviolet (EUV) and X-ray emission. We aim to explore the nature of the so-called micro-flares in CBPs associated with jets in coronal holes and mini coronal mass ejections in the quiet Sun. Methods. Co-observations from the Atmospheric Imaging Assembly (AIA) and Helioseismic Magnetic Imager (HMI) on board the Solar Dynamics Observatory as well as GONG Hα images are used together with a non-linear force free field (NLFFF) relaxation approach, where the latter is based on a time series of HMI line-of-sight magnetograms. Results. A mini-filament (MF) that formed beneath the CBP arcade about 3−4 h before the eruption is seen in the Hα and EUV AIA images to lift up and erupt triggering the formation of an X-ray jet. No significant photospheric magnetic flux concentration displacement (convergence) is observed and neither is magnetic flux cancellation between the two main magnetic polarities forming the CBP in the time period leading to MF lift-off. The CBP micro-flare is associated with three flare kernels that formed shortly after the MF lift-off. No observational signature is found for magnetic reconnection beneath the erupting MF. The applied NLFFF modelling successfully reproduces both the CBP loop complex as well as the magnetic flux rope that hosts the MF during the build-up to the eruption. Conclusions. The applied NLFFF modelling is able to clearly show that an initial potential field can be evolved into a non-potential magnetic field configuration that contains free magnetic energy in the region that observationally hosts the eruption. The comparison of the magnetic field structure shows that the magnetic NLFFF model contains many of the features that can explain the different observational signatures found in the evolution and eruption of the CBP. In the future, it may eventually indicate the location of destabilisation that results in the eruptions of flux ropes.

The large field-of-view of the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) instrument onboard the PRoject for Onboard Autonomy 2 (PROBA2) spacecraft provides a unique opportunity to study extended coronal structures observed in the EUV in conjunction with global coronal magnetic field simulations. A global non-potential magnetic field model is used to simulate the evolution of the global corona from 1 September 2014 to 31 March 2015, driven by newly emerging bipolar active regions determined from Helioseismic and Magnetic Imager (HMI) magnetograms. We compare the large-scale structure of the simulated magnetic field with structures seen off-limb in SWAP EUV observations. In particular, we investigate how successful the model is in reproducing regions of closed and open structures, the scale of structures, and compare the evolution of a coronal fan observed over several rotations. The model is found to accurately reproduce observed large-scale, off-limb structures. When discrepancies do arise they mainly occur off the east solar limb due to active regions emerging on the far side of the Sun, which cannot be incorporated into the model until they are observed on the Earth-facing side. When such “late” active region emergences are incorporated into the model, we find that the simulated corona self-corrects within a few days, so that simulated structures off the west limb more closely match what is observed. Where the model is less successful, we consider how this may be addressed, through model developments or additional observational products.