Ocean Modeling Josephansong
Ocean Modeling Josephansong
MODELING
Joseph K. Ansong
COESSING-2018
OUTLINE
• INTRODUCTION
• MOTIVATION
• EQUATIONS OF MOTION
• EXAMPLES
• CHALLENGES
Introduction: Approaches to oceanic studies
Observational oceanography
Emily/Stephan et al.
Introduction: Approaches to oceanic studies
Satellite oceanography
Ebenezer et al.
Introduction: Approaches to oceanic studies
Laboratory tank
Experiments-
Emily/Stephan
Ocean modeling
Riccardo/Brian/Dimitris/Joseph
What is an ocean model?
It is a representation, in the form of
equations/computer code,
describing physical processes of
our understanding of how the ocean
works.
-Dr. Stephenie Waterman
What is an ocean model?
Physical processes:
a) Ocean movement/dynamics, including
horizontal and vertical advection
www.hycom.org
Motivation: Why model the ocean?
Understanding the 3D dynamics of the ocean on a
GLOBAL scale.
?
Prof. Modeller
Motivation: global ocean
currents
Internal gravity waves
¶u ¶v ¶w
+ + =0
¶x ¶y ¶z
Equation for tracers (Temperature, Salinity, and others),
Advection-diffusion equation:
¶T r
c
+ u ×ÑT = kT c Ñ T
c 2 c
¶t
T c = Tracers
Equation of state (Linear)
r = r0 [1 - a (T - T0 ) + b (S- S0 )]
r0 = 1028 kg / m 3
coefficients of thermal , a ,
T0 = 10 C = 283K
o
and saline contraction, b
S0 = 35 psu
D. Boundary conditions
Discretize equations
• Continuous equations
algebraic equations
(discrete set of operations)
• Discretization methods:
– Finite difference methods
– Finite element methods
– Finite volume methods Example early model grid by
Lewis Fry Richardson (1928)
Model grid: horizontal
• Regular grids: regularly
spaced lines
GHANA
Regional grid
Clever solution: tripolar grid
-circular grid laid over Arctic region with poles on land
Model grid: vertical
• z-coordinate system
based on a series of
depth levels. Easy to
setup. Difficult to locally
increase resolution.
• terrain-following
coordinate system.
Mimics bathymetry and
allows higher resolution
near ocean floor.
Model grid: vertical
• density (isopycnal)-
coordinate system
based on density
layers. Great in the
deep ocean where
there’s less diapycnal
mixing. Poor in
regions with high
vertical mixing.
Model grid: vertical
• hybrid-coordinate
applies the best
suited coordinate
system in different
regions. Gives
improved results but
at a high
computational cost.
Boundary conditions
• Free surface
– Flux exchanges at surface: momentum and tracer
(winds, solar radiation, rainfall, precipitation, etc).
• Ocean bottom
– Topography/bathymetry
– Velocity normal to bottom is zero
– Lateral boundaries (open/closed)
Boundary conditions
• Step-by-step method of
Boundary conditions
solution
(model time stepping) EQUATIONS
Etopo5: http://www.ngdc.noaa.gov/mgg/global/etopo5.HTML
Results: Tides
Results: tidal time series