Abstract. We study a finite volume element discretization of a nonlinear parabolic equation in a convex polygonal domain. We show existence of the discrete solution and derive error estimates in $L_2$- and $H^1$-norms. We also consider a linearized method and provide numerical results to illustrate our theoretical findings.