The first determination of the surface temperature of stars other than the Sun is due to the Hungarian astrophysicist B\'ela Hark\'anyi. Prompted by the recent unprecedented increase in the availability of stellar temperature estimates from Gaia, coinciding with the 150th anniversary of Hark\'anyi's birth, this article presents the life and work of this neglected, yet remarkable figure in the context of the history of stellar astrophysics.

A review of solar cycle prediction methods and their performance is given, including early forecasts for cycle 25. The review focuses on those aspects of the solar cycle prediction problem that have a bearing on dynamo theory. The scope of the review is further restricted to the issue of predicting the amplitude (and optionally the epoch) of an upcoming solar maximum no later than right after the start of the given cycle. In their overall performance during the course of the last few solar cycles, precursor methods have clearly been superior to extrapolation methods. One method that has yielded predictions consistently in the right range during the past few solar cycles is the polar field precursor. Nevertheless, some extrapolation methods may still be worth further study. Model based forecasts are quickly coming into their own, and, despite not having a long proven record, their predictions are received with increasing confidence by the community.

The choice of free parameters in surface flux transport (SFT) models describing the evolution of the large-scale poloidal magnetic field of the Sun is critical for the correct reproduction of the polar magnetic flux built up during a solar cycle, which in turn is known to be a good predictor of the amplitude of the upcoming cycle. For an informed choice of parameters it is important to understand the effect and interplay of the various parameters and to optimize the models for the polar magnetic field. Here we present the results of a large-scale systematic study of the parameter space in an SFT model where the source term representing the net effect of tilted flux emergence was chosen to represent a typical, average solar cycle as described by observations. Comparing the results with observational constraints on the spatiotemporal variation of the polar magnetic field, as seen in magnetograms for the last four solar cycles, we mark allowed and excluded regions in the 3D parameter space defined by the flow amplitude u0, the magnetic diffusivity eta and the decay time scale tau, for three different assumed meridional flow profiles. Without a significant decay term in the SFT equation (i.e., for tau >10 yr) the global dipole moment reverses too late in the cycle for all flow profiles and parameters, providing independent supporting evidence for the need of a decay term, even in the case of identical cycles. An allowed domain is found to exist for tau values in the 5-10 yr range for all flow profiles considered. Generally higher values of eta (500-800 km^2/s) are preferred though some solutions with lower eta are still allowed.

New developments in surface flux transport modeling and theory of flux transport dynamos have given rise to the notion that certain large active regions with anomalous properties (location, tilt angle and/or Hale/non-Hale character) may have a major impact on the course of solar activity in subsequent years, impacting also on the amplitude of the following solar cycles. Here we discuss our current understanding of the role of such "rogue" active regions in cycle-to-cycle variations of solar activity.

After reviewing potential early indicators of an upcoming solar cycle at high latitudes, we focus attention on the rush-to-the-poles (RTTP) phenomenon in coronal green line emission. Considering various correlations between properties of the RTTP with the upcoming solar cycle we find a correlation between the rate of the RTTP and the time delay until the maximum of the next solar cycle. On the basis of this correlation and the known internal regularities of the sunspot number series we predict that, following a minimum in 2019, cycle 25 will peak in late 2024 at an amplitude of about 130 (in terms of smoothed monthly revised sunspot numbers). This slightly exceeds the amplitude of cycle 24 but it would still make cycle 25 a fairly weak cycle.

The origin of cycle-to-cycle variations in solar activity is currently the focus of much interest. It has recently been pointed out that large individual active regions with atypical properties can have a significant impact on the long term behaviour of solar activity. We investigate this possibility in more detail using a recently developed 2$\times$2D dynamo model of the solar magnetic cycle. We find that even a single "rogue" bipolar magnetic region (BMR) in the simulations can have a major effect on the further development of solar activity cycles, boosting or suppressing the amplitude of subsequent cycles. In extreme cases an individual BMR can completely halt the dynamo, triggering a grand minimum. Rogue BMRs also have the potential to induce significant hemispheric asymmetries in the solar cycle. To study the effect of rogue BMRs in a more systematic manner, a series of dynamo simulations were conducted, in which a large test BMR was manually introduced in the model at various phases of cycles of different amplitudes. BMRs emerging in the rising phase of a cycle can modify the amplitude of the ongoing cycle while BMRs emerging in later phases will only impact subsequent cycles. In this model, the strongest impact on the subsequent cycle occurs when the rogue BMR emerges around cycle maximum at low latitudes but the BMR does not need to be strictly cross-equatorial. Active regions emerging as far as $20^{\circ}$ from the equator can still have a significant impact. We demonstrate that the combined effect of the magnetic flux, tilt angle and polarity separation of the BMR on the dynamo is via their contribution to the dipole moment, $\delta D_{\mathrm{BMR}}$. Our results indicate that prediction of the amplitude, starting epoch and duration of a cycle requires an accurate accounting of a broad range of active regions emerging in the previous cycle.

We study the different patterns of interannual magnetic variability in stars on or near the lower main sequence, approximately solar-type (G-K dwarf) stars in time series of 36 years from the Mount Wilson Observatory Ca\,{\sc ii}\,H\&K survey. Our main aim is to search for correlations between cycles, activity measures and ages. Time-frequency analysis has been used to discern and reveal patterns and morphology of stellar activity cycles, including multiple and changing cycles, in the datasets. Both the results from short-term Fourier transform and its refinement using the Choi-Williams distribution, with better frequency resolution, are presented in this study. Rotational periods of the stars were derived using multi-frequency Fourier analysis.From the studied 29 stars we found at least one activity cycle on 28 of them. Twelve stars, with longer rotational periods ($39.7\pm6.0$ days) have simple, smooth cycles, and the rest of the stars, with on-average much faster rotation ($18.1\pm12.2$ days) show complex and sometimes vigorously changing, multiple cycles. The cycles are longer and quite uniform in the first group ($9.7\pm1.9$ years), while are generally shorter and with greater variety in the second one ($7.6\pm4.9$). There is a clear age division between stars with smooth and complex cycles that follows the known separation between the older and younger stars at around 2 to 3~Gyr of age.

Structures seen in idealized numerical experiments on compressible magnetoconvection in an imposed strong vertical magnetic field show important differences from those detected in observations or realistic numerical simulations of sunspot umbrae. To elucidate the origin of these discrepancies, we present a series of idealized 3D compressible magnetoconvection experiments that differ from previous such experiments in several details, bringing them closer to realistic solar conditions. An initially vertical magnetic field B 0 is imposed on a time snapshot of fully developed solar-like turbulent convection in a layer bounded by a stable layer from above. Upon relaxation to a statistically steady state, the structure of the flow field and magnetic field is examined. Instead of the vigorous granular convection (GRC) well known to take place in magnetized or weakly magnetized convection, for high values of B0 heat is transported by small-scale convection (SSC) in the form of narrow, persistent convective columns consisting of slender upflows accompanied by adjacent downflow patches, which are reminiscent of the 'convectons' identified in earlier semianalytic models. For moderate field strengths, flux separation (FXS) is observed: isolated field-free inclusions of GRC are embedded in a strongly magnetized plasma with SSC. Between the SSC and FXS regimes, a transitional regime (F/S) is identified where convectons dynamically evolve into multiply segmented granular inclusions and back. Our results agree in some aspects more closely with observed umbral structures than earlier idealized models, because they do reproduce the strong localized, patchy downflows immediately adjacent to the narrow convective columns. Based on recent observations of umbral dots, we suggest that in some cases the conditions in sunspot umbrae correspond to the newly identified F/S transitional regime.

This article reviews some of the leading results obtained in solar dynamo physics by using temporal oscillator models as a tool to interpret observational data and dynamo model predictions. We discuss how solar observational data such as the sunspot number is used to infer the leading quantities responsible for the solar variability during the last few centuries. Moreover, we discuss the advantages and difficulties of using inversion methods (or backward methods) over forward methods to interpret the solar dynamo data. We argue that this approach could help us to have a better insight about the leading physical processes responsible for solar dynamo, in a similar manner as helioseismology has helped to achieve a better insight on the thermodynamic structure and flow dynamics in the Sun's interior.

We study the behaviour of the van der Pol oscillator when either its damping parameter $\mu$ or its nonlinearity parameter $\xi$ is subject to additive or multiplicative random noise. Assuming various power law exponents for the relation between the oscillating variable and the sunspot number, for each case we map the parameter plane defined by the amplitude and the correlation time of the perturbation and mark the parameter regime where the sunspot number displays solar-like behaviour. Solar-like behaviour is defined here as a good correlation between the rise rate and cycle amplitude {\it and} the lack of a good correlation between the decay rate and amplitude, together with significant ($\ga 10$\,%) r.m.s. variation in cycle lengths and cycle amplitudes. It is found that perturbing $\mu$ alone the perturbed van der Pol oscillator does not show solar-like behaviour. When the perturbed variable is $\xi$, solar-like behaviour is displayed for perturbations with a correlation time of about 3--4 years and significant amplitude. Such studies may provide useful constraints on solar dynamo models and their parameters.

A simple analytical relation of form {\alpha} = 2 {\kappa} -1 between the magnetic energy spectral exponent {\alpha} of the turbulent magnetic field in the solar photosphere and its magnetic flux cancellation exponent {\kappa}, valid under certain restrictive assumptions, is tested and extended outside its range of validity in a series of Monte Carlo simulations. In these numerical tests artificial "magnetograms" are constructed in 1D and 2D by superposing a discrete set of Fourier modes of the magnetic field distribution with amplitudes following a power law spectrum and measuring the cancellation function on these simulated magnetograms. Our results confirm the validity of the analytical relation and extend it to the domain {\alpha} < -1 where {\kappa} ---> 0 as {\alpha} ---> - infinity. The observationally derived upper limit of 0.38 on {\kappa} implies {\alpha} < -0.24 in the granular size range, apparently at odds with a small scale dynamo driven in the inertial range.

The tachocline is important in the solar dynamo for the generation and the storage of the magnetic fields. A most plausible explanation for the confinement of the tachocline is given by the fast tachocline model in which the tachocline is confined by the anisotropic momentum transfer by the Maxwell stress of the dynamo generated magnetic fields. We employ a flux transport dynamo model coupled with the simple feedback formula of this fast tachocline model which basically relates the thickness of the tachocline to the Maxwell stress. We find that this nonlinear coupling not only produces a stable solar-like dynamo solution but also a significant latitudinal variation in the tachocline thickness which is in agreement with the observations.

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 latitude dependence are in excellent agreement with seismic results. In the 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.

A review of solar cycle prediction methods and their performance is given, including forecasts for cycle 24 and focusing on aspects of the solar cycle prediction problem that have a bearing on dynamo theory. The scope of the review is further restricted to the issue of predicting the amplitude (and optionally the epoch) of an upcoming solar maximum no later than right after the start of the given cycle. Prediction methods form three main groups. Precursor methods rely on the value of some measure of solar activity or magnetism at a specified time to predict the amplitude of the following solar maximum. Their implicit assumption is that each numbered solar cycle is a consistent unit in itself, while solar activity seems to consist of a series of much less tightly intercorrelated individual cycles. Extrapolation methods, in contrast, are based on the premise that the physical process giving rise to the sunspot number record is statistically homogeneous, i.e., the mathematical regularities underlying its variations are the same at any point of time, and therefore it lends itself to analysis and forecasting by time series methods. Finally, instead of an analysis of observational data alone, model based predictions use physically (more or less) consistent dynamo models in their attempts to predict solar activity. In their overall performance precursor methods have clearly been superior to extrapolation methods. Nevertheless, some extrapolation methods may still be worth further study. Model based forecasts have not yet have had a chance to prove their skills. One method that has yielded predictions consistently in the right range during the past few solar cycles is that of K. Schatten et al., whose approach is mainly based on the polar field precursor. The incipient cycle 24 will probably mark the end of the Modern Maximum, with the Sun switching to a state of less strong activity.

Parker's analytic Cartesian interface dynamo is generalized to the case of a shear layer of finite thickness and low resistivity ("tachocline"), bounded by a perfect conductor ("radiative zone") on the one side, and by a highly diffusive medium ("convective zone") supporting an $\alpha$-effect on the other side. In the limit of high diffusivity contrast between the shear layer and the diffusive medium, thought to be relevant for the Sun, a pair of exact dispersion relations for the growth rate and frequency of dynamo modes is analytically derived. Graphic solution of the dispersion relations displays a somewhat unexpected, non-monotonic behaviour, the mathematical origin of which is elucidated. The dependence of the results on the parameter values (dynamo number and shear layer thickness) is investigated. The implications of this result for the solar dynamo problem are discussed.

A didactic introduction to current thinking on some aspects of the solar dynamo is given for geophysicists and planetary scientists.

While obviously having a common root, solar and planetary dynamo theory have taken increasingly divergent routes in the last two or three decades, and there are probably few experts now who can claim to be equally versed in both. Characteristically, even in the fine and comprehensive book "The magnetic Universe" (Rudiger & Hollerbach 2004), the chapters on planets and on the Sun were written by different authors. Separate reviews written on the two topics include Petrovay (2000}, Charbonneau (2005), Choudhuri (2008) on the solar dynamo and Glatzmaier (2002), Stevenson (2003) on the planetary dynamo. In the following I will try to make a systematic comparison between solar and planetary dynamos, presenting analogies and differences, and highlighting some interesting recent results.

The chemistry of molecular clouds has been studied for decades, with an increasingly general and sophisticated treatment of the reactions involved. Yet the treatment of turbulent diffusion has remained extremely sketchy, assuming simple Fickian diffusion with a scalar diffusivity D. However, turbulent flows similar to those in the interstellar medium are known to give rise to anomalous diffusion phenomena, more specifically superdiffusion (increase of the diffusivity with the spatial scales involved). This paper considers to what extent and in what sense superdiffusion modifies molecular abundances in interstellar clouds. For this first exploration of the subject we employ a very rough treatment of the chemistry and the effect of non-unifom cloud density on the diffusion equation is also treated in a simplified way. The results nevertheless clearly demonstrate that the effect of superdiffusion is quite significant, abundance values at a given radius being modified by order of unity factors. The sense and character of this influence is highly nontrivial.

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 analoguous 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.

Recently wide publicity has been given to a claim by T. Vachaspati that "black holes do not exist", that the objects known as black holes in astrophysics should rather be called "black stars" and they not only do not have event horizons but actually can be the source of spectacular gamma ray bursts. In this short essay (no flimsier than the original preprint where these extravagant claims appeared) I demonstrate that these ill-considered claims are clearly wrong. Yet they present a good occasion to reflect on some well known but little discussed conceptual difficulties which arise when applying relativistic terminology in an astrophysical context.

Observations indicate that in plage areas (i.e. in active regions outside sunspots) acoustic waves travel faster than in quiet sun, leading to shortened travel times and higher p-mode frequencies. While it is clear that the ultimate cause of any difference between quiet sun and plage is the presence of magnetic fields of order 100 G in the latter, the mechanism by which the magnetic field exerts its influence has not yet been conclusively identified. One possible such mechanism is suggested by the observation that granular motions in plage areas tend to be slightly ``abnormal'', dampened compared to quiet sun. In this paper we consider the effect that abnormal granulation observed in active regions should have on the propagation of acoustic waves. Any such effect is found to be limited to a shallow surface layer where sound waves propagate nearly vertically. The magnetically suppressed turbulence implies higher sound speeds, leading to shorter travel times. This time shift Dt is independent of the travel distance, while it shows a characteristic dependence on the assumed plage field strength. As a consequence of the variation of the acoustic cutoff with height, Dt is expected to be significantly higher for higher frequency waves within the observed regime of 3-5 mHz. The lower group velocity near the upper reflection point further leads to an increased envelope time shift, as compared to the phase shift. $p$-mode frequencies in plage areas are increased by a corresponding amount, Dnu/nu = nu*Dt. These characteristics of the time and frequency shifts are in accordance with observations. The calculated overall amplitude of the time and frequency shifts are comparable to, but still significantly (factor of 2 to 5) less than suggested by measurements.

The mean value of the normalized current helicity in solar active regions is on the order of 1e-8 1/m, negative in the northern hemisphere, positive in the southern hemisphere. Observations indicate that this helicity has a subsurface origin. Possible mechanisms leading to a twist of this amplitude in magnetic flux tubes include the solar dynamo, convective buffeting of rising flux tubes, and the accretion of weak external poloidal flux by a rising toroidal flux tube. After briefly reviewing the observational and theoretical constraints on the origin of helicity, we present a recently developed detailed model for poloidal flux accretion.

One possible scenario for the origin of the solar tachocline, known as the "fast tachocline", assumes that the turbulent diffusivity exceeds eta>10^9 cm^2/s. In this case the dynamics will be governed by the dynamo-generated oscillatory magnetic field on relatively short timescales. Here, for the first time, we present detailed numerical models for the fast solar tachocline with all components of the magnetic field calculated explicitly, assuming axial symmetry and a constant turbulent diffusivity eta and viscosity nu. We find that a sufficiently strong oscillatory poloidal field with dipolar latitude dependence at the tachocline-convective zone boundary is able to confine the tachocline. Exploring the three-dimensional parameter space defined by the viscosity in the range log(nu)=9-11, the magnetic Prandtl number in the range Prm=0.1-10, and the meridional flow amplitude (-3 to +3 cm/s), we also find that the confining field strength B_conf, necessary to reproduce the observed thickness of the tachocline, increases with viscosity nu, with magnetic Prandtl number nu/eta, and with equatorward meridional flow speed. Nevertheless, the resulting B_conf values remain quite reasonable, in the range 10^3-10^4 G, for all parameter combinations considered here. The thickness of the tachocline shows a marked dependence on both time and latitude. A comparison with seismic constraints suggests that best agreement with our models is achieved for the highest values of nu and Prm considered here.

Helioseismic measurements indicate that the solar tachocline is very thin, its full thickness not exceeding 4% of the solar radius. The mechanism that inhibits differential rotation to propagate from the convective zone to deeper into the radiative zone is not known, though several propositions have been made. In this paper we demonstrate by numerical models and analytic estimates that the tachocline can be confined to its observed thickness by a poloidal magnetic field B_p of about one kilogauss, penetrating below the convective zone and oscillating with a period of 22 years, if the tachocline region is turbulent with a diffusivity of eta 10^10 cm^2/s (for a turbulent magnetic Prandtl number of unity). We also show that a similar confinement may be produced for other pairs of the parameter values (B_p, eta). The assumption of the dynamo field penetrating into the tachocline is consistent whenever eta>10^9 cm^2/s.

Following an earlier proposal by Choudhuri (2003) for the origin of twist in the magnetic fields of solar active regions, we model the penetration of a wrapped up background poloidal field into a toroidal magnetic flux tube rising through the solar convective zone. The rise of the straight, cylindrical flux tube is followed by numerically solving the induction equation in a comoving Lagrangian frame, while an external poloidal magnetic field is assumed to be radially advected onto the tube with a speed corresponding to the rise velocity. One prediction of our model is the existence of a ring of reverse current helicity on the periphery of active regions. On the other hand, the amplitude of the resulting twist depends sensitively on the assumed structure (diffuse vs. concentrated/intermittent) of the active region magnetic field right before its emergence, and on the assumed vertical profile of the poloidal field. Nevertheless, in the model with the most plausible choice of assumptions a mean twist comparable to the observations can result rather naturally. Our results indicate that the contribution of this mechanism to the twist can be quite significant, and under favourable circumstances it can potentially account for most of the current helicity observed in active regions.

Parker's interface dynamo is generalized to the case when a homogeneous flow is present in the high-diffusivity (upper) layer in the lateral direction (i.e. perpendicular to the shear flow in the lower layer). This is probably a realistic first representation of the situation near the bottom of the solar convective zone, as the strongly subadiabatic stratification of the tachocline (lower layer in the interface dynamo) imposes a strong upper limit on the speed of any meridional flow there. Analytic solutions to the eigenvalue problem are presented for the cases of vanishing diffusivity contrast and infinite diffusivity contrast, respectively. Unlike the trivial case of a homogeneous system, the ability of the meridional flow to reverse the propagation of the dynamo wave is strongly reduced in the interface dynamo. In particular, in the limit of high diffusivity contrast relevant to the solar case it is found that a meridional flow of realistic amplitude cannot reverse the direction of propagation of the dynamo wave. The implications of this result for the solar dynamo problem are discussed.

The first consistent model for the turbulent tachocline is presented, with the turbulent diffusivity computed within the model instead of being specified arbitrarily. For the origin of the 3D turbulence a new mechanism is proposed. Owing to the strongly stable stratification, the mean radial shear is stable, while the horizontal shear is expected to drive predominantly horizontal, quasi-2D motions in thin slabs. Here I suggest that a major source of 3D overturning turbulent motions in the tachocline is the secondary shear instability due to the strong, random vertical shear arising between the uncorrelated horizontal flows in neighbouring slabs. A formula for the vertical diffusivity due to this turbulence, equation (9), is derived and applied in a simplified 1D model of the tachocline. It is found that Maxwell stresses due to an oscillatory poloidal magnetic field of a few hundred gauss are able to confine the tachocline to a thickness below 5 Mm. The integral scale of the 3D overturning turbulence is the buoyancy scale, on the order of 10 km and its velocity amplitude is a few m/s, yielding a vertical turbulent diffusivity on the order of 10^8 cm^2/s.

We present a model for the lower overshoot layer of the Sun, based on the realistic solar stratification, without the use of a "mixing-length" parameter, by solving the system of Reynolds momentum equations using the closure formalism of Canuto & Dubovikov (1997, 1998). A fixed value of velocity anisotropy is assumed, and the local convection model is assumed to be valid for the convectively unstable layer. In accordance with seismic constraints, overshoot (defined as the amount by which the convectively mixed zone extends beyond its boundary in local theory) is found to be as low as about 6 percent of the pressure scale height, and it is not bounded by a discontinuity from below.

We present a model for torsional oscillations where the inhibiting effect of active region magnetic fields on turbulence locally reduces turbulent viscous torques, leading to a cycle- and latitude-dependent modulation of the differential rotation. The observed depth dependence of torsional oscillations as well as their phase relationship with the sunspot butterfly diagram are reproduced quite naturally in this model. The resulting oscillation amplitudes are significantly smaller than observed, though they depend rather sensitively on model details. Meridional circulation is found to have only a weak effect on the oscillation pattern.

In contrast to the situation with the geodynamo, no breakthrough has been made in the solar dynamo problem for decades. Since the appearance of mean-field electrodynamics in the 1960's, the only really significant advance was in the field of flux tube theory and flux emergence calculations. These new results, together with helioseismic evidence, have led to the realization that the toroidal magnetic flux giving rise to activity phenomena must be stored and presumably generated below the convection zone proper, in what I will call the DOT (Dynamo-Overshoot-Tachoclyne) layer. The only segment of the problem we can claim to basically understand is the transport of flux from this layer to the surface. On the other hand, as reliable models for the DOT layer do not exist we are clueless concerning the precise mechanisms responsible for toroidal/poloidal flux conversion and for characteristic migration patterns (extended butterfly diagram) and periodicities. Even the most basic result of mean-field theory, the interpretation of the butterfly diagram as an alpha-omega dynamo wave, has been questioned. This review therefore will necessarily ask more questions than give answers. Some of these key questions are - Structure of the DOT layer - alpha-quenching and distributed dynamo - High-latitude migration patterns and their interpretation - The ultimate fate of emerged flux

The precise nature of photospheric flows, and of the transport effects they give rise to, has been the subject of intense debate in the last decade. Here we attempt to give a brief review of the subject emphasizing interdisciplinary (solar physics - turbulence theory) aspects, key open questions, and recent developments.

In a statistical analysis of Debrecen Photoheliographic Results sunspot area data we find that the logarithmic deviation (log D)' of the area decay rate D from the parabolic mean decay law (derived in the first paper in this series) follows a Gaussian probability distribution. As a consequence, the actual decay rate D and the its time average are also characterized by approximately lognormal distributions, as found in an earlier work. The correlation time of (log D)' is about 3 days. We find a significant physical anticorrelation between (log D)' and the amount of plage magnetic flux of the same polarity in an annulus around the spot on Kitt Peak magnetograms. The anticorrelation is interpreted in terms of a generalization of the turbulent erosion model of sunspot decay to the case when the flux tube is embedded in a preexisting homogeneous "plage" field. The decay rate is found to depend inversely on the value of this plage field, the relation being very close to logarithmic, i.e. the plage field acts as multiplicative noise in the decay process. A Gaussian probability distribution of the field strength in the surrounding plage will then naturally lead to a lognormal distribution of the decay rates, as observed. It is thus suggested that, beside other multiplicative noise sources, the environmental effect of surrounding plage fields is a major factor in the origin of lognormally distributed large random deviations from the mean law in the sunspot decay rates.

The axisymmetric component of the large-scale solar magnetic fields has a pronounced poleward branch at higher latitudes. In order to clarify the origin of this branch we construct an axisymmetric model of the passive transport of the mean poloidal magnetic field in the convective zone, including meridional circulation, anisotropic diffusivity, turbulent pumping and density pumping. For realistic values of the transport coefficients we find that diffusivity is prevalent, and the latitudinal distribution of the field at the surface simply reflects the conditions at the bottom of the convective zone. Pumping effects concentrate the field to the bottom of the convective zone; a significant part of this pumping occurs in a shallow subsurface layer, normally not resolved in dynamo models. The phase delay of the surface poloidal field relative to the bottom poloidal field is found to be small. These results support the double dynamo wave models, may be compatible with some form of a mixed transport scenario, and exclude the passive transport theory for the origin of the polar branch.

Results from a numerical and analytical investigation of the solution of a nonlinear axially symmetric diffusion equation for the magnetic field are presented for the case when the nonlinear dependence of the diffusivity nu(B) on the magnetic field satisfies basic physical requirements. We find that for sufficiently strong nonlinearity (i.e. for sufficiently strong reduction of nu inside the tube) a current sheet is spontaneously formed around the tube within one diffusion timescale. This sheet propagates inwards with a velocity inversely proportional to the ratio of the field strength just inside the current sheet to the equipartition field strength B0/Be, so the lifetime of a tube with constant internal flux density is increased approximately by a factor not exceeding B0/Be, even for infinitely effective inhibition of turbulence inside the tube. Among the applications of these results we point out that toroidal flux tubes in the solar convective zone are subject to significant flux loss owing to turbulent erosion on a timescale of about 1 month, and that turbulent erosion may be responsible for the formation of a current sheet around a sunspot. It is further proposed that, despite the simplifying assumptions involved, our solutions correctly reflect the essential features of the sunspot decay process.

In a statistical study of the decay of individual sunspots based on DPR data we find that the mean instantaneous area decay rate is related to the spot radius r and the maximal radius r0 as D=C r/r0, C=32.0 +- 0.26 MSH/day. This implies that sunspots on the mean follow a parabolic decay law; the traditional linear decay law is excluded by the data. The validity of the Gnevyshev-Waldmeier relationship between the maximal area A0 and lifetime T of a spot group, A0/T = 10 MSH/day, is also demonstrated for individual sunspots. No evidence is found for a supposed supergranular "quantization" of sunspot areas. Our results strongly support the recent turbulent erosion model of sunspot decay while all other models are excluded.

The degree of linear polarization has recently been found to show wide random variations over the solar disk. These variations are presumably at least partly due to fluctuations in the flux density of turbulent photospheric magnetic fields and associated variations in the degree of Hanle depolarization. In order to understand the origin of such large scale fluctuations of the turbulent magnetic flux density we develop a phenomenological model to calculate the spatial Fourier spectrum of the fluctuations of turbulent magnetic fields in the solar photosphere and convective zone. It is found that if the model parameters are fitted to turbulence closure models and numerical experiments the characteristic scale of the fluctuations is by about an order of magnitude larger than the turbulence scale (the scale of the granulation), owing to the more effective quenching of small-scale fluctuations by turbulent diffusion.

In recent years, our knowledge of photospheric magnetic fields went through a thorough transformation--nearly unnoticed by dynamo theorists. It is now practically certain that the overwhelming majority of the unsigned magnetic flux crossing the solar surface is in turbulent form (intranetwork and hidden fields). Furthermore, there are now observational indications (supported by theoretical arguments discussed in this paper) that the net polarity imbalance of the turbulent field may give a significant or even dominant contribution to the weak large-scale background magnetic fields outside unipolar network areas. This turbulent magnetic field consists of flux tubes with magnetic fluxes below 1e10 Wb (1e18 Mx). The motion of these thin tubes is dominated by the drag of the surrounding flows, so the transport of this component of the solar magnetic field must fully be determined by the kinematics of the turbulence (i.e. it is "passive"), and it can be described by a one-fluid model like mean-field theory (MFT). This paper reviews the theory of passive magnetic field transport using mostly first (and occasionally higher) order smoothing formalism; the most important transport effects are however also independently derived using Lagrangian analysis for a simple two-component flow model. Solar applications of the theory are also presented. Among some other novel findings it is proposed that the observed unsigned magnetic flux density in the photosphere requires a small-scale dynamo effect operating in the convective zone and that the net polarity imbalance in turbulent (and, in particular, hidden) fields may give a major contribution to the weak large-scale background magnetic fields on the Sun.

A model for homogeneous anisotropic incompressible turbulence is proposed. The model generalizes the GISS model of homogeneous isotropic turbulence; the generalization involves the solution of the GISS equations along a set of integration paths in wavenumber (k-) space. In order to make the problem tractable, these integration paths ("cascade lines") must be chosen in such a way that the behaviour of the energy spectral function along different cascade lines should be reasonably similar. In practice this is realized by defining the cascade lines as the streamlines of a cascade flow; in the simplest case the source of this flow may be identified with the source function of the turbulence. Owing to the different approximations involved, the resulting energy spectral function is not exact but is expected to give good approximative values for the bulk quantities characterising the turbulent medium, and for the measure of the anisotropy itself in particular. The model is then applied to the case of low Prandtl number thermal convection. The energy spectral function and the bulk quantities characterizing the flow are derived for different values of the parameter $S=Ra \sigma$. The most important new finding is that unlike the anisotropy of the most unstable mode in linear stability analysis the anisotropy of the turbulence does not grow indefinitely with increasing S but it rather saturates to a relatively moderate finite asymptotic value.