We analyze the collective surface diffusion coefficient, D c , near a first-order phase transition at which two phases coexist and the surface coverage, θ, drops from one single-phase value, θ + , to the other one, θ − . Contrary to other studies, we consider the temperatures that are sufficiently subcritical. Using the local equilibrium approximation, we obtain, both numerically and analytically, the dependence of D c on the coverage and system size, N , near such a transition. In the twophase regime, when θ ranges between θ − and θ + , the diffusion coefficient behaves as a sum of two hyperbolas, D c ≈ A/N |θ − θ − | + B/N |θ − θ + |. The steep hyperbolic increase in D c near θ ± rapidly slows down when the system gets from the two-phase regime to either of the single-phase regimes (when θ gets below θ − or above θ + ), where it approaches a finite value. The crossover behavior of D c between the two-phase and single-phase regimes is described by a rather complex formula involving the Lambert function. We consider a lattice-gas model on a triangular lattice to illustrate these general results, applying them to four specific examples of transitions exhibited by the model.