In 2021, we have developed a toolbox to accelerate the calculation of mass balances according to the Pareto-Soil method from geostatistical pollution models. The modelling results are processed in scripts in which the calculations are automated: they can thus be repeated easily for different variants, by modifying for example the thresholds of applied contents or the extension of the considered geographical domain. The cases treated during the last years have indeed shown us that the choice of the variants of depollution or remediation is iterative, with back and forth on the objectives and the treatment scenarios to be tested. This situation made it necessary to adapt our working methods to better support professionals and managers of polluted sites and soils in this reflection process. Our tools allow us to quickly produce excel files of the percentages of soil volumes and pollutant masses – by content intervals or by threshold exceedances – for all depths or by depth steps. The results are also provided in the form of figures. The scripts can be modified to be adapted to the specificities of the studied sites. Results are provided with estimation uncertainties rigorously quantified by the geostatistical approach. The uncertainties associated with the volume and mass percentages are represented as error bars on the figures.