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.