\documentclass[reqno]{amsart}
\usepackage{hyperref}
\usepackage{amssymb}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2016 (2016), No. 73, pp. 1--11.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2016 Texas State University.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2016/73\hfil 2D Boussinesq equations]
{Blow-up criterion for the 2D Euler-Boussinesq system in terms of temperature}

\author[C. Qian \hfil EJDE-2016/73\hfilneg]
{Chenyin Qian}

\address{Chenyin Qian \newline
Department of Mathematics, 
Zhejiang Normal University, 
Jinhua 321004, China}
\email{qcyjcsx@163.com}

\thanks{Submitted December 18, 20915. Published March 15, 2016.}
\subjclass[2010]{35Q35, 35B35, 35B65, 76D05}
\keywords{2D Boussinesq equation; blow-up criterion;  Besov space; 
\hfill\break\indent transport equation}

\begin{abstract}
 In this article, we study the blow-up slutions for the 2D Euler-Boussinesq
 equation. In particular, it is shown that if
 $$
 \int_0^{T^*} \sup_{r\geq 2}\frac{\|\Lambda^{1-\alpha}
 \theta(t)\|_{L^{r}}} {\sqrt{r\log r}}\,\mathrm{d}t<\infty \quad \text{or}\quad
 \int_0^{T^*} \|\Lambda^{1-\alpha} \theta\|_{\dot{B}^0_{\infty,\infty}}\,\mathrm{d}t
 <\infty,
 $$
 then the local solution can be continued to the global one.
 This is an improvement of  classical Lipschitz-type blow-up criterion
 ($\|\nabla\theta\|_{L^1_tL^{\infty}}$) in terms of the
 temperature  $\theta$.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\allowdisplaybreaks

\section{Introduction}

The 2D incompressible generalized
Boussinesq equations with  the fractional Laplacian dissipation is of the type
\begin{equation} \label{e1.1} 
\begin{gathered}
    \partial_t v+(v\cdot\nabla)v+\nu\Lambda^{\beta}v+\nabla \pi
=\theta \mathrm{e}_2,\quad (x, t)\in\mathbb{R}^2\times(0,\infty),\\
  \partial_t  \theta +(v\cdot\nabla)\theta+\kappa\Lambda^{\alpha}\theta=0,\\
 \nabla\cdot v=0,\\
 v(x, 0)=v^0, \quad \theta(x, 0)=\theta^0,
\end{gathered}
\end{equation}
 where $v=v(x,t)=(v_1(x,t), v_2(x,t))$, $\pi=\pi(x,t)$ and $\theta=\theta(x,t)$
 stand for, respectively, the velocity vector
field,   the pressure   and the temperature. Here, the constants
$\nu\geq 0$ and $\kappa\geq0$ denote viscosity coefficient and
thermal diffusivity coefficient respectively, and $\mathrm{e}_2 = (0, 1)$.
$\Lambda=\sqrt{-\Delta}$ is the Zygmund operator, and
$\Lambda^{\alpha}$ is defined by the Fourier transform,
$$
\widehat{\Lambda^{\alpha}f}(\xi)=|\xi|^{\alpha}\widehat{f}(\xi),\quad
\widehat{f}(\xi)=\int_{\mathbb{R}^2}e^{-ix\cdot\xi}f(x)\,\mathrm{d}x.
$$
The study of the standard 2D incompressible Boussinesq system
(with $\nu>0, \kappa>0$ and $\alpha=\beta=2$ in \eqref{e1.1})
can be traced to 1980s (see \cite{[6]}). Later, there are many works
considering the global well-posedness problem for the standard 2D Boussinesq
system  without the viscosity ($\nu=0$, $\kappa>0$) or without
the thermal diffusivity ($\nu>0$, $\kappa=0$), see
\cite{[11],[8],[21],[23],[22]}.

 Recently, the 2D incompressible Boussinesq system with fractional Laplacian
generalizations have  attracted considerable attention. For example,
 Hmidi and Zerguine  obtained the global well-posedness of Euler-Boussinesq
system \eqref{e1.1} with $1<\alpha\leq2$ in \cite{[13]}.
 Hmidi, Keraani and Rousset  obtained the global
well-posedness  for the Euler-Boussinesq system with critical
dissipation (namely, $\alpha=1$) in \cite{[27]}.
There are many other related results to
the Boussinesq equations system \eqref{e1.1}, we refer the reader to
\cite{[16],[17],[81],[19],[15],[14],[82],[83],[18]}.
 From above mentioned result, we see that it is difficult for us to get
the global well-posedness for the system  \eqref{e1.1} with super-critical
dissipation ($\alpha+\beta<1$). Therefore, people may turn to the  blow-up
criteria in terms of velocity ($v$) and temperature ($\theta$).
As for the blow-up  criterion for Boussinesq equations in terms of $\theta$,
 we refer readers  to \cite{[24]} (see also \cite{[8]}), in
which the authors  Chae and  Nam obtained the blow-up criterion
$\|\nabla\theta\|_{L^1_tL^{\infty}}$ in the framework of
$L^2$.  The purposes of this article  is to
 establish  some  blow-up criteria better than $\|\nabla\theta\|_{L^1_tL^{\infty}}$.


  In view of  \cite{[13]}, one can  establish the local
well-posedness results for the system \eqref{e1.1} with
 $ 0<\alpha\leq 2, \kappa>0 $
and
\begin{equation} \label{C2}
 v^0\in B^{s}_{p,1}(\mathbb{R}^2)\quad \text{with}\quad
\operatorname{div}v^0=0, \quad \theta^0\in B^{s-\alpha}_{p,1}(\mathbb{R}^2),
\end{equation}
where $s\geq1+\frac{2}{p}$ with $p\in]1,\infty[$. Here, we define the
function space of solutions as follows
\begin{equation}\label{dsanz1}
\mathcal{X}_{T}^{s,p}:=
\mathcal{C}([0,T];B^{s}_{p,1}(\mathbb{R}^2))
\times \left(\mathcal{C}([0,T];B^{s-\alpha}_{p,1}(\mathbb{R}^2))\cap
L^{1}([0,T];{B}^{s}_{p,1}(\mathbb{R}^2)\right),
\end{equation}
where $B^{s-\alpha}_{p,1}(\mathbb{R}^2)$ and $B^{s}_{p,1}(\mathbb{R}^2)$
are Besov spaces (see Section 2).
Now, we give the main results of this article.

\begin{theorem}\label{thm1.1}
Let $1/2<\alpha\leq1$  and $(v,\theta)\in \chi_{T^*}^{s,p} $
be the local unique solution of \eqref{e1.1} with  initial data
satisfying \eqref{C2}, where $T^*$ is  the maximal existence time.
If
\begin{equation}\label{4} 
\int_0^{T^*} \sup_{r\geq 2}\frac{\|\Lambda^{1-\alpha} \theta(t)\|_{L^{r}}} 
{\sqrt{r\log r}}\,\mathrm{d}t<\infty,
\end{equation}
 or
\begin{equation}\label{5} 
\int_0^{T^*} \|\Lambda^{1-\alpha}
\theta\|_{\dot{B}^0_{\infty,\infty}}\,\mathrm{d}t<\infty,
\end{equation}
then the solution $(v,\theta)$ can be continued beyond $T^*$.
\end{theorem}

\section{Notation and preliminaries}

 We begin this section with dyadic decomposition. 
Let $\mathcal {S}$ be the Schwartz class of rapidly  decreasing functions. Let
functions $\chi, \varphi \in\mathcal {S}(\mathbb{R}^d)$ supported in
$\mathfrak{B}=\{ \xi\in\mathbb{R}^d: |\xi|\leq 4/3\}$ and
$\mathfrak{C}=\{\xi\in\mathbb{R}^d: 3/4\leq|\xi|\leq 8/3\}$
respectively, such that
\begin{gather*}
\sum_{j\in\mathbb{Z}}\varphi(2^{-j}\xi)=1, \forall
\xi\in\mathbb{R}^d\backslash\{0\}, \\
\chi(\xi)+\sum_{j\geq 0}\varphi(2^{-j}\xi)=1, \forall
\xi\in\mathbb{R}^d.
\end{gather*}
For $u\in\mathcal {S}^{\prime}$, we set
\begin{gather*}
\varDelta_{-1}u=\chi(D)u \quad  \forall q\in\mathbb{N},\quad
\varDelta_{q}u=\varphi(2^{-q}D)u \quad \forall q\in\mathbb{\mathbb{Z}},\\
\dot{\varDelta}_{q}u=\varphi(2^{-q}D)u.
\end{gather*}
The following low-frequency cut-off will be also used:
$$
{\mathcal {S}}_q u=\sum_{-1\leq j\leq q-1}\varDelta_{j}u 
\quad\text{and}\quad \dot{\mathcal {S}}_q u=\sum_{j\leq
q-1}\dot{\varDelta}_{j}u
$$ 
We now recall the definitions of Besov
spaces through the dyadic decomposition.
 Let $s\in \mathbb{R}$, $p, q \in [1,\infty]$, the inhomogeneous Besov space
${B}_{p,q}^{s}(\mathbb{R}^d)$ is the set of tempered distribution
$u$ such that
$$
\|u\|_{{B}_{p,q}^{s}}
:=\left(2^{js}\|\varDelta_j u\|_{L^{p}}\right)_{\ell^{q}}<\infty.
$$
To define the homogeneous Besov spaces, we first denote by 
$\mathcal {S}^{\prime}/\mathcal {P}$ the space of tempered distributions
modulo polynomials. Thus we define the space
$\dot{B}_{p,q}^{s}(\mathbb{R}^d)$ as  the set of  distribution
$u\in\mathcal {S}^{\prime}/\mathcal {P}$ such that
$$
\|u\|_{\dot{B}_{p,q}^{s}}:=\big(2^{js}\|\dot{\varDelta}_j
u\|_{L^{p}}\big)_{\ell^{q}}<\infty.
$$
We point out that if  $s>0$ then we have
${B}_{p,q}^{s}(\mathbb{R}^d)=\dot{B}_{p,q}^{s}(\mathbb{R}^d)\cap
L^{p}(\mathbb{R}^d)$ and
$$
\|u\|_{{B}_{p,q}^{s}}\approx\|u\|_{\dot{B}_{p,q}^{s}}+\|u\|_{L^{p}}.
$$
In our next study we require two kinds of coupled space-time Besov
spaces. The first is defined in the following manner: for $T>0$ and
$q\geq1$, we denote by $L^{r}_{T}\dot{B}_{p,q}^{s}$ the set of all
tempered distribution u satisfying
$$
\|u\|_{L^{r}_{T}\dot{B}_{p,q}^{s}}:=\|\Big(2^{qs}\|\dot{\varDelta}_q
u\|_{L^{p}}\Big)_{\ell^{q}}\|_{L^{r}_{T}}<\infty.
$$
The second mixed space is $\widetilde{L}^{r}_{T}\dot{B}_{p,q}^{s}$
which is the set of tempered distribution u satisfying
$$
\|u\|_{\widetilde{L}^{r}_{T}\dot{B}_{p,q}^{s}}
:=\Big(2^{qs}\|\dot{\varDelta}_q
u\|_{L^{r}_{T}L^{p}}\Big)_{\ell^{q}}<\infty.
$$



\section{Proof of main results}

In this section we use $\Phi_k$ to denote function of the form
$$
\Phi_k(t)=C_0\underbrace{\exp(\dots\exp}_{k\text{ times}}(C_0t)\dots),
$$
where $C_0$ depends on the involved norms of the initial data and
its value may vary from line to line up to some absolute constants.
We will make an intensive use (without mentioning it) of the
following trivial facts
$$
\int_0^t\Phi_k(\tau)d\tau\leq\Phi_k(t)\quad \text{and}\quad
\exp\Big(\int_0^t\Phi_k(\tau)d\tau\Big)\leq \Phi_{k+1}(t).
$$ 
Firstly, we introduce a pseudo-differential operator
$\mathcal {R}_{\alpha}$ defined by 
$\mathcal {R}_{\alpha}:=\Lambda^{-\alpha}\partial_1=\Lambda^{1-\alpha}\mathcal {R},
0<\alpha<1$, where $\mathcal {R}:=\frac{\partial_1}{\Lambda}$ is the
usual Riesz transform. For \eqref{e1.1}, the vorticity  equation is
\begin{equation}\label{4.1}
\partial_t \omega+v\cdot\nabla\omega=\partial_1\theta,
\end{equation}
and the acting of $\mathcal {R}_{\alpha}$ on the temperature
equation, we have
\begin{equation}\label{4.2}
\partial_t \mathcal {R}_{\alpha}\theta+v\cdot
\nabla\mathcal {R}_{\alpha}\theta+\kappa\Lambda^{\alpha}\mathcal {R}_{\alpha}\theta
=-[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta.
\end{equation}
Denote $\digamma=\omega+\mathcal {R}_{\alpha}\theta$. Thus we obtain
\begin{equation}\label{4.3}
\partial_t \digamma+v\cdot\nabla\digamma
=-[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta.
\end{equation}
Firstly, we give the following  crucial  Lemmas  which are  useful
for us to  proof the main results.

\begin{lemma}[{\cite{[38]}}] \label{lem2.2}
Assume that $v$ is divergence-free and that $f$ satisfies the
transport equation on $\mathbb{R}^d$,
\begin{equation} \label{e3.4}
\begin{gathered}
 \partial_t f+v\cdot\nabla f=g, \\
 f|_{t=0}=f_0. 
\end{gathered}
\end{equation} 
There exists a constant $C$, depending only
on $d$, such that for all $1\leq p,r\leq\infty$ and
$t\in\mathbb{R}_{+}$, we have
\begin{equation}\label{24}
\|f\|_{\widetilde{L}^{\infty}_{t}(B^0_{p,r})}\leq
C\big(\|f_0\|_{B^0_{p,r}}+\|g\|_{\widetilde{L}^{1}_{t}(B^0_{p,r})}\big)
\Big(1+\int_{0}^t\|\nabla v(\tau)\|_{L^{\infty}}d\tau\Big).
\end{equation}
\end{lemma}

\begin{lemma}[\cite{[13]}] \label{lem2.4}
Let $v$ be a solution of the incompressible Euler system on $\mathbb{R}^2$
\begin{equation} \label{28}
\begin{gathered}
 \partial_t v+v\cdot\nabla v+\nabla\pi=f, \\
 v(x, 0)=v^0,\\
 \operatorname{div}v=0.
\end{gathered}
\end{equation}
 Then for $s>-1$, $(p,r)\in(1,\infty)\times[1,\infty]$ we have
\begin{equation} \label{29}
 \|v\|_{\widetilde{L}^{\infty}_{t}(B^{s}_{p,r})}\leq
 C\mathrm{e}^{C V(t)}\Big(\|v^0\|_{B^{s}_{p,r}}+  \int_{0}^{t}
 \mathrm{e}^{-C V(\tau)}\|f(\tau)\|_{B^{s}_{p,r}}d\tau\Big),
\end{equation}
with $ V(t):=\int_0^t\|\nabla v(\tau)\|_{L^{\infty}}d\tau$.
\end{lemma}

\begin{lemma}\label{lem4.1}
Let $\alpha\in(0,1)$, $v$ be a smooth divergence-free vector
field.
\begin{itemize}
\item[(i)] for every $(s,p,r)\in
(-1,\alpha)\times[2,\infty]\times[1,\infty]$, there exists a
constant $C>0$ such that
\begin{equation}\label{4.3l}
\|[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta\|_{B^{s}_{p,r}}
\leq C \|\nabla v\|_{L^{p}}\left(\|\theta\|_{B^{s+1-\alpha}_{\infty,r}}
+\|\theta\|_{L^p}\right),
\end{equation}

\item[(ii)] for every $(r,\varrho)\in[1,\infty]\times (1,\infty)$ and
$\epsilon>0 $, there exists a constant $C>0$ such that
\begin{equation}\label{4.3ll}
\|[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta\|_{B^0_{\infty,r}}
\leq C(\|\omega\|_{L^{\infty}}+\|\omega\|_{L^{\varrho}})
\left(\|\theta\|_{B^{\epsilon+1-\alpha}_{\infty,r}}+\|\theta\|_{L^\varrho}\right).
\end{equation}
\end{itemize}
\end{lemma}

\begin{proof}
Part (i) can be found in \cite[Proposition 3.3]{[18]}. For  the
second part, we can imitate the proof of 
\cite[Theorem 3.3(2)]{[27]} to get (ii). We omit the detail  here.
\end{proof}

\begin{lemma}\label{lem2.3}
Let  $\alpha \in[0,2]$, $\kappa>0$ and $v$ be a smooth
divergence-free vector-field of $\mathbb{R}^2$. Let $\theta$ be a
smooth solution of
\begin{equation} \label{26}
\begin{gathered}
 \partial_t\theta+v\cdot\nabla\theta+\kappa\Lambda^{\alpha}\theta=f, \\
 \theta(x, 0)=\theta^0.
\end{gathered}
 \end{equation} 
Then 
(i)  for every $\rho\in[1,\infty]$, $p\in[1,\infty]$, $s>-1$  one
has (see \cite{[13]})
\begin{equation} \label{27}
 \|\theta\|_{\widetilde{L}^{\rho}_{t}(B_{p,1}^{s+\frac{\alpha}{\rho}})}\leq
 C\mathrm{e}^{C V(t)}\Big(\|\theta^0\|_{B^{s}_{p,1}}(1+t^{\frac{1}{\rho}})+
 \int_{0}^{t}
 \mathcal {T}_{s}(\tau)d\tau+\|f\|_{L^{1}_{t}(B^{s}_{p,1})}\Big)
 \end{equation}
with
$$
V(t):=\int_0^t\|\nabla v(\tau)\|_{L^{\infty}}d\tau,\quad
\mathcal {T}_{s}(t):=\|\nabla
\theta(t)\|_{L^{\infty}}\|v(t)\|_{B^{s}_{p,1}}1_{[1,\infty)}(s).
$$

(ii) for every $p\in[1,\infty]$, we have the $L^p$ estimates (see
\cite{[069]})
\begin{equation} \label{27j}
 \|\theta\|_{L^p}
\leq  \|\theta^0\|_{L^{p}}+  \int_{0}^{t}  \|f(\tau)\|_{L^p}d\tau.
 \end{equation}

(iii) for every $p\in[1,\infty]$ (see \cite{[27]})
\begin{equation} \label{2g7}
 \|\theta\|_{\widetilde{L}^{\infty}_{t}B_{p,1}^0}
\leq  C\big(\|\theta^0\|_{B^0_{p,1}}+\|f\|_{L^{1}_{t}B^0_{p,1}}\big)
\Big(1+\int_0^t\|\nabla  v(\tau)\|_{L^{\infty}}d\tau\Big).
\end{equation}

(iv) for $p\in(1,\infty)$, $\rho\in[1,\infty]$ and $f=0$, there
exists a constant $C$ such that (see \cite{[13]})
\begin{equation} \label{9jj1}
 \sup_{q\in \mathbb{N}}2^{q\frac{\alpha}{\rho}}
\|\varDelta_q\theta\|_{L^{\rho}_{t}L^{p}}
\leq C  \|\theta^0\|_{L^{p}}+C  \|\theta^0\|_{L^{\infty}}
\|\omega\|_{L^{1}_{t}L^{p}}.
 \end{equation}
\end{lemma}

\begin{proof}[Proof of Theorem \ref{thm1.1}]
By using the above Lemmas and the equation \eqref{4.3}, we  first
estimate $\|\nabla v(t)\|_{L^{\infty}}$. We do this in several steps.
\smallskip

\noindent\textbf{Step 1:} Estimation of $\|\omega(t)\|_{L^{a}}$, 
for some $\max\{2,p\}<a<\infty$.
From \eqref{4.3}, for  any $t<T^*$, by Lemma \ref{lem2.3} (ii) we have
\begin{equation}\label{4.4}
\|\digamma(t)\|_{L^{a}}\leq\|\digamma(0)\|_{L^{a}}
+\int_0^t\|[\mathcal {R}_{\alpha},
v\cdot\nabla]\theta(\tau)\|_{L^{a}}d\tau.
\end{equation}
Note that (see \cite[p 516]{[069]})
$\int_0^t\Lambda^{\alpha}\mathcal {R}_{\alpha}\theta|\mathcal
{R}_{\alpha}\theta|^{a-2}\mathcal {R}_{\alpha}\theta d\tau\geq 0$,
by \eqref{4.2} we have
\begin{equation}\label{4.5}
\|\mathcal {R}_{\alpha}\theta(t)\|_{L^{a}}\leq\|\mathcal
{R}_{\alpha}\theta^0\|_{L^{a}} +\int_0^t\|[\mathcal
{R}_{\alpha}, v\cdot\nabla]\theta(\tau)\|_{L^{a}}d\tau.
\end{equation}
Since $v^0\in B^{s}_{p,1}$, $s\geq 1+2/p$, we see that
$\omega^0\in L^{p}\cap L^{\infty}$, and for  
$\theta^0\in B^{s-\alpha}_{p,1}$ with $0<\alpha\leq1$,  we have 
$\theta^0\in L^{p}\cap L^{\infty}$ note that $p<a<\infty$, the embedding
$B^{s-\alpha}_{p,1}\hookrightarrow
B^{s-\alpha-\frac{2}{p}+\frac{2}{a}}_{a,1}\hookrightarrow
B^{1-\alpha+\frac{2}{a}}_{a,1}\hookrightarrow B^{1-\alpha}_{a,1}$,
we  have $\mathcal {R}_{\alpha}\theta^0\in L^{a}$.  From \eqref{4.4}
and \eqref{4.5} we obtain
\begin{equation}\label{4.6}
\begin{aligned}
\|\omega\|_{L^a}
&\leq\|\digamma(t)\|_{L^a}+\|\mathcal
{R}_{\alpha}\theta(t)\|_{L^a}\\
&\leq C\big(\|\omega^0\|_{L^{p}\cap
L^{\infty}}+\|\theta^0\|_{B^{s-\alpha}_{p,1}}\big)
+C\int_0^t\|[\mathcal {R}_{\alpha},
v\cdot\nabla]\theta(\tau)\|_{L^{a}}d\tau.
\end{aligned}
\end{equation}
Using the classical embedding $B^0_{a,1}\hookrightarrow L^{a}$,
and Lemma \ref{lem4.1} (i), we have
\begin{equation}\label{4.7}
\begin{aligned}
\|[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta(t)\|_{L^{a}}
&\leq C\|[\mathcal {R}_{\alpha},
v\cdot\nabla]\theta(t)\|_{B^0_{a,1}}\\
&\leq C \|\nabla v(t)\|_{L^a}
\big(\|\theta(t)\|_{B^{1-\alpha}_{\infty,1}}+\|\theta(t)\|_{L^a}\big).
\end{aligned}
\end{equation}
From  $\|\theta(t)\|_{L^a}\leq \|\theta^0\|_{L^a},
\forall t\geq 0$, it follows that
\begin{equation}\label{4.8}
\begin{aligned}
\|\omega(t)\|_{L^a}
&\leq C\left(\|\omega^0\|_{L^{p}\cap
L^{\infty}}+\|\theta^0\|_{B^{s-\alpha}_{p,1}}\right) \\
&\quad +C\int_0^t\|\omega(\tau)\|_{L^a}
\left(\|\theta(\tau)\|_{B^{1-\alpha}_{\infty,1}}+\|\theta^0\|_{L^a}\right)d\tau.
\end{aligned}
\end{equation}
According to Gronwall's inequality we obtain
\begin{equation}\label{4.9}
\|\omega\|_{L^a}\leq
Ce^{Ct}e^{C\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}}.
\end{equation}
Next we estimate $\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}$,
let $N\in\mathbb{N}$, by the Littlewood-Paley decomposition, by
condition \eqref{4} we see that
%\label{4.10b}
\begin{align*}
&\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}} \\
&\leq \|\mathcal{S}_N\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}+\sum_{q\geq N}
2^{q(1-\alpha)}\|\varDelta_q\theta\|_{L^{1}_{t}L^{\infty}}\\
&\leq Ct\|\theta^0\|_{L^{\infty}}+C\int_0^t\sum_{0\leq q< N}
\|\varDelta_q\Lambda^{1-\alpha}\theta(\tau)\|_{L^{\infty}}d\tau
 +\sum_{q\geq N} 2^{q(1-\alpha)}\|\varDelta_q\theta\|_{L^{1}_{t}L^{\infty}}
 \\
&\leq Ct\|\theta^0\|_{L^{\infty}}+C\int_0^t\sum_{0\leq q< N}
2^{q\frac{2}{b}}\|\varDelta_q\Lambda^{1-\alpha}\theta(\tau)\|_{L^{b}}d\tau
+\sum_{q\geq N} 2^{q(1-\alpha)}\|\varDelta_q\theta\|_{L^{1}_{t}L^{\infty}}
 \\
&\leq C\Big(\sum_{0\leq q< N}2^{q\frac{2}{b}}\Big)\sqrt{b\log b}
\int_0^t\sup_{r\geq2}
\frac{\|\Lambda^{1-\alpha}\theta(\tau)\|_{L^{r}}} {\sqrt{r\log r}}d\tau
\\
&\quad  +Ct\|\theta^0\|_{L^{\infty}}+C\sum_{q\geq N}
2^{q(1-\alpha+\frac{2}{a})}\|\varDelta_q\theta\|_{L^{1}_{t}L^{a}}
 \\
&\leq Ct\|\theta^0\|_{L^{\infty}}+C2^{N\frac{2}{b}}
\sqrt{b\log b}+C\sum_{q\geq N}
2^{q(1-\alpha+\frac{2}{a})}\|\varDelta_q\theta\|_{L^{1}_{t}L^{a}}.
\end{align*}
Since $\alpha>1/2$, we choose $a$ large enough such that
$1-2\alpha+4/a<0$, and using Lemma \ref{lem2.3} (iv), we have
\begin{equation}\label{4.11}
\begin{aligned}
&\sum_{q\geq N} 2^{q(1-\alpha+\frac{2}{a})}\|\varDelta_q\theta\|_{L^{1}_{t}L^{a}}\\
& \leq \sum_{q\geq N}
2^{q(1-2\alpha+\frac{2}{a})}\Big(\|\theta^0\|_{L^a}
+\|\theta^0\|_{L^\infty}\int_0^t\|\omega(\tau)\|_{L^a}d\tau\Big)\\
&\leq C\|\theta^0\|_{L^{a}}+
2^{N(1-2\alpha+\frac{2}{a})}\|\theta^0\|_{L^\infty}
\int_0^t\|\omega(\tau)\|_{L^a}d\tau.
\end{aligned}
\end{equation}
 Therefore, 
\begin{align*}
\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}
&\leq C(t+1)\|\theta^0\|_{L^{\infty}\cap
L^{a}}+C2^{N\frac{2}{b}}b^{\frac{1}{2}+\epsilon}\\
&\quad +2^{N(1-2\alpha+\frac{2}{a})}\|\theta^0\|_{L^\infty}
\int_0^t\|\omega(\tau)\|_{L^a}d\tau,
\end{align*}
for any $0<\epsilon<1/2$. Setting $b=N$ and selecting $N$ as follows
$$
N=\Big[\frac{\log\big(e+\int_0^t\|\omega(\tau)\|_{L^a}d\tau\big)}
{(2\alpha-1-2/a)\log2}\Big]+2.
$$
Then 
\begin{equation}\label{4.12}
\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}\leq
C(t+1)\|\theta^0\|_{L^{\infty}\cap L^{a}}
+C\Big[\log\Big(e+\int_0^t\|\omega(\tau)\|_{L^a}d\tau\Big)
\Big]^{\frac{1}{2}+\epsilon}.
\end{equation}
Combining \eqref{4.9}  with \eqref{4.12}, we have
\begin{equation}\label{4.13}
\begin{aligned}
\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}
&\leq C(t+1)\|\theta^0\|_{L^{\infty}\cap L^{a}}
+C\Big[\log\Big(e+\int_0^t\|\omega(\tau)\|_{L^a}d\tau\Big)
 \Big]^{\frac{1}{2}+\epsilon}\\
&\leq C(t+1)\|\theta^0\|_{L^{\infty}\cap L^{a}}
+C\Big[\log\Big(e+tCe^{Ct}e^{C\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}}\Big)
\Big]^{\frac{1}{2}+\epsilon}\\
&\leq C(t^2+t+1)\left(\|\theta^0\|_{L^{\infty}\cap L^{a}}+1\right)+
C\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}^{\frac{1}{2}+\epsilon},
\end{aligned}
\end{equation}
it follows that
$$
\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}\leq C(1+t+t^2).
$$
By \eqref{4.9}, we obtain
\begin{equation}\label{4.14}
\|\omega(t)\|_{L^a}\leq \Phi_1(t),
\end{equation}
where $\Phi_k(t), k=1,2,3,\dots$ is same the as as in subsection 3.3.
\smallskip

\noindent\textbf{Step 2:}
 Estimation of $\|\omega(t)\|_{L^{\infty}}$.
By using the maximum principle for the transport equation
\eqref{4.3} we have
\begin{equation}\label{4.17}
\|\digamma(t)\|_{L^{\infty}}\leq\|\digamma(0)\|_{L^{\infty}}
+\int_0^t\|[\mathcal {R}_{\alpha},
v\cdot\nabla]\theta(\tau)\|_{L^{\infty}}d\tau,
\end{equation}
By using Lemma \ref{lem2.3} (ii) with $p=\infty$, we see that
\eqref{4.2} follows that
\begin{equation}\label{4.18}
\|\mathcal {R}_{\alpha}\theta(t)\|_{L^{\infty}}\leq\|\mathcal
{R}_{\alpha}\theta^0\|_{L^{\infty}} +\int_0^t\|[\mathcal
{R}_{\alpha}, v\cdot\nabla]\theta(\tau)\|_{L^{\infty}}d\tau.
\end{equation}
For  $\theta^0\in B^{s-\alpha}_{p,1}$ with $0<\alpha\leq1$, 
$s\geq 1+2/p$, we have  $B^{s-\alpha}_{p,1}\hookrightarrow
B^{s-\alpha-\frac{2}{p}}_{\infty,1}\hookrightarrow
B^{1-\alpha}_{\infty,1}$, we  have $\mathcal {R}_{\alpha}\theta^0\in
L^{\infty}$.  From \eqref{4.17} and \eqref{4.18} we obtain
\begin{equation}\label{4.19}
\begin{aligned}
\|\omega\|_{L^\infty}
&\leq\|\digamma(t)\|_{L^\infty}+\|\mathcal
{R}_{\alpha}\theta(t)\|_{L^\infty}\\
&\leq\|\omega^0\|_{L^{\infty}}+\|\theta^0\|_{B^{s-\alpha}_{p,1}}
+2\int_0^t\|[\mathcal {R}_{\alpha},
v\cdot\nabla]\theta(\tau)\|_{L^{\infty}}d\tau.
\end{aligned}
\end{equation}
Using the classical  embedding $B^0_{\infty,1}\hookrightarrow
L^{\infty}$, and Lemma \ref{lem4.1} (i), we have
\begin{equation}\label{4.20}
\begin{aligned}
\|[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta\|_{L^{\infty}}
&\leq C\|[\mathcal {R}_{\alpha},
v\cdot\nabla]\theta\|_{B^0_{\infty,1}}\\
&\leq C(\|\omega\|_{L^{\infty}}+\|\omega\|_{L^{a}})
\left(\|\theta\|_{B^{\epsilon+1-\alpha}_{\infty,1}}+\|\theta\|_{L^a}\right).
\end{aligned}
\end{equation}
Combining \eqref{4.19} and \eqref{4.20} we have
\begin{equation}\label{4.21}
\begin{aligned}
&\|\omega(t)\|_{L^\infty} \\
&\leq C+C\int_0^t(\|\omega(\tau)\|_{L^{\infty}}+\|\omega(\tau)\|_{L^{a}})
\left(\|\theta(\tau)\|_{B^{\epsilon+1-\alpha}_{\infty,1}}
 +\|\theta(\tau)\|_{L^a}\right)d\tau\\
&\leq C+\|\omega\|_{L^{\infty}_tL^{a}}
 \left(\|\theta\|_{L^{1}_{t}B^{\epsilon+1-\alpha}_{\infty,1}}+t\|\theta^0\|_{L^a}
 \right)\\
&\quad +C\int_0^t\|\omega(\tau)\|_{L^{\infty}}
\left(\|\theta(\tau)\|_{B^{\epsilon+1-\alpha}_{\infty,1}}
+\|\theta^0\|_{L^a}\right)d\tau.
\end{aligned}
\end{equation}
We claim that
\begin{equation}\label{4.22}
\|\theta\|_{L^{1}_tB^{\epsilon+1-\alpha}_{\infty,1}}\leq
\Phi_1(t),
\end{equation} 
for some suitable $\epsilon>0$.
In fact that, from step 1  and the Lemma \ref{lem2.3} (iv)  we have
$$
\|\theta\|_{\widetilde{L}^{1}_tB^{\alpha}_{a,\infty}}\leq \Phi_1(t),
$$
and then  for every $\sigma<\alpha$, we have
$$
\|\theta\|_{{L}^{1}_tB^{\sigma}_{a,1}}
\leq\|\theta\|_{\widetilde{L}^{1}_tB^{\alpha}_{a,\infty}}
\leq \Phi_1(t).
$$
We  choose that $\sigma=\alpha-1/a$, and we keep in mind that $a$
satisfies $1-2\alpha+4/a<0$, therefore $2\alpha-1-4/a>0$, we choose
$\epsilon$ satisfies $0<\epsilon<2\alpha-1-3/a$, then we have
$\epsilon+1-\alpha<\alpha-3/a=\sigma-2/a$ (here we note that the
selected parameter $a$ is make sure $\alpha>3/a$, and we have
$\sigma-2/a>0$). Therefore, we have the embedding
$B^{\alpha-1/a}_{a,1}=B^{\sigma}_{a,1}\hookrightarrow
B^{\sigma-2/a}_{\infty,1}\hookrightarrow
B^{\epsilon+1-\alpha}_{\infty,1}$, and we finally obtain \eqref{4.22}.
Thus, by \eqref{4.22} and step 1,  one has
\begin{equation}\label{4.23}
\|\omega(t)\|_{L^\infty}\leq
\Phi_{1}(t)+C\int_0^t\|\omega(\tau)\|_{L^{\infty}}
\left(\|\theta(\tau)\|_{B^{\epsilon+1-\alpha}_{\infty,1}}
+\|\theta^0\|_{L^a}\right)d\tau.
\end{equation}
Again using \eqref{4.22} and Gronwall's inequality, we obtain
$$
\|\omega(t)\|_{L^\infty}\leq\Phi_2(t).
$$
\smallskip

\noindent\textbf{Step 3:} Estimation of $\|\nabla v(t)\|_{L^{\infty}}$.
 By using Lemma \ref{lem2.2} and Lemma \ref{lem2.3}(iii),  from \eqref{4.2} 
and \eqref{4.3}  we have
 \begin{equation}\label{4.27}
\|\digamma(t)\|_{B^0_{\infty,1}}+\|\mathcal {R}_{\alpha}
\theta(t)\|_{B^0_{\infty,1}}\leq\left(C+\|[\mathcal {R}_{\alpha},
v\cdot\nabla]\theta\|_{L^{1}_{t}B^0_{\infty,1}}\right) (1+\|\nabla
v\|_{L^{1}_{t}L^{\infty}}),
\end{equation}
Thanks to Lemma \ref{lem4.1}, and by step 1, step 2, we have
\begin{equation}\label{4.28}
\begin{aligned}
&\|[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta\|_{L^{1}_{t}B^0_{\infty,1}}\\
&\leq C\int_{0}^t (\|\omega(\tau)\|_{L^{\infty}}+\|\omega(\tau)\|_{L^{a}})
\left(\|\theta(\tau)\|_{B^{\epsilon+1-\alpha}_{\infty,1}}
+\|\theta(\tau)\|_{L^a}\right)d\tau
\leq \Phi_2(t).
\end{aligned}
\end{equation}
Therefore, we have
$$
\|\omega(t)\|_{B^0_{\infty,1}}\leq\|\digamma(t)\|_{B^0_{\infty,1}}+\|\mathcal
{R}_{\alpha} \theta(t)\|_{B^0_{\infty,1}}\leq\Phi_2(t)
\left(1+\|\nabla v\|_{L^{1}_{t}L^{\infty}}\right).
$$
On the other hand, we have
\begin{equation}\label{4.29}
\begin{aligned}
\|\nabla v(t)\|_{L^{\infty}}
& \leq \|\nabla \varDelta_{-1}v(t)\|_{L^{\infty}}+\sum_{q\in\mathbb{N}}\|
\varDelta_{q}\nabla v(t)\|_{L^{\infty}}\\
&\leq \|\omega(t)\|_{L^a}+\|\omega(t)\|_{B^0_{\infty,1}}\\
&\leq \Phi_1(t)+\|\omega(t)\|_{B^0_{\infty,1}}.
\end{aligned}
\end{equation}
Using \eqref{4.29}, we have
$$
\|\omega(t)\|_{B^0_{\infty,1}}\leq\Phi_2(t)
\Big(1+\int_0^t\|\omega(\tau)\|_{B^0_{\infty,1}}d\tau\Big).
$$
From Gronwall's inequality we obtain
$\|\omega(t)\|_{B^0_{\infty,1}}\leq \Phi_3(t)$.
Going back to \eqref{4.29}, we obtain
$$
\|\nabla v\|_{L^{\infty}}\leq \Phi_3(t).
$$

 Next, we give the estimate of $\|
v(t)\|_{{B}^{s}_{p,1}}, \| \theta(t)\|_{{B}^{s-\alpha}_{p,1}}$.
Since $s\geq1+2/p$,  we have $B^{s-\alpha}_{p,1}\hookrightarrow
B^{1-\alpha}_{\infty,1}$. By Lemma \ref{lem2.3}(i), we have
\begin{equation} \label{2j7jj}
\begin{aligned}
\|\theta(t)\|_{\widetilde{L}^{\infty}_{t}({B}^{1-\alpha}_{\infty,1})}
+\int_0^t\| \theta(\tau)\|_{{B}^{1}_{\infty,1}}d\tau
&\leq  C\mathrm{e}^{C V(t)}(1+t)\|\theta^0\|_{B^{1-\alpha}_{\infty,1}}\\
&\leq  C\mathrm{e}^{C V(t)}(1+t)\|\theta^0\|_{B^{s-\alpha}_{p,1}}.
 \end{aligned}
\end{equation}
On the other hand, by Lemma \ref{lem2.4} we have
$$
\|
v(t)\|_{\widetilde{L}^{\infty}_{t}({B}^{s-\alpha}_{p,1})}\leq
Ce^{CV(t)}\Big( \|
v^0\|_{{B}^{s}_{p,1}}+\int_0^t\|\theta(\tau)\|_{{B}^{s-\alpha}_{p,1}}d\tau\Big).
$$
Therefore, by the embedding $B^0_{\infty,1}\hookrightarrow
L^{\infty}$ we have
\begin{equation} \label{2j7j8j}
\begin{aligned}
&\int_{0}^{t}
 \| \nabla \theta(\tau)\|_{L^{\infty}} \|
v(\tau)\|_{{B}^{s-\alpha}_{p,1}}d\tau\\
&\leq \sup_{0\leq\tau\leq t}\|
v(\tau)\|_{{B}^{s-\alpha}_{p,1}}\int_{0}^{t}
 \|\theta(\tau)\|_{{B}^{1}_{\infty,1}}d\tau \\
&\leq Ce^{CV(t)}(1+t)\|\theta^0\|_{B^{s-\alpha}_{p,1}}
\Big( \|v^0\|_{{B}^{s}_{p,1}}+\int_0^t\|\theta(\tau)\|_{{B}^{s-\alpha}_{p,1}}
d\tau\Big).
\end{aligned}
\end{equation}
Applying  \eqref{2j7j8j} and  \eqref{27}, the estimate of $\theta$  reads 
as follows
\begin{equation} \label{092j7}
\begin{aligned}
& \|\theta(t)\|_{\widetilde{L}^{\infty}_{t}({B}^{s-\alpha}_{p,1})}+\int_0^t\|
\theta(\tau)\|_{{B}^{s}_{p,1}}d\tau\\
&\leq C\mathrm{e}^{C V(t)}\Big(\|\theta^0\|_{B^{s-\alpha}_{p,1}}(1+t)+
 \int_{0}^{t}
 \| \nabla \theta(\tau)\|_{L^{\infty}} \|
v(\tau)\|_{{B}^{s-\alpha}_{p,1}}d\tau\Big)\\
&\leq  C\mathrm{e}^{C V(t)}(1+t)\|\theta^0\|_{B^{s-\alpha}_{p,1}}
\Big(1+\| v^0\|_{{B}^{s}_{p,1}}
+ \int_{0}^{t}\| \theta(\tau)\|_{{B}^{s-\alpha}_{p,1}}d\tau\Big).
 \end{aligned}
\end{equation}
Combining  the estimates of $v$ (applying Lemma \ref{lem2.4} again)
and \eqref{092j7}, we have
\begin{equation}\label{9999}
\begin{aligned}
&\| v(t)\|_{{B}^{s}_{p,1}}+\|
\theta(t)\|_{{B}^{s-\alpha}_{p,1}}+\int_0^t\|
\theta(\tau)\|_{{B}^{s}_{p,1}}d\tau\\
&\leq  C\mathrm{e}^{C V(t)}(1+t)\|\theta^0\|_{B^{s-\alpha}_{p,1}}
\Big(1+\|v^0\|_{{B}^{s}_{p,1}}+  \int_{0}^{t}
 \| \theta(\tau)\|_{{B}^{s-\alpha}_{p,1}}d\tau\Big)\\
&\leq  \Phi_4(t)\Big(1+ \int_{0}^{t}\|
 \theta(\tau)\|_{{B}^{s-\alpha}_{p,1}}d\tau\Big).
\end{aligned}
\end{equation}
Therefore, by Gronwall's inequality, we finally obtain
\begin{equation}\label{9999b}
\| v(t)\|_{{B}^{s}_{p,1}}+\|
\theta(t)\|_{{B}^{s-\alpha}_{p,1}}+\int_0^t\|
\theta(\tau)\|_{{B}^{s}_{p,1}}d\tau\leq \Phi_5(t).
\end{equation}
This completes the first part (when assumption \eqref{4} holds) of this theorem.

 The proof of the second part (when assumption \eqref{5} holds) of this theorem
is quite similar to the one in the first part. The main
difference is the estimates of $\|\omega\|_{L^a}$ in  step 1.
 We begin with \eqref{4.7}, and by the
embedding $B^0_{a,2}\hookrightarrow L^{a}$ and Lemma \ref{lem4.1} (i), we have
\begin{equation}\label{4.7j}
\begin{aligned}
\|[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta(t)\|_{L^{a}}
&\leq C\|[\mathcal {R}_{\alpha}, v\cdot\nabla]\theta(t)\|_{B^0_{a,2}}\\
&\leq C \|\nabla v(t)\|_{L^a}\left(\|\theta(t)\|_{B^{1-\alpha}_{\infty,2}}
 +\|\theta(t)\|_{L^a}\right).
\end{aligned}
\end{equation}
It follows that
\begin{equation}\label{4.8j}
\begin{aligned}
\|\omega(t)\|_{L^a}
&\leq C\left(\|\omega^0\|_{L^{p}\cap
L^{\infty}}+\|\theta^0\|_{B^{s-\alpha}_{p,1}}\right)\\
&\quad +C\int_0^t\|\omega(\tau)\|_{L^a}
\left(\|\theta(\tau)\|_{B^{1-\alpha}_{\infty,2}}+\|\theta^0\|_{L^a}\right)d\tau.
\end{aligned}
\end{equation}
According to Gronwall's inequality we obtain
\begin{equation}\label{4.za9}
\|\omega\|_{L^a}
\leq Ce^{Ct}e^{C\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,2}}}.
\end{equation}
The estimate of $\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,2}}$ is
as follows, let $N\in\mathbb{N}$, by the Littlewood-Paley
decomposition, by condition \eqref{5} we see that
 %\label{4.10}
\begin{align*}
&\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,2}}\\
&\leq \|\mathcal {S}_N\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,2}}
 +\|(\text{Id}-\mathcal{S}_N)\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,1}}\\
&\leq \|\mathcal {S}_N\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,2}}
 +\sum_{q\geq N} 2^{q(1-\alpha)}\|\varDelta_q\theta\|_{L^{1}_{t}L^{\infty}}\\
&\leq Ct\|\theta^0\|_{L^{\infty}}
+C\int_0^t\Big(\sum_{0\leq q< N}
\|\varDelta_q\Lambda^{1-\alpha}\theta(\tau)\|^2_{L^{\infty}}\Big)^{1/2}d\tau
+\sum_{q\geq N} 2^{q(1-\alpha)}\|\varDelta_q\theta\|_{L^{1}_{t}L^{\infty}}
 \\
&\leq C\sqrt{N}\int_0^t\|\Lambda^{1-\alpha}\theta(\tau)
\|_{\dot{B}^0_{\infty,\infty}}d\tau
 +Ct\|\theta^0\|_{L^{\infty}}+C\sum_{q\geq N}
2^{q(1-\alpha+\frac{2}{a})}\|\varDelta_q\theta\|_{L^{1}_{t}L^{a}}\\
&\leq C\sqrt{N}  +Ct\|\theta^0\|_{L^{\infty}}+C\sum_{q\geq N}
2^{q(1-\alpha+\frac{2}{a})}\|\varDelta_q\theta\|_{L^{1}_{t}L^{a}}.
\end{align*}
Then, as for \eqref{4.11}, we can choose suitable $N$ such that
\begin{equation}\label{4.1d2j}
\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,2}}\leq
C(t+1)\|\theta^0\|_{L^{\infty}\cap L^{a}}
+C\Big[\log\Big(e+\int_0^t\|\omega(\tau)\|_{L^a}d\tau\Big)\Big]^{1/2}.
\end{equation}
Combining \eqref{4.za9}  with \eqref{4.1d2j}, we have
$$
\|\theta\|_{L^{1}_{t}B^{1-\alpha}_{\infty,2}}\leq \Phi_1(t).
$$
For the rest, we can follow the same process as above, and complete the proof.
\end{proof}

\subsection*{Acknowledgements}
I would like to  thank the referee for a careful reading of the work 
and many valuable comments. 
The research  is supported  by the Natural Science Foundation of 
Zhejiang Province (LQ16A010001) and NSFC 11501517.

\begin{thebibliography}{00}

\bibitem{[11]}  H. Abidi,  T. Hmidi;
\emph{On the global well-posedness for Boussinesq
system}. J. Differential Equations 233(2007)199-220.

\bibitem{[16]}  D. Adhikari, C. Cao,  J. Wu;
\emph{The 2D Boussinesq equations with
vertical viscosity and vertical diffusivity}. J. Differential
Equations. 249(2010)1078-1088.

\bibitem{[17]}  D. Adhikari, C. Cao, J. Wu;
\emph{Global regularity results for the 2D Boussinesq equations with
vertical dissipation}. J. Differential Equations
251(2011)1637-1655.

\bibitem{[38]} H. Bahouri,  R. Danchin,  J. Y. Chemin;
\emph{Fourier Analysis and Nonlinear Partial Differential Equations}, 
A Series of Comprehensive Studies in Mathematics. DOI 10.1007/978-3-642-16830-7,
Springer Heidelberg Dordrecht London New York. Springer-Verlag
Berlin Heidelberg (2011).

\bibitem{[6]} J. R. Cannon,  E. DiBenedetto;
\emph{The Initial Value Problem for the
Boussinesq Equations with Data in $L^p$}, in: Lecture Notes in
Math., vol. 771(1980), Springer, Berlin, pp. 129-144.

\bibitem{[81]}  C. Cao, J. Wu;
Global regularity for the 2D anisotropic
Boussinesq equations with vertical dissipation, (2011) arXiv:1108.2678v1
[math.AP].

\bibitem{[8]} D. Chae;
\emph{Global regularity for the 2D Boussinesq equations with
partial viscosity terms}, Adv. Math. 203 (2006), 497-513.

\bibitem{[24]} D. Chae;
\emph{Local existence and blow-up criterion for the Euler equations 
in the Besov spaces}.
Asymptotic Analysis. 39 (2004), 339-358.

\bibitem{[19]} D. Chae, J. Wu;
\emph{The 2D Boussinesq equations with logarithmically supercritical
velocities}. {Advances in Mathematics}. 230 (2012), 1618-1645.

\bibitem{[069]} A. C\'ordoba,  D. C\'ordoba;
\emph{A maximum principle applied
to quasi-geostrophic equations}. {Comm. Math. Phys.} 249 (2004), 511-528

\bibitem{[15]} R. Danchin,  M. Paicu;
\emph{Global existence results for the anisotropic Boussinesq system in 
dimension two}. {Math. Models Methods Appl. Sci.} 21 (2011), 421-457.

\bibitem{[82]} R. Danchin,  M.  Paicu, (  Existence and uniqueness results for
the Boussinesq system with data in Lorentz spaces,{ Phys. D}.
237(2008)1444-1460.

\bibitem{[83]} R. Danchin,  M. Paicu,    Global
well-posedness issues for the inviscid Boussinesq system with
Yudovich's type data, {Comm. Math. Phys.} 290(2009)1-14.

\bibitem{[21]}  T. Hmidi,  S. Keraani;
\emph{On the global well-posedness of the
two-dimensional Boussinesq system with a zero diffusivity},
{Adv. Differential Equations}. 12(2007)461-480.

\bibitem{[23]}  T. Hmidi,  S. Keraani;
\emph{ On the global well-posedness of the Boussinesq system with
zero viscosity}, {Indiana Univ. Math. J.} 58(2009)1591-1618.

\bibitem{[27]}  T. Hmidi,  S. Keraani, F. Rousset;
\emph{Global well-posedness for Euler-Boussinesq system with critical dissipation}, 
{Comm. Partial Differential Equations}. 36(2011)420-445.

\bibitem{[13]} T. Hmidi, M. Zerguine;
\emph{On the global well-posedness of the Euler-Boussinesq system with fractional
dissipation}. {Physica D}. 239(2010)1387-1401.

\bibitem{[14]}  T. Hmidi,  S. Keraani, F. Rousset;
\emph{Global well-posedness for a Boussinesq-Navier-Stokes system with 
critical dissipation}, {J. Differential Equations }. 249 (2010), 2147-2174.

\bibitem{[22]} T. Hou,  C. Li;
\emph{Global well-posedness of the viscous Boussinesq
equations}, {Discrete Contin. Dyn. Syst.}  12(2005)1-12.

\bibitem{[9]} A. Majda;
\emph{Introduction to PDEs and Waves for the Atmosphere and
Ocean}, Courant Lecture Notes in Mathematics, vol. 9, (2003) AMS/CIMS.

\bibitem{[18]}  C. Miao, L.  Xue;
\emph{On the global well-posedness of a class of
Boussinesq-Navier-Stokes systems}, NoDEA Nonlinear
Differential Equations Appl. 18(6) (2011), 707-735.

\end{thebibliography}

\end{document}
