Neutral density revisited

The possibilities of defining and computing an approximately neutral densit y variable are reexamined in this paper. There are three desirable properti es that a neutral density variable should possess. Firstly, the isosurfaces of this variable should coincide with (approximately) neutral surfaces. Th is would facilitate the analysis of hydrographic data on the most appropria te mixing and spreading surfaces. Secondly, the horizontal gradients of the neutral density should agree with the gradients of the in situ density, an d thirdly the vertical gradient of the neutral density variable should be p roportional to the static stability of the water column. A density variable that approximates the latter two properties can be used in ocean circulati on models based on layer coordinates, and would reduce substantial errors i n present isopycnal models due to the use of a potential density variable. No variable can possess all the three properties simultaneously. The variab le gamma(n) introduced by Jackett and McDougall (1997: J. Phys. Oceanogr. 2 7, 237-263) satisfies the first of the properties exactly but is not design ed for the use in models. Based on climatological data in the North Atlanti c, an alternative neutral density variable <(nu)over tilde>(S, Theta) is de fined, which is shown to approximate the two gradient criteria much better than any potential density. We suggest that this neutral density variable m ay be useful in isopycnal ocean models as an alternative to potential densi ty, since it could significantly reduce errors in thermal wind relation and vertical stability. (C) 1999 Elsevier Science Ltd. All rights reserved.