Inference for two-stage sampling in spatial surveys

This paper develops a design-based asymptotic theory for two-stage sampling over a continuous spatial domain. The target parameters are integral totals, or smooth functions of such totals, defined over a fixed bounded territory partitioned into an increasingly fine collection of primary sampling units. Within this fixed-area framework, we derive the order of the variance components of the Horvitz-Thompson estimator when secondary sampling units are points selected from continuous sub-regions. We establish design consistency of the estimator, as well as consistency of variance estimators under explicit regularity conditions on inclusion probabilities, sampling densities, and the study variable. The results are extended to plug-in estimators of smooth scaleinvariant functions of totals. For high-entropy first-stage designs, we further show that a H{á}jek-type variance estimator based only on first-order inclusion probabilities is consistent, and we prove asymptotic normality of both total and plug-in estimators. A simulation study inspired by forest inventory applications illustrates the finite-sample performance of the proposed estimators and variance estimators.