Postprocessing and corrections of bathymetry derived from sidescansonar systems: application with SeaMARC II

A procedure for postprocessing bathymetry data provided by a phase-measuring sidescan sonar system is presented. The data were collected with the SeaMARC II system, and are generally characterized by a high level of noise and uneven spatial sampling. Before any spatial filtering is applied, data are selected to remove most of the obvious artifacts and to retain instantaneous depth profiles whose slant ranges increase monotonically from a central location to the edges of the swath. An extrapolation scheme, patterned after a potential field, is proposed to fill gaps in the coverage or to extend the bathymetric swath to that of the corresponding sidescan image when regridding the data to a rectangular frame. To fill the near nadir gap typically found in these data, a specific interpolation methodology is developed that takes into account the slant range of the first bottom return as received by the sidescan sonar itself or by a shipboard echo-sounder. Spatial low-pass filtering is applied through convolutions with parabolic windows whose width is proportional to the footprint of the acoustic beam along track and roughly 1/8 of the swath width across track. Mismatches of contour lines between adjacent tracks are reduced through a statistical method design to correct systematic profile errors     

Toward remote seafloor classification using the angular response ofacoustic backscattering: a case study from multiple overlapping GLORIAdata

While the average seafloor backscatter strength within a narrow range of grazing angles can be used as a first-order classification tool, this technique often fails to distinguish seafloors of known differing geological character. In order to resolve such ambiguities, it is necessary to examine the variation in backscatter strength as a function of grazing angle. For this purpose, a series of multiply overlapping GLORIA sidescan sonar images (6.5 kHz) have been obtained in water depths ranging from 1000 to 2500 m. To constrain the placement of acoustic backscatter measurements and to measure the true impinging angle of the incident wave, the corresponding seafloor was simultaneously surveyed using the Seabeam multibeam system. As a result of the multiple overlap, the angular response of seafloor backscatter strength may be derived for regions much smaller than the swath width. By using the derived angular response of seafloor backscatter strength in regions for which sediment samples exist, an empirical seafloor classification scheme is proposed based on the shape, variance, and magnitude of the angular response. Because of the observed variability in the shape of the angular response with differing seafloor types, routine normalization of single-pass swath data to an equivalent single grazing angle image cannot be achieved. As a result, for the case of single-pass surveys, confident seafloor classification may only be possible for regions approaching the scale of the swath width     

Swath bathymetry with GLORIA

For many years, GLORIA has been producing sonar images of the deep ocean floor. In the mid-1980's, the SeaMARC II system came to prominence producing depth values as well as sonar images. The basic method compares the phases of the signals returning from the seafloor to two rows of transducers. The phase differences are converted into angles of arrival and together with the arrival times converted into range and depth values. This capability has now been added to the GLORIA system. The fact that GLORIA uses a 2s FM pulse means the backscattered reverberation can come from a strip of seafloor up to 1.5 km wide. To accommodate this, overlapping complex FFT's are used to produce a time-frequency matrix for the returning signals. In this matrix, a constant range feature appears as a diagonal. Phases are then calculated using a least-mean-squares estimate along diagonals. The main source of error and bias is due to surface reflection, and this is taken into account. The GLORIA swath bathymetry system was tested on two cruises and it was possible to produce depth contours with a good level of confidence. The total swath width was over eight water depths and would have been greater with a more favorable velocity profile. Comparison with other bathymetry data (such as multibeam systems) showed excellent correlation, having a standard deviation of only 4% of total water depth     

Remote estimation of surficial seafloor properties through the application Angular Range Analysis to multibeam sonar data

The variation of the backscatter strength with the angle of incidence is an intrinsic property of the seafloor, which can be used in methods for acoustic seafloor characterization. Although multibeam sonars acquire backscatter over a wide range of incidence angles, the angular information is normally neglected during standard backscatter processing and mosaicking. An approach called Angular Range Analysis has been developed to preserve the backscatter angular information, and use it for remote estimation of seafloor properties. Angular Range Analysis starts with the beam-by-beam time-series of acoustic backscatter provided by the multibeam sonar and then corrects the backscatter for seafloor slope, beam pattern, time varying and angle varying gains, and area of insonification. Subsequently a series of parameters are calculated from the stacking of consecutive time series over a spatial scale that approximates half of the swath width. Based on these calculated parameters and the inversion of an acoustic backscatter model, we estimate the acoustic impedance and the roughness of the insonified area on the seafloor. In the process of this inversion, the behavior of the model parameters is constrained by established inter-property relationships. The approach has been tested using a 300 kHz Simrad EM3000 multibeam sonar in Little Bay, NH. Impedance estimates are compared to in situ measurements of sound speed. The comparison shows a very good correlation, indicating the potential of this approach for robust seafloor characterization.     

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     

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.     

Benefits of swath mapping for the identification of marine habitats in the New Caledonia Economic Zone

The ZoNéCo programme is devoted to the evaluation of the marine resources of the Economic Zone of New Caledonia. The results are essentially dependent on the quality of the seafloor mapping. From 1993 to 1996, four geological and geophysical surveys using the EM12 DUAL multibeam echosounder provided swath-mapping and acoustic imagery data of the seafloor of selected sites on the northern and southern parts of the Norfolk ridge, the Loyalty basin, around the Loyalty islands and in the westernmost part of the Economic Zone of New Caledonia. The accuracy of these documents shows the morphology of the seafloor in detail and allows rocky substratum to be differentiated from muddy bottom. It allows favorable emplacements of future exploratory fishing surveys to be determined. The benefits of swath mapping are illustrated by the Halipro 2 deep sea trawling cruise (1996) which used the swath mapping data of ZoNéCo 1 cruise (1993) on the southern prolongation of the New Caledonia mainland and Loyalty Islands.     

Sidescan sonar image processing techniques

A four-step processing sequence is described to produce image mosaics from the various segments of a sidescanned acoustic imaging survey of a given seafloor area. Starting with data consisting for each ping of acoustic backscatter levels versus horizontal range across-track, median prefiltering is used first to reduce the influence of outliers on subsequent linear processes. Artifacts that are clearly unrelated to the backscattering properties of the seafloor are then isolated on a ping by ping basis through a spectral analysis that relies on a decomposition using Chebyshev polynomials to filter the low spatial frequency components of the image. Contrast enhancement is then achieved through an original implementation of the classical gray level histogram equalization technique by balancing local versus global histogram contributions. Pixels are mapped on a geographic grid taking due account of the geometry of the measurement and of the spacing between pings to minimize along-track smearing of features. Examples of results obtained with these processing techniques are given for SeaMARC II data recorded during a complete survey of Fieberling Guyot (32°.5 N, 128° W)     

Automatic Registration of TOBI Side-Scan Sonar and Multi-Beam Bathymetry Images for Improved Data Fusion

Deep towed side-scan sonar vehicles such as TOBI acquire high quality imagery of the seafloor with very high spatial resolution but poor locational accuracy. Fusion of the side-scan sonar data with bathymetry data from an independent source is often desirable to reduce ambiguity in geological interpretations, to aid in slant-range correction and to enhance seafloor representation. The main obstacle to fusion is accurate registration of the two datasets.The application of hierarchical chamfer matching to the registration of TOBI side-scan sonar images and multi-beam swath bathymetry is described. This matches low level features such as edges in the TOBI image, with corresponding features in a synthetic TOBI image created by simulating the flight of the TOBI vehicle through the bathymetry. The method is completely automatic, relatively fast and robust, and much easier than manual registration. It allows accurate positioning of the TOBI vehicle, enhancing its usefulness as a research tool. The method is illustrated by automatic registration of TOBI and multi-beam bathymetry data from the Mid-Atlantic Ridge.     

Stochastic modeling of seafloor morphology; Resolution oftopographic parameters by sea beam data

The authors explore the resolving power of an inversion algorithm which estimates five parameters of the seafloor covariance function from a single swath of multibeam echosounding data. The resolving power is evaluated as a function of the swath length, the orientation of ship track with respect to topographic grain, and the response width of the sounding system. The analysis is conducted by inverting sets of synthetic data with known statistics. The mean and standard deviation of the inverted parameters can be directly compared with the input parameters and the standard errors output from the inversion. Experiments show that resolution of the covariance parameters is strongly dependent on the number of characteristic lengths which are sampled. Root-mean-square seafloor height can be estimated to within ~15%, and anisotropic orientation to within ~5% (for a strong lineation), using track lengths as short as three characteristic lengths     

Processing and analysis of Simrad multibeam sonar data

The common approach to analysing data collected with multibeam and sidescan sonars is to visually interpret charts of contoured bathymetry and mosaics of seabed images. However, some of the information content is lost by processing the data into charts because this involves some averaging; the analysis might uncover more information if done on the data at an earlier stage in the processing. Motivated by this potential, I have created a software system which can be used to analyse data collected with Simrad EM1000 (shallow water) and EM12 (deep water) multibeam sonars, as well as to generate bathymetry contour charts and backscatter mosaics. The system includes data preprocessing, such as navigation filtering, depth filtering (removal of outlying values), and amplitude mapping using the multibeam bathymetry to correctly position image pixels across the swath. The data attributes that can be analysed include the orientation and slope of the seafloor, and the mean signal strength for each sounding. To determine bathymetry attributes such as slope, the soundings across a number of beams and across a series of pings are grouped and a least-squares plane fitted to them. Bathymetric curvature is obtained by detrending the grouped data using the least-squares plane and fitting a paraboloid to the residuals. The magnitudes and signs of the paraboloid's coefficients reveal depressions and hills and their orientations. Furthermore, the seafloor geology can be classified using a simple combination of these attributes. For example, flat-lying sediments can be classified where the backscatter, slope and curvature fall below specified values.     

Feasibility of interferometric swath bathymetry using GLORIA, along-range sidescan

The feasibility of adding an interferometric swath bathymetric system to GLORIA, a 6.6 kHz long-range sidescan sonar, is discussed. The size of GLORIA's low-frequency transducer arrays and towfish precludes significant modifications, but even without such changes bathymetric errors could be several tens of metres over a usable swath somewhat smaller than the normal GLORIA swath. A swath bathymetry based on GLORIA will have random errors depending strongly on wind speed, water depth, and swath width. Within the range of these parameters, root-mean-square bathymetry errors in the range of 1-100 m can be expected     

A processing requirement and resolution capability comparison ofside-scan and synthetic-aperture sonars

The processing requirements and resolution capabilities of both side-look sonar (SLS) and synthetic-aperture sonar (SAS) systems are outlined. Side-look sonar is presented as a real-beam imaging technique along with expressions for relevant system- and image-related parameters. Synthetic-aperture sonar is discussed, and the limitations imposed by the speed of sound in the ocean environment are identified. A specific side-look system (SeaMARC I) is presented under two configurations and comparable SAS designs are proposed. Based on the examples provided by the SeaMARC I system and the hypothetical SAS designs, it is shown that single-beam SAS systems can be designed to achieve area coverage rates comparable to single-beam side-scan systems, yet with improved azimuth resolution     

Shadow Enhancement in Synthetic Aperture Sonar Using Fixed Focusing

A shadow cast by an object on the seafloor is important information for target recognition in synthetic aperture sonar (SAS) images. Synthetic aperture imaging causes a fundamental limitation to shadow clarity because the illuminator is moved during the data collection. This leads to a blend of echo and shadow, or geometrical fill-in in the shadow region. The fill-in is most dominant for widebeam synthetic aperture imaging systems. By treating the shadow as a moving target and compensating for the motion during the synthetic aperture imagery, we avoid the geometrical shadow fill-in. We show this to be equivalent to fixing the focus at the range of the shadow caster. This novel technique, referred to as fixed focus shadow enhancement (FFSE) can be used directly as an imaging method on hydrophone data or as a postprocessing technique on the complex SAS image. We demonstrate the FFSE technique on simulated data and on real data from a rail-based SAS, and on two different SAS systems operated on a HUGIN autonomous underwater vehicle.      

Structural development of central North Fiji Basin triple junction

The structural development of the central North Fiji Basin triple junction is revealed in SeaMARC II sidescan imagery. All three limbs have unique morphotectonic characteristics and the adjacent sea-floor fabric is oblique to the limb axes. The triple junction is interpreted to be a recent phenomenon and to have formed by a discrete jump to its present site.     

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     

Acoustic imagery using a multibeam bathymetric system

Historically, measurement and collection of deep‐ocean acoustic imagery are accomplished by towed sidescan systems. Recently, work has been performed to extract acoustic imagery from current hull‐mounted wide‐swath bathymetric sonars with minimal hardware modification. Past work of deriving acoustic imagery from swath sonars has been performed primarily with SeaBeam's sixteen 2^2/3 ° preformed beams. The Navy is investigating the feasibility of extracting an acoustic image from the Sonar Array Survey Systems (SASS), a high‐resolution (1^o beams) wide‐fan (90°) bathymetric system. Due to the large data volume (approximately 1 MB per ping), SASS normally discards the raw acoustic returns once bathymetry is calculated. In early 1991 the Naval Air Development Center (NADC) installed the hardware on board the USNS Maury to capture and record the raw acoustic signal (inphase and quadrature) from the SASS's 144 hydrophones for later inversion to a backscatter image. Preliminary qualitative mosaics of the sidescan images show promising results and warrant further development.     

High-resolution mapping of large gas emitting mud volcanoes on the Egyptian continental margin (Nile Deep Sea Fan) by AUV surveys

Two highly active mud volcanoes located in 990–1,265 m water depths were mapped on the northern Egyptian continental slope during the BIONIL expedition of R/V Meteor in October 2006. High-resolution swath bathymetry and backscatter imagery were acquired with an autonomous underwater vehicle (AUV)-mounted multibeam echosounder, operating at a frequency of 200 kHz. Data allowed for the construction of ~1 m pixel bathymetry and backscatter maps. The newly produced maps provide details of the seabed morphology and texture, and insights into the formation of the two mud volcanoes. They also contain key indicators on the distribution of seepage and its tectonic control. The acquisition of high-resolution seafloor bathymetry and acoustic imagery maps with an AUV-mounted multibeam echosounder fills the gap in spatial scale between conventional multibeam data collected from a surface vessel and in situ video observations made from a manned submersible or a remotely operating vehicle.     

High-resolution seafloor mapping using the DSL-120 sonar system: Quantitative assessment of sidescan and phase-bathymetry data from the Lucky Strike segment of the Mid-Atlantic Ridge

High-resolution, side-looking sonar data collected near the seafloor (100 m altitude) provide important structural and topographic information for defining the geological history and current tectonic framework of seafloor terrains. DSL-120 kHz sonar data collected in the rift valley of the Lucky Strike segment of the Mid-Atlantic Ridge near 37° N provide the ability to quantitatively assess the effective resolution limits of both the sidescan imagery and the computed phase-bathymetry of this sonar system. While the theoretical, vertical and horizontal pixel resolutions of the DSL-120 system are <1 m, statistical analysis of DSL-120 sonar data collected from the Lucky Strike segment indicates that the effective spatial resolution of features is 1–2 m for sidescan imagery and 4 m for phase-bathymetry in the seafloor terrain of the Mid-Atlantic Ridge rift valley. Comparison of multibeam bathymetry data collected at the sea-surface with deep-tow DSL-120 bathymetry indicates that depth differences are on the order of the resolution of the multibeam system (10–30 m). Much of this residual can be accounted for by navigational mismatches and the higher resolving ability of the DSL-120 data, which has a bathymetric footprint on the seafloor that is 20 times smaller than that of hull-mounted multibeam at these seafloor depths (2000 m). Comparison of DSL-120 bathymetry with itself on crossing lines indicates that residual depth values are ±20 m, with much of that variation being accounted for by navigational errors. A DSL-120 survey conducted in 1998 on the Juan de Fuca Ridge with better navigation and less complex seafloor terrain had residual depth values half those of the Lucky Strike survey. The quality of the bathymetry data varies as a function of position within the swath, with poorer data directly beneath the tow vehicle and also towards the swath edges.Variations in sidescan amplitude observed across the rift valley and on Lucky Strike Seamount correlate well with changes in seafloor roughness caused by transitions from sedimented seafloor to bare rock outcrops. Distinct changes in sonar backscatter amplitude were also observed between areas covered with hydrothermal pavement that grade into lava flows and the collapsed surface of the lava lake in the summit depression of Lucky Strike Seamount. Small features on the seafloor, including volcanic constructional features (e.g., small cones, haystacks, fissures and collapse features) and hydrothermal vent chimneys or mounds taller than 2 m and greater than 9 m² in surface area, can easily be resolved and mapped using this system. These features at Lucky Strike have been confirmed visually using the submersible Alvin, the remotely operated vehicle Jason, and the towed optical/acoustic mapping system Argo II.