There is no perfect analogue because gravitational flux is not a (common) feature of GR.
Recall that Gauss' integral law $$\oint_{\partial V} \mathbf{g}\cdot\mathrm{d}\mathbf{A}=-4\pi G\underbrace{\int_V\rho\,\mathrm{d}V}_M$$
is just the differential law $$\operatorname{div}\mathbf{g}=-4\pi G\rho$$
integrated + the divergence theorem. So the two are equivalent.
Now, there is a pretty good analogue for the differential law in GR, namely the Einstein equations. It is well known that the time-time component of the Einstein equations reduces to the Poisson equation $\Delta\phi=4\pi G\rho$ for "nearly Newtonian" spacetimes. (The Poisson equation is equivalent to Gauss' law because $\mathbf{g}=-\operatorname{grad}\phi$.)
But it is not so simple to integrate the Einstein equations and get an integral law.
For asymptotically flat spacetimes, we have the ADM momentum equation
$$P^\mu=\int_\Sigma\sqrt{-g}\star [g^{\mu\sigma}(T_{\sigma\nu}+t_{\text{LL}\sigma\nu})\theta^\nu]$$
(Here $\Sigma$ is a spacelike hypersurface, $\star$ is the Hodge dual, $T$ is the energy-momentum tensor, $t_{\text{LL}}$ is the Landau-Lifsitz pseudotensor and $\theta^\mu$ is a dual tetrad.) The time component of this equation gives the "mass" of the spacetime, like the RHS of Gauss' law above. For for details, cf. e.g. N. Straumann, General Relativity (2013).