Does anyone know where I could find a proof of a variant of a version of the mean-value property for harmonic functions in Riemannian manifolds?  I'm actually more interested in using an elliptic operator applied to a function, u, integrated over geodesic balls, to control the size of u.  Locally is enough; I don't care about issues involving caustics, et cetera.