Large-Scale Structure of the Sunspot Activity and the Background Magnetic Field Formation

Large-scale distribution of the sunspot activity of the Sun has been analyzed by using a technique worked out previously (Erofeev, 1997) to study long-lived, non-axisymmetric magnetic structures with different periods of rotation. Results of the analysis have been compared with those obtained by analyzing both the solar large-scale magnetic field and large-scale magnetic field simulated by means of the well-known flux transport equation using the sunspot groups as a sole source of new magnetic flux in the photosphere. A 21-year period (1964–1985) has been examined.The rotation spectra calculated for the total time interval of two 11-year cycles indicate that sunspot activity consists of a series of discrete components (modes) with different periods of rotation. The largest-scale component of the sunspot activity reveals modes with 27-day and 28-day periods of rotation situated, correspondingly, in the northern and southern hemispheres of the Sun, and two modes with rotation periods of about 29.7 days situated in both hemispheres. Such a modal structure of the sunspot activity agrees well with that of the large-scale solar magnetic field. Moreover, the magnetic field distribution simulated with the flux transport equation also reveals the same modal structure. However, such an agreement between the large-scale solar magnetic field and both the sunspot activity and simulated magnetic field is unstable in time; so, it is absent in the northern hemisphere of the Sun during solar cycle No. 20. Thus the sources of magnetic flux responsible for formation of the large-scale, rigidly rotating magnetic patterns appear to be closely connected, but are not identical with the discrete modes of the sunspot activity.     

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.     

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.     

The decay of the large-scale solar magnetic field

In the absence of new bipolar sources of flux, the large-scale magnetic field at the solar photosphere decays due to differential rotation, meridional flow, and supergranular diffusion. The rotational shear quickly winds up the nonaxisymmetric components of the field, increasing their latitudinal gradients and thus the rates of diffusive mixing of their flux. This process is particularly effective at mid latitudes, where the rotational shear is largest, so that eventually low- and high-latitude remnants of the initial, nonaxisymmetric field pattern survive. In this paper I solve analytically the transport equation describing the evolution of the large-scale photospheric field, to study its time-asymptotic behavior. The solutions are rigidly rotating, uniformly decaying distributions of flux, wound up by differential rotation and localized near either the equator or the poles. A balance between azimuthal transport of flux by the rotational shear and meridional transport by the diffusion gives rise to the rigidly rotating field patterns. The time-scale on which this balance is achieved, and also on which the nonaxisymmetric flux decays away, is the geometric mean of the short time-scale for shearing by differential rotation and the long time-scale for dispersal by supergranular diffusion. A poleward meridional flow alters this balance on its own, intermediate time-scale, accelerating the decay of the nonaxisymmetric flux at low latitudes. Such a flow also hastens the relaxation of the axisymmetric field to a modified dipolar configuration.     

Formation of polar starspots through meridional circulation

To explain the observed intermingling of polarities in the magnetic field distributions of rapidly rotating stars, surface magnetic flux transport models demand the presence of fast meridional flows.We combine simulations of the pre-eruptive and post-eruptive magnetic flux transport in cool stars to investigate the influence of a fast meridional circulation on the latitudinal eruption pattern of magnetic flux tubes and on the polar magnetic field properties. Magnetic flux tubes rising through the convection zone experience an enhanced latitude-dependent poleward deflection through meridional flows, which renders the wings of stellar butterfly diagrams convex. The larger amount of magnetic flux emerging at higher latitudes supports the intermingling of opposite polarities of polar magnetic fields and yields magnetic flux densities in the polar regions about 20% higher than in the case disregarding the pre-eruptive deflection. Taking the pre-eruptive evolution of magnetic flux into account therefore eases the need for the fast meridional flows predicted by previous investigations. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The concentration of the large-scale solar magnetic field by a meridional surface flow

We calculate analytical and numerical solutions to the magnetic flux transport equation in the absence of new bipolar sources of flux, for several meridional flow profiles and a range of peak flow speeds. We find that a poleward flow with a broad profile and a nominal 10 m s^–1 maximum speed concentrates the large-scale field into very small caps of less than 15° half-angle, with average field strengths of several tens of gauss, contrary to observations. A flow which reaches its peak speed at a relatively low latitude and then decreases rapidly to zero at higher latitudes leads to a large-scale field pattern which is consistent with observations. For such a flow, only lower latitude sunspot groups can contribute to interhemispheric flux annihilation and the resulting decay and reversal of the polar magnetic fields.     

The origin of the 28- to 29-day recurrent patterns of the solar magnetic field

Numerical simulations of the Sun's mean line-of-sight magnetic field suggest an origin for the 28-to 29-day recurrent patterns of the field and its associated interplanetary phenomena. The patterns are caused by longitudinal fluctuations in the eruption of new magnetic flux, the transport of this flux to mid latitudes by supergranular diffusion and meridional flow, and the slow rotation of the resulting flux distributions at the 28- to 29-day periods characteristic of those latitudes.E. O. Hulburt Center for Space Research.Laboratory for Computational Physics.     

Large-Scale Transport of Magnetic Flux on the sun

The structure of the large-scale background magnetic field evolves in time and space. The large-scale horizontal transport velocity field of the magnetic flux patterns was inferred over the whole solar photosphere in the course of two solar activity cycles from year 1976 to 1999. The method of velocity determination and the testing procedures of the velocity accuracy are presented. The non-axially symmetric component of the horizontal velocity was found and both zonal and meridional velocity regions were described. The horizontal large-scale transport velocity regions vary in shape and the intensity during different phases of the 11-year solar activity cycle. The total horizontal transport velocity is characterized by the presence of variable amounts of the vector field vortices with symmetric orientation relative to the solar equator. The zonal velocity regions, distributed inside of the zonal belt limited by latitudes ± 35°, are persistent for about 4 Carrington rotations. Recurrent structures of similar velocity distributions are not coherent over the whole solar photosphere.     

The impact of meridional circulation on stellar butterfly diagrams and polar caps

Observations of rapidly rotating solar-like stars show a significant mixture of opposite-polarity magnetic fields within their polar regions. To explain these observations, models describing the surface transport of magnetic flux demand the presence of fast meridional flows. Here, we link subsurface and surface magnetic flux transport simulations to investigate (i) the impact of meridional circulations with peak velocities of  ≤125 m s⁻¹  on the latitudinal eruption pattern of magnetic flux tubes and (ii) the influence of the resulting butterfly diagrams on polar magnetic field properties. Prior to their eruption, magnetic flux tubes with low field strengths and initial cross-sections below  ∼300 km  experience an enhanced poleward deflection through meridional flows (assumed to be polewards at the top of the convection zone and equatorwards at the bottom). In particular, flux tubes which originate between low and intermediate latitudes within the convective overshoot region are strongly affected. This latitude-dependent poleward deflection of erupting magnetic flux renders the wings of stellar butterfly diagrams distinctively convex. The subsequent evolution of the surface magnetic field shows that the increased number of newly emerging bipoles at higher latitudes promotes the intermingling of opposite polarities of polar magnetic fields. The associated magnetic flux densities are about 20 per cent higher than in the case disregarding the pre-eruptive deflection, which eases the necessity for fast meridional flows predicted by previous investigations. In order to reproduce the observed polar field properties, the rate of the meridional circulation has to be of the order of 100 m s⁻¹, and the latitudinal range from which magnetic flux tubes originate at the base of the convective zone (≲50°) must be larger than in the solar case (≲35°).     

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.     

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 Relationship Between Kinematics and Spatial Structure of the Large-Scale Solar Magnetic Field

The rotation of large-scale solar magnetic fields has been investigated by analysing a 20-yr series of synoptic maps of the radial magnetic field. For this purpose, a specially adapted method of spectral analysis was used. We calculated rotation spectra of the magnetic field as functions of the rotation period, heliographic latitude, and longitudinal wave number, k. These spectra reveal the existence of a number of discrete, rigidly rotating components (modes) of the magnetic field, whose rotation periods lie in the wide range from 26.5 to 30.5 days. The significant spectral maxima lie in the (rotation period–latitude) plane close to the curve that represents the differential rotation of small-scale magnetic features. For the first harmonic of the magnetic field (k=1) the properties of the rotation spectra are consistent with those reported by Antonucci, Hoeksema, and Scherrer (1990). However, the distribution of the rigidly rotating modes over rotation period and their latitudinal structure change systematically with the harmonic number k. As k increases, the mean distance P in rotation period between the modes decreases, from 1.2 days for k=1 to 0.3–0.5 days for k=4. This decreasing period separation is accompanied by a decrease of the characteristic latitude separation between the mode maxima. The latitudinal and longitudinal discrete spatial scales of the non-axisymmetric magnetic field appear to be connected with each other, as well as with the temporal scale P.     

The Horizontal Component of Photospheric Plasma Flows During the Emergence of Active Regions on the Sun

The dynamics of horizontal plasma flows during the first hours of the emergence of active region magnetic flux in the solar photosphere have been analyzed using SOHO/MDI data. Four active regions emerging near the solar limb have been considered. It has been found that extended regions of Doppler velocities with different signs are formed in the first hours of the magnetic flux emergence in the horizontal velocity field. The flows observed are directly connected with the emerging magnetic flux; they form at the beginning of the emergence of active regions and are present for a few hours. The Doppler velocities of flows observed increase gradually and reach their peak values 4?–?12 hours after the start of the magnetic flux emergence. The peak values of the mean (inside the ±?500 m?s^?1 isolines) and maximum Doppler velocities are 800?–?970 m?s^?1 and 1410?–?1700 m?s^?1, respectively. The Doppler velocities observed substantially exceed the separation velocities of the photospheric magnetic flux outer boundaries. The asymmetry was detected between velocity structures of leading and following polarities. Doppler velocity structures located in a region of leading magnetic polarity are more powerful and exist longer than those in regions of following polarity. The Doppler velocity asymmetry between the velocity structures of opposite sign reaches its peak values soon after the emergence begins and then gradually drops within 7?–?12 hours. The peak values of asymmetry for the mean and maximal Doppler velocities reach 240?–?460 m?s^?1 and 710?–?940 m?s^?1, respectively. An interpretation of the observable flow of photospheric plasma is given.     

Meridional drift of the large-scale solar magnetic fields in different phases of solar activity

An analysis is made of mean latitudinal profiles of the meridional drift of the large-scale solar magnetic fields. The previously detected equatorward migration of the drift pattern in the course of a cycle is confirmed. Evidence for the existence of a near-equatorial narrow zone ±7° with an equatorward drift with a rate of about 1 m s^-1 is obtained. The study revealed a significant difference in shapes and variations of average drift profiles for the large-scale magnetic field and small magnetic features (Komm, Howard, and Harvey, 1993).     

Quenching of Meridional Circulation in Flux Transport Dynamo Models

Guided by the recent observational result that the meridional circulation of the Sun becomes weaker at the time of the sunspot maximum, we have included a parametric quenching of the meridional circulation in solar dynamo models such that the meridional circulation becomes weaker when the magnetic field at the base of the convection zone is stronger. We find that a flux transport solar dynamo tends to become unstable on including this quenching of meridional circulation if the diffusivity in the convection zone is less than about 2×10¹¹ cm² s⁻¹. The quenching of α, however, has a stabilizing effect and it is possible to stabilize a dynamo with low diffusivity with sufficiently strong α-quenching. For dynamo models with high diffusivity, the quenching of meridional circulation does not produce a large effect and the dynamo remains stable. We present a solar-like solution from a dynamo model with diffusivity 2.8×10¹² cm² s⁻¹ in which the quenching of meridional circulation makes the meridional circulation vary periodically with solar cycle as observed and does not have any other significant effect on the dynamo.     

Simulations of the sun's polar magnetic fields during sunspot cycle 21

Regarding new bipolar magnetic regions as sources of flux, we have simulated the evolution of the radial component of the solar photospheric magnetic field during 1976–1984 and derived the corresponding evolution of the line-of-sight polar fields as seen from Earth. The observed timing and strength of the polar-field reversal during cycle 21 can be accounted for by supergranular diffusion alone, for a diffusion coefficient of 800 km² s^-1. For an assumed 300 km² s^-1 rate of diffusion, on the other hand, a poleward meridional flow with a moderately broad profile and a peak speed of 10 m s^-1 reached at about 5° latitude is required to obtain agreement between the simulated and observed fields. Such a flow accelerates the transport of following-polarity flux to the polar caps, but also inhibits the diffusion of leading-polarity flux across the equator. For flows faster than about 10 m s^-1 the latter effect dominates, and the simulated polar fields reverse increasingly later and more weakly than the observed fields.Laboratory for Computational Physics and Fluid Dynamics.E. O. Hulburt Center for Space Research.     

Sector Structure, Rotation, and Cyclic Evolution of Large-Scale Solar Magnetic Fields

Auto-correlation analysis was performed using digitized synoptic charts of photospheric magnetic fields for the past three solar activity cycles (1965–1994). The obtained correlograms were used to study the rotation and the zonal-sector structure of large-scale solar magnetic fields all over the observable region of heliolatitudes in various phases of solar activity. It is shown that the large-scale system of solar magnetic fields is rather complex and comprises at least three different systems. One is a global rigidly rotating system. It determines the cyclic variation of magnetic fields and is probably responsible for the behavior of magnetic fields in the polar zones. Another is a rigidly rotating 4-sector structure in the central (equatorial and mid-latitude) zone. The third is a differentially rotating system that determines the behavior of the LSSMF structure elements with a size of 30–60° and less. This one is the most noticeable in the central zone and absent in the polar zones. Various cyclic and rotation parameters of the three field structures are discussed.     

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.     

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.     

Stability and equilibrium of emerged magnetic flux

We analyse stability and equilibrium of a unipolar large-scale magnetic field pervading a plane horizontal subphotospheric layer with the possible implications for sunspots in mind. Eddy diffusivity is applied to account for the effects of the small-scale convective turbulence. Diffusivity quenching by magnetic field results in a secondary large-scale instability. A linear stability analysis is performed to define the marginal stability boundary in parametric space and the unstable mode structure. The nonlinear dynamics of the unstable modes are followed numerically. The original state of a uniform vertical magnetic field is transformed via the instability into the nonlinear dynamical equilibrium with a highly intermittant distribution of the magnetic field. Magnetic flux is concentrated in a relatively small area surrounded by an almost field-free region. The role of the fluid motion in the hydromagnetic equilibrium is emphasized. Although the relevance of the instability to the process of sunspot formation is rather questionable, the resulting equilibrium structures are similar to mature spots in their thermal and magnetic properties. Also, the simulated flow structure agrees with helioseismic tomography results.