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.

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-$\alpha$ 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-$\alpha$ 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-$\alpha$ 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Magnetic reconnection, the rearrangement of magnetic field topology, is a fundamental physical process in magnetized plasma systems all over the universe1,2. Its process is difficult to be directly observed. Coronal structures, such as coronal loops and filament spines, often sketch the magnetic field geometry and its changes in the solar corona3. Here we show a highly suggestive observation of magnetic reconnection between an erupting solar filament and its nearby coronal loops, resulting in changes in connection of the filament. X-type structures form when the erupting filament encounters the loops. The filament becomes straight, and bright current sheets form at the interfaces with the loops. Many plasmoids appear in these current sheets and propagate bi-directionally. The filament disconnects from the current sheets, which gradually disperse and disappear, reconnects to the loops, and becomes redirected to the loop footpoints. This evolution of the filament and the loops suggests successive magnetic reconnection predicted by theories1 but rarely detected with such clarity in observations. Our results on the formation, evolution, and disappearance of current sheets, confirm three-dimensional magnetic reconnection theory and have implications for the evolution of dissipation regions and the release of magnetic energy for reconnection in many magnetized plasma systems.

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.

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.

Understanding the 3D structure of coronal mass ejections (CMEs) is crucial for understanding the nature and origin of solar eruptions. To derive information on the 3D structure of CMEs from white-light (total and polarized brightness) images, the polarization ratio technique is widely used. The soon-to-be-launched METIS coronagraph on board Solar Orbiter will use this technique to produce new polarimetric images. We determine the accuracy at which the position of the centre of mass, direction and speed of propagation, and the column density of the CME can be determined along the line of sight. We perform a 3D MHD simulation of a flux rope ejection where a CME is produced. From the simulation we (i) synthesize the corresponding METIS white-light (total and polarized brightness) images and (ii) apply the polarization ratio technique to these synthesized images and compare the results with the known density distribution from the MHD simulation. We find that the polarization ratio technique reproduces with high accuracy the position of the centre of mass along the line of sight. However, some errors are inherently associated with this determination. The polarization ratio technique also allows information to be derived on the real 3D direction of propagation of the CME. In addition, we find that the column density derived from white-light images is accurate and we propose an improved technique where the combined use of the polarization ratio technique and white-light images minimizes the error in the estimation of column densities. Our method allows us to thoroughly test the performance of the polarization ratio technique and allows a determination of the errors associated with it, which means that it can be used to quantify the results from the analysis of the forthcoming METIS observations in white light (total and polarized brightness).

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Recent advances in spectropolarimetry have allowed the reconstruction of stellar coronal magnetic fields. This uses Zeeman-Doppler magnetograms of the surface magnetic field as a lower boundary condition. The ZDI maps, however, suffer from the absence of information about the magnetic field over regions of the surface due to the presence of dark starspots and portions of the surface out of view due to a tilt in the rotation axis. They also suffer from finite resolution which leads to small scale field structures being neglected. This paper explores the effects of this loss of information on the extrapolated coronal fields. For this we use simulated stellar surface magnetic maps for two hypothetical stars. Using the potential field approximation, the coronal fields and emission measures are calculated. This is repeated for the cases of missing information due to, (i) starspots, (ii) a large area of the stellar surface out of view, (iii) a finite resolution. The largest effect on the magnetic field structure arises when a significant portion of the stellar surface remains out of view. This changes the nature of the field lines that connect to this obscured hemisphere. Nonetheless, the field structure in the visible hemisphere is reliably reproduced. Thus the calculation of the locations and surface filling factors of accretion funnels is reasonably well reproduced for the observed hemisphere. The decrease with height of the magnetic pressure, which is important in calculating disc truncation radii for accreting stars, is also largely unaffected in the equatorial plane. The fraction of surface flux that is open and therefore able to supply angular momentum loss in a wind, however, is often overestimated in the presence of missing flux.

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.

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.

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.

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.

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.

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.

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.

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.

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).

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.

Infrared magnitude-redshift relations for the 3CR and 6C samples of radio galaxies are presented for a wide range of plausible cosmological models, including those with non-zero cosmological constant OmegaLambda. Variations in the galaxy formation redshift, metallicity and star formation history are also considered. The results of the modelling are displayed in terms of magnitude differences between the models and no-evolution tracks, illustrating the amount of K-band evolution necessary to account for the observational data. Given a number of plausible assumptions, the results of these analyses suggest that: (i) cosmologies which predict T_0xH_0>1 (where T_0 denotes the current age of the universe) can be excluded; (ii) the star formation redshift should lie in the redshift interval 5<z<20, values towards the lower end of the range being preferred in cosmologies with larger values of T_0xH_0; (iii) the Einstein-de Sitter model provides a reasonable fit to the data; (iv) models with finite values of OmegaLambda can provide good agreement with the observations only if appropriate adjustments of other parameters such as the galaxy metallicities and star-formation histories are made. Without such modifications, even after accounting for stellar evolution, the high redshift radio galaxies are more luminous (ie. more massive) than those nearby in models with finite OmegaLambda, including the favoured model with Omega=0.3, OmegaLambda=0.7. For cosmological models with larger values of T_0xH_0, the conclusions are the same regardless of whether any adjustments are made or not. The implications of these results for cosmology and models of galaxy formation are discussed.