The polarity of the toroidal magnetic field in the solar convection zone periodically reverses in the course of the 11/22-year solar cycle. Among the various processes that contribute to the removal of `old-polarity' toroidal magnetic flux is the emergence of flux loops forming bipolar regions at the solar surface. We quantify the loss of subsurface net toroidal flux by this process. To this end, we determine the contribution of an individual emerging bipolar loop and show that it is unaffected by surface flux transport after emergence. Together with the linearity of the diffusion process this means that the total flux loss can be obtained by adding the contributions of all emerging bipolar magnetic regions. The resulting total loss rate of net toroidal flux amounts to 1.3 X 10^{15} Mx/s during activity maxima and 6.1 X 10^{14} Mx/s during activity minima, to which ephemeral regions contribute about 90% and 97%, respectively. This rate is consistent with the observationally inferred loss rate of toroidal flux into interplanetary space and corresponds to a decay time of the subsurface toroidal flux of about 12 years, also consistent with a simple estimate based on turbulent diffusivity. Consequently, toroidal flux loss by flux emergence is a relevant contribution to the budget of net toroidal flux in the solar convection zone. That the toroidal flux loss rate due to flux emergence is consistent with what is expected from turbulent diffusion, and that the corresponding decay time is similar to the length of the solar cycle are important constraints for understanding the solar cycle and the Sun's internal dynamics.

Power spectra of solar activity based on historical records of sunspot numbers and on cosmogenic isotopes show peaks with enhanced power apart from the dominant 11-year solar cycle, such as the 90-year Gleissberg cycle or the 210-year de Vries cycle. In a previous paper we have shown that the overall shape of the power spectrum is well represented by the results of the generic normal form model for a noisy and weakly nonlinear limit cycle, with parameters all determined by observations. Using this model as a null case, we show here that all local peaks with enhanced power, apart from the 11-year band, are consistent with realisation noise. Even a $3\sigma$ peak is expected to occur with a probability of about 0.25 at least once among the 216 period bins resolved by the cosmogenic isotope data. This casts doubt upon interpretations of such peaks in terms of intrinsic periodicities of the solar dynamo process.

In 1969 Leighton developed a quasi-1D mathematical model of the solar dynamo, building upon the phenomenological scenario of Babcock(1961). Here we present a modification and extension of Leighton's model. Using the axisymmetric component of the magnetic field, we consider the radial field component at the solar surface and the radially integrated toroidal magnetic flux in the convection zone, both as functions of latitude. No assumptions are made with regard to the radial location of the toroidal flux. The model includes the effects of turbulent diffusion at the surface and in the convection zone, poleward meridional flow at the surface and an equatorward return flow affecting the toroidal flux, latitudinal differential rotation and the near-surface layer of radial rotational shear, downward convective pumping of magnetic flux in the shear layer, and flux emergence in the form of tilted bipolar magnetic regions. While the parameters relevant for the transport of the surface field are taken from observations, the model condenses the unknown properties of magnetic field and flow in the convection zone into a few free parameters (turbulent diffusivity, effective return flow, amplitude of the source term, and a parameter describing the effective radial shear). Comparison with the results of two-dimensional flux transport dynamo codes shows that the model captures the essential features of these simulations. We carry out a parameter study over the four-dimensional parameter space and identify the parameter ranges that provide solar-like solutions. Dipole parity is always preferred and solutions with periods around 22 years and a correct phase difference between flux emergence in low latitudes and the strength of the polar fields are found for a return flow speed around 2~m/s, turbulent diffusivity <80~km^2/s, and dynamo excitation not too far above the threshold (linear growth rate <0.1/yr).

Surface flux transport simulations for the descending phase of cycle 24 using random sources (emerging bipolar magnetic regions) with empirically determined scatter of their properties provide a prediction of the axial dipole moment during the upcoming activity minimum together with a realistic uncertainty range. The expectation value for the dipole moment around 2020 $(2.5\pm1.1\,$G) is comparable to that observed at the end of cycle 23 (about $2\,$G). The empirical correlation between the dipole moment during solar minimum and the strength of the subsequent cycle thus suggests that cycle 25 will be of moderate amplitude, not much higher than that of the current cycle. However, the intrinsic uncertainty of such predictions resulting from the random scatter of the source properties is considerable and fundamentally limits the reliability with which such predictions can be made before activity minimum is reached.

In order to match observed properties of the solar cycle, flux-transport dynamo models require the toroidal magnetic flux to be stored in a region of low magnetic diffusivity, typically located at or below the bottom of the convection zone. We infer the turbulent magnetic diffusivity affecting the toroidal field on the basis of empirical data. We consider the time evolution of mean latitude and width of the activity belts of solar cycles 12--23 and their dependence on cycle strength. We interpret the decline phase of the cycles as a diffusion process. The activity level of a given cycle begins to decline when the centers of its equatorward propagating activity belts come within their width (at half maximum) from the equator. This happens earlier for stronger cycles because their activity belts are wider. From that moment on, the activity and the belt width decrease in the same manner for all cycles, independent of their maximum activity level. In terms of diffusive cancellation of opposite-polarity toroidal flux across the equator, we infer the turbulent diffusivity experienced by the toroidal field, wherever it is located, to be in the range 150--450~km$^2$s$^{-1}$. Strong diffusive latitudinal spreading of the toroidal flux underneath the activity belts can be inhibited by an inflow towards the toroidal field bands in the convection zone with a magnitude of several meters per second. The inferred value of the turbulent magnetic diffusivity affecting the toroidal field agrees, to order of magnitude, with estimates based on mixing-length models for the solar convection zone. This is at variance with the requirement of flux-transport dynamo models. The inflows required to keep the toroidal field bands together before they approach the equator are similar to the inflows towards the activity belts observed with local helioseismology.

The ongoing 11-year cycle of solar activity is considerably less vigorous than the three cycles before. It was preceded by a very deep activity minimum with a low polar magnetic flux, the source of the toroidal field responsible for solar magnetic activity in the subsequent cycle. Simulation of the evolution of the solar surface field shows that the weak polar fields and thus the weakness of the present cycle 24 are mainly caused by a number of bigger bipolar regions emerging at low latitudes with a `wrong' (i.e., opposite to the majority for this cycle) orientation of their magnetic polarities in the North-South direction, which impaired the growth of the polar field. These regions had a particularly strong effect since they emerged within $\pm10^\circ$ latitude from the solar equator.

We discuss recent observational, theoretical and modeling progress made in understanding the Sun's internal dynamics, including its rotation, meridional flow, convection and overshoot. Over the past few decades, substantial theoretical and observational effort has gone into appreciating these aspects of solar dynamics. A review of these observations, related helioseismic methodology and inference and computational results in relation to these problems is undertaken here.

Magnetic activity in Sun-like and low-mass stars causes X-ray coronal emission, which is stronger for more rapidly rotating stars. This relation is often interpreted in terms of the Rossby number, i.e., the ratio of rotation period to convective overturn time. We reconsider this interpretation on the basis of the observed X-ray emission and rotation periods of 821 stars with masses below 1.4 Msun. A generalized analysis of the relation between X-ray luminosity normalized by bolometric luminosity, L_X/L_bol, and combinations of rotational period, P, and stellar radius, R, shows that the Rossby formulation does not provide the solution with minimal scatter. Instead, we find that the relation L_X/L_bol ~ P^{-2}R^{-4} optimally describes the non-saturated fraction of the stars. This relation is equivalent to L_X ~ P^{-2}, indicating that the rotation period alone determines the total X-ray emission. Since L_X is directly related to the magnetic flux at the stellar surface, this means that the surface flux is determined solely by the star's rotation and is independent of other stellar parameters. While a formulation in terms of a Rossby number would be consistent with these results if the convective overturn time scales exactly as L_bol^{-1/2}, our generalized approach emphasizes the need to test a broader range of mechanisms for dynamo action in cool stars.

The tilt angles of sunspot groups represent the poloidal field source in Babcock-Leighton-type models of the solar dynamo and are crucial for the build-up and reversals of the polar fields in Surface Flux Transport (SFT) simulations. The evolution of the polar field is a consequence of Hale's polarity rules, together with the tilt angle distribution which has a systematic component (Joy's law) and a random component (tilt-angle scatter). We determine the scatter using the observed tilt angle data and study the effects of this scatter on the evolution of the solar surface field using SFT simulations with flux input based upon the recorded sunspot groups. The tilt angle scatter is described in our simulations by a random component according to the observed distributions for different ranges of sunspot group size (total umbral area). By performing simulations with a number of different realizations of the scatter we study the effect of the tilt angle scatter on the global magnetic field, especially on the evolution of the axial dipole moment. The average axial dipole moment at the end of cycle 17 (a medium-amplitude cycle) from our simulations was 2.73G. The tilt angle scatter leads to an uncertainty of 0.78 G (standard deviation). We also considered cycle 14 (a weak cycle) and cycle 19 (a strong cycle) and show that the standard deviation of the axial dipole moment is similar for all three cycles. The uncertainty mainly results from the big sunspot groups which emerge near the equator. In the framework of Babcock-Leighton dynamo models, the tilt angle scatter therefore constitutes a significant random factor in the cycle-to-cycle amplitude variability, which strongly limits the predictability of solar activity.

It is proposed that the observed near-surface inflows towards the active regions and sunspot zones provide a nonlinear feedback mechanism that limits the amplitude of a Babcock-Leighton-type solar dynamo and determines the variation of the cycle strength. This hypothesis is tested with surface flux transport simulations including converging latitudinal flows that depend on the surface distribution of magnetic flux. The inflows modulate the build-up of polar fields (represented by the axial dipole) by reducing the tilt angles of bipolar magnetic regions and by affecting the cross-equator transport of leading-polarity magnetic flux. With flux input derived from the observed record of sunspot groups, the simulations cover the period between 1874 and 1980 (corresponding to solar cycles 11 to 20). The inclusion of the inflows leads to a strong correlation of the simulated axial dipole strength during activity minimum with the observed amplitude of the subsequent cycle. This in agreement with empirical correlations and in line with what is expected from a Babcock-Leighton-type dynamo. The results provide evidence that the latitudinal inflows are a key ingredient in determining the amplitude of solar cycles.

Magnetic flux is continuously being removed and replenished on the solar surface. To understand the removal process we carried out 3D radiative MHD simulations of the evolution of patches of photospheric magnetic field with equal amounts of positive and negative flux. We find that the flux is removed at a rate corresponding to an effective turbulent diffusivity, of 100--340 km^2/s, depending on the boundary conditions. For average unsigned flux densities above about 70 Gauss, the percentage of surface magnetic energy coming from different field strengths is almost invariant. The overall process is then one where magnetic elements are advected by the horizontal granular motions and occasionally come into contact with opposite-polarity elements. These reconnect above the photosphere on a comparatively short time scale after which the U loops produced rapidly escape through the upper surface while the downward retraction of inverse-U loops is significantly slower, because of the higher inertia and lower plasma beta in the deeper layers.

Numerical simulations indicate that the filamentation of sunspot penumbrae and the associated systematic outflow (the Evershed effect) are due to convectively driven fluid motions constrained by the inclined magnetic field. We investigate whether these motions, in particular the upflows in the bright filaments and the downflows at their edges can be reliably observed with existing instrumentation. We use a snapshot from a sunspot simulation to calculate 2D maps of synthetic line profiles for the spectral lines Fe\sci 7090.4 \AA ~ and C\sci 5380.34 \AA. The maps are spatially and spectrally degraded according to typical instrument properties. Line-of-sight velocities are determined from line bisector shifts. We find that the detectability of the convective flows is strongly affected by spatial smearing, particularly so for the downflows. Furthermore, the line-of-sight velocities are dominated by the Evershed flow unless the observation is made very near to disk center. These problems may have compromised recent attempts to detect overturning penumbral convection. Lines with a low formation height are best suited to detect the convective flows.

To reproduce the weak magnetic field on the polar caps of the Sun observed during the declining phase of cycle 23 poses a challenge to surface flux transport models since this cycle has not been particularly weak. We use a well-calibrated model to evaluate the parameter changes required to obtain simulated polar fields and open flux that are consistent with the observations. We find that the low polar field of cycle 23 could be reproduced by an increase of the meridional flow by 55% in the last cycle. Alternatively, a decrease of the mean tilt angle of sunspot groups by 28% would also lead to a similarly low polar field, but cause a delay of the polar field reversals by 1.5 years in comparison to the observations.

We use the historic record of sunspot groups compiled by the Royal Greenwich Observatory together with the sunspot number to derive the statistical properties of sunspot group emergence in dependence of cycle phase and strength. In particular we discuss the latitude, longitude, area and tilt angle of sunspot groups as functions of the cycle strength and of time during the solar cycle. Using these empirical characteristics the time-latitude diagram of sunspot group emergence (butterfly diagram) is reconstructed from 1700 onward on the basis of the Wolf and group sunspot numbers. This reconstruction will be useful in studies of the long-term evolution of the Sun's magnetic field.

We have used semi-synthetic records of emerging sunspot groups based on sunspot number data as input for a surface flux transport model to reconstruct the evolution of the large-scale solar magnetic field and the open heliospheric flux from the year 1700 onward. The statistical properties of the semi-synthetic sunspot group records reflect those of the observed the Royal Greenwich Observatory photoheliographic results. These include correlations between the sunspot numbers and sunspot group latitudes, longitudes, areas and tilt angles. The reconstruction results for the total surface flux, the polar field, and the heliospheric open flux (determined by a current sheet source surface extrapolation) agree well with the available observational or empirically derived data and reconstructions. We confirm a significant positive correlation between the polar field during activity minimum periods and the strength of the subsequent sunspot cycle, which has implications for flux transport dynamo models for the solar cycle. Just prior to the Dalton minimum, at the end of the 18th century, a long cycle was followed by a weak cycle. We find that introducing a possibly `lost' cycle between 1793 and 1800 leads to a shift of the minimum of the open flux by 15 years which is inconsistent with the cosmogenic isotope record.

Changes in the solar surface area covered by small-scale magnetic elements are thought to cause long-term changes in the solar spectral irradiance, which are important for determining the impact on Earth's climate. To study the effect of small-scale magnetic elements on total and spectral irradiance, we derive their contrasts from 3-D MHD simulations of the solar atmosphere. Such calculations are necessary since measurements of small-scale flux tube contrasts are confined to a few wavelengths and suffer from scattered light and instrument defocus, even for space observations. To test the contrast calculations, we compare rms contrasts from simulations with those obtained with the broad-band filter imager mounted on the Solar Optical Telescope (SOT) onboard the Hinode satellite and also analyse centre-to-limb variations (CLV). The 3-D MHD simulations include the interaction between convection and magnetic flux tubes. They have been run with non-grey radiative transfer using the MURaM code. Simulations have an average vertical magnetic field of 0G, 50G, and 200G. Emergent intensities are calculated with the spectral synthesis code ATLAS9 and are convolved with a theoretical point-spread function to account for the properties of the observations' optical system. We find reasonable agreement for simulated and observed intensity distributions in the visible continuum bands. Agreement is poorer for the CN and G-Bands. The analysis of the simulations exhibits a potentially more realistic centre-to-limb behaviour than calculations based on 1-D model atmospheres. We conclude that starting from 3-D MHD simulations represents a powerful approach to obtaining intensity contrasts for a wide wavelength coverage and for different positions on the solar disk. This also paves the way for future calculations of facular and network contrast as a function of magnetic fluxes.

Until today, the small size of magnetic elements in quiet Sun areas has required the application of indirect methods, such as the line-ratio technique or multi-component inversions, to infer their physical properties. A consistent match to the observed Stokes profiles could only be obtained by introducing a magnetic filling factor that specifies the fraction of the observed pixel filled with magnetic field. Here, we investigate the properties of a small magnetic patch in the quiet Sun observed with the IMaX magnetograph on board the balloon-borne telescope Sunrise with unprecedented spatial resolution and low instrumental stray light. We apply an inversion technique based on the numerical solution of the radiative transfer equation to retrieve the temperature stratification and the field strength in the magnetic patch. The observations can be well reproduced with a one-component, fully magnetized atmosphere with a field strength exceeding 1 kG and a significantly enhanced temperature in the mid- to upper photosphere with respect to its surroundings, consistent with semi-empirical flux tube models for plage regions. We therefore conclude that, within the framework of a simple atmospheric model, the IMaX measurements resolve the observed quiet-Sun flux tube.

We present the first observations of the transverse component of photospheric magnetic field acquired by the imaging magnetograph Sunrise/IMaX. Using an automated detection method, we obtain statistical properties of 4536 features with significant linear polarization signal. Their rate of occurrence is 1-2 orders of magnitude larger than values reported by previous studies. We show that these features have no characteristic size or lifetime. They appear preferentially at granule boundaries with most of them being caught in downflow lanes at some point in their evolution. Only a small percentage are entirely and constantly embedded in upflows (16%) or downflows (8%).

We studied the effect of the perturbation of the meridional flow in the activity belts detected by local helioseismology on the development and strength of the surface magnetic field at the polar caps. We carried out simulations of synthetic solar cycles with a flux transport model, which follows the cyclic evolution of the surface field determined by flux emergence and advective transport by near-surface flows. In each hemisphere, an axisymmetric band of latitudinal flows converging towards the central latitude of the activity belt was superposed onto the background poleward meridional flow. The overall effect of the flow perturbation is to reduce the latitude separation of the magnetic polarities of a bipolar magnetic region and thus diminish its contribution to the polar field. As a result, the polar field maximum reached around cycle activity minimum is weakened by the presence of the meridional flow perturbation. For a flow perturbation consistent with helioseismic observations, the polar field is reduced by about 18% compared to the case without inflows. If the amplitude of the flow perturbation depends on the cycle strength, its effect on the polar field provides a nonlinearity that could contribute to limiting the amplitude of a Babcock-Leighton type dynamo.

We demonstrate that a magneto-convection simulation incorporating essential physical processes governing solar surface convection exhibits turbulent small-scale dynamo action. By presenting a derivation of the energy balance equation and transfer functions for compressible magnetohydrodynamics (MHD), we quantify the source of magnetic energy on a scale-by-scale basis. We rule out the two alternative mechanisms for the generation of small-scale magnetic field in the simulations: the tangling of magnetic field lines associated with the turbulent cascade and Alfvenization of small-scale velocity fluctuations ("turbulent induction"). Instead, we find the dominant source of small-scale magnetic energy is stretching by inertial-range fluid motions of small-scale magnetic field lines against the magnetic tension force to produce (against Ohmic dissipation) more small-scale magnetic field. The scales involved become smaller with increasing Reynolds number, which identifies the dynamo as a small-scale turbulent dynamo.

The existence of a turbulent small-scale solar surface dynamo is likely, considering existing numerical and laboratory experiments, as well as comparisons of a small-scale dynamo in MURaM simulations with Hinode observations. We find the observed peaked probability distribution function (PDF) from Stokes-V magnetograms is consistent with a monotonic PDF of the actual vertical field strength. The cancellation function of the vertical flux density from a Hinode SP observation is found to follow a self-similar power law over two decades in length scales down to the ~200 km resolution limit. This provides observational evidence that the scales of magnetic structuring in the photosphere extend at least down to 20 km. From the power law, we determine a lower bound for the true quiet-Sun mean vertical unsigned flux density of ~43 G, consistent with our numerically-based estimates that 80% or more of the vertical unsigned flux should be invisible to Stokes-V observations at a resolution of 200 km owing to cancellation. Our estimates significantly reduce the order-of-magnitude discrepancy between Zeeman- and Hanle-based estimates.

We obtain information about the magnetic flux present in the quiet Sun by comparing radiative MHD simulations with Hinode/SP observations, with particular emphasis on the role of surface dynamo action. Simulation runs with different magnetic Reynolds numbers (Rm) are used together with observations at different heliocentric angles with different levels of noise. The results show that simulations with an imposed mixed-polarity field and Rm below the threshold for dynamo action reproduce the observed vertical flux density, but do not display a sufficiently high horizontal flux density. Surface dynamo simulations at the highest Rm feasible at the moment yield a ratio of the horizontal and vertical flux density consistent with observational results, but the overall amplitudes are too low. Based on the properties of the local dynamo simulations, a tentative scaling of the magnetic field strength by a factor 2 - 3 reproduces the signal observed in the internetwork regions. We find an agreement with observations at different heliocentric angles. The mean field strength in internetwork, implied by our analysis, is roughly 170 G at the optical depth unity. Our study shows that surface dynamo could be responsible for most of the magnetic flux in the quiet Sun outside the network given that the extrapolation to higher Rm is valid.

By coupling a solar surface flux transport model with an extrapolation of the heliospheric field, we simulate the evolution of the Sun's open magnetic flux and the heliospheric current sheet (HCS) based on observational data of sunspot groups since 1976. The results are consistent with measurements of the interplanetary magnetic field near Earth and with the tilt angle of the HCS as derived from extrapolation of the observed solar surface field. This opens the possibility for an improved reconstruction of the Sun's open flux and the HCS into the past on the basis of empirical sunspot data.

Realistic 3D radiative MHD simulations reveal the magneto-convective processes underlying the formation of the photospheric fine structure of sunspots, including penumbral filaments and umbral dots. Here we provide results from a statistical analysis of simulated umbral dots and compare them with reports from high-resolution observations. A multi-level segmentation and tracking algorithm has been used to isolate the bright structures in synthetic bolometric and continuum brightness images. Areas, brightness, and lifetimes of the resulting set of umbral dots are found to be correlated: larger umbral dots tend to be brighter and live longer. The magnetic field strength and velocity structure of umbral dots on surfaces of constant optical depth in the continuum at 630 nm indicate that the strong field reduction and high velocities in the upper parts of the upflow plumes underlying umbral dots are largely hidden from spectro-polarimetric observations. The properties of the simulated umbral dots are generally consistent with the results of recent high-resolution observations. However, the observed population of small, short-lived umbral dots is not reproduced by the simulations, possibly owing to insufficient spatial resolution.

We study the observational signature of flux emergence in the photosphere using synthetic data from a 3D MHD simulation of the emergence of a twisted flux tube. Several stages in the emergence process are considered. At every stage we compute synthetic Stokes spectra of the two iron lines Fe I 6301.5 {\AA} and Fe I 6302.5 {\AA} and degrade the data to the spatial and spectral resolution of Hinode's SOT/SP. Then, following observational practice, we apply Milne-Eddington-type inversions to the synthetic spectra in order to retrieve various atmospheric parameters and compare the results with recent Hinode observations. During the emergence sequence, the spectral lines sample different parts of the rising flux tube, revealing its twisted structure. The horizontal component of the magnetic field retrieved from the simulations is close to the observed values. The flattening of the flux tube in the photosphere is caused by radiative cooling, which slows down the ascent of the tube to the upper solar atmosphere. Consistent with the observations, the rising magnetized plasma produces a blue shift of the spectral lines during a large part of the emergence sequence.

The structure and dynamics of small vertical photospheric magnetic flux concentrations has been often treated in the framework of an approximation based upon a low-order truncation of the Taylor expansions of all quantities in the horizontal direction, together with the assumption of instantaneous total pressure balance at the boundary to the non-magnetic external medium. Formally, such an approximation is justified if the diameter of the structure (a flux tube or a flux sheet) is small compared to all other relevant length scales (scale height, radius of curvature, wavelength, etc.). The advent of realistic 3D radiative MHD simulations opens the possibility of checking the consistency of the approximation with the properties of the flux concentrations that form in the course of a simulation. We carry out a comparative analysis between the thin flux tube/sheet models and flux concentrations formed in a 3D radiation-MHD simulation. We compare the distribution of the vertical and horizontal components of the magnetic field in a 3D MHD simulation with the field distribution in the case of the thin flux tube/sheet approximation. We also consider the total (gas plus magnetic) pressure in the MHD simulation box. Flux concentrations with super-equipartition fields are reasonably well reproduced by the second-order thin flux tube/sheet approximation. The differences between approximation and simulation are due to the asymmetry and the dynamics of the simulated structures.

Long-term trends in the solar spectral irradiance are important to determine the impact on Earth's climate. These long-term changes are thought to be caused mainly by changes in the surface area covered by small-scale magnetic elements. The direct measurement of the contrast to determine the impact of these small-scale magnetic elements is, however, limited to a few wavelengths, and is, even for space instruments, affected by scattered light and instrument defocus. In this work we calculate emergent intensities from 3-D simulations of solar magneto-convection and validate the outcome by comparing with observations from Hinode/SOT. In this manner we aim to construct the contrast at wavelengths ranging from the NUV to the FIR.

Using turbulent MHD simulations (magnetic Reynolds numbers up to 8000) and Hinode observations, we study effects of turbulence on measuring the solar magnetic field outside active regions. Firstly, from synthetic Stokes V profiles for the FeI lines at 630.1 and 630.2 nm, we show that a peaked probability distribution function (PDF) for observationally-derived field estimates is consistent with a monotonic PDF for actual vertical field strengths. Hence, the prevalence of weak fields is greater than would be naively inferred from observations. Secondly, we employ the fractal self-similar geometry of the turbulent solar magnetic field to derive two estimates (numerical and observational) of the true mean vertical unsigned flux density. We also find observational evidence that the scales of magnetic structuring in the photosphere extend at least down to an order of magnitude smaller than 200 km: the self-similar power-law scaling in the signed measure from a Hinode magnetogram ranges (over two decades in length scales and including the granulation scale) down to the 200 km resolution limit. From the self-similar scaling, we determine a lower bound for the true quiet-Sun mean vertical unsigned flux density of ~50 G. This is consistent with our numerically-based estimates that 80% or more of the vertical unsigned flux should be invisible to Stokes-V observations at a resolution of 200 km owing to the cancellation of signal from opposite magnetic polarities. Our estimates significantly reduce the order-of-magnitude discrepancy between Zeeman- and Hanle-based estimates.

We present results from numerical modeling of emerging flux regions on the solar surface. The modeling was carried out by means of 3D radiative MHD simulations of the rise of buoyant magnetic flux tubes through the convection zone and into the photosphere. Due to the strong stratification of the convection zone, the rise results in a lateral expansion of the tube into a magnetic sheet, which acts as a reservoir for small-scale flux emergence events at the scale of granulation. The interaction of the convective downflows and the rising magnetic flux undulates it to form serpentine field lines emerging into the photosphere. Observational characteristics including the pattern of emerging flux regions, the cancellation of surface flux and associated high speed downflows, the convective collapse of photospheric flux tubes, the appearance of anomalous darkenings, the formation of bright points and the possible existence of transient kilogauss horizontal fields are discussed in the context of new observations from the Hinode Solar Optical Telescope. Implications for the local helioseismology of emerging flux regions are also discussed.

Results of a 3D MHD simulation of a sunspot with a photospheric size of about 20 Mm are presented. The simulation has been carried out with the MURaM code, which includes a realistic equation of state with partial ionization and radiative transfer along many ray directions. The largely relaxed state of the sunspot shows a division in a central dark umbral region with bright dots and a penumbra showing bright filaments of about 2 to 3 Mm length with central dark lanes. By a process similar to the formation of umbral dots, the penumbral filaments result from magneto-convection in the form of upflow plumes, which become elongated by the presence of an inclined magnetic field: the upflow is deflected in the outward direction while the magnetic field is weakened and becomes almost horizontal in the upper part of the plume near the level of optical depth unity. A dark lane forms owing to the piling up of matter near the cusp-shaped top of the rising plume that leads to an upward bulging of the surfaces of constant optical depth. The simulated penumbral structure corresponds well to the observationally inferred interlocking-comb structure of the magnetic field with Evershed outflows along dark-laned filaments with nearly horizontal magnetic field and overturning perpendicular (`twisting') motion, which are embedded in a background of stronger and less inclined field. Photospheric spectral lines are formed at the very top and somewhat above the upflow plumes, so that they do not fully sense the strong flow as well as the large field inclination and significant field strength reduction in the upper part of the plume structures.

Observations with the Hinode spectro-polarimeter have revealed strong horizontal internetwork magnetic fields in the quiet solar photosphere. We aim at interpreting the observations by means of results from numerical simulations. Radiative MHD simulations of dynamo action by near-surface convection are analyzed with respect to the relation between vertical and horizontal magnetic field components. The dynamo-generated fields show a clear dominance of the horizontal field in the height range where the spectral lines used for the observations are formed. The ratio between the averaged horizontal and vertical field components is consistent with the values derived from the observations. This behavior results from the intermittent nature of the dynamo field with polarity mixing on small scales in the surface layers. Our results provide further evidence that local near-surface dynamo action contributes significantly to the solar internetwork fields.

Various methods (or recipes) have been proposed to predict future solar activity levels - with mixed success. Among these, some precursor methods based upon quantities determined around or a few years before solar minimum have provided rather high correlations with the strength of the following cycles. Recently, data assimilation with an advection-dominated (flux-transport) dynamo model has been proposed as a predictive tool, yielding remarkably high correlation coefficients. After discussing the potential implications of these results and the criticism that has been raised, we study the possible physical origin(s) of the predictive skill provided by precursor and other methods. It is found that the combination of the overlap of solar cycles and their amplitude-dependent rise time (Waldmeier's rule) introduces correlations in the sunspot number (or area) record, which account for the predictive skill of many precursor methods. This explanation requires no direct physical relation between the precursor quantity and the dynamo mechanism (in the sense of the Babcock-Leighton scheme or otherwise).

Flow-induced instabilities are relevant for the storage and dynamics of magnetic fields in stellar convection zones and possibly also in other astrophysical contexts. We continue the study started in the first paper of this series by considering the stability properties of longitudinal flows along magnetic flux tubes. A linear stability analysis was carried out to determine criteria for the onset of instability in the framework of the approximation of thin magnetic flux tubes. In the non-dissipative case, we find Kelvin-Helmholtz instability for flow velocities exceeding a critical speed that depends on the Alfv{\'e}n speed and on the ratio of the internal and external densities. Inclusion of a friction term proportional to the relative transverse velocity leads to a friction-driven instability connected with backward (or negative energy) waves. We discuss the physical nature of this instability. In the case of a stratified external medium, the Kelvin-Helmholtz instability and the friction-driven instability can set in for flow speeds significantly lower than the Alfv{\'e}n speed. Dissipative effects can excite flow-driven instability below the thresholds for the Kelvin-Helmholtz and the undulatory (Parker-type) instabilities. This may be important for magnetic flux storage in stellar convection zones and for the stability of astrophysical jets.

We present results of synthetic spectro-polarimetric diagnostics of radiative MHD simulations of solar surface convection with magnetic fields. Stokes profiles of Zeeman-sensitive lines of neutral iron in the visible and infrared spectral ranges emerging from the simulated atmosphere have been calculated in order to study their relation to the relevant physical quantities and compare with observational results. We have analyzed the dependence of the Stokes-I line strength and width as well as of the Stokes-V signal and asymmetries on the magnetic field strength. Furthermore, we have evaluated the correspondence between the actual velocities in the simulation with values determined from the Stokes-I (Doppler shift of the centre of gravity) and Stokes-V profiles (zero-crossing shift). We confirm that the line weakening in strong magnetic fields results from a higher temperature (at equal optical depth) in the magnetic flux concentrations. We also confirm that considerable Stokes-V asymmetries originate in the peripheral parts of strong magnetic flux concentrations, where the line of sight cuts through the magnetopause of the expanding flux concentration into the surrounding convective donwflow.

Context: Observations indicate that the `quiet' solar photosphere outside active regions contains considerable amounts of magnetic energy and magnetic flux, with mixed polarity on small scales. The origin of this flux is unclear. Aims: We test whether local dynamo action of the near-surface convection (granulation) can generate a significant contribution to the observed magnetic flux. Methods: We have carried out MHD simulations of solar surface convection, including the effects of strong stratification, compressibility, partial ionization, radiative transfer, as well as an open lower boundary. Results: Exponential growth of a weak magnetic seed field (with vanishing net flux through the computational box) is found in a simulation run with a magnetic Reynolds number of about 2600. The magnetic energy approaches saturation at a level of a few percent of the total kinetic energy of the convective motions. Near the visible solar surface, the (unsigned) magnetic flux density reaches at least a value of about 25 G. Conclusions: A realistic flow topology of stratified, compressible, non-helical surface convection without enforced recirculation is capable of turbulent local dynamo action near the solar surface.

We study the emergence of magnetic flux from the near-surface layers of the solar convection zone into the photosphere. To model magnetic flux emergence, we carried out a set of numerical radiative magnetohydrodynamics simulations. Our simulations take into account the effects of compressibility, energy exchange via radiative transfer, and partial ionization in the equation of state. All these physical ingredients are essential for a proper treatment of the problem. Furthermore, the inclusion of radiative transfer allows us to directly compare the simulation results with actual observations of emerging flux. We find that the interaction between the magnetic flux tube and the external flow field has an important influence on the emergent morphology of the magnetic field. Depending on the initial properties of the flux tube (e.g. field strength, twist, entropy etc.), the emergence process can also modify the local granulation pattern. The emergence of magnetic flux tubes with a flux of $10^{19}$ Mx disturbs the granulation and leads to the transient appearance of a dark lane, which is coincident with upflowing material. These results are consistent with observed properties of emerging magnetic flux.

We study the origin of the predictive skill of some methods to forecast the strength of solar activity cycles. A simple flux transport model for the azimuthally averaged radial magnetic field at the solar surface is used, which contains a source term describing the emergence of new flux based on observational sunspot data. We consider the magnetic flux diffusing over the equator as a predictor, since this quantity is directly related to the global dipole field from which a Babcock-Leighton dynamo generates the toroidal field for the next activity cycle. If the source is represented schematically by a narrow activity belt drifting with constant speed over a fixed range of latitudes between activity minima, our predictor shows considerable predictive skill with correlation coefficients up to 0.95 for past cycles. However, the predictive skill is completely lost when the actually observed emergence latitudes are used. This result originates from the fact that the precursor amplitude is determined by the sunspot activity a few years before solar minimum. Since stronger cycles tend to rise faster to their maximum activity (known as the Waldmeier effect), the temporal overlapping of cycles leads to a shift of the minimum epochs that depends on the strength of the following cycle. This information is picked up by precursor methods and also by our flux transport model with a schematic source. Therefore, their predictive skill does not require a memory, i.e., a physical connection between the surface manifestations of subsequent activity cycles.

We study the structure and reveal the physical nature of the reversed granulation pattern in the solar photosphere by means of 3-dimensional radiative hydrodynamics simulations. We used the MURaM code to obtain a realistic model of the near-surface layers of the convection zone and the photosphere. The pattern of horizontal temperature fluctuations at the base of the photosphere consists of relatively hot granular cells bounded by the cooler intergranular downflow network. With increasing height in the photosphere, the amplitude of the temperature fluctuations diminishes. At a height of z=130-140 km in the photosphere, the pattern of horizontal temperature fluctuations reverses so that granular regions become relatively cool compared to the intergranular network. Detailed analysis of the trajectories of fluid elements through the photosphere reveal that the motion of the fluid is non-adiabatic, owing to strong radiative cooling when approaching the surface of optical depth unity followed by reheating by the radiation field from below. The temperature structure of the photosphere results from the competition between expansion of rising fluid elements and radiative heating. The former acts to lower the temperature of the fluid whereas the latter acts to increase it towards the radiative equilibrium temperature with a net entropy gain. After the fluid overturns and descends towards the convection zone, radiative energy loss again decreases the entropy of the fluid. Radiative heating and cooling of fluid elements that penetrate into the photosphere and overturn do not occur in equal amounts. The imbalance in the cumulative heating and cooling of these fluid elements is responsible for the reversal of temperature fluctuations with respect to height in the photosphere.

We carry out numerical simulations of the surface evolution of bipolar magnetic regions (BMRs) and magnetic spots on stars, which have radii and surface rotational shears of AB Doradus, the Sun, and the HR 1099 primary. The surface flux transport model is based on the magnetic induction equation for radial fields under the effects of surface differential rotation, meridional flow, and turbulent diffusion due to convective flow patterns. We calculate the flux evolution and the lifetimes of BMRs and unipolar starspots, varying the emergence latitude, surface shear rate, and tilt angle. For BMRs comparable to the largest observed on the Sun, we find that varying the surface flows and the tilt angle modifies the lifetimes over a range of one month. For very large BMRs (area ~10% of the stellar surface) the assumption of a tilt angle increasing with latitude leads to a significant increase of lifetime, as compared to the case without tilt. Such regions can evolve to polar spots that live more than a year. Adopting the observed weak latitudinal shear and the radius of the active subgiant component of HR 1099, we find longer BMR lifetimes as compared to the more strongly sheared AB Dor case. Random emergence of six additional tilted bipoles in an activity belt at 60 degrees latitude enhanced the lifetimes of polar caps up to 7 years. We have also compared the evolution and lifetime of monolithic starspots with those of conglomerates of smaller spots of similar total area. We find similar decay patterns and lifetimes for both configurations.

Results from a realistic simulation of 3D radiative magneto-convection in a strong background magnetic field corresponding to the conditions in sunspot umbrae are shown. The convective energy transport is dominated by narrow upflow plumes with adjacent downflows, which become almost field-free near the surface layers. The strong external magnetic field forces the plumes to assume a cusp-like shape in their top parts, where the upflowing plasma loses its buoyancy. The resulting bright features in intensity images correspond well (in terms of brightness, size, and lifetime) to the observed umbral dots in the central parts of sunspot umbrae. Most of the simulated umbral dots have a horizontally elongated form with a central dark lane. Above the cusp, most plumes show narrow upflow jets, which are driven by the pressure of the piled-up plasma below. The large velocities and low field strengths in the plumes are effectively screened from spectroscopic observation because the surfaces of equal optical depth are locally elevated, so that spectral lines are largely formed above the cusp. Our simulations demonstrate that nearly field-free upflow plumes and umbral dots are a natural result of convection in a strong, initially monolithic magnetic field.

Customary two-dimensional flux transport models for the evolution of the magnetic field at the solar surface do not account for the radial structure and the volume diffusion of the magnetic field. When considering the long-term evolution of magnetic flux, this omission can lead to an unrealistic long-term memory of the system and to the suppression of polar field reversals. In order to avoid such effects, we propose an extension of the flux transport model by a linear decay term derived consistently on the basis of the eigenmodes of the diffusion operator in a spherical shell. A decay rate for each eigenmode of the system is determined and applied to the corresponding surface part of the mode evolved in the flux transport model. The value of the volume diffusivity associated with this decay term can be estimated to be in the range 50--100 km^2/s by considering the reversals of the polar fields in comparison of flux transport simulations with observations. We show that the decay term prohibits a secular drift of the polar field in the case of cycles of varying strength, like those exhibited by the historical sunspot record.

After a dynamically active emergence phase, magnetic flux at the solar surface soon ceases to show strong signs of the subsurface dynamics of its parent magnetic structure. This indicates that some kind of disconnection of the emerged flux from its roots in the deep convection zone should take place. We propose a mechanism for the dynamical disconnection of the surface flux based upon the buoyant upflow of plasma along the field lines. Such flows arise in the upper part of a rising flux loop during the final phases of its buoyant ascent towards the surface. The combination of the pressure buildup by the upflow and the cooling of the upper layers of an emerged flux tube by radiative losses at the surface lead to a progressive weakening of the magnetic field in several Mm depth. When the field strength has become sufficiently low, convective motions and the fluting instability disrupt the flux tube into thin, passively advected flux fragments, thus providing a dynamical disconnection of the emerged part from its roots. We substantiate this scenario by considering the quasi-static evolution of a sunspot model under the effects of radiative cooling, convective energy transport, and pressure buildup by a prescribed inflow at the bottom of the model. For inflow speeds in the range shown by simulations of thin flux tubes, we find that the disconnection takes place in a depth between 2 and 6 Mm for disconnection times up to 3 days.

This paper is a ``Habilitationsschrift'', a second thesis required until recently by universities in Germany and in a few other countries to obtain the right to lecture. It was accepted by the University of Goettingen in 1990 after review by a number of german and international experts. Although the introduction and the references represent the state of research in 1990, most of the remaining content is still relevant and has never been published elsewhere. The most important part is the derivation of a linear stability formalism for thin magnetic flux tubes following an arbitrary path in a gravitationally stratified medium with a stationary velocity. It was found later (Ferriz-Mas & Schuessler, Geophys. Astrophys. Fluid Dyn. vol. 72, 209; 1995) that, for consistency, the inertial term in the equation of motion for the external medium should be included in Eq. (3.24), which leads to an additional term in the stability equations in the case of a spatially varying external velocity. This term is missing in the present text, but can be easily introduced into the formalism. The full abstract can be found at the beginning of the paper.

The spatio-temporal relationship between the sign of the observed radial component of the magnetic field at the solar surface and the sign of the toroidal field as inferred from Hale's polarity rules for sunspots is usually interpreted as signifying the phase relation between the poloidal and the toroidal magnetic field components involved in the solar dynamo process. This has been taken as a constraint for models of the solar dynamo. This note draws attention to the fact that the observed phase relation is naturally and inevitably produced by the emergence of tilted bipolar regions and flux transport through surface flows, without any necessity of recourse to the dynamo process. Consequently, there is no constraint on dynamo models resulting from the observed phase relation.

We address the question whether the properties of the observed latitude-time diagram of sunspot occurence (the butterfly diagram) provide evidence for the operation of a flux-transport dynamo, which explains the migration of the sunspot zones and the period of the solar cycle in terms of a deep equatorward meridional flow. We show that the properties of the butterfly diagram are equally well reproduced by a conventional dynamo model with migrating dynamo waves, but without transport of magnetic flux by a flow. These properties seem to be generic for an oscillatory and migratory field of dipole parity and thus do not permit an observational distinction between different dynamo approaches.

We consider the mechanical equilibrium and stability of matter trapped in the magnetosphere of a rapidly rotating star. Assuming a dipolar magnetic field and arbitrary inclination of the magnetic axis with respect to the axis of rotation we find stable equilibrium positions a) in a (warped) disk roughly aligned with the magnetic equatorial plane and b) at two locations above and below the disk, whose distance from the star increases with decreasing inclination angle between dipole and rotation axis. The distribution of matter is not strongly affected by allowing for a spatial offset of the magnetic dipole. These results provide a possible explanation for some observations of corotating localized mass concentrations in hot magnetic stars.

The extension of the sunspot number series backward in time is of considerable interest for dynamo theory, solar, stellar, and climate research. We have used records of the Be-10 concentration in polar ice to reconstruct the average sunspot activity level for the period between the year 850 to the pr esent. Our method uses physical models for processes connecting the Be-10 concentration with the sunspot number. The reconstruction shows reliably that the period of high solar activity during the last 60 years is unique throughout the past 1150 years. This nearly triples the time interval for which such a statement could be made previously

In a recent paper (Usoskin et al., 2002a), we have reconstructed the concentration of the cosmogenic $^{10}$Be isotope in ice cores from the measured sunspot numbers by using physical models for $^{10}$Be production in the Earth's atmosphere, cosmic ray transport in the heliosphere, and evolution of the Sun's open magnetic flux. Here we take the opposite route: starting from the $^{10}$Be concentration measured in ice cores from Antarctica and Greenland, we invert the models in order to reconstruct the 11-year averaged sunspot numbers since 850 AD. The inversion method is validated by comparing the reconstructed sunspot numbers with the directly observed sunspot record since 1610. The reconstructed sunspot record exhibits a prominent period of about 600 years, in agreement with earlier observations based on cosmogenic isotopes. Also, there is evidence for the century scale Gleissberg cycle and a number of shorter quasi-periodicities whose periods seem to fluctuate in the millennium time scale. This invalidates the earlier extrapolation of multi-harmonic representation of sunspot activity over extended time intervals.

Surface reconstructions of active close binary stars based on photometric and spectroscopic observations reveal non-uniform starspot distributions, which indicate the existence of preferred spot longitudes (with respect to the companion star). We consider the equilibrium and linear stability of toroidal magnetic flux tubes in close binaries to examine whether tidal effects are capable to initiate the formation of rising flux loops at preferred longitudes near the bottom of the stellar convection zone. The tidal force and the deviation of the stellar structure from spherical symmetry are treated in lowest-order perturbation theory assuming synchronised close binaries with orbital periods of a few days. The frequency, growth time, and spatial structure of linear eigenmodes are determined by a stability analysis. We find that, despite their small magnitude, tidal effects can lead to a considerable longitudinal asymmetry in the formation probability of flux loops, since the breaking of the axial symmetry due to the presence of the companion star is reinforced by the sensitive dependence of the stability properties on the stellar stratification and by resonance effects. The orientation of preferred longitudes of loop formation depends on the equilibrium configuration and the wave number of the dominating eigenmode. The change of the growth times of unstable modes with respect to the case of a single star is very small.

Observations of magnetically active close binaries with orbital periods of a few days reveal the existence of starspots at preferred longitudes (with respect to the direction of the companion star). We numerically investigate the non-linear dynamics and evolution of magnetic flux tubes in the convection zoneof a fast-rotating component of a close binary system and explore whether the tidal effects are able to generate non-uniformities in the surface distribution of erupting flux tubes. Assuming a synchronised system with a rotation period of two days and consisting of two solar-type components, both the tidal force and the deviation of the stellar structure from spherical shape are considered in lowest-order perturbation theory. The magnetic field is initially stored in the form of toroidal magnetic flux rings within the stably stratified overshoot region beneath the convection zone. Once the field has grown sufficiently strong, instabilities initiate the formation of rising flux loops, which rise through the convection zone and emerge at the stellar surface. We find that although the magnitude of tidal effects is rather small, they nevertheless lead to the formation of clusters of flux tube eruptions at preferred longitudes on opposite sides of the star, which result from the cumulative and resonant character of the action of tidal effects on rising flux tubes. The longitude distribution of the clusters depends on the initial parameters of flux tubes in the overshoot region like magnetic field strength and latitude, implying that there is no globally unique preferred longitude along a fixed direction.