Direct measurement and matched-field inversion approaches to arrayshape estimation

Accurate knowledge of array shape is essential for carrying out full wavefield (matched-field) processing. Direct approaches to array element localization (AEL) include both nonacoustic (tilt-heading sensors) and acoustic (high-frequency, transponder-based navigation) methods. The low-frequency signature emitted from a distant source also can be used in an inversion approach to determine array shape. The focus of this paper is on a comparison of the array shape results from these three different methods using data from a 120-m aperture vertical array deployed during SWellEx-3 (Shallow Water evaluation cell Experiment 3). Located 2 m above the shallowest array element was a self-recording package equipped with depth, tilt, and direction-of-tilt sensors, thereby permitting AEL to be performed non-acoustically. Direct AEL also was performed acoustically by making use of transponder pings (in the vicinity of 12 kHz) received by high-frequency hydrophones spaced every 7.5 m along the vertical array. In addition to these direct approaches, AEL was carried out using an inversion technique where matched-field processing was performed on a multitone (50-200 Hz), acoustic source at various ranges and azimuths from the array. As shown, the time-evolving array shape estimates generated by all three AEL methods provide a consistent picture of array motion throughout the 6-h period analyzed     

Sensor position estimation and source ranging in a shallow water environment

Underwater acoustic transient signals are generated mechanically at known positions along a wharf. These signals are received by a wide aperture planar array of four underwater acoustic sensors, whose positions relative to the wharf are unknown. A method is described that enables the positions of the sensors to be estimated from accurate differential time-of-arrival measurements (with 0.1 /spl mu/s precision) as the signal wavefronts traverse the array. A comparison of the estimated positions with the nominal positions of the first three sensors, which form a 20-m-wide aperture horizontal line array, reveals a 2-cm displacement of the middle sensor from the line array axis. This slight bowing of the line array results in overranging (bias error of 3%) when the wavefront curvature method is used with the nominal collinear sensor positions to locate a static source of active sonar transmissions at a range of 59.2 m. The use of the spherical intersection method coupled with the estimated sensor positions of the line array provides an order of magnitude improvement in the range estimate (within 0.3% of the actual value). However, systematic ranging errors are observed when the sound propagation medium becomes nonstationary. Next, the differences in the arrival times of the direct path and boundary-reflected path signals at the middle sensor of the wide aperture line array are estimated using the differential phase residue of the analytic signal at the sensor output. These multipath delays are used to estimate the range and depth of the source. Although the average value of the multipath range estimates is within 0.5% of the actual value, the variance of the range estimates is 50 times larger when compared with the results of the spherical intersection and wavefront curvature methods. The multipath delay data are also processed to provide a reliable estimate of the temporal variation in the water depth enabling the tidal variation to be observed.     

Centimeter-Level Positioning of Seafloor Acoustic Transponders from a Deeply-Towed Interrogator

An array of three seafloor transponders was acoustically surveyed to centimeter precision with a deeply-towed interrogator. Measurements of two-way acoustic travel time and hydrostatic pressure made as the interrogator was towed above the array were combined in a least-squares adjustment to estimate the interrogator and transponder positions in two surveys spanning two years. No transponder displacements were expected at this site in the interior of the Juan de Fuca Plate (48^?11′ N, 127^?12′ W) due to the lack of active faults. This was confirmed to a precision of ±2 cm by least-squares adjustment. Marginally detectable blunders in the observations were shown to affect the transponder position estimates by no more than 3 mm, demonstrating the geometric strength of the data set. The accumulation of many hundreds of observations resulted in a significant computational burden on the least-squares inversion procedure. The sparseness of the normal matrix was exploited to reduce by a factor of 1000 the number of calculations. The acoustic survey results suggested that the near-bottom sound speed fields during the two surveys were in better agreement than inferred from yearly single-profile conductivity, temperature, and pressure (CTD) measurements.     

Geoacoustic seafloor exploration with a towed array in a shallowwater area of the Strait of Sicily

Acoustic propagation in shallow water is greatly dependent on the geoacoustic properties of the seabottom. This paper exploits this dependence for estimating geoacoustic sediment properties from the bottom acoustic returns of known signals received on a hydrophone line array. There are two major issues in this approach: one is the feasibility of acoustic inversion with a limited aperture line array, the other is related to the knowledge of the geometry of the experimental configuration. To test the feasibility of this approach, a 40-hydrophone-4-m spaced towed array together with a low-frequency acoustic source, was operated at a shallow water site in the Strait of Sicily. In order to estimate the array deformation in real time, it has been equipped with a set of nonacoustic positioning sensors (compasses, tilt-meters, pressure gauges). The acoustic data were inverted using two complementary approaches: a genetic algorithm (GA) like approach and a radial basis functions (RBF) inversion scheme. More traditional methods, based on core sampling, seismic survey and geophone data, together with Hamilton's regression curves, have also been employed on the same tracks, in order to provide a ground truth reference environment. The results of the experiment, can be summarized as follows: 1) the towed array movement is not negligible for the application considered and the use of positioning sensors are essential for a proper acoustic inversion, 2) the inversion with GA and RBF are in good qualitative agreement with the ground truth model, and 3) the GA scheme tends to have better stability properties. On the other hand, repeated in version of successive field measurements requires much less computational effort with RBF     

Estimates of geoacoustic model parameters from inversions of horizontal and vertical line array data

This paper describes results from geoacoustic inversion of low-frequency acoustic data recorded at a receiving array divided into two sections, a sparse bottom laid horizontal array (HLA) and a vertical array (VLA) deployed in shallow water. The data are from an experiment conducted by the Norwegian Defence Research Establishment (FFI) in the Barents Sea, using broadband explosives (shot) sources. A two-layer range-independent geoacoustic model, consistent with seismic profiles from the area, described the environment. Inversion for geoacoustic model parameters was carried out using a fast implementation of the hybrid adaptive simplex simulated annealing (ASSA) inversion algorithm, with replica fields computed by the ORCA normal mode code. Low-frequency (40-128 Hz) data from six shot sources at ranges 3-9 km from the array were considered. Estimates of sediment and substrate p-wave velocities and sediment thickness were found to be consistent between independent inversions of data from the two sections of the array.     

Geoacoustic model parameter estimation using a bottom-moored hydrophone array

This paper describes results from an experiment carried out to investigate geoacoustic inversion with a bottom-moored hydrophone array located in the shallow waters of the Timor Sea off the northern coast of Australia. The array consisted of two arms in a V shape, horizontally moored at a site that was essentially flat over a large area. Hydrophone positions were estimated using an array element localization (AEL) technique that established relative uncertainties of less than 1 m on the seafloor. The data used for geoacoustic inversion were from experiments with continuous wave (CW) tones in the 80- to 195-Hz band transmitted from a towed projector. A hybrid search algorithm determined the set of geoacoustic model parameters that maximized the Bartlett fit (averaged coherently spatially at each tone and incoherently over frequency) between the measured and modeled data at the array. Due to the long range experimental geometry, the inversion was sensitive to attenuation in the sediment. The inverted geoacoustic profile performed well in a simple test for localizing the sound source at other sites in the vicinity of the array. Range-depth localization performance for the horizontal array was comparable to that for an equivalent vertical array.     

Geoacoustic Information Content of Horizontal Line Array Data

This paper examines the effectiveness of horizontal line arrays (HLAs) for matched-field inversion (MFI) by quantifying geoacoustic information content for a variety of experiment and array factors, including array length and number of sensors, source range and bearing, source-frequency content, and signal-to-noise ratio (SNR). Emphasis is on bottom-moored arrays, while towed arrays are also considered, and a comparison with vertical line array (VLA) performance is made. The geoacoustic information content is quantified in terms of marginal posterior probability distributions (PPDs) for model parameters estimated using a fast Gibbs sampler approach to Bayesian inversion. This produces an absolute, quantitative estimate of the geoacoustic parameter uncertainties which can be directly compared for various experiment and array factors.     

Joint direction-of-arrival and array-shape tracking for multiplemoving targets

The problem of tracking the directions-of-arrival (DOAs) of multiple moving sonar targets with an array of passive sensors is complicated by sensor movement. An algorithm for the joint tracking of source DOAs and sensor positions is presented to address this problem. Initial maximum-likelihood estimates of source DOAs and sensor positions are refined by Kalman filtering. Spatio-temporally correlated array movement is considered. Source angle dynamics are used to achieve correct data association. The new technique is capable of performing well for the difficult cases of sources that cross in angle as well as for fully coherent sources. Computer simulations show that the approach is robust in the presence of array motion modeling uncertainty and effectively reduces dependence on expensive and possibly unreliable hardware     

Accuracy assessment of GPS/Acoustic positioning using a Seafloor Acoustic Transponder System

The accuracy of GPS/Acoustic positioning is crucial for monitoring seafloor crustal deformation. However, the slant range residual is currently the only indicator used to evaluate the precision of positioning seafloor transponders. This study employs a unique Seafloor Acoustic Transponder System (SATS) to evaluate the accuracy of GPS/Acoustic seafloor positioning. The SATS has three transponders and an attitude sensor in a single unit, which provides true lengths of transponder baselines and true attitude of the SATS to ensure assessment reliability and validity. The proposed approach was tested through a GPS/Acoustic experiment, in which an off-the-shelf acoustic system was used to collect range measurements. Using GPS/Acoustic geodetic observations, the positions of three transponders on the SATS were estimated by an optimization technique combined with ray-tracing calculations. The accuracy of the GPS/Acoustic seafloor positioning is assessed by comparing the true baselines and attitude with the results derived from the position estimates of the three transponders. A sensitivity analysis is conducted to investigate the robustness of the GPS/Acoustic positioning results to changes of sound speed. Experimental results demonstrate that the use of the SATS can help to assess the validity of the GPS and acoustic travel time measurements in the GPS/Acoustic seafloor positioning.     

New buoy observation system for tsunami and crustal deformation

We have developed a new system for real-time observation of tsunamis and crustal deformation using a seafloor pressure sensor, an array of seafloor transponders and a Precise Point Positioning (PPP ) system on a buoy. The seafloor pressure sensor and the PPP system detect tsunamis, and the pressure sensor and the transponder array measure crustal deformation. The system is designed to be capable of detecting tsunami and vertical crustal deformation of ±8 m with a resolution of less than 5 mm. A noteworthy innovation in our system is its resistance to disturbance by strong ocean currents. Seismogenic zones near Japan lie in areas of strong currents like the Kuroshio, which reaches speeds of approximately 5.5 kt (2.8 m/s) around the Nankai Trough. Our techniques include slack mooring and new acoustic transmission methods using double pulses for sending tsunami data. The slack ratio can be specified for the environment of the deployment location. We can adjust slack ratios, rope lengths, anchor weights and buoy sizes to control the ability of the buoy system to maintain freeboard. The measured pressure data is converted to time difference of a double pulse and this simple method is effective to save battery to transmit data. The time difference of the double pulse has error due to move of the buoy and fluctuation of the seawater environment. We set a wire-end station 1,000 m beneath the buoy to minimize the error. The crustal deformation data is measured by acoustic ranging between the buoy and six transponders on the seafloor. All pressure and crustal deformation data are sent to land station in real-time using iridium communication.     

Monitoring moored instrument motion by optimal estimation

Many oceanographic applications require the positioning of the underwater sensor at measurement times. We consider here the case of subsurface moored tomographic instruments, where the distance between source and receiver must be known within a few meters. For that purpose, a long baseline array is deployed: this system includes a navigator, attached to the mooring element and an array of three transponders set on the ocean bottom. To process the navigation data collected with such system, we have developed a method based on optimal estimation. The triangulation problem is not a basic spherical constraints one and the specificity of deep underwater positioning, related to the variability of the ocean sound speed profile are pointed out. Correcting terms are proposed and introduced into the system. Simultaneous inversion of all data, defining an overconstrained problem allows to estimate biases and errors. The algorithm is applied here to a dataset collected in the Azores-Canary basin during CAMBIOS experiment.     

Methods of array element localization for a towed underwateracoustic array

The location of the hydrophones on a towed underwater acoustic array as a function of time (array element localization) is needed for signal processing. Methods to perform this localization using least squares polynomial fitting to data from depth sensors, heading sensors, and sensors detecting a ping from a single source are discussed. Arc distance along the array is used as the independent parameter so that all solutions are constrained to be space curves. Examples of application to real data are presented, and techniques to discriminate against bad sensor data are discussed     

MUSIC算法在矢量水听器阵中的应用研究

文中推导了阵元为矢量水听器时的MU S IC算法。对直线阵情况进行了计算机仿真,结果表明,矢量阵MU-S IC算法与声压阵相比具有更高的方位分辨力。     

Modal time-series structure in a shallow-water environment [HudsonCanyon region]

The broad-band acoustic characterization of the Hudson Canyon region off the New Jersey Continental Shelf is studied with an analysis of pressure time series generated by small explosive sources and recorded on a vertical line array (VLA). The average water depth is about 72 m and the average sound-speed profile (SSP) is downward-refracting over the midportions of the water column. The seabed is characterized by sediment layers possessing sand-like characteristics. The sound-speed structure of the water column and the seabed structure create distinguishing modal features in the impulse response in the 250-500-Hz hand. The details of the depth and range dependence of the time series on the VLA are sensitive to small perturbations of the structure of the upper layer of the SSP, the water depth, and the seabed structure. This sensitivity of the acoustic field is investigated using a broad-band range-dependent normal mode model called NAUTILUS. The representation of the spatial and temporal structure of the time series in terms of a modal structure reveals several unique effects of the SSP and the geoacoustic structure of the bottom on the group velocity of the modes over a large bandwidth. Individual modes can be identified in the measured data using direct data-simulation comparisons. Cross-correlation values between data and simulations in a 155-ms time window generally vary from 0.7 to 0.9 for sensors below the thermocline but are much smaller for sensors above the thermocline     

Target localization and tracking from Doppler-shift measurements

It is possible to localize a nonmaneuvering moving source radiating a constant frequency tone from measurements of the Doppler-shifted frequency at several sensors. Due to the nonlinear nature of the problem, it is necessary to find the solution by grid searches. However, if measurements of the rates of frequency changes are available, the search is only in three dimensions instead of the normal five in source frequency, its x-y positions, and speeds. The validity of combining frequency and frequency-rate measurements is confirmed with simulation studies. The overall system includes a least-squares track-sort algorithm to differentiate the true track from the extraneous track, and a Kalman tracker for the prediction of future source positions, thereby reducing the grid search size. An error analysis relating localization accuracy to uncertainties in frequency measurements and sensor positions is also given     

Absolute positioning of an autonomous underwater vehicle using GPS and acoustic measurements

Kinematic global positioning system (GPS) positioning and underwater acoustic ranging can combine to locate an autonomous underwater vehicle (AUV) with an accuracy of /spl plusmn/30cm (2-/spl sigma/) in the global International Terrestrial Reference Frame 2000 (ITRF2000). An array of three precision transponders, separated by approximately 700 m, was established on the seafloor in 300-m-deep waters off San Diego. Each transponder's horizontal position was determined with an accuracy of /spl plusmn/8 cm (2-/spl sigma/) by measuring two-way travel times with microsecond resolution between transponders and a shipboard transducer, positioned to /spl plusmn/10 cm (2-/spl sigma/) in ITRF2000 coordinates with GPS, as the ship circled each seafloor unit. Travel times measured from AUV to ship and from AUV to transponders to ship were differenced and combined with AUV depth from a pressure gauge to estimate ITRF2000 positions of the AUV to /spl plusmn/1 m (2-/spl sigma/). Simulations show that /spl plusmn/30 cm (2-/spl sigma/) absolute positioning of the AUV can be realized by replacing the time-difference approach with directly measured two-way travel times between AUV and seafloor transponders. Submeter absolute positioning of underwater vehicles in water depths up to several thousand meters is practical. The limiting factor is knowledge of near-surface sound speed which degrades the precision to which transponders can be located in the ITRF2000 frame.     

Analysis of a possible sea floor strain measurement system

Abstract

A logical approach to large area sea floor strain measurement is to use a set of precision acoustic transponders interrogated successively from a large number of different positions by a near‐bottom survey vehicle. Limiting errors in such an approach can be segregated into two classes implying two different scales on which averaging should be carried out. First are those arising from individual travel time and vehicle depth measurements. Second are those associated with imperfect knowledge of the sound propagation speed. The first are essentially independent from one measurement to the next; for the second, statistically independent observations must be separated by the order of 100 m. Several thousand observations of successive travel times to individual transponders would thus be smoothed to produce the equivalent of a smaller number (few hundred) of sets of simultaneous range observations, and then these sets would be used to determine the transponder array geometry. Computer simulations using realistic assumptions show that centimeter‐level accuracy can be achieved over areas several kilometers across.     

A Portable Matched-Field Processing System Using Passive Acoustic Time Synchronization

A portable matched-field processing (MFP) system for tracking marine mammals is presented, constructed by attaching a set of autonomous flash-memory acoustic recorders to a rope to form a four-element vertical array, or "insta-array." The acoustic data are initially time-synchronized by performing a matched-field global inversion using acoustic data from an opportunistic source, and then by exploiting the spatial coherence of the ocean ambient noise background to measure and correct for the relative clock drift between the autonomous recorders. The technique is illustrated by using humpback whale song collected off the eastern Australian coast to synchronize the array, which is then used to track the dive profile of the whale using MFP methods. The ability to deploy autonomous instruments into arbitrary "insta-array" geometries with conventional fishing gear may permit nonintrusive array measurements in regions currently too isolated, expensive, or environmentally hostile for standard acoustic equipment     

Acoustic tracking of a freely drifting sonobuoy field: experimental results

This paper describes a regularized acoustic inversion algorithm for tracking individual elements of a freely drifting sonobuoy field using measured acoustic arrival times from a series of impulsive sources. The acoustic experiment involved 11 sonobuoys distributed over an 8/spl times/6-km field, with a total of six sources deployed over 72 min. The inversion solves for an independent track for each sonobuoy (parameterized by the sonobuoy positions at the time of each source transmission), as well as for the source positions and transmission instants. Although this is a strongly under-determined problem, meaningful solutions are obtained by incorporating a priori information consisting of prior estimates (with uncertainties) for the source positions and initial sonobuoy positions and a physical model for sonobuoy motion along preferentially smooth tracks. The inversion results indicate that the sonobuoys move approximately 260-700 m during the source-deployment period. Closely spaced sonobuoys move along similar tracks; however, there is considerable variability in track directions over the entire field. Positioning uncertainties in horizontal coordinates are estimated using a Monte Carlo appraisal procedure to be approximately 100 m in an absolute sense and 65 m in a relative sense. A sensitivity study indicates that the uncertainties of the a priori position estimates are the limiting factor for track accuracy, rather than data uncertainties or source configuration.     

A hybrid simplex genetic algorithm for estimating geoacousticparameters using matched-field inversion

Matched-fieId inversion (MFI) undertakes to estimate the geometric and geoacoustic parameters in an ocean acoustic scenario by matching acoustic field data recorded at hydrophone array with numerical calculations of the field. The model which provides the best fit to the data is the estimate of the actual experimental scenario. MFI provides a comparatively inexpensive method for estimating ocean bottom parameters over an extensive area. The basic components of the inversion process are a sound propagation model and matching (minimization) algorithm. Since a typical MFI problem requires a large number of computationally intensive sound propagation calculations, both of these components have to be efficient. In this study, a hybrid inversion algorithm which uses a parabolic equation propagation model and combines the downhill simplex algorithm with genetic algorithms is introduced. The algorithm is demonstrated on synthetic range-dependent shallow-water data generated using the parabolic equation propagation model. The performance for estimating the model parameters is compared for realistic signal-to-noise ratios in the synthetic data