A method includes using a computing system in generating a land surface or sea floor digital image from post-stack seismic data and determining a clathrate stability zone from the land surface or sea floor digital image. The method also includes generating a laterally extensive subsurface digital image from the seismic data that is within or below the clathrate stability zone and filtering the laterally extensive subsurface digital image to generate a pull-up digital image. The method also includes generating a coarse-grained sediment digital image from the seismic data within the clathrate stability zone. The method also includes spatially superimposing the coarse-grained sediment digital image and the pull-up digital image to generate a combined digital image. The method also includes determining an overlap region from the combined digital image. The overlap region is indicative of a clathrate deposit.