Northern Puerto Rico hosts karst within carbonate units of Oligocene and Miocene age. Lidar-derived elevation data, acquired in 2018, were used to create a mosaic of 3-meter resolution working digital elevation models (DEMs), from which surface depressions were identified using a semi-automated workflow in ArcGIS. Depressions in the automated inventory were systematically checked using aerial imagery, lidar-derived imagery, and 3D viewing of the lidar imagery. Center points were created for each surface depression and were used to create a point density raster. The density raster displays the number of closed depression points per square kilometer.