The rise of twisted magnetic flux tubes

Sunspots are caused by the eruption of magnetic flux tubes through the solar photosphere: current theories of the internal magnetic field of the Sun suggest that such tubes must rise relatively unscathed from the base of the convection zone. In order to understand how the structure of the magnetic field within a buoyant flux tube affects its stability as it rises, we have considered the quasi-two-dimensional rise of isolated magnetic flux tubes through an adiabatically stratified atmosphere. The magnetic field is initially helical; we have investigated a range of initial field configurations, varying the distribution and strength of the twist of the field.     

Development of structure in pores and sunspots: flows around axisymmetric magnetic flux tubes

Flux elements, pores and sunspots form a family of magnetic features observed at the solar surface. As a first step towards developing a fully non-linear model of the structure of these features and of the dynamics of their interaction with solar convection, we conduct numerical experiments on idealized axisymmetric flux tubes in a compressible convecting atmosphere in cylindrical boxes of radius up to 8 times their depth. We find that the magnetic field strength of the flux tubes is roughly independent of both distance from the centre and the total flux content of the flux tube, but that the angle of inclination from the vertical of the field at the edge of the tube increases with flux content. In all our calculations, fluid motion converges on the flux tube at the surface. The results compare favourably with observations of pores; in contrast, large sunspots lie at the centre of an out-flowing moat cell. We conjecture that there is an inflow hidden beneath the penumbrae of large spots, and that this inflow is responsible for the remarkable longevity of such features.     

A New Model for Propagating Parts of EIT Waves: A Current Shell in a CME

EIT waves are observed in EUV as bright fronts. Some of these bright fronts propagate across the solar disk. EIT waves are all associated with a flare and a CME and are commonly interpreted as fast-mode magnetosonic waves. Propagating EIT waves could also be the direct signature of the gradual opening of magnetic field lines during a CME. We quantitatively addressed this alternative interpretation. Using two independent 3D MHD codes, we performed nondimensional numerical simulations of a slowly rotating magnetic bipole, which progressively result in the formation of a twisted magnetic flux tube and its fast expansion, as during a CME. We analyse the origins, the development, and the observability in EUV of the narrow electric currents sheets that appear in the simulations. Both codes give similar results, which we confront with two well-known SOHO/EIT observations of propagating EIT waves (7 April and 12 May 1997), by scaling the vertical magnetic field components of the simulated bipole to the line of sight magnetic field observed by SOHO/MDI and the sign of helicity to the orientation of the soft X-ray sigmoids observed by Yohkoh/SXT. A large-scale and narrow current shell appears around the twisted flux tube in the dynamic phase of its expansion. This current shell is formed by the return currents of the system, which separate the twisted flux tube from the surrounding fields. It intensifies as the flux tube accelerates and it is co-spatial with weak plasma compression. The current density integrated over the altitude has the shape of an ellipse, which expands and rotates when viewed from above, reproducing the generic properties of propagating EIT waves. The timing, orientation, and location of bright and faint patches observed in the two EIT waves are remarkably well reproduced. We conjecture that propagating EIT waves are the observational signature of Joule heating in electric current shells, which separate expanding flux tubes from their surrounding fields during CMEs or plasma compression inside this current shell. We also conjecture that the bright edges of halo CMEs show the plasma compression in these current shells.     

The impact of meridional circulation on stellar butterfly diagrams and polar caps

Observations of rapidly rotating solar-like stars show a significant mixture of opposite-polarity magnetic fields within their polar regions. To explain these observations, models describing the surface transport of magnetic flux demand the presence of fast meridional flows. Here, we link subsurface and surface magnetic flux transport simulations to investigate (i) the impact of meridional circulations with peak velocities of  ≤125 m s⁻¹  on the latitudinal eruption pattern of magnetic flux tubes and (ii) the influence of the resulting butterfly diagrams on polar magnetic field properties. Prior to their eruption, magnetic flux tubes with low field strengths and initial cross-sections below  ∼300 km  experience an enhanced poleward deflection through meridional flows (assumed to be polewards at the top of the convection zone and equatorwards at the bottom). In particular, flux tubes which originate between low and intermediate latitudes within the convective overshoot region are strongly affected. This latitude-dependent poleward deflection of erupting magnetic flux renders the wings of stellar butterfly diagrams distinctively convex. The subsequent evolution of the surface magnetic field shows that the increased number of newly emerging bipoles at higher latitudes promotes the intermingling of opposite polarities of polar magnetic fields. The associated magnetic flux densities are about 20 per cent higher than in the case disregarding the pre-eruptive deflection, which eases the necessity for fast meridional flows predicted by previous investigations. In order to reproduce the observed polar field properties, the rate of the meridional circulation has to be of the order of 100 m s⁻¹, and the latitudinal range from which magnetic flux tubes originate at the base of the convective zone (≲50°) must be larger than in the solar case (≲35°).     

Helicity of solar active regions from a dynamo model

We calculate helicities of solar active regions based on the idea that poloidal flux lines get wrapped around a toroidal flux tube rising through the convection zone, thereby giving rise to the helicity. We use our solar dynamo model based on the Babcock-Leighton α-effect to study how helicity varies with latitude and time.     

The role of emerging bipoles in the formation of a sunspot penumbra

The generation of magnetic flux in the solar interior and its transport from the convection zone into the photosphere, the chromosphere, and the corona will be in the focus of solar physics research for the next decades. With 4 m class telescopes, one plans to measure essential processes of radiative magneto‐hydrodynamics that are needed to understand the nature of solar magnetic fields. One key‐ingredient to understand the behavior of solar magnetic field is the process of flux emergence into the solar photosphere, and how the magnetic flux reorganizes to form the magnetic phenomena of active regions like sunspots and pores. Here, we present a spectropolarimetric and imaging data set from a region of emerging magnetic flux, in which a proto‐spot without penumbra forms a penumbra. During the formation of the penumbra the area and the magnetic flux of the spot increases. First results of our data analysis demonstrate that the additional magnetic flux, which contributes to the increasing area of the penumbra, is supplied by the region of emerging magnetic flux. We observe emerging bipoles that are aligned such that the spot polarity is closer to the spot. As an emerging bipole separates, the pole of the spot polarity migrates towards the spot, and finally merges with it. We speculate that this is a fundamental process, which makes the sunspot accumulate magnetic flux. As more and more flux is accumulated a penumbra forms and transforms the proto‐spot into a full‐fledged sunspot (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

prior to surface appearance of sunspot flux tubes and magneto-thermal pulsation of the Sun

We found an evidence that the luminosity of the Sun systematically decreased about 20 days before sunspot surface appearance by analysing time-lag correlation of time derivatives of running mean time profiles of the data of the Active Cavity Radiometer Irradiance Monitor (ACRIM) I experiment on board of Solar Maximum Mission (SMM) and of the data of the daily sunspot number. This indicates that sunspot flux tube cooling and heat transport blocking by the flux tubes start to take place in the interior of the solar convection zone well before the sunspot surface appearance. From this finding and our previous finding that the luminosity of the Sun systematically increased and the blocked heat appeared on the surface about 50 days after the sunspot surface appearance, a new view of sunspot formation and dynamics and a new view of the luminosity modulation emerged. (i) Sunspots of a solar cycle are formed from clusters of flux tubes which can be seen in the running mean time profile of the sunspot number as a peak with duration on the order of 100 to 200 days. (ii) Heat flow is blocked by the cluster of sunspot flux tubes inside the convection zone to decrease the luminosity about 20 days before the surface emergence of the sunspot cluster. (iii) The blocked heat appears on the surface about 50 days after the surface emergence of the cluster of sunspot flux tubes to heat up the surface. This appears as a thermal pulse in the running mean time profile of the ACRIM dat in between the peaks of the sunspot running mean time profile. This process of heating the surface makes the temperature gradient less steep and weakens the buoyancy of sunspot flux tubes below the surface. (vi) The radiative cooling of the surface layer by the excess heat release steepens the temperature gradient so that the buoyancy of the sub-surface magnetic flux tubes becomes stronger to cause the next surge of emergence of a cluster of sunspots and other magnetic activities, which creates a peak in the time profile of the sunspot number. We call this peak a magnetic pulse of the Sun and the coupled process of alternating pulsed appearance of heat and sunspots the magneto-thermal pulsation of the Sun.     

Spot signatures of a solar‐type dynamo in close binary stars

Magnetic activity signatures in the atmosphere of active stars can be used to place constrains on the underlying processes of flux transport and dynamo operation in its convective envelope. The ‘solar paradigm’ for magnetic activity suggests that the magnetic field is amplified and stored at the base of the convection zone. Once a critical field strength is exceeded, perturbations initiate the onset of instabilities and the growth of magnetic flux loops, which rise through the convection zone, emerge at the stellar surface, and eventually lead to the formation of starspots and active regions. In close binaries, the proximity of the companion star breaks the rotational symmetry. Although the magnitude of tidal distortions is rather small, non‐linear MHD simulations have nevertheless shown in the case of main‐sequence binary components that they can cause non‐uniform surface distributions of flux tube eruptions. The present work extends the investigation to post‐mainsequence components to explore the specific influence of the stellar structure on the surface pattern of erupting flux tubes. In contrast to the case of main‐sequence components, where the consistency between simulation results and observations supports the presumption of a solar‐like dynamo mechanism, the numerical results here do not recover the starspot properties frequently observed on evolved binary components. This aspect points out an insufficiency of the applied flux tube model and leads to the conclusion that additional flux transport and possibly amplification mechanisms have to be taken into account. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stability of magnetic flux tubes in close binary stars

Photometric and Doppler imaging observations of active binaries indicate the existence of starspots at preferred longitudes (position angles with respect to the companion star). We investigate the stability of magnetic flux tubes in the convection zone of close, fast‐rotating binary stars and explore whether the observed preferred longitudes could be caused by tidal forces and the deformation of the active star. We assume a synchronized binary system with spin axes perpendicular to the orbital plane and a rotation period of a few days. The tidal force and the deviation from spherical structure are considered in lowest‐order perturbation theory. The magnetic field is in the form of toroidal magnetic flux rings, which are stored in mechanical equilibrium within the stably stratified overshoot region beneath the convection zone until the field has grown sufficiently strong for the undulatory instability to initiate the formation of rising loops. Frequencies and geometry of stable as well as growth rates of unstable eigenmodes are determined by linear stability analysis. Particular consideration is given to the question whether the effects of tidal forces and perturbations of the stellar structure can force a rising flux loop to enter the convection zone at specific longitudes.     

How to Reach Superequipartition Field Strengths in Solar Magnetic Flux Tubes

A number of independent arguments indicate that the toroidal flux system responsible for the sunspot cycle is stored at the base of the convection zone in the form of flux tubes with field strength close to 10⁵ G. Although the evidence for such strong fields is quite compelling, how such field strength can be reached is still a topic of debate. Flux expulsion by convection should lead to about the equipartition field strength, but the magnetic energy density of a 10⁵-G field is two orders of magnitude larger than the mean kinetic energy density of convective motions. Line stretching by differential rotation (i.e., the “Ω effect” in the classical mean-field dynamo approach) probably plays an important role, but arguments based on energy considerations show that it does not seem feasible that a 10⁵-G field can be produced in this way. An alternative scenario for the intensification of the toroidal flux system in the overshoot layer is related to the explosion of rising, buoyantly unstable magnetic flux tubes, which opens a complementary mechanism for magnetic-field intensification. A parallelism is pointed out with the mechanism of “convective collapse” for the intensification of photospheric magnetic flux tubes up to field strengths well above equipartition; both mechanisms, which are fundamentally thermal processes, are reviewed.     

Meridional flow profile measurements with SOHO/MDI

We present meridional flow measurements of the Sun using a novel helioseismic approach for analyzing SOHO/MDI data in order to push the current limits in radial depth. Analyzing three consecutive months of data during solar minimum, we find that the meridional flow is as expected poleward in the upper convection zone, turns equatorward at a depth of around 40 Mm (∼ 0.95 R_⊙), and possibly changes direction again in the lower convection zone. This may indicate two meridional circulation cells in each hemisphere, one beneath the other. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Formation of polar starspots through meridional circulation

To explain the observed intermingling of polarities in the magnetic field distributions of rapidly rotating stars, surface magnetic flux transport models demand the presence of fast meridional flows.We combine simulations of the pre-eruptive and post-eruptive magnetic flux transport in cool stars to investigate the influence of a fast meridional circulation on the latitudinal eruption pattern of magnetic flux tubes and on the polar magnetic field properties. Magnetic flux tubes rising through the convection zone experience an enhanced latitude-dependent poleward deflection through meridional flows, which renders the wings of stellar butterfly diagrams convex. The larger amount of magnetic flux emerging at higher latitudes supports the intermingling of opposite polarities of polar magnetic fields and yields magnetic flux densities in the polar regions about 20% higher than in the case disregarding the pre-eruptive deflection. Taking the pre-eruptive evolution of magnetic flux into account therefore eases the need for the fast meridional flows predicted by previous investigations. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the Connection Between Mean Field Dynamo Theory and Flux Tubes

Mean field dynamo theory deals with various mean quantities and does not directly throw any light on the question of existence of flux tubes. We can, however, draw important conclusions about flux tubes in the interior of the Sun by combining additional arguments with the insights gained from solar dynamo solutions. The polar magnetic field of the Sun is of order 10 G, whereas the toroidal magnetic field at the bottom of the convection zone has been estimated to be 100000 G. Simple order-of-magnitude estimates show that the shear in the tachocline is not sufficient to stretch a 10 G mean radial field into a 100000 G mean toroidal field. We argue that the polar field of the Sun must get concentrated into intermittent flux tubes before it is advected to the tachocline. We estimate the strengths and filling factors of these flux tubes. Stretching by shear in the tachocline is then expected to produce a highly intermittent magnetic configuration at the bottom of the convection zone. The meridional flow at the bottom of the convection zone should be able to carry this intermittent magnetic field equatorward, as suggested recently by Nandy and Choudhuri (2002). When a flux tube from the bottom of the convection zone rises to a region of pre-existing poloidal field at the surface, we point out that it picks up a twist in accordance with the observations of current helicities at the solar surface.     

Dynamical processes in flux tabes and their role in chromospheric heating

We model the dynamical interaction between magnetic flux tubes and granules in the solar photosphere which leads to the excitation of transverse (kink) and longitudinal (sausage) tube waves. The investigation is motivated by the interpretation of network oscillations in terms of flux tube waves. The calculations show that for magnetic field strengths typical of the network, the energy flux in transverse waves is higher than in longitudinal waves by an order of magnitude. But for weaker fields, such as those that might be found in internetwork regions, the energy fluxes in the two modes are comparable. Using observations of footpoint motions, the energy flux in transverse waves is calculated and the implications for chromospheric heating are pointed out.     

Sausage MHD Waves in Incompressible Flux Tubes with Twisted Magnetic Fields

Twisted magnetic flux tubes are of considerable interest because of their natural occurrence from the Sun’s interior, throughout the solar atmosphere and interplanetary space up to a wide range of applicabilities to astrophysical plasmas. The aim of the present work is to obtain analytically a dispersion equation of linear wave propagation in twisted incompressible cylindrical magnetic waveguides and find appropriate solutions for surface, body and pseudobody sausage modes (i.e. m = 0) of a twisted magnetic flux tube embedded in an incompressible but also magnetically twisted plasma. Asymptotic solutions are derived in long- and short-wavelength approximations. General solutions of the dispersion equation for intermediate wavelengths are obtained numerically. We found, that in case of a constant, but non-zero azimuthal component of the equilibrium magnetic field outside the flux tube the index ν of Bessel functions in the dispersion relation is not integer any more in general. This gives rise to a rich mode-structure of degenerated magneto-acoustic waves in solar flux tubes. In a particular case of a uniform magnetic twist the total pressure is found to be constant across the boundary of the flux tube. Finally, the effect of magnetic twist on oscillation periods is estimated under solar atmospheric conditions. It was found that a magnetic twist will increase, in general, the periods of waves approximately by a few percent when compared to their untwisted counterparts.     

Comparing Simulations of Rising Flux Tubes Through the Solar Convection Zone with Observations of Solar Active Regions: Constraining the Dynamo Field Strength

We study how active-region-scale flux tubes rise buoyantly from the base of the convection zone to near the solar surface by embedding a thin flux tube model in a rotating spherical shell of solar-like turbulent convection. These toroidal flux tubes that we simulate range in magnetic field strength from 15 kG to 100 kG at initial latitudes of 1^° to 40^° in both hemispheres. This article expands upon Weber, Fan, and Miesch (Astrophys. J. 741, 11, 2011) (Article 1) with the inclusion of tubes with magnetic flux of 10²⁰ Mx and 10²¹ Mx, and more simulations of the previously investigated case of 10²² Mx, sampling more convective flows than the previous article, greatly improving statistics. Observed properties of active regions are compared to properties of the simulated emerging flux tubes, including: the tilt of active regions in accordance with Joy’s Law as in Article 1, and in addition the scatter of tilt angles about the Joy’s Law trend, the most commonly occurring tilt angle, the rotation rate of the emerging loops with respect to the surrounding plasma, and the nature of the magnetic field at the flux tube apex. We discuss how these diagnostic properties constrain the initial field strength of the active-region flux tubes at the bottom of the solar convection zone, and suggest that flux tubes of initial magnetic field strengths of ≥?40 kG are good candidates for the progenitors of large (10²¹ Mx to 10²² Mx) solar active regions, which agrees with the results from Article 1 for flux tubes of 10²² Mx. With the addition of more magnetic flux values and more simulations, we find that for all magnetic field strengths, the emerging tubes show a positive Joy’s Law trend, and that this trend does not show a statistically significant dependence on the magnetic flux.     

The effect of Kelvin-Helmholtz instability on rising flux tubes in the convection zone

If the solar dynamo operates at the bottom of the convection zone, then the magnetic flux created there has to rise to the surface. When the convection zone is regarded as passive, the rising flux is deflected by the Coriolis force to emerge at rather high latitudes, poleward of typical sunspot zones (Choudhuri and Gilman, 1987; Choudhuri, 1989). Choudhuri and D'Silva (1990) included the effects of convective turbulence on the rising flux through (a) giant cell drag and (b) momentum exchange by small-scale turbulence. The momentum exchange mechanism could enable flux tubes of radii not more than a few hundred km to emerge radially at low latitudes, but the giant cell drag mechanism required unrealistically small flux tube radii (a few meters for a reasonable giant cell upflow) to counteract the Coriolis force. We now include the additional effect of Kelvin-Helmholtz instability in a symmetrical flux ring caused by the azimuthal flow induced during its rise. The azimuthal flow crosses the threshold for the instability only if there is a giant cell upflow to drag the flux tubes appreciably. In the absence of such a drag, as in the case of a passive convection zone or in the case of momentum exchange by small-scale turbulence, the azimuthal velocity never becomes large enough to cause the instability, leaving the results of the previous calculations unaltered. The giant cell drag, aided by Kelvin-Helmholtz instability, however, becomes now a viable mechanism for curbing the Coriolis force - 10⁴ G flux tubes with radii of a few hundred km being dragged radially by upflows of 70 m s^-1.     

Convective intensification of magnetic fields in the quiet Sun

Kilogauss-strength magnetic fields are often observed in intergranular lanes at the photosphere in the quiet Sun. Such fields are stronger than the equipartition field B _e, corresponding to a magnetic energy density that matches the kinetic energy density of photospheric convection, and comparable with the field B _p that exerts a magnetic pressure equal to the ambient gas pressure. We present an idealized numerical model of three-dimensional compressible magnetoconvection at the photosphere, for a range of values of the magnetic Reynolds number. In the absence of a magnetic field, the convection is highly supercritical and characterized by a pattern of vigorous, time-dependent, 'granular' motions. When a weak magnetic field is imposed upon the convection, magnetic flux is swept into the convective downflows where it forms localized concentrations. Unless this process is significantly inhibited by magnetic diffusion, the resulting fields are often much greater than B _e and the high magnetic pressure in these flux elements leads to their being partially evacuated. Some of these flux elements contains ultraintense magnetic fields that are significantly greater than B _p. Such fields are contained by a combination of the thermal pressure of the gas and the dynamic pressure of the convective motion, and they are constantly evolving. These ultraintense fields develop owing to non-linear interactions between magnetic fields and convection; they cannot be explained in terms of 'convective collapse' within a thin flux tube that remains in overall pressure equilibrium with its surroundings.     

Mechanism of cyclically polarity reversing solar magnetic cycle as a cosmic dynamo

We briefly describe historical development of the concept of solar dynamo mechanism that generates electric current and magnetic field by plasma flows inside the solar convection zone. The dynamo is the driver of the cyclically polarity reversing solar magnetic cycle. The reversal process can easily and visually be understood in terms of magnetic field line stretching and twisting and folding in three-dimensional space by plasma flows of differential rotation and global convection under influence of Coriolis force. This process gives rise to formation of a series of huge magnetic flux tubes that propagate along iso-rotation surfaces inside the convection zone. Each of these flux tubes produces one solar cycle. We discuss general characteristics of any plasma flows that can generate magnetic field and reverse the polarity of the magnetic field in a rotating body in the Universe. We also mention a list of problems which are currently being disputed concerning the solar dynamo mechanism together with observational evidences that are to be constraints as well as verifications of any solar cycle dynamo theories of short and long term behaviors of the Sun, particularly time variations of its magnetic field, plasma flows, and luminosity.