Seafloor pockmarks occur worldwide and may represent millions of m 3 of continental shelf erosion, but few numerical analyses of their morphology and spatial distribution of pockmarks exist. We introduce a quantitative definition of pockmark morphology and, based on this definition, propose a three-step geomorphometric method to identify and extract pockmarks from high-resolution swath bathymetry. We apply this GIS-implemented approach to 25 km 2 of bathymetry collected in the Belfast Bay, Maine USA pockmark field. Our model extracted 1767 pockmarks and found a linear pockmark depth-to-diameter ratio for pockmarks field-wide. Mean pockmark depth is 7.6 m and mean diameter is 84.8 m. Pockmark distribution is non-random, and nearly half of the field's pockmarks occur in chains. The most prominent chains are oriented semi-normal to the steepest gradient in Holocene sediment thickness. A descriptive model yields field-wide spatial statistics indicating that pockmarks are distributed in non-random clusters.Results enable quantitative comparison of pockmarks in fields worldwide as well as similar concave features, such as impact craters, dolines, or salt pools. Published by Elsevier B.V.
IntroductionFirst identified in muddy sediments of the Scotian Shelf (King and MacLean, 1970), pockmarks are seafloor depressions that are found worldwide in a variety of geologic settings (Hovland and Judd, 1988;Judd and Hovland, 2007). These craters can measure hundreds of meters in diameter, may occur in chains kilometers long and, where present in extensive fields, may dominate the seafloor surface (Fader, 1991;Rogers et al., 2006;Pilcher and Argent, 2007). Despite global distribution and general association with seafloor fluid escape, the mechanisms for pockmark formation and evolution remain uncertain (Ussler et al., 2003).Analysis of pockmark morphology and spatial distribution relative to antecedent geology and subsurface fluids (e.g., methane) can provide insight into fluid-migration pathways, pockmark field evolution, and possible mechanisms for pockmark generation and maintenance. In the absence of high-resolution seafloor bathymetry data, previous characterizations of entire pockmark fields relied upon visual interpretation of acoustic backscatter data for pockmark delineation, size statistics and spatial distribution (Fader, 1991;Kelley et al., 1994;Gontz et al., 2002;Rogers et al., 2006). Although acoustic backscatter data were the best available in the cited studies, interpreting size dimensions of concave features, such as pockmarks from these data is often ambiguous (Song, 2007). High-resolution bathymetry data collected by multibeam echosounder and swath sonar technologies enable the study of seafloor morphology to reach scales and resolutions similar to studies in subaerial geomorphology based on digital elevation models (DEMs) (Hughes Clarke et al., 1996). With one exception ) these new technologies have not been applied to an entire pockmark field. Instead, whole-field spatial and morphologic analysis has ebbed, r...