An increase in the coverage dependent adsorption probability with increasing coverage, i.e., an autocatalytic adsorption phenomenon, has frequently been observed. We present a time dependent (kinetic) Monte Carlo simulation scheme which can also account for this phenomenon. The scheme is tested for different adsorption scenarios and, is compared with experimental results of CO adsorption on polar ZnO surfaces [T. Becker et al., J. Chem. Phy. 113, 6334 (2000)]. (C) 2003 American Vacuum Society.