Dynamics of the Tachocline

The solar tachocline is an internal region of the Sun possessing strong radial and latitudinal shears straddling the base of the convective envelope. Based on helioseismic inversions, the tachocline is known to be thin (less than 5% of the solar radius). Since the first theory of the solar tachocline in 1992, this thinness has not ceased to puzzle solar physicists. In this review, we lay out the grounds of our understanding of this fascinating region of the solar interior. We detail the various physical mechanisms at stake in the solar tachocline, and put a particular focus on the mechanisms that have been proposed to explain its thinness. We also examine the full range of MHD processes including waves and instabilities that are likely to occur in the tachocline, as well as their possible connection with active region patterns observed at the surface. We reflect on the most recent findings for each of them, and highlight the physical understanding that is still missing and that would allow the research community to understand, in a generic sense, how the solar tachocline and stellar tachocline are formed, are sustained, and evolve on secular timescales.

1 The solar tachocline 1.1 What is known about the solar tachocline?
After several decades of observations, one of the great achievements of helioseismology remains the inversion of the internal rotational profile of the Sun (Brown et al., 1989).The surface differential rotation prevails through the whole convective zone, as seen in Figure 1 (Thompson et al., 2003).In the radiative interior, the rotation rate is uniform and matches the rotation rate of the convective envelope at mid-latitudes (about 430 nHz near latitude ±35 • ).The transition from a differential to uniform rotation occurs in a transition layer known as the tachocline (Hughes et al., 2007).This region lies just beneath the convection zone and likely has a prolate form, i.e., it is located at ∼ 0.693R ⊙ close the equator, and at ∼ 0.717R ⊙ at higher latitudes (Charbonneau et al., 1999;Basu and Antia, 2003).Its thickness slightly varies according to the definition considered in the helioseismic forward modeling determinations (see, for instance Kosovichev, 1996;Corbard et al., 1999;Elliott and Gough, 1999;Basu and Antia, 2001), yet it is certainly very thin (less than 5% of the solar radius).Recent results considering ∼ 20 years of observations indicate that neither the position or the thickness of the tachocline significantly change with the solar cycle (Basu and Antia, 2019).The tachocline is subject to strong latitudinal and radial shears.The radial shear in the mean azimuthal flow is the strongest (by almost an order of magnitude) and changes sign at latitude ±35 • .Both the radial and the latitudinal shears exhibit significant changes in their amplitudes with the solar cycle.Interestingly, these changes have been found to be different between the cycles 23 and 24 (Basu and Antia, 2019).
This strong large-scale shear drew the attention to the tachocline, and led it to be considered a prominent player in many dynamo models of the Sun known as 'interface dynamos' (Parker, 1993).Indeed, such a strong large-scale shear is a very efficient converter of poloidal fields into toroidal field in mean-field models (see Chapters 1,12 and 15).Nevertheless, the real role played by the tachocline in solar and stellar dynamos is still controversial.On one hand, ZDI observations suggest that the magnetic topology changes from a simple dipolar structure for young TTauri, fully-convective, objects to a complex multipolar topology for objects that have already developed a radiative core (Gregory Fig. 1 Left: Solar differential rotation, as deduced from helioseismology (Schou et al., 1998;Thompson et al., 2003).The surface exhibits a shear layer (often dubbed near-surface shear layer ).The bulk of the convective envelope shows a latitudinal differential rotation with mostly conical iso-contours, except perhaps at the equator where contours are more aligned with the rotation axis.At the base of the convection zone, the rotation becomes progressively solid, creating a sheared layer called the tachocline.Right: Radial cuts of the solar differential rotation profile at latitudes 0 • , 30 • , 45 • , 60 • and 75 • (the dashed line shows the results with a different inversion technique using regularized least-squares, see Thompson et al. (2003) and references therein for more details).et al., 2012).On the other hand, following earlier works from Mohanty and Basri (2003); Reiners et al. (2012), Wright and Drake (2016) unambiguously showed that the X-ray luminosity of fully-convective stars followed the same empirical trend than slightly hotter stars with a radiative core, showing no strong difference in the case a tachocline could be present.In addition, Route (2016) discovered magnetic activity cycles in ultra-cool dwarfs which are fullyconvective, again showing seemingly no strong impact of the existence or nonexistence of a tachocline.In the theoretical field, global numerical simulations (e.g.Strugarek et al. 2017) also showed that a tachocline was actually not necessary to obtain Sun-like cyclic dynamos, and Bice and Toomre (2020) showed that a tachocline was not necessary to generate strong wreaths of toroidal field within the convective core of M-dwarfs.It could nevertheless influence how the dynamo operates.Indeed, numerical experiments show that instabilities in the tachocline could affect the magnetic cycle period (Lawson et al., 2015), by modulating (due to instabilities) the Poynting flux permeating the base of the convective envelope (more on this in Section 3).More recently Guerrero et al. (2019) also reported a similar instability in their numerical endeavors and Brun et al. (2022) found hints in their simulations that taking the tachocline into account could affect the cycle period (see also Beaudoin et al., 2018).

Why is it so thin?
This simple question about the solar tachocline has puzzled and fascinated solar physicists since its discovery.Indeed, its extreme thinness remains one of the key mysteries of our star.Right after its discovery, Spiegel and Zahn (1992) produced the first hydrodynamic theory of the solar tachocline.They analyzed the long-term equilibrium of this thin shear layer, and showed that it should burrow on secular time-scales due to the phenomenon of radiative spreading.In a nutshell, the latitudinal shear with conical iso-contours requires an associated temperature gradient to break the Taylor-Proudman constraint.A meridional circulation then follows, spreading the latitudinal shear downwards on an Eddington-Sweet timescale, defined as where the square of the Brunt-Väisälä frequency is defined as , Ω ⊙ is the mean rotation of the Sun, R the radius of the tachocline and κ the thermal diffusion coefficient.The Eddington-Sweet timescale based on molecular thermal diffusion reaches about 10 20 years for the solar tachocline (Brun and Zahn, 2006).
The radiative spreading (or burrowing) has been revisited multiple times ever since (Garaud and Brummell, 2008;Garaud and Acevedo Arreguin, 2009;Wood andBrummell, 2012, 2018), confirming essentially the robustness of the mechanism highlighted by Spiegel and Zahn (1992).A back-of-the-envelope calculation quickly shows that according to this theory, the solar tachocline should have spread by at least 0.3R ⊙ into the radiative zone at the solar age.Hence, additional mechanisms must be at play to confine the solar tachocline to its extreme thinness smaller than 5% R ⊙ .
Multiple physical scenarios have been proposed to confine the solar tachocline.They can be decomposed into two families: the fast tachocline scenarios, relying on physical mechanisms on timescales of months to years; and the slow tachocline scenarios, relying on physical mechanisms acting on timescales longer than a millennium.The large community debate about the solar tachocline confinement can therefore be largely attributed to the difficulty to adequately model it, because of the multiple time-scales involved in its dynamics and evolution.
The multiple physical mechanism behind confinement scenarios proposed to explain the thinness of the solar tachocline are illustrated in Figure 2 and ordered as a function of their physical timescale (the differentiation between slow and first processes in the tachocline was first described by Gilman 2000, see e.g. the Table 1 there).The various processes can be summarized as follows (we will dive deeper into the last three in Section 2 and Section 3).
• Anisotropic viscosity and turbulent transport.It was originally proposed by Spiegel and Zahn (1992) that turbulent viscosity in the upper tachocline was likely to be anisotropic due to the strong stratification encountered by convective plumes penetrating from above.They modelled it by simply considering a dominant horizontal viscosity, and in this case the tachocline could be confined to its present thickness, provided the effective (turbulent) Prandtl number was large enough.This scenario was opposed by Gough and McIntyre (1998), who argued that turbulence was unlikely to transport angular momentum in a diffusive manner on such a quasi-2D layer, basing their argument on what occurs in the Earth atmosphere.In fact, this debate is not settled today.In a series of analytical papers (Kim, 2005;Leprovost and Kim, 2006;Kim and Leprovost, 2007;Leprovost and Kim, 2009), Kim and Leprovost showed quite clearly that depending on the considered characteristics of the solar tachocline (horizontal shear, radial shear, stratification, magnetism, waves) turbulent transport could actually be either diffusive or anti-diffusive.These different behaviours were confirmed to exist with numerical simulations by Miesch (2003).Indeed, the consideration of low-Pr physics changes the situation when compared to the Earth atmosphere, making the direction and amplitude of the turbulent transport in the tachocline not so obvious to estimate (Garaud, 2020).Furthermore, Tobias et al. (2007) showed that in the presence of large-scale magnetic fields, turbulence could sometimes even not transport angular momentum efficiently at all.This debate was further pushed forward recently by Chen and Diamond (2020), who developed a more complete view of the problem and proposed that the tachocline could be confined by the combined effect of the relaxation of potential vorticity gradient along a resisto-elastic drag originating from tangled magnetic fields in the solar tachocline.Unsurprisingly, turbulent transport therefore still remains one of the most promising and debated origin for the solar tachocline confinement.• Magnetic confinement with a fossil field .Following Spiegel and Zahn (1992), Gough and McIntyre (1998) wrote a very influential paper on the confinement of the solar tachocline, after a seminal idea by Rudiger and Kitchatinov (1997) who proposed a confinement of the tachocline due to the existence of a fossil magnetic field.Arguing that turbulent transport could not behave like a diffusion (see the previous point for a more critical view on this argument), they looked for another possible source for the tachocline confinement.They naturally turned to magnetism.In their seminal paper, they proposed that a magnetic field of fossil origin confined inside the solar radiation zone could balance out the inward radiative spreading of the tachocline through Lorentz forces.Nevertheless, magnetic fields diffuse on Ohmic, secular timescales and an additional mechanism was also required to confine the magnetic field itself, leading to a double-confinement scenario.
The large-scale meridional flow accompanying a priori the burrowing of the tachocline (Wood and Brummell, 2012) was in turn invoked to play this confining role, leading to an overall balance of the solar tachocline.The Gough and McIntyre (1998) scenario possesses an appealing advantage: it offers a truly large-scale theory for the confinement of the solar tachocline, that does not depend on the detailed behavior of complex, turbulent MHD flows at the top of the solar tachocline.Nevertheless, one potential problem with this scenario lies in the possibility for the Sun to possess of a magnetic-field zone between the top of the convection and the confinement layer where the fossil field Lorentz force is at play.This complex scenario was then questioned with various series of numerical simulations (e.g.Strugarek et al. 2011, Acevedo-Arreguin et al. 2013).We will come back to this active debate in Section 2.2, along with the most recent numerical endeavors on this topic.• Magnetic confinement with a dynamo field .The dynamo magnetic field of the Sun visible at its photosphere certainly permeates deep inside the solar interior.At the very least, down to the solar tachocline along with the differential rotation and the turbulent convective motions.Many dynamo theories have historically given the tachocline an important role (see e.g.Chapter 15).If the cyclic magnetic field of the Sun permeates the upper tachocline, could it contribute to confining it to its observed thinness?This interesting question was first addressed by Forgács-dajka and Petrovay (2001) and further developed in a series of papers by the same authors.This scenario historically deeply interested Jean-Paul Zahn, who developed an extension of the initial development by Forgács-dajka and Petrovay (2001).This work was later refined and published in Barnabé et al. (2017) (more on the details in Section 2.1).• Magneto-hydrodynamic instabilities.As a sheared, stratified, magnetized thin layer, one may naturally question the global stability of the solar tachocline, and whether any MHD instabilities could actually help to explain its thinness (e.g. Watson 1981;Gilman and Fox 1997;Charbonneau et al. 1999).A general introduction to this topic can be found in section 8.2 of Miesch (2005).
We now first turn to the magnetic tachocline confinement problem (Section 2) and move on to the detailed description of the magneto-hydrodynamic instabilities in the tachocline and their impact on our understanding of the tachocline of the solar magnetism in Section 3.

Magnetic confinement scenarios of the tachocline
We now give a detailed description of two mechanisms that have been proposed to confine the tachocline to its present thickness, based on the existence of internal magnetic fields either generated by dynamo processes (Section 2.1) or of fossil origin (Section 2.2).

Confinement with a dynamo field
The solar magnetic cycle of 11 years, visible at the photosphere, sustains strong magnetic fields inside our star.But how strong?It is known for a fact that at the surface of the Sun, the dipolar magnetic field has a typical value of a few Gauss over the solar surface (DeRosa et al., 2012).Moreover, magnetic flux concentrations in the form of sunspots can reach magnetic field intensities of hundreds of kG.Deeper down in the solar convective envelope, where the overall dynamo is supposed to be seated, we do not have any precise observational constraints on the field amplitude.Antia et al. (2000) have tried to derive an upper limit for the magnetic field amplitude through helioseismic inversion and found that it should be below 30 T -a somewhat large upper value.Models of flux emergence from the base of the convection zone leading to realistic flux concentration at the photosphere yield a toroidal field of the order of 1 − 10 T (Fan, 2009;Jouve et al., 2018).The possibility of confining the solar tachocline with the oscillating dynamo field penetrating from the upper convective envelope was first suggested by Forgács-dajka and Petrovay (2001), and further explored by Barnabé et al. (2017).A minimal 1D mean-field model for this scenario, which can be mathematically written as where ξ is the depth below the bottom of the convection zone.This model traces the evolution of the azimuthally-averaged poloidal magnetic field a, toroidal magnetic field b, and differential rotation u, over a domain of depth δ below the convection zone.For the sake of simplicity, a is assumed to be constant in latitude, while the latitudinal shape of the toroidal field and of the differential rotation is characterized by a unique wave number k.The fields are coupled through the induction term C A au in Eq. 3, and the Lorentz force C L ab in Eq. 4. This simplified model allows assessing whether the Lorentz force can confine the tachocline by preventing the downward spread of the differential rotation into the underlying radiative interior.In this work, C A , C L and Pm = ν/η (the magnetic Prandtl number) were varied to assess in which parameter regime a confinement was realized.Let's consider first the original scenario proposed by Forgács-dajka and Petrovay (2001) where an oscillating magnetic field is imposed at the top of the domain with a period P cyc = 22 years, with the poloidal and toroidal fields being out of phase by π/2.The two left panels of Figure 3 show the evolution of the toroidal field (upper panel) and of the differential rotation (lower panel).As expected, the toroidal field oscillates due to the boundary forcing and penetrates down to the electro-magnetic skin-depth H skin = 2ηP cyc .The differential rotation u penetrates over about the same depth, as seen in the lower panel (the average over P cyc is shown by the thick black line).In this case, the tachocline is indeed confined, as the initial profile shown by the dashdotted line was actually penetrating deeper down than the converged rotation profile.
The first criticism of the work of Forgács-dajka and Petrovay (2001) was the assumed phase between the toroidal and poloidal magnetic field.Indeed, this phasing can play a critical role to ensure that the Lorentz force is actually strong enough to confine the tachocline.This hypothesis was alleviated in Barnabé et al. (2017) (see Figure 3), who found that the phasing of the two components of the magnetic field actually marginally affect the net differential rotation.In all cases (shown in green), the tachocline remains confined, with a somewhat different radial profile.The conclusion of this work is that this confinement scenario is actually quite robust with respect to the relative phase between the poloidal and toroidal components of the dynamo field.
The second and more important criticism to this scenario lies in its simplicity.Indeed, the original scenario did not consider that the transport in tachocline could be either diffusive or anti-diffusive, or that it could be subject to a radiative spreading (see Section 1.2).Again, this limitation was also alleviated in the work of Barnabé et al. (2017) who identified the region in parameter space where the confinement could be strong enough to oppose such spreading mechanisms.Interestingly, these results show that the dynamo-based confinement scenario seems robust for a large variety of turbulent transport scenarios.
The confinment due to a dynamo generated magnetic field has furthermore been recently confirmed by global simulations (Matilsky et al., 2022).Their model considers a fraction of a stable layer and a convection zone with a thermodynamic stratification that resembles the solar interior up to 0.947R ⊙ and rotates 3 times faster than the Sun which allows to get fast equator and slower poles.If the simulation does not include magnetic field, the angular momentum spreads from the convection zone to the stable interior because of the viscosity.If the simulation includes magnetic field, a dynamo develops generating non-axisymmetric wreaths of toroidal field around the equator (compare left and right panels of Figure 4).This field has a small phase difference with the poloidal magnetic field.This phase difference is instrumental in generating axisymmetric large-scale Maxwell stresses which balances the viscous stresses preventing the inwards transport of angular momentum.There is a skin-depth of penetration of the magnetic field in the stable layers such as in the models of Forgács-dajka and Petrovay (2001); Barnabé et al. (2017).
In such a confinement scenario, the differential rotation exhibits a variability on the magnetic cycle timescale, over the electro-magnetic skin-depth.Such strong variability was not systematically detected in the Sun (Howe et al., 2011).The spatial resolution limit of helioseismic inversions is about 5% of the solar radius (Howe, 2009).If the tachocline is confined by the oscillating dynamo field, this yields an upper limit of about 5 × 10 10 cm 2 /s for the Ohmic diffusivity (Barnabé et al., 2017).
One can summarize the present understanding of the tachocline confinement scenario by an oscillating dynamo magnetic field as follow: • The confinement scenario based on the cyclic dynamo field of the Sun is a robust scenario with respect to the various hypothetical physical ingredients participating in the angular momentum transport • For realistic solar parameters, this confinement scenario works only if the magnetic field is able to permeate sufficiently deep.In other words, if the turbulent Ohmic diffusivity in the upper tachocline exceeds 10 8 cm 2 /s, it is very likely that the cyclic, large scale field of the Sun is able to confine the tachocline to a thickness smaller than 5% of the solar radius • The plausibility of this confinement scenario has been confirmed through global MHD simulations of convective dynamos (Matilsky et al., 2022).Although the parameter regime of these simulations is still far from realistic, it is an encouraging alternative to explore in future simulations.

Confinement with a fossil field
The confinement of the solar tachocline by a fossil field was proposed by Gough and McIntyre (1998).This scenario was based on equilibrium arguments, but lacked the geometrical realism of a 3D modelling.In particular, it relies on the existence of a meridional flow in the upper tachocline which prevents the magnetic field to permeate above.Conveniently, the radiative spreading theory of Spiegel and Zahn (1992) predicts that such a meridional circulation necessarily co-exists as the differential rotation burrows down in the radiative core on secular timescales (Wood and Brummell, 2012).Nevertheless, this scenario has one major difficulty (setting aside the aforementioned inner zone of tachocline assumed to be magnetic field-free).Indeed, such a meridional flow has to go down at some latitude and up at some others to conserve mass.In the upwelling part, the magnetic field is likely to be entrained into the convection zone and therefore does not remain confined.This is what was initially found in numerical simulation of this confinement scenario in Brun and Zahn (2006); Strugarek et al. (2011) using the ASH code.3D visualizations of this phenomenon are shown in Figure 5. Initially (left panel), a fossil magnetic field is assumed to be confined inside the radiative core of the Sun.The model is then evolved, and the magnetic field is twisted in the tachocline (as seen in the middle panel).It eventually finds a way to connect to the convective envelope (third panel) as was just explained.This connection path can prevent the fossil field scenario (Gough and McIntyre, 1998) to establish if the angular momentum exchanged by magnetic torques between the convective and the radiative zones is strong enough.Indeed, the first consequence is that the radiative core is not in solid body rotation anymore, which is of course in disagreement with the observational constraints at our disposal (see Figure 1).The modelling efforts of Strugarek et al. (2011) were revisited by Acevedo-Arreguin et al. ( 2013) because their parameter regime was preventing a true radiative burrowing of Spiegel and Zahn (1992) to occur, by instead favoring viscous burrowing.Indeed, the meridional flow in their case was solely driven by the differential rotation inside the convection zone, and was consequently less energetic and less extended than a meridional circulation accompanying a radiative burrowing (Wood and Brummell, 2018).In this latter work, they unambiguously showed that the differential rotation can spread downward in various ways depending on the adimensional parameter where Pr = ν/κ is the Prandtl number.If σ > 1, the differential rotation spreads viscously.If σ < 1, the spread is controlled by thermal dissipation and a non-negligible meridional circulation develops as the differential rotation burrows.This regime is particularly challenging to model in the solar tachocline, because the Prandtl number is expected to be much less than 10 −3 there.New 3D global simulations of this problem, carried out in the correct σ-parameter regime, indicate that the meridional flow can indeed halt efficiently the outward magnetic spreading in the equatorial region compared to the results of Strugarek et al. (2011).Nevertheless, even in this regime the regions where the meridional flow turns upwards still allow the magnetic field to permeate into the overlying convective envelope and magnetic torques lead to the internal radiative zone to depart from solid-body rotation (Strugarek & Brun, in prep.).Concurrently, another mechanism may contribute to the confinement of the magnetic field in the radiative interior.The turbulent magnetic pumping is a mean-field effect that advects a large-scale field away from the regions where turbulence is strongest (Tobias et al., 1998).According to mean-field theory, the amplitude of the advection is proportional to minus the gradient of the turbulent magnetic diffusivity, γ = − 1 2 ∇η T (Kitchatinov and Rüdiger, 2008).Thus, it is expected that the field is transported from the convection zone downward into the radiative interior.Alternatively, the upwards diffusion of the fossil field could be balanced by its turbulent advection.In mean field models, this effect was verified by Kitchatinov and Rüdiger (2008).They found that total confinement can be found for ∇η T ≳ 10 5 .The recent global simulations of Korre et al. (2021) explore this effect in the Boussinesq approximation and in the non-rotating case.They can progressively increase ∇η T by increasing the Rayleigh (Ra) number in the convection zone, i.e., by making more vigorous convection.In their highest resolution case, with Ra = 10 9 , their simulations evince confinement of the fossil field (see bottom panels in Fig. 5).Yet, a much higher Ra is still needed to achieve the advection values predicted by Kitchatinov and Rüdiger (2008).So far, it is also unclear how rotation may affect these results.In addition, this mechanism relies on the strength of convection in the deep convection zone, which is debated today based on various helioseismic inversions and numerical simulations results (Hanasoge et al., 2015;Hotta et al., 2023).
To conclude, the present understanding of the tachocline confinement scenario by a fossil magnetic field can be summarized as follows: • In the correct burrowing parameter regime of the solar differential rotation (i.e. the radiative spreading), the meridional flow accompanying the burrowing can efficiently halt the outward spread of the hypothetical fossil magnetic field where the flow and the magnetic field are horizontal.• This spread cannot be easily halted in upflow regions of the meridional circulation, and it remains in question to know how the solar convection zone can reach the level of turbulent magnetic pumping needed to compensate the loss of confinement of the fossil field in these regions.• If the spread of the fossil field cannot be halted at some specific latitudes, a magnetic connection between the convective and radiative zone establishes.
The associated angular momentum transport leads to a significant departure of the inner radiative zone from solid-body rotation, that is incompatible with present constraints from helioseismology.
3 Global (M)HD of the solar tachocline After having reviewed the confinement mechanisms of the solar tachocline, we now turn to the more generic MHD phenomenon that could be taking place in this particular region of our Sun, and that can affect the large-scale flows and the magnetism of our star.Differential rotation (DR) inferred from observations in the tachocline (see Fig. 1) can be sustained by the angular momentum transport in latitude and radius, due to the influence of rotation on convection.This differential rotation can itself be unstable to global hydrodynamic disturbances of low longitudinal wave number m.If the tachocline differential rotation is perturbed by some random disturbances or by disturbances that include some longitudinal structures in the form of normal modes, these disturbances either can grow by extracting energy from the differential rotation, or decay away by transferring energy to the differential rotation (Gilman and Cally, 2007;Dikpati, 2012).In the nonlinear regime they can quasiperiodically exchange energies with the differential rotation, and form patterns that have properties of Rossby waves.In particular, the disturbance patterns propagate in the longitudinal direction and have a tilt in latitude.Because of their latitudinal tilts they are able to either extract energy from the DR or can give energy back to the DR.
In addition, a large amount of toroidal magnetic field is expected to exist at the solar tachocline due to its subadiabatic stratification and the strong radial and latitudinal shear.Therefore, toroidal fields in the tachocline are likely closely coupled by MHD processes to the differential rotation and the aforementioned hydrodynamic disturbances, leading to global MHD instabilities.Their behavior depends on the relative energy present in the differential rotation compared to that of the toroidal field (Mestel and Moss, 1983;Charbonneau and MacGregor, 1993;Charbonneau, 2004;Spruit, 2004).The disturbances of toroidal fields, coexisting with differential rotation, can again produce modes which retain some properties of Rossby waves, particularly their phase speed in longitude and latitudinal tilt structures (Gilman and Fox, 1997).But both differential rotation and toroidal fields change these speeds significantly compared to the well known retrograde propagating Rossby waves for a uniformly rotating thin spherical shell of constant thickness, first found by Haurwitz (1940).
Moreover, the presence of magnetic field within the stably stratified radiative interior of stars (the tachocline being subadiabatic qualifies it as such a region) further leads to a large range of MHD instabilities (Spruit, 1999;Jouve et al., 2015).Several MHD instabilities present in stellar radiative interiors have hence been studied in details, for instance that of a large scale poloidal field with closed field lines (Markey and Tayler, 1973;Wright, 1973;Braithwaite, 2007), that of strong toroidal magnetic field (akin to the tachocline situation) (Tayler, 1973;Pitts and Tayler, 1985;Guerrero et al., 2019) and the MRI (Menou and Le Mer, 2006;Masada, 2011;Jouve et al., 2020) to cite only a few emblematic ones.It was for instance shown that the ratio between the poloidal and toroidal magnetic field components or energies is key to maintaining or not the magnetic structure into a stable configuration on long timescales (Wright, 1973;Tayler, 1980;Braithwaite and Spruit, 2004;Braithwaite, 2008Braithwaite, , 2009;;Duez and Mathis, 2010;Monteiro et al., 2023;Fuller and Mathis, 2023).Special attention has been given to the interaction between a large scale field with a large scale mean flow (meridional circulation or differential rotation, Mestel and Moss 1977;Mestel and Weiss 1987).In the situation where a large scale shear is present and maintained, it was shown that a dynamo loop could exist, using the shear as Ω-effect and the generation of non axisymmetric fields and flows through MHD instabilities as α-effect (Spruit, 2002;Cline et al., 2003;Braithwaite, 2006;Zahn et al., 2007;Zaire et al., 2022;Petitdemange et al., 2023).However, such dynamo loop in the radiative interior, with the large scale shear being maintained by external processes such as the convective envelope above or the angular momentum extraction due to stellar wind torque, are not easy to realize as demonstrated in Zahn et al. (2007).Indeed, they required an effective α effect originating from the correlations of the field and flow perturbations, i.e. corresponding to strong electromotive force ϵ = ⟨v ′ × b ′ ⟩.Such correlation can be provided by the MHD instabilities but are not easy to achieve or maintain on long timescales because they can easily be damped when the large scale field aligns with the large scale shear flow or vice and versa, a phenomenon known as Ferraro's law of iso-rotation.Nevertheless, it is remarkable that the joint presence of magnetic fields (either poloidal, toroidal, or both) and large scale shears yields a rich realm of physical mechanisms that are likely acting within stars such as our Sun and in their tachoclines.Since both HD and MHD global processes are likely to operate in the tachocline, we discuss them in more details in the following subsections.

Global HD instabilities
In strictly 2D spherical shells, a DR profile of the solar type1 is unstable to low longitudinal wave numbers only if the equator-to-pole difference in rotation is nearly 30% (Watson, 1981).Nevertheless, the instability can occur for smaller equator-to-pole contrast if (i) Ω is expressed using both the sin2 (θ) and sin 4 (θ) terms, such as Ω = Ω 0 − s 1 sin 2 θ − s 2 sin 4 θ, and/or (ii) 3D effects are considered (e.g. the threshold for the DR instability drops substantially in a shallow-water model when moving from purely 2D to quasi-3D, see Dziembowski and Kosovichev 1987;Charbonneau et al. 1999;Garaud 2001 )In the case (i), the tachocline differential rotation can be still unstable in 2D.In the case (ii), particularly for a layer like the convective overshoot layer of the tachocline that is only slightly subadiabatic on average (thus allowing more vertical motion), differential rotation can be unstable for much smaller amplitude than 30% equator-to-pole differential rotation (Dikpati and Gilman, 2001a).In this latter case, there also appear bulges and depressions in the top surface of the tachocline (corresponding to the tachocline-convection zone interface).Because Coriolis and pressure gradient forces almost balance (as in a geostrophic balance), the perturbation flow becomes clockwise on the bulges, and counterclockwise in the depressions.These also are properties of Rossby waves patterns, whether stable or unstable.
Vertical motions are also correlated with the bulges and depressions, leading to non-zero kinetic helicity.This in turn provides an additional source of alpha effect for solar dynamo models (Dikpati and Gilman, 2001b) that could play a role in parity selection in solar dynamo models (Bonanno et al., 2002;Belucz et al., 2015).

Rossby waves
In the Earth's atmosphere, where Rossby waves were first discovered, the forcing primarily comes from differential heating from the Sun; this is thermal forcing.Coriolis forces are always present in a rotating body, but they do not provide a source of energy to drive the Rossby waves.They only alter the direction of fluid motion in the wave.In the case of the Sun, unlike the Earth's atmosphere, the energy is coming from the mechanical forcing of differential rotation at the top of the tachocline.However, like in the Earth's atmosphere, the Coriolis forces again can alter the direction of plasma flow, but can't provide energy.
Very much like the Earth's Rossby waves, solar Rossby waves were also first predicted theoretically before they were directly observed (Dziembowski and Kosovichev, 1987;Gilman and Fox, 1997;Zaqarashvili et al., 2010).These works found that the disturbances that perturb the tachocline differential rotation form the patterns that have Rossby waves-like properties such as the ones described in Section 3.1.Later on, various observations have indicated the existence of Rossby waves in the Sun (McIntosh et al., 2017;Löptien et al., 2018;Harris et al., 2022).
In fact, even before the role of Rossby waves as perturbation to the global flows and magnetic fields, the idea that the Sun might have an MHD dynamo driven by Rossby waves was explored by Gilman (1968Gilman ( , 1969a,b),b).That model assumed the existence of a 'thermal wind' inside the convection zone, with latitudinal and radial rotation gradients as well as an associated latitudinal specific entropy gradient.Gilman found that Rossby waves and DR can together drive a cyclic dynamo with a period within an order of magnitude of the observed sunspot cycle.The solar tachocline, for which the model assumptions were plausible, was not discovered until almost thirty years later, and observational evidence of Rossby waves in the Sun has been found only over the past five years.
So-called thermal Rossby waves have been discovered theoretically for both the convection zones and radiative interiors (Glatzmaier and Gilman, 1981;Hindman and Jain, 2022).Those may be more like Earth's Rossby waves, in that they contain thermodynamic effects and rely for some of their properties on the outward decline in plasma density, as well as on the spherical geometry of a thick layer as opposed to thin layer that is the tachocline.Thermal Rossby waves are yet to be detected in the Sun.
It is worth noting that in the presence of a latitudinal DR that is stable to perturbations, Rossby waves can theoretically still exist.In that case, their patterns do not have tilts in latitude, and so they cannot extract energy from the reference state.Their phase velocity differs from the uniform rotation case by an amount close to that of the local rotation speed at the latitude where the relative amplitude of the Rossby wave peaks.In an unstable differential rotation, Rossby waves behave similarly, but they have tilts in latitude and a phase speed in longitude close to the rotation speed of the latitude where the unstable mode peaks.These properties are similar for purely 2D waves, and the quasi-3D waves characteristic of hydrodynamic shallow water models (Dikpati and Gilman, 2001a).
With the addition of a toroidal field, the picture changes significantly, depending on the latitudinal profile and strength of the field.We first note that Rossby waves are still nearly hydrodynamic (with retrograde propagation).In addition, a very slow prograde wave starts to appear, and vanishes in the limit of zero toroidal field (Dikpati et al., 2020;Hori et al., 2020).For much stronger fields, magnetic effects overpower Coriolis forces, leading to Alfvén waves propagating with nearly equal speeds in both prograde and retrograde directions.Unstable MHD waves for even moderately strong but latitudinally compact toroidal fields propagate at the rotational speed of the latitude where the toroidal field peaks.Therefore, these waves are likely to be retrograde at high latitudes (poleward of sunspot zones), but could be prograde for a certain range in low latitudes (see, e.g.Dikpati et al., 2018).
It is to be noted that all Rossby waves are inertial waves, however all inertial waves are not Rossby waves.Recently, solar inertial oscillations were studied by Triana et al. (2022) in a 3D spherical shell model of the convection zone containing a homogeneous, incompressible, viscous fluid, using an eigensystem formalism.Retrograde as well as prograde Rossby modes were found, along with other inertial modes, which are distinct from Rossby modes and have radial velocities comparable to horizontal velocities.Rossby modes found by Triana et al. (2022) are similar to those found by Bekki et al. (2022) using a similar model but in a fully compressible regime.

Global MHD instabilities
The addition of a toroidal field to 2D, quasi-3D shallow-water type, and 3D thin-shell type hydrodynamical models has a destabilizing effect, as we noted in the previous section.In that case, instabilities for low longitudinal wave numbers can occur down to very small differential rotations.Which longitudinal wave number is most unstable depends on the strength and profile of the toroidal field.Narrower toroidal bands render higher longitudinal wavenumbers unstable.But very narrow bands of about 2 • latitudinal width are not unstable because they cease to sense the difference in rotation rate across the band.So there is an optimum bandwidth for instability in the context of the solar DR profile (Gilman and Dikpati, 2000).
The existence of unstable toroidal bands at particular latitudes with low longitudinal wave number could provide templates for the longitude and latitude distribution of active regions on the solar surface.The bulges in the tachocline that happen to contain strong toroidal field may indeed be the source for emerging magnetic flux that produces new active regions (Dikpati et al., 2017).
Global HD and MHD instabilities in the tachocline evolve non-linearly due to the interactions among differential rotation, magnetic fields and Rossby waves, which are essentially the disturbance patterns.Nonlinear evolution of these instabilities can produce many interesting features which are discussed in the following two subsections, for pure 2D models and quasi-3D thin-shell shallow-water type models.

Nonlinear evolution of 2D HD and MHD instabilities
In a pure 2D tachocline, nonlinear evolution of the HD instabilities can produce high-latitude jets.Due to the action of Reynolds stresses, the angular momentum can be extracted from the differential rotation specifically from the mid-latitudes, and can be deposited on the poleward sides, creating prograde jets there.Traces of such tachocline jets were observationally detected by helioseismic analysis (Christensen-Dalsgaard et al., 2005).
For a 2D magnetized tachocline, nonlinear evolution of MHD instabilities lead to two interesting features.In the case of broad toroidal fields open into clam-shell pattern (Cally, 2001).Conversely, a narrow toroidal band produces tipping or deformation patterns depending on the field strength.If the toroidal band harbours strong fields, it behaves like a 'steel' ring, and it tips to produce an m = 1 longitudinal pattern (Cally et al., 2003).If the field is weak, the ring deforms and produces m > 1 patterns in longitude.

Nonlinear evolution of quasi-3D instabilities
In quasi-3D thin-shell and shallow-water type models, the formation of highlatitude jets, as well as the clamshell opening of broad toroidal fields and tipping of narrow bands occur too.Nevertheless, additional important nonlinear effects have been found in quasi-3D model simulations, which are not present in 2D models.Due to the presence of potential energy, along with kinetic and magnetic energies in a shallow-water model, nonlinear interactions among differential rotation, toroidal field and Rossby waves can produce the 'Tachocline Nonlinear Oscillations' (TNOs), which arise as a consequence of quasi-periodic exchange of energies among these three reservoirs.These oscillations are produced in both HD and MHD cases, very much like 'nonlinear Orr mechanism' in fluid mechanics.The period of oscillation is typically 6-18 months in the context of the solar interior.This oscillation period interestingly coincides with the observations of Rieger-type periodicity, or 'seasons' of solar activity, characterized by intervals of more intense activity alternating with more quiet periods.This suggests that the origin of such 'seasons' of could be the TNOs (Dikpati et al., 2017).
Figure 6 describes the physics of TNOs in an ideal hydrodynamic system, in which they occur due to nonlinear exchange between kinetic energy of differential rotation and perturbation energy of Rossby waves.We note that the TNO physics in the MHD case is a bit more complex in the MHD case, and defer the interested reader to Figure 6 of Dikpati et al. (2018) for their extended description.An interesting property of TNOs is that they do not damp out in several years' runs in an ideal MHD simulation.Evidence for their existence has been furthermore obtained recently in a fully nonlinear direct numerical simulation of tachocline instabilities using generalized quasilinear approximation by Plummer et al. (2019), where the inclusion of magnetic diffusivity is shown to create an abrupt transition of short-term to long-term oscillations between kinetic and magnetic energies.
Another notable nonlinear effect arises from the action of the magnetic stress (often called the 'mixed stress'), which arises due to the cross-correlation between perturbation magnetic and velocities fields in unstable global MHD modes.This stress can cause extraction of energy from the center of the band, and deposition of energy on the both shoulders of the band.In the nonlinear evolution, this process can continue, and at a point the energy at the bandcenter can be depleted so much as to cause the band to split into two.The poleward part of the split-ring slips fast towards the pole due to curvature stress, which is responsible for causing the famous "poleward-slip" instability of a toroidal ring (Spruit and van Ballegooijen, 1982).A 6 • toroidal ring  a)), because their eastward tilts transport angular momentum away from the equator; consequently the polar region spins up and the equatorial region spins down, resulting into a decrease in pole-to-equator differential rotation (thick blue curve in right frame of panel a), compared to the original (thin dashed curve).Panel (b) shows that tilts become neutral when no more energy can be extracted out of differential rotation.However, the perturbation flow-patterns overshoot from neutral tilt to acquire westward tilts, as shown in panel (c) (middle frame); hence flow-patterns can transport positive angular momentum back from high to low latitudes (see the black arrows on the flow-patterns of the bottom panel), spinning up the equatorial region again.Consequently, differential rotation is restored back, and becomes the source of energy to be extracted by the perturbation flowpatterns.In a nearly dissipationless system this process repeats, very much like nonlinear Orr mechanism in fluid dynamics, and leads to an oscillatory exchange of energies between Rossby waves and tachocline differential rotation (Dikpati et al., 2017).
splitting into two can be seen in this ring-split movie.The toroidal ring is displayed in white arrow-vectors on the colormap, which represents deformation of thin fluid-layer's top-surface (red/orange denotes bulges and blue/dark-blue depressions).The movie shows that the ring is just split into two at t=6.7023 (24.5 days), and the split part of the ring is seen to move poleward at t=7.6556 (30 days).The snapshots are presented in dimensionless time (dimensional time can be obtained by multiplying by 3.65 in the unit of days).
The equatorward part can move slowly to the equator because of the energy getting deposited from the band-center to the equatorward shoulder of the band.This provides a possible physical mechanism for a pair of spots, often emerging at the same longitude but separated in latitude by more than 20degrees (Dikpati et al., 2021).The Figure 7 shows an observational evidence of a pair of spots that could have emerged from split-ring.MHD Rossby waves could have an effect on the solar magnetic activity cycle (Raphaldini and Raupp, 2015) and on space weather (Dikpati and McIntosh, 2020).Here for the general readers, we display a diagram to show how magnetically modified solar Rossby waves can create spatio-temporal patterns on dynamo-generated magnetic fields at the tachocline, very much like the "jet streams" produced by the interaction of Earths' atmospheric Rossby waves with the mean zonal flow.As the "jet streams" steer the terrestrial weather by causing cold wave or draught, respectively transporting cold from high to midlatitudes or heat from low to high latitudes, similarly spatio-temporal magnetic patterns created by the interaction of solar MHD Rossby waves could determine the plausible latitude-longitude locations of flux emergence on the solar surface (see, e.g. Figure 8).We point the reader to the comprehensive review by Zaqarashvili et al. (2021) for more details about Rossby waves and their roles in producing various astrophysical phenomena.

Magnetic buoyancy instability
Concomitantly to global HD and MHD instabilities, perturbations to a tube or layer of magnetic field stored at the tachocline may trigger the buoyancy instability.Over the years, it has been assumed that this instability is responsible for the emergence of magnetic flux from the tachocline to the solar surface, leading perhaps for the formation of sunspots.Furthermore, the turbulence generated by this process may in principle generate a net electromotive force, and therefore an α-effect.
The most relevant aspect of the buoyancy instability is that the conditions for its excitation are easily fulfilled for the magnetic fields developed by the strong radial shear at tachocline levels.The general idea was proposed by Parker (1955) considering a magnetized tube of plasma in pressure equilibrium with its surroundings.If the temperature of the plasma inside the tube is the same as in the exterior, the density inside will be smaller and the tube will become a bubble.This out of equilibrium situation was ever since extensively studied formally as an instability.For instance, Newcomb (1961) used the energy principle of Bernstein et al. (1958) to find the stability condition for interchange modes where the wavenumber along the field tube tends to zero (i.e., there is no deformation of the tube along the direction of the magnetic field).The instability turned out to depend on the vertical gradient of the magnetic field (Thomas and Nye, 1975).The use of the energy principle for an unstable parcel of gas was also performed for non-ideal plasmas (Acheson and Gibbons, 1978).Yet this approach becomes insufficient for more realistic models including rotation and/or shear (Tobias and Hughes, 2004).These aspects, however, were studied through local instability analysis for interchange and, three-dimensional, undular modes (e.g., Gilman, 1970;Acheson and Gibbons, 1978;Hughes, 1985a,b).In general, small magnetic diffusivity and large thermal conductivity have a destabilizing effect (since the instability may develop faster than the magnetic field diffusion, and the positive entropy gradient of a stable layer may decrease, respectively).Rotation was found to have a stabilizing effect, yet the instability depends strongly on the vertical profile of the magnetic field, and rotation-induced shear may therefore contribute to make axisymmetric the easily excited non-axisymmetric buoyant modes.A complete review can be found in (Hughes, 2007).More recently, Gilman (2018b,a) extended the linear stability analysis to the spherical domain, including rotation and vertical shear, and exploring together rotational and buoyancy instabilities.His results indicate that the subadiabatic tachocline is unstable to negative vertical shear, such as it is observed at high latitudes in the Sun.Therefore, the magnetic field diffuses rapidly there and the buoyancy instability becomes inefficient.On the other hand, positive shear, as in the solar lower latitudes, is a stable configuration to rotational instabilities, and the buoyancy instability develops for magnetic fields above ∼ 9 kG.Above this threshold, plenty of unstable modes may become unstable, depending on the amplitude of the magnetic field.
The non-linear evolution of buoyant magnetic fields has been extensively studied through MHD numerical simulations.Initial studies were performed in two-dimensional domains and explored the evolution of single flux tubes (Schuessler, 1979;Moreno-Insertis and Emonet, 1996;Fan et al., 1998), or the development of interchange modes in a magnetic layer (Cattaneo and Hughes, 1988).These studies found the development of the mushroom or umbrella shape in the emerging structures.These structures develop strong vorticity at their edges which ends destroying their coherence.Cattaneo et al. (1990) found in 2D simulations that the addition of twist to the initial magnetic field prevents the development of vorticity and facilitates the emergence of the field.It turned out, however, that in the three-dimensional case the axisymmetric perturbations are unstable only during the first stages of evolution (Matthews et al., 1995;Wissink et al., 2000;Fan, 2001), or remain stable (Fan, 2001), depending on the ambient stratification.The undular 3D modes are dominant, and the bending of the magnetic field lines helps the magnetic structures to remain cohesive during the emergence (Abbett et al., 2001).
Concurrently, the buoyancy instability has been studied under the assumption that below the convection zone the magnetic field is organized in the form of thin flux tubes with radius much smaller than the local pressure scale height.The properties of such tubes change along them but not over their cross section.Thus, the evolution equations may be reduced to one dimension (Spruit, 1981).The simulations commonly consider the full longitudinal domain, and impose non-axisymmetric perturbations, which, after the instability develops, forms arch structures anchored at the base of the tube in the solar convection zone (e.g., Choudhuri, 1989;Ferriz-Mas and Schüssler, 1993) or the subadiabatic tachocline (e.g., Caligari et al., 1995).Because computing the evolution of thin flux tubes is generally inexpensive, a large amount of work has been performed under this approximation, exploring the role of rotation or the conditions for reproducing the tilt angle observed in active regions (e.g., Choudhuri and Gilman, 1987;D'Silva and Choudhuri, 1993).In order to overcome the effect of the Coriolis force, and to emerge with the tilt angles compatible with observations, the initial field strength of the thin flux tubes must be ≥ 10 5 G.A complete description of results and references can be found in Fan (2021).The most recent work on this topic (Weber et al., 2011(Weber et al., , 2013;;Weber and Fan, 2015) explored the evolution of thin flux tubes in the presence of turbulent convection provided by a global simulation (Miesch et al., 2006).The results show an interesting interplay between buoyancy and advection, with initial flux tubes with strength between 40 − 50 kG being the ones that better resemble observational tilt angles, and with rapid rising times.It is worth noting that the thin flux tube approximation is not except of criticism.Different works have shown that the internal dynamics of the tube, neglected Panels (a) and (b): convection pattern and magnetic field lines resulting from a numerical simulation with forced radial shear beneath a convectively unstable layer.If the shear is sufficiently strong, the magnetic elements become buoyant forming Ω-loops.This magnetic field modifies the pattern of convection (adapted from Guerrero and Käpylä, 2011).Panel (c): volume rendering of the buoyant magnetic elements developed from a global dynamo simulation.The iso-surface in blue is in equipartition with the kinetic energy whereas the regions in red have four times more energy (adapted from Nelson et al., 2013).
on the averaged formulation of the thin flux tube approximation, may play a relevant role during the emerging process (Hughes et al. 1998;Hughes andFalle 1998, or more recently Martínez-Sykora et al. 2015).
In the same spirit, Cartesian three-dimensional MHD simulations have explored the evolution of imposed magnetic tubes in the presence of convection (of course this not exactly happening at the stable tachocline but above it, see Fan et al., 2003), the self generation of a magnetic layer by a vertical shear (Vasil andBrummell, 2008, 2009), or even the vertical shear beneath a convectively unstable layer (Guerrero and Käpylä, 2011).Convection strongly affects the evolution of a magnetic structure and creates arched tubes anchored by the downflow lanes with a fat magnetic flux configuration frozen-in with the convectively rising plasma.Under shear conditions more extreme that in the solar tachocline, undular magnetic structures develop and may become buoyancy unstable (Fig. 9(a) and (b)).If the shear is sustained by a mechanical forcing then a self-sustained dynamo can exist because an incoherent, yet different from zero, α-effect can persist (Guerrero and Käpylä, 2011).An oscillatory behavior can also be found if the field is replenished by the boundary condition (Kersalé et al., 2007).
Three dimensional simulations have also been performed in spherical coordinates considering the emergence of magnetic flux tubes in adiabatic (Jouve and Brun, 2007;Fan, 2008;Fournier et al., 2017, e.g.,) and superadiabatic layers (Jouve and Brun, 2009).Even though these works preclude the development of the buoyancy instability in the subadiabatic tachocline, they clearly expose how complex it is for a buoyant structures to cross a strongly stratified layer, specially in the presence of turbulent convection and its associated large scale mean flows (Jouve and Brun, 2009) or of background dynamo magnetic fields (Pinto and Brun, 2013).Reproducing the observed features of active regions, therefore, requires a fine tuning of parameters as the magnitude of the magnetic field, twist, and wave-number of the instability.Note that the origin and magnitude of the twist considered in these models do not have proper justification.Most recently Manek et al. (2018); Manek and Brummell (2021) explored the evolution of twisted magnetic tubes embedded in an adiabatic layer with a background magnetic field.They found that there is a lower limit of the background field that allows the tube emergence, and that magnetic structures where the twist and the background field are aligned rise more easily than tubes with the opposite configuration.Finally, magnetic buoyant elements have been observed in convective global dynamo simulations where the magnetic field is not imposed but develops self-consistently (Nelson et al., 2011(Nelson et al., , 2014;;Fan and Fang, 2014).These works showed that convection may not be a hindrance to magnetic flux emergence but may contribute to the rise of magnetic loops (see Fig. 9(c)).It is not yet clear, however, whether these buoyant magnetic elements may reach the surface levels and form bipolar regions.Another kind of simulations have been performed to explore the conditions of the surface emergence phenomena (e.g., see the review of Cheung and Isobe, 2014, and references therein).These models are far from the scope of this review.

Double band and the extended solar cycle
The examination of torsional oscillations (Guerrero et al. 2016), the coronal green line emission (Tappin and Altrock 2013) and the studies of ephemeral active regions (Martin 2018) provide hints of the existence of extended solar cycles.In these observables, the solar activity cycles also manifest at higher latitudes (around 60-70 degrees) than the classical butterfly diagram for the sunspot cycle.It is generally observed at these latitudes ahead of the start of the new sunspot cycle.This high-latitude component of the solar cycle gradually moves towards the equator, maintaining a separation of about 30 • (see, e.g., the figure 1 of Belucz et al. 2023) with the low-latitude active cycle's branch.When the active cycle's branch reaches near the equator to be ready to annihilate with its opposite-hemisphere counterparts, the high-latitude branch reaches about 30-35 degrees, and could then start producing sunspots of the next cycle.Motivated by these observations, one may consider the possibility that there are double, subsurface magnetic bands with opposite directions in each hemisphere, with a low-latitude toroidal strip and a weak high-latitude  2023) considered a double-band system located within the solar tachocline and consisting of two oppositely-directed toroidal bands in each hemisphere, and study their MHD instability in a quasi-3D MHD shallow-water model.As we make the two oppositely-directed bands in each hemisphere migrate equatorward in the model (Belucz et al., 2023), it is found that the high-and low-latitude bands interact in the same hemisphere when they are at latitude 60 • and 30 • (Figure 10).Here, band-interaction occurs by teleconnection mechanismwe briefly describe the physics below, before elaborating on the MHD of the double-band in the solar tachocline.
Here, band-interaction occurs by a 'teleconnection' mechanism.While teleconnection has long been known in meteorology (Lorenz, 1951) and atmospheric circulation (Blackmon et al., 1984), this concept is relatively new for the Sun (see, e.g., Leamon et al. 2021).Teleconnection is essentially a contemporaneous correlation among various global fluid properties and parameters at remotely separated regions and their influence on each other.Teleconnection was demonstrated in a pioneering paper by Wallace and Gutzler (1981) to be created by planetary waves in the North Atlantic and North Pacific Oscillations.In the Sun, Rossby waves act as perturbations to the unperturbed DR and magnetic fields, nonlinearly exchanging energy among the perturbations and reference states of the system.This long range communication can happen in latitude (even across the equator) as well as in longitude, from one side of the Sun to the other.In the case of dynamo-generated magnetic bands in a solar interior, the basic physics of the global connections was described in figure 4 of Cally et al. (2003).
The main conclusion from the MHD of double-band (Belucz et al., 2023) is that, if an extended solar cycle is represented by the two magnetic bands in the interior, generated by dynamo, namely a magnetic band at sunspot latitudes to represent an active-cycle and an oppositely-directed magnetic band at high-latitude band at 60-degree latitude to represent the extended part of the cycle, the model shows that the high-and low-latitude bands interact and telecommunicate among themselves, but not with their opposite-hemisphere counterparts, until this double-band-system has migrated enough towards the equator.The cross-equatorial communication starts when the low-latitude bands in the north and south hemispheres are closer than 15-degree from the equator, because then the low-latitude bands reach the optimal distance between them to start the teleconnection.This is approximately the phase of the solar cycle when the cycle activity reaches the peak.The low-latitude bands in the North and South start the teleconnection across the equator to indicate the start of the declining phase.

Conclusions and Perspectives
The basic physics of the solar tachocline began to be revealed to solar physicists with the seminal paper of Spiegel and Zahn (1992).Since then, the physics of this intriguing layer has led to many debates in the research community.In particular, its extreme thinness continues to be a source of questioning.The tachocline is a turbulent, doubly-sheared, highly stratified, low-Pr, and magnetized thin layer.It is a promising subject of research for any fluid-dynamics enthusiast.We have seen in the preceding sections that advanced, global 3D models could be developed to address in a major way the dynamics of the tachocline.On the other hands, such models are not resolved enough to model adequately all the small-scale turbulent transport in the solar tachocline.In that respect, several routes can be highlighted to continue deciphering the mysteries of the tachocline: • In the past two years, the study of turbulent transport in the solar tachocline has been again revisited.On one hand, Garaud (2020) has developed in a clear manner the argument that the low-Pr regime of the solar tachocline makes the estimate of the net turbulent transport non-trivial.On the other hand, the effect of magnetic field on the transport is still not fully understood.Following the important work of Tobias et al. (2007), Chen and Diamond (2020) have developed a theoretical argument for a complex tachocline confinement leveraging dynamo-generated tangled magnetic fields.We are today in the position of simulating the highly turbulent solar tachocline in a more and more realistic regime thanks to the massive development of 3D global numerical codes.Dedicated, high-resolution numerical experiment could be built to address frontally this problem, e.g.taking into account the convective penetration in the upper tachocline and the associated magnetic field sustained by dynamo action above.• The role of magnetoshear instabilities in the dynamics of the tachocline and of the solar dynamo is progressing actively recently.These have been found in some global models of the solar interior (Lawson et al., 2015;Guerrero et al., 2019;Jouve et al., 2020), and were interpreted to have different effects on the global dynamo.While the solar dynamo is responsible for generating the spot-producing magnetic fields and their migration in latitude, the global MHD of the tachocline provides a mechanism to create the longitudinal distribution of active regions.As discussed in section 3.3.2,there exists a major opportunity to link via data assimilation the observed spatio-temporal distribution of active regions with the simulated "imprints" of their latitudelongitude locations created by the global tachocline instabilities of the spotproducing magnetic fields (see, e.g.figure 12 of Dikpati and McIntosh 2020).• The tachocline coincides at least in part with the convective-radiative transition in Sun-like stars.Angular momentum can be transported from one region to the other, especially when the star changes its structure or when the external convection zone spins down due to magnetic torques from its wind.In the theory of cool-stars' rotation, the coupling between these two regions plays an important role in the pre-main sequence and early main sequence phases.For instance, Gallet and Bouvier (2013) parameterized such angular momentum exchange timescale to reproduce the rotational distribution observed in open clusters (see also Benbakoura et al. 2019;Ahuir et al. 2020).An assessment of this coupling timescale in realistic models of the tachocline (see e.g.Brun et al., 2011) would therefore be extremely useful today to improve our understanding of the secular evolution of the Sun and solar-type stars.
which is a major facility sponsored by the National Science Foundation under cooperative agreement 1852977, and also acknowledges support from several NASA grants, namely NASA-LWS award number

Fig. 2
Fig.2Schematic of the physical processes at play in the various proposed scenarios to confine the tachocline.The typical timescale of the physical process increases from left to right.The solar differential rotation in the tachocline is illustrated by the background colored zones (see also Fig.1).

Fig. 3
Fig. 3 Left panels: Original confinement scenario of Forgács-dajka and Petrovay (2001) as revisted by Barnabé et al. (2017).The upper panel shows the toroidal field b, and the lower panel the differential rotation u, as a function of depth ξ from left to right.The skin depth H skin = 2ηPcyc is shown by the vertical dashed line.The color and transparency of the line label different phases in one magnetic cycle.Right panel: Study of the influence of the phasing between the two components a and b of the magnetic field.The angle ϕ corresponds to this dephasing.The blue curve is the reference model of Forgács-dajka and Petrovay (2001) that differs from the green models in the boundary condition imposed on a.A reference hydrodynamic model is shown by the red curve.All models considering a dynamo cycle (red and green curve) manage to confine the tachocline regardless of the assumed relative phase ϕ.

Fig. 4
Fig. 4 Magnetic field confinement due to a convective dynamo generated in the global simulations of Matilsky et al. (2022).Panel (a) shows viscous tachocline spreading in a hydrodynamic simulation.Panel (b) shows the formation of a tachocline and a solid body rotating radiative interior in a MHD simulation.

Fig. 5
Fig. 5 Top: 3D visualizations of the tachocline confinement scenario with a fossil field, adapted from Strugarek et al. (2011).The convective motions are shown in the lower cutout of the sphere, with orange-red volumes denoting upward motions and blue volumes downward motions.A few selected field lines are traced from the middle of the radiative zone.The initial confined dipole is shown on the left panel.In the middle panel, the field starts to be sheared in the lower tachocline.In the right panel, the field has permeated through the tachocline and is now connected to the surface.Bottom: confinement of a fossil field through convective turbulent pumping at low (left) and high (right) Rayleigh numbers, adapted from Korre et al. (2021).The black contours show the dipole field lines with values of the streamfunction depicted by the colormap.

Fig. 6
Fig.6Panels a, b, c show schematically how Rossby waves and tachocline latitudinal differential rotation nonlinearly interact.Time goes from top panel to bottom.Energy from unperturbed differential rotation (thick blue curve in left frame of panel(a)) is extracted by perturbation flow-patterns (eastward-oriented red ellipses in the north hemisphere, plotted in middle frame of panel (a)), because their eastward tilts transport angular momentum away from the equator; consequently the polar region spins up and the equatorial region spins down, resulting into a decrease in pole-to-equator differential rotation (thick blue curve in right frame of panel a), compared to the original (thin dashed curve).Panel (b) shows that tilts become neutral when no more energy can be extracted out of differential rotation.However, the perturbation flow-patterns overshoot from neutral tilt to acquire westward tilts, as shown in panel (c) (middle frame); hence flow-patterns can transport positive angular momentum back from high to low latitudes (see the black arrows on the flow-patterns of the bottom panel), spinning up the equatorial region again.Consequently, differential rotation is restored back, and becomes the source of energy to be extracted by the perturbation flowpatterns.In a nearly dissipationless system this process repeats, very much like nonlinear Orr mechanism in fluid dynamics, and leads to an oscillatory exchange of energies between Rossby waves and tachocline differential rotation(Dikpati et al., 2017).

Fig. 7
Fig. 7 Left panel displays active regions' toroid-pattern for July 19, 2000, if all the active regions emerged from one wide toroidal magnetic band from the base of the convection zone; the right panel displays two split toroid-patterns if active regions emerged from spotproducing toroidal rings that have undergone dynamical splitting.North and South toroids are displayed in blue and red respectively (Dikpati et al., 2021).

Fig. 8
Fig. 8 Large-scale Rossby waves propagating around the Sun may play a role in triggering solar storms (yellow and deep red regions indicate where sunspot activity is unlikely and more likely, respectively).These waves are akin to Rossby waves in Earth's atmosphere, which influence weather, including potential heat waves (red arrows) and cold snaps (blue arrows) (schematic figure adapted from SCIENCE NEWS BY AGU, Eos.org 51, vol.101, no. 8, August 2020.)

Fig. 9
Fig. 9Panels (a) and (b): convection pattern and magnetic field lines resulting from a numerical simulation with forced radial shear beneath a convectively unstable layer.If the shear is sufficiently strong, the magnetic elements become buoyant forming Ω-loops.This magnetic field modifies the pattern of convection (adapted fromGuerrero and Käpylä, 2011).Panel (c): volume rendering of the buoyant magnetic elements developed from a global dynamo simulation.The iso-surface in blue is in equipartition with the kinetic energy whereas the regions in red have four times more energy (adapted fromNelson et al., 2013).

Fig. 10
Fig. 10 Snapshots of velocity and magnetic disturbance patterns are shown in the latitudelongitude plane for the asymmetric m=1 modes.Black arrows represent magnetic field vector, white arrows flow vectors and color map the tachocline top surface deformation (red-yellow denote the bulging of top surface, and green-blue the depression).This figure shows the eigenfunctions for the overshoot part of tachocline (G = 0.5) and the bands at 60 • − 30 • (aa), 50 • − 20 • (ab), 40 • − 10 • (ac) and 35 • − 5 • (ad)(Belucz et al., 2023).
80NSSC20K0355, NASA-HSR award number 80NSSC21K1676, subaward from JHU/APL's NASA-HSR grant with award number 80NSSC21K1678 and subaward from Stanford's COFFIES Phase II NASA-DRIVE Center with award number 80NSSC22M0162.Bernadett Belucz's work is supported by Newton International Fellowship of The Royal Society, program number NIF-R1-192417.GG acknowledges support from NASA grants NNX14AB70G, 80NSSC20K0602, and 80NSSC20K1320.Competing Interests.The authors declare they have no conflicts of interest.