Starting with a many-atom master equation of a kinetic, restricted solid-on-solid (KRSOS) model with external material deposition, we investigate nonlinear aspects of the passage to a mesoscale description for a crystal surface in 1+1 dimensions. This latter description focuses on the motion of an atomic line defect (i.e. a step), which is defined via appropriate statistical average over KRSOS microstates. Near thermodynamic equilibrium and for low enough supersaturation, we show that this mesoscale picture is reasonably faithful to the Burton-Cabrera-Frank (BCF) stepflow model. More specifically, we invoke a maximum principle in conjunction with asymptotic error estimates to derive the elements of the BCF model: (i) a diffusion equation for the density of adsorbed adatoms; (ii) a step velocity law; and (iii) a linear relation for the mass flux of adatoms at the step. In this vein, we also provide a criterion by which the adatom flux remains linear in supersaturation, suggesting a range of non-equilibrium conditions under which the BCF model remains valid. Lastly, we make use of kinetic Monte Carlo simulations to numerically study effects that drive the system to higher supersaturations -e.g. deposition of material onto the surface from above. Using these results, we describe empirical corrections to the BCF model that amount to a nonlinear relation for the adatom flux at the step edge.