A recent perturbation model (Votyakov and Bonanos, 2014) is solved with an algebraic approximation. Results are in agreement with independent numerical results of Yang and Garimella (2010). The algebraic model reveals several scaling ratios as well as a non-monotonic behavior of the thermocline thickness as a function of fluid Reynolds number. (C) 2015 Elsevier Ltd. All rights reserved.