On Magnetic Coupling Between the Two Hemispheres in Solar Dynamo Models

By introducing an asymmetry between the two hemispheres, we study whether the solar dynamo solutions in the two hemispheres remain coupled with each other. Our calculations are based on the solar dynamo code SURYA, which incorporates the helioseismically-determined solar-rotation profile, a Babcock–Leighton α effect concentrated near the surface, and a meridional circulation. When the magnetic coupling between the hemispheres is enhanced by either increasing the diffusion or introducing an α effect distributed throughout the convection zone, we find that the solutions in the two hemispheres evolve together with a single period even when we make the meridional circulation or the α effect different in the two hemispheres. On the other hand, when the hemispheric coupling is weaker for other values of parameters, an asymmetry between the hemispheres can make solutions in the two hemispheres evolve independently with different periods.     

On the possibility of a bimodal solar dynamo

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

Modeling a Shallow Solar Dynamo

Photospheric ephemeral regions (EPRs) cover the Sun like a magnetic carpet. From this, we update the Babcock – Leighton solar dynamo. Rather than sunspot fields appearing in the photosphere de novo from eruptions originating in the deep interior, we consider that sunspots form directly in the photosphere by a rapid accumulation of like-sign field from EPRs. This would only occur during special circumstances: locations and times when the temperature structure is highly superadiabatic and contains a large subsurface horizontal magnetic field (only present in the Sun’s lower latitudes). When these conditions are met, superadiabatic percolation occurs, wherein an inflow and downflow of gas scours the surface of EPRs to form active regions. When these conditions are not met, magnetic elements undergo normal percolation, wherein magnetic elements move about the photosphere in Brownian-type motions. Cellular automata (CA) models are developed that allow these processes to be calculated and thereby both small-scale and large-scale models of magnetic motions can be obtained. The small-scale model is compared with active region development and Hinode observations. The large-scale CA model offers a solar dynamo, which suggests that fields from decaying bipolar magnetic regions (BMRs) drift on the photosphere driven by subsurface magnetic forces. These models are related to observations and are shown to support Waldmeier’s findings of an inverse relationship between solar cycle length and cycle size. Evidence for significant amounts of deep magnetic activity could disprove the model presented here, but recent helioseismic observations of “butterfly patterns” at depth are likely just a reflection of surface activity. Their existence seems to support the contention made here that the field and flow separate, allowing cool, relatively field-free downdrafts to descend with little field into the nether worlds of the solar interior. There they heat by compression to form a hot solar-type Santa Ana wind deep below active regions.     

Intensity and Magnetic Field Distribution of Sunspots

We study the relationship between the brightness (I) and magnetic field (B) distributions of sunspots using 272 samples observed at the San Fernando Observatory and the National Solar Observatory, Kitt Peak, whose characteristics varied widely. We find that the I – B relationship has a quadratic form for the spots with magnetic field less than about 2000 G. The slope of the linear part of the I – B curve varies by about a factor of three for different types of spots. In general the slope increases as the spot approaches disk center. The I – B slope does not have a clear dependency on the spot size but the lower limit appears to increase as a function of the ratio of umbra and penumbra area. The I – B slope changes as a function of age of the sunspots. We discuss various sunspot models using these results.     

A Study of Surges: II. On the Relationship between Chromospheric Surges and Coronal Mass Ejections

Liu et al. (Astrophys. J. 628, 1056, 2005a) described one surge – coronal mass ejection (CME) event showing a close relationship between solar chromospheric surge ejection and CME that had not been noted before. In this work, large Hα surges (>72 Mm, or 100 arcsec) are studied. Eight of these were associated with CMEs. According to their distinct morphological features, Hα surges can be classified into three types: jetlike, diffuse, and closed loop. It was found that all of the jetlike surges were associated with jetlike CMEs (with angular widths ≤30 degrees); the diffuse surges were all associated with wide-angle CMEs (e.g., halo); the closed-loop surges were not associated with CMEs. The exclusive relation between Hα surges and CMEs indicates difference in magnetic field configurations. The jetlike surges and related narrow CMEs propagate along coronal fields that are originally open. The unusual transverse mass motions in the diffuse surges are suggested to be due to magnetic reconnections in the corona that produce wide-angle CMEs. For the closed-loop surges, their paths are just outlining stable closed loops close to the solar surface. Thus no CMEs are associated with them.     

On The Torsional Oscillations In Babcock–Leighton Solar Dynamo Models

The torsional oscillations at the solar surface have been interpreted by Schüssler and Yoshimura as being generated by the Lorentz force associated with the solar dynamo. It has been shown recently that they are also present in the upper half of the solar convection zone (SCZ). With the help of a solar dynamo model of the Babcock–Leighton type studied earlier, the longitudinal component of the Lorentz force, L , is calculated, and its sign or isocontours, are plotted vs. time, t, and polar angle, (the horizontal and vertical axis respectively). Two cases are considered, (1) differential rotation differs from zero only in the tachocline, (2) differential rotation as in (1) in the tachocline, and purely latitudinal and independent of depth in the bulk of the SCZ. In the first case the sign of L is roughly independent of latitude (corresponding to vertical bands in the t, plot), whereas in the second case the bands show a pole–equator slope of the correct sign. The pattern of the bands still differs, however, considerably from that of the helioseismic observations, and the values of the Lorentz force are too small at low latitudes. It is all but certain that the toroidal field that lies at the origin of the large bipolar magnetic regions observed at the surface, must be generated in the tachocline by differential rotation; the regeneration of the corresponding poloidal field, B _p has not yet been fully clarified. B _p could be regenerated, for example, at the surface (as in Babcock–Leighton models), or slightly above the tachocline, (as in interface dynamos). In the framework of the Babcock-Leighton models, the following scenario is suggested: the dynamo processes that give rise to the large bipolar magnetic regions are only part of the cyclic solar dynamo (to distinguish it from the turbulent dynamo). The toroidal field generated locally by differential rotation must contribute significantly to the torsional oscillations patterns. As this field becomes buoyant, it should give rise, at the surface, to the smaller bipolar magnetic regions as, e.g., to the ephemeral bipolar magnetic regions. These have a weak non-random orientation of magnetic axis, and must therefore also contribute to the source term for the poloidal field. Not only the ephemeral bipolar regions could be generated in the bulk of the SCZ, but many of the smaller bipolar regions as well (at depths that increase with their flux), all contributing to the source term for the poloidal field. In contrast to the butterfly diagram that provides only a very weak test of dynamo theories, the pattern of torsional oscillations has the potential of critically discriminating between different dynamo models.     

On the Compatibility of a Flux Transport Dynamo with a Fast Tachocline Scenario

The compatibility of the fast-tachocline scenario with a flux-transport dynamo model is explored. We employ a flux-transport dynamo model coupled with simple feedback formulae relating the thickness of the tachocline to the amplitude of the magnetic field or to the Maxwell stress. The dynamo model is found to be robust against the nonlinearity introduced by this simplified fast-tachocline mechanism. Solar-like butterfly diagrams are found to persist and, even without any parameter fitting, the overall thickness of the tachocline is well within the range admitted by helioseismic constraints. In the most realistic case of a time- and latitude-dependent tachocline thickness linked to the value of the Maxwell stress, both the thickness and its latitudinal dependence are in excellent agreement with seismic results. In nonparametric models, cycle-related temporal variations in tachocline thickness are somewhat larger than admitted by helioseismic constraints; we find, however, that introducing a further parameter into our feedback formula readily allows further fine tuning of the thickness variations.     

Solar Cycle Variation of the Interplanetary Forward Shock Drivers Observed at 1 AU

Forecasting space weather more accurately from solar observations requires an understanding of the variations in physical properties of interplanetary (IP) shocks as solar activity changes. We examined the characteristics (occurrence rate, physical parameters, and types of shock driver) of IP shocks. During the period of 1995 – 2001, a total of 249 forward IP shocks were observed. In calculating the shock parameters, we used the solar wind data from Wind at the solar minimum period (1995 – 1997) and from ACE since 1998 including the solar maximum period (1999 – 2001). Most of IP shocks (68%) are concentrated in the solar maximum period. The values of physical quantities of IP shocks, such as the shock speed, the sonic Mach number, and the ratio of plasma density compression, are larger at solar maximum than at solar minimum. However, the ratio of IMF compression is larger at solar minimum. The IP shock drivers are classified into four groups: magnetic clouds (MCs), ejecta, high speed streams (HSSs), and unidentified drivers. The MC is the most dominant and strong shock driver and 150 out of total 249 IP shocks are driven by MCs. The MC is a principal and very effective shock driver not only at solar maximum but also at solar minimum, in contrast to results from previous studies, where the HSS is considered as the dominant IP shock driver.     

The Solar Optical Telescope of Solar-B (Hinode): The Optical Telescope Assembly

The Solar Optical Telescope (SOT) aboard the Solar-B satellite (Hinode) is designed to perform high-precision photometric and polarimetric observations of the Sun in visible light spectra (388 – 668 nm) with a spatial resolution of 0.2 – 0.3 arcsec. The SOT consists of two optically separable components: the Optical Telescope Assembly (OTA), consisting of a 50-cm aperture Gregorian with a collimating lens unit and an active tip-tilt mirror, and an accompanying Focal Plane Package (FPP), housing two filtergraphs and a spectro-polarimeter. The optomechanical and optothermal performance of the OTA is crucial to attain unprecedented high-quality solar observations. We describe in detail the instrument design and expected stable diffraction-limited on-orbit performance of the OTA, the largest state-of-the-art solar telescope yet flown in space.     

Return Meridional Flow in the Convection Zone from Latitudinal Motions of Umbrae of Sunspot Groups

We have derived the velocities of meridional flows by measuring the latitudinal motions (or drifts) of umbrae of spot groups classified into three categories of area: 0 – 5 μ, 5 – 10 μ, and >10 μ (μ area in millionths of the solar hemisphere). The latitudinal drifts (or the meridional flows) in all three categories are directed equatorward in both the northern and southern hemispheres. By sorting the spot groups into three area classes, we are able to relate the respective latitudinal drifts with the three depths in the convection zone where the footpoints of the flux loops of the spot groups of each area class are anchored. We obtain estimates of the anchor depths through a comparison of the rotation rates of the spot groups of each area class with the rotation-rate profiles from helioseismic inversions. The equatorward drifts obtained provide estimates of the meridional flows at the three depths in the convection zone and thereby suggest the presence of return meridional flows as envisaged in the flux-transport dynamo models, which have remained undetected so far. The data sources for this study are measurements of positions and areas of umbrae of sunspots from the photographic white-light images of the Sun of the Kodaikanal Observatory archives for the period 1906 – 1987 and a very similar, but independent, data set from the Mt. Wilson Observatory archives for the period 1917 – 1985.     

A two-layer αω-dynamo model with dynamic feedback on the ω-effect

In an attempt to produce a simple representation of an interface dynamo, I examine a dynamo model composed of two one-dimensional (radially averaged) pseudo-spherical layers, one in the convection zone and possessing an α-effect, and the other in the tachocline and possessing an ω-effect. The two layers communicate by means of an analogue of Newton's law of cooling, and a dynamical back-reaction of the magnetic field on ω is provided. Extensive bifurcation diagrams are calculated for three separate values of η, the ratio of magnetic diffusivities of the two layers. I find recognizable similarities to, but also dramatic differences from, the comparable one-layer model examined by Roald &38; Thomas. In particular, the solar-like dynamo mode found previously is no longer stable in the two-layer version; in its place there is a sequence of periodic, quasi-periodic and chaotic modes probably created in a homoclinic bifurcation. These differences are important enough to provide support for the view that the solar dynamo cannot be meaningfully modelled in one dimension.     

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)     

RHESSI Microflares: II. Implications for Loop Structure and Evolution

We present simple analytic models which predict the peak X-ray emission measure and temperature attained in flares in which the chromospheric evaporation process takes place either in a single ‘monolithic’ loop or in a loop consisting of several filaments that are created successively as the energy release process proceeds in time. As possible mechanisms driving chromospheric evaporation we consider both classical heat conduction from the loop top and non-thermal electron beams. The model predictions are tested for a set of 18 well studied RHESSI microflares. The results suggest beam driven evaporation in filamented loops as being capable of accounting for the observed emission measures and temperatures though there are issues with the very high beam densities needed. On the other hand, estimates of the emission measures achieved by conductive evaporation which are derived by using the Rosner – Tucker – Vaiana (RTV) scaling law are much larger than the observed ones. Possible reasons for this discrepancy are discussed.     

Predicting solar ‘climate’ by assimilating magnetic data into a flux-transport dynamo

Observational and theoretical knowledge about global-scale solar dynamo ingredients have reached the stage that it is possible to calibrate a flux-transport dynamo for the Sun by adjusting only a few tunable parameters. The important ingredients in this class of model are differential rotation (Omega-effect), helical turbulence (alpha-effect), meridional circulation and turbulent diffusion. The meridional circulation works as a conveyor belt and governs the dynamo cycle period. Meridional circulation and magnetic diffusivity together govern the memory of the Sun's past magnetic fields. After describing the physical processes involved in a flux-transport dynamo, we will show that a predictive tool can be built from it to predict mean solar cycle features by assimilating magnetic field data from previous cycles. We will discuss the theoretical and observational connections among various predictors, such as dynamo-generated toroidal flux integral, cross-equatorial flux, polar fields and geomagnetic indices. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Solar activity forecast with a dynamo model

Although systematic measurements of the Sun's polar magnetic field exist only from mid-1970s, other proxies can be used to infer the polar field at earlier times. The observational data indicate a strong correlation between the polar field at a sunspot minimum and the strength of the next cycle, although the strength of the cycle is not correlated well with the polar field produced at its end. This suggests that the Babcock–Leighton mechanism of poloidal field generation from decaying sunspots involves randomness, whereas the other aspects of the dynamo process must be reasonably ordered and deterministic. Only if the magnetic diffusivity within the convection zone is assumed to be high (of order  10¹² cm² s⁻¹  ), we can explain the correlation between the polar field at a minimum and the next cycle. We give several independent arguments that the diffusivity must be of this order. In a dynamo model with diffusivity like this, the poloidal field generated at the mid-latitudes is advected toward the poles by the meridional circulation and simultaneously diffuses towards the tachocline, where the toroidal field for the next cycle is produced. To model actual solar cycles with a dynamo model having such high diffusivity, we have to feed the observational data of the poloidal field at the minimum into the theoretical model. We develop a method of doing this in a systematic way. Our model predicts that cycle 24 will be a very weak cycle. Hemispheric asymmetry of solar activity is also calculated with our model and compared with observational data.     

A Simple Model to Describe Solar Cycle Periodicities below 11 Years

Together with the main 11-year cycle, solar activity also displays intracycle periodicities. A simple nonlinear model that describes the 11-year solar cycle with subperiodicities can be derived from the usual α – ω dynamo theory in the form of a Van der Pol equation with a forcing term. In this paper the results obtained from the Van der Pol oscillator describing the amplitude modulations and periodicities observed from the data set of the global daily coronal emission of the Fe xiv line at 530.3 nm are presented.     

A dynamo model for axisymmetric and non-axisymmetric solar magnetic fields

More and more observations are showing a relatively weak, but persistent, non-axisymmetric magnetic field co-existing with the dominant axisymmetric field on the Sun. Its existence indicates that the non-axisymmetric magnetic field plays an important role in the origin of solar activity. A linear non-axisymmetric  α²– Ω  dynamo model is derived to explore the characteristics of the axisymmetric  ( m = 0)  and the first non-axisymmetric  ( m = 1)  modes and to provide a theoretical basis with which to explain the 'active longitude', 'flip-flop' and other non-axisymmetric phenomena. The model consists of an updated solar internal differential rotation, a turbulent diffusivity varying with depth, and an α-effect working at the tachocline in a rotating spherical system. The difference between the  α²–Ω  and the  α–Ω  models and the conditions that favour the non-axisymmetric modes under solar-like parameters are also presented.     

Simulating Solar Cycles in Northern and Southern Hemispheres by Assimilating Magnetic Data into a Calibrated Flux-Transport Dynamo

We use the flux-transport dynamo prediction scheme introduced by Dikpati, de Toma, and Gilman (Geophys. Res. Lett. 33, L05102, 2006) to make separate simulations and predictions of sunspot cycle peaks for northern and southern hemispheres. Despite the division of the data, the skill level achieved is only slightly lower than that achieved for the sum of both hemispheres. The model shows skill at simulating and predicting the difference in peaks between North and South, provided that difference is more than a few percent. The simulation and prediction skill is achieved without adjustment to any parameters of the model that were used when peaks for the sum of North and South sunspot areas was simulated. The results are also very insensitive to the averaging length applied to the input data, provided the simulations and predictions are for peaks defined by averaging the observations over at least 13 rotations. However, in its present form, the model is not capable of skillfully simulating or predicting short-time-scale features of individual solar cycles.     

Differential rotation and meridional flow in the solar convection zone and beneath

The influence of the basic rotation on anisotropic and inhomogeneous turbulence is discussed in the context of differential rotation theory. An improved representation for the original turbulence leads to a Λ‐effect which complies with the results of 3D numerical simulations. The resulting rotation law and meridional flow agree well with both the surface observations (∂Ω/∂r < 0 and meridional flow towards the poles) and with the findings of helioseismology. The computed equatorward flow at the bottom of convection zone has an amplitude of about 10 m/s and may be significant for the solar dynamo. The depth of the meridional flow penetration into the radiative zone is proportional to ν^0.5_core, where ν_core is the viscosity beneath the convection zone. The penetration is very small if the tachocline is laminar. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)