Showing 1–50 of 50 results for author: Mackay, D H

Exploring the Solar Poles: The Last Great Frontier of the Sun

Authors: Dibyendu Nandy, Dipankar Banerjee, Prantika Bhowmik, Allan Sacha Brun, Robert H. Cameron, S. E. Gibson, Shravan Hanasoge, Louise Harra, Donald M. Hassler, Rekha Jain, Jie Jiang, Laurène Jouve, Duncan H. Mackay, Sushant S. Mahajan, Cristina H. Mandrini, Mathew Owens, Shaonwita Pal, Rui F. Pinto, Chitradeep Saha, Xudong Sun, Durgesh Tripathi, Ilya G. Usoskin

Abstract: Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodyna… ▽ More Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere. △ Less

Submitted 30 December, 2022; originally announced January 2023.

Comments: This White Paper was submitted in 2022 to the United States National Academies Solar and Space Physics (Heliophysics) Decadal Survey

A Statistical Comparison of EUV Brightenings Observed by SO/EUI with Simulated Brightenings in Non-potential Simulations

Authors: Krzysztof Barczynski, Karen A. Meyer, Louise K. Harra, Duncan H. Mackay, Frederic Auchere, David Berghmans

Abstract: The High Resolution Imager (HRI_EUV) telescope of the Extreme Ultraviolet Imager (EUI) instrument onboard Solar Orbiter has observed EUV brightenings, so-called campfires, as fine-scale structures at coronal temperatures. The goal of this paper is to compare the basic geometrical (size, orientation) and physical (intensity, lifetime) properties of the EUV brightenings with regions of energy dissip… ▽ More The High Resolution Imager (HRI_EUV) telescope of the Extreme Ultraviolet Imager (EUI) instrument onboard Solar Orbiter has observed EUV brightenings, so-called campfires, as fine-scale structures at coronal temperatures. The goal of this paper is to compare the basic geometrical (size, orientation) and physical (intensity, lifetime) properties of the EUV brightenings with regions of energy dissipation in a non-potential coronal magnetic field simulation. In the simulation, HMI line-of-sight magnetograms are used as input to drive the evolution of solar coronal magnetic fields and energy dissipation. We applied an automatic EUV brightening detection method to EUV images obtained on 30 May 2020 by the HRI_EUV telescope. We applied the same detection method to the simulated energy dissipation maps from the non-potential simulation to detect simulated brightenings. We detected EUV brightenings with density of 1.41x10^{-3} brightenings/Mm^2 in the EUI observations and simulated brightenings between 2.76x10^{-2} - 4.14x10^{-2} brightenings/Mm^2 in the simulation, for the same time range. Although significantly more brightenings were produced in the simulations, the results show similar distributions of the key geometrical and physical properties of the observed and simulated brightenings. We conclude that the non-potential simulation can successfully reproduce statistically the characteristic properties of the EUV brightenings (typically with more than 85% similarity); only the duration of the events is significantly different between observations and simulation. Further investigations based on high-cadence and high-resolution magnetograms from Solar Orbiter are under consideration to improve the agreement between observation and simulation. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 23 pages, 8 figures; Paper accepted to publication in Solar Physics

Eruptions from coronal bright points: A spectroscopic view by IRIS of a mini-filament eruption, QSL reconnection, and reconnection-driven outflows

Authors: Maria S. Madjarska, Duncan H. Mackay, Klaus Galsgaard, Thomas Wiegelmann, Haixia Xie

Abstract: The present 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 IRIS. We analyse IRIS spectroscopic and slit-jaw imaging observations as well as images taken in the extreme-ul… ▽ More The present 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 IRIS. We analyse IRIS spectroscopic and slit-jaw imaging observations as well as images taken in the extreme-ultraviolet channels of AIA, and line-of-sight magnetic-field data from HMI onboard the SDO. We also employ an NLFFF relaxation approach based on the HMI magnetogram time series. We identify a strong small-scale brightening as a micro-flare in a CBP. The mini-eruption manifests with 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 A 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 detect in both the spectral and imaging IRIS data and AIA data, strong flows along and at the edges of the elongated feature which we believe 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. The QSL reconnection site has the same spectral appearance as the so-called explosive events identified by strong blue- and red-shifted emission, thus answering a long outstanding question about the true nature of this spectral phenomenon. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: accepted for publication in A&A, the abstract is reduced in size, for the movies contact MM

Journal ref: A&A 660, A45 (2022)

Determining the source and eruption dynamics of a stealth CME using NLFFF modelling and MHD simulations

Authors: Stephanie L. Yardley, Paolo Pagano, Duncan H. Mackay, Lisa A. Upton

Abstract: 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. 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… ▽ More 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. 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. 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. 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 A 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. 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. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: Accepted in A&A, 17 pages, 12 figures

Journal ref: A&A 652, A160 (2021)

Simulating the Coronal Evolution of Bipolar Active Regions to Investigate the Formation of Flux Ropes

Authors: Stephanie L. Yardley, Duncan H. Mackay, Lucie M. Green

Abstract: 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 et al. A time sequence of cleaned photospheric line-of-sight magnetograms, that covers the entire evolution of each AR, is used to drive the simulation. A comparison of the simulated coronal magnetic field with the 1… ▽ More 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 et al. A time sequence of cleaned photospheric line-of-sight magnetograms, that 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 A 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. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 32 pages, 8 figures, accepted for publication

Identifying solar-like magnetic cycles with Zeeman-Doppler-Imaging (ZDI)

Authors: Lisa T. Lehmann, Gaitee A. J. Hussain, Aline A. Vidotto, Moira M. Jardine, Duncan H. Mackay

Abstract: We are reaching the point where spectropolarimetric surveys have run for long enough to reveal solar-like magnetic activity cycles. In this paper we investigate what would be the best strategy to identify solar-like magnetic cycles and ask which large-scale magnetic field parameters best follow a solar-type magnetic cycle and are observable with the Zeeman-Doppler-Imaging (ZDI) technique. We appro… ▽ More We are reaching the point where spectropolarimetric surveys have run for long enough to reveal solar-like magnetic activity cycles. In this paper we investigate what would be the best strategy to identify solar-like magnetic cycles and ask which large-scale magnetic field parameters best follow a solar-type magnetic cycle and are observable with the Zeeman-Doppler-Imaging (ZDI) technique. We approach these questions using the 3D non-potential flux transport simulations of \cite{Yeates2012} modelling the solar vector magnetic field over 15 years (centred on solar cycle 23). The flux emergence profile was extracted from solar synoptic maps and used as input for a photospheric flux transport model in combination with a non-potential coronal evolution model. We synthesise spectropolarimetric data from the simulated maps and reconstruct them using ZDI. The ZDI observed solar cycle is set into the context of other cool star observations and we present observable trends of the magnetic field topology with time, sunspot number and S-index. We find that the axisymmetric energy fraction is the best parameter of the ZDI detectable large-scale field to trace solar-like cycles. Neither the surface averaged large-scale field or the total magnetic energy is appropriate. ZDI seems also to be able to recover the increase of the toroidal energy with S-index. We see further that ZDI might unveil hints of the dynamo modes that are operating and of the global properties of the small-scale flux emergence like active latitudes. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: 21 pages, 19 figures, accepted for publication in MNRAS

Eruptions from coronal hole bright points: observations and non-potential modelling

Authors: Maria S. Madjarska, Klaus Galsgaard, Duncan H. Mackay, Kostadinka Koleva, Momchil Dechev

Abstract: A single case study of a CBP in an equatorial coronal hole with an exceptionally large size is investigated to extend 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 the extreme-ultraviolet (EUV) and X-ray emission. We aim to explore the nature of the so-called micro-flares in CBPs… ▽ More A single case study of a CBP in an equatorial coronal hole with an exceptionally large size is investigated to extend 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 the 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. Co-observations from the Atmospheric Imaging Assembly (AIA) and Helioseismic Magnetic Imager (HMI) on board the Solar Dynamics Observatory, and GONG Halpha 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. A mini-filament (MF) that formed beneath the CBP arcade around 3-4 h before the eruption is seen in the Halpha 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 the MF liftoff. The CBP micro-flare is associated with three flare kernels that formed shortly after the MF liftoff. No observational signature is found for 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. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 13 pages, 11 figures, the abstract has been reduced from its original length, accepted by A&A

Journal ref: A&A 643, A19 (2020)

Investigation of the Middle Corona with SWAP and a Data-Driven Non-Potential Coronal Magnetic Field Model

Authors: Karen A. Meyer, Duncan H. Mackay, Dana-Camelia Talpeanu, Lisa A. Upton, Matthew J. West

Abstract: The large field-of-view of the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) instrument on board the PRoject for Onboard Autonomy 2 (PROBA2) spacecraft provides a unique opportunity to study extended coronal structures observed in EUV in conjunction with global coronal magnetic field simulations. A global non-potential magnetic field model is used to simulate the evolu… ▽ More The large field-of-view of the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) instrument on board the PRoject for Onboard Autonomy 2 (PROBA2) spacecraft provides a unique opportunity to study extended coronal structures observed in 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. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 26 pages, 11 figures

Hydrogen Non-Equilibrium Ionisation Effects in Coronal Mass Ejections

Authors: P. Pagano, A. Bemporad, D. H. Mackay

Abstract: A new generation of coronagraphs to study the solar wind and CMEs are being developed and launched. These coronagraphs will heavily rely on multi-channel observations where visible light (VL) and UV-EUV observations provide new plasma diagnostics. One of these instruments, Metis on board ESA-Solar Orbiter, will simultaneously observe VL and the UV Lyman-$α$ line. The number of neutral Hydrogen ato… ▽ More A new generation of coronagraphs to study the solar wind and CMEs are being developed and launched. These coronagraphs will heavily rely on multi-channel observations where visible light (VL) and UV-EUV observations provide new plasma diagnostics. One of these instruments, Metis on board ESA-Solar Orbiter, will simultaneously observe VL and the UV Lyman-$α$ line. The number of neutral Hydrogen atoms (a small fraction of coronal protons) is a key parameter for deriving plasma properties such as temperature from the observed Lyman-$α$ line intensity. However, these measurements are significantly affected if non-equilibrium ionisation effects occur. The aim of this work is to determine if non-equilibrium ionisation effects are relevant in CMEs and in particular when and in which regions of the CME. We use a magneto-hydrodynamic simulation of a magnetic flux rope ejection to generate a CME. From this we then reconstruct the ionisation state of Hydrogen atoms in the CME by evaluating both the advection of neutral and ionised Hydrogen atoms and the ionisation and recombination rates in the MHD simulation. We find that the equilibrium ionisation assumption holds mostly in the core of the CME. In contrast non-equilibrium ionisation effects are significant at the CME front, where we find about 100 times more neutral Hydrogen atoms than prescribed by ionisation equilibrium conditions, even if this neutral Hydrogen excess might be difficult to identify due to projection effects. This work provides key information for the development of a new generation of diagnostic techniques that aim at combining visible light and Lyman-$α$ line emissions. The results show that non-ionisation equilibrium effects need to be considered when we analyse CME fronts. To incorrectly assume equilibrium ionisation in these regions would lead to a systematic underestimate of plasma temperatures. △ Less

Submitted 27 March, 2020; originally announced March 2020.

Measuring stellar magnetic helicity density

Authors: K. Lund, M. Jardine, L. T. Lehmann, D. H. Mackay, V. See, A. A. Vidotto, J. -F. Donati, R. Fares, C. P. Folsom, S. V. Jeffers, S. C. Marsden, J. Morin, P. Petit

Abstract: Helicity is a fundamental property of a magnetic field but to date it has only been possible to observe its evolution in one star - the Sun. In this paper we provide a simple technique for mapping the large-scale helicity density across the surface of any star using only observable quantities: the poloidal and toroidal magnetic field components (which can be determined from Zeeman-Doppler imaging)… ▽ More Helicity is a fundamental property of a magnetic field but to date it has only been possible to observe its evolution in one star - the Sun. In this paper we provide a simple technique for mapping the large-scale helicity density across the surface of any star using only observable quantities: the poloidal and toroidal magnetic field components (which can be determined from Zeeman-Doppler imaging) and the stellar radius. We use a sample of 51 stars across a mass range of 0.1-1.34 M$_\odot$ to show how the helicity density relates to stellar mass, Rossby number, magnetic energy and age. We find that the large-scale helicity density increases with decreasing Rossby number $R_o$, peaking at $R_o \simeq 0.1$, with a saturation or decrease below that. For both fully- and partially-convective stars we find that the mean absolute helicity density scales with the mean squared toroidal magnetic flux density according to the power law: $|\langle{h\,}\rangle|$ $\propto$ $\langle{\rm{B_{tor}}^2_{}\,\rangle}^{0.86\,\pm\,0.04}$. The scatter in this relation is consistent with the variation across a solar cycle, which we compute using simulations and observations across solar cycles 23 and 24 respectively. We find a significant decrease in helicity density with age. △ Less

Submitted 31 January, 2020; originally announced January 2020.

Comments: 11 pages, 7 figures

Understanding the plasma and magnetic field evolution of a filament using observations and Nonlinear force-free field modelling

Authors: Stephanie L. Yardley, Antonia Savcheva, Lucie M. Green, Lidia van Driel-Gesztelyi, David Long, David R. Williams, Duncan H. Mackay

Abstract: We present observations and magnetic field models of an intermediate filament present on the Sun in August 2012, associated with a polarity inversion line that extends from AR 11541 in the east into the quiet sun at its western end. A combination of SDO/AIA, SDO/HMI, and GONG H alpha data allow us to analyse the structure and evolution of the filament from 2012 August 4 23:00 UT to 2012 August 6 0… ▽ More We present observations and magnetic field models of an intermediate filament present on the Sun in August 2012, associated with a polarity inversion line that extends from AR 11541 in the east into the quiet sun at its western end. A combination of SDO/AIA, SDO/HMI, and GONG H alpha data allow us to analyse the structure and evolution of the filament from 2012 August 4 23:00 UT to 2012 August 6 08:00 UT when the filament was in equilibrium. By applying the flux rope insertion method, nonlinear force-free field models of the filament are constructed using SDO/HMI line-of-sight magnetograms as the boundary condition at the two times given above. Guided by observed filament barbs, both modelled flux ropes are split into three sections each with a different value of axial flux to represent the non-uniform photospheric field distribution. The flux in the eastern section of the rope increases by 4$\times$10$^{20}$ Mx between the two models, which is in good agreement with the amount of flux cancelled along the internal PIL of AR 11541, calculated to be 3.2$\times$10$^{20}$ Mx. This suggests that flux cancellation builds flux into the filament's magnetic structure. Additionally, the number of field line dips increases between the two models in the locations where flux cancellation, the formation of new filament threads and growth of the filament is observed. This suggests that flux cancellation associated with magnetic reconnection forms concave-up magnetic field that lifts plasma into the filament. During this time, the free magnetic energy in the models increases by 0.2$\times$10$^{31}$ ergs. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 17 pages, 10 figures

A New Space Weather Tool for Identifying Eruptive Active Regions

Authors: P. Pagano, D. H. Mackay, S. L. Yardley

Abstract: One of the main goals of solar physics is the timely identification of eruptive active regions. Space missions such as Solar Orbiter or future Space Weather forecasting missions would largely benefit from this achievement. Our aim is to produce a relatively simple technique that can provide real time indications or predictions that an active region will produce an eruption. We expand on the theore… ▽ More One of the main goals of solar physics is the timely identification of eruptive active regions. Space missions such as Solar Orbiter or future Space Weather forecasting missions would largely benefit from this achievement. Our aim is to produce a relatively simple technique that can provide real time indications or predictions that an active region will produce an eruption. We expand on the theoretical work of \citet{Pagano2019fp} that was able to distinguish eruptive from non-eruptive active regions. From this we introduce a new operational metric that uses a combination of observed line-of-sight magnetograms, 3D data-driven simulations and the projection of the 3D simulations forward in time. Results show that the new metric correctly distinguishes active regions as eruptive when observable signatures of eruption have been identified and as non-eruptive when there are no observable signatures of eruption. After successfully distinguishing eruptive from non-eruptive active regions we illustrate how this metric may be used in a "real-time" operational sense were three levels of warning are categorised. These categories are: high risk (red), medium risk (amber) and low risk (green) of eruption. Through considering individual cases we find that the separation into eruptive and non-eruptive active regions is more robust the longer the time series of observed magnetograms used to simulate the build up of magnetic stress and free magnetic energy within the active region. Finally, we conclude that this proof of concept study delivers promising results where the ability to categorise the risk of an eruption is a major achievement. △ Less

Submitted 9 October, 2019; originally announced October 2019.

A Prospective New Diagnostic Technique for Distinguishing Eruptive and Non-Eruptive Active Regions

Authors: P. Pagano, D. H. Mackay, S. L. Yardley

Abstract: Active regions are the source of the majority of magnetic flux rope ejections that become Coronal Mass Ejections (CMEs). To identify in advance which active regions will produce an ejection is key for both space weather prediction tools and future science missions such as Solar Orbiter. The aim of this study is to develop a new technique to identify which active regions are more likely to generate… ▽ More Active regions are the source of the majority of magnetic flux rope ejections that become Coronal Mass Ejections (CMEs). To identify in advance which active regions will produce an ejection is key for both space weather prediction tools and future science missions such as Solar Orbiter. The aim of this study is to develop a new technique to identify which active regions are more likely to generate magnetic flux rope ejections. The new technique will aim to: (i) produce timely space weather warnings and (ii) open the way to a qualified selection of observational targets for space-borne instruments. We use a data-driven Non-linear Force-Free Field (NLFFF) model to describe the 3D evolution of the magnetic field of a set of active regions. We determine a metric to distinguish eruptive from non-eruptive active regions based on the Lorentz force. Furthermore, using a subset of the observed magnetograms, we run a series of simulations to test whether the time evolution of the metric can be predicted. The identified metric successfully differentiates active regions observed to produce eruptions from the non-eruptive ones in our data sample. A meaningful prediction of the metric can be made between 6 to 16 hours in advance. This initial study presents an interesting first step in the prediction of CME onset using only LOS magnetogram observations combined with NLFFF modelling. Future studies will address how to generalise the model such that it can be used in a more operational sense and for a variety of simulation approaches. △ Less

Submitted 24 August, 2019; originally announced August 2019.

Eruptions from Quiet Sun Coronal Bright Points. II. Non-Potential Modeling

Authors: Klaus Galsgaard, Maria S. Madjarska, Duncan H. Mackay, Chaozhou Mou

Abstract: Our recent observational study shows that the majority of coronal bright points (CBPs) in the quiet Sun are sources of one or more eruptions during their lifetime. Here, we investigate the non-potential time dependent structure of the magnetic field of the CBP regions with special emphasis on the time evolving magnetic structure at the spatial locations where the eruptions are initiated. The magne… ▽ More Our recent observational study shows that the majority of coronal bright points (CBPs) in the quiet Sun are sources of one or more eruptions during their lifetime. Here, we investigate the non-potential time dependent structure of the magnetic field of the CBP regions with special emphasis on the time evolving magnetic structure at the spatial locations where the eruptions are initiated. The magnetic structure is evolved in time using a Non-Linear Force Free Field (NLFFF) relaxation approach, based on a time series of Helioseismic and Magnetic Imager (HMI) longitudinal magnetograms. This results in a continuous time series of NLFFFs. The time series is initiated with a potential field extrapolation based on a magnetogram taken well before the time of the eruptions. This initial field is then evolved in time in response to the observed changes in the magnetic field distribution at the photosphere. The local and global magnetic field structures from the time series of NLFFF field solutions are analysed in the vicinity of the eruption sites at the approximate times of the eruptions. The analysis shows that many of the CBP eruptions reported in Mou etl. (2018) contain twisted flux tube located at the sites of eruptions. The presence of flux ropes at these locations provides in many cases a direct link between the magnetic field structure, their eruption and the observation of mini coronal mass ejections (mini-CMEs). It is found that all repetitive eruptions are homologous. The NLFFF simulations show that twisted magnetic field structures are created at the locations hosting eruptions in CBPs. These twisted structures are produced by footpoint motions imposed by changes in the photospheric magnetic field observations. The true nature of the micro-flares remains unknown. Further 3D data-driven MHD modelling is required to show how these twisted regions become unstable and erupt. △ Less

Submitted 28 January, 2019; originally announced January 2019.

Comments: 16 pages, 16 figures, in press in A&A

Journal ref: A&A 623, A78 (2019)

Observing the simulations: Applying ZDI to 3D non-potential magnetic field simulations

Authors: L. T. Lehmann, G. A. J. Hussain, M. M. Jardine, D. H. Mackay, A. A. Vidotto

Abstract: The large-scale magnetic fields of stars can be obtained with the Zeeman-Doppler-Imaging (ZDI) technique, but their interpretation is still challenging as the contribution of the small-scale field or the reliability of the reconstructed field properties is still not fully understood. To quantify this, we use 3D non-potential magnetic field simulations for slowly rotating solar-like stars as inputs… ▽ More The large-scale magnetic fields of stars can be obtained with the Zeeman-Doppler-Imaging (ZDI) technique, but their interpretation is still challenging as the contribution of the small-scale field or the reliability of the reconstructed field properties is still not fully understood. To quantify this, we use 3D non-potential magnetic field simulations for slowly rotating solar-like stars as inputs to test the capabilities of ZDI. These simulations are based on a flux transport model connected to a non-potential coronal evolution model using the observed solar flux emergence pattern. We first compare four field prescriptions regarding their reconstruction capabilities and investigate the influence of the spatial resolution of the input maps on the corresponding circularly polarised profiles. We then generate circularly polarised spectra based on our high resolution simulations of three stellar models with different activity levels, and reconstruct their large-scale magnetic fields using a non-potential ZDI code assuming two different stellar inclination angles. Our results show that the ZDI technique reconstructs the main features of slowly rotating solar-like stars but with $\sim\,$one order of magnitude less magnetic energy. The large-scale field morphologies are recovered up to harmonic modes $\ell \sim 5$, especially after averaging over several maps for each stellar model. While ZDI is not able to reproduce the input magnetic energy distributions across individual harmonic modes, the fractional energies across the modes are generally within $20\,\%$ agreement. The fraction of axisymmetric and toroidal field tends to be overestimated for stars with solar flux emergence patterns for more pole-on inclination angles. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 26 pages, 16 figures, accepted for publication in MNRAS

The role of flux cancellation in eruptions from bipolar active regions

Authors: S. L. Yardley, L. M. Green, L. Van Driel-Gesztelyi, D. R. Williams, D. H. Mackay

Abstract: The physical processes or trigger mechanisms that lead to the eruption of coronal mass ejections (CMEs), the largest eruptive phenomenon in the heliosphere, are still undetermined. Low-altitude magnetic reconnection associated with flux cancellation appears to play an important role in CME occurrence as it can form an eruptive configuration and reduce the magnetic flux that contributes to the over… ▽ More The physical processes or trigger mechanisms that lead to the eruption of coronal mass ejections (CMEs), the largest eruptive phenomenon in the heliosphere, are still undetermined. Low-altitude magnetic reconnection associated with flux cancellation appears to play an important role in CME occurrence as it can form an eruptive configuration and reduce the magnetic flux that contributes to the overlying, stabilising field. We conduct the first comprehensive study of 20 small bipolar active regions in order to probe the role of flux cancellation as an eruption trigger mechanism. We categorise eruptions from the bipolar regions into three types related to location and find that the type of eruption produced depends on the evolutionary stage of the active region. In addition we find that active regions that form eruptive structures by flux cancellation (low-altitude reconnection) had, on average, lower flux cancellation rates than the active region sample as a whole. Therefore, while flux cancellation plays a key role, by itself it is insufficient for the production of an eruption. The results support that although flux cancellation in a sheared arcade may be able to build an eruptive configuration, a successful eruption depends upon the removal of sufficient overlying and stabilising field. Convergence of the bipole polarities also appears to be present in regions that produce an eruption. These findings have important implications for understanding the physical processes that occur on our Sun in relation to CMEs and for space weather forecasting. △ Less

Submitted 31 August, 2018; originally announced August 2018.

Comments: 16 pages, 4 figures

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

Authors: A. R. Yeates, T. Amari, I. Contopoulos, X. Feng, D. H. Mackay, Z. Mikić, T. Wiegelmann, J. Hutton, C. A. Lowder, H. Morgan, G. Petrie, L. A. Rachmeler, L. A. Upton, A. Canou, P. Chopin, C. Downs, M. Druckmüller, J. A. Linker, D. B. Seaton, T. Török

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

Submitted 2 August, 2018; originally announced August 2018.

Comments: 29 pages, 11 figures, accepted for publication in Space Science Reviews

Connecting the large- and the small-scale magnetic fields of solar-like stars

Authors: L. T. Lehmann, M. M. Jardine, D. H. Mackay, A. A. Vidotto

Abstract: A key question in understanding the observed magnetic field topologies of cool stars is the link between the small- and the large-scale magnetic field and the influence of the stellar parameters on the magnetic field topology. We examine various simulated stars to connect the small-scale with the observable large-scale field. The highly resolved 3D simulations we used couple a flux transport model… ▽ More A key question in understanding the observed magnetic field topologies of cool stars is the link between the small- and the large-scale magnetic field and the influence of the stellar parameters on the magnetic field topology. We examine various simulated stars to connect the small-scale with the observable large-scale field. The highly resolved 3D simulations we used couple a flux transport model with a non-potential coronal model using a magnetofrictional technique. The surface magnetic field of these simulations is decomposed into spherical harmonics which enables us to analyse the magnetic field topologies on a wide range of length scales and to filter the large-scale magnetic field for a direct comparison with the observations. We show that the large-scale field of the self-consistent simulations fits the observed solar-like stars and is mainly set up by the global dipolar field and the large-scale properties of the flux pattern, e.g. the averaged latitudinal position of the emerging small-scale field and its global polarity pattern. The stellar parameters flux emergence rate, differential rotation and meridional flow affect the large-scale magnetic field topology. An increased flux emergence rate increases the magnetic flux in all field components and an increased differential rotation increases the toroidal field fraction by decreasing the poloidal field. The meridional flow affects the distribution of the magnetic energy across the spherical harmonic modes. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Comments: 20 pages, 16 figures, accepted for publication in MNRAS

Photospheric observations of surface and body modes in solar magnetic pores

Authors: Peter H. Keys, Richard J. Morton, David B. Jess, Gary Verth, Samuel D. T. Grant, Mihalis Mathioudakis, Duncan H. Mackay, John G. Doyle, Damian J. Christian, Francis P. Keenan, Robertus Erdelyi

Abstract: Over the past number of years, great strides have been made in identifying the various low-order magnetohydrodynamic wave modes observable in a number of magnetic structures found within the solar atmosphere. However, one aspect of these modes that has remained elusive, until now, is their designation as either surface or body modes. This property has significant implications on how these modes tr… ▽ More Over the past number of years, great strides have been made in identifying the various low-order magnetohydrodynamic wave modes observable in a number of magnetic structures found within the solar atmosphere. However, one aspect of these modes that has remained elusive, until now, is their designation as either surface or body modes. This property has significant implications on how these modes transfer energy from the waveguide to the surrounding plasma. Here, for the first time to our knowledge, we present conclusive, direct evidence of these wave characteristics in numerous pores which were observed to support sausage modes. As well as outlining methods to detect these modes in observations, we make estimates of the energies associated with each mode. We find surface modes more frequently in the data, and also that surface modes appear to carry more energy than those displaying signatures of body modes. We find frequencies in the range of ~2 to 12 mHz with body modes as high as 11 mHz, but we do not find surface modes above 10 mHz. It is expected that the techniques we have applied will help researchers search for surface and body signatures in other modes and in differing structures to those presented here. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 20 pages, 9 figures, Accepted for publication in ApJ

A new technique for observationally derived boundary conditions for space weather

Authors: P. Pagano, D. H. Mackay, A. R. Yeates

Abstract: In recent years, space weather research has focused on developing modelling techniques to predict the arrival time and properties of coronal mass ejections (CMEs) at the Earth. The aim of this paper is to propose a new modelling technique suitable for the next generation of Space Weather predictive tools that is both efficient and accurate. The aim of the new approach is to provide interplanetary… ▽ More In recent years, space weather research has focused on developing modelling techniques to predict the arrival time and properties of coronal mass ejections (CMEs) at the Earth. The aim of this paper is to propose a new modelling technique suitable for the next generation of Space Weather predictive tools that is both efficient and accurate. The aim of the new approach is to provide interplanetary space weather forecasting models with accurate time dependent boundary conditions of erupting magnetic flux ropes in the upper solar corona. To produce boundary conditions, we couple two different modelling techniques, MHD simulations and a quasi-static non-potential evolution model. Both are applied on a spatial domain that covers the entire solar surface. The non-potential model uses a time series of observed synoptic magnetograms to drive the non-potential quasi-static evolution of the coronal magnetic field. This allows us to follow the formation and loss of equilibrium of magnetic flux ropes. Following this a MHD simulation captures the dynamic evolution of the erupting flux rope. The present paper focuses on the MHD simulations that follow the ejection of magnetic flux ropes to 4$R_\odot$. We first propose a technique for specifying the pre-eruptive plasma properties in the corona. Next, time dependent MHD simulations describe the ejection of two magnetic flux ropes, that produce time dependent boundary conditions for the magnetic field and plasma at 4$R_{\odot}$. In the present paper, we show that the dual use of quasi-static non-potential magnetic field simulations and full time dependent MHD simulations can produce realistic inhomogeneous boundary conditions for space weather forecasting tools. Before a fully operational model can be produced there are a number of technical and scientific challenges that still need to be addressed. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Simulating the Coronal Evolution of AR 11437 using SDO/HMI Magnetograms

Authors: Stephanie L. Yardley, Duncan H. Mackay, Lucie M. Green

Abstract: The coronal magnetic field evolution of AR 11437 is simulated by applying the magnetofrictional relaxation technique of Mackay et al. (2011). A sequence of photospheric line-of-sight magnetograms produced by SDO/HMI are used to drive the simulation and continuously evolve the coronal magnetic field of the active region through a series of non-linear force-free equilibria. The simulation is started… ▽ More The coronal magnetic field evolution of AR 11437 is simulated by applying the magnetofrictional relaxation technique of Mackay et al. (2011). A sequence of photospheric line-of-sight magnetograms produced by SDO/HMI are used to drive the simulation and continuously evolve the coronal magnetic field of the active region through a series of non-linear force-free equilibria. The simulation is started during the first stages of the active region emergence so that its full evolution from emergence to decay can be simulated. A comparison of the simulation results with SDO/AIA observations show that many aspects of the active region's observed coronal evolution are reproduced. In particular, it shows the presence of a flux rope, which forms at the same location as sheared coronal loops in the observations. The observations show that eruptions occur on 2012 March 17 at 05:09 UT and 10:45 UT and on 2012 March 20 at 14:31 UT. The simulation reproduces the first and third eruption, with the simulated flux rope erupting roughly 1 and 10 hours before the observed ejections, respectively. A parameter study is conducted where the boundary and initial conditions are varied along with the physical effects of Ohmic diffusion, hyperdiffusion and an additional injection of helicity. When comparing the simulations, the evolution of the magnetic field, free magnetic energy, relative helicity and flux rope eruption timings do not change significantly. This indicates that the key element in reproducing the coronal evolution of AR 11437 is the use of line-of-sight magnetograms to drive the evolution of the coronal magnetic field. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Comments: 18 pages, 13 figures, accepted in ApJ

Observations and Modelling of the Pre-Flare Period of the 29 March 2014 X1 Flare

Authors: M. M. Woods, L. K. Harra, S. A. Matthews, D. H. Mackay, S. Dacie, D. M. Long

Abstract: On the 29 March 2014 NOAA active region (AR) 12017 produced an X1 flare which was simultaneously observed by an unprecedented number of observatories. We have investigated the pre-flare period of this flare from 14:00 UT until 19:00 UT using joint observations made by the Interface Region Imaging Spectrometer (IRIS) and the Hinode Extreme Ultraviolet Imaging Spectrometer (EIS). Spectral lines prov… ▽ More On the 29 March 2014 NOAA active region (AR) 12017 produced an X1 flare which was simultaneously observed by an unprecedented number of observatories. We have investigated the pre-flare period of this flare from 14:00 UT until 19:00 UT using joint observations made by the Interface Region Imaging Spectrometer (IRIS) and the Hinode Extreme Ultraviolet Imaging Spectrometer (EIS). Spectral lines providing coverage of the solar atmosphere from chromosphere to the corona were analysed to investigate pre-flare activity within the AR. The results of the investigation have revealed evidence of strongly blue-shifted plasma flows, with velocities up to 200 km/s, being observed 40 minutes prior to flaring. These flows are located along the filament present in the active region and are both spatially discrete and transient. In order to constrain the possible explanations for this activity, we undertake non-potential magnetic field modelling of the active region. This modelling indicates the existence of a weakly twisted flux rope along the polarity inversion line in the region where a filament and the strong pre-flare flows are observed. We then discuss how these observations relate to the current models of flare triggering. We conclude that the most likely drivers of the observed activity are internal reconnection in the flux rope, early onset of the flare reconnection, or tether cutting reconnection along the filament. △ Less

Submitted 23 January, 2017; originally announced January 2017.

Comments: 27 pages, 15 figures, accepted for publication in Solar Physics

The energy budget of stellar magnetic fields: comparing non-potential simulations and observations

Authors: L. T. Lehmann, M. M. Jardine, A. A. Vidotto, D. H. Mackay, V. See, J. -F. Donati, C. P. Folsom, S. V. Jeffers, S. C. Marsden, J. Morin, P. Petit

Abstract: The magnetic geometry of the surface magnetic fields of more than 55 cool stars have now been mapped using spectropolarimetry. In order to better understand these observations, we compare the magnetic field topology at different surface scale sizes of observed and simulated cool stars. For ease of comparison between the high-resolution non-potential magnetofrictional simulations and the relatively… ▽ More The magnetic geometry of the surface magnetic fields of more than 55 cool stars have now been mapped using spectropolarimetry. In order to better understand these observations, we compare the magnetic field topology at different surface scale sizes of observed and simulated cool stars. For ease of comparison between the high-resolution non-potential magnetofrictional simulations and the relatively low-resolution observations, we filter out the small-scale field in the simulations using a spherical harmonics decomposition. We show that the large-scale field topologies of the solar-based simulations produce values of poloidal/toroidal fields and fractions of energy in axisymmetric modes that are similar to the observations. These global non-potential evolution model simulations capture key magnetic features of the observed solar-like stars through the processes of surface flux transport and magnetic flux emergence. They do not, however, reproduce the magnetic field of M-dwarfs or stars with dominantly toroidal field. Furthermore, we analyse the magnetic field topologies of individual spherical harmonics for the simulations and discover that the dipole is predominately poloidal, while the quadrupole shows the highest fraction of toroidal fields. Magnetic field structures smaller than a quadrupole display a fixed ratio between the poloidal and toroidal magnetic energies. △ Less

Submitted 26 October, 2016; originally announced October 2016.

Comments: 6 pages, 3 figures, accepted for publication in MNRAS Letters

Solar coronal magnetic fields derived using seismology techniques applied to omnipresent sunspot waves

Authors: D. B. Jess, V. E. Reznikova, R. S. I. Ryans, D. J. Christian, P. H. Keys, M. Mathioudakis, D. H. Mackay, S. Krishna Prasad, D. Banerjee, S. D. T. Grant, S. Yau, C. Diamond

Abstract: Sunspots on the surface of the Sun are the observational signatures of intense manifestations of tightly packed magnetic field lines, with near-vertical field strengths exceeding 6,000 G in extreme cases. It is well accepted that both the plasma density and the magnitude of the magnetic field strength decrease rapidly away from the solar surface, making high-cadence coronal measurements through tr… ▽ More Sunspots on the surface of the Sun are the observational signatures of intense manifestations of tightly packed magnetic field lines, with near-vertical field strengths exceeding 6,000 G in extreme cases. It is well accepted that both the plasma density and the magnitude of the magnetic field strength decrease rapidly away from the solar surface, making high-cadence coronal measurements through traditional Zeeman and Hanle effects difficult since the observational signatures are fraught with low-amplitude signals that can become swamped with instrumental noise. Magneto-hydrodynamic (MHD) techniques have previously been applied to coronal structures, with single and spatially isolated magnetic field strengths estimated as 9-55 G. A drawback with previous MHD approaches is that they rely on particular wave modes alongside the detectability of harmonic overtones. Here we show, for the first time, how omnipresent magneto-acoustic waves, originating from within the underlying sunspot and propagating radially outwards, allow the spatial variation of the local coronal magnetic field to be mapped with high precision. We find coronal magnetic field strengths of 32 +/- 5 G above the sunspot, which decrease rapidly to values of approximately 1 G over a lateral distance of 7000 km, consistent with previous isolated and unresolved estimations. Our results demonstrate a new, powerful technique that harnesses the omnipresent nature of sunspot oscillations to provide magnetic field mapping capabilities close to a magnetic source in the solar corona. △ Less

Submitted 19 May, 2016; originally announced May 2016.

Comments: Published in Nature Physics, 12, 179 (2016; DOI:10.1038/nphys3544). 22 pages (Letter + Supplementary Material), including 7 figures

Journal ref: Nature Physics, 12, 179 (2016)

Stellar Coronal Response to Differential Rotation and Flux Emergence

Authors: G. P. S. Gibb, D. H. Mackay, M. M. Jardine, A. R. Yeates

Abstract: We perform a numerical parameter study to determine what effect varying differential rotation and flux emergence has on a star's non-potential coronal magnetic field. In particular we consider the effects on the star's surface magnetic flux, open magnetic flux, mean azimuthal field strength, coronal free magnetic energy, coronal heating and flux rope eruptions. To do this, we apply a magnetic flux… ▽ More We perform a numerical parameter study to determine what effect varying differential rotation and flux emergence has on a star's non-potential coronal magnetic field. In particular we consider the effects on the star's surface magnetic flux, open magnetic flux, mean azimuthal field strength, coronal free magnetic energy, coronal heating and flux rope eruptions. To do this, we apply a magnetic flux transport model to describe the photospheric evolution, and couple this to the non-potential coronal evolution using a magnetofrictional technique. A flux emergence model is applied to add new magnetic flux onto the photosphere and into the corona. The parameters of this flux emergence model are derived from the solar flux emergence profile, however the rate of emergence can be increased to represent higher flux emergence rates than the Sun's. Overall we find that flux emergence has a greater effect on the non-potential coronal properties compared to differential rotation, with all the aforementioned properties increasing with increasing flux emergence rate. Although differential rotation has a lesser effect on the overall coronal properties compared to flux emergence, varying differential rotation does alter the coronal structure. As the differential rotation rate increases, the corona becomes more open, and more non-potential. △ Less

Submitted 10 March, 2016; originally announced March 2016.

Comments: 15 pages, 18 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 456, Issue 4, p.3624-3637, 03/2016

Influence of Non-Potential Coronal Magnetic Topology on Solar-Wind Models

Authors: S. J. Edwards, A. R. Yeates, F. -X. Bocquet, D. H. Mackay

Abstract: By comparing a magneto-frictional model of the low coronal magnetic field to a potential-field source-surface model, we investigate the possible impact of non-potential magnetic structure on empirical solar-wind models. These empirical models (such as Wang-Sheeley-Arge) estimate the distribution of solar-wind speed solely from the magnetic-field structure in the low corona. Our models are computed… ▽ More By comparing a magneto-frictional model of the low coronal magnetic field to a potential-field source-surface model, we investigate the possible impact of non-potential magnetic structure on empirical solar-wind models. These empirical models (such as Wang-Sheeley-Arge) estimate the distribution of solar-wind speed solely from the magnetic-field structure in the low corona. Our models are computed in a domain between the solar surface and 2.5 solar radii, and are extended to 0.1 AU using a Schatten current-sheet model. The non-potential field has a more complex magnetic skeleton and quasi-separatrix structures than the potential field, leading to different sub-structure in the solar-wind speed proxies. It contains twisted magnetic structures that can perturb the separatrix surfaces traced down from the base of the heliospheric current sheet. A significant difference between the models is the greater amount of open magnetic flux in the non-potential model. Using existing empirical formulae this leads to higher predicted wind speeds for two reasons: partly because magnetic flux tubes expand less rapidly with height, but more importantly because more open field lines are further from coronal-hole boundaries. △ Less

Submitted 2 November, 2015; originally announced November 2015.

Comments: 21 pages, 14 figures, accepted for publication in Solar Physics

Simulating AIA observations of a flux rope ejection

Authors: Paolo Pagano, Duncan H. Mackay, Stefaan Poedts

Abstract: Extreme ultraviolet (EUV) images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) are providing new insights into the early phase of CME evolution. Observations now show the ejection of magnetic flux ropes from the solar corona and how they evolve into CMEs. These observations are difficult to interpret in terms of basic physical mechanisms and quantities. T… ▽ More Extreme ultraviolet (EUV) images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) are providing new insights into the early phase of CME evolution. Observations now show the ejection of magnetic flux ropes from the solar corona and how they evolve into CMEs. These observations are difficult to interpret in terms of basic physical mechanisms and quantities. To fully understand CMEs we need to compare equivalent quantities derived from both observations and theoretical models. To this end we aim to produce synthesised AIA observations from simluations of a flux rope ejection. To carry this out we include the role of thermal conduction and radiative losses, both of which are important for determining the temperature distribution of the solar corona during a CME. We perform a simulation where a flux rope is ejected from the solar corona. From the density and temperature of the plasma in the simulation we synthesise AIA observations. The emission is then integrated along the line of sight using the instrumental response function of AIA. We sythesise observations of AIA in the channels at 304 A, 171 A, 335 A, and 94 A. The synthesised observations show a number of features similar to actual observations and in particular reproduce the general development of CMEs in the low corona as observed by AIA. In particular we reproduce an erupting and expanding arcade in the 304 A and 171 A channels with a high density core. The ejection of a flux rope reproduces many of the features found in the AIA observations. This work is therefore a step forward in bridging the gap between observations and models, and can lead to more direct interpretations of EUV observations in terms of flux rope ejections. We plan to improve the model in future studies in order to perform a more quantitative comparison. △ Less

Submitted 31 July, 2014; originally announced July 2014.

Journal ref: A&A 568, A120 (2014)

Stellar Differential Rotation and Coronal Timescales

Authors: G. P. S. Gibb, M. M. Jardine, D. H. Mackay

Abstract: We investigate the timescales of evolution of stellar coronae in response to surface differential rotation and diffusion. To quantify this we study both the formation time and lifetime of a magnetic flux rope in a decaying bipolar active region. We apply a magnetic flux transport model to prescribe the evolution of the stellar photospheric field, and use this to drive the evolution of the coronal… ▽ More We investigate the timescales of evolution of stellar coronae in response to surface differential rotation and diffusion. To quantify this we study both the formation time and lifetime of a magnetic flux rope in a decaying bipolar active region. We apply a magnetic flux transport model to prescribe the evolution of the stellar photospheric field, and use this to drive the evolution of the coronal magnetic field via a magnetofrictional technique. Increasing the differential rotation (i.e. decreasing the equator-pole lap time) decreases the flux rope formation time. We find that the formation time is dependent upon the geometric mean of the lap time and the surface diffusion timescale. In contrast, the lifetime of flux ropes are proportional to the lap time. With this, flux ropes on stars with a differential rotation of more than eight times the solar value have a lifetime of less than two days. As a consequence, we propose that features such as solar-like quiescent prominences may not be easily observable on such stars, as the lifetimes of the flux ropes which host the cool plasma are very short. We conclude that such high differential rotation stars may have very dynamical coronae. △ Less

Submitted 12 July, 2014; originally announced July 2014.

Journal ref: MNRAS (October 1, 2014) 443 (4): 3251-3259

The Solar Cycle Variation of Topological Structures in the Global Solar Corona

Authors: Sarah J. Platten, Clare E. Parnell, Andrew L. Haynes, Eric R. Priest, Duncan H. Mackay

Abstract: The distribution of magnetic flux across the solar photosphere results in a complex web of coronal magnetic field structures. To understand this complexity, the magnetic skeleton of the coronal field can be calculated. The skeleton highlights the separatrix surfaces that divide the field into topologically distinct regions, allowing open-field regions on the solar surface to be located. Furthermor… ▽ More The distribution of magnetic flux across the solar photosphere results in a complex web of coronal magnetic field structures. To understand this complexity, the magnetic skeleton of the coronal field can be calculated. The skeleton highlights the separatrix surfaces that divide the field into topologically distinct regions, allowing open-field regions on the solar surface to be located. Furthermore, separatrix surfaces and their intersections (separators) are important likely energy release sites. This paper investigates, throughout the solar cycle, the nature of coronal magnetic-field topologies that arise under the potential-field source-surface approximation. In particular, we characterise the typical global fields at solar maximum and minimum. Global magnetic fields are extrapolated from observed Kitt Peak and SOLIS synoptic magnetograms, from Carrington rotations 1645 to 2144, using the potential-field source-surface model. Hence, variations in the coronal skeleton are studied over three solar cycles. Key building blocks which make up magnetic fields are identified and classified according to the nature of their separatrix surfaces. The magnetic skeleton reveals that, at solar maximum, the global coronal field involves a multitude of topological structures at all latitudes. Many open-field regions exist originating anywhere on the photosphere. At solar minimum, the coronal topology is heavily influenced by the solar dipole. A strong dipole results in a simple large-scale structure involving just two large polar open-field regions, but, at short radial distances between plus or minus 60 deg latitude, the small-scale topology is complex. If the solar dipole if weak, as in the recent minimum, then the low-latitude quiet-sun magnetic fields may be globally significant enough to create many disconnected open-field regions at low latitudes, in addition to the two polar open-field regions. △ Less

Submitted 20 June, 2014; originally announced June 2014.

Journal ref: Astronomy & Astrophysics, Volume 565, id.A44, 14 pp. 2014

Dynamic properties of bright points in an active region

Authors: Peter H. Keys, Mihalis Mathioudakis, David B. Jess, Duncan H. Mackay, Francis P. Keenan

Abstract: Context. Bright points (BPs) are small-scale, magnetic features ubiquitous across the solar surface. Previously, we have observed and noted their properties for quiet Sun regions. Here, we determine the dynamic properties of BPs using simultaneous quiet Sun and active region data. Methods. High spatial and temporal resolution G-band observations of active region AR11372 were obtained with the Rapi… ▽ More Context. Bright points (BPs) are small-scale, magnetic features ubiquitous across the solar surface. Previously, we have observed and noted their properties for quiet Sun regions. Here, we determine the dynamic properties of BPs using simultaneous quiet Sun and active region data. Methods. High spatial and temporal resolution G-band observations of active region AR11372 were obtained with the Rapid Oscillations in the Solar Atmosphere instrument at the Dunn Solar Telescope. Three subfields of varying polarity and magnetic flux density were selected with the aid of magnetograms obtained from the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. Bright points within these subfields were subsequently tracked and analysed. Results. It is found that BPs within active regions display attenuated velocity distributions with an average horizontal velocity of ~0.6 km/s, compared to the quiet region which had an average velocity of 0.9 km/s. Active region BPs are also ~21% larger than quiet region BPs and have longer average lifetimes (~132s) than their quiet region counterparts (88 s). No preferential flow directions are observed within the active region subfields. The diffusion index (gamma) is estimated at ~1.2 for the three regions. Conclusions. We confirm that the dynamic properties of BPs arise predominately from convective motions. The presence of stronger field strengths within active regions is the likely reason behind the varying properties observed. We believe that larger amounts of magnetic flux will attenuate BP velocities by a combination of restricting motion within the intergranular lanes and by increasing the number of stagnation points produced by inhibited convection. Larger BPs are found in regions of higher magnetic flux density and we believe that lifetimes increase in active regions as the magnetic flux stabilises the BPs. △ Less

Submitted 15 May, 2014; originally announced May 2014.

Comments: Accepted for publication in A&A, 9 pages, 3 figures

The Influence of the Magnetic Field on Running Penumbral Waves in the Solar Chromosphere

Authors: D. B. Jess, V. E. Reznikova, T. Van Doorsselaere, P. H. Keys, D. H. Mackay

Abstract: We use images of high spatial and temporal resolution, obtained using both ground- and space-based instrumentation, to investigate the role magnetic field inclination angles play in the propagation characteristics of running penumbral waves in the solar chromosphere. Analysis of a near-circular sunspot, close to the center of the solar disk, reveals a smooth rise in oscillatory period as a functio… ▽ More We use images of high spatial and temporal resolution, obtained using both ground- and space-based instrumentation, to investigate the role magnetic field inclination angles play in the propagation characteristics of running penumbral waves in the solar chromosphere. Analysis of a near-circular sunspot, close to the center of the solar disk, reveals a smooth rise in oscillatory period as a function of distance from the umbral barycenter. However, in one directional quadrant, corresponding to the north direction, a pronounced kink in the period-distance diagram is found. Utilizing a combination of the inversion of magnetic Stokes vectors and force-free field extrapolations, we attribute this behaviour to the cut-off frequency imposed by the magnetic field geometry in this location. A rapid, localised inclination of the magnetic field lines in the north direction results in a faster increase in the dominant periodicity due to an accelerated reduction in the cut-off frequency. For the first time we reveal how the spatial distribution of dominant wave periods, obtained with one of the highest resolution solar instruments currently available, directly reflects the magnetic geometry of the underlying sunspot, thus opening up a wealth of possibilities in future magneto-hydrodynamic seismology studies. In addition, the intrinsic relationships we find between the underlying magnetic field geometries connecting the photosphere to the chromosphere, and the characteristics of running penumbral waves observed in the upper chromosphere, directly supports the interpretation that running penumbral wave phenomena are the chromospheric signature of upwardly-propagating magneto-acoustic waves generated in the photosphere. △ Less

Submitted 29 October, 2013; originally announced October 2013.

Comments: Accepted for publication in ApJ (11 pages, 5 figures)

Effect of gravitational stratification on the propagation of a CME

Authors: Paolo Pagano, Duncan H. Mackay, Stefaan Poedts

Abstract: Our aim is to study the role of gravitational stratification on the propagation of CMEs. In particular, we assess how it influences the speed and shape of CMEs and under what conditions the flux rope ejection becomes a CME or when it is quenched. We ran a set of MHD simulations that adopt an eruptive initial magnetic configuration that has already been shown to be suitable for a flux rope ejection… ▽ More Our aim is to study the role of gravitational stratification on the propagation of CMEs. In particular, we assess how it influences the speed and shape of CMEs and under what conditions the flux rope ejection becomes a CME or when it is quenched. We ran a set of MHD simulations that adopt an eruptive initial magnetic configuration that has already been shown to be suitable for a flux rope ejection. We varied the temperature of the backgroud corona and the intensity of the initial magnetic field to tune the gravitational stratification and the amount of ejected magnetic flux. We used an automatic technique to track the expansion and the propagation of the magnetic flux rope in the MHD simulations. From the analysis of the parameter space, we evaluate the role of gravitational stratification on the CME speed and expansion. Our study shows that gravitational stratification plays a significant role in determining whether the flux rope ejection will turn into a full CME or whether the magnetic flux rope will stop in the corona. The CME speed is affected by the background corona where it travels faster when the corona is colder and when the initial magnetic field is more intense. The fastest CME we reproduce in our parameter space travels at 850 km/s. Moreover, the background gravitational stratification plays a role in the side expansion of the CME, and we find that when the background temperature is higher, the resulting shape of the CME is flattened more. Our study shows that although the initiation mechanisms of the CME are purely magnetic, the background coronal plasma plays a key role in the CME propagation, and full MHD models should be applied when one focusses especially on the production of a CME from a flux rope ejection. △ Less

Submitted 25 October, 2013; originally announced October 2013.

Comments: accepted on A&A

Solar Magnetic Carpet III: Coronal Modelling of Synthetic Magnetograms

Authors: K. A. Meyer, D. H. Mackay, A. A. van Ballegooijen, C. E. Parnell

Abstract: This paper is the third in a series of papers working towards the construction of a realistic, evolving, non-linear force-free coronal field model for the solar magnetic carpet. Here, we present preliminary results of 3D time-dependent simulations of the small-scale coronal field of the magnetic carpet. Four simulations are considered, each with the same evolving photospheric boundary condition: a… ▽ More This paper is the third in a series of papers working towards the construction of a realistic, evolving, non-linear force-free coronal field model for the solar magnetic carpet. Here, we present preliminary results of 3D time-dependent simulations of the small-scale coronal field of the magnetic carpet. Four simulations are considered, each with the same evolving photospheric boundary condition: a 48 hr time series of synthetic magnetograms produced from the model of Meyer, Mackay, van Ballegooijen and Parnell, 2011, Solar Phys., 272, 29. Three simulations include a uniform, overlying coronal magnetic field of differing strength, the fourth simulation includes no overlying field. The build-up, storage and dissipation of magnetic energy within the simulations is studied. In particular, we study their dependence upon the evolution of the photospheric magnetic field and the strength of the overlying coronal field. We also consider where energy is stored and dissipated within the coronal field. The free magnetic energy built up is found to be more than sufficient to power small-scale, transient phenomena such as nanoflares and X-ray bright points, with the bulk of the free energy found to be stored low down, between 0.5-0.8 Mm. The energy dissipated is presently found to be too small to account for the heating of the entire quiet Sun corona. However, the form and location of energy dissipation regions are in qualitative agreement with what is observed on small scales on the Sun. Future MHD modelling using the same synthetic magnetograms may lead to a higher energy release. △ Less

Submitted 6 March, 2013; originally announced March 2013.

The Sun's Global Photospheric and Coronal Magnetic Fields: Observations and Models

Authors: D. H. Mackay, A. R. Yeates

Abstract: In this review, our present day understanding of the Sun's global photospheric and coronal magnetic fields is discussed from both observational and theoretical viewpoints. Firstly, the large-scale properties of photospheric magnetic fields are described, along with recent advances in photospheric magnetic flux transport models. Following this, the wide variety of theoretical models used to simulat… ▽ More In this review, our present day understanding of the Sun's global photospheric and coronal magnetic fields is discussed from both observational and theoretical viewpoints. Firstly, the large-scale properties of photospheric magnetic fields are described, along with recent advances in photospheric magnetic flux transport models. Following this, the wide variety of theoretical models used to simulate global coronal magnetic fields are described. From this, the combined application of both magnetic flux transport simulations and coronal modeling techniques to describe the phenomena of coronal holes, the Sun's open magnetic flux and the hemispheric pattern of solar filaments is discussed. Finally, recent advances in non-eruptive global MHD models are described. While the review focuses mainly on solar magnetic fields, recent advances in measuring and modeling stellar magnetic fields are described where appropriate. In the final section key areas of future research are identified. △ Less

Submitted 28 November, 2012; originally announced November 2012.

Comments: 61 pages, 18 figures. Published version available via http://www.livingreviews.org/lrsp-2012-6

Journal ref: Living Rev. Solar Phys., 9 (2012) 6

Solar Magnetic Carpet II: Coronal Interactions of Small-Scale Magnetic Fields

Authors: K. A. Meyer, D. H. Mackay, A. A. van Ballegooijen

Abstract: This paper is the second in a series of studies working towards constructing a realistic, evolving, non-potential coronal model for the solar magnetic carpet. In the present study, the interaction of two magnetic elements is considered. Our objectives are to study magnetic energy build up, storage and dissipation as a result of emergence, cancellation, and flyby of these magnetic elements. In the… ▽ More This paper is the second in a series of studies working towards constructing a realistic, evolving, non-potential coronal model for the solar magnetic carpet. In the present study, the interaction of two magnetic elements is considered. Our objectives are to study magnetic energy build up, storage and dissipation as a result of emergence, cancellation, and flyby of these magnetic elements. In the future these interactions will be the basic building blocks of more complicated simulations involving hundreds of elements. Each interaction is simulated in the presence of an overlying uniform magnetic field, which lies at various orientations with respect to the evolving magnetic elements. For these three small-scale interactions, the free energy stored in the field at the end of the simulation ranges from $0.2-2.1\times 10^{26}$ ergs, while the total energy dissipated ranges from $1.3-6.3\times 10^{26}$ ergs. For all cases, a stronger overlying field results in higher energy storage and dissipation. For the cancellation and emergence simulations, motion perpendicular to the overlying field results in the highest values. For the flyby simulations, motion parallel to the overlying field gives the highest values. In all cases, the free energy built up is sufficient to explain small-scale phenomena such as X-ray bright points or nanoflares. In addition, if scaled for the correct number of magnetic elements for the volume considered, the energy continually dissipated provides a significant fraction of the quiet Sun coronal heating budget. △ Less

Submitted 16 November, 2012; originally announced November 2012.

Journal ref: Solar Physics (2012), Volume 278, Issue 1, pp.149-175

Chirality of High Latitude Filaments over Solar Cycle 23

Authors: A. R. Yeates, D. H. Mackay

Abstract: A non-potential quasi-static evolution model coupling the Sun's photospheric and coronal magnetic fields is applied to the problem of filament chirality at high latitudes. For the first time, we run a continuous 15 year simulation, using bipolar active regions determined from US National Solar Observatory, Kitt Peak magnetograms between 1996 and 2011. Using this simulation, we are able to address… ▽ More A non-potential quasi-static evolution model coupling the Sun's photospheric and coronal magnetic fields is applied to the problem of filament chirality at high latitudes. For the first time, we run a continuous 15 year simulation, using bipolar active regions determined from US National Solar Observatory, Kitt Peak magnetograms between 1996 and 2011. Using this simulation, we are able to address the outstanding question of whether magnetic helicity transport from active latitudes can overcome the effect of differential rotation at higher latitudes. Acting alone, differential rotation would produce high latitude filaments with opposite chirality to the majority type in each hemisphere. We find that differential rotation can indeed lead to opposite chirality at high latitudes, but only for around 5 years of the solar cycle following the polar field reversal. At other times, including the rising phase, transport of magnetic helicity from lower latitudes overcomes the effect of in situ differential rotation, producing the majority chirality even on the polar crowns at polar field reversal. These simulation predictions will allow for future testing of the non-potential coronal model. The results indicate the importance of long-term memory and helicity transport from active latitudes when modeling the structure and topology of the coronal magnetic field at higher latitudes. △ Less

Submitted 11 June, 2012; originally announced June 2012.

Comments: 15 pages, 5 figures, accepted for publication in ApJ Letters

Using Kepler transit observations to measure stellar spot belt migration rates

Authors: J. Llama, M. Jardine, D. H. Mackay, R. Fares

Abstract: Planetary transits provide a unique opportunity to investigate the surface distributions of star spots. Our aim is to determine if, with continuous observation (such as the data that will be provided by the Kepler mission), we can in addition measure the rate of drift of the spot belts. We begin by simulating magnetic cycles suitable for the Sun and more active stars, incorporating both flux emerg… ▽ More Planetary transits provide a unique opportunity to investigate the surface distributions of star spots. Our aim is to determine if, with continuous observation (such as the data that will be provided by the Kepler mission), we can in addition measure the rate of drift of the spot belts. We begin by simulating magnetic cycles suitable for the Sun and more active stars, incorporating both flux emergence and surface transport. This provides the radial magnetic field distribution on the stellar surface as a function of time. We then model the transit of a planet whose orbital axis is misaligned with the stellar rotation axis. Such a planet could occult spots at a range of latitudes. This allows us to complete the forward modelling of the shape of the transit lightcurve. We then attempt the inverse problem of recovering spot locations from the transit alone. From this we determine if transit lightcurves can be used to measure spot belt locations as a function of time. We find that for low-activity stars such as the Sun, the 3.5 year Kepler window is insufficient to determine this drift rate. For more active stars, it may be difficult to distinguish subtle differences in the nature of flux emergence, such as the degree of overlap of the "butterfly wings". The rate and direction of drift of the spot belts can however be determined for these stars. This would provide a critical test of dynamo theory. △ Less

Submitted 16 February, 2012; originally announced February 2012.

Comments: 5 pages. Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

Solar Magnetic Carpet I: Simulation of Synthetic Magnetograms

Authors: Karen A. Meyer, Duncan H. Mackay, Adriaan A. van Ballegooijen, Clare E. Parnell

Abstract: This paper describes a new 2D model for the photospheric evolution of the magnetic carpet. It is the first in a series of papers working towards constructing a realistic 3D non-potential model for the interaction of small-scale solar magnetic fields. In the model, the basic evolution of the magnetic elements is governed by a supergranular flow profile. In addition, magnetic elements may evolve thr… ▽ More This paper describes a new 2D model for the photospheric evolution of the magnetic carpet. It is the first in a series of papers working towards constructing a realistic 3D non-potential model for the interaction of small-scale solar magnetic fields. In the model, the basic evolution of the magnetic elements is governed by a supergranular flow profile. In addition, magnetic elements may evolve through the processes of emergence, cancellation, coalescence and fragmentation. Model parameters for the emergence of bipoles are based upon the results of observational studies. Using this model, several simulations are considered, where the range of flux with which bipoles may emerge is varied. In all cases the model quickly reaches a steady state where the rates of emergence and cancellation balance. Analysis of the resulting magnetic field shows that we reproduce observed quantities such as the flux distribution, mean field, cancellation rates, photospheric recycle time and a magnetic network. As expected, the simulation matches observations more closely when a larger, and consequently more realistic, range of emerging flux values is allowed (4e16 - 1e19 Mx). The model best reproduces the current observed properties of the magnetic carpet when we take the minimum absolute flux for emerging bipoles to be 4e16 Mx. In future, this 2D model will be used as an evolving photospheric boundary condition for 3D non-potential modeling. △ Less

Submitted 4 August, 2011; originally announced August 2011.

Comments: 33 pages, 16 figures, 5 gif movies included: movies may be viewed at http://www-solar.mcs.st-and.ac.uk/~karen/movies_paper1/

Journal ref: Solar Phys., 272, 29-58 2011

Modeling the Dispersal of an Active Region: Quantifying Energy Input into the Corona

Authors: Duncan H. Mackay, L. M. Green, Aad van Ballegooijen

Abstract: In this paper a new technique for modeling non-linear force-free fields directly from line of sight magnetogram observations is presented. The technique uses sequences of magnetograms directly as lower boundary conditions to drive the evolution of coronal magnetic fields between successive force-free equilibria over long periods of time. It is illustrated by applying it to MDI observations of a de… ▽ More In this paper a new technique for modeling non-linear force-free fields directly from line of sight magnetogram observations is presented. The technique uses sequences of magnetograms directly as lower boundary conditions to drive the evolution of coronal magnetic fields between successive force-free equilibria over long periods of time. It is illustrated by applying it to MDI observations of a decaying active region, NOAA AR 8005. The active region is modeled during a 4 day period around its central meridian passage. Over this time, the dispersal of the active region is dominated by random motions due to small scale convective cells. Through studying the build up of magnetic energy in the model, it is found that such small scale motions may inject anywhere from $2.5-3 \times 10^{25}$ erg s$^{-1}$ of free magnetic energy into the coronal field. Most of this energy is stored within the center of the active region in the low corona, below 30 Mm. After 4 days the build-up of free energy is 10% that of the corresponding potential field. This energy buildup, is sufficient to explain the radiative losses at coronal temperatures within the active region. Small scale convective motions therefore play an integral part in the energy balance of the corona. This new technique has wide ranging applications with the new high resolution, high cadence observations from the SDO:HMI and SDO:AIA instruments. △ Less

Submitted 25 February, 2011; originally announced February 2011.

Comments: 13 Figures, 2 movies

Journal ref: ApJ, 729, 97 (2011)

A Non-potential Model for the Sun's Open Magnetic Flux

Authors: A. R. Yeates, D. H. Mackay, A. A. van Ballegooijen, J. A. Constable

Abstract: Measurements of the interplanetary magnetic field (IMF) over several solar cycles do not agree with computed values of open magnetic flux from potential field extrapolations. The discrepancy becomes greater around solar maximum in each cycle, when the IMF can be twice as strong as predicted by the potential field model. Here we demonstrate that this discrepancy may be resolved by allowing for elec… ▽ More Measurements of the interplanetary magnetic field (IMF) over several solar cycles do not agree with computed values of open magnetic flux from potential field extrapolations. The discrepancy becomes greater around solar maximum in each cycle, when the IMF can be twice as strong as predicted by the potential field model. Here we demonstrate that this discrepancy may be resolved by allowing for electric currents in the low corona (below 2.5 solar radii). We present a quasi-static numerical model of the large-scale coronal magnetic evolution, which systematically produces these currents through flux emergence and shearing by surface motions. The open flux is increased by 75%-85% at solar maximum, but only 25% at solar minimum, bringing it in line with estimates from IMF measurements. The additional open flux in the non-potential model arises through inflation of the magnetic field by electric currents, with super-imposed fluctuations due to coronal mass ejections. The latter are modelled by the self-consistent ejection of twisted magnetic flux ropes. △ Less

Submitted 21 June, 2010; originally announced June 2010.

Comments: 10 pages, 8 figures, accepted for publication in J. Geophys. Res. A

Physics of Solar Prominences: II - Magnetic Structure and Dynamics

Authors: D. H. Mackay, J. T. Karpen, J. L. Ballester, B. Schmieder, G. Aulanier

Abstract: Observations and models of solar prominences are reviewed. We focus on non-eruptive prominences, and describe recent progress in four areas of prominence research: (1) magnetic structure deduced from observations and models, (2) the dynamics of prominence plasmas (formation and flows), (3) Magneto-hydrodynamic (MHD) waves in prominences and (4) the formation and large-scale patterns of the filam… ▽ More Observations and models of solar prominences are reviewed. We focus on non-eruptive prominences, and describe recent progress in four areas of prominence research: (1) magnetic structure deduced from observations and models, (2) the dynamics of prominence plasmas (formation and flows), (3) Magneto-hydrodynamic (MHD) waves in prominences and (4) the formation and large-scale patterns of the filament channels in which prominences are located. Finally, several outstanding issues in prominence research are discussed, along with observations and models required to resolve them. △ Less

Submitted 11 January, 2010; originally announced January 2010.

Comments: 75 pages, 31 pictures, review paper

Comparison of a Global Magnetic Evolution Model with Observations of Coronal Mass Ejections

Authors: A. R. Yeates, G. D. R. Attrill, Dibyendu Nandy, D. H. Mackay, P. C. H. Martens, A. A. van Ballegooijen

Abstract: The relative importance of different initiation mechanisms for coronal mass ejections (CMEs) on the Sun is uncertain. One possible mechanism is the loss of equilibrium of coronal magnetic flux ropes formed gradually by large-scale surface motions. In this paper, the locations of flux rope ejections in a recently-developed quasi-static global evolution model are compared with observed CME source… ▽ More The relative importance of different initiation mechanisms for coronal mass ejections (CMEs) on the Sun is uncertain. One possible mechanism is the loss of equilibrium of coronal magnetic flux ropes formed gradually by large-scale surface motions. In this paper, the locations of flux rope ejections in a recently-developed quasi-static global evolution model are compared with observed CME source locations over a 4.5-month period in 1999. Using EUV data, the low-coronal source locations are determined unambiguously for 98 out of 330 CMEs. Despite the incomplete observations, positive correlation (with coefficient up to 0.49) is found between the distributions of observed and simulated ejections, but only when binned into periods of one month or longer. This binning timescale corresponds to the time interval at which magnetogram data are assimilated into the coronal simulations, and the correlation arises primarily from the large-scale surface magnetic field distribution; only a weak dependence is found on the magnetic helicity imparted to the emerging active regions. The simulations are limited in two main ways: they produce fewer ejections, and they do not reproduce the strong clustering of observed CME sources into active regions. Due to this clustering, the horizontal gradient of radial photospheric magnetic field is better correlated with the observed CME source distribution (coefficient 0.67). Our results suggest that, while the gradual formation of magnetic flux ropes over weeks can account for many observed CMEs, especially at higher latitudes, there exists a second class of CMEs (at least half) for which dynamic active region flux emergence on shorter timescales must be the dominant factor. △ Less

Submitted 17 December, 2009; originally announced December 2009.

Comments: 27 pages, 7 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.709:1238-1248,2010

Initiation of Coronal Mass Ejections in a Global Evolution Model

Authors: A. R. Yeates, D. H. Mackay

Abstract: Loss of equilibrium of magnetic flux ropes is a leading candidate for the origin of solar coronal mass ejections (CMEs). The aim of this paper is to explore to what extent this mechanism can account for the initiation of CMEs in the global context. A simplified MHD model for the global coronal magnetic field evolution in response to flux emergence and shearing by large-scale surface motions is d… ▽ More Loss of equilibrium of magnetic flux ropes is a leading candidate for the origin of solar coronal mass ejections (CMEs). The aim of this paper is to explore to what extent this mechanism can account for the initiation of CMEs in the global context. A simplified MHD model for the global coronal magnetic field evolution in response to flux emergence and shearing by large-scale surface motions is described and motivated. Using automated algorithms for detecting flux ropes and ejections in the global magnetic model, the effects of key simulation parameters on the formation of flux ropes and the number of ejections are considered, over a 177-day period in 1999. These key parameters include the magnitude and sign of magnetic helicity emerging in active regions, and coronal diffusion. The number of flux ropes found in the simulation at any one time fluctuates between about 28 and 48, sustained by the emergence of new bipolar regions, but with no systematic dependence on the helicity of these regions. However, the emerging helicity does affect the rate of flux rope ejections, which doubles from 0.67 per day if the bipoles emerge untwisted to 1.28 per day in the run with greatest emerging twist. The number of ejections in the simulation is also increased by 20%-30% by choosing the majority sign of emerging bipole helicity in each hemisphere, or by halving the turbulent diffusivity in the corona. For reasonable parameter choices, the model produces approximately 50% of the observed CME rate. This indicates that the formation and loss of equilibrium of flux ropes may be a key element in explaining a significant fraction of observed CMEs. △ Less

Submitted 28 April, 2009; originally announced April 2009.

Comments: 15 pages, 9 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.699:1024-1037,2009

Modelling the Global Solar Corona III: Origin of the Hemispheric Pattern of Filaments

Authors: A. R. Yeates, D. H. Mackay

Abstract: We consider the physical origin of the hemispheric pattern of filament chirality on the Sun. Our 3D simulations of the coronal field evolution over a period of 6 months, based on photospheric magnetic measurements, were previously shown to be highly successful at reproducing observed filament chiralities. In this paper we identify and describe the physical mechanisms responsible for this success… ▽ More We consider the physical origin of the hemispheric pattern of filament chirality on the Sun. Our 3D simulations of the coronal field evolution over a period of 6 months, based on photospheric magnetic measurements, were previously shown to be highly successful at reproducing observed filament chiralities. In this paper we identify and describe the physical mechanisms responsible for this success. The key mechanisms are found to be (1) differential rotation of north-south polarity inversion lines, (2) the shape of bipolar active regions, and (3) evolution of skew over a period of many days. As on the real Sun, the hemispheric pattern in our simulations holds in a statistical sense. Exceptions arise naturally for filaments in certain locations relative to bipolar active regions, or from interactions between a number of active regions. △ Less

Submitted 2 October, 2008; originally announced October 2008.

Comments: 14 pages, 5 figures, accepted for publication in Solar Physics

Evolution and Distribution of Current Helicity in Full-Sun Simulations

Authors: A. R. Yeates, D. H. Mackay, A. A. van Ballegooijen

Abstract: Current helicity quantifies the location of twisted and sheared non-potential structures in a magnetic field. We simulate the evolution of magnetic fields in the solar atmosphere in response to flux emergence and shearing by photospheric motions. In our global-scale simulation over many solar rotations the latitudinal distribution of current helicity develops a clear statistical pattern, matchin… ▽ More Current helicity quantifies the location of twisted and sheared non-potential structures in a magnetic field. We simulate the evolution of magnetic fields in the solar atmosphere in response to flux emergence and shearing by photospheric motions. In our global-scale simulation over many solar rotations the latitudinal distribution of current helicity develops a clear statistical pattern, matching the observed hemispheric sign at active latitudes. In agreement with observations there is significant scatter and intermixing of both signs of helicity, where we find local values of current helicity density that are much higher than those predicted by linear force-free extrapolations. Forthcoming full-disk vector magnetograms from Solar Dynamics Observatory will provide an ideal opportunity to test our theoretical results on the evolution and distribution of current helicity, both globally and in single active regions. △ Less

Submitted 13 May, 2008; originally announced May 2008.

Comments: 4 pages, 3 figures, accepted for publication in ApJ Letters

Modelling the Global Solar Corona II: Coronal Evolution and Filament Chirality Comparison

Authors: A. R. Yeates, D. H. Mackay, A. A. van Ballegooijen

Abstract: The hemispheric pattern of solar filaments is considered using newly-developed simulations of the real photospheric and 3D coronal magnetic fields over a 6-month period, on a global scale. The magnetic field direction in the simulation is compared directly with the chirality of observed filaments, at their observed locations. In our model the coronal field evolves through a continuous sequence o… ▽ More The hemispheric pattern of solar filaments is considered using newly-developed simulations of the real photospheric and 3D coronal magnetic fields over a 6-month period, on a global scale. The magnetic field direction in the simulation is compared directly with the chirality of observed filaments, at their observed locations. In our model the coronal field evolves through a continuous sequence of nonlinear force-free equilibria, in response to the changing photospheric boundary conditions and the emergence of new magnetic flux. In total 119 magnetic bipoles with properties matching observed active regions are inserted. These bipoles emerge twisted and inject magnetic helicity into the solar atmosphere. When we choose the sign of this active-region helicity to match that observed in each hemisphere, the model produces the correct chirality for up to 96% of filaments, including exceptions to the hemispheric pattern. If the emerging bipoles have zero helicity, or helicity of the opposite sign, then this percentage is much reduced. In addition, the simulation produces a higher proportion of filaments with the correct chirality after longer times. This indicates that a key element in the evolution of the coronal field is its long-term memory, and the build-up and transport of helicity from low to high latitudes over many months. It highlights the importance of continuous evolution of the coronal field, rather than independent extrapolations at different times. This has significant consequences for future modelling such as that related to the origin and development of coronal mass ejections. △ Less

Submitted 19 November, 2007; originally announced November 2007.

Comments: 21 pages, 6 figures, accepted for publication in Solar Physics (Springer)

Exploring the Physical Basis of Solar Cycle Predictions: Flux Transport Dynamics and Persistence of Memory in Advection versus Diffusion Dominated Solar Convection Zones

Authors: A. R. Yeates, D. Nandy, D. H. Mackay

Abstract: The predictability, or lack thereof, of the solar cycle is governed by numerous separate physical processes that act in unison in the interior of the Sun. Magnetic flux transport and the finite time delay it introduces, specifically in the so-called Babcock-Leighton models of the solar cycle with spatially segregated source regions for the alpha and omega effects, play a crucial rule in this pre… ▽ More The predictability, or lack thereof, of the solar cycle is governed by numerous separate physical processes that act in unison in the interior of the Sun. Magnetic flux transport and the finite time delay it introduces, specifically in the so-called Babcock-Leighton models of the solar cycle with spatially segregated source regions for the alpha and omega effects, play a crucial rule in this predictability. Through dynamo simulations with such a model, we study the physical basis of solar cycle predictions by examining two contrasting regimes, one dominated by diffusive magnetic flux transport in the solar convection zone, the other dominated by advective flux transport by meridional circulation. Our analysis shows that diffusion plays an important role in flux transport, even when the solar cycle period is governed by the meridional flow speed. We further examine the persistence of memory of past cycles in the advection and diffusion dominated regimes through stochastically forced dynamo simulations. We find that in the advection-dominated regime, this memory persists for up to three cycles, whereas in the diffusion-dominated regime, this memory persists for mainly one cycle. This indicates that solar cycle predictions based on these two different regimes would have to rely on fundamentally different inputs - which may be the cause of conflicting predictions. Our simulations also show that the observed solar cycle amplitude-period relationship arises more naturally in the diffusion dominated regime, thereby supporting those dynamo models in which diffusive flux transport plays a dominant role in the solar convection zone. △ Less

Submitted 7 September, 2007; originally announced September 2007.

Comments: 35 pages, 12 figures, submitted to ApJ

Modelling the Global Solar Corona: Filament Chirality Observations and Surface Simulations

Authors: A. R. Yeates, D. H. Mackay, A. A. van Ballegooijen

Abstract: The hemispheric pattern of solar filaments is considered in the context of the global magnetic field of the solar corona. In recent work Mackay and van Ballegooijen have shown how, for a pair of interacting magnetic bipoles, the observed chirality pattern could be explained by the dominant range of bipole tilt angles and helicity in each hemisphere. This study aims to test this earlier result th… ▽ More The hemispheric pattern of solar filaments is considered in the context of the global magnetic field of the solar corona. In recent work Mackay and van Ballegooijen have shown how, for a pair of interacting magnetic bipoles, the observed chirality pattern could be explained by the dominant range of bipole tilt angles and helicity in each hemisphere. This study aims to test this earlier result through a direct comparison between theory and observations, using newly-developed simulations of the actual surface and 3D coronal magnetic fields over a 6-month period, on a global scale. In this paper we consider two key components of the study; firstly the observations of filament chirality for the sample of 255 filaments, and secondly our new simulations of the large-scale surface magnetic field. Based on a flux-transport model, these will be used as the lower boundary condition for the future 3D coronal simulations. Our technique differs significantly from those of other authors, where the coronal field is either assumed to be purely potential, or has to be reset back to potential every 27 days in order that the photospheric field remain accurate. In our case we ensure accuracy by the insertion of newly-emerging bipolar active regions, based on observed photospheric synoptic magnetograms. The large-scale surface field is shown to remain accurate over the 6-month period, without any resetting. This new technique will enable future simulations to consider the long-term build-up and transport of helicity and shear in the coronal magnetic field, over many months or years. △ Less

Submitted 22 July, 2007; originally announced July 2007.

Comments: 26 pages, 11 figures, accepted for publication in Solar Physics (Springer)

The impact of meridional circulation on stellar butterfly diagrams and polar caps

Authors: V. Holzwarth, D. H. Mackay, M. Jardine

Abstract: Observations of rapidly rotating solar-like stars show a significant mixture of opposite-polarity magnetic fields within their polar regions. To explain these observations, models describing the surface transport of magnetic flux demand the presence of fast meridional flows. Here, we link sub-surface and surface magnetic flux transport simulations to investigate (i) the impact of meridional circ… ▽ More Observations of rapidly rotating solar-like stars show a significant mixture of opposite-polarity magnetic fields within their polar regions. To explain these observations, models describing the surface transport of magnetic flux demand the presence of fast meridional flows. Here, we link sub-surface and surface magnetic flux transport simulations to investigate (i) the impact of meridional circulations with peak velocities of <125m/s on the latitudinal eruption pattern of magnetic flux tubes and (ii) the influence of the resulting butterfly diagrams on polar magnetic field properties. Prior to their eruption, magnetic flux tubes with low field strengths and initial cross sections below about 300km experience an enhanced poleward deflection through meridional flows. In particular flux tubes which originate between low and intermediate latitudes within the convective overshoot region are strongly affected. This latitude-dependent poleward deflection of erupting magnetic flux renders the wings of stellar butterfly diagrams distinctively convex. The subsequent evolution of the surface magnetic field shows that the increased number of newly emerging bipoles at higher latitudes promotes the intermingling of opposite polarities of polar magnetic fields. The associated magnetic flux densities are about 20% higher than in the case disregarding the pre-eruptive deflection, which eases the necessity for fast meridional flows predicted by previous investigations. In order to reproduce the observed polar field properties, the rate of the meridional circulation has to be on the order of 100m/s, and the latitudinal range from which magnetic flux tubes originate at the base of the convective zone (<50degrees) must be larger than in the solar case (<35degrees). △ Less

Submitted 5 April, 2006; originally announced April 2006.

Comments: 15 pages, 12 figures, accepted for publication by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.369:1703-1718,2006

Inferring X-ray coronal structures from Zeeman-Doppler images

Authors: M. Jardine, K. Wood, A. Collier Cameron, J. -F. Donati, D. H. Mackay

Abstract: We have modelled the X-ray emission from the young rapid rotator AB Doradus (P = 0.514 days) using as a basis Zeeman-Doppler maps of the surface magnetic field. This allows us to reconcile the apparently conflicting observations of a high X-ray emission measure and coronal density with a low rotational modulation in X-rays. The technique is to extrapolate the coronal field from the surface maps… ▽ More We have modelled the X-ray emission from the young rapid rotator AB Doradus (P = 0.514 days) using as a basis Zeeman-Doppler maps of the surface magnetic field. This allows us to reconcile the apparently conflicting observations of a high X-ray emission measure and coronal density with a low rotational modulation in X-rays. The technique is to extrapolate the coronal field from the surface maps by assuming the field to be potential. We then determine the coronal density for an isothermal corona by solving hydrostatic equilibrium along each field line and scaling the plasma pressure at the loop footpoints with the magnetic pressure. We set the density to zero along those field lines that are open and those where at any point along their length the plasma pressure exceeds the magnetic pressure. We then calculate the optically thin X-ray emission measure and rotational modulation for models with a range of coronal densities. Although the corona can be very extended, much of the emission comes from high-latitude regions close to the stellar surface. Since these are always in view as the star rotates, there is little rotational modulation. We find that emission measures in the observed range of 10^52.8 - 10^53.3 per cubic centimetre can be reproduced with densities in the range 10^9 - 10^10.7 per cubic centimetre for coronae at temperatures of 10^6-10^7 K. △ Less

Submitted 24 July, 2002; originally announced July 2002.

Comments: 7 pages, 9 figures, to be published in Monthly Notices of the Royal Astronomical Society

Journal ref: Mon.Not.Roy.Astron.Soc. 336 (2002) 1364