I'm a bit confused about the tensor density weight of Poisson brackets in general relativity and their covariance. It's perhaps related to being unclear as to what happens when I integrate a scalar density of some weight other than 1. Let's say I have the Poisson bracket of General Relativity in the 3+1 ADM formalism acting on some local scalar $f(x)$ on a space-time slice and some scalar quantity $G$. ($G$ could be the Hamiltonian, and $f(x)$ could be a scalar, but it could also be a scalar density, e.g .$\sqrt{g}$ which changes things but not the essence of what I'm asking). The Poisson bracket is given by \begin{align} \{f(x),G\}&=\int d^3y \Big[ \frac{\delta f(x)} {\delta g_{ab}(y)}\frac{\delta G}{\delta \pi^{ab}(y)} - \frac{\delta f(x)} {\delta \pi^{ab}(y)}\frac{\delta G}{\delta g_{ab}(y)}\Big] \\ &\stackrel{?}{=} \frac{\delta f(x)} {\delta g_{ab}(x)}\frac{\delta G}{\delta \pi^{ab}(x)} - \frac{\delta f(x)} {\delta \pi^{ab}(x)}\frac{\delta G}{\delta g_{ab}(x)} \end{align} with $g_{ab}$ the 3-metric and using the convention of taking its conjugate momenta $\pi^{ab}$ a tensor density of weight one (since we derive it from the Lagrangian density). 2 questions: The first is that the tensor weight of the first expression seems to be -2 (plus whatever comes with $f$, since I have the $d^3y$ on the top and the $\delta \pi^{ab}$ on the bottom. Since the left hand side is usually something like $\partial_tf(x)$, I would have expected it to have a tensor weight of 1. And this expression doesn't look like it will give diffeomorphism invariance, although I accept that it must (I guess one needs to consider how the 3-manifold sits in the 4-manifold for this).
There is some discussion of the invariance properties of the Poisson bracket here: Poisson brackets in curved spacetime, but I don't find it particularly enlightening. Anyone have a simple explanation?