Does anyone know where I could find a proof of a variant of a version of the meanvalue 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.

Schoen & Yau, Lectures on Differential Geometry (International Press, 1994), Chapter II, Section 6: Mean value inequality for subharmonic functions. This is what you want. 


You may try: (1) Grigorʹyan, A. A. Stochastically complete manifolds and summable harmonic functions. (Russian) Izv. Akad. Nauk SSSR Ser. Mat. 52 (1988), no. 5, 11021108, 1120; translation in Math. USSRIzv. 33 (1989), no. 2, 425–432 (2) Leon Karp, Subharmonic functions, harmonic mappings and isometric immersions (pp. 133142); in Seminar on Differential Geometry. Papers presented at seminars held during the academic year 1979–1980. Edited by Shing Tung Yau. Annals of Mathematics Studies, 102. Princeton University Press, Princeton, N.J.; University of Tokyo Press, Tokyo, 1982. 

