A new look at the relationship between activity, dynamo number and Rossby number in late-type stars

The correlation between stellar activity, as measured by the indicator Δ R _HK, and the Rossby number Ro in late-type stars is revisited in light of recent developments in solar dynamo theory. Different stellar interior models, based on both mixing-length theory and the full spectrum of turbulence, are used in order to see to what extent the correlation of activity with Rossby number is model dependent, or otherwise can be considered universal. Although we find some modest model dependence, we find that the correlation of activity with Rossby number is significantly better than with rotation period alone for all the models we consider. Dynamo theory suggests that activity should scale with the dynamo number. A current model of the solar dynamo, the so-called interface dynamo, proposes that the amplification of the toroidal magnetic field by differential rotation (the ω -effect) and the production of the poloidal magnetic field from toroidal by helical turbulence (the α -effect) take place in different, adjacent layers near the base of the convection zone. A new scale analysis based on the interface dynamo shows that the appropriate dynamo number does not depend on the Rossby number alone, but also depends on an additional dimensionless factor related to the differential rotation. This leads to a new interpretation of the correlation between activity and Rossby number, which in turn leads to some conclusions about the magnitude of differential rotation in the dynamo layers of late-type main-sequence stars.     

The slender solar tachocline: a magnetic model

A model for the sharp transition from differential rotation in the solar convection zone to rigid rotation in the radiative interior is presented. Differential rotation in the radiative zone is shown to be quenched efficiently by an internal magnetic field. The poloidal field amplitude, B₀, is the input parameter for our model which determines the transition layer thickness and the toroidal field strength. It is illustrated analytically and confirmed numerically that a rather small field, B₀ = 10⁻⁴ Gauss, suffices to satisfy the helioseismological restrictions on the depth of the differential rotation penetration below the convection zone. The transition layer thickness decreases further with increasing B₀. The toroidal field amplitude B ≃ 200 Gauss is almost independent of B₀.     

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)     

A joined model for solar dynamo and differential rotation

A model for the solar dynamo, consistent in global flow and numerical method employed with the differential rotation model, is developed. The magnetic turbulent diffusivity is expressed in terms of the entropy gradient, which is controlled by the model equations. The magnetic Prandtl number and latitudinal profile of the alpha-effect are specified by fitting the computed period of the activity cycle and the equatorial symmetry of magnetic fields to observations. Then, the instants of polar field reversals and time-latitude diagrams of the fields also come into agreement with observations. The poloidal field has a maximum amplitude of about 10 Gs in the polar regions. The toroidal field of several thousand Gauss concentrates near the base of the convection zone and is transported towards the equator by the meridional flow. The model predicts a value of about 10³⁷ erg for the total magnetic energy of large-scale fields in the solar convection zone.     

Reflection and Ducting of Gravity Waves Inside the Sun

Internal gravity waves excited by overshoot at the bottom of the convection zone can be influenced by rotation and by the strong toroidal magnetic field that is likely to be present in the solar tachocline. Using a simple Cartesian model, we show how waves with a vertical component of propagation can be reflected when traveling through a layer containing a horizontal magnetic field with a strength that varies with depth. This interaction can prevent a portion of the downward traveling wave energy flux from reaching the deep solar interior. If a highly reflecting magnetized layer is located some distance below the convection zone base, a duct or wave guide can be set up, wherein vertical propagation is restricted by successive reflections at the upper and lower boundaries. The presence of both upward and downward traveling disturbances inside the duct leads to the existence of a set of horizontally propagating modes that have significantly enhanced amplitudes. We point out that the helical structure of these waves makes them capable of generating an α-effect, and briefly consider the possibility that propagation in a shear of sufficient strength could lead to instability, the result of wave growth due to over-reflection.     

Hydrodynamic Source of the 11-year Solar Variations

In the traditional axisymmetric models of the 11-year solar cycle, oscillations of the magnetic fields appear in the background of nonoscillating (over time scale considered) turbulent velocity fields and differential rotation. In this paper, an alternative approach is developed: The excitation of magnetic oscillations with the 22-year period is the consequence of hydrodynamic oscillations with the 11-year period. In the excitation of hydrodynamic oscillations, two processes taking place in high latitudes near the interface between the convective and radiative zones play a key role. One is forcing of the westerly zonal flow, the conditions for which are due to deformation of the interfacial surface. The other process is the excitation of a shear instability of zonal flow as a consequence of a strong radial gradient of angular velocity. The development of a shear instability at some stage brings about the disruption of the forcing of differential rotation. In the first (hydrodynamic) part of the paper, the dynamics of axisymmetric flows near the bottom of the convection zone is numerically simulated. Forcing of differential rotation having velocity shear in latitude and the existence of solutions in the form of torsional waves with the 11-year oscillation period are shown. In the second part the dynamics of the magnetic field is studied. The most pronounced peculiarities of the solutions are the existence of forced oscillations with the 22-year period and the drift of the toroidal magnetic field component from the mid latitudes to the equator. In high and low latitudes after cycle maximum, the toroidal component is of opposite sign in accordance with observations. In the third part, the transport of momentum from the bottom of the convection zone to the outer surface by virtue of diffusivity is considered. The existence of some sources of differential rotation in the convection zone is not implied. A qualitative correspondence of the differential rotation profile in the bulk of the convection zone and on its outer surface to experimental data is shown. The time correspondence between torsional and magnetic oscillations is also in accordance with observations.     

The Alpha Effect and the Observed Twist and Current Helicity of Solar Magnetic Fields

We present a straightforward comparison of model calculations for the α-effect, helicities, and magnetic field line twist in the solar convection zone with magnetic field observations at atmospheric levels. The model calculations are carried out in a mixing-length approximation for the turbulence with a profile of the solar internal rotation rate obtained from helioseismic inversions. The magnetic field data consist of photospheric vector magnetograms of 422 active regions for which spatially-averaged values of the force-free twist parameter and of the current helicity density are calculated, which are then used to determine latitudinal profiles of these quantities. The comparison of the model calculations with the observations suggests that the observed twist and helicity are generated in the bulk of the convection zone, rather than in a layer close to the bottom. This supports two-layer dynamo models where the large-scale toroidal field is generated by differential rotation in a thin layer at the bottom while the α-effect is operating in the bulk of the convection zone. Our previous observational finding was that the moduli of the twist factor and of the current helicity density increase rather steeply from zero at the equator towards higher latitudes and attain a certain saturation at about 12 – 15^∘. In our dynamo model with algebraic nonlinearity, the increase continues, however, to higher latitudes and is more gradual. This could be due to the neglect of the coupling between small-scale and large-scale current and magnetic helicities and of the latitudinal drift of the activity belts in the model.     

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)     

Nonlinear Evolution of 2d Tachocline Instabilities

A spectral method is used to explore the nonlinear evolution of known linear instabilities in a 2D differentially rotating magneto-hydrodynamic shell, representing the solar tachocline. Several simulations are presented, with a range of outcomes for the magnetic field configuration. Most spectacularly, the `clam instability', which occurs for solar differential rotation and a strong broad toroidal magnetic field structure, results in the field tipping over by 90° and reconnecting. A common characteristic of all the simulations though is that the nonlinear instabilities produce a strong angular momentum mixing effect which pushes the rotation towards a solid body form. It is argued that this may be the mechanism required by the model of Spiegel and Zahn to limit the tachocline's thickness.     

Large-scale flows excited by magnetic fields in the solar convective zone

Observations demonstrate a nearly 22-year periodic zonal flow superimposed on general solar differential rotation (LaBonte and Howard, 1982) and some meridional motions (e.g., Tuominen, Tuominen, and Kyrolänen, 1983). Such flows can be excited by the magnetic wave generated by the dynamo in the solar convective zone.An approximate analytical solution for the zonal and meridional flows for a given magnetic wave is constructed. This approach is justified by the fact that the magnetic field is generated by differential rotation and mean helicity, and the magnetic field in the time interval under consideration does not affect much this main flow; it can, however, strongly influence the perturbations of this flow.The density gradient in the convective zone is taken into account as an essential point in the solution construction. The solution agreed well with observational features and, in particular, it gives a phase shift between the rotational (zonal) wave and solar activity. A polar branch of the rotational wave can be described as an effect created by a poleward moving dynamo wave.Secular variations in the symmetrical part of the differential rotation and in the asymmetry between the north and south hemispheres are predicted.The alternative approaches to the explanation of the origin of the observed large-scale flows are discussed.     

Some recent developments in solar dynamo theory

We discuss the current status of solar dynamo theory and describe the dynamo model developed by our group. The toroidal magnetic field is generated in the tachocline by the strong differential rotation and rises to the solar surface due to magnetic buoyancy to create active regions. The decay of these active regions at the surface gives rise to the poloidal magnetic field by the Babcock-Leighton mechanism. This poloidal field is advected by the meridional circulation first to high latitudes and then down below to the tachocline. Dynamo models based on these ideas match different aspects of observational data reasonably well.     

On the solar rotation and activity

The interaction between differential rotation and magnetic fields in the solar convection zone was recently modelled by Brun (2004). One consequence of that model is that the Maxwell stresses can oppose the Reynolds stresses, and thus contribute to the transport of the angular momentum towards the solar poles, leading to a reduced differential rotation. So, when magnetic fields are weaker, a more pronounced differential rotation can be expected, yielding a higher rotation velocity at low latitudes taken on the average. This hypothesis is consistent with the behaviour of the solar rotation during the Maunder minimum. In this work we search for similar signatures of the relationship between the solar activity and rotation determined tracing sunspot groups and coronal bright points. We use the extended Greenwich data set (1878–1981) and a series of full-disc solar images taken at 28.4 nm with the EIT instrument on the SOHO spacecraft (1998–2000). We investigate the dependence of the solar rotation on the solar activity (described by the relative sunspot number) and the interplanetary magnetic field (calculated from the interdiurnal variability index). Possible rotational signatures of two weak solar activity cycles at the beginning of the 20th century (Gleissberg minimum) are discussed. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The rotation of magnetic loop systems in the solar atmosphere

Eclipse photographs indicate that large regions of the inner solar corona are confined in various types of closed magnetic configurations and, as a result, do not participate in the general solar wind expansion. In this paper, the rotation of initially poloidal loop configurations of this type, as influenced by differential rotation of the footpoints, is investigated. The analysis is restricted to axially symmetric fields and it is assumed that the toroidal magnetic field induced by differential rotation is small as compared to the initial poloidal field. This restricts the validity of the analysis to times less than about one month.The most interesting physical situation is that of flux tubes existing in one solar hemisphere only, one end of the tube being fixed in the photosphere at a higher latitude than the other. As a consequence, the lower end of the tube rotates at a faster rate than the upper end. Solution of the pertinent equations reveals that the angular velocity measured along a field line increases monotonically from its value at the poleward footpoint to that at the lower footpoint. The variation of angular velocity along the field depends upon the field geometry only and is not directly related to the variation of angular velocity along the solar surface between the footpoints. Depending upon the field configuration, both outward radial increases and decreases are possible. Using the Newton and Nunn model for the surface differential rotation rate, the angular velocity distribution on two particularly simple types of closed magnetic loop systems is determined analytically. It is shown that the angular velocity increases outward in the polar regions but decreases outward near the equator - leading to a decrease in differential rotation with height.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Complexes of activity of the solar cycle and very large scale convection

The magnetic field pattern associated with large scale convective motions, which are much larger than the supergranules and have been conceived as a source of maintenance of the solar differential rotation, is calculated in the framework of a slowly and differentially rotating thin spherical shell, including the effects of thermal conductivity and viscosity. The approximations of Boussinesq are used and the initial state of the magnetic field is assumed to be purely toroidal.The resulting magnetic field pattern rotates rigidly on the differentially rotating Sun with some phase delay to the convective pattern, if it is assumed that only the predominant mode with the maximum growth rate is actually realized in the solar convection zone. The obtained magnetic and convective patterns and their properties seem to explain naturally the various aspects of large scale ordering of solar activity such as the existence and behavior of complexes of activity, the rigid body rotation of proton flare active longitudes, their association with UMR's, the existence of ghost and mirror image of UMR's themselves and the fact that the rotational period derived from sunspot data is shorter than that derived spectroscopically from fluid velocity.     

On the interaction between differential rotation and magnetic fields in the Sun

We have performed 3-D numerical simulations of compressible convection under the influence of rotation and magnetic fields in spherical shells. They aim at understanding the subtle coupling between convection, rotation and magnetic fields in the solar convection zone. We show that as the magnetic Reynolds number is increased in the simulations, the magnetic energy saturates via nonlinear dynamo action, to a value smaller but comparable to the kinetic energy contained in the shell, leading to increasingly strong Maxwell stresses that tend to weaken the differential rotation driven by the convection. These simulations also indicate that the mean toroidal and poloidal magnetic fields are small compared to their fluctuating counterparts, most of the magnetic energy being contained in the non-axisymmetric fields. The intermittent nature of the magnetic fields generated by such a turbulent convective dynamo confirms that in the Sun the large-scale ordered dynamo responsible for the 22-year cycle of activity can hardly be located in the solar convective envelope.     

The solar dynamo

The basic features of the solar activity mechanism are explained in terms of the dynamo theory of mean magnetic fields. The field generation sources are the differential rotation and the mean helicity of turbulent motions in the convective zone. A nonlinear effect of the magnetic field upon the mean helicity results in stabilizing the amplitude of the 22-year oscillations and forming a basic limiting cycle. When two magnetic modes (with dipole and quadrupole symmetry) are excited nonlinear beats appear, which may be related to the secular cycle modulation.The torsional waves observed may be explained as a result of the magnetic field effect upon rotation. The magnetic field evokes also meriodional flows.Adctual variations of the solar activity are nonperiodic since there are recurrent random periods of low activity of the Maunder minimum type. A regime of such a magnetic hydrodynamic chaos may be revealed even in rather simple nonlinear solar dynamo models.The solar dynamo gives rise also to three-dimensional, non-axisymmetric magnetic fields which may be related to a sector structure of the solar field.     

Latitudinal migration of sunspots based on the ESAI database

The latitudinal migration of sunspots toward the equator,which implies there is propagation of the toroidal magnetic flux wave at the base of the solar convection zone,is one of the crucial observational bases for the solar dynamo to generate a magnetic field by shearing of the pre-existing poloidal magnetic field through differential rotation.The Extended time series of Solar Activity Indices(ESAI)elongated the Greenwich observation record of sunspots by several decades in the past.In this study,ESAI's yearly mean latitude of sunspots in the northern and southern hemispheres during the years 1854 to 1985 is utilized to statistically test whether hemispherical latitudinal migration of sunspots in a solar cycle is linear or nonlinear.It is found that a quadratic function is statistically significantly better at describing hemispherical latitudinal migration of sunspots in a solar cycle than a linear function.In addition,the latitude migration velocity of sunspots in a solar cycle decreases as the cycle progresses,providing a particular constraint for solar dynamo models.Indeed,the butterfly wing pattern with a faster latitudinal migration rate should present stronger solar activity with a shorter cycle period,and it is located at higher latitudinal position,giving evidence to support the Babcock-Leighton dynamo mechanism.     

Double maxima of 11-year solar cycles

We propose a scenario to explain the observed phenomenon of double maxima of sunspot cycles, including the generation of a magnetic field near the bottom of the solar convection zone (SCZ) and the subsequent rise of the field from the deep layers to the surface in the royal zone. Five processes are involved in the restructuring of the magnetic field: the Ω-effect, magnetic buoyancy, macroscopic turbulent diamagnetism, rotary ?ρ-effect, and meridional circulation. It is found that the restructuring of magnetism develops differently in high-latitude and equatorial domains of the SCZ. A key role in the proposed mechanism of the double maxima is played by two waves of toroidal fields from the lower base of the SCZ to the solar surface in the equatorial domain. The deep toroidal fields are excited by the Ω-effect near the tachocline at the beginning of the cycle. Then these fields are transported to the surface due to the combined effect of magnetic buoyancy, macroscopic turbulent diamagnetism, and the rotary magnetic ?ρ-flux in the equatorial domain. After a while, these magnetic fragments can be observed as bipolar sunspot groups at the middle latitudes in the royal zone. This first, upward-directed wave of toroidal fields produces the main maximum of sunspot activity. However, the underlying toroidal fields in the high-latitude polar domains are blocked at the beginning of the cycle near the SCZ bottom by two antibuoyancy effects — the downward turbulent diamagnetic transfer and the magnetic ?ρ-pumping. In approximately 1 or 2 years, a deep equatorward meridional flow transfers these fields to low-latitude parts of the equatorial domain (where there are favorable conditions for magnetic buoyancy), and the belated magnetic fields (the second wave of toroidal fields) rise to the surface. When this second batch of toroidal fields comes to the solar surface at low latitudes, it leads to the second sunspot maximum.     

Turbulent dynamo near tachocline and reconstruction of azimuthal magnetic field in the solar convection zone

In order to extend the abilities of the αΩ dynamo model to explain the observed regularities and anomalies of the solar magnetic activity, the negative buoyancy phenomenon and the magnetic quenching of the α effect were included in the model, as well as newest helioseismically determined inner rotation of the Sun were used. Magnetic buoyancy constrains the magnitude of toroidal field produced by the Ω effect near the bottom of the solar convection zone (SCZ). Therefore, we examined two “antibuoyancy” effects: i) macroscopic turbulent diamagnetism and ii) magnetic advection caused by vertical inhomogeneity of fluid density in the SCZ, which we call the ∇ρ effect. The Sun's rotation substantially modifies the ∇ρ effect. The reconstruction of the toroidal field was examined assuming the balance between mean‐field magnetic buoyancy, turbulent diamagnetism and the rotationally modified ∇ρ effect. It is shown that at high latitudes antibuoyancy effects block the magnetic fields in the deep layers of the SCZ, and so the most likely these deep‐rooted fields could not become apparent at the surface as sunspots. In the near‐equatorial region, however, the upward ∇ρ effect can facilitate magnetic fields of about 3000 – 4000 G to emerge through the surface at the sunspot belt. Allowance for the radial inhomogeneity of turbulent velocity in derivations of the helicity parameter resulted in a change of sign of the α effect from positive to negative in the northern hemisphere near the bottom of the SCZ. The change of sign is very important for direction of the Parker's dynamo‐waves propagation and for parity of excited magnetic fields. The period of the dynamo‐wave calculated with allowance for the magnetic quenching is about seven years, that agrees by order of magnitude with the observed mean duration of the sunspot cycles. Using the modern helioseismology data to define dynamo‐parameters, we conclude that north‐south asymmetry should exist in the meridional field. At low latitudes in deep layers of the SCZ, the αΩ dynamo excites most efficiency the dipolar mode of the meridional field. Meanwhile, in high‐latitude regions a quadrupolar mode dominates in the meridional field. The obtained configuration of the net meridional field is likely to explain the magnetic anomaly of polar fields (the apparent magnetic “monopole”) observed near the maxima of solar cycles. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic Helicity as a Predictor of the Solar Cycle

It is well known that the polar magnetic field is at its maximum during solar minima, and that the behaviour during this time acts as a strong predictor of the strength of the following solar cycle. This relationship relies on the action of differential rotation (the Omega effect) on the poloidal field, which generates the toroidal flux observed in sunspots and active regions. We measure the helicity flux into both the northern and the southern hemispheres using a model that takes account of the Omega effect, which we apply to data sets covering a total of 60 years. We find that the helicity flux offers a strong prediction of solar activity up to five years in advance of the next solar cycle. We also hazard an early guess as to the strength of Solar Cycle 25, which we believe will be of similar amplitude and strength to Cycle 24.