Source author record

D. H. Mackay

D. H. Mackay appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2020arXiv

Hydrogen Non-Equilibrium Ionisation Effects in Coronal Mass Ejections

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.

preprint2020arXiv

Measuring stellar magnetic helicity density

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.

preprint2019arXiv

A New Space Weather Tool for Identifying Eruptive Active Regions

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Stellar Coronal Response to Differential Rotation and Flux Emergence

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.

preprint2015arXiv

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

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.

preprint2014arXiv

Stellar Differential Rotation and Coronal Timescales

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.

preprint2013arXiv

Solar Magnetic Carpet III: Coronal Modelling of Synthetic Magnetograms

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.

preprint2013arXiv

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

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.

preprint2012arXiv

Chirality of High Latitude Filaments over Solar Cycle 23

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Using Kepler transit observations to measure stellar spot belt migration rates

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.

preprint2010arXiv

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

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.

preprint2010arXiv

Physics of Solar Prominences: II - Magnetic Structure and Dynamics

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.

preprint2009arXiv

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

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.

preprint2009arXiv

Initiation of Coronal Mass Ejections in a Global Evolution Model

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.

preprint2007arXiv

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

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.