A Proposed Paradigm for Solar Cycle Dynamics Mediated via Turbulent Pumping of Magnetic Flux in Babcock-Leighton type Solar Dynamos
At present, Babcock-Leighton flux transport solar dynamo models appear as the most promising model for explaining diverse observational aspects of the sunspot cycle. The success of these flux transport dynamo models is largely dependent upon a single-cell meridional circulation with a deep equatorward component at the base of the Sun’s convection zone. However, recent observations suggest that the meridional flow may in fact be very shallow (confined to the top 10% of the Sun) and more complex than previously thought. Taken together these observations raise serious concerns on the validity of the flux transport paradigm. By accounting for the turbulent pumping of magnetic flux as evidenced in magnetohydrodynamic simulations of solar convection, we demonstrate that flux transport dynamo models can generate solar-like magnetic cycles even if the meridional flow is shallow. Solar-like periodic reversals is recovered even when meridional circulation is altogether absent, however, in this case the solar surface magnetic field dynamics does not extend all the way to the polar regions. Very importantly, our results demonstrate that the Parker-Yoshimura sign rule for dynamo wave propagation can be circumvented in Babcock-Leighton dynamo models by the latitudinal component of turbulent pumping – which can generate equatorward propagating sunspot belts in the absence of a deep, equatorward meridional flow. We also show that variations in turbulent pumping coefficients can modulate the solar cycle amplitude and periodicity. Our results suggest the viability of an alternate magnetic flux transport paradigm – mediated via turbulent pumping – for sustaining solar-stellar dynamo action.
Department of Physical Sciences, Indian Institute of Science Education and Research, Kolkata \move@AU\move@AF\@affiliation Center of Excellence in Space Sciences India, IISER Kolkata, Mohanpur 741246, West Bengal, India
The cycle of sunspots involves the generation and recycling of the Sun’s toroidal and poloidal magnetic field components. The magnetohydrodynamic (MHD) dynamo mechanism that achieves this is sustained by the energy of solar internal plasma motions such as differential rotation, turbulent convection and meridional circulation. The toroidal field is generated through stretching of the poloidal component by differential rotation [Parker 1955] and is believed to be stored and amplified at the overshoot layer [Moreno-Insertis et al. 1992] beneath the base of the solar convection zone (SCZ). Strong toroidal flux tubes are unstable to magnetic buoyancy and erupt through the surface producing sunspots, which are strongly magnetized and have a systematic tilt [Hale 1908, Hale et al. 1919]. The poloidal field is believed to be regenerated through a combination of helical turbulent convection (traditionally known as the mean-field -effect; [Parker 1955]) in the main body of the SCZ and the redistribution of the magnetic flux of tilted bipolar sunspot pairs (the Babcock-Leighton process; [Babcock 1961, Leighton 1969]).
Despite early, pioneering attempts to self-consistently model the interactions of turbulent plasma flows and magnetic fields in the context of the solar cycle [Gilman 1983, Glatzmaier 1985] such full MHD simulations are still not successful in yielding solutions that can match solar cycle observations. This task is indeed difficult, for the range of density and pressure scale heights, scale of turbulence and high Reynolds number that characterize the SCZ is difficult to capture even in the most powerful supercomputers. An alternative approach to modelling the solar cycle is based on solving the magnetic induction equation in the SCZ with observed plasma flows as inputs and with additional physics gleaned from simulations of convection and flux tube dynamics. These so called flux transport dynamo models have shown great promise in recent years in addressing a wide variety of solar cycle problems [Charbonneau 2010, Ossendrijver 2003].
In particular, solar dynamo models based on the Babcock-Leighton mechanism for poloidal field generation have been more successful in explaining diverse observational features of the solar cycle [Dikpati & Charbonneau 1999, Nandy & Choudhuri 2002, Chatterjee et al. 2004, Choudhuri et al. 2004, Guerrero & de Gouveia Dal Pino 2007, Nandy et al. 2011, Choudhuri & Karak 2012, Hazra et al. 2014, Passos et al. 2014]. Recent observations also strongly favor the Babcock-Leighton mechanism as a major source for poloidal field generation [Dasi-Espuig et al. 2010, Muñoz- Jaramillo et al. 2013]. In this scenario, the poloidal field generation is essentially predominantly confined to near-surface layers. For the dynamo to function efficiently, the toroidal field that presumably resides deep in the interior has to reach the near-surface layers for the Babcock-Leighton poloidal source to be effective. This is achieved by the buoyant transport of magnetic flux from the Sun’s interior to its surface (through sunspot eruptions). Subsequent to this the poloidal field so generated at near-surface layers must be transported back to the solar interior, where differential rotation can generate the toroidal field. The deep meridional flow assumed in such models (See Fig. 1, left-hemisphere) plays a significant role in this flux transport process and is thought to govern the period of the sunspot cycle [Charbonneau & Dikpati 2000, Hathaway et al. 2003, Yeates et al. 2008, Hazra et al. 2014]. Moreover, a fundamentally crucial role attributed to the deep equatorward meridional flow is that it allows the Parker-Yoshimura sign rule [Parker 1955, Yoshimura 1975] to be overcome, which would otherwise result in poleward propagating dynamo waves in contradiction to observations that the sunspot belt migrates equatorwards with the progress of the cycle [Choudhuri et al. 1995, Hazra et al. 2014, Passos et al. 2015, Belusz et al. 2015].
While the poleward meridional flow at the solar surface is well observed (Hathaway & Rightmire 2010; 2011) the internal meridional flow profile has remained largely unconstrained. A recent study utilizing solar supergranules [Hathaway 2012] suggests that the meridional flow is confined to within the top 10% of the Sun (Fig. 1, right-hemisphere) – much shallower than previously thought. Independent studies utilizing helioseismic inversions are also indicative that the equatorward meridional counterflow may be located at shallow depths [Mitra-Kraev & Thompson 2007, Zhao et al. 2013]. The latter also infer the flow to be multi-cellular and more complex. These studies motivate exploring alternative paradigms for flux transport dynamics in Babcock-Leighton type models of the solar cycle which are crucially dependent on meridional circulation linking the two segregated dynamo source regions in the SCZ. This leads us to consider the role of turbulent pumping.
Magnetoconvection simulations supported by theoretical considerations have established that turbulent pumping preferentially transports magnetic fields vertically downwards [Brandenburg et al. 1996, Tobias et al. 2001, Ossendrijver et al. 2002, Dorch & Nordlund 2001, Käpylä et al. 2006, Pipin & Seehafer 2009, Racine et al. 2011, Rogachevskii et al. 2011, Augustson et al. 2015, Warnecke et al. 2016, Simard et al. 2016] – likely mediated via strong downward convective plumes which are particularly effective on weak magnetic fields (such as the poloidal component). In strong rotation regimes, there is also a significant latitudinal component of turbulent pumping. In particular, two studies, one utilizing mean-field dynamo simulations [Brandenburg et al. 1992] and the other utilizing turbulent three dimensional magnetoconvection simulations [Ossendrijver et al. 2002] recognized the possibility that turbulent pumping may contribute to the equatorward propagation of the toroidal field belt. We note that most Babcock-Leighton kinematic flux transport solar dynamo models do not include the process of turbulent pumping of magnetic flux. The few studies that exist on the impact of turbulent pumping in the context of flux transport dynamo models show it to be dynamically important in flux transport dynamics, the maintenance of solar-like parity and solar-cycle memory [Guerrero & de Gouveia Dal Pino 2008, Karak & Nandy 2012, Jiang et al. 2013]. In their model with turbulent pumping, Guerrero & de Gouveia Dal Pino (2008) used a spatially distributed -coefficient in the near-surface layers to model the Babcock-Leighton poloidal source and a meridional circulation whose equatorward component penetrated up to , i.e., more than half the depth of the SCZ; therefore, from this modelling it is not possible to segregate the contributions of turbulent pumping and meridional flow (the peak latitudinal component of the former coincides with the equatorward component of the latter) to the toroidal field migration.
Here, utilizing a newly developed state-of-the-art flux transport dynamo model where a double-ring algorithm is utilized to model the Babcock-Leighton process, we explore the impact of turbulent pumping in flux transport dynamo models with nonexistent, or shallow meridional circulation. Our results indicate the possibility of an alternative flux transport paradigm for the solar cycle in which turbulent pumping of magnetic flux resolves the problems posed by a shallow (or inconsequential) meridional flow.
Our flux transport solar dynamo model solves for the coupled, evolution equation for the axisymmetric toroidal and poloidal components of the solar magnetic fields:
where, is the toroidal component of magnetic field and is the vector potential for the poloidal component of magnetic field. is the meridional flow, is the differential rotation, is the turbulent magnetic diffusivity and . For the differential rotation and diffusivity profile, we use an analytic fit to the observed solar differential rotation (the near-surface shear layer is not included) and a two-step turbulent diffusivity profile (which ensures a smooth transition to low levels of diffusivity beneath the base of the convection zone) (For detailed profile, see Hazra & Nandy 2013). We use the same meridional flow profile as defined in Hazra & Nandy (2013). Our flow profile has penetration depth of to represent deep meridional flow situation, and to represent shallow meridional flow situation. We set the peak speed of the meridional flow to be 15 ms (near mid-latitudes). The second term on the RHS of the toroidal field evolution equation acts as the source term for the toroidal field (rotational shear), while in the poloidal field evolution equation, the source term, , is due to the Babcock-Leighton mechanism. Here we use a double-ring algorithm for buoyant sunspot eruptions that best captures the Babcock-Leighton mechanism for poloidal field generation [Durney 1997, Nandy & Choudhuri 2001, Muñoz-Jaramillo et al. 2010, Hazra & Nandy 2013] and which has been tested thoroughly in other contexts. Specifics about our double ring algorithm can be found in Hazra & Nandy (2013) and Hazra (2016; PhD Thesis).
To bring out the significance of the recent observations, we first consider a single cell, shallow meridional flow, confined only to the top 10% of the convection zone (Fig. 1, right-hemisphere). In the first scenario we seek to answer the following question: Can solar-like cycles be sustained through magnetic field dynamics completely confined to the top 10% of the Sun?
In these simulations initialized with antisymmetric toroidal field condition (with initial B 100 kG), we first allow magnetic flux tubes to buoyantly erupt from 0.90 (i.e., the depth to which the shallow flow is confined) when they exceed a buoyancy threshold of Gauss (G). In this case we find that the simulated fields fall and remain below this threshold (at all latitudes at 0.90 ) with no buoyant eruptions, implying that a Babcock-Leighton type solar dynamo cannot operate in this case. Dikpati et al. (2002) considered the contribution of the near-surface shear layer in their simulations (which we have not) and concluded that this near-surface layer contributes only about 1 kG to the total toroidal field production and hence insufficient to drive a large-scale dynamo. Guerrero & de Gouveia Dal Pino (2008) also utilized a near-surface shear layer with radial pumping and found solar-like solutions only under special circumstances; however, given that for this particular case they utilized a local -effect for the latter simulations (with a spatially distributed -effect in the near-surface layer) it is not evident that these simulations are relatable to the Babcock-Leighton solar dynamo concept. The upper layers of the SCZ is highly turbulent and storage and amplification of strong magnetic flux tubes may not be possible in these layers [Parker 1975, Moreno-Insertis 1983] and therefore this result is not unexpected. While Brandenburg (2005) has conjectured that the near-surface shear layer may be able to power a large-scale dynamo, this remains to be convincingly demonstrated in the context of a Babcock-Leighton dynamo.
In the second scenario with a shallow meridional flow, we allow magnetic flux tubes to buoyantly erupt from 0.71 , i.e. from base of the convection zone. In this case we get periodic solutions but analysis of the butterfly diagrams (taken both at the base of SCZ and near solar surface) shows that the toroidal field belts have almost symmetrical poleward and equatorward branches with no significant equatorward migration (see Fig. 2). Moreover, as already noted by Guerrero, G. & de Gouveia Dal Pino (2008), the solutions with shallow meridional flow always display quadrupolar parity in contradiction with solar cycle observations. Clearly, a shallow flow poses a serious problem for solar cycle models.
We now introduce both radial and latitudinal turbulent pumping in our dynamo model to explore whether a Babcock Leighton flux transport dynamo can operate with meridional flow which is much shallower than previously assumed; we also extend this study to the scenario where meridional flow is altogether absent.
The turbulent pumping profile is determined from independent MHD simulations of solar magnetoconvection [Ossendrijver et al. 2002, Käpylä et al. 2006]. Profiles for radial and latitudinal turbulent pumping ( and ) are:
The value of and determines the amplitude of and respectively. Fig. 3 (top and bottom plot) shows that radial pumping speed (dashed lines) is negative throughout the convection zone corresponding to downward advective transport and vanishes below . The radial pumping speed is maximum near the poles and decreases towards the equator. Fig. 3 (top and bottom plot) shows that the latitudinal pumping speed (solid lines) is positive (negative) in the convection zone in the northern (southern) hemisphere and vanishes below the overshoot layer. This corresponds to equatorward latitudinal pumping throughout the convection zone.
Dynamo simulations with turbulent pumping generate solar-like magnetic cycles (Fig. 4 and Fig. 5). Now the toroidal field belt migrates equatorward, the solution exhibits solar-like parity and the correct phase relationship between the toroidal and poloidal components of the magnetic field (see Fig. 5). Evidently, the coupling between the poloidal source at the near-surface layers with the deeper layers of the convection zone where the toroidal field is stored and amplified, the equatorward migration of the sunspot-forming toroidal field belt and correct solar-like parity is due to the important role played by turbulent pumping. We note if the speed of the latitudinal pumping in on order of 1.0 ms the solutions are always of dipolar parity irrespective of whether one initializes the model with dipolar or quadrupolar parity. Interestingly, the latitudinal migration rate of the sunspot belt as observed is of the same order.
The above result begs the question whether flux transport solar dynamo models based on the Babcock-Leighton mechanism that include turbulent pumping can operate without any meridional plasma flow. To test this, we remove meridional circulation completely from our model and perform simulations with turbulent pumping included. We find that this model generates solar-like sunspot cycles with periodic reversals (see Fig. 6) which are qualitatively similar to the earlier solution with both pumping and shallow meridional flow. However, we find that the surface magnetic field dynamics related to polar field reversal is limited to within 60 degrees latitudes in both the hemispheres. At higher latitudes (near the poles) the field is very weak and almost non-varying over solar cycle timescales. This is expected if the surface magnetic field dynamics is governed primarily by diffusion. Based on this result, we argue that this scenario of non-existent meridional circulation is not supported by current observations of surface dynamics which seem to suggest that the fields do migrate all the way to the poles.
Two important characteristics associated with the solar magnetic cycle are its amplitude and periodicity. While the periodicity of the cycle predominantly depends on the recycling time between toroidal and poloidal field, its amplitude depends on a variety of factors including dynamo source strengths and relative efficacy of transport timescales with respect to the turbulent diffusion timescale. We explore the dependency of the solar cycle period and amplitude to variations in the transport coefficients to explore the subtleties of the interplay between diverse flux transport processes. Figure 7 shows the dependency of cycle amplitude and periodicity on different velocity components like turbulent pumping and (shallow) meridional flow. A parametric analysis of this dependency yields the following relationships for cycle period () and cycle amplitude (Amp):
which is gleaned from simulations within the following ranges: , and ; and are radial and latitudinal turbulent pumping speeds, and is the shallow meridional flow speed.
This analysis shows that cycle period and amplitude are both governed by diverse transport coefficients such as meridional flow speed, and radial and latitudinal components of turbulent pumping. As radial turbulent pumping carries the flux directly to the base of the convection zone where toroidal field is amplified, increase in the radial turbulent pumping speed leads to a decrease in cycle period. Increasing latitudinal pumping also has a similar effect on period which is similar to what is achieved by increasing meridional flow speed, namely a faster transport through the shear layer and thus shorter cycle periods. The cycle amplitude decreases on increasing the latitudinal pumping or meridional flow speed and this is due to the fact that less time is available for toroidal field induction when it is swept at a faster rate through the rotational shear layers. In surface flux transport models, a similar effect is found but due to a different reason – wherein a faster meridional flow reduces the polar field strength because it takes flux of both polarity and deposits this at the poles (in effect carrying less net flux to the poles); in these simulations with a shallow meridional flow and the double-ring algorithm a similar mechanism could also be contributing to an overall reduction of the field strength. What is interesting to note though is the positive dependence of cycle amplitude on the radial pumping speed. We believe that a faster radial pumping moves the poloidal field down to the generating layers of the toroidal field in the deeper parts of the convection zone faster, thus allowing for less turbulent decay in the poloidal field strength; this eventually results in a stronger poloidal field in the SCZ which generates a stronger toroidal component.
We note that the derived exponents for the cycle period above differ from that determined by Guerrero & de Gouveia Dal Pino (2008). The cycle period in our simulations is more strongly dependent on the latitudinal speed of turbulent pumping and less so on meridional circulation, whereas in Guerrero & de Gouveia Dal Pino (2008) it is the exact reverse. In our model the meridional flow is very shallow and limited to only the top of the SCZ, whereas in the model setup of Guerrero & de Gouveia Dal Pino (2008), the meridional flow penetrates down to about 0.8 ; this we believe makes their dynamo cycle periods more sensitive to meridional flow as compared to latitudinal pumping.
Generally, we find solar-like solutions in a modest turbulent pumping speed range on the order of 1 ms. This parameter study shows that our result are robust to reasonable variations in turbulent pumping coefficients and also points to how the latter may determine solar cycle strength and periodicity.
As there is some uncertainty regarding the exact details of turbulent pumping profiles, we have tested an alternative turbulent pumping profile based on Warnecke et al. (2016). Recent magnetoconvection simulations performed by Warnecke et al. (2016) suggest that radial pumping is downward throughout the convection zone below and upward above , while latitudinal pumping is poleward at the surface and equatorward at the base of the convection zone. Our generated turbulent pumping profiles in the northern hemisphere (defined within ) based on the suggestions of Warnecke et al. (2006) are:
where i.e. the penetration depth of the latitudinal pumping. The amplitudes of and are determined by the value of and respectively. Turbulent pumping profiles in the southern hemisphere are generated by replacing colatitude by . Fig. 8 (the top and middle panels) show that our generated turbulent pumping profiles capture the basic essence of the suggestions made by Warnecke et al. (2016). Our simulations (with shallow meridional flow) and the more complex turbulent pumping profile gleaned from Warnecke et al. (2016) reproduce broad features of the solar cycle and are qualitatively similar to those detailed earlier.
In summary, we have demonstrated that flux transport dynamo models of the solar cycle based on the Babcock-Leighton mechanism for poloidal field generation does not require a deep equatorward meridional plasma flow to function effectively. In fact, our results indicate that when turbulent pumping of magnetic flux is taken in to consideration, dynamo models can generate solar-like magnetic cycles even without any meridional circulation although the surface magnetic field dynamics does not reach all the way to the polar regions in this case. Our conclusions are robust across a modest range of plausible parameter space for turbulent pumping coefficients and also indicate some tolerance for diverse pumping profiles.
These findings have significant implications for our understanding of the solar cycle. First of all, the serious challenges that were apparently posed by observations of a shallow (and perhaps complex, multi-cellular) meridional flow on the very premise of flux transport dynamo models stands resolved. Turbulent pumping essentially takes over the role of meridional circulation by transporting magnetic fields from the near-surface solar layers to the deep interior, ensuring that efficient recycling of toroidal and poloidal field components across the SCZ is not compromised. While these findings augur well for dynamo models of the solar cycle, they also imply that we need to revisit many aspects of our current understanding if indeed meridional circulation is not as effective as previously thought. For example, our simulations indicate that variations in turbulent pumping speeds can be an effective means for the modulation of solar cycle periodicity and amplitude.
It has been argued earlier that the interplay between competing flux transport processes determine the dynamical memory of the solar cycle governing solar cycle predictability [Yeates et al. 2008]. If turbulent pumping is the dominant flux transport process as seems plausible based on the simulations presented herein, the cycle memory would be short and this is indeed supported by independent studies [Karak & Nandy 2012] and solar cycle observations [Muñoz- Jaramillo et al. 2013]. It is noteworthy that on the other hand, if meridional circulation were to be the dominant flux transport process, the solar cycle memory would be relatively longer and last over several cycles. This is not borne out by observations.
Previous results in the context of the maintenance of solar-like dipolar parity have relied on a strong turbulent diffusion to couple the Northern and Southern hemispheres of the Sun [Chatterjee et al. 2004], or a dynamo -effect which is co-spatial with the deep equatorward counterflow in the meridional circulation assumed in most flux transport dynamo models [Dikpati & Gilman 2001]. However, our results indicate that turbulent pumping is equally capable of coupling the Northern and Southern solar hemispheres and aid in the maintenance of solar-like dipolar parity. This is in keeping with earlier, independent simulations based on a somewhat different dynamo model [Guerrero & de Gouveia Dal Pino 2008].
Most importantly, our results point out an alternative to circumventing the Parker-Yoshimura sign rule constraint [Parker 1955, Yoshimura 1975] in Babcock-Leighton type solar dynamos that would otherwise imply poleward propagating sunspot belts in conflict with observations. Brandenburg et al. (1992) and Ossendrijver et al. (2002) had already pointed towards this possibility in the context of mean-field dynamo models. While a deep meridional counterflow is currently thought to circumvent this constraint and force the toroidal field belt equatorward, our results convincingly demonstrate that the latitudinal component of turbulent pumping provides a viable alternative to overcoming the Parker-Yoshimura sign rule in Babcock-Leighton models of the solar cycle (even in the absence of meridional circulation).
We note however that our theoretical results should not be taken as support for the existence of a shallow meridional flow, rather we point out that flux transport dynamo models of the solar cycle are equally capable for working with a shallow or non-existent meridional flow, as long as the turbulent pumping of magnetic flux is accounted for; this is particularly viable when turbulent pumping has a dynamically important latitudinal component. Taken together, these insights suggest a plausible new paradigm for dynamo models of the solar cycle, wherein, turbulent pumping of magnetic flux effectively replaces the important roles that are currently thought to be mediated via a deep meridional circulation within the Sun’s interior. Since the dynamical memory and thus predictability of the solar cycle depends on the dominant mode of magnetic flux transport in the Sun’s interior, this would also imply that physics-based prediction models of long-term space weather need to adequately include the physics of turbulent pumping of magnetic fields.
We acknowledge the referee of this manuscript for useful suggestions. We thank Jörn Warnecke for helpful discussions related to the adaptation of the turbulent pumping profile from Warnecke et al. 2016. We are grateful to the Ministry of Human Resource Development, Council for Scientific and Industrial Research, University Grants Commission of the Government of India and a NASA Heliophysics Grand Challenge Grant for supporting this research.
- Augustson et al. 2015 Augustson, K., Brun, A.S., Miesch, M. & Toomre, J. 2015, ApJ, 809, 149
- Babcock 1961 Babcock, H.W. 1961, ApJ, 133, 572
- Brandenburg et al. 1992 Brandenburg, A., Moss, D. & Tuominen, I. 1992, in The Solar Cycle, ed. K. L., Harvey, ASP Conf. Ser., 27, 542
- Brandenburg et al. 1996 Brandenburg, A., Jennings, R.L., Nordlund, A., Rieutord, M., Stein, R.M. & Tuominen, I. 1996, JFM, 306, 325
- Brandenburg 2005 Brandenburg, A. 2005,ApJ, 625, 539
- Belusz et al. 2015 Belusz, B., Dikpati, M. & Forgacs-Dajka 2015, ApJ, 806, 169
- Charbonneau & Dikpati 2000 Charbonneau, P., Dikpati, M. 2000, ApJ, 543, 1027
- Charbonneau 2010 Charbonneau, P. 2010, Living Rev. Sol. Phys., 7, 3
- Chatterjee et al. 2004 Chatterjee, P., Nandy, D., & Choudhuri, A.R. 2004, A&A, 427, 1019
- Choudhuri et al. 1995 Choudhuri, A. R., Schüssler, M., & Dikpati, M. 1995, A&A, 303, L29
- Choudhuri et al. 2004 Choudhuri, A.R., Chatterjee, P., & Nandy, D. 2004, ApJL, 615, L57
- Choudhuri & Karak 2012 Choudhuri, A.R. & Karak, B.B. 2012, PhRvL, 109, 171103
- Dasi-Espuig et al. 2010 Dasi-Espuig, M. et al. 2010, A&A, 518, A7
- Dikpati & Charbonneau 1999 Dikpati, M., & Charbonneau, P. 1999, ApJ, 518, 508
- Dikpati & Gilman 2001 Dikpati, M., & Gilman, P.A. 2001, ApJ, 559, 428
- Dikpati et al. 2002 Dikpati, M., Corbard, T., Thompson, M.J. & Gilman, P.A. 2002, ApJ, 575, L41
- Dorch & Nordlund 2001 Dorch, S.B.F. & Nordlund, A. 2001, A&A, 365, 562
- Durney 1997 Durney, B.R. 1997, ApJ, 486, 1065
- Gilman 1983 Gilman, P.A. 1983, ApJ, 53, 243
- Glatzmaier 1985 Glatzmaier, G.A. 1985, ApJ, 291, 300
- Guerrero & de Gouveia Dal Pino 2007 Guerrero, G. & de Gouveia Dal Pino, E.M. 2007, A&A, 464,341
- Guerrero & de Gouveia Dal Pino 2008 Guerrero, G. & de Gouveia Dal Pino, E.M. 2008, A&A, 485, 267
- Hale 1908 Hale, G.E. 1908, ApJ, 133, 572
- Hale et al. 1919 Hale, G.E., Ellerman, F., Nicholson, F.R. & Joy, A.B. 1919, ApJ, 49, 153
- Hathaway et al. 2003 Hathaway, D. H., Nandy, D., Wilson, R. M., & Reichmann, E. J. 2003, ApJ, 589, 665.
- Hathaway & Rightmire 2010 Hathaway, D.H. & Rightmire, L. 2010, Science 327, 1350
- Hathaway & Rightmire 2011 Hathaway, D.H. & Rightmire, L. 2011, ApJ 729, 80
- Hathaway 2012 Hathaway, D.H. , 2012, ApJ, 760, 84
- Hazra & Nandy 2013 Hazra, S. & Nandy, D. 2013, In ASI Conf. Ser. 10, International Symposium on Solar Terrestrial Physics, ed. N. Gopalswamy, S.S. Hasan, P.B. Rao & P. Subramanian, 115
- Hazra et al. 2014 Hazra, G., Karak, B.B. & Choudhuri, A.R. 2014, ApJ, 782, 93
- Hazra et al. 2014 Hazra, S., Passos, D. & Nandy, D. 2014, ApJ, 789, 5
- Hazra 2016 Hazra, S. 2016, PhD Thesis, IISER Kolkata, arXiv:1604.00563
- Jiang et al. 2013 Jiang, J., Cameron, R. H., Schmitt, D. & Isik, E. 2013, A&A, 553, A128
- Käpylä et al. 2006 Käpylä, P.J., Korpi, M.J., Tuominen, I. 2006, Astron Nachr, 327, 884
- Karak & Nandy 2012 Karak, B.B., & Nandy, D. 2012, ApJL, 296, 1671
- Leighton 1969 Leighton, R.B. 1969, ApJ, 156, 1
- Mitra-Kraev & Thompson 2007 Mitra-Kraev, U., & Thompson, M.J. 2007, Astronomische Nachrichten, 328,1009
- Moreno-Insertis 1983 Moreno-Insertis, F. 1983, A&A, 122, 241
- Moreno-Insertis et al. 1992 Moreno-Insertis, F., Schüssler, M., Ferriz-Mas, A., A&A, 264, 686
- Muñoz-Jaramillo et al. 2009 Muñoz-Jaramillo, A., Nandy, D., & Martens, P. C. H. 2009, ApJ, 698, 461
- Muñoz-Jaramillo et al. 2010 Muñoz-Jaramillo, A., Nandy, D., Martens, P.C.H. & Yeates, A.R. 2010, ApJL, 720, L20
- Muñoz- Jaramillo et al. 2013 Muñoz-Jaramillo, A., Dasi-Espuig, M., Balmceda, L.A., Kirova, N.A., & Deluca, E.E. 2013, ApJL, 767, L25
- Nandy & Choudhuri 2001 Nandy, D., & Choudhuri, A. R. 2001, ApJ, 551, 576
- Nandy & Choudhuri 2002 Nandy, D., & Choudhuri, A. R. 2002, Science, 296, 1671
- Nandy et al. 2011 Nandy, D., Muñoz-Jaramillo, A., & Martens, P. C. H. 2011, Nature, 471, 80
- Ossendrijver et al. 2002 Ossendrijver, M., Stix, M., Brandenburg, A., & Rüdiger, G. 2002, A&A, 394, 735
- Ossendrijver 2003 Ossendrijver M 2003, A&A Rev, 11, 287
- Parker 1955 Parker, E.N. 1955, ApJ, 122, 293
- Parker 1975 Parker, E.N. 1975, ApJ, 202, 523
- Passos et al. 2014 Passos, D., Nandy, D., Hazra, S., & Lopes, I. 2014, A&A, 563, A18
- Passos et al. 2015 Passos, D., Charbonneau, P. & Miesch, M. 2015, ApJL, 800, L18
- Pipin & Seehafer 2009 Pipin, V.V. & Seehafer, N. 2009, ApJ, 493, 819
- Racine et al. 2011 Racine, E., Charbonneau, P., Bouchat, A. & Smolarkiewicz, P.K. 2011, ApJ, 735, 46
- Rogachevskii et al. 2011 Rogachevskii, I., Kleeorin, N., Käpylä, P.J. & Brandenburg, A. 2011, Phys Rev E, 84, 056314
- Simard et al. 2016 Simard, C., Charbonneau, P. & Dube, C. 2016, ArXiv 1604.01533v1
- Tobias et al. 2001 Tobias, S.M., Brummell, N.H., Clune, T.L. & Toomre, J. 2001, ApJ, 549, 1183
- Warnecke et al. 2016 Warnecke, J., Rheinhardt, M., Käpylä, M.J, Käpylä, P.J. & Brandenburg, A. 2016, ArXiv 1601.03730
- Yoshimura 1975 Yoshimura, H. 1975, ApJ, 201, 740
- Yeates et al. 2008 Yeates, A.R., Nandy, D. & Mackay, D.H. 2008, ApJ, 673, 544
- Zhao et al. 2013 Zhao, J., Bogart, R.S., Kosovichev, A.G., Duvall, Jr., T.L. & Hartlep, T. 2013, ApJL, 774, L29