We introduce variational approximations for curve evolutions in two-dimensional Riemannian manifolds that are conformally flat, i.e. conformally equivalent to the Euclidean space. Examples include the hyperbolic plane, the hyperbolic disk, the elliptic plane as well as any conformal parameterization of a two-dimensional surface in R d , d ≥ 3. In these spaces we introduce stable numerical schemes for curvature flow and curve diffusion, and we also formulate a scheme for elastic flow. Variants of the schemes can also be applied to geometric evolution equations for axisymmetric hypersurfaces in R d . Some of the schemes have very good properties with respect to the distribution of mesh points, which is demonstrated with the help of several numerical computations. 1 2 .With the help of the splitting, we propose in Section 3 a semi-implicit scheme for which stability can be shown.The outline of this paper is as follows. In Section 2 we derive the governing equations for curvature flow, curve diffusion and elastic flow, provide weak formulations and relate the introduced flows to geometric evolution equations for axisymmetric hypersurfaces. In Section 3 we introduce finite element approximations and show existence and uniqueness as well as stability results. Section 4 is devoted to several numerical results, which