In this paper we generalize the algebraic density property to not necessarily smooth affine varieties relative to some closed subvariety containing the singular locus. This property implies the remarkable approximation results for holomorphic automorphisms of the Andersen Lempert theory. We show that an affine toric variety X satisfies this algebraic density property relative to a closed T-invariant subvariety Y if and only if X \ Y not equal T. For toric surfaces we are able to classify those which possess a strong version of the algebraic density property (relative to the singular locus). The main ingredient in this classification is our proof of an equivariant version of Brunella's famous classification of complete algebraic vector fields in the affine plane. (C) 2014 Elsevier B.V. All rights reserved.