This paper presents an integrative hierarchical stepwise sampling (IHS) method and two case studies to compare it with stratified random sampling (SRS) and conditioned Latin hypercube sampling (cLHS). The first comparison between IHS and SRS was conducted for mapping sand content of two soil layers in a study area in Anhui Province, China. Two sample sets of the same sample size were collected in the field based on IHS and SRS. The second case study is a simulation study, where we compared IHS and cLHS for mapping soil series in the Raffelson watershed in Wisconsin (USA). The study used an accurate and detailed soil series map produced previously as a proxy of the real soil distribution. Virtual samples with nine sample sizes designed by IHS and cLHS were collected on the soil map. For both case studies, an individual predictive soil mapping method was employed and independent validation samples were used to evaluate the mapping accuracies. Results indicate that IHS generally performs better than SRS for capturing distributions of the environmental variables. It obtained higher mapping accuracies than SRS at different sample sizes. On the other hand, cLHS appears to provide a better representation for distributions of the environmental variables than IHS, but the mapping accuracies with IHS are higher than those with cLHS at nearly all sample sizes. Finally, both case studies showed that IHS provides valuable information on representativeness of the samples.Abbreviations: FCM, fuzzy c-means; IHS, integrative hierarchical stepwise sampling; cLHS, conditioned Latin hypercube sampling; SRS, stratified random sampling. S oil sampling plays an essential role in soil survey. A sampling design specifies where and how many to sample in the field. Different sampling design strategies will produce different soil distribution maps and greatly affect the mapping accuracies (Brus and de Gruijter, 1997;van Groenigen et al., 2000;de Gruijter et al., 2006;Brus et al., 2006;Gregoire and Valentine, 2007;Heim et al., 2009). Because the collection of soil samples in the field is often time-consuming and resource-intensive, it is important to apply efficient sampling design strategies.When prior knowledge on an area's local soil variation is limited, simple random sampling has been a commonly used sampling strategy (Oliver and Webster, 1986;Grinand et al., 2008;Walvoort et al., 2010 • Evaluation of an integrative hierarchical stepwise (IHS) sampling method by comparing it with stratified random sampling (SRS) and conditioned Latin hypercube sampling (cLHS) through two case studies.• IHS obtained higher mapping accuracies than SRS and cLHS at nearly all sample sizes.• IHS provides valuable information on the representativeness of samples.• SRS and cLHS were found to generate unstable results on sample sets and soil maps.