A reflected energy prediction model for long-range hydroacoustic reflection in the oceans

Acoustic energy from underwater earthquakes and explosions can propagate over long distances with very little attenuation in the deep ocean. When this sound encounters a seamount, island, or continental margin, it can scatter and again propagate over long distances. Hydrophones in the deep sound channel can detect these reflections tens of minutes or hours after arrivals from the direct source-to-receiver path. This paper presents the Reflected Energy Prediction (REP) model, a model for predicting these reflected arrivals. For a given source and receiver, the REP model uses a detailed knowledge of the underwater environment and components of the Hydroacoustic Coverage Assessment Model, HydroCAM, to predict the impulse response of the ocean. When this impulse response is convolved with a source function, a waveform envelope prediction is made that can be compared with recorded data. In this paper we present the model and a few applications of the model using data recorded from earthquakes and explosions in the Atlantic and Indian Oceans. These examples illustrate the use of the model and initial steps toward model calibration.     

Locating animals from their sounds and tomography of the atmosphere: experimental demonstration.

Calling animals are located using widely distributed receivers, and the sounds from the animals are used to map the sound speed and wind fields by means of tomography. In particular, two Red-Winged Blackbirds Agelaius phoeniceus are correctly located within a meter using recordings from five receivers spread over a 20 by 30 m region. The demonstration hinges on two new developments. First, a new algorithm for blindly estimating the impulse response of the channel is shown capable of estimating the differences in the time of first arrivals at two receivers. Since it is known that the first arrivals travel along nearly straight paths, the difference in time constrains the animal's location to a hyperboloid, and the animal is located by intersecting hyperboloids from many pairs of receivers. Second, in order to accurately find the intersection point and map the sound speed and wind fields using tomography, a nonlinear equation is solved. The new algorithm for blindly estimating the impulse response of a channel offers a new way for locating sounds and making tomographic maps of the environment without any requirement for a model for the propagation of sound such as is needed for focalization and matched field processing.     

Modal analysis of broadband acoustic receptions at 3515-km range in the North Pacific using short-time Fourier techniques

In 1995-1996 the Acoustic Thermometry of Ocean Climate (ATOC) experiment provided an opportunity to study long-range broadband transmissions over a series of months using mode-resolving vertical arrays. A 75-Hz source off the California coast transmitted broadband pulses to receiving arrays in the North Pacific, located at ranges of 3515 and 5171 km. This paper develops a short-time Fourier transform (STFT) processor for estimating the signals propagating in the lowest modes of the ocean waveguide and applies it to analyze data from the ATOC experiment. The STFT provides a convenient framework for examining processing issues associated with broadband signals. In particular, this paper discusses the required frequency resolution for mode estimation, analyzes the broadband performance of two standard modal beamforming algorithms, and explores the time/frequency tradeoffs inherent in broadband mode processing. Short-time Fourier analysis of the ATOC receptions at 3515 km reveals a complicated arrival structure in modes 1-10. This structure is characterized by frequency-selective fading and a high degree of temporal variability. At this range the first ten modes have equal average powers, and the magnitude-squared coherence between the modes is effectively zero. The coherence times of the peaks in the STFT mode estimates are on the order of 5.5 min. An analysis of mean arrival times yields modal dispersion curves and indicates that there are statistically significant shifts in travel time over 5 months of ATOC transmissions.     

Comparison between ocean-acoustic fluctuations in parabolic-equation simulations and estimates from integral approximations

Line-integral approximations to the acoustic path integral have been used to estimate the magnitude of the fluctuations in an acoustic signal traveling through an ocean filled with internal waves. These approximations for the root-mean-square (rms) fluctuation and the bias of travel time, rms fluctuation in a vertical arrival angle, and the spreading of the acoustic pulse are compared here to estimates from simulations that use the parabolic equation (PE). PE propagations at 250 Hz with a maximum range of 1000 km were performed. The model environment consisted of one of two sound-speed profiles perturbed by internal waves conforming to the Garrett-Munk (GM) spectral model with strengths of 0.5, 1, and 2 times the GM reference energy level. Integral-approximation (IA) estimates of rms travel-time fluctuations were within statistical uncertainty at 1000 km for the SLICE89 profile, and in disagreement by between 20% and 60% for the Canonical profile. Bias estimates were accurate for the first few hundred kilometers of propagation, but became a strong function of time front ID beyond, with some agreeing with the PE results and others very much larger. The IA structure functions of travel time with depth are predicted to be quadratic with the form theta(2)vc0(-2)deltaz(2), where deltaz is vertical separation, c0 is a reference sound speed, and thetav is the rms fluctuation in an arrival angle. At 1000 km, the PE results were close to quadratic at small deltaz, with values of thetav in disagreement with those of the integral approximation by factors of order 2. Pulse spreads in the PE results were much smaller than predicted by the IA estimates. Results imply that acoustic tomography of internal waves at ranges up to 1000 km can use the IA estimate of travel-time variance with reasonable reliability.     

一种基于序贯估计的直达声区水面舰船 被动测距方法 *

该文利用基于射线的盲解卷积方法，从直达声区的水面舰船噪声中提取出船和锚系于海底的垂直接收阵之间的时域信道响应，并利用直达波在不同阵元相对于参考阵元的到达时间差，通过序贯方法，利用射线模型和声速剖面信息，对水面舰船距接收阵的距离进行了估计。通过处理海深约为580 m的2016年美国圣巴巴拉海峡的实验数据，对1.6~3.5 km直达声区范围内Anna Maersk商船与垂直阵之间的距离进行了估计，验证了测距方法的有效性，并将结果与系统测量值和几何方法的估计值进行了比较。由于该方法不需要对海底参数进行估计，所以在海底参数未知时要优于传统匹配场方法；在声速剖面存在跃层且海底为多层分布的复杂信道条件下，该方法仍能对距离进行有效估计，且与测量值的相对误差在6%以内，小于几何方法的估计误差，测距结果精度较高。     

Mapping of the ocean surface wind by ocean acoustic interferometers

Measurements of marine surface winds are crucial to understanding mechanical and thermodynamic forces on the ocean. Satellite measurements of surface winds provide global coverage but are problematic at high wind speeds. Acoustic techniques of wind speed retrieval, and even for tracking hurricanes, have been suggested as an alternative since wind is a strong source of ambient noise in the ocean. Such approaches involve near-local measurements with bottom-mounted hydrophones located close to the area of interest. This paper suggests a complementary approach: measuring directivity of low-frequency ambient noise in the horizontal plane. These measurements would employ long vertical line arrays (VLAs) spanning a significant portion of the ocean waveguide. Two VLAs separated by a distance of some tens of kilometers and coherently measuring acoustic pressure form a single ocean interferometer. By sampling the area of interest from different perspectives with at least two interferometers, marine surface winds might be mapped over horizontal scales of the order of 1000 km with about 10 km resolution (more specifically, the 10 km resolution here means that contribution from the basis functions representing surface wind field with the scale of spatial variations of the order of 10 km can be resolved; independent retrieval of the wind within 10(4) cells of a corresponding grid is hardly possible). An averaging time required to overcome statistical variability in the noise field is estimated to be about 3 h. Numerical simulations of propagation conditions typical for the North Atlantic Ocean are presented.     

Temporal coherence of acoustic signals in a fluctuating ocean

Temporal coherence of acoustic signals propagating in a fluctuating ocean is important for many practical applications and has been studied intensively experimentally. However, only a few theoretical formulations of temporal coherence exist. In the present paper, a three-dimensional (3D) modal theory of sound propagation in a fluctuating ocean is used to derive closed-form equations for the spatial-temporal coherence function of a broadband signal. The theory is applied to the analysis of the temporal coherence of a monochromatic signal propagating in an ocean perturbed by linear internal waves obeying the Garrett-Munk (G-M) spectral model. In particular, the temporal coherence function is calculated for propagation ranges up to 10(4) km and for five sound frequencies: 12, 25, 50, 75, and 100 Hz. Then, the dependence of the coherence time (i.e., the value of the time lag at which the temporal coherence decreases by a factor of e) on range and frequency is studied. The results obtained are compared with experimental data and predictions of the path-integral theory.     

A note on the proper frequency resolution for the room transfer function in the phased beam tracing method

In the simulation of acoustic response of a room, one cannot satisfactorily estimate the transfer function or the impulse response with a sparse frequency resolution. We have studied about the frequency resolution for a satisfactory estimation of room transfer function or the reverberation time using the phased beam tracing method. In order to investigate what happens with the frequency resolution, a number of frequency resolution (0.01-1 Hz), which are usually employed in the numerical analysis, was tested for a room model. It was found that the total accumulated phase in a transfer function and the late part of an impulse response of a room are influenced by the frequency resolution. Main reason for the difference is the correct detection of non-minimum phase zeros, depending on the frequency resolution. A criterion for a proper frequency resolution was suggested by considering the modal density of three-dimensional space. When a sparse dataset was initially given, we showed that cubic spline interpolation can be used to enhance the precision of detection for non-minimum phase zeros.     

深海表面层声速剖面变化对声传播的影响

在我国部分海域，表面声速常常在较短的时间内从弱正梯度变为弱负梯度，这将对表面附近的声传播产生巨大影响。本文通过数值模拟的方法，讨论了表面层声速梯度变化对声传播的影响，同时解释了某次海上实验出现的异常现象：在同一地点同一航向，前后两天海况相同，水文相差不大，第一天的实验可以在50km之外接收到信号，而第二天只能在30km处才能收到信号。     

Range dependence of the vertical structure of the sound field in the ocean

In the ocean without fluctuations, the sound field is calculated by the method of geometrical acoustics with allowance for purely water-path rays in a sound channel of canonical shape with a thickness of 4 km for distances of 500 and 2000 km. The sound field is determined as a sum of individual rays arriving at a given point with their own amplitudes and phases. It is shown that the vertical structure of the sound field consists of a number of caustics separated by regions with a quasi-random distribution of the field whose amplitude is much smaller than that in the caustics. At a fixed distance, the number of caustics is equal to the difference between the numbers of the ray turning points at the boundaries of the departure angle range. As the distance from the source increases, the number of caustics increases proportionally to distance.     

The use of in-situ test method EN 1793-6 for measuring the airborne sound insulation of noise barriers

The in situ measurement of the airborne sound insulation, as outlined in EN 1793-6:2012, is becoming a common means of quantifying the performance of road traffic noise reducing devices. Newly installed products can be tested to reveal any construction defects and periodic testing can help to identify long term weaknesses in a design. The method permits measurements to be conducted in the presence of background noise from traffic, through the use of impulse response measurement techniques, and is sensitive to sound leakage. Factors influencing the measured airborne sound insulation are discussed, with reference to measurements conducted on a range of traffic noise barriers located around Auckland, New Zealand. These include the influence of sound leakage in the form of hidden defects and visible air gaps, signal-to-noise ratio, and noise barrier height. The measurement results are found to be influenced by the presence of hidden defects and small air gaps, with larger air gaps making the choice of measurement position critical. A signal-to-noise ratio calculation method is proposed, and is used to show how the calculated airborne sound insulation varies with signal-to-noise ratio. It is shown that the measurement results are influenced by barrier height, through the need for reduced length Adrienne temporal windows to remove the diffraction components, prohibiting the direct comparison of results from noise barriers with differing heights.     

On the annoyance caused by impulse sounds produced by small, medium-large, and large firearms

A laboratory study was designed in which the annoyance was investigated for 14 different impulse sound types produced by various firearms ranging in caliber from 7.62 to 155 mm. Sixteen subjects rated the annoyance for the simulated conditions of (1) being outdoors, and (2) being indoors with the windows closed. In the latter case, a representative outdoor-to-indoor reduction in sound level was applied. It was anticipated that the presumed additional annoyance caused by the "heaviness" of the impulse sounds might be predicted from the difference between the C-weighted sound exposure level (CSEL; LCE) and the A-weighted sound exposure level (ASEL; LAE). In the outdoor rating conditions, the annoyance was almost entirely determined by ASEL. The explained variance, r2, in the mean ratings by ASEL was 0.95. In the indoor rating conditions, however, the explained variance in the annoyance ratings by (outdoor) ASEL was significantly increased from r2 = 0.87 to r2= 0.97 by adding the product (LCE-LAE)(LAE-alpha) as a second variable. In combination with a 12-dB adjustment for small firearms, the present results showed that for the entire set of impulse sounds rated indoors with windows closed, the rating sound level, Lr, is given by Lr=LAE +12dB+beta(LCE-LAE)(LAE-alpha), with alpha=45dB and beta=0.015dB(-1). For the outdoor rating condition, the optimal parameter values were equal to alpha=57 dB and, again, beta=0.015 dB(-1). In validation studies, in which the effects of the present rating procedure will be compared to field data, it has to be determined to what extent the constants alpha and beta have to be adjusted.     

Neural time and movement time in choice of whistle or pulse burst responses to different auditory stimuli by dolphins

Echolocating dolphins emit trains of clicks and receive echoes from ocean targets. They often emit each successive ranging click about 20 ms after arrival of the target echo. In echolocation, decisions must be made about the target--fish or fowl, predator or food. In the first test of dolphin auditory decision speed, three bottlenose dolphins (Tursiops truncatus) chose whistle or pulse burst responses to different auditory stimuli randomly presented without warning in rapid succession under computer control. The animals were trained to hold pressure catheters in the nasal cavity so that pressure increases required for sound production could be used to split response time (RT) into neural time and movement time. Mean RT in the youngest and fastest dolphin ranged from 175 to 213 ms when responding to tones and from 213 to 275 ms responding to pulse trains. The fastest neural times and movement times were around 60 ms. The results suggest that echolocating dolphins tune to a rhythm so that succeeding pulses in a train are produced about 20 ms over target round-trip travel time. The dolphin nervous system has evolved for rapid processing of acoustic stimuli to accommodate for the more rapid sound speed in water compared to air.     

Conditioning-induced protection from impulse noise in female and male chinchillas

Sound conditioning (pre-exposure to a moderate-level acoustic stimulus) can induce resistance to hearing loss from a subsequent traumatic exposure. Most sound conditioning experiments have utilized long-duration tones and noise at levels below 110 dB SPL as traumatic stimuli. It is important to know if sound conditioning can also provide protection from brief, high-level stimuli such as impulses produced by gunfire, and whether there are differences between females and males in the response of the ear to noise. In the present study, chinchillas were exposed to 95 dB SPL octave band noise centered at 0.5 kHz for 6 h/day for 5 days. After 5 days of recovery, they were exposed to simulated M16 rifle fire at a level of 150 dB peak SPL. Animals that were sound conditioned showed less hearing loss and smaller hair cell lesions than controls. Females showed significantly less hearing loss than males at low frequencies, but more hearing loss at 16 kHz. Cochleograms showed slightly less hair cell loss in females than in males. The results show that significant protection from impulse noise can be achieved with a 5-day conditioning regimen, and that there are consistent differences between female and male chinchillas in the response of the cochlea to impulse noise.     

海洋信道随机性引起的船舶噪声线谱衰减

阐述准周期性随机声脉冲序列作为船舶噪声模型的普适性,侧重讨论船舶噪声准周期性随机声脉冲序列所产生的线谱经随机海洋信道传输后的变化。指出在大多数工况下,船舶噪声并不是以几种不同频率的加性正弦波直接叠加在宽带连续谱噪声上,而是一类本身并不含有任何有限强度加性正弦波分量的信号。研究结果表明:海洋信道多途效应引起船舶噪声脉冲序列信号时域上时延扩展;海洋信道的随机性使船舶噪声随机声脉冲序列信号准周期性下降,导致船舶噪声线谱相对幅度“额外”快速衰减,对于存在大范围不均匀水团的海洋信道,50 km传播距离上其介质随机起伏引起的船舶噪声线谱“额外”衰减最高可达11 dB。海洋信道随机性引起的船舶噪声线谱衰减不可忽视。      

Model-oriented ocean tomography using higher frequency, bottom-mounted hydrophones

A tomographic scheme is presented that ingests ocean acoustic measurements into an ocean model using data from bottom-mounted hydrophones. The short distances between source-receiver pairs (1-10 km) means arrival times at frequencies of 8-11 kHz are readily detectable and often distinguishable. The influence of ocean surface motion causes considerable variability in acoustic travel times. Techniques are presented for measuring travel times and removing the variability due to surface waves. An assimilation technique is investigated that uses differences in measured and modeled acoustic travel times to impose corrections on the oceanographic model. Equations relating travel time differences to oceanographic variables are derived, and techniques are presented for estimating the acoustic and ocean model error covariance matrices. One test case using a single source-receiver pair shows that the tomographic information can have an impact on constraining the solution of the ocean circulation model but can also introduce biases in the predictions. A second test case utilizes knowledge of a bias in a model-predicted variable to limit grid cells that are impacted by the tomographic data. In this case, using the tomographic data results in significant improvements in the model predictions without introducing any biases.     

Blast noise characteristics as a function of distance for temperate and desert climates

Variability in received sound levels were investigated at distances ranging from 4 m to 16 km from a typical blast source in two locations with different climates and terrain. Four experiments were conducted, two in a temperate climate with a hilly terrain and two in a desert climate with a flat terrain, under a variety of meteorological conditions. Sound levels were recorded in three different directions around the source during the summer and winter seasons in each location. Testing occurred over the course of several days for each experiment during all 24 h of the day, and meteorological data were gathered throughout each experiment. The peak levels (L(Pk)), C-weighted sound exposure levels (CSEL), and spectral characteristics of the received sound pressure levels were analyzed. The results show high variability in L(Pk) and CSEL at distances beyond 2 km from the source for each experiment, which was not clearly explained by the time of day the blasts occurred. Also, as expected, higher frequency energy is attenuated more drastically than the lower frequency energy as the distance from the source increases. These data serve as a reference for long-distance blast sound propagation.     

Closed form solutions for the acoustical impulse response over a masslike or an absorbing plane

The transient sound field caused by a Dirac delta impulse function above an infinite locally reacting plane can be calculated by applying the inverse Fourier transform of the corresponding half-space Green's function in frequency domain. As a starting point, the representation given by Ochmann [J. Acoust. Soc. Am. 116(6), 3304-3311 (2004)] is used, which consists of discrete and continuous superposition of point sources. For a locally reacting plane with masslike character and also with pure absorbing behavior, it is possible to express the resulting impulse response in closed form. Such a result is surprising, because corresponding formulations in the frequency domain are not available yet. Hence, the first main result is the closed form solution Eq. (22) for an impulse response over an infinite plane with a pure imaginary impedance. The second main result is the closed form solution Eq. (53) for an impulse response over an infinite plane with a pure real impedance. As a particular application of both main results, a convolution technique is used for deriving formulas for point sources with a general time dependency. For special signals like an exponentially decaying time signal or a triangular shaped impulse, the resulting sound field can be presented in terms of elementary functions.     

Internal-wave time evolution effect on ocean acoustic rays

A range-dependent field of sound speed in the ocean, c(x,z), caused by internal waves, can give rise to instabilities in acoustic ray paths. Past work has shown the importance of the background, range-independent, sound-speed profile; the ray initial conditions; the source-receiver geometry (depths and range); and the strength of the internal waves. However, in the past the time evolution of the internal waves has been ignored on the grounds that the speed of internal waves is much slower than the speed of the acoustic wave. It is shown here by numerical simulation that two rays with identical initial conditions, traveling through an ocean with the same background profile and the same random realization of internal waves, but with the internal waves frozen in one case and evolving in the other, travel significantly different trajectories. The dependence of this "frozen-unfrozen" difference on the initial ray launch angle, the background profile, and the strength of the internal-wave spectrum, is investigated. The launch-angle difference that generates similar arrival-depth differences to those induced by internal-wave time evolution is on the order of 100 microrad. The pattern of differences is measured here by the arrival depth at the final range of 1000 km. The observed pattern as a function of launch angle, change in the background profile, and change in internal-wave strength is found to be nearly the same for "frozen-unfrozen" change as for a slight change in launch angle.