1 / 47

Representing Gravity Current Entrainment in Large-Scale Ocean Models

Representing Gravity Current Entrainment in Large-Scale Ocean Models Robert Hallberg (NOAA/GFDL & Princeton U.) With significant contributions from Laura Jackson, Sonya Legg, and the Gravity Current Entrainment Climate Process Team:. NOAA/GFDL: S. Griffies, R. Hallberg, S. Legg, L. Jackson*

shaun
Download Presentation

Representing Gravity Current Entrainment in Large-Scale Ocean Models

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Representing Gravity Current Entrainment in Large-Scale Ocean Models Robert Hallberg (NOAA/GFDL & Princeton U.) With significant contributions from Laura Jackson, Sonya Legg, and the Gravity Current Entrainment Climate Process Team: NOAA/GFDL: S. Griffies, R. Hallberg, S. Legg, L. Jackson* NCAR: G. Danabasoglu, P. Gent, W. Large, W. Wu* U. Miami: E. Chassingnet, T. Ozgokmen, H. Peters, Y. Chang* WHOI: J. Price, J. Yang, U. Riemenschneider* Lamont Doherty: A. Gordon George Mason: P. Schopf Princeton U.: T. Ezer (Plus ~12 active collaborators) *Postdocs funded by the CPT http://www.cpt-gce.org

  2. An Idealized Rotating OverflowDOME Test Case 1 (Legg, et al., Ocean Modelling, 2006) Near-bottom tracer concentration with contours of buoyancy Tracer concentration just west of the inflow Dx=500m, Dz=30m MITgcm Simulation

  3. An Idealized Rotating OverflowDOME Test Case 1 (Legg, et al., Ocean Modelling, 2006) Near-bottom tracer concentration with contours of buoyancy Tracer concentration just west of the inflow Dx=500m, Dz=30m MITgcm Simulation

  4. Resolvable by large-scale models Hydraulic control at sill Geostrophic adjustment of plume along slope Downslope transport of dense water (some model types?) Some geostrophic eddy effects? Detrainment at neutral density Require Parameterization Exchange through subgridscale straits Shear instability and entrainment (TURBULENCE!!!) Bottom boundary layer mixing and drag processes (TURBULENCE!!!) Some eddy effects? Flow down narrow channels? z y x Important Processes in Overflows Shear instability & entrainment Hydraulic control at sill Geostrophic eddies Bottom friction Bottom-stress mixing Downslope descent Detrainment Physical processes in overflows

  5. Overview • A tour of overflows • Oceanic Gravity Currents are important in the formation and transformation of the majority of deep water masses. • Important Processes in Typical Oceanic Dense Gravity Currents: • Hydraulic or tidal control of source water flows, often in narrow straits • Downslope descent (gravitational, Ekman driven, and eddy induced) • Shear-driven mixing at the plume top • Bottom boundary layer mechanical stirring within the plume • Thermobaric influences of the ocean’s nonlinear equation of state • Detrainment at the neutral depth • Challenges for representing overflows in large-scale models: • Avoiding inherent problems with excessive numerical entrainment • Source water supply (representing the unresolved) • Studies of equilibrium stratified shear instability. • A new shear-driven turbulence mixing parameterization • A new bottom-turbulence mixing parameterization

  6. Mediterranean Outflow Plume • Without its 3-fold entrainment, Mediterranean Outflow water would fill the bottom of the Atlantic • Gibraltar itself exhibits rectified tidal exchange in conjunction with hydraulic control • Because of thermobaricity, salty Mediterranean water has a greater density at lower pressures, contributing to shallow detrainment. Gibraltar Velocities over the Tidal Cycle (CANIGO cruises Send & Baschek, JGR 2001)

  7. Climatological Salinity at 1000 m Depth

  8. Faroe Bank Channel and Denmark Strait Outflows Density along axis of Faroe Bank Channel Denmark Strait Sea Surface Temperature Denmark Strait: J. Girton; FBC: C. Mauritzen, J. Price

  9. Abyssal Overflows – the Romanche Fracture Zone Potential Temperature along Romanche Fracture Zone Potential Temperature at 5000 m Depth Ferron et al., JPO 1998

  10. z y x Steps in Adequately Representing Gravity Currents • Supply source water to the plume with the right rate and properties. • Model must be able to represent downslope flow without excessive numerical entrainment. • Parameterize entrainment & mixing to the right extent. • Parameterize subgridscale circulations? (e.g. eddies, flow in small channels). Shear instability & entrainment Hydraulic control at sill Geostrophic eddies Bottom friction Bottom-stress mixing Downslope descent Detrainment Physical processes in overflows

  11. Source Water properties depend on the right large-scale circulation and properties. Several important source waters enter through very narrow channels! Gibraltar is ~12 km wide. Red Sea outflow channel is ~5 km wide. Faroe Bank channel is ~15 km wide at depths that matter. Channels that are much smaller than the model grid require special treatment – e.g. partial barriers. Source water supply The topography around Gibraltar, with a 1° grid (black), and the coastline (blue) that GFDL’s 1° global isopycnal model uses.

  12. i , j æ ö h A 1 1 ç ÷ = + d - d q f v u i j ç ÷ i , j q q D D h è ø A h x y Representing Straits with Partially Open Faces(Work with A. Adcroft, GFDL) Partially open faces can dramatically improve simulations of overflows that pass through narrow straits. The model equations need to be modified to be energetically consistent. E.g. Sadourny’s 1975 Energy conserving discretization of the shallow water equations: Terms underlined in red are affected directly by using the partially open faces. Terms underlined in blue are affected indirectly (i.e. no code changes).

  13. Resolution requirements for avoiding numerical entrainment in descending gravity currents. Z-coordinate: Require that AND to avoid numerical entrainment. (Winton, et al., JPO 1998) Suggested solutions for Z-coordinate models: • "Plumbing" parameterization of downslope flow: Beckman & Doscher (JPO 1997), Campin & Goose (Tellus 1999). • Adding a separate, resolved, terrain-following boundary layer: Gnanadesikan (~1998), Killworth & Edwards (JPO 1999),Song & Chao (JAOT 2000). • Add a nested high-resolution model in key locations? • No existing scheme is entirely satisfactory! Sigma-coordinate:Avoiding entrainment requires that Isopycnal-coordinate:Numerical entrainment is not an issue - BUT • If resolution is inadequate, no entrainment can occur. Need

  14. Diapycnal Mixing Equations in Isopycnic Coordinates • In isopycnic coordinates, diapycnal diffusion is nonlinear • The discrete form leads to a coupled set of nonlinear differential equations These can be solved implicitly and iteratively, with an arbitrary distribution of diffusivities to avoid the impossible time-step limit (Hallberg, MWR 2000) • The work-diffusivity relationship is exact in density coordinates. • Entrainment can also be parameterized directly, based upon resolved shear Richardson numbers and a reinterpretation of the Ellis & Turner (1959) bulk Richardson number parameterization (Hallberg, MWR 2000). This parameterization gives entraining gravity currents that are qualitatively similar to observations, but has subsequently been improved upon.

  15. Constant Diffusivity Richardson Number Mixing

  16. DOME Model Intercomparisons and Resolution Dependence (Legg et al., Ocean Modelling 2006) MITgcm (Z-coordinate) with Convective Adjustment 2.5 km x 60 m 10 km x 144 m 50 km x 144 m HIM (isopycnal coordinate) with shear Ri# param. 10 km x 25 Layer 50 km x 25 Layer

  17. Plume Entrainment as a Function of Resolutionfor 6 DOME Test Cases Entrainment Rate Near Source (nondim.) Final Plume Buoyancy (m s-2) Horizontal Grid Spacing (km) Horizontal Grid Spacing (km) Solid lines: MITgcm (Z-coordinate) Dashed lines: HIM (Isopycnal coordinate) For full details, see Legg et al., Ocean Modelling 2006.

  18. Parameterizing Overflow Entrainment:Observations of Bulk Entrainment in Oceanicand Laboratory Gravity Currents (J. Price) A bulk entrainment law applies, provided the Reynolds number is not quite small.

  19. Examples of Gravity Current Mixing Parameterizations: • Generic shear parameterizations – e.g. KPP (Large et al., 1994): Typically calibrated for the Equatorial Undercurrent. • Two-equation turbulence closures (e.g. Mellor-Yamada; k-e; k-w). • Plume-specific parameterizations – e.g. Ellison & Turner (1959) bulk Ri# parameterization reinterpreted for shear Ri# (Hallberg, 2000): This can be cast as a diffusivity, D is over an unstable region: May Need Resolution Dependence!

  20. Simulated Mediterranean Outflow Plume(Papadakis et al., Ocean Modelling 2003) Zonal Velocity Salinity in 3 Isopycnal Layers Salinity

  21. A Non-rotating Overflow Entering a Stratified Environment (Courtesy T. Özgökmen)

  22. LES and Parameterized Overflow Entrainment(Xu, Chang, Peters, Özgökmen, and Chassignet, Ocean Modelling in press)

  23. Failure and Success of Existing Parameterizations • A universal parameterization can have no dimensional “constants”. • KPP’s interior shear mixing (Large et al., 1994) and Pacanowski and Philander (1982) both use dimensional diffusivities. • The same parameterization should work for all significant shear-mixing. • In GFDL’s HIM-based coupled model, Hallberg (2000) gives too much mixing in the Pacific Equatorial Undercurrent or too little in the plumes with the same settings. • To be affordable in climate models, must accommodate time steps of hours. • Longer than the evolution of turbulence. • Longer than the timescale for turbulence to alter its environment. • 2-equation (e.g. Mellor-Yamada, k-, or k-w) closure models may be adequate. • The TKE equations are well-understood, but the second equation (length-scale, or dissipation rate, or vorticity) tend to be ad-hoc (but fitted to observations) • Need to solve the vertical columns implicitly in time for: • TKE • Dissipation/vorticity • Stratification (T & S) • (and 5.) Shear (u & v) • Simpler sets of equations may be preferable. • Many use boundary-layer length scales (e.g. Mellor-Yamada) and are not obvious appropriate for interior shear instability. However, sensible results are often obtained by any scheme that mixes rapidly until the Richardson number exceeds some critical value.

  24. z x z Temperature (°C) x 3-DNS of Shear Instability(L. Jackson, R. Hallberg, & S. Legg in prep.) Temperature during initial development of Kelvin-Helmholtz instabilities Kelvin-Helmholtz instability 3D stratified turbulence Representative instantaneous along-channel Cross-section in statistical steady state Temperature (°C)

  25. Considerations for a Parameterization of StratifiedShear Instability S = ||∂U/∂z|| [s-1] Velocity shear N2 = -g/r∂r/∂z [s-2] Buoyancy Frequency H[m] Vertical extent of small Ri Q [m2 s-2] Turbulent kinetic energy per unit mass u* = (t/r)1/2 [m s-1] Friction velocity (for boundary turbulence) z* [m] Distance from boundary (for boundary turbulence) • Mixing should vanish if the shear Richardson number (Ri = N2/S2) exceeds ~1/4 everywhere • Vigorous mixing may extend past the region of small Ri. • Homogeneous stratified turbulence is often characterized by the buoyancy length scale • Kelvin-Helmholtz (K-H) saturation velocity scales are ~ H S. • K-H instabilities span the region of small Ri, i.e. length scales of ~ H. • Mixing-length arguments suggest peak K-H-type diffusivities scaling as ~H2S. • Near solid boundaries, length scales are proportional to the distance from the boundaries, and diffusivities are ~ 0.4u*z*.

  26. Translating “Entrainment Rate” parameterizations into diffusive parameterizations (L. Jackson) The diffusion of density can be linked to entrainment parameterizations by combining the density conservation equation: with the continuity equation in density coordinates: The latter equality is ill-behaved when ∂r/∂z=0, but with constant stratification it reduces to ET parameterisation (Hallberg, 2000)

  27. S=||Uz|| Entrainment-law derived theory for Shear-driven mixing k = 0 at solid boundaries • Properties: • Uses a length scale which is a combination of the width of the low Ri region (where F(Ri)>0) and the buoyancy length scale LBuoy = Q1/2/N. • Decays exponentially away from low Ri region • Vertically uniform, unbounded limit: • Ellison and Turner limit (large Q) reduces to form similar to ET parameterisation • Unstratified limit: similar to law-of-the-wall theories of parabolic diffusivity between two boundaries and log-like profiles of velocity near the boundaries.

  28. TKE Budget to Complement Proposed Diffusivity Equation • Assumptions: • Q reaches steady state faster than background flow is evolving so no DQ/Dt term • Assume Pr = 1 (for now) • Q0 needed to avoid singularity in diffusivity equation (solution not sensitive to Q0 and 0) • Parameterization of dissipation as c(Q-Q0)N • Q intended for use in diffusivity equation is due to turbulent kinetic energy only - difficult to compare to results from DNS because of internal waves.

  29. Equilibrium DNS of Shear-driven Stratified Turbulence • Non-hydrostatic direct numerical simulations (MITgcm) • 2m x 2m x 2.5m with grid size ~ 2.5mm in centre • Molecular viscosity and diffusivity, Kolmogorov scale mostly resolved. • Cyclic domain in x,y • Shear and jet profiles • Statistically steady state reached • Force average velocity profiles to evolve to given profile • Initially constant stratification and relaxed to initial density profile • All profiles are spatially averaged in x and y and time averaged

  30. z x z Temperature (°C) x 3-DNS of Shear Instability(L. Jackson, R. Hallberg, & S. Legg in prep.) Temperature during initial development of Kelvin-Helmholtz instabilities Kelvin-Helmholtz instability 3D stratified turbulence Representative instantaneous along-channel Cross-section in statistical steady state Temperature (°C)

  31. Buoyancy flux (m2/s3) DNS Shear-Instability Results and the Proposed Parameterization DNS data New parameterisation (Jackson et al.) ET parameterisation (Hallberg 2000) F(Ri) = 0.15*(1-Ri/0.25)/(1-0.9*Ri/0.25), c=1.9 F(Ri) = 0.15*(1-Ri/0.8)/(1+1.0*Ri/0.8), c=1.7 F(Ri) = 0.12*(1-Ri/0.25)/(1-0.9*Ri/0.25), c=1.24

  32. Buoyancy flux (m2/s3) DNS of Shear Instability and Existing 2-equation closures Black: DNS Results Green: GOTM k-w Blue: GOTM k-e Red: Mellor-Yamata 2.5 Jackson et al., proposed parameterization:

  33. Buoyancy flux (m2/s3) DNS Jet results DNS data New parameterisation (Jackson et al.) ET parameterisation (Hallberg 2000) F(Ri) = 0.15*(1-Ri/0.25)/(1-0.9*Ri/0.25), c=1.9 F(Ri) = 0.15*(1-Ri/0.8)/(1+1.0*Ri/0.8), c=1.7 F(Ri) = 0.12*(1-Ri/0.25)/(1-0.9*Ri/0.25), c=1.24

  34. Buoyancy flux (m2/s3) Existing 2-equation Closures Compared to DNS Jet Black: DNS Results Green: GOTM k-w Blue: GOTM k-e Red: Mellor-Yamata 2.5 Jackson et al., proposed parameterizations:

  35. Diapycnal Diffusivities Diagnosed from 3-D DNS Diagnosed diapycnal diffusivity (m2 s-1) Gradient Richardson number

  36. Buoyancy flux (m2/s3) Shear Instability with a Larger Ri#Not yet equilibrated?

  37. 500 m x 30 m MITgcm Illustrating the power of the CPT paradigm At the start of the CPT, with thick, nonrotating plumes entering ambient stratification, GFDL’s Isopycnal coordinate model (HIM) would give plumes that split in two. Such split plumes do not occur in nonhydrostatic “truth” simulations. 10 km x 25 layer HIM Ellison & Turner Mixing Only

  38. Observed profiles from Red Sea plume from RedSOX (H. Peters) Actively mixing Interfacial Layer Shear Ri# Param. Appropriate Here. Well-mixed Bottom Boundary Layer

  39. Bottom Boundary Layer Mixing • Diapycnal mixing of density requires work. • The rate at which bottom drag extracts energy from the resolved flow is straightforward to calculate. • Assumptions: • 20%? of the extracted energy is available to drive mixing. • Available work decays away from the bottom with e-folding scale of • Mixing completely homogenizes the near bottom water until the energy source is exhausted. Legg, Hallberg, & Girton, Ocean Modelling, 2006

  40. 500 m x 30 m MITgcm With thick plumes, both Interfacial and and Drag-induced Mixing are needed. (Legg et al., Ocean Modelling, 2006) 10 km x 25 layer HIM Ellison & Turner + Drag Mixing Ellison & Turner Mixing Only

  41. Year 5 salinity along 38.5°N in GFDL’s 1° Global Isopycnal Model Adding the Legg et al. bottom-drag mixing parameterization leads to dramatic improvements in an IPCC-class ocean model. Double Mediterranean plumes without bottom-drag mixing

  42. Summary • Overflows are critical in the formation of most deep-ocean water masses. • Turbulent mixing with the right rate is critical for models to obtain the right properties. (Otherwise in a stratified ambient environment, the plunging plume entrains the wrong water.) • Large-scale models require parameterizations of such mixing that capture both the equilibrium turbulence and (sometimes) its equilibrium modification of the resolved flow. • Mature Kelvin-Helmholtz-like mixing is significant in the interfacial layers atop gravity currents. • Existing parameterizations do not appear to work very well in detail based on comparisons with DNS (although they may work well enough for some overflows). • Laura Jackson (Princeton/GFDL CPT postdoc) has a new 2-equation (diffusivity – TKE) shear-driven turbulent mixing parameterization that looks very promising. • Bottom-stress driven turbulence is significant for homogenizing the bottom boundary layer, and must be parameterized. http://www.cpt-gce.org

More Related