In this work, we used molecular dynamics (MD) to perform trajectory simulations of ice-like argon and amorphous silica aggregates on atomically smooth highly ordered pyrolytic graphite (HOPG) and a comparatively rougher quartz surface. It was found that at all incidence velocities, the quartz surface was stickier than the HOPG surface. The sticking probabilities and elastic moduli obtained from MD were then used to model surface evolution at a micron length scale using kinetic Monte Carlo (kMC) simulations. Rules were derived to control the number of sites available for the process execution in kMC to accurately model erosion of HOPG by atomic oxygen (AO) attack and ice-nucleation on surfaces. It was observed that the effect of defects was to increase the material erosion rate, while that of aggregate nucleation was to lower it. Similarly, simulations were performed to study the effects of AO attack and N2 adsorption-desorption on surface evolution and it was found that N2 adsorption-desorption limits the surface available for erosion by AO attack.