“…By a 'temperature patch solution', we mean a solution of (1.1) where θ is an indicator function for each time t ≥ 0. These are similar to the sharp front solutions of the SQG equation and their variants which have been and are extensively studied [47], [46], [20], [19], [18], [21], [13], [5], [7], [38], [39], [34], [35], [32], [33], and the more classical theory of vortex patches for the 2D Euler equation as laid out in [41], and also [10], [11], and [37]. The Boussinesq system also supports initial data of 'Yudovich' type.…”