A Rigorous Derivation of the Water Mass Transformation Framework, the Relation between Mixing and Diasurface Exchange Flow, and Links to Recent Theories in Estuarine Research

Knut Klingbeil aLeibniz Institute for Baltic Sea Research Warnemünde (IOW), Rostock, Germany

Search for other papers by Knut Klingbeil in
Current site
Google Scholar
PubMed
Close
https://orcid.org/0000-0002-6736-1260
and
Erika Henell aLeibniz Institute for Baltic Sea Research Warnemünde (IOW), Rostock, Germany

Search for other papers by Erika Henell in
Current site
Google Scholar
PubMed
Close
https://orcid.org/0000-0003-4381-3100
Free access

Abstract

In this paper we present the analytical derivation of a local water mass transformation (WMT) framework for an individual water column. We exactly formulate the mapping of the governing equations from geopotential coordinates to an arbitrary tracer space. Unique definitions for the local effective vertical diasurface fluxes are given. In tracer space we derive new relations between the local diatracer fluxes and the mixing per tracer class. The key relation between the effective vertical diatracer velocity and the mixing per tracer class directly formulates how the overturning circulation is linked to local tracer variance dissipation. Horizontal integration of the governing equations in tracer space and the relations between the diatracer quantities finally recovers the well-known integral WMT formulations.

© 2023 American Meteorological Society. This published article is licensed under the terms of the default AMS reuse license. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Knut Klingbeil, knut.klingbeil@io-warnemuende.de

Abstract

In this paper we present the analytical derivation of a local water mass transformation (WMT) framework for an individual water column. We exactly formulate the mapping of the governing equations from geopotential coordinates to an arbitrary tracer space. Unique definitions for the local effective vertical diasurface fluxes are given. In tracer space we derive new relations between the local diatracer fluxes and the mixing per tracer class. The key relation between the effective vertical diatracer velocity and the mixing per tracer class directly formulates how the overturning circulation is linked to local tracer variance dissipation. Horizontal integration of the governing equations in tracer space and the relations between the diatracer quantities finally recovers the well-known integral WMT formulations.

© 2023 American Meteorological Society. This published article is licensed under the terms of the default AMS reuse license. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Knut Klingbeil, knut.klingbeil@io-warnemuende.de

1. Introduction

In his seminal work, Walin (1977) presented integral budget equations for the volume V>(S) and salt content S>(S) of all water within a semienclosed basin with salinities larger than a prescribed salinity S,
V>(S)tQ>(S)Qdia(S)=0,
S>(S)tF>(S)Fdia(S)=0.
The volume V>(S) and its salt content S>(S) are illustrated in Fig. 1. The reference to the original equations in Walin (1977), and a comparison to the original notation, is given in Table 1. The integrated volume and salt fluxes entering across the open ocean boundary with salinities larger than S are given by Q>(S) and F>(S). Across the isohaline the integrated diahaline fluxes Qdia(S) and Fdia(S) are entering, with the diahaline salt flux consisting of advective (SQdia) and diffusive contributions (Jdia):
Fdia(S)=SQdia(S)+Jdia(S).
Walin (1977) showed that for no diffusive fluxes across the open boundary, i.e., F>(S)=SSq(S)dS with q(S) = −∂Q>/∂S being the volume flux per salinity class, the integrated diahaline volume flux is determined by the integrated diahaline diffusive salt flux
Qdia(S)=Jdia(S)S.
In analogy, Walin (1982) presented budgets for volumes V(T) of all water with temperatures smaller or equal to T, and additionally considered integrated volume and temperature fluxes across the ocean surface, Qη,≤(T) and Fη,(T), respectively. For Fη,(T)=TTqη(T)dT+Jη,(T), with qη(T) = ∂Qη,≤/∂T being the volume flux per temperature class, and Jη,(T) being the integrated nonadvective temperature flux across the ocean surface, relation (3) was extended to
Qdia(T)=Jdia(T)T+Jη,(T)T.
Walin’s concept of analyzing budgets of water masses defined by a distinct property, e.g., salinity or temperature, is the foundation of the water mass transformation (WMT) framework. The transformation of water masses due to interior and boundary processes leads to mass fluxes (volume fluxes under the Boussinesq approximation) between water masses, i.e., across the isosurface of the defining property, and is measured by Qdia, see (3) and (4). A review of the integral WMT framework, its derivation from differential volumes in tracer space and references for the numerous applications can be found in Groeskamp et al. (2019).
Fig. 1.
Fig. 1.

Idealized illustration of control volumes bounded by an isohaline with salinity S and a boundary transect toward the open ocean. For simplicity a stable salinity stratification is assumed, such that the volume V>(S) containing water with salinities larger than S is located below the isohaline, and the remaining volume V(S) containing water with salinities less or equal to S above the isohaline. The integrated diahaline volume flux Qdia is a measure for the water mass transformation. Integrated surface boundary and open boundary volume fluxes are depicted as Qη,≤, Q, and Q>, defined based on how the salinity of the transported water relates to S. Integrated salt fluxes are denoted by F and alter the salt contents S(S) and S>(S) inside the corresponding volumes. Surface boundary fluxes Qη,> and Fη,>, transporting or exchanging water with salinities larger than S, are not sketched for simplicity. Additional freshwater river discharge is depicted as Qr. The arrows indicate the direction of positive fluxes. Walin (1977) presented the budget equations for V>(S) and S>(S); see (1a) and (1b). In many studies on estuaries the budgets for V(S) and S(S), (5a) and (5b) are used.

Citation: Journal of Physical Oceanography 53, 12; 10.1175/JPO-D-23-0130.1

Table 1.

Overview of main quantities in an isohaline framework with salinities larger than the specific isohaline, and comparison to earlier publications. Wang et al. (2017) and Walin (1977) presented the equivalent equations horizontally integrated over a single water column and an ocean basin with vertical transects, respectively. In the integrated equations relation (23) holds and there is no need to distinguish between diahaline and effective vertical diahaline fluxes.

Table 1.
With the focus on integrated diahaline fluxes in classical estuaries, MacCready et al. (2002) presented the budgets for V(S) and S(S), neglecting surface boundary fluxes but including explicit freshwater river discharge Qr,
V(S)tQ(S)+Qdia(S)=Qr,
S(S)tF(S)+Fdia(S)=0.
The budgets (5a) and (5b) offer to study the water masses inside a classical estuary. The special case of an estuarine volume completely bounded by an isohaline, i.e., without open boundary fluxes (Q,F), was considered by MacCready and Geyer (2001) and Hetland (2005). Contrary, if the estuarine volume is completely bounded by the open boundary, see, e.g., MacCready (2011) or Burchard et al. (2019), the diahaline fluxes (Qdia,Fdia) drop out. For this case, MacCready (2011) developed the concept of total exchange flow (TEF): based on the volume flux per salinity class q(S) = ∂Q/∂S = −∂Q>/∂S, the time-averaged open boundary fluxes can be decomposed into net in- and outflow contributions, and associated salinities for the in- and outflowing water masses can be calculated, finally recovering the Knudsen relations (Knudsen 1900; Burchard et al. 2018).
Burchard (2020) supplemented (5a) and (5b) by a budget equation for the content of squared salinity. Its sink term M(S)=V(S)χdV, being the volume-integrated dissipation of salinity variance χ (Burchard and Rennau 2008), is a measure for mixing inside V(S). With
𝕞(S)=MS(integratedmixingpersalinityclass),
Burchard (2020) derived the linear relation
𝕞(S)=2SQr
as the universal law of estuarine mixing for estuarine volumes completely bounded by the isohaline with salinity S, and assuming sufficiently long time averaging (denoted by angle brackets). Burchard et al. [2021, their (19)] showed the general relation between the integrated diahaline diffusive salt flux and the integrated mixing per salinity class
Jdia(S)=12𝕞(S).
Although the integrated diahaline fluxes and mixing (Qdia, Jdia, and 𝕞) provide valuable information about interior processes, they do not provide information on how these quantities are distributed on the isohaline within the basin. Separating the domain into individual water columns facilitates a detailed view of the horizontal distribution. Volume integrated WMT budget equations for individual water columns have been presented in, e.g., Wang et al. (2017) and Holmes et al. (2021). However, in those works no analytical derivation and definition of the local diahaline quantities was done. Under the assumption of a stable salinity stratification, Li et al. (2022) derived budget equations for the thickness D>(S) of all water with salinities larger than S as the distance of the isohaline from the bottom, and the corresponding salt content S>(S). They also introduced the local diahaline volume flux per unit horizontal area udia,z and the local mixing per salinity class m(S). Differences in notations are again compared in Tables 1 and 2. However, because of the conceptual limitation due to the assumed stable salinity stratification, their formulated relations are not generally valid.
Table 2.

Overview of main integrated quantities in an isohaline framework with salinities smaller or equal to the specific isohaline, and comparison to earlier publications.

Table 2.
In this manuscript we present the analytical derivation of a local WMT framework and generally valid mathematical definitions for the corresponding local quantities. In contrast to Groeskamp et al. (2019), we start our derivations in classical geopotential coordinates and exactly formulate the mapping to tracer space. Our new unified presentation easily links the various existing integrated formulations. One major new result of our derivations is an analytical equation that shows how the overturning circulation is linked to local small-scale tracer variance dissipation. Until now, this can only be anticipated from the integrated formulation
Qdia(S)=12𝕞(S)S,
that we obtain by combining (3) and (8).

The paper starts with the governing equations in Cartesian coordinates (section 2a) and the mathematical formulation for the mapping from vertical coordinates to a general coordinate space of an arbitrary property s in section 2b. The mapping of the local budgets for volume, tracer content and squared tracer content, as well as the definitions for the local diasurface quantities, are presented in the following sections 2c2e. The relations between these diasurface quantities in tracer coordinates are derived in section 2f. In section 2g we show how the local budgets and relations recover the well-known integral formulations presented in the introduction. The new local relations are demonstrated in section 3. A summary and conclusions are given in section 4. In the appendixes, details can be found about the unified notation (appendix A), the relation between the general mapping and coordinate transformations (appendix B), the mapping of the budget equations (appendix C), integral relations in tracer space (appendix D), and the used Leibniz rule (appendix E).

2. Derivation of a local WMT formulation

a. Governing equations in Cartesian coordinates

Under the Boussinesq approximation, mass conservation is turned into volume conservation which is locally described by the incompressibility constraint
ux+υy+wz=0.
In (10) the flux divergence is written in Cartesian coordinates x, y, z, and the velocity components in the corresponding directions are given by u, υ, and w, respectively. Kinematic boundary conditions at the ocean surface and the bottom are given by
ηt+u(η)ηx+υ(η)ηyw(η)=uη,z,
u(H)Hx+υ(H)Hy+w(H)=0,
with uη,z being the volume flux across the ocean surface per unit horizontal area (positive into the water column), e.g., due to precipitation, evaporation, ice melting, or freezing.
The local budget equation for an arbitrary tracer with concentration cz(x, y, z, t) reads
czt+x(uczKhczx=jx)=fx+y(υczKhczy=jy)=fy+z(wczKυczz=jz)=fz=π.
The notation cz is used to indicate the functional dependence on z coordinates. The tracer flux in the x direction is given by fx and consists of advective (ucz) and diffusive (jx) contributions. The fluxes into the other directions are defined accordingly. Diffusive fluxes are given in terms of horizontal and vertical diffusivities Kh and Kυ, respectively, but can also be extended to consider more advanced diffusivity tensors. The source term π on the right-hand side of (12) accounts for local tracer production, e.g., due to shortwave radiation for temperature. Flux boundary conditions at the ocean surface and bottom are given by
cz(η)ηt+fx(η)ηx+fy(η)ηyfz(η)=fη,z,
fx(H)Hx+fy(H)Hy+fz(H)=0,
with fη,z being the surface tracer flux per unit horizontal area (positive into the water column). Bottom tracer fluxes, e.g., geothermal heating, are neglected but can easily be added to the equations.

b. Mapping to an arbitrary s space

In this work we follow Walin’s construction of control volumes bounded by a distinct water mass property. But, in contrast to Walin, we will apply the construction to a water column. Let the function sz(z) describe the vertical distribution of an arbitrary property (e.g., temperature, salinity or density) within the water column between the bottom at z = −H(x, y) and the ocean surface at z = η(x, y, t). For any prescribed property value s we need to identify the parts of the water column, i.e., the vertical positions z, where water masses with sz(z) ≤ s or sz(z) > s are located. This is straight forward for strictly monotonic functions sz(z), e.g., for stably stratified temperature or salinity distributions as sketched in Figs. 2a and 2c. Only in this monotonic case there exists a unique inverse function zs(s), which directly provides the vertical position of the prescribed s isosurface that separates the two control “volumes” in the water column. This inverse function does not exist for nonmonotonic functions sz(z), e.g., for tracer distributions with inversions as sketched in Fig. 2e. In this nonmonotonic case the isosurface exists at multiple positions within the water column and the control volumes consist of noncontiguous parts. Following Wolfe (2014) and Klingbeil et al. (2019), a general mathematical formulation that correctly considers also these noncontiguous parts, is possible by means of the Heaviside step function H{x} and the filter functions
H{x}=1H{x}={1,forx0,0,forx>0,
H>{x}=H{x}={0,forx0,1,forx>0.
With these filter functions the cumulative thickness of the control volumes, i.e., all parts of the water column where sz(z) ≤ s or sz(z) > s, respectively, can be calculated as
D(s)=HηH{sz(z)s}dz,
D>(s)=HηH>{sz(z)s}dz.
Due to the filter functions in the integrands of (15a) and (15b), the integrals only count those parts of the water column where sz(z) ≤ s or sz(z) > s, respectively. Figure 2 illustrates the functioning of the filter functions, as well as the resulting cumulative thicknesses D and D> for strictly monotonic and nonmonotonic functions sz(z). For strictly monotonic functions sz(z), and depending on the sign of ∂sz/∂z, the cumulative thickness represents the distance of the isosurface to the ocean surface or bottom. For nonmonotonic functions the cumulative thickness consists of the thicknesses from the individual noncontiguous parts.
Fig. 2.
Fig. 2.

Control volumes within a water column separated by an isosurface of property s. The vertical distribution of the property is given by the function sz(z). (a),(c),(e) Examples for strictly monotonic and nonmonotonic distributions are illustrated. For a prescribed value s the filter functions H and H>, defined in (14a) and (14b), identify the parts of the water column with sz(z) ≤ s or sz(z) > s, respectively. (b),(d),(f) Examples for the location of the isosurface in xz space in the vicinity of the water column are sketched. The local volume budget (18) describes the relation between the cumulative thickness D and horizontal transport U, as well as the diasurface velocities udiaz for the two control volumes.

Citation: Journal of Physical Oceanography 53, 12; 10.1175/JPO-D-23-0130.1

In the following we will present the derivations only for the “≤” case for simplicity. The unified derivation for both cases can be found in the appendixes. The differential thickness, i.e., the thickness per s class, is given by
D(s)s=(15a)Hηδ{sz(z)s}dz,
and ∂D/∂sds represents the thickness of the infinitesimal ds band around the s isosurface. The Dirac delta distribution δ{x}=H{x}/x is the derivative of the Heaviside step function. It only lets vertical positions contribute to the integral, where the specific isosurface s is present. In analogy to D and ∂D/∂s, the cumulative and differential formulations for any quantity ψz can be obtained by
Ψ(s)=Hηψz(z)H{sz(z)s}dz,
Ψ(s)s=Hηψz(z)δ{sz(z)s}dz.
The cumulative and differential quantities in (17a) and (17b) are defined for any value of s. Therefore, these relations can be used to map data from z coordinates to an arbitrary s space, similar to a classical coordinate transformation, but applicable also for nonmonotonic functions sz(z). See appendix B for more details.

c. Local volume budget

Following the derivation of (C15) in appendix C, vertical integration of (10) over the parts of the water column where sz(z) ≤ s yields the cumulative volume budget in s space
D(s)t+U(s)x+V(s)y+udia,zs(s)=uη,zH{sηs},
with the cumulative thickness D from (15a), the cumulative transports
U(s)=Hηu(z)H{sz(z)s}dz,
V(s)=Hηυ(z)H{sz(z)s}dz,
and the effective vertical diasurface velocity (positive into the direction of increasing s)
udia,zs(s)=Hηudia,zz(z)|szz|δ{sz(z)s}dz(effectiveverticaldiasurfacevelocity).
The cumulative transports in (19a) and (19b) are the basis for generalized streamfunctions in s space (Döös and Webb 1994; Nurser and Lee 2004; MacCready 2011; Wolfe 2014). The effective vertical diasurface velocity udia,zs in s space, defined in (20), represents the volume flux across the given s isosurface per unit horizontal area.1 It is a key quantity in this paper. The term effective stresses that the integrand can contribute with compensating values from multiple positions in the water column where the specific isosurface s is present (see Fig. 2f). At each vertical position in the water column, the local diasurface volume flux per unit horizontal area, udia,zz in z space in the integrand of (20), is given by
udia,zz(z)=udiaz(z)|n(z)z|1=udiaz(z)|szz|1|sz|,
in terms of the diasurface velocity udiaz. The relation between udiaz and udia,zz in (21) is illustrated in Fig. 3. The diasurface velocity udiaz is the component of the relative velocity normal to the isosurface that is currently located at the specific position z and moving with velocity us,
udiaz(z)=[u(z)us(z)]n(z)=1|sz|DszDt(diasurfacevelocity),
with n(z) = ∇sz/|∇sz| and z denoting the unit normal vectors at the isosurface and in the vertical, respectively. The term D/Dt is the material derivative. As sketched by the geometric relations in Fig. 3, the diasurface transport across an infinitesimal isosurface area element dA can be reformulated in terms of the local diasurface volume flux per unit horizontal area udia,zz and the horizontal projection of the isosurface area dAz,
udiazdA=(21)udia,zz|zn|dA=dAz=udiaz=udia,zzdAz.
Relation (23) is well known from the treatment of boundary fluxes at the ocean surface and bottom, e.g., (C4a) and (C5a), because the mathematical formulation is identical for any surface (Griffies 2004, section 6.7).
Fig. 3.
Fig. 3.

Sketch for relation between area elements and diasurface velocities (udiaz,udia,zz). The infinitesimal area element dA is located on the isosurface with the local unit normal vector n. Its horizontal projection is given by dAz, with z being the vertical unit vector. By construction udiazdA=udia,zzdAz holds; see (23). The same relations are valid for the total and diffusive tracer fluxes (fdiaz,fdia,zz) and (jdiaz,jdia,zz), respectively.

Citation: Journal of Physical Oceanography 53, 12; 10.1175/JPO-D-23-0130.1

The local diasurface volume flux per unit horizontal area, udia,zz in (21), is formally not defined for vertical isosurface elements, i.e., for ∂sz/∂z = 0. In our framework for a vertical water column this does not pose a problem, because the formally infinite udia,zz is multiplied with |∂sz/∂z| = 0 in (20) and an associated horizontal area dAz = 0 in (23). Instead, the diasurface transport through a vertical isosurface is already accounted for as part of the cumulative horizontal transports (19a) and (19b). Moreover, due to the associated modified horizontal transport divergence in (18), this “missing” contribution to udia,zs in (20) will simply appear in an adjacent water column. Therefore, the globally integrated diasurface volume flux
Qdia(s)=udia,zs(s)dAz(integrateddiasurfacevolumeflux)
again correctly accounts for all contributions. In case of a nonmonotonic coordinate distribution, with an isosurface occurring at multiple positions within the water column, the area dAz is still only counted once for the specific s (and not for every position z) in (24), because the effect of multiple isosurfaces is already accounted for in the multiple contributions of udia,zz to udia,zs in (20).
For strictly monotonic coordinate distributions sz(z), and depending on the sign of ∂sz/∂z, the cumulative transports (19a) and (19b) represent the transports below or above the isosurface, respectively (see Figs. 2b,d). In this monotonic case a unique inverse function zs(s) exists, and the effective vertical diasurface velocity in (20) only contains the contribution from the one vertical position of the isosurface,
udia,zs(s)=(20),(B3)Hηudia,zz[zs(s)]δ{ss}ds=(B3)udia,zz[zs(s)]forstrictlymonotonicsz(z).
Finally, the term on the right hand side of (18) with uη,z considers the volume fluxes across the ocean surface into water column volumes that contain the isosurface sη. In case of geometric coordinates sη refers to the coordinate value at the ocean surface sz(η). In case of s being a water mass property sη refers to the property value that is transported with the volume flux.

d. Local tracer budget

Following the derivation of (C16) in appendix C, vertical integration of (12) over the parts of the water column where sz(z) ≤ s yields the cumulative tracer budget in s space
C(s)t+Fx,(s)x+Fy,(s)y+fdia,zs(s)=Π(s)+fη,zH{sηs},
with the cumulative tracer content
C(s)=Hηcz(z)H{sz(z)s}dz,
the cumulative tracer fluxes
Fx,(s)=Hηfx(z)H{sz(z)s}dz,
Fy,(s)=Hηfy(z)H{sz(z)s}dz,
the cumulative tracer production
Π(s)=Hηπ(z)H{sz(z)s}dz,
and the effective vertical diasurface tracer flux (positive into the direction of increasing s)
fdia,zs(s)=Hηfdia,zz(z)|szz|δ{sz(z)s}dz(effectiveverticaldiasurfacetracerflux).
For the interpretation of the effective vertical diasurface tracer flux, and the quantities in the following relations, we refer to the analogous explanations for the effective diasurface velocity udia,zs in (20). The local diasurface tracer flux per unit horizontal area in the integrand of (30) is defined as
fdia,zz(z)=fdiaz(z)|n(z)z|1=fdiaz(z)|szz|1|sz|,
with the diasurface tracer flux fdiaz being the component of the tracer flux normal to the moving isosurface that is currently located at the specific position z,
fdiaz(z)=cz(z)udiaz(z)+j(z)n(z)=jdiaz(z)=(22)1|sz|(czDszDt+jxszx+jyszy+jzszz)=1|sz|(czszt+fxszx+fyszy+fzszz)(diasurfacetracerflux).
For strictly monotonic coordinate distributions sz(z) the effective vertical diasurface tracer flux in (30) only contains the contribution from the one vertical position of the isosurface, analogously to the relation for the diasurface velocity (25), such that
fdia,zs(s)=cs(s)udia,zs(s)+jdia,zs(s)forstrictlymonotonicsz(z)
holds, with the diffusive contribution jdia,zs defined equivalently to fdia,zs in (30),
jdia,zs(s)=Hηjdia,zz(z)|szz|δ{sz(z)s}dz
=Hη(KhczxszxKhczyszyKυczzszz)×δ{sz(z)s}dz.
Finally, the last term on the right hand side of (26) with fη,z considers the tracer flux across the ocean surface into water column volumes that contain the isosurface sη. In case of salinity, that is for c = S, there is no salt flux across the atmosphere–ocean interface, but potentially across ice–ocean interfaces.

e. Local squared tracer budget

Following Burchard and Rennau (2008), and according to (C3), an equation for the squared tracer can be derived from the tracer equation (12),
(cz)2t+fx(2)x+fy(2)y+fz(2)z=2czπχ,
with the fluxes of squared tracer defined in (C1) and the tracer variance dissipation
χ=2[Kh(czx)2+Kh(czy)2+Kυ(czz)2].
Flux boundary conditions at the ocean surface and the bottom are given by (C4d) and (C5d). Note, that the boundary fluxes of squared tracer are already determined by the volume and tracer fluxes according to (C4e) and (C5e).
Following the derivation of (C14) in appendix C, vertical integration of (35) over the parts of the water column where sz(z) ≤ s yields a budget equation for the cumulative squared tracer content in s space. Its sink term is given by the cumulative content of variance dissipation and quantifies the mixing inside D:
M(s)=Hηχ(z)H{sz(z)s}dz(cumulativetracermixing).
From (37) the mixing per s class can be defined as
m(s)=Ms=Hηχ(z)δ{sz(z)s}dz(mixingpersclass).

f. Mapping to tracer coordinates

In the previous sections we derived the governing equations in a general coordinate space by mapping to an arbitrary property s. In the following we will demonstrate, that for the special choice of tracer coordinates, i.e., s = c with c being the tracer in the tracer equation (12), useful relations between the diasurface (now diatracer) quantities can be derived. These relations in tracer space are the local equivalents of the integrated formulations presented in the introduction.

1) Relation between diatracer diffusive tracer fluxes and mixing

If the tracer is mapped to its own coordinate system, the effective vertical diatracer tracer flux (30) can be written as
fdia,zc(c)=cudia,zc(c)+jdia,zc(c),
which, in contrast to (33), does not require any monotonicity constraints on the tracer distribution. The effective vertical diatracer diffusive tracer flux jdia,zc obeys
jdia,zc(c)=(34b)Hη[Kh(czx)2+Kh(czy)2+Kυ(czz)2]=(36)χ/2δ{cz(z)c}dz=(37)12(Mc)=(38)m(c),
with the cumulative mixing M defined in (37) and m(c) being the mixing per tracer class. Relation (40) is the local equivalent of the integrated formulation (8) and the correct generalization of (9) in Li et al. (2022).

2) Ocean surface boundary fluxes

The volume and tracer fluxes across the ocean surface in (18) and (26) are formally considered in water column volumes that contain the isosurface sη. In tracer space the individual fluxes enter different tracer classes and thus are considered in different water column volumes:
uη,zH{sηs}uη,zH{cηc},
fη,zH{sηs}cηuη,zH{cηc}+jη,zH{cz(η)c}.
Each advective contribution enters the tracer class that corresponds to the tracer concentration cη transported by uη,z. The nonadvective tracer flux jη,z is assumed to enter the tracer class that is currently located at the ocean surface, given by cz(η). The consideration of advective contributions in tracer classes with concentrations different from the actual sea surface tracer concentration is mathematically valid, because tracer coordinates are not required to be monotonic and layers of cη can be constructed anywhere in the water column. For example, incoming freshwater volume fluxes always enter the zero salinity class, independent of the actual sea surface salinity.

3) Budget for internal tracer content

With (39), (41a), and (41b) a budget equation for the internal tracer content CI,≤ can be derived by combining (18) and (26),
t(CcD)=CI,(c)+x(Fx,cU)FIx,(c)+y(Fy,cV)FIy,(c)+jdia,zc(c)=Π(c)+(cηc)uη,zH{cηc}+jη,zH{cz(η)c},
with internal tracer fluxes FIx,≤ and FIy,≤. As motivated by Holmes et al. (2019) and Bladwell et al. (2021) for the integrated internal heat and salt content, respectively, CI,≤ quantifies the tracer content relative to the bounding tracer concentration c and the so-called external tracer content cD. As shown below, in this work the budget for the internal tracer content will be used to derive a relation between the diatracer velocity and the diatracer diffusive tracer fluxes. Alternative expressions for the internal tracer content and the associated internal tracer fluxes can be derived by means of (D1a) and partial integration:
CI,(c)=cD(c)dc,
FIx,(c)=cU(c)dc+Jx,(c),
FIy,(c)=cV(c)dc+Jy,(c).

4) Relation between diatracer velocity and diffusive tracer fluxes

With (43a)(43c) and (18), the budget for the internal tracer content (42) can be reformulated2 to
cudia,zc(c)dc+jdia,zc(c)+Jx,(c)x+Jy,(c)y=Π(c)+jη,zH{cz(η)c}.
Taking the derivative with respect to c yields the differential formulation
udia,zc(c)+jdia,zcc+x(Jx,c)+y(Jy,c)=Πc+jη,zδ{cz(η)c}.
Neglecting horizontal diffusive tracer fluxes and tracer production simplifies the general relation to
udia,zc(c)=jdia,zcc+jη,zδ{cz(η)c}.
Following (46), the local diatracer velocity (positive into the direction of increasing tracer concentration) across isosurfaces in the interior of the water column is always into the direction of larger jdia,zc, i.e., into the direction of stronger diasurface diffusive tracer fluxes. Of course, the net diffusive tracer transport is always into the direction of decreasing tracer concentration, but it can be either diverging (jdia,zc/c>0, i.e., dominated by the transport from the present tracer concentration toward lower one), or converging (jdia,zc/c<0, i.e., dominated by the transport from higher tracer concentration toward the present one). Divergence is associated with a decreasing volume of high tracer concentrations and an increasing volume of lower ones, i.e., a volume transport from high to low tracer concentrations (udia,zc<0). Vice versa, convergence is associated with water mass being transformed from low to high tracer concentrations (udia,zc>0). Relation (46) is the local equivalent3 of the integrated formulation (4) and the correct generalization of (18) in Li et al. (2022).

5) Relation between diatracer velocity and mixing

The combination of (40) and (46) directly formulates how the diatracer exchange flow is linked to local tracer variance dissipation:
udia,zc(c)=12mc+jη,zδ{cz(η)c}.
Following (47), the local diatracer velocity across isosurfaces in the interior of the water column is always into the direction of increased mixing. Relation (47) is the local equivalent of the integrated formulation (9), deduced in the introduction for salinity without salt fluxes across the ocean surface.

g. Integrated budgets and relations

Finally, we will show that integration of the local budgets and relations, derived in the previous sections, recovers the integral formulations presented in the introduction.

1) Integrated volume and tracer content budgets

Integration of the local volume and tracer content budgets (18) and (26) in tracer space over a horizontal area Az, and combination with (41a) and (41b) yields prognostic equations for volumes V(c) of all water with tracer concentrations ≤ c and their tracer content C(c):
tD(c)dAz=V(c)+[U(c)x+V(c)y]dAz=Q(c)Qr,(c)+udia,zc(c)dAz=Qdia(c)=uη,zH{cηc}dAz=Qη,(c),
tC(c)dAz=C(c)+[Fx,(c)x+Fy,(c)y]dAz=F(c)Fr,(c)+fdia,zc(c)dAz=Fdia(c)=Π(c)dAz=P(c)+[cηuη,z{cηc}+jη,z{cz(η)c}]dAz=Fη,(c).
The volumes are bounded by an open boundary and the isosurface with tracer concentration c. Figure 1 illustrates volumes for the special case of a stable salinity stratification. The equations derived here are also valid for nonmonotonic tracer distributions with inversions and noncontiguous volumes. The evolution of the volumes and their tracer content depends on the integrated fluxes entering across the open boundary (Q,F), integrated diatracer fluxes (Qdia,Fdia) and integrated fluxes entering across the ocean surface (Qη,Fη). Incoming integrated fluxes due to rivers (Qr,,Fr,) are treated separately from the integrated open boundary fluxes. The integrated tracer production is denoted by P. The integrated diatracer tracer flux is given by
Fdia(c)=(39)cQdia(c)+jdia,zc(c)dAz=Jdia(c),
and recovers (2).
The general formulation of integrated fluxes across the ocean surface and due to rivers in (48a) and (48b) supports loads with different tracer concentrations. For the special case of freshwater fluxes in salinity space the formulation simplifies, because the associated volume fluxes per salinity class are given by qη(S) = ∂Qη,≤/∂S = Qηδ{SηS} and qr(S) = ∂Qr,≤/∂S = Qrδ{SrS} with Sη = Sr = 0, respectively. Then the integrated fluxes become
Qη,(S)=Qη,Jη,=Fη,=0(integratedfreshwatersurfacefluxes),
Qr,(S)=Qr,Jr,=Fr,=0(integratedfreshwaterriverdischarge).
With these integrated freshwater fluxes the integral budgets (48a) and (48b) recover (5a) and (5b) used by, e.g., MacCready et al. (2002). Analogously, horizontal integration of (C15) and (C16) also recover Walin’s (1a) and (1b).

2) Relation between integrated diatracer diffusive tracer flux and mixing

The integrated diatracer diffusive tracer flux Jdia obeys
Jdia(c)=(40)12m(c)dAz=𝕞(c)=12cM(c)dAz=M(c),
which recovers (8). The integrated mixing per tracer class 𝕞 can be either obtained by horizontal integration of the local mixing per tracer class m, or following (6) as the derivative of the integrated tracer mixing M.

3) Integrated internal tracer content budget

By combination of the integral budgets (48a) and (48b), or directly by horizontal integration of the local formulation (42), the budget for the integrated internal tracer content can be obtained as
t(CcV)=CI,(c)(FcQ)=FI,(c)(Fr,cQr,)=FIr,(c)+Jdia(c)=P(c)+(Fη,cQη,)=FIη,(c),
which recovers the internal heat content budget (12) in Holmes et al. (2019) and the internal salt content budget (10) in Bladwell et al. (2021). A comparison to the notation of Bladwell et al. (2021) is given in Table 2.
With the budget for the internal salt content, the mixing per salinity class inside a classical estuarine volume can be diagnosed,
𝕞(S)=2{t[S(S)SV(S)][F(S)SQ(S)]+S(Qr+Qη)},
where (50a), (50b), and (51) has been used. This formulation generalizes (7) in Li et al. (2022) by including integrated fluxes across the ocean surface and the open boundary. For an estuarine volume completely bounded by the isohaline the integrated open boundary fluxes Q(S) and F(S) do not exist (MacCready and Geyer 2001; Burchard 2020). In this case, and after sufficiently long time averaging, (53) recovers the universal law of estuarine mixing (7).

4) Relation between integrated diatracer volume and diffusive tracer fluxes

The integrated formulation of relation (46) reads
Qdia(c)=Jdiac+cjη,zH{cz(η)c}dAz=Jη,,
and recovers (4). In contrast to the local relation (46), by means of the divergence theorem the integral relation (54) does not require to neglect horizontal diffusive tracer fluxes in general, but only across the open boundary. For a global volume such an open boundary does not even exist, and relation (54) holds for any horizontal diffusive tracer fluxes.

5) Relation between integrated diatracer volume flux and mixing

Combination of (51) and (54), or direct integration of the local relation (47), yields
Qdia=12𝕞c+Jη,c,
and generalizes (9) deduced in the introduction.

3. Demonstration

a. Idealized steady state and a word of caution

In a steady-state balance between vertical advection and vertical diffusion the tracer equation (12) reduces to wcz/∂z = ∂(Kυ ∂cz/∂z)/∂z. Assuming the resulting horizontally homogeneous tracer distribution is strictly monotonic with ∂cz/∂z > 0, the mixing inside an arbitrary dz interval is given by dM = χdz = 2Kυ(∂cz/∂z)2dz = 2Kυ ∂cz/∂zdcmdc. Combination of both equations yields w∂cz/∂z = ∂(m/2)/∂z = ∂(m/2)/∂c ∂cz/∂z and thus w = ∂(m/2)/∂c, which resembles the newly derived relation between the diatracer velocity and mixing (47). Indeed, for this special setting udia,zc=(25)udia,zz=(21),(22)w holds, and the simplified notation (omitting the superscripts for the actual coordinate space) is appropriate as long as the quantities are evaluated in the correct space. However, in more general settings we strongly encourage the use of the precise notation developed during the course of the manuscript, clearly indicating the coordinate space and preventing misunderstandings about derivatives. In this context, it is also important to remember that derivatives of the local quantities with respect to a tracer are evaluated at constant horizontal position (x, y), i.e., should not be misinterpreted as three-dimensional gradients in tracer direction.

b. Model results for an idealized tidal estuary

To demonstrate the derived relations in salinity space, the diagnostic calculation of diahaline mixing, volume and salt fluxes has been implemented into the coastal ocean model GETM (Burchard and Bolding 2002; Klingbeil et al. 2018). A simulation has been carried out for an idealized tidal estuary, which already represents realistic features of water mass transformation and overturning circulation, while still being a simple benchmark test case. Detailed results from realistic ocean modeling applications are presented in Henell et al. (2023) and Reese et al. (2023).

The estuary has a length of 100 km and a width of 500 m. The water depth decreases linearly from 15 m at the open ocean boundary to 5 m at the river end (see Fig. 4a). The estuary is forced by a harmonic semidiurnal tide with an elevation amplitude of 0.6 m and a constant salinity of 30 g kg−1 at the ocean boundary. At the river end a freshwater discharge of 50 m3 s−1 is prescribed. Horizontal diffusive fluxes are switched off. The estuarine circulation is discussed in detail in Burchard et al. (2019) and Klingbeil et al. (2019), such that we focus on the new diahaline analyses here. Figures 4b and 4c show the diagnosed diahaline diffusive fluxes and mixing between incoming high-saline ocean water and estuarine water of lower salinities. The resulting diahaline exchange flow in Figs. 4d–f indicates the transformation of the incoming high-saline water masses into lower-salinity water all the way up to the reversal point 60 km up-estuary, and the transformation of the outflowing brackish water masses into higher-salinity water again while flowing back toward the open ocean. This overturning circulation in salinity space complements the classical picture of estuarine circulation (e.g., Hansen and Rattray 1965).

Fig. 4.
Fig. 4.

Model results for an idealized tidal estuary in a periodic state. (a) Tidally averaged salinity distribution (blue contours) and salinity variance dissipation χ (colored) are shown. (b)–(d) The tidal mean diahaline diffusive salt flux jdia,zS, mixing per salinity class m, and diahaline velocity udia,zS have been diagnosed during model runtime, independently from (26), (38), and (18), respectively. In (b) and (d), the diahaline fluxes from higher to lower salinity are negative. (e),(f) The quantities are calculated from (b) and (c), respectively. The almost perfect agreement between the diahaline quantities in (b) and (c), as well as in (d)–(f) demonstrates the analytically derived relations (40), (46), and (47).

Citation: Journal of Physical Oceanography 53, 12; 10.1175/JPO-D-23-0130.1

4. Summary and conclusions

In this paper we analytically derived the local water mass transformation (WMT) framework for an individual water column. The derivation is based on the mapping of the governing equations from geopotential coordinates to an arbitrary tracer space. Exact local definitions for the effective vertical diasurface velocity, tracer flux and its diffusive contribution are given in (20), (30), and (34a), respectively. Mixing is defined in terms of the dissipation of tracer variance according to (37) and (38). In tracer space new relations between the local diatracer quantities and the mixing per tracer class have been derived as (46), (40), and (47). The key relation (47) between the effective vertical diatracer velocity and the mixing per tracer class directly formulates how the overturning circulation is linked to local tracer variance dissipation. It was demonstrated for different settings. The relations could also be used to study the diapycnal mixing and circulation in bottom boundary layers in the abyssal ocean described by Ferrari et al. (2016), and consistently link it to the large-scale overturning circulation. The local quantities and relations in tracer space represent an intermediate formulation between the direct ones in geopotential coordinates, for example presented by Garrett et al. (1995) and Ferrari et al. (2016), and the ones from the integral WMT framework presented in the introduction. The effective vertical diasurface quantities offer a convenient calculation of integrated diasurface fluxes by simple integration over the horizontal area instead of the complexly shaped area of the isosurface. By horizontally integrating the budget equations in tracer space and the relations between the diatracer quantities we finally recovered the well-known integral WMT formulations. Similar to the extension of the integral formulation by Hieronymus et al. (2014), an extension of our local formulation to multidimensional tracer space is possible. A summary of local and integral quantities and equations presented in this paper, as well as a comparison with published formulations, is given in Tables 1 and 2. The local diatracer fluxes and the local diatracer mixing provide new and valuable information on the spatial distribution of these quantities, and will enable new insights for coastal and large-scale oceanography.

1

In this paper all diasurface fluxes measured as “per unit horizontal area” are denoted by subscript “dia,z”.

2

Due to cuη,zH{cηc}dc=(cηc)uη,zH{cηc} the advective flux across the ocean surface in (42) cancels with the one from (18).

3

Here, jη,zδ{cz(η)c}=(jη,zH{cz(η)c})/c holds.

Acknowledgments.

This paper is a contribution to the project M5 of the Collaborative Research Centre TRR 181 “Energy Transfers in Atmosphere and Ocean” (project 274762653), funded by the German Research Foundation (DFG). We are grateful to Hans Burchard (IOW, Germany) for inspiring discussions and valuable comments. We also acknowledge the comments from two anonymous reviewers.

Data availability statement.

Model data shown in Fig. 4 have been obtained from a simulation with the coastal ocean model GETM (www.getm.eu). The model setup, simulation data, and scripts are archived at https://doi.org/10.5281/zenodo.8136528 (Klingbeil 2023a). The used model source code is archived at https://doi.org/10.5281/zenodo.8137091 (Klingbeil 2023b).

APPENDIX A

Unified Notation

For brevity we will use
D,>(s)=HηH,>{sz(z)s}dz
as a combined notation for (15a) and (15b) throughout the appendixes. In both cases the differential thickness is given by
±D,>(s)(A1)s=Hηδ{sz(z)s}dz.
In (A2) and in the following the upper sign (here “+”) always refers to the “≤” case and the lower sign (here “−”) to the “>” case. In analogy to D≤,> and ±∂D≤,>/∂s, for any quantity ψz the cumulative and differential formulations in s space can be obtained by
Ψ,>(s)=Hηψz(z)H,>{sz(z)s}dz,
±Ψ,>(s)s=Hηψz(z)δ{sz(z)s}dz.

APPENDIX B

Mapping as a Generalized Coordinate Transformation

In this paper the vertical coordinate z is replaced by an arbitrary coordinate s. For geometric coordinates, e.g., geopotential (s = z) or sigma coordinates with s = σ = (zη)/(η + H), the coordinate distribution sz(z) is strictly monotonic between the bottom at z = −H(x, y) and the ocean surface at z = η(x, y, t). The same holds for stably stratified temperature or salinity distributions (see Figs. 2a,c). Under this monotonicity constraint there exists a well-defined coordinate transformation to a unique inverse function zs(s) such that an arbitrary quantity ψ can be uniquely defined in s coordinates by
ψs(s)=ψz[zs(s)]forstrictlymonotonicsz(z),
with superscripts indicating the coordinate space the different functions are defined for. In contrast, for nonmonotonic distributions sz(z), e.g., for tracer coordinates where single tracer values may occur at multiple positions in the water column (e.g., Fig. 2e), there exists no function zs(s) and the usual transformation procedure (B1) cannot be applied. In this case a general mapping to s space can formally be defined by thickness-weighted quantities in infinitesimally small bins according to (B2a),
ψs(s)=Hηψz(z)δ{sz(z)s}dzHηδ{sz(z)s}dz,
=±Ψ,>(s)s±D,>(s)s.
The equivalent formulation (B2b) is of more practical use and is based on the cumulative thickness D≤,> and the cumulative quantity Ψ≤,> defined in (A1) and (A3). By construction (B2a) and (B2b) map data from z coordinates to s levels, and do not require any monotonicity constraints on the coordinate distribution sz(z). For strictly monotonic functions the integrals in (B2a) can be transformed by means of
Hηψz(z)δ{sz(z)s}dz=+ψz[zs(s)]|zss|δ{ss}ds=ψz[zs(s)]|zss|forstrictlymonotonicsz(z),
where in the last step the sampling property of the δ distribution was used. Analogously, the differential thickness ±∂D≤,>/∂s simplifies to the magnitude of the Jacobian of the coordinate transformation,
±D,>s=|zss|=|szz|1forstrictlymonotonicsz(z),
and (B2a) recovers (B1). The formulations (B2a) and (B2b) are motivated by (13) and (14) in Klingbeil et al. (2019) from the context of generalized thickness-weighted averaging.

APPENDIX C

Derivation of the p-Generalized Tracer Budget in s Space

Because the integrations for volume, tracer and tracer-square content are similar, the following derivations consider the tracer concentration c to the power of an arbitrary integer p ≥ 0. The flux of cp in the x direction (the fluxes into the other directions are defined accordingly) is given by
fx(p)=u(cz)pKh(cz)px=jx(p),=p(cz)p1fx(p1)(cz)pu,
with fx(0)=u and fx(1)=fx being the velocity and tracer flux, respectively. According to (C1), for p ≥ 2, fx(p) is determined by u and fx. With the incompressibility constraint (10) it can be shown that
fx(p)x+fy(p)y+fz(p)z=p(cz)p1(fxx+fyy+fzz)p(p1)2(cz)p2×2[Kh(czx)2+Kh(czy)2+Kυ(czz)2]=χ.
Multiplication of the tracer equation (12) by p(cz)p1 and reformulating the fluxes with (C2) yields the prognostic equation for (cz)p:
(cz)pt+fx(p)x+fy(p)y+fz(p)z=p(cz)p1πp(p1)2(cz)p2χ.
Boundary fluxes across the ocean surface and the bottom per unit horizontal area (positive into the water column) are defined by
fη,z(p)=fη(p)|nηz|1
=fη(p)|(ηz)|
={[cz(η)]p[u(η)uη]+j(p)(η)}(ηz)
=[cz(η)]pηt+fx(p)(η)ηx+fy(p)(η)ηyfz(p)(η)
=(C1)p[cz(η)]p1fη,z(p1)[cz(η)]puη,z, and
fH,z(p)=fH(p)|nHz|1
=fH(p)|(z+H)|
={[cz(H)]pu(H)+j(p)(H)}(z+H)
=fx(p)(H)Hx+fy(p)(H)Hy+fz(p)(H)
=(C1)p[cz(H)]p1fH,z(p1)[cz(H)]puH,z,
with the unit normal vectors nη = ∇(ηz)/|∇(ηz)| and nH = ∇(z + H)/|∇(z + H)|, and the vertical unit vector z. For p = 0 the relations (C4d) and (C5d) recover the well-known kinematic boundary conditions (11a) and (11b), and for p = 1 the tracer flux boundary conditions (13a) and (13b). According to (C4e) and (C5e), for p ≥ 2 the boundary fluxes are determined by the ones for volume and tracer.
Following (A3) the p-generalized cumulative tracer content is defined as
C,>(p)(s)=Hη[cz(z)]pH,>{sz(z)s}dz,
with C,>(0)=D,> and C,>(1)=C,>. In analogy the p-generalized cumulative tracer fluxes are defined as
Fx,,>(p)(s)=Hηfx(p)(z)H,>{sz(z)s}dz,
Fy,,>(p)(s)=Hηfy(p)(z)H,>{sz(z)s}dz,
with Fx,,>(0)=U,>,Fy,,>(0)=V,>, Fx,,>(1)=Fx,,> and Fy,,>(1)=Fy,,>. For (C6)(C7b) Leibniz rule (E1) yields
C,>(p)(s)t=Hη(cz)psztδ{sz(z)s}dz+Hη(cz)ptH,>{sz(z)s}dz+[cz(η)]pηtH,>{sηs},
Fx,,>(p)(s)x=Hηfx(p)szxδ{sz(z)s}dz+Hηfx(p)xH,>{sz(z)s}dz+fx(p)(η)ηxH,>{sηs}+fx(p)(H)HxH,>{sHs},
Fy,,>(p)(s)y=Hηfy(p)szyδ{sz(z)s}dz+Hηfy(p)yH,>{sz(z)s}dz+fy(p)(η)ηyH,>{sηs}+fy(p)(H)HyH,>{sHs}
Application of (E2) for the vertical flux gives
Hηz[H,>{sz(z)s}fz(p)]dz=Hηδ{sz(z)s}fz(p)szzdz+HηH,>{sz(z)s}fz(p)zdzfz(p)(η)H,>{sηs}fz(p)(H)H,>{sHs}
With (C4d), (C5d), and (C9), summation of (C8a)(C8c) yields
C,>(p)(s)t+Fx,,>(p)(s)x+Fy,,>(p)(s)y=Hη[(cz)pszt+fx(p)szx+fy(p)szy+fz(p)szz]δ{sz(z)s}dz+Hη[(cz)pt+fx(p)x+fy(p)y+fz(p)z]H,>{sz(z)s}dz+fη,z(p)H,>{sηs}+fH,z(p)H,>{sHs}
The integrand of the first integral on the right hand side of (C10) can be further reformulated in terms of the local diasurface flux fdia(p),z, defined as the component of the relative flux normal to the isosurface that is currently located at the specific position z and moving with velocity us:
fdia(p),z=(cz)p(uus)n+j(p)n=1|sz|[(cz)pDszDt+jx(p)szx+jy(p)szy+jz(p)szz]
=1|sz|[(cz)pszt+fx(p)szx+fy(p)szy+fz(p)szz]
=(C1)p(cz)p1fdiaz(p1)(cz)pudiaz.
The local diasurface velocity and the local diasurface tracer flux are given by fdia(0),z=udiaz in (22) and fdia(1),z=fdiaz in (32), respectively. With the unit normal vector n = ∇sz/|∇sz|, the diasurface fluxes are positive into the direction of increasing s. In analogy to the formulations of the fluxes across the ocean surface and bottom (C4a) and (C5a), the associated local diasurface tracer flux per unit horizontal area formally can be defined as
fdia,z(p),z=fdia(p),z|nz|1=fdia(p),z|szz|1|sz|.
Using (C11b) and (C12), the first integral on the right hand side of (C10) can be defined as the effective vertical diasurface flux
fdia,z(p),s(s)=Hηfdia,z(p),z(z)|szz|δ{sz(z)s}dz,
The integrand of the second integral on the right hand side of (C10) can be replaced by inserting the prognostic equation for (cz)p, (C3), such that the following budget equation for the p-generalized cumulative tracer content can be derived:
C,>(p)(s)t+Fx,,>(p)(s)x+Fy,,>(p)(s)y±fdia,(z)(p),s(s)=Hη[p(cz)p1πp(p1)2(cz)p2χ]H,>{sz(z)s}dz+fη,z(p)H,>{sηs}+fH,z(p)H,>{sHs}
Neglecting bottom fluxes, for p = 0 the volume budget
D,>(s)t+U,>(s)x+V,>(s)y±udia,zs(s)=uη,zH,>{sηs},
and for p = 1 the tracer budget
C,>(s)t+Fx,,>(s)x+Fy,,>(s)y±fdia,zs(s)=Π,>(s)+fη,zH,>{sηs},
with the cumulative tracer production
Π,>(s)=Hηπ(z)H,>{sz(z)s}dz
are obtained, respectively. Taking the derivative of (C14) with respect to s and using (B2b) yields the differential formulation
t(±D,>sc(p),s)+x(±D,>sfx(p),s)+y(±D,>sfy(p),s)+fdia,z(p),ss=Hη[p(cz)p1πp(p1)2(cz)p2χ]δ{sz(z)s}dz+fη,z(p)δ{sηs}+fH,z(p)δ{sHs},
with the thickness-weighted quantities in s space, c(p),s=(C,>(p)/s)/(D,>/s), fx(p),s=(Fx,,>(p)/s)/(D,>/s) and fy(p),s accordingly, from (B2b). For strictly monotonic coordinates, i.e., with (B4), the budget (C18) with p = 0 and p = 1 recovers the budget equations in generalized vertical coordinates which are the basis for the layer-integrated equations in ocean models, see, e.g., Griffies (2004) and Klingbeil et al. (2018).

APPENDIX D

Further Relations in Tracer Space

For tracer coordinates, i.e., for s = c, and by means of the Newton–Leibniz theorem the following useful transformation relation holds for an arbitrary quantity ψ:
Hη(cz)pψzH,>{ cz(z)c }dz=Hη(cz)pψzH,>{ cz(z)() }dz=0+c[ cHη(cz)pψzH,>{ cz(z)c }dz ]dc=±c[ Hη(cz)pψzδ{ cz(z)c }dz ]dc=±c(c)p[ Hηψzδ{ cz(z)c }dz ]dc=±c(c)p[ ±cHηψzH,>{ cz(z)c }dz ]dc=(A3)±c(c)p(±Ψ,>c)dc=+(c)p(±Ψ,>c)H,>{ cz(z)c }dc=(B2b)+(c)p[ ±D,>cψc(c) ]H,>{ cz(z)c }dc.
By means of (D1a), the integration on the right hand side of (C14) can be carried out in tracer space,
C,>(p)(c)t+Fx,,>(p)(c)x+Fy,,>(p)(c)y±fdia,z(p),c(c)=±c[p(c)p1(±Π,>c)p(p1)2(c)p2(±M,>c)]dc+fη,z(p)H,>{cηc}+fH,z(p)H,>{cHc},
with the cumulative tracer production Π≤,> defined in (C17), and the cumulative tracer mixing M,>(s)=Hηχ(z)H,>{sz(z)s}dz. Subsequent derivation with respect to c yields
t(±D,>ccp)+x[±D,>cfx(p),c]+y[±D,>cfy(p),c]+fdia,z(p),cc=pcp1(±Π,>c)p(p1)2cp2(±M,>c)+fη,z(p)δ{cηc}+fH,z(p)δ{cHc},
which simplifies the first term and the right hand side from the general formulation (C18).

APPENDIX E

Leibniz Rule

For ξ ∈{x, y, t} the Leibniz rule reads
ξH(ξ)η(ξ)ψz(ξ,z)H,>{sz(ξ,z)s}dz=Hηψzszξδ{sz(ξ,z)s}dz+HηψzξH,>{sz(z)s}dz+ψz(η)ηξH,>{sz(η)s}+ψz(H)HξH,>{sz(H)s}.
But note that the bottom topography H will be considered as constant in time t in this manuscript.
For the vertical direction it can be shown that
Hηz[ψzH,>{sz(z)s}]dz=Hηψzszzδ{sz(z)s}dz+HηψzzH,>{sz(z)s}dz,ψz(η)H,>{sz(η)s}ψz(H)H,>{sz(H)s},
where the third line follows directly from the first one, using the Newton–Leibniz theorem.

In the rest of the manuscript we replace sz(η) and sz(−H) in the arguments of the Heaviside functions in (E1) and (E2) by sη and s−H.

REFERENCES

  • Bladwell, C., R. M. Holmes, and J. D. Zika, 2021: Internal salt content: A useful framework for understanding the oceanic branch of the water cycle. J. Phys. Oceanogr., 51, 21672179, https://doi.org/10.1175/JPO-D-20-0212.1.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., 2020: A universal law of estuarine mixing. J. Phys. Oceanogr., 50, 8193, https://doi.org/10.1175/JPO-D-19-0014.1.

  • Burchard, H., and K. Bolding, 2002: GETM, a general estuarine transport model: Scientific documentation. Tech. Rep. EUR 20253 EN, 159 pp., https://publications.jrc.ec.europa.eu/repository/handle/JRC23237.

  • Burchard, H., and H. Rennau, 2008: Comparative quantification of physically and numerically induced mixing in ocean models. Ocean Modell., 20, 293311, https://doi.org/10.1016/j.ocemod.2007.10.003.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., and Coauthors, 2018: The Knudsen theorem and the total exchange flow analysis framework applied to the Baltic Sea. Prog. Oceanogr., 165, 268286, https://doi.org/10.1016/j.pocean.2018.04.004.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., X. Lange, K. Klingbeil, and P. MacCready, 2019: Mixing estimates for estuaries. J. Phys. Oceanogr., 49, 631648, https://doi.org/10.1175/JPO-D-18-0147.1.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., U. Gräwe, K. Klingbeil, N. Koganti, X. Lange, and M. Lorenz, 2021: Effective diahaline diffusivities in estuaries. J. Adv. Model. Earth Syst., 13, e2020MS002307, https://doi.org/10.1029/2020MS002307.

    • Search Google Scholar
    • Export Citation
  • Döös, K., and D. J. Webb, 1994: The deacon cell and the other meridional cells of the southern ocean. J. Phys. Oceanogr., 24, 429442, https://doi.org/10.1175/1520-0485(1994)024<0429:TDCATO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ferrari, R., A. Mashayek, T. J. McDougall, M. Nikurashin, and J.-M. Campin, 2016: Turning ocean mixing upside down. J. Phys. Oceanogr., 46, 22392261, https://doi.org/10.1175/JPO-D-15-0244.1.

    • Search Google Scholar
    • Export Citation
  • Garrett, C., K. Speer, and E. Tragou, 1995: The relationship between water mass formation and the surface buoyancy flux, with application to Phillips’ Red Sea model. J. Phys. Oceanogr., 25, 16961705, https://doi.org/10.1175/1520-0485(1995)025<1696:TRBWMF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Griffies, S. M., 2004: Fundamentals of Ocean Climate Models. Princeton University Press, 518 pp., https://doi.org/10.2307/j.ctv301gzg.

  • Groeskamp, S., S. M. Griffies, D. Iudicone, R. Marsh, A. J. G. Nurser, and J. D. Zika, 2019: The water mass transformation framework for ocean physics and biogeochemistry. Annu. Rev. Mar. Sci., 11, 271305, https://doi.org/10.1146/annurev-marine-010318-095421.

    • Search Google Scholar
    • Export Citation
  • Hansen, D. V., and M. Rattray Jr., 1965: Gravitational circulation in straits and estuaries. J. Mar. Res., 23, 104122.

  • Henell, E., H. Burchard, U. Gräwe, and K. Klingbeil, 2023: Spatial composition of the diahaline overturning circulation in a fjord-type, non-tidal estuarine system. J. Geophys. Res. Oceans, https://doi.org/10.1029/2023JC019862, in press.

    • Search Google Scholar
    • Export Citation
  • Hetland, R. D., 2005: Relating river plume structure to vertical mixing. J. Phys. Oceanogr., 35, 16671688, https://doi.org/10.1175/JPO2774.1.

    • Search Google Scholar
    • Export Citation
  • Hieronymus, M., J. Nilsson, and J. Nycander, 2014: Water mass transformation in salinity–temperature space. J. Phys. Oceanogr., 44, 25472568, https://doi.org/10.1175/JPO-D-13-0257.1.

    • Search Google Scholar
    • Export Citation
  • Holmes, R. M., J. D. Zika, and M. H. England, 2019: Diathermal heat transport in a global ocean model. J. Phys. Oceanogr., 49, 141161, https://doi.org/10.1175/JPO-D-18-0098.1.

    • Search Google Scholar
    • Export Citation
  • Holmes, R. M., J. D. Zika, S. M. Griffies, A. M. Hogg, A. E. Kiss, and M. H. England, 2021: The geography of numerical mixing in a suite of global ocean models. J. Adv. Model. Earth Syst., 13, e2020MS002333, https://doi.org/10.1029/2020MS002333.

    • Search Google Scholar
    • Export Citation
  • Klingbeil, K., 2023a: Diahaline warner estuary (GETM model results). Zenodo, https://doi.org/10.5281/zenodo.8136528.

  • Klingbeil, K., 2023b: Source code for the coastal ocean model GETM (tef branch). Zenodo, https://doi.org/10.5281/zenodo.8137091.

  • Klingbeil, K., F. Lemarié, L. Debreu, and H. Burchard, 2018: The numerics of hydrostatic structured-grid coastal ocean models: State of the art and future perspectives. Ocean Modell., 125, 80105, https://doi.org/10.1016/j.ocemod.2018.01.007.

    • Search Google Scholar
    • Export Citation
  • Klingbeil, K., J. Becherer, E. Schulz, H. E. de Swart, H. M. Schuttelaars, A. Valle-Levinson, and H. Burchard, 2019: Thickness-weighted averaging in tidal estuaries and the vertical distribution of the Eulerian residual transport. J. Phys. Oceanogr., 49, 18091826, https://doi.org/10.1175/JPO-D-18-0083.1.

    • Search Google Scholar
    • Export Citation
  • Knudsen, M., 1900: Ein hydrographischer Lehrsatz. Ann. Hydrogr. Maritimen Meteor., 28, 316320.

  • Li, X., M. Lorenz, K. Klingbeil, E. Chrysagi, U. Gräwe, J. Wu, and H. Burchard, 2022: Salinity mixing and diahaline exchange flow in a large multi-outlet estuary with islands. J. Phys. Oceanogr., 52, 21112127, https://doi.org/10.1175/JPO-D-21-0292.1.

    • Search Google Scholar
    • Export Citation
  • MacCready, P., 2011: Calculating estuarine exchange flow using isohaline coordinates. J. Phys. Oceanogr., 41, 11161124, https://doi.org/10.1175/2011JPO4517.1.

    • Search Google Scholar
    • Export Citation
  • MacCready, P., and W. R. Geyer, 2001: Estuarine salt flux through an isohaline surface. J. Geophys. Res., 106, 11 62911 637, https://doi.org/10.1029/2001JC900006.

    • Search Google Scholar
    • Export Citation
  • MacCready, P., R. D. Hetland, and W. R. Geyer, 2002: Long-term isohaline salt balance in an estuary. Cont. Shelf Res., 22, 15911601, https://doi.org/10.1016/S0278-4343(02)00023-7.

    • Search Google Scholar
    • Export Citation
  • Nurser, A. J. G., and M.-M. Lee, 2004: Isopycnal averaging at constant height. Part I: The formulation and a case study. J. Phys. Oceanogr., 34, 27212739, https://doi.org/10.1175/JPO2649.1.

    • Search Google Scholar
    • Export Citation
  • Reese, L., U. Gräwe, K. Klingbeil, X. Li, M. Lorenz, and H. Burchard, 2023: Local mixing determines spatial structure of diahaline exchange flow in a mesotidal estuary: A study of extreme runoff conditions. J. Phys. Oceanogr., https://doi.org/10.1175/JPO-D-23-0052.1, in press.

    • Search Google Scholar
    • Export Citation
  • Walin, G., 1977: A theoretical framework for the description of estuaries. Tellus, 29A, 128136, https://doi.org/10.3402/tellusa.v29i2.11337.

    • Search Google Scholar
    • Export Citation
  • Walin, G., 1982: On the relation between sea-surface heat flow and thermal circulation in the ocean. Tellus, 34A, 187195, https://doi.org/10.3402/tellusa.v34i2.10801.

    • Search Google Scholar
    • Export Citation
  • Wang, T., W. R. Geyer, and P. MacCready, 2017: Total exchange flow, entrainment, and diffusive salt flux in estuaries. J. Phys. Oceanogr., 47, 12051220, https://doi.org/10.1175/JPO-D-16-0258.1.

    • Search Google Scholar
    • Export Citation
  • Wolfe, C. L., 2014: Approximations to the ocean’s residual circulation in arbitrary tracer coordinates. Ocean Modell., 75, 2035, https://doi.org/10.1016/j.ocemod.2013.12.004.

    • Search Google Scholar
    • Export Citation
Save
  • Bladwell, C., R. M. Holmes, and J. D. Zika, 2021: Internal salt content: A useful framework for understanding the oceanic branch of the water cycle. J. Phys. Oceanogr., 51, 21672179, https://doi.org/10.1175/JPO-D-20-0212.1.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., 2020: A universal law of estuarine mixing. J. Phys. Oceanogr., 50, 8193, https://doi.org/10.1175/JPO-D-19-0014.1.

  • Burchard, H., and K. Bolding, 2002: GETM, a general estuarine transport model: Scientific documentation. Tech. Rep. EUR 20253 EN, 159 pp., https://publications.jrc.ec.europa.eu/repository/handle/JRC23237.

  • Burchard, H., and H. Rennau, 2008: Comparative quantification of physically and numerically induced mixing in ocean models. Ocean Modell., 20, 293311, https://doi.org/10.1016/j.ocemod.2007.10.003.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., and Coauthors, 2018: The Knudsen theorem and the total exchange flow analysis framework applied to the Baltic Sea. Prog. Oceanogr., 165, 268286, https://doi.org/10.1016/j.pocean.2018.04.004.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., X. Lange, K. Klingbeil, and P. MacCready, 2019: Mixing estimates for estuaries. J. Phys. Oceanogr., 49, 631648, https://doi.org/10.1175/JPO-D-18-0147.1.

    • Search Google Scholar
    • Export Citation
  • Burchard, H., U. Gräwe, K. Klingbeil, N. Koganti, X. Lange, and M. Lorenz, 2021: Effective diahaline diffusivities in estuaries. J. Adv. Model. Earth Syst., 13, e2020MS002307, https://doi.org/10.1029/2020MS002307.

    • Search Google Scholar
    • Export Citation
  • Döös, K., and D. J. Webb, 1994: The deacon cell and the other meridional cells of the southern ocean. J. Phys. Oceanogr., 24, 429442, https://doi.org/10.1175/1520-0485(1994)024<0429:TDCATO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ferrari, R., A. Mashayek, T. J. McDougall, M. Nikurashin, and J.-M. Campin, 2016: Turning ocean mixing upside down. J. Phys. Oceanogr., 46, 22392261, https://doi.org/10.1175/JPO-D-15-0244.1.

    • Search Google Scholar
    • Export Citation
  • Garrett, C., K. Speer, and E. Tragou, 1995: The relationship between water mass formation and the surface buoyancy flux, with application to Phillips’ Red Sea model. J. Phys. Oceanogr., 25, 16961705, https://doi.org/10.1175/1520-0485(1995)025<1696:TRBWMF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Griffies, S. M., 2004: Fundamentals of Ocean Climate Models. Princeton University Press, 518 pp., https://doi.org/10.2307/j.ctv301gzg.

  • Groeskamp, S., S. M. Griffies, D. Iudicone, R. Marsh, A. J. G. Nurser, and J. D. Zika, 2019: The water mass transformation framework for ocean physics and biogeochemistry. Annu. Rev. Mar. Sci., 11, 271305, https://doi.org/10.1146/annurev-marine-010318-095421.

    • Search Google Scholar
    • Export Citation
  • Hansen, D. V., and M. Rattray Jr., 1965: Gravitational circulation in straits and estuaries. J. Mar. Res., 23, 104122.

  • Henell, E., H. Burchard, U. Gräwe, and K. Klingbeil, 2023: Spatial composition of the diahaline overturning circulation in a fjord-type, non-tidal estuarine system. J. Geophys. Res. Oceans, https://doi.org/10.1029/2023JC019862, in press.

    • Search Google Scholar
    • Export Citation
  • Hetland, R. D., 2005: Relating river plume structure to vertical mixing. J. Phys. Oceanogr., 35, 16671688, https://doi.org/10.1175/JPO2774.1.

    • Search Google Scholar
    • Export Citation
  • Hieronymus, M., J. Nilsson, and J. Nycander, 2014: Water mass transformation in salinity–temperature space. J. Phys. Oceanogr., 44, 25472568, https://doi.org/10.1175/JPO-D-13-0257.1.

    • Search Google Scholar
    • Export Citation
  • Holmes, R. M., J. D. Zika, and M. H. England, 2019: Diathermal heat transport in a global ocean model. J. Phys. Oceanogr., 49, 141161, https://doi.org/10.1175/JPO-D-18-0098.1.

    • Search Google Scholar
    • Export Citation
  • Holmes, R. M., J. D. Zika, S. M. Griffies, A. M. Hogg, A. E. Kiss, and M. H. England, 2021: The geography of numerical mixing in a suite of global ocean models. J. Adv. Model. Earth Syst., 13, e2020MS002333, https://doi.org/10.1029/2020MS002333.

    • Search Google Scholar
    • Export Citation
  • Klingbeil, K., 2023a: Diahaline warner estuary (GETM model results). Zenodo, https://doi.org/10.5281/zenodo.8136528.

  • Klingbeil, K., 2023b: Source code for the coastal ocean model GETM (tef branch). Zenodo, https://doi.org/10.5281/zenodo.8137091.

  • Klingbeil, K., F. Lemarié, L. Debreu, and H. Burchard, 2018: The numerics of hydrostatic structured-grid coastal ocean models: State of the art and future perspectives. Ocean Modell., 125, 80105, https://doi.org/10.1016/j.ocemod.2018.01.007.

    • Search Google Scholar
    • Export Citation
  • Klingbeil, K., J. Becherer, E. Schulz, H. E. de Swart, H. M. Schuttelaars, A. Valle-Levinson, and H. Burchard, 2019: Thickness-weighted averaging in tidal estuaries and the vertical distribution of the Eulerian residual transport. J. Phys. Oceanogr., 49, 18091826, https://doi.org/10.1175/JPO-D-18-0083.1.

    • Search Google Scholar
    • Export Citation
  • Knudsen, M., 1900: Ein hydrographischer Lehrsatz. Ann. Hydrogr. Maritimen Meteor., 28, 316320.

  • Li, X., M. Lorenz, K. Klingbeil, E. Chrysagi, U. Gräwe, J. Wu, and H. Burchard, 2022: Salinity mixing and diahaline exchange flow in a large multi-outlet estuary with islands. J. Phys. Oceanogr., 52, 21112127, https://doi.org/10.1175/JPO-D-21-0292.1.

    • Search Google Scholar
    • Export Citation
  • MacCready, P., 2011: Calculating estuarine exchange flow using isohaline coordinates. J. Phys. Oceanogr., 41, 11161124, https://doi.org/10.1175/2011JPO4517.1.

    • Search Google Scholar
    • Export Citation
  • MacCready, P., and W. R. Geyer, 2001: Estuarine salt flux through an isohaline surface. J. Geophys. Res., 106, 11 62911 637, https://doi.org/10.1029/2001JC900006.

    • Search Google Scholar
    • Export Citation
  • MacCready, P., R. D. Hetland, and W. R. Geyer, 2002: Long-term isohaline salt balance in an estuary. Cont. Shelf Res., 22, 15911601, https://doi.org/10.1016/S0278-4343(02)00023-7.

    • Search Google Scholar
    • Export Citation
  • Nurser, A. J. G., and M.-M. Lee, 2004: Isopycnal averaging at constant height. Part I: The formulation and a case study. J. Phys. Oceanogr., 34, 27212739, https://doi.org/10.1175/JPO2649.1.

    • Search Google Scholar
    • Export Citation
  • Reese, L., U. Gräwe, K. Klingbeil, X. Li, M. Lorenz, and H. Burchard, 2023: Local mixing determines spatial structure of diahaline exchange flow in a mesotidal estuary: A study of extreme runoff conditions. J. Phys. Oceanogr., https://doi.org/10.1175/JPO-D-23-0052.1, in press.

    • Search Google Scholar
    • Export Citation
  • Walin, G., 1977: A theoretical framework for the description of estuaries. Tellus, 29A, 128136, https://doi.org/10.3402/tellusa.v29i2.11337.

    • Search Google Scholar
    • Export Citation
  • Walin, G., 1982: On the relation between sea-surface heat flow and thermal circulation in the ocean. Tellus, 34A, 187195, https://doi.org/10.3402/tellusa.v34i2.10801.

    • Search Google Scholar
    • Export Citation
  • Wang, T., W. R. Geyer, and P. MacCready, 2017: Total exchange flow, entrainment, and diffusive salt flux in estuaries. J. Phys. Oceanogr., 47, 12051220, https://doi.org/10.1175/JPO-D-16-0258.1.

    • Search Google Scholar
    • Export Citation
  • Wolfe, C. L., 2014: Approximations to the ocean’s residual circulation in arbitrary tracer coordinates. Ocean Modell., 75, 2035, https://doi.org/10.1016/j.ocemod.2013.12.004.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    Idealized illustration of control volumes bounded by an isohaline with salinity S and a boundary transect toward the open ocean. For simplicity a stable salinity stratification is assumed, such that the volume V>(S) containing water with salinities larger than S is located below the isohaline, and the remaining volume V(S) containing water with salinities less or equal to S above the isohaline. The integrated diahaline volume flux Qdia is a measure for the water mass transformation. Integrated surface boundary and open boundary volume fluxes are depicted as Qη,≤, Q, and Q>, defined based on how the salinity of the transported water relates to S. Integrated salt fluxes are denoted by F and alter the salt contents S(S) and S>(S) inside the corresponding volumes. Surface boundary fluxes Qη,> and Fη,>, transporting or exchanging water with salinities larger than S, are not sketched for simplicity. Additional freshwater river discharge is depicted as Qr. The arrows indicate the direction of positive fluxes. Walin (1977) presented the budget equations for V>(S) and S>(S); see (1a) and (1b). In many studies on estuaries the budgets for V(S) and S(S), (5a) and (5b) are used.

  • Fig. 2.

    Control volumes within a water column separated by an isosurface of property s. The vertical distribution of the property is given by the function sz(z). (a),(c),(e) Examples for strictly monotonic and nonmonotonic distributions are illustrated. For a prescribed value s the filter functions H and H>, defined in (14a) and (14b), identify the parts of the water column with sz(z) ≤ s or sz(z) > s, respectively. (b),(d),(f) Examples for the location of the isosurface in xz space in the vicinity of the water column are sketched. The local volume budget (18) describes the relation between the cumulative thickness D and horizontal transport U, as well as the diasurface velocities udiaz for the two control volumes.

  • Fig. 3.

    Sketch for relation between area elements and diasurface velocities (udiaz,udia,zz). The infinitesimal area element dA is located on the isosurface with the local unit normal vector n. Its horizontal projection is given by dAz, with z being the vertical unit vector. By construction udiazdA=udia,zzdAz holds; see (23). The same relations are valid for the total and diffusive tracer fluxes (fdiaz,fdia,zz) and (jdiaz,jdia,zz), respectively.

  • Fig. 4.

    Model results for an idealized tidal estuary in a periodic state. (a) Tidally averaged salinity distribution (blue contours) and salinity variance dissipation χ (colored) are shown. (b)–(d) The tidal mean diahaline diffusive salt flux jdia,zS, mixing per salinity class m, and diahaline velocity udia,zS have been diagnosed during model runtime, independently from (26), (38), and (18), respectively. In (b) and (d), the diahaline fluxes from higher to lower salinity are negative. (e),(f) The quantities are calculated from (b) and (c), respectively. The almost perfect agreement between the diahaline quantities in (b) and (c), as well as in (d)–(f) demonstrates the analytically derived relations (40), (46), and (47).

All Time Past Year Past 30 Days
Abstract Views 1402 837 0
Full Text Views 1794 1521 119
PDF Downloads 539 249 29