Geostatistics meets a growing interest for the remediation forecast of potentially contaminated sites, by providing adapted methods to perform both chemical and radiological pollution mapping, to estimate contaminated volumes, potentially integrating auxiliary information, and to set up adaptive sampling strategies. As part of demonstration studies carried out for GeoSiPol (Geostatistics for Polluted Sites), geostatistics has been applied for the detailed diagnosis of a former oil depot in France. The ability within the geostatistical framework to generate pessimistic / probable / optimistic scenarios for the contaminated volumes allows a quantification of the risks associated to the remediation process: e.g. the financial risk to excavate clean soils, the sanitary risk to leave contaminated soils in place. After a first mapping, an iterative approach leads to collect additional samples in areas previously identified as highly uncertain. Estimated volumes are then updated and compared to the volumes actually excavated. This benchmarking therefore provides a practical feedback on the performance of the geostatistical methodology.