We study the dynamics of the geodesic flow for a class of non-complete Riemannian metrics on a negatively curved surface M . These finite area surfaces are composed of finitely many “singular” surfaces of revolutions of the form y = x, r > 1 for 0 ≤ x ≤ 1 (thin pieces), together with connecting surfaces of bounded negative curvature (thick pieces). The curvature at every point is negative and bounded away from zero, but is unbounded in the “cusps.” The geodesic flow is non-complete because geodesics corresponding to singular unit tangent vectors pointing into the origin hit the cusps in finite time and then cease to be defined. Such singular unit tangent vectors are dense in the unit tangent bundle. Our main result is ergodicity of the geodesic flow (Theorem 5.1). It immediately follows that these metrics have dense geodesics in the unit tangent bundle SM . We also prove that the closed geodesics are dense in the unit tangent bundle (Theorem 6.1) The motivation for considering these metrics arises from studying the geodesic flow for the Weil-Petersson (WP) metric on two dimensional moduli spaces of Riemann surfaces., e.g., the moduli space for the once punctured torus. In several fundamental ways the geodesic flow we consider provides a good model for the WP geodesic flow. For example, Wolpert showed that the WP metric is also non-complete, has finite area, has negative curvature bounded away from zero, and is unbounded in the ”cusps.” Also, various authors have shown that the cusps can be approximated by singular surface of revolutions. At present, there is only a single technical obstruction (additional estimates on the derivatives of the GF in the cusp) that prevent us from applying this general method to prove ergodicity of the WP geodesic flow. See Section 7. The geometric properties of both types of flows imply that the geodesic flow is a uniformly hyperbolic dynamical system with singularities. The non-completeness causes pathologies in the stable and unstable manifolds, as for many billiard flows. For example, stable and unstable manifolds at a point, if they exist, may intersect a cusp point, and thus have only finite length. One needs extensions of “Pesin theory” to systems with singularities to study these geodesic flows. Another challenging aspect is that these surfaces may be simply connected, making the use of any of the traditional arguments involving boundaries of the covering spaces inapplicable. Studying the dynamics requires the development of new approches. Our strategy to prove ergodicity of the geodesic flow is to first establish non-uniform hyperbolicity, i.e., the Lyapunov exponents are non-zero at almost every point. We then prove that there