## 1. Introduction

Studies of coastal environments often require high-resolution velocity measurements throughout the water column over a broad area. The need for high resolution in space and time suggests use of a high-frequency acoustic Doppler current profiler (ADCP) that can detect finescale variations in currents (Hickey et al. 1998; Moum et al. 2003). However, such instruments have a limited vertical range and can sample only a small part of a typical ∼100-m coastal water column from a fixed depth. The need for broad horizontal coverage and the inherent contradiction between range and resolution suggest deployment of a 1200-kHz ADCP on a towed, undulating vehicle. Use of the enhanced sampling rate via RD Instruments’ (RDI’s) mode 12 feature is useful to improve accuracy. Mode 12 uses smaller depth cells, making it possible to measure highly sheared profiles of fast flows. In mode 1, small depth cells can result in noisy data unless long time averages are performed. Mode 12 gets around this by increasing the ping rate, allowing the depth cell size to be small without increasing the averaging time and data noise. In mode 12, each ping consists of a sequence of subpings, which are summed to ping averages.

The River Influence on the Shelf Ecosystems (RISE) project focuses on sampling the highly mobile Columbia River plume environment. The small time and space scales of the plume require rapid, high-resolution (in time and space) sampling. To meet these sampling criteria, we mounted an RDI 1200-kHz ADCP on a TRIAXUS towfish, developed by MacArtney Underwater Technology Group Company, which was steerable such that it surfaced well outside the visible ship wake (MacArtney Underwater Technology Group 2005). The TRIAXUS towfish is steerable in 3D and, in a configuration installed on the R/V *Point Sur*, can reach depths of ∼200 m. The ADCP on the TRIAXUS measures ocean current velocities relative to the TRIAXUS vehicle. To obtain the current velocities in the earth coordinates, we have to determine the TRIAXUS movements relative to the ship or the earth coordinates. The primary difficulty in carrying out this procedure is to precisely determine the TRIAXUS vehicle movements. Visbeck (2002) used a linear inverse method to obtain depth profiles of the velocities from lowered ADCP (LADCP) sampling data. In this note, we introduce a simple method to process the TRIAXUS ADCP measured velocity data and recover the ocean current velocities in the earth coordinates.

## 2. Methodology

*û*is the velocity removing its mean;

*S*1 is the number of velocity profile cells;

*S*2 is the number of velocity profiles; subscript

*p*represents the

*p*th velocity profile cell; and superscripts

*A*and

*B*denote different velocity profiles. The other is the rms of the detrended velocities,

*û*′ is the detrended velocity. The rms

_{1}represents the difference between velocity profiles and rms

_{2}shows the data uncertainty within velocity profiles. For the case of Fig. 1, rms

_{1}= 0.08 m s

^{−1}, and rms

_{1}= 0.05 m s

^{−1}. Here, we have rms

_{1}< 2 rms

_{2}.

To correct the TRIAXUS ADCP-measured velocities, we had another 1200-kHz ADCP mounted on a pole fixedly attached to the ship. The ship-mounted ADCP had the same configurations as that on the TRIAXUS, and measured a velocity layer from 3.0- to 20-m depth. The vessel velocities were determined using a shipboard GPS, and therefore, we can obtain absolute current velocities in the earth coordinates, which are employed as references for the TRIAXUS ADCP-measured velocities. However, the TRIAXUS vehicle sometimes cycled to a depth beyond the range of ship-mounted ADCP measurements, so that the TRIAXUS and ship-mounted pole ADCPs did not overlap. In this situation, the following method is used to solve the problem. Figure 2 is a schematic map of the two ADCP measurements. In the figure, *u _{i,j}* and

*u*represent the velocities measured by the TRIAXUS ADCP. The subscripts

_{k,j}*i*and

*k*denote the

*i*th and

*k*th ensemble, and the subscript

*j*represents the

*j*th bin, so

*j*is a function of depth;

*U*and

_{i,j}*U*are the ship-mounted ADCP measurements and have already been corrected to be in the earth coordinates by using the ship’s GPS navigation data. There are offsets between

_{k,j}*u*and

_{i,j}*U*and between

_{i,j}*u*and

_{k,j}*U*, resulting from TRIAXUS vehicle motions relative to the vessel. For the method, there are two steps as below.

_{k,j}- For the ensembles with overlapping bins between the two ADCPs (ensembles
*i*,*i*+ 1,*k*− 1,*k*in Fig. 2), the TRIAXUS vehicle velocities (*u*and^{t}_{i}*u*) relative to the vessel (at ensembles^{t}_{k}*i*and*k*, respectively) satisfywhere*N*is the bin number (46),*M*is the overlapping bin number, and the superscript*t*represents TRIAXUS parameters. Using Eq. (3), we can determine TRIAXUS vehicle velocities*u*and^{t}_{i}*u*, as well as^{t}_{k}*u*^{t}_{i+}_{1}and*u*^{t}_{k−}_{1}in the ensembles with overlapping bins between TRIAXUS and ship-mounted ADCPs. - For the ensembles without overlapping bins between the two ADCPs (ensembles
*i*+ 2,*i*+ 3, . . . ,*k*− 3,*k*− 2 in Fig. 2), we determine the TRIAXUS vehicle velocities by comparing the neighboring ensembles. In Fig. 2, one can see that*u*_{i+}_{2}and_{,j}*u*_{i+}_{1}have vertical overlapping bins. We assume the two neighboring ensembles of velocity measurements have a negligible difference, which gives_{,j}where*L*_{1}represents the first bin of ensemble*i*+ 1 overlapping with*i*+ 2;*u*^{c}_{i}_{+1}denotes the corrected TRIAXUS ADCP-measured current velocity. Using Eq. (4), we get TRIAXUS vehicle velocity*u*^{t}_{i}_{+2}, and in same manner, we solve the TRIAXUS velocities at ensembles*i*+ 3,*i*+ 4, . . . ,*k*− 3, and*k*− 2 step by step. However, this procedure may magnify errors from*i*+ 3 to*k*− 2, since we do not consider the difference between two neighbor ensembles. The errors could be accumulated from*i*+ 3 to*k*− 2, so the accuracies of*u*and_{c}*u*decrease step by step. For ensemble_{t}*k*− 1, if we neglect the difference between*k*− 2 and*k*− 1, we havewhere*L*_{2}represents the first bin of ensemble*k*− 1 overlapping with*k*− 2. Using Eq. (5), we can obtain the TRIAXUS velocity at time*k*− 2 (*u*^{t}_{k−}_{2}). Furthermore,*u*^{t}_{k−}_{3},*u*^{t}_{k−}_{4}, . . . , and*u*^{t}_{i+2}are solved. The derivation procedure is from*k*− 2 to*i*+ 2, and therefore, the accuracies decrease from*k*− 2 to*i*+ 2, which are opposite to those of first group of the solution. Let the first solution group be referred to as*u*^{t,i}_{m,j}, and the second as*u*^{t,k}_{m,j}, where*m*represents ensembles*i*+ 2,*i*+ 3, . . . ,*k*− 3,*k*− 2. To make optimal estimations of the TRIAXUS velocities between*i*+ 2 and*k*− 2, we use a weighted average to combine two groups of solutions, namelywhere*α*_{1}and*α*_{2}are weighted coefficients related to the distances between*i*and*m*and between*m*and*k*, respectively. The corrected TRIAXUS ADCP-measured velocity can be derived by

## 3. Results

Figures 3a and 3b show a segment of TRIAXUS ADCP-measured *u* and *υ* velocities, respectively, from which the ship navigation velocities were removed. The segment is between 46.9990°N, 125.3376°W and 47.0000°N, 125.2711°W, and is a part of the measurement transect in Fig. 3c. Figures 3a and 3b exhibit great undulations in the raw velocities, revealing that the TRIAXUS motion relative to the vessel contaminates the velocity observations.

We corrected the TRIAXUS ADCP-measured velocities using the method introduced here. The results are shown in Fig. 4. Figures 4a and 4b give the ship-mounted ADCP-measured velocity *u* component and the corrected TRIAXUS ADCP-measured *u*, respectively. Figures 4c and 4d specify the *υ* components of the ship-mounted ADCP and corrected TRIAXUS-measured velocities, respectively. In the ship transect between 47°N, 125.2741°W and 46.9986°N, 125.2030°W (Figs. 4a and 4c on Julian day of 2005 between 150.5879 and 150.6079), one can see that a northeastward current occupied the surface layer with a depth of 5 m. Between 5 and 8 m, this current turned eastward. This feature is also shown in the TRIAXUS-measured velocities in Figs. 4b and 4d. In Fig. 4a, there was a westward current appearing in the ship transect between 46.9970°N, 125.1436°W and 47.0186°N, 125.0885°W (Fig. 4a on Julian day of 2005 between 150.6216 and 150.6391) below 5 m, which was also captured by the TRIAXUS ADCP (Fig. 4b). At the 10-m depth in the same position, Fig. 4c shows a southward current that also appeared in Fig. 4d. In addition, Figs. 4b and 4d indicate that there was a southwestward current occupying a depth between 30 and 40 m in this subtransect. In Fig. 4b, one can see a positive velocity region, which reveals that an eastward current existed in the ship transect between 47.0351°N, 125.0610°W and 47.0349°N, 125.0166°W (on Julian day of 2005 between 150.6489 and 150.6607), and strengthened at a depth of 50–60 m. These results show that the TRIAXUS ADCP measurements can detect deeper current structures than the ship-mounted ADCP.

Comparisons between the corrected TRIAXUS ADCP-measured velocities and the ship-mounted ADCP measurements are given in Fig. 5. Figures 5a and 5b show the *u* and *υ* components of the ship-mounted and TRIAXUS ADCP-measured velocities at the 8-m depth, respectively. The rms of the velocity difference between the two ADCP measurements is 0.024 m s^{−1} for *u* and 0.025 m s^{−1} for *υ*. Figures 5c and 5d, 5e and 5f and 5g and 5h show the time series of the velocity *u* and *υ* components at depths of 12, 14, and 16 m, respectively, whereas at these depths, the rms difference between two ADCP measurements for *u* and *υ* is calculated as 0.026 and 0.023 m s^{−1}, 0.042 and 0.025 m s^{−1}, and 0.042 and 0.023 m s^{−1}, respectively.

The TRIAXUS vehicle velocities relative to the vessel are derived and shown in Fig. 6. Figure 6a illustrates the *u* component and Fig. 6c displays the *υ* component. One can see that the velocities oscillate between −0.8 and 0.8 m s^{−1} in a certain frequency. This oscillation is consistent with the TRIAXUS depth undulation cycle. We calculate the power spectra of the TRIAXUS motion velocities, which are shown in Figs. 6b (*u* component) and 6d (*υ* component). There are three spectrum peaks in the figures; the largest is at the frequency 0.0068 Hz (period: 2.45 min), the second is at 0.0132 Hz (period: 1.26 min), and the third is at 0.0196 Hz (period: 51 s). Obviously, the period 51 s corresponds to the TRIAXUS vehicle depth oscillation cycle. The second and third could be the harmonics of the low-frequency peak because the TRIAXUS is not perfectly sinusoidal.

## 4. Discussion and summary

In Visbeck’s method (VM), a least squares equation was developed by using all the ADCP bins including the overlapping ones. To solve the barotropic velocity profile from the least squares equation, he used three constraints: zero current assumption, bottom tracking, and ship navigation. These constraints could construct an invertible model matrix ( ) and make the least squares equation solvable. To recover all depths of the ADCP velocity, we also use overlapping bins as a connection between neighboring ensembles. However, the three constraints proposed by Visbeck (2002) are not available for our cases. Thus, we use the ship-mounted ADCP velocity measurements for constraints rather than the zero current assumption, bottom tracking, and ship navigation. This is like building a new model matrix for the VM. The method in this note does not need the matrix inversion, but uses simple, step-by-step calculations to obtain the velocity profile, and therefore takes less computer time. In addition, this method uses the ship ADCP measurements as the references, which allows the variations of ocean current velocities between ensembles. In general, this method provides an alternative way to recover the undulated ADCP velocity profiles.

The TRIAXUS towed vehicle can oscillate in a water depth and make the onboard high-frequency ADCP acquire a deeper velocity profile. In this note, we introduce a method to remove the TRIAXUS ADCP vehicle velocities from the measurements and obtain the TRIAXUS vehicle motions relative to the ship in the earth coordinates. The results show that the corrected TRIAXUS ADCP-measured velocities could capture not only surface layer current features as measured by a ship-mounted ADCP but also the features in a deep depth up to 60 m. The comparison between the two ADCP measurements shows that rms differences between the velocities measured by ship-mounted and TRIAXUS ADCPs at depths of 8, 12, 14, and 16 m are less than 0.042 m s^{−1}. An analysis of the derived TRIAXUS vehicle motion velocities reveals that the vehicle horizontal velocities relative to the vessel have an oscillation consistent with the vehicle depth undulation.

## Acknowledgments

The study is supported by the National Science Foundation RISE (River Influences on Ecosystems) Project OCE 0239072. We thank Captain Ron L. Short of the R/V *Point Sur* and Marine Technicians Stewart Lamberdin and Ben Jokinen for their superb support of in situ data collection. The authors are also grateful to the anonymous reviewers for their valuable suggestions and comments.

## REFERENCES

Hickey, B. M., Pietrafesa L. J. , Jay D. A. , and Boicourt W. C. , 1998: The Columbia River plume study, subtidal variability in the velocity and salinity fields.

,*J. Geophys. Res.***103****,**10339–10368.MacArtney Underwater Technology Group, 2005: TRIAXUS towed undulator. Doc. 2005-001E. [Available online at http://www.triaxus.com/filer/1021/TRIAXUS-Towed-Undulator.625266203754055.pdf.].

Moum, J. N., Farmer D. M. , Smyth W. D. , Armi L. , and Vagle S. , 2003: Structure and generation of turbulence at interfaces strained by internal solitary waves propagating shoreward over the continental shelf.

,*J. Phys. Oceanogr.***33****,**2093–2112.Visbeck, M., 2002: Deep velocity profiling using lowered acoustic Doppler current profiles: Bottom tracking and inverse solution.

,*J. Atmos. Oceanic Technol.***19****,**794–807.

Fig. 2. Schematic diagram of the TRIAXUS and ship-mounted ADCP measurements. In fact, the ship-mounted and TRIAXUS ADCPs had forty-six 0.5-m bins.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 2. Schematic diagram of the TRIAXUS and ship-mounted ADCP measurements. In fact, the ship-mounted and TRIAXUS ADCPs had forty-six 0.5-m bins.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 2. Schematic diagram of the TRIAXUS and ship-mounted ADCP measurements. In fact, the ship-mounted and TRIAXUS ADCPs had forty-six 0.5-m bins.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 3. Uncorrected TRIAXUS ADCP-measured velocities, (a) *u* component and (b) *υ* component, in a (c) segment of a cruise transect between 46.9990°N, 125.3376°W and 47.0000°N, 125.2711°W.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 3. Uncorrected TRIAXUS ADCP-measured velocities, (a) *u* component and (b) *υ* component, in a (c) segment of a cruise transect between 46.9990°N, 125.3376°W and 47.0000°N, 125.2711°W.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 3. Uncorrected TRIAXUS ADCP-measured velocities, (a) *u* component and (b) *υ* component, in a (c) segment of a cruise transect between 46.9990°N, 125.3376°W and 47.0000°N, 125.2711°W.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 4. The ship-mounted ADCP-measured velocities, (a) *u* and (c) *υ* components, and the corrected TRIAXUS ADCP-measured velocities, (b) *u* and (d) *υ* components.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 4. The ship-mounted ADCP-measured velocities, (a) *u* and (c) *υ* components, and the corrected TRIAXUS ADCP-measured velocities, (b) *u* and (d) *υ* components.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 4. The ship-mounted ADCP-measured velocities, (a) *u* and (c) *υ* components, and the corrected TRIAXUS ADCP-measured velocities, (b) *u* and (d) *υ* components.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 5. Comparisons between the ship-mounted ADCP-measured velocities and the corrected TRIAXUS ADCP-measured velocities at depths of (a),(b) 8, (c),(d) 12, (e),(f) 14, and (g),(h) 16 m. The (left) *u* component and (right) *υ* component.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 5. Comparisons between the ship-mounted ADCP-measured velocities and the corrected TRIAXUS ADCP-measured velocities at depths of (a),(b) 8, (c),(d) 12, (e),(f) 14, and (g),(h) 16 m. The (left) *u* component and (right) *υ* component.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 5. Comparisons between the ship-mounted ADCP-measured velocities and the corrected TRIAXUS ADCP-measured velocities at depths of (a),(b) 8, (c),(d) 12, (e),(f) 14, and (g),(h) 16 m. The (left) *u* component and (right) *υ* component.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 6. The (a) *u* and (b) *υ* components of the TRIAXUS vehicle velocities and the power spectra of the velocity (c) *u* and (d) *υ* components.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 6. The (a) *u* and (b) *υ* components of the TRIAXUS vehicle velocities and the power spectra of the velocity (c) *u* and (d) *υ* components.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1

Fig. 6. The (a) *u* and (b) *υ* components of the TRIAXUS vehicle velocities and the power spectra of the velocity (c) *u* and (d) *υ* components.

Citation: Journal of Atmospheric and Oceanic Technology 25, 9; 10.1175/2008JTECHO534.1