Systematic microstructure variability in double-diffusively stable coastal waters of nonuniform density gradient

Author Posting. © American Geophysical Union, 2002. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 107, C10 (2002): 3144, doi:10.1029/2001JC000844.Conductivity microstructure, water velocity, and stratification were measured during a tow-yo transect near the New England shelf/slope front in early August 1997. Velocity data were collected with an acoustic Doppler profiler on the ship. The other data were collected with a towed platform. Estimates of χ, the rate of dissipation of temperature variance, were computed from the conductivity data with vertical resolution of 0.3 m. Relationships between χ and shear, temperature gradient, buoyancy frequency (N), and gradient Richardson number (Ri) were explored, with special focus on measurements taken in waters stable to double-diffusive processes (to avoid ambiguity of interpretation) and exhibiting variable density gradient (N ranging from 5 to 40 cph). For this subset of data, χ computed from data grouped into classes of local mean temperature gradient (inline equation/dz) was proportional to inline equation/dz to the 0.7 power, which is consistent with diapycnal thermal eddy diffusivity K being proportional to (inline equation/dz)−1.3 within the framework of the Osborn-Cox model that relates χ and the mean temperature gradient to the heat flux. No correlation between K and Ri was observed, with Ri computed at 4 m vertical scale, so that systematic inhomogeneous large-scale forcing is not responsible for a false correlation of K and inline equation/dz. Water mass salinity characteristics in the area caused N2 to be proportional to (inline equation/dz)4/5, rather than to inline equation/dz as in the isohaline case, giving rise to the steep inverse relation K = 10−10N−3.3, with N in radians/s and diffusivity in m2/s. The fit K = 2 × 10−9N−2.5 results if one questionable data ensemble is disregarded. These relations are comparable to results obtained previously from the near-bottom tow data. They are not intended to be universal formulae but are meant to describe the conditions we encountered. They are not expected to hold at high and low values of N outside of our measurement range. An interpretation is that under these conditions the less strongly stratified (lower N) layers in this shelf area are more prone to instability of the larger-scale shear than the intervening interfaces, with the subsequent greater energy dissipation in the layers leading to higher buoyancy flux KN2 in the layers than in the interfaces.This work was supported by the U.S. Office of Naval Research (Grants N-00014-95-1-0633 and N00014-95-1-1064) and by a WHOI postdoctoral scholarship award to C.R.R