Geometric corrections on sidescan sonar images based on bathymetry. Application with SeaMARC II and Sea Beam data

Acoustic backscatter images of the seafloor obtained with sidescan sonar systems are displayed most often using a flat bottom assumption. Whenever this assumption is not valid, pixels are mapped incorrectly in the image frame, yielding distorted representations of the seafloor. Here, such distortions are corrected by using an appropriate representation of the relief, as measured by the sonar that collected the acoustic backscatter information. In addition, all spatial filtering operations required in the pixel relocation process take the sonar geometry into account. Examples of the process are provided by data collected in the Northeastern Pacific over Fieberling Guyot with the SeaMARC II bathymetric sidescan sonar system and the Sea Beam multibeam echo-sounder. The nearly complete (90%) Sea Beam bathymetry coverage of the Guyot serves as a reference to quantify the distortions found in the backscatter images and to evaluate the accuracy of the corrections performed with SeaMARC II bathymetry. As a byproduct, the processed SeaMARC II bathymetry and the Sea Beam bathymetry adapted to the SeaMARC II sonar geometry exhibit a 35m mean-square difference over the entire area surveyed.On leave at the Naval Research Laboratory, Code 7420, Washington D.C. 20375-5350.     

Simultaneous operation of the Sea Beam multibeam echo-sounder andthe SeaMARC II bathymetric sidescan sonar system

An experiment aboard the Scripps Institution of Oceanography's RV Thomas Washington has demonstrated the seafloor mapping advantages to be derived from combining the high-resolution bathymetry of a multibeam echo-sounder with the sidescan acoustic imaging plus wide-swath bathymetry of a shallow-towed bathymetric sidescan sonar. To a void acoustic interference between the ship's 12-kHz Sea Beam multibeam echo-sounder and the 11-12-kHz SeaMARC II bathymetric sidescan sonar system during simultaneous operations, Sea Beam transmit cycles were scheduled around SeaMARC II timing events with a sound source synchronization unit originally developed for concurrent single-channel seismic, Sea Beam, and 3.5-kHz profile operations. The scheduling algorithm implemented for Sea Beam plus SeaMARC II operations is discussed, and the initial results showing their combined seafloor mapping capabilities are presented     

Seafloor acoustic remote sensing with multibeam echo-sounders and bathymetric sidescan sonar systems

This paper examines the potential for remote classification of seafloor terrains using a combination of quantitative acoustic backscatter measurements and high resolution bathymetry derived from two classes of sonar systems currently used by the marine research community: multibeam echo-sounders and bathymetric sidescans sonar systems. The high-resolution bathymetry is important, not only to determine the topography of the area surveyed, but to provide accurate bottom slope corrections needed to convert the arrival angles of the seafloor echoes received by the sonars into true angles of incidence. An angular dependence of seafloor acoustic backscatter can then be derived for each region surveyed, making it possible to construct maps of acoustic backscattering strength in geographic coordinates over the areas of interest. Such maps, when combined with the high-resolution bathymetric maps normally compiled from the data output by the above sonar systems, could be very effective tools to quantify bottom types on a regional basis, and to develop automatic seafloor classification routines.     

Recommendations for improved and coherent acquisition and processing of backscatter data from seafloor-mapping sonars

Multibeam echosounders are becoming widespread for the purposes of seafloor bathymetry mapping, but the acquisition and the use of seafloor backscatter measurements, acquired simultaneously with the bathymetric data, are still insufficiently understood, controlled and standardized. This presents an obstacle to well-accepted, standardized analysis and application by end users. The Marine Geological and Biological Habitat Mapping group (Geohab.org) has long recognized the need for better coherence and common agreement on acquisition, processing and interpretation of seafloor backscatter data, and established the Backscatter Working Group (BSWG) in May 2013. This paper presents an overview of this initiative, the mandate, structure and program of the working group, and a synopsis of the BSWG Guidelines and Recommendations to date. The paper includes (1) an overview of the current status in sensors and techniques available in seafloor backscatter data from multibeam sonars; (2) the presentation of the BSWG structure and results; (3) recommendations to operators, end-users, sonar manufacturers, and software developers using sonar backscatter for seafloor-mapping applications, for best practice methods and approaches for data acquisition and processing; and (4) a discussion on the development needs for future systems and data processing. We propose for the first time a nomenclature of backscatter processing levels that affords a means to accurately and efficiently describe the data processing status, and to facilitate comparisons of final products from various origins.     

基于回波时间的多波束测深与声纳数据匹配方法

以多波束精确的水深数据为参照源,采用原始回波时间对多波束测深数据与其同源声纳数据进行匹配,从而获得高精度和高分辨率的海底影像数据,并避免了传统声纳图像处理过程中斜距改正所带来的几何形变。匹配结果采用光照图输出,并与三维水深图、原始声纳图像和CARIS处理后的声纳图像进行比较分析。该方法有效地提高了多波束数据的利用率,增强了对海底地形的探测分辨率。     

多波束反向散射强度数据处理研究

在探讨多波束测深系统反向散射强度与海底底质类型的关系基础上,研究影响反向散射强度的各种因素,主要分析了海底地形起伏、中央波束区反射信号对反向散射强度的影响,并给出了消除这些影响的方法;将处理后的“纯”反向散射强度数据镶嵌生成海底声像图,为海底底质类型划分以及地貌解译提供了基础数据和辅助判读依据.     

Constrained iterative deconvolution applied to SeaMARC I sidescansonar imagery

Images collected by any sidescan sonar system represent the convolution of the acoustic beam pattern of the instrument with the true echo amplitude distribution over the seafloor. At typical low speeds, the 1.7° beam width of SeaMARC I (seafloor mapping and remote characterization) results in multiple insonification of individual targets, particularly at the outside of the swath. A nonlinearly constrained iterative deconvolution technique developed for radar applications can be applied to SeaMARC I imagery to reduce the effect of the beam pattern and equalize the spectral content of the image across the swath. Since the deconvolution is implemented in the along-track direction, the registration of individual scan lines must be precisely corrected before the operator is applied. The deconvolution operator must be modeled to account for beam shape, vehicle speed, swath width, slant range, and ping rate. The method is numerically stable and increases the effective resolution of the image, but results in some loss of dynamic range. The technique is applied to target recognition and imagery from volcanic terrains of the central Juan de Fuca Ridge     

Processing Multibeam Backscatter Data

A new highly precise source of data has recently become available using multibeam sonar systems in hydrography. Multibeam sonar systems can provide hydrographic quality depth data as well as high-resolution seafloor sonar images. We utilize the seafloor backscatter strength data of each beam from multibeam sonar and the automatic classification technology so that we can get the seafloor type identification maps. In this article, analyzing all kinds of error effects in backscatter strength, data are based on the relationship between backscatter strength and seafloor types. We emphasize particularly analyzing the influences of local bottom slope and near nadir reflection in backscatter strength data. We also give the correction algorithms and results of these two influent factors. After processing the raw backscatter strength data and correcting error effects, we can get processed backscatter strength data which reflect the features of seafloor types only. Applying the processed backscatter strength data and mosaicked seafloor sonar images, we engage in seafloor classification and geomorphy interpretation in future research.     

Processing Multibeam Backscatter Data

A new highly precise source of data has recently become available using multibeam sonar systems in hydrography. Multibeam sonar systems can provide hydrographic quality depth data as well as high-resolution seafloor sonar images. We utilize the seafloor backscatter strength data of each beam from multibeam sonar and the automatic classification technology so that we can get the seafloor type identification maps. In this article, analyzing all kinds of error effects in backscatter strength, data are based on the relationship between backscatter strength and seafloor types. We emphasize particularly analyzing the influences of local bottom slope and near nadir reflection in backscatter strength data. We also give the correction algorithms and results of these two influent factors. After processing the raw backscatter strength data and correcting error effects, we can get processed backscatter strength data which reflect the features of seafloor types only. Applying the processed backscatter strength data and mosaicked seafloor sonar images, we engage in seafloor classification and geomorphy interpretation in future research.     

Differential phase estimation with the SeaMARCII bathymetricsidescan sonar system

A maximum-likelihood estimator is used to extract differential phase measurements from noisy seafloor echoes received at pairs of transducers mounted on either side of the SeaMARC II bathymetric sidescan sonar system. Carrier frequencies for each side are about 1 kHz apart, and echoes from a transmitted pulse 2 ms long are analyzed. For each side, phase difference sequences are derived from the full complex data consisting of base-banded and digitized quadrature components of the received echoes. With less bias and a lower variance, this method is shown to be more efficient than a uniform mean estimator. It also does not exhibit the angular or time ambiguities commonly found in the histogram method used in the SeaMARC II system. A figure for the estimation uncertainty of the phase difference is presented, and results are obtained for both real and simulated data. Based on this error estimate and an empirical verification derived through coherent ping stacking, a single filter length of 100 ms is chosen for data processing applications     

Fault scarp identification in side-scan sonar and bathymetry images from the Mid-Atlantic Ridge using wavelet-based digital filters

Digital filters designed using wavelet theory are applied to high resolution deep-towed side-scan sonar data from the median valley walls, crestal mountains, and flanks of the Mid-Atlantic Ridge at 29°10 N. With proper tuning, the digital filters are able to identify the location, orientation, length, and width of highly reflective linear features in sonar images. These features are presumed to represent the acoustic backscatter from axis-facing normal faults. The fault locations obtained from the digital filters are well correlated with visual geologic interpretation of the images. The side-scan sonar images are also compared with swath bathymetry from the same area. The digitally filtered bathymetry images contain nine of the eleven faults identified by eye in the detailed geologic interpretation of the side-scan data. Faults with widths (measured perpendicular to their strike) of less than about 150 m are missed in the bathymetry analysis due to the coarser resolution of these data. This digital image processing technique demonstrates the potential of wavelet-based analysis to reduce subjectivity and labor involved in mapping and analyzing topographic features in side-scan sonar and bathymetric image data.     

Inhomogeneous substrate analysis using EM300 backscatter imagery

Backscatter reflectivity from multibeam echo-sounders provides a powerful tool to efficiently characterize seafloor substrates. A comprehensive EM300 bathymetric and backscatter survey has been completed of Cook Strait, in central New Zealand. This paper presents a detailed analysis of the realtime corrections applied to the raw EM300 multibeam data and additional corrections required to compute angular variations of the backscatter strength. The corrections, including the local absorption coefficient, the influence of seafloor topography and sound refraction in the water column, are determined for different Cook Strait seafloor substrates. Modifying MB-System software code, we extracted the backscatter signal parameters in order to quantify the raw backscatter strength and apply additional processing. Profiles of backscatter strength versus incidence angle were computed for a variety of sites characterized by flat seafloor and homogeneous substrates, and for which ground-truth data were available. For each homogeneous site, different but characteristic backscatter profiles are observed that can be interpreted in terms of sediment facies. To analyze heterogeneous substrates, we present a statistical technique, based on a 3-dimensional distribution of (incidence angle, backscatter strength) couples that preserves geological information of the substrate components. This analysis, using backscatter data acquired on a submarine volcano, north of New Zealand, clearly differentiates soft sediments and lava flows within a heterogeneous substrate.     

Multisource multibeam backscatter data: developing a strategy for the production of benthic habitat maps using semi-automated seafloor classification methods

The establishment of multibeam echosounders (MBES) as a mainstream tool in ocean mapping has facilitated integrative approaches towards nautical charting, benthic habitat mapping, and seafloor geotechnical surveys. The bathymetric and backscatter information generated by MBES enables marine scientists to present highly accurate bathymetric data with a spatial resolution closely matching that of terrestrial mapping, and can generate customized thematic seafloor maps to meet multiple ocean management needs. However, when a variety of MBES systems are used, the creation of objective habitat maps can be hindered by the lack of backscatter calibration, due for example, to system-specific settings, yielding relative rather than absolute values. Here, we describe an approach using object-based image analysis to combine 4 non-overlapping and uncalibrated (backscatter) MBES coverages to form a seamless habitat map on St. Anns Bank (Atlantic Canada), a marine protected area hosting a diversity of benthic habitats. The benthoscape map was produced by analysing each coverage independently with supervised classification (k-nearest neighbor) of image-objects based on a common suite of 7 benthoscapes (determined with 4214 ground-truthing photographs at 61 stations, and characterized with backscatter, bathymetry, and bathymetric position index). Manual re-classification based on uncertainty in membership values to individual classes—especially at the boundaries between coverages—was used to build the final benthoscape map. Given the costs and scarcity of MBES surveys in offshore marine ecosystems—particularly in large ecosystems in need of adequate conservation strategies, such as in Canadian waters—developing approaches to synthesize multiple datasets to meet management needs is warranted.     

Mapping marine habitats with high resolution sidescan sonar

The application of marine geophysics and GIS techniques to the characterization of benthic habitats has increased the ability of fisheries managers to assess distribution and habitat types beyond common practices. We report upon a 150 kHz sidescan sonar survey offshore of Kruzof Island, Alaska undertaken to characterize rockfish (Sebastes) habitat. Using GIS, MapGrafix and Map¹Factory we determined the percentage of seafloor cover that exists in our survey area. Bathymetry in the study area was determined with sidescan interferometry. All XYZ data were gridded using Surfer and plotted in shaded relief, bathymetric contour, and 3-dimensional formats. Contoured bathymetry was used as an over-lay in MapGrafix. Small sub-areas were extracted from the bathymetric data for closer study, and gridded in Surfer. Areas of the mosaic where backscatter patterns were not distinct were verified with hand samples and video collected with the submersible Delta. The use of submersibles for verification of interpreted lithologies and surface textures enables a high degree of accuracy for the interpretations. Lithotypes were lumped into larger groups based on morphology and fish associations with different morphologies verified using the submersible. The accuracy of digital maps from high-resolution sidescan sonar data allows a close quantification of the areal extents of these important features, directing the application of management strategies to critical areas.     

A SeaMARC II survey of recent submarine volcanism near Easter Island

SeaMARC II side-scan sonar data reveal that a large area of seafloor north and west of Easter Island has been disrupted by recent submarine volcanism. A large volcanic area begins approximately 60 km WNW of the island and extends for over 130 km to the west. The volcanic field is characterized by high backscatter intensity in the side-scan sonar records and is elevated 400–1000 m above the N-S seafloor fabric that surrounds it. This field, the Abu Volcanic Field, covers at least 2500 km² and appears to consist of recent lava flows and small volcanoes. Backscatter intensity of the Abu Volcanic Field is similar to that of the adjacent ridge flank which is less than 0.4 Ma, suggesting a similar age for its formation. Two additional areas of high backscatter immediately north of Easter Island cover a combined area of over 300 km². The sidescan sonar records show that these features are clearly of volcanic origin and are not debris flows from the nearby island. The flows are nearly 300 m thick and are morphologically similar to subaerial pahoehoe lava shields. Their high backscatter indicates that they are also the products of relatively recent submarine volcanic activity. The presence of these large areas of recent volcanism in the vicinity of Easter Island has important implications for the various models that have been proposed to explain the origin of the Easter Seamount Chain. In addition, the similar ages of Easter Island and the Easter Microplate suggest that the presence of a hotspot near or beneath this fast-spreading portion of the East Pacific Rise about 4.5 m.y. ago may have initiated the large-scale rift propagation that created the microplate.     

Processing of high-frequency multibeam echo sounder data for seafloor characterization

Processing simultaneous bathymetry and backscatter data, multibeam echosounders (MBESs) show promising abilities for remote seafloor characterization. High-frequency MBESs provide a good horizontal resolution, making it possible to distinguish fine details at the water-seafloor interface. However, in order to accurately measure the seafloor influence on the backscattered energy, the recorded sonar data must first be processed and cleared of various artifacts generated by the sonar system itself. Such a preprocessing correction procedure along with the assessment of its validity limits is presented and applied to a 95-kHz MBES (Simrad EM 1000) data set. Beam pattern effects, uneven array sensitivities, and inaccurate normalization of the ensonified area are removed to make possible further quantitative analysis of the corrected backscatter images. Unlike low-frequency data where the average backscattered energy proves to be the only relevant feature for discriminating the nature of the seafloor, high-frequency MBES backscatter images exhibit visible texture patterns. This additional information involves different statistical distributions of the backscattered amplitudes obtained from various seafloor types. Non-Rayleigh statistics such as K-distributions are shown to fit correctly the skewed distributions of experimental high-frequency data. Apart from the effect of the seafloor micro-roughness, a statistical model makes clear a correlation between the amplitude statistical distributions and the signal incidence angle made available by MBES bathymetric abilities. Moreover, the model enhances the effect of the first derivative of the seafloor backscattering strength upon statistical distributions near the nadir and at high incidence angles. The whole correction and analysis process is finally applied to a Simrad EM 1000 data set.     

A framework to quantify uncertainties of seafloor backscatter from swath mapping echosounders

Multibeam echosounders (MBES) have become a widely used acoustic remote sensing tool to map and study the seafloor, providing co-located bathymetry and seafloor backscatter. Although the uncertainty associated with MBES-derived bathymetric data has been studied extensively, the question of backscatter uncertainty has been addressed only minimally and hinders the quantitative use of MBES seafloor backscatter. This paper explores approaches to identifying uncertainty sources associated with MBES-derived backscatter measurements. The major sources of uncertainty are catalogued and the magnitudes of their relative contributions to the backscatter uncertainty budget are evaluated. These major uncertainty sources include seafloor insonified area (1–3 dB), absorption coefficient (up to >?6 dB), random fluctuations in echo level (5.5 dB for a Rayleigh distribution), and sonar calibration (device dependent). The magnitudes of these uncertainty sources vary based on how these effects are compensated for during data acquisition and processing. Various cases (no compensation, partial compensation and full compensation) for seafloor insonified area, transmission losses and random fluctuations were modeled to estimate their uncertainties in different scenarios. Uncertainty related to the seafloor insonified area can be reduced significantly by accounting for seafloor slope during backscatter processing while transmission losses can be constrained by collecting full water column absorption coefficient profiles (temperature and salinity profiles). To reduce random fluctuations to below 1 dB, at least 20 samples are recommended to be used while computing mean values. The estimation of uncertainty in backscatter measurements is constrained by the fact that not all instrumental components are characterized and documented sufficiently for commercially available MBES. Further involvement from manufacturers in providing this essential information is critically required.     

Analysis of seafloor backscatter strength dependence on the survey azimuth using multibeam echosounder data

The sediment backscatter strength measured by multibeam echosounders is a key feature for seafloor mapping either qualitative (image mosaics) or quantitative (extraction of classifying features). An important phenomenon, often underestimated, is the dependence of the backscatter level on the azimuth angle imposed by the survey line directions: strong level differences at varying azimuth can be observed in case of organized roughness of the seabed, usually caused by tide currents over sandy sediments. This paper presents a number of experimental results obtained from shallow-water cruises using a 300-kHz multibeam echosounder and specially dedicated to the study of this azimuthal effect, with a specific configuration of the survey strategy involving a systematic coverage of reference areas following “compass rose” patterns. The results show for some areas a very strong dependence of the backscatter level, up to about 10-dB differences at intermediate oblique angles, although the presence of these ripples cannot be observed directly—neither from the bathymetry data nor from the sonar image, due to the insufficient resolution capability of the sonar. An elementary modeling of backscattering from rippled interfaces explains and comforts these observations. The consequences of this backscatter dependence upon survey azimuth on the current strategies of backscatter data acquisition and exploitation are discussed.     

Mapping benthic habitat in regions of gradational substrata: An automated approach utilising geophysical,geological, and biological relationships

The relationship between acoustic backscatter, sediment characteristics and benthic habitat is examined using high-resolution sidescan sonar data collected at the Loch Linnhe artificial reef site on the west coast of Scotland. The site is typical for the continental shelf of NW Europe, with a mix of seabed environments from muddy to coarse, stony substrata on a 10–100 m length scale. A sidescan sonar mosaic was produced and classified according to derived backscatter parameters (mean, median and standard deviation of the backscatter values) using an unsupervised classification procedure. The accuracy of the final classified map was assessed by comparison with a ground-truthing survey in which the biological habitat was derived from underwater video footage. The sidescan correctly predicted seabed surface characteristics of observed biological habitat with 78% accuracy. A second, and more challenging test of the acoustic data to correctly predict biological habitat was made by comparing it with data from 21 grab sampling stations. These stations were divided into three groups using multivariate statistical techniques based on their backscatter properties. Benthic assemblage structure was found to be significantly distinct between the high and low, and the medium and low backscatter stations. There was a low to moderate but significant correlation between the multivariate patterns of acoustic backscatter, benthic assemblage structure, and particle size distribution. The work shows that even in areas with subtle and gradational changes in substratum, the sidescan was able to predict biological community with an acceptable accuracy.     

Geomorphology,acoustic backscatter,and processes in Santa Monica Bay from multibeam mapping

Santa Monica Bay was mapped in 1996 using a high-resolution multibeam system, providing the first substantial update of the submarine geomorphology since the initial compilation by Shepard and Emery [(1941) Geol. Soc. Amer. Spec. Paper 31]. The multibeam mapping generated not only high-resolution bathymetry, but also coregistered, calibrated acoustic backscatter at 95 kHz. The geomorphology has been subdivided into six provinces; shelf, marginal plateau, submarine canyon, basin slope, apron, and basin. The dimensions, gradients, and backscatter characteristics of each province is described and related to a combination of tectonics, climate, sea level, and sediment supply. Fluctuations of eustatic sea level have had a profound effect on the area; by periodically eroding the surface of Santa Monica plateau, extending the mouth of the Los Angeles River to various locations along the shelf break, and by connecting submarine canyons to rivers. A wetter glacial climate undoubtedly generated more sediment to the rivers that then transported the increased sediment load to the low-stand coastline and canyon heads. The trends of Santa Monica Canyon and several bathymetric highs suggest a complex tectonic stress field that has controlled the various segments. There is no geomorphic evidence to suggest Redondo Canyon is fault controlled. The San Pedro fault can be extended more than 30 km to the northwest by the alignment of a series of bathymetric highs and abrupt changes in direction of channel thalwegs.