  In this article we prove that a classical $XY$ model subjected to weak i.i.d.
random field pointing in a fixed direction exhibits residual magnetic order in
$\mathbb{Z}^2$ and aligns perpendicular to the random field direction. The
paper is a sequel to \cite{NC} where the three-dimensional case was treated.
Our approach is based on a multi-scale Peierls contour argument developed in
\cite{NC}. On the microscopic scale we extract energetic costs from the
occurrence of contours, which themselves are defined on a macroscopic scale.
The technical challenges in $\mathbb{Z}^2$ stem from difficulties controlling
the size and roughness of the fluctuation fields which model the short
length-scale oscillations of near-optimizers of the random field Hamiltonian.
