Abstract. We construct a convolution-thresholding approximation scheme for the geometric surface evolution in the case when the velocity of the surface at each point is a given function of the mean curvature. Conditions for the monotonicity of the scheme are found and the convergence of the approximations to the corresponding viscosity solution is proved. We also discuss some aspects of the numerical implementation of such schemes and present several numerical results. 1. Introduction. The topic of curvature flows of different types was popular during the last 20 years and is still popular in both pure and applied mathematics. By curvature flow we mean a family {Γ t } t≥0 of hypersurfaces in R n depending on time t with local normal velocity equal to the mean curvature or a function of it for generalized curvature flows. The mean curvature in turn denotes here the sum of principal curvatures.In the three-dimensional case a smooth initial surface can develop singularities after some finite time. There have been several successful attempts to deal with singularities and topological complications: the varifold approach [7], [2], the phase field method [14], [8], and the level-set method. This approach was suggested in the physical literature [26] and was extensively developed for numerical purposes by Osher and Sethian [27]. The main idea of this method is to evolve some continuous function u : [0, ∞) × R n → R in such a way that Γ t ⊂ R n would always be a level-set of u (x, t), i.e., Γ t = {x ∈ R n : u (x, t) = 0} for all t ≥ 0. In the case of the mean curvature flow, the evolution equation for u turns out to be