An advantage of geostatistics for mapping soil properties is that, given a statistical model of the variable of interest, we can make a rational decision about how densely to sample so that the map is sufficiently precise. However, uncertainty about the statistical model affects this process. In this paper we show how Bayesian methods can be used to support decision making on sampling with an uncertain model, ensuring that the probability of meeting certain levels of precision is high enough.