The Sun’s Large-Scale Magnetic Field and Its Long-Term Evolution

The Sun’s large-scale external field is formed through the emergence of magnetic flux in active regions and its subsequent dispersal over the solar surface by differential rotation, supergranular convection, and meridional flow. The observed evolution of the polar fields and open flux (or interplanetary field) during recent solar cycles can be reproduced by assuming a supergranular diffusion rate of 500 – 600 km² s⁻¹ and a poleward flow speed of 10 –20 m s⁻¹. The nonaxisymmetric component of the large-scale field decays on the flow timescale of ∼1 yr and must be continually regenerated by new sunspot activity. Stochastic fluctuations in the longitudinal distribution of active regions can produce large peaks in the Sun’s equatorial dipole moment and in the interplanetary field strength during the declining phase of the cycle; by the same token, they can lead to sudden weakenings of the large-scale field near sunspot maximum (Gnevyshev gaps). Flux transport simulations over many solar cycles suggest that the meridional flow speed is correlated with cycle amplitude, with the flow being slower during less active cycles.     

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.     

Systematically Asymmetric Heliospheric Magnetic Field: Evidence for a Quadrupole Mode and Non-Axisymmetry with Polarity Flip-Flops

Recent studies of the heliospheric magnetic field (HMF) have detected interesting, systematic hemispherical and longitudinal asymmetries which have a profound significance for the understanding of solar magnetic fields. The in situ HMF measurements since the 1960s show that the heliospheric current sheet (HCS) is systematically shifted (coned) southward during solar minimum times, leading to the concept of a bashful ballerina. While temporary shifts can be considerably larger, the average HCS shift (coning) angle is a few degrees, less than the 7.2^∘ tilt of the solar rotation axis. Recent solar observations during the last two solar cycles verify these results and show that the magnetic areas in the northern solar hemisphere are larger and their intensity weaker than in the south during long intervals in the late declining to minimum phase. The multipole expansion reveals a strong quadrupole term which is oppositely directed to the dipole term. These results imply that the Sun has a symmetric quadrupole S0 dynamo mode that oscillates in phase with the dominant dipole A0 mode. Moreover, the heliospheric magnetic field has a strong tendency to produce solar tilts that are roughly opposite in longitudinal phase. This implies is a systematic longitudinal asymmetry and leads to a “flip-flop” type behaviour in the dominant HMF sector whose period is about 3.2 years. This agrees very well with the similar flip-flop period found recently in sunspots, as well as with the observed ratio of three between the activity cycle period and the flip-flop period of sun-like stars. Accordingly, these results require that the solar dynamo includes three modes, A0, S0 and a non-axisymmetric mode. Obviously, these results have a great impact on solar modelling.     

The decay of the mean solar magnetic field

We have analyzed the effects that differential rotation and a hypothetical meridional flow would have on the evolution of the Sun's mean line-of-sight magnetic field as seen from Earth. By winding the large-scale field into strips of alternating positive and negative polarity, differential rotation causes the mean-field amplitude to decay and the mean-field rotation period to acquire the value corresponding to the latitude of the surviving unwound magnetic flux. For a latitudinally broad two-sector initial field such as a horizontal dipole, the decay is rapid for about 5 rotations and slow with a t ^–1/2 dependence thereafter. If a poleward meridional flow is present, it will accelerate the decay by carrying the residual flux to high latitudes where the line-of-sight components are small. The resulting decay is exponential with an e-folding time of 0.75 yr (10 rotations) for an assumed 15 m s^–1 peak meridional flow speed.E.O. Hulburt Center for Space Research.Laboratory for Computational Physics.     

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)     

Sector-structured interplanetary magnetic field associated with the fast plasma streams in 1985–1996

An analysis of 373 well-defined high-speed solar-wind streams observed at 1 AU during the years 1985–1996 is outlined. The distribution of the occurrence of these streams as a function of Bartels rotation days using the dominant polarity of the interplanetary magnetic field (IMF) associated with the referred fast streams shows that a four-sector pattern for the positive IMF polarity and a two-sector pattern for the negative IMF polarity are the dominant features in the investigated period. The high-speed streams seem to occur at preferred Bartels days: positive polarity streams are most frequent near Bartels days 5 and 18, while negative polarity streams are most frequent in days 14 and 23. Moreover, the corotating streams with positive IMF polarity prefer to occur in days 5 and 18 of the Bartels rotation period, whereas flare-generated streams with negative IMF polarity occur in days 14 and 23. The observed distribution of Bartels days is probably related to the distribution of the solar sources of high-speed solar wind streams as the solar wind carries with it the photospheric magnetic polarity of the solar source region. In addition, the distribution of the streams reveals a similar behaviour during the ascending and the declining phase of the last solar cycle (22nd) in contrast to the previous one where it has an opposite appearance. Determined differences in the characteristics of the sector structured IMF associated with the fast streams of the last cycle with the previous one (21st) and some similarities with the alternate solar cycle (20th) seem to be attributed to the 22-year magnetic cycle and to the polarity reversals of the polar magnetic field of the Sun. As the magnetic sectors are due to multiple crossings of the solar equatorial plane by a large-scale, warped heliospheric current sheet, it is suggested that the two-sector pattern arises from a tilted solar magnetic dipole component and the more commonly observed four-sector pattern from a quadrupole component of the solar interplanetary magnetic field.     

Solar-Cycle-Related Variations in the Solar Differential Rotation and Meridional Flow: A Comparison

We have analysed a large set of sunspot group data (1874 – 2004) and find that the meridional flow strongly varies with the phase of the solar cycle, and the variation is quite different in the northern and the southern hemispheres. We also find the existence of considerable cycle-to-cycle variation in the meridional velocity, and about a 11-year difference between the phases of the corresponding variations in the northern and the southern hemispheres. In addition, our analysis also indicates the following: (i) the existence of a considerable difference (about 180^°) between the phases of the solar-cycle variations in the latitude-gradient terms of the northern and the southern hemispheres’ rotations; (ii) the existence of correlation (good in the northern hemisphere and weak in the southern hemisphere) between the mean solar-cycle variations of meridional flow and the latitude-gradient term of solar rotation; (iii) in the northern hemisphere, the cycle-to-cycle variation of the mean meridional velocity leads that of the equatorial rotation rate by about 11 years, and the corresponding variations have approximately the same phase in the southern hemisphere; and (iv) the directions of the mean meridional velocity is largely toward the pole in the longer sunspot cycles and largely toward the equator in the shorter cycles.     

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.     

Large-scale motions in the sun

Large-scale solar motions comprise differential rotation (with latitudinal, and perhaps radial gradients), axially symmetric meridional motions, and possible asymmetric motions (giant convective cells or Rossby-type waves or both). These motions must be basic in any satisfactory theory of the changing pattern of solar magnetic fields and of the 22-yr cycle. In the present paper available data are discussed and, as far as possible, evaluated and explained.Rotational measurements are based on the changing positions of discrete features such as sunspots, on Doppler shifts, on geophysical changes and on statistical evaluation of the motions of diffuse objects. The first mentioned, comprising faculae, sunspots, K-corona (to latitudes 45°) and filaments, show agreement better than 0.7 %. A new formula for surface rotation _s , based on faculae and sunspot data, is _s = 14.52 – 2.48 sin² b – 2.51 sin⁶ b deg day^–1, where b is latitude, and validity may extend to about 70°. Errors in Doppler shift measurements and statistical treatments are discussed. There is evidence of a much slower coronal rate at high latitudes, and of a slower sub-surface rate at lower latitudes.Ordered meridional motions have been revealed by statistical investigations of the positions of spot groups, of spots and of filaments. All these results seem explicable in terms of an oscillating hydro-magnetic circulation in each hemisphere. These have both 11-yr and 22-yr components, and these periods are provided by a general dipole field of about one gauss, together with a pair of toroidal fields centred at latitudes ±16° and of average strength of order 10 G.Evidence of large-scale (perhaps 3 × 10⁵ km), irregular surface motions is provided by the distribution of surface magnetic flux, the motions of sunspots, and Doppler-shift observations; it is supported by Ward's theory of the equatorial acceleration. The possibility is suggested that these asymmetric motions also drive the oscillatory meridional motions.     

Coronal general magnetic field evolution as a new parameter of the solar cycle

The magnetic field lines of the corona associated with the solar-cycle surface general magnetic field are calculated by a potential-field approximation to study the solar-cycle evolution of the geometry of the coronal field. The surface field evolution used here is the radial field evolution, predicted by a model of the solar cycle driven by the dynamo action of the global convection, and justified observationally using Mount Wilson magnetic synoptic chart data. The evolution of the calculated coronal general field is now good for comparison with observations and shows the following. (i) The field of the polar and high-altitude corona has dipolar structure in almost all phases of the solar cycle except in a short time interval around maximum phase despite the quadrupolar structure of the general magnetic field at the surface; quadrupolar field forms loop-like structure in the lower corona. The almost-dipolar structure of the polar and high-latitude corona and the loop formation of the equatorial lower corona explain the appearance of the undisturbed minimum corona observed at eclipses. (ii) The polar field lines are directed almost radially at the minimum phase, which should be responsible for polar plumes. The field lines slowly open up to participate in the loop-like structure of the equatorial lower corona, and rapidly change their structure and polarity at the maximum phase, to resume the almost radial configuration slowly, (iii) During the rapidly changing maximum phase, the field lines do not penetrate deep into the interplanetary space resulting in the absence of polar plumes and the appearance of the circular corona- the maximum corona. In this phase, the coronal field should not be approximated by a dipole field. The surface field evolution which can explain such behaviors of the corona is characteristic of the solar-cycle process dominated by the latitudinal gradient of the differential rotation. If the radial gradient dominated in the subsurface process, the coronal evolution would look quite different and would show latitudinal propagation of enhancement of activity. Although nonaxisymmetric features should be superposed on the axisymmetric general field to express the real corona, the general field can be a basic coronal field in studying long-term interaction between the convection zone and the interstellar space especially in studying the magnetic braking of the solar rotation.     

An Exploration of Non-kinematic Effects in Flux Transport Dynamos

Recent global magnetohydrodynamical simulations of solar convection producing a large-scale magnetic field undergoing regular, solar-like polarity reversals also present related cyclic modulations of large-scale flows developing in the convecting layers. Examination of these simulations reveal that the meridional flow, a crucial element in flux transport dynamos, is driven at least in part by the Lorentz force associated with the cycling large-scale magnetic field. This suggests that the backreaction of the field onto the flow may have a pronounced influence on the long-term evolution of the dynamo. We explore some of the associated dynamics using a low-order dynamo model that includes this Lorentz force feedback. We identify several characteristic solutions which include single period cycles, period doubling and chaos. To emulate the role of turbulence in the backreaction process we subject the model to stochastic fluctuations in the parameter that controls the Lorentz force amplitude. We find that short term fluctuations produce long-term modulations of the solar cycle and, in some cases, grand minima episodes where the amplitude of the magnetic field decays to near zero. The chain of events that triggers these quiescent phases is identified. A subsequent analysis of the energy transfer between large-scale fields and flows in the global magnetohydrodynamical simulation of solar convection shows that the magnetic field extracts energy from the solar differential rotation and deposits part of that energy into the meridional flow. The potential consequences of this marked departure from the kinematic regime are discussed in the context of current solar cycle modeling efforts based on flux transport dynamos.     

The Evolution of the Sun's Open Magnetic Flux – II. Full Solar Cycle Simulations

In this paper the origin and evolution of the Sun's open magnetic flux is considered by conducting magnetic flux transport simulations over many solar cycles. The simulations include the effects of differential rotation, meridional flow and supergranular diffusion on the radial magnetic field at the surface of the Sun as new magnetic bipoles emerge and are transported poleward. In each cycle the emergence of roughly 2100 bipoles is considered. The net open flux produced by the surface distribution is calculated by constructing potential coronal fields with a source surface from the surface distribution at regular intervals. In the simulations the net open magnetic flux closely follows the total dipole component at the source surface and evolves independently from the surface flux. The behaviour of the open flux is highly dependent on meridional flow and many observed features are reproduced by the model. However, when meridional flow is present at observed values the maximum value of the open flux occurs at cycle minimum when the polar caps it helps produce are the strongest. This is inconsistent with observations by Lockwood, Stamper and Wild (1999) and Wang, Sheeley, and Lean (2000) who find the open flux peaking 1–2 years after cycle maximum. Only in unrealistic simulations where meridional flow is much smaller than diffusion does a maximum in open flux consistent with observations occur. It is therefore deduced that there is no realistic parameter range of the flux transport variables that can produce the correct magnitude variation in open flux under the present approximations. As a result the present standard model does not contain the correct physics to describe the evolution of the Sun's open magnetic flux over an entire solar cycle. Future possible improvements in modeling are suggested.     

An equatorially symmetric rotator model for magnetic stars

Suggested by theoretical and empirical investigations, a very simple magnetic field distribution at the surface of magnetic stars is proposed. This model corresponds to the superposition of a dipole and a quadrupole, both lying in the equatorial plane with their axes in the same direction. The observed variations of the effective magnetic field strength of the sixteen stars investigated can be represented surprisingly well with this model. When the quadrupole contribution is large, the variation of the magnetic field vector over the surface of the star differs characteristically from that of the decentred dipole model; thus the difficulty with the latter, regarding the representation of both the effective magnetic field and the mean surface field with the same parameters may perhaps be avoided.     

Nonlinear boussinesq convective model for large scale solar circulations

We present extensive numerical calculations for a model of thermal convection of a Boussinesq fluid in an equatorial annulus of a rotating spherical shell. The convection induces and maintains differential rotation and meridian circulation. The model is solved for an effective Prandtl number P = 1, with effective Taylor number T in the range 10² <T <10⁶, and effective Rayleigh number R between the critical value for onset of convection, and a few times that value. With = 2.6 × 10^–6 s^–1, d = 1.4 × 10¹⁰ cm (roughly the depth of the solar convection zone) the range of Taylor number is equivalent to kinematic viscosities between 10¹⁴ and 10¹² cm² s^–1, which encompasses eddy viscosities estimated from mixing length theory applied to the Sun.The convection does generally make equatorial regions rotate faster, the more so as T is increased, but local equatorial deceleration near the surface is also produced at intermediate T for large enough R above critical. The differential rotation is maintained primarily through momentum transport in the cells up the gradient, rather than by meridian circulation. Differential rotation energy increases relative to cell energy with increasing T, surpassing it near T = 3 × 10⁴. The differential rotation tends to stretch out the convective cells, analogously to what is thought to happen to solar magnetic regions. Differential rotation and meridian circulation energies are nearly equal for T = 10³, but the meridian circulation energy falls off relative to differential rotation like T ^–1 for larger T. The meridian circulation is always toward the poles near the surface, contrary to models of Kippenhahn, Cocke, Köhler, and Durney and Roxburgh. The radial shear produced in the differential rotation is almost always positive, as in the Köhler model, but contrary to the assumptions made by Leighton for his random walk solar cycle model.Solutions in the neighborhood of T = 3 × 10⁴ seem to compare best with various solar observations including differential rotation amplitude, cell wavelength, tilted structure, horizontal momentum transport, and weak meridian circulation. The local equatorial deceleration (equatorward of 10–15°) has not been observed, although the techniques of data analysis may not have been sensitive to it. The most important deficiency of the model is that all the solutions with T 10³ show the vertical heat transport a rather strong function of latitude, with a maximum at the equator, no evidence of which is seen at the solar surface.The National Center for Atmospheric Research is sponsored by The National Science Foundation.     

Differential rotation and sector structure of solar magnetic fields

The differential rotation and sector structure of solar magnetic fields has been studied using digitized data on photospheric magnetic fields recorded at the Mount Wilson Observatory during the period August 1959–May 1970. The power spectra show considerable power in high-frequency peaks, corresponding to harmonic components with wavelengths less than 1/10 solar rotation. Calculations for a series of shorter time intervals show how the distribution of power over the various harmonic components in the sector pattern varies strongly with the solar cycle. The equatorial rotation rate of solar magnetic fields is about 0.1 km s^-1 faster than that of the photospheric plasma determined from Doppler shifts. It is shown that the Doppler measurements mainly refer to the non-network regions. The differential flow of 0.1 km s^-1 forms streamlines around the magnetic fine structures. The different rotation rates of various solar features can be explained in terms of the rotation rates of magnetic and non-magnetic regions. The rotation rates of the magnetic fields in active and quiet regions agree at the equator. At higher latitudes, however, the background fields deviate less from solid-body rotation, indicating that their source is below the deepest layers to which the sunspot magnetic fields penetrate. This suggests that turbulent diffusion of the field in old active regions may not be the main source for the background magnetic field, but that the source is located close to a rigidly rotating solar core with a synodic rotation period of 26.87 days.     

Self-excitation of motions in near-convection zones and the generation of solar magnetic field

In view of the recently discovered time variations in rotation velocity within the solar differentially rotating tachocline (Howe et al. 2000), we study conditions for the equilibrium and excitation of motions in nonrigidly rotating magnetized layers of the radiative zones located near the boundaries of the convection zone. The emphasis is on the possible relationship between quasi-periodic tachocline pulsations and the generation of a nonaxisymmetric magnetic field in the convection zone. This field generation is studied under the assumption that it results from a reduction in the expenditure of energy on convective heat transport. The (antisymmetric about the equator) field is shown to increase in strength if there are both a radial gradient in angular velocity and steady-state axisymmetric meridional circulation of matter. The sense of circulation is assumed to change (causing the sign of the generated field to change) after the maximum permissible field strength is reached. This is apparently attributable to the excitation of the corresponding turbulent viscosity of the medium. It is also important that the cyclic field variations under discussion are accompanied by variations in solar-type dipole magnetic field.     

Magnetic helicity in stellar dynamos: new numerical experiments

The theory of large scale dynamos is reviewed with particular emphasis on the magnetic helicity constraint in the presence of closed and open boundaries. In the presence of closed or periodic boundaries, helical dynamos respond to the helicity constraint by developing small scale separation in the kinematic regime, and by showing long time scales in the nonlinear regime where the scale separation has grown to the maximum possible value. A resistively limited evolution towards saturation is also found at intermediate scales before the largest scale of the system is reached. Larger aspect ratios can give rise to different structures of the mean field which are obtained at early times, but the final saturation field strength is still decreasing with decreasing resistivity. In the presence of shear, cyclic magnetic fields are found whose period is increasing with decreasing resistivity, but the saturation energy of the mean field is in strong super‐equipartition with the turbulent energy. It is shown that artificially induced losses of small scale field of opposite sign of magnetic helicity as the large scale field can, at least in principle, accelerate the production of large scale (poloidal) field. Based on mean field models with an outer potential field boundary condition in spherical geometry, we verify that the sign of the magnetic helicity flux from the large scale field agrees with the sign of α. For solar parameters, typical magnetic helicity fluxes lie around 10⁴⁷ Mx² per cycle.     

Meridional drift in the large-scale solar magnetic field pattern

A method of two-dimensional correlation functions has been applied to a sequence of synoptic maps of the large-scale magnetic field to obtain the meridional drift pattern of field structures. The meridional drift profile obtained is antisymmetric about the equator. The meridional drift is directed from the equator to the poles at latitudes below 45°. A maximum drift velocity of 11–13 m s^–1 is attained in the latitude range 30°. A picture of the space-time distribution of meridional drift is also obtained, which may be interpreted as resulting from the effect of azimuthal convective rolls (3 rolls per hemisphere) on the large-scale magnetic field. Rolls originate at high latitudes following the cycle maximum, and migrate equatorwards until the minimum of the next cycle. The picture in the equatorial region can correspond to convective rolls with lifetimes of about two years, or to the process of interaction of rolls from two hemispheres.     

The sun's magnetic sector structure

The synoptic appearance of solar magnetic sectors is studied using 454 sector boundaries observed at Earth during 1959–1973. The sectors are clearly visible in the photospheric magnetic field. Sector boundaries can be clearly identified as north-south running demarcation lines between regions of persistent magnetic polarity imbalances. These regions extend up to about 35 ° of latitude on both sides of the equator. They generally do not extend into the polar caps. The polar cap boundary can be identified as an east-west demarcation line marking the poleward limit of the sectors. The typical flux imbalance for a magnetic sector is about 4 × 10²¹ Mx.     

Time-Distance Helioseismology with f Modes as a Method for Measurement of Near-Surface Flows

Travel times measured for the f mode have been used to study flows near the solar surface in conjunction with simultaneous measurements of the magnetic field. Previous flow measurements of Doppler surface rotation, small magnetic feature rotation, supergranular pattern rotation, and surface meridional circulation have been confirmed. In addition, the flow in supergranules due to Coriolis forces has been measured. The spatial and temporal power spectra for a six-day observing sequence have been measured.