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.     

Geoacoustic inversion of broad-band acoustic data in shallow water on a single hydrophone

The inversion of broad-band low-frequency acoustic signals received on sparse arrays can lead to robust and efficient estimations of sea-bed properties. This paper describes a shallow-water geoacoustic inversion scheme based on the use of a model-based matched-impulse response on a single hydrophone. Results from the INTIMATE'96 experiment on the Portuguese shelf break are reviewed. In order to minimize the effects of strong time variability due to internal tides, only the time-stable waterborne bottom-surface reflected arrivals are exploited. A quasi-linear inversion algorithm is first applied to refine the geometry of the experiment. Then, inversion of bottom parameters is performed with an objective function that only makes use of the bottom-surface reflected arrivals' amplitudes. The experimental results show that broad-band transmissions (300-800 Hz) received on a single hydrophone, combined with the use of a simple eigenray code, are sufficient to correctly resolve geometrical parameters and bottom features. The analysis of the reflection coefficients both on simulated and real data helps to understand the validity of the inverted parameters and to derive the basis of an equivalent medium concept for geoacoustic inversion based on a "through-the-sensor" approach.     

声相关海流剖面仪测速研究

声相关海流剖面仪(ACCP)适合深海测速,能够测得相对海底的绝对速度。文章在前人工作基础上,以声纳方程为基础,构造海底混响回波仿真模型,并对基于K irchhoff近似和利用微扰法的R ay le igh-R ice近似相结合的散射模型赋予新的物理意义,并将其应用到声相关测速仿真研究过程中,使仿真更加接近于海洋的实际环境。文章利用H ilbert变换进行包络检波来减小测速误差,并且根据接收信噪比得到了不同航速时能够达到最大底跟踪的距离,这与参考文献给出的结果是一致的。     

Estimation of bottom scattering strength from measured and modeledmid-frequency sonar reverberation levels

Hamilton-type geoacoustic models were developed for Area Foxtrot, a shallow water test bed south of Long Island, for emerging active sonar systems where the surface sediment type is highly spatially variable. Reverberation levels (RL) were modeled using the finite-element parabolic equation (FEPE) propagation model to augment the generic sonar model (GSM) propagation model because the bottom loss model in GSM did not estimate transmission loss (TL) accurately in shallow water. FEPE estimates reveal that there is a greater than 15 dB difference between TL for sand and that for silt-day sediments in Area Foxtrot. The comparison between modeled RL and measured RL (from a 1991 active sonar exercise) enabled bottom scattering strength kernels to be developed for Area Foxtrot. Bottom scattering strength was found to be a function of sediment type. Hard sand sediment has a bottom scattering strength which obeys Lambert's law (sin² &thetas;) while that of silt-clay sediment is consistent with sub-bottom volume scattering (sine). The RLs in Area Foxtrot are azimuth-dependent and are a function of TL and bottom scattering strength (and hence bottom sediment type). Sonar beams steered towards the hard sand show higher RLs than for silt-clay, and knowledge of the sediment type and its spatial variation must be known to model RL accurately. A method to determine sediment type using measured RLs and RL slopes is given     

Accuracy of the spatial representation of the seafloor with bathymetric sidescan sonars

When isobath maps of the seafloor are constructed with a bathymetric sidescan sonar system the position of each sounding is derived from estimates of range and elevation. The location of each pixel forming the acoustic backscatter image is calculated from the same estimates. The accuracy of the resulting maps depends on the acoustic array geometry, on the performances of the acoustic signal processing, and on knowledge of other parameters including: the platform's navigation, the sonar transducer's attitude, and the sound rays' trajectory between the sonar and the seafloor. The relative importance of these factors in the estimation of target location is assesed. The effects of the platform motions (e.g. roll, pitch, yaw, sway, surge and heave) and of the uncertainties in the elevation angle measurements are analyzed in detail. The variances associated with the representation (orientation and depth) of a plane, rectangular patch of the seafloor are evaluated, depending on the geometry of the patch. The inverse problem is addressed. Its solution gives the lateral dimensions of the spatial filter that must be applied to the bathymetric data to obtain specified accuracies of the slopes and depths. The uncertainty in the estimate of elevation angle, mostly due to the acoustic noise, is found to bring the main error contribution in across-track slope estimates. It can also be critical for along-track slope estimates, overshadowing error contributions due to the platform's attitude. Numerical examples are presented.On leave at the Naval Research Laboratory, Code 7420, Washington D.C. 20375-5350, U.S.A.     

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.     

Seafloor profiling by a wideband sonar: simulation,frequency-response optimization, and results of a brief sea test

An ahead-looking probe of some kind, optical or acoustic, is critical when one is attempting seafloor exploration from a mobile platform. A single-frequency, split aperture sonar system can be used for this purpose, but a wideband monopulse sonar offers many advantages. It computes a running estimate of the vertical directional cosine of the source of the echo, and can thus reveal the positions of multiple wave scatterers as long as their echoes can still be time resolved. Theoretical studies of its performance have been made previously, but were directly applicable only to extremely simple seafloor geometries. A new time-domain digital simulation that largely circumvents this limitation has been developed. The simulation also provides a means for testing the theory and optimizing system parameters. The reverberation model does not account for some features of acoustic backscattering such as diffraction, but it is believed to be adequate for the investigation of most signal processing aspects of the sonar system. The theory of the simulation is developed and several examples are presented and discussed. In addition, some preliminary results are presented from a sea test that used the air-sea interface as a surrogate seafloor     

Simulating Realistic-Looking Sediment Ripple Fields

Sandy sediment ripples impact sonar performance in coastal waters through Bragg scattering. Observations from data suggest that sandy ripple elevation relative to the mean seafloor as a function of the horizontal coordinates is not Gaussian distributed; specifically, peak amplitude fading over space associated with a random Gaussian process is largely absent. Such a non-Gaussian nature has implications for modeling acoustic scattering from, and penetration into, sediments. An algorithm is developed to generate ripple fields with a given power spectrum; these fields have non-Gaussian statistics and are visually consistent with data. Higher order statistics of these ripple fields and their implications to sonar detection are discussed.      

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 Geomorphology and Processes on the Western Continental Margin of India based on Multibeam Acoustic Backscatter Data Analysis

An assessment of the multibeam sonar data of the central Western Continental Margins of India has been carried out to evaluate the seafloor geomorphology and processes by examining the geomorphological attributes e.g., slope, sediments, structures, etc. associated with geomorphic features. The variation in relief and the features located in the region have been mapped and interpreted collectively by utilizing several geospatial mapping tools. The backscatter strength across the area, apparently congruent with the local relief, has helped to examine the sediment movement on the seafloor. The prominent features found in the region include faults, pockmarks, mounds, submarine terraces, and submerged fossil reefs. Several areas with varying topography engender comparable fractal dimension at short scale breaks, and the probability density functions (PDFs) utilizing backscatter data depicting overlapping classes. The present study highlights how fractals and scale break parameters can be utilized to determine the seafloor processes and associated sedimentological dynamics in a complex geographical environment with strong bottom currents, seasonal upwelling, and faulted structure. The role and impact of the various geomorphic processes on the reworking of sediment movement and the overall progression of the seafloor morphology has been revealed for the first time in this part of the ocean bottom.     

Study on the high precision acoustic measurement techniques for determining temperature field around seafloor hydrothermal vent

This paper presents the basis of acoustic method used for temperature field measurement of seafloor hydrothermal vent and two techniques of the parabolic interpolation and the bending compensation of propagation paths of acoustic signal are introduced. Experimental research is performed to exactly rebuild the temperature field around hot springs on the floor of Qiezishan Lake, Yunnan, China. The accuracy of the travel time estimation has been improved based on the aforementioned technique and method. At the same time, by comparison of the results of temperature field with different means, the max absolute error, the maximum relative error and the root mean square error are given. It shows that the technique and the method presented in the paper can be applied to the temperature field measurement detector around the seafloor hydrothermal vent. It also has a good accuracy.     

Measurement of sound transmission and signal gain in the complex Strait of Korea

The experiment, The Acoustic Characterization Test III, was conducted in the oceanographically complex Strait of Korea to accurately measure the sound transmission under known environmental conditions. Geoacoustic profiles derived from geophysical measurements, measured bathymetry, and sound-speed profiles were the basis for range dependent parabolic equation (PE) calculations. Agreement between measured and calculated transmission loss was obtained with an attenuation profile in the near water-sediment interface layer with a dependence on frequency to the 1.8 power consistent with measurements in other sand-silt areas. Since the environment was oceanographically complex and the shipping noise levels were high, the coherency of the sound transmission was estimated using relative signal gain (RSG). RSG was taken as the difference between the gain calculated with PE and measured with the array and at longer ranges and higher frequencies was found to be approximately -2 dB with a signal gain coefficient of variation of 5%. This RSG degradation, attributed to the random signal phase fluctuations resulting from scattering from the surfaces and volume of the waveguide, yielded using a Gaussian coherence function a spatial coherence length of 30/spl lambda/ @ 400 Hz-40 km. In addition, high resolution imaging of five targets with two bottom mounted arrays illustrate the achievable performance of low-to-mid frequency active sonar in this environment.     

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.     

一种分层海底反向散射模型

In order to predict the bottom backscattering strength more accurately, the stratified structure of the seafloor is considered. The seafloor is viewed as an elastic half-space basement covered by a fluid sediment layer with finite thickness. On the basis of calculating acoustic field in the water, the sediment layer, and the basement, four kinds of scattering mechanisms are taken into account, including roughness scattering from the water-sediment interface, volume scattering from the sediment layer, roughness scattering from the sediment-basement interface,and volume scattering from the basement. Then a backscattering model for a stratified seafloor applying to low frequency(0.1–10 kHz) is established. The simulation results show that the roughness scattering from the sediment-basement interface and the volume scattering from the basement are more prominent at relative low frequency(below 1.0 kHz). While with the increase of the frequency, the contribution of them to total bottom scattering gradually becomes weak. And the results ultimately approach to the predictions of the high-frequency(10–100 kHz) bottom scattering model. When the sound speed and attenuation of the shear wave in the basement gradually decrease, the prediction of the model tends to that of the full fluid model, which validates the backscattering model for the stratified seafloor in another aspect.     

Overview of SAX99: environmental considerations

A 1-km² area located 2 km off the Florida Panhandle (30°22.6'N; 86°38.7'W) was selected as the site to conduct high-frequency acoustic seafloor penetration, sediment propagation, and bottom scattering experiments. Side scan, multibeam, and normal incidence chirp acoustic surveys as well as subsequent video surveys, diver observations, and vibra coring, indicate a uniform distribution of surficial and subbottom seafloor characteristics within the area. The site, in 18-19 m of water, is characterized by 1-2-m-thick fine-to-medium clean sand and meets the logistic and scientific requirements specified for the acoustic experiments. This paper provides a preliminary summary of the meteorological, oceanographic, and seafloor conditions found during the experiments and describes the important physical and biological processes that control the spatial distribution and temporal changes in these characteristics     

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.     

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