The two-dimensional (2D) incompressible Boussinesq system is not only an important model in geophysics, but also retains some key features of the 3D Euler and Navier-Stokes equations such as the vortex stretching mechanism. ...