Development of structure in pores and sunspots: flows around axisymmetric magnetic flux tubes

Flux elements, pores and sunspots form a family of magnetic features observed at the solar surface. As a first step towards developing a fully non-linear model of the structure of these features and of the dynamics of their interaction with solar convection, we conduct numerical experiments on idealized axisymmetric flux tubes in a compressible convecting atmosphere in cylindrical boxes of radius up to 8 times their depth. We find that the magnetic field strength of the flux tubes is roughly independent of both distance from the centre and the total flux content of the flux tube, but that the angle of inclination from the vertical of the field at the edge of the tube increases with flux content. In all our calculations, fluid motion converges on the flux tube at the surface. The results compare favourably with observations of pores; in contrast, large sunspots lie at the centre of an out-flowing moat cell. We conjecture that there is an inflow hidden beneath the penumbrae of large spots, and that this inflow is responsible for the remarkable longevity of such features.     

Numerical simulations of rotating axisymmetric sunspots

A numerical model of axisymmetric convection in the presence of a vertical magnetic flux bundle and rotation about the axis is presented. The model contains a compressible plasma described by the non-linear MHD equations, with density and temperature gradients simulating the upper layer of the Sun's convection zone. The solutions exhibit a central magnetic flux tube in a cylindrical numerical domain, with convection cells forming collar flows around the tube. When the numerical domain is rotated with a constant angular velocity, the plasma forms a Rankine vortex, with the plasma rotating as a rigid body where the magnetic field is strong, as in the flux tube, while experiencing sheared azimuthal flow in the surrounding convection cells, forming a free vortex. As a result, the azimuthal velocity component has its maximum value close to the outer edge of the flux tube. The azimuthal flow inside the magnetic flux tube and the vortex flow is prograde relative to the rotating cylindrical reference frame. A retrograde flow appears at the outer wall. The most significant convection cell outside the flux tube is the location for the maximum value of the azimuthal magnetic field component. The azimuthal flow and magnetic structure are not generated spontaneously, but decay exponentially in the absence of any imposed rotation of the cylindrical domain.     

Magnetic flux separation in photospheric convection

Three-dimensional non-linear magnetoconvection in a strongly stratified compressible layer exhibits different patterns as the strength of the imposed magnetic field is reduced. There is a transition from a magnetically dominated regime, with small-scale convection in slender hexagonal cells, to a convectively dominated regime, with clusters of broad rising plumes that confine the magnetic flux to narrow lanes where fields are locally intense. Both patterns can coexist for intermediate field strengths, giving rise to flux separation: clumps of vigorously convecting plumes, from which magnetic flux has been excluded, are segregated from regions with strong fields and small-scale convection. A systematic numerical investigation of these different states shows that flux separation can occur over a significant parameter range and that there is also hysteresis. The results are related to the fine structure of magnetic fields in sunspots and in the quiet Sun.     

Modelling photospheric magnetoconvection

The increasing power of computers makes it possible to model the non-linear interaction between magnetic fields and convection at the surfaces of solar-type stars in ever greater detail. We present the results of idealized numerical experiments on two-dimensional magnetoconvection in a fully compressible perfect gas. We first vary the aspect ratio λ of the computational box and show that the system runs through a sequence of convective patterns, and that it is only for a sufficiently wide box (λ ≥ 6) that the flow becomes insensitive to further increases in λ. Next, setting λ = 6, we decrease the field strength from a value strong enough to halt convection and find transitions to small-scale steady convection, next to spatially modulated oscillations (first periodic, then chaotic) and then to a new regime of flux separation, with regions of strong field (where convection is almost completely suppressed) separated by broad convective plumes. We also explore the effects of altering the boundary conditions and show that this sequence of transitions is robust. Finally, we relate these model calculations to recent high-resolution observations of solar magnetoconvection, in plage regions as well as in light bridges and the umbrae of sunspots.     

The emergence of magnetic flux

This paper first summarizes the morphology and dynamics of emerging flux regions and arch filament systems and then discusses detailed observations of a particular active region with emerging magnetic flux.The central part of the growing active region shows abnormal granulation and a weak magnetic field that, locally, is transverse. In the border zone, strong downward flows occur in the chromopshere and photosphere (small features with strong magnetic fields (faculae, pores) are formed here.) Near the leading and following edge, sunspots are formed by the coalescence of such small magnetic elements.The observational data are interpreted by means of a heuristic model of an emergent magnetic loop-shaped bundle consisting of many flux tubes. In this model we incorporate the theory of convective collapse and the buoyancy of flux tubes. The observed complexity in the structure and dynamics, including strong transverse fields and velocity shear, is attributed to the emergence of several flux regions within the active region at different orientations.     

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.     

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)     

An accretion disc model with a magnetic wind and turbulent viscosity

A model is presented for an accretion disc with turbulent viscosity and a magnetically influenced wind. The magnetic field is generated by a dynamo in the disc, involving the turbulence and radial shear. Disc-wind solutions are found for which the wind mass flux is sufficient to play a major part in driving an imposed steady inflow, but small enough for most material to be accreted on to the central object. Constraints arise for the magnetic Reynolds and Prandtl numbers in terms of the turbulent Mach number and vertical length-scale of the disc's horizontal magnetic field. It is shown that the imposition of a stellar boundary condition enhances the wind mass flux in the very inner region of the disc and may result in jet formation.     

Turbulent effects of sunspot magnetic field reconstruction

We study the effects of two-dimensional turbulence generated in sunspot umbra due to strong magnetic fields and Alfven oscillations excited in sunspots due to relatively weak magnetic fields on the evolution of sunspots. Two phases of sunspot magnetic field decaying are shown to exist. The initial rapid phase of magnetic field dissipation is due to two-dimensional turbulence. The subsequent slow phase of magnetic field decaying is associated with Alfven oscillations. Our results correspond to observed data that provide evidence for two types of sunspot evolution. The effect of macroscopic diamagnetic expulsion of magnetic field from the convective zone or photosphere toward sunspots is essential in supporting the long-term stability and equilibrium of vertical magnetic flux tubes in sunspots.     

The stability of a velocity shear in the presence of an inhomogeneous magnetic field

The stability of a velocity shear in the presence of a parallel but non-uniform magnetic field is considered in general terms. Two special cases are then investigated; (i) the well known case of a plane interface at which a discontinuity in the magnetic field coincides with the velocity shear; (ii) an axially symmetric flow in which discontinuities in the magnetic and velocity fields occur at a cylindrical surface whose axis is parallel to the flow. In the first case the flow is stabilized if the rms Alfvén velocity of the magnetic field exceeds the shear velocity; a result consistent with that obtained by other writers. In the second case it is shown that the discontinuity in the magnetic field increases the stability of the system. The significance of this result for the stability of the flux ropes associated with sunspots in the solar convection zone is considered.     

On the appearance of magnetic flux in the solar photosphere

Ideas and models for the appearance of photospheric magnetic structure are confronted with observational data. Some findings are: The magnetic flux emerging in an active region consists of a bundle of flux tubes which were already concentrated before penetrating into the photosphere. A model of a rising bunch of flux tubes joining into a few strands at larger depths describes the coalescence of spots near the leading and following edges of the active region while more flux may surface near the center of the region. There is no observational evidence for appreciable helical twists in the flux bundles.Throughout the region's lifetime the magnetic elements move coherently, the whole magnetic structure rotates faster than the quiet photosphere. In active regions the convective flow at scales larger than the granulation is arrested by the magnetic structure. The long-lived supergranular cells around spots and in the enhanced network in turn determine the decay properties of spots and facular clusters. The modulation of the convective flow by the magnetic structure explains the slow dispersal of faculae.The hierarchy of magnetic elements (sunspots-pores-knots-facular clusters-facular points) may be explained by a set of magnetostatic flux tube models in the top of the convection zone. The underlying assumptions are that the heat flow along the magnetic field is reduced and that there is no heat exchange across the field except by radiation.A tentative model is proposed to account for the amplification, ascent and emergence of intense flux bundles. The assumptions are: (i) the field is concentrated in toroidal bundles by differential rotation, (ii) in the deep convection zone flux bundles are contained by the external turbulent pressure, and (iii) for field strengths up to the equipartition value efficient lateral heat exchange is possible. After a loop has surfaced radiative cooling and subsequent convective downflow reduce the temperature in the top of the flux tubes which then contract to field strengths well above the local equipartition value. There the heat flow is channelled along the field, which creates the conditions for the magnetostatic flux tube models without requiring a blocking of the heat flow somewhere within the tubes.The paper contains a brief review on the evolution of the magnetic field from the emergence in active regions up to the enigmatic disappearance, and a list of topics for further observational investigation.     

Local Helioseismology of Sunspots: Current Status and Perspectives

Mechanisms of the formation and stability of sunspots are among the longest-standing and intriguing puzzles of solar physics and astrophysics. Sunspots are controlled by subsurface dynamics, hidden from direct observations. Recently, substantial progress in our understanding of the physics of the turbulent magnetized plasma in strong-field regions has been made by using numerical simulations and local helioseismology. Both the simulations and helioseismic measurements are extremely challenging, but it is becoming clear that the key to understanding the enigma of sunspots is a synergy between models and observations. Recent observations and radiative MHD numerical models have provided a convincing explanation for the Evershed flows in sunspot penumbrae. Also, they lead to the understanding of sunspots as self-organized magnetic structures in the turbulent plasma of the upper convection zone, which are maintained by a large-scale dynamics. Local helioseismic diagnostics of sunspots still have many uncertainties, some of which are discussed in this review. However, there have been significant achievements in resolving these uncertainties, verifying the basic results by new high-resolution observations, testing the helioseismic techniques by numerical simulations, and comparing results obtained by different methods. For instance, a recent analysis of helioseismology data from the Hinode space mission has successfully resolved several uncertainties and concerns (such as the inclined-field and phase-speed filtering effects) that might affect the inferences of the subsurface wave-speed structure of sunspots and the flow pattern. It is becoming clear that for the understanding of the phenomenon of sunspots it is important to further improve the helioseismology methods and investigate the whole life cycle of active regions, from magnetic flux emergence to dissipation. The Solar Dynamics Observatory mission has started to provide data for such investigations.     

A new model for flux emergence and the evolution of sunspots and the large-scale fields

Existing models for the evolution of sunspots and sunspot groups, describing the subsurface structure of the magnetic fields and their interactions with the convective motions, are briefly reviewed. It is shown that they are generally unable to account for the most recent data concerning the relationship between the large-scale solar magnetic field structures and the magnetic fields of active regions. In particular, it is shown that the former do not arise directly from the decay of the latter, as required by the Babcock model and all other models based on it. Other observations which are not adequately explained by current models are also cited.A new model is put forward based on the expulsion of toroidal magnetic flux by the dominant (i.e. giant) cells of the convection zone. The flux expelled above these cells forms the large-scale field and thus the configuration of this field provides a clue to the structure of the giant cell patterns. The flux expelled below the cells becomes twisted into a rope as in the Babcock model but a loop or stitch forms only in the region of upflow of the giant cells. The interaction of this loop with intermediate-sized cells as it rises to the surface determines the configuration and extent of the active region which appears at the surface. The compatibility of the model with other observations is discussed and its implications for theories of the solar cycle are noted.     

Total resonant absorption of acoustic oscillations in sunspots

The question of total resonant absorption of acoustic oscillations in sunspots is studied for cylindrical 1-D flux tubes that are stratified only in the radial direction and surrounded by a uniform, non-magnetic plasma. The numerical investigation of Goossens and Poedts (1992) in linear resistive MHD is taken further by increasing the strength of the azimuthal magnetic field in the equilibrium flux tubes. For relatively strong azimuthal magnetic fields, total absorption is found over a relatively wide range of spot radii.     

Temporal Evolution of Sunspot Areas and Estimation of Related Plasma Flows

The increased amount of information provided by ongoing missions such as the Solar Dynamics Observatory (SDO) represents a great challenge for the understanding of basic questions such as the internal structure of sunspots and how they evolve with time. Here, we contribute with the exploitation of new data, to provide a better understanding of the separate growth and decay of sunspots, umbra, and penumbra. Using fuzzy sets to compute separately the areas of sunspot umbra and penumbra, the growth and decay rates for active regions NOAA 11117, NOAA 11428, NOAA 11429, and NOAA 11430 are computed from the analysis of intensitygrams obtained by the Helioseismic and Magnetic Imager onboard SDO. A simplified numerical model is proposed for the decay phase, whereby an empirical irrotational and uniformly convergent horizontal velocity field interacting with an axially symmetric and height-invariant magnetic field reproduces the large-scale features of the much more complex convection observed inside sunspots.     

Flows around sunspots and pores

We report on three sequences of high-resolution white-light and magnetogram observations obtained in the summer of 1989. The duration of sub-arcsecond seeing was three to four hours on each day. Study of the white-light and magnetogram data yields the following results:

1. For all but one of the sunspots we have observed, both dark fibrils and bright grains in the inner part of the penumbra of sunspots move toward the umbra with a speed of about 0.5 km s^-1. In the outer part of the penumbra, movement is away from the umbra. The one exception is a newly formed spot, which has inflow only in its penumbra.
2. Granular flows converge toward almost every pore, even before its formation. Pores are observed to form by the concentration of magnetic flux already existing in the photosphere. The pores (or small sunspots), in turn, then move and concentrate to form bigger sunspot.
3. We followed an emerging flux region (EFR) from 29 to 31 July, 1989 that was composed of a large number of bipoles with magnetic polarities mixed over a large area in the first day of its birth. As time went on, polarities sorted out: the leading polarity elements moved in one direction; the following, the opposite. During the process a large number of cancellations occurred, with some sub-flares and surges observed simultaneously. After about 24 hours, the positive and negative fluxes were essentially separated.
4. We find two kinds of photospheric dark alignments in the region of new flux emergence: (a) alignments connecting two poles of opposite magnetic polarity form the tops of rising flux tubes; (b) alignments corresponding to the magnetic flux of one polarity, which we call elongated pores.

     

Convective intensification of magnetic fields in the quiet Sun

Kilogauss-strength magnetic fields are often observed in intergranular lanes at the photosphere in the quiet Sun. Such fields are stronger than the equipartition field B _e, corresponding to a magnetic energy density that matches the kinetic energy density of photospheric convection, and comparable with the field B _p that exerts a magnetic pressure equal to the ambient gas pressure. We present an idealized numerical model of three-dimensional compressible magnetoconvection at the photosphere, for a range of values of the magnetic Reynolds number. In the absence of a magnetic field, the convection is highly supercritical and characterized by a pattern of vigorous, time-dependent, 'granular' motions. When a weak magnetic field is imposed upon the convection, magnetic flux is swept into the convective downflows where it forms localized concentrations. Unless this process is significantly inhibited by magnetic diffusion, the resulting fields are often much greater than B _e and the high magnetic pressure in these flux elements leads to their being partially evacuated. Some of these flux elements contains ultraintense magnetic fields that are significantly greater than B _p. Such fields are contained by a combination of the thermal pressure of the gas and the dynamic pressure of the convective motion, and they are constantly evolving. These ultraintense fields develop owing to non-linear interactions between magnetic fields and convection; they cannot be explained in terms of 'convective collapse' within a thin flux tube that remains in overall pressure equilibrium with its surroundings.     

On the out of phase appearance of large-scale diffuse magnetic field of the Sun with respect to sunspots

We assume the large-scale diffuse magnetic field of the Sun to originate from the poloidal component of a dynamo operating at the base of the convection zone, whereas the sunspots are due to the toroidal component. The evolution of the poloidal component is studied to model the poleward migration of the diffuse field seen on the solar surface and the polar reversal at the time of sunspots maxima (Dikpati and Choudhuri 1994, 1995).     

Equilibrium and decay of sunspots

We suggest a quantitative sunspot model developed in terms of mean-field magnetohydrodynamics (MHD). The model consistently describes the distributions of magnetic field, fluid velocity, and thermodynamic parameters in a sunspot and the surrounding matter. Two versions of the model allow the MHD equilibrium in sunspots and their slow decay to be analyzed. The baroclinic flow converging to the sunspot plays an important role in the equilibrium. Several calculated characteristics—almost uniform distributions of brightness and magnetic field inside sunspots, their abrupt changes at the boundary, and nearly linear decreases in the area and magnetic flux of decaying sunspots with time—qualitatively agree with the observations.