A Comparison of Sparse and Non-sparse Techniques for Electric-Field Inversion from Normal-Component Magnetograms

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

Chirality of high-latitude filaments over Solar Cycle 23

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

Kinematic active region formation in a three-dimensional solar dynamo model

We propose a phenomenological technique for modelling the emergence of active regions within a three-dimensional, kinematic dynamo framework. By imposing localized 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 three-dimensional 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

Optimal unstirred state of a passive scalar

Given a passive tracer distribution, what is the simplest unstirred pattern that may be reached under incompressible advection? This question is partially motivated by recent studies of three-dimensional (3-D) magnetic reconnection, in which the patterns of a topological invariant called the field line helicity greatly simplify until reaching a relaxed state. We test two approaches: a variational method with minimal constraints, and a magnetic relaxation scheme where the velocity is determined explicitly by the pattern of. Both methods achieve similar convergence for simple test cases. However, the magnetic relaxation method guarantees a monotonic decrease in the Dirichlet seminorm of, and is numerically more robust. We therefore apply the latter method to two complex mixed patterns modelled on the field line helicity of 3-D magnetic braids. The unstirring separates into a small number of large-scale regions determined by the initial topology, which is well preserved during the computation. Interestingly, the velocity field is found to have the same large-scale topology as. Similarity to the simplification found empirically in 3-D magnetic reconnection simulations supports the idea that advection is an important principle for field line helicity evolution.</p

How many active regions are necessary to predict the solar dipole moment?

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.Comment: 15 pages, 10 figures, accepted for publication in Ap

Pathways of Large-scale Magnetic Couplings Between Solar Coronal Events

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 timescale of the large-scale corona, estimated from the extreme-ultraviolet afterglow, is on average longer than the mean time between coronal mass ejections (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

Estimating the Rate of Field Line Braiding in the Solar Corona by Photospheric Flows

In this paper, we seek to understand the timescale in 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 (FTTE), 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 FTTE, the values obtained from the observations from the plage region are around an order of magnitude smaller. This behavior is carried over to the FTTE. Nevertheless, the results suggest that the photospheric motions induce complex tangling of the coronal field on a timescale of hours

A Model for Coronal Hole Bright Points and Jets Due to Moving Magnetic Elements

Coronal jets and bright points occur prolifically in predominantly unipolar magnetic regions, such as coronal holes (CHs), 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 "flyby" 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 polarity inversion line was strongly sheared, on the other hand, filament channels repeatedly formed and erupted via the breakout mechanism, explosively increasing the interchange reconnection and generating nonhelical jets. The flyby produced even more energetic and explosive jets. Our results explain several key aspects of CH bright points and jets, and the relationships between them

Surface Flux Transport on the Sun

We review the surface flux transport model for the evolution of magnetic flux patterns on the Sun’s surface. Our underlying motivation is to understand the model’s prediction of the polar field (or axial dipole) strength at the end of the solar cycle. The main focus is on the “classical” model: namely, steady axisymmetric profiles for differential rotation and meridional flow, and uniform supergranular diffusion. Nevertheless, the review concentrates on recent advances, notably in understanding the roles of transport parameters and – in particular – the source term. We also discuss the physical justification for the surface flux transport model, along with efforts to incorporate radial diffusion, and conclude by summarizing the main directions where researchers have moved beyond the classical model

Coronal Magnetic Field Evolution from 1996 to 2012: Continuous Non-potential Simulations

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