I am working with the Laplacian on a Riemannian manifold $(M,g)$ (compact, without boundary). In spherical geodesic coordinates $(r, \sigma)$ around some arbitrary $x \in M$ (where $\sigma$ denotes the angular coorinates taken together), the Laplacian looks like $\Delta = \frac {\partial ^2} {\partial r^2} + H(x,r) \frac \partial {\partial r} + \Delta_{S(x,r)}$ where $S(x,r)$ is the geodesic sphere of radius $r$ around $x$ and $H(x,r)$ is the total mean curvature of $S(x,r)$.
I need to do some calculations at the end of which I must make $r \to 0$ (i.e. evaluate things in $x$). To this end, it would be very helpful to know the first few terms of the development of $H(x,r)$ with respect to $r$. So far, I only know that $H(x,r) = \frac {n-1} r + O(r)$, which is good but not enough. Does anybody know the coefficient of $r$ (and possibly of $r^2$) in the above expansion?