In the geothermal Euganean area (Veneto region, NE Italy) water temperatures range from 60 to 86 degreesC. The aquifer considered is rocky and the production wells in this study have a depth ranging from 300 to 500 m. For exploitation purposes, it is important to identify zones with a high probability that the temperature is more than 80 degreesC and zones with a high probability that the temperature is less than 70 degreesC First, variographic analysis was conducted from 186 temperature data of thermal ground waters. This analysis gave results that are consistent with the main regional tectonic structure, the NW-SE trending "Schio-Vicenza" fault system. Then indicator variograms of the second, fifth, and eighth decile were compared to identify the spatial continuity at different thresholds. The Unacceptability of a multigaussian hypothesis of the random function and the necessity to know the cumulative distribution function in any location, suggested the use of a nonparametric geostatistical procedure such as indicator kriging. Thus, indicator variograms at the cutoffs of 65, 70, 73, 75, 78, 80, 82, and 84 degreesC were analyzed, fitted, and used during the indicator kriging procedure. Finally. probability maps were derived from postprocessing indicator kriging results. These maps identified scarcely exploited areas with a high probability of the temperature being higher than 80 degreesC, between 70 and 80'C and areas with high probability of the temperature being below 70 degreesC.