Daniel Rouseff

Underwater acoustic communication by passive-phase conjugation: theory and experimental results

A new method for coherent underwater acoustic communication called passive phase conjugation is evaluated. The method is so named because of conceptual similarities to active phase conjugation methods that have been demonstrated in the ocean. In contrast to active techniques, however, the array in passive phase conjugation needs only receive. The procedure begins with a source transmitting a single probe pulse. After waiting for the multipathed arrivals to clear, the source then transmits the data stream. At each element in the distant receiving array, the received probe is cross-correlated with the received data stream. This cross-correlation is done in parallel at each array element and the results are summed across the array to achieve the final communication signal suitable for demodulation. As the ocean changes, it becomes necessary to break up the data stream and insert new probe pulses. Results from an experiment conducted in Puget Sound near Seattle are reported. Measurements were made at multiple ranges and water depths in range-dependent environments.

Multichannel equalization by decision-directed passive phase conjugation: experimental results

An adaptive technique for underwater acoustic communication using passive phase conjugation (PPC) is developed. Multipath channel-parameter identification is accomplished by decision-directed model building and finite-window block-updated least squares computed by LSQR (an iterative linear systems solver). The resulting channel estimates are then used by the PPC processor to generate decisions for use in the next processing block. This architecture effectively accomplishes array equalization with low computation cost in shallow-water environments that exhibit rapidly fluctuating multipath scattering. The performance on shallow-water acoustic communications channels is demonstrated at ranges of 0.9-4.6 km under windy surface conditions and shipping noise, using measured wide-band telemetry data with binary phase-shift keying signaling. The algorithm is evaluated with sparse receiver apertures using subsets of a 14-element array.

Detection of high-intensity focused ultrasound liver lesions using dynamic elastometry.

A novel ultrasound technique was developed for detecting the distribution of stiffness in biological tissue. The method, which we call ‘dynamic elastometry,’ involves applying a low-frequency vibration (≤5 Hz) to the tissue and measuring the resulting velocity pattern within the sample using Doppler spectral analysis. Based upon the velocity differences, an elastically stiff region can be differentiated from surrounding soft tissue. Dynamic elastometry was used to both detect and quantify lesions produced by high-intensity focused ultrasound (HIFU) in porcine livers. Measurements of the lesion position and length agreed well with independent geometric measurements. The mean and standard deviation of the differences between the two types of measurement were −0.01cm and 0.10 cm for lesion position, and −0.05 cm and 0.12 cm for lesion length, respectively. The relative stiffness between lesions and normal liver tissue was estimated by the velocity gradient ratio. Results were compared with the Youngs modulus ratios between lesion and normal liver tissue obtained from mechanical measurement. The dynamic elastometric estimates had a strong linear correlation with the mechanical measurements (r = 0.93) but were smaller than the latter by 26%.

Effect of Shallow Water Internal Waves on Ocean Acoustic Striation Patterns

Abstract Contour plots of underwater acoustic intensity, mapped in range and frequency, often exhibit striations. It has been claimed that a scalar parameter ‘beta’, defined in terms of the slope of the striations, is invariant to the details of the acoustic waveguide. In shallow water, the canonical value is β=1. In the present paper, the waveguide invariant is modelled as a distribution rather than a scalar. The effects of shallow water internal waves on the distribution are studied by numerical simulation. Realizations of time-evolving shallow water internal wave fields are synthesized and acoustic propagation simulated using the parabolic equation method. The waveguide invariant distribution is tracked as the internal wave field evolves in time. Both random background internal waves and more event-like solitary internal waves are considered.

Modeling the Waveguide Invariant as a Distribution

The “invariant parameter” called “beta” is often useful for describing the acoustic interference pattern in a waveguide. For some shallow water waveguides, the measured acoustic intensity might contain contributions from several propagating acoustic modes. For each pair of these modes, a different value for the waveguide invariant might apply. If the acoustic intensity is measured over some distributed aperture and finite bandwidth, it may become difficult to assign a single value to beta. In the present work, the waveguide invariant is treated as a distribution. An algorithm for estimating this distribution for a general measurement geometry is developed. The algorithm is exercised for different classes of shallow water waveguides. When the propagation is dominated by modes interacting with the sea surface, the distribution can be sharply peaked. For cases where the sound speed profile creates a duct, the distribution is more diffuse. The effects of source/receiver depth, range, bandwidth and bottom atte...

Coherence of acoustic modes propagating through shallow water internal waves

The 1995 Shallow Water Acoustics in a Random Medium (SWARM) experiment [Apel et al., IEEE J. Ocean. Eng. 22, 445-464 (1997)] was conducted off the New Jersey coast. The experiment featured two well-populated vertical receiving arrays, which permitted the measured acoustic field to be decomposed into its normal modes. The decomposition was repeated for successive transmissions allowing the amplitude of each mode to be tracked. The modal amplitudes were observed to decorrelate with time scales on the order of 100 s [Headrick et al., J. Acoust. Soc. Am. 107(1), 201-220 (2000)]. In the present work, a theoretical model is proposed to explain the observed decorrelation. Packets of intense internal waves are modeled as coherent structures moving along the acoustic propagation path without changing shape. The packets cause mode coupling and their motion results in a changing acoustic interference pattern. The model is consistent with the rapid decorrelation observed in SWARM. The model also predicts the observed partial recorrelation of the field at longer time scales. The model is first tested in simple continuous-wave simulations using canonical representations for the internal waves. More detailed time-domain simulations are presented mimicking the situation in SWARM. Modeling results are compared to experimental data.

Tomographic reconstruction of stratified fluid flow

A method for imaging a moving fluid is proposed and evaluated by numerical simulation. A cross-section of a three-dimensional fluid is probed by high-frequency acoustic waves from several different directions. Assuming straight-ray geometric acoustics, the time of flight depends on both the scaler sound speed and the vector fluid velocity. By appropriately combining travel times, projections of both the sound speed and the velocity are isolated. The sound speed is reconstructed using the standard filtered backprojection algorithm. Though complete inversion of velocity is not possible, sufficient information is available to recover the component of fluid vorticity transverse to the plane of insonification. A new filtered backprojection algorithm for vorticity is developed and implemented. To demonstrate the inversion procedure, a 3-D stratified fluid is simulated and travel time data are calculated by path integration. These data are then inverted to recover both the scaler sound speed and the vorticity of the evolving flow. >

Blind deconvolution for robust signal estimation and approximate source localization.

Synthetic time reversal (STR) is a technique for blind deconvolution in an unknown multipath environment that relies on generic features (rays or modes) of multipath sound propagation. This paper describes how ray-based STR signal estimates may be improved and how ray-based STR sound-channel impulse-response estimates may be exploited for approximate source localization in underwater environments. Findings are based on simulations and underwater experiments involving source-array ranges from 100 m to 1 km in 60 -m-deep water and chirp signals with a bandwidth of 1.5-4.0 kHz. Signal estimation performance is quantified by the correlation coefficient between the source-broadcast and the STR-estimated signals for a variable number N of array elements, 2 ≤ N ≤ 32, and a range of signal-to-noise ratio (SNR), -5 dB ≤ SNR ≤ 30 dB. At high SNR, STR-estimated signals are found to have cross-correlation coefficients of ∼90% with as few as four array elements, and similar performance may be achieved at a SNR of nearly 0 dB with 32 array elements. When the broadband STR-estimated impulse response is used for source localization via a simple ray-based backpropagation scheme, the results are less ambiguous than those obtained from conventional broadband matched field processing.

Intersymbol interference in underwater acoustic communications using time-reversal signal processing.

Coherent underwater communication is hampered by the time spread inherent to acoustic propagation in the ocean. Because time-reversal signal processing produces pulse compression, communications has been suggested as a natural application of the technique. Passive versions of time-reversal processing use a receive-only array to do combined temporal and spatial matched filtering. It can be shown, however, that the pulse compression it achieves is not perfect and that an equalizer that relies solely on time-reversal processing will have an error floor caused by uncompensated intersymbol interference (ISI). In the present paper, a physics-based model is developed for the uncompensated ISI in a passive time-reversal equalizer. The model makes use of a normal-mode expansion for the acoustic field. The matched-filtering integral is approximated and the intermediate result interpreted using the waveguide invariant. After combining across the array and sampling, formal statistical averages of the soft demodulation output are calculated. The results show how performance scales with bandwidth, with the number and position of array elements, and with the length of the finite impulse response matched filters. Good agreement is obtained between the predicted scaling and that observed in field experiments.

Striation-based beamforming for estimating the waveguide invariant with passive sonar

The waveguide invariant summarizes the pattern of constructive and destructive interference between acoustic modes propagating in the ocean waveguide. For many sonar signal-processing schemes, it is essential to know the correct numerical value for the waveguide invariant. While conventional beamforming can estimate the ratio between the waveguide invariant and the range to the source, it cannot unambiguously separate the two terms. In the present work, striation-based beamforming is developed. It is shown that the striation-based beamformer can be used to produce an estimate for the waveguide invariant that is independent of the range. Simulation results are presented.