Helicity of solar active regions from a dynamo model

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

Towards a model for the solar dynamo

We study a mean field model of the solar dynamo, in which the non-linearity is the action of the azimuthal component of the Lorentz force of the dynamo-generated magnetic field on the angular velocity. The underlying zero-order angular velocity is consistent with recent determinations of the solar rotation law, and the form of the alpha effect is chosen so as to give a plausible butterfly diagram. For small Prandtl numbers we find regular, intermittent and apparently chaotic behaviour, depending on the size of the alpha coefficient. For certain parameters, the intermittency displays some of the characteristics believed to be associated with the Maunder minimum. We thus believe that we are capturing some features of the solar dynamo.     

Solar cycle according to mean magnetic field data

To investigate the shape of the solar cycle, we have performed a wavelet analysis of the large–scale magnetic field data for 1960–2000 for several latitudinal belts and have isolated the following quasi-periodic components: ∼22, 7 and 2 yr. The main 22-yr oscillation dominates all latitudinal belts except the latitudes of ±30° from the equator. The butterfly diagram for the nominal 22-yr oscillation shows a standing dipole wave in the low-latitude domain  (∣θ∣≤ 30°)  and another wave in the sub-polar domain  (∣θ∣≥ 35°)  , which migrates slowly polewards. The phase shift between these waves is about π. The nominal 7-yr oscillation yields a butterfly diagram with two domains. In the low-latitude domain  (∣θ∣≤ 35°)  , the dipole wave propagates equatorwards and in the sub-polar region, polewards. The nominal 2-yr oscillation is much more chaotic than the other two modes; however the waves propagate polewards whenever they can be isolated.
We conclude that the shape of the solar cycle inferred from the large-scale magnetic field data differs significantly from that inferred from sunspot data. Obviously, the dynamo models for a solar cycle must be generalized to include large-scale magnetic field data. We believe that sunspot data give adequate information concerning the magnetic field configuration deep inside the convection zone (say, in overshoot later), while the large-scale magnetic field is strongly affected by meridional circulation in its upper layer. This interpretation suggests that the poloidal magnetic field is affected by the polewards meridional circulation, whose velocity is comparable with that of the dynamo wave in the overshoot layer. The 7- and 2-yr oscillations could be explained as a contribution of two sub-critical dynamo modes with the corresponding frequencies.     

Evidence for a long-term variation of the dynamo action responsible for the solar magnetic cycle

By using the sunspot time series as a proxy, we have made a detailed analysis of the mean solar magnetic field over the last two and half centuries, by means of a reconstruction of its phase space. We find evidence of a long-term trend variation of some of the solar physical processes (over a few decades) that might be responsible for the apparent erratic behaviour of the solar magnetic cycle. The analysis is done by means of a careful study of the axisymmetric dynamo model equations, where we show that the temporal counterpart of the magnetic field can be described by a self-regulated two-dimensional dynamic system, usually known as a Van der Pol–Duffing oscillator. Our results suggest that during the last two and half centuries, the velocity of the meridional flow, v _p, and the efficiency of the α mechanism responsible for the conversion of toroidal magnetic field into poloidal magnetic field might have suffered variations that can explain the observed variability in the solar cycle.     

The solar dynamo in the light of the distribution of various sunspot magnetic classes over butterfly diagram

We present the data concerning the distribution of various sunspot magnetic classes over the solar butterfly diagram and discuss how this data can inform solar dynamo models. We use the statistics of sunspots that violate the Hale polarity law to estimate the ratio of the fluctuating and mean components of the toroidal magnetic field inside the solar convective zone. An analysis of the spatial distribution of bipolar, unipolar and complex sunspot groups in the context of simple dynamo models results in the conclusion that the mean toroidal field is relatively simple and maintains its shape during the course of the solar cycle (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Flux-dominated solar dynamo model with a thin shear layer

Flux-dominated solar dynamo models have demonstrated to reproduce the main features of the large scale solar magnetic cycle, however the use of a solar like differential rotation profile implies in the the formation of strong toroidal magnetic fields at high latitudes where they are not observed. In this work, we invoke the hypothesis of a thin-width tachocline in order to confine the high-latitude toroidal magnetic fields to a small area below the overshoot layer, thus avoiding its influence on a Babcock-Leighton type dynamo process. Our results favor a dynamo operating inside the convection zone with a tachocline that essentially works as a storage region when it coincides with the overshoot layer. (© 2007 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.     

Mode enslavement in a two-layer stellar dynamo

We consider a conventional stellar  α²ω  -dynamo with dynamo generators localized in two spherical shells separated by a passive layer. The signs of the α-effect as well as rotational shear in the dynamo active layers can be chosen to give dynamo waves that propagate in opposite directions (poleward and equatorward) if the layers are considered separately in the framework of the Parker migratory dynamo. In a sequence of numerical experiments we show that the variety of dynamo-generated magnetic configurations in the system under discussion is quite rich. We identify the possibility of almost independent dynamo waves existing in the two layers as well as enslavement of one layer by the other, and of activity waves generated by a joint action of the two layers. We suggest some qualitative explanations of the behaviour and discuss also the limited nature of these explanations. This variety of phenomena suggests previously underexploited freedoms to understand how predictions of dynamo theory may accommodate the observed solar and stellar activity phenomenology.     

On the possibility of a bimodal solar dynamo

A simple way to couple an interface dynamo model to a fast tachocline model is presented, under the assumption that the dynamo saturation is due to a quadratic process and that the effect of finite shear layer thickness on the dynamo wave frequency is analogous to the effect of finite water depth on surface gravity waves. The model contains one free parameter which is fixed by the requirement that a solution should reproduce the helioseismically determined thickness of the tachocline. In this case it is found that, in addition to this solution, another steady solution exists, characterized by a four times thicker tachocline and 4–5 times weaker magnetic fields. It is tempting to relate the existence of this second solution to the occurrence of grand minima in solar activity. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Calibration of Vector Magnetogram with the Nonlinear Least-squares Fitting Technique

To acquire Stokes profiles from observations of a simple sunspot with the Video Vector Magnetograph at Huairou Solar Observing Station (HSOS), we scanned the FeI λ5324.19A line over the wavelength interval from 150mA redward of the line center to 150 mA blueward, in steps of 10 mA. With the technique of analytic inversion of Stokes profiles via nonlinear least-squares, we present the calibration coefficients for the HSOS vector magnetic magnetogram. We obtained the theoretical calibration error with linear expressions derived from the Unno-Becker equation under weak-field approximation.     

Cluster analysis for pattern recognition in solar butterfly diagrams

We investigate to what extent the wings of solar butterfly diagrams can be separated without an explicit usage of Hale's polarity law as well as the location of the solar equator. We apply two algorithms of cluster analysis for this purpose, namely DBSCAN and C‐means, and demonstrate their ability to separate the wings of contemporary butterfly diagrams based on the sunspot group density in the diagram only. Then we apply the method to historical data concerning the solar activity in the 18th century (Staudacher data). The method separates the two wings for Cycle 2, but fails to separate them for Cycle 1. In our opinion, this finding supports the interpretation of the Staudacher data as an indication of the unusual nature of the solar cycle in the 18th century (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Search for relationship between duration of the extended solar cycles and amplitude of sunspot cycle

Duration of the extended solar cycles is taken into the consideration. The beginning of cycles is counted from the moment of polarity reversal of large-scale magnetic field in high latitudes, occurring in the sunspot cycle n till the minimum of the cycle n + 2. The connection between cycle duration and its amplitude is established. Duration of the “latent” period of evolution of extended cycle between reversals and a minimum of the current sunspot cycle is entered. It is shown, that the latent period of cycles evolution is connected with the next sunspot cycle amplitude and can be used for the prognosis of a level and time of a sunspot maximum. The 24th activity cycle prognosis is made. The found dependences correspond to transport dynamo model of generation of solar cyclicity, it is possible with various speed of meridional circulation. Long-term behavior of extended cycle's lengths and connection with change of a climate of the Earth is considered. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Looking inside the butterfly diagram

The suitability of Maunder's butterfly diagram to give a realistic picture of the photospheric magnetic flux large scale distribution is discussed. The evolution of the sunspot zone in cycle 20 through 23 is described. To reduce the noise which covers any structure in the diagram, a smoothing algorithm has been applied to the sunspot data. This operation has eliminated any short period fluctuation, and given visibility to long duration phenomena. One of these phenomena is the fact that the equatorward drift of the spot zone center of mass results from the alternation of several prograde (namely, equatorward) segments with other stationary or poleward segments. The long duration of the stationary/retrograde phases as well as the similarities among the spot zone alternating paths in the cycles under examination prevent us from considering these features as meaningless fluctuations, randomly superimposed on the continuous equatorward migration. On the contrary, these features should be considered physically meaningful phenomena, requiring adequate explanations. Moreover, even the smoothed spotted area markedly oscillates. The compared examination of area and spot zone evolution allows us to infer details about the spotted area distribution inside the butterfly diagram. Links between the changing structure of the spot zone and the tachocline rotation rate oscillations are proposed. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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

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

Solar dynamo models with α ‐effect and turbulent pumping from local 3D convection calculations

Results from kinematic solar dynamo models employing α ‐effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2–3 m s^–1, having scaled the local quantities with the convective velocity at the bottom of the convection zone from a solar mixing‐length model. Rotation profile of the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface shear layer on the dynamo solutions. With these choices of the small‐ and large‐scale velocity fields, we obtain estimate of the ratio of the two induction effects, C _α /C _Ω ≈ 10^–3, which we keep fixed in all models. We also include a one‐cell meridional circulation pattern having a magnitude of 10–20 m s^–1 near the surface and 1–2 m s^–1 at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the α ‐effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of latitude. Turbulent pumping of the mean fields is predominantly down‐ and equatorward. The anisotropies in the turbulent diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B _r B _ϕ . Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the α ‐effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5° … 60°). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar‐type solutions. The role of the surface shear layer is found to be important only in the case where the α ‐effect has an appreciable magnitude near the surface. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)