The Lagrangian-Averaged Navier-Stokes alpha (LANS-α) Turbulence Model in Primitive Equation Ocean Modeling
|
|
- Mercy Marion Merritt
- 5 years ago
- Views:
Transcription
1 The Lagrangian-Averaged Navier-Stokes alpha (LANS-α) Turbulence Model in Primitive Equation Ocean Modeling Mark R. Petersen with Matthew W. Hecht, Darryl D. Holm, and Beth A. Wingate Los Alamos National Laboratory Outline POP ocean model & climate change assessment LANS-α implementation in POP Idealized test case: the channel domain The real thing: the North Atlantic D2Hfest LA-UR July 24, 2007
2 CO2 ice core record currently: 380 ppm Homo sapiens neanderthalensis Homo sapiens archaic Homo sapiens sapiens Homo erectus domesticated plants & animals Siegenthaler et.al. Science 2005
3 Community Climate System Model Collaboration of: National Center for Atmospheric Research (NCAR) in Boulder, CO Los Alamos National Laboratory (LANL) Atmosphere (NCAR) Land Surface (NCAR) Flux Coupler Ocean POP (LANL) Sea Ice (LANL)
4 IPCC - Intergovernmental Panel on Climate Change Created in 1988 by World Meteorological Organization (WMO) and United Nations Environment Programme (UNEP) Role of IPCC: assess on a comprehensive, objective, open and transparent basis the scientific, technical and socio-economic information relevant to understanding: the scientific basis of risk of human-induced climate change its potential impacts and options for adaptation and mitigation. Main activity: Assessment reports Third Assessment Report: 2001 Fourth Assessment Report: 2007 Fifth: planned for 2013
5 IPCC scenarios of future emissions 1: global population levels off, declines after : continuously increasing population A: slower conversion to clean & efficient technologies A1FI: fossil intensive A1T: non-fossil intensive A1B: balance of F&T A2 B: faster conversion to clean & efficient technologies B1 B2 economic models IS92a: business as usual (extrapolation from current rates of increase) carbon cycle models
6 scenarios CO 2 emissions CO 2 concentration A: slower conversion to clean & efficient technologies B: faster conversion to clean & efficient technologies 1: global population levels off, declines after : continuously increasing population A1FI: fossil intensive A1T: non-fossil intensive A1B: balance of F&T A2 B1 B2 economic models carbon cycle models Final product average of ensemble Temperature change scenario A2 ensemble of climate models
7 Ensemble mean: temperature changes
8 IPCC: Estimates of confidence TAR p.82
9 Parallel Ocean Program (POP) Bryan-Cox type model, z-level vertical grid, finite difference model conservation of momentum " t u + u# $u % f & u = %' 0 %1 $p + A M $ h 2 u + " z µ" z u advection conservation of mass for incompressible fluid advection Coriolis " h # u + $ z w = 0 pressure gradient conservation of tracers (temperature, salinity) " t # + u$ %# = A H % h 2 # + " z &" z # + Q diffusion diffusion source/ sink grid u w ϕ t p ρ 0 T S hor. velocity vertical velocity tracer time pressure density temperature salinity hydrostatic in the vertical "p "z = #$g equation of state " = "( T,S, p)
10 POP: 0.1º resolution, speed
11 Parallel Ocean Program (POP) Resolution is costly, but critical to the physics Climate simulations low resolution: 1 deg (100 km) long duration: 100s of years fully coupled to atmosphere, etc. Eddy-resolving sim. high resolution: 0.1 deg (10 km) short duration: years ocean only Rossby Radius of deformation Potential temperature 0.8º x 0.8º grid Potential temperature 0.1º x 0.1º grid
12 low resolution: 0.8º What do you get with higher resolution? Small-scale turbulence and eddies transport energy and heat. cost of doubling horizontal grid is factor of º Reynolds decomposition: u = u + u " total 0.2º time average perturbation These become more realistic with higher resolution: eddy heat transport: u " 2 T " eddy kinetic energy: ( 1 2 u " u " ) 2 feedback of small-scale features on the largescale mean flow - important for oceanic jets vertical temperature profile high resolution: 0.1º
13 Outline POP ocean model & climate change assessment LANS-α implementation in POP Idealized test case: the channel domain The real thing: the North Atlantic
14 Lagrangian-Averaged Navier-Stokes Equation (LANS-α) Two ways to take averages: Lagrangian and Eulerian Lagrangian averaged velocity actual velocity v Lagrangian averaged trajectory particle trajectory Eulerian averaged velocity, u, is average at x o v Lagrangian averaged velocity u Eulerian averaged velocity smooth u = ( 1"# 2 $ 2 ) "1 v Helmholtz operator rough
15 Lagrangian-Averaged Navier-Stokes Equation (LANS-α) v rough u smooth larger α smooths more v Lagrangian averaged velocity u Eulerian averaged velocity smooth u = ( 1"# 2 $ 2 ) "1 v Helmholtz operator rough "v dt + u# $v + v j$u j % f & u = % 1 ' 0 $p + ($ 2 v + F advection extra nonlinear term Coriolis pressure gradient diffusion
16 Lagrangian-Averaged Navier-Stokes Equation (LANS-α) "v dt + u# $v + v j$u j % f & u = % 1 ' 0 $p + ($ 2 v + F advection extra nonlinear term Coriolis pressure gradient diffusion $ v 1 # x u 1 + v 2 # x u 2 + v 3 # x u 3 ' & ) v j "u j = v 1 # y u 1 + v 2 # y u 2 + v 3 # y u 3 & ) % v 1 # z u 1 + v 2 # z u 2 + v 3 # z u 3 ( Extra nonlinear term is in alpha model, but not in Leray model. This term is required for conservation of PV.
17 tracer equation momentum equation tracer equation momentum equation advection Standard POP "# dt + u$ %# = D H (#) + D V (#) advection extra nonlinear term advection Coriolis advection Coriolis diffusion "u dt + u# $u % f & u + metric(u) = %' %1 0 $p + F H (u) + F V (u) Helmholtz inversion u = ( 1"# 2 $ 2 ) "1 v e.g. centrifugal POP-alpha "# dt + u$ %# = D H (#) + D V (#) diffusion e.g. centrifugal pressure gradient pressure gradient rough velocity, diffusion smooth velocity, "v dt + u# $v + u j$v j % f & u + metric(u) = %' %1 0 $p + F H (v) + F V (v) v u diffusion
18 The POP-alpha model Issues: 1. How do we implement the alpha model within the barotropic/baroclinic splitting of POP? full ocean (vertical section) level 1 level 2 level 3 level K fast surface gravity waves vertically integrated barotropic part single layer implicit time step slower internal gravity waves baroclinic part multiple levels explicit time step vertical mean = 0 level 1 level 2 level 3 level K
19 Outline of algorithm - Standard POP Start: variables known at step n+1 temperature, salinity, density, pressure at n n baroclinic velocities u k barotropic velocity U n Baroclinic computations levels k =1Kkm leap frog step: n u ˆ +1 k = u n"1 n k + 2#tRHS k where RHS contains: advection metric (centrifugal) Coriolis diffusion Barotropic computations implicit solve for the free surface height " n +1 new barotropic velocity: U n +1 = U n"1 + 2#tRHS where RHS contains: vertically integrated forcing "# n$1,"# n,"# n +1 subtract depth average: n u +1 n k = u ˆ +1 k " 1 H km # k=1 n u ˆ +1 k dz k add barotropic velocity back to baroclinic: n u +1 n k = u +1 k + U n +1 step n+1 complete
20 Outline of algorithm - POP-alpha Start: variables known at step n+1 temperature, salinity, density, pressure at n baroclinic smooth and rough velocities u n n k,v k barotropic smooth and rough velocities U n,v n Baroclinic computations levels k =1Kkm leap frog step: n v ˆ +1 k = v n"1 n k + 2#tRHS k where RHS contains: advection metric (centrifugal) Coriolis diffusion extra nonlinear term, subtract depth average: n v +1 n k = v ˆ +1 k " 1 H km # k=1 n u +1 n k = smooth(ˆ v +1 k ) n u +1 n k = u ˆ +1 k " 1 H km # k=1 "u T # v n v ˆ +1 k dz k n u ˆ +1 k dz k Barotropic computations implicit solve for the free surface height " n +1 new barotropic velocity: V n +1 = V n"1 + 2#tRHS where RHS contains: vertically integrated forcing "# n$1,"# n,"# n +1 U n +1 = smooth(v n +1 ) add barotropic velocity back to baroclinic: n u +1 n k = u +1 k + U n +1 n v +1 n k = v +1 k + V n +1 step n+1 complete
21 Barotropic Algorithm - Pop-alpha Simultaneously solve for: free surface height " n +1 and both velocities, U n +1,V n +1 Invert using iterative CG routine smoothing within each iteration is too costly! (" # H(1$% 2 " 2 ) $1 " $ 2 &g' )( n +1 = 2 1 $ % "#g HUn&1 & 2 ' n + 1 $ % #g" 2 #g H(1&(2 $ 2 ) &1 G n & BU n & #g$ ' n&1 + ' n V n +1 = V n"1 + 2#t[ G n " BU n " $g% (& n"1 + & n + & n +1 )] [ ( )] momentum forcing terms U n +1 = (1"# 2 $ 2 ) "1 V n +1 momentum forcing terms resolution:
22 (" # H(1$% 2 " 2 ) $1 " $ 2 &g' )( n +1 = 2 Barotropic Algorithm - Pop-alpha Simultaneously solve for: free surface height Invert using iterative CG routine 1 $ % "#g HUn&1 & 2 ' n + 1 $ % #g" 2 #g H(1&(2 $ 2 ) &1 G n & BU n & #g$ ' n&1 + ' n V n +1 = V n"1 + 2#t[ G n " BU n " $g% (& n"1 + & n + & n +1 )] " n +1 and both velocities, U n +1,V n +1 smoothing within each iteration is too costly! What if we eliminate just this one smoothing step? [ ( )] momentum forcing terms U n +1 = (1"# 2 $ 2 ) "1 V n +1 momentum forcing terms resolution:
23 Outline POP ocean model & climate change assessment LANS-α implementation in POP Idealized test case: the channel domain The real thing: the North Atlantic
24 up N E The test problem: Idealization of Antarctic Circumpolar Current surface 12ºC thermal forcing 2ºC periodic bndry solid boundary solid boundary zonal wind periodic bndry Surface temperature deep-sea ridge isotherms Baroclinic Instability: 1.Eastward zonal wind causes northward Ekman transport 2.Circulation tilts isotherms (lines of constant temperature) 3.Potential energy of baroclinic instability converted to kinetic energy. 4.Small scale turbulence (eddies) transport heat and kinetic energy Thermocline depth is determined by the eddy transport quantities.
25 The Baroclinic Instability Surface temperature Kinetic energy POP POP POP low resolution-0.8º standard POP high resolution-0.1º standard POP resolution: 0.2 (low) (high) Eddy kinetic energy POP Potential temperature - vertical cross section POP POP resolution: 0.2 (low) (high) Depth of 6C isotherm high res. depth low resolution-0.8º standard POP depth high resolution-0.1º standard POP med. res. S N S N low res.
26 POP-alpha Results Kinetic energy Eddy kinetic energy POP-α POP-α POP POP-α POP POP POP POP-α POP POP resolution: 0.8º 0.8º 0.4º 0.4º 0.2º Depth of 6C isotherm 0.2 POP resolution: 0.8º 0.8º 0.4º 0.4º 0.2º Computation time (log scale) POP 0.4 POP-α POP POP-α 0.4 POP 0.8 POP-α POP POP-α 0.8 POP resolution: 0.8º 0.8º 0.4º 0.4º 0.2º
27 Can you see more eddies with LANS-alpha? POP, 0.4º resolution Pot. Temperature POP-α, 0.4º resolution Pot. Temperature POP, 0.2º resolution Pot. Temperature Velocity Velocity Velocity All sections at a depth of 1600m rough velocity: red, smooth: black
28 Dispersion Relation for LANS-α using linearized shallow water equations normally: with LANS-α: Gravity waves " 2 = k 2 gh " 2 = k 2 gh 1+ # 2 k 2 frequency ω, wavenumber k, gravity g, height H Dispersion relation normal with LANS-α Rossby waves #k$ normally: " = k 2 +1/R 2 LANS-α: " = normal with LANS-α #k$ ( ) +1/R 2 k 2 1+ % 2 k 2 Rossby radius R = gh / f 0 beta " = # y f Dispersion relation LANS-α slows down gravity and Rossby waves at high wave number.
29 What does LANS-α do to the Rossby Radius? Solve for k R, the wavenumber of the Rossby Radius: Use that to find R*, the effective Rossby Radius, as a function of α: Dispersion relation k R k R * normal LANS-α, small α LANS-α, larger α LANS-α makes the Rossby Radius effectively larger.
30 The POP-alpha model Issues: 3. How do we smooth the velocity in an Ocean General Circulation Model? v rough u smooth Helmholtz inversion is costly! ( ) "1 v u = 1"# 2 $ 2 or: use a filter u(x) = # G(r)v(x " r)dr for example, a top-hat filter
31 1D Filter Instabilities filter, G(r) u i = bv i"1 + v i + bv i b If b = 0.5, smoothing filters out the Nyquist frequency. gridpoints v rough u smooth u i = bv i"1 + v i + bv i b = 1 ("1) ("1) = 0 The smooth velocity u is blind to this oscillation. Therefore, the free surface height " cannot counter it!
32 1D Filter Instabilities filter, G(r) u i = bv i"1 + v i + bv i b If b = 0.5, smoothing filters out the Nyquist frequency. v rough 370 m/s u smooth 30 m/s Condition for stability is b < 0.5
33 Filter: Conditions for stability rough velocity v filter: b < 1/2 stable b stable b stable b b < 1/2+c b+c < 1 c < 1/2 stable b,c b+d < 1/2+c b+c < 1+2d c < 1/2 d < 1/2+b b+d < 1/2+c+e b+c+e < 1+2d c < 1/2+e d+2e < 1+b
34 Filter analysis: Helmholtz inversion Green s function Take v to be a point source: Then compute u = (1" # $! Can use this to understand filter near boundaries: 2 2 "1 ) v
35 Filters filter width 3 filter width 5 Wider filters result in: Stronger smoothing Effects are like larger α More computation More ghostcells Temperature - hor. mean filter width 7 high res POP filter width 9 med. res POP
36 Filters filter width 3 filter width 5 Wider filters result in: Stronger smoothing Effects are like larger α More computation More ghostcells Temperature - hor. mean filter width 7 high res POP filter width 9 med. res POP med. res POP-α width 9 width 7 width 5 width 3
37 Helmholtz inversion: vary alpha Kinetic energy Filters: vary the filter width Kinetic energy Eddy kinetic energy Eddy kinetic energy
38 Adding LANS-α increases computation time by <30% We can take smaller timesteps with LANS-α
39 How does LANS-α compare with other turbulence models? POP Hypervisc. only POP-α Hypervisc, LANS- α GM Hypervisc, Gent-McW. GM-α Hypervisc, Gent-McW, LANS-α 2ν Hypervisc, 2x viscosity coef. 4ν Hypervisc, 4x viscosity coef. Kinetic Energy POP POP-α POP GM-α 2ν 4ν GM Depth of 6C isotherm POP, higher res. GM GM-α POP-α POP POP-α 0.4 resolution Eddy Kinetic Energy 0.2 res. POP 2ν 4ν POP GM-α 2ν 4ν GM 0.4 resolution 0.2 res.
40 Leray Model is about half as strong as LANS-α Kinetic Energy Depth of 6C isotherm Eddy Kinetic Energy
41 Outline POP ocean model & climate change assessment LANS-α implementation in POP Idealized test case: the channel domain The real thing: the North Atlantic
42 POP simulations of the North Atlantic POP, 0.28º resolution Sea surface height POP, 0.1º resolution Sea surface height More realistic Gulf Stream and Northwest corner at high resolution Smith et. al jpo
43 POP simulations of the North Atlantic POP, 0.28º resolution Eddy kinetic energy POP, 0.1º resolution Eddy kinetic energy Higher eddy kinetic energy at high res. Smith et. al jpo
44 LANS-α in POP: North Atlantic simulations The punch lines: 1. Due to rough boundaries and fast jets in realistic domains, both POP-α and Leray have numerical instabilities in the North Atlantic, particularly at boundaries. 2. The Leray Model runs longer in the North Atlantic, and has higher EKE and visibly more vortices (but still boundary issues). 3. The key is to use the right boundary conditions for the smoothing step (Helmholtz or filter). spurious oscilations in rough vel. filter weights: filter weights:
45 LANS-α in POP: North Atlantic simulations Even with lower filter weights, numerical instabilities at boundary stop code after 6000 time steps (70 days):
46 Leray Model runs stably in North Atlantic "v dt + u# $v + v j$u j % f & u = % 1 extra nonlinear term not in Leray model ' 0 $p + ($ 2 v + F POP 0.2º 3yr mean temp. & velocity POP-Leray 0.2º 3yr mean temp. & velocity Leray Model has visibly more vortices, but still boundary issues.
47 POP-Leray has higher KE and EKE than POP POP 0.2º Eddy kinetic energy (3yr mean) POP-Leray 0.2º Eddy kinetic energy (3yr mean) cross section through North Atlantic current Globally averaged EKE (3 yr mean): POP 0.2º: 11.1 POP-Leray 0.2º: 13.1 POP 0.1º: 29.4 (Smith et. al. 2000) POP-Leray
48 POP-Leray has higher KE and EKE that POP POP 0.2º Eddy kinetic energy POP 0.1º EKE, varying viscosity high visc. POP-Leray 0.2º Eddy kinetic energy more realistic low visc. (Bryan et. al. 2007)
49
50 What are my boundary conditions? B.C. v = 0 u = 0 I ve tried many! Equation " t v + u# $v + v j $u j % f & u = % 1 ' 0 $p + ($ 2 v + F u = ( 1"# 2 $ 2 ) "1 v u = filter(v) Option 1: shrink filter at boundary Option 2: shrink filter near boundary Option 3: make filter weights=0 on land water land water land water land
51 A Possibility: Use variable alpha " t v + u# $v + $u 2 $% 2 (x) + v j $u j & f ' u = & 1 ( 0 $p + )$ 2 v + F u = ( 1" 2#(x)$#(x) % $ "# 2 $ 2 ) "1 v We are thinking about this
52 Summary Higher resolution can solve all of your problems. You can t possibly have high enough resolution to solve your problems. The LANS-alpha model captures higher-resolution effects in our test problem, where eddies near the grid-scale are important. POP-Leray runs longer than POP- α in the North Atlantic domain, and shows promising signs, like higher eddy activity Both POP-α and POP-Leray have problems with the rough boundaries and topography of the North Atlantic. Further work on boundary conditions for LANS-alpha, with Helmholtz or filter smoothing, is required.
2. Baroclinic Instability and Midlatitude Dynamics
2. Baroclinic Instability and Midlatitude Dynamics Midlatitude Jet Stream Climatology (Atlantic and Pacific) Copyright 26 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without
More informationGoals of this Chapter
Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature in the presence of positive static stability internal gravity waves Conservation of potential vorticity in the presence
More informationATMOSPHERIC AND OCEANIC FLUID DYNAMICS
ATMOSPHERIC AND OCEANIC FLUID DYNAMICS Fundamentals and Large-scale Circulation G E O F F R E Y K. V A L L I S Princeton University, New Jersey CAMBRIDGE UNIVERSITY PRESS An asterisk indicates more advanced
More information1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must
Lecture 5: Waves in Atmosphere Perturbation Method With this method, all filed variables are separated into two parts: (a) a basic state part and (b) a deviation from the basic state: Perturbation Method
More informationDaniel J. Jacob, Models of Atmospheric Transport and Chemistry, 2007.
1 0. CHEMICAL TRACER MODELS: AN INTRODUCTION Concentrations of chemicals in the atmosphere are affected by four general types of processes: transport, chemistry, emissions, and deposition. 3-D numerical
More information6 Two-layer shallow water theory.
6 Two-layer shallow water theory. Wewillnowgoontolookatashallowwatersystemthathastwolayersofdifferent density. This is the next level of complexity and a simple starting point for understanding the behaviour
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationRotating stratified turbulence in the Earth s atmosphere
Rotating stratified turbulence in the Earth s atmosphere Peter Haynes, Centre for Atmospheric Science, DAMTP, University of Cambridge. Outline 1. Introduction 2. Momentum transport in the atmosphere 3.
More informationJet Formation in the Equatorial Oceans Through Barotropic and Inertial Instabilities. Mark Fruman
p. 1/24 Jet Formation in the Equatorial Oceans Through Barotropic and Inertial Instabilities Mark Fruman Bach Lien Hua, Richard Schopp, Marc d Orgeville, Claire Ménesguen LPO IFREMER, Brest, France IAU
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More informationResolution-dependent Eddy Parameterizations for Large-scale Ocean Models
Resolution-dependent Eddy Parameterizations for Large-scale Ocean Models Robert Hallberg with contributions from Alistair Adcroft NOAA / GFDL and Princeton University Mercator/Tripolar Resolution Required
More informationNote that Rossby waves are tranverse waves, that is the particles move perpendicular to the direction of propagation. f up, down (clockwise)
Ocean 423 Rossby waves 1 Rossby waves: Restoring force is the north-south gradient of background potential vorticity (f/h). That gradient can be due to either the variation in f with latitude, or to a
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationGravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017
Lecture 5: Waves in Atmosphere Perturbation Method Properties of Wave Shallow Water Model Gravity Waves Rossby Waves Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature
More information2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves
Chapter. The continuous equations φ=gh Φ=gH φ s =gh s Fig..5: Schematic of the shallow water model, a hydrostatic, incompressible fluid with a rigid bottom h s (x,y), a free surface h(x,y,t), and horizontal
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationR. Hallberg, A. Adcroft, J. P. Dunne, J. P. Krasting and R. J. Stouffer NOAA/GFDL & Princeton University
Sensitivity of 21st Century Steric Sea Level Rise to Ocean Model Formulation R. Hallberg, A. Adcroft, J. P. Dunne, J. P. Krasting and R. J. Stouffer NOAA/GFDL & Princeton University Hallberg, R., A. Adcroft,
More informationThe general circulation: midlatitude storms
The general circulation: midlatitude storms Motivation for this class Provide understanding basic motions of the atmosphere: Ability to diagnose individual weather systems, and predict how they will change
More informationNorth Atlantic circulation in three simulations of 1/12, 1/25, and 1/50
North Atlantic circulation in three simulations of 1/12, 1/2, and 1/ Xiaobiao Xu and Eric Chassignet Center for ocean-atmospheric prediction studies Florida State University Motivation Numerical models
More informationIsland Wakes in Shallow Water
Island Wakes in Shallow Water Changming Dong, James C. McWilliams, et al Institute of Geophysics and Planetary Physics, University of California, Los Angeles 1 ABSTRACT As a follow-up work of Dong et al
More informationAtmosphere, Ocean and Climate Dynamics Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 12.003 Atmosphere, Ocean and Climate Dynamics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Problem
More information10 Shallow Water Models
10 Shallow Water Models So far, we have studied the effects due to rotation and stratification in isolation. We then looked at the effects of rotation in a barotropic model, but what about if we add stratification
More informationA comparison of two formulations of barotropic-baroclinic splitting for layered models of ocean circulation
A comparison of two formulations of barotropic-baroclinic splitting for layered models of ocean circulation Robert L. Higdon Department of Mathematics Oregon State University Corvallis, Oregon 97331-4605,
More information8 3D transport formulation
8 3D transport formulation 8.1 The QG case We consider the average (x; y; z; t) of a 3D, QG system 1. The important distinction from the zonal mean case is that the mean now varies with x, or (more generally)
More informationTorben Königk Rossby Centre/ SMHI
Fundamentals of Climate Modelling Torben Königk Rossby Centre/ SMHI Outline Introduction Why do we need models? Basic processes Radiation Atmospheric/Oceanic circulation Model basics Resolution Parameterizations
More informationisopycnal outcrop w < 0 (downwelling), v < 0 L.I. V. P.
Ocean 423 Vertical circulation 1 When we are thinking about how the density, temperature and salinity structure is set in the ocean, there are different processes at work depending on where in the water
More informationDynamics Rotating Tank
Institute for Atmospheric and Climate Science - IACETH Atmospheric Physics Lab Work Dynamics Rotating Tank Large scale flows on different latitudes of the rotating Earth Abstract The large scale atmospheric
More informationModeling Large-Scale Atmospheric and Oceanic Flows 2
Modeling Large-Scale Atmospheric and Oceanic Flows V. Zeitlin known s of the Laboratoire de Météorologie Dynamique, Univ. P. and M. Curie, Paris Mathematics of the Oceans, Fields Institute, Toronto, 3
More informationOCEAN MODELING II. Parameterizations
OCEAN MODELING II Parameterizations Gokhan Danabasoglu Oceanography Section Climate and Global Dynamics Division National Center for Atmospheric Research NCAR is sponsored by the National Science Foundation
More informationFrictional Damping of Baroclinic Waves
Frictional Damping of Baroclinic Waves HHH, 26th May 2006 Bob Plant With thanks to Stephen Belcher, Brian Hoskins, Dan Adamson Motivation Control simulation, T+60 Simulation with no boundary layer turbulence,
More informationA Two-Level Time-Stepping Method for Layered Ocean Circulation Models: Further Development and Testing
A Two-Level Time-Stepping Method for Layered Ocean Circulation Models: Further Development and Testing Robert L. Higdon Department of Mathematics Oregon State University Corvallis, Oregon 97331-4605 This
More informationPAPER 333 FLUID DYNAMICS OF CLIMATE
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 1:30 pm to 4:30 pm Draft 21 June, 2016 PAPER 333 FLUID DYNAMICS OF CLIMATE Attempt no more than THREE questions. There are FOUR questions in total.
More informationOcean General Circulation Modeling
Lecture 8 Ocean General Circulation Modeling 9. The equations of motion: Navier-Stokes The governing equations for a real fluid are the Navier-Stokes equations (conservation of linear momentum and mass
More informationthat individual/local amplitudes of Ro can reach O(1).
Supplementary Figure. (a)-(b) As Figures c-d but for Rossby number Ro at the surface, defined as the relative vorticity ζ divided by the Coriolis frequency f. The equatorial band (os-on) is not shown due
More informationd v 2 v = d v d t i n where "in" and "rot" denote the inertial (absolute) and rotating frames. Equation of motion F =
Governing equations of fluid dynamics under the influence of Earth rotation (Navier-Stokes Equations in rotating frame) Recap: From kinematic consideration, d v i n d t i n = d v rot d t r o t 2 v rot
More informationROSSBY WAVE PROPAGATION
ROSSBY WAVE PROPAGATION (PHH lecture 4) The presence of a gradient of PV (or q.-g. p.v.) allows slow wave motions generally called Rossby waves These waves arise through the Rossby restoration mechanism,
More informationAtmosphere, Ocean and Climate Dynamics Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 12.003 Atmosphere, Ocean and Climate Dynamics Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Contents
More informationLecture 1. Equations of motion - Newton s second law in three dimensions. Pressure gradient + force force
Lecture 3 Lecture 1 Basic dynamics Equations of motion - Newton s second law in three dimensions Acceleration = Pressure Coriolis + gravity + friction gradient + force force This set of equations is the
More informationSIO 210 Introduction to Physical Oceanography Mid-term examination November 3, 2014; 1 hour 20 minutes
NAME: SIO 210 Introduction to Physical Oceanography Mid-term examination November 3, 2014; 1 hour 20 minutes Closed book; one sheet of your own notes is allowed. A calculator is allowed. (100 total points.)
More informationInfluence of forced near-inertial motion on the kinetic energy of a nearly-geostrophic flow
Abstract Influence of forced near-inertial motion on the kinetic energy of a nearly-geostrophic flow Stephanne Taylor and David Straub McGill University stephanne.taylor@mail.mcgill.ca The effect of forced
More informationTurbulence. 2. Reynolds number is an indicator for turbulence in a fluid stream
Turbulence injection of a water jet into a water tank Reynolds number EF$ 1. There is no clear definition and range of turbulence (multi-scale phenomena) 2. Reynolds number is an indicator for turbulence
More information3. Midlatitude Storm Tracks and the North Atlantic Oscillation
3. Midlatitude Storm Tracks and the North Atlantic Oscillation Copyright 2006 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without permission. EFS 3/1 Review of key results
More informationpsio 210 Introduction to Physical Oceanography Mid-term examination November 3, 2014; 1 hour 20 minutes Answer key
NAME: psio 210 Introduction to Physical Oceanography Mid-term examination November 3, 2014; 1 hour 20 minutes Answer key Closed book; one sheet of your own notes is allowed. A calculator is allowed. (100
More informationAFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES
AFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES Victoria Dollar RTG Seminar Research - Spring 2018 April 16, 2018 Victoria Dollar ASU April 16, 2018 1 / 26 Overview Introduction Rossby waves and
More informationTracer Based Ages, Transit Time Distributions, and Water Mass Composition: Observational and Computational Examples
Tracer Based Ages, Transit Time Distributions, and Water Mass Composition: Observational and Computational Examples Frank Bryan Climate and Global Dynamics Division National Center for Atmospheric Research
More informationChapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationNCAR Global Atmospheric Core Workshop, Boulder, June 2008
NCAR Global Atmospheric Core Workshop, Boulder, June 2008 D. Majewski based on Christiane Jablonowski (email: cjablono@umich.edu) University of Michigan Goals of the Test Suite NASA/GFDL Test cases should
More informationThe General Circulation of the Atmosphere: A Numerical Experiment
The General Circulation of the Atmosphere: A Numerical Experiment Norman A. Phillips (1956) Presentation by Lukas Strebel and Fabian Thüring Goal of the Model Numerically predict the mean state of the
More informationInternal Wave Modeling in Oceanic General Circulation Models
Internal Wave Modeling in Oceanic General Circulation Models Flavien Gouillon, Eric P. Chassignet Center for Ocean - Atmospheric Prediction Studies The Florida State University Introduction One of the
More informationDiagnosis of a Quasi-Geostrophic 2-Layer Model Aaron Adams, David Zermeño, Eunsil Jung, Hosmay Lopez, Ronald Gordon, Ting-Chi Wu
Diagnosis of a Quasi-Geostrophic 2-Layer Model Aaron Adams, David Zermeño, Eunsil Jung, Hosmay Lopez, Ronald Gordon, Ting-Chi Wu Introduction For this project we use a simple two layer model, which is
More information( ) (9.1.1) Chapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation. 9.1 Geostrophy and scaling.
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationPhysical Oceanography, MSCI 3001 Oceanographic Processes, MSCI Dr. Katrin Meissner Week 5.
Physical Oceanography, MSCI 3001 Oceanographic Processes, MSCI 5004 Dr. Katrin Meissner k.meissner@unsw.e.au Week 5 Ocean Dynamics Transport of Volume, Heat & Salt Flux: Amount of heat, salt or volume
More informationLecture 8. Lecture 1. Wind-driven gyres. Ekman transport and Ekman pumping in a typical ocean basin. VEk
Lecture 8 Lecture 1 Wind-driven gyres Ekman transport and Ekman pumping in a typical ocean basin. VEk wek > 0 VEk wek < 0 VEk 1 8.1 Vorticity and circulation The vorticity of a parcel is a measure of its
More information8/21/08. Modeling the General Circulation of the Atmosphere. Topic 4: Equatorial Wave Dynamics. Moisture and Equatorial Waves
Modeling the General Circulation of the Atmosphere. Topic 4: Equatorial Wave Dynamics D A R G A N M. W. F R I E R S O N U N I V E R S I T Y O F W A S H I N G T O N, D E P A R T M E N T O F A T M O S P
More informationNonlinear baroclinic dynamics of surface cyclones crossing a zonal jet
Nonlinear baroclinic dynamics of surface cyclones crossing a zonal jet Jean-Baptiste GILET, Matthieu Plu and Gwendal Rivière CNRM/GAME (Météo-France, CNRS) 3rd THORPEX International Science Symposium Monterey,
More informationSome remarks on climate modeling
Some remarks on climate modeling A. Gettelman & J. J. Hack National Center for Atmospheric Research Boulder, Colorado USA Selected overheads by Doug Nychka Outline Hierarchy of atmospheric modeling strategies
More informationDownscaling Global Warming with a Regional Ocean- Atmosphere Model over the Tropical Atlantic
Downscaling Global Warming with a Regional Ocean- Atmosphere Model over the Tropical Atlantic Role of equatorial ocean dynamics: equatorial upwelling and ocean mesoscale variability Hyodae Seo and Shang-Ping
More information9 Rossby Waves. 9.1 Non-divergent barotropic vorticity equation. CSU ATS601 Fall (Holton Chapter 7, Vallis Chapter 5)
9 Rossby Waves (Holton Chapter 7, Vallis Chapter 5) 9.1 Non-divergent barotropic vorticity equation We are now at a point that we can discuss our first fundamental application of the equations of motion:
More informationDonald Slinn, Murray D. Levine
2 Donald Slinn, Murray D. Levine 2 Department of Civil and Coastal Engineering, University of Florida, Gainesville, Florida College of Oceanic and Atmospheric Sciences, Oregon State University, Corvallis,
More informationTransformed Eulerian Mean
Chapter 15 Transformed Eulerian Mean In the last few lectures we introduced some fundamental ideas on 1) the properties of turbulent flows in rotating stratified environments, like the ocean and the atmosphere,
More informationATS 421/521. Climate Modeling. Spring 2015
ATS 421/521 Climate Modeling Spring 2015 Lecture 9 Hadley Circulation (Held and Hou, 1980) General Circulation Models (tetbook chapter 3.2.3; course notes chapter 5.3) The Primitive Equations (tetbook
More informationDynamics and Kinematics
Geophysics Fluid Dynamics () Syllabus Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3,
More informationThe Planetary Circulation System
12 The Planetary Circulation System Learning Goals After studying this chapter, students should be able to: 1. describe and account for the global patterns of pressure, wind patterns and ocean currents
More informationAtmosphere, Ocean, Climate Dynamics: the Ocean Circulation EESS 146B/246B
Atmosphere, Ocean, Climate Dynamics: the Ocean Circulation EESS 146B/246B Instructor: Leif Thomas TA: Gonçalo Zo Zo Gil http://pangea.stanford.edu/courses/eess146bweb/ Course Objectives Identify and characterize
More informationESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves
ESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves Reference: An Introduction to Dynamic Meteorology (4 rd edition), J.R. Holton Atmosphere-Ocean Dynamics, A.E. Gill Fundamentals of Atmospheric
More informationDynamics of wind-forced coherent anticyclones in the open ocean
Dynamics of wind-forced coherent anticyclones in the open ocean Inga Koszalka, ISAC-CNR, Torino, Italy Annalisa Bracco, EAS-CNS, Georgia Tech Antonello Provenzale, ISAC-CNR, Torino, Italy Idealized set-up
More informationChapter 2. Quasi-Geostrophic Theory: Formulation (review) ε =U f o L <<1, β = 2Ω cosθ o R. 2.1 Introduction
Chapter 2. Quasi-Geostrophic Theory: Formulation (review) 2.1 Introduction For most of the course we will be concerned with instabilities that an be analyzed by the quasi-geostrophic equations. These are
More informationWhy build a climate model
Climate Modeling Why build a climate model Atmosphere H2O vapor and Clouds Absorbing gases CO2 Aerosol Land/Biota Surface vegetation Ice Sea ice Ice sheets (glaciers) Ocean Box Model (0 D) E IN = E OUT
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationA 1/10th Degree Global Ocean Simulation Using the Parallel Ocean Program
A 1/10th Degree Global Ocean Simulation Using the Parallel Ocean Program Mathew E Maltrud Fluid Dynamics Group MS B216 Los Alamos National Laboratory Los Alamos, NM 87545 phone: (505) 667-9097 fax: (505)
More informationLecture 25: Ocean circulation: inferences from geostrophic and thermal wind balance
Lecture 25: Ocean circulation: inferences from geostrophic and thermal wind balance November 5, 2003 Today we are going to study vertical sections through the ocean and discuss what we can learn about
More informationSusan Bates Ocean Model Working Group Science Liaison
Susan Bates Ocean Model Working Group Science Liaison Climate Simulation Laboratory (CSL) Accelerated Scientific Discovery (ASD) NCAR Strategic Capability (NSC) Climate Process Teams (CPTs) NSF Earth System
More informationLETTERS. Influence of the Thermohaline Circulation on Projected Sea Level Rise
VOLUME 13 JOURNAL OF CLIMATE 15 JUNE 2000 LETTERS Influence of the Thermohaline Circulation on Projected Sea Level Rise RETO KNUTTI AND THOMAS F. STOCKER Climate and Environmental Physics, Physics Institute,
More informationBALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity
BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere Need to introduce a new measure of the buoyancy Potential temperature θ In a compressible fluid, the relevant measure
More informationEquilibrium surface elevation distribution in a baroclinic ocean with realistic bottom topography
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 08, NO. C0, 330, doi:0.029/2002jc00579, 2003 Equilibrium surface elevation distribution in a baroclinic ocean with realistic bottom topography Srdjan Dobricic Institute
More informationLecture 1. Amplitude of the seasonal cycle in temperature
Lecture 6 Lecture 1 Ocean circulation Forcing and large-scale features Amplitude of the seasonal cycle in temperature 1 Atmosphere and ocean heat transport Trenberth and Caron (2001) False-colour satellite
More informationWind Gyres. curl[τ s τ b ]. (1) We choose the simple, linear bottom stress law derived by linear Ekman theory with constant κ v, viz.
Wind Gyres Here we derive the simplest (and oldest; Stommel, 1948) theory to explain western boundary currents like the Gulf Stream, and then discuss the relation of the theory to more realistic gyres.
More informationLFV in a mid-latitude coupled model:
LFV in a mid-latitude coupled model: Gulf Stream influence on the NAO M. Ghil 1,2, E. Simonnet 3, Y. Feliks 1,4 1 Dept. of Atmospheric and Oceanic Sciences and IGPP, UCLA, USA. 2 Geosciences Dept. and
More informationQuasi-geostrophic ocean models
Quasi-geostrophic ocean models March 19, 2002 1 Introduction The starting point for theoretical and numerical study of the three dimensional large-scale circulation of the atmosphere and ocean is a vorticity
More informationEART164: PLANETARY ATMOSPHERES
EART164: PLANETARY ATMOSPHERES Francis Nimmo Last Week Radiative Transfer Black body radiation, Planck function, Wien s law Absorption, emission, opacity, optical depth Intensity, flux Radiative diffusion,
More informationClimate of an Earth- like Aquaplanet: the high- obliquity case and the <dally- locked case
Climate of an Earth- like Aquaplanet: the high- obliquity case and the
More informationCourse , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 9: Numerical Models
Course 1.81, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 9: Numerical Models Phillips Numerical Experiment: We are going to discuss Phillips (1956) original numerical experiment
More informationMERIDIONAL OVERTURNING CIRCULATION: SOME BASICS AND ITS MULTI-DECADAL VARIABILITY
MERIDIONAL OVERTURNING CIRCULATION: SOME BASICS AND ITS MULTI-DECADAL VARIABILITY Gokhan Danabasoglu National Center for Atmospheric Research OUTLINE: - Describe thermohaline and meridional overturning
More informationConvection Induced by Cooling at One Side Wall in Two-Dimensional Non-Rotating Fluid Applicability to the Deep Pacific Circulation
Journal of Oceanography Vol. 52, pp. 617 to 632. 1996 Convection Induced by Cooling at One Side Wall in Two-Dimensional Non-Rotating Fluid Applicability to the Deep Pacific Circulation ICHIRO ISHIKAWA
More informationBarotropic geophysical flows and two-dimensional fluid flows: Conserved Quantities
Barotropic geophysical flows and two-dimensional fluid flows: Conserved Quantities Di Qi, and Andrew J. Majda Courant Institute of Mathematical Sciences Fall 2016 Advanced Topics in Applied Math Di Qi,
More informationOcean Modeling II. Parameterized Physics. Peter Gent. Oceanography Section NCAR
Ocean Modeling II Parameterized Physics Peter Gent Oceanography Section NCAR Ocean Modelling Challenges LO-RES (3 o ) O(100+ years/day) WORKHORSE (1 o ) O(10-100 years/day) HI-RES (0.1 o ) O(1 year/day)
More informationGeophysics Fluid Dynamics (ESS228)
Geophysics Fluid Dynamics (ESS228) Course Time Lectures: Tu, Th 09:30-10:50 Discussion: 3315 Croul Hall Text Book J. R. Holton, "An introduction to Dynamic Meteorology", Academic Press (Ch. 1, 2, 3, 4,
More informationEddy PV Fluxes in a One Dimensional Model of Quasi-Geostrophic Turbulence
Eddy PV Fluxes in a One Dimensional Model of Quasi-Geostrophic Turbulence Christos M.Mitas Introduction. Motivation Understanding eddy transport of heat and momentum is crucial to developing closure schemes
More informationPhysical Oceanography, MSCI 3001 Oceanographic Processes, MSCI Dr. Katrin Meissner Ocean Dynamics.
Physical Oceanography, MSCI 3001 Oceanographic Processes, MSCI 5004 Dr. Katrin Meissner k.meissner@unsw.e.au Ocean Dynamics The Equations of Motion d u dt = 1 ρ Σ F dt = 1 ρ ΣF x dt = 1 ρ ΣF y dw dt =
More informationDeveloping a nonhydrostatic isopycnalcoordinate
Developing a nonhydrostatic isopycnalcoordinate ocean model Oliver Fringer Associate Professor The Bob and Norma Street Environmental Fluid Mechanics Laboratory Dept. of Civil and Environmental Engineering
More informationFully-Coupled and Forced Ocean Sea-Ice Simulations Towards CESM2 (since mini-breck) Defining the Ocean Component of CESM2
Fully-Coupled and Forced Ocean Sea-Ice Simulations Towards CESM2 (since mini-breck) Defining the Ocean Component of CESM2 OMWG PLAN FOR DEFINING THE CESM2 OCEAN COMPONENT Timeline 01 March 2016 22 April
More informationModel equations for planetary and synoptic scale atmospheric motions associated with different background stratification
Model equations for planetary and synoptic scale atmospheric motions associated with different background stratification Stamen Dolaptchiev & Rupert Klein Potsdam Institute for Climate Impact Research
More informationIntroduction to Atmospheric Circulation
Introduction to Atmospheric Circulation Start rotating table Cloud Fraction Dice Results from http://eos.atmos.washington.edu/erbe/ from http://eos.atmos.washington.edu/erbe/ from http://eos.atmos.washington.edu/erbe/
More informationRoughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman
Roughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman 1. Characteristics of the Roughness Sub layer With well understood caveats, the time averaged statistics of flow in the atmospheric
More informationINTERNAL GRAVITY WAVES
INTERNAL GRAVITY WAVES B. R. Sutherland Departments of Physics and of Earth&Atmospheric Sciences University of Alberta Contents Preface List of Tables vii xi 1 Stratified Fluids and Waves 1 1.1 Introduction
More information2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven Ocean Fronts
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. 2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven
More information