Geostrophic adjustment
Question: How does the adjustment to geostrophy take place?
H(1 + )
H(1 - )
0 x
h x( ) H H o sgn( )x
sgn( ) ,
x for x . for x
R S
T
1 1 0 0
u
t fv c
2 x
v
t fu c
2 y
t
u x
v
y 0
c gH
2
2 2 2
2
2
2 0
t c
x y f
F
H G I
K J
u-momentum
v-momentum continuity
(1)
(2)
(3)
Take the divergence of (1) and (2) and use (3) to eliminate the divergence
The nonrotating case (f = 0)
( , )x t 12 o sgn(x ct ) sgn( x ct ) u x t( , ) 12 co sgn(x ct ) sgn( x ct )
( , )x t 12 F x ct( ) G x ct( )
u x t( , ) c o
• the motion remains irrotational ( = 0) then
2 2
2 2
2
2
t c 2
x y
F
H G I K J
• assume two-dimensional motion / y 0
for | |x ct u x t( , ) 0 for | |x ct
time
The rotating case (f nonzero)
• The two moving steps will experience a Coriolis force acting to the right (NH case)
a flow in the -ve y-direction
this flow will experience a Coriolis force in the
-ve x-direction, opposing the flow in the x-direction
• The final adjusted state may be calculated without consideration of the transients
t f u
x
v
F
yH G IKJ
0The vorticity equation is
t
u x
v
y 0
The divergence equation is
t f
F
H G IKJ
0 the linearized form of the potential vorticity equation Q' f is the perturbation potential vorticity
t f
F
H G IKJ
0 Q x y t' ( , , ) Q x y' ( , , )0a parcel retains its original Q’
/ f o sgn( )x
2
2 2 2
2 2 2
t c
x f f o x
sgn( )
c d
dx f f o x
2 2
2 2 2
sgn( )
( ) exp( / ) exp( / )
x x L for x
x L for x
o R
R
R S
T
11 00LR c f/ gH f/ Here, at t = 0, u = v = 0 and = osgn(x)
2
2 2 2
2
2
2 0
t c
x y f
F
H G I
K J
in the steady state
Solution is:
Solution is continuous at and symmetric about x = 0
-3 -2 -1 1 2 3
1 2 3
-1 -2
-3
v
x
x
gH
( ) exp( / ) exp( / )
x x L for x
x L for x
o R
R
R S
T
11 00The transient problem
2 1
2 2 2
21 2
1 0
t c
x f H
1
0
o steady
o R
x
x L x at t
sgn( )
exp( / )sgn( ) .
1 exp[ (i kx t)]
2 f2 k c2 2
2
2 2 2
2 2 2
t c
x f f o x
sgn( )
The equation for is now
Put
Then Put
1 0 2 2
2
z
o k kx
ksin dk
Can show that at t = 0
LR1
sin(kx t) sin( kx t) 2sin kx cost
1 0 2 2
2
z
o k kx t
k dk
sin cos
/ y 0 1 sin kxcost For t > 0
Note that
When
u kx t v f
kx t
cos sin , cos cos
u gH c
kx t
k dk
o
z
2
0 2 2
cos sin
u gH c J f t x c for x ct
for x ct
o
R S
| T
|
( / ) ( (0 , / ) ) ,2 2 2 12
0
x
x
u u
x
v v
Energy considerations
t H u v gH
x gHu
y gHv
12 2 2 1
2 2 2 0
( ) ( ) ( )
E / t F X t( , ) F(X t, ) 0
E H u v gH dx
X
X
z
12 ( 2 2 ) 12 22
uH u
t fv c
F
x H G I
K J
2 0
vH v
t fu c
F
y H G I
K J
2 0
gH
t
u x
v
F
y H G I K J
0 +
+
/ y 0 x X
F x t( , ) gHu
E H u v gH dx
X
X
z
12 ( 2 2 ) 12 22F x t( , ) gHu
Total perturbation wave energy per unit length in the y-direction
KE PE
is the rate of energy transfer per unit length in the y-direction at x
Nonrotating situation
at t = 0: K E 0, PE 12 gH2o2
after the passage of the wave K E 12 H u2 12 H c2o2, P E 0 c gH
Rotating situation
Note: the Coriolis forces do not directly appear in the energy equations
However: the energy changes in the adjustment problem are profoundly affected by rotation
PE gH dx gH e dx
Hc L
o steady
X X
o
x L
o R
R
z
12 2 2 2 12 2 2z
023
2 2
1 1
( ) ( / )
KE Hv dx gH e dx
Hc L
steady X
X
o
x L
o R
R
z
12 2 12 2 2z
0 212
2 2
2
/ )
= one third of the PE released!
The remaining energy is radiated away as inertia-gravity waves in the transient part of the solution.
Notes:
• Energy is hard to extract from the available potential energy in a rotating fluid. The geostrophic equilibrium state which is established retains a certain
amount of the initial available potential energy.
• The steady solution is not one of rest, but is a geostrophically-balanced flow.
• The steady solution is degenerate - any velocity field in geostrophic balance satisfies the continuity equation exactly. Therefore the steady solution cannot be found by looking for a solution of the steady-state equations - some other information is needed.
• The required information is supplied by the principle that potential vorticity is materially conserved (or locally conserved in the linear problem).
• The equation for the steady solution contains a length scale, LR = c/|f|, where c = (gH)1/2 is the wave speed in the absence of rotation.
Rossby Length
LR is a fundamental length scale in atmosphere-ocean dynamics. It is the horizontal scale at which the gross effects of rotation are of comparable importance with gravitational (or buoyancy) effects
Early in the adjustment stage in the transient problem, the change in level is confined to a small distance and the horizontal pressure gradient is
comparatively large. Accordingly, gravity dominates the flow behaviour.
Thus at scales small compared with the Rossby length, the adjustment is approximately the same as in a non-rotating system.
Later, as the change in surface elevation is spread over a distance
comparable with the Rossby length, the Coriolis acceleration becomes just as important as the pressure gradient term and thus rotation leads to a response that is radically different from that in the non-rotating case.
LR is also an important scale for the geostrophic equilibrium solution as well.
In the problem analyzed, the discontinuity did not spread out indefinitely, but only over a distance of the order of LR.
For geostrophic, or quasi-geostrophic flow, LR is the scale for which the two contributing terms to the perturbation potential vorticity Q′ are of the same order.
For a sinusoidal variation of surface elevation with wavenumber k, the ratio of the vorticity term to the gravitational term in Q′ is .
Therefore, for short waves (1/k << LR) the vorticity term dominates while for long waves (1/k >> LR) gravitational effects associated with the free surface perturbation dominate.
/ f L k2R 2
For quasi-geostrophic wave motions, the ratio characterizes not only
the partition of perturbation potential vorticity, but also the partition of energy.
L k2R 2
Q' f
12
2 2
gH
z
wave dx2 2 R
KE L k PE
The perturbation PE integrated over one wavelength is
2 2
1
2 Hwave (u v )dx
The perturbation KE integrated over one wavelength is
For quasi-geostrophic waves periodic in x, one can show that (Ex. 9.2)
Thus short-wavelength geostrophic flow contains mainly
kinetic energy, whereas long-wavelength geostrophic flow has most of its energy in the potential form.
The situation is different for inertia-gravity waves (see Ex. 9.2).
Changes in are associated with changes in the mass field, whereas changes in are associated with changes in the velocity field.
For large scales, the potential vorticity perturbation is mainly associated with perturbations in the mass field and energy
changes are mainly a result of potential energy changes.
For small scales, the potential vorticity perturbation is
mainly associated with perturbations in the velocity field and energy changes are mainly a result of kinetic energy changes.
A distinction can be made between the adjustment processes at different scales.
At large scales (1/k >> LR) , it is the mass field that is
determined by the initial potential vorticity, and the velocity field is merely that which is in geostrophic equilibrium with the mass field.
In other words, the large-scale velocity field adjusts to be in equilibrium with the large-scale mass field.
In contrast, at small scales (1/k << LR) , it is the velocity field that is determined by the initial potential vorticity, and the mass field is merely that which is in geostrophic equilibrium with the velocity field.
In other words, the mass field adjusts to be in equilibrium with the velocity field.
Balanced adjustment
u u + u
x x
Sx
v
y
h + h
h f
Balanced adjustment
-a a F
initial state
a later state
u
-a a
x x Special cases
-a a initial state
0 u
x . u
addition of mass
H
S 0
v v
Horizontal cross-section through a field of cumulus clouds
heating function:
δ(y)sin (πz/H)
The response to an idealized slab-symmetric narrow sinusoidal heat source in a uniformly-stratified fluid layer.
Horizontal section at mid-height of the buoyancy perturbation at three times due to a periodic array of heat sources turned
the Gaussian heat source