Stability of toroidal magnetic fields in the solar tachocline and beneath

Stability of toroidal magnetic field in a stellar radiation zone is considered for the cases of uniform and differential rotation. In the rigidly rotating radiative core shortly below the tachocline, the critical magnetic field for instability is about 600 G. The unstable disturbances for slightly supercritical fields have short radial scales ∼1 Mm. Radial mixing produced by the instability is estimated to conclude that the internal field of the sun can exceed the critical value of 600 G only marginally. Otherwise, the mixing is too strong and not compatible with the observed lithium abundance. Analysis of joint instability of differential rotation and toroidal field leads to the conclusion that axisymmetric models of the laminar solar tachocline are stable to nonaxisymmetric disturbances. The question of whether sun-like stars can posses tachoclines is addressed with positive answer for stars with rotation periods shorter than about two months. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-axisymmetric solar magnetic fields

A simple non-linear, non-axisymmetric mean field dynamo model is applied to a differentially rotating spherical shell. Two approximations are used for the angular velocity, to represent what is now believed to be the solar rotation law. In each case, stable solutions are found which possess a small non-axisymmetric field component. Although the model has a number of obvious shortcomings, it may be relevant to the problem of the solar active longitudes.     

Force-free magnetic fields in the solar atmosphere

Reliable measurements of the solar magnetic field are restricted to the level of the photosphere. For about half a century attempts have been made to calculate the field in the layers above the photosphere, i.e. in the chromosphere and in the corona, from the measured photospheric field. The procedure is known as magnetic field extrapolation. In the superphotospheric parts of active regions the magnetic field is approximately force-free, i.e. electric currents are aligned with the magnetic field. The practical application to solar active regions has been largely confined to constant-α or linear force-free fields, with a spatially constant ratio, α, between the electric current and the magnetic field. We review results obtained from extrapolations with constant-α force-free fields, in particular on magnetic topologies favourable for flares and on magnetic and current helicities. Presently, different methods are being developed to calculate non-constant-α or nonlinear force-free fields from photospheric vector magnetograms. We also briefly discuss these methods and present a comparison of a linear and a nonlinear force-free magnetic field extrapolation applied to the same photospheric boundary data. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamics of the solar tachocline – I. An incompressible study

Gough & McIntyre have suggested that the dynamics of the solar tachocline are dominated by the advection–diffusion balance between the differential rotation, a large-scale primordial field and baroclinicly driven meridional motions. This paper presents the first part of a study of the tachocline, in which a model of the rotation profile below the convection zone is constructed along the lines suggested by Gough & McIntyre and solved numerically. In this first part, a reduced model of the tachocline is derived in which the effects of compressibility and energy transport on the system are neglected; the meridional motions are driven instead by Ekman–Hartmann pumping. Through this simplification, the interaction of the fluid flow and the magnetic field can be isolated and is studied through non-linear numerical analysis for various field strengths and diffusivities. It is shown that there exists only a narrow range of magnetic field strengths for which the system can achieve a nearly uniform rotation. The results are discussed with respect to observations and to the limitations of this initial approach. A following paper combines the effects of realistic baroclinic driving and stratification with a model that closely follows the lines of work of Gough & McIntyre.     

Background magnetic fields during last three cycles of solar activity

This paper describes our studies of evolution of the solar magnetic field with different sign and field strength in the range from –100 G to 100 G. The structure and evolution of large‐scale magnetic fields on the Sun during the last 3 cycles of solar activity is investigated using magnetograph data from the Kitt Peak Solar Observatory. This analysis reveals two groups of the large‐scale magnetic fields evolving differently during the cycles. The first group is represented by relatively weak background fields, and is best observed in the range of 3–10 Gauss. The second group is represented by stronger fields of 75–100 Gauss. The spatial and temporal properties of these groups are described and compared with the total magnetic flux. It is shown that the anomalous behaviour of the total flux during the last cycle can be found only in the second group. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A dynamo model for axisymmetric and non-axisymmetric solar magnetic fields

More and more observations are showing a relatively weak, but persistent, non-axisymmetric magnetic field co-existing with the dominant axisymmetric field on the Sun. Its existence indicates that the non-axisymmetric magnetic field plays an important role in the origin of solar activity. A linear non-axisymmetric  α²– Ω  dynamo model is derived to explore the characteristics of the axisymmetric  ( m = 0)  and the first non-axisymmetric  ( m = 1)  modes and to provide a theoretical basis with which to explain the 'active longitude', 'flip-flop' and other non-axisymmetric phenomena. The model consists of an updated solar internal differential rotation, a turbulent diffusivity varying with depth, and an α-effect working at the tachocline in a rotating spherical system. The difference between the  α²–Ω  and the  α–Ω  models and the conditions that favour the non-axisymmetric modes under solar-like parameters are also presented.     

Latitudinal shear instability in the solar tachocline

In an attempt to explain the observed rotation profile in the solar radiative zone and the tachocline, Spiegel & Zahn proposed a model based on anisotropic turbulent angular momentum transport. Although very successful in reproducing some of the features of the solar tachocline, their model assumes without verification that the origin of the turbulence could be caused by latitudinal shear instability. This paper studies the weakly non-linear evolution of two-dimensional shear instability, in which the interaction between the global rotation profile and the Reynolds stresses can be described self-consistently. Provided that the initial rotation profile is sufficiently close to marginal stability (which is the case of the solar tachocline), the instability is shown to saturate and to relax to a marginally stable state, which differs very little from the observed rotation profile. It is therefore likely that the tachocline is in a state of marginal stability with respect to latitudinal shear instability, and shows that angular momentum transport in the tachocline is unlikely to be caused by shear-induced turbulence.     

Differential rotation of solar magnetic fields

The connection of the differential rotation of solar magnetic fields with the field sign and strength is studied. The synoptic maps of magnetic fields over the last three solar cycles taken at the Kitt Peak Observatory served as input data for the study. The algorithm of magnetic field filtering over 14 chosen strengt intervals and successive 5-degree latitude zones was applied to these data. The Fourier transform of the time series obtained was then used. Analysis of the power spectra led to the conclusion that there are two types of magnetic fields. These differ in strength (0–50 and 50–700 G) and rotation characteristics. The rotation differentiality for strong magnetic field is almost twice as large as that for weak magnetic fields.     

Dynamo action in simulations of penetrative solar convection with an imposed tachocline

We summarize new and continuing three-dimensional spherical shell simulations of dynamo action by convection allowed to penetrate downward into a tachocline of rotational shear. The inclusion of an imposed tachocline allows us to examine several processes believed to be essential in the operation of the global solar dynamo, including differential rotation, magnetic pumping, and the stretching and organization of fields within the tachocline. In the stably stratified core, our simulations reveal that strong axisymmetric magnetic fields (of ∼ 3000 G strength) can be built, and that those fields generally exhibit a striking antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. In the convection zone above, fluctuating fields dominate over weaker mean fields. New calculations indicate that the tendency toward toroidal fields of antisymmetric parity is relatively insensitive to initial magnetic field configurations; they also reveal that on decade-long timescales, the magnetic fields can briefly enter (and subsequently emerge from) states of symmetric parity.We have not yet observed any overall reversals of the field polarity, nor systematic latitudinal propagation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Fractal properties of solar magnetic fields

We study the spatial properties of solar magnetic fields using data from the Solar Vector Magnetograph of the Marshall Space Flight Center (MSFC) (FeI 5250.2 Å) and SOHO/MDI longitudinal magnetic field measurements (Ni 6767.8 Å) (96-min full-disk maps). Our study is focused on two objects: the fractal properties of sunspots and the fractal properties of the spatial magnetic field distribution of active and quiet regions considered as global structures. To study the spatial structure of sunspots, we use a well-known method of determining the fractal dimension based on an analysis of the perimeter—area relation. To analyze the fractal properties of the spatial magnetic field distribution over the solar surface, we use a technique developed by Higuchi. We have revealed the existence of three families of self-similar contours corresponding to the sunspot umbra, penumbra, and adjacent photosphere. The fractal coefficient has maxima near the umbra—penumbra and penumbra—photosphere boundaries. The fractal dependences of the longitudinal and transverse magnetic field distributions are similar, but the fractal numbers themselves for the transverse fields are larger than those for the longitudinal fields approximately by a factor of 1.5. The fractal numbers decrease with increasing mean magnetic field strength, implying that the magnetic field distribution is more regular in active regions.     

Constraints on angular momentum transport in the Sun from simulations of the tachocline

I present a new numerical tool for studying the interaction of meridional flows and magnetic fields, and study their role in establishing angular-momentum balance in the solar radiative zone. Quantitative comparisons with helioseismic observations provide stringent constraints on existing models of the dynamics of the solar interior. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An interpretation of rapid changes in the magnetic field associated with solar flares

The energy source of a flare is the magnetic field in the corona. A topological model of the magnetic field is used here for interpreting the recently discovered drastic changes in magnetic field associated with solar flares. The following observational results are self‐consistently explained: (1) the transverse field strength decreases at outer part of active regions and increases significantly in their centers; (2) the center‐of‐mass positions of opposite magnetic polarities converge towards the magnetic neutral line just after flares onset; (3) the magnetic flux of active regions decreases steadily during the course of flares. For X‐class flares, almost 50% events show such changes. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Filaments and large-scale magnetic fields

We have used Stanford magnetic field maps to construct distributions of longitudinal magnetic field gradients in the neighbourhood of polarity inversion lines. The distributions were constructed with proper account of the type of the polarity inversion lines and of the existence or absence of dark filaments above them. It is shown that for polarity inversion lines that pass inside active regions or on their boundary, grad B_II distributions for portions of the lines with persisting filament are shifted toward lower values of gradient as compared with grad B_II distributions for portions of the lines without filaments. The influence of the spatial resolution of the magnetograms upon polarity inversion line characteristics is discussed.     

Nonlinear simulations of magnetic instabilities in stellar radiation zones: The role of rotation and shear

Using the 3-dimensional ASH code, we have studied numerically the instabilities that occur in stellar radiation zones in presence of large-scale magnetic fields, rotation and large-scale shear. We confirm that some configurations are linearly unstable, as predicted by Tayler and collaborators, and we determine the saturation level of the instability. We find that rotation modifies the peak of the most unstable wave number of the poloidal instability but not its growth rate as much as in the case of the m = 1 toroidal instability for which it is changed to σ = /Ω. Further in the case with rotation and shear, we found no sign of the dynamo mechanism suggested recently by Spruit even though we possess the essential ingredients (Tayler's m = 1 instability and a large scale shear) supposedly at work. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Changes in the photospheric magnetic field associated with solar flares

The flare-related, persistent and abrupt changes in the photospheric magnetic field have been reported by many authors during recent years. These bewildering observational results pose a challenge to the current flare theories in which the photospheric magnetic field usually remains unchanged in the eruption. In this paper, changes in the photosphere magnetic field during the solar eruption are investigated based on the catastrophe model. The results indicate that the projection effect is an important source that yields the change in the observed photospheric magnetic field in the line-of-sight. Furthermore one may observe the change in the normal component of magnetic field if the spectrum line used to measure the photospheric magnetic field does not exactly come from the photospheric surface. Our results also show that the significance of selecting the correct spectral lines to study the photospheric field becomes more apparent for the magnetic configurations with complex boundary condition (or background field).     

A synoptic program for large solar telescopes: Cyclic variation of turbulent magnetic fields

Upcoming large solar telescopes will offer the possibility of unprecedented high resolution observations. However, during periods of non‐ideal seeing such measurements are impossible and alternative programs should be considered to best use the available observing time. We present a synoptic program, currently carried out at the Istituto Ricerche Solari Locarno (IRSOL), to monitor turbulent magnetic fields employing the differential Hanle effect in atomic and molecular lines. This program can be easily adapted for the use at large telescopes exploring new science goals, nowadays impossible to achieve with smaller telescopes. The current, interesting scientific results prove that such programs are worthwhile to be continued and expanded in the future. We calculate the approximately achievable spatial resolution at a large telescope like ATST for polarimetric measurements with a noise level below 5 × 10^‐5 and a temporal resolution which is sufficient to explore variations on the granular scale. We show that it would be important to optimize the system for maximal photon throughput and to install a high‐speed camera system to be able to study turbulent magnetic fields with unprecedented accuracy (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)