A&A 508, 979992 (2009)
Flow instabilities of magnetic flux tubes
IV. Flux storage in the solar overshoot region^{}
E. Isik^{1,3}  V. Holzwarth^{2,3}
1 
Department of Mathematics and Computer
Science, Istanbul Kültür University, Ataköy Campus,
Bakirköy 34156, Istanbul, Turkey
2 
FraunhoferInstitut für Kurzzeitdynamik, ErnstMachInstitut,
Eckerstraße 4, 79104 Freiburg im Breisgau, Germany
3 
MaxPlanckInstitut für Sonnensystemforschung,
MaxPlanckStraße 2, 37191 KatlenburgLindau, Germany
Received 3 July 2009 / Accepted 22 September 2009
Abstract
Context. Flowinduced instabilities of magnetic flux tubes
are relevant to the storage of magnetic flux in the interiors of stars
with outer convection zones. The stability of magnetic fields in
stellar interiors is of importance to the generation and transport of
solar and stellar magnetic fields.
Aims. We consider the effects of material flows on the dynamics
of toroidal magnetic flux tubes located close to the base of the solar
convection zone, initially within the overshoot region. The problem is
to find the physical conditions in which magnetic flux can be stored
for periods comparable to the dynamo amplification time, which is of
the order of a few years.
Methods. We carry out nonlinear numerical simulations to
investigate the stability and dynamics of thin flux tubes subject to
perpendicular and longitudinal flows. We compare the simulations with
the results of simplified analytical approximations.
Results. The longitudinal flow instability induced by the
aerodynamic drag force is nonlinear in the sense that the growth rate
depends on the perturbation amplitude. This result is consistent with
the predictions of linear theory. Numerical simulations without
friction show that nonlinear Parker instability can be triggered below
the linear threshold of the field strength, when the difference in
superadiabaticity along the tube is sufficiently large. A localised
downflow acting on a toroidal tube in the overshoot region leads to
instability depending on the parameters describing the flow, as well as
the magnetic field strength. We determined ranges of the flow
parameters for which a linearly Parkerstable magnetic flux tube is
stored in the middle of the overshoot region for a period comparable to
the dynamo amplification time.
Conclusions. The longitudinal flow instability driven by
frictional interaction of a flux tube with its surroundings is relevant
to determining the storage time of magnetic flux in the solar overshoot
region. The residence time for magnetic flux tubes with
Mx
in the convective overshoot layer can be comparable to the dynamo
amplification time, provided that the average speed and the duration of
an external downflow do not exceed about 50 m s ^{1} and 100 days, respectively, and that the lateral extension of the flow is smaller than about .
Key words: Sun: interior  Sun: magnetic fields  magnetohydrodynamics (MHD)
1 Introduction
Observations of large solar active regions are indicative of an organised subsurface magnetic field along the azimuthal (eastwest) direction, with opposite polarity orientations (from east to west or vice versa) in the northern and southern hemispheres. Emerging magnetic structures are in a filamentary state, in the form of magnetic flux concentrations (flux tubes) of various sizes (e.g., sunspots, pores). Flux tubes rise in the convection zone, emerge at the surface, and form bipolar magnetic regions, which follow polarity rules (Hale's law) and systematic tilt angles (Joy's law).Theoretical studies indicate that weak magnetic fields are transported by flux expulsion (Schüssler 1984) and convective pumping (Tobias et al. 2001) to the lower boundary of the convection zone, where the toroidal magnetic field is amplified by radial and latitudinal velocity shear in the solar tachocline (see e.g., Solanki et al. 2006). The stably stratified lower convection zone, in particular the convective overshoot layer, is a likely location for the generation and storage of the largescale azimuthal magnetic flux. Numerical studies simulating a layer of horizontal magnetic field in the bottom of the solar convection zone indicate that an initially uniform field underlying a fieldfree layer leads to the formation of magnetic flux tubes by magnetic RayleighTaylor instability (e.g., Fan 2001). Following its formation, a toroidal magnetic flux tube can reach a mechanical equilibrium close to the bottom of the convection zone (MorenoInsertis et al. 1992). The emergence of magnetic flux tubes driven by magnetic buoyancy and the properties of active regions (lowlatitude emergence, tilt angles, proper motions of sunspots) require azimuthal flux densities of the order of 10^{5} G in the overshoot region (Schüssler 1996; Fan et al. 1994; D'Silva & Choudhuri 1993; Caligari et al. 1995). The possibilities for the generation of these flux densities have been reviewed by Schüssler & Rempel (2002), Schüssler & FerrizMas (2003), and FerrizMas & Steiner (2007).
The average duration of the solar activity cycle is about 11 years and the total magnetic flux emerging within one cycle ranges between orders of 10^{24}10^{25} Mx. The amplification of the toroidal magnetic field to its maximum strength in about half an activity cycle raises the following question: how can the toroidal flux be stored stably for at least a few years in the dynamo amplification region? A related problem concerns the feedback of magnetic flux loss on the rate of toroidal flux generation. To more clearly understand magnetic flux generation and storage in the convective overshoot region, it is important to determine and constrain the effects of flows on the stability and dynamics of magnetic flux tubes. Here, we focus on: (a) the nonlinear development of flowinduced flux tube instabilities, which can be dynamically significant in the course of toroidal field amplification; and (b) effects of perpendicular flows on the storage of toroidal flux tubes. Consequently, some of the flow properties prevailing in the overshoot region can also be constrained, by requiring that flux tubes with field strengths of up to a few times 10^{4} G are stored in the overshoot region for about a few years.
We consider the possibility that the frictioninduced instability (Holzwarth 2008; Holzwarth et al. 2007) leads to flux loss from the overshoot layer for field strengths below the Parker instability limit. We refer to such Parkerstable flux tubes as ``subcritical'' throughout the paper. We investigate effects of flows on magnetic flux tubes in mechanical equilibrium to obtain quantitative estimates, and to answer the following question: how strongly do (i) finiteamplitude perturbations of flux tubes by perpendicular flows (see also Schüssler & Ferriz Mas 2007, Paper I); and (ii) the frictional interaction of the tube with its surroundings limit the residence time of flux tubes in the overshoot region?
We approach the problem by considering the nonlinearity of frictioninduced instability (Sect. 2), including a discussion of the effects of finite perturbations on the magnetic buoyancy instability (Sect. 2.3), and determine the ranges of flow parameters that allow the storage of subcritical flux tubes subject to radial flows (Sect. 3).
The analyses and simulations were carried out in the framework of the thin flux tube approximation (Spruit 1981), which at present is the only existing approach that can deal with the small magnetically induced variations in density, pressure, and temperature corresponding to the high plasma (>10^{5}) of the deep solar convection zone.
2 Frictioninduced instabilities
The mechanical equilibrium of a toroidal flux tube in the solar convection zone requires that the plasma within the tube rotates faster than the surrounding medium. MorenoInsertis et al. (1992) demonstrated that a toroidal flux tube rotating initially at the same rate as the external medium can reach this equilibrium by developing an internal prograde flow. The speed of this ``equilibrium'' flow depends on the magnetic field strength, depth, and latitude. A flux tube subject to nonaxisymmetric perturbations (varying in azimuth) will develop MHD waves that propagate in prograde and retrograde directions. If the internal flow speed (in the restframe of the external medium) is greater than the phase speed of the slowest transversal retrograde wave mode, this wave mode will be advected in the prograde direction, while being amplified (Holzwarth et al. 2007, hereafter Paper II). The critical flux density for the instability is lower than that of the Parker instability. Holzwarth (2008, hereafter Paper III) extended the analysis in Paper II to toroidal flux tubes, thus including the effect of magnetic curvature force and rotation.
The linear analyses in Papers II and III considered Stokestype friction, which is linear in perpendicular velocity. To test the results of the linear stability analysis in the nonlinear context, we shall consider the drag force, which is quadratic in velocity. This test can be implemented by carrying out fully nonlinear numerical simulations of flux tubes subject to nonaxisymmetric perturbations.
2.1 Linear results
The hydrodynamic drag force exerted by a steady flow with velocity
perpendicular to a straight cylinder is given by
where is the drag coefficient, is the density of the external medium neighbouring the equilibrium flux tube, and is the crosssectional radius of the equilibrium tube.
The Stokestype friction force per unit volume,
where the constant is the friction coefficient, was considered in Papers II and III to include the frictiontype deceleration in the linearised equation of motion. In Paper II, it was shown that the critical field strength for frictioninduced instability hardly depends on the value of . Comparing Eq. (2) with Eq. (1), one finds that the friction coefficient is of the order
where the drag coefficient, , is set to be unity, because of the cylindrical crosssection and high Reynolds number (Batchelor 1967).
We consider a toroidal flux tube located in the convective overshoot region, parallel to the equatorial plane. For the stratification of the ambient medium, we use a model convection zone developed by Skaley & Stix (1991), which uses a nonlocal treatment of convection as described by Shaviv & Salpeter (1973). In the model, the overshoot region extends about 10^{4} km below the base of the convection zone, which is defined as the depth at which the convective energy flux changes its sign, at about r=512 Mm. The thickness of the overshoot layer corresponds to about 20% of the local pressure scale height. Throughout the paper, the terms bottom, middle, and top (levels) of the overshoot region are used to refer to the radial positions at, respectively, 2000, 5000, and 8000 km above the lower boundary of the overshoot region, which is at a radius of r=502 Mm ( ).
The linear stability analysis of Paper III has provided growth rates of the frictioninduced instability for toroidal flux tubes. In Fig. 1, growth rates are shown as a function of the field strength (B_{0}) and latitude () located in the middle of the overshoot region, for s^{1} and rigid solar rotation^{}. In Paper III, the value of was chosen by determining the average perpendicular velocity of the mass elements of the flux tube during the initial stages of nonlinear numerical simulations for a flux tube with Mx. The frictioninduced instability sets in for G. The onset of Parker instability is in the interval G at low latitudes.
Figure 1: Stability diagram for toroidal flux tubes located in the middle of the solar overshoot region. The shaded areas indicate unstable configurations with the grey scale representing the growth rate. Dotted lines outline regions of Parker instability. 

Open with DEXTER 
2.2 Nonlinear simulations
We carried out numerical simulations of a toroidal flux tube
using a semiimplicit finitedifference scheme developed
by MorenoInsertis (1986), that was
extended to three dimensions and spherical geometry by Caligari et al. (1995).
The numerical procedure is based on
the equations of ideal magnetohydrodynamics in the framework of
the thin flux tube approximation (Spruit 1981), in the form
given by FerrizMas & Schüssler (1995,1993).
In the numerical scheme,
the flux tube is described by a string of Lagrangian mass elements,
which move in three dimensions under the effects of various body
forces.
The equation of motion for the material inside the toroidal flux tube,
in a reference frame rotating with the angular velocity of the tube,
,
is written as
=  
(4) 
where D/D is the Lagrangian derivative, the subscript i denotes quantities inside the flux tube (the subscript e denotes external quantities in the following). The terms on the right hand side of the equation are, respectively, the total pressure force (gas and magnetic), magnetic tension force, effective gravity (including the centrifugal force), Coriolis force, and the hydrodynamic drag force, which is given by Eq. (1). We assume rigid rotation in the solar interior. Although rotational velocity shear can be incorporated into the framework of our model (FerrizMas & Schüssler 1995; Caligari et al. 1995), we disregard it for the sake of clarity, to focus on the fundamental effects of radial flows on flux tube dynamics. Caligari et al. (1995) found that external shear does not have a significant effect on the stability properties and dynamics of flux tubes in the Sun.
Figure 2: Comparison of the analytical solutions for the growth rate of the frictioninduced instability (lines) with the numerical simulations (symbols) for flux tubes at the middle of the overshoot region and latitude. The initial perturbation amplitude (), and the Stokes drag coefficient () are indicated in the legend. 

Open with DEXTER 
The initial value for the crosssectional radius of each flux tube is taken to be km for all simulations. For B_{0}=10^{5} G, this corresponds to a magnetic flux of about Mx, which is typical of a bipolar magnetic region of moderate size on the solar surface.
The flux tube in mechanical equilibrium is initially perturbed by
small displacements in three dimensions.
The azimuthal dependencies of the initial perturbations in
each coordinate k were taken to be in the form
where the amplitudes of all modes are equal. In cases of instability, the growth rate of the perturbation was determined by an exponential fit to the radial location of the top of the growing loop. Each fit was limited to the time interval corresponding to the initial exponential growth phase of the instability.
We compare the linear results for a range of friction parameters , with the results of nonlinear simulations in Fig. 2, which shows the growth rates of the frictioninduced instability as a function of the magnetic field strength, for latitude in the middle of the overshoot region, for chosen values of . The two plateaus where the curves for the linear solutions converge correspond to Parkerunstable regions according to the linear analysis (regions enclosed by dots in Fig. 1). Thus the curve for in Fig. 2 can be seen as a horizontal cut through Fig. 1 at . Each set of numerical simulations were performed with a fixed value of the total initial perturbation amplitude in the radial direction, , ranging between 51 km and 2539 km. The values of were chosen such that the curves correspond to the range of growth rates found in the simulations. The simulations exhibit an overall similarity with the linear results in terms of the field strength dependence of the growth rate. The linear growth rate of the instability based on a fixed value of corresponds approximately to a certain perturbation amplitude, , as a function of B_{0} in the linearly Parkerstable and frictionally unstable intermediate regime.
The coefficient is a measure of the strength of frictional coupling of the oscillating flux tube with the surrounding medium. For , the growth rate is proportional to , because the frictional coupling facilitates the amplification of perturbations. For , on the other hand, the growth rates begin to decrease with increasing , because too large friction impedes perpendicular movements of the flux tube and thus the development of overstability. The values chosen in Fig. 2 correspond to the regime in which friction has a destabilising effect. The reason for the correspondance between and is the quadratic dependence of the drag force on the perpendicular velocity: a larger initial perturbation leads to a larger perpendicular velocity and thus to a larger , meaning a stronger Stokesfrictional coupling. The nonlinear growth of the instability is shown in a supplemented animation (see the online appendix, Fig. C.1), for a flux tube with G, , and km.
The variation in growth time as a function of perturbation amplitude for G (see Fig. 2) is shown in Fig. 3. It shows that the frictioninduced instability grows faster for larger perturbations, following a power law.
Figure 3: Growth time as a function of the displacement amplitude, G, , at the middle of the overshoot region. The upper horizontal axis shows the flow velocity amplitude that leads to the displacement according to the analytical approximation discussed in Sect. 3.1, for m=5 and m s^{1}. 

Open with DEXTER 
There are two reasons for the nonlinear behaviour of the instability, i.e., that the growth rate is proportional to the initial perturbation amplitude: (i) is proportional to , i.e, the drag force is nonlinear; (ii) the spatial variation in the superadiabaticity, , of the external medium along the flux tube becomes increasingly effective with increasing , such that the excessive magnetic buoyancy at the highest location of the tube provides an additional upward acceleration. However, the latter ``effect'' is not likely to be the dominant source of nonlinearity, because otherwise there would be far poorer agreement between the linear and nonlinear results in Fig. 2. This point will be investigated further in Sect. 2.3.
We carried out additional simulations for linearly Parkerstable configurations with a perturbation of m=1, by varying the amplitude. The results are given in Table 1. In all cases but one for G and , friction was taken into account. For G, the flux tube is linearly stable (see Fig. 1), and is also stable to the finite perturbations applied in the simulations. For stronger fields the frictioninduced instability sets in, the growth rate then being proportional to the perturbation amplitude. For a given latitude and perturbation amplitude, a higher B_{0} leads to a more rapid growth, as predicted by the linear approach. Owing to nonlinear effects, the initial m=1 perturbation initiates perturbations of higherorder modes, in particular if the initial amplitude is large. In the case of unstable flux tubes, the eigenmode with the shortest growth time dominates the evolution and causes an increase in the top position of the flux tube. During its early evolution, this rise is approximately exponential and can be characterised by the efolding times given in Table 1. In those cases marked with asterisks, the top position of the perturbed flux tube decreases with time. Since each eigenmode decays on its own individual damping time and modes may interact nonlinearly, a unique efolding timescale of the overall decay of the top position is not possible. After all, in these cases we find that the top position decreases to about one tenth of its original value within 10^{3}  10^{4} days. The configuration G and is very close to the Parker instability boundary (see Fig. 1), and the result for a sufficiently large perturbation is a rapid growth of an instability induced by the steep gradient of the external medium along the tube, which we consider in the next section.
Table 1: Growth times for various finite perturbations.
2.3 Effect of stratification
Figure 4: Radial profile of the superadiabaticity in the outer 40 per cent of the Sun, according to the solar model used (Skaley & Stix 1991). The arrows and the corresponding circles mark the boundaries of the overshoot region. 

Open with DEXTER 
The most important external quantity in determining the stability of a toroidal flux tube in the overshoot region is the superadiabaticity, . Figure 4 shows the radial profile of in the outer 40% of the solar interior, according to the model under consideration (Skaley & Stix 1991). From the upper radiative zone to the lower convection zone, the superadiabaticity increases from negative values and changes its sign at about . It increases by a factor of about through the overshoot layer. We assume that a toroidal flux tube is deformed in such a way that it extends between the boundaries of the overshoot region. The parts of the tube that extend to higher layers (smaller ) will experience larger buoyancy force and thus be destabilised, whereas the parts that extend to deeper layers will be stabilised owing to sufficiently negative superadiabaticity. These nonlinear effects may lead to the formation of a rising flux loop.
Figure 5 shows the initial superadiabaticity difference, , between the surroundings of the top and bottom parts of the flux tube as a function of the initial perturbation amplitude, , for the case presented in Fig. 3 (Sect. 2.2). We note that the relative difference is lower than unity ( ).
To check whether a variation in is responsible for the nonlinear dependence of the growth rate on the perturbation amplitude, we carried out simulations with the same initial conditions as given above (Fig. 5), but without the drag force, to avoid a mixing with the frictional instability, thus isolating the effect. No instability was found within the considered range of perturbation amplitudes. This indicates that the effect is not responsible for the nonlinearity of the frictional instability for .
Figure 5: The difference in superadiabaticity between the surroundings of the top and bottom parts of perturbed flux tubes, which are initially located at the middle overshoot zone with and G, the same initial conditions as in Fig. 3. The drag force has been taken into account. The numbers indicate growth times of the instability for each simulation. 

Open with DEXTER 
In order to constrain the difference of superadiabaticity required for the effect to have a significant role in destabilising a flux tube, we carried out numerical simulations at the bottom of the overshoot layer, where the radial gradient of is steeper than in the upper overshoot region. Two sets of simulations were carried out for latitudes and and G, without the drag force. The two cases correspond to linearly Parkerstable configurations, very close to the instability boundary, which is about 2 and 4 kG larger for the high and lowlatitude cases, respectively. For both cases, we find instability for km. Figure 6 shows as a function of for the highlatitude case. A rising loop is formed for km, which corresponds to . The growth rates are about 60 days for km and 20 days for km, which are comparable with the rise time in the convection zone proper. The rising loop is triggered by nonlinear effects owing to a sufficiently large magnetic buoyancy difference along the flux tube, rather than by overstability of interacting wave modes.
Figure 6: Same as Fig. 5, for G, , at the bottom of the overshoot region, without the drag force. When the initial perturbation amplitude exceeds 10^{3} km, the flux tube becomes unstable owing to buoyancy differences along the tube. 

Open with DEXTER 
To summarise the results obtained in this section, we find that a superadiabaticity difference of the order of 10^{5} is required for a subcritical toroidal flux tube to become unstable in the overshoot layer, due to the effect.
2.4 The slingshot case
One may expect that rapidly rising flux loops can be formed by localised upflows at the top of the overshoot region. An interesting question is whether such an eruption of a small part of the flux tube can originate close to the top of the overshoot region. To test this possibility, we carried out a numerical simulation for a flux tube at the top of the overshoot region with G and , which is a linearly Parkerunstable configuration. The perturbation was applied with km and m=15, which corresponds to an azimuthal extension of for each lifted portion. In the course of the simulation, the wave energy in high azimuthal modes is gradually transferred into lowerm modes. The tube enters into the convection zone with m=2 mode, forming two largescale loops. This is consistent with the prediction of the linear stability analysis.
Smallscale loops can also originate in the bottom of the overshoot region, where the effect (Sect. 2.3) is significant. A slingshot effect can occur at the radiative zone boundary, if a sufficiently strong localised downflow pushes a small part of the tube downward. Subsequently, the submerged part would be ejected upwards owing to strong buoyancy. We defined initial conditions describing localised downward perturbations of Gaussian shape with the azimuthal extension ranging from about to ( ) at the bottom of the overshoot region. At the beginning of the simulation, the submerged part of the tube rises rapidly from the radiative zone boundary. However, because of the drag force, which is proportional to the square of the perpendicular velocity, it has already been rapidly decelerated within the overshoot region and its translational kinetic energy has been partly transferred into MHD waves propagating along the tube.
3 Displacement of a toroidal flux tube by radial flows
We have so far considered the dynamics of toroidal flux tubes based on the assumption of spatially perturbed initial configurations, without an explicit external driving force. In the following sections, we investigate the effects of external flows perpendicular to the tube axis, which displace the tube owing to the drag force, given by Eq. (1). The purposes are (1) to quantify the effects of external flows on the subsequent displacement; and (2) to test the possibility of storing a toroidal magnetic flux tube with Mx in the solar overshoot region for a few years, which is comparable to the dynamo amplification time. We examine the effects of spatially periodic radial flows in Sect. 3.1, and of localised downflows in Sect. 3.2.
3.1 Azimuthally periodic flow: linear analysis
We assume that a toroidal flux tube is deformed by azimuthally periodic radial flows to such an extent that the drag force is balanced by buoyancy and magnetic tension. By solving the linearised equations of motion for the perturbations of a thin flux tube, we can derive an analytical expression relating flux tube parameters to flow parameters. This relation allows us to estimate the conditions in which the deformation becomes so large that parts of the flux tube enter radiative and/or convection zones. These deformations can destabilise a flux tube in the overshoot region, e.g., by the effect.
3.1.1 Solution procedure
We adopt the linearised equations of thin magnetic flux
tubes as given by
FerrizMas & Schüssler (1995,1993) and apply the drag force exerted by an
azimuthally periodic flow of the form
The azimuthal wavenumber, m, measures the azimuthal extension of the perpendicular flow, . We assume that the flux tube reaches a stationary equilibrium state, where the drag force is balanced by buoyancy and magnetic tension forces, so that the time derivative in the momentum equation vanishes. The components of the resulting linearised equations of motion for perturbations in steady state are given in Appendix A, Eqs. (A.1)(A.3). For an azimuthally periodic flow, the components of the displacement as a function of azimuth in cylindrical coordinates are
where the complex amplitudes are given by Eqs. (A.9)(A.11). The phase difference between the azimuthal displacement and the spherical radial displacement is . The azimuthal perturbation leads to diverging flows around tube crests and converging flows around the troughs of the tube, to restore hydrostatic equilibrium (see Fig. 7, middle panel). Such a flow pattern increases the density deficit in the tube crests, thus it has a destabilising effect. The azimuthal perturbation does not have a significant effect on the radial position of the tube crests, owing to (i) the relatively small ratio ; and (ii) the phase relation between and .
In the present analysis, we assume that any displacement of mass elements along the tube is negligible before the stationary equilibrium is reached. This process involves nonlinear variations in density, internal flow speed, and magnetic field strength as functions of azimuth and time, which are outside the scope of our linear approach. We present the justification of this assumption and discuss its limitations in Appendix B.
Figure 7: Geometry of a flux tube deformed by azimuthally periodic flows, at three initial locations (solid lines): (from top to bottom of the figure) the top, the middle, and the bottom of the overshoot region. The latitude is and the field strength is G. The distance from the middle of the overshoot region in kilometres is shown on the left axis, and in units of local pressure scale height (at ) on the right axis. The azimuthal wavenumber is chosen such that the troughs or crests of the tube partially enter the convection or the radiative zones. The radius of the tube is 1000 km. Arrows denote the relative strength of the resultant perturbation, which is led by the external flow. 

Open with DEXTER 
3.1.2 Parameter study
Figure 8: Variation in the maximal displacement of flux tubes by azimuthally periodic perpendicular flows with four different azimuthal wavenumbers, m, as a function of flow speed ( G, R_{t}=1000 km). The vertical axes denote the radial distance with respect to the middle of the overshoot region (on the right axes in units of the local pressure scale height at the midovershoot region). Solid lines show the maximal downward displacement of the troughs of a tube starting in the bottom of the overshoot region. Dashed lines show the maximal displacement for the troughs and crests of tubes starting in the middle of the overshoot region, and the dashdotted lines show the displacement of the crests of tubes starting near the top of the overshoot region. For clarity, the effect of upward (downward) flows on flux tubes near the bottom ( top) of the overshoot region is not shown. The shaded regions denote the radiative zone below and the convection zone above. 

Open with DEXTER 
Dependence on the flow velocity and the azimuthal wavenumber.
Figure 8 shows the variation in , which is the radial distance of the crest of the perturbed tube from the middle of the overshoot region, as a function of the velocity amplitude of the radial flow for various azimuthal wavenumbers. The minimum flow velocity required to advect the flux tube up to the convection zone or down to the radiative zone boundaries corresponds to the intersections of the curves with the boundaries of the overshoot layer. For a fixed flow speed and increasing m, it becomes more difficult to displace the tube crests from the equilibrium configuration, because the magnetic curvature force increases. For a given m, faster flows lead to larger perturbations. In the case for which the tube extends between the boundaries of the overshoot layer, the variation in buoyancy along the tube can be significant, because the superadiabaticity changes by about three orders of magnitude between the upper and lower boundaries (cf. Sect. 2.3, Fig. 4).Figure 9: Displacement of flux tubes by radial external flows that are periodic in azimuth, as a function of field strength. The left panel is for m=8, and the right panel for m=20, to show the contrast between the effects of flows with large and small wavelengths. The radius, the latitude, the initial depth of the flux tubes, the vertical axes, and the line styles in the plots are the same as in Fig 8. The tube is subjected to an azimuthally periodic flow with an amplitude of 10, 26, 31 m s^{1} for the bottom, middle, and top of the overshoot region, respectively. The long dashed lines indicate the boundaries of the radiation and convection zones. 

Open with DEXTER 
Dependence on the magnetic field strength.
Figure 9 shows the variation in the displacement amplitude as a function of its field strength, for km, and the amplitude of the azimuthally periodic perpendicular flow is assumed to be 14 m s^{1} at the bottom, 26 m s^{1} at the middle, and 31 m s^{1} at the top of the overshoot region, as in Fig. 7. As the azimuthal wavenumber of the flow is increased, the tension force resists the deformation of the tube more strongly, so that the perturbation weakens. For G and , the deformation of the tube by the flow becomes smaller than 200 km. For m=20, a perpendicular flow with a speed of 26 m s^{1}, applied to a tube with B_{0}=10^{5} G, located at the middle of the overshoot region, displaces it by an extent of about (300 km). The corresponding relative change in superadiabaticity between the depths of the crests and the troughs is (see Fig. 5), so the effect is ineffective in this case.
Dependence on the tube radius.
We now set the field strength to G and vary the tube radius, . The resulting functional dependence is shown in Fig. 10 for m=8 and m=20, for tubes starting at the three sets of depths and convective speeds chosen above for Fig. 9. Because the drag force is inversely proportional to R_{t}, thicker tubes are less affected by the flow. For m=20, tubes thicker than about 2000 km in diameter do not reach the boundaries of the overshoot region. At first glance, this indicates that thicker tubes may be stored in the overshoot region for longer times than thinner tubes. However, diameters larger than about a few thousand kilometres are not relevant to the present context, because the thin flux tube approximation is not valid if the tube diameter is comparable to the thickness of the overshoot layer: the radial variation in superadiabaticity across the tube becomes nonnegligible, inducing a strong variation in the buoyancy over the tube crosssection (see Fig. 4).The upper horizontal axis in Fig. 3 shows the perpendicular velocity amplitude of the external flow, leading to a given maximum radial displacement of the flux tube, , for m=5. Using the correspondence between and (Sect. 3.1.1), we can make the following estimation, using Fig. 3. A toroidal flux tube with , G, and km can be stored in the middle of the overshoot region for about 3 years, provided that the external flow velocity does not exceed about . This semianalytical estimate concerning the storage of a magnetic flux tube in the overshoot layer is tested using nonlinear numerical simulations involving a localised downflow, in Sect. 3.2.2.
Figure 10: Same as Fig. 9, for G and varying tube radius. 

Open with DEXTER 
3.2 Localised radial flow
We next consider the effects of both a localised radial flow, which has a finite extension in azimuth and latitude, and a longitudinal flow. The former describes, e.g., a convective downdraft penetrating into the overshoot layer, and the latter is required to define the initial mechanical equilibrium state.
3.2.1 Linear analysis
Before treating the nonlinear evolution of flux tubes under the combined effects of a localised downflow and the longitudinal flow, we develop an analytical approach to help us to understand the basic physics. We consider a toroidal flux tube subject to a localised downflow in the middle of the overshoot region, and describe the downward flow speed by a Gaussian function of the azimuth and obtain the stationary solution in a way similar to the one in Sect. 3.1. The description of the flow field and the stationary solution for the displacement are given in Appendix A.2. Figure 11 shows the radiusazimuth diagrams in which the dashed curves represent the stationary solutions for a flow with m s^{1} (minus sign means that the flow is directed downward) and various field strengths. For simplicity, the flux tube is located in the equatorial plane. For B_{0}=10^{4} G, the downflow leads to a valleyshaped deformation in the flux tube. In this case, the shape of the flux tube is mostly determined by the azimuthal dependence of the perpendicular flow speed. For G and G, magnetic tension increasingly affects the stationary equilibrium shape of the flux tube. For a sufficiently strong magnetic field, a localised downflow leads to a deformation with a larger azimuthal extension than that of the downflow itself, owing to magnetic tension.
Figure 11: The radial position of an equatorial flux tube as a function of the azimuth, as determined by linear analysis ( dashed curves, m s ^{1}), and by nonlinear simulation of the SF50 case, after 8 months ( solid curves). The internal flow (in the rest frame of the external medium) is from left to right. The initial field strength is different in each panel. The initial radial position is cm (the middle of the overshoot region). The dashdotted horizontal lines indicate the boundaries of the overshoot region. 

Open with DEXTER 
3.2.2 Nonlinear simulations
We assume a Gaussian profile for the external flow field,where a(t) describes the time variation. We set cm and km, , , so that the flow speed has its maximum value, , at the upper boundary of the overshoot region, and falls to about one fifth of its maximum in the middle of the overshoot region. The lateral extensions and are free parameters. For all the simulations, the flux tube is located in the middle of the overshoot layer, and its crosssectional radius is 10^{3} km.
Stationary flow (SF)
To test whether the stationary equilibria predicted in Sect. 3.2.1 can occur in the nonlinear case, we carried out numerical simulations by applying a localised downflow, which becomes stationary after a given time, days. We consider a toroidal flux tube, which is initially in mechanical equilibrium. We assume that an external downflow gradually develops around the point , such that the time profile in Eq. (8) is assumed to beWe present two subsets of simulations in this section: for the subset SF50, m s^{1}, and for SF10 10 m s^{1}. For each subset, the free parameters are the latitude, magnetic field strength, and the azimuthal and latitudinal extensions of the downflow, with the condition . We define the ``rise time'' of an unstable flux tube as the time between t=0 and the end of the simulation when the tip of the rising loop arrives at about , where the thin flux tube approximation breaks down.
Snapshots from the set SF50 at t=8 months for various initial field strengths are shown with solid curves in Fig. 11, in comparison with the linear results for stationary equilibrium (Sect. 3.2.1). The direction of the internal flow along the flux tube is from left to right. The stationary solution is nearly identical to the nonlinear simulation snapshot for B_{0}=10^{4} G. The broadening of the submerged part for higher B_{0} (see Sect. 3.2.1) is visible in the numerical solutions. However, the simulations deviate from the linear results substantially on the left side of the downflow region, to an increasing extent for higher B_{0}. The azimuthally symmetric external downflow deforms the flux tube in an asymmetric form, owing to the internal flow. The central part of the tube reaches a dynamical equilibrium after about 100 days, whereas the tube as a whole never comes to equilibrium in any case: for low field strengths ( G), the upward portion creates a transverse wave, which propagates leftward (retrograde). In the course of its subsequent evolution, the wave energy is transferred to lower azimuthal modes, while the submerged portion of the tube persists as long as the downflow speed is kept constant. For higher values of B_{0}, the magnetic tension force limits the extent of the downward displacement, whereas the excess magnetic buoyancy in the left wing is enhanced. This leads to a Parkerunstable loop for G, because the top of the loop reaches layers of sufficiently high superadiabaticity, and the density deficit in the loop top increases with time.
Figure 12: Physical quantities illustrating the effect of the external downflow, as a function of azimuth. The snapshots are taken at t=5 days of the SF50 simulation run for G. The external downflow is centred around , which is indicated by a vertical line. We plot a) the radial displacement of the tube from its initial position, b) the internal flow speed; c) the difference between the internal and external densities; d) the deviation of field strength from B_{0}; e) the internal pressure (dashed line shows the external pressure), and f) the internal density. 

Open with DEXTER 
To explain the physical reasons for the asymmetric deformation, we show in Fig. 12 the azimuthal profiles of various physical variables at an early stage of the simulation (t=5 days), for the SF50 case, where G. Once the external downflow begins to advect a part of the tube downward to a smaller radius (Fig. 12a), angular momentum conservation leads to a local increase in the internal flow speed (b). The maximum internal flow speed is reached at the left wing of the submerged portion of the tube, because of the additional acceleration by the component of gravity along the tube. Through the right side of the submerged part, the internal flow decelerates back to its initial equilibrium speed. This asymmetric flow leads to a density deficit at the left part and a density excess at the right side (c). The resulting positive buoyancy in the left part increases with time and forms an upward moving loop, which is driven further by the accelerating internal flow on the right wing of the loop and the higher superadiabaticity at the top location (see Fig. 11, the bottom panel). Starting from the initial phases, the accelerating flow also leads to a slightly lower internal pressure at the left part (e). To balance the external pressure, the magnetic field strength increases on the left part (d), and this contributes to the subsequent growth of the magnetic buoyancy instability. For G, the instability does not set in, because the internal flow required for the initial mechanical equilibrium is too slow for a sufficient density deficit to develop on the left wing of the submerged part. For G, the transversal wave propagates leftward and quits the downflow region, well before a sufficient density deficit on the left wing develops.
Table 2 presents the rise times (in cases of instability) of flux tubes at different latitudes, field strengths, and lateral extensions of the downflow, for SF50 and SF10 cases. Linearly Parkerunstable cases are shown in italics.
Table 2: Rise times for stationary flow simulations.
For the set SF50, the cases with field strengths up to G are stable. For the linearly Parkerstable cases ( G), greater lateral extensions of the downflow would destabilise the tube, because for larger , the magnetic curvature force is smaller, thus the perturbed part of the tube can descend deeper, increasing further the internal flow speed in the left wing. In this regime, a stronger magnetic field also leads to faster growth of the loop, because the internal equilibrium flow becomes faster, which develops the density deficit on the left wing more rapidly. For the linearly Parkerunstable cases (denoted in italics in Table 2), the rise times are generally longer, because in most cases the initial central deformation of the tube excites standing waves with m=5. This limits the growth of the loop on the left side of the submerged part, and the flux tube enters the convection zone only after the growth of an oscillatory instability with m=1. For SF10 simulations, the downflow is too weak to trigger nonlinear buoyancy instability, so the linearly Parkerstable cases are stable in the nonlinear case, whereas for G Parker instability sets in, in accordance with the prediction of linear perturbation theory.Transient flow (TF)
Downflows penetrating into the solar
overshoot layer are probably transient, i.e., they are decelerated by
the increasingly stable stratification on their way towards the
radiative zone.
To test the consequences of such a flow, we assume a Gaussian
profile for the time variation in the downflow speed in
Eq. (8), such that
where is the time when the maximum flow speed is reached. The flow sets in at t=0 with . The initial effect of the downflow on the flux tube is similar to that in the SF case. However, once the downflow ceases, it leaves behind transversal tube waves propagating along the tube, which interfere with each other. For unstable cases, the relative longitudinal flow leads to instability by frictional coupling, with the growth rate depending on the internal flow velocity (thus the field strength) and the perturbation amplitude (cf. Sect. 2). The difference between SF and TF is that in the former case the stationary flow had continued for a sufficiently long time to allow the formation of a buoyant loop on one side of the downflow region (Fig. 12).
We carried out a systematic survey of numerical simulations by varying the flow duration, , and velocity amplitude, in the ranges m s^{1} and d, for . The constants and were chosen such that the initial speed is of the order of . For G, the growth time of the frictioninduced instability is longer than the upper time limit, which we set to be 7 years. Instabilities with growth times longer than a few years are not relevant in the present context, because within that time the toroidal magnetic field must be amplified in the solar tachocline to its full strength. For G, which is about 20 kG lower than the Parker instability threshold for linear perturbations, the flux tube develops an unstable loop within about 4 years for some values of .
The evolution of the tube for m s^{1} and downflow durations of 60 days (TF60) and 180 days (TF180) is shown in supplementary animations (see the online appendix, Figs. C.2 and C.3). The tube rises over about 10 years for TF60, and about 3.3 years for TF180. In the initial stages of TF180, the evolution of the tube resembles the strongfield behaviour for the SF50 case (cf. Fig. 11, bottom panel).
Figure 13 shows the rise times as a function of and . Only the final month of the rise time is spent in the convection zone proper, the rest being spent in the overshoot region. A relatively fast downflow ( 40 m s^{1}) lasting less than about 2 months, or a slow downflow ( 20 m s^{1}) lasting less than about 6 months allow the storage of the flux tube in the overshoot region for more than about 3 years. As a rough estimate of the flow duration in the overshoot region, we also calculated the convective turnover time, , in the middle of the overshoot layer. To lead to the formation of an unstable flux loop out of a tube with G, either the flow should last much longer than the turnover time or the flow should reach much higher speeds than the range assumed here.
Figure 13: The variation in the rise time for a linearly Parkerstable flux tube in the midovershoot region, as a function of the duration and maximum speed of the localised external downflow with , , Mm (convection zone  overshoot layer boundary) for the TF simulations. The field strength and latitude of the initial tube are G and (cf. Fig. 1). The diamonds represent the ``stable'' cases, for which no emergence takes place before 7 years, and the size of the filled circles is inversely scaled with the rise time. The solid curve shows the convective turnover time as an indicator for the flow duration, estimated as . 

Open with DEXTER 
4 Discussion
We have investigated flowinduced instabilities of toroidal magnetic
flux tubes in the solar overshoot region, to extend the results
of Paper III to the nonlinear regime and to quantitatively
estimate the effects of convective flows on the
storage of magnetic flux in the solar overshoot region.
Our simulations confirm the main results of the linear stability
analysis of
Paper III (Sect. 2.2).
The perpendicular velocity component can be approximated generally
in the form
,
where
is the amplitude of the radial perturbation and
is the eigenfrequency of the fastest growing
unstable wave mode. Substituting this expression into Eq. (3),
we obtain
For a given value of , which is specified for a given set of (depth), , and B_{0}, increases with increasing perturbation amplitude, as shown by numerical simulations in Sect. 2.2 (Fig. 2). Considering the relations obtained in Sect. 3.1, a perpendicular flow with and m=5 in the middle of the overshoot region can displace a flux tube with G to km (cross symbols in Fig. 2). Substituting into Eq. (11), this amounts to , which is of the same order as the value chosen for the full curve in Fig. 2. Considering the numerical results in Sect. 2.2 (see Fig. 3), we conclude that an external flow with a velocity amplitude of 110 m s^{1} perpendicular to the tube axis leads to instability with a growth time of the order of 1000 days (2.74 years) or longer. This range of velocities is consistent with the estimates of van Ballegooijen (1982) for the convective velocities in the overshoot region, supporting the possibility of storing magnetic flux tubes that contain fluxes of the order of 10^{21} Mx. In a more recent study, Brummell et al. (2002) carried out 3D numerical simulations of penetrative convection, which infer penetration depths between and , when they extrapolate the dimensionless numbers in their simulations to solar values. These values are also comparable to the radial perturbation amplitudes that we have assumed here.
The numerical experiments presented in Sect. 2.3 have shown us that the nonlinear instability occurring in the linearly Parkerstable regime is induced mainly by the frictional coupling of the flux tube with its surroundings. We have also found unstable flux tube configurations, for which the effect plays a significant role in the dynamics. These numerical experiments were made without considering the drag force and in the bottom of the overshoot layer, where the radial gradient of superadiabaticity is steeper than in the remainder of the overshoot region. This result may be relevant during the final phases of the decay of the largescale toroidal field in the Sun: if we assume that the toroidal flux at the upper layers of the overshoot region has already been removed to a large extent at this phase, a flux tube that forms near the bottom of the overshoot region can be destabilised rapidly by strong (possibly rare) convective downflows, leading to a few active regions during a solar minimum. We conjecture that the effect is not the main source of flux loss from the overshoot layer, with a possible exception in the bottom of the overshoot region, provided that overshooting convective flows are sufficiently strong.
In test simulations, we have found that for perturbations larger than about (55 km), the growth rates start to deviate from the predictions of linear stability analysis (FerrizMas & Schüssler 1995), because nonlinear effects govern the dynamics and determine the growth rate of instability, through either the frictioninduced instability (Sect. 2.2) or the effect (Sect. 2.3), depending on the depth and the longitudinal flow speed.
In seeking possibilities of a slingshot effect that leads to smallscale loops in thin flux tubes, we have found that this effect is inefficient in removing magnetic flux from the overshoot region. The common result of the simulations is that a loop driven by a transient downflow hits the radiative zone and bounces back rapidly. However, in its way through the overshoot region it is strongly decelerated mainly by friction.
In a parameter study surveying analytical estimates of the displacement amplitude as a function of field strength (Sect. 3.1.2), we have found that a flow pattern with m=20 and m s^{1} acting on a tube with B_{0}=10^{5} G leads to a perturbation of about 300 km. The corresponding difference in the superadiabaticity of the stratification is (see Fig. 5). Therefore, we do not expect azimuthally periodic flows with short wavelengths ( ) to trigger flux tube instabilities for G, provided that the displacement is less than about 300 km, in other words, m s^{1}.
After understanding the fundamental effects of azimuthally periodic flows on flux tubes, we have considered the case of a localised downflow in Sect. 3.2. For a given flow pattern, we have calculated the evolved states of flux tubes with various field strengths corresponding to linearly Parkerstable cases, using the steadystate approximation in conjunction with numerical simulations. The experiments presented in Fig. 11 can also be interpreted in terms of a toroidal field strength increase with time, in the rising phase of solar activity. As the toroidal field is amplified, localised downflows will have stronger disruptive effects on flux tubes, owing to the increasing internal flow speed, which is determined by the mechanical equilibrium condition. For a stationary downflow, we thus suggest that the lower limit to the field strength for which flux loops start entering the convection zone is of the order of G for the middle of the overshoot region.
Proceeding to nonstationary, transient flows, we have set up a survey of simulations for G in the middle of the overshoot region for Gaussian time profiles (TF case). Depending on the dynamical properties of the downflow, we have calculated the evolution of a flux tube with an upper time limit of 7 years. However, in the first 7 years of an activity cycle, the amplification of the largescale toroidal field in the tachocline cannot be neglected. If a flux tube with a subcritical field strength, say G, remains in the overshoot region for a few years, it is likely that it will be amplified and eventually form a Parkerunstable loop for G, which can lead to the formation of active regions with the observed properties in the photosphere.
We have found that the storage of a toroidal magnetic flux tube with Mx for times comparable to the dynamo amplification time is possible. Within the scope of the thin flux tube approximation and nonlocal mixing length models of the solar interior, we have not found any significant (hydrodynamically or magnetically induced) stability problem that impedes the construction of Parkerunstable tubes with fluxes of the order of 10^{22} Mx.
We have considered the problem of storing a toroidal flux tube field rather than the problem of building up the magnetic field in the first place. We have assumed that the field is amplified in the overshoot region, e.g., by rotational shear, on a timescale of a few years. The stationary equilibria of subcritical flux tubes ( G) advected by radial flows indicate deformations comparable to the size of the overshoot region (see Fig. 9). On the one hand, it may be argued that these deformations can destabilise the tube, e.g., by effect, at times comparable to the duration of overshooting convective flows (see Sect. 3.1). However, nonlinear simulations (Fig. 11) show deviations from linear estimates, as the field strength increases towards G: the tube is destabilised as the field strength is increased. These simulations indicate that toroidal magnetic fields can be stored within the overshoot region in the course of their amplification up to the critical field strength (10^{5} G) required to explain general properties of active regions. Our results can be crosschecked when the spatiotemporal structure of penetrative convection near the bottom of the convection zone is observed by helioseismology, or when 3D hydrodynamic simulations with realistic Reynolds numbers and sufficiently high resolution become available.
5 Conclusions
Based on analytical steadystate approximations and nonlinear simulations of toroidal flux tubes in the solar convective overshoot region, we reach the following conclusions: The flow instability driven by the frictional coupling of transversal MHD waves with external flows is nonlinear in the sense that the growth rate is a function of the initial perturbation amplitude. This is consistent with the results of the linear stability analysis of Paper III, in which the parameter is proportional to the speed of the perpendicular flow and thus to the amplitude of the subsequent perturbation. Therefore, the perturbation amplitude determines the parameter in the linear analysis.
 Significant buoyancy variations (effect) along magnetic flux tubes can lead to a nonlinear buoyancy instability. The effect most likely occurs in regions of large radial gradients of superadiabaticity such as the bottom of the overshoot layer, where a sufficient variation along the tube ( ) is easier to attain than in the upper layers.
 For flux tubes in the solar overshoot layer, we have established links between the perpendicular flow velocity amplitude, its spatial (azimuthal) extension, and the resulting radial displacement and instabilities induced by effect or longitudinal flows.
 To store magnetic flux with flux densities between 10^{4}10^{5} G for times comparable to the dynamo amplification time in the convective overshoot layer, the average flow speed and the flow duration must not exceed about 50 m s ^{1} and 100 days, respectively, and that the azimuthal extension of the flow is lower than about . If we assume that these conditions prevail in the solar overshoot region, then magnetic fluxes of up to 10^{22} Mx can be stored within thin flux tubes during the dynamo amplification phase.
The authors are grateful to Manfred Schüssler for useful discussions and suggestions in the course of the study, and acknowledge the referee for the suggestions, which helped in improving the manuscript.
Appendix A: Flux tube subject to radial flows
The equations for a flux ring at an arbitrary latitude , in the limit of , are given by FerrizMas & Schüssler (1995) (Sect. 4.1). Consider that a flow along the spherical radial direction with an amplitude , is applied to a flux tube with radius located at latitude . This flow will exert a drag force perpendicular to the tube axis, , where is the unperturbed external density. After a certain time, the drag force and the restoring forces of magnetic tension, Coriolis, and buoyancy will balance each other and a stationary equilibrium will be reached. Adopting the equation of motion for linearised perturbations (FerrizMas & Schüssler 1995) for the stationary force equilibrium yields the inhomogeneous system of equations
where
(A.4)  
(A.5) 
Here we have assumed that the magnitude of the drag force is comparable to a firstorder perturbation.
A.1 Azimuthally periodic flow
We define the flow velocity in the form
for a given
azimuthal wavenumber m. The resulting drag force is given by
where is the maximum value of the velocity. The resulting displacement is of the form
Substituting Eqs. (A.6 and A.7) into Eqs. (A.1)(A.3), we obtain
This equation can be solved in a straightforward way using Cramer's rule. The solution for the complex amplitudes yields
Substituting each component (A.9)(A.11) into Eq. (A.7), we obtain
A.2 Localised flow
In the case of a localised perpendicular flow, we define the velocity
and the drag force in the following forms:
where is the Fourier transform of the function . For the resulting displacement, , we follow the method described in Sect. A.1 and write the components in the Fourier series representation
where the complex constants are given by Eqs. (A.9)(A.11). The final solution for the components of the perturbation resulting from a localised perpendicular flow has the form
Appendix B: The stationary equilibrium approximation
B.1 Comparison of azimuthal and radial effects
Suppose that a toroidal flux tube is in mechanical equilibrium.
Then apply a localised external flow in the spherical radial direction.
The parts of the tube affected by the flow are advected
in the direction of the flow, until the drag force is balanced by the
restoring forces of buoyancy and magnetic tension.
To justify ignoring any displacement of mass elements
along the tube, a necessary condition is
that the time for the radial displacement, ,
is short compared
to that of the displacement along the tube axis
(tangential perturbation, ), so that
where r_{0} is the radial coordinate of the equilibrium tube, is the propagation velocity of the perturbation along the tube, and is the maximum flow speed, and can be determined from the solution of the eigenvector problem using the linearised equation of motion for firstorder perturbations (see Sect. B.2). The lefthand side of Eq. (B.1) represents a timescale for the reestablishment of hydrostatic equilibrium, and the righthand side represents the time it takes for a dynamical balance between the drag force and the restoring forces is reached.
As an example, we consider an equatorial flux tube, in the middle of the overshoot layer. The external flow is assumed to be localised (Eq. (A.14)) and downwards with a maximum speed of 10 m s^{1}. This is the same flow configuration as that described in Sects. 3.2.1 and A.2. For simplicity, we disregard the effects of longitudinal wave modes and the presence of azimuthal flow, which is required by the mechanical equilibrium condition. Table B.1 shows the values for both sides of the inequality (B.1), as a function of the azimuthal wavenumber of the external flow, for G. We also give the values of the downward displacement amplitude, , produced by the external downflow. For comparison, numerical simulations were performed by setting up azimuthally periodic flow configurations for various azimuthal wavenumbers. In each case, the flow sets in at t=0 and reaches its full strength ( m s^{1}) in 30 days. The tabulated values are the maximum displacement at t=240 days, for which the drag force exerted by the external flow is already balanced by buoyancy and magnetic tension. The variation in displacement amplitude as a function of m for azimuthally periodic flows shows general agreement with that for the linearly estimated amplitudes of each azimuthal wave mode corresponding to a localised downflow.
The ratio decreases for smaller m. When becomes comparable to (as the ratio falls below about 10), the azimuthal perturbation might have a significant effect on the density distribution along the tube, such that our assumption of a stationary equilibrium between the normal forces becomes inaccurate.
Table B.1: Characteristic times for radial and azimuthal displacements.
B.2 Flux tube in the equatorial plane: eigenvector problem
For simplicity, we consider a flux tube in the equatorial plane, to obtain eigenvectors for the wavelike solutions of the perturbed flux tube. In this case, the momentum equation for linear perturbations is simpler than for a nonequatorial flux tube, and therefore it is simpler to make a physical interpretation.
To obtain relations between the flow velocity and the resulting
displacement, we consider the linearised equations of motion for
perturbations
,
where R is the distance from the rotation
axis,
is the longitude, and
is the latitude. Here, only radial
and azimuthal components (perturbations in the plane) are considered,
because the equation for the latitudinal component is decoupled in the equatorial
plane. Following Ferriz Mas & Schüssler (1993), the linearised equations
of motion for the perturbations in the plane read
where
The meanings of the various symbols in the above equations are as follows: is the local pressure scale height, is the superadiabaticity, p_{i0} and are the equilibrium values for the pressure and density inside the tube, respectively, r_{0} is the spherical radial position of the equilibrium tube, and are, respectively, the angular velocities of the plasma in the tube and the external medium, is the Alfvén velocity, B_{0} is the magnetic field strength, and g_{0} the gravitational acceleration with a radial dependence , where is taken to be 2, which means thatthe small contribution of the mass inside the convection zone is neglected. The substitution of the ansatz
into Eqs. (B.2)(B.3) leads to the system of equations
For the existence of nontrivial solutions, the characteristic
determinant must vanish. This yields the dispersion relation,
which is a 4thorder polynomial in terms of the eigenfrequency.
The eigenvector in the azimuthal direction is determined by one of the
equations in the system
(B.9B.10):
from which we obtain, by using the ansatz Eq. (B.8), the velocity
of the azimuthal perturbation, as an estimate of the rate at which the hydrostatic equilibrium along the perturbed flux tube is established. We use this result in evaluating Eq. (B.1), to check for the validity of the assumption of stationary equilibrium of a flux tube being displaced by an external flow.
References
 Batchelor, G. K. 1967, Fluid dynamics (Cambridge, England: Cambridge University Press)
 Brummell, N. H., Clune, T. L., & Toomre, J. 2002, ApJ, 570, 825 [NASA ADS] [CrossRef]
 Caligari, P., MorenoInsertis, F., & Schüssler, M. 1995, ApJ, 441, 886 [NASA ADS] [CrossRef]
 D'Silva, S., & Choudhuri, A. R. 1993, A&A, 272, 621 [NASA ADS]
 Fan, Y. 2001, ApJ, 546, 509 [NASA ADS] [CrossRef]
 Fan, Y., Fisher, G. H., & McClymont, A. N. 1994, ApJ, 436, 907 [NASA ADS] [CrossRef]
 FerrizMas, A., & Steiner, O. 2007, Sol. Phys., 246, 31 [NASA ADS] [CrossRef]
 FerrizMas, A., & Schüssler, M. 1993, Geophysical and Astrophysical Fluid Dynamics, 72, 209 [NASA ADS] [CrossRef]
 FerrizMas, A., & Schüssler, M. 1995, Geophysical and Astrophysical Fluid Dynamics, 81, 233 [NASA ADS] [CrossRef]
 Holzwarth, V. 2008, A&A, 485, 351 [NASA ADS] [EDP Sciences] [CrossRef]
 Holzwarth, V., Schmitt, D., & Schüssler, M. 2007, A&A, 469, 11 [NASA ADS] [EDP Sciences] [CrossRef]
 MorenoInsertis, F. 1986, A&A, 166, 291 [NASA ADS]
 MorenoInsertis, F., Schüssler, M., & FerrizMas, A. 1992, A&A, 264, 686 [NASA ADS]
 Schüssler, M. 1984, A&A, 140, 453 [NASA ADS]
 Schüssler, M. 1996, in Solar and Astrophysical Magnetohydrodynamic Flows, 1737
 Schüssler, M., & FerrizMas, A. 2003, in Advances in Nonlinear Dynamos, ed. A. FerrizMas, & M. Núñez (London/New York: Taylor & Francis), 123
 Schüssler, M., & Ferriz Mas, A. 2007, A&A, 463, 23 [NASA ADS] [EDP Sciences] [CrossRef]
 Schüssler, M., & Rempel, M. 2002, in From Solar Min to Max: Half a Solar Cycle with SOHO, ed. A. Wilson, ESA SP, 508, 499
 Shaviv, G., & Salpeter, E. E. 1973, ApJ, 184, 191 [NASA ADS] [CrossRef]
 Skaley, D., & Stix, M. 1991, A&A, 241, 227 [NASA ADS]
 Solanki, S. K., Inhester, B., & Schüssler, M. 2006, Rep. Progr. Phys., 69, 563 [NASA ADS] [CrossRef]
 Spruit, H. C. 1981, A&A, 98, 155 [NASA ADS]
 Tobias, S. M., Brummell, N. H., Clune, T. L., & Toomre, J. 2001, ApJ, 549, 1183 [NASA ADS] [CrossRef]
 van Ballegooijen, A. A. 1982, A&A, 113, 99 [NASA ADS]
Online Material
Appendix C: Animations
We provide three animated GIF files showing the simulated evolution
of perturbed flux tubes, available online:
frict_inst.gif
TF180.gif
TF60.gif
The animation frict_inst.gif shows the final phases of the development of the frictioninduced instability in the overshoot region, between t=700 d and t=1100 d (Sect. 2.2). A snapshot from the animation is shown in Fig. C.1. The asterisk signs represent selected mass elements, which all move rightwards in the initial phases, owing to the internal equilibrium flow. We assume that G, , and the initial radial perturbation amplitude is km. The initial location is the middle of the overshoot region.
Figure C.1: The shape of a flux tube with G, and km at t=1030 d, in the plane. The horizontal lines mark the boundaries of the overshoot region. 

Open with DEXTER 
The animation files TF60.gif and TF180.gif show the evolution of a flux tube subject to a radial downflow with a duration of 60 and 120 days, respectively, until t=300 days (Sect. 3.2.2). Two snapshots are shown in Figs. C.2 and C.3, corresponding to the time when the downflow ceases in each case, i.e., at t=60 d and t=180 d. Note that for the longerduration flow, the flux tube is radially more disturbed when the action of the downflow is finished.
Figure C.2: The shape of a flux tube with G, , at t=60 d, when the transient downflow ( m s^{1}) ceases. The horizontal lines mark the boundaries of the overshoot region and the vertical line denote the azimuthal location of the centre of the downflow. 

Open with DEXTER 
Figure C.3: Same as Fig. C.2, at t=180 d, when the transient downflow ceases. 

Open with DEXTER 
Footnotes
 ... region^{}
 Appendix C and movies are only available in electronic form at http://www.aanda.org
 ... rotation^{}
 Throughout the paper, the index ``0'' refers to quantities pertaining to the mechanical equilibrium state.
All Tables
Table 1: Growth times for various finite perturbations.
Table 2: Rise times for stationary flow simulations.
Table B.1: Characteristic times for radial and azimuthal displacements.
All Figures
Figure 1: Stability diagram for toroidal flux tubes located in the middle of the solar overshoot region. The shaded areas indicate unstable configurations with the grey scale representing the growth rate. Dotted lines outline regions of Parker instability. 

Open with DEXTER  
In the text 
Figure 2: Comparison of the analytical solutions for the growth rate of the frictioninduced instability (lines) with the numerical simulations (symbols) for flux tubes at the middle of the overshoot region and latitude. The initial perturbation amplitude (), and the Stokes drag coefficient () are indicated in the legend. 

Open with DEXTER  
In the text 
Figure 3: Growth time as a function of the displacement amplitude, G, , at the middle of the overshoot region. The upper horizontal axis shows the flow velocity amplitude that leads to the displacement according to the analytical approximation discussed in Sect. 3.1, for m=5 and m s^{1}. 

Open with DEXTER  
In the text 
Figure 4: Radial profile of the superadiabaticity in the outer 40 per cent of the Sun, according to the solar model used (Skaley & Stix 1991). The arrows and the corresponding circles mark the boundaries of the overshoot region. 

Open with DEXTER  
In the text 
Figure 5: The difference in superadiabaticity between the surroundings of the top and bottom parts of perturbed flux tubes, which are initially located at the middle overshoot zone with and G, the same initial conditions as in Fig. 3. The drag force has been taken into account. The numbers indicate growth times of the instability for each simulation. 

Open with DEXTER  
In the text 
Figure 6: Same as Fig. 5, for G, , at the bottom of the overshoot region, without the drag force. When the initial perturbation amplitude exceeds 10^{3} km, the flux tube becomes unstable owing to buoyancy differences along the tube. 

Open with DEXTER  
In the text 
Figure 7: Geometry of a flux tube deformed by azimuthally periodic flows, at three initial locations (solid lines): (from top to bottom of the figure) the top, the middle, and the bottom of the overshoot region. The latitude is and the field strength is G. The distance from the middle of the overshoot region in kilometres is shown on the left axis, and in units of local pressure scale height (at ) on the right axis. The azimuthal wavenumber is chosen such that the troughs or crests of the tube partially enter the convection or the radiative zones. The radius of the tube is 1000 km. Arrows denote the relative strength of the resultant perturbation, which is led by the external flow. 

Open with DEXTER  
In the text 
Figure 8: Variation in the maximal displacement of flux tubes by azimuthally periodic perpendicular flows with four different azimuthal wavenumbers, m, as a function of flow speed ( G, R_{t}=1000 km). The vertical axes denote the radial distance with respect to the middle of the overshoot region (on the right axes in units of the local pressure scale height at the midovershoot region). Solid lines show the maximal downward displacement of the troughs of a tube starting in the bottom of the overshoot region. Dashed lines show the maximal displacement for the troughs and crests of tubes starting in the middle of the overshoot region, and the dashdotted lines show the displacement of the crests of tubes starting near the top of the overshoot region. For clarity, the effect of upward (downward) flows on flux tubes near the bottom ( top) of the overshoot region is not shown. The shaded regions denote the radiative zone below and the convection zone above. 

Open with DEXTER  
In the text 
Figure 9: Displacement of flux tubes by radial external flows that are periodic in azimuth, as a function of field strength. The left panel is for m=8, and the right panel for m=20, to show the contrast between the effects of flows with large and small wavelengths. The radius, the latitude, the initial depth of the flux tubes, the vertical axes, and the line styles in the plots are the same as in Fig 8. The tube is subjected to an azimuthally periodic flow with an amplitude of 10, 26, 31 m s^{1} for the bottom, middle, and top of the overshoot region, respectively. The long dashed lines indicate the boundaries of the radiation and convection zones. 

Open with DEXTER  
In the text 
Figure 10: Same as Fig. 9, for G and varying tube radius. 

Open with DEXTER  
In the text 
Figure 11: The radial position of an equatorial flux tube as a function of the azimuth, as determined by linear analysis ( dashed curves, m s ^{1}), and by nonlinear simulation of the SF50 case, after 8 months ( solid curves). The internal flow (in the rest frame of the external medium) is from left to right. The initial field strength is different in each panel. The initial radial position is cm (the middle of the overshoot region). The dashdotted horizontal lines indicate the boundaries of the overshoot region. 

Open with DEXTER  
In the text 
Figure 12: Physical quantities illustrating the effect of the external downflow, as a function of azimuth. The snapshots are taken at t=5 days of the SF50 simulation run for G. The external downflow is centred around , which is indicated by a vertical line. We plot a) the radial displacement of the tube from its initial position, b) the internal flow speed; c) the difference between the internal and external densities; d) the deviation of field strength from B_{0}; e) the internal pressure (dashed line shows the external pressure), and f) the internal density. 

Open with DEXTER  
In the text 
Figure 13: The variation in the rise time for a linearly Parkerstable flux tube in the midovershoot region, as a function of the duration and maximum speed of the localised external downflow with , , Mm (convection zone  overshoot layer boundary) for the TF simulations. The field strength and latitude of the initial tube are G and (cf. Fig. 1). The diamonds represent the ``stable'' cases, for which no emergence takes place before 7 years, and the size of the filled circles is inversely scaled with the rise time. The solid curve shows the convective turnover time as an indicator for the flow duration, estimated as . 

Open with DEXTER  
In the text 
Figure C.1: The shape of a flux tube with G, and km at t=1030 d, in the plane. The horizontal lines mark the boundaries of the overshoot region. 

Open with DEXTER  
In the text 
Figure C.2: The shape of a flux tube with G, , at t=60 d, when the transient downflow ( m s^{1}) ceases. The horizontal lines mark the boundaries of the overshoot region and the vertical line denote the azimuthal location of the centre of the downflow. 

Open with DEXTER  
In the text 
Figure C.3: Same as Fig. C.2, at t=180 d, when the transient downflow ceases. 

Open with DEXTER  
In the text 
Copyright ESO 2009