The exact area calculation of irregularly distributed data is in the focus of all territorial geochemical balancing methods or definition of protection zones. Especially in the deep sea environment the interpolation of measurements into surfaces represents an important gain of information, because of cost- and time-intensive data acquisition.The geostatistical interpolation method indicator kriging therefore is applied for an accurate mapping of the spatial distribution of benthic communities following a categorical classification scheme at the deep-sea submarine Håkon Mosby Mud Volcano. Georeferenced video mosaics were obtained during several dives by the Remotely Operated Vehicle Victor6000 in a water depth of 1260 m. Mud volcanoes are considered as significant source locations for methane indicated by unique chemoautotrophic communities as Beggiatoa mats and pogonophoran tube worms. For the detection and quantification of their spatial distribution 2840 georeferencedvideo mosaics were analysed by visual inspection. Polygons, digitised on the georeferenced images within a GIS, build the data basis for geostatistically interpolated mono-parametric surface maps. Indicator kriging is applied to the centroids of the polygons calculating surfacemaps.The quality assessment of the surface maps is conducted by leave-one-out cross-validation evaluating the fit of the indicator kriging variograms by using statistical mean values. Furthermore, the estimate was evaluated by a validation dataset of the visual inspection of 530 video mosaics not included within the interpolation, thus, proving the interpolated surfaces independently. With regard to both validating mechanisms, we attained satisfying results and we provided each category applied for the identification of biogeochemical habitats with apercentage probability value of occurrence.