Daniele Del Santo
email: delsanto@units.it

Francesco Fanelli
email: fanelli@math.univ-lyon1.fr

Gabriele Sbaiz
email: gabriele.sbaiz@phd.units.it

Aneta Wróblewska-Kamińska

On the influence of gravity in the dynamics of geophysical flows

Keywords: 2020 Mathematics Subject Classification: 35Q86 (primary), 35B25, 76U60, 35B40, 76M45 (secondary) barotropic Navier-Stokes-Coriolis system, gravity, stratification effects, low Mach, Froude and Rossby numbers, multiscale limit

In the present paper, we study a multiscale limit for the barotropic Navier-Stokes system with Coriolis and gravitational forces, for vanishing values of the Mach, Rossby and Froude numbers (Ma, Ro and Fr, respectively). The focus here is on the effects of gravity: albeit remaining in a low stratification regime Ma/Fr → 0, we consider scaling for the Froude number which go beyond the "critical" value Fr = √ Ma. The rigorous derivation of suitable limiting systems for the various choices of the scaling is shown by means of a compensated compactness argument. Exploiting the precise structure of the gravitational force is the key to get the convergence.

Introduction

In this paper we continue the investigation we began in [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF], about multiscale analysis of mathematical models for geophysical flows. Our focus here is on the effect of gravity in regimes of low stratification, but which go beyond a choice of the scaling that, in light of previous results, we call "critical".

1

In order to explain better all this, let us introduce some physics about the problem we are interested in, and give an overview of related studies. We present in Section 2 the precise system we will work on, and the statements of our main results.

Some physical considerations

By definition (see e.g. [START_REF] Cushman-Roisin | Introduction to Geophysical Fluid Dynamics[END_REF], [START_REF] Pedlosky | Geophysical fluid dynamics[END_REF]), geophysical flows are flows whose dynamics is characterised by large time and space scales. Typical examples are currents in the atmosphere and the ocean, but of course there are many other cases where such fluids occur out of the Earth, like flows on stars or other celestial bodies.

At those scales, the effects of the rotation of the ambient space (which in the previous examples is the Earth) are no more negligible, and the fluid motion undergoes the action of a strong Coriolis force. A simplistic assumption, which is however often adopted in physical and mathematical studies, consists in restricting the attention to flows at mid-latitudes, i.e. flows which take place far enough from the poles and the equatorial zone. Thus, if we denote by ≥ 0 the density of the fluid and by u ∈ R 3 its velocity field, the Coriolis force may be represented in the following form:

C( , u) := 1 Ro e 3 × u , (1.1)
where e 3 = (0, 0, 1), the symbol × denotes the classical external product of vectors in R 3 and Ro > 0 is the so-called Rossby number, a physical adimensional parameter linked to the speed of rotation of the Earth. In particular, the previous definition implies that the rotation is approximated to take place around the vertical axis, and its strength does not depend on the latitude. We point out that, despite all these simplifications, the obtained model is already able to give a quite accurate description of several physically relevant phenomena occurring in the dynamics of geophysical fluids (see e.g. [START_REF] Pedlosky | Geophysical fluid dynamics[END_REF], [START_REF] Chemin | Mathematical geophysics. An introduction to rotating fluids and the Navier-Stokes equations[END_REF]).

In geophysical fluid dynamics, effects of the fast rotation are predominant; this translates into the fact that the value of Ro is very small. As a matter of fact, the Rossby number Ro is defined as the ratio between the nonlinear acceleration to the Coriolis parameter term, namely

Ro := U ref f ref L ref ,
where U ref , L ref and f ref are respectively the horizontal velocity scale, the horizontal length scale and the reference Coriolis frequency (see e.g. [START_REF] Kundu | Fluid mechanics. Fifth edition[END_REF] for more details). For instance, for a typical atmospheric value of U ref ∼ 10 m/s, f ref ∼ 10 -4 s -1 and L ref ∼ 1000 km, the Rossby number turns out to be 0.1; its value is even smaller for many flows in the oceans. As established by the Taylor-Proudman theorem in geophysics, the fast rotation imposes a certain rigidity/stability, as it undresses the motion of any vertical variation, and forces it to take place on planes orthogonal to the rotation axis. Thus, the dynamics becomes purely two-dimensional and horizontal, and the fluid tends to move along vertical columns. However, such an ideal configuration is hinder by another fundamental force acting at geophysical scales, the gravity, which works to restore vertical stratification of the density. The gravitational force may be represented by the term G( ) := -1 Fr 2 e 3 , where Fr > 0 is the Froude number, another physical adimensional parameter, which measures the importance of the stratification effects in the dynamics. In geophysics (see again [START_REF] Kundu | Fluid mechanics. Fifth edition[END_REF] for details), Fr represents the square root of the ratio between inertia and gravity, namely

Fr := U ref √ g L ref ,
where g is the acceleration of gravity.

As it happens for the Rossby number, at large scales also the Froude parameter is typically very small. Thus, the competition between the stabilisation effect of the Coriolis force and the vertical stratification due to gravity, is translated in the model into the competition between the orders of magnitude of the two parameters Ro and Fr.

Actually, it turns out that the gravity G acts in combination with pressure forces. Restricting from now on our attention to the case of compressible fluids, like currents in the atmosphere for instance, and neglecting for a while heat transfer processes, the pressure term arising in the mathematical model takes the form P( ) := 1 Ma 2 ∇p( ) , where p is a known smooth function of the density (and, in the general case, of the temperature of the fluid) and Ma > 0 is the so-called Mach number, a third fundamental adimensional parameter which sets the size of isentropic departures from incompressible flow: the more Ma is small, the more compressibility effects are low. Also the value of Ma is very small for geophysical flows, since it is defined as

Ma := U ref c
with c being the sound speed (for instance, in the oceans the typical sound speed is c ∼ 1520 m/s).

As it is customary in physical studies, because of the complexity of the model, one would like to derive reduced models for geophysical flows, which however are able to retain most of the properties of the original system. The problem is that the three terms C, G and P enter into play in the model with a very large prefactor in front of them, owing to the smallness of the values of Ro, Fr and Ma respectively. For assessing their relative importance and their influence in the dynamics, one fixes a choice of their orders of magnitude. Actually (see e.g. the discussion in Section 1.4 of [START_REF] Cushman-Roisin | Introduction to Geophysical Fluid Dynamics[END_REF]), there is some arbitrariness in doing so, depending on the specific properties of the physical system and on the processes one wants to put the accent on.

In general, geophysical fluid dynamics is a multiscale process, meaning that Earth's rotation, gravity and pressure forces act, and interact, at different scales in the system. In other words, Ro, Fr and Ma have different orders of magnitude, and only for some specific choices, all (or some) of them are in balance.

Multiscale analysis: an overview of previous results

At the mathematical level, in the last 30 years there has been a huge amount of works devoted to the rigorous justification, in various functional frameworks, of the reduced models considered in geophysics. Studies have been carried out in various contexts: for instance, focusing only on the effects of the low Mach number, or on its interplay with a low Froude number regime.

Reviewing the whole literature about this subject goes far beyond the scopes of this introduction, therefore we make the choice to report only on works which deal with the presence of the Coriolis force (1.1). We also decide to leave aside from the discussion the case of incompressible models, because (owing to the rigidity imposed by the divergence-free constraint on the velocity field of the fluid) less pertinent for multiscale analysis. We refer to book [START_REF] Chemin | Mathematical geophysics. An introduction to rotating fluids and the Navier-Stokes equations[END_REF] and the references therein for a panorama of known results for incompressible homogeneous fluids, even though more recent developments have been made (see e.g. [START_REF] Scrobogna | Highly rotating fluids with vertical stratification for periodic data and vanishing vertical viscosity[END_REF] for a case where stratification is considered). Notice that there are also a few recent works [START_REF] Fanelli | Asymptotics of fast rotating density-dependent fluids in two space dimensions[END_REF], [START_REF] Cobb | On the fast rotation asymptotics of a non-homogeneous incompressible MHD system[END_REF], [START_REF] Sbaiz | Fast rotation limit for the 2-D non-homogeneous incompressible Euler equations[END_REF], dealing with incompressible non-homogeneous fluids, but results in that direction are only partial and the general picture still remains poorly understood at present.

The framework of compressible fluid models, instead, provides a much richer setting for the multiscale analysis of geophysical flows. In what we are going to say, we make the choice of focusing on works which deal with viscous flows and which perform the asymptotic study for general ill-prepared initial data. However, the literature about the subject is more ample than that.

First results in that direction were presented in [START_REF] Feireisl | A singular limit for compressible rotating fluids[END_REF], [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF] for the barotropic Navier-Stokes system. There, the authors investigated the combined effect of a strong Coriolis force (low Rossby number limit) and of weak compressibility of the fluid (low Mach number limit), under the scaling

Ma = ε m and Ro = ε , with m ≥ 1 , (1.2)
where ε ∈ ]0, 1] is a small parameter, which one wants to let go to 0 in order to derive an asymptotic model. Notice that in [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF] the effects due to the centrifugal force were considered as well, but this imposed the severe restriction m > 10. In the case m = 1 in (1.2), the system presents an isotropic scaling, since the Rossby and Mach numbers act at the same order of magnitude and they keep in balance in the limit process. This balance takes the name of quasi-geostrophic balance, and the limit system is identified as the so-called quasi-geostrophic equation for the stream function of the target velocity field. When m > 1, instead, the pressure and Coriolis forces act at different scales, the former one having a predominant effect on the dynamics of the fluid. At the mathematical level, the anisotropy of scaling generates some complications in the analysis; in [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF] this issue was handled by the use of dispersive estimates, which allowed to show convergence to a 2-D incompressible Navier-Stokes system. We refer to [START_REF] Fanelli | Highly rotating viscous compressible fluids in presence of capillarity effects[END_REF] for a similar study in the context of capillary models. There, the choice m = 1 was made, but the anisotropy was given by the scaling fixed for the internal forces term (the so-called Korteweg stress tensor). In addition, we refer to [START_REF] Fanelli | Incompressible and fast rotation limit for barotropic Navier-Stokes equations at large Mach numbers[END_REF] for the case of large Mach numbers, namely for the case when 0 ≤ m < 1 in (1.2). Since, in that instance, the pressure gradient is not strong enough to compensate the Coriolis force, in order to find some interesting dynamics in the limit one has to introduce a penalisation of the bulk viscosity coefficient.

In [START_REF] Feireisl | Scale interactions in compressible rotating fluids[END_REF], [START_REF] Feireisl | Multiple scales and singular limits for compressible rotating fluids with general initial data[END_REF] the effects of gravity were added, under the scaling

Fr = ε n , with 1 ≤ n < m 2 . (1.3)
In particular, in those works one had m > 2. As before, a planar incompressible Navier-Stokes system was identified as the limiting system, but, as already mentioned, the anisotropy of scaling created several difficulties in the analysis. We refer to [START_REF] Kwon | Multiscale analysis in the compressible rotating and heat conducting fluids[END_REF] and [START_REF] Kwon | Derivation of geostrophic equations as a rigorous limit of compressible rotating and heat conducting fluids with the general initial data[END_REF] for related studies in the context of the full Navier-Stokes-Fourier system, under the same choices of the scaling (notice that, in [START_REF] Kwon | Derivation of geostrophic equations as a rigorous limit of compressible rotating and heat conducting fluids with the general initial data[END_REF], the case m = 1 was considered, but the gravitational force was not penalised at all). The asymptotic results of [START_REF] Feireisl | Scale interactions in compressible rotating fluids[END_REF], [START_REF] Feireisl | Multiple scales and singular limits for compressible rotating fluids with general initial data[END_REF], [START_REF] Kwon | Multiscale analysis in the compressible rotating and heat conducting fluids[END_REF] and [START_REF] Kwon | Derivation of geostrophic equations as a rigorous limit of compressible rotating and heat conducting fluids with the general initial data[END_REF] are all based on a fine combination of the relative entropy/relative energy method with dispersive estimates derived from oscillatory integrals, and a strong argument which allows to handle the ill-preparation of the data (typically, the use of relative energy estimates requires to consider well-prepared initial data). In all those works, a vanishing viscosity regime was also considered.

In our recent work [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF], devoted to the full Navier-Stokes-Fourier system in presence of stratification, we were able to improve the choice of the scaling (1.3) and take the endpoint case n = m/2, with m ≥ 1 as in (1.2). In passing, we mention that also effects of the centrifugal force were considered in [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF], but this imposed the additional constraint m ≥ 2 on the order of the Mach number (which, besides, refined the restriction in [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF]). The improvement on the orders of the scaling was possible, essentially due to a different technique employed for proving convergence, based on compensated compactness arguments. We refer to [START_REF] Gallagher | Weak convergence results for inhomogeneous rotating fluid equations[END_REF] for the first implementation of that method in the context of fast rotating fluids, to [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF], [START_REF] Fanelli | A singular limit problem for rotating capillary fluids with variable rotation axis[END_REF], [START_REF] Fanelli | Incompressible and fast rotation limit for barotropic Navier-Stokes equations at large Mach numbers[END_REF] for other applications in the case of non-homogeneous flows. In particular, the convergence is not quantitative at all, but only qualitative. This technique is purely based on the algebraic structure of the system, which allows to find smallness (and vanishing to the limit) of suitable non-linear quantities, and fundamental compactness properties for other quantities (linked to the vorticity of the fluid and to the variations of the density); such compactness properties were already put in evidence in [START_REF] Fanelli | Asymptotics of fast rotating density-dependent fluids in two space dimensions[END_REF] (see also [START_REF] Cobb | On the fast rotation asymptotics of a non-homogeneous incompressible MHD system[END_REF]) in the context of non-homogeneous incompressible fluids in fast rotation. All these features were enough to pass to the limit in the primitive system, and derive the limiting dynamics: a 2-D incompressible Navier-Stokes system when m > 1, a quasi-geostrophic equation for the stream function of the limit velocity when m = 1.

An important point of the study performed in [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF] is that the scaling n = m/2 allowed to deduce some stratification effect in the limit. More precisely, although the limit dynamics was horizontal and two-dimensional, as dictated by the Taylor-Proudman theorem, stratification appeared in the functions representing departures of the density and temperature from the respective equilibria. On the contrary, in previous works like [START_REF] Feireisl | Scale interactions in compressible rotating fluids[END_REF], [START_REF] Feireisl | Multiple scales and singular limits for compressible rotating fluids with general initial data[END_REF], [START_REF] Kwon | Multiscale analysis in the compressible rotating and heat conducting fluids[END_REF], [START_REF] Kwon | Derivation of geostrophic equations as a rigorous limit of compressible rotating and heat conducting fluids with the general initial data[END_REF], based on the scaling (1.3), stratification effects were completely absent. In this sense, we call the endpoint case n = m/2 "critical".

To conclude this part, we mention that all the results quoted so far concern various regimes of low stratification, meaning that, according to the scaling in (1.2), (1.3), one has Ma Fr -→ 0 when ε → 0 + .

The strong stratification regime, namely when the ratio Ma/Fr is of order O(1), is particularly delicate for fast rotating fluids. This is in stark contrast with the results available about the derivation of the anelastic approximation, where rotation is neglected: we refer e.g. to [START_REF] Masmoudi | Rigorous derivation of the anelastic approximation[END_REF], [START_REF] Bresch | An example of low Mach (Froude) number effects for compressible flows with nonconstant density (height) limit[END_REF], [START_REF] Feireisl | On singular limits arising in the scale analysis of stratified fluid flows[END_REF] and, more recently, [START_REF] Fanelli | Low Mach number limit for degenerate Navier-Stokes equations in presence of strong stratification[END_REF] (see also [START_REF] Feireisl | Singular limits in thermodynamics of viscous fluids[END_REF] and references therein for a more detailed account of previous works). The reason for that has to be ascribed exactly to the competition between vertical stratification (due to gravity) and horizontal stability (which the Coriolis force tends to impose): in the strong stratification regime, vertical oscillations of the solution (seem to) persist in the limit, and the available techniques do not allow at present to deal with this problem in its full generality. Nonetheless, partial results have been obtained in the case of well-prepared initial data, by means of a relative entropy method: we refer to [START_REF] Feireisl | Rotating compressible fluids under strong stratification[END_REF] for the first result, where the mean motion is derived, and to [START_REF] Bocchi | Anisotropy and stratification effects in the dynamics of fast rotating compressible fluids[END_REF] for an analysis of Ekman boundary layers in that framework.

A short overview of the contents of the paper

In the present work, we continue our investigation from [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF], devoted to the multiscale analysis of systems for geophysical fluids and the derivation of reduced models. For clarity of exposition, we neglect here heat transfer processes in the fluid, and focus on the classical barotropic Navier-Stokes system as the primitive system; the more general case of the Navier-Stokes-Fourier system can be handled at the price of some additional technicalities (as done in [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF]). Also, we simplify the model by neglecting effects due to the centrifugal force. On the one hand, this choice is not dramatic from the physical viewpoint (see the discussion in [START_REF] Cushman-Roisin | Introduction to Geophysical Fluid Dynamics[END_REF], for instance); on the other hand, we could include the presence of the centrifugal force, after imposing some additional restrictions on the order of magnitude of the Mach number. We refer to Section 2 below for the presentation of the precise equations we are going to consider in this paper.

We work in the context of global in time finite energy weak solutions to the barotropic Navier-Stokes system with Coriolis force, which provides a good setting for studying singular limits for that system. We consider the general case of ill-prepared initial data.

Our goal here is to go beyond the "critical" choice Fr = √ Ma performed in [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF], and investigate other regimes where the stratification has an even more important effect. More precisely, we fix the following choice for the parameters m and n appearing in (1.2) and (1. The restriction n < 1 when m = 1 is imposed in order to avoid a strong stratification regime: as already mentioned before, it is not clear at present how to deal with this case for general ill-prepared initial data, as all the available techniques seem to break down in that case. On the other hand, the restriction 2 n ≤ m + 1 (for m > 1) looks to be of technical nature. However, it comes out naturally in at least two points of our analysis, and it is not clear to us if, and how, it is possible to bypass it and consider the remaining range of values (m + 1)/2 < n < m. Let us point out that, in our considerations, the relation n < m holds always true, so we will always work in a low stratification regime. At the qualitative level, our main results will be quite similar to the ones presented in [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF], in particular the limit dynamics will be the same (after distinguishing the two cases m > 1 and m = 1). We refer again to Section 2 for the precise statements. In this paper, the main point we put the accent on is how using in a fine way not only the structure of the system, but also the precise structure of each term in order to pass to the limit. To be more precise, the fact of considering values of n going above the threshold 2n = m is made possible thanks to special algebraic cancellations involving the gravity term in the system of wave equations. Such cancellations owe very much to the peculiar form of the gravitational force, which depends on the vertical variable only, and they do not appear, in general, if one wants to consider the action of different forces on the system. As one may easily guess, the case 2n = m + 1 is more involved: indeed, this choice of the scaling implies the presence of an additional bilinear term of order O(1) in the computations; in turn, this term might not vanish in the limit, differently to what happens in the case 2n < m + 1. In order to see that this does not occur, and that this term indeed disappears in the limit process, one has to use more thoroughly the structure of the system to control the oscillations.

Before moving on, let us give an outline of the paper. In Section 2 we collect our assumptions and we state our main results. In Section 3 we show the main consequences of the finite energy condition on the family of weak solutions we are going to consider. Namely, we derive uniform bounds in suitable norms, which allow us to extract weak-limit points, and we explore the constraints those limit points have to satisfy. In Sections 4 and 5, we complete the proof of our main results, showing convergence in the weak formulation of the equations in the cases m > 1 and m = 1, respectively, via a compensated compactness argument. We conclude the paper with Appendix A, where we present some tools from Littlewood-Paley decomposition, which have been needed in our analysis. Given p ∈ [1, +∞], by L p (B) we mean the classical space of Lebesgue measurable functions g, where |g| p is integrable over B (with the usual modifications for the case p = +∞). Sometimes, given T > 0 and (p, q) ∈ [1, +∞] 2 , we use the symbol L p T (L q ) to denote the space L p 0, T ; L q (B) . Given k ≥ 0, we denote by W k,p (B) the Sobolev space of functions which belongs to L p (B) together with all their derivatives up to order k. When p = 2, we set W k,2 (B) = H k (B). For the sake of simplicity, we will often omit from the notation the set B, that we will explicitly point out if needed.

In the whole paper, the symbols c and C will denote generic multiplicative constants, which may change from line to line, and which do not depend on the small parameter ε. Sometimes, we will explicitly point out the quantities on which these constants depend, by putting them inside brackets.

Let f ε 0<ε≤1 be a sequence of functions in a normed space Y . If this sequence is bounded in Y , we use the notation

f ε ε ⊂ Y .
Next, let us introduce some notation specific to fluids in fast rotation.

If B is a domain in R 3 , we decompose x ∈ B into x = (x h , x 3 ), with x h ∈ R 2 denoting its horizontal component. Analogously, for a vector-field v = (v 1 , v 2 , v 3 ) ∈ R 3 , we set v h = (v 1 , v 2
) and we define the differential operators ∇ h and div h as the usual operators, but acting just with respect to x h . In addition, we define the operator

∇ ⊥ h := -∂ 2 , ∂ 1 .
Finally, the symbol H denotes the Helmholtz projector onto the space of solenoidal vector fields in B, while H h denotes the Helmholtz projection on R 2 . Observe that, in the sense of Fourier multipliers, one has

H h f = -∇ ⊥ h (-∆ h ) -1 curl h f
. Moreover, since we will deal with a periodic problem in the x 3 -variable, we also introduce the following decomposition: for a vector-field X, we write

X(x) = X (x h ) + X(x) with X (x h ) := 1 |T 1 | T 1 X(x h , x 3 ) dx 3 , (1.5)
where T 1 := [-1, 1]/ ∼ is the one-dimensional flat torus (here ∼ denotes the equivalence relation which identifies -1 and 1) and T 1 denotes its Lebesgue measure. Notice that X has zero vertical average, and therefore we can write X(x) = ∂ 3 Z(x) with Z having zero vertical average as well.

The primitive system

As already said in the introduction, in this paper we assumed that the motion of the fluid is described by a rescaled version of the barotropic Navier-Stokes system with Coriolis and gravitational forces. Thus, given a small parameter ε ∈ ]0, 1], the system reads as follows:

∂ t ε + div ( ε u ε ) = 0 (NSC 1 ε ) ∂ t ( ε u ε ) + div ( ε u ε ⊗ u ε ) + 1 ε e 3 × ε u ε + 1 ε 2m ∇ x p( ε ) = div S(∇ x u ε ) + ε ε 2n ∇ x G , (NSC 2 ε )
where we recall that m and n are taken according to (1.4). The unknowns in the previous equations are the density ε = ε (t, x) ≥ 0 of the fluid and its velocity field

u ε = u ε (t, x) ∈ R 3 , where t ∈ R + and x ∈ Ω := R 2 × ]0, 1[. The viscous stress tensor in (NSC 2 ε
) is given by Newton's rheological law

S(∇ x u ε ) = µ ∇ x u ε + ∇ T x u ε - 2 3 div u ε Id Id Id + η div u ε Id Id Id , (2.1)
where µ > 0 is the shear viscosity and η ≥ 0 represents the bulk viscosity. The term e 3 × ε u ε takes into account the (strong) Coriolis force acting on the fluid. As for the gravitational force, it is physically relevant to assume that

G(x) = -x 3 . (2.2)
The precise expression of G will be useful in some computations below, although some generalisations are certainly possible.

The system is supplemented with complete slip boundary conditions, namely

u ε • n |∂Ω = 0 and [S(∇ x u ε )n] × n |∂Ω = 0 , (2.3)
where n denotes the outer normal to the boundary ∂Ω = {x 3 = 0} ∪ {x 3 = 1}. Notice that this is a true simplification, because it avoid complications due to the presence of Ekman boundary layers, when passing to the limit ε → 0 + .

Remark 2.1

As is well-known (see e.g. [START_REF] Ebin | Viscous fluids in a domain with frictionless boundary[END_REF]), equations (NSC 1 ε )-(NSC 2 ε ), supplemented by the complete slip boundary boundary conditions from (2.3), can be recasted as a periodic problem with respect to the vertical variable, in the new domain

Ω = R 2 × T 1 , with T 1 := [-1, 1]/ ∼ ,
where ∼ denotes the equivalence relation which identifies -1 and 1. Indeed, the equations are invariant if we extend ρ and u h as even functions with respect to x 3 , and u 3 as an odd function.

In what follows, we will always assume that such modifications have been performed on the initial data, and that the respective solutions keep the same symmetry properties. Now we need to impose structural restrictions on the pressure function p. We assume that

p ∈ C 1 [0, ∞) ∩ C 2 (0, ∞), p(0) = 0, p ( ) > 0 for all ≥ 0 . (2.4)
Additionally to (2.4), we require that

exists γ > 3 2 such that lim →+∞ p ( ) γ-1 = p ∞ > 0 . (2.5)
Without loss of generality, we can suppose that p has been renormalised so that p (1) = 1.

Equilibrium states

Next, we focus our attention on the so-called equilibrium states. For each value of ε ∈ ]0, 1] fixed, the equilibria of system (NSC

1 ε )-(NSC 2 ε ) consist of static densities ε satisfying ∇ x p( ε ) = ε 2(m-n) ε ∇ x G in Ω . (2.6)
Equation (2.6) identifies ε up to an additive constant: taking the target density to be 1, we get

H ( ε ) = ε 2(m-n) G + H (1) , where H( ) = 1 p(z) z 2 dz . (2.7)
Notice that relation (2.7) implies that

H ( ) = p ( ) and H (1) = 1 .
Therefore, we infer that, whenever m ≥ 1 and m > n as in the present paper, for any x ∈ Ω one has ε (x) -→ 1 in the limit ε → 0 + . More precisely, the next statement collects all the necessary properties of the static states. It corresponds to Lemma 2.3 and Proposition 2.5 of [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF].

Proposition 2.2 Let the gravitational force G be given by (2.2). Let ε 0<ε≤1 be a family of static solutions to equation (2.6)

in Ω = R 2 × ]0, 1[.
Then, there exist an ε 0 > 0 and a 0 < ρ * < 1 such that ε ≥ ρ * for all ε ∈ ]0, ε 0 ] and all x ∈ Ω. In addition, for any ε ∈ ]0, ε 0 ], one has:

| ε (x) -1| ≤ C ε 2(m-n) , for a constant C > 0 which is uniform in x ∈ Ω and in ε ∈ ]0, 1].
Without loss of any generality, we can assume that ε 0 = 1 in Proposition 2.2.

In light of this analysis, it is natural to try to solve system (NSC

Initial data and finite energy weak solutions

In view of the boundary conditions (2.8) "at infinity", we assume that the initial data are close (in a suitable sense) to the equilibrium states ε that we have just identified. Namely, we consider initial densities of the following form:

0,ε = ε + ε m (1) 0,ε . (2.9)
For later use, let us introduce also the following decomposition of the initial densities:

0,ε = 1 + ε 2(m-n) R 0,ε with R 0,ε = r ε + ε 2n-m (1) 0,ε , r ε := ε -1 ε 2(m-n) .
(2.10)

Notice that the r ε 's are in fact data of the system, since they only depend on p and G.

We suppose the density perturbations

(1) 0,ε to be measurable functions and satisfy the control

sup ε∈ ]0,1] (1) 0,ε (L 2 ∩L ∞ )(Ω) ≤ c , (2.11)
together with the "mean-free condition"

Ω (1) 0,ε dx = 0 .
As for the initial velocity fields, we assume the following uniform bound:

sup ε∈ ]0,1] ε u 0,ε L 2 (Ω) ≤ c .
(2.12)

Remark 2.3 In view of Proposition 2.2, the condition in (2.12) immediately implies that

sup ε∈ ]0,1] u 0,ε L 2 (Ω) ≤ c .
Thanks to the previous uniform estimates, up to extraction, we can identify the limit points

(1) 0

:= lim ε→0 (1) 0,ε weakly- * in L ∞ (Ω) ∩ L 2 (Ω) (2.13)
u 0 := lim ε→0 u 0,ε weakly in L 2 (Ω) . (2.14)
At this point, let us specify better what we mean for finite energy weak solution (see [START_REF] Feireisl | Singular limits in thermodynamics of viscous fluids[END_REF] for details).

Definition 2.4

Let Ω = R 2 × ]0, 1[ . Fix T > 0 and ε > 0. Let ( 0,ε , u 0,ε ) be an initial datum satisfying (2.9)-(2.12). We say that the couple ( ε , u ε ) is a finite energy weak solution of the system (NSC 1 ε )-(NSC 2 ε ) in ]0, T [ ×Ω, supplemented with the boundary conditions (2.3) and far field conditions (2.8), related to the initial datum ( 0,ε , u 0,ε ), if the following conditions hold:

(i) the functions ε and u ε belong to the class

ε ≥ 0 , ε -ε ∈ L ∞ 0, T ; L 2 + L γ (Ω) , u ε ∈ L 2 0, T ; H 1 (Ω) , u ε • n |∂Ω = 0 ;
(ii) the equations have to be satisfied in a distributional sense:

- T 0 Ω ( ε ∂ t ϕ + ε u ε • ∇ x ϕ) dxdt = Ω 0,ε ϕ(0, •) dx (2.15) for any ϕ ∈ C ∞ c ([0, T [ ×Ω) and T 0 Ω -ε u ε • ∂ t ψ -ε [u ε ⊗ u ε ] : ∇ x ψ + 1 ε e 3 × ( ε u ε ) • ψ - 1 ε 2m p( ε )div ψ dxdt (2.16) = T 0 Ω -S(∇ x u ε ) : ∇ x ψ + 1 ε 2n ε ∇ x G • ψ dxdt + Ω 0,ε u 0,ε • ψ(0, •) dx for any test function ψ ∈ C ∞ c ([0, T [ ×Ω; R 3 ) such that ψ • n |∂Ω = 0;
(iii) the energy inequality holds for almost every t ∈ (0, T ):

Ω 1 2 ε |u ε | 2 (t) dx + 1 ε 2m Ω E ( ε , ε ) (t) dx + t 0 Ω S(∇ x u ε ) : ∇ x u ε dxdτ (2.17) ≤ Ω 1 2 0,ε |u 0,ε | 2 dx + 1 ε 2m Ω E ( 0,ε , ε ) dx ,
where the function

E (ρ, ε ) := H(ρ) -(ρ -ε ) H ( ε ) -H( ε ) (2.18)
is the relative internal energy of the fluid, with H given by (2.7).

The solutions is global if the previous conditions hold for all T > 0.

Under the assumptions fixed above, for any fixed value of the parameter ε ∈ ]0, 1], the existence of a global in time finite energy weak solution ( ε , u ε ) to system (NSC 1 ε )-(NSC 2 ε ), related to the initial datum ( 0,ε , u 0,ε ), in the sense of the previous definition, can be proved as in the classical case, see e.g. [START_REF] Lions | Mathematical topics in Fluid Mechanics[END_REF], [START_REF] Feireisl | Dynamics of viscous compressible fluids[END_REF]. Notice that the mapping t → ( ε u ε )(t, •) is weakly continuous, and one has ( ε

) |t=0 = 0,ε together with ( ε u ε ) |t=0 = 0,ε u 0,ε .
We remark also that, in view of (NSC 1 ε ), the total mass is conserved in time, in the following sense: for almost every t ∈ [0, +∞[ , one has

Ω ε (t) -ε dx = 0 .
To conclude, we point out that, in our framework of finite energy weak solutions to the primitive system, inequality (2.17) will be the only tool to derive uniform estimates for the family of weak solutions we are going to consider.

Main results

We can now state our main results. We point out that, due to the scaling (1.4), the relation m > n is always true, so we will always be in a low stratification regime.

The first statement concerns the case when the effects linked to the pressure term are predominant with respect to the fast rotation, i.e. m > 1.

Theorem 2.5 Let Ω = R 2 × ]0, 1[ and G ∈ W 1,∞ (Ω) be as in (2.2). Take m > 1 and m + 1 ≥ 2n > m.
For any fixed value of ε ∈ ]0, 1], let initial data ( 0,ε , u 0,ε ) verify the hypotheses fixed in Paragraph 2.1.2, and let ( ε , u ε ) be a corresponding weak solution to system (NSC 1 ε )-(NSC 2 ε ), supplemented with the structural hypotheses (2.1) on S(∇ x u ε ) and with boundary conditions (2.3) and far field conditions (2.8). Let u 0 be defined as in (2.14).

Then, for any T > 0 one has the convergence properties

ε → 1 strongly in L ∞ 0, T ; L min{2,γ} loc (Ω) u ε U weakly in L 2 0, T ; H 1 (Ω) ,
where U = (U h , 0), with U h = U h (t, x h ) such that div h U h = 0. In addition, the vector field U h is a weak solution to the following homogeneous incompressible Navier-Stokes system in R + × R 2 ,

∂ t U h + div h U h ⊗ U h + ∇ h Γ -µ∆ h U h = 0 , (2.19)
for a suitable pressure function Γ ∈ D (R + × R 2 ) and related to the initial condition

U |t=0 = H h u h 0 .
When m = 1, the Mach and Rossby numbers have the same order of magnitude, and they keep in balance in the whole asymptotic process, realising in this way the so-called quasi-geostrophic balance in the limit. The next statement is devoted to this case.

Theorem 2.6 Let Ω = R 2 × ]0, 1[ and let G ∈ W 1,∞ (Ω) be as in (2.2). Take m = 1 and 1/2 < n < 1.
For any fixed value of ε ∈ ]0, 1], let initial data ( 0,ε , u 0,ε ) verify the hypotheses fixed in Paragraph 2.1.2, and let ( ε , u ε ) be a corresponding weak solution to system (NSC 1 ε )-(NSC 2 ε ), supplemented with the structural hypotheses (2.1) on S(∇ x u ε ) and with boundary conditions (2.3) and far field conditions (2.8).

Let

(1) 0 , u 0 be defined as in (2.13) and (2.14). Then, for any T > 0 one has the following convergence properties:

ε → 1 strongly in L ∞ 0, T ; L min{2,γ} loc (Ω) (1)
ε := ε -ε ε * (1)
weakly- * in

L ∞ 0, T ; L 2 + L γ (Ω) u ε U weakly in L 2 0, T ; H 1 (Ω) ,
where, as above, U = (U h , 0), with U h = U h (t, x h ) such that div h U h = 0. Moreover, one has the relation 1) , and (1) satisfies (in the weak sense) the quasi-geostrophic equation

U h = ∇ ⊥ h (
∂ t (1) -∆ h (1) -∇ ⊥ h (1) • ∇ h ∆ h (1) + µ∆ 2 h (1) = 0 , (2.20)
supplemented with the initial condition

(1) -∆ h (1) |t=0 =
(1) 0

-curl h u h 0 .

Consequences of the energy inequality

In Definition 2.4, we have postulated that the family of weak solutions ε , u ε ε considered in Theorems 2.5 and 2.6 satisfies the energy inequality (2.17). In this section we take advantage of that fact to infer uniform bounds for ε , u ε ε , see Section 3.1. Thanks to those bounds, we can extract weak-limit points of the sequence of solutions and deduce some properties those limit points have to satisfy, see Section 3.2.

Uniform bounds and weak limits

This section is devoted to establish uniform bounds on the sequence ε , u ε ε . This can be done as in the classical case (see e.g. [START_REF] Feireisl | Singular limits in thermodynamics of viscous fluids[END_REF] for details), since the Coriolis term does not contribute to the total energy balance of the system. However, for the reader's convenience, let us present some details.

To begin with, let us introduce a partition of the space domain Ω into the so-called "essential" and "residual" sets. For this, for t > 0 and for all ε ∈ ]0, 1], we define the sets

Ω ε ess (t) := x ∈ Ω ε (t, x) ∈ [1/2 ρ * , 2] , Ω ε res (t) := Ω \ Ω ε ess (t) ,
where the positive constant ρ * > 0 has been defined in Proposition 2.2. Then, given a function h, we write

h = [h] ess + [h] res , where [h] ess := h 1 Ω ε ess (t) .
Here above, 1 A denotes the characteristic function of a set A ⊂ Ω.

Next, we observe that

E ρ(t, x), ε (x) ess ∼ [ρ -ε (x)] 2 ess and E ρ(t, x), ε (x) res ≥ C 1 + ρ(t, x) γ res ,
where ε is the static density state identified in Section 2.1.1 and E is given by (2.18). Here above, the multiplicative constants are all strictly positive and may depend on ρ * and we agree to write A ∼ B whenever there exists a "universal" constant c > 0 such that

(1/c) B ≤ A ≤ c B.
Thanks to the previous observations, we easily see that, under the assumptions fixed in Section 2 on the initial data, the right-hand side of (2.17) is uniformly bounded for all ε ∈ ]0, 1]. Specifically, we have

Ω 1 2 0,ε |u 0,ε | 2 dx + 1 ε 2m Ω E ( 0,ε , ε ) dx ≤ C .
Owing to the previous inequalities and the finite energy condition (2.17) on the family of weak solutions, it is quite standard to derive, for any time T > 0 fixed and any ε ∈ ]0, 1], the following estimates:

sup t∈[0,T ] √ ε u ε L 2 (Ω; R 3 ) ≤ c (3.1) sup t∈[0,T ] ε -ε ε m ess (t) L 2 (Ω) ≤ c (3.2) sup t∈[0,T ] Ω 1 M ε res [t] dx ≤ c ε 2m (3.3) sup t∈[0,T ] Ω [ ε ] γ res (t) dx ≤ c ε 2m (3.4) T 0 ∇ x u ε + ∇ T x u ε - 2 3 div u ε Id Id Id 2 L 2 (Ω; R 3×3 ) dt ≤ c . (3.5)
We refer to [START_REF] Feireisl | Singular limits in thermodynamics of viscous fluids[END_REF] (see also [START_REF] Feireisl | A singular limit for compressible rotating fluids[END_REF], [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF], [START_REF] Fanelli | Incompressible and fast rotation limit for barotropic Navier-Stokes equations at large Mach numbers[END_REF] and [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF]) for the details of the computations.

Owing to (3.5) and a generalisation of the Korn inequality (see e.g. Chapter 10 of [START_REF] Feireisl | Singular limits in thermodynamics of viscous fluids[END_REF]), we gather that ∇u ε ε ⊂ L 2 T (L 2 ). On the other hand, by arguing as in [START_REF] Feireisl | A singular limit for compressible rotating fluids[END_REF], we can use (3.1), (3.3) and (3.4) to deduce that also u ε ε ⊂ L 2 T (L 2 ). Putting those bounds together, we finally infer that

T 0 u ε 2 H 1 (Ω; R 3 ) dt ≤ c . (3.6)
In particular, there exist U ∈ L 2 loc R + ; H 1 (Ω; R 3 ) such that, up to a suitable extraction (not relabelled here), we have

u ε U in L 2 loc R + ; H 1 (Ω; R 3 ) . (3.7)
Let us move further and consider the density functions. The previous estimates on the density tell us that we must find a finer decomposition for the densities. As a matter of fact, for any time T > 0 fixed, we have

ε -1 L ∞ T (L 2 +L γ +L ∞ ) ≤ c ε 2(m-n) . (3.8)
In order to see (3.8) we write

| ε -1| ≤ | ε -ε | + | ε -1| . (3.9) From (3.2), we infer that ε -ε ess is of order O(ε m ) in L ∞ T (L 2 ).
For the residual part of the same term, we can use (3.4) to discover that it is of order O(ε 2m/γ ). Observe that, if 1 < γ < 2, the higher order is O(ε m ), whereas, in the case γ ≥ 2, by use of (3.4) and (3.3) again, it is easy to get

[ ε -ε ] res 2 L ∞ T (L 2 ) ≤ C ε 2m . (3.10)
Finally, we apply Proposition 2.2 to control the last term in the right-hand side of (3.9). In the end, estimate (3.8) is proved.

This having been established, and keeping in mind the notation introduced in (2.9) and (2.10), we can introduce the density oscillation functions

R ε := ε -1 ε 2(m-n) = r ε + ε 2n-m (1)
ε ,

where we have defined

ε (t, x) := ε -ε ε m and r ε (x) := ε -1 ε 2(m-n) .
(3.11)

Thanks again to (3.2), (3.4) and Proposition 2.2, we see that the above quantities verify the following uniform bounds, for any time T > 0 fixed:

sup ε∈ ]0,1] (1) ε L ∞ T (L 2 +L γ (Ω)) ≤ c and sup ε∈ ]0,1] r ε L ∞ (Ω) ≤ c . (3.12)
In view of the previous properties, there exist (1) ∈ L ∞ loc (R + ; L 2 + L γ ) and r ∈ L ∞ such that (up to the extraction of a new suitable subsequence), for any T > 0 we have

(1) ε * (1) weakly- * in L ∞ (0, T ; L 2 + L γ (Ω)) and r ε * r weakly- * in L ∞ (Ω). (3.13)
In particular, we get

R ε * r weakly- * in L ∞ 0, T ; L min{γ,2} loc (Ω) .
Remark 3.1 Observe that, owing to (3.10), when γ ≥ 2 we get

sup ε∈ ]0,1] (1) ε L ∞ T (L 2 ) ≤ c .
Therefore, in that case we actually have that (1) ∈ L ∞ T (L 2 ) and that

(1) ε * (1) in L ∞ T (L 2 ). Analogously, when γ ≥ 2 we also get ε -1 L ∞ T (L 2 +L ∞ ) ≤ c ε 2(m-n) .

Constraints on the limit

In this section, we establish some properties that the limit points of the family ε , u ε ε , which have been identified here above, have to satisfy. We first need a preliminary result about the decomposition of the pressure function, which will be useful in the following computations.

Lemma 3.2 Let (m, n) ∈ R 2 verify the condition m + 1 ≥ 2n > m ≥ 1.
Let p be the pressure term satisfying the structural hypotheses (2.4)-(2.5). Then, for any ε ∈ ]0, 1], one has

1 ε 2m ∇ x p( ε ) -p( ε ) = 1 ε m ∇ x p (1) (1) ε + 1 ε 2n-m ∇ x Π ε , (3.14)
where the functions

(1) ε have been introduced in (3.11) and, for all T > 0, the family Π ε ε verifies the uniform bound

Π ε L ∞ T (L 1 +L 2 +L γ ) ≤ C . (3.15)
When γ ≥ 2, one can dispense of the space L γ in the above control of Π ε ε .

Proof: We start by writing simple algebraic computations:

1 ε 2m ∇ x p( ε ) -p( ε ) = 1 ε 2m ∇ x p( ε ) -p( ε ) -p ( ε )( ε -ε ) + 1 ε m ∇ x p ( ε ) -p (1) (1) ε + 1 ε m ∇ x p (1) (1) ε . (3.16)
We start by analysing the first term on the right-hand side of (3.16). For the essential part, we can employ a Taylor expansion to write

[p( ε ) -p( ε ) -p ( ε )( ε -ε )] ess = p (z ε )( ε -ε ) 2
ess , where z ε is a suitable point between ε and ε . Thanks to the uniform bound (3.2), we have that this term is of order O(ε 2m ) in L ∞ T (L 1 ), for any T > 0 fixed. For the residual part, we can use (3.3) and (3.4), together with the boundedness of the profiles ε (keep in mind Proposition 2.2), to deduce that

[p( ε ) -p( ε ) -p ( ε )( ε -ε )] res L ∞ T (L 1 ) ≤ C ε 2m .
We refer to e.g. Lemma 4.1 of [START_REF] Fanelli | Incompressible and fast rotation limit for barotropic Navier-Stokes equations at large Mach numbers[END_REF] for details.

In a similar way, a Taylor expansion for the second term on the right-hand side of (3.16) gives

p ( ε ) -p (1) (1) ε = p (η ε )( ε -1) (1)
ε where η ε is a suitable point between ε and 1. Owing to Proposition 2.2 again and to bound (3.12), we infer that this term is of order O(ε 2(m-n) ) in L ∞ T (L 2 + L γ ), for any time T > 0 fixed. Then, defining

Π ε := 1 ε 2(m-n) p( ε ) -p( ε ) ε m -p (1) (1)
ε we have the control (3.15). The final statement concerning the case γ ≥ 2 easily follows from Remark 3.1. This completes the proof of the lemma. Remark 3.3 Notice that the last term appearing in (3.14) is singular in ε. This is in stark contrast with the situation considered in previous works, see e.g. [START_REF] Feireisl | A singular limit for compressible rotating fluids[END_REF], [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF], [START_REF] Feireisl | Multiple scales and singular limits for compressible rotating fluids with general initial data[END_REF] and [START_REF] Fanelli | Incompressible and fast rotation limit for barotropic Navier-Stokes equations at large Mach numbers[END_REF]. However, its gradient structure will play a fundamental role in the computations below.

This having been pointed out, we can now analyse the constraints on the weak-limit points (1) , U , identified in relations (3.7) and (3.13) above.

The case of large values of the Mach number: m > 1

We start by considering the case of anisotropic scaling, namely m > 1 and m + 1 ≥ 2n > m. Notice that, in particular, one has m > n.

Proposition 3.4 Let m > 1 and m + 1 ≥ 2n > m in (NSC 1 ε )-(NSC 2 ε ). Let ( ε , u ε )
ε be a family of weak solutions, related to initial data ( 0,ε , u 0,ε ) ε verifying the hypotheses of Section 2.1.2. Let ( (1) , U ) be a limit point of the sequence

(1) ε , u ε ε , as identified in Section 3.1. Then U = U h , 0 , with U h = U h (t, x h ) and div h U h = 0 , (3.17)
∇ x (1) = 0 in D (R + × Ω) . (3.18)
Proof: First of all, let us consider the weak formulation of the mass equation (NSC

1 ε ). Take a test function ϕ ∈ C ∞ c R + × Ω and denote [0, T ] × K := supp ϕ.
Then by (2.15) we have

- T 0 K ε -1 ∂ t ϕ dxdt - T 0 K ε u ε • ∇ x ϕ dxdt = K 0,ε -1 ϕ(0, • ) dx .
We can easily pass to the limit in this equation, thanks to the strong convergence ε -→ 1, provided by (3.8), and the weak convergence of u ε in L 2 T L 6 loc , provided by (3.7) and Sobolev embeddings. Notice that one always has 1/γ + 1/6 ≤ 1. In this way, we find

- T 0 K U • ∇ x ϕ dxdt = 0 (3.19)
for ϕ taken as above. Since the choice of ϕ is arbitrary, we obtain that div U = 0 a.e. in R + × Ω .

(3.20)

Next, we test the momentum equation (NSC 2 ε ) on ε m φ, for a smooth compactly supported φ. Using the uniform bounds established in Section 3.1, it is easy to see that the term presenting the time derivative, the viscosity term and the convective term all converges to 0, in the limit ε → 0 + . Since m > 1, also the Coriolis term vanishes when ε → 0 + . It remains us to consider the pressure and gravity terms in the weak formulation (2.16) of the momentum equation: using relation (2.6), we see that we can couple them to write

1 ε 2m ∇ x p( ε ) - 1 ε 2n ε ∇ x G = 1 ε 2m ∇ x p( ε ) -p( ε ) -ε m-2n (1) ε ∇ x G . (3.21)
By (3.12) and the fact that m > n, we readily see that the last term in the right-hand side of (3.21) converges to 0, when tested against any smooth compactly supported ε m φ. At this point, we use Lemma 3.2 to treat the first term on the right-hand side of (3.21). So, taking φ ∈ C ∞ c ([0, T [ ×Ω) (for some T > 0), we test the momentum equation against ε m φ: using (3.13), in the limit ε → 0 + we find that T 0 Ω p (1) (1) div φ dxdt = 0 .

Recalling that p (1) = 1, the previous relation implies (3.18) for (1) . In particular, that relation implies that (1) (t, x) = c(t) for almost all (t, x) ∈ R + × Ω, for a suitable function c = c(t) depending only on time. Now, in order to see effects due to the fast rotation in the limit, we need to "filter out" the contribution coming from the low Mach number. To this end, we test (NSC 2 ε ) on ε φ, where this time we take φ = curl ψ, for some smooth compactly supported ψ ∈ C ∞ c [0, T [ ×Ω , with T > 0. Once again, by uniform bounds we infer that the ∂ t term, the convective term and the viscosity term all converge to 0 when ε → 0 + . As for the pressure and the gravitational force, we argue as in (3.21). Since the structure of φ kills any gradient term, we are left with the convergence of the integral

T 0 Ω ε m-2n+1 (1) ε ∇ x G • φ dx dt -→ δ 0 (m -2n + 1) T 0 Ω (1) ∇ x G • φ dx dt ,
where δ 0 (ζ) = 1 if ζ = 0, δ 0 (ζ) = 0 otherwise. Finally, arguing as done for the mass equation, we see that the Coriolis term converges to the integral T 0 Ω e 3 × U • φ. Consider the case m+1 > 2n for a while. Passing to the limit for ε → 0 + , we find that H (e 3 × U ) = 0, which implies that e 3 × U = ∇ x Φ, for some potential function Φ. From this relation, one easily deduces that Φ = Φ(t, x h ), i.e. Φ does not depend on x 3 , and that the same property is inherited by U h = U 1 , U 2 , i.e. one has U h = U h (t, x h ). Furthermore, since U h = -∇ ⊥ Φ, we get that div h U h = 0. At this point, we combine this fact with (3.20) to infer that ∂ 3 U 3 = 0; but, thanks to the boundary condition (2.3), we must have U • n |∂Ω = 0, which implies that U 3 has to vanish at the boundary of Ω. Thus, we finally deduce that U 3 ≡ 0, whence (3.17) follows.

Now, let us focus on the case when m + 1 = 2n. The previous computations show that, when ε → 0 + , we get

e 3 × U + (1) ∇ x G = ∇ x Φ in D (R + × Ω) , (3.22)
for a new suitable function Φ. However, owing to (3.18), we see that (1) 1) G ; hence, the previous relations can be recasted as e 3 × U = ∇ x Φ, for a new scalar function Φ. Therefore, the same analysis as above applies, allowing us to gather (3.17) also in the case m + 1 = 2n.

∇ x G = ∇ x (

The case m = 1

Now we focus on the case m = 1. In this case, the fast rotation and weak compressibility effects are of the same order: this allows to reach the so-called quasi-geostrophic balance in the limit.

Proposition 3.5 Take m = 1 and 1/2 < n <

1 in system (NSC 1 ε )-(NSC 2 ε ). Let ( ε , u ε ) ε be a family of weak solutions to (NSC 1 ε )-(NSC 2 ε )
, associated with initial data ( 0,ε , u 0,ε ) verifying the hypotheses fixed in Section 2.1.2. Let ( (1) , U ) be a limit point of the sequence

(1)
ε , u ε ε , as identified in Section 3.1. Then,

(1) = (1) (t, x h ) and U = U h , 0 , with U h = ∇ ⊥ h (1) a.e. in R + × R 2 . (3.23)
In particular, one has U h = U h (t, x h ) and div h U h = 0.

Proof: Arguing as in the proof of Proposition 3.4, it is easy to pass to the limit in the continuity equation. In particular, we obtain again relation (3.20) for U .

Only the analysis of the momentum equation changes a bit with respect to the previous case m > 1. Now, since the most singular terms are of order ε -1 (keep in mind Lemma 3.2), we test the weak formulation (2.16) of the momentum equation against ε φ, where φ is a smooth compactly supported function. Similarly to what done above, the uniform bounds of Section 3.1 allow us to infer that the only quantity which does not vanish in the limit is the sum of the terms involving the Coriolis force, the pressure and the gravitational force: more precisely, using also Lemma 3.2 and (2.6), we have

e 3 × ε u ε + ∇ x p( ε ) -p( ε ) ε -ε 2(1-n) (1) ε ∇ x G = O(ε)
in the sense of D (R + × Ω). Following the same computations performed in the proof of Proposition 3.4, in the limit ε → 0 + it is easy to get that

e 3 × U + ∇ x p (1) (1) = 0 in D R + × Ω .
After recalling that p (1) = 1, this equality can be equivalently written as

e 3 × U + ∇ x (1) = 0 a.e. in R + × Ω .
Notice that U is in fact in L 2 loc (R + ; L 2 ), therefore so is ∇ x (1) ; hence the previous relation is in fact satisfied almost everywhere in R + × Ω.

At this point, we can repeat the same argument used in the proof of Proposition 3.4 to deduce (3.23). The proposition is thus proved.

Convergence in the case m > 1

In this section, we complete the proof of Theorem 2.5. Namely, we show convergence in the weak formulation of the primitive system, in the case when m > 1 and m + 1 ≥ 2n > m.

In Proposition 3.4, we have already seen how passing to the limit in the mass equation. However, problems arise when tackling the convergence in the momentum equation. Indeed, the analysis carried out so far is not enough to identify the weak limit of the convective term ε u ε ⊗u ε , which is highly non-linear. For proving that this term converges to the expected limit U ⊗ U , the key point is to control the strong oscillations in time of the solutions, generated by the singular terms in the momentum equation. For this, we will use a compensated compactness argument and exploit the algebraic structure of the wave system underlying the primitive equations (NSC 1 ε )-(NSC 2 ε ).

In Section 4.1, we start by giving a quite accurate description of those fast oscillations. Then, using that description, we are able, in Section 4.2, to establish two fundamental properties: on the one hand, strong convergence of a suitable quantity related to the velocity fields; on the other hand, that the other terms which do not involve that quantity tend to vanish when ε → 0 + . In turn, this allows us to complete, in Section 4.3, the proof of the convergence.

Analysis of the acoustic waves

The goal of the present subsection is to describe the fast time oscillations of the solutions. First of all, we recast our equations into a wave system. Then, we establish uniform bounds for the quantities appearing in the wave system. Finally, we apply a regularisation in space procedure for all the quantities, which is preparatory in view of the computations of Section 4.2.

Formulation of the acoustic equation

We introduce the quantity

V ε := ε u ε .
Then, straightforward computations show that we can recast the continuity equation in the form

ε m ∂ t (1) ε + div V ε = 0 , (4.1)
where

(1) ε is defined in (3.11). Next, thanks to Lemma 3.2 and the static relation (2.6), we can derive the following form of the momentum equation:

ε m ∂ t V ε + ε m-1 e 3 × V ε + p (1) ∇ x (1) ε = ε 2(m-n) (1) ε ∇ x G -∇ x Π ε (4.2) + ε m div S(∇ x u ε ) -div ( ε u ε ⊗ u ε ) .
Then, if we define

f ε := div S(∇ x u ε ) -ε u ε ⊗ u ε and g ε := (1) ε ∇ x G -∇ x Π ε , (4.3)
recalling that we have normalised the pressure function so that p (1) = 1, we can recast the primitive system (NSC 1 ε )-(NSC 2 ε ) in the following form:

ε m ∂ t (1) ε + div V ε = 0 ε m ∂ t V ε + ∇ x (1) ε + ε m-1 e 3 × V ε = ε m f ε + ε 2(m-n) g ε . (4.4)
We remark that system (4.4) has to be read in the weak sense:

for any ϕ ∈ C ∞ c [0, T [ ×Ω , one has -ε m T 0 Ω (1) ε ∂ t ϕ - T 0 Ω V ε • ∇ x ϕ = ε m Ω (1) 0,ε ϕ(0) ,
and also, for any

ψ ∈ C ∞ c [0, T [ ×Ω; R 3 such that (ψ • n) |∂Ω = 0, one has -ε m T 0 Ω V ε • ∂ t ψ - T 0 Ω (1)
ε div ψ + ε m-1 T 0 Ω e 3 × V ε • ψ = ε m Ω 0,ε u 0,ε • ψ(0) + ε m T 0 Ω f ε • ψ + ε 2(m-n) T 0 Ω g ε • ψ .
Here we use estimates of Section 3.1 in order to show uniform bounds for the solutions and the data in the wave equation (4.4). We start by dealing with the "unknown" V ε . Splitting the term into essential and residual parts, one can obtain for all T > 0:

V ε L ∞ T (L 2 +L 2γ/(γ+1) ) ≤ c . (4.5)
In the next lemma, we establish bounds for the source terms in the system of acoustic waves (4.4).

Lemma 4.1 Write f ε = div f ε and g ε = g 1 ε -∇ x Π ε
, where we have defined g 1 ε :=

ε ∇ x G and the functions Π ε have been introduced in Lemma 3.2.

For any T > 0 fixed, one has the uniform embedding properties

f ε ε ⊂ L 2 T (L 2 + L 1 ) and g 1 ε ε ⊂ L 2 T (L 2 + L γ ) .
In the case γ ≥ 2, we may dispense with the space L γ in the control of g 1 ε ε . In particular, the sequences f ε ε and g ε ε , defined in system (4.4), are uniformly bounded in the space L 2 [0, T ]; H -s (Ω) , for all s > 5/2. Proof: From (3.1), (3.5) and (3.6), we immediately infer the uniform bound for the family

f ε ε in L 2 T (L 1 + L 2 )
, from which we deduce also the uniform boundedness of f ε ε in L 2 T (H -s ), for any s > 5/2. Next, for bounding g 1 ε ε we simply use (3.12), together with Remark 3.1 when γ ≥ 2. Keeping in mind the bounds established in Lemma 3.2, the uniform estimate for g ε ε follows.

Regularization and description of the oscillations

As already mentioned in Remark 2.1, in order to apply the Littlewood-Paley theory, it is convenient to reformulate problem (NSC 1 ε )-(NSC 2 ε ) in the new domain (which we keep calling Ω, with a little abuse of notation)

Ω := R 2 × T 1 , with T 1 := [-1, 1]/ ∼ .
In addition, to avoid the appearing of (irrelevant) multiplicative constants in the computations, we suppose that the torus T 1 has been renormalised so that its Lebesgue measure is equal to 1. Now, for any M ∈ N we consider the low-frequency cut-off operator S M of a Littlewood-Paley decomposition, as introduced in equation (A.1) of the Appendix. Then, we define

(1) ε,M = S M (1) ε and V ε,M = S M V ε . (4.6)
The previous regularised quantities satisfy the following properties.

Proposition 4.2 For any T > 0, we have the following convergence properties, in the limit M → +∞:

sup 0<ε≤1 (1)
ε -

(1) ε,M L ∞ ([0,T ];H -s ) -→ 0 ∀ s > max 0, 3 1 γ - 1 2 sup 0<ε≤1 V ε -V ε,M L ∞ ([0,T ];H -s ) -→ 0 ∀ s > 3 2 γ . (4.7)
Moreover, for any M > 0, the couple (

ε,M , V ε,M ) satisfies the approximate wave equations

ε m ∂ t (1)
ε,M + div V ε,M = 0 ε m ∂ t V ε,M + ε m-1 e 3 × V ε,M + ∇ x (1) ε,M = ε m f ε,M + ε 2(m-n) g ε,M (4.8)
where (f ε,M ) ε and (g ε,M ) ε are families of smooth (in the space variables) functions satisfying, for any s ≥ 0, the uniform bounds

sup 0<ε≤1 f ε,M L 2 ([0,T ];H s ) + sup 0<ε≤1 g ε,M L ∞ ([0,T ];H s ) ≤ C(s, M ) , (4.9)
where the constant C(s, M ) depends on the fixed values of s ≥ 0 and M > 0, but not on ε > 0.

Proof: Thanks to characterization (A.2) of H s , properties (4.7) are straightforward consequences of the uniform bounds establish in Section 3.1. For instance, let us consider the functions [START_REF] Bahouri | Fourier Analysis and Nonlinear Partial Differential Equations[END_REF] ε : when γ ≥ 2, owing to Remark 3.1 one has

(1) ε ε ⊂ L ∞ T (L 2
), and then we use estimate (A.3) from the Appendix. When 1 < γ < 2, instead, we first apply the dual Sobolev embedding to infer that

(1) ε ε ⊂ L ∞ T (H -σ ), with σ = σ(γ) = 3 1/γ -1/2 ,
and then we use (A.3) again. The bounds for the momentum V ε ε can be deduced by a similar argument, after observing that 2γ/(γ + 1) < 2 always.

Next, applying the operator S M to (4.4) immediately gives us system (4.8), where we have set

f ε,M := S M f ε and g ε,M := S M g ε .
Thanks to Lemma 4.1 and (A.2), it is easy to verify inequality (4.9).

We will need also the following important decomposition for the momentum vector fields V ε,M and their curl . Proposition 4.3 For any M > 0 and any ε ∈ ]0, 1], the following decompositions hold true:

V ε,M = ε 2(m-n) t 1 ε,M + t 2 ε,M and curl x V ε,M = ε 2(m-n) T 1 ε,M + T 2 ε,M ,
where, for any T > 0 and s ≥ 0, one has

t 1 ε,M L 2 ([0,T ];H s ) + T 1 ε,M L 2 ([0,T ];H s ) ≤ C(s, M ) t 2 ε,M L 2 ([0,T ];H 1 ) + T 2 ε,M L 2 ([0,T ];L 2 ) ≤ C ,
for suitable positive constants C(s, M ) and C, which are uniform with respect to ε ∈ ]0, 1].

Proof: We decompose V ε,M = ε 2(m-n) t 1 ε,M + t 2 ε,M
, where we define

t 1 ε,M := S M ε -1 ε 2(m-n) u ε and t 2 ε,M := S M (u ε ) . (4.10)
The decomposition of curl x V ε,M follows after setting T j ε,M := curl x t j ε,M , for j = 1, 2. We have to prove uniform bounds for all those terms, by using the estimates established in Section 3.1 above. First of all, we have that u ε ε ⊂ L 2

T (H 1 ), for any T > 0 fixed. Then, we immediately gather the sought bounds for the vector fields t 2 ε,M and T 2 ε,M .

For the families of t 1 ε,M and T 1 ε,M , instead, we have to use the bounds provided by (3.8) and (when γ ≥ 2) Remark 3.1. In turn, we see that for any T > 0:

ε -1 ε 2(m-n) u ε ⊂ L 2 T (L 1 + L 2 + L 6γ/(γ+6) ) → L 2 T (H -σ ) ,
for some σ > 0 large enough. Therefore, the claimed bounds follow thanks to the regularising effect of the operators S M . The proof of the proposition is thus completed.

Convergence of the convective term

In this subsection we show the convergence of the convective term. The first step is to reduce its analysis to the case of smooth vector fields V ε,M .

Lemma 4.4 Let T > 0. For any ψ ∈ C ∞ c [0, T [ ×Ω; R 3 , we have lim M →+∞ lim sup ε→0 + T 0 Ω ε u ε ⊗ u ε : ∇ x ψ dx dt - T 0 Ω V ε,M ⊗ V ε,M : ∇ x ψ dx dt = 0 .
Proof: The proof is very similar to the one of Lemma 4.5 from [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF], for this reason we just outline it. One starts by using the decomposition ε = 1 + ε 2(m-n) R ε to reduce (owing to the uniform bounds of Section 3.1) the convective term to the "homogeneous counterpart": for any test function

ψ ∈ C ∞ c R + × Ω; R 3 , one has lim ε→0 + T 0 Ω ε u ε ⊗ u ε : ∇ x ψ dx dt - T 0 Ω u ε ⊗ u ε : ∇ x ψ dx dt = 0 .
Notice that, here, one has to use that γ ≥ 3/2. After that, we write Notice that, in order to pass to the limit in the weak formulation of the momentum equation and derive the limit system, it is enough to use test functions ψ as above.

u ε = S M (u ε ) + (Id -S M )u ε = t 2 ε,M + (Id -S M )u ε . Using Proposition 4.3 and the fact that (Id -S M ) u ε L 2 T (L 2 ) ≤ C 2 -M ∇ x u ε L 2 T (L 2 ) ≤ C 2 -M ,
Thus, for ψ as in (4.12), we have to pass to the limit in the term

- T 0 Ω V ε,M ⊗ V ε,M : ∇ x ψ = T 0 Ω div (V ε,M ⊗ V ε,M ) • ψ .
Notice that the integration by parts above is well-justified, since all the quantities inside the integrals are smooth with respect to the space variable. Owing to the structure of the test function, and resorting to the notation introduced in (1.5), we remark that we can write

T 0 Ω div (V ε,M ⊗ V ε,M ) • ψ = T 0 R 2 T 1 ε,M + T 2 ε,M • ψ h ,
where we have defined the terms

T 1 ε,M := div h V h ε,M ⊗ V h ε,M and T 2 ε,M := div h V h ε,M ⊗ V h ε,M . (4.13)
In the next two paragraphs, we will deal with each one of those terms separately. We borrow most of the arguments from [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF] (see also [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF], [START_REF] Fanelli | Incompressible and fast rotation limit for barotropic Navier-Stokes equations at large Mach numbers[END_REF] for a similar approach). However, the special structure of the gravity force will play a key role here, in order (loosely speaking) to compensate the stronger singularity due to our scaling 2n > m.

Convergence of the vertical averages

We start by dealing with T 1 ε,M . It is standard to write

T 1 ε,M = div h V h ε,M ⊗ V h ε,M = div h V h ε,M V h ε,M + V h ε,M • ∇ h V h ε,M (4.14) = div h V h ε,M V h ε,M + 1 2 ∇ h V h ε,M 2 + curl h V h ε,M V h ε,M ⊥ .
Notice that the second term is a perfect gradient, so it vanishes when tested against divergence-free test functions. Hence, we can treat it as a remainder, in the sense of (4.11).

For the first term in the second line of (4.14), instead, we take advantage of system (4.8): averaging the first equation with respect to x 3 and multiplying it by V h ε,M , we arrive at

div h V h ε,M V h ε,M = -ε m ∂ t (1) ε,M V h ε,M = ε m (1) ε,M ∂ t V h ε,M + R ε,M ,
in the sense of distributions. We remark that the term presenting the total derivative in time is in fact a remainder, thanks to the factor ε m in front of it. Now, we use the horizontal part of (4.8), where we first take the vertical average and then multiply by

(1)
ε,M : since m > 1, we gather

ε m (1) ε,M ∂ t V h ε,M = - (1) ε,M ∇ h (1) ε,M -ε m-1 (1) ε,M V h ε,M ⊥ + ε m (1) ε,M f h ε,M + ε 2(m-n) (1) ε,M g h ε,M = -ε m-1 (1) ε,M V h ε,M ⊥ + R ε,M ,
where we have repeatedly exploited the properties proved in Proposition 4.2 and we have included in the remainder term also the perfect gradient. Inserting this relation into (4.14) yields

T 1 ε,M = curl h V h ε,M -ε m-1 (1) ε,M V h ε,M ⊥ + R ε,M .
Observe that the first term appearing in the right-hand side of the previous relation is bilinear. Thus, in order to pass to the limit in it, one needs some strong convergence property. As a matter of fact, in the next computations we will work on the regularised wave system (4.8) to show that the quantity

γ ε,M := curl h V h ε,M -ε m-1 (1) ε,M
is compact in some suitable space. In particular, as m > 1, also curl h V h ε,M is compact. In order to see this, we write the vertical average of the first equation in (4.8) as

ε 2m-1 ∂ t (1) ε,M + ε m-1 div h V h ε,M = 0 .
Next, we take the vertical average of the horizontal components of the second equation in (4.8) and then apply curl h : we obtain

ε m ∂ t curl h V h ε,M + ε m-1 div h V h ε,M = ε m curl h f h ε,M + ε 2(m-n) curl h g h ε,M .
At this point, we recall the definition (4.3) of g ε , and we see that curl h g h ε,M ≡ 0. This property is absolutely fundamental, since it allows to erase the last term in the previous relation, which otherwise would have represented an obstacle to get compactness of the γ ε,M 's. Indeed, thanks to this observation, we can sum up the last two equations to get

∂ t γ ε,M = curl h f h ε,M . (4.15)
Using estimate (4.9) in Proposition 4.2, we discover that, for any M > 0 fixed, the family (∂ t γ ε,M ) ε is uniformly bounded (with respect to ε) in e.g. the space L 2 T (L 2 ). On the other hand, we have that, again for any M > 0 fixed, the sequence (γ ε,M ) ε is uniformly bounded (with respect to ε) e.g. in the space L 2

T (H 1 ). Since the embedding H 1 loc → L 2 loc is compact, the Aubin-Lions Lemma implies that, for any M > 0 fixed, the family (γ ε,M ) ε is compact in L 2 T (L 2 loc ). Then, up to extraction of a suitable subsequence (not relabelled here), that family converges strongly to a tempered distribution γ M in the same space. Now, we have that γ ε,M converges strongly to

γ M in L 2 T (L 2 loc ) and V h ε,M converges weakly to V h M in L 2 T (L 2
loc ) (owing to Proposition 4.3, for instance). Then, we deduce that

γ ε,M V h ε,M ⊥ -→ γ M V h M ⊥ in D R + × R 2 .
Observe that, by definition of γ ε,M , we must have γ M = curl h V h M . On the other hand, owing to Proposition 4.3 and (4.10), we know that V h M = S M U h . Therefore, in the end we have proved that, for m > 1 and m + 1 ≥ 2n > m, one has the convergence (at any M ∈ N fixed, when ε → 0

+ ) T 0 R 2 T 1 ε,M • ψ h dx h dt -→ T 0 R 2 curl h S M U h S M (U h ) ⊥ • ψ h dx h dt , (4.16)
for any T > 0 and for any test-function ψ as in (4.12).

Vanishing of the oscillations

We now focus on the term T 2 ε,M , defined in (4.13). Recall that m > 1. In what follows, we consider separately the two cases m + 1 > 2n and m + 1 = 2n. As a matter of fact, in the case when m + 1 = 2n, a bilinear term involving g ε,M has no power of ε in front of it, so it is not clear that it converges to 0 and, in fact, it might persist in the limit, giving rise to an additional term in the target system. To overcome this issue and show that this actually does not happen, we deeply exploit the structure of the wave system to recover a quantitative smallness for that term (namely, in terms of positive powers of ε).

The case m + 1 > 2n

Starting from the definition of T 2 ε,M , the same computations as above yield

T 2 ε,M = div h ( V h ε,M ) V h ε,M + 1 2 ∇ h | V h ε,M | 2 + curl h V h ε,M V h ε,M ⊥ . (4.17)
Let us now introduce the quantities

Φ h ε,M := ( V h ε,M ) ⊥ -∂ -1 3 ∇ ⊥ h V 3 ε,M and ω 3 ε,M := curl h V h ε,M . Then we can write curl V ε,M h = ∂ 3 Φ h ε,M and curl V ε,M 3 = ω 3 ε,M .
In addition, from the momentum equation in (4.8), where we take the mean-free part and then the curl , we deduce the equations

   ε m ∂ t Φ h ε,M -ε m-1 V h ε,M = ε m ∂ -1 3 curl f ε,M h + ε 2(m-n) ∂ -1 3 curl g ε,M h ε m ∂ t ω 3 ε,M + ε m-1 div h V h ε,M = ε m curl h f h ε,M . (4.18)
Making use of the relations above, recalling the definitions in (4.3), and thanks to Propositions 4.2 and 4.3, we can write

curl h V h ε,M V h ε,M ⊥ = ω 3 ε,M V h ε,M ⊥ = ε∂ t Φ h ε,M ⊥ ω 3 ε,M -ε ω 3 ε,M ∂ -1 3 curl f ε,M h ⊥ -ε m+1-2n ω 3 ε,M ∂ -1 3 curl g ε,M h ⊥ = -ε Φ h ε,M ⊥ ∂ t ω 3 ε,M + R ε,M = Φ h ε,M ⊥ div h V h ε,M + R ε,M , (4.19)
where, again, the equalities hold in the sense of distributions. We point out that, thanks to the scaling m + 1 > 2n, we could include in the remainder also the last term appearing in the second equality, which was of order O(ε m+1-2n ).

Hence, putting the gradient term into R ε,M , from (4.17) we arrive at

T 2 ε,M = div h V h ε,M V h ε,M + Φ h ε,M ⊥ + R ε,M = div V ε,M V h ε,M + Φ h ε,M ⊥ -∂ 3 V 3 ε,M V h ε,M + Φ h ε,M ⊥ + R ε,M . (4.20)
At this point, the computations mainly follow the same lines of [START_REF] Feireisl | Multi-scale analysis of compressible viscous and rotating fluids[END_REF] (see also [START_REF] Fanelli | Incompressible and fast rotation limit for barotropic Navier-Stokes equations at large Mach numbers[END_REF], [START_REF] Del Santo | Wróblewska-Kamińska: A multiscale problem for viscous heat-conducting fluids in fast rotation[END_REF]). First of all, we notice that (in the last line) the second term on the right-hand side is another remainder. Indeed, using the definition of the function Φ h ε,M and the fact that the test function ψ does not depend on x 3 , one has

∂ 3 V 3 ε,M V h ε,M + Φ h ε,M ⊥ = ∂ 3 V 3 ε,M V h ε,M + Φ h ε,M ⊥ -V 3 ε,M ∂ 3 V h ε,M + Φ h ε,M ⊥ = R ε,M - 1 2 ∇ h V 3 ε,M 2 = R ε,M .
Next, in order to deal with the first term on the right-hand side of (4.20), we use the first equation in (4.8) to obtain (in the sense of distributions)

div V ε,M V h ε,M + Φ h ε,M ⊥ = -ε m ∂ t (1) ε,M V h ε,M + Φ h ε,M ⊥ + R ε,M = ε m (1) ε,M ∂ t V h ε,M + Φ h ε,M ⊥ + R ε,M .
Now, equations (4.8) and (4.18) immediately yield that

ε m (1) ε,M ∂ t V h ε,M + Φ h ε,M ⊥ = R ε,M - (1)
ε,M ∇ h

ε,M = R ε,M - 1 2 ∇ h (1) ε,M 2 = R ε,M . (1)
This relation finally implies that T 2 ε,M = R ε,M is a remainder, in the sense of relation (4.11): for any T > 0 and any test-function ψ as in (4.12), one has the convergence (at any M ∈ N fixed, when ε → 0

+ ) T 0 R 2 T 2 ε,M • ψ h dx h dt -→ 0 . (4.21)
The case m + 1 = 2n

In the case m + 1 = 2n, most of the previous computations may be reproduced exactly in the same way. The only (fundamental) change concerns relation (4.19): since now m + 1 -2n = 0, that equation now reads

curl h V h ε,M V h ε,M ⊥ = Φ h ε,M ⊥ div h V h ε,M -ω 3 ε,M ∂ -1 3 curl g ε,M h ⊥ + R ε,M , (4.22)
and, repeating the same computations performed for T 2 ε,M in the previous paragraph, we have

T 2 ε,M = R ε,M -ω 3 ε,M ∂ -1 3 curl g ε,M h ⊥ .
Hence, the main difference with respect to the previous case is that, now, we have to take care of the

term ω 3 ε,M ∂ -1 3 curl g ε,M h ⊥
, which is non-linear and of order O(1), so it may potentially give rise to oscillations which persist in the limit.

In order to show that this does not happen, we make use of definition (4.3) of g ε to compute

curl g ε,M h,⊥ = curl (1) ε,M ∇ x G -∇ x Π ε,M h,⊥ =    -∂ 2 (1) ε,M ∂ 1 (1) ε,M 0    h,⊥ = -∇ h (1)
ε,M .

From this relation, in turn we get

T 2 ε,M = R ε,M + ω 3 ε,M ∂ -1 3 ∇ h (1) ε,M . (4.23)
Now, we have to employ the potential part of the momentum equation in (4.8), which has not been used so far. Taking the oscillating component of the solutions, we obtain

∇ h (1) ε
,M = -ε m ∂ t V h ε,M -ε m-1 ( V h ε,M ) ⊥ + ε m f h ε,M + ε 2(m-n) g h ε,M = -ε m ∂ t V h ε,M + R ε,M .
Inserting this relation into (4.23) and using (4.18), we finally gather

T 2 ε,M = -ε m ω 3 ε,M ∂ t ∂ -1 3 V h ε,M + R ε,M = ε m ∂ t ω 3 ε,M ∂ -1 3 V h ε,M + R ε,M = R ε,M ,
because we have taken m > 1. This relation finally implies that, also in the case when m + 1 = 2n, T 2 ε,M is a remainder: for any T > 0 and any test-function ψ as in (4.12), one has the convergence (4.21).

The limit system

Thanks to the computations of the previous subsection, we can now pass to the limit in equation (2.16). Recall that m > 1 and m + 1 ≥ 2n > m here.

To begin with, we take a divergence-free test function ψ as in (4.12), specifically

ψ = ∇ ⊥ h φ, 0 , with φ ∈ C ∞ c [0, T [ ×R 2 , φ = φ(t, x h ) . (4.24)
We point out that, since all the integrals will be made on R 2 (in view of the choice of the test functions in (4.24) above), we can safely work on the domain Ω = R 2 × ]0, 1[ . In addition, for ψ as in (4.24), all the gradient terms vanish identically, as well as all the contributions due to the vertical component of the equation. In particular, we do not see any contribution of the pressure and gravity terms: equation (2.16) becomes Making use of the uniform bounds of Section 3.1, we can pass to the limit in the ∂ t term and in the viscosity term. Moreover, our assumptions imply that 0,ε u 0,ε u 0 in e.g. L 2 loc . Next, the Coriolis term can be dealt with in a standard way: using the structure of ψ and the mass equation (2.15), we can write

T 0 Ω -ε u h ε • ∂ t ψ h -ε u h ε ⊗ u h ε : ∇ h ψ h + 1 ε ε u h ε ⊥ • ψ h dx dt
T 0 Ω 1 ε ε u h ε ⊥ • ψ h = T 0 R 2 1 ε ε u h ε • ∇ h φ = -ε m-1 T 0 R 2 (1) ε ∂ t φ -ε m-1 R 2
(1) 0,ε φ(0, •) , which of course converges to 0 when ε → 0 + . It remains us to tackle the convective term ε u h ε ⊗ u h ε . For it, we take advantage of Lemma 4.4 and relations (4.16) and (4.21), but we still have to take care of the convergence for M → +∞ in (4.16). We start by performing equalities (4.14) backwards in the term on the right-hand side of (4.16): thus, we have to pass to the limit for M → +∞ in

T 0 R 2 U h M ⊗ U h M : ∇ h ψ h dx h dt .
Now, we remark that, since U h ∈ L 2 T (H 1 ) by (3.7), from (A.2) we gather the strong convergence

S M U h -→ U h in L 2
T (H s ) for any s < 1, in the limit for M → +∞. Then, passing to the limit for M → +∞ in the previous relation is an easy task: we finally get that, for ε → 0 + , one has

T 0 Ω ε u h ε ⊗ u h ε : ∇ h ψ h -→ T 0 R 2 U h ⊗ U h : ∇ h ψ h .
In the end, we have shown that, letting ε → 0 + in (4.25), one obtains

T 0 R 2 U h • ∂ t ψ h + U h ⊗ U h : ∇ h ψ h dx h dt = T 0 R 2 µ∇ h U h : ∇ h ψ h dx h dt - R 2 u h 0 • ψ h (0, •) dx h ,
for any test function ψ as in (4.12). This implies (2.19), concluding the proof of Theorem 2.5.

5 Proof of the convergence for m = 1

In the present section, we complete the proof of the convergence in the case m = 1 and 1/2 < n < 1. We will use again the compensated compactness argument depicted in Section 4.2, and in fact most of the computations apply also in this case.

Analysis of the acoustic-Poincaré waves

When m = 1, the wave system (4.4) takes the form

ε ∂ t (1) ε + div V ε = 0 ε ∂ t V ε + ∇ x (1) ε + e 3 × V ε = ε f ε + ε 2(1-n) g ε , (5.1)
where

(1) ε ε and V ε ε are defined as in Section 4.1.1. This system is supplemented with the initial datum (1) 0,ε , 0,ε u 0,ε .

Next, we regularise all the quantities, by applying the Littlewood-Paley cut-off operator S M to (5.1): we deduce that [START_REF] Bahouri | Fourier Analysis and Nonlinear Partial Differential Equations[END_REF] ε,M and V ε,M , defined as in (4.6), satisfy the regularised wave system

ε ∂ t (1) ε,M + div V ε,M = 0 ε ∂ t V ε,M + ∇ x (1)
ε,M + e 3 × V ε,M = ε f ε,M + ε 2(1-n) g ε,M , (5.2) in the domain R + × Ω, where we recall that f ε,M := S M f ε and g ε,M := S M g ε . It goes without saying that a result similar to Proposition 4.2 holds true also in this case.

As it is apparent from the wave system (5.1) and its regularised version, when m = 1 the pressure term and the Coriolis term are in balance, since they are of the same order. This represents the main change with respect to the case m > 1, and it comes into play in the compensated compactness argument. Therefore, despite most of the computations may be repeated identical as in the previous section, let us present the main points of the argument.

Handling the convective term when m = 1

Let us take care of the convergence of the convective term in the case when m = 1.

First of all, it is easy to see that Lemma 4.4 still holds true. Therefore, given a test function ψ ∈ C ∞ c [0, T [ ×Ω; R 3 such that div ψ = 0 and ∂ 3 ψ = 0, we have to pass to the limit in the term

- T 0 Ω V ε,M ⊗ V ε,M : ∇ x ψ = T 0 Ω div (V ε,M ⊗ V ε,M ) • ψ = T 0 R 2 T 1 ε,M + T 2 ε,M • ψ h ,
where we agree again that the torus T 1 has been normalised so that its Lebesgue measure is equal to 1 and we have adopted the same notation as in (4.13).

At this point, we notice that the analysis of T 2 ε,M can be performed as in Section 4.2.2, because we have m + 1 > 2n, i.e. n < 1. Mutatis mutandis, we find relation (4.21) also in the case m = 1.

Let us now deal with the term T 1 ε,M . Arguing as in Section 4.2.1, we may write it as

T 1 ε,M = curl h V h ε,M - (1) ε,M V h ε,M ⊥ + R ε,M .
Now we use the horizontal part of (5.2): averaging it with respect to the vertical variable and applying the operator curl h , we find

ε ∂ t curl h V h ε,M + div h V h ε,M = ε curl h f h ε,M .
Taking the difference of this equation with the first one in (5.2), we discover that

∂ t γ ε,M = curl h f h ε,M , where γ ε,M := curl h V h ε,M - (1)
ε,M .

An argument analogous to the one used after (4.15) above, based on Aubin-Lions Lemma, shows that γ ε,M ε is compact in e.g. L 2 T (L 2 loc ). Then, this sequence converges strongly (up to extraction of a suitable subsequence, not relabelled here) to a tempered distribution γ M in the same space.

Using the previous property, we may deduce that

γ ε,M V h ε,M ⊥ -→ γ M V h M ⊥ in D R + × R 2 ,

  3): we assume that either m > 1 and m < 2 n ≤ m + 1 , or m = 1 and 1 2 < n < 1 . (1.4)

  Some notation and conventions. Let B ⊂ R n . The symbol 1 B denotes the characteristic function of B. The notation C ∞ c (B) stands for the space of C ∞ functions on R n and having compact support in B. The dual space D (B) is the space of distributions on B.

  which holds in view of estimate (A.3) from the Appendix and the uniform bound (3.6), one can conclude. From now on, for notational convenience, we generically denote by R ε,M any remainder term, that is any term satisfying the property lim M • ψ dx dt = 0 (4.11) for all test functions ψ ∈ C ∞ c [0, T [ ×Ω; R 3 lying in the kernel of the singular perturbation operator, namely (in view of Proposition 3.4) such that ψ ∈ C ∞ c [0, T [ ×Ω; R 3 such that div ψ = 0 and ∂ 3 ψ = 0 . (4.12)

  x u ε ) : ∇ x ψ dx dt + Ω 0,ε u 0,ε • ψ(0, •) dx .

In this section, we introduce the primitive system and formulate our working hypotheses (see Section 2.1), then we state our main results (in Section 2.2).

We agree that f (D) stands for the pseudo-differential operator u → F -1 [f (ξ) u(ξ)].

Acknowledgements

The work of the second and third authors has been partially supported by the project CRISIS (ANR-20-CE40-0020-01), operated by the French National Research Agency (ANR). The last author is supported by (Polish) National Center of Science grant 2020/38/E/ST1/00469.

The first and the third authors are members of the Italian Institute for Advanced Mathematics (IN-dAM) group.

where we have

M . Owing to the regularity of the target velocity U h , we can pass to the limit also for M → +∞, as detailed in Section 4.3 above. Thus, we find

for all test functions ψ such that div ψ = 0 and ∂ 3 ψ = 0. Recall the convention 1) when m = 1 (keep in mind Proposition 3.5), the last term in the integral on the right-hand side is actually zero.

End of the proof

Thanks to the previous analysis, we are now ready to pass to the limit in equation (2.16). For this, we take a test-function ψ as in (4.24); notice in particular that div ψ = 0 and ∂ 3 ψ = 0. Then, once again all the gradient terms and all the contributions coming from the vertical component of the momentum equation vanish identically, when tested against such a ψ. Recall that all the integrals will be performed in R 2 . So, equation (2.16) reduces to

For the rotation term, we can test the first equation in (5.1) against φ to get

whence we deduce that

0,ε φ(0, •) .

In addition, the convergence of the convective term has been performed in (5.3). As for other terms, we can argue as in Section 4.3. Hence, letting ε → 0 + in the equation above, we get

which is the weak formulation of equation (2.20). In the end, also Theorem 2.6 is proved.

A Appendix -A few tools from Littlewood-Paley theory

Let us present some tools from Littlewood-Paley theory, which we have exploited in our analysis. We refer e.g. to Chapter 2 of [START_REF] Bahouri | Fourier Analysis and Nonlinear Partial Differential Equations[END_REF] for details. For simplicity of exposition, we deal with the R d case, with d ≥ 1;

however, the whole construction can be adapted also to the d-dimensional torus T d , and to the "hybrid" case R d1 × T d2 . First of all, we introduce the so-called Littlewood-Paley decomposition. We fix a smooth radial function χ such that Supp χ ⊂ B(0, 2), χ ≡ 1 in a neighborhood of B(0, 1) and the map r → χ(r e) is non-increasing over R + for all unitary vectors e ∈ R d . Set ϕ (ξ) = χ (ξ) -χ (2ξ) and ϕ j (ξ) := ϕ(2 -j ξ) for all j ≥ 0. The dyadic blocks (∆ j ) j∈Z are defined by 1

For any j ≥ 0 fixed, we also introduce the low frequency cut-off operator

Note that S j is a convolution operator. More precisely, after defining

for all j ∈ N and all tempered distributions u ∈ S we have that S j u = K j * u. Thus the L 1 norm of K j is independent of j ≥ 0. This implies that S j maps continuously L p into itself, for any 1 ≤ p ≤ +∞. Moreover, the following property holds true: for any u ∈ S , then one has the equality u = j ∆ j u in the sense of S . Let us also recall the so-called Bernstein inequalities.

Lemma A.1 Let 0 < r < R. A constant C exists so that, for any non-negative integer k, any couple (p, q) in [1, +∞] 2 , with p ≤ q, and any function u ∈ L p , we have, for all λ > 0,

By use of Littlewood-Paley decomposition, we can define the class of Besov spaces.

Definition A.2 Let s ∈ R and 1 ≤ p, r ≤ +∞. The non-homogeneous Besov space B s p,r is defined as the subset of tempered distributions u for which u B s p,r := 2 js ∆ j u L p j≥-1 r < +∞ . Besov spaces are interpolation spaces between Sobolev spaces. In fact, for any k ∈ N and p ∈ [1, +∞] we have the chain of continuous embeddings B k p,1 → W k,p → B k p,∞ , which, in the case when 1 < p < +∞, can be refined to B k p,min(p,2) → W k,p → B k p,max(p,2) . In particular, for all s ∈ R we deduce that B s 2,2 ≡ H s , with equivalence of norms:

Observe that, from that equivalence, we easily get the following property: for any f ∈ H s and any j ∈ N, one has Id -S j f H σ ≤ C ∇f H s-1 2 -j(s-σ) for all σ ≤ s , (A.3)

where C > 0 is a "universal" constant, independent of f , j, s and σ. This inequality has been repeatedly used in our computations.