A Remarkable Recent Transition in the Solar Dynamo

We summarize the major aspects of the remarkable, fairly long lasting period (∼2005\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sim 2005$\end{document} to ∼2010\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sim 2010$\end{document}) of low solar activity, that we will call the Transition. It is the transitional stage between the Grand Maximum of the 20th century and a forthcoming (most probably Regular) episode of solar activity. The various kinds of activity in the functioning of the equatorial components of the solar dynamo before and during the Transition are summarized. While the behavior of unipolar magnetic regions and their rest-latitudes already gave very early indications – mid 20th century – of the forthcoming Transition, more such indications became available around 1995 and the main part of it occurred between 2005 and 2010. Some of the inferences are discussed. We submit the hypothesis that the solar tachocline undergoes pulsations and we present some helioseismic evidences. In that scenario we find that its equatorial part has moved downward over a fairly small semi-amplitude (∼0.03\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sim 0.03$\end{document} solar radii) during the time of the Transition. There are several indications, apart from this ‘pulsation’, that the tachocline may even be pulsating with still smaller amplitudes in more modes. We speculate about the physical mechanism(s).


Introduction
When, around the year 2005, the past Schwabe cycle (#23) in solar activity had come down to its minimum, the expected rise to the maximum of cycle #24 appeared to be considerably delayed. Contrary to initial expectations, the usual solar activity, such as the Total Solar Irradiation (TSI), sunspots, Coronal Mass Ejections, flares . . . only started to increase in intensity, number and degree of activity in the period between 2008 and 2010.
This relatively long-lasting minimum had several unusual aspects. Earlier, Duhau and de Jager (2008) distinguished three kinds of relatively large episodes of solar activity, each of them containing about a dozen or more Schwabe cycles, viz. Grand Maxima, Grand Minima and episodes of Regular variations. Usoskin et al. (2014), who studied solar activity over a period of three millennia, confirmed this subdivision. The reality of the recent (20th century) Grand Maximum is confirmed in critical studies by Solanki et al. (2004), Lockwood et al. (2009Lockwood et al. ( , 2010, Lockwood (2012) and Clette et al. (2014), although, after the introduction of the new system of sunspot counting (Clette et al. 2014;Svalgaard and Schatten 2016), the maximum is less extreme than was thought before. Clette et al. (2014) call it the 'Modern Maximum'. In order to avoid confusion we will continue to call it the '20th century Grand Maximum' or briefly the 'Grand Maximum'.
The Grand Maximum started in 1924 (Duhau and Chen 2002; while Lockwood et al. (2009) place the beginning in 1920. For the end of it  find 2008 (cf. our Table 1 later in this paper, and the corresponding text), while Lockwood et al. (2009) The period that we are discussing here is the time interval between the past Grand Maximum and a forthcoming Episode of which  and de Jager and Duhau (2012) give arguments that it will be of the Regular type. In this paper we will call this transitional period the Transition. Feynman and Ruzmaikin (2014) give examples of more such transitions that may have existed in the past.
Other authors have already mentioned this Transition. We refer to papers and reviews by Russell et al. (2010), Lockwood et al. (2010), Lockwood (2012), Nandy et al. (2011), Solanki and Krivova (2011), de Jager and Duhau (2011), . Quoting Russell et al. (2010): "Previous solar minima had occurred in 1996, 1986, 1976 and 1966, so  Uncertainty remained at that time, as appears from a quotation by Weiss (2010): Is this just an abnormal fluctuation or are we about to experience a Maunder-like Grand Minimum? The answer can be given, based on a look a the sunspot curve. This recent extended minimum is the longest of the past century. More of them have occurred earlier in the past (some of them listed by Feynman and Ruzmaikin 2014). Examples of sunspot cycles with extreme minima are cycle #13 (maximum in 1894; minimum sunspot number 1.7), and cycle #14 (maximum in 1906; minimum sunspot number 1.5). We compare these with the last fully observed cycle #23: its minimum sunspot number was 1.7 -just comparable. We note that these sunspot numbers are given according to the 'classical' counting system -cf. the remark at the end of this section. The difference between these extreme minima and the present Transition phase is that this last one is the first one that could be observed by modern observational techniques. This period of solar relative inactivity had several aspects: these are partly related to the polar and partly to the equatorial part of solar activity, the two main solar magnetic regions. In this paper we first summarize and analyse existing observations. Thereafter these observations are synthesized into a scenario that links them to the two magnetically active solar regions and we examine the implications for our current understanding of part of the functioning of the solar dynamo. We will forward a new suggestion, viz. that the solar tachocline is pulsating. Quantitatively, we find that these 'pulsations' need to have a remarkably small amplitude in order to explain the observations. Finally, we present some reflections on the possible mechanism.

Remark on Sunspot Numbers and Numbering System
Around the beginning of this century the Zürich sunspot numbering system introduced by Wolf and developed further by Waldmeier and successors (The Wolf Sunspot Numbers, generally called International Sunspot Numbers, ISN) was replaced by a new system of sunspot numbering (the Modified International Sunspot Numbers (MISN; Clette et al. 2014), while another new counting system is called the Group Sunspot Numbers (GSN; Hoyt and Schatten 1998;Hathaway et al. 2002;Svalgaard and Schatten 2016). The main reason for the change to the MISN was to overcome the 'Waldmeier discontinuity' (also described by Aparicio et al. 2012), while the introduction of the Group Sunspot Numbers was bases on the -likely -assumption that sunspot groups represent solar equatorial magnetic variability better than the Zürich system.
After these modifications it is inevitable that in this paper -which has partly the character of a review paper -diagrams or formulas appear that are based on one of the different systems of sunspot numbering. We realize that this is a confusing situation and try, in each particular case, to inform the reader about the system of numbering that we are dealing with.

The Role of the Polar and Equatorial Magnetic Regions in Controlling Solar Activity Variations; Proxies for the Main Variables
The fact that the solar surface presents two magnetically active regions is often overlooked, particularly by climatologists who study sun-climate relationships. In their research they nearly invariably restrict themselves to searching for the relation between climate and the magnetic activity in the sunspot region, where the latter is characterized by the number of sunspots on the disk. But the situation that the polar and equatorial areas of the sun are both dominated by magnetic fields and structures with approximately equally strong total magnetic fluxes is already known since the time that  introduced technical innovations in their spectral observations of magnetic fields that permitted them to observe fields as small as 0.5 Gauss. Babcock and Livingston (1958) studied the polar magnetic fields. According to their measurements the northern field had at that time an average strength of 0.2 G, while the southern average field strength was 0.6 G. They found that, on the average, the time of maximum of the polar field was delayed by 3 years after sunspot minimum. Earlier work by Grotrian and Künzel (1950) had already given values for the total magnetic flux of the equatorial (sunspot) region and somewhat later work (review by de Jager 1959, p. 340 and following) showed that the polar and equatorial fluxes are approximately comparable in magnitude. The same result was obtained later by Makarov et al. (2004): the total flux from the 'large unipolar fields' (of order of a few gauss) is according to them by about a factor 2 larger than the flux from the polar faculae (about a few hundred Gauss each) and the latter is about a factor 2 larger than the total flux due to the sunspots. However, for the corresponding local magnetic energies the order is the opposite. In connection to this, we should add that the main issue is not how much flux the poles have in comparison with active lower latitudes, but the fact that they seem to be the clear precursors of the subsequent Schwabe cycle.
In view of the comparable total fluxes of the polar and equatorial magnetic fields it is necessary to consider both elements in discussions of solar variability. The role of the large unipolar fields is another item that has, so far, hardly been studied (cf. Sect. 6 of this paper).
The polar regions have only been subject of regular observations during the past half century and for a study of the long-term variation of polar activity the introduction of a proxy is necessary. The same applies to the equatorial region. The magnetic fields of the sunspots have been measured since the early part of the 20th century but not systematically so. Systematic observations of the facular magnetic field strengths around spots are known only since the 2nd part of that century. Therefore, we need proxies in order to describe the time dependence of these two variables.
The sunspot number is related to the magnetic flux contained in the (equatorial) sunspot belt. That activity shows itself in different ways: sunspots, faculae, flares, Coronal Mass Ejections, just to mention the most important components. Since that ensemble of activities is only well known since the 2nd half of the 20th century, its values are usually simulated by either the International Sunspot Number or the Group Sunspot Number. Both are generally taken as suitable proxies. A better proxy for the total magnetic energy contained in active regions would have been the integrated magnetic flux of all spots present, but the latter quantity has only been measured during part of the 20th century. Therefore, in the absence of a better proxy for the long-term total equatorial activity, one generally chooses R max , the maximum international sunspot numbers, or else, the maximum value of the recently introduced Group Number for subsequent Schwabe cycles. The R max data are known from 1610 onward.
A direct measure for the total polar activity would be the sum of the magnetic fields of the polar faculae. There have been measured since the beginning of the 20th century (Sheeley 1966(Sheeley , 1976(Sheeley , 1991(Sheeley , 2008 and their calibrated values over the past 100 years are clearly correlated with the heliospheric magnetic fields. The measured values of the polar faculae are also in good agreement with both the polar fields and the polar flux estimates. This makes these data excellent proxies for studying the polar magnetic fields (Muñoz-Jaramillo et al. 2012, 2013a, 2013b. A drawback is that data of the polar facular flux are only known since the beginning of the 20th century. An alternative proxy for the total polar magnetic field strength is therefore the commonly used value of aa min , the minimum value of the geomagnetic aa index for subsequent Schwabe cycles (Russell 1975;Russell and Mulligan 1995;Hathaway et al. 1999;Duhau and Chen 2002) The aa min data are known from 1870 onward (Mayaud 1972) and the series has been extended down to 1844 by Nevalinna and Kataja (1993). The aa min data from 1856 onward were homogenized by M. Lockwood (priv. comm.) and published by de Jager and Duhau (2012). Still earlier aa min data (before 1844) do exist in the literature (cf. Nagovitsyn 2005) but these are based on an assumed correlation between R max and aa min , while an essential part of this paper is that such a correlation should not be assumed and to the best has to be studied. Therefore such extrapolated data will not be used in this paper. (An exception is made for Figs. 12a and part of 12b, but we have to realize that for these diagrams the parts from before 1844 are unreliable.) The fact that aa min is a suitable proxy for the polar magnetic fields is also shown by diagrams such as Fig. 1 which presents the relation between aa min and the total polar magnetic flux DM max . The correlation between these two quantities (aa min and polar flux) is excellent. Fig. 1 The relation between aa min and the total polar magnetic flux DM max . The aa min data are homogenized by M. Lockwood from 1856 onward and published by de Jager and Duhau (2012) The aa index as such carries information on both the toroidal and the poloidal components of the magnetic field. The first is related to the (equatorial) activity centers, the latter rather to the corotating interaction regions in the solar wind. When the aa index reaches its minimum value the former contribution is very small. For that reason the use of the aa min proxy is advocated.

The Main Aspects of Solar Activity During the Transition
Turning to the theme of this paper, the Transition of 2005-2010, we reproduce as Fig. 2a a diagram due to Lockwood, published as an addendum to a review by de Jager (2012) after a study week at the Berne International Space Science Institute (Bengtsson et al. 2012). It shows nine different aspects of solar variability. Arrows point to the time of minimum or maximum, where the maxima refer to those variables that reach their maximum values at times of minimum solar activity, e.g. cosmic ray maximum appears at times of minimal solar wind flux.
As appears from the text in the frame of the diagram there are, roughly, two periods of solar extreme activity (in the sense as described above) during the period of the Transition. One of them is between 2008.5 and 2009.5 and a second one is between 2009.5 and 2010.5. The first time of extreme values is reached by the sunspot number, the middle UV spectral irradiance (SSI) and the MgII spectral lines in the near UV part of the spectrum. These variables are all directly related to the activity regions in the sunspot belt and to their medium-energy radiation fluxes. They are, hence, in one way or the other related to the Total Solar Irradiance (TSI), and in turn to the sunspot number R.
The late group of extreme values, those around 2010, is mainly related to the corpuscular solar emission of which the open solar flux (Lockwood et al. 2009;Lockwood 2013) is an important component. The magnetized part of the corpuscular radiation is related to the cosmic ray flux (cf. Sect. 5), because the latter is modulated by the interplanetary magnetized solar plasma. Maximum magnetized solar plasma corresponds to minimum cosmic ray flux and reversely so.
We compare these observations with those of the X-ray emitting jets emanating from coronal holes (Savcheva et al. 2009), while realizing that the coronal holes are chiefly related to the polar fields. The jets show a motion transverse to the jet outflow direction. This motion reflects the magnetic polarity orientation of the ephemeral active region at the base of the jets. The direction of the motion changed between 2006-2007 and 2008-2009. Remarkably enough, the preferred motion in the first period reflected the active regions of the coming sunspot cycle (#24) while those of the latter period reflected the cycle #25. These observations confirm that new cycles of solar activity start at high latitudes.
The data from Fig. 2a find further support in Fig. 2b, which shows the time variation of several variables that are strongly associated with the sunspot number, hence with the equatorial part of solar activity. The data from this latter diagram have been derived from several sources, that are listed in a diagram by Balogh et al. (2014). They show that all these sunspot-associated variables reach their minimum values in 2008.5-2009.0 with an uncertainty of ±0.5 year.
To further illustrate the variations relative to the sunspot belt activity we reproduce, as Fig. 3, a diagram due to Yeo et al. (2014). It gives an overview of the variation of the TSI from 1972 through 2011, hence including the years of the Transition. Among other things it also shows the increased accuracy of measurements with time. In the earlier part of the measurements, around 1986 the spreading is large, e.g. the ACRIM data being lower than all others, while in 2009 they are all nearly equal. The average value of the four curves confirms Lockwood's time of minimum TSI and Fig. 3 also shows that the (average of the) minimum TSI value of 2008.5 is the lowest of the three cycles shown. The difference is insignificantly small: it falls within the error bounds. Hence, this information cannot be used for the present study but it seems worth to keep it in mind for later investigations (Yeo et al. 2015). Anyway, realizing that the TSI reflects the concentration of magnetic elements in the chromosphere and specifically the magnetic flux of the facular fields around the sunspots, it shows that during the Transition the magnetic activity in the sunspot belt was very low.
De  and de Jager and Nieuwenhuijzen (2013) have found that 34 % or 24 % respectively, of the sun-related long-term average terrestrial ground temperature variations (as measured from 1610 till 1950 and till 1912, respectively) are proportional to variations of the polar magnetic flux. They suggested that these latter variations are linked to the variations in the solar corpuscular flux, most probably the flux emanating from the coronal holes. It is significant that during the century before the year 2000 the average interplanetary magnetic field had doubled (Solanki et al. 2000) but that the solar wind velocity decreased steadily since about 2000. It reached a minimum value after 2009 (Lockwood et al. 2009;reproduced as Fig. 11 by Russell et al. 2010). This is one of the bases of part of the conclusions summarized in our Figs. 2a and 2b.
This particle flux, being a major part of the open solar flux (cf. Lockwood 2013) is most active at times near sunspot minimum and is associated with the polar magnetic regions and particularly with the coronal holes. These, in turn, are related to the poloidal field of the solar dynamo. Therefore, these observations confirm the above suggestion that the late (∼ 2010) minimum activity should be associated with the poloidal (polar) part of the sun's magnetic field while, as we noted before, the earlier TSI minimum (of 2008.5) is related to the equatorial (toroidal) part.
In conclusion, the extreme values of solar activity variation in the early part of the Transition (∼ 2008.5) are related to variations in the equatorial magnetic fluxes, while those that regulate the later part (∼ 2009-2010) are connected to that in the polar magnetic fields.

Umbral Magnetic Fields
A prominent aspect of solar variability during and partly preceding the Transition is an observed relatively stronger variation of the average maximum sunspot umbral field strength  The total solar irradiance (measurements from four monitors). Between cycles #21 and #24 the average TSI went through a minimum that may have been deeper than the average value of most of that of the foregoing ones (Yeo et al. 2014). Vertical lines mark the time of sunspot cycle minima data were obtained by measuring the Zeeman splitting of this near-infrared Fe I line from spectra taken at the McMath-Pierce telescope. In spot umbrae this line is always completely split and it thus directly yields field strength and not flux (as given by other magnetographs). The spatial umbral resolution was 2.5 arcsec, to which seeing effects should be added.
In papers initiated by Livingston (2006, 2009), cf. also, , Watson et al. (2009Watson et al. ( , 2014, Livingston and Watson (2015), Rezaie et al. (2015), there was initial evidence that the average maximum umbral magnetic field strength decreased steadily from 1999-2005. This decrease is now understood to have been biased by a then observational interest in large spots. After 2005, and through to the present (2016), all spots on the visible disk were observed (in the allotted telescope time of a few days per month) and after 2010 the average approached a daily value of slightly above 2000 G, as shown in Fig. 4a. This result was confirmed in broad outline by Watson et al. (2009Watson et al. ( , 2014; cf. The initial belief was that the decrease to around 2000 G during the Transition was a unique property of the Transition, but that conclusion was reformulated on basis of observations published by  and ; cf. Fig. 4d. They found, from observations covering the past half-century, that the average maximum umbral field strength is virtually never constant, contrary to what had always implicitly been assumed by the solar community. An extended series of field flux measurements over the past five decades yields evidence that the cyclic diminishing effect is part of a periodic phenomenon: in successive Schwabe cycles the umbral maximum field flux goes up-and downwards, in agreement with the sunspot number. But during the Transition the decrease was more clearly evident. Figure 4d   "pores" whether or not we see penumbrae. Pores are treated herein the same way as sunspots. We should add that pores can be missed if the seeing is inadequate.
The data of Fig. 4a show a flat trend of about 2200 G after 2005. A comparison of these data with those of Fig. 4d may suggest a difference between the two series of measurements but that is hard to prove because the one series of measurements virtually ends where the other starts. In the too short period of overlap between Figs. 4a and 4d there are only small differences. That makes one wonder about the precise values in the years of overlap. It is only to be regretted that there is no better overlap between the two series of measurements.
We conclude that the image of Fig. 4d is confirmed in broad outline by Figs. 4a, 4c, and 4e.
Before going to an interpretation of the data we need to add a few lines on their delicate character. The data of Fig. 4c are based on observations in sunspot umbrae, approximately 5 to 6 days per month. One observation is made in each spot that was visible, at the darkest position in their umbrae. Therefore, as said above, these observations are not compatible with those shown in Figs. 4a, 4b and 4d. Figure 4c is interesting by itself in that it displays the day-to-day means with their overall average (the constant dashed line).
After these remarks we conclude: Based on available information, as summarized in Figs The intriguing question is how this will continue in coming years. Summarizing these conclusions: Eleven-year variations of the average umbral magnetic field values are apparently a common phenomenon but the minimum value was more pronounced and lasted longer during the recent Transition, which shows in Figs. 4 as an extended minimum.  give evidence that the Livingston-Penn effect is regulated by the relative contribution of small and large sunspots in the ensemble of sunspots, in such a way that the relative number of smaller spots was larger during the extended minimum. They stress that there is no monotonous centennial decline and that the umbral magnetic field measurements just exhibit cyclic variations.
This conclusion is confirmed in the various components of Fig. 4.  Potgieter (2008). Hence, the overall variation of the intensity of cosmic rays at the Earth depends on various aspects of solar, interplanetary and heliospheric activity. For details cf. the review by Potgieter (2013). The modeling of the global heliosphere with comprehensive MHD models indicates that the heliosphere is highly asymmetric in terms of a nose-tail and even a north-south (meridional) direction and that conditions in the interstellar medium could influence the extent of the heliosphere (e.g. Borrmann and Fichtner 2005), and therefore also cosmic rays over millennia. For a comprehensive review on several aspects of space climate related to cosmic rays, cf. Scherer et al. (2006) and for a comprehensive review of solar activity over millennia, Usoskin (2013). Cosmogenic radionuclides also provide records of the paleocosmic radiation that extend to more than 103 centuries into the past and are as such indicators of changes in space climate (e.g. . Studies of cosmic ray variations over the past 9400 years, and their implications regarding solar activity and the solar dynamo are also reported (e.g. ).

Cosmic Ray Observations
The long-term variation as actually observed by cosmic ray detectors on ground level (called neutron monitors) is shown in Fig. 5a. This detector counts all cosmic rays with energies above about 7 GeV. The plots are shown together with sunspot numbers and the tilt angle of the wavy heliospheric current sheet (HCS). The latter is an excellent proxy for solar activity from charged particle point of view and leads to the distinct profiles which characterize the 22-year cycle seen in the upper panel of Fig 5a. Note in Fig. 5a how on top of the larger 11-year cycle, the cosmic ray intensity typically peaks in so-called A < 0 solar magnetic field epochs (e.g. 1960s, 1980s, etc.) but being flattish during A > 0 epochs (1970s, 1990s, etc.). The dominant modulation process in this regard is particle drifts in addition to the conventional convection, diffusion and adiabatic energy changes. During A > 0 polarity cycles all positively particles drift toward the Earth mainly through the polar regions of the heliosphere while during A < 0, when they drift inward mainly through the equatorial region, they encounter the wavy HCS. Consequently, this strong relation to the changing HCS waviness (tilt angles) before and after a solar minimum shaped the cosmic ray time profiles into forming peaks. This process was repeated during the 2009 solar minimum (A < 0 cycle) but the time profile was less sharp as depicted in Fig. 5a, indicative that the last solar minimum was indeed different from previous A < 0 cycles.
Another consequence of drifts is charge-sign dependence in cosmic ray modulation because negative particles will drift differently towards the Earth. This is a main feature in the study of how galactic cosmic rays respond to what happens to the global magnetic field of the Sun; for details, see the review by Potgieter (2014). The full reversal of the polarity of the global solar magnetic field during periods of extreme solar activity is another important aspect of how the 22-year cycle is apparent in the heliosphere.
An important aspect for the present study is the rise of the cosmic ray flux to reach in December 2009 the highest maximum that have been observed since the beginning of the space research era. This is shown in Fig. 5b for the PAMELA observations (Adriani et al. 2013;. Instead of grouping with the previously observed A < 0 spectra (in red) the December 2009 spectrum was even higher than the A > 0 spectra (in blue) which, according to drift modulation models, should always be higher than A < 0 cycles at energies below 500 MeV. This was even measured at ground level as shown in the top panel of Fig. 5a, although at these higher energies the effect, described here, is diminished. As noted, the shape of the intensity-time profile was far less sharply peaked around 2009 than during previous A < 0 polarity cycles. This is a reflection of the exceptional transparency of the interplanetary medium and heliosphere for galactic cosmic rays at that time, because of the very low level of solar activity in 2009 and the preceding years. The date, December 2009, confirms the conclusion proton spectra observed over the past five decades during solar minimum activity. The A < 0 magnetic polarity cycles are shown in red and the A > 0 cycles in blue. The PAMELA spectrum for the end of 2009 is indicated by the star-symbols and is the highest cosmic ray spectrum observed since the beginning of the space research era. References to these observations as shown in the legend are listed by Strauss and Potgieter (2014) of Figs. 2a and 2b that the solar particle flux reached a minimum in 2009 ± 0.5. In addition to PAMELA observations, Mewaldt et al. (2010) reported that cosmic rays with high rigidity reached record setting intensities during this time, while Heber et al. (2009) andMcDonald et al. (2010) also noted this extraordinary event from space bound observations. It followed from these cosmic ray observations for this period that the delay between the time for minimum sunspot numbers and maximum cosmic ray intensities was at least three times longer than during previous even numbered solar cycles (see also e.g., Kane 2011; Aslam and Badruddin 2012;Mewaldt 2013).
It can also be noted from Fig. 5a that sunspot numbers (middle panel) have been very low during the present period of maximum solar activity, indicating that the Sun has indeed been quieter than during recent previous cycles, even during solar maximum periods and that cosmic rays are responding accordingly. The time-profile of the tilt angles, on the other hand, does not reveal a quieter Sun suggesting that although the Sun becomes overall quieter, with global solar wind turbulence significantly less than the average value, the solar magnetic field still follows the usual time behavior. This could be interpreted that even during an epoch such as the Maunder Minimum the polarity of the solar magnetic field had kept changing in a 22-year fashion.
However, the magnitude B of the solar magnetic field observed at the Earth decreased significantly from 2006 to 2009 to become the lowest observed so far. This has contributed significantly to the observed higher cosmic rays during this period because cosmic ray diffusion and drifts scale essentially as 1/B. Svalgaard and Cliver (2007) proposed that B at Earth has a 'floor' value of ≈ 46 microgauss in yearly averages, which should be obtained during solar minima. They attributed the floor to a constant baseline solar open flux. In 2008 to 2009, the notion of such a floor was undercut by annual averages of ∼ 40 microgauss. Cliver and Ling (2011) revised both the level and the concept of such a floor and concluded that B ∼ 28 microgauss could be a 'floor' value. They suggested that the source of the floor in B is the slow solar wind; the Earth was in slow solar-wind flows ∼ 70 % of the time during 2009. The decrease in the tilt angle was, in contrast, rather ordinary: it had not decreased as rapidly as the magnitude of the heliomagnetic field at Earth during this period, but eventually also reached a minimum value at the end of 2009.
The large decreasing and increasing steps over relatively short periods shown in the upper panel of Fig. 5a is indicative of large transients originating from the Sun and that, when propagating throughout the heliosphere, subsequently disturb the quiet solar minimum conditions and eventually lead to a minimum cosmic ray flux at the Earth during periods of extreme solar activity. In this regard a prominent role is played by Coronal Mass Ejections (CME) as indicated by Cliver and Ling (2001a). These CME's are mostly correlated to the equatorial solar activity but can cause, through merging and interaction, cosmic ray regions to extend towards higher heliolatitudes and throughout the heliosphere (see e.g. Potgieter 1993). A significant decrease in the number of CMEs was observed in -2009(Vourlidas et al. 2010) which means that these cosmic rays encountered fewer barriers into the inner heliosphere so that significantly fewer transients such as Forbush decreases had been observed during the present solar cycle.

Unipolar Magnetic Regions; Their Significance for Solar Variability and Their Relation to Coronal Holes
In this Section we consider the origin and development of the solar large unipolar magnetic regions, and their relation to the coronal holes, their Hale boundaries (Svalgaard and Wilcox 1976) and their rest latitudes (Makarov et al. 2004) during the past Grand Maximum of solar activity. Distinction should be made between apparent and real unipolar regions. There are many seemingly unipolar regions on the sun, that are so small that their unipolar character can be negated nor confirmed. Lamb et al. (2008) find that the majority of the particle flux of the quiet sun originates from areas of coalescence of previously existing but unresolved flux, that show up as arcsecond structures. They also estimate that their rate of coalescence, integrated over the sun, is 7 × 10 21 Maxwell/hr. They further suggest that "currently unresolved emergences may be at least as important as ephemeral regions emergences to the total overall flux budget".
Next to these quasi-unipolar regions there are the larger real unipolar regions. They have field strengths of a few gauss. Their latitude bounds are fairly similar to those of the active regions and they encompass between one third and one tenth of the solar surface. Their fundamental role in governing solar variability has been grossly underestimated for a long time and was only realized quite recently.
We describe their discovery and significance. With the increased observational possibilities of the fifties of the past century, , (cf. also Simpson et al. 1955) detected large areas on the sun that were characterized by only one magnetic polarity, the unipolar magnetic field regions. Since their average field strengths are only of the order of a few Gauss, their observations demand highly developed instrumental facilities, and that was the reason why their study progressed only slowly. Initially this even led to a wrong interpretation. Leighton (1964) suggested that they originate out of bipolar active regions when these are decaying and that idea found initial support in observations by Bumba and Howard (1965). These authors found that their latitudinal distribution coincided fairly well with that of the bipolar active regions. Further support of that suggestion was presented by Giovanelli (1982).
It was at about the same time that McIntosh and Wilson (1985), also guided by earlier work from Svalgaard and Wilcox (1976) found that the above scenario cannot be correct, if only for the reason that unipolar magnetic field regions do appear at the surface before the disappearance of spot groups. Later investigations gradually showed the fundamental importance of these regions for solar variability, particularly for the origin and developments of sunspot pairs, cf. Solanki (2003), Akasofu (2014Akasofu ( , 2015. The present views are represented by a schematic diagram published by Akasofu (2015) that is partly based on work by McIntosh (1981). It is reproduced here as Fig. 6. Unipolar magnetic regions occur on the solar surface in an alternating positive-negative sequence. These two parts are separated by a neutral line, the Hale boundary. The spots of the two opposed polarities of sunspot groups occur on the two sides of the Hale boundary. Single spots of the same polarity as the unipolar magnetic regions may occur inside the unipolar magnetic regions. The difference between the even-and odd-numbered solar cycles (the Hale law) and also the different behavior on either side of the equator are illustrated by the two parts of Fig. 6. A further comment is that only the latitudinal dependence is shown in the diagrams. With regard to this, the examination of a large number of magnetic maps, depending on the polarity of the dipole field (sunspot cycle) shows that one of the unipolar regions is often, but not always, connected to the polar field. Because the relationship between the polar and unipolar fields is badly known (but for the polar migration of unipolar field), the diagram does not include the polar field. Also for that reason Fig. 6 is a schematic representation.
The systematic equatorward drift of the region of occurrence of sunspots (the 'butterfly diagram') may find an explanation in a study by Nakagawa (1971). In his model largescale deep-seated whirling circulation patterns, oriented along the meridians, occur on either  Akasofu (2015) hemisphere. While the surface part of the whirl is magnetically positive, the inverse is the case with the deeper part of the whirl and inversely so. The surface part shows as magnetically neutral regions. It is these deeper whirls that are suggested to be responsible for the butterfly effect. Makarov et al. (2004) found that new unipolar regions are generated during each cycle in a more or less symmetrical way near the equator (in the beginning: very irregularly). They have opposite signs in the two opposite hemispheres. The borderline between two large scale unipolar regions is an apparently neutral line (neutral in the line of sight). Zwaan (1987) was one of the first to emphasize the role of unipolar magnetic regions in the formation of coronal holes. These latter structures had already been discovered by Waldmeier (1957) in HeI observations but they only came out clearly when X ray observations became available. Their magnetic fields are unobservable or weak and in the latter case they are predominantly unipolar. The field lines stretch out away from the sun, to result in diverging open field lines (Zwaan 1987). They are present during the major part of the solar cycle but are most pronounced during sunspot minimum. As a fairly general rule they are absent around maximum sunspot number. Bohlin (1977) showed that a new coronal hole is formed after the emergence of a nearby active region. Later (quotation from Zwaan 1987) it became apparent that "a large open magnetic structure which may show a coronal hole during part of its existence is usually associated with a succession of different active regions", and that "long-lived coronal holes are caused by a specific pattern in the emergence of active regions". Transient jets occur in polar coronal holes; some of them have a motion that is transverse to the outflow direction. In the early phase of the Transition the jet-base ephemeral active regions in the coronal holes had a preferred east-west direction, that changes into the opposite in the later part, around late 2008 (Savcheva et al. 2009). This demonstrates that coronal holes find their source in the unipolar magnetic regions. The importance of coronal holes versus Coronal Mass Ejections as sources for cosmic ray modulation was studied by Cliver and Ling (2001b). The latter are more important than the former, but the latitude of the coronal hole plays a strong role in the sense that the open flux from high latitudes is strongest correlated with cosmic ray intensities.
A remarkable effect, most probably related to the present theme of the Transition was mentioned by Makarov et al. (2004). The position of the low-latitude termination of the Hale boundary came gradually closer to the equator in the course of the past 20th century Grand Maximum. During cycle 23, hence close to the time of the Transition, it was already nearing the equator (Fig. 7).
Since we now know that these regions are the areas of birth of coronal holes we conclude that during cycle 24 (and later?) coronal holes will extend down to very low latitudes, even to close to the equator. The coronal holes are the sources of (mainly weakly magnetic) solar winds and this generates the conclusion that the solar wind must be strong and mainly with a low magnetic flux in (the) coming cycle(s). This has indeed been observed (Lockwood et al. 2009, Lockwood 2013, and that finding is a fine confirmation of the Makarov-Tlatov-Callebaut observations. This phenomenon may, perhaps only partly, be related to the latitude drift of sunspots (Richardson and Schwarzschild 1953;Tuominen 1954): at latitudes below 26 • the spots move equatorward, the faster when closer to the equator, with an extreme value of 0.017 degrees per day at 2 • latitude. At higher latitudes they move poleward, but with slower speed.
During the past few Schwabe cycles the lower rest latitudes have arrived at latitudes that are fairly close to the solar equator, leaving only little space for a large unipolar region between it and the equator. The opposite situation applies in the other hemisphere. This observation too may be relevant for understanding the recent phase transition.

Decreasing Equatorial Activity Since ∼ 1970
In line with the foregoing is the general decrease of equatorial activity in the second half of the past century, as illustrated in Figs. 8 and 9.
The maximum field strength of the toroidal field as measured by R max was steadily decreasing, at least during the past two or three solar cycles. This decrease is not only seen in R max but also in those other variables of which the variability is driven by the poloidal field variations of which the strength is measured by aa min . The Total Solar Irradiance (TSI) has been nearly constant during that time and its value in 2011 was still well above the extrapolated value of the maximum of solar cycle #12, a cycle that is quite similar to #24. We note that the TSI has been measured only during the last three cycles, but it has been extrapolated downward on basis of its good correlation with the sunspot number. The number of flares (blue curve) decreased strongly after sunspot cycle 22.
The progressive decline of sunspot activity prior to the Transition is not restricted to the present transition between the 20th century Grand Maximum and the subsequent Episode. There was also "a progressive decline of solar activity before the onset of the Maunder Minimum" (Clette et al. 2014).
A year before the maximum of cycle 24 the TSI was still at a level comparable with that at solar maximum of cycle 20 and well above its value at solar maximum #12. The above differences in the variables may have relevance to the matter discussed in the preceding Section, where we described the gradual equatorward movement of the rest latitudes of the large-scale unipolar regions. Later in this paper (Sect. 11 and subsequent sections) we investigate to what degree these movements of solar field-related parameters may be interpreted as forerunners of the Transition.

Polar Variability in the Framework of the Total Solar Variation
In the later part of the polar cycle the unipolar magnetic regions gradually move poleward following and replacing the previously present unipolar large region (of opposite sign), each in its own hemisphere. Ultimately (in the subsequent cycle) the large scale unipolar regions disappear near their poles.
The poleward motion of the unipolar regions can be visualized in Fig. 9a, a diagram due to Ananthakrishnan (as early as 1954; cf. also the modern version by Hathaway in Fig. 9b). These diagrams show the distribution of prominences as a function of latitude and time. The rationale of this statement is that many prominences mark the boundaries of opposite magnetic fields (Fig. 9). This diagram also makes clear how and why the unipolar magnetic areas change into polar holes when arriving at higher latitudes. "The poleward motion of Fig. 9 (a) Isopleths of mean daily prominence area from 1905 till 1952. Remarkable are the equatorward motion of the major prominence area and the poleward motion of the high-latitude prominence zones. From Ananthakrishnan (1954). (b) A modern version of (a) that includes the recent transition; a magnetic butterfly diagram from the longitudinal averaged magnetic fields, based on observations at Kitt Peak and from SOHO. Copied from Hathaway (2010) the solar prominence areas (and the green line emission) is connected with polarity reversal and the ensuing development of the polar fields" (Cliver 2014).
Just like the equatorial region, the polar area is a seat of various manifestations of magnetic fields. These are the polar bright points and the associated ephemeral magnetic regions and polar faculae. An important feature is the coronal holes. These are non-magnetic regions or at best regions with a small magnetic field that are mainly situated in polar areas but that do sometimes extend down to fairly low latitudes as already mentioned in the preceding Section (Zwaan 1987;Makarov et al. 2004;Lockwood et al. 2009;Lockwood 2013). The behavior of the polar and equatorial activity during and prior to the phase transition is shown in Fig. 10. The polar magnetic flux DM has been measured directly at the Wilcox Solar Observatory while the equatorial flux is given by its proxy, the sunspot number R. To mark the difference between successive Schwabe cycles in the framework of their Hale cycles we have allocated alternating positive and negative signs to the data from consecutive cycles, in accordance with the polarity of the polar field. From Figs. 9 and 10 we therefore conclude that the polar field strength at solar sunspot minimum is a precursor of the subsequent sunspot maximum; cf. also Schatten et al. (1978) and Svalgaard et al. (2005).
Remarkable are the differences between the present Schwabe cycle (#24) and the two preceding ones in the period between ∼ 2000 and 2010. The maximum of sunspot cycle #23 is lower than that of the previous ones and the value of R flattens down to close to zero after about 2007. At the same time the polar cycle after ∼ 2002 is quite lower than the previous ones and seems to stay at about a constant value after ∼ 2004.
Noteworthy is also the gradual decline of the total solar magnetic field over the past few decades. Its values during cycle 23 were discussed by Hoeksema (2010); a related study of the flux emergence is by Schmieder et al. (2014). Field variation during the past three cycles is summarized by Lee et al. (2011) from whom we copied our Fig. 11. That diagram shows that the start of the gradual decline of the total photospheric magnetic flux is already observable around 1980. From Lee et al. (2011) we quote: The solar field is reduced over the entire surface of the sun, while the polar fields are about half those observed during the previous minimum period. Duhau (2003) remarked that solar activity was in a declining episode, at least as early as since 1993.
Another aspect of Fig. 10 is the following: The usually accepted scenario is that polar and equatorial maxima and minima alternate: maximum of the one coincides more or less, but very roughly, with the other's minimum and inversely so. Figure 10 shows that such is the case for cycles 21 and 22: there polar maxima and equatorial minima occur simultaneously in 1976 and 1986, but that is not the case for cycle 23. For that cycle polar maximum and equatorial minimum occur in 1994 and 1997 respectively. For cycle 24 the difference in time is still larger: they occur in 2004 and ∼ 2010 respectively. We think that this difference may be helpful for understanding the functioning of the solar dynamo during the Transition.
Summarizing this section: We illustrate the gradual decline of the total solar magnetic field over the past few decades. At the same time the large scale unipolar regions disappeared near the poles. Figure 10 shows that the time difference between polar maximum and equatorial minimum increased in the course of cycles 22 to 24. The first signs of the gradual decline of the photospheric magnetic flux are already recognizable as early as around 1980.

Summarizing Observational Data on the Transition
The various observations, summarized in the first part of this paper show the remarkable behavior of the sun during the period of the Transition. We defined it as the period from ∼ 2005 through ∼ 2010, but there are other observations indicating an earlier gradual move towards the Transition. It is even possible to state that the first indications of a precursor towards the Transition were already discernible right after the maximum of the Grand Maximum of the past century, i.e. from about 1960 onward (cf. Figs. 5,7 and 8). This shows that the Transition is not just a short-period and impulsive feature but that it is being prepared during a longer period.
These considerations force us to investigate the various periodicities in the solar dynamo, and how they manifest themselves in the photospheric and near-photospheric aspects of solar activity. Ultimately these observations ask for an examination of the stability of the tachocline and the way it directs the functioning of the solar dynamo.
We summarize the main features of the Transition: the sunspot cycle and during the Transition it took its smallest value since ∼1960. It is ascribed to a larger number of weaker spots (Fig. 4a). -Rest latitudes of unipolar magnetic fields -these being the places of birth of coronal holes -moved gradually equatorward during the past Grand Maximum and arrived close to the equator during or after cycle 24 (Fig. 7).
-The start of the gradual decline of the photospheric magnetic field was already observable around 1980 (Fig. 11). -Maximum sunspot numbers of successive cycles steadily decrease after the maximum (1960) of the 20th century Grand Maximum (Fig. 8). -As a general rule the high-latitude prominence areas move poleward after the equatorward motion of the lower latitude prominences (Fig. 9). -The polar magnetic maximum occurred during the Transition some 6 years after equatorial magnetic minimum (Fig. 10).   Other information on the starting time of the Transition can be derived from observations of the solar microwave radiation, of which the 10.7 cm observations form a continuous sequence since around 1955. The variations in this flux are due to those in the equatorial facular fields, just as the Total Solar Irradiance (TSI) is. Hence, the 10.7 cm radiation is fully comparable to the TSI (cf. also our Fig. 3), but the amplitude of the variations in 10.7 cm is much more pronounced. The variations in this radiation match excellently with the frequency shift of solar internal oscillations derived from helioseismic observations (Broomhall et al. 2009), but Basu et al. (2012;cf. also Balogh et al. 2014, p. 11 and Fig. 5) found deviations from that match that were most significant for the observations at lowest frequencies (1860-2400 μHz). At these frequencies deviations are first discernible in 1995 and only later at higher frequencies.

Starting Time of the Transition
From these data we conclude that (a) first signs of a forthcoming change of the kind of Episode were visible in the eighties of the last century, that (b) very first precursors to the Transition were already visible around 1995 and that (c) the full start came some ten years later. The implications for the driving mechanism will be discussed in the following Sections.

The Phase Diagram and the Predictability of Long-Term Solar Variability
Phase transitions occur between Grand Episodes of solar variability. A very well-known Grand Episode is the Maunder Minimum that lasted from about 1610 till about 1740. Another episode, one of the highest of the past ten thousand years (Usoskin and Mursula 2003;Solanki et al. 2004) is the recent 20th Century (Grand) Maximum, this being the period of large solar activity during the main part of the 20th century. It was characterized by strong Schwabe cycles with their associated phenomena such as an abundant number of solar flares and also by related features such as the increased interplanetary magnetic flux (Solanki et al. 2000). The polar flux too reached at that time higher values than ever measured before. The question is relevant if there is anything that is systematic or predictable in the phase transitions between successive Episodes. To that end Duhau (2003) introduced a phase diagram, consisting of a plot of the proxy of the equatorial magnetic field R max against that of the polar field aa min (for the definition of the latter cf. Fig. 1 and the corresponding text). The resulting phase diagram is shown in the two components of Fig. 12. Here, some comments are in order.
Empirically, Duhau found that there exist a unique point in that diagram, with coordinates R max = 93.4 (defined in the earlier system of sunspot counting) and aa min = 10.3, that she baptized the Transition Point.
Actually, the two coordinates of this point were initially found as the average values of R max resp. aa min over the past few centuries. After having derived these two values (Duhau 2003), it was a fascinating addition to dynamo theory to find that in the past four centuries new solar activity episodes tend to start when the (R max , aa min ) curve in the phase diagram crossed this Transition Point or passed near to it . This empirical finding of the special character of the Transition Point may have implications for our understanding of some of the mechanisms active in the solar dynamo.
Similar phase transitions took place around 1628 (onset of the Maunder Minimum), in 1740 (end of Maunder Minimum, and transition to a period of Regular Oscillations; our Table 1 gives 1719, on basis of the phase diagram of Fig. 12) in 1923/24 (onset of 20th century Grand Maximum), and the most recent one started around 2008. Other possible transitions may have occurred in 1230, 1349 and 1405 (Duhau and de Jager 2008), but these are unverifiable, being based on extrapolated, hence uncertain values of the proxy data of the magnetic fluxes: First, the aa min data are only known after 1844 and earlier aa min data have been derived by Nagovitsyn (2005) by extrapolation. But this extrapolation rests for an important part on the high degree of correlation between R max and aa min , while the essential aspects of the phase diagram are based on the subtle differences between both. Therefore only the late (1800-present) diagram is reliable. We publish the earlier one too, for illustrative purposes.
Another comment is the following: the R max data in the phase diagrams that have been published so far were always derived by using the 'classical' Zürich international sunspot numbers (ISN). We have decided to derive the present phase diagrams on basis of the new Group Numbers (GSN; Svalgaard and Schatten 2016; cf. end paragraph of Sect. 1 for the definitions). In order to eliminate relatively shortlived variations such as those due to volcanoes, el Nino, the Schwabe cycle and the like, the Svalgaard-Schatten GSN data have been smoothed with a Lowess smoothing technique, by using an 18 years smoothing interval, which is the optimal value (de Jager and Usoskin 2006). By this smoothing these fluctuations are not eliminated but they are no longer dominating.
Thus, the data figuring in the phase diagrams of Fig. 12 were derived. The Modified International Sunspot Numbers have values of the order of many tens to above one hundred while the classical International Sunspot Numbers and the associated Group Sunspot Number values are about ten times smaller. Therefore, the R max value of the Transition Point, from here onward applied to the GSN data, is in the new diagrams roughly some ten times smaller than the value that was used previously when it was based on the ISN data.
The diagrams shows close passages to the transition point to occur in 1628 (onset of the Maunder Minimum) and again in 1719 (onset of the period of Regular Oscillations). Hence, according to this diagram, the Maunder Minimum terminated around 1719, although direct observations rather favor the time around 1740. And indeed, the next following nearest year  de Jager and Nieuwenhuijzen (2013). In the present copy use is made of Lockwood's homogenized aa min data (Lockwood, priv. com.), smoothed by an 18 years Lowess technique and the newly introduced Sunspot Group Numbers (Svalgaard and Schatten 2016), smoothed likewise is 1746, but its distance to the Transition Point has then grown to about 1.5 times that of 1719. This deviating result may be due to the fact that the aa min data from before 1844 are not based on direct measurements but on extrapolations (as stated above) and may, hence, not be fully reliable.
A later close passage occurred in 1924. That year does not deviate much from the one given by Lockwood et al. (2009)  We next examine if there is some system in the 'close passages' to the Transition Point. To that end we define the 'distance of closest passage' (de Jager and Nieuwenhuijzen 2013) as the root of the quadratic sum of the differential aa min and R max values, both with respect to the coordinates of the Transition Point and both normalized by dividing them by the values of these coordinates. These new coordinates are based on the Group Numbers; they are 7.51 for R max while we confirm the value 10.32 for aa min , which is equal to the values used in earlier studies. We calculate the distance to the Transition Point as distance = ( R max /7.51) 2 + ( aa min /10.32) 2 1/2 . (1) The resulting data are specified in Table 1. Table 1 gives the year 1924 as the starting date for the 20th century Grand Maximum. That year is derived from the phase diagram (Fig. 12) in the way described above. We mentioned already that Lockwood et al. (2009)  The present authors were hoping that the value of the 'distance' might, in some way be related to the character of the next grand episode but there are only weak such indications: the start of the Maunder Minimum, being the most extreme of the recent Grand Episodes, is related to the shortest distance. Clearly, a few more centuries have to pass before reliable conclusions can be taken. An essential result is that the extended transition phase, that lasted from ∼ 2005 till ∼ 2010, is related to the start of another Grand Episode. Hence, it marks the end of the 20th century Grand Maximum (Duhau andde Jager 2008, 2010; cf. also the review de Jager and Duhau 2012).
The next question is that of the character of the next Grand Episode. A Grand Minimum does not seem probable. Steinhilber et al. (2009Steinhilber et al. ( , 2010 have shown that such type minima occur only during negative phases of the Hallstatt period of solar activity and that statement was confirmed by de Jager and Duhau (2012). The Hallstatt cycle (Clilverd et al. 2003(Clilverd et al. , 2006Muscheler et al. 2004) has a length of 2300 years, and it is, as observations show (cf. Fig. 10 of Versteegh 2005), based on the regular accumulation of Grand Minima. These accumulations repeat themselves with a periodicity of ∼ 2300 years. Since the Hallstatt Fig. 13 The recent Hallstatt oscillation (green; right ordinate), and the Gleissberg cycle (blue; left ordinate) and their respective predictions till the year 2200 (black lines). The zero line was crossed at the beginning of the of the 21th century. These results are based on the Modified International Sunspot Numbers (MISN; Clette et al. 2014Clette et al. , 2015, extended backward from 1705 till 290 by the Schove (1955) time series cycle passed through its zero line early in the 21th century (our Fig. 13) one may not expect a Grand Minimum to come. A Grand Maximum is not very probable either, on statistical grounds, for the simple reason that a succession of two such maxima has not been observed before. Therefore, another period of Regular Oscillations, such as occurred between 1740 and 1923 seems the most probable one (cf. also de Jager and Duhau 2012). This conclusion is confirmed in a study presently underway. We interpret the recent phase transition as the crucial step towards that Episode.
Sunspot cycle 23 occurred close to the beginning of the new Episode. The next following cycle, 24, with the observed smoothed maximum value of 68 sunspot numbers -the smallest of the past century -was reached after a small but significant Gnevyshev gap (for its definition cf. Basilevskaya et al. 2006). From the more successful early forecasts of that maximum we mention e.g. Duhau (2003) who came already at an early time to an expected sunspot maximum of 87.5 ± 23.5 Another forecast based on similar points of departure is that of  and that of  who expected for the maximum sunspot number of cycle 24 values that were 30 % and 35 % respectively smaller than that of cycle 23, of which the maximum sunspot number was 121. These forecasts too were only slightly higher that the observed value. Later forecasts too came closer to the observed values, e.g. Solanki and Krivova (2011) foresaw "between 60 and 100".

The Tachocline; The Origin of Sunspots
A dynamo is a physical system that transforms kinetic into magnetic energy. In this Section we restrict our discussion of the solar dynamo to the aspect that is important for the origin of sunspots. Essentially, we refer to a few properties of the Babcock-Leighton phenomenological model. In choosing that model we base ourselves also on Muñoz-Jaramillo et al. (2013a, 2013b, Cameron and Schüssler (2015) and Petrie (2015) who find that this model well covers the observational data.
The solar dynamo is situated in the outer convection zone of the sun. That zone has a depth of ∼ 2 × 10 5 km (Christensen-Dalsgaard 1991;Basu 1997;Brun and Toomre 2002;Cally 2003). In order to allow for the observed time scale of solar activity, which is of the order of 10 years, a skin depth of the relevant oscillator is needed of the order of 10 5 km. This is one of the reasons for identifying the dynamo with parts of the outer convection zone (review by Rüdiger 2008). At the bottom of that zone the correlation time of the turbulent cells is about equal to the value of the solar rotation period, and that is one of the reasons for situating the essential part of the dynamo in that region. Part of the dynamo is the tachocline which consists of the shearing layers at and below the base of the outer convection zone. The tachocline is a relatively thin, but as many observations show, not perfectly spherical layer (Elliot and Gough 1999;Antia and Basu 2011). It is 27,000 km thick at the equator and 30,000 km at latitude 60 • (Charbonneau et al. 1999). The outer boundary is prolate while the inner boundary is nearly spherical (Antia and Basu 2011). It is situated at an average depth of 215,000 km below the solar surface, just below the base of the outer solar convection zone, where the combination of convective motions and differential rotation powers the magnetic fields. Detached elements of the internal fields originate by kink-instability and may thereupon ascend to appear at the surface as the various elements of solar variability. These are, among other phenomena, the sunspot belt, the polar magnetic areas and the large unipolar areas.
At the tachocline level the interplay of convective motions and differential solar rotation leads to the origin of a complicated nonlinear system of magnetic fields, based on associated electric currents (Charbonneau et al. 1999;Basu and Schou 2000;Käpylä et al. 2006). The tachocline, hence, is the result of an interplay between the motions generated by the large-scale differential solar rotation and the small-scale stochastic motions due to the anisotropic turbulence generation by the convective motions in the outer solar convection zone. This interplay between motions at very different scales greatly complicates the theory of the dynamo (Brandenburg et al. 2012), which is (Weiss et al. 1984;Weiss 1988) an example of deterministic chaos. In the traditional theoretical approach the high conductivity hampers the development of large-scale flow and associated large-scale fields, but in a different approach by Tobias and Cattaneo (2013) high-resolution dynamo simulations succeed in producing large-scale structures and flows at high magnetic Reynolds numbers. In their solution propagating dynamo waves are a dominant feature of the solar cycle.
Due to the internal solar shearing motions the magnetic fields are stretched into a direction parallel to the solar equator. Thus, ultimately the system of magnetic fields becomes mainly toroidal. By continued stretching the toroidal field strength increases in the course of subsequent solar rotations. When these field strengths reach a critical value, the resulting toroidal field becomes unstable. Smaller closed loops detach from it, due to kink instability. They have smaller specific weights and hence stronger buoyancy than the surrounding solar gas, because of their high magnetic field strengths. Therefore these loops rise, to reach the solar photosphere after one or a few months. There they appear as a pair of sunspots (Solanki et al. 2004). This hypothetical scenario for the origin of sunspots is supported by helioseismic observation of sunspots on their way upward; at the time of the observations they were still as deep as 65,000 km below the solar surface (Ilonides et al. 2011).
For a quantitative analysis of the evolution of the solar magnetic field it is important to know the above mentioned critical value for which sunspot forming loops originate by kink instability of the tachocline. Early determinations are 100 kilogauss (D' Silva and Choudhuri 1993;Caligari et al. 1995 andRuzmaikin 1998) and 45 kilogauss (Dikpati 2005). On basis of these data de Jager (2005) suggested a critical value of 80 kG. More recent work show that a minimum field > 40 kG is needed to form solar active regions (Weber et al. 2013). Quoting from Weber: At 100 kG, the development of Omega-shaped rising loops is mainly controlled by the growth of the magnetic buoyancy instability. However, at low field strengths of 15 kG, the development of rising Omega-shaped loops is largely controlled by convective flows, and properties of the emerging loops are significantly changed compared to previous results in the absence of convection. With convection, rise times are drastically reduced (from years to a few months), loops are able to emerge at low latitudes, and tilt angles of emerging loops are consistent with Joy's Law for initial field strengths of greater than or equal to 40 kG. In a more recent study in which the effects of radiative heating in turbulent solar-like convection are included, Weber and Fan (2015), found that a minimum value is needed that should in any case be stronger than 60 kG. On basis of the above data we suggest to use in our further calculations a critical value of 60 kG. The uncertainty of this value is obvious; let us assume a value of ±20 kG.
Extended descriptions and analyses of the dynamo are found in many review papers. We refer to the papers by Küker et al. (2001), Käpylä et al. (2006) and the reviews by Fisher et al. (2000); Tobias (2002);Ossendrijver (2003); Charbonneau (2010); de Jager and Duhau (2011);Schmieder et al. (2014);Fletcher et al. (2015). Many more papers could be mentioned but we do not consider that necessary because a treatment of dynamo theories is outside the scope and purposes of the present paper. Only one essential aspect of the dynamo theories is considered in this paper. It is the matter of rising loops that appear as sunspots and activity centers after having detached from the tachocline by kink-instability, viz. the abovedescribed mechanism that is, along with other authors, dealt with by D' Silva and Choudhuri (1993), Caligari et al. (1995), Ruzmaikin (1998), Dikpati (2005, Weber et al. (2013) and Weber and Fan (2015).
A related additional aspect is that the dynamo is a nonlinear system governed by deterministic chaos (Weiss 1988(Weiss , 2010Mininni and Gomez 2002;Charbonneau et al. 2004;Weiss and Thompson 2009). But this should be specified: Quoting Weiss (2010): While modulation leading to grand maxima and grand minima does most probably correspond to the non-linear pattern revealed in simplistic dynamo models, the small-scale turbulent motion cannot be modeled in detail. Based on an earlier study by ,  could show that the occurrence of Grand Minima could be imitated in numerical experiments on the basis of a simple dynamo model with random fluctuations of the dynamo governing parameters. They also found that during the Grand Minimum (reference: the Maunder Minimum) the Schwabe cycle continues to exist, though at a very reduced scale. These studies lead us to conclude that the solar dynamo is a non-linear rather than a chaotic system.
The essential questions are here: what values of the dynamo parameters are critical in driving the grand minima?  found that grand minima occur for specific values of the parameters in the dynamo equations leading to random fluctuations in the dynamo governing equations. And why does another Grand Episode start when the proxies for the polar and the equatorial magnetic fluxes arrive simultaneously at their average values (= the Transition Point)? This mere fact forces us to consider the fundamental question: which dynamo parameters were critical in driving the solar dynamo during or around the phase Transition?

The Tachocline and the Transition; An Attempt to Synthesize
We want to use the observations described earlier in this paper (cf. Sects 9 and 10), together with the considerations described in Sect. 12, to improve our understanding of the structure and functioning of the solar equatorial activity during and preceding the Transition.

Decrease of Umbral Magnetic Field Strengths
Sunspots originate in kink-unstable loops that detach from the toroidal field in the tachocline. These loops thereupon rise to appear at the solar surface as sunspots after a couple of months. They can only detach from the toroidal field in the tachocline when the toroidal field in the tachocline has reached a minimum magnetic field strength. An average value of that value, derived from data given by various authors (summarized in the preceding section) is 60 kilogauss. That value does not depend markedly on the position of the tachocline in the sun. The umbral field strength at the time they arrive at the surface depends on the difference in gas density between the tachocline and the solar surface regions. The components of Fig. 4 shows that the average umbral field strength decreased during the 10 years of the Transition from ∼ 2700 to ∼ 2000 Gauss. Since the surface magnetic field in the sunspots depends on the distance traversed from the tachocline to the surface, this must imply that this distance has increased during the years of the Transition. In other words: the equatorial diameter of the tachocline has gradually decreased between the years 2000 and 2010. There are as yet (2016) no indications that it started to increase, as happened before during earlier Schwabe cycles. The emerging picture is that a least during the past Grand Maximum the average umbral field strength was oscillating, and keeping pace with the Schwabe cycles. Figure 7 shows that the so-called rest latitudes of the neutral lines between the large scale unipolar regions decreased gradually over the past 11 Schwabe cycles. The lower one of the rest latitudes gradually approached the solar equator during the 20th century Grand Maximum, leaving a small band of unipolar region growing from around the equator to the (lower) rest latitude. This is another signal of the Transition. This may be an indication of a gradual evolution of the meridional circulation, which may, in turn, be related to some kind of relatively longperiod oscillation of the tachocline.

A Pulsating Tachocline?
In this final Section we introduce the hypothesis that the tachocline undergoes systematic motions with relatively small amplitudes, such that its equatorial parts show a component that moves periodically up-and downward. These motions could be called pulsations but they might, upon further research, appear to have another character.
The magnetic loops that will appear as sunspots at the photospheric level, have detached themselves from the tachocline when the latters field strength had grown to 6 × 10 4 Gauss, a value that is the average of the various determinations that have so far been published, as described earlier in this paper. We add at this point that this value of 6 × 10 4 Gauss is not critical for the discussion below, as will appear also from the following: While rising, the loops expand and their field strengths decrease, in relation to the (decreasing) surrounding gas density; this causes the loops to expand. The deeper the tachocline source is situated, the smaller is therefore the resulting sunspot's photospheric magnetic field. Hence, the actual problem is that of the relation between the depth variation and related amplitude of the pulsating tachocline and the observed photospheric sunspot field strengths.
Assume an average ('standard') tachocline depth z(0), where z is the vertical ordinate. We place z = 0 at this assumed 'standard' tachocline depth. The consequent umbral magnetic field corresponding with the standard tachocline depth will be called B 0 . We now assume that the tachocline moves up or down, above or below the standard level. When the tachocline is situated at another depth than the standard value the consequent umbral magnetic field also assumes another value, that we will call B phot . Hence the ratio between the observed photospheric field strength B phot and the standard strength B o is where P is the gas pressure. The next problem is that of the relation between the depth of the tachocline and the value of B phot , in short that of the distance between the standard tachocline depth and the actual depth, in relation to the photospheric values. In still other words: we want to find the semi-amplitude of the assumed pulsation.
Combining Eqs. (2) and (3) yields The various components of Fig. 4 show that in the period between 1960 and 2000 the average umbral magnetic field had an average strength of ∼ 2450 G (cf. Sect. 4) with a weakly defined semi-amplitude of some 300 G. According to Eq. (4) this implies a tachocline pulsation with a semi-amplitude of ∼ 12,000 km. The full pulsation, twice this value, may be compared with the tachocline thickness, about 25,000 to 30,000 km. Hence, we are dealing with a fairly small full amplitude, still inside the average ('standard') thickness of the tachocline.
With regard to the Transition we find that a decrease of the sunspot field from 2500 to 2000 Gauss, as happened between the years 2000 and 2010 (cf. Figs. 4a and 4b), corresponds, according to Eq. (4), with a semi-amplitude of 21,000 km. Hence, the recent Transition can be explained by a downward pulsation of the tachocline over that distance. Again: a distance smaller than the thickness of the tachocline. The small amplitude of the pulsations, where the tachocline remains in broad outline within its boundaries may explain why this hypothetical pulsation has not been detected by helioseismic techniques. Yet, Antia and Basu (2011) find from their helioseismic analysis that some parameters do show changes with time, while Basu et al. (2012) mention the thinning of the relevant layers during solar cycle 23.

Pulsations?
The question is relevant whether the tachocline is really 'pulsating'. Another viewpoint may be that the tachocline, as observed by helioseismology, is undergoing a differential rotation with respect to the solar core that leads to torsional oscillation of the whole convective layer (cf. Charbonneau 2010). In that line,  wrote ". . . we conclude that torsional oscillations stabilize the tachocline-convective layer system motions, leading to transfer of angular momentum in the latitudinal direction from the core spin to the tachocline. . . " Flux transport dynamo models (for a review see Dikpati and Gilman 2009) seem to contain all the necessary ingredients for realistic simulations of solar dynamo action. However, they predict a strong sunspot cycle 24 (Dikpati et al. 2006), comparable with those that occurred in the recent Grand Maximum. Observations instead indicate that this episode has ended (see . While these models assume a constant poloidal source at the tachocline-convective layer boundary these recent results indicate that the failure of that model is due to the fact that the conditions at the tachoclineconvective layer boundary suddenly changed during the transition(s), a possibility that is not incorporated in the flux transport models. These changing boundary conditions during the chaotic transitions are likely due to a sudden change of orientation of the tachoclineconvective motions with respect to the relic field seated at the core. In fact, for starting the dynamo loop, flux transport dynamo models need a seed field at the tachocline-convective layer boundary. Also, a relic field appears to be essentially necessary for having the radiative interior rotating as a rigid body and evidences indicating that his field is pointing southward have been provided since Cowling (1945) suggested its existence (reviews in Mursula et al. 2001, Duhau and. The persistence of a negative polar field component at the transition, as found here, gives further support to these findings.

Non-linear Pulsations?
These foregoing considerations seem to support the hypothesis that the tachocline may present motions of different kind that could be interpreted as pulsations with a small, hardly detectable amplitude. The most striking component is the pulsation of the equatorial diameter with an eleven years period and with a semi-amplitude of 12,000 km (0.017R 0 ). Another conclusion is that the Transition corresponds to a downward displacement of the tachocline over 21,000 km.
There is also evidence that there must be more periodicities acting than this sole one. This evidence is given by the four additional major periodicities: the two Gleissberg periods (the short and the long one), the Suess and the Hallstatt periodicities. Their analysis can be used to yield information on the variation of the dynamo (cf. also Fig. 13): During the negative phases of the Hallstatt periodicity of some 2000 years, Grand Minima cluster together (Steinhilber et al. 2010). They do not occur during the positive phase. It is known that the minimum umbral field strength that is needed for a spot to be visible as such, is ∼ 1500 G. In the very extreme case of a hypothetical Grand Minimum with zero sunspots, the umbral field strength must, hence, be smaller than or at most equal to 1500 G. That would be the case for a downward movement of the tachocline over 45,000 km (0.06R 0 ). It seems probable that this fairly large shift of the tachocline depth is the reason why such extreme Maunder-type minima have not yet been observed.

Origin of Pulsations
The origin of the tachocline 'pulsations' could most probably be found in the frame of the meridional circulation because the physical conditions in the tachocline area are at least partly related to the meridional circulation. By this circulation, magnetic flux and angular momentum are transported. The circulation consists of large cells that show up near the surface as an equatorward flow of ∼ 20 m/s. The return flow, deeper inside, is some ten times slower. The total circulation time is ∼ 10 years. The situation that Grand Minima occur solely during negative phases of the Hallstatt oscillation may be explained by the assumption that the Hallstatt oscillation is due to a change of symmetry of the tachocline-convective motions with respect to the relic field seated at the core.
But the situation may be more complicated. Noteworthy is the finding of a double (or even triple) cell structure (Haber et al. 2002;Zhao et al. 2013). Haber et al. found that the equatorward streaming covers depths between 0.83 and 0.92R sun and the poleward one between 0.75-0.83R sun and 0.92-1.0R sun . These boundaries are no fixed quantities, they change pseudo-randomly in the course of the solar cycle. Variation of the flow during the Transition are noteworthy: during solar cycle 23 the meridional flow extended up to the poles while in cycle 22 the circulation had switched equatorward near latitude 60 • (Dikpati et al. 2010).The authors suggested that this may have been at the origin of the longer duration of cycle 23. Nandy et al. (2011) introduced a fast meridional flow in the first half of the solar polar cycle followed by a slower flow in the second half. It appears that this introduction is able to reproduce the main characteristics of the sunspot cycle 23. Zhao et al. (2013) even find indications that the meridional flow does not extend as deep as to the bottom of the convection zone. Various computer simulations by other authors arrive at similar results: the velocity of the meridional circulation and its changes seem an essential aspect for explaining the Transition. More stochastic phenomena may occur too: a case is described in which a submerged cell, originating in the northern hemisphere disrupts the poleward flow and the symmetry in the streaming (Haber et al. 2002).
For understanding the assumed tachocline 'pulsations', we note that the radial distance of the tachocline to the solar center depends, among other things, on the amount of downward overshoot from the convective region. A sheer layer of radial differential rotation defines the conditions at the bottom of the convection zone (Choudhuri et al. 1995). This limits the meridional circulation to a lower bound of ∼ 0.7R sun (Gilman and Miesch 2004) and we suggest that this, in turn, influences or even determines the lowest position that the tachocline can take and the consequent umbral field strength. The physical origin of the tachocline pulsations may be sought in this direction. Numerical model calculations (Sule et al. 2005) indeed show how the meridional circulation changes the shape and position of the tachocline. Near the pole the tachocline is thinnest. This may find its origin in the counter-rotating meridional circulation that appears to develop in the polar regions MacGregor 2010, 2011).
It is obvious that the above considerations must be worked out and developed on base of an MHD analysis. A promising development that runs more or less parallel to this line of thought is found in the work of Mininni and Gomez (2002) and in that of Hiremath (2010). In the first quoted paper the MHD equations are developed and described. In the second paper, Hiremath (2010), mentions two possible ways of approach in understanding the solar dynamo: the turbulent dynamo and the MHD oscillatory mechanisms, and discusses the latter on basis of the two equations for an axisymmetric magnetic field structure in cylindrical geometry, as forwarded by Parker (1955). This theme has more recently been studied by  and .These equations couple the differential rotation to the poloidal and toroidal fields. Their application forecasts Alfvén waves that originate in the convective layer. They do explain the 22-years periodicity and it is obvious from the used equations, that there must be more such quasi-periodicities. It is not improbable that these are at the basis of the 'pulsations' of the tachocline.
The crucial question, yet to be studied, is if this MHD approach makes it possible to obtain quantitative agreement with the postulated up-and downward motions of the tachocline. Answering this question goes beyond the present review but is an obvious follow-up.