Yes, Birman-Series is true in variable negative curvature. Here is a sketch.
Suppose that $S$ is a closed, connected, oriented surface of genus $g > 1$. Fix a metric $g$ of variable negative curvature on $S$. Let $\Lambda$ be the Birman-Series set - the union of all simple geodesics.
The theory of train tracks is "soft", and so goes through in the variable negative curvature case. (For example, see Bonahon's 1996 paper Shearing hyperbolic surfaces, bending pleated surfaces and Thurston’s symplectic form.) So there are a finite number of almost geodesic train tracks in $S$ that, between them, carry all simple geodesics (with very small, and very well-spaced, "errors"). In any one track $\tau$, the set of carried curves of length at most $L$ has size at most $O(L^k)$ for some $k \leq 6 g - 6$. Let $\Lambda_\tau$ be the subset of $\Lambda$ obtained by taking the union of simple geodesics carried by $\tau$. Let $\Lambda_\tau(L)$ be the union of simple closed geodesics carried by $\tau$ with length at most $L$.
In variable negative curvature we still have "exponential convergence" of geodesics. So, let $b$ be a branch of $\tau$, one of the tracks given above. Let $I$ be a small arc dual to $b$. We examine the intersection of $I$ with $\Lambda_\tau$. These cluster about the intersections of $I$ with $\Lambda_\tau(L)$. We deduce that we can cover the intersection of $I$ with $\Lambda_\tau$ with polynomially many (in $L$) exponentially small (in $L$) intervals. Thus the cross-sectional Cantor set (to $\Lambda_\tau$) has Hausdorff dimension zero.
Here is a related example that confused me for a little while.
Suppose that $S$ is the surface of genus two. Consider a sequence $(S, g_t)_{t \to 0}$ of variable negative curvature metrics that Gromov-Hausdorff converge to the octagon surface $(S, g_0)$: the translation surface obtained by gluing opposite sides of a regular octagon. Let $\Lambda_t$ be the Berman-Series set for $(S, g_t)$ - the union of all simple geodesics (or equivalently, the closure of the union of all simple closed geodesics). Then the sets $\Lambda_t$ Hausdorff converge to $\Lambda_0$. However, by the above, if $t > 0$ then $\Lambda_t$ has Hausdorff dimension one while $\Lambda_0 = S$ has Hausdorff dimension two.
So we have a natural example of Hausdorff dimension jumping under taking Hausdorff limits. However, this sort of thing happens often enough that we should not be surprised. (Or at least, we should not be as surprised as I was when I realised the above was not an obstruction to the train track theory going through!)