Design-based consistency of the Horvitz–Thompson estimator under spatial sampling with applications to environmental surveys