In geostatistics, Gaussian random fields are often used to model heterogeneities of soil or subsurface parameters. To give spatial approximations these fields, they discretized. Then, different techniques geostatistical inversion condition them on measurement data. Among techniques, Markov chain Monte Carlo (MCMC) stand out, because yield asymptotically unbiased conditional realizations. Howeve...