“…Moreover, it will be very hard to get the high‐order estimates because we need a higher order estimate when we have obtained a high‐order estimate (e.g., when we estimate ∥∇ u ∥, we need to estimate ∥Δ u ∥); it is an endless iteration. However, in the case of the full dissipation ( ν = ν ( θ ) > 0, κ = κ ( θ ) > 0), this problem will not happen (e.g., ), because the higher order terms on the right can be controlled by the left terms. In order to overcome them, we are going to utilize the vanishing diffusivity method to prove the existence of global quasi‐strong solutions and get some higher order uniform estimates, and furthermore prove the global well‐posedness of the Boussinesq system () and ().…”