“…The choice of the plurigaussian approach is motivated by the fact that it enables the incorporation of geological knowledge, through the definition of the contact relationships between the alteration domains, the modeling of spatially-varying domain proportions, and the reproduction of conditioning data from exploration drill holes. In the last decade, this approach has proven to be versatile and has been applied to the modeling of ore deposits with varied geological settings, such as granitehosted deposits [19], roll-front uranium deposits [8,16], kimberlite pipes [6], copper porphyries [4,17,21], nickel laterites [14], or carbonate-hosted lead-zinc deposits [22].…”