Mapping climatic temperature changes in the ocean with acoustictomography: navigational requirements

In eddy-resolving hydrodynamic models, first-mode baroclinic Rossby waves linked to El Nino/Southern Oscillation are the dominant features which change basin-wide temperatures below the seasonal thermocline in the northeast Pacific at periods less than a decade. Simulations are carried out in which Rossby waves are mapped using acoustic tomography. Based on the model which propagated these waves, a Kalman filter is used to map temperature signals for a year. The modeled data are taken from a dense network of acoustic tomography sections. At 300-m depth, where the temperature perturbations associated with Rossby waves are about ±1°C, 80% to 90% of the model variance is accounted for with tomographic estimates. The corresponding standard deviations of the estimates are less than 0.1°C at 400-km resolution. About 80% of the model variance is accounted for with tomography when the navigational errors of the sources and receivers are as poor as one kilometer. Consequently, it may be unnecessary to accurately navigate actual tomographic instruments to map climate change. Modeling results are insensitive to: 1) a reduction in data due to a significant number of instruments which fail; 2) whether the instruments are mobile or fixed; 3) the detailed trajectories of mobile receivers; 4) the shape of the a priori spectrum of ocean fluctuations; 5) the corrections to the acoustic travel-time biases; and 6) the errors in the sound-speed algorithm. In basin-scale arrays, the modeled variance of acoustic travel time depends on the horizontal wavenumber of temperature as k^-5.5. Because sound has little sensitivity to small wavelengths, modeled Rossby waves can be mapped in a day from a few sources and of order ten receivers. The results only depend on the model having large scales in space and time     

Estimating climatic temperature change in the ocean with syntheticacoustic apertures

An acoustic tomography simulation is carried out in the eastern North Pacific ocean to assess whether climate trends are better detected and mapped with mobile or fixed receivers. In both cases, acoustic signals from two stationary sources are transmitted to ten receivers. Natural variability of the sound-speed field is simulated with the Naval Research Laboratory (NRL) layered-ocean model. A sequential Kalman-Bucy filter is used to estimate the sound speed field, where the a priori error covariance matrix of the parameters is estimated from the NRL model. A spatially homogeneous climate trend is added to the NRL fluctuations of sound speed, but the trend is not parameterized in the Kalman filter. Acoustic travel times are computed between the sources and receivers by combining sound speeds from the NRL model with those from the unparameterized climate trend. The effects of the unparameterized climate trend are projected onto parameters which eventually drift beyond acceptable limits. At that time, the unparameterized trend is detected. Mobile and fixed receivers detect the trend at about the same time. At detection time, however, maps from fixed receivers are less accurate because some of the unparameterized climate trend is projected onto tile spatially varying harmonics of the sound-speed field. With mobile receivers, the synthetic apertures suppress the projection onto these harmonics. Instead, the unparametrized trend is correctly projected onto the spatially homogeneous portion of the parameterized sound-speed field     

Refraction of Sound in a Horizontally Inhomogeneous, Time-Dependent Ocean

Measurements of the three-dimensional (3-D) structure of a sound-speed field in the ocean with the spatial and temporal resolution required for prediction of acoustic fields are extremely demanding in terms of experimental assets, and they are rarely available in practice. In this study, a simple analytic technique is developed within the ray approximation to quantify the uncertainty in acoustic travel time and propagation direction that results from an incomplete knowledge or purely statistical characterization of sound-speed variability in the horizontal plane. Variation of frequency of an acoustic wave emitted by a narrowband source due to a temporal variation of environmental parameters is considered for deterministic and random media. In a random medium with locally statistically homogeneous, time-dependent 3-D fluctuations of the sound speed, calculation of the signal frequency and bearing angle variances as well as the travel-time bias due to horizontal refraction is approximately reduced to integration of respective statistical parameters of the environmental fluctuations along a ray in a background, range-dependent, deterministic medium. The technique is applied to acoustic transmissions in a coastal ocean, where tidally generated nonlinear internal waves are the prevailing source of sound-speed fluctuations, and in a deep ocean, where the fluctuations are primarily due to spatially diffuse internal waves with the Garrett–Munk spectrum. The significance of 3-D and four-dimensional (4-D) acoustic effects in deep and shallow water is discussed.     

Optimum design of the ocean acoustic tomography system for the Sea of Japan

An ocean acoustic tomography system covering the region of 800×1000 km with the spatial resolution of eddy-resolving scales has been designed on the basis of computer experiments using the hydrographic data collected in the Sea of Japan. The optimum number of acoustic sources required for 20 acoustic receivers was estimated as 13 by changing the source number. The spatial resolution for the optimum system was 41 km smaller than the dominant size of meso-scale eddies in the Sea of Japan. The effect of travel-time errors on tomographic maps is also quantified.     

Multimegameter-range acoustic data obtained by bottom-mountedhydrophone arrays for measurement of ocean temperature

Acoustic signals transmitted from the ATOC source on Pioneer Seamount off the coast of California have been received at various sites around the Pacific Basin since January 1996. We describe data obtained using bottom-mounted receivers, including US Navy Sound Surveillance System arrays, at ranges up to 5 Mm from the Pioneer Seamount source. Stable identifiable ray arrivals are observed in several cases, but some receiving arrays are not well suited to detecting the direct ray arrivals. At 5-Mm range, travel-time variations at tidal frequencies (about 50 ms peak to peak) agree well with predicted values, providing verification of the acoustic measurements as well as the tidal model. On the longest and northernmost acoustic paths, the time series of resolved ray travel times show an annual cycle peak-to-peak variation of about 1 s and other fluctuations caused by natural oceanic variability. An annual cycle is not evident in travel times from shorter acoustic paths in the eastern Pacific, though only one realization of the annual cycle is available. The low-pass-filtered travel times are estimated to an accuracy of about 10 ms. This travel-time uncertainty corresponds to errors in range- and depth-averaged temperature of only a few millidegrees, while the annual peak-to-peak variation in temperature averaged horizontally over the acoustic path and vertically over the upper 1 km of ocean is up to 0.5°C     

The coupled oceanographic–tomographic analysis and prediction system: A practical means for utilizing high frequency tomography in the world's oceans

Details are presented of a methodology that utilizes acoustic travel time information in an ocean circulation model. Recent developments of this model-oriented tomography are discussed, representing some significant improvements over earlier formulations. More accurate means of determining the arrival times of specific ray paths are detailed, along with a means of estimating possible errors in the calculated travel times. The assimilation of the observed arrival time information into an ocean model is achieved using a Kalman gain, and more advanced expressions for calculating the Kalman gain are presented. A formulation to account for errors in the stated positions of a source and receiver is also presented. It is shown that the methodology performs fairly well in reproducing observed travel time anomalies. However, the model-predicted anomalies along a specific ray path may not always track the observed anomalies for that path when assimilating multiple ray path data. Results indicate that additional work is required to determine a means of handling observed arrival time data without having prior knowledge of the magnitude of errors in the observations. Results from simulation experiments provide estimates of: (1) potential errors when the travel times for ray paths are only sampled at discreet intervals as opposed to continuously and (2) to what degree acoustic data can be expected to “correct” model-predicted fields.     

Comparison of Kalman and least squares filters for locatingautonomous very low frequency acoustic sensors

The Marine Physical Laboratory has designed, fabricated, and taken to sea self-contained, freely drifting acoustic sensors which can measure signal propagation and ambient ocean noise in the 1-20-Hz band for up to 25-hour periods. The deployment of several freely drifting floats forms an array of sensors whose outputs can be combined after the experiment with a beamformer. A Kalman filter and a least-squares estimator have been developed to estimate float positions from travel-time measurements. Computer simulation is used to compare filter performance-under several deployment scenarios. Results show that the Kalman filter performs better than the least-squares filter when the floats are subjected to small-magnitude accelerations between measurements. Neither filter was sensitive to relatively major changes in deployment geometry as long as the sound-speed profile is known exactly     

海洋声层析观测技术和方法

叙述海洋声层析观测系统,以声线传播时间层析为重点概括了海洋声层析的基本原理和其他主要方法,共6个方法。对运用海洋声层析观测来反演海洋状态问题的建立、求解及其误差来源作了分析和讨论。以测量声线传播时间为例介绍了海洋声层析观测系统主要设计技术。     

Reciprocal acoustic transmissions: Instrumentation for Mesoscale monitoring of ocean currents

By simultaneously transmitting acoustic pulses in opposite directions between two points in midocean, one can separate the effects of ocean currents on acoustic propagation from the effects of sound-speed structure. Reciprocal acoustic transmissions can therefore be used to measure ocean currents. Acoustic transceivers have been designed and built to measure the mean currents between two points separated by 300 km. The equipment functioned satisfactorily during a sbort test conducted during 1983. Preliminary analysis of that experiment has yielded differential travel times that appear reasonable, but more work is required to relate the differential travel times to meaningful ocean-current estimates.     

Consideration of fine-scale coastal oceanography and 3-D acoustics effects for the ESME sound exposure model

Results and recommendations for evaluating the effects of fine-scale oceanographic scattering and three-dimensional (3-D) acoustic propagation variability on the Effects of Sound on the Marine Environment (ESME) acoustic exposure model are presented. Pertinent acoustic scattering theory is briefly reviewed and ocean sound-speed fluctuation models are discussed. Particular attention is given to the nonlinear and linear components of the ocean internal wave field as a source of sound-speed inhomogeneities. Sound scattering through the mainly isotropic linear internal wave field is presented and new results relating to acoustic scattering by the nonlinear internal wave field in both along and across internal wave wavefront orientations are examined. In many cases, there are noteworthy fine-scale induced intensity biases and fluctuations of order 5-20 dB.     

Temporal resolutions of time-reversal and passive-phase conjugation for underwater acoustic communications

In this paper, we study the temporal resolution of a time-reversal or passive-phase conjugation process as applied to underwater acoustic communications. Specifically, we address 1) the time resolution or the pulse width of a back-propagated time-compressed pulse as compared with the original transmitted pulse; 2) the effectiveness of temporal focusing as measured by the peak-to-sidelobe ratio of the back-propagated or phase-conjugated pulse (both pulse elongation and sidelobe leakages are causes of intersymbol interference and bit errors for communications); 3) the duration of temporal focusing or the temporal coherence time of the underwater acoustic channel; and 4) the stability of temporal focusing as measured by the phase fluctuations of successive pulses (symbols). Binary phase-shift keying signals collected at sea from a fixed source to a fixed receiver are used to extract the above four parameters and are compared with simulated results. Mid-frequency (3-4-kHz) data were collected in a dynamic shallow-water environment, exhibiting high temporal fluctuations over a scale of minutes. Despite this, the channel is found to be highly coherent over a length of 17 s. As a result, only one probe signal is used for 17 s of data. The bit error rate and variance of the symbol phase fluctuations are measured as a function of the number of receivers. They are of the same order as that calculated from the simulated data. The agreement suggests that these two quantities could be modeled for a communication channel with high coherence time. The phase variance can be used to determine the maximum data rate for a phase-shift keying signal for a given signal bandwidth and a given number of receivers.     

Ocean acoustic tomography: estimating the acoustic travel time withphase

Continuous acoustic transmission (133 Hz, 60-ms resolution) between a bottom-mounted source near Oahu, Hawaii, and a bottom-mounted receiver at 4000-km range near the coast of northern California was recorded to learn how to measure precisely the travel time so that basin-scale fluctuations in the Pacific can be detected. Daily incoherent averages of some of the multipaths exhibited stability during this period. The standard deviation of the travel time of the resolved peaks in the daily incoherent averages is about 30 ms. An acoustic method, based on cross-correlation, is derived to estimate the change in the average acoustic phase (travel time) to a precision of about 0.018 cycles (135 μs) every 2 min. Travel-time estimates based on the cross-correlator reduce the aberrations due to internal waves by about 19 dB in comparison with CW transmissions. The new travel-time estimator is applied to the measurements to examine some of the fluctuations of the Pacific     

Differences between passive-phase conjugation and decision-feedback equalizer for underwater acoustic communications

Passive-phase conjugation (PPC) uses passive time reversal to remove intersymbol interferences (ISIs) for acoustic communications in a multipath environment. It is based on the theory of signal propagation in a waveguide, which says that the Green's function (or the impulse-response function) convolved with its time-reversed conjugate, summed over a (large-aperture) vertical array of receivers (denoted as the Q function) is approximately a delta function in space and time. A decision feedback equalizer (DFE) uses a nonlinear filter to remove ISI based on the minimum mean-square errors (mmse) between the estimated symbols and the true (or decision) symbols. These two approaches are motivated by different principles. In this paper, we analyze both using a common framework. We note the commonality and differences, and pros and cons, between the two methods and compare their performance in realistic ocean environments, using simulated and at-sea data. The performance measures are mean-square error (mse), output signal-to-noise ratio (SNR), and bit-error rate (BER) as a function of the number of receivers. For a small number of receivers, the DFE outperforms PPC in all measures. The reason for poor PPC performance is that, for a small number of receivers, the Q function has nonnegligible sidelobes, resulting in nonzero ISI. As the number of receivers increases, the BER for both processors approaches zero, but at a different rate. The modeled performance differences (in mse and SNR) between PPC and DFE are in general agreement with the measured values from at-sea data, providing a basis for performance prediction.     

基于最小二乘和卡尔曼滤波方法进行原子时预报的研究

介绍和分析了最小二乘和卡尔曼滤波方法在时间预报中的应用。通过IGS站提供的钟差数据,分别运用这两种方法对其中六个原子钟进行了时间预报的实验。通过对预报结果进行分析,结论证明为了取得较好的预报效果,不同的预报方法和钟参数的模型对于观测数据的要求有一些差别。从实验数据所成图形来看,当采用一天的观测数据进行模型预报时,最小二乘法的预报精度比卡尔曼滤波法稍高一些。     

Inversion for Time-Evolving Sound-Speed Field in a Shallow Ocean by Ensemble Kalman Filtering

In the context of the recent Maritime Rapid Environmental Assessment/Blue Planet 2007 sea trial (MREA/BP07), this paper presents a range-resolving tomography method based on ensemble Kalman filtering of full-field acoustic measurements, dedicated to the monitoring of environmental parameters in coastal waters. The inverse problem is formulated in a state–space form wherein the time-varying sound-speed field (SSF) is assumed to follow a random walk with known statistics and the acoustic measurements are a nonlinear function of the SSF and the bottom properties. The state–space form enables a straightforward implementation of a nonlinear Kalman filter, leading to a data assimilation problem. Surface measurements augment the measurement vector to constrain the range-dependent structure of the SSF. Realistic scenarios of vertical slice shallow-water tomography experiments are simulated with an oceanic model, based on the MREA/BP07 experiment. Prior geoacoustic inversion on the same location gives the bottom acoustic properties that are input to the propagation model. Simulation results show that the proposed scheme enables the continuous tracking of the range-dependent SSF parameters and their associated uncertainties assimilating new measurements each hour. It is shown that ensemble methods are required to properly manage the nonlinearity of the model. The problem of the sensitivity to the vertical array (VA) configuration is also addressed.      

海洋声层析的基本原理和应用

从射线声学和简正波声学的角度，概述了海洋声层析的基本理论，包括射线走时反演、简正波走时反演、简正波相位反演和简正波水平折射层析。海洋声层析以反演海水温度和流速为基础。还总结了声层析在海洋学研究中的应用。     

Shallow-Water Acoustic Tomography Performed From a Double-Beamforming Algorithm: Simulation Results

Recent shallow-water experiments in sea channels have been performed using two vertical coplanar densely sampled source and receive arrays. Applying a double-beamforming algorithm on the two arrays both on synthetic numerical simulations and on experimental data sets, we extract efficiently source and receive angles as well as travel times for a large number of acoustic rays that propagate and bounce in the shallow-water waveguide. We then investigate how well sound-speed variations in the waveguide are reconstructed using a ray time-delay tomography based on a Bayesian inversion formulation. We introduce both data and model covariance matrices and we discuss on the synthetic numerical example how to choose the a priori information on the sound-speed covariance matrix. We attribute the partial sound-speed reconstruction to the ray-based tomography and we suggest that finite-frequency effects should be considered as the vertical and horizontal size of the Fresnel zone significantly spreads in the waveguide. Finally, the contribution of a different set of ray angles for tomography goal is also presented.      

海洋资料同化对气候季节-年际预测技巧及初始场的影响试验

集合卡尔曼滤波(Ensemble Kalman filter, EnKF)是一种国内外广泛使用的海洋资料同化方案, 用集合成员的状态集合表征模式的背景误差协方差, 结合观测误差协方差, 计算卡尔曼增益矩阵, 有效地将观测信息添加到模式初始场中。由于季节、年际预测很大程度上受到初始场的影响, 因此资料同化可以提高模式的预测性能。本文在NUIST-CFS1.0预测系统逐日SST nudging的初始化方案上, 利用EnKF在每个月末将全场(full field)海表温度(sea surface temperature, SST)、温盐廓线(in-situ temperature and salinity profiles, T-S profiles)以及卫星观测海平面高度异常(sea level anomalies, SLA)观测资料同化到模式初始场中, 对比分析了无海洋资料同化以及加入同化后初始场的区别、加入海洋资料同化后模式提前1~24个月预测性能的差异以及对于厄尔尼诺-南方涛动(El Niño-southern oscillation, ENSO)预测技巧的影响。结果表明, 加入海洋资料同化能有效地改进初始场, 并且呈现随深度增加初始场改进越显著的特征。加入同化后, 对全球SST、次表层海水温度的平均预测技巧均有一定的提高, 也表现出随深度增加预测技巧改进越明显的特征。但加入海洋资料同化后, 模式对ENSO的预测技巧有所下降, 可能是由于模式误差的存在, 使得同化后的预测初始场从接近观测的状态又逐渐恢复到与模式动力相匹配的状态, 加剧了赤道太平洋冷舌偏西、中东部偏暖的气候平均态漂移。     

Multiple-source localization using GPS technology and receivedarrival time structure analysis in an air-deployed system

An efficient and robust method has been developed to locate multiple impulsive sources in an ocean environment. Global position system (GPS) receivers were installed on sonobuoys to obtain their locations within a few meters of accuracy. A sonobuoy field was deployed in a ring-type pattern. Charges were then set off at arbitrary locations within the ring, High-resolution plots were used to obtain direct path and/or first bottom bounce arrivals on each buoy. A model grid of arrival times was constructed, corresponding to the dimensions of the buoy field. A ray model previously developed here at the Applied Research Laboratories at the University of Texas at Austin (ARL:UT) was used to obtain model travel times. The minimum value of the least-square-type error between the real arrival times and the modeled travel times resulted in an unambiguous location of the source, within the limits of the grid spacing chosen. This value was calculated by picking one receiver as the reference and then summing the timing errors of the remaining receivers relative to the reference. Successive iterations with finer grid spacings result in source localization within the accuracy of the buoy locations. The localization routine was extended by allowing permutations of the pulse arrivals on each buoy to account for multiple sources closely separated in time and/or space. An automated correlation technique is presented as an alternative to the leading edge-detection method used here for obtaining relative arrival times. Two proof-of-concept experiments were performed and some results of data obtained at Lake Travis and the Gulf of Mexico are presented     

Coastal acoustic tomography system and its field application

The coastal acoustic tomography system (CATS), composed of five moored acoustic stations, has been constructed to measure current fields. The system is developed with special considerations in mind, including the use of Global Positioning System clock signals in the synchronization of the system clock timing among the multiple acoustic stations, and the use of the differently coded Gold sequences to identify the acoustic signals corresponding to individual stations from a received signal. The CATS was successfully applied to map the structure of strongly nonlinear tidal currents in the coastal sea. In spite of the limited spatial resolution caused by inadequate sound transmission data, the two-dimensional tidal vortices features of growth, translation, and decay processes are reconstructed through an inverse analysis of the acoustic travel time obtained among the station pairs. It is evident that the CATS is a powerful tool for measuring variable current fields generated in the coastal seas