After the field project, all 19 sonics were operated in the 89-cm-diameter, low-turbulence wind tunnel of the NCAR Sensor Calibration Laboratory. The working section of this wind tunnel is too small to accurately calibrate the sonics, but the wind tunnel is sufficiently repeatable to examine the consistency of the sonics. This was done by operating the wind tunnel over the range from 0 to 20 m/sec and calculating the linear regression of the wind components measured by each sonic with respect to the wind tunnel speed measured by a pitot-static tube. The orthogonal wind components were transformed into path-parallel wind components using the nominal angles of the paths:

`U _{a} = (-U + sqrt(3)*W)/2
U_{b} = (U + sqrt(3)*V)/4 + sqrt(3)*W/2
U_{c} = (U - sqrt(3)*V)/4 + sqrt(3)*W/2 `

The following tables list the intercept `a` and slope `b` of a linear regression of the sonic wind components onto the tunnel speed measured by the pitot tube, e.g.

`U _{a} = a + b*U_{pitot} `

Sonic Position | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
---|---|---|---|---|---|---|---|---|---|

CSAT3 S/N | 246 | 247 | 364 | 366 | 367 | 368 | 369 | 370 | 373 |

U_{a}, a (cm/s) |
1.8 | 0.6 | 8.4 | -0.3 | 1.5 | 1.2 | -0.7 | 0.1 | 3.8 |

U_{a}, b |
-0.472 | -0.470 | -0.464 | -0.471 | -0.478 | -0.471 | -0.471 | -0.475 | -0.469 |

U_{b}, a (cm/s) |
3.2 | 1.8 | 2.3 | 1.7 | 0.4 | -1.0 | 0.2 | 0.6 | 1.3 |

U_{b}, b |
0.267 | 0.264 | 0.252 | 0.264 | 0.261 | 0.267 | 0.260 | 0.262 | 0.259 |

U_{c}, a (cm/s) |
-0.3 | 1.0 | 0.2 | 0.4 | 3.1 | 0.2 | 1.6 | 3.2 | 3.5 |

U_{c}, b |
0.253 | 0.269 | 0.277 | 0.263 | 0.274 | 0.272 | 0.267 | 0.274 | 0.268 |

Sonic Position | 1 | 2 | 3 | 4 | 5 | 6 | Spare | 7 | 8 | 9 |
---|---|---|---|---|---|---|---|---|---|---|

CSAT3 S/N | 374 | 375 | 376 | 377 | 378 | 379 | 800 | 853 | 855 | 856 |

U_{a}, a (cm/s) |
-0.5 | 8.2 | 10.6 | -2.0 | 4.4 | -1.2 | -0.1 | -1.7 | -1.2 | -5.2 |

U_{a}, b |
-0.479 | -0.499 | -0.473 | -0.468 | -0.468 | -0.470 | -0.466 | -0.462 | -0.463 | -0.475 |

U_{b}, a (cm/s) |
2.8 | -7.9 | 0.6 | -1.9 | 1.7 | -4.2 | 1.8 | 1.7 | 5.6 | 0.7 |

U_{b}, b |
0.264 | 0.267 | 0.277 | 0.269 | 0.263 | 0.266 | 0.264 | 0.270 | 0.271 | 0.260 |

U_{c}, a (cm/s) |
2.8 | 2.2 | 6.1 | 2.5 | 2.9 | 1.4 | -0.1 | 3.9 | 0.6 | -13.7 |

U_{c}, b |
0.271 | 0.257 | 0.270 | 0.252 | 0.261 | 0.267 | 0.263 | 0.264 | 0.260 | 0.253 |

If the sonics were perfectly aligned in the wind tunnel and there was no flow distortion caused by the sonic, then the values of the slopes should be 0.5, 0.25 and 0.25. The calculated slopes for `U _{a}` are slightly smaller and those for

Most troubling, however, are the intercepts for sonics 364 (`U _{a}`), 375 (

The calibration by Campbell has been completed and the following results received from Ed Swiatek (2/8/05):

s/n 0375 - The head was broken at one of the 45 degree joints. This joint was repaired once before in the past. We are not sure if the join simply failed again or the head was mishandled. The head was repaired and recalibrated using embedded code 3.0 and 4.0.

s/n 0856 - There was a significant amount of water inside the transducers. This water changed the response of the transducers. The response further changed when the water froze. The water was "baked" out of all the transducers. This CSAT3 was never calibrated using embedded code 3.0, thus it was only calibrated using 4.0. (TWH: This was our spare sonic and to the best recollection of Jim Edson, Steve Semmer and me it never left the laboratory at WHOI. It was exposed to rain during the pre-project flow distortion test, but to the best of my knowledge this was its only use outdoors.)

s/n 0364 - It appears that the calibration drifted on this CSAT3. Larry Jacobsen is analyzing the before and after data in an attempt to understand why that calibration drifted. This CSAT3 was calibrated with embedded code 3.0 and 4.0.

s/n 0376 - It appears that the calibration drifted on this CSAT3 as well. Larry Jacobsen is also analyzing the before and after data in an attempt to understand why that calibration drifted. This CSAT3 was calibrated with embedded code 3.0 and 4.0.

*J. Edson*

Figure 1. OHATS sonic array, 8/6/04 (sonics at each height are numbered 1-9 from left to right).

*J. Edson*

Figure 2. OHATS sonic array, 10/26/04.

During previous ISFF projects, sonic data have been corrected for tilt of the anemometer using the planar fit technique, which determines the sonic pitch and roll angles by fitting a plane to the measured vertical velocity as a function of wind direction and assuming that deviations from a horizontal plane are caused by sonic tilt. The fit is applied for a restricted range of wind directions that does not include directions downwind of known obstacles such as the tower on which the sonic is mounted or, in the case of the horizontal array technique used in OHATS, directions downwind of adjacent sonics. However even for a narrow range of wind directions, say within ± 45° of normal to the horizontal sonic array, it is evident that there is flow distortion upwind of the ASIT tower which causes a true non-zero vertical velocity.

Figure 3 shows an example planar fit for the center sonic at a height of 5 m. The 30-minute-averaged elevation angle of the wind is plotted versus wind direction, where zero azimuth angle is a wind normal to the array. A tilted sonic in a planar wind field (zero 30-minute-average true vertical velocity) maps out a sinusoidal curve as shown on the plot. The distorted flow field upwind of ASIT does not fall on the best-fit sine curve. This occurs for |azimuth| > 90° because the sonics are then downwind of other sonics and the ASIT legs. More seriously, however, the flow is also distorted upwind of the ASIT tower, i.e. for -45° < azimuth < 45°.

*T. Horst*

Figure 4. MotionPak (small gray box) mounted on upper (5.6 m) sonic array.

As a consequence of the flow distortion upwind of the ASIT tower, a less objective method is used to rotate the sonic data to a nominally vertical orientation. First the MotionPak data are used to (objectively) remove the pitch and roll angles of the sonic array support structure. The array pitch and roll angles are determined from the 5-minute-average accelerations measured by the MotionPak along three orthogonal axes,

`pitch = -sin ^{-1}(a_{x}/g) `

`roll = sin ^{-1}(a_{y}/g cos(pitch)) `

where the acceleration of gravity is calculated as

`g = (a _{x}^{2} + a_{y}^{2} + a_{z}^{2})^{1/2} `,

`a _{x}` is the horizontal acceleration normal to the array (positive in a direction from the array to the ASIT tower),

Figure 7 shows the pitch and roll angles of the array after correction for the roll offset. Note that beginning around noon on September 18, a wind of 10-15 m/s from the back of the array permanently reduced the pitch angle by about 0.5°. Another strong wind event from the opposite wind direction occurred in the third week of October.

Finally, the sonic pitch associated with sag of the sonic booms is determined by (subjectively) assuming that the vertical velocity at the location of each sonic is equal to zero for a wind direction normal to the array. These derived pitch angles range up to ± 2°, but 60% are within the range ± 1°. Although the assumption used in this final step is arbitrary, the pitch angles determined in this manner are qualitatively consistent with the relative pitch of adjacent sonics as seen in photographs of the array, Figures 1-2. (Note in particular that the relative pitch angles of the sonics appear to be unchanged over the deployment period.) These angles, shown in the following table, are *subtracted* from the pitch measured by the MotionPak, e.g. a sonic pitch offset of 2.25° (sonic 1 at 5 m) subtracted from an array pitch of -3° produces a net sonic pitch of -5.25°.

Sonic Position: | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
---|---|---|---|---|---|---|---|---|---|

Height=5.6 m | 0.05 | 0.8 | -0.2 | 0.5 | 0.2 | 1.3 | -2.0 /-1.35 | -0.85 | 0.7 /1.15 |

Height=5.0 m | 2.25 | 0.1 | 0.1 | 1.7 | -1.15 | 0.5 | 2.05 | 0.25 | 0.85 |

The individual tilt plots can be seen by means of the links in the table. The sonic position in the table is its ordinal location when facing into the array toward the ASIT tower and counting from left to right. The 7th and 9th sonics at 5.6 m appear to have changed their pitch during a high wind event from the back of the array near noon on September 18. The double entries for those sonics refer to the pitch angles prior to and following that event.

The planar fit technique is not applicable to the OHATS sonic data because of flow distortion upwind of the ASIT tower. As a consequence the absolute values of the pitch and roll angles of the sonics cannot be determined, although the relative pitch angles of the sonics are probably accurate to better than 1°. Because we do not know the relative alignment of the MotionPak and the sonics, the uncertainty in sonic pitch and roll will cause errors when attempting to remove the array movement from the sonic data. It may be that these errors are acceptable, but this remains to be determined.