Single hydrophone source localization

The method presented in this paper assumes that the received signal is a linear combination of delayed and attenuated uncorrelated replicas of the source emitted waveform. The set of delays and attenuations, together with the channel environmental conditions, provide sufficient information for determining the source location. If the transmission channel is assumed known, the source location can be estimated by matching the data with the acoustic field predicted by the model conditioned on the estimated delay set. This paper presents alternative techniques that do not directly attempt to estimate time delays from the data but, instead, estimate the subspace spanned by the delayed source signal paths. Source Localization is then done using a family of measures of the distance between that subspace and the subspace spanned by the replicas provided by the model. Results obtained on the INTIMATE'96 data set, in a shallow-water acoustic channel off the coast of Portugal, show that a sound source emitting a 300-800-Hz LFM sweep could effectively be localized in range or depth over an entire day     

VLF source localization with a freely drifting acoustic sensorarray

The source localization and tracking capability of the freely drifting Swallow float volumetric array is demonstrated with the matched-field processing (MFP) technique using the 14-Hz CW data collected during a 1989 float experiment conducted in the northeast Pacific. Initial MFP of the experimental data revealed difficulties in estimating the source depth and range while the source azimuth estimate was quite successful. The main cause of the MFP performance degradation was incomplete knowledge of the environment. An environment adaptation technique using a global optimization algorithm was developed to alleviate the environmental mismatch problem, allowing the ocean-acoustic environment to be adapted to the acoustic data in a matched-field sense. Using the adapted environment, the 14-Hz source was successfully localized and tracked in azimuth and range within a region of interest using the MFP technique at a later time interval. Two types of environmental parameters were considered, i.e., sound speed and modal wave number. While both approaches yield similar results, the modal wave number adaptation implementation is more computationally efficient     

Matched-field processing for broad-band source localization

In the Hudson Canyon experiment, a sound source moved at a constant depth in 73 m of water while transmitting four tonals. The signal was received on a vertical array of hydrophones that spanned the water column. The data set from this experiment has become a standard test case for studying source tracking using matched field processing. As part of that process it was important to first determine a suitable environment model and demonstrate the feasibility of matched-field processing. In this paper, we provide the background on the original data processing that was done to accomplish this. Several interesting results emerged from that study. Frequency averaging was demonstrated to be extremely beneficial when used with the Bartlett processor. However, the popular Minimum Variance processor performed poorly. Finally we discuss a very simple approach to combining the energy coherently that provided significantly improved results.     

A study on propeller noise source localization in a cavitation tunnel

Localizing noise sources in cavitation experiments is an important research subject along with predicting noise levels. A cavitation tunnel propeller noise localization method is presented. Propeller noise measurement experiments were performed in the MOERI cavitation tunnel. To create cavitating conditions, a wake-generating dummy body was devised. In addition, 10 hydrophones were put inside a wing-shaped casing to minimize the unexpected flow inducing noise around the hydrophones. After measuring both of the noises of the rotating propeller behind the dummy body and acoustic signals transmitted by a virtual source, the data were processed via three objective functions based on the ideas of matched field processing and source strength estimation to localize noises on the propeller plane. In this paper, the measured noise analysis and the localization results are presented. Through the experiments and the analysis, it was found that the source localization methods that have been used in shallow water applications could be successfully adapted to the cavitation tunnel experiments.     

The spherical interpolation method of source localization

A closed-form least squares approximate maximum likelihood method for localization of broad-band emitters from time-difference-of-arrival (TDOA) measurements, called the spherical interpolation (SI) method, is presented. The localization formula is derived from least squares "equation-error" minimization. Computer simulation results show that the SI method has variance approaching the Cramer-Rap lower bound.     

A skewed eddy of Batchelor-modon type

A Batchelor-modon eddy is a highly specialized nonlinear vortex pair, whose potential vorticity depends linearly on the stream function viewed from the coordinates moving with the translation velocity of the eddy. To generalize it, a skewed model is developed by introducing a cubic nonlinearity in addition to the linear term. A perturbation analysis shows that the eddy region is no longer a circle but is elongated longitudinally or transversely according as the sign of the cubic term. Moreover, the eddy is slightly flattened or steepened. The cubic term increases or decreases the translation velocity, if the average radius and the amplitude are fixed. A numerical experiment on anf-plane is carried out to show that these skewed eddies retain their initial forms even after they turn a corner of the basin; they are as stable as (first-mode) standard Batchelor-modon eddies. The present skewed model gives a reasonable qualitative interpretation of deformed eddies which result from merging of two eddies or from initially Gaussian eddies near the boundary.     

Broad-band matched-field source localization in the east China Sea

An incoherent broad-band frequency (100-200 Hz) domain Bartlett processor is applied to the wide-band source (WBS) signals for source localization. The coupled normal mode-parabolic equation theory based on the WKBZ approximation is used to calculate the replica fields in the sloping bottom environment. The experimental analysis shows that the accuracy of the source localization is largely improved with the consideration of the slope of bottom. The range estimates of the majority of signals by localized by matched-field processing in the range from 30.0 to 50.0 km are consistent with the global positioning system measured ranges.     

基于密度峰值聚类的中尺度涡轨迹自动追踪方法

中尺度涡信息的提取包括涡旋的识别和轨迹追踪,其自动识别与追踪对于基于海量数据的中尺度涡分析十分重要。传统涡旋轨迹自动追踪方法一般需要预先设定搜索半径的阈值,存在一定的主观性。针对传统中尺度涡轨迹追踪方法存在的问题,论文从聚类的角度出发,提出基于密度峰值聚类算法实现对涡旋轨迹的自动追踪,并以南海中尺度涡追踪为例,将基于聚类的追踪算法与传统的相似度追踪算法进行比较分析。结果表明:（1）基于密度峰值聚类算法,可实现对海洋中尺度涡的自动追踪,该算法涡旋追踪准确率优于传统相似度算法;（2）该涡旋追踪算法对资料的完整性依赖度较低,特别是对于存在部分缺损数据的情况仍能较准确追踪;（3）该追踪算法克服了传统涡旋追踪算法需要预先设定搜索半径阈值的问题,自适应性更强。     

不确定海洋环境中基于贝叶斯理论的多声源定位算法

环境参数失配导致定位性能大幅度下降是匹配场定位所面临的难题之一。应用贝叶斯理论对环境聚焦,是当前解决该难题的研究热点。环境聚焦方法的实质是将未知环境参数和声源位置联合优化估计,当出现多个目标时,估计的参数会随着声源个数成倍增加,因此不得不利用有限的观测信息来实现众多参数的估计。本文采用最大似然比方法,获得信号源谱和误差项的最大似然估计,实现这些敏感性较弱参数的间接反演,有效降低了反演参数维数和定位算法复杂度。针对遗传算法的早熟和稳定性差的问题,改进了似然函数的经验表达式。将多维后验概率密度在参数起伏变化范围内积分,得到反演参数的一维边缘概率分布,求解最优值的同时进行反演结果的不确定性分析。本文仿真了位于相同距离、不同深度的两个声源,使用仿真实验验证了提出算法的有效性。     

Robust source localization in shallow water based on vector optimization

Owing to the multipath effect, the source localization in shallow water has been an area of active interest. However, most methods for source localization in shallow water are sensitive to the assumed model of the underwater environment and have poor robustness against the underwater channel uncertainty, which limit their further application in practical engineering. In this paper, a new method of source localization in shallow water, based on vector optimization concept, is described, which is highly robust against environmental factors affecting the localization, such as the channel depth, the bottom reflection coefficients, and so on. Through constructing the uncertainty set of the source vector errors and extracting the multi-path sound rays from the sea surface and bottom, the proposed method can accurately localize one or more sources in shallow water dominated by multipath propagation. It turns out that the natural formulation of our approach involves minimization of two quadratic functions subject to infinitely many nonconvex quadratic constraints. It shows that this problem (originally intractable) can be reformulated in a convex form as the so-called second-order cone program (SOCP) and solved efficiently by using the well-established interior point method, such as the software tool, SeDuMi. Computer simulations show better performance of the proposed method as compared with existing algorithms and establish a theoretical foundation for the practical engineering application.     

黑潮的涡分辨率数值模拟

本文利用涡分辨率的HYCOM模式,以NCEP月平均再分析资料(1979—1993)为驱动场,并采用单向数值嵌套的方式对黑潮流域进行数值模拟,成功模拟了黑潮流域的高分辨率流场特征。模拟结果显示:黑潮路径符合前人对黑潮的认识;在地形和流量的共同作用下,黑潮对吕宋海峡的入侵呈现多平衡态的特征;日本以南的黑潮路径发生多种时间尺度的摆动(从季节内到年际)。黑潮在PN断面上流速跟同期观测十分相符,流轴集中在陆架破折处,季节变化较弱。台湾岛以东黑潮,东海黑潮以及吐噶喇海峡黑潮的流量符合对应时期观测,并且各自呈现出不同的变化特点。     

A boundary-value problem for the parameterized mesoscale eddy transport

We present a physically and numerically motivated boundary-value problem for each vertical ocean column, whose solution yields a parameterized mesoscale eddy-induced transport streamfunction. The new streamfunction is a nonlocal function of the properties of the fluid column. It is constructed to have a low baroclinic mode vertical structure and to smoothly transition through regions of weak stratification such as boundary layers or mode waters. It requires no matching conditions or regularization in unstratified regions; it satisfies boundary conditions of zero transport at the ocean surface and bottom; and it provides a sink of available potential energy for each vertical seawater column, but not necessarily at each location within the column. Numerical implementation of the methodology requires the solution of a one-dimensional tridiagonal problem for each vertical column. To illustrate the approach, we present an analytical example based on the nonlinear Eady problem and two numerical simulations.     

Interaction of waves with non-colinear currents

The present paper reports on a study of the interaction of a current-free monochromatic surface wave field with a wave-free uniform current field in a three-dimensional flow frame. The wave and the current fields are not necessarily collinear with each other. The formulation of the wave-current field is done under the assumption of irrotational and inviscid flow. We have developed the three dimensional expressions describing the characteristics of the combined flow in terms of mass, momentum and energy transport conservation. These equations are found efficient to describe the sought-for combined wave-current field. The parameters describing the wave-current field after the interaction are the surface disturbance amplitude and length, mean water depth, mean current-like parameter and direction of the combined flow, which would be calculated from a set of equations that satisfy conservation of mean mass, momentum and energy flux and a dispersion relation on the free surface before and after the interaction. The results are shown in terms of relative changes in wave heights and lengths, current-like parameters and final directions obtained for the combined wave-current field with respect to current-free wave and wave-free current pre-interaction parameters.     

含强水流高阶Boussinesq水波方程

采用摄动法并利用已建立的纯波情况下高阶Boussinesq方程,建立了可以考虑强水流与波浪相互作用的高阶Boussinesq方程.水流速度与波浪群速具有相同量级,且随时间和空间的变化尺度远大于波浪周期和波长.方程色散性近似到[4/4]阶Pade展开,对浅水情况方程可以是完全非线性的,可适用于波流相互作用的强非线性问题.通过将水流存在时波长和波幅的结果与一阶斯托克斯波结果对比,讨论了具有不同近似程度的3种含波流相互作用的Boussinesq方程的适用性.     

A North Atlantic deep-water eddy in the Agulhas Current system

One hundred and twelve stations of CTDO2 and LADCP were collected in the Agulhas Current system as part of the Agulhas Undercurrent experiment (AUCE) in March 2003. Along an offshore section, at approximately 35.6°S and 27.3°E to the northwest of the tip of the Agulhas Plateau, an unusual feature was revealed between 2200 and 3500 m depth, imbedded in the northward moving NADW layer. An anomalously high salinity of 34.83, 0.03 saltier than the surrounding water, was observed. Maximums in the potential temperature and oxygen were also found, with isotherms dropping by about 250 m over 50 km and a doming of the oxygen layers. From the convex lens structure of the neutral surfaces, we conclude that we sampled an anticyclonic eddy of NADW. Since the LADCP data reveal deep velocities up to 20 cm s⁻¹, yet no anticyclonic circulation, whereas the geostrophic velocity referenced to the bottom shows a weak anticyclonic circulation, we inferred that we sampled the outer edge of the eddy and not its core. From an analysis of the water properties within the eddy and a comparison with known properties in the SE Atlantic Ocean and SW Indian Ocean, we conclude that the eddy was formed in the Agulhas Retroflection region. We speculate that the eddy was the result of an instability in the NADW slope current, which flows from the SE Atlantic around the Agulhas Bank. A deeply penetrating Agulhas Ring spun up the deep waters, pinching off an eddy, which later detached from the slope current and was carried southward. Once offshore, it coupled with the surface Agulhas Return Current, whose meandering path advected the eddy northeastward and ejected it over the Agulhas Plateau.     

A sheetflow sediment transport model for skewed-asymmetric waves combined with strong opposite currents

Prototype scale physical model tests were conducted to investigate the sheetflow sediment transport of uniform sand under different skewed-asymmetric oscillatory flows with and without the presence of relatively strong currents in the opposite direction against wave propagation. Experiments show that in most cases with fine sands, the “cancelling effect” which balances the on-/off-shore net transport under pure asymmetric/skewed oscillatory flows and results a moderate net transport was developed for combined skewed-asymmetric shaped oscillations. However, under certain conditions (T > 5 s) with coarse sands, the onshore sediment transport was enhanced for combined skewed-asymmetric flows. Additionally, the new experimental data under collinear oscillatory flows and strong currents show that offshore net transport rates increase with decreasing velocity skewness and acceleration skewness. Sediment movement behaviors were investigated through analysis of experimental data obtained from the image analysis technique and attempts were made to estimate and formulate the sheetflow layer thickness. Accordingly, sediment transport under oscillatory sheetflow conditions was studied and successfully explained by comparing the bed shear stress and the phase lag parameter at each half cycle. Consequently, these parameters were incorporated in an improved Dibajinia and Watanabe's type sediment transport model. The formula is calibrated against a comprehensive experimental data (331 in total). Good agreement obtained between predictions and measurements shows that the new formula is fulfilled for practical purposes.     

Tracing nearshore bottom currents with sea-bed drifters

The experience acquired during 36 Woodhead drifter releases performed along the Tuscan coast is briefly reported. Without a costly apparatus, residual current patterns may be easily obtained with this method. Care should be taken, however, when sediment dispersal patterns are required, since weak currents which drag drifters slowly, displace only finer sediments.     

A Lagrangian biogeochemical study of an eddy in the Northeast Atlantic

We report the results of an experiment in the Northeast Atlantic in which sulphur hexafluoride (SF₆) was released within an eddy and the behaviour of trace gases, nutrients and productivity followed within a Lagrangian framework over a period of 24 days. Measurements were also made in the air above the eddy in order to estimate air-sea exchange rates for some components. The physical, biological and biogeochemical properties of the eddy resemble those of other eddies studied in this area, suggesting that the results we report may be applicable beyond the specific eddy studied. During a period of low wind speed at the start of the experiment, we are able to quantitatively describe and balance the nutrient and carbon budgets for the eddy. We also report concentrations of various trace gases in the region which are similar to those observed in other studies and we estimate exchange rates for several trace gases. We show that the importance of gas exchange over other loss terms varies with time and also varies for the different gases. We show that the various trace gases considered (CO₂, dimethyl sulphide (DMS), N₂O, CH₄, non-methane-hydrocarbons, methyl bromide, methyl iodide and volatile selenium species) are all influenced by physical and biological processes, but the overall distribution and temporal variability of individual gases are different to one another. A storm disrupted the stratification in the eddy during the experiment, resulting in enhanced nutrient supply to surface waters, enhanced gas exchange rates and a change in plankton community, which we quantify, although overall productivity was little changed. Emphasis is placed on the regularity of storms in the temperate ocean and the importance of these stochastic processes in such systems.