The Importance of Photospheric Intense Flux Tubes for Coronal Heating

Dissipation of magnetic energy in the corona requires the creation of very fine scale-lengths because of the high magnetic Reynolds number of the plasma. The formation of current sheets is a natural possible solution to this problem and it is now known that a magnetic field that is stressed by continous photospheric motions through a series of equilibria can easily form such sheets. Furthermore, in a large class of 3D magnetic fields without null points there are locations, called quasi-separatrix layers (QSLs), where the field-line linkage changes drastically. They are the relevant generalisation of normal separatrices to configurations without nulls: along them concentrated electric currents are formed by smooth boundary motions and 3D magnetic reconnection takes place when the layers are thin enough. With a homogenous normal magnetic field component at the boundaries, the existence of thin enough QSL to dissipate magnetic energy rapidly requires that the field is formed by flux tubes that are twisted by a few turns. However, the photospheric field is not homogeneous but is fragmented into a large number of thin flux tubes. We show that such thin tubes imply the presence of a large number of very thin QSLs in the corona. The main parameter on which their presence depends is the ratio between the magnetic flux located outside the flux tubes to the flux inside. The thickness of the QSLs is approximately given by the distance between neighbouring flux tubes multiplied by the ratio of fluxes to a power between two and three (depending on the density of flux tubes). Because most of the photospheric magnetic flux is confined in thin flux tubes, very thin QSLs are present in the corona with a thickness much smaller than the flux tube size. We suggest that a turbulent resistivity is triggered in a QSL, which then rapidly evolves into a dynamic current sheet that releases energy by fast reconnection at a rate that we estimate to be sufficient to heat the corona. We conclude that the fragmentation of the photospheric magnetic field stimulates the dissipation of magnetic energy in the corona.     

Magnetic-field structure associated with coronal streamers

An interesting coronal structure was observed during the solar eclipse of May 30, 1965. This comprised a series of bright arches centered approximately on a quiescent prominence. A bright ray originated near the top of one of the arches and pointed almost radially away from the photosphere. The ray could be followed for 1.5 solar radii and was deflected towards a direction parallel to the equatorial plane. By comparing the photographs with Fraunhofer maps and magnetograms, the following interpretation of the structure was obtained. The prominence lies above the neutral line of an extended bipolar magnetic region. The bright arches coincide with flux tubes arising from small photospheric regions of enhanced magnetic-field strength. The ray represents a projection view of a thin region of enhanced plasma density in the neighborhood of a current sheet which separates two flux tubes of opposite polarity. The ray is interpreted as a coronal streamer, and it is suggested that all streamers are related to current sheets.     

Statistical Flux Tube Properties of 3D Magnetic Carpet Fields

The quiet-Sun photosphere consists of numerous magnetic flux fragments of both polarities that evolve with granular and supergranular flow fields. These concentrations give rise to a web of intermingled magnetic flux tubes which characterise the coronal magnetic field. Here, the nature of these flux tubes is studied. The photosphere is taken to be the source plane and each photospheric fragment is represented by a series of point sources. By analysing the potential field produced by these sources, it is found that the distribution of flux tube lengths obtained by (i) integrating forward from positive sources and (ii) tracing back from negative sources is highly dependent on the total flux imbalance within the region of interest. It is established that the relation between the footpoint separation of a flux tube and its height cannot be assumed to be linear. Where there is a significant imbalance of flux within a region, it is found that fragments of the dominant polarity will have noticeably more connections, on average, than the minority polarity fragments. Despite this difference, the flux from a single fragment of either polarity is typically divided such that (i) 60–70% connects to one opposite-polarity fragment, (ii) 25–30% goes to a further 1 to 2 opposite-polarity fragments, and (iii) any remaining flux may connect to as many as another 50 or more other opposite-polarity fragments. This is true regardless of any flux imbalance within the region. It is found that fragments connect preferentially to their nearest neighbours, with, on average, around 60–70% of flux closing down within 10 Mm of a typical fragment. Only 50% of the flux in a quiet region extends higher than 2.5 Mm above the solar surface and 5–10% extends higher than 25 Mm. The fragments that contribute to the field above this height cover a range of sizes, with even the smallest of fragments contributing to the field at heights of over 50 Mm.     

Topology and current ribbons: A model for current,reconnection and flaring in a complex,evolving corona

Magnetic field enters the corona from the interior of the Sun through isolated magnetic features on the solar surface. These features correspond to the tops of submerged magnetic flux tubes, and coronal field lines often connect one flux tube to another, defining a pattern of inter-linkage. Using a model field, in which flux tubes are represented as point magnetic charges, it is possible to quantify this inter-linkage. If the coronal field were current-free then motions of the magnetic features would change the inter-linkage through implicit (vacuum) magnetic reconnection. Without reconnection the conductive corona develops currents to avoid changing the flux linkage. This current forms singular layers (ribbons) flowing along topologically significant field lines called separators. Current ribbons store magnetic energy as internal stress in the field: the amount of energy stored is a function of the flux tube displacement. To explore this process we develop a model called the minimum-current corona (MCC) which approximates the current arising on a separator in response to displacement of photospheric flux. This permits a model of the quasi-static evolution of the corona above a complex active region. We also introduce flaring to rapidly change the flux inter-linkage between magnetic features when the internal stress on a separator becomes too large. This eliminates the separator current and releases the energy stored by it. Implementation of the MCC in two examples reveals repeated flaring during the evolution of simple active regions, releasing anywhere from 10²⁷–10²⁹ ergs, at intervals of hours. Combining the energy and frequency gives a general expression for heat deposition due to flaring (i.e., reconnection).     

On the appearance of magnetic flux in the solar photosphere

Ideas and models for the appearance of photospheric magnetic structure are confronted with observational data. Some findings are: The magnetic flux emerging in an active region consists of a bundle of flux tubes which were already concentrated before penetrating into the photosphere. A model of a rising bunch of flux tubes joining into a few strands at larger depths describes the coalescence of spots near the leading and following edges of the active region while more flux may surface near the center of the region. There is no observational evidence for appreciable helical twists in the flux bundles.Throughout the region's lifetime the magnetic elements move coherently, the whole magnetic structure rotates faster than the quiet photosphere. In active regions the convective flow at scales larger than the granulation is arrested by the magnetic structure. The long-lived supergranular cells around spots and in the enhanced network in turn determine the decay properties of spots and facular clusters. The modulation of the convective flow by the magnetic structure explains the slow dispersal of faculae.The hierarchy of magnetic elements (sunspots-pores-knots-facular clusters-facular points) may be explained by a set of magnetostatic flux tube models in the top of the convection zone. The underlying assumptions are that the heat flow along the magnetic field is reduced and that there is no heat exchange across the field except by radiation.A tentative model is proposed to account for the amplification, ascent and emergence of intense flux bundles. The assumptions are: (i) the field is concentrated in toroidal bundles by differential rotation, (ii) in the deep convection zone flux bundles are contained by the external turbulent pressure, and (iii) for field strengths up to the equipartition value efficient lateral heat exchange is possible. After a loop has surfaced radiative cooling and subsequent convective downflow reduce the temperature in the top of the flux tubes which then contract to field strengths well above the local equipartition value. There the heat flow is channelled along the field, which creates the conditions for the magnetostatic flux tube models without requiring a blocking of the heat flow somewhere within the tubes.The paper contains a brief review on the evolution of the magnetic field from the emergence in active regions up to the enigmatic disappearance, and a list of topics for further observational investigation.     

Time-dependent interaction of granules with magnetic flux tubes

The time-dependent interaction of the granulation velocity field with a magnetic flux tube is investigated here. It is seen that when a magnetic field line is displaced normal to itself so as to simulate thebuffeting action of granules, a flow of gas is initiated along the field. By choosing a lateral velocity field which is consistent with observations of granules, it is found that the resulting gas motion is a downward flow with a velocity compatible with the observed downflow in isolated photospheric flux tubes. It is therefore proposed that the observed photospheric downflow is a manifestation of the interaction of granules with flux tubes.     

Topological Departures from Translational Invariance along a Filament Observed by THEMIS

We study the topology of the 3D magnetic field in a filament channel to address the following questions: Is a filament always formed in a single flux tube? How does the photospheric magnetic field lead to filament interruptions and to feet formation? What is the relation between feet-related field lines and the parasitic polarities? What can topological analyses teach us about EUV filament channels? To do so, we consider a filament observed on 6 October 2004 with THEMIS/MTR, in Hα with the full line profile simultaneously and cospatially with its photospheric vector magnetic field. The coronal magnetic field was calculated from a “linear magnetohydrostatic” extrapolation of a composite THEMIS-MDI magnetogram. Its free parameters were adjusted to get the best match possible between the distribution of modeled plasma-supporting dips and the Hα filament morphology. The model results in moderate plasma β≤1 at low altitudes in the filament, in conjunction with non-negligible departures from force-freeness measured by various metrics. The filament here is formed by a split flux tube. One part of the flux tube is rooted in the photosphere aside an observed interruption in the filament. This splitted topology is due to strong network polarities on the edge of the filament channel, not to flux concentrations closer to the filament. We focus our study to the northwest portion of the filament. The related flux tube is highly fragmented at low altitudes. This fragmentation is due to small flux concentrations of two types. First, some locally distort the tube, leading to noticeable thickness variations along the filament body. Second, parasitic polarities, associated with filament feet, result in secondary dips above the related local inversion line. These dips belong to long field lines that pass below the flux tube. Many of these field lines are not rooted near the related foot. Finally, the present model shows that the coronal void interpretation cannot be ruled out to interpret the wideness of EUV filament channels.     

Global Simulations of Differentially Rotating Magnetized Disks: Formation of Low-beta Filaments and Structured Coronae

We examine the excitation of oscillations in the magnetic network of the Sun through the footpoint motion of photospheric magnetic flux tubes located in intergranular lanes. The motion is derived from a time series of high-resolution G-band and continuum filtergrams using an object-tracking technique. We model the response of the flux tube to the footpoint motion in terms of the Klein-Gordon equation, which is solved analytically as an initial value problem for transverse (kink) waves. We compute the wave energy flux in upward-propagating transverse waves. In general we find that the injection of energy into the chromosphere occurs in short-duration pulses, which would lead to a time variability in chromospheric emission that is incompatible with observations. Therefore, we consider the effects of turbulent convective flows on flux tubes in intergranular lanes. The turbulent flows are simulated by adding high-frequency motions (periods 5-50 s) with an amplitude of 1 km s(-1). The latter are simulated by adding random velocity fluctuations to the observationally determined velocities. In this case, we find that the energy flux is much less intermittent and can in principle carry adequate energy for chromospheric heating.     

Fip Enhancement by Alfvén Ionization

Alfvén ionization is offered as a possible mechanism underlying the enhanced population of low first ionization potential (FIP) species in the solar corona. In this process, the photospheric flow impinging on the magnetic structure of a coronal flux tube collides with, and displaces, ions in the magnetised plasma within the flux tube. This leads to pockets of charge imbalance that persist due to the impeded electron transport perpendicular to the magnetic field. The localised electric field then energises electrons to the impact ionization energy threshold of low-FIP components in the surface flow. Such species remain trapped in the plasma, and drift up the magnetic structure, causing a localised population enhancement compared to photospheric levels. We find that this mechanism successfully accounts for observed biases for flow speeds known to exist in the photosphere, and moreover explains certain anomalous abundances which do not fit into existing theories.     

On Signatures of Twisted Magnetic Flux Tube Emergence

Recent studies of NOAA active region 10953, by Okamoto et al. (Astrophys. J. Lett. 673, 215, 2008; Astrophys. J. 697, 913, 2009), have interpreted photospheric observations of changing widths of the polarities and reversal of the horizontal magnetic field component as signatures of the emergence of a twisted flux tube within the active region and along its internal polarity inversion line (PIL). A filament is observed along the PIL and the active region is assumed to have an arcade structure. To investigate this scenario, MacTaggart and Hood (Astrophys. J. Lett. 716, 219, 2010) constructed a dynamic flux emergence model of a twisted cylinder emerging into an overlying arcade. The photospheric signatures observed by Okamoto et al. (2008, 2009) are present in the model although their underlying physical mechanisms differ. The model also produces two additional signatures that can be verified by the observations. The first is an increase in the unsigned magnetic flux in the photosphere at either side of the PIL. The second is the behaviour of characteristic photospheric flow profiles associated with twisted flux tube emergence. We look for these two signatures in AR 10953 and find negative results for the emergence of a twisted flux tube along the PIL. Instead, we interpret the photospheric behaviour along the PIL to be indicative of photospheric magnetic cancellation driven by flows from the dominant sunspot. Although we argue against flux emergence within this particular region, the work demonstrates the important relationship between theory and observations for the successful discovery and interpretation of signatures of flux emergence.     

Interaction and Dynamics of the Photospheric Network Magnetic Elements

Small-scale magnetic elements in the quiet photospheric network are believed to play a key role in the energy flow from the solar surface to upper layers of atmosphere. Their intense hydro-magnetic activity includes merging and fragmentation of same polarity fluxes, `total' or partial cancellation of neighboring flux elements of opposite polarity, dynamic appearance and disappearance of compact bipoles, etc. We study the general features of these processes, and show that non-collinearity of flux tubes, sharp stratification of low atmosphere and finite plasma beta lead to several specific effects in the interacting flux tubes that may explain the morphological properties of network magnetic field and also provide a mechanism for the energy build up and release in the nearby chromosphere and transition region. We show that during the collision of flux tubes in the photosphere reconnection occurs regardless of whether the flux tubes are of opposite or of the same polarity. But the dynamics of reconnection products are significantly different and lead to different macroscopic effects that can be observed.     

RENDERING THREE-DIMENSIONAL SOLAR CORONAL STRUCTURES

An X-ray or EUV image of the corona or chromosphere is a 2D representation of an extended 3D complex for which a general inversion process is impossible. A specific model must be incorporated in order to understand the full 3D structure. We approach this problem by modeling a set of optically-thin 3D plasma flux tubes which we render these as synthetic images. The resulting images allow the interpretation of the X-ray/EUV observations to obtain information on (1) the 3D structure of X-ray images, i.e., the geometric structure of the flux tubes, and on (2) the internal structure using specific plasma characteristics, i.e., the physical structure of the flux tubes. The data-analysis technique uses magnetograms to characterize photospheric magnetic fields and extrapolation techniques to form the field lines. Using a new set of software tools, we have generated 3D flux tube structures around these field lines and integrated the plasma emission along the line of sight to obtain a rendered image. A set of individual flux-tube images is selected by a non-negative least-squares technique to provide a match with an observed X-ray image. The scheme minimizes the squares of the differences between the synthesized image and the observed image with a non-negative constraint on the coefficients of the brightness of the individual flux-tube loops. The derived images are used to determine the specific photospheric foot points and physical data, i.e., scaling laws for densities and loop lengths. The development has led to computer efficient integration and display software that is compatible for comparison with observations (e.g., Yohkoh SXT data, NIXT, or EIT). This analysis is important in determining directly the magnetic field configuration, which provides the structure of coronal loops, and indirectly the electric currents or waves, which provide the energy for the heating of the plasma. We have used very simple assumptions (i.e., potential magnetic fields and isothermal corona) to provide an initial test of the techniques before complex models are introduced. We have separated the physical and geometric contributions of the emission for a set of flux tubes and concentrated, in this initial study, on the geometric contributions by making approximations to the physical contributions. The initial results are consistent with the scaling laws derived from the Yohkoh SXT data.     

Formation of a Flare-Productive Active Region: Observation and Numerical Simulation of NOAA AR 11158

We present a comparison of the Solar Dynamics Observatory (SDO) analysis of NOAA Active Region (AR) 11158 and numerical simulations of flux-tube emergence, aiming to investigate the formation process of this flare-productive AR. First, we use SDO/Helioseismic and Magnetic Imager (HMI) magnetograms to investigate the photospheric evolution and Atmospheric Imaging Assembly (AIA) data to analyze the relevant coronal structures. Key features of this quadrupolar region are a long sheared polarity inversion line (PIL) in the central δ-sunspots and a coronal arcade above the PIL. We find that these features are responsible for the production of intense flares, including an X2.2-class event. Based on the observations, we then propose two possible models for the creation of AR 11158 and conduct flux-emergence simulations of the two cases to reproduce this AR. Case 1 is the emergence of a single flux tube, which is split into two in the convection zone and emerges at two locations, while Case 2 is the emergence of two isolated but neighboring tubes. We find that, in Case 1, a sheared PIL and a coronal arcade are created in the middle of the region, which agrees with the AR 11158 observation. However, Case 2 never builds a clear PIL, which deviates from the observation. Therefore, we conclude that the flare-productive AR 11158 is, between the two cases, more likely to be created from a single split emerging flux than from two independent flux bundles.     

On the origin of strong downdrafts associated with the birth of sunspots

During the emerging phase of sunspot pores strong downdrafts are observed in photospheric lines, whereas the upward flow is not detected within the observational accuracy. For the study of the origin of these downdrafts, we present the numerical solutions of non-steady hydrodynamic equations in one dimension along the rising magnetic flux tubes. In these solutions, it is assumed that the solar gas inside the tube is initially in the hydrostatic HSRA and then the tube rises with a given velocity and configuration. The results reproduce well the observed radial velocity distributions, hence it is concluded that the observed strong downdrafts originate from the sliding matter along the rising magnetic flux tubes.     

Linear MHD Sausage Waves in Compressible Magnetically Twisted Flux Tubes

Oscillations of magnetic flux tubes are of great importance as they contain information about the geometry and fine structure of the flux tubes. Here we derive and analytically solve in terms of Kummer’s functions the linear governing equations of wave propagation for sausage surface and body modes (m=0) of a magnetically twisted compressible flux tube embedded in a compressible uniformly magnetized plasma environment in cylindrical geometry. A general dispersion relation is obtained for such flux tubes. Numerical solutions for the phase velocity are obtained for a wide range of wavenumbers and for varying magnetic twist. The effect of magnetic twist on the period of oscillations of sausage surface modes for different values of the wavenumber and vertical magnetic field strength is calculated for representative photospheric and coronal conditions. These results generalize and extend previous studies of MHD waves obtained for incompressible or for compressible but nontwisted flux tubes. It is found that magnetic twist may change the period of sausage surface waves of the order of a few percent when compared to counterparts in straight nontwisted flux tubes. This information will be most relevant when high-resolution observations are used for diagnostic exploration of MHD wave guides in analogy to solar-interior studies by means of global eigenoscillations in helioseismology.     

Non-equilibrium of Magnetic Flux Tubes emerging into the Solar Corona

A lack of equilibrium of twisted magnetic flux tubes emerging from the photosphere into the corona is considered. Assuming mass and flux conservation in the tube and an isothermal tube that is in thermal equilibrium with the surrounding plasma, it is shown that a sufficently rapid temperature increase through the transition zone may lead to the loss of magnetohydrostatic equilibrium of the emerging flux tube due to the enhancement of the plasma pressure inside the tube. The non-equilibrium leads to a rapid expansion of the tube to reach a new equilibrium state. The rise and expansion of the tube before and after the non-equilibrium are accompanied by an increase in the twist of the magnetic field. This may lead to the field exceeding the threshold for the onset of the kink instability and a subsequent explosive release of magnetic energy.     

Thin flux tube models for star spots

The possibility of understanding stellar activity as an up‐scaled version of the activity of our Sun is investigated. A theoretical model to explain properties of sunspots is used for explaining observed latitudes of star spots. The model is based on thin‐flux‐tube simulations that study the path of magnetic flux tubes from their origin in a stellar overshoot layer to photospheric layers. A direct comparison of the simulation results with individual stars is given. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The magnetic structure of arch filament systems

Photographic-type magnetograms are used in conjunction with H filtergrams to study the structure and evolution of magnetic fields associated with arch filament systems. The magnetograms show that the opposite ends of the arch filaments are indeed rooted in photospheric magnetic fields of opposite polarity. Furthermore, these magnetic field systems are in every case new magnetic flux appearing at the solar surface. Time lapse studies show the detailed process by which the flux tubes emerge through the surface. First, supergranules bring individual strands of magnetic flux to the surface and sweep the two feet of the flux tube to opposite sides of the supergranule. Then, the flux tube rises through the chromosphere, creating a visible arch filament. It is also shown that the observed rotation of the axis of an arch filament system in the plane of the solar surface is caused by the emergence of successive flux loops, each possessing different axial tilts.     

Alfvén waves in the solar atmosphere

The nonlinear propagation of Alfvén waves on open solar magnetic flux tubes is considered. The flux tubes are taken to be vertical and axisymmetric, and they are initially untwisted. The Alfvén waves are time-dependent axisymmetric twists. Their propagation into the chromosphere and corona is investigated by solving numerically a set of nonlinear time-dependent equations, which couple the Alfvén waves into motions parallel to the initial magnetic field (motion in the third coordinate direction is artificially suppressed). The principal conclusions are: (1) Alfvén waves can steepen into fast shocks in the chromosphere. These shocks can pass through the transition region into the corona, and heat the corona. (2) As the fast shocks pass through the transition region, they produce large-velocity pulses in the direction transverse to B _o. The pulses typically have amplitudes of 60 km s^–1 or so and durations of a few tens of seconds. Such features may have been observed, suggesting that the corona is in fact heated by fast shocks. (3) Alfvén waves exhibit a strong tendency to drive upward flows, with many of the properties of spicules. Spicules, and the observed corrugated nature of the transition region, may therefore be by-products of magnetic heating of the corona. (4) It is qualitatively suggested that Alfvén waves may heat the upper chromosphere indirectly by exerting time-dependent forces on the plasma, rather than by directly depositing heat into the plasma.