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

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

On the helicity of open magnetic fields

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. However, 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

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

On the compatibility of a flux transport dynamo with a fast tachocline scenario

The compatibility of the fast tachocline scenario with a flux transport dynamo model is explored. We employ a flux transport dynamo model coupled with simple feedback formulae relating the thickness of the tachocline to the amplitude of the magnetic field or to the Maxwell stress. The dynamo model is found to be robust against the nonlinearity introduced by this simplified fast tachocline mechanism. Solar-like butterfly diagrams are found to persist and, even without any parameter fitting, the overall thickness of the tachocline is well within the range admitted by helioseismic constraints. In the most realistic case of a time and latitude dependent tachocline thickness linked to the value of the Maxwell stress, both the thickness and its latitude dependence are in excellent agreement with seismic results. In the nonparametric models, cycle related temporal variations in tachocline thickness are somewhat larger than admitted by helioseismic constraints; we find, however, that introducing a further parameter into our feedback formula readily allows further fine tuning of the thickness variations.Comment: Accepted in Solar Physic

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

Solar Magnetic Carpet I: Simulation of Synthetic Magnetograms

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.Comment: 33 pages, 16 figures, 5 gif movies included: movies may be viewed at http://www-solar.mcs.st-and.ac.uk/~karen/movies_paper1

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