We calculate Sorkin's manifestly covariant entanglement entropy S for a massive and massless minimally coupled free Gaussian scalar field for the de Sitter horizon and Schwarzschild de Sitter horizons respectively in d > 2. In de Sitter spacetime we restrict the Bunch-Davies vacuum in the conformal patch to the static patch to obtain a mixed state. The finiteness of the spatial L 2 norm in the static patch implies that S is well defined for each mode. We find that S for this mixed state is independent of the effective mass of the scalar field, and matches that of [1], where, a spatial density matrix was used to calculate the horizon entanglement entropy. Using a cut-off in the angular modes we show that S ∝ Ac, where Ac is the area of the de Sitter cosmological horizon. Our analysis can be carried over to the black hole and cosmological horizon in Schwarzschild de Sitter spacetime, which also has finite spatial L 2 norm in the static regions. Although the explicit form of the modes is not known in this case, we use the boundary conditions of [2] for a massless minimally coupled scalar field, to find the mode-wise S b,c , where b, c denote the black hole and de Sitter cosmological horizons, respectively. As in the de Sitter calculation we see that S b,c ∝ A b,c after taking a cut-off in the angular modes.