Two-way time spreading and path loss in shallow water at 20-40 kHz

Two-way time spreading and path-loss measurements were collected in water 100 m deep, off the coast of Nova Scotia. Data were collected at frequencies of 20-22 kHz, 27-29 kHz, and 35-37 kHz using linear FM pulses 0.160 s in duration. The source-receiver was an anchored, high-frequency active sonar, and the target was a free-drifting echo repeater. Sonar and target positions were recorded using a portable tracking range. In the paper, two-way time spreading and path loss measurements are compared with modeled estimates obtained using an enhanced version of the generic sonar model (GSM). The GSM estimates of time spreading due to multipath propagation compare favorably with the experimental data. The model indicates that the path loss for individual eigenrays was extremely sensitive to fluctuations in the sound-speed profile. This led to substantial variation in the model output depending on the choice of profile. In place of the model, an empirical estimate of path loss was computed from the data. We obtained a two-way spreading loss of 2[18.4log₁₀(R)] where R is the range from sonar to target. The data were also used to compute the standard deviation of the received echo intensity at each frequency. The standard deviation was computed two different ways. First it was computed using the peak echo level from each of the pulses at a given frequency. Then, it was computed from the total energy received from each of the pings. At all frequencies, the standard deviation was 1-2 dB lower when computed from the total received energy     

A model for numerical simulation of nonstationary sonar reverberation using linear spectral prediction

An innovative approach to the numerical generation of nonstationery reverberation time series is presented and demonstrated. The computer simulated reverberation time series are of high quality, in that they are accurate representations of those which would result from an actual sonar system (transmit/receive and horizontal/ vertical beampatterns; pulse type, shape, length, and power; frequency and sampling rate), platform (speed and depth), and environment (wind speed and direction, backscattering strengths, and propagation loss). Volume, surface, and/or bottom reverberation as seen by a multiple beam sonar on a moving platform is generated. The approach utilizes recent developments in linear spectral prediction research in which the spectra of stochastic processes are modeled as rational functions and algorithms are used to efficiently compute optimal estimates of coefficients which specify the spectra. A two-fold sequence is formulated; first, the expected reverberation spectra for all beams are predicted and, second, the stochastic time series are generated from the expected spectra. The expected spectra are predicted using a numerical implementation, referred to as the REVSPEC (reverberation spectrum) model, of a general formulation of Faure, Ol'shevskii, and Middleton. Given the spectra, the Levinson-Durbin method is used to solve the Yule-Walker equations of the autoregressive formulation of linear spectral prediction. The numerical implementation of the approach, referred to as the REVSIM (reverberation simulation) model, produces nonstationary coherent multiple-beam reverberation time series. The formulation of the REVSIM model is presented and typical results given. A comparison is made between the simulation outputs of the REVSIM model and those of the REVGEN (reverberation generator) model, a standard well-accepted time series simulation model, to demonstrate the validity of the new approach.     

Do equilibrium conditions exist during sediment transport studies on continental margins? An example from the East China Sea

Analysis of high-resolution (<1 m) GeoPulse seismic reflection and Chirp sonar data from the continental margin of the East China Sea reveals the presence of a thin (1- to 5-m-thick) sediment veneer over a highly reflective transgressive surface. Due to reverberation of the output signal and the effect of the bubble pulse, this veneer may erroneously be interpreted as a laterally continuous feature. Sonar data, which does not suffer from these ghosting effects, reveals that the sea bottom is often an erosional surface that exposes older, relict sediments, which were deposited under a radically different hydrodynamic regime than exists in these locations today.     

Theoretical accuracy of Doppler navigation sonars and acoustic Doppler current profilers

This paper addresses the problem of Doppler shift estimation in Doppler sonar systems. The analysis focuses on the single-beam geometry formed by a circular planar array and considers both narrow-band (or so-called incoherent) and wide-band (or coherent) Doppler sonars, transmitting, respectively, one long continuous-wave pulse and a train of short continuous-wave pulses. The correlation function of the reverberation signal at the beam output is derived for volume reverberation. Directive transmission or reception and a combination of both is considered. Estimation theory is applied to derive the Cramer-Rao bound of the Doppler parameter estimate. The effect of pulse duration, sonar geometry, beamwidth, and signal-to-noise ratio are discussed. The accuracy of coherent and incoherent systems is compared for a specific case.     

Statistical characterization of active sonar reverberation using extreme value theory

The statistics of reverberation in active sonar are characterized by non-Rayleigh distributed amplitudes in the normalized matched filter output. Unaccounted for, this property can lead to high false-alarm rates in fixed-threshold detectors. A new approach to modeling threshold-crossing statistics based on extreme value theory is proposed, which uses the generalized Pareto distribution as the unique asymptotic model of the tail distribution, valid at large thresholds. Methods of parameter estimation are discussed and applied to active sonar reverberation collected on a hull-mounted sonar system. The statistics of reverberation in active sonar are found to generally have a power-law behavior in the tails with a shape parameter that is persistent in time and bandwidth dependent. The threshold needed for accurate parameter estimation is generally found to be well below that of typical fixed-threshold detectors.     

Simulation of non-Rayleigh reverberation and clutter

The simulation of active sonar reverberation time series has traditionally been done using either a computationally intensive point-scatterer model or a Rayleigh-distributed reverberation-envelope model with a time-varying power level. Although adequate in scenarios where reverberation arises from a multitude of scatterers, the Rayleigh model is not representative of the target-like non-Rayleigh reverberation or clutter commonly observed with modern high-resolution sonar systems operating in shallow-water environments. In this paper, techniques for simulating non-Rayleigh reverberation are developed within the context of the finite-number-of-scatterers representation of K-distributed reverberation, which allows control of the reverberation-envelope statistics as a function of system (beamwidth and bandwidth) and environmental (scatterer density and size) parameters. To avoid the high computational effort of the point-scatterer model, reverberation is simulated at the output of the matched filter and is generated using efficient approximate methods for forming K-distributed random variables. Finite impulse response filters are used to introduce the effects of multipath propagation and the shape of the reverberation power spectrum, the latter of which requires the development of a prewarping of the K distribution parameters to control the reverberation-envelope statistics. The simulation methods presented in this paper will be useful in the testing and evaluation of active sonar signal processing algorithms, as well as for simulation-based research on the effects of the sonar system and environment on the reverberation-envelope probability density function.     

Shallow-water reverberation level: measurement technique and initial reference values

A quality database of reverberation is absolutely essential if one is to understand the shallow-water reverberation problem. However, to get wideband reverberation levels (RL) simultaneously for both short and long ranges at low- and mid-frequencies is a delicate task that can be subject to errors. This paper introduces a simple method to get RL for the Asian Sea International Acoustics Experiment in the East China Sea (ASIAEX01). Special attention is paid to the measurements of the RL at short- and mid-ranges. With this method, one does not need to accurately calibrate hydrophones and measurement systems, or to measure absolute source level (SL). It can avoid signal overflow and saturation problems caused by powerful sound sources. The RL (relative to SL) at 1 s (or at 2 s) after an explosive source is detonated is defined as the initial reference reverberation level (IRRL). The IRRLs from four sites with different sandy sediments and different water depths have been given as a function of frequency in the 150-2500 Hz range. A mathematical model gives a physical explanation of the measured IRRL data. The resultant RL and IRRL may offer some reference values for the design of reverberation measurements or numerical simulations of shallow-water reverberation and bottom scattering.     

Data-model comparisons of reverberation at three shallow-watersites

Reverberation measurements made by the SACLANT Undersea Research Centre at three shallow-water sites (130-190-m depth) are compared with each other and with estimates from the DREA normal-mode reverberation model OGOPOGO. The experiments over silt-clay and sand seabeds were conducted at slightly bistatic geometries (0.7-6.0-km source-receiver separation), using explosive sources detonated at mid-water depths. The signals were received on hydrophones of either a vertical or horizontal array and analyzed in one-tenth-decade frequency bands from 25 to 1000 Hz. The data are compared with each other to investigate the site differences and frequency dependencies, and with the estimates from the reverberation model OGOPOGO to interpret the data and to obtain a qualitative measure of the scattering. For modeling purposes, geoacoustic models of the seabed were assumed, and the reverberation data were fitted by adjusting the Lambert bottom scattering coefficients. Good model agreement was obtained with both individual hydrophone and data. Though somewhat sensitive to the geoacoustic the Lambert coefficients give a measure of the frequency dependence of the scattering. For the silt-clay bottom, the scattering is weak but is independent of frequency; for the sand bottoms, the scattering is stronger and increases with frequency. These results are compared with estimates from other experiments     

On shallow-water bottom reverberation frequency dependence

Determinations of bottom scattering strength in the decade below 1 kHz under downward refracting conditions have been made using acoustic reverberation and transmission data from the 2001 East China Sea Asian Seas International Acoustic Experiment (ASIAEX). The measurements were performed using explosive sources and receiving hydrophones in ship-suspended vertical-line arrays. The focus of this paper has been the dependence of bottom scattering strength on the frequency and characterization of the uncertainties associated with the extraction of scattering strength from reverberation. The derived bottom scattering strength gradually rises with frequency from 100-300 Hz and then more rapidly above 300 Hz. A potential explanation suggests that the frequency variation results from two scattering mechanisms, rough layer scattering at the low end of the band and sediment near-surface volume scattering at the high end. The spatial extrapolation of these results is explored by comparing them with similarly derived scattering strengths using data obtained under the Navy's Harsh Environments Program at a somewhat separated site (56 km) under environmental conditions similar to those during ASIAEX. In the ASIAEX analysis, it has been found that the largest source of uncertainty in the scattering-strength frequency dependence arises from persistence of finite-amplitude effects associated with the source signal.     

Estimation of Transmission Loss in the Presence of Geoacoustic Inversion Uncertainty

A common problem in sonar system prediction is that the ocean environment is not well known. Utilizing probabilistic based results from geoacoustic inversions we characterize parameters relevant to sonar performance. This paper describes the estimation of transmission loss and its statistical properties based on posterior parameter probabilities obtained from inversion of ocean acoustic array data. This problem is solved by first finding an ensemble of relevant environmental model parameters and the associated posterior probability using a likelihood based inversion of the acoustic array data. In a second step, each realization of these model parameters is weighted with their posterior probability to map into the transmission loss domain. This approach is illustrated using vertical-array data from a recent benchmark data set and from data acquired during the Asian Seas International Acoustics Experiment (ASIAEX) 2001 in the East China Sea. The environmental parameters are first estimated using a probabilistic-based geoacoustic inversion technique. Based on the posterior probability that each of these environmental models fits the ocean acoustic array data, each model is mapped into transmission loss. This enables us to compute a full probability distribution for the transmission loss at selected frequencies, ranges, and depths, which potentially could be used for sonar performance prediction.     

"Principal component inverse" algorithm for detection in thepresence of reverberation

Detection in the presence of reverberation is often difficult in active sonar, due to the reflection/diffusion/diffraction of the transmitted signal by the ocean surface, ground, and volume. A modelization of reverberation is often used to improve detection because classical algorithms are inefficient. A commonly used reverberation model is colored and nonstationary noise. This model leads to elaborate detection algorithms which normalize and whiten reverberation. In this paper, we focus on a more deterministic model which considers reverberation as a sum of echoes issued from the transmitted signal. The Principal Component Inverse (PCI) algorithm is used with this model to estimate and delete the reverberation echoes. A rank analysis of the observation matrix shows that PCI is efficient in this configuration under some conditions, such as when the transmitted signal is Frequency Modulated. Both methods are validated with real sonar surface reverberation noise. We show that whitening has poor performance when reverberation and target echo have the same properties, while PCI maintains the same performance whatever the reverberation characteristics. Further, we extend the algorithms to spatio-temporal data. We propose a new algorithm for PCI which allows better echo separation. This new method is shown to be more efficient on real spatio-temporal data     

A statistical study of acoustic signals backscattered from the seabottom

Backscattering measurements were performed in shallow water on sand, gravel, and clay bottoms. The equipment included a parametric array that emitted pulses of differential frequencies (8 to 40 kHz) with a 3° directivity. The ranges did not exceed 50 m. The grazing angles varied from 4° to 90°. The bottom backscattering strength does not depend on the emitted pulse type (frequency and length). If one fits a Lambert law to the variations of the backscattering strength versus the grazing angle, the value at the origin fluctuates between-15 and-22 dB without any clear effect from the different bottom types. Statistical tests show that under the experimental measurement conditions: (1) the alternative received signal does not generally follow a normal distribution; (2) among five classical distributions in sonar and radar that have been fitted to the detected-integrated signal (exponential, Weibull, chi-2, log-normal, Rice), the best-fitted law is the log-normal; (3) signals backscattered by separated areas of the same bottom can hardly be regarded as stationary and, even less, homogeneous; and (4) with an anisotropic bottom topography the statistical properties depend on the aspect under which this topography is seen     

Shallow water reverberation: normal-mode model predictions comparedwith bistatic towed-array measurements

A normal-mode model for calculating reverberation in shallow water is presented. Some illustrative calculations are given for the bistatic case and for vertical and horizontal line-array receivers. Emphasis is on comparison with measurements of bistatic reverberation obtained at a shallow-water area in the Mediterranean. The data are from explosive sources received by a towed array, analyzed in one-tenth-decade frequency bands at subkilohertz frequencies. Model calculations for a flat-bottomed environment indicate a strong dependence on propagation conditions and a weak dependence on beam steering direction. Preliminary comparisons give quite good agreement between measured reverberation and model predictions, but point to the need for extending modeling efforts to handle range-dependent environments     

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

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

多波束声呐图像条带中央和边缘残差处理方法

多波束声呐图像是进行海底底质分类的主要数据源之一,由于受海洋噪声、声波散射和混响、仪器设备等因素影响,其经各项常规改正后仍存在明显残差,突出表现在中央波束区和条带重叠区,难以形成高质量的声呐图像。文中分析了多波束声呐图像残差的成因及影响,提出了一种基于多条带最小二乘拟合的多波束声呐图像残差处理方法。首先,得到相邻声脉冲(ping)信号中央区域、重叠区域以及整体趋势的拟合函数;然后,通过拟合函数计算得到中央和重叠区域的残差改正系数;最后,通过改正系数进行残差改正。实验分析表明,该方法在保留原始细节的基础上,有效削弱了残差对声呐图像的影响,对多波束声呐图像处理具有参考和应用价值。     

Mine avoidance techniques for underwater vehicles

It is shown that by implementing certain mine avoidance techniques, an underwater vehicle equipped with an obstacle avoidance sonar (OAS) and a navigation system can safely navigate an unknown minefield. The mine avoidance techniques take into account the physical limitations of the sonar and the navigation system, the maneuverability constraints on the underwater vehicle, and the required safe standoff distance from all mines. Extensive computer simulations have verified the mine avoidance capability in more than 50 different minefields. In all 50 simulations the vehicle reached a predetermined end point and maintained at least the specified, minimum safe standoff distance from each mine. The simulation accurately models the major difficulties associated with the sonar, the navigation system, and the vehicle dynamics. The sonar model includes surface, bottom, and volume reverberation; thermal, ambient, and flow noises; actual receiver and projector beam patterns; and false alarms and missed detections. The navigation system model contains the effects of biases, random noises, and scale factor errors. The vehicle dynamic model simulates angular velocities and accelerations associated with underwater vehicles     

Analysis of shallow-water experimental acoustic data including acomparison with a broad-band normal-mode-propagation model

Channel temporal variability, resulting from fluctuations in oceanographic parameters, is an important issue for reliable communications in shallow-water-long-range acoustic propagation. As part of an acoustic model validation exercise, audio-band acoustic data and oceanographic data were collected from shallow waters off the West Coast of Scotland. These data have been analyzed for temporal effects. The average impulse response for this channel has been compared with simulations using a fast broad-band normal-mode propagation model. In this paper, we also introduce a novel technique for estimating and removing the bistatic reverberation contribution from the data. As propagation models do not necessarily account for reverberation, it has to be extracted from the signals when comparing measured and modeled transmission loss     

General motion estimation from correlation sonar

The displacement of a sonar array can be estimated accurately using the correlation of bottom reverberation signals, at successive sweeps. This is the principle of the pulsed correlation log. In this paper, the conditions and accuracy of array horizontal translation estimation are analyzed. The case of a three-dimensional array is considered first, to show that in this case arbitrary translation and rotation can be estimated. The case of a plane array is then analyzed and it is shown that such an array allows estimation of horizontal translation. The derivation relies on modeling the space-time correlation function of bottom reverberation, which is assumed isotropic. Both directive and omnidirectional transmissions are considered. Accuracy of displacement estimates are derived, showing the influence of wavelength, grazing angle, bandwidth, number of overlapping hydrophones, and reverberation-to-noise ratio     

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.     

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.