Let $(M,g)$ be a negatively curved, closed Riemannian manifold. I'll ask the question first, then explain the involved players. This data defines the Patterson-Sullivan length function, \begin{align} L_{g}:\mathcal{M}_{-}\rightarrow \mathbb{R}, \end{align} where $\mathcal{M}_{-}$ is the space of Riemannian metrics on $M$ which have negative sectional curvature. The question is the following, suppose $g_{t}$ is a $C^k$-one parameter family of negatively curved Riemannian metrics defined on some open interval $(a,b).$ What is the regularity of $L_g:(a,b)\rightarrow \mathbb{R}$ defined by $t\mapsto L_{g}(g_t)?$ In an ideal world, if the variation in metrics is real analytic, this function would also vary real analytically. There are some known instances of this phenomenon. For instance, the paper of Bridgeman-Taylor shows the veracity of this result for the slightly different case of convex-cocompact hyperbolic metrics (which are infinite volume, but it's a very similar setup). https://www2.bc.edu/~bridgem/papers/wpextension.pdf
Unfortunately, their proof relies crucially on some properties of the deformation theory of convex-cocompact hyperbolic metrics, and so cannot be carried over to this context.
Some explanation now. For a pair of negatively curved metrics $g,g',$ there is a conjugacy (for the $\pi_1(M)$-action) of the space of geodesics in the universal cover $\widetilde{M},$ \begin{align} f:\mathcal{G}(g)\rightarrow \mathcal{G}(g'). \end{align} The Patterson-Sullivan geodesic current $\sigma_g$ (for g) is a measure on the space of geodesics $\mathcal{G}(g)$ which is invariant under the action of $\pi_1(M).$ Using $f,$ one obtains the push-forward measure $f_{*}(\sigma_g).$ Combining this measure with arc-length measure for $g'$ along geodesics gives a measure on the total space of an $\mathbb{R}$-bundle over the space of geodesics, \begin{align} T^{1}(\widetilde{M})\rightarrow \mathcal{G}(g'). \end{align} The map here sends a unit vector to the unique geodesic to which it is tangent. This picture is invariant under the $\pi_1(M)$-action, and taking the quotient by the action yields a measure on the unit tangent bundle $T^{1}(M).$ The function $L_g(g')$ records the total mass of this measure.
I didn't define the Patterson-Sullivan geodesic current here, but I'm banking on the fact that anyone who answers will already know the construction of this measure.
If anyone would like more detail, I'll be happy to write it up. Let me know.