News Column

A Bistatic Phased-Array Doppler Sonar for Wave-Current Research

July 1, 2014

Smith, Jerome A


Wave breaking and wave-forced flows are important to air-sea interactions and to the transport and dispersal of materials at sea. But recent measurements have shown a discrepancy in the Eulerian response to wave groups compared to scientists' current theoretical understanding of wave-current interactions. Flow structures on scales of centimeters to meters occur underneath breaking waves, and larger-scale flows are driven by wave-current interactions (e.g., Langmuir circulation, alongshore flows). Such details of the vertically resolved flow are just beginning to be modeled, and observational guidance is needed. Here a new instrument is described that is intended to measure waves and currents over a 2D vertical plane underwater, resolving two components of velocity on this plane. Initial observations were made near the Scripps Pier (La Jolla, California), where steep waves and strong currents can be reliably found, yet logistics are not too burdensome. To get the spatial resolution desired using 200-kHz sound, ping-to-ping ''coherent processing'' would have be used for Doppler estimation; however, near shore the reverberations remain too strong for far too long to get any coherence, unlike the previous experience in deep water. In view of this, using much higher frequencies (>1 MHz) with ''incoherent processing'' is suggested; the increased attenuation at higher frequencies then would subdue the reverberation problem, but with comparable space-time resolution.

(ProQuest: ... denotes formulae omitted.)

1. Introduction

Surface waves are central to air-sea interactions (e.g., Thorpe 1982; Edson and Fairall 1994; Pattison and Belcher 1999; Veron et al. 2008; Fan et al. 2009). For example, breaking waves directly inject momentum, gas, and turbulent energy into the mixed layer of the sea and have similar dramatic effects in the air. Although the implications of the intermittent stress transfer due to breaking have historically been the subject of much discussion, they have only recently been simulated and studied in detail (Sullivan et al. 2004, 2007). These simulations indicate that vortical structures resulting from breaking events influence the motion to a depth several times greater than the wave amplitude, deeper than previously thought.

Waves also transport mass and momentum (Stokes 1847; Longuet-Higgins 1953). Longuet-Higgins and Stewart (1962, 1964, p. 530) identified and evaluated the ''excess flow of momentum due to the presence of the waves'' and, in analogy to optics, named it the ''radia- tion stress.'' Changes in the radiation stress (momen- tum flux) of the waves are compensated for by changes in the mean field, so that overall momentum is con- served. For example, Longuet-Higgins and Stewart (1960, 1961) described the generation of group-bound forced long waves. Field measurements adequate for comparison were not obtained until decades later (Herbers and Guza 1991; Herbers et al. 1994; Smith 2006a). In this last case, a discrepancy was identified: the measured response at the surface was found to be more than twice what theory suggests, as discussed further below. For more details than given below, see Smith (2006a);foradiscussionofthewavegroupsre- sponsible, see Smith and Brulefert (2010).

Wave-current and wave-topography interactions are important in a variety of scenarios, including wave- induced setdown and setup near shore (Longuet-Higgins and Stewart 1964; Bowen et al. 1968); generation of long- shore currents (Bowen 1969; Longuet-Higgins 1970a,b); the interaction of freely propagating long and short sur- face waves (Longuet-Higgins 1969; Hasselmann 1971; Garrett and Smith 1976; and many others); and the gen- eration of Langmuir circulation (LC), a prominent form of motion in the wind-driven surface mixed layer (e.g., Langmuir1938; Craik1977;Leibovich 1980; Lietal. 1995; Skyllingstad and Denbo 1995; McWilliams et al. 1997; McWilliams and Sullivan 2000; Phillips 2002; Sullivan et al. 2007; Kukulka et al. 2009). A key term in the LC case involves the bending of vortex lines by the vertically varying Stokes drift of the waves, called the ''vortex force'' by Craik and Leibovich (1976).Initsvertically integrated form, this term corresponds to a mean mo- mentum change to compensate for the change in wave momentum as the waves are refracted (Garrett 1976; Leibovich 1980; Smith 1980), which is a part of the radi- ation stress (Smith 2006b), identified first (but not used) by Hasselmann (1971). Much of the work to date on wave-current interactions still depends on depth- integrated balances, as used in the original radiation stress formulation of the interaction. For depth-resolving simulations, the depth-dependent variations of the Stokes drift, including direction as well as magnitude, are needed to evaluate this vortex force and to properly distribute the remaining radiation stress terms over depth.

As mentioned above, the theory as developed to date may not be complete. Referring to comparisons of the predicted and measured surface current response to wave groups, Smith (2006a, p. 1382) states, ''It is found that Eulerian counterflows occur that cancel the esti- mated Stokes drift variations at the surface completely. This is a stronger surface response than expected from the group-bound forced wave analysis'' (mentioned above). Also, similar sized Eulerian counterflows have been observed in laboratory experiments (Kemp and Simons 1982, 1983; Jiang and Street 1991; Nepf 1992; Groeneweg and Klopman 1998; Swan et al. 2001; Monismith et al. 2007). ''The changes in the Eulerian current profile from case to case (no waves, down-current waves, up-current waves) are comparable to (and oppose) the Stokes drift profile calculated from the wave param- eters as given (in both magnitude and shape)'' (Smith 2006a, p. 1382). While the Coriolis force would eventually cause such depth-by-depth cancellation (McWilliams and Restrepo 1999), in these experiments the cancellation is observed immediately, over times small compared to the passing of a group, far too fast for vorticity to diffuse downward from the surface. In contrast, pressure fields measured in intermediate depth water have been found consistent with second-order theory (Herbers and Guza 1991; Herbers et al. 1994). However, the expected re- sponse at intermediate depths differs markedly from that in deep water, and it is not clear how pressure measure- ments made close to the bottom should correspond to currents throughout the water column.

Measured velocities in both the laboratory and field show persistent differences from theory that indicate the need for further work and understanding, especially re- garding the vertical structure of the forcing and response. This is particularly vital and timely because some ob- servational guidance is needed to aid development of models that resolve the vertical structure of the flows, which is especially important near shore, where wave- current interactions can be the dominant driving force (e.g., Groeneweg and Battjes 2003; Mellor 2003; Özkan- Haller and Li 2003; McWilliams et al. 2004; Smith 2006b; Newberger and Allen 2007; Ardhuin et al. 2008; Kumar et al. 2012). The objective here is to obtain measurements of at least two velocity components over a vertical plane in order to examine the time/space characteristics of both the wave motions and the resulting response as wave groups pass and, in particular, to explore the vertical structure of the response.

Another great challenge in the study of the air-sea interface is to measure and understand natural breaking waves and the resulting subsurface flow structures. It is exceedingly hard to get quantitative measurements closer than a few wave amplitudes to the surface with actively breaking waves, although many promising techniques are being developed (e.g., Vagle and Farmer 1992; Lamarre and Melville 1994b; Farmer et al. 1998; Veron and Melville 1999; Terrill and Melville 2000; Melville et al. 2002; Wu and Nepf 2002; Gemmrich and Farmer 2004; Gemmrich et al. 2008). The objective here is to char- acterize the subsurface flow structures under naturally occurring breaking waves. Emphasis is not on the small- scale turbulent kinetic energy (TKE) or dissipation (cf. Agrawal et al. 1992; Terray et al. 1996; Gemmrich and Farmer 2004), but on the slightly larger-scale vortical structures (cf. Melville et al. 2002, p. 221). As noted there, ''The dominant feature of the post-breaking ve- locity field is the coherent eddy that drifts slowly down- stream at the same speed as the patch of smaller-scale turbulence.'' It has also been shown that breaker shape has an effect on the flow and resulting vortex (cf. Nepf et al. 1998). These ''coherent structures'' should be readily measurable by the proposed technique, although we cannot expect the same level of detail as is possible in the laboratory. Observations to date can help broadly characterize the scales and frequency of breaking waves (Gemmrich and Farmer 1999; Melville and Matusov 2002), and laboratory experiments can reproduce the flow structures under a few select breaker types (cf. Melville et al. 2002). However, there remain two other problems requiring observations from nature for guid- ance: 1) natural wave groups generally persist longer than the ''point break'' subjects of laboratory studies, and in addition to breaking repeatedly, they should be accom- panied by the group-forced response referred to above (Smith 2006a,b); and 2) waves in nature seldom align perfectly with the background vertical shear, as they do in the laboratory. Thus, the wave groups, breaking, and mean flow, including the response, should be studied to- gether, both to facilitate separating their effects and to include any joint effects resulting from their (nonlinear) superposition, for example, in the mean-flow response.

Here we describe our recent efforts to develop the technology to make continuous multicomponent veloc- ity estimates over a 2D ''vertical slice'' area.

2. Approach

The approach evolved from a prototype deployment of an up-looking coherent phased-array system done in 1999. That prototype provided roughly 16 independent beams covering a 228-wide pie-shaped wedge, extending upward from a deployed depth of 13.5 m. The wedge was sampled up to 52 times per second, a rate suitable for examination of the rapid evolution of the flow under breaking conditions (Fig. 1). The surface was located in each frame using a multibeam technique developed to minimize the effect of dropouts caused by dense bubble clouds. The prototype data analysis showed that image correlation techniques are effective for estimating the subsurface vertical velocities associated with the bubble clouds as well as for locating the surface (Smith 2001). It was also found that, under the actively breaking crest, a phase shift occurs in the vertical motion with depth, such that the deeper motion (1-3 m below the surface) is delayed by up to a second or so relative to the surface itself (Smith 2005). With multiple beams, the shielding of the surface by dense bubble clouds is rarely com- plete. The bubble clouds are sufficiently nonuniform that some fraction of the beams can almost always ''sneak through'' to the surface and back. Because at- tenuation is also correlated with the bubble-cloud den- sity, one may worry about a corresponding reduction in sound speed, which would affect estimates of both the range and velocity; however, for acoustic frequencies greater than 40 kHz or so, the sound speed changes duetobubbledensityarenegligible(Lamarre and Melville 1994a; Terrill and Melville 1997). Thus, from the portions observed between the densest bubble clouds, it should be feasible to reconstruct the larger- scale motions such as the ''primary vortex'' resulting from breaking.

In addition to preliminary results on wave breaking (Smith 2001, 2005), the prototype deployment has pro- vided many ideas for improvements, in terms of array and deployment geometry, most of which were implemented in the newer instrument described here. Chief among these are 1) place the instrument closer to the surface to reduce the lag between pings and facilitate coherent Doppler processing; 2) increase the number of channels to improve horizontal resolution and permit better digital beamforming; 3) expand to bistatic or multistatic transmit/receive geometry to get more than one com- ponent of velocity; and 4) extend the beamformed angle out to a full 1808 aperture to reduce spatial aliasing.

The basic concept for acoustic measurement of veloc- ities is straightforward, and the performance character- istics for the simple case of a monostatic narrow-beam Doppler sonar are fairly well known for both coherent (e.g., Lohrmann et al. 1990 and references therein) and incoherent (e.g., Pinkel and Smith 1992 and references therein) operation. Extension to a 2D area of measure- ments is also straightforward-for example, using a linear array of receivers with digital beamforming (cf. Smith 2002). As demonstrated by Smith (2002),performance along any one of many digitally formed beams is in line with that of a single simple beam. Bistatic configurations have also been used, for example, with a narrow-beam transmitter and several fan-beam receivers offset in var- ious directions to obtain multiple velocity components along the line defined by the transmitted beam (e.g., Rolland and Lemmin 1997; Hurther and Lemmin 1998, 2008; Hay et al. 2012a,b, and references therein), and have also been used extensively in atmospheric radar applications.

3. Development of the bistatic phased-array Doppler sonar system (BiPADS)

The basic approach was to build on prior proven de- signs, extending the capabilities: 64 receive channels; purely resistive tuning (relaxing the restriction to a nar- row frequency range); and two or more transmitters located at a distance in the view plane to provide mul- tiple components of velocity estimated from Doppler shifts, using a bistatic configuration (see Fig. 2).

Use of a bistatic configuration with a digital multi- beam receiver (noncollocated transmitters and an array of receivers) introduces two complexities to processing: 1) the location from which backscatter is received is no longer a simple linear function of the time delay since transmission; and 2) the Doppler shift corresponds to the velocity component along the direction that bisects the angle between transmitter, target location, and receiver (and hence is a function of both receive angle and range or time delay). While the geometry guiding this is straightforward (and is illustrated below, after defining the proposed array geometry), the high data rate contemplated here means that efficient algorithms are needed. In the following subsections, we review 1) the deployment geometry, 2) acoustic properties and limitations, and 3) bandwidth and processing re- quirements to obtain estimation errors within tolera- ble bounds.

There are several technological obstacles to be dealt with (aside from the electromechanical complexity of the system itself), including the effects of strong advec- tion on the two acoustic paths and variations in the speed of sound in the presence of bubbles. For sound at frequencies well over 40 kHz, the latter is a weak effect, with variations limited to under 0.33% or so at 40 kHz and decreasing above that (Lamarre and Melville 1994a). For round-trip distances of up to 15 m, this could account for misidentification of the scatterer location by up to 5 cm, about the resolution hoped for; however, it would require the bubble cloud to fill the entire water column over the whole acoustic path. For bubble clouds on the order a few meters in size, the error in location is pro- portionally less and can likely be ignored (conversely, the sound speed anomaly field likely cannot be estimated from the data with any useful accuracy). Similarly, a 0.33% error in the velocity estimate would likely not be detectable.

a. Bistatic geometry

Interpretation of the backscatter data can be consid- ered a two-step process: first, the time-angle input data has to be interpolated to a regular grid; second, the di- rection and errors of each velocity estimate have to be assessed.

We start by defining a grid of ''output locations'' G 5 (gx, gy); for example, a regular rectangular grid with 5-cm spacing, relative to the center of the receive array, defined to be at R 5 (0, 0) (see Fig. 3). The transmitter locations are defined as T 5 (tx, ty) on one side and S 5 (sx, sy) on the other. In the following, the method is described for a particular grid point G and transmitter T; extension to other grid points and to the other trans- mitter is straightforward.

The angle from the center of the receive array to G is

... (3.1)

(where the two-argument version of arctan is used to eliminate the 1808 ambiguity of the one-argument ver- sion). The time delay from transmitting to receiving the signal from G is just the time it takes for sound to travel the distance from T to G plus the distance from G to R: Let

... (3.2)

so the time delay for the signal from G is written as

... (3.3)

where c is the speed of sound. Standard processing of the receiver array data yields a 2D array of backscatter re- turns versus angle and delay time, so the above equation enables us to locate the returns coming from point G (etc.).

While much can be learned from just acoustic inten- sity maps, we would also like to estimate the Doppler shift and to know to what directional component of scatterer velocity this corresponds. The Doppler shift arises from the rate of change of the total pathlength. First, define the (radian) acoustic wavenumbers along each path: kT (along LT from T to G)andkR (along LR from G to R). Then, the Doppler shift from transmitter T, DvT, can be written as

... (3.4)

(note the minus sign on kR due to its definition as going from G to R). We evaluate the vector k (kT and kR) from the center frequency v0, the speed of sound c, and the known path directions:

... (3.5)

where ?R is as given above, and (similarly)

... (3.6)

Then we find

... (3.7)

where we have written V(gx, gy) 5 (Vx, Vy) and used trigonometric angle-sum and angle-difference identities. As anticipated at the start of the section, the direction of the Doppler estimate component bisects the angle TGR (see Fig. 3). Also, the size of the Doppler shift is reduced relative to the monostatic case, approaching zero response as the angles approach opposing directions. Graphically, in Fig. 3, the Doppler shift estimates correspond to ve- locity components perpendicular to the ellipses of con- stant time delay, and where the two sets of ellipses are more nearly perpendicular to each other there is more nearly equal resolution of both components of velocity.

Naming the corresponding velocity component ''VT,'' it would be estimated as

... (3.8)

Note that the paths, angles, Doppler estimates, and sound vector ks are all functions of the gridpoint loca- tion G as well as which transmit is active. This in effect yields a set of transforms to go from the beamformed angle-delay matrix for the backscatter from each trans- mitter to output arrays of intensity and Doppler shift (or VT and VS) versus location on a uniform Cartesian grid.

b. Combining the transmits

The two offset transmitters can provide Doppler esti- mates of velocity components in two different directions, generally not orthogonal but hopefully not parallel ei- ther. Since the Doppler estimates inherently contain noise (as discussed below), simple Gaussian elimina- tion is unstable, and an estimate explicitly accounting for the noise is preferred. Here we employ the ''optimal estimate'' derived by Smith (2002). We define the di- rection of the two estimated components at each grid point, as evaluated above for the transmitter T at grid point G:let

... (3.9)

be the direction of the estimated velocity component VT (as above), and for the ''second transmit'' component VS, let

... (3.10)

Now denote the RMS Doppler velocity errors (to be discussed below) associated with these components as ęT and ęS , and normalize them by the (assumed iso- tropic) RMS true velocity variances: let

... (3.11)

Then we write estimates for the Cartesian components of velocity in the form

... (3.12)

Smith's (2002) solution corresponds to

... (3.13a)

... (3.13b)

... (3.13c)

... (3.13d)

(note that the denominators are all the same). In the (unrealistic) limit of no noise, this reduces to Gaussian elimination. In the limit of large noise on one of the transmitters (S, thus NS /'), that estimate is ignored, and the other is just reduced by its own error, consistent with minimizing the rms error of the result, for example,

... (3.14)

4. Acoustic noise and Doppler error

There are two kinds of noise to deal with: acoustic ambient noise and Doppler estimation error. Because there is a finite amount of ambient noise in the ocean and the backscattered signal attenuates with distance through the water, the signal will eventually fade beyond detection and the estimate becomes purely noise. How- ever, at close range, the acoustic signal-to-noise ratio is quite large, and the ambient noise can be neglected. Even so, there remains a finite amount of Doppler error.

For ping-to-ping coherent processing, returns from the same time delay for consecutive ping pairs are cor- related to estimate the Doppler shift (Lhermitte and Serafin 1984; Zedel 2008). In this case, the Doppler noise depends on the inherent decorrelation time of the scat- terers, the speed of advection of scatterers into and out of the sample volume, and on the reverberation character- istics of the earlier transmissions over the time period of the subsequent ones. A potential advantage of sampling a 2D planar field is that the correlations (or covariances) can be estimated from averages on a moving trajectory on that plane, reducing the advective decorrelation effect. Presuming the carrier (center) frequency is previously remixed to zero, there should be no problem with phase stability on a moving trajectory.

In deep water, as in the Research Platform (R/P) Floating Instrument Platform (FLIP) trials mentioned above, the pings decayed rapidly enough that the time between pings was mainly limited by the desired range: with 50 pings per second, ranges up to 15 m were sampled, but the correlations remaining at the 200-kHz frequencies employed were too small and erratic to encourage Doppler estimation. In shallow water, in contrast, we found that reverberation persists for too long compared to the decorrelation time scale to get useful data.

To estimate the net noise variance s2 relative to the signal strength, we make use of the coherence:

... (4.1)

where Coh is the ratio of the time-lagged covariance hCti to the received variance hs2i, based on an average, denoted by the angle brackets (hi), over M sample pairs:

... (4.2)

where the asterisk (*) denotes the complex conjugate (see appendix A).

Given the relative noise level s2, the error associated with the covariances upon which the Doppler velocities are based can be estimated. This turns out to depend also on how the M pairs to average are arranged. One way is to average over successive samples (pings) in time:

... (4.3)

where the time between samples (pings) is t . In contrast, the other way involves independent sample pairs (e.g., in a coherent processing scheme, range averaging the products from a single ping pair, assuming each range bin is independent), so with r as the range index and t as the time between pings (again):

... (4.4)

When the true signal covariance is all real (say 5 1), the Doppler error shows up in the imaginary part of the var- iance. For independent sample averages, the covariance error is split equally between the real and imaginary parts:

... (4.5)

where the ''1'' (real) corresponds to the true covariance, and the term on the right is the net Doppler error vari- ance (see appendix B). For sequential sample averages, the total noise variance is the same (i.e., double that term), but the repeated entries in successive sample pairs have the effect of shifting a portion of the error from the imaginary to the real part:

... (4.6a)

... (4.6b)

(see appendix B; Fig. 4). This means the Doppler error is less for sequential sample averaging than for independent averages with the same number of sample pairs, unless the real-part variance is large enough to cause many er- roneous negative real-part values. For example, assuming normally distributed noise, if s2 . M/6, then more than 1% of the real-part error values would exceed the signal strength (51), and half of these would be on the negative side. Note that the above equation can be rewritten with the true covariance represented as Aeiu , but both the math and discussion would be more cumbersome.

5. Prototype bistatic deployment

Trials with a prototype BiPADS were done near the end of the Scripps Pier (La Jolla, California; see Fig. 5) to assess the potential for estimating two in-plane com- ponents of velocity, using transmitters offset in the same plane as the receive array (see Fig. 3). The prototype BiPADS was deployed off the side of the pier, on a frame with long arms that unfold to hold the transmitters at fixed locations ;4 m to either side of the receive array. The base configuration is 200-kHz center frequency. Beamforming (to 63 beams, dropping the outermost ''Nyquist angle'') was done via time delay; the angular resolution at the center angle is about 1.38,increasing toward the edges. In these trials, we tested the instru- ment's capabilities, first trying a variety of sample and ping rates, and then a variety of center frequencies. Ini- tially. we ran it at the 12.5-kHz sample rate, but we found we could reliably increase this to a maximum ;14.7-kHz rate, limited by the bandwidth of the total throughput from all 64 channels to the host computer. Next we tested the ping interval, starting with a conservative 12.5 pings per second and noting that it took 40 ms before the re- turns decayed substantially. So, we increased the repeat rate to 25 pings per second (1/40 ms) to see if we could get any detectable coherence. While we got fairly satisfactory imaging of the intensity returns (see Fig. 6), 40 ms is too long to obtain usable coherence between pings: they were uniformly statistically indistinguishable from zero.

For the sample frame illustrated in Fig. 6, the pings alternated between transmitters (at 40-ms intervals), and a pair of pings was averaged, including transmits from both sides. This was mainly to establish whether we could alternate between transmitters and image the same locations (we did not even attempt ping pairs or sequences from each side). The surface is about 6 m up on average, the waves (swell) were ;1-m peak to trough, and the wind was dead calm. The surface is dimmer than usual in such images, both because of the calm condi- tions and because we put sound-absorbing panels above and below the transmitter bars, limiting the aperture to ;308, reducing specular hits on the surface above the receivers as well as bottom reflections. There are discrete intensity features at the surface and in the water that move with the waves; the criterion for ''satisfactory im- aging'' here is that we could successfully navigate most of these bright spots from the two transmitters to coincide (there will always be some ''outliers'' due to directional aliasing, slightly different sample volumes, and aniso- tropic scatterers). It is reassuring that these spots are mostly in the same place no matter which side the sound comes from; that is, that our mapping from the data- stream to physical coordinates is successful: proper lo- cation of these intensity features is sensitive to the detailed geometry of deployment, and to the precision of timing. Unfortunately, while the mapping from raw data to a Cartesian grid was successful, the reverberations remained too strong for too long to be able to get any ping-to-ping coherence for Doppler estimation.

More interesting in some ways were the results of our tests of the system's ability to run off resonance (the transducers are all designed for 200 kHz). We success- fullyranat 170,180, 190,200, 225,250, 275,300, 325,350, 400, 450, 500, 550, and even 600 (kHz). Some aspects of note: at 170 kHz, the receive array is ''endfire'' capable (i.e., can beamform over 1808). At 400 kHz, there was interference from an unidentified source; the ability to shift operating frequency is critical for avoiding such nuisances. At all frequencies, the signal was still strong at the maximum distances (;10 m) for the same power setting, with little change in intensity. As the frequency increases, attenuation should increase, but the decor- relation time decreases at about the same rate through this range.

To summarize the results from these deployments, we achieved most of the improvements suggested by prior work but failed at the final step:

1) Move closer to surface. Check: 6 m.

2) Increase the number of elements for better beam- forming. Check: 64.

3) Expand to bistatic geometry. Check.

4) Extend beamforming to 1808. Check, at 170 kHz.

5) Use purely resistive tuning to broaden range of usable frequencies. Check.

6) Ping fast enough for coherence. Fail (at least in shallow water).

6. Discussion

The logical approach is to return to reliable and robust incoherent Doppler processing at a much higher fre- quency to get the resolution desired (say, 10-cm spatial resolution and 10 cm s21 velocity resolution). How high a frequency and what bandwidth are required to obtain this? Using repeat-sequence coding (Pinkel and Smith 1992), the single-ping range-velocity uncertainty prod- uct DRDV is estimated as (Pinkel and Smith 1992)

... (6.1)

where c is the speed of sound (;1500 m s21), f is the center frequency (Hz), and L is the number of ''bits'' (samples) in the code. The maximum number of bits, or ''code length,'' that can be used depends on the fre- quency, bandwidth (Df ), and a ''wrap velocity'' (Vw) chosen to exceed 99% of actual velocity values antici- pated (we typically choose Vw ' 2ms21). Setting the covariance phase anomaly to p at a time lag of one code length yields the requirement

... (6.2)

So, for example, if we assume we can achieve a bandwidth of 10% of the carrier frequency (and with the values quoted above for c and Vw), then we find L # 18.75, independent of frequency. The code length must be an integer, and it must perform well. There is a good ''Barker-like'' code of length 19; the next shorter good performer is 13.

Putting (6.2) into (6.1) yields the net resolution versus frequency and bandwidth:

... (6.3)

With the above values of c, Vw, and bandwidth ratio, we get DRDV 5 (25 321 m2 s22)f21. At 2 MHz, this yields an rms error of 0.012 66 m2 s21, corresponding to roughly (11 cm)(11 cm s21), close to the target resolutions. With a 200-kHz sample rate, the 13-bit code would corre- spond to just under 5 cm in range, so the desired range resolution would correspond to two code lengths; in contrast, the 19-bit code would correspond to just over 7 cm. However, both codes could accommodate ; 14-cm range bins (with a correspondingly smaller velocity er- ror): 3 3 13 bits or 2 3 19 bits are close.

The other limitation is the maximum attainable range: as frequency goes up, so does attenuation. The beam geometry considered here is close enough to that used in Smith (1989) that we can examine the range limit at maximum usable power (limited by nonlinear satura- tion) as estimated there. Unfortunately, those calcula- tions stop just shy of 1 MHz; but above 200 kHz or so, the maximum attainable range goes very nearly like 1/f,so we can extrapolate: this yields the (surely optimistic) estimates of 80 m at 1 MHz, or 40 m at 2 MHz. In prac- tice, there is also the matter of less efficient scatterers (no bubbles are small enough to resonate at megahertz frequencies). However, experience with commercial sonars operated near 1 MHz in incoherent mode in- dicates ranges of 11 m are achievable, which is in line with our desired working range (on the order of 10 m). Limiting the range to 10-15 m also means that there is time for 50 pings per second or more, enabling signifi- cant additional noise reduction by time averaging; thus, it may also be feasible to operate down to 1 MHz, if need be, to get more range. However, it looks likely that 2 MHz will be workable, maybe even optimal.

With incoherent processing, there is a third source of noise, due to ''self-clutter'' from the nonoverlapping portion of each sample pair. The analysis of section 4 still applies to the other sources of noise (i.e., to both ambient noise and anomalous scatterers); however, the roles of range averaging versus ping averaging are reversed: for incoherent processing, range averaging is ''sequential,'' while each ping is independent. This could be helpful in extending the usable range, where the ambient noise is the limiting factor: for example, there can be progres- sively more range averaging at farther ranges.

Acknowledgments. The early prototype field work (1999) and initial analysis were supported by ONR (Grant N00014-02-1-0855), with additional analysis and modeling supported by NSF (Grant OCE 06-23679). The new instrument development was supported by NSF (Grant OCE-06-23679). Thanks are due to the OPG en- gineering team of M. Goldin, M. Bui, Tyler Hughens, and A. Aja for their tireless efforts to design, construct, de- ploy, and operate the BiPADS, and to help develop the software to acquire and analyze the data. Thanks are also due to many others for the many interesting and educa- tional conversations.


Agrawal, Y. C., E. A. Terray, M. A. Donelan, P. A. Hwang, A. J. Williams III, W. M. Drennan, K. K. Kahma, and S. A. Krtaigorodski, 1992: Enhanced dissipation of kinetic energy beneath surface waves. Nature, 359, 219-220, doi:10.1038/ 359219a0.

Ardhuin, F., N. Rascle, and K. A. Belibassakis, 2008: Explicit wave-averaged primitive equations using a generalized La- grangian mean. Ocean Modell., 20, 35-60, doi:10.1016/ j.ocemod.2007.07.001.

Bowen, A. J., 1969: The generation of longshore currents on a plane beach. J. Mar. Res., 27, 206-215.

_____, D. L. Inman, and V. P. Simmons, 1968: Wave 'setdown' and wave setup. J. Geophys. Res., 73, 2569-2577, doi:10.1029/ JB073i008p02569.

Craik, A. D. D., 1977: The generation of Langmuir circulation by an instability mechanism. J. Fluid Mech., 81, 209-223, doi:10.1017/S0022112077001980.

_____, and S. Leibovich, 1976: A rational model for Langmuir circulations. J. Fluid Mech., 73, 401-426, doi:10.1017/ S0022112076001420.

Edson, J. B., and C. W. Fairall, 1994: Spray droplet modeling. 1. Lagrangian model simulation of the turbulent transport of evaporating droplets. J. Geophys. Res., 99, 25 295-25 311, doi:10.1029/94JC01883.

Fan, Y., I. Ginis, and T. Hara, 2009: The effect of wind-wave- current interaction on air-sea momentum fluxes and ocean response in tropical cyclones. J. Phys. Oceanogr., 39, 1019- 1034, doi:10.1175/2008JPO4066.1.

Farmer, D. M., S. Vagle, and A. D. Booth, 1998: A free-flooding acoustical resonator for measurement of bubble size distributions. J. Atmos. Oceanic Technol., 15, 1132-1146, doi:10.1175/1520-0426(1998)015,1132:AFFARF.2.0.CO;2.

Garrett, C., 1976: Generation of Langmuir circulations by surface waves-A feedback mechanism. J. Mar. Res., 34, 117-130.

_____, and J. Smith, 1976: On the interaction between long and short surface waves. J. Phys. Oceanogr., 6, 925-930, doi:10.1175/ 1520-0485(1976)006,0925:OTIBLA.2.0.CO;2.

Gemmrich, J. R., and D. M. Farmer, 1999: Observations of the scale and occurrence of breaking surface waves. J. Phys. Oceanogr., 29, 2595-2606, doi:10.1175/1520-0485(1999)029,2595: OOTSAO.2.0.CO;2.

_____, and _____, 2004: Near-surface turbulence in the presence of breaking waves. J. Phys. Oceanogr., 34, 1067-1086, doi:10.1175/ 1520-0485(2004)034,1067:NTITPO.2.0.CO;2.

_____, M. L. Banner, and C. Garrett, 2008: Spectrally resolved energy dissipation rate and momentum flux of breaking waves. J. Phys. Oceanogr., 38, 1296-1312, doi:10.1175/ 2007JPO3762.1.

Groeneweg, J., and G. Klopman, 1998: Changes of the mean ve- locity profiles in the combined wave-current motion de- scribed in a GLM formulation. J. Fluid Mech., 370, 271-296, doi:10.1017/S0022112098002018.

_____, and J. A. Battjes, 2003: Three-dimensional wave effects on a steady current. J. Fluid Mech., 478, 325-343, doi:10.1017/ S0022112002003476.

Hasselmann, K., 1971: On the mass and momentum transfer be- tween short gravity waves and larger-scale motions. J. Fluid Mech., 50, 189-205, doi:10.1017/S0022112071002520.

Hay, A. E., L. Zedel, R. Cheel, and J. Dillon, 2012a: Observations of the vertical structure of turbulent oscillatory boundary layers above fixed roughness beds using a prototype wideband co- herent Doppler profiler: 1. The oscillatory component of the flow. J. Geophys. Res., 117, C03005, doi:10.1029/2011JC007113.

_____, _____, _____, and _____, 2012b: Observations of the vertical structure of turbulent oscillatory boundary layers above fixed roughness using a prototype wideband coherent Doppler profiler: 2. Turbulence and stress. J. Geophys. Res., 117, C03006, doi:10.1029/2011JC007114.

Herbers, T. H. C., and R. T. Guza, 1991: Wind-wave nonlinearity observed at the sea floor. Part I: Forced-wave energy. J. Phys. Oceanogr.,21,1740-1761,doi:10.1175/1520-0485(1991)021,1740: WWNOAT.2.0.CO;2.

_____, S. Elgar, and R. T. Guza, 1994: Infragravity-frequency (0.005- 0.05 Hz) motions on the shelf. Part I: Forced waves. J. Phys. Oceanogr., 24,917-927, doi:10.1175/1520-0485(1994)024,0917: IFHMOT.2.0.CO;2.

Hurther, D., and U. Lemmin, 1998: A constant-beam-width transducer for 3D acoustic Doppler profile measurements in open-channel flows. Meas. Sci. Technol., 9, 1706-1714, doi:10.1088/0957-0233/9/10/010.

_____, and _____2008: Improved turbulence profiling with field- adapted acoustic Doppler velocimeters using a bifrequency Doppler noise suppression method. J. Atmos. Oceanic Technol., 25, 452-463, doi:10.1175/2007JTECHO512.1.

Jiang, J.-Y., and R. L. Street, 1991: Modulated flows beneath wind- ruffled, mechanically generated water waves. J. Geophys. Res., 96, 2711-2721, doi:10.1029/90JC02259.

Kemp, P. H., and R. R. Simons, 1982: The interaction of waves and a turbulent current; waves propagating with the current. J. Fluid Mech., 116, 227-250, doi:10.1017/S0022112082000445.

_____, and _____, 1983: The interaction of waves and a turbulent current: Waves propagating against the current. J. Fluid Mech., 130, 73-89, doi:10.1017/S0022112083000981.

Kukulka,T.,A.J.Plueddemann,J.H.Trowbridge,andP.P.Sullivan, 2009: Significance of Langmuir circulation in upper ocean mixing: Comparison of observations and simulations. Geo- phys. Res. Lett., 36, L10603, doi:10.1029/2009GL037620.

Kumar, N., G. Voulgaris, J. C. Warner, and M. Olabarrieta, 2012: Implementation of the vortex force formalism in the coupled ocean-atmosphere-wave-sediment transport (COAWST) mod- eling system for inner shelf and surf zone applications. Ocean Modell., 47, 65-95, doi:10.1016/j.ocemod.2012.01.003.

Lamarre, E., and W. K. Melville, 1994a: Sound-speed measure- ments near the ocean surface. J. Acoust. Soc. Amer., 96, 3605- 3616, doi:10.1121/1.410578.

_____, and _____, 1994b: Void-fraction measurements and sound- speed fields in bubble plumes generated by breaking waves. J. Acoust. Soc. Amer., 95, 1317-1328, doi:10.1121/1.408572.

Langmuir, I., 1938: Surface motion of water induced by wind. Science, 87, 119-123, doi:10.1126/science.87.2250.119.

Leibovich, S., 1980: On wave-current interaction theories of Langmuir circulations. J. Fluid Mech., 99, 715-724, doi:10.1017/S0022112080000857.

Lhermitte, R., and R. Serafin, 1984: Pulse-to-pulse coherent Doppler sonar signal processing techniques. J. Atmos. Oceanic Technol., 1, 293-308, doi:10.1175/1520-0426(1984)001,0293: PTPCDS.2.0.CO;2.

Li, M., K. Zahariev, and C. Garrett, 1995: Role of Langmuir cir- culation in the deepening of the ocean surface mixed layer. Science, 270, 1955-1957, doi:10.1126/science.270.5244.1955.

Lohrmann, A., B. Hackett, and L. P. RØed, 1990: High resolution measurements of turbulence, velocity and stress using a pulse- to-pulse coherent sonar. J. Atmos. Oceanic Technol., 7, 19-37, doi:10.1175/1520-0426(1990)007,0019:HRMOTV.2.0.CO;2.

Longuet-Higgins, M. S., 1953: Mass transport in water waves. Philos. Trans. Roy. Soc. London, A245, 535-581, doi:10.1098/ rsta.1953.0006.

_____, 1969: A nonlinear mechanism for the generation of sea waves. Proc. Roy. Soc. London, A311, 371-389, doi:10.1098/ rspa.1969.0123.

_____, 1970a: Longshore currents generated by obliquely incident sea waves: 1. J. Geophys. Res., 75, 6778-6789, doi:10.1029/ JC075i033p06778.

_____, 1970b: Longshore currents generated by obliquely incident sea waves: 2. J. Geophys. Res., 75, 6790-6801, doi:10.1029/ JC075i033p06790.

_____, and R. W. Stewart, 1960: Changes in the form of short gravity waves on long waves and tidal currents. J. Fluid Mech., 8, 565- 583, doi:10.1017/S0022112060000803.

_____, and _____, 1961: The changes in amplitude of short gravity waveson steadynon-uniformcurrents.J. FluidMech., 10,529- 549, doi:10.1017/S0022112061000342.

_____, and _____, 1962: Radiation stress and mass transport in gravity waves, with application to 'surf-beats.' J. Fluid Mech., 13, 481- 504, doi:10.1017/S0022112062000877.

_____, and _____, 1964: Radiation stresses in water waves; A physical discussion, with applications. Deep-Sea Res. Oceanogr. Abstr., 11, 529-562, doi:10.1016/0011-7471(64)90001-4.

McWilliams, J. C., and J. M. Restrepo, 1999: The wave-driven ocean circulation. J. Phys. Oceanogr., 29, 2523-2540, doi:10.1175/ 1520-0485(1999)029,2523:TWDOC.2.0.CO;2.

_____, and P. P. Sullivan, 2000: Vertical mixing by Langmuir cir- culations. Spill Sci. Technol. Bull., 6, 225-237, doi:10.1016/ S1353-2561(01)00041-X.

_____, _____, and C. H. Moeng, 1997: Langmuir turbulence in the ocean. J. Fluid Mech., 334, 1-30, doi:10.1017/S0022112096004375.

_____, J. M. Restrepo, and E. M. Lane, 2004: An asymptotic theory for the interaction of waves and currents in coastal waters. J. Fluid Mech., 511, 135-178, doi:10.1017/S0022112004009358.

Mellor, G., 2003: The three-dimensional current and surface wave equations. J. Phys. Oceanogr., 33, 1978-1989, doi:10.1175/ 1520-0485(2003)033,1978:TTCASW.2.0.CO;2.

Melville, W. K., and P. Matusov, 2002: Distribution of breaking waves at the ocean surface. Nature, 417, 58-63, doi:10.1038/ 417058a.

_____, F. Veron, and C. J. White, 2002: The velocity field under breaking waves: Coherent structures and turbulence. J. Fluid Mech., 454, 203-233, doi:10.1017/S0022112001007078.

Monismith, S. G., E. A. Cowen, H. M. Nepf, J. Magnaudet, and L. Thais, 2007: Laboratory observations of mean flows under surface gravity waves. J. Fluid Mech., 573, 131-147, doi:10.1017/ S0022112006003594.

Nepf, H. M., 1992: The production and mixing effects of Langmuir circulations. Ph.D. thesis, Stanford University, 268 pp.

_____, C.H.Wu,andE.S.Chan,1998:Acomparisonoftwo-andthree- dimensional wave breaking. J. Phys. Oceanogr., 28, 1496-1510, doi:10.1175/1520-0485(1998)028,1496:ACOTAT.2.0.CO;2.

Newberger, P.A., and J. S. Allen,2007: Forcing a three-dimensional, hydrostatic, primitive-equation model for application in the surf zone: 2. Application to DUCK94. J. Geophys. Res., 112, C08019, doi:10.1029/2006JC003474.

Özkan-Haller, H. T., and Y. Li, 2003: Effects of wave-current in- teraction on shear instabilities of longshore currents. J. Geo- phys. Res., 108, 3139, doi:10.1029/2001JC001287.

Pattison,M. J., and S.E. Belcher, 1999: Production rates ofsea-spray droplets. J. Geophys. Res., 104, 18 397-18 407, doi:10.1029/ 1999JC900090.

Phillips, W. R. C., 2002: Langmuir circulations beneath growing or decaying surface waves. J. Fluid Mech., 469, 317-342, doi:10.1017/S0022112002001908.

Pinkel, R., and J. A. Smith, 1992: Repeat-sequence coding for im- provedprecisionofDoppler sonar andsodar. J.Atmos.Oceanic Technol., 9, 149-163, doi:10.1175/1520-0426(1992)009,0149: RSCFIP.2.0.CO;2.

Rolland, T., and U. Lemmin, 1997: A two-component acoustic velocity profiler for use in turbulent open-channel flow. J. Hydraul.Res., 35, 545-561, doi:10.1080/00221689709498410.

Skyllingstad, E. D., and D. W. Denbo, 1995: An ocean large-eddy simulation of Langmuir circulations and convection in the sur- face mixed layer. J. Geophys. Res., 100, 8501-8522, doi:10.1029/ 94JC03202.

Smith, J. A., 1980: Waves, currents, and Langmuir circulation. Ph.D. thesis, Dalhousie University, 242 pp.

_____, 1989: Doppler sonar and surface waves: Range and reso- lution. J. Atmos. Oceanic Technol., 6, 680-696, doi:10.1175/ 1520-0426(1989)006,0680:DSASWR.2.0.CO;2.

_____, 2001: Phased-array Doppler sonar measurements of near- surface motion: Langmuir circulation, surface waves, and breaking. Acoustical Oceanography, T. G. Leighton et al., Ed., Vol. 23, Part 2, Proceedings of the Institute of Acoustics Conference on Acoustical Oceanography, Bath University Press, 163-170.

_____, 2002: Continuous time-space sampling of near-surface veloc- ities using sound. J. Atmos. Oceanic Technol., 19, 1860-1872, doi:10.1175/1520-0426(2002)019,1860:CTSSON.2.0.CO;2.

_____, 2005: Acoustically observed bubble cloud evolution under deep-water breaking waves. Proc. Underwater Acoustic Mea- surements: Technologies and Results, Heraklion, Crete, Greece, Acoustical Society of America, 6 pp. [Available online at (UAM).pdf.]

_____, 2006a: Observed variability of ocean wave stokes drift, and the Eulerian response to passing groups. J. Phys. Oceanogr., 36, 1381-1402, doi:10.1175/JPO2910.1.

_____, 2006b: Wave-current interactions in finite depth. J. Phys. Oceanogr., 36, 1403-1419, doi:10.1175/JPO2911.1.

_____, and C. Brulefert, 2010: Evolution of persistent wave groups. J. Phys. Oceanogr., 40, 67-84, doi:10.1175/2009JPO4149.1.

Stokes, G. G., 1847: On the theory of oscillatory waves. Trans. Cambridge Philos. Soc., 8, 441-455.

Sullivan, P. P., J. C. McWilliams, and W. K. Melville, 2004: The oceanic boundary layer driven bywave breaking with stochastic variability. Part 1. Direct numerical simulations. J. Fluid Mech., 507, 143-174, doi:10.1017/S0022112004008882.

_____, _____, and _____, 2007: Surface gravity wave effects in the oceanic boundary layer: Large-eddy simulation with vortex force and stochastic breakers. J. Fluid Mech., 593, 405-452, doi:10.1017/S002211200700897X.

Swan, C., I. P. Cummins, and R. L. James, 2001: An experimental study of two-dimensional surface water waves propagating on depth-varying currents. Part 1. Regular waves. J. Fluid Mech., 428, 273-304, doi:10.1017/S0022112000002457.

Terray, E. A., M. A. Donelan, Y. C. Agrawal, W. M. Drennan, K. K. Kahma, A. J. Williams, P. A. Hwang, and S. A. Kitaigorodskii, 1996: Estimates of kinetic energy dissipation under break- ing waves. J. Phys. Oceanogr., 26, 792-807, doi:10.1175/ 1520-0485(1996)026,0792:EOKEDU.2.0.CO;2.

Terrill, E. J., and W. K. Melville, 1997: Sound-speed measurements in the surface-wave layer. J. Acoust. Soc. Amer., 102, 2607- 2625, doi:10.1121/1.420315.

_____, and _____, 2000: A broadband acoustic technique for mea- suring bubble size distributions: Laboratory and shallow water measurements. J. Atmos. Oceanic Technol., 17, 220-239, doi:10.1175/1520-0426(2000)017,0220:ABATFM.2.0.CO;2.

Thorpe, S. A., 1982: On the clouds of bubbles formed by breaking wind-waves in deep water, and their role in air-sea gas transfer. Philos. Trans. Roy. Soc. London, A304, 155-210, doi:10.1098/rsta.1982.0011.

Vagle, S., and D. M. Farmer, 1992: The measurement of bubble- size distributions by acoustical backscatter. J. Atmos. Oceanic Technol., 9, 630-644, doi:10.1175/1520-0426(1992)009,0630: TMOBSD.2.0.CO;2.

Veron, F., and W. K. Melville, 1999: Pulse-to-pulse coherent Doppler measurements of waves and turbulence. J. Atmos. Oceanic Technol., 16, 1580-1597, doi:10.1175/1520-0426(1999)016,1580: PTPCDM.2.0.CO;2.

_____, _____, and L. Lenain, 2008: Wave-coherent air-sea heat flux. J. Phys. Oceanogr., 38, 788-802, doi:10.1175/ 2007JPO3682.1.

Wu, C. H., and H. M. Nepf, 2002: Breaking criteria and energy losses for three-dimensional wave breaking. J. Geophys. Res., 107, 3177, doi:10.1029/2001JC001077.

Zedel, L., 2008: Modeling pulse-to-pulse coherent Doppler sonar. J. Atmos. Oceanic Technol., 25, 1834-1844, doi:10.1175/ 2008JTECHO585.1.


Scripps Institution of Oceanography, La Jolla, California

(Manuscript received 17 September 2013, in final form 27 March 2014)

Corresponding author address: Jerome A. Smith, Marine Phys- ical Laboratory, Scripps Institution of Oceanography, University of California, San Diego, 9500 Gilman Drive, Mail Code 0213, La Jolla, CA 92093-0213.


DOI: 10.1175/JTECH-D-13-00187.1

(ProQuest: Appendix omitted.)

For more stories covering the world of technology, please see HispanicBusiness' Tech Channel

Source: Journal of Atmospheric and Oceanic Technology

Story Tools Facebook Linkedin Twitter RSS Feed Email Alerts & Newsletters