In this paper we address a discrepancy between the surface flux evolution in a 3D kinematic dynamo model and a 2D surface flux transport model that has been closely calibrated to the real Sun. We demonstrate that the difference is due to the connectivity of active regions to the toroidal field at the base of the convection zone, which is not accounted for in the surface-only model. Initially, we consider the decay of a single active region, firstly in a simplified Cartesian 2D model and subsequently the full 3D model. By varying the turbulent diffusivity profile in the convection zone, we find that increasing the diffusivity - so that active regions are more rapidly disconnected from the base of the convection zone - improves the evolution of the surface field. However, if we simulate a full solar cycle, we find that the dynamo is unable to sustain itself under such an enhanced diffusivity. This suggests that in order to accurately model the solar cycle, we must find an alternative way to disconnect emerging active regions, whilst conserving magnetic flux.

Models for astrophysical plasmas often have magnetic field lines that leave the boundary rather than closing within the computational domain. Thus, the relative magnetic helicity is frequently used in place of the usual magnetic helicity, so as to restore gauge invariance. We show how to decompose the relative helicity into a relative field-line helicity that is an ideal-magnetohydrodynamic invariant for each individual magnetic field line, and vanishes along any field line where the original field matches the reference field. Physically, this relative field-line helicity is a magnetic flux, whose specific definition depends on the gauge of the reference vector potential on the boundary. We propose a particular `minimal' gauge that depends only on the reference field and minimises this boundary contribution, so as to reveal topological information about the original magnetic field. We illustrate the effect of different gauge choices using the Low-Lou and Titov-Demoulin models of solar active regions. Our numerical code to compute appropriate vector potentials and relative field-line helicity in Cartesian domains is open source and freely available.

Coronal jets and bright points occur prolifically in predominantly unipolar magnetic regions, such as coronal holes, where they appear above minority-polarity intrusions. Intermittent low-level reconnection and explosive, high-energy-release reconnection above these intrusions are thought to generate bright points and jets, respectively. The magnetic field above the intrusions possesses a spine-fan topology with a coronal null point. The movement of magnetic flux by surface convection adds free energy to this field, forming current sheets and inducing reconnection. We conducted three-dimensional magnetohydrodynamic simulations of moving magnetic elements as a model for coronal jets and bright points. A single minority-polarity concentration was subjected to three different experiments: a large-scale surface flow that sheared part of the separatrix surface only, a large-scale surface flow that also sheared part of the polarity inversion line surrounding the minority flux, and the latter flow setup plus a "fly-by" of a majority-polarity concentration past the moving minority-polarity element. We found that different bright-point morphologies, from simple loops to sigmoids, were created. When only the field near the separatrix was sheared, steady interchange reconnection modulated by quasi-periodic, low-intensity bursts of reconnection occurred, suggestive of a bright point with periodically varying intensity. When the field near the PIL was strongly sheared, on the other hand, filament channels repeatedly formed and erupted via the breakout mechanism, explosively increasing the interchange reconnection and generating non-helical jets. The fly-by produced even more energetic and explosive jets. Our results explain several key aspects of coronal-hole bright points and jets, and the relationships between them.

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.

We test recent claims that the polar field at the end of Cycle 23 was weakened by a small number of large, abnormally oriented regions, and investigate what this means for solar cycle prediction. We isolate the contribution of individual regions from magnetograms for Cycles 21, 22 and 23 using a 2D surface flux transport model, and find that although the top ~10% of contributors tend to define sudden large variations in the axial dipole moment, the cumulative contribution of many weaker regions cannot be ignored. In order to recreate the axial dipole moment to a reasonable degree, many more regions are required in Cycle 23 than in Cycles 21 and 22 when ordered by contribution. We suggest that the negative contribution of the most significant regions of Cycle 23 could indeed be a cause of the weak polar field at the following cycle minimum and the low-amplitude Cycle 24. We also examine the relationship between a region's axial dipole moment contribution and its emergence latitude, flux, and initial axial dipole moment. We find that once the initial dipole moment of a given region has been measured, we can predict the long-term dipole moment contribution using emergence latitude alone.

The topology of magnetic fields near the open-closed flux boundary in the Sun's corona is an important influencing factor in the process of interchange reconnection, whereby plasma is exchanged between open and closed flux domains. Maps of the magnetic squashing factor at the radial outer boundary in coronal field models reveal the presence of the so-called `S-web', and suggest that interchange reconnection could potentially deposit closed coronal material into high-latitude regions far from the heliospheric current sheet. Here we demonstrate that certain features of the S-web reveal the underlying topological structure of the magnetic field. Specifically, in order for the arcing bands of highly squashed magnetic flux of the S-web to terminate or intersect away from the helmet streamer apex, there must be a null spine line that maps a finite segment of the photospheric open-closed boundary up to a singular point in the open flux domain. We propose that this association between null spine lines and arc termination points may be used to identify locations in the heliosphere that are preferential for the appearance of solar energetic particles or slow solar wind plasma with certain characteristics.

In this paper we seek to understand the timescale on which the photospheric motions on the Sun braid coronal magnetic field lines. This is a crucial ingredient for determining the viability of the braiding mechanism for explaining the high temperatures observed in the corona. We study the topological complexity induced in the coronal magnetic field, primarily using plasma motions extracted from magneto-convection simulations. This topological complexity is quantified using the field line winding, finite time topological entropy and passive scalar mixing. With these measures we contrast mixing efficiencies of the magneto-convection simulation, a benchmark flow known as a ``blinking vortex', and finally photospheric flows inferred from sequences of observed magnetograms using local correlation tracking. While the highly resolved magneto-convection simulations induce a strong degree of field line winding and finite time topological entropy, the values obtained from the observations from the plage region are around an order of magnitude smaller. This behavior is carried over to the finite time topological entropy. Nevertheless, the results suggest that the photospheric motions induce complex tangling of the coronal field on a timescale of hours.

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.

On 21 August 2017 a total solar eclipse swept across the contiguous United States providing excellent opportunities for diagnostics of the Sun's corona. The Sun's coronal structure is notoriously difficult to observe except during solar eclipses; thus theoretical models must be relied upon for inferring the underlying magnetic structure of the Sun's outer atmosphere. These models are necessary for understanding the role of magnetic fields in the heating of the corona to a million degrees and generation of severe space weather. Here we present a methodology for predicting the structure of the coronal field based on long-term surface flux transport simulations whose output is utilized to extrapolate the coronal magnetic field structures. This prescription was applied to the 21 August 2017 solar eclipse. Post-eclipse analysis shows good agreement between model simulated and observed coronal structures and their locations on the limb. We demonstrate that slow changes in the Sun's surface magnetic field distribution driven by long-term flux emergence and evolution govern large-scale coronal structures with a (plausibly cycle-phase dependent) dynamical memory timescale on the order of few solar rotations -- opening up the possibility of large-scale, global corona predictions at least a month in advance.

Predictions of the solar wind at Earth are a central aspect of space weather prediction. The outcome of such a prediction, however, is highly sensitive to the method used for computing the magnetic field in the corona. We analyze the impact of replacing the potential field coronal boundary conditions, as used in operational space weather prediction tools, by non-potential conditions. For this, we compare the predicted solar wind plasma parameters with observations at 1 AU for two six-months intervals, one at solar maximum and one in the descending phase of the current cycle. As a baseline, we compare with the operational Wang-Sheeley-Arge model. We find that for solar maximum, the non-potential coronal model and an adapted solar wind speed formula lead to the best solar wind predictions in a statistical sense. For the descending phase, the potential coronal model performs best. The Wang-Sheeley-Arge model outperforms the others in predicting high speed enhancements and streamer interactions. A better parameter fitting for the adapted wind speed formula is expected to improve the performance of the non-potential model here.

The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field regions in the model should approximately correspond to coronal holes observed in emission, and (2) the magnitude of the open magnetic flux in the model should match that inferred from in situ spacecraft measurements. In this study, we calculate both MHD and PFSS solutions using fourteen different magnetic maps produced from five different types of observatory magnetograms, for the time period surrounding July, 2010. We have found that for all of the model/map combinations, models that have coronal hole areas close to observations underestimate the interplanetary magnetic flux, or, conversely, for models to match the interplanetary flux, the modeled open field regions are larger than coronal holes observed in EUV emission. In an alternative approach, we estimate the open magnetic flux entirely from solar observations by combining automatically detected coronal holes for Carrington rotation 2098 with observatory synoptic magnetic maps. This approach also underestimates the interplanetary magnetic flux. Our results imply that either typical observatory maps underestimate the Sun's magnetic flux, or a significant portion of the open magnetic flux is not rooted in regions that are obviously dark in EUV and X-ray emission.

Formed through magnetic field shearing and reconnection in the solar corona, magnetic flux ropes are structures of twisted magnetic field, threaded along an axis. Their evolution and potential eruption are of great importance for space weather. Here we describe a new methodology for the automated detection of flux ropes in simulated magnetic fields, based on fieldline helicity. More robust than previous methods, it also allows the full volume extent of each flux rope structure to be determined. Our Flux Rope Detection and Organization (FRoDO) code is publicly available, and measures the magnetic flux and helicity content of pre-erupting flux ropes over time, as well as detecting eruptions. As a first demonstration the code is applied to the output from a time-dependent magnetofrictional model, spanning 1996 June 15 - 2014 February 10. Over this period, 1561 erupting and 2099 non-erupting magnetic flux ropes are detected, tracked, and characterized. For this particular model data, erupting flux ropes have a mean net helicity magnitude of $(2.66 \pm 6.82) \times 10^{43}$ Mx$^2$, while non-erupting flux ropes have a significantly lower mean of $(4.04 \pm 9.25) \times 10^{42}$ Mx$^2$, although there is overlap between the two distributions. Similarly, the mean unsigned magnetic flux for erupting flux ropes is $(4.04 \pm 6.17) \times 10^{21}$ Mx, significantly higher than the mean value of $(7.05 \pm 16.8) \times 10^{20}$ Mx for non-erupting ropes. These values for erupting flux ropes are within the broad range expected from observational and theoretical estimates, although the eruption rate in this particular model is lower than that of observed coronal mass ejections. In future the FRoDO code will prove a valuable tool for assessing the performance of different non-potential coronal simulations and comparing them with observations.

Accurate prediction of solar activity calls for precise calibration of solar cycle models. Consequently we aim to find optimal parameters for models which describe the physical processes on the solar surface, which in turn act as proxies for what occurs in the interior and provide source terms for coronal models. We use a genetic algorithm to optimize surface flux transport models using National Solar Observatory (NSO) magnetogram data for Solar Cycle 23. This is applied to both a 1D model that inserts new magnetic flux in the form of idealized bipolar magnetic regions, and also to a 2D model that assimilates specific shapes of real active regions. The genetic algorithm searches for parameter sets (meridional flow speed and profile, supergranular diffusivity, initial magnetic field, and radial decay time) that produce the best fit between observed and simulated butterfly diagrams, weighted by a latitude-dependent error structure which reflects uncertainty in observations. Due to the easily adaptable nature of the 2D model, the optimization process is repeated for Cycles 21, 22, and 24 in order to analyse cycle-to-cycle variation of the optimal solution. We find that the ranges and optimal solutions for the various regimes are in reasonable agreement with results from the literature, both theoretical and observational. The optimal meridional flow profiles for each regime are almost entirely within observational bounds determined by magnetic feature tracking, with the 2D model being able to accommodate the mean observed profile more successfully. Differences between models appear to be important in deciding values for the diffusive and decay terms. In like fashion, differences in the behaviours of different solar cycles lead to contrasts in parameters defining the meridional flow and initial field strength.

By defining an appropriate field line helicity, we apply the powerful concept of magnetic helicity to the problem of global magnetic field evolution in the Sun's corona. As an ideal-magnetohydrodynamic invariant, the field line helicity is a meaningful measure of how magnetic helicity is distributed within the coronal volume. It may be interpreted, for each magnetic field line, as a magnetic flux linking with that field line. Using magneto-frictional simulations, we investigate how field line helicity evolves in the non-potential corona as a result of shearing by large-scale motions on the solar surface. On open magnetic field lines, the helicity injected by the Sun is largely output to the solar wind, provided that the coronal relaxation is sufficiently fast. But on closed magnetic field lines, helicity is able to build up. We find that the field line helicity is non-uniformly distributed, and is highly concentrated in twisted magnetic flux ropes. Eruption of these flux ropes is shown to lead to sudden bursts of helicity output, in contrast to the steady flux along the open magnetic field lines.

Accurate estimates of the horizontal electric field on the Sun's visible surface are important not only for estimating the Poynting flux of magnetic energy into the corona but also for driving time-dependent magnetohydrodynamic models of the corona. In this paper, a method is developed for estimating the horizontal electric field from a sequence of radial-component magnetic field maps. This problem of inverting Faraday's law has no unique solution. Unfortunately, the simplest solution (a divergence-free electric field) is not realistically localized in regions of non-zero magnetic field, as would be expected from Ohm's law. Our new method generates instead a localized solution, using a basis pursuit algorithm to find a sparse solution for the electric field. The method is shown to perform well on test cases where the input magnetic maps are flux balanced, in both Cartesian and spherical geometries. However, we show that if the input maps have a significant imbalance of flux - usually arising from data assimilation - then it is not possible to find a localized, realistic, electric field solution. This is the main obstacle to driving coronal models from time sequences of solar surface magnetic maps.

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.

Predicting the final state of turbulent plasma relaxation is an important challenge, both in astrophysical plasmas such as the Sun's corona and in controlled thermonuclear fusion. Recent numerical simulations of plasma relaxation with braided magnetic fields identified the possibility of a novel constraint, arising from the topological degree of the magnetic field-line mapping. This constraint implies that the final relaxed state is drastically different for an initial configuration with topological degree 1 (which allows a Taylor relaxation) and one with degree 2 (which does not reach a Taylor state). Here we test this transition in numerical resistive-magnetohydrodynamic simulations, by embedding a braided magnetic field in a linear force-free background. Varying the background force-free field parameter generates a sequence of initial conditions with a transition between topological degree 1 and 2. For degree 1, the relaxation produces a single twisted flux tube, while for degree 2 we obtain two flux tubes. For predicting the exact point of transition, it is not the topological degree of the whole domain that is relevant, but only that of the turbulent region.

As an observational case study, we consider the origin of a prominent poleward surge of leading polarity, visible in the magnetic butterfly diagram during Solar Cycle 24. A new technique is developed for assimilating individual regions of strong magnetic flux into a surface flux transport model. By isolating the contribution of each of these regions, the model shows the surge to originate primarily in a single high-latitude activity group consisting of a bipolar active region present in Carrington Rotations 2104-05 (November 2010-January 2011) and a multipolar active region in Rotations 2107-08 (February-April 2011). This group had a strong axial dipole moment opposed to Joy's law. On the other hand, the modelling suggests that the transient influence of this group on the butterfly diagram will not be matched by a large long-term contribution to the polar field, because of its location at high latitude. This is in accordance with previous flux transport models.

We investigate the evolution of field line helicity for magnetic fields that connect two boundaries without null points, with emphasis on localized finite-B magnetic reconnection. Total (relative) magnetic helicity is already recognized as an important topological constraint on magnetohydrodynamic processes. Field line helicity offers further advantages because it preserves all topological information and can distinguish between different magnetic fields with the same total helicity. Magnetic reconnection changes field connectivity and field line helicity reflects these changes; the goal of this paper is to characterize that evolution. We start by deriving the evolution equation for field line helicity and examining its terms, also obtaining a simplified form for cases where dynamics are localized within the domain. The main result, which we support using kinematic examples, is that during localized reconnection in a complex magnetic field, the evolution of field line helicity is dominated by a work-like term that is evaluated at the field line endpoints, namely the scalar product of the generalized field line velocity and the vector potential. Furthermore, the flux integral of this term over certain areas is very small compared to the integral of the unsigned quantity, which indicates that changes of field line helicity happen in a well-organized pairwise manner. It follows that reconnection is very efficient at redistributing helicity in complex magnetic fields despite having little effect on the total helicity.

In this work we take advantage of eleven different sunspot group, sunspot, and active region databases to characterize the area and flux distributions of photospheric magnetic structures. We find that, when taken separately, different databases are better fitted by different distributions (as has been reported previously in the literature). However, we find that all our databases can be reconciled by the simple application of a proportionality constant, and that, in reality, different databases are sampling different parts of a composite distribution. This composite distribution is made up by linear combination of Weibull and log-normal distributions -- where a pure Weibull (log-normal) characterizes the distribution of structures with fluxes below (above) $10^{21}$Mx ($10^{22}$Mx). We propose that this is evidence of two separate mechanisms giving rise to visible structures on the photosphere: one directly connected to the global component of the dynamo (and the generation of bipolar active regions), and the other with the small-scale component of the dynamo (and the fragmentation of magnetic structures due to their interaction with turbulent convection). Additionally, we demonstrate that the Weibull distribution shows the expected linear behavior of a power-law distribution (when extended into smaller fluxes), making our results compatible with the results of Parnell et al. (2009).

We reconsider the topological interpretation of magnetic helicity for magnetic fields in open domains, and relate this to the relative helicity. Specifically, our domains stretch between two parallel planes, and each of these ends may be magnetically open. It is demonstrated that, while the magnetic helicity is gauge-dependent, its value in any gauge may be physically interpreted as the average winding number among all pairs of field lines with respect to some orthonormal frame field. In fact, the choice of gauge is equivalent to the choice of reference field in the relative helicity, meaning that the magnetic helicity is no less physically meaningful. We prove that a particular gauge always measures the winding with respect to a fixed frame, and propose that this is normally the best choice. For periodic fields, this choice is equivalent to measuring relative helicity with respect to a potential reference field. But for aperiodic fields, we show that the potential field can be twisted. We prove by construction that there always exists a possible untwisted reference field.

We report on the first observation of a single hybrid magnetic structure that contains both a pseudostreamer and a double helmet streamer. This structure was originally observed by the SWAP instrument aboard the PROBA2 satellite between 5 and 10~May~2013. It consists of a pair of filament channels near the south pole of the sun. On the western edge of the structure, the magnetic morphology above the filaments is that of a side-by-side double helmet streamer, with open field between the two channels. On the eastern edge, the magnetic morphology is that of a coronal pseudostreamer without the central open field. We investigated this structure with multiple observations and modelling techniques. We describe the topology and dynamic consequences of such a unified structure.

We estimate the energy input into the solar corona from photospheric footpoint motions, using observations of a plage region by the Hinode Solar Optical Telescope. Assuming a perfectly ideal coronal evolution, two alternative lower bounds for the Poynting flux are computed based on field line footpoint trajectories, without requiring horizontal magnetic field data. When applied to the observed velocities, a bound based solely on displacements between the two footpoints of each field line is tighter than a bound based on relative twist between field lines. Depending on the assumed length of coronal magnetic field lines, the higher bound is found to be reasonably tight compared with a Poynting flux estimate using an available vector magnetogram. It is also close to the energy input required to explain conductive and radiative losses in the active region corona. Based on similar analysis of a numerical convection simulation, we suggest that observations with higher spatial resolution are likely to bring the bound based on relative twist closer to the first bound, but not to increase the first bound substantially. Finally, we put an approximate upper bound on the magnetic energy by constructing a hypothetical ``unrelaxed'' magnetic field with the correct field line connectivity.

We study coherent structures in solar photospheric flows in a plage in the vicinity of the active region AR 10930 using the horizontal velocity data derived from Hinode/SOT magnetograms. Eulerian and Lagrangian coherent structures are detected by computing the Q-criterion and the finite-time Lyapunov exponents of the velocity field, respectively. Our analysis indicates that, on average, the deformation Eulerian coherent structures dominate over the vortical Eulerian coherent structures in the plage region. We demonstrate the correspondence of the network of high magnetic flux concentration to the attracting Lagrangian coherent structures (a-LCS) in the photospheric velocity based on both observations and numerical simulations. In addition, the computation of a-LCS provides a measure of the local rate of contraction/expansion of the flow.

We propose a phenomenological technique for modelling the emergence of active regions within a three-dimensional, kinematic dynamo framework. By imposing localised velocity perturbations, we create emergent flux-tubes out of toroidal magnetic field at the base of the convection zone, leading to the eruption of active regions at the solar surface. The velocity perturbations are calibrated to reproduce observed active region properties (including the size and flux of active regions, and the distribution of tilt angle with latitude), resulting in a more consistent treatment of flux-tube emergence in kinematic dynamo models than artificial flux deposition. We demonstrate how this technique can be used to assimilate observations and drive a kinematic 3D model, and use it to study the characteristics of active region emergence and decay as a source of poloidal field. We find that the poloidal components are strongest not at the solar surface, but in the middle convection zone, in contrast with the common assumption that the poloidal source is located near the solar surface. We also find that, while most of the energy is contained in the lower convection zone, there is a good correlation between the evolution of the surface and interior magnetic fields.

The high-cadence, comprehensive view of the solar corona by SDO/AIA shows many events that are widely separated in space while occurring close together in time. In some cases, sets of coronal events are evidently causally related, while in many other instances indirect evidence can be found. We present case studies to highlight a variety of coupling processes involved in coronal events. We find that physical linkages between events do occur, but concur with earlier studies that these couplings appear to be crucial to understanding the initiation of major eruptive or explosive phenomena relatively infrequently. We note that the post-eruption reconfiguration time scale of the large-scale corona, estimated from the EUV afterglow, is on average longer than the mean time between CMEs, so that many CMEs originate from a corona that is still adjusting from a previous event. We argue that the coronal field is intrinsically global: current systems build up over days to months, the relaxation after eruptions continues over many hours, and evolving connections easily span much of a hemisphere. This needs to be reflected in our modeling of the connections from the solar surface into the heliosphere to properly model the solar wind, its perturbations, and the generation and propagation of solar energetic particles. However, the large-scale field cannot be constructed reliably by currently available observational resources. We assess the potential of high-quality observations from beyond Earth's perspective and advanced global modeling to understand the couplings between coronal events in the context of CMEs and solar energetic particle events.

A topological flux function is introduced to quantify the topology of magnetic braids: non-zero line-tied magnetic fields whose field lines all connect between two boundaries. This scalar function is an ideal invariant defined on a cross-section of the magnetic field, whose integral over the cross-section yields the relative magnetic helicity. Recognising that the topological flux function is an action in the Hamiltonian formulation of the field line equations, a simple formula for its differential is obtained. We use this to prove that the topological flux function uniquely characterises the field line mapping and hence the magnetic topology. A simple example is presented.

Coupled flux transport and magneto-frictional simulations are extended to simulate the continuous magnetic field evolution in the global solar corona for over 15 years, from the start of Solar Cycle 23 in 1996. By simplifying the dynamics, our model follows the build-up and transport of electric currents and free magnetic energy in the corona, offering an insight into the magnetic structure and topology that extrapolation-based models can not. To enable these extended simulations, we have implemented a more efficient numerical grid, and have carefully calibrated the surface flux transport model to reproduce the observed large-scale photospheric radial magnetic field, using emerging active regions determined from observed line-of-sight magnetograms. This calibration is described in some detail. In agreement with previous authors, we find that the standard flux transport model is insufficient to simultaneously reproduce the observed polar fields and butterfly diagram during Cycle 23, and that additional effects must be added. For the best-fit model, we use automated techniques to detect the latitude-time profile of flux ropes and their ejections over the full solar cycle. Overall, flux ropes are more prevalent outside of active latitudes but those at active latitudes are more frequently ejected. Future possibilities for space-weather prediction with this approach are briefly assessed.

We introduce a topological flux function to quantify the topology of magnetic braids: non-zero, line-tied magnetic fields whose field lines all connect between two boundaries. This scalar function is an ideal invariant defined on a cross-section of the magnetic field, and measures the average poloidal magnetic flux around any given field line, or the average pairwise crossing number between a given field line and all others. Moreover, its integral over the cross-section yields the relative magnetic helicity. Using the fact that the flux function is also an action in the Hamiltonian formulation of the field line equations, we prove that it uniquely characterizes the field line mapping and hence the magnetic topology.

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.

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.

Aims. We show how the build-up of magnetic gradients in the Sun's corona may be inferred directly from photospheric velocity data. This enables computation of magnetic connectivity measures such as the squashing factor without recourse to magnetic field extrapolation. Methods.Assuming an ideal evolution in the corona, and an initially uniform magnetic field, the subsequent field line mapping is computed by integrating trajectories of the (time-dependent) horizontal photospheric velocity field. The method is applied to a 12 hour high-resolution sequence of photospheric flows derived from Hinode/SOT magnetograms. Results. We find the generation of a network of quasi-separatrix layers in the magnetic field, which correspond to Lagrangian coherent structures in the photospheric velocity. The visual pattern of these structures arises primarily from the diverging part of the photospheric flow, hiding the effect of the rotational flow component: this is demonstrated by a simple analytical model of photospheric convection. We separate the diverging and rotational components from the observed flow and show qualitative agreement with purely diverging and rotational models respectively. Increasing the flow speeds in the model suggests that our observational results are likely to give a lower bound for the rate at which magnetic gradients are built up by real photospheric flows. Finally, we construct a hypothetical magnetic field with the inferred topology, that can be used for future investigations of reconnection and energy release.

We investigate the relaxation of braided magnetic loops in order to find out how the type of braiding via footpoint motions affects resultant heating of the loop. Two magnetic loops, braided in different ways, are used as initial conditions in resistive MHD simulations and their subsequent evolution is studied. The fields both undergo a resistive relaxation in which current sheets form and fragment and the system evolves towards a state of lower energy. In one case this relaxation is very efficient with current sheets filling the volume and homogeneous heating of the loop occurring. In the other case fewer current sheets develop, less magnetic energy is released in the process and a patchy heating of the loop results. The two cases, although very similar in their setup, can be distinguished by the mixing properties of the photospheric driver. The mixing can be measured by the topological entropy of the plasma flow, an observable quantity.

The definition and measurement of magnetic reconnection in three-dimensional magnetic fields with multiple reconnection sites is a challenging problem, particularly in fields lacking null points. We propose a generalization of the familiar two-dimensional concept of a magnetic flux function to the case of a three-dimensional field connecting two planar boundaries. Using hyperbolic fixed points of the field line mapping, and their global stable and unstable manifolds, we define a unique flux partition of the magnetic field. This partition is more complicated than the corresponding (well-known) construction in a two-dimensional field, owing to the possibility of heteroclinic points and chaotic magnetic regions. Nevertheless, we show how the partition reconnection rate is readily measured with the generalized flux function. We relate our partition reconnection rate to the common definition of three-dimensional reconnection in terms of integrated parallel electric field. An analytical example demonstrates the theory, and shows how the flux partition responds to an isolated reconnection event.

A topological constraint on the dynamics of a magnetic field in a flux tube arises from the fixed point indices of its field line mapping. This can explain unexpected behaviour in recent resistive-magnetohydrodynamic simulations of magnetic relaxation. Here we present the theory for a general periodic flux tube, representing, for example, a toroidal confinement device or a solar coronal loop. We show how an ideal dynamics on the side boundary of the tube implies that the sum of indices over all interior fixed points is invariant. This constraint applies to any continuous evolution inside the tube, which may be turbulent and/or dissipative. We also consider the analogous invariants obtained from periodic points (fixed points of the iterated mapping). Although there is a countably infinite family of invariants, we show that they lead to at most two independent dynamical constraints. The second constraint applies only in certain magnetic configurations. Several examples illustrate the theory.

The emergence of tilted bipolar active regions and the dispersal of their flux, mediated via processes such as diffusion, differential rotation and meridional circulation is believed to be responsible for the reversal of the Sun's polar field. This process (commonly known as the Babcock-Leighton mechanism) is usually modeled as a near-surface, spatially distributed $\alpha$-effect in kinematic mean-field dynamo models. However, this formulation leads to a relationship between polar field strength and meridional flow speed which is opposite to that suggested by physical insight and predicted by surface flux-transport simulations. With this in mind, we present an improved double-ring algorithm for modeling the Babcock-Leighton mechanism based on active region eruption, within the framework of an axisymmetric dynamo model. Using surface flux-transport simulations we first show that an axisymmetric formulation -- which is usually invoked in kinematic dynamo models -- can reasonably approximate the surface flux dynamics. Finally, we demonstrate that our treatment of the Babcock-Leighton mechanism through double-ring eruption leads to an inverse relationship between polar field strength and meridional flow speed as expected, reconciling the discrepancy between surface flux-transport simulations and kinematic dynamo models.

The final state of turbulent magnetic relaxation in a reversed field pinch is well explained by Taylor's hypothesis. However, recent resistive-magnetohydrodynamic simulations of the relaxation of braided solar coronal loops have led to relaxed fields far from the Taylor state, despite the conservation of helicity. We point out the existence of an additional topological invariant in any flux tube with non-zero field: the topological degree of the field line mapping. We conjecture that this constrains the relaxation, explaining why only one of three example simulations reaches the Taylor state.

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.

The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free "potential field" extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helicity due to flux emergence and large-scale photospheric motions. This helicity concentrates into twisted magnetic flux ropes, which may lose equilibrium and be ejected. Here, we consider how the magnetic structure predicted by this model-in particular the flux ropes-varies over the solar activity cycle, based on photospheric input data from six periods of cycle 23. The number of flux ropes doubles from minimum to maximum, following the total length of photospheric polarity inversion lines. However, the number of flux rope ejections increases by a factor of eight, following the emergence rate of active regions. This is broadly consistent with the observed cycle modulation of coronal mass ejections, although the actual rate of ejections in the simulation is about a fifth of the rate of observed events. The model predicts that, even at minimum, differential rotation will produce sheared, non-potential, magnetic structure at all latitudes.

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.