In order to efficiently analyse the vast amount of data generated by solar space missions and ground-base instruments, modern machine learning techniques such as decision trees, support vector machines (SVMs) and neural networks can be very useful. In this paper we present initial results from using a convolutional neural network (CNN) to analyse observations from the Atmospheric Imaging Assembly (AIA) in the 1600A wavelength. The data is pre-processed to locate flaring regions where flare ribbons are visible in the observations. The CNN is created and trained to automatically analyse the shape and position of the flare ribbons, by identifying whether each image belongs into one of four classes: two-ribbon flare, compact/circular ribbon flare, limb flare or quiet Sun, with the final class acting as a control for any data included in the training or test sets where flaring regions are not present. The network created can classify flare ribbon observations into any of the four classes with a final accuracy of 94%. Initial results show that most of the images are correctly classified with the compact flare class being the only class where accuracy drops below 90% ad some observations are wrongly classified as belonging to the limb class.

Understanding how the release of stored magnetic energy contributes to the generation of non-thermal high energy particles during solar flares is an important open problem in solar physics. Magnetic reconnection plays a fundamental role in the energy release and conversion processes taking place during flares. A common approach for investigating particle acceleration is to use test particles in fields derived from magnetohydrodynamic (MHD) simulations of reconnection. These MHD simulations use anomalous resistivities that are much larger than the Spitzer resistivity based on Coulomb collisions. The processes leading to enhanced resistivity should also affect the test particles. We explore the link between resistivity and particle orbits building on a previous study using a 2D MHD simulation of magnetic reconnection. This paper extends the previous investigation to a 3D magnetic reconnection configuration and to study the effect on test particle orbits. We carried out orbit calculations using a 3D MHD simulation of separator reconnection. We use the relativistic guiding centre approximation including stochastic pitch angle scattering. The effects of varying the resistivity and the models for pitch angle scattering on particle orbit trajectories, final positions, energy spectra, final pitch angle distribution, and orbit duration are all studied in detail. Pitch angle scattering widens collimated beams of orbit trajectories, allowing orbits to access previously unaccessible field lines; this causes final positions to spread to topological structures that were previously inaccessible. Scattered orbit energy spectra are found to be predominantly affected by the level of anomalous resistivity, with the pitch angle scattering model only playing a role in isolated cases. Scattering is found to play a crucial role in determining the pitch angle and orbit duration distributions.

Sudden destabilisations of the magnetic field, such as those caused by spontaneous reconnection, will produce waves and/or flows. Here, we investigate the nature of the plasma motions resulting from spontaneous reconnection at a 3D separator. In order to clearly see these perturbations, we start from a magnetohydrostatic equilibrium containing two oppositely-signed null points joined by a generic separator along which lies a twisted current layer. The nature of the magnetic reconnection initiated in this equilibrium as a result of an anomalous diffusivity is discussed in detail in \cite{Stevenson15_jgra}. The resulting sudden loss of force balance inevitably generates waves that propagate away from the diffusion region carrying the dissipated current. In their wake a twisting stagnation-flow, in planes perpendicular to the separator, feeds flux back into the original diffusion site (the separator) in order to try to regain equilibrium. This flow drives a phase of slow weak impulsive-bursty reconnection that follows on after the initial fast-reconnection phase.

Magnetic separators, which lie on the boundary between four topologically-distinct flux domains, are prime locations in three-dimensional magnetic fields for reconnection, especially in the magnetosphere between the planetary and interplanetary magnetic field and also in the solar atmosphere. Little is known about the details of separator reconnection and so the aim of this paper, which is the first of two, is to study the properties of magnetic reconnection at a single separator. Three-dimensional, resistive magnetohydrodynamic numerical experiments are run to study separator reconnection starting from a magnetohydrostatic equilibrium which contains a twisted current layer along a single separator linking a pair of opposite-polarity null points. The resulting reconnection occurs in two phases. The first is short involving rapid-reconnection in which the current at the separator is reduced by a factor of around 2.3. Most ($75\%$) of the magnetic energy is converted during this phase, via Ohmic dissipation, directly into internal energy, with just $0.1\%$ going into kinetic energy. During this phase the reconnection occurs along most of the separator away from its ends (the nulls), but in an asymmetric manner which changes both spatially and temporally over time. The second phase is much longer and involves slow impulsive-bursty reconnection. Again Ohmic heating dominates over viscous damping. Here, the reconnection occurs in small localised bursts at random anywhere along the separator.

The aim of this work is to investigate and characterise particle behaviour in an (observationally-driven) 3D MHD model of the solar atmosphere above a slowly evolving, non-flaring active region. We use a relativistic guiding-centre particle code to investigate the behaviour of selected particle orbits, distributed throughout a single snapshot of the 3D MHD simulation. Two distinct particle acceleration behaviours are recovered, which affect both electrons and protons: (i) direct acceleration along field lines and (ii) tangential drifting of guiding centres with respect to local magnetic field. However, up to 40\% of all particles actually experience a form of (high energy) particle trap, because of changes in the direction of the electric field and unrelated to the strength of the magnetic field; such particles are included in the first category. Additionally, category (i) electron and proton orbits undergo surprisingly strong acceleration to non-thermal energies ($\lesssim42$MeV), because of the strength and extent of super-Dreicer electric fields created by the MHD simulation. Reducing the electric field strength of the MHD model does not significantly affect the efficiency of the (electric field-based) trapping mechanism, but does reduce the peak energies gained by orbits. We discuss the implications for future experiments, which aim to simulate non-flaring active region heating and reconnection.

The aim of this work is to investigate and characterise particle behaviour in a 3D magnetohydrodynamic (MHD) model of a reconnecting magnetic separator. We use a relativistic guiding-centre test-particle code to investigate electron and proton acceleration in snapshots from 3D MHD separator reconnection experiments, and compare the results with findings from an analytical separator reconnection model studied in a previous investigation. The behaviour and acceleration of large distributions of particles are examined in detail for both analytical and numerical separator reconnection models. Accelerated particle orbit trajectories are shown to follow the separator before leaving the system along the separatrix surface of one of the nulls (determined by particle species) in the system of both models. A sufficiently localised electric field about the separator causes the orbits to appear to follow the spine bounding the separatrix surface field lines instead. We analyse and discuss the locations and spread of accelerated particle orbit final positions, which are seen to change over time in the numerical separator reconnection model. We deduce a simple relationship between the final energy range of particle orbits and the model dimensions, and discuss its implications for observed magnetic separators in the solar corona.

While theoretical models and simulations of magnetic reconnection often assume symmetry such that the magnetic null point when present is co-located with a flow stagnation point, the introduction of asymmetry typically leads to non-ideal flows across the null point. To understand this behavior, we present exact expressions for the motion of three-dimensional linear null points. The most general expression shows that linear null points move in the direction along which the vector field and its time derivative are antiparallel. Null point motion in resistive magnetohydrodynamics results from advection by the bulk plasma flow and resistive diffusion of the magnetic field, which allows non-ideal flows across topological boundaries. Null point motion is described intrinsically by parameters evaluated locally; however, global dynamics help set the local conditions at the null point. During a bifurcation of a degenerate null point into a null-null pair or the reverse, the instantaneous velocity of separation or convergence of the null-null pair will typically be infinite along the null space of the Jacobian matrix of the magnetic field, but with finite components in the directions orthogonal to the null space. Not all bifurcating null-null pairs are connected by a separator. Furthermore, except under special circumstances, there will not exist a straight line separator connecting a bifurcating null-null pair. The motion of separators cannot be described using solely local parameters, because the identification of a particular field line as a separator may change as a result of non-ideal behavior elsewhere along the field line.

Magnetic reconnection, a change of magnetic field connectivity, is a fundamental physical process in which magnetic energy is released explosively. It is responsible for various eruptive phenomena in the universe. However, this process is difficult to observe directly. Here, the magnetic topology associated with a solar reconnection event is studied in three dimensions (3D) using the combined perspectives of two spacecraft. The sequence of extreme ultraviolet (EUV) images clearly shows that two groups of oppositely directed and non-coplanar magnetic loops gradually approach each other, forming a separator or quasi-separator and then reconnecting. The plasma near the reconnection site is subsequently heated from $\sim$1 to $\ge$5 MK. Shortly afterwards, warm flare loops ($\sim$3 MK) appear underneath the hot plasma. Other observational signatures of reconnection, including plasma inflows and downflows, are unambiguously revealed and quantitatively measured. These observations provide direct evidence of magnetic reconnection in a 3D configuration and reveal its origin.

Magnetic fields permeate the entire solar atmosphere weaving an extremely complex pattern on both local and global scales. In order to understand the nature of this tangled web of magnetic fields, its magnetic skeleton, which forms the boundaries between topologically distinct flux domains, may be determined. The magnetic skeleton consists of null points, separatrix surfaces, spines and separators. The skeleton is often used to clearly visualize key elements of the magnetic configuration, but parts of the skeleton are also locations where currents and waves may collect and dissipate. In this review, the nature of the magnetic skeleton on both global and local scales, over solar cycle time scales, is explained. The behaviour of wave pulses in the vicinity of both nulls and separators is discussed and so too is the formation of current layers and reconnection at the same features. Each of these processes leads to heating of the solar atmosphere, but collectively do they provide enough heat, spread over a wide enough area, to explain the energy losses throughout the solar atmosphere? Here, we consider this question for the three different solar regions: active regions, open-field regions and the quiet Sun. We find that the heating of active regions and open-field regions is highly unlikely to be due to reconnection or wave dissipation at topological features, but it is possible that these may play a role in the heating of the quiet Sun. In active regions, the absence of a complex topology may play an important role in allowing large energies to build up and then, subsequently, be explosively released in the form of a solar flare. Additionally, knowledge of the intricate boundaries of open-field regions (which the magnetic skeleton provides) could be very important in determining the main acceleration mechanism(s) of the solar wind.

Separators, which are in many ways the three-dimensional equivalent to two-dimensional nulls, are important sites for magnetic reconnection. Magnetic reconnection occurs in strong current layers which have very short length scales. The aim of this work is to explore the nature of current layers around separators. A separator is a special field line which lies along the intersection of two separatrix surfaces and forms the boundary between four topologically distinct flux domains. In particular, here the current layer about a separator that joins two 3D nulls and lies along the intersection of their separatrix surfaces is investigated. A magnetic configuration containing a single separator embedded in a uniform plasma with a uniform electric current parallel to the separator is considered. This initial magnetic setup, which is not in equilibrium, relaxes in a non-resistive manner to form an equilibrium. The relaxation is achieved using the 3D MHD code, Lare3d, with resistivity set to zero. A series of experiments with varying initial current are run to investigate the characteristics of the resulting current layers present in the final (quasi-) equilibrium states. In each experiment, the separator collapses and a current layer forms along it. The dimensions and strength of the current layer increase with initial current. It is found that separator current layers formed from current parallel to the separator are twisted. Also the collapse of the separator is a process that evolves like an infinite-time singularity where the length, width and peak current in the layer grow slowly whilst the depth of the current layer decreases.

While the exact acceleration mechanism of energetic particles during solar flares is (as yet) unknown, magnetic reconnection plays a key role both in the release of stored magnetic energy of the solar corona and the magnetic restructuring during a flare. Recent work has shown that special field lines, called separators, are common sites of reconnection in 3D numerical experiments. To date, 3D separator reconnection sites have received little attention as particle accelerators. We investigate the effectiveness of separator reconnection as a particle acceleration mechanism for electrons and protons. We study the particle acceleration using a relativistic guiding-centre particle code in a time-dependent kinematic model of magnetic reconnection at a separator. The effect upon particle behaviour of initial position, pitch angle and initial kinetic energy are examined in detail, both for specific (single) particle examples and for large distributions of initial conditions. The separator reconnection model contains several free parameters and we study the effect of changing these parameters upon particle acceleration, in particular in view of the final particle energy ranges which agree with observed energy spectra.

The distribution of magnetic flux across the solar photosphere results in a complex web of coronal magnetic field structures. To understand this complexity, the magnetic skeleton of the coronal field can be calculated. The skeleton highlights the separatrix surfaces that divide the field into topologically distinct regions, allowing open-field regions on the solar surface to be located. Furthermore, separatrix surfaces and their intersections (separators) are important likely energy release sites. This paper investigates, throughout the solar cycle, the nature of coronal magnetic-field topologies that arise under the potential-field source-surface approximation. In particular, we characterise the typical global fields at solar maximum and minimum. Global magnetic fields are extrapolated from observed Kitt Peak and SOLIS synoptic magnetograms, from Carrington rotations 1645 to 2144, using the potential-field source-surface model. Hence, variations in the coronal skeleton are studied over three solar cycles. Key building blocks which make up magnetic fields are identified and classified according to the nature of their separatrix surfaces. The magnetic skeleton reveals that, at solar maximum, the global coronal field involves a multitude of topological structures at all latitudes. Many open-field regions exist originating anywhere on the photosphere. At solar minimum, the coronal topology is heavily influenced by the solar dipole. A strong dipole results in a simple large-scale structure involving just two large polar open-field regions, but, at short radial distances between plus or minus 60 deg latitude, the small-scale topology is complex. If the solar dipole if weak, as in the recent minimum, then the low-latitude quiet-sun magnetic fields may be globally significant enough to create many disconnected open-field regions at low latitudes, in addition to the two polar open-field regions.

In this paper we study current accumulations in 3D "tilted" nulls formed by a folding of the spine and fan. A non-zero component of current parallel to the fan is required such that the null's fan plane and spine are not perpendicular. Our aims are to provide valid magnetohydrostatic equilibria and to describe the current accumulations in various cases involving finite plasma pressure.To create our equilibrium current structures we use a full, non-resistive, magnetohydrodynamic (MHD) code so that no reconnection is allowed. A series of experiments are performed in which a perturbed 3D tilted null relaxes towards an equilibrium via real, viscous damping forces. Changes to the initial plasma pressure and to magnetic parameters are investigated systematically.An initially tilted fan is associated with a non-zero Lorentz force that drives the fan and spine to collapse towards each other, in a similar manner to the collapse of a 2D X-point. In the final equilibrium state for an initially radial null with only the current perpendicular to the spine, the current concentrates along the tilt axis of the fan and in a layer about the null point with a sharp peak at the null itself. The continued growth of this peak indicates that the system is in an asymptotic regime involving an infinite time singularity at the null. When the initial tilt disturbance (current perpendicular to the spine) is combined with a spiral-type disturbance (current parallel to the spine), the final current density concentrates in three regions: one on the fan along its tilt axis and two around the spine, above and below the fan. The increased area of current accumulation leads to a weakening of the singularity formed at the null. The 3D spine-fan collapse with generic current studied here provides the ideal setup for non-steady reconnection studies.

We provide a valid magnetohydrostatic equilibrium from the collapse of a 2D X-point in the presence of a finite plasma pressure, in which the current density is not simply concentrated in an infinitesimally thin, one-dimensional current sheet, as found in force-free solutions. In particular, we wish to determine if a finite pressure current sheet will still involve a singular current, and if so, what is the nature of the singularity. We use a full MHD code, with the resistivity set to zero, so that reconnection is not allowed, to run a series of experiments in which an X-point is perturbed and then is allowed to relax towards an equilibrium, via real, viscous damping forces. Changes to the magnitude of the perturbation and the initial plasma pressure are investigated systematically. The final state found in our experiments is a "quasi-static" equilibrium where the viscous relaxation has completely ended, but the peak current density at the null increases very slowly following an asymptotic regime towards an infinite time singularity. Using a high grid resolution allows us to resolve the current structures in this state both in width and length. In comparison with the well known pressureless studies, the system does not evolve towards a thin current sheet, but concentrates the current at the null and the separatrices. The growth rate of the singularity is found to be tD, with 0 < D < 1. This rate depends directly on the initial plasma pressure, and decreases as the pressure is increased. At the end of our study, we present an analytical description of the system in a quasi-static non-singular equilibrium at a given time, in which a finite thick current layer has formed at the null.

The removal of magnetic flux from the quiet-sun photosphere is important for maintaining the statistical steady-state of the magnetic field there, for determining the magnetic flux budget of the Sun, and for estimating the rate of energy injected into the upper solar atmosphere. Magnetic feature death is a measurable proxy for the removal of detectable flux. We used the SWAMIS feature tracking code to understand how nearly 20000 detected magnetic features die in an hour-long sequence of Hinode/SOT/NFI magnetograms of a region of quiet Sun. Of the feature deaths that remove visible magnetic flux from the photosphere, the vast majority do so by a process that merely disperses the previously-detected flux so that it is too small and too weak to be detected. The behavior of the ensemble average of these dispersals is not consistent with a model of simple planar diffusion, suggesting that the dispersal is constrained by the evolving photospheric velocity field. We introduce the concept of the partial lifetime of magnetic features, and show that the partial lifetime due to Cancellation of magnetic flux, 22 h, is 3 times slower than previous measurements of the flux turnover time. This indicates that prior feature-based estimates of the flux replacement time may be too short, in contrast with the tendency for this quantity to decrease as resolution and instrumentation have improved. This suggests that dispersal of flux to smaller scales is more important for the replacement of magnetic fields in the quiet Sun than observed bipolar cancellation. We conclude that processes on spatial scales smaller than those visible to Hinode dominate the processes of flux emergence and cancellation, and therefore also the quantity of magnetic flux that threads the photosphere.

This paper is the third in a series of papers working towards the construction of a realistic, evolving, non-linear force-free coronal field model for the solar magnetic carpet. Here, we present preliminary results of 3D time-dependent simulations of the small-scale coronal field of the magnetic carpet. Four simulations are considered, each with the same evolving photospheric boundary condition: a 48 hr time series of synthetic magnetograms produced from the model of Meyer, Mackay, van Ballegooijen and Parnell, 2011, Solar Phys., 272, 29. Three simulations include a uniform, overlying coronal magnetic field of differing strength, the fourth simulation includes no overlying field. The build-up, storage and dissipation of magnetic energy within the simulations is studied. In particular, we study their dependence upon the evolution of the photospheric magnetic field and the strength of the overlying coronal field. We also consider where energy is stored and dissipated within the coronal field. The free magnetic energy built up is found to be more than sufficient to power small-scale, transient phenomena such as nanoflares and X-ray bright points, with the bulk of the free energy found to be stored low down, between 0.5-0.8 Mm. The energy dissipated is presently found to be too small to account for the heating of the entire quiet Sun corona. However, the form and location of energy dissipation regions are in qualitative agreement with what is observed on small scales on the Sun. Future MHD modelling using the same synthetic magnetograms may lead to a higher energy release.

Determining the heating mechanism (or mechanisms) that causes the outer atmosphere of the Sun, and many other stars, to reach temperatures orders of magnitude higher than their surface temperatures has long been a key problem. For decades the problem has been known as the coronal heating problem, but it is now clear that `coronal heating' cannot be treated or explained in isolation and that the heating of the whole solar atmosphere must be studied as a highly coupled system. The magnetic field of the star is known to play a key role, but, despite significant advancements in solar telescopes, computing power and much greater understanding of theoretical mechanisms, the question of which mechanism or mechanisms are the dominant supplier of energy to the chromosphere and corona is still open. Following substantial recent progress, we consider the most likely contenders and discuss the key factors that have made, and still make, determining the actual (coronal) heating mechanism (or mechanisms) so difficult.

Context: The majority of studies on stressed 3D magnetic null points consider magnetic reconnection driven by an external perturbation, but the formation of a genuine current sheet equilibrium remains poorly understood. This problem has been considered more extensively in two-dimensions, but lacks a generalization into 3D fields. Aims: 3D magnetic nulls are more complex than 2D nulls and the field can take a greater range of magnetic geometries local to the null. Here, we focus on one type and consider the dynamical non-resistive relaxation of 3D spiral nulls with initial spine-aligned current. We aim to provide a valid magnetohydrostatic equilibrium, and describe the electric current accumulations in various cases, involving a finite plasma pressure. Methods: A full MHD code is used, with the resistivity set to zero so that reconnection is not allowed, to run a series of experiments in which a perturbed spiral 3D null point is allowed to relax towards an equilibrium, via real, viscous damping forces. Changes to the initial plasma pressure and other magnetic parameters are investigated systematically. Results: For the axi-symmetric case, the evolution of the field and the plasma is such that it concentrates the current density in two cone-shaped regions along the spine, thus concentrating the twist of the magnetic field around the spine, leaving a radial configuration in the fan plane. The plasma pressure redistributes in order to maintain the current density accumulations. However, it is found that changes in the initial plasma pressure do not modify the final state significantly. In the cases where the initial magnetic field is not axi-symmetric, a infinite-time singularity of current perpendicular to the fan is found at the location of the null.

The highly dynamical, complex nature of the solar atmosphere naturally implies the presence of waves in a topologically varied magnetic environment. Here, the interaction of waves with topological features such as null points is inevitable and potentially important for energetics. The low resistivity of the solar coronal plasma implies that non-MHD effects should be considered in studies of magnetic energy release in this environment. This paper investigates the role of the Hall term in the propagation and dissipation of waves, their interaction with 2D magnetic X-points and the nature of the resulting reconnection. A Lagrangian remap shock-capturing code (Lare2d) is used to study the evolution of an initial fast magnetoacoustic wave annulus for a range of values of the ion skin depth in resistive Hall MHD. A magnetic null-point finding algorithm is also used to locate and track the evolution of the multiple null-points that are formed in the system. Depending on the ratio of ion skin depth to system size, our model demonstrates that Hall effects can play a key role in the wave-null interaction. In particular, the initial fast-wave pulse now consists of whistler and ion-cyclotron components; the dispersive nature of the whistler wave leads to (i) earlier interaction with the null, (ii) the creation of multiple additional, transient nulls and, hence, an increased number of energy release sites. In the Hall regime, the relevant timescales (such as the onset of reconnection and the period of the oscillatory relaxation) of the system are reduced significantly, and the reconnection rate is enhanced.

The sudden reconnection of a non-force free 2D current layer, embedded in a low-beta plasma, triggered by the onset of an anomalous resistivity, is studied in detail. The resulting behaviour consists of two main phases. Firstly, a transient reconnection phase, in which the current in the layer is rapidly dispersed and some flux is reconnected. This dispersal of current launches a family of small amplitude magnetic and plasma perturbations, which propagate away from the null at the local fast and slow magnetosonic speeds. The vast majority of the magnetic energy released in this phase goes into internal energy of the plasma, and only a tiny amount is converted into kinetic energy. In the wake of the outwards propagating pulses, an imbalance of Lorentz and pressure forces creates a stagnation flow which drives a regime of impulsive bursty reconnection, in which fast reconnection is turned on and off in a turbulent manner as the current density exceeds and falls below a critical value. During this phase, the null current density is continuously built up above a certain critical level, then dissipated very rapidly, and built up again, in a stochastic manner. Interestingly, the magnetic energy converted during this quasi-steady phase is greater than that converted during the initial transient reconnection phase. Again essentially all the energy converted during this phase goes directly to internal energy. These results are of potential importance for solar flares and coronal heating, and set a conceptually important reference for future 3D studies.

Magnetic neutral points, where the magnitude of the magnetic field vanishes locally, are potential locations for energy conversion in the solar corona. The fact that the magnetic field is identically zero at these points suggests that for the study of current sheet formation and of any subsequent resistive dissipation phase, a finite beta plasma should be considered, rather than neglecting the plasma pressure as has often been the case in the past. The rapid dissipation of a finite current layer in non-force-free equilibrium is investigated numerically, after the sudden onset of an anomalous resistivity. The aim of this study is to determine how the energy is redistributed during the initial diffusion phase, and what is the nature of the outward transmission of information and energy. The resistivity rapidly diffuses the current at the null point. The presence of a plasma pressure allows the vast majority of the free energy to be transferred into internal energy. Most of the converted energy is used in direct heating of the surrounding plasma, and only about 3% is converted into kinetic energy, causing a perturbation in the magnetic field and the plasma which propagates away from the null at the local fast magnetoacoustic speed. The propagating pulses show a complex structure due to the highly non-uniform initial state. It is shown that this perturbation carries no net current as it propagates away from the null. The fact that, under the assumptions taken in this paper, most of the magnetic energy released in the reconnection converts internal energy of the plasma, may be highly important for the chromospheric and coronal heating problem.

Context. For the last thirty years, most of the studies on the relaxation of stressed magnetic fields in the solar environment have onlyconsidered the Lorentz force, neglecting plasma contributions, and therefore, limiting every equilibrium to that of a force-free field. Aims. Here we begin a study of the non-resistive evolution of finite beta plasmas and their relaxation to magnetohydrostatic states, where magnetic forces are balanced by plasma-pressure gradients, by using a simple 2D scenario involving a hydromagnetic disturbance to a uniform magnetic field. The final equilibrium state is predicted as a function of the initial disturbances, with aims to demonstrate what happens to the plasma during the relaxation process and to see what effects it has on the final equilibrium state. Methods. A set of numerical experiments are run using a full MHD code, with the relaxation driven by magnetoacoustic waves damped by viscous effects. The numerical results are compared with analytical calculations made within the linear regime, in which the whole process must remain adiabatic. Particular attention is paid to the thermodynamic behaviour of the plasma during the relaxation. Results. The analytical predictions for the final non force-free equilibrium depend only on the initial perturbations and the total pressure of the system. It is found that these predictions hold surprisingly well even for amplitudes of the perturbation far outside the linear regime. Conclusions. Including the effects of a finite plasma beta in relaxation experiments leads to significant differences from the force-free case.

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.

The role of null-point reconnection in a 3D numerical MHD model of solar emerging flux is investigated. The model consists of a twisted magnetic flux tube rising through a stratified convection zone and atmosphere to interact and reconnect with a horizontal overlying magnetic field in the atmosphere. Null points appear as the reconnection begins and persist throughout the rest of the emergence, where they can be found mostly in the model photosphere and transition region, forming two loose clusters on either side of the emerging flux tube. Up to 26 nulls are present at any one time, and tracking in time shows that there is a total of 305 overall, despite the initial simplicity of the magnetic field configuration. We find evidence for the reality of the nulls in terms of their methods of creation and destruction, their balance of signs, their long lifetimes, and their geometrical stability. We then show that due to the low parallel electric fields associated with the nulls, null-point reconnection is not the main type of magnetic reconnection involved in the interaction of the newly emerged flux with the overlying field. However, the large number of nulls implies that the topological structure of the magnetic field must be very complex and the importance of reconnection along separators or separatrix surfaces for flux emergence cannot be ruled out.

The importance of magnetic reconnection as an energy release mechanism in many solar, stellar, magnetospheric and astrophysical phenomena has long been recognised. Reconnection is the only mechanism by which magnetic fields can globally restructure, enabling them to access a lower energy state. Over the past decade, there have been some major advances in our understanding of three-dimensional reconnection. In particular, the key characteristics of 3D magnetohydrodynamic (MHD) reconnection have been determined. For instance, 3D reconnection (i) occurs with or without nulls, (ii) occurs continuously and continually throughout a diffusion region and (iii) is driven by counter rotating flows. Furthermore, analysis of resistive 3D MHD magnetic experiments have revealed some intriguing effects relating to where and how reconnection occurs. To illustrate these new features, a series of constant-resistivity experiments, involving the interaction of two opposite-polarity magnetic sources in an overlying field, are considered. Such a simple interaction represents a typical building block of the Sun's magnetic atmosphere. By following the evolution of the magnetic topology, we are able to explain where, how and at what rate the reconnection occurs. Remarkably there can be up to five energy release sites at anyone time (compared to one in the potential case) and the duration of the interaction increases (more than doubles) as the resistivity decreases (by a factor of 16). The decreased resistivity also leads to a higher peak ohmic dissipation and more energy being released in total, as a result of a greater injection of Poynting flux.

Magnetic null points can be located numerically in a potential field extrapolation or their average density can be estimated from the Fourier spectrum of a magnetogram. We use both methods to compute the null point density from a quiet Sun magnetogram made with Hinode's NFI and from magnetograms from SOHO's MDI in both its high-resolution and low-resolution modes. All estimates of the super-chromospheric column density (z>1.5 Mm) agree with one another and with the previous measurements: 0.003 null points per square Mm of solar surface.

The coronal magnetic field above a particular photospheric region will vanish at a certain number of points, called null points. These points can be found directly in a potential field extrapolation or their density can be estimated from Fourier spectrum of the magnetogram. The spectral estimate, which assumes that the extrapolated field is random, homogeneous and has Gaussian statistics, is found here to be relatively accurate for quiet Sun magnetograms from SOHO's MDI. The majority of null points occur at low altitudes, and their distribution is dictated by high wavenumbers in the Fourier spectrum. This portion of the spectrum is affected by Poisson noise, and as many as five-sixths of null points identified from a direct extrapolation can be attributed to noise. The null distribution above 1500 km is found to depend on wavelengths that are reliably measured by MDI in either its low-resolution or high-resolution mode. After correcting the spectrum to remove white noise and compensate for the modulation transfer function we find that a potential field extrapolation contains, on average, one magnetic null point, with altitude greater than 1.5 Mm, above every 322 square Mm patch of quiet Sun. Analysis of 562 quiet Sun magnetograms spanning the two latest solar minimum shows that the null point density is relatively constant with roughly 10% day-to-day variation. At heights above 1.5 Mm, the null point density decreases approximately as the inverse cube of height. The photospheric field in the quiet Sun is well approximated as that from discrete elements with mean flux 1.0e19 Mx distributed randomly with density n=0.007 per square Mm.

With the recent launch of the Hinode satellite our view of the nature and evolution of quiet-Sun regions has been improved. In light of the new high resolution observations, we revisit the study of the quiet Sun's topological nature. Topology is a tool to explain the complexity of the magnetic field, the occurrence of reconnection processes, and the heating of the corona. This Letter aims to give new insights to these different topics. Using a high-resolution Hinode/SOT observation of the line-of-sight magnetic field on the photosphere, we calculate the three dimensional magnetic field in the region above assuming a potential field. From the 3D field, we determine the existence of null points in the magnetic configuration. From this model of a continuous field, we find that the distribution of null points with height is significantly different from that reported in previous studies. In particular, the null points are mainly located above the bottom boundary layer in the photosphere (54%) and in the chromosphere (44%) with only a few null points in the corona (2%). The density of null points (expressed as the ratio of the number of null points to the number of photospheric magnetic fragments) in the solar atmosphere is estimated to be between 3% and 8% depending on the method used to identify the number of magnetic fragments in the observed photosphere. This study reveals that the heating of the corona by magnetic reconnection at coronal null points is unlikely. Our findings do not rule out the heating of the corona at other topological features. We also report the topological complexity of the chromosphere as strongly suggested by recent observations from Hinode/SOT.

Null points are important locations in vector fields, such as a magnetic field. A new technique (a trilinear method for finding null points) is presented for finding null points over a large grid of points, such as those derived from a numerical experiment. The method was designed so that the null points found would agree with any fieldlines traced using the commonly used trilinear interpolation. It is split into three parts: reduction, analysis and positioning, which, when combined, provide an efficient means of locating null points to a user-defined sub-grid accuracy. We compare the results of the trilinear method with that of a method based on the Poincare index, and discuss the accuracy and limitations of both methods.

Feature tracking and recognition are increasingly common tools for data analysis, but are typically implemented on an ad-hoc basis by individual research groups, limiting the usefulness of derived results when selection effects and algorithmic differences are not controlled. Specific results that are affected include the solar magnetic turnover time, the distributions of sizes, strengths, and lifetimes of magnetic features, and the physics of both small scale flux emergence and the small-scale dynamo. In this paper, we present the results of a detailed comparison between four tracking codes applied to a single set of data from SOHO/MDI, describe the interplay between desired tracking behavior and parameterization of tracking algorithms, and make recommendations for feature selection and tracking practice in future work.

The heating of the solar corona is likely to be due to reconnection of the highly complex magnetic field that threads throughout its volume. We have run a numerical experiment of an elementary interaction between the magnetic field of two photospheric sources in an overlying field that represents a fundamental building block of the coronal heating process. The key to explaining where, how and how much energy is released during such an interaction is to calculate the resulting evolution of the magnetic skeleton. A skeleton is essentially the web of magnetic flux surfaces (called separatrix surfaces) that separate the coronal volume into topologically distinct parts. For the first time the skeleton of the magnetic field in a 3D numerical MHD experiment is calculated and carefully analysed, as are the ways in which it bifurcates into different topologies. A change in topology normally changes the number of magnetic reconnection sites. In our experiment, the magnetic field evolves through a total of six distinct topologies. Initially, no magnetic flux joins the two sources. Then a new type of bifurcation, called a global double-separator bifurcation, takes place: this bifurcation is likely to be one of the main ways in which new separators are created in the corona (separators are field lines at which 3D reconnection takes place). This is the first of five bifurcations in which the skeleton becomes progressively more complex before simplifying. Surprisingly, for such a simple initial state, at the peak of complexity there are five separators and eight flux domains present.

The magnetic field plays a crucial role in heating the solar corona, but the exact energy release mechanism(s) is(are) still unknown. Here, we investigate in detail, the process of magnetic energy release in a situation where two initially independent flux systems are forced into each other. Work done by the foot point motions goes in to building a current sheet in which magnetic reconnection takes place. The scaling relations of the energy input and output are determined as functions of the driving velocity and the strength of fluxes in the independent flux systems. In particular, it is found that the energy injected into the system is proportional to the distance travelled not the rate of travel. Similarly, the rate of Joule dissipation is related to the distance travelled. Hence, rapidly driven foot points lead to bright, intense, but short-lived events, whilst slowly driven foot points produce weaker, but longer-lived brightenings. Integrated over the lifetime of the events both would produce the same heating if all other factors were the same. A strong overlying field has the affect of creating compact flux lobes from the sources. These appear to lead to a more rapid injection of energy, as well as a more rapid release of energy. Thus, the stronger the overlying field the more compact and more intense the heating. This means observers must know the rate of movement of the magnetic fragments involved in an events, as well as determine the strength and orientation of the surrounding field to be able to predict anything about the energy dissipated.

In this paper, we investigate a simple model where two, initially unconnected, flux systems are forced to interact in response to the imposed boundary driving by solving the non-ideal 3D MHD equations numerically. The reconnection rate of the dynamical process is determined and compared with the corresponding rate for the potential evolution of the magnetic field. This shows that the dynamic reconnection rate is about a factor of two smaller than the potential (perfect, instantaneous) rate for realistic solar driving velocities demonstrating that this three-dimensional magnetic reconnection process is fast. The energy input for a fixed advection distance is found to be independent of the driving velocity. The Joule dissipation associated with the reconnection process is also found to be basically dependent on the advection distance rather than driving velocity. This implies that the timescale for the event determines the effect the heating has on the temperature increase. Finally, the numerical experiments indicate that the observational structure of the reconnection site changes dramatically depending on the phase of the evolution of the passage of the two flux sources. In the initial phase, where the sources become connected, the heating is confined to a compact region. For the disconnecting phase the energy gets distributed over a larger area due to the reconnected field line connectivity.