MAR650: Lecture 2: Layer Oceans and Mixed Layer Dynamics z y Wind stress o x Mixed Layer and Wind-induced Ekman layer Interior geostrophic layer Bottom Ekman layer Bottom log-boundary layer 1 1) Wind-induced surface Ekman layer ∂2u − fv = K E m 2 ∂z Coriolis force = Vertical diffusion ∂2v fu = K E m ∂z2 a) Current profile: τ v € The Ekman velocity decreases and rotates clockwise o E with depth Ekman spiral. 45 b) The Ekman layer thickness (depth): € € 2K Directly proportional to turbulent h = m E f viscosity coefficient and inversely hE UE (Ekman volume transport) proportional to the Coriolis parameter. c) The direction of the surface Ekman current: € v The angle between the wind stress and E o o tanα = = 1 , α = 45 surface Ekman current is 45 . On the uE τ v τ northern hemisphere, the surface E o Ekman current is 45 on the right of o o the wind stress. 45 90 d) The total volume transport: € € € U The volume transport is always 90o to the direction of the E wind stress. On the northern hemisphere, it is on the right of 2 the wind stress. Ekman Pumping: 0 ∂u ∂v ∂w ∂U ∂V 1 ∂τ y ( + + )dz = 0 ⇒ w = + = ∫ ∂x ∂y ∂z z=−hE ∂x ∂y ρf ∂x −hE If we extend it to a general case with x and y components of the wind, we get ∂U ∂V 1 ∂τ y ∂τ 1 w = + = ( − x ) = k ⋅ curlτ z= − hE ∂x ∂y ρf ∂x ∂y ρo f €The vertical velocity at the bottom of the surface Ekman layer is proportional to the vorticity of the surface wind stress and it is independent of the detailed structure of the Ekman flow. the vertical velocity is positive for the positive vorticity of the wind stress, and negative for the negative vorticity of the wind stress. curlτ > 0 € hE Ekman pumping 2) Interior geostrophic layer Turbulent mixing is weak, the motion is approximately geostrophic 1 ∂P 1 ∂P u = , v = − Coriolis force = Pressure gradient force g ρf ∂y g ρf ∂x FP FP (The pressure gradient force) t Geostrophic 0 t Vg 1 adjustment t P0 2 € € € Vg P1 Fc Fc P 2 F € € c Fc (The Coriolis force) € Where does the energy come here? Vorticity of the wind€ stress (Ekman pumping) € 3) The Bottom Ekman boundary layer 1 ∂P ∂ 2u fv K − = − + m 2 ρ ∂x ∂z 2 Coriolis force=Pressure gradient force+Vertical diffusion 1 ∂P ∂ v fu = − + K m 2 ρ ∂y ∂z Surface wind stress drives the surface Ekman layer and also transport the energy the interior layer. The energy driving the bottom boundary layer could be come from the interior geostrophic motion. a) Total velocity and Ekman velocity: Total velocity (a sum of geostrophic and Ekman velocities) rotates counterclockwise with depth. It reaches the geostrophic current at the top of the Ekman layer and vanishes at the bottom. Ekman velocity decreases and rotates clockwise with the height (from the bottom. It is equal in magnitude and opposite in direction to the geostrophic velocity at the bottom and vanish at the top of the bottom Ekman layer. v Total velocity Ekman velocity with depth with height 0 u -u 0 uE ug g (top) (bottom) d) Total volume transport: The volume transport is always 90o on the right of the bottom frictional force (on the northern hemisphere). Unlike the volume transport in the surface Ekman layer, the transport in the bottom boundary is related to the Ekman layer thickness and magnitude of the geostrophic current in the interior. 2 τ by K m hE 1 U E = − = − u g = − u g = − u g hE f fhE 2hE 2 Note: τby is the bottom stress exerting from the fluid on the bottom. The bottom frictional force is equal in magnitude and opposite to the bottom stress. τby (transport) UE -τby (bottom frictional force) 4) The bottom log boundary layer In a very thin layer close to the bottom, the turbulent viscosity-induced stress is much larger than the Coriolis and pressure gradient forces. In this layer, the stress does not change with depth: ∂τ This is also called the constant = 0 ∂z stress boundary layer According to the definition of the stress, we have € 2 ∂u ∂u 1 τ u* τ = ρK m ⇒ = ( ) = u* : the frictional velocity ∂z ∂z K m ρ K m Assume that mixing length l = κ z , where κ is von Karman constant , then we have ∂u u u u = * = * ⇒ u = * (ln z + C) ∂z l κz κ Define zo is the roughness height, at which turbulence-induced stress equals to zero, then u* z u = ln zo : 0.01 mm ~ 60 mm κ zo Surface Mixed Layer General features of vertical structure of the water temperature: T ∂T = 0 Mixed layer ∂z Solar radiation, wind mixing/cooling ∂T Thermoclines large Sharp vertical gradient of the water temperature Depth ∂z Deep weak ∂T stratified layer small Weak vertical gradient of the water temperature ∂z Question 1: Why do the thermoclines exist in the ocean? 8 An example: Assume that solar radiation leads to a vertically linear distribution of the water temperate: T = To + γz Wind mixing and cooling cause a vertical mixing, resulting in a mixed layer of h and a temperature as 1 0 h T (T z)dz T = ∫ o + γ = o − γ An averaged temperature from 0 to -h h −h 2 h h T = To − γ Mixed layer 2 Thermoclines A transition zone between the mixed T = To + γz layer to the deep ocean For density: pycnoclines; for salinity: thermohalines 9 Question 2: Is the oceanic current is a turbulent motion? Reynold number: VL V: velocity of the motion R = e γ L: typical scale of the motion γ: fluid’s viscosity Laboratory experiment: Re > 2000 or 4000 (depending on the property of the fluid), the motion become turbulent! In the ocean, 11 -6 2 Re > 10 (for example, in the Gulf Stream, V~ 1 m/s, L ~100 km, and γ ~10 m /s) The oceanic current is a turbulent motion! 10 Question 3: What causes turbulence and how is it dissipated? Assume that Shear Instability u2: the typical scale of the kinetic energy of the large-scale eddies; l: the typical size scale of the large-scale eddies; Large-scale eddies The time scale required for the energy transfer of large-scale eddies to small-scale eddies is l u Small-scale eddies The energy transfer rate from large-scale eddies to small-scale€ eddies equals: turbulent kinetic energy u 2 u 3 ~ = Dissipation time scale l / u l If these small-scale eddy’s energy is dissipated at a €rate of ε, then u 3 Turbulent dissipation rate ε ~ l 11 € Question 4: How does shear instability cause turbulence? Moving a blob of water in a stable density gradient leads to vertical ρ1 oscillations. ^ The buoyancy force defined as its “weight” equals to: ρ2 ^ Fg = −g(ρb − ρ s ) ρ3 ^ ρb is the density of the blob; ρs is the density of the surrounding water ρ4 The frequency of the blob’s oscillation is related to how fast it changes direction, It is controlled by the buoyancy force. In fact, the frequency of the oscillation is equal to g ∂ρ N = − Brunt-Väisälä frequency ρo ∂z Internal wave frequency: f ≤ ω ≤ N How could turbulence mixing occur in a stable stratified fluid? 12 One of the most important forces to cause vertical mixing is the vertical shear of the horizontal velocity, i.e,: Why do the shear of ∂u ∂v , the velocity is so ∂z ∂z important for mixing? Example: w > 0 w = 0 no wind w = 0 Smoke plume cigarette Does any shear of the velocity generate Light In the ocean: turbulence? Heavy 13 Richardson number The occurrence of the turbulence mixing is dampened by strong density gradient, but strengthened by strong shear of the velocity. Whether or not turbulence mixing could occur depends on a ratio of vertical stratification to the vertical shear of the velocity. Richardson number is defined as a ratio of the square of the buoyancy frequency to the square of the vertical shear of the horizontal velocity: 2 2 ∂u Ri = N / ∂z In general, turbulence develops as Ri < 0.25, but recent microstructure measurement shows that critical Richardson number is bigger for the interior ocean mixing or over the slope. Question 5: How is a mixed layer developed? What are the key balance of turbulence energy during the development of a mixed layer? 14 A mixed layer is developed with 4 stages: t (time) t h1 t1/3 h h2 t2 t1/3 h3 Monin-Obukov Depth h4 Stage 1: t << 1/f (time scale << the inertial period), wind stress is limited in a thin layer close to the surface, h ~ t Stage 2: the mixed layer continues to deepen, but cooler water below the mixed layer is entrained, the deepening speed of the mixed layer decreases: h ~ t1/3 Stage 3: a strong vertical shear of the velocity is established, shear instability speed up mixing. In this stage, h ~ t2 Stage 4: Energy generated by the vertical shear of the velocity is balanced by the heat-induced buoyancy force, the mixed layer stops deepening. This is a slowly process during which h ~ t1/3. u* = τ ο / f : frictional velocity; ρ : reference density ~ 2 m u 3 o ρo o * : surface wind stress, Q: the surface heat flux, Monin-Obukov Depth: h = τ ο 15 gαQ mo: turbulent dissipation coefficient € € The most important properties of a mixed layer is ∂V ∂T ∂ S ∂ρ = 0; = 0, = 0, = 0 ∂ z ∂ z ∂ z ∂ z Question 6: For a given same wind stress, the wind could drive an Ekman flow, but also wind mixing leads €to a mixed layer.