arxiv: v1 [math.ap] 13 Dec 2012

Transcription

1 THE D INCOMPRESSIBLE BOUSSINESQ EQUATIONS WITH GENERAL CRITICAL DISSIPATION QUANSEN JIU 1, CHANGXING MIAO, JIAHONG WU 3 AND ZHIFEI ZHANG 4 arxiv:11.37v1 [math.ap] 13 Dec 01 Abstract. This paper aims at the global regularity problem concerning the D incompressible Boussinesq equations with general critical dissipation. The critical dissipation refers to α+β = 1 when Λ α ( ) α and Λ β represent the fractional Laplacian dissipation in the velocity and the temperature equations, respectively. We establish the global regularity for the general case with α +β = 1 and α 0 < α < 1. The cases when α = 1 and when α = 0 were previously resolved by Hmidi, Keraani and Rousset [7, 8]. The global existence and uniqueness is achieved here by exploiting the global regularity of a generalized critical surface quasi-gesotrophic equation as well as the regularity of a combined quantity of the vorticity and the temperature. 1. Introduction This paper studies the global(in time) regularity of solutions to the D incompressible Boussinesq equations with a general critical dissipation t u+u u+νλ α u = p+θe, x R, t > 0, u = 0, x R, t > 0, (1.1) t θ +u θ +κλ β θ = 0, x R, t > 0, u(x,0) = u 0 (x), θ(x,0) = θ 0 (x), x R, where u = u(x, t) denotes the D velocity, p = p(x, t) the pressure, θ = θ(x, t) the temperature, e the unit vector in the vertical direction, and ν > 0, κ > 0, 0 < α < 1 and 0 < β < 1 are real parameters. Here Λ = represents the Zygmund operator with Λ α being defined through the Fourier transform, namely Λ α f(ξ) = ξ α f(ξ), where the Fourier transform is given by f(ξ) = e ix ξ f(x)dx. R (1.1) generalizes the standard D Boussinesq equations in which the dissipation is given by the Laplacian operator. The D Boussinesq equations and their fractional Laplacian generalizations have attracted considerable attention recently due to their physical applications and mathematical significance. The Boussinesq equations model geophysical flows such as atmospheric fronts and oceanic circulation, and play an important role in the study of Raleigh-Bernard convection (see, e.g., [13, 3, 35, 41]). Mathematically the Date: October 18, Mathematics Subject Classification. 35Q35, 35B35, 35B65, 76D03. Key words and phrases. Critical Boussinesq equations, global regularity. 1

2 Q. JIU, C. MIAO, J. WU AND Z. ZHANG D Boussinesq equations serve as a lower dimensional model of the 3D hydrodynamics equations. In fact, the Boussinesq equations retain some key features of the 3D Navier- Stokes and the Euler equations such as the vortex stretching mechanism. As pointed out in [36], the inviscid Boussinesq equations can be identified with the 3D Euler equations for axisymmetric flows. One main focus of recent research on the D Boussinesq equations has been on the global regularity issue when the dissipation is given by a fractional Laplacian or is present only in one direction (see, e.g., [1,, 3, 6, 7, 9, 10, 16, 19, 0, 1,, 4, 5, 6, 7, 8, 9, 30, 31, 33, 34, 37, 39, 40, 46, 47]). Most recent work targets the critical and the supercritical cases. The critical case refers to α+β = 1 in (1.1) while the supercritical case corresponds to α + β < 1. In [7, 8] Hmidi, Keraani and Rousset were able to establish the global regularity for two critical cases: (1.1) with α = 1 and κ = 0 and (1.1) with ν = 0 and β = 1. Miao and Xue in [37] obtained the global regularity for (1.1) with ν > 0, κ > 0 and α ( 6 6, 1 4 ), β ( 1 α, min { 7+ }) 6 α, α(1 α), α. 5 6 α In addition, Constantin and Vicol [16] verified the global regularity of (1.1) with ν > 0, κ > 0, α (0,), β (0,), β > +α. We briefly mention that the logarithmically supercritical case has also been dealt with and global regularity has been established ([10, 4, 30, 31]). This paper aims at the global regularity of (1.1) with a general critical dissipation 0 < α < 1, 0 < β < 1, α+β = 1 and we succeed in the case when α 0 < α < 1, where α 0 is given in (1.) below. More precisely we are able to prove the following theorem. Theorem 1.1. Let α 0 < α < 1 and α+β = 1, where α 0 = (1.) Assume that u 0 B σ,1 (R ) with σ 5 and θ 0 B,1 (R ). Then (1.1) has a unique global solution (u,θ) satisfying, for any 0 < T <, u C([0,T];B σ,1 (R )) L 1 ([0,T];B σ+α,1 (R )), θ C([0,T];B,1(R )) L 1 ([0,T];B +β,1 (R )). (1.3) Here B s q,r with s R and q,r [1, ] denotes an inhomogeneous Besov space and its precise definition is given in Section 4. The key component in the proof of Theorem 1.1 is to establish the global a priori bounds in the class defined in (1.3). This does not appear to be trivial and the energy methods are not sufficient for this purpose. Although the global bounds for u in L ([0,T];L ) and θ in L ([0,T];L q ) with q [, ] can be easily obtained, the global bounds for the derivatives are not evident. To avoid the

3 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 3 pressure term, we resort to the vorticity formulation { t ω +u ω +νλ α ω = 1 θ, u = ψ, ψ = ω or u = 1 ω. (1.4) However, the vortex stretching term 1 θ appears toprevent usfromproving anyglobal bound for ω. A natural idea would be to eliminate 1 θ from the vorticity equation. For notational convenience, we set ν = κ = 1 in (1.1) throughout the rest of this paper. Realizing that Λ α ω 1 θ = Λ α (ω Λ α 1 θ), we can hide 1 θ by considering the new quantity G = ω R α θ with R α = Λ α 1, which satisfies t G+u G+Λ α G = [R α,u ]θ+λ β α 1 θ. (1.5) Here we have used the standard commutator notation [R α,u ]θ = R α (u θ) u R α θ. (1.5) can be obtained by taking the difference of the vorticity equation and the resulting equation after applying R α to the temperature equation. The quantity G was first introduced in [7] to deal with the critical case when α = 1 and β = 0. Although (1.5) appears to be more complicated than the vorticity equation, but the commutator term [R α,u ]θ is less singular than 1 θ in the vorticity equation. By obtaining a suitable bound for [R α,u ]θ, we are able to obtain a global bound for G L when α > 4. In addition, by fully exploiting the dissipation with α > 4, a global bound is 5 5 also established for G L q when q is in the range < q < q 0 8 4α 8 7α. (1.6) This global bound for G L q enables us to gain further regularity for G. In fact, we establish that, for α 0 < α and s 3α, the Besov norm G B s q, obeys, for any T > 0 and t T, G(t) B s q, C, (1.7) where C is a constant depending on T and the norms of the initial data. In contrast to the critical case with α = 1 and κ = 0 dealt with in [7], the general critical case appears to be more difficult. The regularity of G here does not translate to the regularity on the vorticity ω since the corresponding regularity of R α θ is not known for α < 1. This paper offers a different approach by gaining further regularity through the temperature equation t θ+u θ+λ β θ = 0. (1.8) Since u is determined by ω through the Biot-Savart law, or u = 1 ω and ω = G+R α θ, we can decompose u into two parts, u = 1 ω = 1 G+ 1 R α θ ũ+v. (1.9) For α 0 < α and as a consequence of (1.7), ũ is regular in the sense that ũ L = 1 G L C G B s q, C.

4 4 Q. JIU, C. MIAO, J. WU AND Z. ZHANG In addition, when α+β = 1, v in terms of θ can be written as v = 1 Λ (1 β) 1 θ. Therefore, (1.8) is almost a generalized critical surface quasi-geostrophic (SQG) type equation first studied in [14] except that u here contains an extra regular velocity ũ. We remark that there is a large literature on the SQG equation and interested readers may consult [8, 11, 1, 14, 15, 16, 17, 18, 3] and the references therein. Since energy estimates do not appear to yield the desired global a priori bounds, we employ the approach of Constantin and Vicol [16] to establish the global regularity of (1.8) and (1.9). Different from the Schwartz class setting in [16], the initial data θ 0 here is in H 1 with θ 0 L <. The precise global existence and uniqueness of (1.8) and (1.9) obtained here can be stated as follows. Theorem 1.. Let β (0,1) and 0 < T <. Let ũ be a D vector field satisfying ũ = 0 and M max { ũ L (0,T;L (R )), ũ L (0,T;L (R ))} <. Consider the generalized critical SQG type equation t θ+u θ+λ β θ = 0, x R, t > 0, u = ũ+v, v = Λ 3+β 1 θ, x R, t > 0, θ(x,0) = θ 0 (x), x R. Assume that θ 0 H 1 (R ) with θ 0 L <. Then (1.10) has a unique solution θ C([0,T];H 1 (R )) satisfying θ L (0,T;L (R )) C(M, θ 0 H 1, θ 0 L,T). (1.10) In order to prove this theorem, we need to convert the operator relating v and 1 θ, namely Λ 3+β into an integral form. Since β (0,1), the standard Riesz potential formula does not appear to apply here (see, e.g., [43]). Nevertheless, Λ 3+β can be represented through an integral kernel by making use of the inverse Fourier transform of functions of the form P k(ξ), where P ξ k+ β k is a harmonic polynomial of degree k ([43, p.73]). As a special consequence, the symmetric part of v can be represented as S( v) 1 ( v +( v)t ) = C σ(x y) R x y ( 1θ(x) 1+β 1 θ(y))dy with σ(z) = 1 [ ] z1 z z1 z z z1 z. z 1 z More details can be found in Section. The gained regularity in θ via Theorem 1. allows us to assert the desired regularity in the velocity and the vorticity. Especially, ω L is bounded on any time interval [0, T]. Further regularity leading to (1.3) is established through energy estimates in Besov space settings. With these global bounds at our disposal, the global existence part of Theorem 1.1 follows from a local existence through a standard procedure such as the successive approximation and an extension of the local solution into a global one

5 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 5 with the aid of the global a priori bounds. The uniqueness of solutions in the class (1.3) is clear. The rest of this paper is organized as follows. The second section represents the relation v = 1 Λ (1 β) 1 θ as an integral. Integral formulas for v and its symmetric part are also derived in this section. The third section proves Theorem 1.. The fourth section provides the definitions of functional spaces such as the Besov spaces and related facts. In addition, a commutator estimate in the Besov space setting is also proven in this section. This commutator estimate will be used extensively in the sections that follow. Sections 5 and 6 establish global a priori bounds for G L and for G L q, where q satisfies (1.6). The last section proves Theorem 1.1. To do so, we first obtain a global bound of G in the Besov space Bq, s with any s 3α. This global bound and Theorem 1. yield the proof of Theorem Representing v = Λ 3 β 1 θ and v as integrals In this section we represent v = Λ 3 β 1 θ, its gradient v and the symmetric part of v as integrals. These integral representations will be used in the next section. More precisely, we prove the following lemma. Lemma.1. Let θ be a smooth function of R which is sufficiently rapidly decreasing at. Let β (0,1) and v be given by Then v and v can be written as v = Λ 3 β 1 θ. (.1) R (x y) v(x) = C(β) x y 1θ(y)dy, [ ] 1+β v(x) = C(β) 1 0 R x y 1θ(y)dy 1+β (x y) (x y) (1+β)C(β) R x y 3+β 1 θ(y)dy, (.) where C(β) is a constant depending on β only, and a b denotes the tensor product of two vectors a and b, namely a b = (a i b j ). Especially the symmetric part of v, denoted by S( v), is given by S( v) 1 ( v+( v)t ) = C where C is constant depending on β only and σ(z) = 1 z σ(x y) R x y ( 1θ(x) 1+β 1 θ(y))dy, (.3) [ ] z1 z z1 z z1. z z 1 z To derive the formulas in Lemma.1, we need a fact stated in the following lemma, which can be found in Stein s book ([43, p.73]). We note that the Fourier transform

6 6 Q. JIU, C. MIAO, J. WU AND Z. ZHANG f(ξ) in Stein s book ([43, p.46]) is defined by f(ξ) = e πix ξ f(x)dx, R d which is normally the definition of the inverse Fourier transform. This explains why the Fourier transform f in Stein s book is here changed to the inverse Fourier transform f. Lemma.. Let P k (x) with x R d be a homogeneous harmonic polynomial of degree k, where k 1 is an integer. Let 0 < α < d. Then, for a constant C(d,k,α) depending on d,k and α only ( ) Pk (ξ) = C(d,k,α) P k(x) ξ d+k α x k+α in the sense that R d P k (ξ) P k (x) φ(ξ)dξ = C(d,k,α) ξ d+k α R x d k+α φ (x)dx for every φ which is sufficiently rapidly decreasing at, and whose inverse Fourier transform has the same property. With this lemma at our disposal, we can now prove Lemma.1. Proof of Lemma.1. By (.1), According to Lemma., v(ξ) = iξ ξ 3+β 1 θ(ξ). v(x) = i ( ξ ξ 3+β) 1 θ = C(β) x 1θ = C(β) 1+β x R (x y) x y 1+β 1θ(y)dy. Since β (0,1), the kernel in the representation of v is not singular and we have [ ] (x y) v = C(β) x R x y 1+β 1 θ(y)dy ( ) = C(β) 1 0 R x y 1θ(y)dy 1+β (1+β)C(β) (x y) (x y) R x y 3+β 1 θ(y)dy. (.3) is obtained by taking the symmetric part of (.) and inserting 1 θ(x) in the resulting integral. The inserted term does not contribute to the integral thanks to the fact that, for any r > 0, σ(z)dz = 0. z =r This completes the proof of Lemma.1.

7 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 7 3. Global regularity for an active scalar with critical dissipation The goal of this section is to prove Theorem 1., which states the global regularity of a generalized SQG type equation with critical dissipation. The proof is obtained by modifying the approach of Constantin and Vicol [16]. Different from [16], the functional setting here is weaker. To prove Theorem 1., we first establish the global existence and uniqueness of (1.10) when ũ is smooth and θ 0 is smooth and decays sufficiently fast at infinity. Theorem 3.1. Let β (0,1] and 0 < T <. Let θ 0 C (R ) and decays sufficiently fast at. Let ũ be a smooth D vector field satisfying ũ = 0 and M max { ũ L (0,T;L (R )), ũ L (0,T;L (R ))} <. Then (1.10) has a unique global smooth solution θ on [0,T]. In addition, θ L (0,T;L (R )) C(M, θ 0 L, θ 0 L,T), (3.1) where the constant C in the inequality above depends only on the quantities inside the parenthesis. We will be more specific about the decay requirement on θ 0 in the proof of Theorem 3.1. The rest of this section proves Theorem 3.1 and then Theorem 1.. The proof of Theorem 3.1 employs the method of Constantin and Vicol [16]. We recall a basic concept. For a given δ > 0, a function f is said to have only small shocks with a parameter δ or simply f OSS δ if there is L > 0 such that f(x) f(y) δ whenever x y < L. (3.) TheproofofTheorem 3.1consists oftwo mainparts. The first partshows thatθ OSS δ for a suitable δ = δ( θ 0 L ) implies (3.1). L in (3.) is not required to be big in order for this part of the result to hold. What we really need here is a smoothness property on θ and it suffices for θ OSS δ with a small L, say L < 1. The second part proves that θ 0 OSSδ for some L > 0 implies that θ OSS δ for the same L. For the sake of 4 clarity, we present each part as a proposition. Proposition 3.. Let β (0,1] and 0 < T <. Let ũ be a vector field satisfying ũ = 0 and M max { ũ L (0,T;L (R )), ũ L (0,T;L (R ))} <. If θ is in OSS δ uniformly on [0,T] with δ given by C δ = β (3.3) β θ 0 L for a suitable pure constant C > 0 (independent of θ 0 ), then θ is bounded as in (3.1). Proposition 3.3. Let β (0,1] and 0 < T <. Let ũ be a vector field satisfying ũ = 0 and M max { ũ L (0,T;L (R )), ũ L (0,T;L (R ))} <. Assume that θ 0 OSS δ/4 for some δ > 0. Then θ OSS δ uniformly on [0,T].

8 8 Q. JIU, C. MIAO, J. WU AND Z. ZHANG With the two propositions above in our disposal, we can now prove Theorem 3.1. Proof of Theorem 3.1. Since θ 0 C (R ) and decays sufficiently fast at, it is easy to check θ 0 OSS δ/4 for δ given by (3.3). There are two different ways to achieve this. The first is to use the simple inequality θ 0 (x) θ 0 (y) θ 0 L x y and take L = δ/(4 θ 0 L ). Alternatively, we first take R > 0 such that θ 0 (x) < δ/8 for any x R. Then, by the uniform continuity of θ 0 in the disk x R, there is L 1 > 0 such that θ 0 (x) θ 0 (y) δ/4 when x y L 1. Taking L = min{r,l 1 }, we obtain θ 0 OSS δ/4. By Proposition 3.3, θ OSS δ uniformly on [0,T]. By Proposition 3., θ satisfies (3.1). The global existence and uniqueness follows from a local well-posedness and an extension to a global solution through the global bound for θ L in (3.1). We omit further details. This completes the proof of Theorem 3.1. Proof of Theorem 1.. We first regularize ũ and the initial data. For ǫ > 0, we define ρ ǫ to be the standard mollifier, namely x ) ρ(x) = ρ( x ) C0 (R ), ρ 0, ρ(x)dx = 1, ρ ǫ (x) = ǫ ρ(. R ǫ In addition, let χ ǫ be the standard smooth cutoff, namely χ ǫ (x) = χ(ǫx),χ C0 (R ) and χ(x) = 1 in x. Now we define ũ ǫ = ρ ǫ ũ, θ ǫ 0 = ρ ǫ (χ ǫ θ 0 ). and consider the following regularized initial-value problem t θ ǫ +u ǫ θ ǫ +Λ β θ ǫ = 0, x R, t > 0, u ǫ = ũ ǫ +v ǫ, v ǫ = Λ 3+β 1 θ ǫ, x R, t > 0, θ ǫ (x,0) = θ0(x), ǫ x R. (3.4) Then θ ǫ 0 C (R ) and has the decay properties required in Theorem 3.1. Therefore, by Theorem 3.1, (3.4) has a unique global smooth solution θ ǫ satisfying θ ǫ L (0,T;L (R )) C(M, θ ǫ 0 L, θǫ 0 L,T). Since θ0 ǫ L θ 0 L 3 6 θ L θ L θ L and θ0 ǫ L θ 0 L, θ ǫ admits a global bound that is uniform with respect to ǫ, θ ǫ L (0,T;L (R )) C(M, θ 0 H 1, θ 0 L,T). In addition, a simple energy estimate shows that In fact, (3.5) follows from the energy inequality θ ǫ H 1 C(M, θ 0 H 1, θ 0 L,T). (3.5) 1 d dt θǫ β L + Λ1+ θ ǫ L ( ũǫ L + v ǫ L ) θ ǫ L,

9 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 9 and the global bounds ũ ǫ L ũ L and v ǫ L = Λ 3+β 1 θ ǫ L C( θ ǫ L + θ ǫ L ) (3.6) C(M, θ 0 H 1, θ 0 L,T). The interpolation inequality in (3.6) can be proven through the Littlewood-Paley decomposition and the details are deferred to the end of Section 4. Then θ ǫ has a weak limit θ H 1, namely θ ǫ θ weakly in H 1 as ǫ 0. (3.7) In addition, it can be shown that θ ǫ θ strongly in L as ǫ 0. (3.8) This can be achieved by estimating the difference θ ǫ 1 θ ǫ through energy estimates. Let u ǫ 1 and u ǫ be the corresponding velocities. Then θ = θ ǫ 1 θ ǫ and u = u ǫ 1 u ǫ satisfy Therefore, Noticing that t θ +u θ ǫ 1 +u ǫ θ +Λ β θ = 0. d dt θ L + Λβ/ θ L = u θ ǫ 1 θdx. R u = u ǫ 1 u ǫ = ũ ǫ 1 ũ ǫ +v ǫ 1 v ǫ, we obtain, by Hölder s inequality, u θ ǫ 1 θdx ũ ǫ 1 ũ ǫ L θ ǫ 1 L θ L R + v ǫ 1 v ǫ L β θ ǫ 1 L 1 β θ L. Due to ũ ǫ 1 ũ ǫ L C ǫ 1 ǫ M and recalling (3.5), we have ũ ǫ 1 ũ ǫ L θ ǫ 1 L θ L C ǫ 1 ǫ θ L. where C = C(M, θ 0 H 1, θ 0 L,T). Since v ǫ 1 v ǫ = 1 Λ (1 β) 1 θ, and by Hardy-Littlewood-Sobolev inequality, By a simple interpolation inequality, Therefore, v ǫ 1 v ǫ L β C Λ (1 β) θ L β C θ L. θ ǫ 1 L 1 β θǫ 1 1 β L θ ǫ 1 β L C(M, θ 0 H 1, θ 0 L,T). d dt θ L + Λβ/ θ L C ǫ 1 ǫ θ L + C θ L.

10 10 Q. JIU, C. MIAO, J. WU AND Z. ZHANG Consequently θ(t) L θ ǫ 1 (t) θ ǫ (t) L C( ǫ 1 ǫ + θ ǫ 1 0 θǫ 0 L )ect C( ǫ 1 ǫ + θ ǫ 1 0 θ 0 L + θ 0 θ ǫ 0 L )e Ct C max{ǫ 1,ǫ }e Ct, where C = C(M, θ 0 H 1, θ 0 L,T). This proves (3.8). By the interpolation inequality, for 0 < s < 1, θ ǫ θ H s C θǫ θ 1 s /s L θ ǫ θ s /s H s, (3.7) and (3.8) imply that, for 0 < s < 1, θ ǫ θ strongly in H s as ǫ 0. This strong convergence allows us to show that θ is the corresponding global solution of (1.10) satisfying θ L ([0,T];H 1 ) and θ L C(M, θ 0 H 1, θ 0 L,T). The continuity of θ in time, namely θ C([0,T];H 1 ), follows from a standard approach (see, e.g., [36, p.111] or [4, p.138]) and we omit the details. Finally, it is not hard to verify that such solutions are unique. This completes the proof of Theorem 1.. We now turn to the proofs of Propositions 3. and 3.3. To prove Proposition 3., we will make use of the following lower bound obtained in [16]. This lower bound improves a pointwise inequality of Córdoba and Córdoba [17]. Lemma 3.4. Let β (0,) and q [1, ]. Let f C 1 (R d ) decay sufficiently fast at. Then the pointwise lower bound holds, f Λ β ( f) 1 Λβ( f ) + 1 f + D( f)+ C 0 f βq q+d βq, q+d L q where C 0 = C 0 (d,β,q) is a constant and D is given by the principal value integral D(g) = C(d,β)P.V. Proof of Proposition 3.. We show that, if Rd g(x) g(y) x y d+β dy 0. θ OSS δ with δ = C β β θ 0 L for a suitable C, then θ satisfies (3.1). Clearly, θ satisfies t ( θ)+u ( θ)+λ β ( θ) = u( θ). Dotting with θ and applying Lemma 3.4 yield 1 ( t +u +Λ β ) θ + 1 θ +β D( θ)+ C 0 θ 0 β L θ S( u) θ S( u) L θ(x), (3.9)

11 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 11 where S( u) denotes the symmetric part of u, or S( u) = 1 (( u)+( u)t ). Clearly S( u) L ũ L + S( v) L M + S( v) L. Now we bound S( v) L. According to Lemma.1, S( v(x)) C σ(x y) R x y ( 1θ(x) 1+β 1 θ(y))dy. To further the estimate, we split the integral over R into three parts: x y ρ, ρ < x y L, x y > L for 0 < ρ < L to be specified later. By Hölder s inequality, the integral over x y ρ is bounded by 1 x y 1θ(x) 1+β 1 θ(y) dy D( θ)ρ 1 β. x y ρ Through integration by parts, the integral over ρ < x y L is bounded by 1 θ(x) θ(y) dy x y +β δ 1 dy x y +β ρ< x y L = Cδρ β. ρ< x y L We remark that a rigorous justification of the estimate above can be carried out through smooth cutoff functions. Again through integration by parts, the integral over x y > L is bounded by Therefore, x y >L 1 x y θ(x) θ(y) dy +β C θ 0 L L β. S( u) L θ(x) C D( θ)ρ 1 β θ(x) +Cδρ β θ(x) +(M +C θ 0 L L β ) θ(x) 1 D( θ)+cρ β θ(x) 4 +Cδρ β θ(x) +(M +C θ 0 L L β ) θ(x), where C s are pure constants depending on β only. If we choose ρ and δ as then 1 ρ = and δ = (4CC 0 θ 0 β L ) 1 β θ(x) 1 C 1+β (4C 0 θ 0 β L ) β S( u) L θ(x) 1 θ +β D( θ)+ +(M +C θ C 0 θ 0 β 0 L L β ) θ(x). L Inserting this bound in (3.9) yields 1 ( t +u +Λ β ) θ + θ(x) +β C 0 θ 0 β L (M +C θ 0 L L β ) θ(x). (3.10),

12 1 Q. JIU, C. MIAO, J. WU AND Z. ZHANG This differential inequality allows us to conclude that where C is given by θ(x,t) max{ θ 0 L, C(M, θ 0 L )}, (3.11) C = ( C 0 (M +C θ 0 L L β) 1 β θ 0 L. (3.1) To show (3.11), we define Γ : [0, ) [0, ) to be a nondecreasing C convex function satisfying Γ(ρ) = 0 for 0 ρ max{ θ 0 L, C} and Γ(ρ) > 0 for ρ > max{ θ 0 L, C}. Multiplying (3.10) by Γ ( θ ) and applying the lower bound Λ β (Γ(f)) Γ (f)λ β f, (3.13) where is valid for any β (0,) and any convex function Γ (see [1, 17]), we obtain 1 ( t +u +Λ β )Γ( θ ) [ ] (M +C θ 0 L L β ) θ(x) 1 θ(x,t) β Γ ( θ ). C β It is easy to verify that the right-hand side of the inequality above is always less than or equal to zero due to the definition of Γ. Therefore, ( t +u +Λ β )Γ( θ ) 0. Thanks to u = 0 and the positivity of integral R f(x) f(x) p Λ β f(x)dx 0, we obtain Γ( θ ) L p Γ( θ 0 ) L p for any 1 p <. Letting p, we have Γ( θ(t) ) L Γ( θ 0 ) L 0, which implies (3.11). This completes the proof of Proposition 3.. To prove Proposition 3.3, we need the lower bound in the following lemma. This lemma can be proven by following the lines of that for Lemma 3.4 (see [16]). Lemma 3.5. Let β (0,). Let h R and δ h θ(x,t) θ(x+h,t) θ(x,t). Then the following pointwise lower bound holds, δ h θλ β δ h θ 1 Λβ (δ h θ) +D h (δ h θ)+c δ hθ(x,t) +β θ β L h β, where C is a constant depending on β only, and R (δ h θ(x,t) δ h θ(y,t)) D h (δ h θ) = cp.v. dy. x y +β Proof of Proposition 3.3. Let h R and consider the evolution of where Φ(h) = e Ψ(h) with Ψ satisfying g(x,t;h) = (δ h θ(x,t)) Φ(h) (3.14) Ψ(h) = Ψ( h ) is nondecreasing, Ψ(0) = 0, Ψ(h) + as h. (3.15)

13 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 13 An explicit form of Ψ will be specified later. First of all, δ h θ satisfies ( t +u +δ h u h +Λ β )δ h θ = 0, (3.16) where δ h u(x,t) = u(x+h,t) u(x,t). Multiplying (3.16) by δ h θ(x,t)φ(h) and applying Lemma 3.5, we obtain ( t +u +δ h u h +Λ β )g +D h (δ h θ)φ(h)+cφ(h) δ hθ(x,t) +β (δ h θ(x,t)) δ h u h Φ(h). Noticing that h Φ(h) = Φ(h) h Ψ(h), we have (δ h θ(x,t)) δ h u h Φ(h) = (δ h θ(x,t)) Φ(h)δ h u h Ψ(h). Therefore, if we write L t +u +δ h u h +Λ β, we obtain g 1+β θ β L h β Lg +D h (δ h θ)φ(h)+c Φ β (h) θ g δ hu Ψ ( h ). (3.17) β L h β Thanks to the assumptions on ũ, Therefore, u = ũ+v satisfies According to Lemma.1, δ h ũ M min{m h, u L (0,T;L )} <. v(x,t) = C δ h u M + δ h v. (3.18) R (x y) x y 1+β 1θ(y)dy. Since β (0,1) and θ is bounded and decays sufficiently fast at, we obtain after integration by parts, (x 1 y 1 )(x y ) v 1 (x,t) = (1+β)C θ(y)dy, R x y 3+β [ 1 v (x,y) = C x y + (1+β)(x ] 1 y 1 ) θ(y) dy. 1+β x y 3+β Therefore, δ h v C C R 1 R x y 1+β δ hθ(y,t) dy 1 x y 1+β δ hθ(y,t) dy+c x y 1 x y >1 1 x y 1+β δ hθ(y,t) dy C θ 0 L L. (3.19) Inserting the estimates of (3.18) and (3.19) in (3.17) and using Φ(h) 1, we obtain Lg +C 1 g 1+ β θ 0 β L h β C g( M + θ 0 L L )Ψ ( h ), (3.0)

14 14 Q. JIU, C. MIAO, J. WU AND Z. ZHANG where C 1 and C are constants depending on β only. Now we choose Ψ( h ) = C 1 (δ 0 ) β C (1 β) θ 0 β L ( M + θ 0 L L ) h 1 β. Certainly this choice satisfies (3.15). Then (3.0) becomes Lg We can then conclude from (3.1) that C 1 g h β θ 0 β L (δ β 0 g β ). (3.1) g(x,t;h) δ 0 for h L. This can be proven by following a similar argument as in the proof of Proposition 3.. Moreprecisely, we multiply (3.1) by Γ (g)with Γ being a non-decreasing smoothconvex function on [0, ) satisfying Γ(ρ) = 0 for 0 ρ δ0 and positive for ρ δ 0. Thanks to (3.13), ) L ( Γ(g) C 1 g h β θ 0 β L (δ β 0 g β ) Γ (g) 0. By first estimating the L q -norm of Γ(g) and then sending q, we obtain that the L -norm of Γ(g) is bounded by its initial L -norm. Since initially g δ 0 and thus the initial L -norm of Γ(g) is zero, we have L -norm of Γ(g) is zero for all time. Therefore g δ 0 and δ h θ(x,t) δ 0 (Φ(h)) 1 δ0 e 1 Ψ(L) 4δ 0. This completes the proof of Proposition Besov spaces and a commutator estimate This section provides the definitions of some of the functional spaces and related facts to be used in the subsequent sections. More details can be found in several books and many papers (see, e.g., [4, 5, 38, 4, 44]). In addition, we prove a commutator estimate to be used extensively in the sections that follow. To introduce the Besov spaces, we start with a few notation. S denotes the usual Schwarz class ands itsdual, thespace oftempered distributions. S 0 denotes a subspace of S defined by { } S 0 = φ S : φ(x)x γ dx = 0, γ = 0,1,, R d and S 0 denotes its dual. S 0 can be identified as S 0 = S /S 0 = S /P where P denotes the space of multinomials. For each j Z, we write A j = { ξ R d : j 1 ξ < j+1}. (4.1)

15 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 15 The Littlewood-Paley decomposition asserts the existence of a sequence of functions {Φ j } j Z S such that and supp Φ j A j, j= Φj (ξ) = Φ 0 ( j ξ) or Φ j (x) = jd Φ 0 ( j x), Φ j (ξ) = { 1, if ξ R d \{0}, 0, if ξ = 0. Therefore, for a general function ψ S, we have Φ j (ξ) ψ(ξ) = ψ(ξ) for ξ R d \{0}. j= In addition, if ψ S 0, then That is, for ψ S 0, and hence j= Φ j (ξ) ψ(ξ) = ψ(ξ) for any ξ R d. j= j= Φ j ψ = ψ Φ j f = f, f S 0 in the sense of weak- topology of S 0. For notational convenience, we define j f = Φ j f, j Z. (4.) Definition 4.1. For s R and 1 p,q, the homogeneous Besov space B s p,q consists of f S 0 satisfying We now choose Ψ S such that Then, for any ψ S, and hence f B s p,q js j f L p l q <. Ψ(ξ) = 1 Ψ ψ + Ψ f + Φ j (ξ), ξ R d. j=0 Φ j ψ = ψ j=0 Φ j f = f (4.3) j=0

16 16 Q. JIU, C. MIAO, J. WU AND Z. ZHANG in S for any f S. To define the inhomogeneous Besov space, we set 0, if j, j f = Ψ f, if j = 1, Φ j f, if j = 0,1,,. (4.4) Definition 4.. The inhomogeneous Besov space B s p,q with 1 p,q and s R consists of functions f S satisfying f B s p,q js j f L p l q <. The Besov spaces B p,q s and Bp,q s with s (0,1) and 1 p,q can be equivalently defined by the norms (Rd ) ( f(x+t) f(x) L p) q 1/q f Bs p,q = dt, t d+sq ) (Rd ( f(x+t) f(x) L p) q 1/q f B s p,q = f L p + dt. t d+sq When q =, the expressions are interpreted in the normal way. Besides the Fourier localization operators j, the partial sum S j is also a useful notation. For an integer j, S j j 1 k= 1 where k is given by (4.4). For any f S, the Fourier transform of S j f is supported on the ball of radius j. Bernstein s inequality is a useful tool on Fourier localized functions and these inequalities trade integrability for derivatives. The following proposition provides Bernstein type inequalities for fractional derivatives. k, Proposition 4.3. Let α 0. Let 1 p q. 1) If f satisfies supp f {ξ R d : ξ K j }, for some integer j and a constant K > 0, then ) If f satisfies ( ) α f L q (R d ) C 1 αj+jd(1 p 1 q ) f L p (R d ). supp f {ξ R d : K 1 j ξ K j } for some integer j and constants 0 < K 1 K, then C 1 αj f L q (R d ) ( ) α f L q (R d ) C αj+jd(1 p 1 q ) f L p (R d ), where C 1 and C are constants depending on α,p and q only. The rest of this section is devoted to the proof of a commutator estimate. We need the following lemma.

17 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 17 Lemma 4.4. Let δ (0,1). Let q, q 1, q,r 1,r [1, ] satisfying = 1 q 1 q q, = 1. r 1 r If f B δ q 1,r 1 (R d ), x δ+ d r 1 φ L r (R d ), and g L q (R d ), then φ (fg) fφ g L q C x δ+ d r 1 φ L r f B δ q 1,r 1 g L q, (4.5) where C is a constant independent of f,g and h. In the borderline case when δ = 1, (4.5) still holds if f Bδ q1,r 1 is replaced by f L q 1 and r = 1. The case when δ = 1, q 1 = q and q = was obtained in [7, p.153]. The fractional case with δ (0,1), q 1 = q and q = was obtained by Chae and Wu in [10]. Here we need to make use of the more general case when q 1 > q. Proof of Lemma 4.4. By Minkowski s inequality followed by Hölder s inequality, for any q [1, ], [ q ] 1/q φ (fg) f(φ g) L q = φ(z)(f(x) f(x z))g(x z)dz dx φ(z) f( ) f( z)) L q 1 g L q dz f( ) f( z)) L q 1 g L q z δ+ d r 1 φ(z) z δ+ d r 1 L r (4.5) then follows from the definition of B δ q 1,r 1. Proposition 4.5. Let α (0,1). Let s (0,1) and δ (0,1) satisfy s+1 α δ < 0. Let q [, ), r [1, ], q 1,q [, ] satisfy 1 q = 1 q q. Then, for a constant C, L r 1 [R α,f]g B s q,r C f g B δ q1, B. (4.6) q s+1 α δ,r When δ = 1, (4.6) is still valid if f B δ q1, is replaced by f L p 1. Proof. Let k 1 be an integer. By the notion of paraproducts, we write k [R α,f]g = J 1 +J +J 3, where J 1 = k (R α (S j 1 f j g) S j 1 f R α j g), j k J = k (R α ( j f S j 1 g) j f R α S j 1 g), j k J 3 = j k 1 k (R α ( j f ) j g) j f R α j g

18 18 Q. JIU, C. MIAO, J. WU AND Z. ZHANG with j = j 1 + j + j+1. Wefirstnotethat, ifthefouriertransformoff issupported in the annulus around radius j, then R α F can be represented as a convolution, R α F = h j F, h j (x) = (d+1 α)j h 0 ( j x) (4.7) for a function h 0 in the Schwartz class S. This can be obtained by simply examining the Fourier transform of R α F. By the definition of Bq,r s, [R α,f]g B s q,r = sk k [R α,f]g L q l r sk J 1 L q l r + sk J L q l r + sk J 3 L q l r. (4.8) Applying Lemma 4.4, we find Thus, J 1 L q C (1 α)k x δ dk h 0 ( k x) L 1 S j 1 f Bδ q1, jg L q By Bernstein s inequality, we have C (1 α δ)k f B δ q 1, jg L q. sk J 1 L q l r C f B g q δ 1, B. (4.9) q s+1 α δ,r J L q C (1 α δ)j f S B δ q1, j 1g L q C (1 α δ)j f B δ q1, m g L q C f B δ q1, m k m k (1 α δ)(k m) (1 α δ)m m g L q. Since s+1 α δ < 0, we obtain, by applying Young s inequality for series, sk J L q l r C f g B δ q1, B. (4.10) q s+1 α δ,r Similarly, we have J 3 L q C j k 1 (1 α δ)j f B δ q1, jg L q. Therefore, for s > 0, by Young s inequality for series, sk J 3 L q l r C sk j k 1 (1 α δ)j f B δ q1, jg L q l r C f g B δ q1, B. (4.11) q s+1 α δ,r Combining (4.8), (4.9), (4.10) and (4.11), we obtain the desired bound in (4.6). This completes the proof of Proposition 4.5. We have used a simple inequality, namely (3.6) in the proof of Theorem 1. and we now prove it. Its proof is deferred to this section since it involves the Littlewood-Paley

19 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 19 decomposition and Bernstein s inequality described above. By Bernstein s inequality, Λ 3+β 1 θ ǫ L 1 Λ 3+β 1 θ ǫ L + j Λ 3+β 1 θ ǫ L j=0 C 1 θ ǫ L + βj j θ ǫ L C θ ǫ L + j=0 (β 1)j j θ ǫ L, j=0 which reduces to (3.6) by recalling that β (0,1). 5. Global L -bound for G This section establishes a global a priori bound for G L. Recall that G = ω R α θ with R α = Λ α 1 (5.1) and G satisfies t G+u G+Λ α G = [R α,u ]θ+λ 1 α 1 θ. (5.) This global bound is valid for any 4 < α < 1. 5 Theorem 5.1. Assume (u 0,θ 0 ) satisfies the assumptions stated in Theorem 1.1. Let (u,θ) be the corresponding solution of (1.1). If 4 < α < 1, then G defined in (5.1) has 5 a global L bound, namely for any T > 0 and t T, G(t) L + t 0 Λ α G(τ) L dτ C(T,u 0,θ 0 ), where C(T,u 0,θ 0 ) is a constant depending on T and u 0 L and θ 0 L L only. To prove Theorem 5.1, the following elementary global a priori bounds will be used. Lemma 5.. Assume that u 0 and θ 0 satisfy the conditions stated in Theorem 1.1. Then the corresponding solution (u,θ) of (1.1) obey the global bounds, for any t > 0, θ(t) L r θ 0 L r for any r, u(t) L + t 0 Λ α u(τ) L dτ ( u 0 L +t θ 0 L ). Proof of Theorem 5.1. Taking the inner product of (5.) with G, we obtain, after integration by parts, 1 d dt G L + Λα G L = J 1 +J, (5.3) where J 1 = G[R α,u ]θdx = G [R α,u]θdx, J = GΛ 1 α 1 θdx.

20 0 Q. JIU, C. MIAO, J. WU AND Z. ZHANG Let q 1,q [, ] satisfying 1 q q = 1. By Hölder s inequality and Proposition 4.5, Setting q 1 = α J 1 Λ α G L Λ 1 α Λ 1 [R α,u]θ L and noticing that C Λ α G L u B α θ q1., B 5 α q, u = ( ) 1 ω = ( ) 1 (G+Λ α 1 θ), we obtain, by the Hardy-Littlewood-Sobolev inequality, u B α q 1, C Λα u L q 1 C Λ α 1 ω L q 1 So, we get C Λ α 1 G L q 1 +C θ L q 1 C G L +C θ 0 L q 1. (5.4) u B α q1, u Bα q1, +C u L C G L +C θ 0 L q 1 +C u L. Thanks to the embedding L q B 5 α q, Therefore, by Young s inequality θ B 5 α q, J Λα G L +C θ 0 L q G L +C θ 0 L q where C s are pure constants. Since 4 5 < α < 1, J θ 0 L Λ α G L when α > 4, we have 5 C θ 0 L q. (5.5) θ 0 L Λ α 4( G 1 α 1) L G 5 4 α L 1 4 Λα G L +C θ 0 α 5α 4 L G L. ( θ0 L q 1 + u 0 L +t θ 0 L ), Inserting the bounds for J 1 and J in (5.3) and applying Gronwall s inequality yield the desired bound. This completes the proof of Theorem Global L q -bound for G with < q < q 0 This section proves a global bound for G L q with < q < q 0, where q 0 is specified in (6.1). This global bound is valid for any 4 5 < α < 1. Theorem 6.1. Let 4 5 < α < 1 and α + β = 1. Consider (1.1) with u 0 H and θ 0 B,1. Let (u,θ) be the corresponding smooth solution. Assume that < q < q 0 with q 0 given by q 0 8 4α 8 7α. (6.1) Then, for any T > 0 and t T, G(t) L q C(T,u 0,θ 0 ),

21 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 1 where C(T,u 0,θ 0 ) is a constant depending on T, u 0 and θ 0 only. Proof. Taking the inner product of (5.) with G G q, we have 1 d qdt G q L + G G q Λ α Gdx = K q 1 +K, (6.) where K 1 = K = G G q [R α,u ]θdx, G G q Λ 1 α 1 θdx. By a pointwise inequality for fractional Laplacians (see [17]) and a Sobolev embedding inequality, G G q Λ α Gdx C Λ α q G dx C 0 G q, L α q where C,C 0 > 0 are constants. To bound K 1, we have for any s (0,1), K 1 Λ s (G G q )Λ 1 s (Λ 1 [R α,u]θ)dx Λ s (G G q ) L Λ 1 s ([R α,u]θ) L. By Proposition 4.5, (5.4) and (5.5), we obtain Λ 1 s [R α,u]θ L = [R α,u]θ H1 s C u B α θ q1, B s α +C u q, L θ 0 L C( G L + θ 0 L q 1 + u L ) θ 0 L q +C u L θ 0 L B(t), where α < s, q 1 =, 1 α q q = 1 and B(t) is a smooth function of t depending on u 0 L and θ 0 L L q 1 L q only. Choosing s satisfying +s α ( α) q and applying Lemma 6. below, we obtain Λ s (G G q ) L C G q L α q G H +s α ( α) q Combining the estimates, we obtain, for any t T, = α K 1 B(t) G q L α q G H α C 0 G q L α q = C G q L α q α G H. +C G q H α, wherec isaconstant depending ont andthenormsoftheinitial data. Since α < s, K (1+Λ) s (G G q ) L (1+Λ) α s θ L C( Λ s (G G q ) L + G q 1 L ) θ 0 L C 0 G q L q α +C G q H α +C G q 1 L θ 0 L.

22 Q. JIU, C. MIAO, J. WU AND Z. ZHANG Inserting the estimates for K 1 and K in (6.) and using the fact in Theorem 5.1 that, for q and any T > 0, T 0 G q H α dt <, we then obtain the desired the global bound. The conditions on the indices are fulfilled if q, s α < 0, +s α ( α) q α > 4 5, This completes the proof of Theorem 6.1. q < q 0 8 4α 8 7α. = α We have used the following lemma in the proof of Theorem 6.1. This lemma generalizes a previous inequality of [7, p.170]. Lemma 6.. Let s (0,1), α (0,1) and q [, ). Then, Λ s (G G q ) L C G q L q α G B s. q ( α) q(1 α), Especially, Λ s (G G q ) L C G q L q α ( α) G H+s α q. (6.3) Proof. Identifying H s with the Besov space B, s, we have Λ s (G G q ) L = G G q ( z) G G q ( ) L z +s dz. Using the simple inequality, for q, and Hölder s inequality, we have a a q b b q C a b ( a q + b q ) G G q ( z) G G q ( ) L C G (q ) L α q Therefore, Λ s (G G q ) L C G (q ) L α q G( z) G( ) L G Bs q ( α) q(1 α), q ( α) q(1 α). (6.3) holds due to the embedding proof of Lemma 6.. +s α ( α) q B, B s q,. This completes the ( α) q(1 α)

23 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 3 7. Proof of Theorem 1.1 This section proves Theorem 1.1. To do so, we need more regularity for G. By fully exploiting the dissipation, we are able to show in Proposition 7.1 below that G in B s q, with any s 3α and q < q 0 is bounded for all time. Proposition 7.1. Consider the IVP (1.1) with α 0 < α < 1. Assume the initial data (u 0,θ 0 ) satisfies the conditions stated in Theorem 1.1 and let (u,θ) be the corresponding solution. If the indices s and q satisfy 0 < s 3α, α 1 < q < q 0 8 4α 8 7α, (7.1) then G obeys the global a priori bound, for any T > 0, G L (0,T;B s q, ) C, where C is a constant depending on T and the initial norms only. This proposition will be proven at the end of this section. With this proposition at our disposal, we are ready to prove Theorem 1.1. Proof of Theorem 1.1. Thekeypartoftheproofistoestablishtheglobala prioribounds of the solution in the functional setting (1.3). According to Proposition 7.1, for any T > 0, and s and q satisfying (7.1), we have Since u = 1 ω, we have G ω Λ α 1 θ L (0,T;Bq, s ), (7.) u = 1 G+ 1 Λ α 1 θ ũ+v. Since α 0 < α < 1, we can choose q satisfying (7.1) such that (3α ) < 0. q In fact, thanks to α 0 < α, we have 3α < 8 4α 8 7α and any q satisfies 3α < q < q 0 8 4α 8 7α wouldwork. Bytheembedding B,1 0 L, Bernstein s inequalityandtheboundedness of Riesz transforms on L q, ũ L 1 j G L C j= 1 j q j G L q j= 1 G B s q, j= 1 j( q s).

24 4 Q. JIU, C. MIAO, J. WU AND Z. ZHANG Taking s 3α but close to 3α, we have s < 0 and thus q M ũ L (0,T;L ) <. Applying Theorem 1. to the equation for θ, { t θ +u θ +Λ 1 α θ = 0, u = ũ+v, v 1 Λ α 1 θ, we obtain θ C([0,T];H 1 ), θ L (0,T;L ) C(M, θ 0 L, θ 0 L,T) <. (7.3) According to the vorticity equation t ω +u ω +Λ α ω = 1 θ, (7.4) we have, thanks to ω 0 B,1 σ 1 (R ) with σ 5 and B1,1 (R ) L (R ), ω(t) L ω 0 L + T 0 1 θ(τ) L dτ. (7.5) (7.3) and(7.5) will lead us to the global a priori bounds in (1.3). To this goal, we employ the following estimate, for any integer j 1 and any t (0,T], αj t ( t ) j ω L dτ C ω 0 L + 1 θ(τ) L dτ j ( 1+t+ u L + ω L 1 t L ) h(t), (7.6) where C isaconstant independent ofj. (7.6) canbeestablished by applying Besov space type estimates on (7.4) and more details can be found in [38, p.134]. As a consequence, we have, for any 0 < α < α, t t α j j ω L dτ = α j j ω L dτ 0 j= 1 j= 1 h(t) 0 (α α)j /(j +1) Ch(t). (7.7) Therefore, by Bernstein s inequality and j u L C j ω L for j 0, u L 1 u L + j u L C u L +C j ω L. (7.8) j=0 j= 1 Combining (7.7) and (7.8) yields, for any 0 < t T, t 0 j=0 u(τ) L dτ <. (7.9) Through standard Besov space type estimates, we can show that Y(t) u(t) B σ,1 + θ(t) B,1 obeys the differential inequality d dt Y(t)+C( u B σ+α + θ,1 B +β) C(1+ u(t) L + θ(t) L )Y(t).,1

25 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 5 Gronwall s inequality together with (7.3) and (7.9) yields the desired property in (1.3). The continuity of u and θ in time, namely u C([0,T];B,1 σ ) and θ C([0,T];B,1 ), follows from a standard approach (see, e.g., [36, p.111] or [4, p.138]). We thus have obtained all the global a priori bounds in (1.3). The global existence part can then be obtained by first establishing the local existence through a standard process such as the successive approximation and then extending to all time via the global a priori bounds. The uniqueness of solutions in the functional setting (1.3) is clear and we thus omit the details. This completes the proof of Theorem 1.1. We finally prove Proposition 7.1. Proof of Proposition 7.1. Recall that G satisfies t G+u G+Λ α G = [R α,u ]θ+λ 1 α 1 θ. (7.10) Let j 1. Applying j to (7.10) and taking the inner product with j G j G q, we obtain 1 d qdt jg q L + q j G j G q Λ α j Gdx = K (j) 1 +K (j) +K (j) 3, (7.11) where K (j) 1 = j G j G q j [R α,u ]θdx, K (j) = j G j G q j (u G)dx, (7.1) K (j) 3 = j G j G q j (Λ 1 α 1 θ)dx. For j 0, the Fourier transform of j G is supported away from the origin and the dissipative part admits a lower bound, j G j G q Λ α j Gdx C αj j G q L q (7.13) for a constant C that depends on q and α only (see, e.g. [11, 45]). For j = 1, the dissipative part is still nonnegative and can be neglected. By Hölder s inequality, K (j) 1 j G q 1 L q j[r α,u ]θ L q. Furthermore, j [R α,u ]θ L q can be estimated in a similar fashion as in the proof of Proposition 4.5 and is bounded by In addition, as in (5.4), j [R α,u ]θ L q C ( α)j u Bα q, θ L + C u L θ L. u Bα q, C Λ α u L q C Λ α 1 G L q +C θ L q C G L q + C θ 0 L q, (7.14) where 1 q = 1 q α + 1. For α (α 0,1), we have < q < q andthus G L q < C. Therefore, K (j) 1 B(t) ( α)j j G q 1 L q, (7.15)

26 6 Q. JIU, C. MIAO, J. WU AND Z. ZHANG where B 1 (t) isasmoothfunctionof tthat depends onthe initial normsonly. To estimate K (j), we apply the notion of paraproducts to write j (u G) = J 1 +J +J 3 +J 4 +J 5, (7.16) where J 1 = ( j (S k 1 u k G) S k 1 u j k G), j k J = (S k 1 u S j u) j k G, j k J 3 = S j u j G, J 4 = j ( k u S k 1 G), j k J 5 = j ( k u k G) k j 1 with k = k 1 + k + k+1. Inserting (7.16) in (7.1) naturally splits the integral in K (j) into five parts K (j) = K (j) 1 +K (j) +K (j) 3 +K (j) 4 +K (j) 5. By Hölder s inequality, K (j) 1 j G q 1 L ( q j (S k 1 u k G) S k 1 u j k G) L q. j k Since the summation above is for k satisfying j k and can be replaced by a constant multiple of the representative term with k = j, we obtain, by a standard commutator estimate and Bernstein s inequality, K (j) 1 C jg q q 1 (1 α)j Λ α S j 1 u L q j G L C j G q q 1 (1 α+ q )j Λ α u L q j G L q C (1 α+ q )j Λ α u L q j G q L q. The second term K (j) can be bounded by K (j) C jg q 1 L q ju L q j j G L. C (1 α+ q )j Λ α j u L q j G q L q, where we have used the lower bound part of the Bernstein inequality αj j u L q C Λ α j u L q. This inequality is valid for j 0. In the case when j = 1, this inequality is not needed and it suffices to apply the upper bound part of the Bernstein inequality 1 u L q C 1 u L.

27 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 7 Due to the divergence-free condition u = 0, we have K (j) 3 and Bernstein s inequality K (j) 4 C j G q 1 L ju S q j 1 G L q C αj j G q 1 L q Λα j u L q m j 1 C (1 α+ q )j j G q 1 L q Λα j u L q K (j) 5 C j G q 1 L q k j 1 C (1 α+ q )j j G q 1 L q = 0. By Hölder s inequality (1+ q )m m G L q m j 1 j j ( k u k G) L q k j 1 (1+ q )(m j) m G L q. (j k)(α q ) Λ α k u L q k G L q. Collecting the estimates for K (j) and bounding Λ α u L q as in (7.14), we have [ K (j) C(1 α+ q )j j G q 1 L q j G L q + (1+ q )(m j) m G L q + k j 1 (j k)(α q ) k G L q By Hölder s inequality and Bernstein s inequality, m j 1 ]. (7.17) K (j) 3 (1 α)j j G q 1 L q jθ L q (1 α)j j G q 1 L q θ 0 L q. (7.18) Combining (7.11), (7.13), (7.15), (7.17) and (7.18), we obtain d dt jg L q +C αj j G L q CB(t) ( α)j + C (1 α+ q )j L(t), where B(t) is a smooth function of t depending on the initial norms only and L = j G L q + (1+ q )(m j) m G L q + (j k)(α q ) k G L q. m j 1 Integrating in time leads to k j 1 j G(t) L q e Cαjt j G(0) L q +C B(t) ( 3α)j t +C (1 α+ q )j e Cαj (t τ) L(τ)dτ, 0 where B(t) is a smooth function of t depending on the initial norms only. Multiplying each side by sj with 0 < s 3α and taking sup with respect to j, we find that G(t) B s q, G(0) B s q, +C B(t)+L 1 +L +L 3, (7.19)

28 8 Q. JIU, C. MIAO, J. WU AND Z. ZHANG where t L 1 = sup C (1 α+ q )j e Cαj (t τ) sj j G(τ) L q dτ, j 1 0 t L = sup C (1 α+ q )j e Cαj (t τ) sj (1+ q )(m j) m G L q dτ, j 1 0 m j 1 t L 3 = sup C (1 α+ q )j e Cαj (t τ) sj (j k)(α q ) k G L q dτ. j 1 0 k j 1 Since 1 α+ q < 0, we choose an integer j 0 such that, for j j 0, C (1 α+ q )j 1 8. Then L 1 sup C (1 α+ 1 j j 0 + sup j 0 +1 j< q )j t 0 e Cαj (t τ) sj j G(τ) L q dτ t C (1 α+ q )j e Cαj (t τ) sj j G(τ) L q dτ 0 C G L (0,T;L q ) (1 α+ q +s)j G(t) L (0,T;B s q, ) L and L 3 obey the same bound. Inserting these bounds in (7.19), we find G L (0,T;B s q, ) C(1+ G(0) B s q, )+C G L (0,T;L q ), where we have written C for max 0 t T B(t). Since u0 B σ,1 with σ 5 and θ 0 B,1, ω 0 B σ 1,1 (R ) B s q, (R ), G(0) = ω 0 Λ α 1 θ 0 B s q,. This finishes the proof of Proposition 7.1. Acknowledgements Jiu was supported by NSF of China under grants No and No Miao was supported by NSF of China under grants No and No Wu was partially supported by NSF grant DMS Zhang was supported by NSF of China under grants No and No , Program for New Century University Talents and Fok Ying Tung Education Foundation. In addition, Jiu and Wu were partially supported by a special grant from NSF of China under No References [1] H. Abidi and T. Hmidi, On the global well-posedness for Boussinesq system, J. Differential Equations 33 (007), [] D. Adhikari, C. Cao and J. Wu, The D Boussinesq equations with vertical viscosity and vertical diffusivity, J. Differential Equations 49 (010),

29 THE D BOUSSINESQ EQUATIONS WITH CRITICAL DISSIPATION 9 [3] D. Adhikari, C. Cao and J. Wu, Global regularity results for the D Boussinesq equations with vertical dissipation, J. Differential Equations 51 (011), [4] H. Bahouri, J.-Y. Chemin and R. Danchin, Fourier Analysis and Nonlinear Partial Differential Equations, Springer, 011. [5] J. Bergh and J. Löfström, Interpolation Spaces, An Introduction, Springer-Verlag, Berlin- Heidelberg-New York, [6] J.R. Cannon and E. DiBenedetto, The initial value problem for the Boussinesq equations with data in L p, Lecture Notes in Math., Vol Springer, Berlin, 1980, pp [7] C. Cao and J. Wu, Global regularity for the D anisotropic Boussinesq equations with vertical dissipation, arxiv: v1 [math.ap] 1 Aug 011. [8] L. Caffarelli and A. Vasseur, Drift diffusion equations with fractional diffusion and the quasigeostrophic equation, Ann. Math. 171 (010), [9] D. Chae, Global regularity for the D Boussinesq equations with partial viscosity terms, Adv. Math. 03 (006), [10] D. Chae and J. Wu, The D Boussinesq equations with logarithmically supercritical velocities, Adv. Math. 30 (01), [11] Q. Chen, C. Miao and Z. Zhang, A new Bernstein s inequality and the D dissipative quasigeostrophic equation, Comm. Math. Phys. 71 (007), [1] P. Constantin, Euler equations, Navier-Stokes equations and turbulence, Mathematical foundation of turbulent viscous flows, 143, Lecture Notes in Math., 1871, Springer, Berlin, 006. [13] P. Constantin and C.R. Doering, Infinite Prandtl number convection, J. Statistical Physics 94 (1999), [14] P. Constantin, G. Iyer and J. Wu, Global regularity for a modified critical dissipative quasigeostrophic equation, Indiana Univ. Math. J. 57 (008), [15] P. Constantin, A. Majda and E. Tabak, Formation of strong fronts in the -D quasi-geostrophic thermal active scalar, Nonlinearity 7 (1994), [16] P. Constantin and V. Vicol, Nonlinear maximum principles for dissipative linear nonlocal operators and applications, Geom. Funct. Anal. (01), DOI: /s [17] A. Córdoba and D. Córdoba, A maximum princple applied to quasi-geostroohhic equations; Comm. Math. Phys. 49 (004), [18] D. Córdoba and C. Fefferman, Scalars convected by a two-dimensional incompressible flow, Commun. Pure Appl. Math. 55 (00), [19] X. Cui, C. Dou and Q. Jiu, Local well-posedness and blow up criterion for the inviscid Boussinesq system in Hölder spaces, J. Partial Differential Equations 5 (01), [0] R. Danchin and M. Paicu, Global well-posedness issues for the inviscid Boussinesq system with Yudovich s type data, Comm. Math. Phys. 90 (009), [1] R. Danchin and M. Paicu, Global existence results for the anisotropic Boussinesq system in dimension two, Math. Models Methods Appl. Sci. 1 (011), [] W. E and C. Shu, Samll-scale structures in Boussinesq convection, Phys. Fluids 6 (1994), [3] A.E. Gill, Atmosphere-Ocean Dynamics, Academic Press (London), 198. [4] T. Hmidi, On a maximum principle and its application to the logarithmically critical Boussinesq system, Anal. Partial Differential Equations 4 (011), [5] T. Hmidi and S. Keraani, On the global well-posedness of the two-dimensional Boussinesq system with a zero diffusivity, Adv. Differential Equations 1 (007), [6] T. Hmidi and S. Keraani, On the global well-posedness of the Boussinesq system with zero viscosity, Indiana Univ. Math. J. 58 (009), [7] T. Hmidi, S. Keraani and F. Rousset, Global well-posedness for a Boussinesq-Navier-Stokes system with critical dissipation, J. Differential Equations 49 (010), [8] T. Hmidi, S. Keraani and F. Rousset, Global well-posedness for Euler-Boussinesq system with critical dissipation, Comm. Partial Differential Equations 36 (011), [9] T. Hou and C. Li, Global well-posedness of the viscous Boussinesq equations, Discrete and Cont. Dyn. Syst. 1 (005), 1-1.

30 30 Q. JIU, C. MIAO, J. WU AND Z. ZHANG [30] D. KC, D. Regmi, L. Tao and J. Wu, The D Euler-Boussinesq equations with a logarithemically supercritical velocity, submitted for publication. [31] D. KC, D. Regmi, L. Tao and J. Wu, The D Boussinesq-Navier-Stokes equations with logarithmically supercritical dissipation, submitted for publication. [3] A. Kiselev, F. Nazarov and A. Volberg, Global well-posedness for the critical D dissipative quasigeostrophic equation, Invent. Math. 167 (007), [33] M. Lai, R. Pan and K. Zhao, Initial boundary value problem for two-dimensional viscous Boussinesq equations, Arch. Ration. Mech. Anal. 199 (011), [34] A. Larios, E. Lunasin and E.S. Titi, Global well-posedness for the D Boussinesq system without heat diffusion and with either anisotropic viscosity or inviscid Voigt-a regularization, arxiv: v1 [math.ap] 5 Oct 010. [35] A.J. Majda, Introduction to PDEs and Waves for the Atmosphere and Ocean, Courant Lecture Notes in Mathematics 9, AMS/CIMS, 003. [36] A.J. Majda and A.L. Bertozzi, Vorticity and Incompressible Flow, Cambridge University Press, 001. [37] C. Miao and L. Xue, On the global well-posedness of a class of Boussinesq- Navier-Stokes systems, NoDEA Nonlinear Differential Equations Appl. 18 (011), [38] C. Miao, J. Wu and Z. Zhang, Littlewood-Paley Theory and its Applications in Partial Differential Equations of Fluid Dynamics, Science Press, Beijing, China, 01 (in Chinese). [39] H.K. Moffatt, Some remarks on topological fluid mechanics, in: R.L. Ricca (Ed.), An Introduction to the Geometry and Topology of Fluid Flows, Kluwer Academic Publishers, Dordrecht, The Netherlands, 001, [40] K. Ohkitani, Comparison between the Boussinesq and coupled Euler equations in two dimensions. Tosio Kato s method and principle for evolution equations in mathematical physics (Sapporo, 001). Surikaisekikenkyusho Kokyuroku No. 134 (001), [41] J. Pedlosky, Geophysical Fluid Dyanmics, Springer-Verlag, New York, [4] T. Runst and W. Sickel, Sobolev Spaces of fractional order, Nemytskij operators and Nonlinear Partial Differential Equations, Walter de Gruyter, Berlin, New York, [43] E. Stein, Singular Integrals and Differentiability Properties of Functions, Princeton Unviersity Press, Princeton, NJ, [44] H. Triebel, Theory of Function Spaces II, Birkhauser Verlag, 199. [45] J. Wu, Lower bounds for an integral involving fractional Laplacians and the generalized Navier- Stokes equations in Besov spaces, Comm. Math. Phys. 63 (006), [46] J. Wu and X. Xu, Well-posedness and inviscid limits of the Boussinesq equations with fractional Laplacian dissipation, submitted for publication. [47] X. Xu, Global regularity of solutions of D Boussinesq equations with fractional diffusion, Nonlinear Analysis: TMA 7 (010), School of Mathematical Sciences, Capital Normal University, Beijing , People s Republic of China address: jiuqs@mail.cnu.edu.cn Institute of Applied Physics and Computational Mathematics, P.O. Box 8009, Beijing , People s Republic of China address: miao changxing@iapcm.ac.cn 3 Department of Mathematics, Oklahoma State University, 401 Mathematical Sciences, Stillwater, OK 74078, USA address: jiahong@math.okstate.edu 4 School of Mathematical Sciences, Peking University, Beijing , People s Republic of China address: zfzhang@math.pku.edu.cn