  The Thorndike et al., (\emph{J. Geophys. Res.} {\bf 80} 4501, 1975) theory of
the ice thickness distribution, $g(h)$, treats the dynamic and thermodynamic
aggregate properties of the ice pack in a novel and physically self-consistent
manner. Therefore, it has provided the conceptual basis of the treatment of
sea-ice thickness categories in climate models. The approach, however, is not
mathematically closed due to the treatment of mechanical deformation using the
redistribution function $\psi$, the authors noting ``The present theory suffers
from a burdensome and arbitrary redistribution function $\psi .$'' Toppaladoddi
and Wettlaufer (\emph{Phys. Rev. Lett.} {\bf 115} 148501, 2015) showed how
$\psi$ can be written in terms of $g(h)$, thereby solving the mathematical
closure problem and writing the theory in terms of a Fokker-Planck equation,
which they solved analytically to quantitatively reproduce the observed winter
$g(h)$. Here, we extend this approach to include open water by formulating a
new boundary condition for their Fokker-Planck equation, which is then coupled
to the observationally consistent sea-ice growth model of Semtner (\emph{J.
Phys. Oceanogr.} {\bf 6}(3), 379, 1976) to study the seasonal evolution of
$g(h)$. We find that as the ice thins, $g(h)$ transitions from a single- to a
double-peaked distribution, which is in agreement with observations. To
understand the cause of this transition, we construct a simpler description of
the system using the equivalent Langevin equation formulation and solve the
resulting stochastic ordinary differential equation numerically. Finally, we
solve the Fokker-Planck equation for $g(h)$ under different climatological
conditions to study the evolution of the open-water fraction.
