Trajectories of fast (0.55 MeV) protons impinging under planar surface channeling conditions with an angle of incidence of several mrad are calculated on KBr(001) surfaces decorated by electron stimulated desorption (ESD). Two-dimensional angular distributions of protons scattered on the surface are obtained. In order to calculate trajectories, a computational technique is needed to form a set of morphologies of the surfaces according to the electron ‰uence. In our method, an assumption is made that the topmost atoms on the low coordinate sites are easily emitted by the ESD processes 5) . By choosing parameters introduced in the method, the oscillatory behavior of the desorption yields demonstrated by other groups is reproduced 2) . The resulting peak intensity of the re‰ected protons from the trajectory on the ensemble of surfaces with varying electron ‰uence oscillates as the electron ‰uence increases. Corresponding angular distributions of re‰ected protons have been measured on the electron irradiated KBr(001) surfaces, and a similar oscillation can be seen on the scattered yields with weak amplitudes.