  In this paper we establish the asymptotic stability of steady solutions for
the Boussinesq systems in the framework of Cartesian product of critical
weak-Morrey spaces on $\mathbb{R}^n$, where $n \geqslant 3$. In our strategy,
we first establish the continuity for the long time of the bilinear terms
associated with the mild solutions of the Boussinesq systems, i.e., the
bilinear estimates by using only the norm of the present spaces. As a direct
consequence, we obtain the existence of global small mild solutions and
asymptotic stability of steady solutions of the Boussinesq systems in the class
of continuous functions from $[0, \infty)$ to the Cartesian product of critical
weak-Morrey spaces. Our techniques consist interpolation of operators, duality,
heat semigroup estimates , Holder and Young inequalities in block spaces (based
on Lorentz spaces) that are preduals of Morrey-Lorentz spaces. Our results are
generalized the previous ones of the steady Boussinesq systems in weak-$L^p$
spaces obtained by Hishida [T. Hishida, {\it On a class of Stable Steady flow
to the Exterior Convection Problem}, Journal of Differential Equations, Vol.
141, Iss. 1 (1997), pages 54-85] and Ferreira et al. [L.C.F. Ferreira and E.J.
Villamizar-Roa, {\it On the stability problem for the Boussinesq equations in
weak-$L^p$ spaces}, Commun. Pure Appl. Anal. (2010), Vol. 9, No. 3, pages
667-684] and of the Navier-Stokes equations in Morrey spaces obtained by Kozono
et al. [H. Kozono and M. Yamazaki, {\it The stability of small stationary
solutions in Morrey Spaces of the Navier-Stokes equation}, Indiana University
Mathematics Journal, Vol. 44, No. 4 (1995), pages 1307-1336].
