Aubin's book - construction of Green's function on compact manifold

Your have two nice texts of Frederic Robert (a descendant of Aubin) about the construction and estimate of Green function

one in french: http://www.iecl.univ-lorraine.fr/~Frederic.Robert/ConstrucGreen.pdf

an other in english where you will find an answer to (b) when assuming Neumann boundary condition: http://www.iecl.univ-lorraine.fr/~Frederic.Robert/NotesGreenNeumannRobert.pdf

Else you can also have a look to appendix A of the book: Blow-up theory for elliptic PDEs in Riemannian geometry, Olivier Druet, Emmanuel Hebey and Frederic Robert, Mathematical Notes, Princeton University Press, Volume 45.

The answer to (a) is easier than you think. As you remark, Aubin shows the estimate holds near the diagonal, so for some $\epsilon$, for all $P$, $Q$ with $\rho(P,Q)<\epsilon$,

$G(P,Q) \leq k \rho(P,Q)^{2-n}.$

However, since the manifold is compact, the continuous function $G$ is bounded (say, by the real number $A$) on the complement of this epsilon-neighbourhood of the diagonal. Thus for all $P$, $Q$

$G(P,Q) \leq max(k, A \rho(P,Q)^{n-2})\rho(P,Q)^{2-n}\leq \max(k, A \ diam(M)^{n-2})\rho(P,Q)^{2-n}$.