A result by Elton(6) states that an iterated function system $$M_n = F_n (M_{n - 1} ),{\text{ }}n \geqslant 1,$$ of i.i.d. random Lipschitz maps F1,F2,... on a locally compact, complete separable metric space \((\mathbb{X},d)\) converges weakly to its unique stationary distribution π if the pertinent Liapunov exponent is a.s. negative and \(\mathbb{E}\log ^ + d(F_1 (x_0 ),x_0 ) 0, where L1 denotes the Lipschitz constant of F1. The same and also polynomial rates have been recently obtained in Alsmeyer and Fuh(1) by different methods. In this article, necessary and sufficient conditions are given for the positive Harris recurrence of (Mn)n≥0 on some absorbing subset \(\mathbb{H}{\text{ of }}\mathbb{X}\). If \(\mathbb{H} = \mathbb{X}\) and the support of π has nonempty interior, we further show that the same respective moment conditions ensuring the weak convergence rate results mentioned above now lead to polynomial, respectively geometric rate results for the convergence to π in total variation ∥⋅∥ or f-norm ∥⋅∥f, f(x)=1+d(x,x0)η for some η∈(0,p]. The results are applied to various examples that have been discussed in the literature, including the Beta walk, multivariate ARMA models and matrix recursions.