Neutral Density

Mathematical expression

A neutral density surface is the surface along which a given water mass will move, remaining neutrally buoyant.

McDougall and Jackett demonstrated that the normal to the neutral surfaces is in the direction of β ∇ S − α ∇ θ , where S is the salinity, θ is the potential temperature, α the thermal expansion coefficient and β the saline concentration coefficient. Thus, neutral surfaces are defined as the surfaces everywhere perpendicular to the vector ρ ( β ∇ S − α ∇ θ ) . For such a surface to exist, its helicity H must be zero; if this condition is respected, a scalar γ n exists and it is the one which satisfies the following formula:

where b is an integrating scalar factor, which is function of space, and R a residual term that is usually different from zero.

This formula represents a coupled system of first-order partial differential equations, that has to be solved to obtain the desired value of γ n . The solutions of (1 ) can be obtained by using numerical techniques.

In the real ocean, the condition of helicity equal to zero is not generally satisfied exactly. Therefore, and because of the non-linear terms in the equation of state, it is impossible to create analytically a Well-defined neutral density surface. There will always be flow through the calculated surfaces, because of the presence of a neutral helicity.

Therefore, it is possible to obtain only a best-fit approximate neutral surface, through which there is no flow of major proportions and along which it is generally accepted that flow takes place. γ n is a Well-defined function and Jackett and McDougall demonstrated that the inaccuracy due to the not exact neutrality is below the present instrumentation error in density. Neutral density surfaces stay within a few tens meters of an ideal surface anywhere in the world.

For how γ n has been defined, neutral density surfaces can be considered the continuous analog of the commonly used potential density surfaces, which are defined over various discrete values of pressures (see for example and ).

Spatial dependence

Given the spatial dependence of the neutral density, its calculation requires the knowledge of the spatial distribution of temperature and salinity in the ocean. Therefore, the definition of γ n has to be linked with a global hydrographic dataset, based on the climatology of the world’s ocean (see World Ocean Atlas and ). In this way, the solution of (1 ) provides values of γ n for a referenced global dataset. The solution of the system for a high resolution dataset would be computationally very expensive. In this case, the original dataset can be sub-sampled and (1 ) can be solved over a more limited set of data.

Algorithm for the computation of neutral surfaces using γ n {\displaystyle \gamma ^{n}\,}

Jackett and McDougall constructed the variable γ n using the data in the “Levitus dataset”. As this dataset consist of measurements of S and T at 33 standard depth levels at a 1° resolution, the solution of (1 ) for such a large dataset would be computationally very expensive. Therefore, they sub-sampled the data of the original dataset onto a 4°x4° grid and solved (1 ) on the nodes of this grid. The authors suggested to solve this system by using a combination of the method of characteristics in nearly 85% of the ocean (the characteristic surfaces of (1 ) are neutral surfaces along which γ n is constant) and the finite differences method in the remaining 15%. The output of these calculations is a global dataset labeled with values of γ n . The field of γ n values resulting from the solution of the differential system (1 ) satisfies (1 ) an order of magnitude better (on average) than the present instrumentation error in density.

The labeled dataset is then used to assign γ n values to any arbitrary hydrographic data at new locations, where values are measured as a function of depth by interpolation to the four closest points in the Levitus atlas.

Practical computation of γ n {\displaystyle \gamma ^{n}\,}

The formation of neutral density surfaces from a given hydrographic observation requires only a call to a computational code that contains the algorithm developed by Jackett and McDougall.

The Neutral Density code comes as a package of Matlab or as a Fortran routine. It enables the user to fit neutral density surfaces to arbitrary hydrographic data and just 2 MBytes of storage are required to obtain an accurately pre-labelled world ocean.

Then, the code permits to interpolate the labeled data in terms of spatial location and hydrography. By taking a weighted average of the four closest casts from the labeled data set, it enables to assign γ n values to any arbitrary hydrographic data.

Another function provided in the code, given a vertical profile of labeled data and γ n surfaces, finds the positions of the specified γ n surfaces within the water column, together with error bars.

The complete code is available through the World Wide Web at http://www.teos-10.org/preteos10_software/ . The code comes with documentation in the form of Readme files.

Advantages of using the neutral density variable

Comparisons between the approximated neutral surfaces obtained by using the variable γ n and the previous commonly used methods to obtain discretely referenced neutral surfaces (see for example Reid (1994), that proposed to approximate neutral surfaces by a linked sequence of potential density surfaces referred to a discrete set of reference pressures) have shown an improvement of accuracy (by a factor of about 5) and an easier and computationally less expensive algorithm to form neutral surfaces. A neutral surface defined using γ n differs only slightly from an ideal neutral surface. In fact, if a parcel moves around a gyre on the neutral surface and returns to its starting location, its depth at the end will differ by around 10m from the depth at the start. If potential density surfaces are used, the difference can be hundreds of meters, a far larger error.