We derive new integral estimates on substatic manifolds with boundary of horizon type, naturally arising in General Relativity. In particular, we generalize to this setting an identity due to Magnanini-Poggesi [24] leading to the Alexandrov Theorem in R-n and improve on a Heintze-Karcher type inequality due to Li-Xia [22]. Our method relies on the introduction of a new vector field with nonnegative divergence, generalizing to this setting the P-function technique of Weinberger [36]. (C) 2022 Elsevier Masson SAS. All rights reserved.