I am using turning band simulations to produce estimates of a quantity distributed across a 2-dimensional field with uncertainty, as described in ICES' "Handbook of Geostatistics in R for Fisheries and Marine Ecology" (Petitgas et al. 2017). I am unclear as to how variability arises in the simulated realisations (i.e. how do different realisations arise from a single dataset and a deterministic variogram). Is there any literature you could point to which would give a thorough explanation of this?
Many thanks.