This process was repeated with the benthic habitat data. Lastly, the Map Algebra function was used to add the two reclassified rasters into a single layer in which larger values reflected both a habitat type and depth that correlated with greater occurrence of debris removed. This ‘‘hot spots’’ map layer was created twice: first, using all benthic habitat types including seagrasses for the full study area (332 km2), and second, restricting the analysis to non-seagrass habitats (i.e. hard
bottom) to target coral reef areas targeted by the managementplan of BNP (80 km2).