BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

A. R. Yeates

A. R. Yeates contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRphysics.plasm-phnlin.CDphysics.flu-dynphysics.space-phastro-phastro-ph.EPmath-phmath.MP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

26 published item(s)

preprint2022arXiv

On the limitations of magneto-frictional relaxation

The magneto-frictional method is used in solar physics to compute both static and quasi-static models of the Sun's coronal magnetic field. Here, we examine how accurately magneto-friction (without fluid pressure) is able to predict the relaxed state in a one-dimensional test case containing two magnetic null points. Firstly, we show that relaxation under the full ideal magnetohydrodynamic equations in the presence of nulls leads necessarily to a non-force-free state, which could not be reached exactly by magneto-friction. Secondly, the magneto-frictional solutions are shown to lead to breakdown of magnetic flux conservation, whether or not the friction coefficient is scaled with magnetic field strength. When this coefficient is constant, flux is initially conserved, but only until discontinuous current sheets form at the null points. In the ensuing weak solution, we show that magnetic flux is dissipated at these current sheets. The breakdown of flux conservation does not occur for an alternative viscous relaxation scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How good is the bipolar approximation of active regions for surface flux transport?

We investigate how representing active regions with bipolar magnetic regions (BMRs) affects the end-of-cycle polar field predicted by the surface flux transport model. Our study is based on a new database of BMRs derived from the SDO/HMI active region patch data between 2010 and 2020. An automated code is developed for fitting each active region patch with a BMR, matching both the magnetic flux and axial dipole moment of the region and removing repeat observations of the same region. By comparing the predicted evolution of each of the 1090 BMRs with the predicted evolution of their original active region patches, we show that the bipolar approximation leads to a 24% overestimate of the net axial dipole moment, given the same flow parameters. This is caused by neglecting the more complex multipolar and/or asymmetric magnetic structures of many of the real active regions, and may explain why previous flux transport models had to reduce BMR tilt angles to obtain realistic polar fields. Our BMR database and the Python code to extract it are freely available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Minimal Helicity of Solar Coronal Magnetic Fields

Potential field extrapolations are widely used as minimum-energy models for the Sun's coronal magnetic field. As the reference to which other magnetic fields are compared, they have -- by any reasonable definition -- no global (signed) magnetic helicity. Here we investigate the internal topological structure that is not captured by the global helicity integral, by splitting it into individual field line helicities. These are computed using potential field extrapolations from magnetogram observations over Solar Cycle 24, as well as for a simple illustrative model of a single bipolar region in a dipolar background. We find that localised patches of field line helicity arise primarily from linking between strong active regions and their overlying field, so that the total unsigned helicity correlates with the product of photospheric and open fluxes. Within each active region, positive and negative helicity may be unbalanced, but the signed helicity is only around a tenth of the unsigned helicity. Interestingly, in Cycle 24, there is a notable peak in unsigned helicity caused by a single large active region. On average, the total unsigned helicity at the resolution considered is approximately twice the typical signed helicity of a single real active region, according to non-potential models in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stellar Coronal Response to Differential Rotation and Flux Emergence

We perform a numerical parameter study to determine what effect varying differential rotation and flux emergence has on a star's non-potential coronal magnetic field. In particular we consider the effects on the star's surface magnetic flux, open magnetic flux, mean azimuthal field strength, coronal free magnetic energy, coronal heating and flux rope eruptions. To do this, we apply a magnetic flux transport model to describe the photospheric evolution, and couple this to the non-potential coronal evolution using a magnetofrictional technique. A flux emergence model is applied to add new magnetic flux onto the photosphere and into the corona. The parameters of this flux emergence model are derived from the solar flux emergence profile, however the rate of emergence can be increased to represent higher flux emergence rates than the Sun's. Overall we find that flux emergence has a greater effect on the non-potential coronal properties compared to differential rotation, with all the aforementioned properties increasing with increasing flux emergence rate. Although differential rotation has a lesser effect on the overall coronal properties compared to flux emergence, varying differential rotation does alter the coronal structure. As the differential rotation rate increases, the corona becomes more open, and more non-potential.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

By comparing a magneto-frictional model of the low coronal magnetic field to a potential-field source-surface model, we investigate the possible impact of non-potential magnetic structure on empirical solar-wind models. These empirical models (such as Wang-Sheeley-Arge) estimate the distribution of solar-wind speed solely from the magnetic-field structure in the low corona. Our models are computed in a domain between the solar surface and 2.5 solar radii, and are extended to 0.1 AU using a Schatten current-sheet model. The non-potential field has a more complex magnetic skeleton and quasi-separatrix structures than the potential field, leading to different sub-structure in the solar-wind speed proxies. It contains twisted magnetic structures that can perturb the separatrix surfaces traced down from the base of the heliospheric current sheet. A significant difference between the models is the greater amount of open magnetic flux in the non-potential model. Using existing empirical formulae this leads to higher predicted wind speeds for two reasons: partly because magnetic flux tubes expand less rapidly with height, but more importantly because more open field lines are further from coronal-hole boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Physical role of topological constraints in localised magnetic relaxation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Source of a Prominent Poleward Surge During Solar Cycle 24

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The coronal energy input from magnetic braiding

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A complete topological invariant for braided magnetic fields

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chirality of High Latitude Filaments over Solar Cycle 23

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lagrangian coherent structures in photospheric flows and their implications for coronal magnetic structure

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

In this review, our present day understanding of the Sun's global photospheric and coronal magnetic fields is discussed from both observational and theoretical viewpoints. Firstly, the large-scale properties of photospheric magnetic fields are described, along with recent advances in photospheric magnetic flux transport models. Following this, the wide variety of theoretical models used to simulate global coronal magnetic fields are described. From this, the combined application of both magnetic flux transport simulations and coronal modeling techniques to describe the phenomena of coronal holes, the Sun's open magnetic flux and the hemispheric pattern of solar filaments is discussed. Finally, recent advances in non-eruptive global MHD models are described. While the review focuses mainly on solar magnetic fields, recent advances in measuring and modeling stellar magnetic fields are described where appropriate. In the final section key areas of future research are identified.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unique Topological Characterization of Braided Magnetic Fields

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Generalized Flux Function for Three-dimensional Magnetic Reconnection

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical constraints from field line topology in magnetic flux tubes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Heating of braided coronal loops

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

Measurements of the interplanetary magnetic field (IMF) over several solar cycles do not agree with computed values of open magnetic flux from potential field extrapolations. The discrepancy becomes greater around solar maximum in each cycle, when the IMF can be twice as strong as predicted by the potential field model. Here we demonstrate that this discrepancy may be resolved by allowing for electric currents in the low corona (below 2.5 solar radii). We present a quasi-static numerical model of the large-scale coronal magnetic evolution, which systematically produces these currents through flux emergence and shearing by surface motions. The open flux is increased by 75%-85% at solar maximum, but only 25% at solar minimum, bringing it in line with estimates from IMF measurements. The additional open flux in the non-potential model arises through inflation of the magnetic field by electric currents, with super-imposed fluctuations due to coronal mass ejections. The latter are modelled by the self-consistent ejection of twisted magnetic flux ropes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solar Cycle Variation of Magnetic Flux Ropes in a Quasi-Static Coronal Evolution Model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Topological Constraints on Magnetic Relaxation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

The relative importance of different initiation mechanisms for coronal mass ejections (CMEs) on the Sun is uncertain. One possible mechanism is the loss of equilibrium of coronal magnetic flux ropes formed gradually by large-scale surface motions. In this paper, the locations of flux rope ejections in a recently-developed quasi-static global evolution model are compared with observed CME source locations over a 4.5-month period in 1999. Using EUV data, the low-coronal source locations are determined unambiguously for 98 out of 330 CMEs. Despite the incomplete observations, positive correlation (with coefficient up to 0.49) is found between the distributions of observed and simulated ejections, but only when binned into periods of one month or longer. This binning timescale corresponds to the time interval at which magnetogram data are assimilated into the coronal simulations, and the correlation arises primarily from the large-scale surface magnetic field distribution; only a weak dependence is found on the magnetic helicity imparted to the emerging active regions. The simulations are limited in two main ways: they produce fewer ejections, and they do not reproduce the strong clustering of observed CME sources into active regions. Due to this clustering, the horizontal gradient of radial photospheric magnetic field is better correlated with the observed CME source distribution (coefficient 0.67). Our results suggest that, while the gradual formation of magnetic flux ropes over weeks can account for many observed CMEs, especially at higher latitudes, there exists a second class of CMEs (at least half) for which dynamic active region flux emergence on shorter timescales must be the dominant factor.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Initiation of Coronal Mass Ejections in a Global Evolution Model

Loss of equilibrium of magnetic flux ropes is a leading candidate for the origin of solar coronal mass ejections (CMEs). The aim of this paper is to explore to what extent this mechanism can account for the initiation of CMEs in the global context. A simplified MHD model for the global coronal magnetic field evolution in response to flux emergence and shearing by large-scale surface motions is described and motivated. Using automated algorithms for detecting flux ropes and ejections in the global magnetic model, the effects of key simulation parameters on the formation of flux ropes and the number of ejections are considered, over a 177-day period in 1999. These key parameters include the magnitude and sign of magnetic helicity emerging in active regions, and coronal diffusion. The number of flux ropes found in the simulation at any one time fluctuates between about 28 and 48, sustained by the emergence of new bipolar regions, but with no systematic dependence on the helicity of these regions. However, the emerging helicity does affect the rate of flux rope ejections, which doubles from 0.67 per day if the bipoles emerge untwisted to 1.28 per day in the run with greatest emerging twist. The number of ejections in the simulation is also increased by 20%-30% by choosing the majority sign of emerging bipole helicity in each hemisphere, or by halving the turbulent diffusivity in the corona. For reasonable parameter choices, the model produces approximately 50% of the observed CME rate. This indicates that the formation and loss of equilibrium of flux ropes may be a key element in explaining a significant fraction of observed CMEs.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

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

The hemispheric pattern of solar filaments is considered using newly-developed simulations of the real photospheric and 3D coronal magnetic fields over a 6-month period, on a global scale. The magnetic field direction in the simulation is compared directly with the chirality of observed filaments, at their observed locations. In our model the coronal field evolves through a continuous sequence of nonlinear force-free equilibria, in response to the changing photospheric boundary conditions and the emergence of new magnetic flux. In total 119 magnetic bipoles with properties matching observed active regions are inserted. These bipoles emerge twisted and inject magnetic helicity into the solar atmosphere. When we choose the sign of this active-region helicity to match that observed in each hemisphere, the model produces the correct chirality for up to 96% of filaments, including exceptions to the hemispheric pattern. If the emerging bipoles have zero helicity, or helicity of the opposite sign, then this percentage is much reduced. In addition, the simulation produces a higher proportion of filaments with the correct chirality after longer times. This indicates that a key element in the evolution of the coronal field is its long-term memory, and the build-up and transport of helicity from low to high latitudes over many months. It highlights the importance of continuous evolution of the coronal field, rather than independent extrapolations at different times. This has significant consequences for future modelling such as that related to the origin and development of coronal mass ejections.

0 citations0 reviewsNo code linkedOpen access