On the twist of emerging flux loops in the solar convection zone

Abstract

We perform numerical simulations of emerging flux loops in the solar convective envelope based on a weakly twisted thin flux tube model recently derived by Longcope and Klapper, generalizing the original formulation for the dynamics of untwisted thin flux tubes by Spruit. The generalized formulation includes the description of torsional Alfvén waves and takes into account the coupling of the writhing motion of the tube axis to the change in the flux tube twist based on the requirement of global helicity conservation of the closed thin flux tube. In this model, the twist of the thin flux tube is described by a quantity q defined as the angular rate of field-line rotation about the tube axis per unit length of the tube. We examine the evolution of twist along Ω-shaped emerging flux loops which are formed as a result of the non-linear growth of the Parker instability of toroidal magnetic flux tubes at the base of the solar convection zone. We find that: (1) In the northern hemisphere, a left-handed twist is generated in the flux tubes as a result of the right-handed tilt or writhe of the emerging loops induced by the Coriolis force. The generated twist increases with the latitude of emergence over the range from 0° to about 38° latitude, but then decreases when the emerging latitude exceeds 38°, because of a change in the preferred eruption pattern. The magnitude of the generated twist q is very small, ≲ 2 × 10-4 rad Mm-1, more than an order of magnitude smaller than the observed amplitude of twist (∼ 0.01 rad Mm-1) in solar active regions. (2) For a toroidal flux ring with a uniform initial twist q0 along the ring, the twist amplitude |q| at the apex of the emerging loop decrease by a factor of about 0.67 because of the stretching of the loop, as it rises from the base of the convection zone to about 20 Mm below the photosphere, at which depth the flux tube can no longer be considered thin. However, because of the more rapid increase of the tube cross-sectional radius a with height, |qa|, which corresponds to the ratio between the azimuthal field to the axial field Bθ/Bl of the tube, increases by a factor of about 2.5 at the apex of the loop, as it rises over the same distance. (3) Because of the effect of the Coriolis force, the distribution of twist along the emerging loop is asymmetric between the leading (in the direction of rotation) and the following sides of the loop. Both |q| and |qa| are greater at the following side than the leading at any depth. Based on the evolution of twist along emerging flux loops, we discuss possible constraints on the twist q0 of initial toroidal flux tubes at the base of the convection zone.