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)     

Interpreting Helioseismic Structure Inversion Results of Solar Active Regions

Helioseismic techniques such as ring-diagram analysis have often been used to determine the subsurface structural differences between solar active and quiet regions. Results obtained by inverting the frequency differences between the regions are usually interpreted as the sound-speed differences between them. These in turn are used as a measure of temperature and magnetic-field strength differences between the two regions. In this paper we first show that the “sound-speed” difference obtained from inversions is actually a combination of sound-speed difference and a magnetic component. Hence, the inversion result is not directly related to the thermal structure. Next, using solar models that include magnetic fields, we develop a formulation to use the inversion results to infer the differences in the magnetic and thermal structures between active and quiet regions. We then apply our technique to existing structure inversion results for different pairs of active and quiet regions. We find that the effect of magnetic fields is strongest in a shallow region above 0.985R _⊙ and that the strengths of magnetic-field effects at the surface and in the deeper (r<0.98R _⊙) layers are inversely related (i.e., the stronger the surface magnetic field the smaller the magnetic effects in the deeper layers, and vice versa). We also find that the magnetic effects in the deeper layers are the strongest in the quiet regions, consistent with the fact that these are basically regions with weakest magnetic fields at the surface. Because the quiet regions were selected to precede or follow their companion active regions, the results could have implications about the evolution of magnetic fields under active regions.     

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)     

Probing signatures of the alpha-effect in the solar convection zone

An attempt to extract maximum information on signatures of the alpha-effect from current helicity and twist density calculations in the solar photosphere is carried out. A possible interpretation of the results for developing the dynamo theory is discussed. The analysis shows that the surface magnetic current helicity is mainly negative/positive in the northern/southern hemispheres of the Sun. This indicates the actual alpha-effect at the photospheric level to be positive/negative, respectively. However, at the bottom of the convection zone, we may assume this effect to change the sign to negative/positive. We reveal some quantities related to the alpha-effect and discuss its spatial and temporal distribution. It is also found that there are a small number of active regions where the sign of the alpha-effect is opposite to that in most active regions. Such exceptional active regions seem to localize at certain active longitudes. We compare the determined regularities with theoretical predictions of the alpha-effect distribution in the solar convection zone.     

Imaging bright-spots in the accretion flow near the black hole horizon of Sgr A*

Solar active regions are distinguished by their strong magnetic fields. Modern local helioseismology seeks to probe them by observing waves which emerge at the solar surface having passed through their interiors. We address the question of how an acoustic wave from below is partially converted to magnetic waves as it passes through a vertical magnetic field layer where the sound and Alfvén speeds coincide (the equipartition level), and find that (i) there is no associated reflection at this depth, either acoustic or magnetic, only transmission and conversion to an ongoing magnetic wave; and (ii) conversion in active regions is likely to be strong, though not total, at frequencies typically used in local helioseismology, with lower frequencies less strongly converted. A simple analytical formula is presented for the acoustic-to-magnetic conversion coefficient.     

On the Relationship between Flux Emergence and CME Initiation

We present in this paper a statistical study aimed at understanding the possible relationship between surface magnetic field variation and CME initiation. The three samples studied comprise 189 CME-source regions, 46 active regions, and 15 newly emerging active regions. Both large-scale and small-scale variations of longitudinal magnetic fields of these regions are studied. To quantitatively study these variations, three physical quantities are calculated: the average total magnetic flux (ATF), the flux variation rate (FVR), and the normalized flux variation rate (NFVR). Our results show that 60% of the CME-source regions are found to have magnetic flux increases during 12 hours before CME eruptions and 40% are found to have magnetic flux decreases. The NFVR of CME-source regions are found to be statistically identical to those of active regions, averaged over 111 hours, and significantly smaller than those of newly emerging active regions. In addition 91% of the CME-source regions are found to have small-scale flux emergence, whereas small-scale flux emergences are also easily identified in active regions during periods with no solar surface activity. Our study suggests that the relationship between flux emergence and CME eruption is complex and the appearance of flux emergence alone is not unique for the initiation of CME eruption.     

On the solar rotation and activity

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

太阳活动区纵向电流密度的计算

本文利用太阳活动区光球横向磁场观测资料推算纵向电流密度分布,论述了具体的计算方法和取得的结果,并简要讨论了太阳活动区电流计算在太阳物理研究中的应用。     

The 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)     

Force-free magnetic fields in the solar atmosphere

Reliable measurements of the solar magnetic field are restricted to the level of the photosphere. For about half a century attempts have been made to calculate the field in the layers above the photosphere, i.e. in the chromosphere and in the corona, from the measured photospheric field. The procedure is known as magnetic field extrapolation. In the superphotospheric parts of active regions the magnetic field is approximately force-free, i.e. electric currents are aligned with the magnetic field. The practical application to solar active regions has been largely confined to constant-α or linear force-free fields, with a spatially constant ratio, α, between the electric current and the magnetic field. We review results obtained from extrapolations with constant-α force-free fields, in particular on magnetic topologies favourable for flares and on magnetic and current helicities. Presently, different methods are being developed to calculate non-constant-α or nonlinear force-free fields from photospheric vector magnetograms. We also briefly discuss these methods and present a comparison of a linear and a nonlinear force-free magnetic field extrapolation applied to the same photospheric boundary data. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Study of polar faculae by means of spectro‐polarimetric observations and radiative transfer calculations

Polar faculae are of special interest for solar physics because of their close relationship to the global magnetic field of the Sun and to solar activity, and because of the recently found kilogauss magnetic fields, which are very unusual for the structures outside active regions at high latitudes of the Sun. The idea is that polar faculae can be represented by bundles of unresolved small‐scale magnetic flux tubes, which are characterized by sizes of about 100 km and strong magnetic fields. High resolution spectro‐polarimetric observations of the considered structures were performed and complemented by the radiation transfer calculations with oblique rays passing through an inhomogeneous magnetic medium. The recent results of observations and numerical calculations are presented.     

3-D non-linear evolution of a magnetic flux tube in a spherical shell: The isentropic case

We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code.We seek to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, especially when the tube is introduced at low latitudes. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Observation of solar differential rotation with the aid of magnetic tracers

We tried to search for the manifestation of differential rotation in the distribution of weak remnants of magnetic fields measured with a very low resolution. We found that, during the periods of low solar activity and in parts of the solar photosphere with smaller density of new magnetic flux sources, it was possible to observe the distribution of magnetic tracers in the form of differential rotation parabolas which increase their curvature from one rotation to the next. The obtained differential rotation rates are not far from those given by highly averaged sunspot data or by the daily magnetic fields. The characteristic differential rotation parabolas as well as specific cellular-like features disturbing their smooth patterns are always formed from fields of one main polarity, the sign of which depends on the phase of the activity cycle.Solar Cycle Workshop Paper.     

Formation of Solar Delta Active Regions： Twist and Writhe of Magnetic Ropes

We analyze the process of formation of delta configuration in some well-known super active regions based on photospheric vector magnetogram observations. It is found that the magnetic field in the initial developing stage of some delta active regions shows a potential-like configuration in the solar atmosphere, the magnetic shear develops mainly near the magnetic neutral line with magnetic islands of opposite polarities, and the large-scale photospheric twisted field forming gradually later. Some results are obtained: (1) The analysis of magnetic writhe of whole active regions cannot be limited in the strong field of sunspots, because the contribution of the fraction of decayed magnetic field is non-negligible. (2) The magnetic model of kink magnetic ropes, supposed to be generated in the subatmosphere, is not consistent with the evolution of large-scale twisted photospheric transverse magnetic field and not entirely consistent with the relationship with magnetic shear in some delta active regions. (3) T     

Signature of differential rotation in solar disk‐integrated chromospheric line emission

UARS SOLSTICE data have been subjected to Fourier and wavelet analyses in order to search for the signature of the solar rotation law in the disk‐integrated irradiance of UV lines. Lyman‐α, Mg II, and Ca II data show a different behaviour. In the SOLSTICE data there are significant temporal variations of the rotation rate of the UV tracers over 5—6 years. Often several distinct rotation periods appear almost simultaneously. Beside the basic period around 27 days there are signals at 32—35 days corresponding to the rotation rate at very high latitudes. For more than 5 years during another period of the solar cycle the rotational behaviour is quite different; there is an indication of differential rotation of active regions in these Ca II ground‐based data. The data contain a wealth of information about the solar differential rotation, but it proves difficult to disentangle the effects of the different emitting sources.     

Distribution of magnetic flux on the solar surface and low-degree p-modes

The frequencies of solar p-modes are known to change over the solar cycle. There is also recent evidence that the relation between frequency shift of low-degree modes and magnetic flux or other activity indicators differs between the rising and falling phases of the solar cycle, leading to a hysteresis in such diagrams. We consider the influence of the changing large-scale surface distribution of the magnetic flux on low-degree ( l ≤3) p-mode frequencies. To that end, we use time-dependent models of the magnetic flux distribution and study the ensuing frequency shifts of modes with different order and degree as a function of time. The resulting curves are periodic functions (in simple cases just sine curves) shifted in time by different amounts for the different modes. We show how this may easily lead to hysteresis cycles comparable to those observed. Our models suggest that high-latitude fields are necessary to produce a significant difference in hysteresis between odd- and even-degree modes. Only magnetic field distributions within a small parameter range are consistent with the observations by Jiménez-Reyes et al. Observations of p-mode frequency shifts are therefore capable of providing an additional diagnostic of the magnetic field near the solar poles. The magnetic distribution that is consistent with the p-mode observations also appears reasonable compared with direct measurements of the magnetic field.     

Rotation of solar structures in the upper chromosphere. II. Time variations in the latitudinal distribution of the rotation of active regions and coronal holes

The time variations in the latitudinal distribution of the rotation of active regions and coronal holes are investigated. The synoptic maps obtained from observations in the He I 1083 nm line at Kitt Peak Observatory over almost three solar cycles are used as observational data. A Fourier analysis of the time series constructed from synoptic maps has yielded the following results. The rotation of active regions differs significantly from the rotation of coronal holes in all parameters: the set of the most significant rotation periods, their latitudinal distribution, and time variations. The rotation of active regions and coronal holes is characterized by variations from cycle to cycle, a time-varying north-south asymmetry. The power spectra for consecutive cycles of solar activity differ significantly for both epochs of high activity and minima. Analysis of the total power of the spectra within four selected intervals of periods from 21 to 33 days has shown that the total power is highest in the intervals of periods 24–27 and 27–30 days. This is valid for both active regions and coronal holes. The correlation between the total powers in the above intervals of periods changes noticeably with time. Long-lived or successively appearing active regions with rotation periods in the range 24–30 days are typical of the time of a sharp decrease in the total equivalent width of active regions. This includes not only the decline time of the 11-year cycles, but also the minima between recurrent activity maxima during one cycle. A predominance of long-lived coronal holes as their total equivalent width decreases is noticeable for coronal holes with rotation periods in the interval 30–33 days. All of the above results suggest that the rotation of solar structures is determined mainly by the subphotospheric sources of specific structures, not by the rotation of the main volumes of solar plasma of the quiet Sun.     

Three Super Active Regions in the Descending Phase of Solar Cycle 23

We analyze the magnetic configurations of three super active regions, NOAA 10484, 10486 and 10488, observed by the Huairou Multi-Channel Solar Telescope (MCST) from 2003 October 18 to November 4. Many energetic phenomena, such as flares (including a X-28 flare) and coronal mass ejections (CMEs), occurred during this period. We think that strong shear and fast emergence of magnetic flux are the main causes of these events. The question is also of great interest why these dramatic eruptions occurred so close together in the descending phase of the solar cycle.     

Fractal properties of solar magnetic fields

We study the spatial properties of solar magnetic fields using data from the Solar Vector Magnetograph of the Marshall Space Flight Center (MSFC) (FeI 5250.2 Å) and SOHO/MDI longitudinal magnetic field measurements (Ni 6767.8 Å) (96-min full-disk maps). Our study is focused on two objects: the fractal properties of sunspots and the fractal properties of the spatial magnetic field distribution of active and quiet regions considered as global structures. To study the spatial structure of sunspots, we use a well-known method of determining the fractal dimension based on an analysis of the perimeter—area relation. To analyze the fractal properties of the spatial magnetic field distribution over the solar surface, we use a technique developed by Higuchi. We have revealed the existence of three families of self-similar contours corresponding to the sunspot umbra, penumbra, and adjacent photosphere. The fractal coefficient has maxima near the umbra—penumbra and penumbra—photosphere boundaries. The fractal dependences of the longitudinal and transverse magnetic field distributions are similar, but the fractal numbers themselves for the transverse fields are larger than those for the longitudinal fields approximately by a factor of 1.5. The fractal numbers decrease with increasing mean magnetic field strength, implying that the magnetic field distribution is more regular in active regions.     

Night‐time science with large solar telescopes: The magnetic Sun through time

Today the Sun has a regular magnetic cycle driven by a dynamo action. But how did this regular cycle develop? How do basic parameters such as rotation rate, age, and differential rotation affect the generation of magnetic fields? Zeeman Doppler imaging (ZDI) is a technique that uses high‐resolution observations in circularly polarised light to map the surface magnetic topology on stars. Utilising the spectropolarimetric capabilities of future large solar telescopes it will be possible to study the evolution and morphology of the magnetic fields on a range of Sun‐like stars from solar twins through to rapidly‐rotating active young Suns and thus study the solar magnetic dynamo through time. In this article I discuss recent results from ZDI of Sun‐like stars and how we can use night‐time observations from future solar telescopes to solve unanswered questions about the origin and evolution of the Sun's magnetic dynamo (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)