1 Background theory
1.1 The PV flux
The flux form of the PV equation has been discussed at length in Haynes and
McIntyre (1987), Marshall and Nurser (1992), Rhines (1993), Schar (1993) and
Bretherton and Schar (1993). Here we very briefly outline the approach,
drawing out issues that are particularly pertinent to the ocean.
First let's clarify some issues related to the flux-form of the PV equation.
Following Bretherton and Schar (1993) we can write, making use of the
definition of PV and that Ñ·w = 0:
|
¶
¶t
|
( rQ ) = |
¶
¶t
|
(-w ·Ñs) = - |
¶
¶t
|
Ñ·(ws) = -Ñ·J |
| (1) |
with J = [(¶)/( ¶t)](ws).
This shows that one can always write a conservation law in flux form
from the definition of PV. This is true for any scalar field s and
non-divergent vector w. If one adds any non-divergent vector
to J, the flux-form equation will still be satisfied. The problem
is then to set the gauge and this choice has to be made on physical grounds.
Expanding the partial derivative in J, one can write:
Now, using the definition w = 2W+Ñ×u and some vector identity, one can write:
J = w |
¶s
¶t
|
+ |
¶u
¶t
|
×Ñs+ Ñ× |
æ ç
è
|
s |
¶u
¶t
|
|
ö ÷
ø
|
|
| (3) |
Note that the third term above is non-divergent and hence does not contribute to the PV flux-form equation. The problem is then to determine the 'gauge' X such that:
J = w |
¶s
¶t
|
+ |
¶u
¶t
|
×Ñs+X |
| (4) |
has physical significance.
X = 0 is not a satisfying choice because that would imply, as noted by BC, that J = 0 in the steady state. Instead, we require that J reduces to the advective flux rQu in the absence of diabatic and frictional forcings. This choice will uniquely determine X and hence J. Using the previous expression for J, X can be written:
X = rQu - w |
¶s
¶t
|
- |
¶u
¶t
|
×Ñs |
| (5) |
Now we use the definitions of PV and the thermodynamic and momentum equations:
|
¶u
¶t
|
= -w×u-Ñp+ |
F
ro
|
Ñr¢ |
| (8) |
In (8), the density r of the Boussinesq fluid is given by r = ro+r¢ where ro is a constant reference density. F is the geopotential and p is the Bernoulli function written in the Boussinesq approximation:
p = |
| u| 2
2
|
+ |
p
ro
|
+ |
r¢
ro
|
F |
| (9) |
In (9), p is the deviation of the pressure from that of a resting, hydrostatically balanced ocean. Replacing in (5), we have:
X = -( w·Ñs)u+w( u·Ñs)+( w×u) ×Ñs+Ñp×Ñs |
| (10) |
(omitting the term [(F)/( ro)]Ñr¢×Ñs). The first three terms on the right hand side cancel leaving:
X sets the gauge for J. Replacing in (4), we find:
J = w |
¶s
¶t
|
+ |
æ ç
è
|
|
¶u
¶t
|
+Ñp |
ö ÷
ø
|
×Ñs |
| (12) |
expression written by Schar (1993) in the atmospheric context.
Note that Eq.(12):
- as noted by Haynes and McIntyre (1987), reveals the `impermeability theorem' in a transparent way - the first
term on the lhs, when projected in the direction normal to the s
surface, is equal to vsrQ where vs = -| Ñs| -1[(¶s)/( ¶t)] is the velocity of
the s surface normal to itself. The remaining terms represent a flux
that is always parallel to the s surface.
- shows that J can be evaluated without explicit reference
to frictional (F) and buoyancy ([(Ds)/( Dt)]) sources, if the rates of change of s and u are known. Of course, as seen above, part of this result originates from the definition of PV. Note however the crucial role played by the 'gauge' term Ñp×Ñs.
- shows that in the steady state p is the streamfunction for the J vector on s surfaces, a very general result first noted
by Schar (1993) and Bretherton and Schar (1993).
To obtain an expression for J that includes explicitely the frictional and buoyancy sources, we follow the same procedure as above, replacing in (12) the rates of change of s and u from the thermodynamic and momentum equations but retaining the terms [(Ds)/( Dt)] and F, respectively. We find:
J = rQu+w |
Ds
Dt
|
+F×Ñs+ |
F
ro
|
Ñr¢×Ñs |
| (13) |
which is the PV flux written down in Marshall and Nurser (MN, Eqs.1b and 1c), but
modified by the non-advective thermobaric term [(F)/( ro)]Ñr¢×Ñs. The latter term (neglected in
MN) is zero if s = s( r¢) - for example
if we ignore the pressure dependence of r¢ on p. However,
if Ñr¢×Ñs ¹ 0, Lagrangian
conservation of PV no longer pertains even in adiabatic, frictionless flow -
see McDougall (1988). But more importantly in the present context, the flux
form of the PV equation can never be compromised.
In the remaining of this paper, unless noted otherwise, we will make use of the form (12) for J.
File translated from TEX by TTH, version 1.95.
On 19 Nov 2002, 16:02.