statistics of broadband transmissions through a range-dependent fluctuating ocean waveguide mark...

1
Statistics of broadband transmissions through a range-dependent fluctuating ocean waveguide Mark Andrews and Purnima Ratilal; Northeastern University, Boston, MA 02115 Source and Tracks Site 1 Site 2 o A SL 10 log 10 ) 2 exp( ) , | ( 2 r f r r G A I o o Likelihood Function: N i i i I W W I W P 1 1 ) ( exp ) / ( ) ( range-dependent correction coefficient Range- dependent Expected intensity : Normalized Likelihood Functions 10*log |P(W)| Parseval Sum Matched filter , range dependent correction to TL (dB/km) Source Level (dB) Range (km) Range (km) dB Parseval Sum Matched Filter ) 2 exp( ) | ( 2 r r r G A I o o Source Level : Acoustic signals transmitted through an ocean waveguide are temporally and spatially varying due to multimodal interference and effects such as internal waves. Estimating parameters of sonar, the environment and scatterers requires a statistical approach that incorporates medium uncertainties into the signal analysis. Short duration broadband pulses were transmitted from a moored source array and measured by a towed horizontal receiving array at varying ranges from the source. The ping-to-ping fluctuations in the measured acoustic data are a result of changes in the waveguide modal intereference structure due to both motion of the array and presence of time dependent random internal waves. Our analysis shows that the Parseval sum and match filter outputs of the broadband transmissions for the direct arrival have significantly smaller standard deviations compared to the instantaneous single frequency transmission. We account for these observations by modeling the broadband acoustic field intensity in both a static waveguide and with Monte Carlo simulations in a fluctuating waveguide. A maximum likelihood estimator is implemented to provide a global inversion of the data for source level and a range-dependent expected intensity, as well as quantifying the match filter degradation in the multi-modal ocean waveguide. ONR Main Acoustic Clutter Experiment (2003), Geoclutter Program Long range SONAR such as OAWRS can be used to instantaneously image schools of fish or other targets. Propagation and statistics of broadband signals must be thoroughly understood to invert for scattering strength Center frequency: 415 Hz Bandwidth: 50 Hz Duration: 1.0 sec Broadband transmissions propagated through an ocean waveguide are attenuated and dispersed as they propagate out in range. dB Center frequency 4.98 1.18 Parseval Sum 2.99 2.56 Match Filter 3.64 1.88 m r C 500 hours C 3 Bistatic experimental setup used for Ocean Acoustic Wave Remote Sensing (OAWRS) Broadband Transmissions LFM transmission Range = 2.5km Range = 13km One Way Propagation Data 2 2 ) , | ( ) ( ) ( f r r G f Q f o f i f i f f o f f T o df f r r G f Q df f dt t E 2 2 2 ) , | ( ) ( ) ( ) ( ) ( ) ( * t t Kq t h M ) 2 exp( ) ( ) 2 exp( ) ( ) ( * M ft j f KQ dt ft j t h f H f i f i f f M o f f df t t f j f r r G f Q K df ft j f H f t MF )] ( 2 exp[ ) , | ( ) ( ) 2 exp( ) ( ) ( ) ( 2 df ft j f r r G f Q t o ) 2 exp( ) , | ( ) ( ) ( Received Signal: Parseval Sum: Match Filter: Match Filter Field: Single Frequency Field: Parseval Sum is the sum of all of the energy in the bandwidth of the signal Match filter extracts the energy that “matches” the original signal Q(f) ) , | ( ) ( ) 2 exp( ) ( ) ( f r r G f Q dt ft j t f o Single Frequency Intensity: Match Filter Intensity: 2 2 2 )] ( 2 exp[ ) , | ( ) ( ) ( f i f f M o df t t f j f r r G f Q t MF 1 2 ) ( f i f f df f Q K Broadband Processing Ocean Waveguide Computational Modeling Broadband signals transmitted directly from the source to receiver traveling along a straight track were collected and processed using a Parseval sum intensity, match filter intensity, and compared to the intensity at the center frequency. The statistics using a running window over each track is shown below. Sound Speed (m/s) Parseval Sum Match Filter The Range Dependent Acoustic Model (RAM) accounts for the environmental parameters of the ocean, such as the sea floor geoacoustic parameters, the sound speed profile, and the changing bathymetry over range (above). It is applied here to investigate the statistics of the processed signal intensity in both space and time, using Monte Carlo simulations. (below) Center Frequency intensity Depth (m) Depth (m) Internal wave correlation time and range for Monte Carlo simulations: Sound speed profiles taken througho ut experime nt Standard Deviations Mean s Statistics of data over 20 pings Calibrating for Source Level and Expected Intensity from one-way Data An approach is demonstrated to calibrate for the source level as well as the range-dependent expected intensity using a maximum likelihood estimator. Statistical spatial and temporal fluctuations of the broadband intenisties, as well as match filter degradation are accounted for using a range-dependent computational model. It is demonstrated that modal interference and modal dispersion are primary causes of these effects. Statistics of dynamic model over 20 pings:: The statistics of the one-way propagated data is used to calibrate for the expected intensity versus range by calibrating both the Source level, as well as a range dependent correction, , to the expected transmission loss. This is accomplished using a maximum likelihood estimator, which incorporates all of the one-way data and the expected transmission loss using a depth average of the RAM model. Expected transmission loss using RAM model Wi = Intensity measurement = time bandwidth product (see below) Time (s) Time (s) Time (s) Frequency (Hz) Frequency (Hz) Frequency (Hz) May 14 2003: 09:32:25 EDT May 15, 2003: 08:49:55 EDT Abstract At depth of 50 m Range (km) Range (km) dB dB dB Calibration Results: Parsev al Sum Matche d Filte r SL (dB) 218.3 217.6 (dB/km ) 0.38 0.30 MF degradation = 0.7dB Range-dependent deg = 0.08dB/km Expected intensities and Data Time Bandwidth Product The time bandwidth product, , is a measurement of the number of independent measurement samples. ) , 2 ( ) log 10 ( e I Intensity standard deviation relates to time-bandwidth product: Conclusions Tianrun Chen, Sunwoong Lee, and Nicholas Makris; Massachusetts Institute of Technology, Cambridge, MA 02139 References N. C. Makris, “The effect of saturated transmission scintillation on ocean acoustic intensity measurements,” J. Acoust. Soc. Am. 100(2), 769-783 (1996). Goodman, Joseph W., Statistical Optics, Wiley, New York, (1985) Acknowledgements This work was funded by the Office of Naval Research and the Sloan Foundation

Upload: stephen-woods

Post on 13-Jan-2016

216 views

Category:

Documents


2 download

TRANSCRIPT

Page 1: Statistics of broadband transmissions through a range-dependent fluctuating ocean waveguide Mark Andrews and Purnima Ratilal; Northeastern University,

Statistics of broadband transmissions through a range-dependent fluctuating ocean waveguide

Mark Andrews and Purnima Ratilal; Northeastern University, Boston, MA 02115

Source and Tracks

Site 1

Site 2oASL 10log10

)2exp(),|(2

rfrrGAI oo

Likelihood Function:

N

i

ii I

WWI

WP1

1

)(

exp)/()(

range-dependent correction coefficient

Range-dependent Expected intensity :

Normalized Likelihood Functions 10*log |P(W)|

Parseval Sum Matched filter

, range dependent correction to TL (dB/km)

So

urce

Le

vel (

dB)

Range (km) Range (km)

dB

Parseval Sum Matched Filter

)2exp()|(2

rrrGAI oo

Source Level :

Acoustic signals transmitted through an ocean waveguide are temporally and spatially varying due to multimodal interference and effects such as internal waves. Estimating parameters of sonar, the environment and scatterers requires a statistical approach that incorporates medium uncertainties into the signal analysis. Short duration broadband pulses were transmitted from a moored source array and measured by a towed horizontal receiving array at varying ranges from the source. The ping-to-ping fluctuations in the measured acoustic data are a result of changes in the waveguide modal intereference structure due to both motion of the array and presence of time dependent random internal waves. Our analysis shows that the Parseval sum and match filter outputs of the broadband transmissions for the direct arrival have significantly smaller standard deviations compared to the instantaneous single frequency transmission. We account for these observations by modeling the broadband acoustic field intensity in both a static waveguide and with Monte Carlo simulations in a fluctuating waveguide. A maximum likelihood estimator is implemented to provide a global inversion of the data for source level and a range-dependent expected intensity, as well as quantifying the match filter degradation in the multi-modal ocean waveguide.

ONR Main Acoustic Clutter Experiment (2003), Geoclutter Program

●Long range SONAR such as OAWRS can be used to instantaneously image schools of fish or other targets.●Propagation and statistics of broadband signals must be thoroughly understood to invert for scattering strength

Center frequency: 415 Hz

Bandwidth: 50 Hz

Duration: 1.0 sec

Broadband transmissions propagated through an ocean waveguide are attenuated and dispersed as they propagate out in range.

dB

Center frequency 4.98 1.18

Parseval Sum 2.99 2.56

Match Filter 3.64 1.88

mrC 500

hoursC 3

Bistatic experimental setup used for Ocean Acoustic Wave Remote Sensing (OAWRS)

Broadband Transmissions

LFM transmission Range = 2.5km Range = 13km

One Way Propagation Data

22),|()()( frrGfQf o

f

i

f

i

f

f

o

f

fTo dffrrGfQdffdttE

222),|()()()(

)()( * ttKqth M )2exp()()2exp()()( *MftjfKQdtftjthfH

f

i

f

i

f

f

Mo

f

f

dfttfjfrrGfQKdfftjfHftMF )](2exp[),|()()2exp()()()(2

dfftjfrrGfQt o )2exp(),|()()(

Received Signal:

Parseval Sum:

Match Filter:

Match Filter Field:

Single Frequency Field:

Parseval Sum is the sum of all of the energy in the bandwidth of the signal

Match filter extracts the energy that “matches” the original signal Q(f)

),|()()2exp()()( frrGfQdtftjtf o

Single Frequency Intensity:

Match Filter Intensity:

2

22)](2exp[),|()()(

f

i

f

f

Mo dfttfjfrrGfQtMF 1

2)(

f

i

f

f

dffQK

Broadband Processing

Ocean Waveguide Computational Modeling

Broadband signals transmitted directly from the source to receiver traveling along a straight track were collected and processed using a Parseval sum intensity, match filter intensity, and compared to the intensity at the center frequency. The statistics using a running window over each track is shown below.

Sound Speed (m/s)Parseval Sum Match Filter

The Range Dependent Acoustic Model (RAM) accounts for the environmental parameters of the ocean, such as the sea floor geoacoustic parameters, the sound speed profile, and the changing bathymetry over range (above). It is applied here to investigate the statistics of the processed signal intensity in both space and time, using Monte Carlo simulations. (below)

Center Frequency intensity

Dep

th (

m)

Dep

th (

m)

Internal wave correlation time and range for Monte Carlo simulations:

Sound speed profiles taken throughout experiment

Standard DeviationsMeans

Sta

tistic

s of

dat

a ov

er 2

0 pi

ngs

Calibrating for Source Level and Expected Intensity from one-way Data

An approach is demonstrated to calibrate for the source level as well as the range-dependent expected intensity using a maximum likelihood estimator.Statistical spatial and temporal fluctuations of the broadband intenisties, as well as match filter degradation are accounted for using a range-dependent computational model. It is demonstrated that modal interference and modal dispersion are primary causes of these effects.

Sta

tistic

s of

dyn

amic

m

odel

ove

r 20

pin

gs::

The statistics of the one-way propagated data is used to calibrate for the expected intensity versus range by calibrating both the Source level, as well as a range dependent correction, , to the expected transmission loss. This is accomplished using a maximum likelihood estimator, which incorporates all of the one-way data and the expected transmission loss using a depth average of the RAM model.

Expected transmission loss using RAM model

Wi = Intensity measurement = time bandwidth product (see below)

Time (s) Time (s)Time (s)

Frequency (Hz) Frequency (Hz)Frequency (Hz)

May 14 2003: 09:32:25 EDT May 15, 2003: 08:49:55 EDT

Abstract

At depth of 50 m

Range (km) Range (km)

dBdB

dB

Calibration Results:

Parseval Sum

Matched Filter

SL (dB) 218.3 217.6

(dB/km) 0.38 0.30

MF degradation = 0.7dB

Range-dependent deg = 0.08dB/km

Expected intensities and Data

Time Bandwidth Product

The time bandwidth product, , is a measurement of the number of independent measurement samples.

),2()log10( eI Intensity standard deviation relates to time-bandwidth product:

Conclusions

Tianrun Chen, Sunwoong Lee, and Nicholas Makris; Massachusetts Institute of Technology, Cambridge, MA 02139

ReferencesN. C. Makris, “The effect of saturated transmission scintillation on ocean acoustic intensity measurements,” J. Acoust. Soc. Am. 100(2), 769-783 (1996).

Goodman, Joseph W., Statistical Optics, Wiley, New York, (1985)

Acknowledgements

This work was funded by the Office of Naval Research and the Sloan Foundation