A 2019 comment by Hochet and Tailleux and the corresponding reply by Holmes et al. discuss the volume and mass balance on a control volume bounded by a given isotherm and the ocean free surface. This note partly reconciles the terms of discrepancy on volume or mass transport appearing in these publications by proving that the integral expressions in the comment of Hochet and Tailleux, using a particular parameterization of the moving surfaces in Cartesian coordinates, correspond to the mass transport across the moving surfaces, as long as the mass density is included, as given in direct vector notation by Holmes et al. in their reply.
Two recent notes on the article by Holmes et al. (2019a, hereinafter HZE), namely, the comment by Hochet and Tailleux (2019, hereinafter HT) and the reply by Holmes et al. (2019b, hereinafter HZEr), discuss the volume and mass balance on a control volume bounded by a given isotherm and the ocean free surface. The main purpose of HT was to provide explicit mathematical expressions linking the volume change of HZE to the diabatic sources and sinks of heat. HZE correctly point out that the control volume is neither steady nor material and express the volume balance in terms of the volume flux across the appropriate moving boundary surfaces in direct vector notation [Eq. (3) in HZE complemented with the usual definition of volume flux across a moving surface given in Eq. (1) in HZEr]. On the other hand HT express the mass balance in terms of horizontal area integrals in Cartesian coordinates involving the fluid flow and the geometry of the pertinent surfaces [Eq. (15) in HT]. Since the relation between both mathematical expressions is not given, the reader may remain confused about which approach is in fact correct. This note reconciles both approaches by showing that, as long as the volume flux in HZEr is used instead of the mass flux, both mathematical expressions are exactly equal.
2. Mass transport across moving surfaces
For managing balance equations in integral form in arbitrary control volumes (i.e., not necessarily material or steady), it is useful to consider the kinematic identity (see Viúdez 2000) expressing the rate of change of mass in an arbitrary volume ,
The above expression means that the rate of change of the extensive quantity with spatial density f(x, t) in an arbitrary volume , which is material with respect to the velocity w(x, t), is equal to the rate of change of the extensive quantity in the fluid particles moving with velocity u(x, t) located in the arbitrary volume minus the flux or transport of f relative to the moving boundary ∂. In the particular case f = 1 Eq. (1) represents the volume balance, and in the case that f is equal to the mass density f(x, t) = ρ(x, t) Eq. (1) represents the mass balance
where we have assumed mass conservation ρ + ρ∇ ⋅ u = 0. The term in the right-hand side of Eq. (2) is equivalent to the term defined by HZEr [their Eq. (1)] when mass transport is considered instead of volume transport. We prove now that this balance is equal to the first equality in HT [their Eq. (15)], which for clarity is reproduced here as
For the upper free surface z = η(x, y, t), we define the function s(x, y, z, t) ≡ z − η(x, y, t) so that the implicit equation for the free surface is s(x, y, z, t) = 0. The velocity of displacement ws(x, y, t) of this free surface satisfies the equation
where the unit vector normal to the free surface
so that, at the upper surface, the right-hand side of Eq. (4) includes the term
Note that, since w is the velocity of displacement of the integration surface, w and dS are tangent vectors.
In a similar way, for the lower, bottom surface, z = −h(x, y, θ, t), and for a constant temperature value θ we define the function b(x, y, z, t) ≡ −z − h(x, y, θ, t) so that the implicit equation for the isothermal surface z = −h(x, y, θ, t) is b(x, y, z, t) = 0. Following the same steps carried out above with the free surface we obtain the velocity of displacement wb(x, y, t) of the isothermal surface:
where the unit normal vector
and therefore at the bottom isothermal surface
Using the above relations for the speed of displacement w of a moving surface defined by a relation of the form F(x, y, z, t) ≡ z − f(x, y, t) = 0, and for a differential area element dS on this surface, it follows that
where u[x, y, f(x, y, t)] = uf(x, y, t) = ufh + wfk is the fluid velocity, and ρ[x, y, f(x, y, t)] = ρf(x, y, t) is the mass density, on the surface z = f(x, y, t). Thus, the term on the right-hand side of Eq. (14) may be kinematically interpreted as the differential flux of mass across the moving surface z = f(x, y, t).
The general theory of mass and volume transport across isosurfaces of a balanced fluid property, that is, a fluid property obeying a balance equation, may be used to derive other balance configurations in the ocean, such as that in Walin (1982) for the heat balance or that in Garrett et al. (1995) for the average buoyancy balance (see Viúdez 2000, his appendix B).
In this note we have partly reconciled the terms of discrepancy on volume or mass transport appearing in two recent notes by proving that the integral expressions in Hochet and Tailleux (2019), using a particular parameterization of the moving surfaces in Cartesian coordinates, correspond to the mass transport across the moving surfaces, as long as the mass density is included, as given in direct vector notation by Holmes et al. (2019b).
I thank one anonymous reviewer for his/her comments. Partial support for this study was obtained through project RTI2018-100844-B-C33 (Spanish Ministry of Science, Innovation and Universities).