Title: Marginal Thermobaric Stability in the Weddell Sea
1Marginal Thermobaric Stability in the Weddell Sea
- Miles McPhee
- McPhee Research Company
2Thermobaric Instability Following Loyning and
Weber, JGR, 102, p. 27875
2-layer system, upper layer colder and less saline
Linearized equation of state
Thermal expansion coefficient increases with depth
ambient
3Weddell Sea
Greenland Sea
4Strength of thermobaric tendency must exceed the
background stratification
Marginal stability line
5(No Transcript)
6(No Transcript)
7(No Transcript)
8(No Transcript)
9(No Transcript)
10ANZFLUX Ship CTD station 50, linearized
112.4 hour average of turbulence measurements
centered at time 206.35 (Warm Regime drift).
Circles are averages lines are twice the std dev
of the 15-min samples.
12(No Transcript)
13Two-layer (Type II) stability diagram following
Akitomo (1999) for idealized Ship Station 50
Increase Sml
14- The thermobaric barrier calculation
- Calculate the actual density (pressure included),
subtract density of a water column with mixed
layer properties. Determine the level (zmax) of
the maximum difference drmax. - Determine the sensible heat that must be vented
to reduce water temperature above zmax to ?ml. - Add the latent heat loss required to increase
salinity (by freezing) enough to eliminate ?? at
zml - Htot is the total heat loss.
15(No Transcript)
16Pentagrams indicate Hto t lt 100 MJ/m2
17(No Transcript)
18(No Transcript)
19(No Transcript)
20(No Transcript)
21(No Transcript)
2227.4 W m-2
McPhee, Kottmeier and Morison, JPO, 1999
23Ice temperatures from the AWI Buoy thermistor
string
24Mean values almost identical 30 W m-2
25Friction velocity prescribed from buoy results.
Heat flux also from buoy ice measurements but
scaled by the calculated ice thickness
Ocean heat flux calculated prognostically, ice
thickness determined by enthalpy balance at the
interface.
Dynamic mixed layer depth based on buoyancy
frequency (pot density). Scalar based on
difference from near surface value.
The model neglects thermobaric effects but
calculates thermobaric barrier parameters at each
time step.
26 The 1-D model forced with buoy data and
initialized with YU075 No thermobaricity effect
considered.
27 The 1-D model forced with buoy data and
initialized with YU075 Eddy viscosity set to 2000
cm2/s across vertical domain after 217.75.
28Horizontally Homogeneous Model Results
- Initialize model with every profile with Htot lt
100 MJ m-2 forced by buoy time series (38) - 27 1-D profiles became unstable by the end of
August
29Is there a simple way of getting a handle on eddy
viscosity and scalar diffusivity when
thermobaric mixing is occurring?
- Parameterize entrainment process in terms of
conversion of PE to TKE - Base the mixing length on a fraction (k) of the
entrained layer depth - Then
30(No Transcript)
31(No Transcript)
32(No Transcript)
33Summary
- Thermobaric instability was not observed directly
during the ANZFLUX 94 project but there is a
strong inference that it occurred shortly after. - About ¼ of the profiles observed with the yoyo
CTD system during the Maud Rise drift went
thermobarically unstable by the end of winter in
a simple 1-D model forced with drifting buoy
data. - In the model, preconditioning of the initial
density profile to include distinct step-like
structure in the upper pyncnocline was necessary
for instability. - Steps were found mostly in the halo region
surrounding Maud Rise (2500-3000 m isobaths)
34Summary (cont)
- There may be a good chance of encountering
episodes of Type II convection near Maud Rise in
late winter. - Measuring turbulent dissipation rates and
turbulent fluxes directly during a Type II
episode is feasible based on ANZFLUX experience.
Such data would be of great value in evaluating
and guiding numerical model development. - Even in the absence of direct Type II convection,
studying processes that maintain the step
structure and pycnocline weather in the Weddell
would add significantly to our understanding of
the system.