physicists have long realized the importance of the earth s magnetic field and that this field originates in the earth s fluid core .
the geodynamo is the process by which the rotating , convecting , electrically conducting molten iron in the earth s fluid core maintains the geomagnetic field against ohmic decay .
the convective processes in the core that produce the velocity fields required for dynamo action are a combination of thermal and compositional convection .
the full dynamo problem requires the examination of the full three dimensional partial differential equations governing convective , incompressible magnetohydrodyamics ( mhd ) . in the past decades
computer models have been used to simulate the actual geodynamo .
see , for example , the review articles by roberts et al on modeling the geodynamo  @xcite and references therein .
however , as glatzmaier et al  @xcite remark , no three dimensional dynamo model has been run at the spatial resolution required to simulate the broad spectrum of turbulence which surely exists in the earth s fluid core .
current computers and numerical methods require the imposition of diffusivities that are several orders of magnitude larger than those which are realistic .
it is therefore reasonable to attempt to gain some insight into the geodynamo by considering a reduction of the full mhd equations to a system that is more tractable , but one that retains many of the essential features of the problem .
the magnetogeostrophic equation proposed by moffatt and loper  @xcite is one such model
. the physical postulates of this model are the following : slow cooling of the earth leads to slow solidification of the liquid metal core onto the solid inner core and releases latent heat of solidification that drives compositional convection in the fluid core .
the arguments for the appropriate ranges of the characteristic length , velocity , and perturbation density are based on these physical postulates .
we first present the full coupled three - dimensional mhd equations for the evolution of the velocity vector @xmath0 , the magnetic field vector @xmath1 and the buoyancy field @xmath2 in the boussinesq approximation and written in the frame of reference rotating with angular velocity @xmath3 .
the physical forces governing this system are coriolis forces , lorentz forces , and gravity . following the notation of moffatt and loper  @xcite we write the equations in terms of dimensionless variables .
the orders of magnitude of the resulting non - dimensional parameters are motivated by the physical postulates of the moffatt and loper model : @xmath4 = - \nabla p +   { e}_2 \cdot \nabla b + r_m b\cdot\nabla b + n^2 \theta   { e}_3 + \eps \delta u \label{eq:1.1 } \\ & r_m \left [ \partial_t b + u \cdot \nabla b - b\cdot \nabla u\right ] =   { e}_2 \cdot \nabla u + \delta b \label{eq:1.2 } \\ & \partial_t \theta + u \cdot \nabla \theta = \eps_\kappa
\delta \theta + s \label{eq:1.3 } \\ & \nabla \cdot u = 0 , \nabla \cdot b = 0 .
\label{eq:1.4}\end{aligned}\ ] ] here @xmath5 is a given smooth function that represents the forcing of the system .
the mathematical statement of the geodynamo problem asks whether there are initial data for the mhd system for which the evolution of the perturbation of the magnetic field @xmath6 grows for sufficiently long time .
this can be interpreted as a question of the existence of instabilities of .    the notation in  is the following .
the cartesian unit vectors are denoted by @xmath7 , and @xmath8 . for simplicity
, we have assumed that the axis of rotation and the gravity @xmath9 are aligned in the direction of @xmath8 .
we have assumed that the magnetic field @xmath1 consists of the sum of an underlying constant field @xmath10 and a perturbation @xmath6 .
our choice of @xmath11 as the underlying magnetic field is consistent with the models where the magnetic field is believed to be predominantly toroidal ( cf .
the dimensionless parameters are @xmath12 where the typical velocity scale @xmath13 is chosen to be @xmath14 here @xmath15 is a typical amplitude of the buoyancy @xmath16 and @xmath17 is a characteristic length scale .
@xmath18 is the sum of the fluid and magnetic pressures , @xmath19 is the kinematic viscosity , @xmath20 is the magnetic diffusivity , and @xmath21 is the molecular diffusivity of the compositional variation that creates an ambient density @xmath22 .
the usual estimate for @xmath20 in the fluid core is @xmath23 .
the values of @xmath19 and @xmath21 are speculative , but likely to be extremely small . for a more detailed discussion of plausible ranges of the physical parameters that are appropriate for the geodynamo
, we refer the reader to the book of ghil and childress  @xcite .
for the regions in the earth s fluid core modeled by moffatt and loper , it is argued that the dimensionless parameters are the following orders of magnitude .
the inverse elsasser number @xmath24 will be of order unity , or somewhat less in the condition of the fluid core , while the rossby number @xmath25 is of order @xmath26 .
the value of @xmath27 is somewhat uncertain since the viscosity @xmath19 is a matter of guess work , but plausibly @xmath28 in the core conditions .
the inverse peclet number @xmath29 is also likely to be extremely small , perhaps of the order of @xmath30 .
according to moffatt and loper  @xcite , the magnetic reynolds number @xmath31 is relatively small in the range of @xmath17 considered . with these estimates of the sizes of the parameters the dominant balance of the leading order terms in equations and give the following reduced system : @xmath32 along with the incompressibility condition @xmath33 the _ linear _ system of equations  determine the differential operators that relate the vector fields @xmath0 and @xmath6 with the scalar buoyancy @xmath2 .
these operators encode the vestiges of the physics in the problem , namely coriolis force , lorentz force , and gravity .
vector manipulations of  give the expression @xmath34 u = n^4 ( { e}_3 \cdot \nabla ) \delta ( { e}_3 \times \nabla \theta ) + n^2 ( { e}_2 \cdot \nabla)^2 \nabla \times ( { e}_3 \times \nabla \theta ) \label{eq:1.8}.\end{aligned}\ ] ] the sole remaining nonlinearity in the system comes from the coupling of with the equation for the time evolution of the @xmath2 .
we call this nonlinear active scalar equation the _ magnetogeostrophic _ ( mg ) equation .
it is the subject of this article .
we study the properties and behavior of the active scalar equation in @xmath35 @xmath36 \label{eq : mg:2}\end{aligned}\ ] ] via an examination of the fourier multiplier operator @xmath37 obtained from , which relates @xmath38 and @xmath16 .
more precisely , @xmath39 = ( { \widehat}{m}_j { \widehat}{\theta})^{\vee}\ ] ] for @xmath40 , and the explicit expressions for the components of @xmath41 as functions of the fourier variable @xmath42 are @xmath43 where @xmath44 on @xmath45 we define @xmath46 , for all @xmath47 , since for self - consistency of the model we assume that @xmath16 and @xmath38 have zero vertical mean , i.e. @xmath48 . several important properties of the @xmath49 s are immediately obvious .
firstly , the functions are strongly anisotropic with respect to the dependence on the integers @xmath50 , and @xmath51 .
this is a consequence of the interplay of the three above mentioned physical forces .
secondly , the function @xmath52 is an even function of @xmath50 , and @xmath51 .
thirdly , we have @xmath53 for any @xmath54 ( we use the standard summation convention on repeated indices ) , which is built into the system since the velocity field is divergence free @xmath55
lastly , although the symbols @xmath49 are bounded ( by a constant that depends only on @xmath56 ) in regions of fourier space where @xmath57 , this is not the case in the `` curved '' regions where @xmath58 , and @xmath59 with @xmath60 . in such regions the symbols are unbounded as @xmath61 with @xmath62 that _ is independent of @xmath56 _ , and this bound is sharp .
hence the mg equation falls in the class of abstract active scalar equations ( see  below ) where the relation between the divergence - free velocity field @xmath38 and the scalar @xmath16 is given as a singular operator of order @xmath63 .
these properties of @xmath41 make the mg equation interesting and challenging mathematically , as well as having a clear physical basis in its derivation from the mhd equations .
, @xmath64 , and @xmath65 , with @xmath66 .
the range of the frequencies is @xmath67 $ ] and @xmath68$].,title="fig:",scaledwidth=30.0% ] , @xmath64 , and @xmath65 , with @xmath66 .
the range of the frequencies is @xmath67 $ ] and @xmath68$].,title="fig:",scaledwidth=30.0% ] , @xmath64 , and @xmath65 , with @xmath66 .
the range of the frequencies is @xmath67 $ ] and @xmath68$].,title="fig:",scaledwidth=30.0% ]    in section  [ sec : diffusive ] we state results proved in @xcite for the diffusive mg equation , i.e. , the case when @xmath69 .
we recall that in the fluid core @xmath29 is positive but probably extremely small . in this case
the mg equation is global well posed .
to prove this fact we use the control on @xmath70 , a quantity that is non - increasing along solutions , and which gives a bound on the drift velocity @xmath38 in the space @xmath71 .
the latter information on the drift is `` critical '' for the natural scaling of the _ linear _ drift - diffusion equation @xmath72 indeed , for @xmath73 divergence - free , we proved in @xcite that the solution @xmath16 of becomes hlder continuous in finite time , by appealing to the classical ideas of de giorgi  @xcite , in the spirit of  @xcite .
this result was also established independently in  @xcite .
see also  @xcite , the more recent  @xcite , and references therein , for regularity issues of linear parabolic equations with singular coefficients and divergence - free drift .
the hlder continuity of solutions to the linear equation   may then be bootstrapped to show that the solutions of the dissipative mg equations become smooth in arbitrarily short finite time  @xcite .
having established the global well - posedness of  for any @xmath69 , we can proceed to examine the stability of equilibria .
we show that the existence of an unstable eigenvalue for the linearized equations implies nonlinear ( lyapunov ) instability of the nonlinear mg equation ( see also  @xcite ) .
we then exhibit the existence of an unstable eigenvalue in a particular example , using the method of continued fractions cf .
the magnitude of this eigenvalue is @xmath74 which , in the context of the fluid core , is extremely large .
thus the mg equation can be used to model very rapid exponential growth in time of @xmath16 , and hence , through  , of the perturbation velocity and magnetic fields , i.e. a dynamo type instability .
we emphasize that this `` mg - dynamo '' instability requires the presence of all three forces  coriolis , lorentz , and gravity  because the dimensionless parameter @xmath24 ( defined above in ) must be strictly in the range @xmath75 .    in section  [ sec : inviscid ]
we examine the non diffusive mg equation and set @xmath76 ( cf .
the situation is then dramatically different from the case @xmath69 . without the laplacian to control the unbounded operator @xmath37 , the mg equation is _ ill - posed _ in the sense of hadamard in sobolev spaces .
such a result utilizes in an essential way both the unboundedness of @xmath52 and the fact that it is an even function of @xmath77 .
the proof again uses the existence of an unstable eigenvalue of the linearized equation , however when @xmath76 it is possible to construct such an eigenvalue that is arbitrarily large and it then follows that there is no lipschitz solution map at time @xmath78 ( see also  @xcite ) .    in section  [ sec : frac ]
we describe results for the fractionally diffusive mg equation ( cf .
@xcite ) , i.e. replace the laplacian by @xmath79 , for @xmath80 .
this situation is non - physical but mathematically interesting in the context of the tools required to treat the fractional laplacian .
we now have a dichotomy across the value @xmath81 . in the range @xmath82
the equation is _ locally _ well - posed , while it is hadamard _ lipschitz ill - posed _ for @xmath83 . at the critical value @xmath84
the problem is globally well - posed for suitably small initial data ( and source term ) , but is ill - posed for sufficiently large initial data . in terms of the non - dimensional parameters introduced above in ,
this dichotomy can be seen as global well - posedness for @xmath85 , and lipschitz ill - posedness for @xmath86 .
a further feature of interest is that the anisotropy of the symbol can be explored to obtain an improvement in the regularity of the solutions when the initial data and source are _ supported on a plane in fourier space_. for such well - prepared data the local existence and uniqueness of solutions can be obtained for all values @xmath80 , and the global existence holds for all initial data when @xmath87 .
we note that some of the theorems discussed in this article are valid for certain classes of partial differential equations cf .
below , and in any spatial dimension @xmath88 . for such theorems ,
the three dimensional mg equation is a particular example of an active scalar equation for which these results hold .
we conclude the manuscript with section  [ sec : concl ] , in which we present a few remaining open problems regarding the mathematical analysis of the mg equations .
in this section we discuss the global regularity of solutions to  and the lyapunov stability of steady states , in the case that @xmath69 .      in order to show that there is no finite time loss of regularity of solutions to  with @xmath69 let us first perform a rough derivative count .
the dissipation is smoothing by two derivatives , while the nonlinearity is also losing two derivatives : the first coming from @xmath89 , while the second one derivative loss is coming from the constitutive law @xmath90 $ ] , in view of . in this sense
we may call the problem `` critical '' ( we make this more precise below ) .
since the operator @xmath37 relating @xmath38 and @xmath16 is an order one pseudo - differenential operator , one may view the mg equation as an example of a more general class of of active scalar equations : @xmath91 with @xmath92 a matrix of zero - order pseudo - differential operators , and we use the summation convention on repeated indices .
we have dropped the smooth source term @xmath93 from since it is not an obstacle to the well - posedness theory , have re - scaled the time variable so that the diffusivity coefficient is @xmath63 instead of @xmath29 ( this is not an issue since we will construct global in time solutions ) , and the @xmath92 corresponding to are given by @xmath94 .
we call this class of active scalar equation `` singular '' since the drift velocity is given in terms of the advected scalar by a constitutive law which is losing derivatives .
we refer the interested reader in to the work of chae , constantin , cordoba , gancedo , and wu  @xcite and also  @xcite for further discussions about the well - posedness theory of singular active scalar equations .
to make the idea of criticality more precise , we note that the active scalar equation  has a natural scaling invariance built into it : if @xmath95 is a solution , then @xmath96 is _ also _ a solution , with corresponding drift given by @xmath97 .
hence @xmath98 is the critical lebesgue space with respect to the natural scaling of the equation .
additionally is a scalar equation , so that @xmath99 norms are non - increasing in time for @xmath100 , and hence the @xmath101 norm is the `` strongest '' a priori controlled norm along solutions . in this sense
we say that  is a `` critical '' equation .
another active scalar equation derived from the navier - stokes equations for a rapidly rotating fluid in a thin shell is the surface quasi - geostrophic equation , introduced by constantin , majda , and tabak  @xcite ( see also  @xcite and references therein ) , which in the critically dissipative case reads @xmath102 and @xmath103 represents the @xmath104 riesz transform
. there are strong analogies between the systems  and . the velocity field is divergence free in both cases , and the @xmath101 norm is critical with respect to the natural scaling of the equation .
the main difference is that in for sqg the dissipative operator is non - local , but the drift velocity is not more singular than the advected scalar .    inspired by the de giorgi - type proof of global regularity for the critically dissipative sqg equation of cafarelli and vasseur  @xcite ( see also  @xcite for different proofs ) , to show that solutions of  do not blow up in finite time , we first look at the associated linear problem , where the drift velocity is _ given _ and obeys the bounds that are a priori available for weak solutions of the nonlinear problem : @xmath105 . in
@xcite we use a de giorgi iteration technique to prove that a weak solutions of the _ linear _ drift - diffusion equation @xmath106 instantly become hlder continuous .
more precisely , the result is :    [ thm : degiorgi ] let the initial data @xmath107 and @xmath108 be a weak solution of , where the drift velocity @xmath109 is given by for a suitable matrix @xmath110 .
then , for any @xmath111 , there exists @xmath112 such that @xmath113 .
the above result was obtained independently by sergin , silvestre , sverak , and zlatos in  @xcite . in view of the constitutive a direct corollary of therem  [ thm : degiorgi ]
is that :    weak solutions @xmath16 of  are hlder continuous in positive time , i.e. they obey @xmath114 for any @xmath115 , and a suitable @xmath112 .    in view of the above discussion about scaling
, hlder regularity is sub - critical information for - and one may bootstrap to prove that the solution is classical :    [ thm : higher ] let @xmath116 be a weak solution of the initial value problem associated to
, with @xmath117 .
then @xmath118 , for any @xmath119 , for some @xmath120 .
the proof of theorem  [ thm : higher ] relies on a suitable bootstrap procedure in besov spaces , in the spirit of  @xcite . repeating the argument for the derivative of @xmath16 one may further obtain that the solution of  is in fact @xmath121 smooth in finite time .
we conclude this subsection with a very brief discussion of the proof of theorem  [ thm : degiorgi ] .
we refer the reader to  @xcite for all details . along the lines of caffarelli and vasseur  (
* theorem 3 ) for the critical sqg equation , the first step consists of showing that a weak solution is bounded in the spacial variable ( see  ( * ? ? ? * appendix a ) for a construction of weak solutions which obey the energy inequality ) .
this is achieved by proving that the function @xmath122 has zero energy if @xmath123 is chosen large enough .
the key ingredient is the inequality @xmath124 which can be shown to hold for weak solutions . using that @xmath109 is divergence
free one obtains the energy inequality @xmath125 for all @xmath126 and @xmath127 . for @xmath128 , and @xmath129 to be chosen sufficiently large , letting @xmath130 , @xmath131 , and @xmath132 it follows from that @xmath133 for some constant @xmath134 .
letting @xmath135 be sufficiently large we thus obtain implies that @xmath136 exponentially as @xmath137 , and hence @xmath138 .
having established the boundedness of the weak solutions , the next step is to prove hlder regularity , which requires a more delicate argument .
the main obstruction to applying the classical parabolic de giorgi estimates via an @xmath99-based caccioppoli inequality ( @xmath139 ) , is that @xmath140 is the derivative of a @xmath141 function . by essentially using the divergence - free nature of @xmath109 , and appealing to the john - nirenberg inequality
, we may however prove a suitable local energy inequality : for any @xmath142 and @xmath143 , we have @xmath144 where @xmath145 is a fixed positive constant . here
we have denoted by @xmath146 the parabolic cylinder @xmath147\times b_\rho(x_0)$ ] for @xmath148 and @xmath149 .
moreover , estimate also holds with @xmath16 replaced by @xmath150 .
the bound may be in turn used , using a suitable iteration argument to obtain an improvement of oscillation result .
assume that @xmath151 , where @xmath152 is arbitrary .
then we have @xmath153 for some positive constant @xmath145 .
we note that the prefactor of @xmath154 of the right side of . besides estimate , the main ingredient in the proof of hlder continuity is the following estimate on the growth of level sets of the solution .
there exist @xmath155 , @xmath156 , and @xmath157 such that for @xmath158 if we know that @xmath159 $ ] we have @xmath160 , @xmath161 , @xmath162 , @xmath163 , and @xmath164 . to prove the above estimate one uses a suitable version of the local parabolic energy inequality .
bounds and are the main ingredients in the de giorgi iteration procedure , which at the end of the day yields the improvement in oscillation @xmath165 for some @xmath80 , where @xmath166 \times b_4 $ ] , @xmath167 \times b_r$ ] , @xmath158 are arbitrary , @xmath168 , and @xmath169 is as above . since @xmath170 , is then classical to show that implies the solution @xmath16 is hlder continuous , with exponent that may be computed explicitly .      in this subsection
we consider the stability of steady states to . first we prove for steady states of the the more general class of equations , with drift velocity given by , that linear implies nonlinear instability .
we then turn to the specific mg equations
, for which we construct an explicit unstable eigenvalue of the linear operator .
let @xmath171 be a steady solution of
, i.e. @xmath172 where @xmath173 is a time independent force , @xmath69 , and @xmath92 are zero order pseudo differential operators . writing the solution of  as @xmath174
we obtain that the perturbation temperature @xmath16 obeys @xmath175 where the _ dissipative linear _
operator @xmath17 is defined as @xmath176 with @xmath177 , and the _ nonlinear _
operator @xmath56 is given by @xmath178 recall the following notion of lyapunov stability :    [ def : nonlinear : instability ] let @xmath179 be a pair of banach spaces .
a steady state @xmath180 is called @xmath179 nonlinearly stable if for any @xmath148 , there exists @xmath181 so that if @xmath182 and @xmath183 , then    = 1em    there exists a global in time solution @xmath16 to the initial value problem  with @xmath184 ;    and we have the bound @xmath185 for a.e .  @xmath186 .
an equilibrium @xmath187 that is not stable in the above sense is called lyapunov unstable .    in the above definition ,
the banach space @xmath188 is the space where a local existence theorem for the nonlinear equations is available ( e.g.  @xmath189 ) , while @xmath190 is the space where the spectrum of the linear operator is analyzed , and where the instability is measured ( e.g.  @xmath191 ) .
the following result states that the linear lyapunov instability of a steady state @xmath180 implies the nonlinear one .
[ thm : linear : nonlinear ] let @xmath180 be a smooth , mean zero solution of . if the associated linear operator @xmath17 , as defined in , has spectrum in the unstable region , then the steady state is @xmath192 lyapunov nonlinearly unstable , for arbitrary @xmath193 .
the proof relies on a bootstrap argument that is by now well - known for fluids equations  @xcite , and an a priori time independent bound on the @xmath194 norm of the solution , which we obtain from the de giorgi @xmath195 estimate .
we only briefly sketch the argument given in full detail in @xcite , to which refer the interested reader .
first , we fix a linearly unstable smooth eigenfunction @xmath196 , with eigenvalue of _ maximal _ real part @xmath197 . in order to prove theorem  [ thm : linear : nonlinear ] , for @xmath198 $ ] we consider solutions @xmath199 of @xmath200 where @xmath201 and @xmath196 are as above , and
show that there exits a constant @xmath202 such that for any @xmath198 $ ] we may find a time @xmath203 with @xmath204 by definition  [ def : nonlinear : instability ] , the bound shows that the trivial solution @xmath205 of  is @xmath192 lyapunov unstable : indeed , an @xmath206 perturbation leads in finite time to an @xmath207 size solution , in finite time .
in order to prove we use the following classical idea : the initial data being an eigenfunction of @xmath17 , will grow exponentially in time , while the nonlinear term may be suitably controlled iso that it can not cancel this exponential growth . to make this idea concrete ,
it is convenient to define @xmath208 , where @xmath209 . due to the shift
, the resolvent of @xmath210 contains the full right - half of the complex plane , and being a bounded perturbation of the laplacian we also have that @xmath211 generates an analytic semigroup over @xmath212 , and @xmath213 is smoothing by two derivatives . then from the duhamel formula and the interpolation bound @xmath214 , we obtain the estimate @xmath215 the key step which allows us to close the estimates is proving the existence of a constant @xmath216 such that @xmath217 for all sufficiently small @xmath27 .
the proof of estimate relies on the @xmath195 bound which may be obtained from theorem  [ thm : degiorgi ] and a suitable bootstrap procedure to reach the @xmath218 regularity level ( see ( * ? ? ?
* appendix a ) ) .
lastly , let @xmath219 and @xmath220 be the maximal time such that @xmath221 from and one may then obtain @xmath222 and then in turn that @xmath223 a lower bound on @xmath224 , similar to but with the linear term minus the @xmath225 power term , may be obtained from the duhamel formula . in turn
, this implies that @xmath226 and concludes the proof of theorem  [ thm : linear : nonlinear ] .      in order to apply theorem  [ thm : linear : nonlinear ] to the ( nonlinear ) mg equations  , we need to find a steady state @xmath180 for which the corresponding linear operator @xmath17 defined in , has unstable spectrum .
to achieve this we use an idea that dates back at least to  @xcite for navier - stokes and  @xcite for euler .
first , note that @xmath227 , with @xmath228 smooth , is automatically a steady state , with associated velocity @xmath229 , if we assume that the source term is given by @xmath230 .
we thus consider the source term to equal @xmath231 , for which we have a steady state solution of given by @xmath232 for the above choice of the steady states , the equation @xmath233 becomes @xmath234 , \label{eq : eigenvalue}\end{aligned}\ ] ] where @xmath235 is the fourier multiplier defined earlier in . next , we make the _ ansatz _ that the temperature @xmath16 which solves has the form @xmath236 for some fixed @xmath50 in @xmath237 , to be chosen precisely later . using
we may compute from the third component of the velocity field @xmath38 associated with @xmath16 given by , insert this into and obtain that @xmath238 and @xmath239 must solve @xmath240 where we have denoted @xmath241 and @xmath242 with @xmath56 being the non - dimensional constant in .
an essential feature of the coefficients @xmath243 is that they form an increasing sequence which grows rapidly , as @xmath244 when @xmath245 .
equation gives a recurrence relation for the sequence @xmath239 ( see e.g. @xcite ) : @xmath246 we emphasize that given any @xmath247 ( which then uniquely defines @xmath248 for all @xmath249 ) and given any @xmath250 , one can solve for all @xmath251 using the recursion relations . the key is however that only for suitable values of @xmath238 do the @xmath252 vanish sufficiently fast so that @xmath16 defined by is sufficiently smooth .
solving for such @xmath238 amounts to solving the infinite continued fraction equation @xmath253 solving is the main difficulty in this argument , since we need not just show there is a solution , but also give bounds on it in terms of the @xmath243 .
we omit these details and summarize the results obtained in ( * ? ? ?
* lemma 2.6 ) , ( * ? ? ?
* lemma 4.2 ) , and ( * ? ? ? * theorem 2.2 ) as follows .
[ lem : contd : frac ] let @xmath254 be defined by and respectively , and let @xmath255 .
there exits a real eigenvalue @xmath256 which solves , and we have the estimate @xmath257 moreover , the corresponding sequence @xmath239 decays exponentially in @xmath258 .
a direct consequence of estimate is that the eigenvalue @xmath259 obeys the upper and lower bounds @xmath260 for @xmath261 fixed , we hence obtain from that there exists at most finitely many frequency pairs @xmath262 for which @xmath263 .
moreover , one may show that in the parameter range of the geodynamo , i.e. @xmath264 , for suitable @xmath265 the left side of is positive ( see below ) , thereby proving :    associated to the steady state @xmath180 defined in , the linearized mg operator @xmath17 defined in  , has unstable spectrum .
lastly , note that one may optimize in @xmath266 , and compute a lower bound for the largest eigenvalue of @xmath17 in terms of the physical parameters @xmath267 , and of @xmath268 .
more precisely , we have that @xmath269 with the lower bound being attained e.g. for ( of course , we need to take integer part ) @xmath270 in particular , we note that as @xmath271 the lower bound given by the right hand side of blows up ( this will be important in sections  [ sec : inviscid ] and  [ sec : frac ] ) .
moreover , for the above choice of @xmath272 and @xmath273 , we have @xmath274 as @xmath271 .
this is in precise accordance to the fact that the absolute value of the fourier multiplier associated to @xmath37 grows at the fastest rate when @xmath275 as @xmath61 .
see .    combining the results in this subsection with theorem  [ thm : linear : nonlinear ] , we have thus constructed an example of solutions to the mg equations which displays strong ( exponential ) growth in @xmath16 and hence in the velocity field @xmath38 by and also in the magnetic field @xmath276 through .
the rate of this growth can be shown to be proportional to @xmath277 which is consistent with the dynamo instability scenario .
as discussed in the introduction , in the fluid core the value of @xmath29 , given by the quotient of @xmath21 the molecular diffusivity coefficient , and @xmath278 , the product of the typical length scale and typical velocity in the fluid , is very small .
therefore , from the physical point of view it is natural to study the non - diffusive ( unforced ) mg equation @xmath279 , \quad \nabla \cdot u = 0 \label{eq : inv:2}\end{aligned}\ ] ] or the corresponding non - diffusive singular active scalar equation when is replaced by , i.e. @xmath280 . as we shall see below
, the system  is _ ill - posed _ in sobolev spaces , in the sense that there is no operator which is lipschitz in time at @xmath78 .      before we discuss the issue of ill - posedness ,
let us briefly discuss why the _ evenness _ of the operator @xmath37 breaks the classical proof of local existence in sobolev spaces for the class of singular active scalar equation . the standard energy - approach to local existence in @xmath281 would entail applying the @xmath282 to , and then take @xmath212 inner products with @xmath283 . morally speaking ( this is rigorous when @xmath284 is an integer in view of the leibniz rule ) we have only one term which seem to prevent closing the estimates at the @xmath281 level : @xmath285 indeed , the other term where more than @xmath286 derivatives lie on a single function is @xmath287 , but this term vanishes upon integrating by parts since @xmath288 . the only hope to treat the term @xmath289 would be to discover a commutator structure .
however , since @xmath37 is _ not anti - symmetric _ ,
i.e. _ even _ in fourier space , we _ can not _
write @xmath290 , where @xmath291 ( -\delta)^{s/2 } \theta \ ; ( -\delta)^{s/2 } \theta dx = { \mathcal t }   + { \mathcal c } \label{eq : s : good}\end{aligned}\ ] ] and @xmath292 $ ] is the usual commutator bracket .
instead for the mg , we have @xmath293 , and we _ can not _ explore the fact that @xmath294 may be bounded in terms of @xmath295 ( if @xmath296 ) .
this is the main reason why we do not seem to be able to close estimates at the sobolev level .    on the other hand ,
it was observed in  @xcite , see also @xcite , that if we replace @xmath90 $ ] ( cf .  ) by a certain class of singular constitutive laws @xmath297 ( cf .  ) , and the matrix @xmath92 is such that @xmath298 for any @xmath299 and @xmath300 ( here @xmath301 is the usual @xmath212 inner product ) , then the corresponding singular active scalar equation  is locally well posed in @xmath281 for @xmath296 .
the main reason is that under these assumptions the operator @xmath302 is anti - symmetric , and we are able to write @xmath303 , using the notation of and .
the commutator term @xmath294 would enable one estimates to close at the @xmath281 level .
lastly , we mention that for analytic initial data , the non - diffusive mg equation is indeed locally well - posed in the class of real - analytic functions , since each term in the equation loses at most one derivative .
this was proven in ( * ? ? ?
* theorem  2.1 ) by using a cauchy - kowalewski - type argument , in the spirit of the gevrey - class regularity proof of  @xcite ( see also  @xcite ) .
in particular , this result justifies the short time numerical simulations of the mg equation are justified , as they only keep track of finitely many fourier modes ( hence the associated functions are real analytic ) .
a cauchy problem for a certain partial differential equation is called _ well - posed in @xmath188 in the sense of hadamard _ , if for any initial data in @xmath188 , the problem has a unique solution in @xmath304 , with @xmath305 depending only on the initial data in the @xmath188-norm , and the solution map @xmath306 satisfies some continuity properties , for a suitable space @xmath307 .
if one of these properties fail , the cauchy problem is called _ ill - posed_. in this section we discuss the failure of the solution map for the non - diffusive mg equation to be _ lipschitz continuous _ with respect to perturbations in the initial data ( in the topology of a sobolev space @xmath188 ) around a specific steady profile @xmath187 .
such results have been recently obtained for several models in fluid dynamics , see e.g.  @xcite and references therein .
the proof consists of a linear and a nonlinear step .
we construct a steady state and a sequence of eigenfunctions , with arbitrarily large eigenvalues , for the linearized mg equation around this steady state .
once these eigenvalues are exhibited , one may use a fairly robust argument ( see , e.g.  @xcite ) to show that this severe linear ill - posedness implies the lipschitz ill - posedness for the nonlinear problem .
we emphasize that the key feature of the mg constitutive law which gives this ill - posednes is the evenness and unboundedness of the symbol associated to @xmath37 .
in fact , the evenness also gives non - uniqueness for for bounded weak solutions of the non - diffusive mg equation , as was recently established by shvydkoy  @xcite .
let @xmath308 be a positive integer .
we use the steady state @xmath309 with associated velocity field @xmath310 . similarly
to , we consider the associated linearized ( non - diffusive ) mg operator defined by @xmath311 \partial_3 \th_0 \label{eq : steady : l}\end{aligned}\ ] ] where the @xmath235 operator is defined in terms of its fourier symbol . with this notation , we have    [ thm : ill : lin ] the cauchy problem associated to the linearized @xmath312 equation @xmath313 is ill - posed in the sense of hadamard over @xmath212 .
that is , for any @xmath314 , there exits a real - analytic initial datum @xmath315 , such that the ivp associated to does not have a solution @xmath316 , which obeys @xmath317 } \| \theta(t,\cdot ) \|_{l^2 } \leq k \|\theta(0,\cdot ) \|_{y } \label{eq : linear : well}\end{aligned}\ ] ] where @xmath318 is any sobolev space that embeds compactly in @xmath212 .
the proof of the above theorem builds on the construction of unstable eigenvalues in the critically diffusive case , more specifically on lemma  [ lem : contd : frac ] .    as in section  [ sec : unstable :
eigenvalue : dissipative ] , cf .
we make the ansatz that @xmath16 is given in terms of a fourier sine series @xmath319 with @xmath272 and @xmath273 to be chosen precisely later .
note that since we are in the non - diffusive case , we need to set @xmath76 , and hence cf .
we have @xmath320 for any @xmath321 .
the value of @xmath243 remain however unchanged . using lemma  [ lem : contd : frac ]
, we obtain an eigenvalue @xmath322 , and an exponentially decaying sequence @xmath239 , such that cf .   the bound @xmath323 holds .
we set @xmath324 and @xmath325 in , where @xmath326 is a positive integer , and obtain the lower bound @xmath327 on the eigenvalue , associated to this choice of @xmath265 . here
we also denote by @xmath329 the eigenfunction constructed in this way . using an additional argument , which requires suitably bounding the exponentially decaying sequence @xmath239 which defines our eigenfunction , one may in fact additionally normalize @xmath16 and prove ( cf .
* lemma 2.7 ) , ( * ? ? ?
* lemma 4.3 ) ) that @xmath330 where @xmath331 is a sufficiently large constant .
for example , if @xmath332 with @xmath333 , we may explicitly compute a suitable such constant .
we conclude the proof of theorem  [ thm : ill : lin ] by using a contradiction argument .
if such a solution would exist , for a fixed @xmath305 and @xmath334 , we pick @xmath335 large enough so that @xmath336 holds .
the linear problem   can be shown to have a unique @xmath212 solution , and hence by construction we have @xmath337 the assumed estimate is contradicted by our choice of @xmath338 and the lower bound .
this concludes the proof of theorem  [ thm : ill : lin ] .      before stating our main ill - posedness result ,
let us recall the definition of lipschitz well - posedness .
[ def : well ] let @xmath339 be sobolev spaces .
the cauchy problem for the non - diffusive mg equation is locally lipschitz @xmath340 _ well - posed _
, if there exist continuous functions @xmath341 , so that for every pair of initial data @xmath342 there exist unique solutions @xmath343 of the initial value problem associated to  , that additionally satisfy @xmath344 for every @xmath345 $ ] , where @xmath346 and @xmath347 .
we note that defintion  [ def : well ] allows the solution map to lose regularity .
since the mg equation is first order , the typical pairs of spaces @xmath340 that we have in mind here are @xmath348 , and @xmath349 , with @xmath350 .
the main result of this section is then :    [ thm : nonlinear : ill ] the non - diffusive mg equations  , are locally lipschitz @xmath340 _ ill - posed _ in sobolev spaces @xmath351 embedded in @xmath352 , in the sense of definition  [ def : well ] .
the proof of the above result is based on the severe linear ill - posedness obtained in theroem  [ thm : ill : lin ] , and a standard perturbative argument  ( * ? ? ? *
theorem 2 ) or  @xcite , which we breifly sketch below .    let @xmath353 be the steady state defined above in , let @xmath354 \partial_r \th_0 $ ] be its associated linearized mg operator , and denote the mg nonlinearity by @xmath355 \partial_j \theta$ ] . by our assumptions on @xmath188 and @xmath318 , we have @xmath356 for some constant @xmath357 .
also , let @xmath358 .    in definition , we chose @xmath359 be the steady state , and let @xmath360 , for some smooth function @xmath361 to be chosen later .
assume by contradiction that the problem is lipschitz locally well - posed in @xmath340 .
one may then show that there exist @xmath362 and @xmath363 such that @xmath364 } \|\theta^{(1)}_\delta(t ) - \th_0\|_{x } \leq k \delta\ ] ] for _ any _ @xmath365 $ ] .
writing the solution @xmath366 as an order @xmath367 perturbation around the steady state , @xmath368 we thus have @xmath369 } \|\psi_\delta(t)\|_x \leq k \label{eq : psi : bound:1}\end{aligned}\ ] ] and @xmath370 solves the cauchy problem @xmath371 the key observation is that when combined with , estimate yields @xmath372 for any @xmath373 $ ] .
therefore , when passing @xmath374 , the nonlinearity in disappears , and we should obtain the linear limiting equation . indeed , from , , and the bounds , we obtain the necessary estimates to apply the aubin - lions compactness argument and show that @xmath375 strongly in @xmath376 and on this time interval the limiting function obeys @xmath377 but @xmath378 also inherits the bound from the sequence @xmath370 , which contradicts thereom  [ thm : ill : lin ]
choose @xmath361 an eigenfunction of @xmath17 as constructed in section  [ sec : ill : linear ] , so that in time @xmath379 the exponential growth of the @xmath212 norm caused the solution to exceed the level @xmath380 .
this concludes our sketch of the proof of theorem  [ thm : nonlinear : ill ] .
the fractionally diffusive mg equation , although non - physical , presents a very challenging and interesting mathematical problem . replacing the laplacian in ( [ eq : mg:1 ] ) by @xmath381 , for @xmath80
, we obtain the active scalar equation @xmath382 \label{eq : mgf:2},\\      & \theta(0,\cdot ) = \theta_0 \label{eq : mgf:3}\end{aligned}\ ] ] where the symbol of the fourier multiplier @xmath37 is given by ( [ eq : m1])([eq : d ] ) .
note that here we work with the non - dimensional variables , so that the initial data and the source term are _ normalized to be of size unity_.    the main feature of this situation is that the structure of @xmath383 ( anisotropy , evenness , and unboundedness ) and a relatively weak smoothing effect of the nonlocal operator @xmath384 produce a sharp dichotomy across the value @xmath84 , as follows . if @xmath82 the equation is _ locally _ well - posed , while it is hadamard _ lipschitz ill - posed _ for @xmath83 . at the threshold value @xmath84
the problem is globally well - posed for @xmath85 , respectively ill - posed for @xmath385 .
furthermore , the anisotropy of the symbols @xmath383 can be used to improve the regularity of the solutions when the initial data and source have a suitably chosen fourier support . for such _ well - prepared _ data the local existence and uniqueness of solutions can be obtained for all values @xmath80 , and the global existence holds for all initial data when @xmath87 .
we remark that in all above mentioned regimes , the problem is _ super - critical _ in the sense described in section [ sec : diffusive ] .
in section [ sec : diffusive ] we stated the well - posedness result when the diffusion in the equation is given by the full laplacian .
if the diffusion in the equation is introduced by the fractional operator @xmath384 , @xmath87 , then we have the following result .
[ thm : fractional_local ] let @xmath87 , and fix @xmath386 .
assume that @xmath387 and @xmath388 have zero - mean on @xmath35 .
then there exists a time @xmath389 and a unique smooth solution @xmath390 of the cauchy problem ( [ eq : mgf:1])([eq : mgf:2 ] ) .    under an additional assumption on the size of @xmath29 ( equivalently , the size of the initial data in the unscaled equations ) ,
we obtain a global existence result .
[ thm : fractional_global ] let @xmath391 and @xmath93 be as in the statement of theorem  [ thm : fractional_local ] , and let @xmath392 have zero - mean on @xmath35 , where @xmath393 .
there exists a large enough number @xmath394 , such that if @xmath395 , then the unique smooth solution @xmath16 of the cauchy problem  is global in time , i.e. @xmath396 .    [ rem : rescale ] note that in  @xcite we did not use the non - dimensional variables with which we are working here .
therefore , the global well - posedness condition stated in  @xcite is that the initial data and source , have to be _ sufficiently small _ , compared to the diffusivity coefficient @xmath21 . in this paper
however , the data and source have been normalized to be of order unity , and the only coefficient remaining is @xmath29 .
it is hence natural that the small data condition turns to a _
large enough _ @xmath29 condition .
this remains valid throughout this section .
the detailed proofs of theorems  [ thm : fractional_local ] and  [ thm : fractional_global ] are given in @xcite , and rely mainly on energy methods and the use of suitable commutator estimates .
we omit these details here .      in this subsection
we consider the stability of steady states to . the construction follows the outline described in section  [ sec : unstable : eigenvalue : dissipative ] . in this case
, we choose @xmath397 and @xmath398 , for some integer @xmath399 .
as before , the first step in proving that the cauchy problem  is lipschitz ill - posed is to consider the linearized problem around the given steady state , and exhibit the existence of unstable eigenvalues with arbitrarily large real part .
thus , we consider @xmath400 where the linear operator @xmath17 is given by @xmath401 we seek the temperature @xmath16 which solves ( [ eq : fractional_inst1 ] ) and has the form @xmath402 for some fixed @xmath403 to be chosen later . inserting ( [ eq : fractional_theta ] ) into ( [ eq : fractional_inst1 ] ) and using the definition of @xmath235 , we obtain that holds , but with @xmath404 denoted by @xmath405 and with @xmath243 given by above . as before ,
these assumptions lead to solving the continued fraction equation for some real @xmath406 .
this is achieved by suitably modifying lemma  [ lem : contd : frac ] to obtain :    [ lemma : sn:5.2 ] let @xmath407 be defined by and respectively , and let @xmath408 .
there exists a real eigenvalue @xmath409 , and a rapidly decaying sequence @xmath410 so that the function @xmath16 defined by solves
. furthermore , we have the lower bound @xmath411    in particular , setting for instance @xmath412 and @xmath413 in , for each @xmath326 we obtain a real eigenvalue @xmath414 of @xmath17 , which is bounded from below as @xmath415 moreover , given a sobolev space @xmath318 that is compactly embedded in @xmath212 , we may ensure the existence of a constant @xmath416 such that the smooth eigenfunction @xmath329 corresponding to @xmath414 obeys .
the immediate consequence of estimate is that due to @xmath417 , we have @xmath418 the above construction proves that the linearized problem is ill - posed in @xmath212 , as we have exhibited eigenvalues of arbitrarily large real part . following the argument outlined in section  [ sec : ill : nonlinear ] , this type of strong ill - posedness for the linearized problem , automatically implies that the nonlinear problem is lipschitz ill - posed in sobolev spaces ( recall defintion  [ def : well ] of lipschitz well - posedness ) .
thus we obtain :    [ thm : nonlinear : ill : frac ] the cauchy problem  , with @xmath419 is locally lipschitz @xmath340 _ ill - posed _ in sobolev spaces @xmath351 embedded in @xmath420 .
the value @xmath84 marks the threshold for proving the local well - posedness in sobolev spaces for . if the initial datum is small with respect to @xmath21 , or equivalently , if @xmath29 is sufficiently large ( see remark  [ rem : rescale ] above ) , then one may use energy estimates to show that there exists a unique global smooth solution evolving from this data . in this direction
we have :    [ thm : gamma=1/2:small ] let @xmath421 and assume that @xmath387 and @xmath388 have zero - mean on @xmath422 .
there exists a sufficiently large number @xmath394 , such that if @xmath395 , then the unique smooth solution @xmath423 of the cauchy problem  is global in time .    as in the case of @xmath87 , the proof is based on the classical kato - ponce commutator estimate .
the proof however does not apply for the case of `` small '' @xmath29 ( large data ) , not even to prove the local existence of solutions .    in the case where @xmath29 is suitably small
, we may construct a steady state such that solutions evolving from initial data which is arbitrarily close to this steady state diverge from it at an arbitrarily large rate , i.e. the problem is hadamard ill - posed .
more precisely , we have :    [ thm : ill:1/2 ] let @xmath424 be sobolev spaces , let @xmath81 , and fix the integer @xmath425 . consider the steady state solution of given by @xmath426 , with source term @xmath427 , and assume that @xmath428 then , given any @xmath389 and any @xmath429 , there exists an initial condition @xmath430 and a corresponding solution @xmath431 of the cauchy problem  , such that @xmath432 } \vert   \theta^{(2)}(t,\cdot ) - \theta^{(1)}(t,\cdot ) \vert_{x } \geq 2 k \vert   \theta^{(2)}_0 - \theta^{(1)}_0 \vert_{y}\end{aligned}\ ] ] where @xmath433 .
that is , the problem  is locally lipschitz ( x , y ) ill - posed when @xmath86 .    to prove the above result , we follow the ill - posedness result in the case @xmath419 .
the main difference is that while for @xmath419 by we have real eigenvalues of arbitrarily large real part , when @xmath434 the lower bound only blows up for @xmath29 sufficiently small , depending on @xmath24 and our choice of @xmath268 .
this condition is rather unsatisfactory , since it implies @xmath271 as @xmath435 .
in fact , we can establish a more robust smallness condition on @xmath29 , in the sense that it is independent of @xmath24 ( and @xmath268 ) . to see this , insert @xmath436 and @xmath437 in the lower bound , where @xmath438 .
we obtain the lower bound @xmath439 therefore , under the assumption , which is independent of @xmath268 or @xmath56 , the right side of diverges to @xmath440 as @xmath441 and hence we again have a sequence of unstable eigenvalues of the linear problem , with arbitrarily large real part .
combined with the nonlinear argument presented in section  [ sec : ill : nonlinear ] this concludes the sketch of the proof of theorem  [ thm : ill:1/2 ] .      in this subsection
we present results which take advantage of the curved nature ( paraboloid - type ) of the regions in frequency space where the symbol @xmath442 grows like @xmath443 .
we observe that if the support of @xmath444 lies on a plane , then the operator @xmath37 behaves like an operator of order zero acting on this @xmath16 , and hence the velocity @xmath38 is as smooth as @xmath16 .
this allows us to improve the well - posedness results in the range @xmath69 and @xmath445 , over the generic case when no conditions on the fourier transforms of the initial data ( and source term ) are imposed .    to be more precise
let us introduce some notation .
for a rational number @xmath446 with @xmath447 , define frequency plane @xmath448 as @xmath449 one may show that whenever @xmath447 , latexmath:[\[\begin{aligned }   \label{eq : m : plane : bound }    this observation , combined with the fact that when the data @xmath15 and source @xmath93 have frequency support on @xmath451 , then so does the solution @xmath452 of the cauchy problem  at all later times @xmath453 , means that the smoothing properties of the fractional diffusion term are stronger than they are for the generic data situation .
hence it is possible to prove stronger regularity results , and in particular the _ local existence _ and uniqueness of sobolev solutions holds _ for all _ values of @xmath80 .
[ thm : well - prepared : local ] let @xmath80 , and fix @xmath454 .
assume that @xmath387 and @xmath388 have zero - mean on @xmath422 and are such that @xmath455 for some fixed @xmath456 and all @xmath457 .
then there exists a @xmath389 and a unique smooth solution @xmath458 of the cauchy problem  , such that @xmath459 for all @xmath460 .
we note that the above theorem holds even if @xmath419 , a setting in which we know that for _ generic _ initial data the problem is ill - posed .
moreover , one may also prove that if the initial data is small with respect to @xmath21 , or equivalently , @xmath85 in the non - dimensional setting of this manuscript , the solution obtained in theorem  [ thm : well - prepared : local ] is global in time .    the main difficulty in the proof of theorem  [ thm : well - prepared : local ] is the construction and uniqueness of solutions , which live on the frequency plane @xmath451 .
this relies on constructing of a suitable sequence of approximations supported on the fourier space on @xmath451 and showing its convergence using energy methods .
the main obstacle one has to overcome is designing an iteration scheme which is both suitable for energy estimates , and is such that in each iteration step the frequency support of the approximation lies on @xmath451 .
more explicitly , a scheme in which @xmath461 is given in terms of @xmath462 preserves the frequency support on @xmath451 , however is not good for closing the estimates at the level of sobolev spaces . on the other hand
if we define @xmath461 in terms of @xmath463 , the energy estimates are now clear but it seems difficult to inductively prove that @xmath464 is supported on @xmath451 .
therefore , we work with a regularized problem , where we introduce an artificial term @xmath465 .
the advantage of considering a stronger dissipative problem is that we can _ construct _ smooth solutions which have frequency support lying on @xmath451 .
for details of the construction we refer the interested reader to @xcite .
lastly , we pass to the limit @xmath466 and obtain a solution of  with desired properties .    if the initial data @xmath15 and the source term @xmath93 have fourier support in @xmath451 , we show that for @xmath467 the local in time solution may be continued for all time .
[ thm : well - prepared : global ] let @xmath468 and @xmath93 be as in the statement of theorem  [ thm : well - prepared : local ] . if @xmath469 , the unique smooth solution @xmath16 of the cauchy problem  is global in time .
the proof hinges on the observation that the solution constructed in theorem [ thm : well - prepared : local ] satisfies @xmath470 .
then , in view of and the hrmander - mikhlin theorem we have @xmath471 for all @xmath472 , and so the velocity @xmath38 is as smooth as the scalar @xmath16 .
hence , the case @xmath473 is `` sub - critical '' for the a priori controlled quantities , and energy estimates in the spirit of those used to prove the global regularity of the sub - critical sqg equation  @xcite , may be used to prove theorem [ thm : well - prepared : global ] .
the case @xmath84 in theorem  [ thm : well - prepared : global ] is more subtle , and is analogous with the problem of global regularity for the critical sqg equation  @xcite . the anisotropic nature of the symbol @xmath37 , makes it more convenient to apply the de giorgi - inspired proof of @xcite in the context of the fractional mg equation with @xmath81 .
the main ingredient is that if @xmath474 and @xmath475 , then @xmath476 .
hence , considering first the linear problem with divergence - free @xmath141 drift we obtain hlder continuity of the solution .
the later norm is sub - critical and one may bootstrap to higher regularity .
we refer the interested reader to  @xcite for details .
we have proven that the magneto - geostrophic model proposed by moffatt and loper is capable of sustaining dynamo action .
the mg - dynamo instability exhibited here relies crucially on the fact that the dimensionless parameter @xmath24   is neither @xmath477 , nor @xmath440 .
the exponential growth is of order @xmath478 , for all @xmath479 , and @xmath480 . in terms of the physics in the model ,
this condition on @xmath24 requires that _ all three forces _ , coriolis , lorentz , and gravity , _ are non - zero _ , making the mg - instability somewhat different in nature from the usual rayleigh - taylor instability .    at the level of mathematical analysis several problems
remain to be addressed , including some well / ill - posedness issues . for the non - diffusive equation ,
_ is there a `` stronger '' form of ill - posedness such as norm inflation or even finite time blowup ? _ can we refine the well - posedness results for the fractionally diffusive problem in the range @xmath481 ?
that is , _ can one establish the global existence of smooth solutions or on the contrary , prove finite time blowup ? _    several , interesting , mathematical issues arise when returning to the derivation of the mg equation from the full three dimensional mhd equations in the spirit of  @xcite .
in particular , we note that moffatt s and loper s initial discussion includes the kinematic viscosity @xmath19 , which is directly proportional to @xmath27 in .
the values of @xmath27 and @xmath29 are highly conjectural , but both are believed to be _ extremely _ small .
therefore , as _ opposed _ to the mg system , if we keep @xmath482 we derive from  the system @xmath483 the presence of the laplacian term in makes the constitutive law relating @xmath38 and @xmath16 `` better '' , in the sense that the @xmath37 operator corresponding to  , now depends both on @xmath56 and @xmath27 , and is in fact smoothing of degree @xmath484 ( see equations ( 29)(30 ) in @xcite ) .
hence , in this case the constitutive law is even better than that for the two - dimensional euler equation in vorticity form , implying that : _ in the range @xmath482 , the system  is globally well - posed , even if @xmath485_.    this observation is interesting in view of the much more challenging problem of _ establishing properties of the solutions to  in the vanishing parameters limits @xmath486 and @xmath271 _ , as proposed by moffatt and loper . the later question is directly related to the turbulent scaling regime for @xmath16 conjectured by moffatt in  @xcite
. `` turbulence '' in the framework of vanishing diffusivity @xmath29 for the forced system  , and the presence of anomalous dissipation of energy for @xmath16 , or `` thetargy '' , may be plausible for very weak solutions ( below onsager critical ) of the non - diffusive equation .
p.  constantin , d.  cordoba , and j.  wu . on the critical dissipative quasi - geostrophic equation .
, 50(special issue):97107 , 2001 .
dedicated to professors ciprian foias and roger temam ( bloomington , in , 2000 ) .
magnetostrophic turbulence and the geodynamo . in y.
kaneda , editor , _ iutam symposium on computational physics and new perspectives in turbulence , nagoya , japan , september , 1114 , 2006 _ , volume  4 of _ iutam bookser .
_ , pages 339346 .
springer , dordrecht , 2008 .