[
Abstract
Spherical shell dynamo models based on rotating convection show that the flow within the tangent cylinder is dominated by an offaxis plume that extends from the inner core boundary to high latitudes and drifts westward. Earlier studies explained the formation of such a plume in terms of the effect of a uniform axial magnetic field that significantly increases the lengthscale of convection in a rotating plane layer. However, rapidly rotating dynamo simulations show that the magnetic field within the tangent cylinder has severe lateral inhomogeneities that may influence the onset of an isolated plume. Increasing the rotation rate in our dynamo simulations (by decreasing the Ekman number ) produces progressively thinner plumes that appear to seek out the location where the field is strongest. Motivated by this result, we examine the linear onset of convection in a rapidly rotating fluid layer subject to a laterally varying axial magnetic field. A cartesian geometry is chosen where the finite dimensions mimic in cylindrical coordinates. The lateral inhomogeneity of the field gives rise to a unique mode of instability where convection is entirely confined to the peakfield region. The localization of the flow by the magnetic field occurs even when the field strength (measured by the Elsasser number ) is small and viscosity controls the smallest lengthscale of convection. The lowest Rayleigh number at which an isolated plume appears within the tangent cylinder in spherical shell dynamo simulations agrees closely with the viscousmode Rayleigh number in the plane layer linear magnetoconvection model. The lowest Elsasser number for plume formation in the simulations is significantly higher than the onset values in linear magnetoconvection, which indicates that the viscous–magnetic mode transition point with spatially varying fields is displaced to much higher Elsasser numbers.
The localized excitation of viscousmode convection by a laterally varying magnetic field provides a mechanism for the formation of isolated plumes within the Earth’s tangent cylinder. The polar vortices in the Earth’s core can therefore be nonaxisymmetric. More generally, this study shows that a spatially varying magnetic field strongly controls the structure of rotating convection at a Rayleigh number not much different from its nonmagnetic value.
Magnetic confinement of rotating convection]Confinement of rotating convection by a laterally varying magnetic field B. Sreenivasan and V. Gopinath]Binod Sreenivasan and Venkatesh Gopinath
1 Introduction
The Earth’s dynamo is powered by thermochemical convection occurring in its liquid iron outer core. The rapid rotation of the Earth’s core divides convection into two regions, inside and outside the tangent cylinder. The tangent cylinder is an imaginary cylinder that touches the solid inner core and cuts the core surface at approximately latitude 70. The tangent cylinder may be approximated by a rotating plane layer in which convection takes place under a predominantly axial () magnetic field and gravity pointing in the downward direction. Strongly ageostrophic motions are needed to transport heat from the inner core boundary to the core–mantle boundary inside the tangent cylinder (Jones 2015), which implies that nonmagnetic convection inside the tangent cylinder starts at a Rayleigh number much higher than the threshold value for convection outside it. At onset, thin upwellings and downwellings aligned with the axis develop along which the vorticity changes sign, in line with the classical picture of rotating RayleighBénard convection in a plane layer.
Observations of secular variation of the Earth’s magnetic field suggest that there are anticyclonic polar vortices in the core (Olson & Aurnou 1999; Hulot et al. 2002). Whereas core flow inversion models support the presence of axisymmetric toroidal motions, it is not clear that relatively smallscale, nonaxisymmetric motions would be dominant (see Holme 2015, and references therein). Nonmagnetic laboratory experiments that simulate the tangent cylinder region (Aurnou et al. 2003) show an ensemble of thin helical plumes extending from the inner core boundary to high latitudes. A largescale anticyclonic zonal flow in the polar regions is suggested, the likely cause of which is a thermal wind (Pedlosky 1987; Sreenivasan & Jones 2006a):
(1.0) 
where is the angular velocity about the rotation axis , is the acceleration due to gravity, is the thermal expansion coefficient and is the temperature perturbation. Equation (1) is obtained by taking the curl of the momentum equation in the inertiafree, inviscid limit. If the polar regions are slightly warmer than the equatorial regions due to a buildup of light material, (1) predicts an axisymmetric anticyclonic circulation near the poles. It remains to be seen whether magnetic laboratory experiments (Aujogue et al. 2016) would support the presence of smallscale, nonaxisymmetric polar circulation.
Numerical simulations of the geodynamo (e.g. Sreenivasan & Jones 2005) present a different picture from nonmagnetic experiments in that the structure of convection within the tangent cylinder is often dominated by an offaxis plume that carries warm fluid from the inner core surface to highlatitude regions (greater than latitude 70). This type of convection also produces a polar vortex because the radially outward flow at the top of the plume interacts with the background rotation (via the Coriolis force) to generate a nonaxisymmetric, anticyclonic flow patch. For supercritical convection in the Earth’s tangent cylinder, one or more strong plumes may be produced which continuously expel magnetic flux from high latitudes, a process that may be inferred from observation of the rather weak flux in this region (Jackson et al. 2000) or the location of the persistent magnetic flux patches just outside the tangent cylinder (Gubbins et al. 2007). To understand the physical origin of the isolated plumes within the tangent cylinder, §2 focuses on their onset; that is, the regime of their first appearance.
The linear theory of magnetoconvection (Chandrasekhar 1961) predicts that onset in a rotating plane layer occurs either as thin viscously controlled columns or largescale magnetic rolls (see, for example, the structures in figure 5(b) and (d) in §3). Sreenivasan & Jones (2006a) equate the critical Rayleigh numbers for the viscous and magnetic branches of onset to obtain the transition point Elsasser number , where is the Ekman number. (Here, measures the uniform magnetic field strength and is the ratio of viscous to Coriolis forces). If the momentum diffusivity is given a ‘turbulent’ value of the order of the magnetic diffusivity, then , so the viscous–magnetic crossover value is . As this is much less than the observed dipole field at the Earth’s core–mantle boundary, Sreenivasan & Jones (2006a) propose that the offaxis plumes within the tangent cylinder may be in the largescale magnetic mode. However, these arguments rely on the assumption of a uniform axial magnetic field permeating the fluid layer, whereas rapidly rotating dynamo simulations show that the magnetic field has severe axial and lateral inhomogeneities. An important aim of our study is to see whether isolated plumes can form via confinement of viscousmode convection by the naturally occurring, laterally varying magnetic field distribution within the tangent cylinder. This necessitates a comparative study across of plume onset in dynamo simulations (§2).
The onset of convection in threedimensional physical systems has been well understood from onedimensional linear onset theory. Early experiments on the onset of convection in a rotating cylinder containing mercury heated from below and placed in a uniform axial magnetic field (Nakagawa 1957) show that the measured critical Rayleigh number agrees closely with that predicted by onedimensional plane layer onset theory (Chandrasekhar 1961). Subsequently, MHD instabilities have been extensively studied using spatially varying imposed fields of the form in cylindrical coordinates () (Malkus 1967; Soward 1979; Jones et al. 2003) or more complex fields thought to be relevant to rotating dynamos (Fearn & Proctor 1983; Kuang & Roberts 1990; Zhang 1995; Longbottom et al. 1995; Tucker & Jones 1997; Sreenivasan & Jones 2011). In these studies, the backreaction of the mean field on convection via the linearized Lorentz force is the main point of interest, while the generation of the mean field itself is decoupled from this process. Although an incomplete representation of the nonlinear dynamo, linear magnetoconvection provides crucial insights into how the field changes the structure of the flow at onset. For a field that is either uniform or of a lengthscale comparable to the depth of the fluid layer, largescale magnetically controlled convection sets in at small Elsasser numbers (e.g. Zhang 1995; Jones et al. 2003). On the other hand, if the lengthscale of the field is small compared to the layer depth as rapidly rotating dynamo models suggest, the viscous–magnetic mode transition point is displaced to Elsasser numbers or higher (Gopinath & Sreenivasan 2015). The fact that smallscale convection is possible for a wide range of suggests that convection in the Earth’s core may operate in the viscous mode.
Linear stability models that consider variation of the basic state variables along two coordinate axes (Theofilis 2011) resolve the perturbations in two finite dimensions, while the third dimension is of infinite extent. Recent examples of linear onset models where perturbations are resolved in more than one direction include that of doublediffusive convection in a rectangular duct with or without a longitudinal flow (Hu et al. 2012), and quasigeostrophic convection in a cylindrical annulus with the gravity pointing radially outward (Calkins et al. 2013). Models of rotating convection subject to laterally varying magnetic fields are not available. Motivated by the onset of localized convection within the tangent cylinder in nonlinear dynamos, §3 examines onset in a rotating plane layer subject to a laterally varying magnetic field. The finite vertical () dimension and one horizontal () dimension in cartesian coordinates mimic the axial () and azimuthal () dimensions respectively in cylindrical coordinates.
For the classical case of convection under a uniform field of , the scale of convection perpendicular to the rotation axis is significantly increased, and this reduces the Ohmic and viscous dissipation rates. As the work done by the buoyancy force need not be high in order to maintain convection, the critical Rayleigh number is much lower than for nonmagnetic convection (Kono & Roberts 2002). On the other hand, if convection under a spatially inhomogeneous field is viscously controlled so that is much smaller than the axial lengthscale of columns, would be comparable to its nonmagnetic value. The lengthscale of convection thus has implications for the power requirement of a rotating dynamo. An obvious counterpoint to this argument is that of subcritical behaviour, wherein saturated (strongfield) numerical dynamos survive at a Rayleigh number lower than that required for a seed field to grow (e.g. Kuang et al. 2008; Sreenivasan & Jones 2011; Hori & Wicht 2013). The role of the selfgenerated magnetic field in lowering the threshold for convection appears to be consistent with the classical theory of convective onset under a uniform magnetic field (Chandrasekhar 1961) that predicts a significant decrease in critical Rayleigh number from its nonmagnetic value. Numerical dynamo simulations at , however, show that subcritical behaviour is preferred for relatively small magnetic Prandtl numbers rather than for (Morin & Dormy 2009; Sreenivasan & Jones 2011), which indicates that a relatively large ratio of the inertial to Coriolis forces in the equation of motion (measured by the Rossby number , where is the magnetic Reynolds number) may promote subcriticality. Furthermore, Sreenivasan & Jones (2011) show that the depth of subcriticality in rotating spherical dynamos is strongly influenced by the kinematic boundary condition. Noslip boundaries produce dominant columnar convection via Ekman pumping, but give a value that is much smaller than for stressfree boundaries where largescale zonal flows dominate even in slightly supercritical convection. Dynamo calculations at lower would help ascertain whether remains relatively constant or decreases with decreasing Ekman number. Another point of relevance here is that the backreaction of the magnetic field on the columnar flow need not drastically change the transverse lengthscale of convection . Sreenivasan et al. (2014) find that the magnetic field enhances the relative kinetic helicity between cyclones and anticyclones, a process that is essentially independent of . Indeed, saturated spherical dynamo models show that the magnetic field does not appreciably increase from its nonmagnetic value (see, for example, Gopinath & Sreenivasan 2015). In short, the magnetic field can enhance helical fluid motion while preserving the smallscale structure produced by rapid rotation.
Presentday dynamo models mostly operate in parameter regimes where the viscous and Ohmic dissipation rates are comparable in magnitude. If Ohmic dissipation at small lengthscales must dominate over viscous dissipation as in liquid metal magnetohydrodynamic turbulence (Davidson 2001), the magnetic diffusivity must far exceed the momentum diffusivity , so that . Dynamos operating in this regime are very likely turbulent, with a welldefined energy cascade from the energy injection scale to the Ohmic dissipation scale. Geodynamo models typically operate at (e.g. Christensen & Wicht 2007) where the turbulent value of is assumed to match . Low, low models are rare because of the computational effort involved in solving them, but linear magnetoconvection models with spatially varying fields are possible at these parameters. Apart from predicting whether convection in the Earth’s core operates in small scales, these models also give the peak local Elsasser numbers in the core that would still yield a volumeaveraged Elsasser number of order unity. The analysis of rapidly rotating convection under a spatially varying magnetic field is partly motivated by these ideas.
In this study, it is shown that a laterally inhomogeneous magnetic field gives rise to isolated columnar vortices in a rotating plane layer at the onset of convection. This mode of onset is linked to the formation of isolated plumes within the tangent cylinder in convectiondriven dynamos. §2 presents nonlinear dynamo simulations where strongly localized convection appears within the tangent cylinder. Since the critical Rayleigh number is much higher within the tangent cylinder than outside it, supercritical dynamo simulations present the opportunity to visualize the onset of isolated plumes within the tangent cylinder. §3 considers the linear onset of convection in a rotating plane layer of finite aspect ratio subject to a laterally varying axial magnetic field. The onset of localized convection within the tangent cylinder is then interpreted in the light of the linear magnetoconvection results. The main results of this paper are summarized in §4.
2 Nonlinear dynamo simulations
The aim of the spherical shell dynamo simulations is to obtain the regime for onset of localized convection within the tangent cylinder. In the Boussinesq approximation (Kono & Roberts 2002), we consider the dynamics an electrically conducting fluid confined between two concentric, corotating spherical surfaces whose radius ratio is 0.35. The main body forces acting on the fluid are the thermal buoyancy force, the Coriolis force originating from the background rotation of the system and the Lorentz force arising from the interaction between the induced electric currents and the magnetic fields. The nondimensional magnetohydrodynamic (MHD) equations for the velocity , magnetic field and temperature are
(2.0)  
(2.0)  
(2.0)  
(2.0) 
The modified pressure in equation (2.0) is given by , where is the fluid pressure. The velocity satisfies the noslip condition at the boundaries and the magnetic field matches a potential field at the outer boundary. Convection is set up in the shell by imposing a temperature difference between the boundaries. The basic state temperature distribution is given by , where . Equations (2.0)–(2.0) are solved by a dynamo code that uses spherical harmonic expansions in () and finite difference discretization in (Willis et al. 2007). The radial grid points are located at the zeros of a Chebyshev polynomial and are clustered near the boundaries.
The dimensionless parameters in equations (2.0)–(2.0) are the Ekman number , the modified Rayleigh number , Elsasser number , Prandtl number and magnetic Prandtl number , which are defined as follows:
(2.0) 
where is the spherical shell thickness, is the kinematic viscosity, is the density, is the thermal diffusivity, is the magnetic diffusivity, is the gravitational acceleration, is the coefficient of thermal expansion, is the superadiabatic temperature difference between the boundaries, is the angular velocity of background rotation and is the magnetic permeability. The ratio is also called the Roberts number, . The Elsasser number is an output that measures the volumeaveraged strength of the selfgenerated magnetic field in the model. In addition, the Elsasser number based on the measured peak axial () magnetic field within the tangent cylinder is also defined.
Two parameter regimes are considered in this study: (a) , , and (b) , . The Roberts number in both regimes, but at the higher the choice of the larger keeps nonlinear inertia small in the simulation (Sreenivasan & Jones 2006b). Runs for are done with 96 finite difference grid points in radius and a maximum spherical harmonic degree . For , 192 radial grid points and a spectral cutoff of are used. Simulations in both parameter regimes produce strongly dipoledominated magnetic fields.
The focus of attention in this study is on the onset of localized convection within the tangent cylinder. For each dynamo calculation, an equivalent nonmagnetic calculation is done for which only the momentum and temperature equations are stepped forward in time. For , convection starts in the tangent cylinder at in both dynamo and nonmagnetic runs, which indicates that the magnetic field does not alter the critical Rayleigh number for onset. (Convection outside the tangent cylinder sets in at a much lower value of ). At , the tangent cylinder is filled with upwellings and downwellings, although the effect of the magnetic field is visible in the enhanced velocity in plumes (compare figures 1b and c). The magnetic field appears to have mostly diffused in from outside, where dynamo action via columnar convection occurs at much lower Rayleigh number (figure 1a). Close to onset, the magnetic field is not affected much by the plumes, which is why this diffused field is largely homogeneous in the azimuthal () direction. At , convection is strong enough to cause some lateral inhomogeneity in the magnetic field. Patches of form at the base of the convection zone (figure 1d) because of convergent flow at the base of plumes. Dominant upwellings (in red) form over the flux patches while weak convection exists in other areas (figure 1e). At , the highly inhomogeneous field patch that develops at the bottom concentrates convection over it and wipes out convection in the rest of the fluid layer (figure 1g,h). A progressive enhancement of occurs until a threshold field strength is reached, upon which convection is supported only in the strongfield region. At subsequent times, the flow follows the path of the peak magnetic field. The nonmagnetic runs at and show a uniformly distributed axial flow structure (figure 1f,i), which suggests that the confinement of convection in the dynamo is due to the laterally inhomogeneous magnetic field that forms within the tangent cylinder.
For , convection outside the tangent cylinder starts at . Figure 2(b) shows that at , smallscale convection is uniformly distributed inside the tangent cylinder. A comparison of the velocities in the dynamo and nonmagnetic runs (figure 2b,c) shows that the magnetic field intensifies the flow even as its smallscale structure is preserved. The scale of the lateral variation of seen in figure 2(d) () is fixed by the preexisting smallscale velocity field interacting with the field diffusing from outside the tangent cylinder, and this can explain why the transverse lengthscale of is appreciably smaller compared to that at the higher Ekman number. The smallscale patches of in turn concentrate smallscale over them, although convection is still active in other regions. The formation of isolated plumes causes a skewness in , with the peak upwelling velocity being approximately twice the downwelling velocity (figure 2e). As is increased to 438, is strong enough to concentrate a smallscale plume over it and suppress convection elsewhere (figure 2g,h). The marked decrease in the azimuthal lengthscale of the plume with decreasing Ekman number (figures 1h and 2h) suggests that, while the plume is magnetically confined, its width (lengthscale perpendicular to the rotation axis, ) may be controlled by the fluid viscosity. As with the higher Ekman number, the nonmagnetic simulations retain the uniformly distributed axial flow structure from the onset of convection (figure 2c,f,i).
Figure 3 shows horizontal () section plots within the tangent cylinder of and for the two Ekman numbers at onset of the offaxis plume. The nonmagnetic is provided for comparison. As is concentrated near the base of the convection zone, the strong correlation between the magnetic field and the flow is clearly visible by looking at two different sections, for and for . The decrease in plume width at the lower Ekman number is evident by comparing the section plots of at the same (figure 3b,e). It is plausible that the magnetic field locally reduces the Rayleigh number for convection from its nonmagnetic value, upon which the plume finds the location where the field is strongest. A strongly supercritical () dynamo simulation suggests that this idea deserves consideration, as the dominant upwelling in the tangent cylinder continuously migrates to the location of the peak magnetic field in a period of less than 0.1 magnetic diffusion time. Further studies at higher are necessary to obtain the regime where the – correlation within the tangent cylinder completely breaks down.
Table 1 presents the parameters and some key properties of the dynamo simulations performed for the two Ekman numbers. The volumeaveraged Elsasser number ( in all runs) does not give any insight into the onset of the localized plume in the tangent cylinder; on the other hand, the Elsasser number calculated based on the peak value in the tangent cylinder shows a clear increase at plume onset. The field components and are a factor lower than .
180  66.19  0.9  1.34  

186  68.11  0.92  2.94  
190  68.58  0.97  3.49  
385  164.95  1.25  6.76  
415  181.81  1.39  7.92  
438  186.05  1.46  16.81 
A key issue that arises from the nonlinear dynamo simulations is whether the isolated plumes that form within the tangent cylinder are viscously or magnetically controlled. Although it may appear from the simulations at that the magnetic field increases the scale of convection at plume onset (see figure 1b and h), a comparison across Ekman numbers shows that the plume width decreases with decreasing Ekman number (figure 3b,e). In addition, the Rayleigh number for plume onset within the tangent cylinder increases with decreasing Ekman number. These findings suggest that the onset of isolated plumes within the tangent cylinder is controlled by the fluid viscosity.
As the sloping boundaries (top and bottom caps) of the tangent cylinder themselves prevent perfect geostrophy, the critical Rayleigh number and wavenumber at which nonmagnetic convection sets in may not be faithfully reproduced by a plane layer linear onset model. On the other hand, if a laterally varying magnetic field strongly localizes convection in the tangent cylinder, a plane layer magnetoconvection model could be a good approximation for the onset of magnetic convection in the tangent cylinder because the change of boundary curvature across a thin plume is small. Therefore, a study of convective onset in a plane layer under a laterally varying magnetic field is justified. This study is presented in the following section.
3 Linear magnetoconvection model
3.1 Problem setup and governing equations
We consider an electrically conducting fluid in a plane layer of finite aspect ratio, where the vertical () and a horizontal () lengthscale are known and the third direction () is of infinite horizontal extent. The and directions mimic the axial () and azimuthal () directions in cylindrical polar coordinates . The basic state temperature gradient across the layer sets up convection under gravity that acts in the negative (downward) direction. The system rotates about the axis. The fluid layer is permeated by a laterally varying magnetic field of the form
(3.0) 
where is a reference magnetic field strength, , and are constants and is the horizontal lengthscale of the magnetic field. The problem setup is shown in figure 4.
In the Boussinesq approximation, the following linearized MHD equations govern the system:
(3.0)  
(3.0)  
(3.0)  
(3.0) 
The dimensionless parameters , , and in equations (3.0)–(3.0) have the same definitions as in (2.0), except that the spherical shell thickness is replaced by the plane layer depth . The Elsasser number is defined based on the reference magnetic field strength.
By applying the operators and to the momentum equation (3.0) and to the induction equation (3.0) and taking the components of the equations, the behaviour of the five perturbation variables – velocity, vorticity, magnetic field, electric current density and temperature – can be obtained. As the Roberts number is set to unity throughout this study, the onset of convection with an axial magnetic field is expected to be stationary for a wide range of Ekman numbers (Aujogue et al. 2015). Furthermore, this study aims to investigate the structure of convection at onset and seek comparisons with the longtime convection pattern within the tangent cylinder in saturated (quasisteady) nonlinear dynamos. The time dependence of the perturbations is therefore not considered, and solutions are sought in the following form:
(3.0) 
where is the wave number in the direction. After introducing this solution into the governing equations, the following system of differential equations is obtained:
(3.0)  
(3.0)  
(3.0)  
(3.0)  
(3.0) 
where and . The variables , , , are related to the eigenfunctions , , , by the identities
(3.0)  
(3.0) 
where is the horizontal Laplacian.
The stability calculations are performed with both stressfree and noslip boundaries on . Electromagnetic conditions are insulating at the top and bottom, although one set of calculations with mixed (bottom perfectly conducting and top insulating) conditions is done to show that the nature of convective onset is not different from that for insulating walls. As isothermal conditions are maintained for the basic state, the temperature perturbation vanishes at the top and bottom. As the horizontal () direction mimics the azimuthal () direction in cylindrical polar coordinates, periodic conditions are set at the side walls. The boundary conditions on are implemented as follows:
(3.0)  
(3.0)  
(3.0)  
(3.0)  
(3.0) 
3.2 Method of solution and benchmarks
The stationary onset of magnetconvection with a laterally varying field is studied for the parameters , and . The generalized eigenvalue problem , where is solved using Matlab. For the set of equations (3.0)–(3.0), the matrices and their elements are presented in Appendix A. A spectral collocation method that uses Chebyshev differentiation in and Fourier differentiation in is used to resolve the eigenfunctions in two dimensions. For problems with variable coefficient terms (as in this study), the spectral collocation method uses simple matrix multiplication in physical space to treat the terms, whereas a pure spectral method would have resulted in convolution sums for such terms that are algebraically complex (Peyret 2002). The drawback of the collocation method, however, is that the differentiation matrices are dense, making computations memoryintensive (Muite 2010; Hu et al. 2012). The construction of the Fourier and Chebyshev differentiation matrices follows a standard approach, and is given in Appendix A for completeness. For with stressfree boundaries, grid independence is secured with points in and points in , so the nonzero elements of and are of size .
Figure 5 shows the existence of a finite number of equally unstable wavenumbers () at the onset of convection in a plane layer of finite aspect ratio. The exact number and values of the unstable wavenumbers are predictable for a given horizontal lengthscale (Appendix B). As , the number of unstable wavenumbers would become infinite. For a given ,
where is the wavenumber and is the resultant wavenumber. Consequently, the last unstable wavenumber coincides with the critical wavenumber for the classical onedimensional plane layer of infinite horizontal extent (Chandrasekhar 1961). For , nonmagnetic convection gives ; and for , the axial velocity at the unstable wavenumber marked A in figure 5(a) shows 11 pairs of rolls (figure 5b), consistent with the fact that . In a similar way, convection under a uniform axial magnetic field for and at the point B in figure 5(c) produces 2 pairs of rolls (figure 5d) because and .
3.3 Onset of convection under a laterally varying magnetic field
We investigate marginalstate convection in a plane layer of depth and horizontal lengthscale permeated by an inhomogeneous axial () magnetic field of the form (3.0) giving strong localization of the field (see figure 4b). The background value of the field is small compared to its peak value but not zero, in line with the axial field distribution within the tangent cylinder in rapidly rotating dynamo simulations. Figures 6 and 7 give the critical Rayleigh number () and wavenumber () diagrams for this field distribution, with the reference states for nonmagnetic convection and homogeneous magnetic field provided for comparison. The critical wavenumber for the homogeneous field is not shown because its value is not unique, as noted earlier in §3.2 (although the critical resultant wavenumber is unique). For the laterally varying field, follows the same trend as for the homogeneous field, being approximately constant in the largewavenumber viscous branch and then falling steeply in the smallwavenumber magnetic branch. The field inhomogeneity, however, displaces the viscous–magnetic mode transition point to a higher Elsasser number. Changing the mechanical and electromagnetic boundary conditions does not alter the basic properties of the regime diagrams, although the numerical values of and differ from one condition to the other. While the viscous branches for insulating and mixed (top insulating and bottom perfectly conducting) electromagnetic conditions largely overlap, the use of mixed conditions moves the viscous–magnetic transition further to the right (compare the blue and magenta lines in figure 6a and b). Table 2 (for stressfree conditions) and table 3 (for noslip conditions) present selected values of the critical parameters spanning the two branches of instability. A notable property of onset in the magnetic mode is that and are nearly independent of the Ekman number , in agreement with the classical picture of onset in an infinite plane layer under a uniform magnetic field (see, for example, figure 3 in Aujogue et al. 2015). In this regime, the critical temperature gradient for convection is independent of viscosity, and it is the magnetic field via the Lorentz force that breaks the TaylorProudman constraint to set up convection in the fluid layer.
(Stressfree)  (Stressfree)  (Stressfree)  



0  237.91  35.35  0  509.41  76.27  0  1096  164.38  
0.001  237.83  35.3  0.001  509.37  76.25  0.001  1095.97  164.3  
0.01  237.82  35.3  0.01  509.30  76.1  0.01  1095.90  164.3  
0.04  237.57  35.1  0.04  509.01  76.1  0.04  1095.72  164.1  
0.1  237.03  35  0.1  508.43  75.9  0.05  1095.61  164.1  
0.2  236.10  34.7  0.11*  507.73  3.1; 75.9  0.051  1095.22  3.06  
0.22  235.90  34.6  0.12  465.91  3.08  0.1  560.33  3.07  
0.238*  235.70  3.2; 34.6  0.15  374.14  3.09  0.15  375.23  3.08  
0.239  233.44  3.16  0.2  282.61  3.11  0.2  283.37  3.1  
0.3  189.01  3.2  0.3  191.92  3.17  0.3  192.20  3.16  
0.5  120.58  3.35  0.5  121.62  3.35  0.5  121.73  3.35  
0.6  104.29  3.5  0.6  105.01  3.48  0.6  105.08  3.48  
0.8  85.28  3.8  0.8  85.68  3.78  0.8  85.72  3.78  
1.0  75.33  4.15  1.0  75.57  4.13  1.0  75.60  4.13 
(Noslip)  (Noslip)  

0  194.65  31.3  0  440.82  70.5 
0.001  194.61  31.3  0.001  440.81  70.4 
0.01  194.56  31.3  0.01  440.80  70.4 
0.04  193.80  31.1  0.04  439.73  70 
0.1  192.21  30.7  0.1  437.45  69.6 
0.15  190.87  30.4  0.107*  436.91  3.3; 69.6 
0.2  189.50  30  0.11  424.03  3.3 
0.225*  188.78  3.6; 30  0.2  256  3.2 
0.23  185.95  3.6  0.3  179.78  3.3 
0.3  154.15  3.5  0.4  140.42  3.3 
0.5  106.54  3.5  0.5  116.83  3.4 
0.6  93.99  3.6  0.6  101.41  3.5 
Figure 8 presents the neutral stability curves extracted from various points on the regime diagrams for two Ekman numbers. For , the laterally varying magnetic field of small strength forces a unique mode of instability (), even as the vestige of the multiplewavenumber, nonmagnetic solution is visible in the oscillations of the curve (figure 8(a), red line). The amplitude of the oscillations decreases with increasing Elsasser number, and for (blue line in figure 8b), the magnetic mode of onset () appears at the same Rayleigh number as the viscous mode (). For , the magnetic mode overtakes the viscous mode as the most unstable. As is increased further, progressively decreases but remains approximately constant (figure 8c). The viscous–magnetic mode transition at takes place over a very narrow range of Elsasser numbers, with showing viscous onset () and showing magnetic onset (; blue line in figure 8d). (The logarithmic axis scale of figure 8(d) shows the large scale separation between the viscous and magnetic modes clearly). The largewavenumber oscillations still exist, although with the laterally varying magnetic field these are never the most unstable modes.
Figures 9 and 10 show the axial velocity at convective onset for two Ekman numbers. (Both and have identical structures.) The main idea that comes out of this study is that convection takes the form of isolated plumes under a laterally varying magnetic field even when the smallest lengthscale of the flow is controlled by viscosity. For a small field strength , a unique mode of instability develops where convection is concentrated in the neighbourhood of the peak magnetic field at (figures 9a,b). (It has been confirmed that moving the peak location of the imposed field by changing the constant in equation (3.0) also moves the location of convection). The large number of convection cells stacked in the plane points to the viscous mode of onset. Convection here is magnetically confined, yet its smallest lengthscale is viscously controlled. The magnetic field can therefore help overcome the TaylorProudman constraint and set up localized convection while not significantly changing the wavenumber of convection from its nonmagnetic value. It is notable that the fieldinduced localization is more pronounced at the lower Ekman number: for and , the rolls at are appreciably thinner than at , although the imposed magnetic field profile is the same in both cases (figures 9a and b; 9c and d). The formation of thin, yet isolated plumes has direct relevance to convection within the tangent cylinder in rapidly rotating spherical dynamo simulations (§2) where similar structures are noted. As the field strength is increased further to , convection at is still in the largewavenumber viscous branch, whereas convection at has crossed over to the smallwavenumber magnetic branch (figures 10(a) and (b) and table 2). From figures 10(c) and (d) (), it is noted that the smallwavenumber convection at both Ekman numbers is almost identical in structure, consistent with the fact that the critical parameters (, ) are independent of Ekman number beyond the viscous–magnetic transition point (e.g. Aujogue et al. 2015).
Figure 11 shows the components of the Coriolis, Lorentz and viscous forces (denoted by subscripts , and respectively) in the momentum equation (3.0). (In this model, the pressure gradient is not solved for). Here , for which gives onset in the viscous branch (table 2). The component of gives . From the plots of the Lorentz and viscous forces shown on the same colour scale (figure 11b,c), it is inferred that the Lorentz force, whose magnitude is times that of the viscous force, is influential in overcoming the TaylorProudman constraint and setting up convection. Interestingly, makes the dominant contribution to the Lorentz force while is slightly smaller in magnitude than the viscous force (see equation 3.0). At onset in the magnetic branch (), the peak value of is three orders of magnitude larger than that of and only one order of magnitude smaller than that of , which emphasizes the wellknown role of the magnetic field in overcoming the rotational constraint.
The mode of convective onset is reflected in the principal balance of terms in the vorticity equation (3.0). For , closely matches (figure 12a), whereas for , closely matches (figure 12b). The Lorentz force term has a negligible contribution to the balance in the former case while the viscous force term has a negligible effect in the latter. The axis range for the two cases is chosen such that the difference in the lengthscale of convection is clear. Although the largescale magnetic mode of onset in figure 12(b) is well understood, the role of the Lorentz force in setting up convection at the small viscous lengthscale (figure 12a, upper panel) has not received much attention in the literature.
Increasing the background magnetic field intensity relative to the peak value reduces the lateral inhomogeneity of the field. For moderate lateral variation, obtained by progressively increasing in the mean field profile (3.0), convection at onset retains the structure of isolated rolls centered at the location of the peak field. A weak lateral variation (), however, gives rise to a cluster of rolls in the viscous mode whose intensity decays from the centre () towards the periphery. As the lateral variation goes to zero, the solution would tend to that for a homogeneous magnetic field (§3.2).
In summary, a laterally varying magnetic field acting on a rotating fluid layer gives rise to a unique mode of instability where convection follows the path of the peak field. This localized excitation of convection is consistent with the idea that the magnetic field generates helical fluid motion in regions that are otherwise quiescent (Sreenivasan & Jones 2011), although here it is shown that the flow lengthscale at onset could be viscously or magnetically controlled. The critical Rayleigh number for magnetic convection increases with decreasing Ekman number in the viscous branch of onset, whereas it is nearly independent of Ekman number in the magnetic branch. The width of the convection zone decreases with decreasing Ekman number in the viscous branch, whereas it is nearly independent of Ekman number in the magnetic branch.
3.4 Implications for onset of localized convection within the tangent cylinder
From a comparison between the plane layer linear magnetoconvection model and the spherical shell dynamo simulations, the following points are noted:

The Rayleigh number for onset of localized convection (in the form of an isolated plume) within the tangent cylinder ( for and for , with noslip boundaries) lies on the viscous branch of onset in the plane layer magnetoconvection model (– for and – for , with noslip boundaries). This agreement between the plane layer and the spherical dynamo Rayleigh numbers rests on two factors: (a) The dominant magnetic field within the tangent cylinder is axial; and (b) the formation of localized convection within the tangent cylinder is practically unaffected by the curvature of the bounding walls. It is notable that the plane layer model does not predict the critical Rayleigh number for nonmagnetic convection in the tangent cylinder, as the sloping walls allow uniformly distributed convection at a lower Rayleigh number (e.g. for ).

The width of an isolated plume within the tangent cylinder markedly decreases with decreasing Ekman number (figure 3b,e), an effect that is noted only in the viscous branch of onset in the magnetoconvection model (figure 9c,d). Had the plume formed in the magnetic branch of onset, its width should have been nearly independent of the Ekman number (figure 10c,d).
The close agreement between the dynamo and viscousbranch magnetoconvection Rayleigh numbers notwithstanding, the critical Elsasser number in the simulation is much higher than that in the linear magnetoconvection model. For example, in the dynamo simulation at , plume onset occurs for while does not depart much from its nonmagnetic value with noslip boundaries, (see tables 1 and 3). This striking difference in between the spherical shell tangent cylinder and the plane layer model is because of the naturally occurring axial variation of (e.g. figure 2c) whose effect is not considered in the layer model. A recent study (Gopinath & Sreenivasan 2015) shows that a horizontal magnetic field of small axial lengthscale (more applicable to convection outside the tangent cylinder than inside) shifts the viscous–magnetic mode transition point from the classical value for a uniform field, to a much higher value , without a drastic change in the critical Rayleigh number. It is therefore possible that this transition point is displaced to large for an axially varying , allowing an intense, spatially varying magnetic field to exist for . Linear magnetoconvection with field variation along the axial coordinate (in addition to , or both) brings additional complexities owing to the presence of a horizontal field component required to satisfy the divergencefree condition of the mean field, and the presence of a mean flow arising from the Magnetic–Coriolis force balance in the vorticity equation. Nevertheless, such a model is useful in predicting the magnetic field intensity required to produce isolated plumes within the tangent cylinder.
4 Concluding remarks
The new results that have come out of this study are summarized below in points (i)–(iv), together with what was known from earlier studies:

A comparison across Ekman numbers of the onset of an isolated plume within the tangent cylinder in dynamo simulations reveals that (a) the Rayleigh number for plume onset increases with decreasing , and (b) the plume width markedly decreases with decreasing . These results bear the hallmark of viscousmode convection. In addition, the strong – correlation suggests that the plume may seek out the location where the field is strongest.
Earlier studies (§1) proposed that the tangent cylinder plume is in the largescale magnetic mode, in which case its onset Rayleigh number and width should be independent of . However, these studies did not examine plume onset across Ekman numbers.

A laterally varying axial magnetic field localizes convection in a rotating plane layer. The onset of convection takes the form of isolated plumes in regions where the magnetic field is strong. Of particular interest is the onset of localized, smallscale convection (e.g. figure 9(c), for and ), in which case the critical Rayleigh number is not significantly different from that for nonmagnetic convection (figure 6(a), blue line).
Earlier onset models (see §1) did not consider the possibility of a laterally varying mean field locally exciting convection in a rotating layer. These models predicted uniformly distributed convection either in the smallscale viscous mode or in the largescale magnetic mode, with the viscous–magnetic transition occurring at .

The Rayleigh number for plume onset within the tangent cylinder agrees closely with the viscousmode Rayleigh number in the plane layer magnetoconvection model (§3.4). This result suggests that the localized convection within the tangent cylinder is in the viscous mode.

It follows from (iii) that the onset of an isolated plume within the tangent cylinder is approximately linear, even as nonlinear dynamo action exists outside the tangent cylinder.
While it is already known that the onset of pure (nonmagnetic) convection inside the tangent cylinder requires a Rayleigh number much higher than the critical Rayleigh number outside it, our study provides an analogous result for magnetic convection.
Since the confinement of convection occurs in both viscous and magnetic modes of onset and the plume width increases at the mode crossover point (figures 9d and 10b), it might appear that a strong magnetic field within the tangent cylinder would give rise to a plume in the magnetic mode. Notably, however, the plume width does not increase with increasing (and ) in the dynamo regime of relatively strong rotation (; figure 2). This indicates that the effect of rotation on the plume width prevails over the effect of the magnetic field, so that the viscous–magnetic mode transition does not occur. It is hence reasonable to suppose that the laterally varying field within the Earth’s tangent cylinder would strongly localize plumes in the smallscale viscous mode, in turn producing nonaxisymmetric polar vortices. The width of plumes is likely determined by the smallest scale that can be supported against magnetic diffusion in the core; it is plausible that this scale has magnetic Reynolds number .
The nonlinear dynamo simulations in this study are far from the low, low regime thought to exist in the Earth’s core. Simulations in which magnetic diffusion is significantly higher than thermal (and viscous) diffusion would help ascertain whether the critical Rayleigh number for plume formation progressively increases or tends to an asymptotic value as is decreased.
The linear stability analysis makes the simplifying assumption that the imposed field is invariant along one of the horizontal directions () that is chosen to be infinite in extent. The threedimensional linear simulation of the case in figure 9(c) with the field varying in both and shows that the confinement in the direction is merely replicated in the direction with no change in the critical Rayleigh number (). Whereas decomposing the perturbations as waves along involves no loss of generality, it offers two distinct advantages: The critical wavenumber () readily confirms whether convection is viscously or magnetically controlled; and the calculations are far less expensive than threedimensional onset simulations. Calculations for are memoryintensive with the spectral collocation method, but the evolution of pure spectral methods may eventually overcome this limitation.
The confinement of rotating convection at small Elsasser number does not imply that the mean magnetic field strength within the Earth’s tangent cylinder should be small. Rather, a field strength of or higher is plausible (table 1). Consideration of the axial inhomogeneity of the magnetic field would likely displace the viscous–magnetic mode crossover point to much higher , which makes smallscale convection a reality for intense, spatially varying fields. Despite the naturally occurring variation of the axial field inside the tangent cylinder, the Rayleigh number for plume onset matches well with the approximately constant Rayleigh number for viscous magnetoconvection. This indicates that the main effect of the variation is to extend the viscous regime to higher Elsasser numbers.
Binod Sreenivasan thanks the Department of Science and Technology (Government of India) for the award of a SwarnaJayanti Fellowship.
A Matrices for linear magnetoconvection in two dimensions
Here and are identity matrices of size and respectively ( and being the number of points in and ), so that has size .
The differential operator matrices in (A 0) are given by
(A 0)  
(A 0)  
(A 0) 
The abbreviated elements of matrix are as follows:
(A 0)  
(A 0)  
(A 0)  
(A 0)  
(A 0)  
(A 0)  
(A 0) 
A standard approach is followed in the construction of the differentiation matrices (Trefethen 2000; Huang et al. 2006). For , the first order Fourier differentiation matrix for even is,
(A 0) 
and for odd ,
(A 0) 
where .
The transformed GaussLobatto points for in the range are given by
(A 0) 
and the first order Chebyshev differentiation matrix is given by
(A 0) 
B Multiple unstable modes in a plane layer of finite aspect ratio
For stationary convection in a plane layer with periodic boundaries spaced a length apart and stressfree boundaries spaced unit distance apart, the axial velocity has the functional form
(B 0) 
Following Chandrasekhar (1961), this solution is introduced into (3.0)–(3.0) to give the characteristic equation
(B 0) 
where
Since , for marginal state (critical) convection we obtain
(B 0) 
For , onset of convection occurs at and . For , can take 9 integer values: . The corresponding critical wavenumbers are
These 9 modes appear at the onset of convection (figure 5(a), black line).
In the presence of a uniform axial () magnetic field, the form of the function in (B 0) gives the following characteristic equation (Chandrasekhar 1961):
(B 0) 
For and , onset of magnetoconvection occurs at and . For , (B 0) gives . The critical wavenumbers for are therefore,
These 3 modes appear at the onset of magnetoconvection (figure 5(c), blue line).
 Aujogue et al. (2016) Aujogue, K., Pothérat, A., Bates, I., Debray, F. & Sreenivasan, B. 2016 Little Earth Experiment: An instrument to model planetary cores. Rev. Sci. Instrum. 87, 084502.
 Aujogue et al. (2015) Aujogue, K., Pothérat, A. & Sreenivasan, B. 2015 Onset of plane layer magnetoconvection at low Ekman number. Phys. Fluids 27, 106602.
 Aurnou et al. (2003) Aurnou, J., Andreadis, S., Zhu, L. & Olson, P. 2003 Experiments on convection in Earth’s core tangent cylinder. Earth Planet. Sci. Lett. 212, 119–134.
 Calkins et al. (2013) Calkins, M. A., Julien, K. & Marti, P. 2013 Threedimensional quasigeostrophic convection in the rotating cylindrical annulus with steeply sloping endwalls. J. Fluid Mech. 732, 214–244.
 Chandrasekhar (1961) Chandrasekhar, S. 1961 Hydrodynamic and hydromagnetic stability. Clarendon Press, Oxford.
 Christensen & Wicht (2007) Christensen, U. R. & Wicht, J. 2007 Numerical dynamo models. In Treatise on Geophysics (ed. P. Olson), , vol. 8, pp. 245–282. Elsevier.
 Davidson (2001) Davidson, P. A. 2001 An Introduction to Magnetohydrodynamics. Cambridge University Press.
 Fearn & Proctor (1983) Fearn, D. R. & Proctor, M. R. E. 1983 Hydromagnetic waves in a differentially rotating sphere. J. Fluid Mech. 128, 1–20.
 Gopinath & Sreenivasan (2015) Gopinath, V. & Sreenivasan, B. 2015 On the control of rapidly rotating convection by an axially varying magnetic field. Geophys. Astrophys. Fluid Dyn. 109, 567–586.
 Gubbins et al. (2007) Gubbins, D., Willis, A. P. & Sreenivasan, B. 2007 Correlation of Earth’s magnetic field with lower mantle thermal and seismic structure. Phys. Earth Planet. Inter. 358, 957–990.
 Holme (2015) Holme, R. 2015 Largescale flow in the core. In Core Dynamics (ed. P. Olson), Treatise on Geophysics, vol. 8, pp. 91–113. Elsevier B. V.
 Hori & Wicht (2013) Hori, K. & Wicht, J. 2013 Subcritical dynamos in the early Mars’ core: Implications for cessation of the past Martian dynamo. Phys. Earth Planet. Int. 219, 21–33.
 Hu et al. (2012) Hu, J., Henry, D., Yin, XY. & BenHadid, H. 2012 Linear biglobal analysis of RayleighBénard instabilities in binary fluids with and without throughflow. J. Fluid Mech. 713, 216–242.
 Huang et al. (2006) Huang, L., Ng, CO. & Chwang, A. T. 2006 A FourierChebyshev collocation method for the mass transport in a layer of powerlaw fluid mud. Comput. Methods Appl. Mech. Engrg. 195, 1136–1153.
 Hulot et al. (2002) Hulot, G., Eymin, C., Langlais, B., Mandea, M. & Olsen, N. 2002 Smallscale structure of the geodynamo inferred from Oersted and Magsat satellite data. Nature 416, 620–623.
 Jackson et al. (2000) Jackson, A., Jonkers, A. R. T. & Walker, M. R. 2000 Four centuries of geomagnetic secular variation from historical records. Phil. Trans. R. Soc. Lond. A 358, 957–990.
 Jones (2015) Jones, C. A. 2015 Thermal and compositional convection in the outer core. In Core Dynamics (ed. P. Olson), Treatise on Geophysics, vol. 8, pp. 115–159. Elsevier B. V.
 Jones et al. (2003) Jones, C. A., Mussa, A. I. & Worland, S. J. 2003 Magnetoconvection in a rapidly rotating sphere: the weakfield case. Proc. R. Soc. Lond. A 459, 773–797.
 Kono & Roberts (2002) Kono, M. & Roberts, P. H. 2002 Recent geodynamo simulations and observations of the geomagnetic field. Rev. Geophys. 40 (4), 1013.
 Kuang et al. (2008) Kuang, W., Jiang, W. & Wang, T. 2008 Sudden termination of Martian dynamo?: implications from subcritical dynamo simulations. Geophys. Res. Lett. 35, L14284.
 Kuang & Roberts (1990) Kuang, W. & Roberts, P. H. 1990 Resistive instabilities in rapidly rotating fluids: Linear theory of the tearing mode. Geophys. Astrophys. Fluid Dyn. 55, 199–239.
 Longbottom et al. (1995) Longbottom, A. W., Jones, C. A. & Hollerbach, R. 1995 Linear magnetoconvection in a rotating spherical shell, incorporating a finitely conducting inner core. Geophys. Astrophys. Fluid Dyn. 80, 205–227.
 Malkus (1967) Malkus, W. V. R. 1967 Hydromagnetic planetary waves. J. Fluid Mech. 28, 793–802.
 Morin & Dormy (2009) Morin, V. & Dormy, E. 2009 The dynamo bifurcation in rotating spherical shells. Intl J. Mod. Phys. B 23, 5467–5482.
 Muite (2010) Muite, B. K. 2010 A numerical comparison of chebyshev methods for solving fourth order semilinear initial boundary value problems. J. Comput. Appl. Math. 234, 317–342.
 Nakagawa (1957) Nakagawa, Y. 1957 Experiments on the instability of a layer of mercury heated from below and subject to the simultaneous action of a magnetic field and rotation. Proc. R. Soc. A 242, 81–88.
 Olson & Aurnou (1999) Olson, P. & Aurnou, J. 1999 A polar vortex in the Earth’s core. Nature 402, 170–173.
 Pedlosky (1987) Pedlosky, J. 1987 Geophysical Fluid Dynamics. New York, USA: SpringerVerlag.
 Peyret (2002) Peyret, R. 2002 Spectral Methods for Incompressible Viscous Flow. New York: SpringerVerlag.
 Soward (1979) Soward, A. M. 1979 Thermal and magnetically driven convection in a rapidly rotating fluid layer. J. Fluid Mech. 90, 669–684.
 Sreenivasan & Jones (2005) Sreenivasan, B. & Jones, C. A. 2005 Structure and dynamics of the polar vortex in the Earth’s core. Geophys. Res. Lett. 32, L20301.
 Sreenivasan & Jones (2006a) Sreenivasan, B. & Jones, C. A. 2006a Azimuthal winds, convection and dynamo action in the polar regions of planetary cores. Geophys. Astrophys. Fluid Dyn. 100, 319–339.
 Sreenivasan & Jones (2006b) Sreenivasan, B. & Jones, C. A. 2006b The role of inertia in the evolution of spherical dynamos. Geophys. J. Int. 164, 467–476.
 Sreenivasan & Jones (2011) Sreenivasan, B. & Jones, C. A. 2011 Helicity generation and subcritical behaviour in rapidly rotating dynamos. J. Fluid Mech. 688, 5–30.
 Sreenivasan et al. (2014) Sreenivasan, B., Sahoo, S. & Dhama, G. 2014 The role of buoyancy in polarity reversals of the geodynamo. Geophys. J. Int. 199, 1698–1708.
 Theofilis (2011) Theofilis, V. 2011 Global Linear Instability. Annu. Rev. Fluid Mech. 43, 319–352.
 Trefethen (2000) Trefethen, L. N. 2000 Spectral Methods in Matlab, 1st edn. Society for Industrial and Applied Mathematics (SIAM).
 Tucker & Jones (1997) Tucker, P. J. Y. & Jones, C. A. 1997 Magnetic and thermal instabilities in a plane layer: I. Geophys. Astrophys. Fluid Dyn. 86, 201–227.
 Willis et al. (2007) Willis, A. P., Sreenivasan, B. & Gubbins, D. 2007 Thermal coremantle interaction: Exploring regimes for ‘locked’ dynamo action. Phys. Earth Planet. Inter. 165, 83–92.
 Zhang (1995) Zhang, K. 1995 Spherical shell rotating convection in the presence of a toroidal magnetic field. Proc. R. Soc. Lond. A 448, 243–268.