The great majority of the world population resides nowadays in urban areas. Understanding their physical and social structure, and especially their urban land use pattern dynamics throughout time, becomes crucial for successful, effective management of such areas. This study is committed to simulate and predict urban land use change in a pilot city belonging to the São Paulo Metropolitan Region, southeast of Brazil, by means of a cellular automata model associated with the Markov chain. This model is driven by data derived from orbital and airborne remotely sensed images and is parameterized by the Bayesian weights of evidence method. Several layers related to infrastructure and biophysical aspects of the pilot city, São Caetano do Sul, were used as evidence in the simulation process. Alternative non-stationary scenarios were generated for the short-run, and the results obtained from past simulations were statistically validated using a multiresolution “goodness-of-fit” metric relying on fuzzy logic. The best simulations reached fuzzy similarity indices around 0.25–0.58 for small neighborhood windows when an exponential decay approach was employed for the analysis, and approximately 0.65–0.95 when a constant decay and larger windows were considered. The adopted Bayesian inference method proved to be a good parameterization approach for simulating processes of urban land use change involving multiple classes and transitions.