The profile of the differential rotation together with the sign of the alpha-effect determine the dynamo wave direction. In early models of the solar dynamo the dynamo wave often leads to a poleward migration of the activity belts. Flux transport by the meridional flow or the effect of the surface shear layer are possible solutions. In a model including the corona, we show that various migrations can be obtained by varying the properties of the corona. A new dynamo of Babcock-Leighton type also leads to the correct equatorward migration by the non-linear relation between flux density and rise time of the flux.

As we are heading towards the next solar cycle, presumably with a relatively small amplitude, it is of significant interest to reconstruct and describe the past grand minima on the basis of actual observations of the time. The Dalton Minimum is often considered one of the grand minima captured in the coverage of telescopic observations. Nevertheless, the reconstructions of the sunspot group number vary significantly, and the existing butterfly diagrams have a large data gap during the period. This is partially because most long-term observations have remained unexplored in historical archives. Therefore, to improve our understanding on the Dalton Minimum, we have located two series of Thadd\"aus Derfflinger's observational records (a summary manuscript and logbooks) as well as his Brander's 5.5-feet azimuthal-quadrant preserved in the Kremsm\"unster Observatory. We have revised the existing Derfflinger's sunspot group number with Waldmeier classification and eliminated all the existing 'spotless days' to remove contaminations from solar meridian observations. We have reconstructed the butterfly diagram on the basis of his observations and illustrated sunspot distributions in both solar hemispheres. Our article aims to revise the trend of Derfflinger's sunspot group number and to bridge a data gap of the existing butterfly diagrams around the Dalton Minimum. Our results confirm that the Dalton Minimum is significantly different from the Maunder Minimum, both in terms of cycle amplitudes and sunspot distributions. Therefore, the Dalton Minimum is more likely a secular minimum in the long-term solar activity, while further investigations for the observations at that time are required.

Between 1761 and 1776, Christian Horrebow made regular observations of sunspots from Rundetaarn in Copenhagen. Based on these observations he writes in 1775 that "it appears that after the course of a certain number of years, the appearance of the Sun repeats itself with respect to the number and size of the spots". Thus, Horrebow hypothesized the idea of a cyclic Sun several decades before Heinrich Schwabe discovered the solar cycle and estimated its period. This proves the ability of Horrebow as a sunspot observer. In this article, we present a general overview of the work of Christian Horrebow, including a brief biography and a complete bibliography. We also present a translation from Danish to English of his writings on sunspots in the Dansk Historisk Almanak. These writings include tables of daily sunspot measurements of which we discuss the completeness.

The number of spots on the surface of the Sun is one of the best tracers of solar variability we have. The sunspot number is not only known to change in phase with the 11-year solar cycles, but also to show variability on longer time scales. It is, however, not only the sunspot number that changes in connection with solar variability. The location of the spots on the solar surface is also known to change in phase with the 11-year solar cycle. This has traditionally been visualised in the so-called butterfly diagram, but this is only well constrained from the beginning of the 19th century. This is unfortunate, as knowledge about the butterfly diagram could aid our understanding of the variability and the Sun-Earth connection. As part of a larger review of the work done on sunspots by the Danish astronomer Christian Horrebow, we here present a reanalysis of Christian Horrebow's notebooks covering the years 1761 and 1764 - 1777. These notebooks have been analysed in at least three earlier studies by Thiele (Astron. Nachr. 50, 257, 1859), d'Arrest (published in Wolf, Astron. Mitt. Eidgenoss. Sternwarte Zur. 4, 77, 1873) and Hoyt and Schatten (Solar Phys. 160, 387, 1995). In this article, we construct a complete record of sunspot positions covering the years 1761 and 1764 - 1777. The resulting butterfly diagram shows the characteristic structure known from observations in the 19th and 20th century. We do see some indications of equatorial sunspots in the observations we have from Cycle 1. However, in Cycle 2, which has much better coverage, we do not see such indications.

Contex: Solar magnetic activity exhibits variations with periods between 1.5--4 years, the so-called quasi-biennial oscillations (QBOs), in addition to the well-known 11-year Schwabe cycles. Solar dynamo is thought to be the responsible mechanism for generation of the QBOs. Aims: In this work, we analyse sunspot areas to investigate the spatial and temporal behaviour of the QBO signal and study the responsible physical mechanisms using simulations from fully nonlinear mean-field flux-transport dynamos. Methods: We investigated the behaviour of the QBOs in the sunspot area data in full disk, and northern and southern hemispheres, using wavelet and Fourier analyses. We also ran solar dynamos with two different approaches to generating a poloidal field from an existing toroidal field, Babcock-Leighton and turbulent $\alpha$ mechanisms. We then studied the simulated magnetic field strengths as well as meridional circulation and differential rotation rates using the same methods. Results: The results from the sunspot areas show that the QBOs are present in the full disk and hemispheric sunspot areas and they show slightly different spatial and temporal behaviours, indicating a slightly decoupled solar hemispheres. The QBO signal is generally intermittent and in-phase with the sunspot area data, surfacing when the solar activity is in maximum. The results from the BL-dynamos showed that they are neither capable of generating the slightly decoupled behaviour of solar hemispheres nor can they generate QBO-like signals. The turbulent $\alpha$-dynamos, on the other hand, generated decoupled hemispheres and some QBO-like shorter cycles. Conclusions: In conclusion, our simulations show that the turbulent $\alpha$-dynamos with the Lorentz force seems more efficient in generating the observed temporal and spatial behaviour of the QBO signal compared with those from the BL-dynamos.

Although the occurrence of solar irradiance variations induced by magnetic surface features (e.g., sunspots, faculae, magnetic network) is generally accepted, the existence of intrinsic luminosity changes due to the internal magnetic fields is still controversial. This additional contribution is expected to be accompanied by radius variations, and to be potentially significant for the climate of the Earth. We aim to constrain theoretically the radius and luminosity variations of the Sun that are due to the effect of the variable magnetic fields in its interior associated with the dynamo cycle. We have extended a one-dimensional stellar evolution code to include several effects of the magnetic fields on the interior structure. We investigate different magnetic configurations, based on both observational constraints and on the output of state-of-the-art mean field dynamo models. We explore both step-like and simply periodic time dependences of the magnetic field peak strength. We find that the luminosity and radius variations are in anti-phase and in phase, respectively, with the magnetic field strength. For peak magnetic field strengths of the order of tens of kilogauss, luminosity variations ranging between 10^{-6} and 10^{-3} (in modulus) and radius variations between 10^{-6} and 10^{-5} are obtained. Modest but significant radius variations (up to 10^{-5} in relative terms) are obtained for magnetic fields of realistic strength and geometry, providing a potentially observable signature of the intrinsic variations. Establishing their existence in addition to the accepted surface effects would have very important implications for the understanding of solar-induced long-term trends on climate.

Context. Solar dynamo models of Babcock-Leighton type typically assume the rise of magnetic flux tubes to be instantaneous. Solutions with high-magnetic-diffusivity have too short periods and a wrong migration of their active belts. Only the low-diffusivity regime with advective meridional flows is usually considered. Aims. In the present paper we discuss these assumptions and applied a time delay in the source term of the azimuthally averaged induction equation. This delay is set to be the rise time of magnetic flux tubes which supposedly form at the tachocline. We study the effect of the delay, which adds to the spacial non-locality a non-linear temporal one, in the advective but particularly in the diffusive regime. Methods. Fournier et al. (2017) obtained the rise time according to stellar parameters such as rotation, and the magnetic field strength at the bottom of the convection zone. These results allowed us to constrain the delay in the mean-field model used in a parameter study. Results. We identify an unknown family of solutions. These solutions self-quench, and exhibit longer periods than their non-delayed counterparts. Additionally, we demonstrate that the non-linear delay is responsible for the recover of the equatorward migration of the active belts at high turbulent diffusivities. Conclusions. By introducing a non-linear temporal non-locality (the delay) in a Babcock-Leighton dynamo model, we could obtain solutions quantitatively comparable to the solar butterfly diagram in the diffusion-dominated regime.

We present reconstructed sunspot positions based on observations reported in letters between Gottfried Kirch and other contemporary astronomers from AD 1680 to 1709, i.e. in the last decades of the Maunder Minimum. The letters were compiled and edited by Herbst (2006). The letters (and observations) from Kirch are mostly by Gottfried Kirch, but some also by his 2nd wife Maria M. Kirch (married 1692) and their son Christfried Kirch (born 1694). Using excerpts from the letters, some with drawings, we found some 35 sunspot groups (often for several days in a row or with interruptions) by Kirch and/or his letter partners (in three cases, only the month is given: 1704 Jan, Feb, 1707 Mar, otherwise always the exact dates) - usually one group at a time. We also found 17 explicit spotless days, several of them new (previously without any known observations). We could constrain the heliographic latitude by Bayesian inference for 19 sunspot groups - five of them completely new (one group 1680 May 20-22 from Kirch and Ihle, one to two groups 1680 Jun 15-23 for Kirch, one group 1684 May 6 from Ihle, and one group 1688 Dec 14-15 from Kirch), while the others mostly agree (within 2sigma) with previously published values for those dates by others. With these data, we then amend the butterfly diagram for the Maunder Minimum. By comparison of our data with the sunspot group catalog in Hoyt & Schatten (1998), we noticed a number of discrepancies, e.g. that dates for British observers in the Maunder Minimum (Flamsteed, Caswell, Derham, Stannyan, Gray, and Sharp) as listed in Hoyt & Schatten (1998) are their original Julian dates, not converted to the Gregorian calendar (10-11 day offset in Hoyt & Schatten). Most of these modifications also apply to the modified sunspot group catalog in Vaquero et al. (2016). We also present two aurorae observed by the Kirchs in 1707 and 1716.

We present and analyse the sunspot observations performed by Franz I. C. Hallaschka in 1814 and 1816. These solar observations were carried out during the so-called Dalton minimum, around the maximum phase of the Solar Cycle 6. These records are very valuable because they allow us to complete observational gaps in the collection of sunspot group numbers, improving its coverage for this epoch. We have analysed and compared the observations made by Hallaschka with the records made by other contemporary observers. Unfortunately, the analysis of the sunspot areas and positions showed that they are too inaccurate for scientific use. But, we conclude that sunspot counts made by Hallaschka are similar to those made by other astronomers of that time. The observations by Hallaschka confirm a low level of the solar activity during the Dalton minimum.

We aim to investigate the nature and occurrence characteristics of grand solar minimum and maximum periods, which are observed in the solar proxy records such as 10Be and 14C, using a fully non-linear Babcock-Leighton type flux-transport dynamo including momentum and entropy equations. The differential rotation and meridional circulation are generated from the effect of turbulent Reynolds stress and are subjected to back-reaction from the magnetic field. To generate grand minimum and maximum-like periods in our simulations, we used random fluctuations in the angular momentum transport process, namely the Lambda-mechanism, and in the Babcock-Leighton mechanism. To characterise the nature and occurrences of the identified grand minima and maxima in our simulations, we used the waiting time distribution analyses, which reflects whether the underlying distribution arises from a random or a memory-bearing process. The results show that, in majority of the cases, the distributions of grand minima and maxima reveal that the nature of these events originates from memoryless processes. We also found that in our simulations the meridional circulation speed tends to be smaller during grand maximum, while it is faster during grand minimum periods. The radial differential rotation tend to be larger during grand maxima, while it is smaller during grand minima. The latitudinal differential rotation on the other hand is found to be larger during grand minima.

We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ.

Context: Long-term variability in solar cycles represents a challenging constraint for theoretical models. Mean-field Babcock-Leighton dynamos that consider non-instantaneous rising flux tubes have been shown to exhibit long-term variability in their magnetic cycle. However a relation that parameterizes the rise-time of non-axisymmetric magnetic flux tubes in terms of stellar parameters is still missing. Aims: We aim to find a general parameterization of the rise-time of magnetic flux tubes for solar-like stars. Methods: By considering the influence of magnetic tension on the rise of non-axisymmetric flux tubes, we predict the existence of a control parameter referred as $\Gamma_{\alpha_1}^{\alpha_2}$. This parameter is a measure of the balance between rotational effects and magnetic effects (buoyancy and tension) acting on the magnetic flux tube. We carry out two series of numerical experiments (one for axisymmetric rise and one for non-axisymmetric rise) and demonstrate that $\Gamma_{\alpha_1}^{\alpha_2}$ indeed controls the rise-time of magnetic flux tubes. Results: We find that the rise-time follows a power law of $\Gamma_{\alpha_1}^{\alpha_2}$ with an exponent that depends on the azimuthal wavenumber of the magnetic flux loop. Conclusions: Compressibility does not impact the rise of magnetic flux tubes, while non-axisymmetry does. In the case of non-axisymmetric rise, the tension force modifies the force balance acting on the magnetic flux tube. We identified the three independent parameters required to predict the rise-time of magnetic flux tubes, that is, the stellar rotation rate, the magnetic flux density of the flux tube, and its azimuthal wavenumber. We combined these into one single relation that is valid for any solar-like star. We suggest using this generalized relation to constrain the rise-time of magnetic flux tubes in Babcock-Leighton dynamo models.

We describe a revised collection of the number of sunspot groups from 1610 to the present. This new collection is based on the work of Hoyt and Schatten (Solar Phys. 179, 189, 1998). The main changes are the elimination of a considerable number of observations during the Maunder Minimum (hereafter, MM) and the inclusion of several long series of observations. Numerous minor changes are also described. Moreover, we have calculated the active-day percentage during the MM from this new collection as a reliable index of the solar activity. Thus, the level of solar activity obtained in this work is greater than the level obtained using the original Hoyt and Schatten data, although it remains compatible with a grand minimum of solar activity. The new collection is available in digital format.

Digital images of the observations printed in the books "Rosa Ursina sive solis" and "Prodromus pro sole mobili" by Christoph Scheiner as well as the drawings from Scheiner's letters to Marcus Welser are analysed in order to obtain information on positions and sizes of sunspots that appeared before the Maunder minimum. In most cases, the given orientation of the ecliptic is used to set up the heliographic coordinate system for the drawings. Positions and sizes are measured manually on the screen. Very early drawings have no indication of their orientation. A rotational matching using common spots of adjacent days is used in some cases, while in other cases, the assumption of images being aligned with a zenith-horizon coordinate system appeared to be the most probable. In total, 8167 sunspots were measured. A distribution of sunspot latitudes versus time (butterfly diagram) is obtained for Scheiner's observations. The observations of 1611 are very inaccurate, the drawings of 1612 have at least an indication of their orientation, while the remaining part of the spot positions from 1618-1631 have good to very good accuracy. We also computed 697 tilt angles of apparently bipolar sunspot groups observed in the period 1618-1631. We find that the average tilt angle of nearly 4 degrees is not significantly different from 20th-century values.

Sunspot positions from various historical sets of solar drawings are analysed with respect to the tilt angles of bipolar sunspot groups. Data by Scheiner, Hevelius, Staudacher, Zucconi, Schwabe, and Spoerer deliver a series of average tilt angles spanning a period of 270 years, additional to previously found values for 20th-century data obtained by other authors. We find that the average tilt angles before the Maunder minimum were not significantly different from the modern values. However, the average tilt angles of a period 50 years after the Maunder minimum, namely for cycles 0 and 1, were much lower and near zero. The normal tilt angles before the Maunder minimum suggest that it was not abnormally low tilt angles which drove the solar cycle into a grand minimum.

Context. Using asteroseismic techniques, it has recently become possible to probe the internal rotation profile of low-mass (~1.1-1.5 Msun) subgiant and red giant stars. Under the assumption of local angular momentum conservation, the core contraction and envelope expansion occurring at the end of the main sequence would result in a much larger internal differential rotation than observed. This suggests that angular momentum redistribution must be taking place in the interior of these stars. Aims. We investigate the physical nature of the angular momentum redistribution mechanisms operating in stellar interiors by constraining the efficiency of post-main sequence rotational coupling. Methods. We model the rotational evolution of a 1.25 Msun star using the Yale Rotational stellar Evolution Code. Our models take into account the magnetic wind braking occurring at the surface of the star and the angular momentum transport in the interior, with an efficiency dependent on the degree of internal differential rotation. Results. We find that models including a dependence of the angular momentum transport efficiency on the radial rotational shear reproduce very well the observations. The best fit of the data is obtained with an angular momentum transport coefficient scaling with the ratio of the rotation rate of the radiative interior over that of the convective envelope of the star as a power law of exponent ~3. This scaling is consistent with the predictions of recent numerical simulations of the Azimuthal Magneto-Rotational Instability. Conclusions. We show that an angular momentum transport process whose efficiency varies during the stellar evolution through a dependence on the level of internal differential rotation is required to explain the observed post-main sequence rotational evolution of low-mass stars.

Extending the knowledge about the properties of solar cycles into the past is essential for understanding the solar dynamo. This paper aims at estimating areas of sunspots observed by Schwabe in 1825-1867 and at calculating the tilt angles of sunspot groups. The sunspot sizes in Schwabe's drawings are not to scale and need to be converted into physical sunspot areas. We employed a statistical approach assuming that the area distribution of sunspots was the same in the 19th century as it was in the 20th century. Umbral areas for about 130,000 sunspots observed by Schwabe were obtained, as well as the tilt angles of sunspot groups assuming them to be bipolar. There is, of course, no polarity information in the observations. The annually averaged sunspot areas correlate reasonably with sunspot number. We derived an average tilt angle by attempting to exclude unipolar groups with a minimum separation of the two alleged polarities and an outlier rejection method which follows the evolution of each group and detects the moment it turns unipolar at its decay. As a result, the tilt angles, although displaying considerable scatter, place the leading polarity on average 5.85+-0.25 closer to the equator, in good agreement with tilt angles obtained from 20th-century data sets. Sources of uncertainties in the tilt angle determination are discussed and need to be addressed whenever different data sets are combined. The sunspot area and tilt angle data are provided online.

We present a few newly found old sunspot observations from the years AD 1708, 1709, and 1710, which were obtained by Peter Becker from Rostock, Germany. For 1709, Becker gave a detailed drawing: he observed a sunspot group made up of two spots on Jan 5, 6, and 7, and just one of the two spots was observed on Jan 8 and 9. We present his drawing and his explanatory text. We can measure the latitude and longitude of these two spots and estimate their sizes for all five days. While the spots and groups in 1708 and the spot on four of the five days in January 1709 were known before from other observers (e.g. Hoyt & Schatten 1998), the location of the spots in early January 1709 were not known before, so that they can now be considered in reconstructed butterfly diagrams. The sunspots detected by Becker on 1709 Jan 5 and 1710 Sep 10 were not known before at all, as the only observer known for those two dates, La Hire, did not detect that spot (group). We estimate new group sunspot numbers for the relevant days, months, and years. The time around 1708--1710 is important, because it documents the recovery of solar activity towards the end of the Maunder Grand Minimum. We also show two new spot observations from G. Kirch for 1708 Sep 13 \& 14 as described in his letter to Wurzelbaur (dated Berlin AD 1708 Dec 19).

Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods: We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results: The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets. Conclusions: We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level.

Surface differential rotation (DR) is one major ingredient of the magnetic field generation process in the Sun and likely in other stars. The term solar-like differential rotation describes the observation that solar equatorial regions rotate faster than polar ones. The opposite effect of polar regions rotating faster than equatorial ones (termed as antisolar DR) has only been observed in a few stars, although there is evidence from theoretical dynamo models. We present a new method to detect the sign of DR (i.e. solar-like or antisolar DR) by analyzing long-term high-precision light curves with the Lomb-Scargle periodogram.We compute the Lomb-Scargle periodogram and identify a set of significant periods $P_k$, which we associate with active regions located at different latitudes on the the stellar surface. If detectable, the first harmonics ($P_k'$) of these periods were identified to compute their peak-height-ratios $r_k:=h(P_k')/h(P_k)$. Spots rotating at lower latitudes generate less sine-shaped light curves, which requires additional power in the harmonics, and results in larger ratios $r_k$. Comparing different ratios $r_k$ and the associated periods $P_k$ yields information about the spot latitudes, and reveals the sign of DR. We tested our method on different sets of synthetic light curves all exhibiting solar-like DR. The number of cases where our method detects antisolar DR is the false-positive rate of our method. Depending on the set of light curves, the noise level, the required minimum peak separation, and the presence or absence of spot evolution, our method fails to detect the correct sign in at most 20%. We applied our method to 50 Kepler G stars and found 21-34 stars with solar-like DR and 5-10 stars with antisolar DR, depending on the minimum peak separation.

We investigate numerically the self-sustained dynamo action in a spinning sphere whose sense of rotation reverses periodically. This system serves as a simple model of a dynamo in small bodies powered by frequent collisions. It is found that dynamo action is possible in some intervals of collision rates. At high Ekman numbers the laminar spin-up flow is helical in the boundary layers and the Ekman circulation together with the azimuthal shear powers the dynamo action. At low Ekman number a non-axisymmetric instability helps the dynamo action. The intermittency of magnetic field occurs at low Ekman number.

The (Zeeman-)Doppler imaging studies of solar-type stars very often reveal large high-latitude spots. This also includes F stars that possess relatively shallow convection zones, indicating that the dynamo operating in these stars differs from the solar dynamo. We aim to determine whether mean-field dynamo models of late-F type dwarf stars can reproduce the surface features recovered in Doppler maps. In particular, we wish to test whether the models can reproduce the high-latitude spots observed on some F dwarfs. The photometric inversions and the surface temperature maps of AF Lep were obtained using the Occamian-approach inversion technique. Low signal-to-noise spectroscopic data were improved by applying the least-squares deconvolution method. The locations of strong magnetic flux in the stellar tachocline as well as the surface fields obtained from mean-field dynamo solutions were compared with the observed surface temperature maps. The photometric record of AF Lep reveals both long- and short-term variability. However, the current data set is too short for cycle-length estimates. From the photometry, we have determined the rotation period of the star to be 0.9660+-0.0023 days. The surface temperature maps show a dominant, but evolving, high-latitude (around +65 degrees) spot. Detailed study of the photometry reveals that sometimes the spot coverage varies only marginally over a long time, and at other times it varies rapidly. Of a suite of dynamo models, the model with a radiative interior rotating as fast as the convection zone at the equator delivered the highest compatibility with the obtained Doppler images.

Most of our knowledge about the Sun's activity cycle arises from sunspot observations over the last centuries since telescopes have been used for astronomy. The German astronomer Gustav Sp\"orer observed almost daily the Sun from 1861 until the beginning of 1894 and assembled a 33-year collection of sunspot data covering a total of 445 solar rotation periods. These sunspot drawings were carefully placed on an equidistant grid of heliographic longitude and latitude for each rotation period, which were then copied to copper plates for a lithographic reproduction of the drawings in astronomical journals. In this article, we describe in detail the process of capturing these data as digital images, correcting for various effects of the aging print materials, and preparing the data for contemporary scientific analysis based on advanced image processing techniques. With the processed data we create a butterfly diagram aggregating sunspot areas, and we present methods to measure the size of sunspots (umbra and penumbra) and to determine tilt angles of active regions. A probability density function of the sunspot area is computed, which conforms to contemporary data after rescaling.

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

The record of solar activity is reviewed here with emphasis on peculiarities. Since sunspot positions tell us a lot more about the solar dynamo than the various global sunspot numbers, we first focus on the records of telescopic observations of sunspots leading to positional information. Then we turn to the proxy record from cosmogenic isotope abundances, which shows recurrent grand minima over the last 9500 years. The apparent distinction between episodes of strong modulation, and intervening episodes with milder modulation and weaker overall activity, hints at the solar dynamo following a variety of solutions, with different symmetries, over the course of millennia.

Aims. Sunspot number is a benchmark series in many studies, but may still contain inhomogeneities and inconsistencies. In particular, an essential discrepancy exists between the two main sunspot number series, Wolf (WSN) and group (GSN) sunspot numbers, before 1848. The source of this discrepancy has so far remained unresolved. However, the recently digitized series of solar observations in 1825-1867 by Samuel Heinrich Schwabe, who was the primary observer of the WSN before 1848, makes such an assessment possible. Methods. We construct sunspot series, similar to WSN and GSN, but using only Schwabe's data. These series, called WSN-S and GSN-S, respectively, were compared with the original WSN and GSN series for the period 1835-1867 to look for possible inhomogeneities. Results. We show that: (1) The GSN series is homogeneous and consistent with the Schwabe data throughout the entire studied period; (2) The WSN series decreases by roughly ~20% around 1848 caused by the change of the primary observer from Schwabe to Wolf and an inappropriate individual correction factor used for Schwabe in the WSN; (3) This implies a major inhomogeneity in the WSN, which needs to be corrected by reducing its values by 20% before 1848; (4) The corrected WSN series is in good agreement with the GSN series. This study supports the earlier conclusions that the GSN series is more consistent and homogeneous in the earlier part than the WSN series.

While the presence of magnetic fields on low-mass stars is attributed to a dynamo process essentially driven by convective motions, the existence of magnetic fields on intermediate-mass stars has very likely other reasons. Presuming that the fields we see are nearly constant in time, the paper focuses on the generation of stable magnetic configurations at the early stages of stellar evolution. The convective processing of an initial magnetic field during the pre-main-sequence phase is studied in a very simple model star. Azimuthal magnetic fields are found to be typical remnants in the upcoming radiative envelope after the convection has receded.

Samuel Heinrich Schwabe made 8486 drawings of the solar disk with sunspots in the period from November 5, 1825 to December 29, 1867. We have measured sunspot sizes and heliographic positions on digitized images of these drawings. A total of about 135,000 measurements of individual sunspots are available in a data base. Positions are accurate to about 5% of the solar radius or to about three degrees in heliographic coordinates in the solar disk center. Sizes were given in 12 classes as estimated visually with circular cursor shapes on the screen. Most of the drawings show a coordinate grid aligned with the celestial coordinate system. A subset of 1168 drawings have no indication of their orientation. We have used a Bayesian estimator to infer the orientations of the drawings as well as the average heliographic spot positions from a chain of drawings of several days, using the rotation profile of the present Sun. The data base also includes all information available from Schwabe on spotless days.

Much of our knowledge about the solar dynamo is based on sunspot observations. It is thus desirable to extend the set of positional and morphological data of sunspots into the past. Gustav Sp\"orer observed in Germany from Anklam (1861-1873) and Potsdam (1874-1894). He left detailed prints of sunspot groups, which we digitized and processed to mitigate artifacts left in the print by the passage of time. After careful geometrical correction, the sunspot data are now available as synoptic charts for almost 450 solar rotation periods. Individual sunspot positions can thus be precisely determined and spot areas can be accurately measured using morphological image processing techniques. These methods also allow us to determine tilt angles of active regions (Joy's law) and to assess the complexity of an active region.

The frequent presence of weak magnetic fields on the surface of spotted late-B stars with HgMn peculiarity in binary systems has been controversial during the two last decades. We re-analyse available spectropolarimetric material by applying the moment technique on spectral lines of inhomogeneously distributed elements separately. Furthermore, we present new determinations of the mean longitudinal magnetic field for the HgMn star HD65949 and the hotter analog of HgMn stars, the PGa star HD19400, using FORS2 installed at the VLT. We also give new measurements of the eclipsing system AR Aur with a primary star of HgMn peculiarity which were obtained with the SOFIN spectropolarimeter installed at the Nordic Optical Telescope. We downloaded from the ESO archive the publically available HARPS spectra for eight HgMn stars and one normal and one superficially normal B-type star obtained in 2010. The application of the moment technique to the HARPS and SOFIN spectra allowed us to study the presence of the longitudinal magnetic field, the crossover effect, and quadratic magnetic fields. Results for the HgMn star HD65949 and the PGa star HD19400 are based on a linear regression analysis of low-resolution spectra obtained with FORS2 in spectropolarimetric mode. Our measurements of the magnetic field with the moment technique using spectral lines of several elements separately reveal the presence of a weak longitudinal magnetic field, a quadratic magnetic field, and the crossover effect on the surface of several HgMn stars as well as normal and superficially normal B-type stars. Furthermore, our analysis suggests the existence of intriguing correlations between the strength of the magnetic field, abundance anomalies, and binary properties.

The sunspot drawings of Johann Staudacher of 1749--1799 were used to determine the solar differential rotation in that period. These drawings of the full disk lack any indication of their orientation. We used a Bayesian estimator to obtain the position angles of the drawings, the corresponding heliographic spot positions, a time offset between the drawings and the differential rotation parameter \delta\Omega, assuming the equatorial rotation period is the same as today. The drawings are grouped in pairs, and the resulting marginal distributions for \delta\Omega were multiplied. We obtain \delta\Omega=-0.048 \pm 0.025 d^-1 (-2.75^o/d) for the entire period. There is no significant difference to the value of the present Sun. We find an (insignificant) indication for a change of \delta\Omega throughout the observing period from strong differential rotation, \delta\Omega\approx -0.07 d^-1, to weaker differential rotation, \delta\Omega\approx-0.04 d^-1.

We investigate the stability of the Hall-MHD system and determine its importance for neutron stars at their birth, when they still consist of differentially rotating plasma permeated by extremely strong magnetic fields. We solve the linearised Hall-MHD equations in a spherical shell threaded by a homogeneous magnetic field. With the fluid/flow coupling and the Hall effect included, the magnetorotational instability and the Hall effect are both acting together. Results differ for magnetic fields aligned with the rotation axis and anti-parallel magnetic fields. For a positive alignment of the magnetic field the instability grows on a rotational time-scale for any sufficiently large magnetic Reynolds number. Even the magnetic fields which are stable against the MRI due to the magnetic diffusion are now susceptible to the shear-Hall instability. In contrast, the negative alignment places strong restrictions on the growth and the magnitude of the fields, hindering the effectiveness of the Hall-MRI. While non-axisymmetric modes of the MRI can be suppressed by strong enough rotation, there is no such restriction when the Hall effect is present. The implications for the magnitude and the topology of the magnetic field of a young neutron star may be significant.

A method for estimating the true meteor rate \lambda\ from a small number of observed meteors n is derived. We employ Bayesian inference with a Poissonian likelihood function. We discuss the choice of a suitable prior and propose the adoption of Jeffreys prior, P(\lambda)=\lambda^{-0.5}, which yields an expectation value E(\lambda) = n+0.5 for any n \geq 0. We update the ZHR meteor activity formula accordingly, and explain how 68%- and 95%-confidence intervals can be computed.

The discovery of exotic abundances, chemical inhomogeneities, and weak magnetic fields on the surface of late B-type primaries in spectroscopic binaries has important implications not only for our understanding of the formation mechanisms of stars with Hg and Mn peculiarities themselves, but also for the general understanding of B-type star formation in binary systems. The origin of the abundance anomalies observed in late B-type stars with HgMn peculiarity is still poorly understood. The connection between HgMn peculiarity and membership in binary and multiple systems is supported by our observations during the last decade. The important result achieved in our studies of a large sample of HgMn stars is the finding that most HgMn stars exhibit spectral variability of various chemical elements, proving that the presence of an inhomogeneous distribution on the surface of these stars is a rather common characteristic and not a rare phenomenon. Further, in the studied systems, we found that all components are chemically peculiar with different abundance patterns. Generally, He and Si variable Bp stars possess large-scale organised magnetic fields that in many cases appear to occur essentially in the form of a single large dipole located close to the centre of the star. The presence of magnetic fields in the atmospheres of HgMn stars has been demonstrated in several studies. In addition to the measurements of longitudinal and quadratic magnetic fields, this work also showed evidence for a relative magnetic intensification of Fe II lines produced by different magnetic desaturations induced by different Zeeman-split components.

A long time-series of sunspot observations is preserved from Samuel Heinrich Schwabe who made notes and drawings of sunspots from 1825-1867. Schwabe's observing records are preserved in the manuscript archives of the Royal Astronomical Society, London. The drawings have now been digitized for future measurements of sunspot positions and sizes. The present work gives an inventory and evaluation of the images obtained from the log books of Schwabe. The total number of full-disk drawings of the sun with spots is 8486, the number of additional verbal reports on sunspots is 3699. There are also 31 reports about possible aurorae.

It is shown that the magnetic current-driven (`kink-type') instability produces flow and field patterns with helicity and even with \alpha-effect but only if the magnetic background field possesses non-vanishing current helicity \bar{\vec{B}}\cdot curl \bar{\vec{B}} by itself. Fields with positive large-scale current helicity lead to negative small-scale kinetic helicity. The resulting \alpha-effect is positive. These results are very strict for cylindric setups without <i>z/I>-dependence of the background fields. The sign rules also hold for the more complicated cases in spheres where the toroidal fields are the result of the action of differential rotation (induced from fossil poloidal fields) at least for the case that the global rotation is switched off after the onset of the instability.

Ap star magnetism is often attributed to fossil magnetic fields which have not changed much since the pre-main-sequence epoch of the stars. Stable magnetic field configurations are known which could persist probably for the entire main-sequence life of the star, but they may not show the complexity and diversity exhibited by the Ap stars observed. We suggest that the Ap star magnetism is not a result of stable configurations, but is the result of an instability based on strong toroidal magnetic fields buried in the stars. The highly nonaxisymmetric remainders of the instability are reminiscent of the diversity of fields seen on Ap stars. The strengths of these remnant magnetic fields is actually between a few per cent up to considerable fractions of the internal toroidal field; this means field strengths of the order of kGauss being compatible with what is observed. The magnetic fields emerge at the surface rather quickly; rough estimates deliver time-scales of the order of a few years. Since rotation stabilizes the instability, normal A stars may still host considerable, invisible toroidal magnetic fields.

This study is concerned with the early evolution of magnetic fields and differential rotation of intermediate-mass stars which may evolve into Ap stars. We report on simulations of the interplay of differential rotation and magnetic fields, the stability limits and non-linear evolution of such configurations, and the prospects of dynamo action from the unstable cases. The axisymmetric problem delivers a balance between field amplification and back-reaction of the magnetic field on the differential rotation. The non-axisymmetric case involves also the Tayler instability of the amplified toroidal fields. We consider limits for field amplification and apply these to young A stars. Apart from its application to Ap stars, the instability is scrutinized for the fundamental possibility of a dynamo. We are not looking for a dynamo as an explanation for the Ap star phenomenon. The kinetic helicity is concentrated near the tangent cylinder of the inner sphere of the computational domain and is negative in the northern hemisphere. This appears to be a ubiquitous effect not special to the Tayler instability. The latter is actually connected with a positive current helicity in the bulk of the spherical shell giving rise to a small, but non-vanishing alpha-effect in non-linear evolution of the instability.

The original sunspot observations by Heinrich Samuel Schwabe of 1825-1867 were digitized and a first subset of spots was measured. In this initial project, we determined more than 14 000 sunspot positions and areas comprising about 11% of the total amount of spots available from that period. The resulting butterfly diagram has a typical appearance, but with evident north-south asymmetries.

We investigate the phase difference of the sunspot cycles in the two hemispheres and compare it with the latitudinal sunspot distribution. If the north-south phase difference exhibits a long-term tendency, it should not be regarded as a stochastic phenomenon. We use datasets of historical sunspot records and drawings made by Staudacher, Hamilton, Gimingham, Carrington, Spouml;rer, and Greenwich observers, as well as the sunspot activity during the Maunder minimum reconstructed by Ribes and Nesme-Ribes. We employ cross-recurrence plots to analyse north-south phase differences. We show that during the last 300 years, the persistence of phase-leading in one of the hemispheres exhibits a secular variation. Changes from one hemisphere to the other leading in phase were registered near 1928 and 1968 as well as two historical ones near 1783 and 1875. A long-term anticorrelation between the hemispheric phase differences in the sunspot cycles and the latitudinal distribution of sunspots was traced since 1750.

We obtained high-resolution, high signal-to-noise UVES and a few lower quality HARPS spectra revealing the presence of resolved magnetically split lines. HD101412 is the first Herbig Ae star for which the rotational Doppler effect was found to be small in comparison to the magnetic splitting. The measured mean magnetic field modulus varies from 2.5 to 3.5kG, while the mean quadratic field was found to vary in the range of 3.5 to 4.8kG. To determine the period of variations, we used radial velocity, equivalent width, line width, and line asymmetry measurements of variable spectral lines of several elements, as well as magnetic field measurements. The most pronounced variability was detected for spectral lines of He I and the iron peak elements, whereas the spectral lines of CNO elements are only slightly variable. From spectral variations and magnetic field measurements we derived a potential rotation period P_rot=13.86d, which has to be proven in future studies with a larger number of observations. It is the first time that the presence of element spots is detected on the surface of a Herbig Ae/Be star. Our previous study of Herbig Ae stars revealed a trend towards stronger magnetic fields for younger Herbig Ae stars, confirmed by statistical tests. This is in contrast to a few other (non-statistical) studies claiming that magnetic Herbig Ae stars are progenitors of the magnetic Ap stars. New developments in MHD theory show that the measured magnetic field strengths are compatible with a current-driven instability of toroidal fields generated by differential rotation in the stellar interior. This explanation for magnetic intermediate-mass stars could be an alternative to a frozen-in fossil field.

Because of the lack of reliable sunspot observation, the quality of sunspot number series is poor in the late 18th century, leading to the abnormally long solar cycle (1784--1799) before the Dalton minimum. Using the newly recovered solar drawings by the 18--19th century observers Staudacher and Hamilton, we construct the solar butterfly diagram, i.e. the latitudinal distribution of sunspots in the 1790's. The sudden, systematic occurrence of sunspots at high solar latitudes in 1793--1796 unambiguously shows that a new cycle started in 1793, which was lost in traditional Wolf's sunspot series. This finally confirms the existence of the lost cycle that has been proposed earlier, thus resolving an old mystery. This letter brings the attention of the scientific community to the need of revising the sunspot series in the 18th century. The presence of a new short, asymmetric cycle implies changes and constraints to sunspot cycle statistics, solar activity predictions, solar dynamo theories as well as for solar-terrestrial relations.

Digitized images of the drawings by J.C. Staudacher were used to determine sunspot positions for the period of 1749-1796. From the entire set of drawings, 6285 sunspot positions were obtained for a total of 999 days. Various methods have been applied to find the orientation of the solar disk which is not given for the vast majority of the drawings by Staudacher. Heliographic latitudes and longitudes in the Carrington rotation frame were determined. The resulting butterfly diagram shows a highly populated equator during the first two cycles (Cycles 0 and 1 in the usual counting since 1749). An intermediate period is Cycle 2, whereas Cycles 3 and 4 show a typical butterfly shape. A tentative explanation may be the transient dominance of a quadrupolar magnetic field during the first two cycles.

A variety of magnetohydrodynamic mechanisms that may play a role in magnetic, chemically peculiar (mCP) stars is reviewed. These involve dynamo mechanisms in laminar flows as well as turbulent environments, and magnetic instabilities of poloidal and toroidal fields as well as combinations of the two. While the proto-stellar phase makes the survival of primordial fields difficult, the variety of magnetic field configurations on mCP stars may be an indication for that they are instability remnants, but there is no process which is clearly superior in explaining the strong fields.

The important result achieved in our recent study of a large sample of HgMn stars using UVES at the VLT and FEROS at the ESO 2.2m telescope is the finding that most HgMn stars exhibit spectral variability of various chemical elements, proving that the presence of an inhomogeneous distribution on the surface of these stars is a rather common characteristics and not a rare phenomenon.

The three-dimensional, hydrodynamic stability of the solar tachocline is investigated based on a rotation profile as a function of both latitude and radius. By varying the amplitude of the latitudinal differential rotation, we find linear stability limits at various Reynolds numbers by numerical computations. We repeated the computations with different latitudinal and radial dependences of the angular velocity. The stability limits are all higher than those previously found from two-dimensional approximations and higher than the shear expected in the Sun. It is concluded that any part of the tachocline which is radiative is hydrodynamically stable against small perturbations.

There are successful approaches to explain the formation of the tachocline by a poloidal magnetic field in the solar core. We present here the first MHD simulations of the solar tachocline which self-consistently include the meridional circulation. We show that the meridional flow significantly changes the shape and the characteristics of the tachocline. We find that after the inclusion of the meridional circulation, a tachocline can be formed even when the poloidal field lines are crossing the boundary between the radiative zone and the convection zone. We also discuss the effects of the magnetic Prandtl number as well as of the magnetic Reynolds number on the properties of the tachocline. The tachocline is much thinner for higher magnetic Reynolds numbers and/or lower magnetic Prandtl numbers. We expect that a poloidal magnetic seed field of around 1 G will be sufficient to produce the tachocline of the Sun. However, the model requires the initial magnetic field to be in a narrow range for satisfying tachocline solutions. The simulations including a stable temperature gradient produce a shallower as well as slower meridional circulation than the ones without it, as desired from the Lithium abundance at the solar surface.

Stars of spectral classes A and late B are almost entirely radiative. CP stars are a slowly rotating subgroup of these stars. It is possible that they possessed long-lived accretion disks in their T Tauri phase. Magnetic coupling of disk and star leads to rotational braking at the surface of the star. Microscopic viscosities are extremely small and will not be able to reduce the rotation rate of the core of the star. We investigate the question whether magneto-rotational instability can provide turbulent angular momentum transport. We illuminate the question whether or not differential rotation is present in CP stars. Numerical MHD simulations of thick stellar shells are performed. An initial differential rotation law is subject to the influence of a magnetic field. The configuration gives indeed rise to magneto-rotational instability. The emerging flows and magnetic fields transport efficiently angular momentum outwards. Weak dependence on the magnetic Prandtl number (~0.01 in stars) is found from the simulations. Since the estimated time-scale of decay of differential rotation is 10^7-10^8 yr and comparable to the life-time of A stars, we find the braking of the core to be an ongoing process in many CP stars. The evolution of the surface rotation of CP stars with age will be an observational challenge and of much value for verifying the simulations.

The hydrodynamic stability of accretion disks is considered. The particular question is whether the combined action of a (stable) vertical density stratification and a (stable) radial differential rotation gives rise to a new instability for nonaxisymmetric modes of disturbances. The existence of such an instability is not suggested by the well-known Solberg-Hoiland criterion. It is also not suggested by a local analysis for disturbances in general stratifications of entropy and angular momentum which is presented in our Section 2 confirming the results of the Solberg-Hoiland criterion also for nonaxisymmetric modes within the frame of ideal hydrodynamics but only in the frame of a short-wave approximation for small m. As a necessary condition for stability we find that only conservative external forces are allowed to influence the stable disk. As magnetic forces are never conservative, linear disk instabilities should only exist in the magnetohydrodynamical regime which indeed contains the magnetorotational instability as a much-promising candidate. To overcome some of the used approximations in a numerical approach,the equations of the compressible adiabatic hydrodynamics are integrated imposing initial nonaxisymmetric velocity perturbations with m=1 to m=200. Only solutions with decaying kinetic energy are found. The system always settles in a vertical equilibrium stratification according to pressure balance with the gravitational potential of the central object. keywords: accretion disks -- hydrodynamic instabilities -- turbulence

The theory of the solar/stellar activity cycles is presented, based on the mean-field concept in magnetohydrodynamics. A new approach to the formulation of the electromotive force and the theory of differential rotation and meridional circulation is described. Dynamo cycles in the overshoot layer and distributed dynamos are compared, with the latter including the influence of meridional flow. The overshoot layer dynamo reproduces the solar cycle periods and the butterfly diagram only if alpha=0 in the convection zone (CZ). The distributed dynamo including meridional flows shows the observed butterfly diagram even with a positive dynamo-alpha in CZ. The nonlinear feedback of strong magnetic fields on differential rotation leads to grand minima in the cyclic activity similar to those observed. Our 2D model contains the large- and small-scale feedback of magnetic fields on diff. rotation and induction in a mean-field formulation (Lambda-, alpha-quenching). Grand minima may also occur if a dynamo occasionally falls below its critical eigenvalue. We never found any indication that such an on-off dynamo collapses by this effect after being excited. The full quenching of turbulence by strong magnetic fields as reduced induction (alpha) and reduced turbulent diffusivity (eta_T) is studied in 1D. We find a stronger dependence of cycle period on dynamo number compared with a pure alpha-quenching model giving a very weak cycle period dependence. Also the temporal fluctuations of alpha and eta_T from a random-vortex simulation were applied to a dynamo. Then the low `quality' of the solar cycle can be explained with a small number of giant cells as dynamo-active turbulence. The transition from almost regular magnetic oscillations (many vortices) to a more or less chaotic time series (very few vortices) is shown.

The possibility of non-helical large scale dynamo action is investigated using three-dimensional simulations of global accretion discs as well as idealized local simulations without rotation and only shear. Particular emphasis is placed on a certain correlation between vorticity and azimuthal velocity gradient which has been predicted to drive large scale dynamo action, independent of the presence or absence of kinetic helicity. In the global disc simulations two types of behaviours are found: those which do show this type of velocity correlation and those which do not. The former ones are typically also the cases where the resistivity is larger. The latter ones show signs typical of dynamo action based on the usual helicity effect. In the idealized simulations without rotation and just shear the above correlation is found to be particularly strong. In both cases there is, as expected, a systematic flux of magnetic helicity through the midplane. However, very little helicity flux leaves the domain through the top and bottom boundaries. The idealized simulations reveal that much of this systematic flux comes from the rotational component of the helicity flux and does not contribute to its divergence.

We perform global three-dimensional simulations of accretion disks integrating the compressible, non-viscous, but diffusive MHD equations. The disk is supposed to be isothermal. We make use of the ZEUS-3D code integrating the MHD equations and added magnetic diffusivity. We measure the efficiency of the angular-momentum transport. Various model simulations delivered transport parameters of alpha_SS=0.01 to 0.05 which are consistent with several local numerical investigations. Two of the models reach a highly turbulent state at which alpha_SS is of the order of 0.1. After a certain stage of saturating of the turbulence, Reynolds stress is found to be negative (inward transport) in many of the models, whereas Maxwell stresses dominate and deliver a positive (outward) total transport. Several of the models yield strongly fluctuating Reynolds stresses, while Maxwell stresses are smooth and always transport outwards. Dynamo action is found in the accretion disk simulations. A positive dynamo-alpha is indicated in the northern hemisphere of the most prominent run, coming along with negative kinetic and current helicities (all having the opposite sign on the southern side). The dipolar structure of the magnetic field is maintained throughout the simulations, although indication for a decay of antisymmetry is found. The simulations covered relatively thick disks, and results of thin-disk dynamo models showing quadrupolar fields may not be compatible with the results presented here.