Proceedings of the Industrial Simulation Conference 2007, 64-68 A new Cellular Automata model for the dynamics of sandpile surfaces is presented in this work. Following the BCRE model (Bouchaud et al., 1994) main ideas, we propose the use of two variables, instead of only the one that has been used in Cellular Automata models so far. The model predicts sandpile properties observed in other models, and in real systems, such as slope roughness decay or uphill propagation of avalanches. In addition, the introduction of the second variable allows the prediction of characteristics, such as the appearance of a second critical angle, that have been predicted by continuous models. The main advantage of the proposed model is that its formulation is simpler and easier to interpret than that of the continuous models. Also, the resulting numerical scheme is simpler and more efficient, allowing for a wider range of applications, including interactive simulation.