In the present study, we co-simulate hydrofacies and piezometric data in order to construct geostatistical realizations of underground geology in an area of the West Thessaly basin. This basin is of great importance in terms of sustainable water management and environmental perspective in Greece. Through Plurigaussian modeling, the hydrofacies are first transformed into Gaussian Random Fields. …