The purpose of the thesis is to analyze the properties of acoustic sound field in the Sizihwan Bay Marine Test Field (SBMTF), and to conduct geoacoustic inversion using the measured data. Two experiments were carried out during February, 2009, and March, 2010. The source was UW350, and the receivers were ITC 6050 hydrophones. The transmitted frequencies lie between 350 Hz and 1250 Hz. On the analysis of acoustic propagation, the study calls for the application of OASES for the analysis of transmission loss and effects of range dependency. Moreover, the uncertainty analysis due to environmental factors was carried out based upon probabilistic approach. For the inversion analysis, the sensitivity of each environmental parameter was first analyzed, and those parameters with high sensitivity were chosen for inversion. Inversion was conducted by the application of SAGA. The results have shown that, due to the fact that the SBMTF is an very shallow water environment with water depth less than 30 m, the acoustic sound field is strongly affected by the boundaries. From the uncertainty analysis, it shows that the transmission loss has high degree of uncertainty resulting from seabed property due to its lack of accurate measurement. The inversion results for water depth and seabed sound speed were obtained, and compared with the measured data. This research calls for experimental design, data processing, software application, and result analysis, offering an overall understanding of the properties of SBMTF that is valuable for future study.

Originally developed in the 1970s for the optimal control of systems governed by partial differential equations, the adjoint method has found several successful applications, e.g., in meteorology with large-scale 3D or 4D atmospheric data assimilation schemes, for carbon cycle data assimilation in biogeochemistry and climate research, or in oceanographic modelling with efficient adjoint codes of ocean general circulation models. Despite the variety of applications in these research fields, adjoint methods have only very recently drawn attention from the ocean acoustics community. In ocean acoustic tomography and geoacoustic inversion, where the inverse problem is to recover unknown acoustic properties of the water column and the seabed from acoustic transmission data, the solution approaches are typically based on travel time inversion or standard matched-field processing in combination with metaheuristics for global optimization. In order to complement the adjoint schemes already in use in meteorology and oceanography with an ocean acoustic component, this thesis is concerned with the development of the adjoint of a full-field acoustic propagation model for shallow water environments. In view of the increasing importance of global ocean observing systems such as the European Seas Observatory Network, the Arctic Ocean Observing System and Maritime Rapid Environmental Assessment (MREA) systems for defence and security applications, the adjoint of an ocean acoustic propagation model can become an integral part of a coupled oceanographic and acoustic data assimilation scheme in the future. Given the acoustic pressure field measured on a vertical hydrophone array and a modelled replica field that is calculated for a specific parametrization of the environment, the developed adjoint model backpropagates the mismatch (residual) between the measured and predicted field from the receiver array towards the source. The backpropagated error field is then converted into an estimate of the exact gradient of the objective function with respect to any of the relevant physical parameters of the environment including the sound speed structure in the water column and densities, compressional/shear sound speeds, and attenuations of the sediment layers and the sub-bottom halfspace. The resulting environmental gradients can be used in combination with gradient descent methods such as conjugate gradient, or Newton-type optimization methods tolocate the error surface minimum via a series of iterations. This is particularly attractive for monitoring slowly varying environments, where the gradient information can be used to track the environmental parameters continuously over time and space. In shallow water environments, where an accurate treatment of the acoustic interaction with the bottom is of outmost importance for a correct prediction of the sound field, and field data are often recorded on non-fully populated arrays, there is an inherent need for observation over a broad range of frequencies. For this purpose, the adjoint-based approach is generalized for a joint optimization across multiple frequencies and special attention is devoted to regularization methods that incorporate additional information about the desired solution in order to stabilize the optimization process. Starting with an analytical formulation of the multiple-frequency adjoint approach for parabolic-type approximations, the adjoint method is progressively tailored in the course of the thesis towards a realistic wide-angle parabolic equation propagation model and the treatment of fully nonlocal impedance boundary conditions. A semi-automatic adjoint generation via modular graph approach enables the direct inversion of both the geoacoustic parameters embedded in the discrete nonlocal boundary condition and the acoustic properties of the water column. Several case studies based on environmental data obtained in Mediterranean shallow waters are used in the thesis to assess the capabilities of adjoint-based acoustic inversion for different experimental configurations, particularly taking into account sparse array geometries and partial depth coverage of the water column. The numerical implementation of the approach is found to be robust, provided that the initial guesses are not too far from the desired solution, and accurate, and converges in a small number of iterations. During the multi-frequency optimization process, the evolution of the control parameters displays a parameter hierarchy which clearly relates to the relative sensitivity of the acoustic pressure field to the physical parameters. The actual validation of the adjoint-generated environmental gradients for acoustic monitoring of a shallow water environment is based on acoustic and oceanographic data from the Yellow Shark '94 and the MREA '07 sea trials, conducted in the Tyrrhenian Sea, south of the island of Elba. Starting from an initial guess of the environmental control parameters, either obtained through acoustic inversion with global search or supported by archival in-situ data, the adjoint method provides an efficient means to adjust local changes with a couple of iterations and monitor the environmental properties over a series of inversions. In this thesis the adjoint-based approach is used, e.g., to fine-tune up to eight bottom geoacoustic parameters of a shallow-water environment and to track the time-varying sound speed profile in the water column. In the same way the approach can be extended to track the spatial water column and bottom structure using a mobile network of sparse arrays. Work is currently being focused on the inclusion of the adjoint approach into hybrid optimization schemes or ensemble predictions, as an essential building block in a combined ocean acoustic data assimilation framework and the subsequent validation of the acoustic monitoring capabilities with long-term experimental data in shallow water environments.

Single-beam echosounders are an inexpensive, practical and non-invasive means of remote sensing the seabed. Ideally, the common single-beam echosounder should be able to tell fishers, navigators, engineers and scientists what the seabed consists of in addition to water depth. Low-frequency underwater acoustic systems (<10 kHz) can do this in some circumstances, but are expensive, offer limited resolution and potentially hazardous to marine mammals. High-frequency systems, such as single and multibeam echosounders, are very effective at mapping bathymetry, but do not characterize the seabed directly. Instead, these systems divide the seabed into self-similar segments or classes, and then rely on ground-truth data (usually sediment grab samples) to assign seabed-type labels such as sand, etc., to the classes. However, inadequate and inaccurate ground-truth is a major problem. Single-beam seabed classification methods also suffer from a lack of discriminatory power and from artefacts such as water depth and seabed slope. The cause of these problems is that the methods lack a basis in physics and are mainly statistical. Then, the central objective in this dissertation is to develop physics-based methods to improve classification and to address the problem of ground-truth by inferring seabed characteristics directly from the acoustics. An overview of current methods is presented along with case studies of single-beam surveys to introduce the current seabed classification method called QTC VIEW™ and to identify specific problems. A physical basis is established in scattering and geometrical theories and observations of field and model data. This leads to new classification and characterization methods that overcome the shortcomings of current seabed classification methods. Advancements also include new physical models of echosounding. The new methods are presented, implemented and evaluated. Highlights of experimental results include a new testbed located in Patricia Bay, British Columbia. The testbed consists of exhaustive ground-truth, surveys and novel controlled experiments with various single-beam echosounders, ranging in frequency from 12 to 200 kHz. Simulated echo time series data from the numerical BORIS model and a new analytic model are used to augment the testbed. Evaluation of experimental results shows the new physics-based methodology improves seabed classification significantly and enables seabed characterization by an uncalibrated single-beam echosounder. / Graduate

In shallow water environments, sound propagation experiences multiple interactions with the surface/bottom interfaces, with hydrodynamic disturbances such as internal waves, and with tides and fronts. It is thus very difficult to make satisfactory predictions of sound propagation in shallow water. Given that many of the ocean characteristics can be modeled as stochastic processes, the statistical measure, spatial coherence, is consequently an important quantity. Spatial coherence provides valuable information for array performance predictions. However, for the case of long-range, low frequency propagation, studies of spatial coherence influenced by various environmental parameters are limited insofar as having the appropriate environmental data with which to model and interpret the results. The comprehensive Asian Seas International Experiment 2001 (ASIAEX01) examined acoustic propagation and scattering in shallow water. Environmental oceanographic data were taken simultaneously with the acoustic data. ASIAEX01 provided a unique data set which enabled separate study of the characteristics of the oceanographic features and their influence on long range sound propagation. In this thesis, the environmental descriptors considered include sediment sound speed and attenuation, background internal waves, episodic non-linear internal waves, and air-sea interface conditions. Using this environmental data, the acoustic data are analyzed to show the characteristics of spatial coherence in a shallow water waveguide. It is shown that spatial coherence can be used as an inversion parameter to extract geoacoustic information for the seabed. Environmental phenomena including internal waves and wind-generated surface waves are also studied. The spatial and temporal variations in the sound field induced by them are presented. In addition, a tank experiment is presented which simulates propagation in a shallow water waveguide over a short range. Based on the data model comparison results, the model proposed here is effective in addressing the major environmental effects on sound propagation in shallow water.

The objective of this thesis is to develop and implement an algorithm for inverting the sound speed profile via estimation of the parameters embedded in a geoacoustic model. The environmental model inscribes a continuously-varying marine sediment layer with density and sound speed distributions represented by the generalized-exponential and inverse-square functions, respectively. Based upon a forward problem of plane-wave reflection from a non-uniform sediment layer overlying a uniform elastic basement, an inversion procedure for estimating the sound speed profile from the reflected sound field under the influence of noise is established and numerically implemented. The inversion invokes a probabilistic approach quantified by the posterior probability density for measuring the uncertainties of the estimated parameters from synthetic noisy data. Preliminary analysis on the solution of the forward problem and the sensitivity of the model parameters is first conducted, leading to a determination of the parameters chosen for inversion in the ensuing study. The parameter uncertainties referenced 1-D and 2-D marginal posterior probability densities are then examined, followed by the statistical estimation for the sound speed profile in terms of 99 % credibility interval. The effects of, the signal-to-noise ratio (SNR), the dimension of data vector, the region in which the data sampled, on the statistical estimation of sound speed profile are demonstrated and discussed.

The overall goal of this thesis is to develop non-linear Bayesian methods for three-dimensional tracking of a moving acoustic source in shallow water despite environmental uncertainty, with application to data from a horizontal line array (HLA) of hydrophones. As a precursor, Bayesian geoacoustic inversion is applied to estimate seabed model parameters and their uncertainties. A simulation study examines the effect of source and array factors on geoacoustic information content in matched-field inversion of HLA data, as quantified in terms of model parameter uncertainties. Bayesian geoacoustic inversion is applied to both controlled-source and ship-noise data from a HLA deployed on the seafloor in a shallow-water experiment conducted in the Barents Sea. A new approach is introduced to account for data error reduction due to averaging data over time-series subsegments (snapshots), based on empirically apportioning measurement and theory error, with effects on inversion results compared to those of existing approaches. It is further demonstrated that combining data from multiple, independent time-series segments (for a moving source) in the inversion can significantly reduce geoacoustic parameter uncertainties. Geoacoustic uncertainties are also shown to depend on ship range and orientation, with lowest uncertainties for short ranges and for the ship stern/propeller oriented toward the array. Sediment sound-speed profile and density estimates from controlled-source and ship-noise data inversions are found to be in good agreement with values from geophysical measurements. Two non-linear Bayesian matched-field inversion approaches are developed for three-dimensional source tracking despite environmental uncertainty. Focalization-tracking maximizes the posterior probability density (PPD) over track and environmental parameters. Synthetic test cases show that the algorithm substantially outperforms tracking with poor environmental estimates and generally obtains results close to those achieved with exact environmental knowledge. Marginalization-tracking integrates the PPD over environmental parameters to obtain joint marginal distributions over source coordinates, from which track uncertainty estimates and the most probable track are extracted. Both approaches are applied to data from the Barents Sea experiment. Focalization-tracking successfully estimates the tracks of the towed source and a surface ship in cases where simpler tracking algorithms fail. Marginalization-tracking generally outperforms focalization-tracking and gives uncertainty estimates that encompass the true tracks.

La connaissance de l'environnement marin est nécessaire pour un grand nombre d'applications dans le domaine de l'acoustique sous-marine comme la communication, la localisation et détection sonar et la surveillance des mammifères marins. Il constitue le moyen principal pour éviter les interférences néfastes entre le milieu naturel et les actions industriels et militaires conduites en zones côtières.Notre travail de thèse se place dans un contexte de sonar actif avec des fréquences allant de 1 kHz à 10 kHz pour des distances de propagations allant de 1 km à plusieurs dizaines de kilomètres. Nous nous intéressons particulièrement aux environnements de propagation grands fonds, à l'utilisation des antennes industrielles comme les antennes de flancs, les antennes cylindriques et les antennes linéaires remorquées, et à l'utilisation de signaux large bande afin de travailler avec des résolutions en distance et en vitesse très élevées. Le travail de recherche présenté dans ce mémoire est dédié à la recherche de nouveaux paramètres discriminants pour la classification de cible sous-marine en sonar actif et notamment à l'estimation de l'immersion instantanée.Cette étude présente : (1) les calculs de nouvelles bornes de Cramer-Rao pour la position d'une cible en distance en et en profondeur, (2) l'estimation conjointe de la distance et de l'immersion d'une cible à partir de la mesure des temps d'arrivées et des angles d'élévations sur une antenne surfacique et (3) l'estimation conjointe de la distance, de l'immersion et du gisement d'une cible à partir de la mesure des temps d'arrivées et des pseudo-gisements sur une antenne linéaire remorquée.Les méthodes développées lors de cette étude ont été validées sur des simulations, des données expérimentales à petite échelle et des données réelles en mer. / The knowledge of the marine environment is required for many underwater applications such as communications, sonar localization and detection, and marine mammals monitoring. It enables preventing harmful interference between the natural environment and industrial and military actions in coastal areas.This thesis work concentrates upton the context of active sonar with frequencies from 1 kHz to 10 kHz and long propagation ranges from 1 km to several tens of kilometers. We also concentrates upon deep water environment, the use of industrial arrays such as cylindrical arrays, flank arrays and linear towed arrays, and the use of large time-bandwidth signals in order to obtain high distance and speed resolutions. This research work is dedicated to the research of new features for the underwater target classification in active sonar, and specifically to the instantaneous target-depth estimation.This thesis presents: (1) calculations of new Cramer-Rao bounds for the target-position in range and in depth, (2) the joint estimation of the target-depth and the target-range from the arrival time and elevation angle measures with a surface array, (3) the joint estimation of the target-depth, the target-range and the target-bearing from the arrival time and pseudo-bearing angle measures with a linear towed array.The methods presented in this manuscript have been benchmarked on simulation, on reduced-scale experimental data and real marine data.

Waterborne acoustic signals carry information about the ocean environment. Ocean geoacoustic inversion is the task of estimating environmental parameters from received acoustic signals by matching the measured sound with the predictions of a physics-based model. A lower bound on the uncertainty associated with environmental parameter estimates, the Cramér-Rao bound, can be calculated from the Fisher information, which is dependent on derivatives of a physics-based model. Physics-based preconditioners circumvent the need for variable step sizes when computing numerical derivatives. This work explores the feasibility of using a neural network to perform geoacoustic inversion for environmental parameters and their associated uncertainties from ship noise spectrogram data. To train neural networks, a synthetic dataset is generated and tested for generalizability against 31 measurements taken during the SBCEX2017 study of the New England Mud Patch.

In this thesis, acoustic inversion methods are employed to estimate array element locations and the geoacoustic properties of the seabed using measured acoustic data consisting of noise from a surface ship in the Gulf of Mexico. The array element localization utilizes relative travel-time information obtained by cross-correlating the recorded time series of ship noise received at spatially separated hydrophones. The relative travel-time data are used in an inversion, based on the regularized least-squares method and the acoustic ray tracing equations, to obtain improved estimates of the receiver and source positions and their uncertainties. Optimization and Bayesian matched-field inversion methods are employed to estimate seabed geoacoustic properties and their uncertainties in the vicinity of a bottom-moored vertical line array using the recorded surface ship noise. This study is used to test the feasibility of matched-field methods to detect temporal changes in the geoacoustic properties of the seabed near a known gas hydrate mound in the Gulf of Mexico. Finally, a synthetic study is performed that demonstrates how ignoring environmental range dependence of seabed sound speed and water depth in matched-field inversion can lead to biases in the estimated geoacoustic parameters. The study considers the distributions of optimal parameter estimates obtained from a large number of range-independent inversions of synthetic data generated for random range-dependent environments. Range-independent Bayesian inversions are also performed on selected data sets and the marginal parameter distributions are examined. Both hard- and soft-bottom environments are examined at a number of scales of variability in sound speed and water depth.

Ce travail de thèse propose d’analyser les signaux issus d’une source omnidirectionnelle sphérique réfléchis par un milieu sédimentaire stratifié et enregistré par une antenne d’hydrophones, en vue de caractériser quantitativement les sédiments marins aux moyennes fréquences, i.e. comprises entre 1 et 10 kHz. La recherche développée dans ce manuscrit propose une méthodologie facilitant la recherche des paramètres géoacoustiques du milieu avec la méthode des sources images, ainsi qu’un ensemble de solutions techniques appropriées afin d’améliorer cette méthode d’inversion récemment développée. La méthode des sources images repose sur une modélisation physique de la réflexion des ondes émises par une source sur un milieu stratifié sous l’approximation de Born. Par conséquent, la réflexion de l’onde sur le milieu stratifié peut être représentée par une collection de sources images, symétriques de la source réelle par rapport aux interfaces, dont les positions spatiales sont liées à la vitesse des ondes acoustiques et aux épaisseurs des couches. L’étude se décline en deux volets : traitements des signaux et inversion des paramètres géoacoustiques. La première partie du travail est focalisée sur le développement de la méthode des sources images. La méthode originelle se basait sur la construction de cartes de migration et de semblance de signaux pour déterminer les paramètres d’entrée de l’algorithme d’inversion qui sont des temps de trajet et des angles d’arrivée. Afin d’éviter cette étape, nous détectons les temps d’arrivée avec l’opérateur d’énergie de Teager-Kaiser (TKEO) et nous trouvons les angles par une méthode de triangulation. Le modèle d’inversion a été ensuite intégré en prenant en compte la possibilité de déformation de l’antenne. Cette partie se termine par une nouvelle approche qui combine TKEO et des méthodes temps fréquence afin d’avoir une bonne détection du temps d’arrivée dans le cas de signaux fortement bruités. Sur le plan du modèle et de l’inversion géoacoustique, nous proposons tout d’abord une description précise du modèle direct en introduisant le concept de sources images virtuelles. Cette étape permet de mieux comprendre l’approche développée. Ensuite, nous proposons une extension de la méthode des sources image pour l’inversion de paramètres géoacoustiques supplémentaires : la densité, l’atténuation et la vitesse des ondes de cisaillement. Cette extension est basée sur les résultats de l’inversion originelle (estimation du nombre de strates, de leur épaisseur, et de la vitesse des ondes de compression) ainsi que sur l’utilisation de l’amplitude des signaux réfléchis. Ces améliorations et extensions de la méthode des sources images sont illustrées par leur application sur des signaux synthétiques et des signaux réels issus d’expérimentations en cuve et à la mer. Les résultats obtenus sont très satisfaisants, tant au niveau des performances de calcul que de la qualité des estimations fournies. / This thesis aims at analyzing the signals emitted from a spherical omnidirectional source reflected by a stratified sedimentary environment and recorded by a hydrophone array in order to characterize quantitatively the marine sediments at medium frequencies, i.e. between 1 and 10 kHz. The research developed in this manuscript provides a methodology to facilitate the estimation of medium geoacoustic parameters with the image source method, and some appropriate technical solutions to improve this recently developed inversion method. The image source method is based on a physical modeling of the wave reflection emitted from a source by a stratified medium under the Born approximation. As result, the reflection of the wave on the layered medium can be represented by a set of image sources, symmetrical to the real source with respect to the interfaces, whose spatial positions are related to the sound speeds and the thicknesses of the layers. The study consists of two parts : signal processing and inversion of geoacoustic parameters. The first part of the work is focused on the development of the image source method. The original method was based on migration and semblance maps of the recorded signals to determine the input parameters of the inversion algorithm which are travel times and arrival angles. To avoid this step, we propose to determine the travel times with the Teager-Kaiser energy operator (TKEO) and the arrival angles are estimate with a triangulation approach. The inversion model is then integrated, taking into account the possible deformation of the antenna. This part concludes with a new approach that combines TKEO and time-frequency representations in order to have a good estimation of the travel times in the case of noisy signals. For the modeling and geoacoustic inversion part, we propose first an accurate description of the forward model by introducing the concept of virtual image sources. This idea provides a deeper understanding of the developed approach. Then, we propose an extension of the image sources method to the estimation of supplementary geoacoustic parameters : the density, the absorption coefficient, and the shear wave sound speed. This extension is based on the results of the original inversion (estimation of the number of layers, their thicknesses, and the pressure sound speeds) and on the use of the amplitudes of the reflected signals. These improvements and extents of the image source method are illustrated by their applications on both synthetic and real signals, the latter coming from tank and at-sea measurements. The obtained results are very satisfactory, from a computational point of view as well as for the quality of the provided estimations.