 Solution

The action contains neither $\dot{N}$ nor ${\dot{N}}^i$, so their conjugate momentum vanishes. They impose constraints rather than independent evolution equations: after the Legendre transform in mechanics, the lapse function and shift vector are Lagrange multipliers in the Hamiltonian formulation of general relativity.

Use precisely the paper's sign convention for extrinsic curvature, $K_{ij}=({\dot{h}}_{ij}-D_i{N}_j-D_j{N}_i)/(2N)$. At fixed $h_{ij},N,N^i$,Consequently the canonical momentum of the spatial metric isThis is a weight-one tensor density, with indices raised and lowered using $h_{ij}$. Its trace and inverse relation areThe sign of $K_{ij}$ used here is opposite to the convention in extrinsic curvature of a spatial hypersurface; mixing the two conventions would reverse the momentum formula.