3-Physics Dudhia
3-Physics Dudhia
3-Physics Dudhia
Jimy
Dudhia
NCAR
WRF
Physics
• Radia;on
– Longwave
(ra_lw_physics)
– Shortwave
(ra_sw_physics)
• Surface
– Surface
layer
(sf_sfclay_physics)
– Land/water
surface
(sf_surface_physics)
• PBL
(bl_pbl_physics)
• Cumulus
parameteriza;on
(cu_physics)
• Microphysics
(mp_physics)
• Turbulence/Diffusion
(diff_opt,
km_opt)
Radia;on
Provides
Atmospheric
temperature
tendency
profile
Surface
radia;ve
fluxes
Radia;on
parameteriza;ons
• Shortwave
(solar)
–
use
astronomical
equa;ons
to
calculate
the
sun’s
posi;on
as
a
func;on
of
;me
of
day
and
day
of
the
year.
Radia;on
as
part
of
the
en;re
model
energy
budget
WRF
Longwave
Radia;on
Schemes
(ra_lw_physics)
• Compute
clear-‐sky
and
cloud
upward
and
downward
radia;on
fluxes
– Consider
IR
emission
from
layers
– Surface
emissivity
based
on
land-‐type
– Flux
divergence
leads
to
cooling
in
a
layer
– Downward
flux
at
surface
important
in
land
energy
budget
– IR
radia;on
generally
leads
to
cooling
in
clear
air
(~2K/day),
stronger
cooling
at
cloud
tops
and
warming
at
cloud
base
Longwave
Radia;on
schemes
ra_lw_physics Scheme Reference Added
1 RRTM Mlawer et al. (1997, JGR) 2000
5 New Goddard Chou and Suarez (2001, NASA Tech Memo) 2011
31 Held-Suarez 2000
Land-surface fluxes
of heat and moisture
PBL
Surface
Fluxes
• Heat,
moisture
and
momentum
H = ρc p ChsΔθ
It is related to the roughness length
and u* by
ku*
Chs =
⎛ z ⎞
€ ln⎜ ⎟ − ψ h
⎝ z0 ⎠
WRF
Surface
Layer
Op;ons
(sf_sfclay_physics)
• Use
similarity
theory
to
determine
exchange
coefficients
and
diagnos;cs
of
2m
T
and
q
and
10
m
winds
• Provide
exchange
coefficient
to
land-‐surface
models
• Provide
fric;on
velocity
to
PBL
scheme
• Provide
surface
fluxes
over
water
points
• Schemes
have
varia;ons
in
stability
func;ons,
roughness
lengths
ARW only
Hurricane
Op;ons
• Ocean
Mixed
Layer
Model
(omlcall=1)
– 1-‐d
slab
ocean
mixed
layer
(specified
ini;al
depth)
– Includes
wind-‐driven
ocean
mixing
for
SST
cooling
feedback
• Alterna;ve
surface-‐layer
op;ons
for
high-‐wind
ocean
surface
(isncflx=1,2)
– Use
with
sf_sfclay_physics=1
– Modifies
Charnock
rela;on
to
give
less
surface
fric;on
at
high
winds
(lower
Cd)
– Modifies
surface
enthalpy
(Ck,
heat/moisture)
either
with
constant
z0q
(isncflx=1),
Garraf
formula;on
(op;on
2)
Frac;onal
Sea
Ice
• frac;onal_seaice=1
-‐
with
input
sea-‐ice
frac;on
data
can
par;;on
land/water
fluxes
within
a
grid
box
WRF
Land-‐Surface
Model
Op;ons
(sf_surface_physics)
• Simple
5-‐layer
soil
model
– No
vegeta;on
or
snow
cover
predic;on,
just
thermal
diffusion
in
soil
layers
• Noah
LSM,
RUC
LSM,
PX
LSM
land-‐surface
models
– Sophis;cated
vegeta;on
model
and
snow
cover
predic;on
Land-‐Surface
Model
Processes
Land-‐Surface
Model
• Predicts
soil
temperature
and
soil
moisture
in
layers
(4
for
Noah,
6
for
RUC,
2
for
PX)
• Predicts
snow
water
equivalent
on
ground.
May
be
in
layers
(RUC)
• May
predict
canopy
moisture
(Noah)
Land-‐Surface
Op;ons
• 5-‐layer
thermal
diffusion
• Noah
LSM
• RUC
LSM
• Pleim-‐Xiu
LSM
Vegeta;on
and
Soil
• Processes
include
evapotranspira;on,
root
zone
and
leaf
effects
• Vegeta;on
frac;on
varies
seasonally
• Considers
vegeta;on
categories
(e.g.
cropland,
forest
types,
etc.)
• Considers
soil
categories
(e.g.
sandy,
clay,
etc.)
for
drainage
and
thermal
conduc;vity
Urban
Effects
• Urban
category
in
LSM
is
usually
adequate
for
larger-‐scale
studies
• Or
can
use
an
urban
model
(sf_urban_physics)
with
Noah
LSM
– Urban
Canopy
Model
– Building
Environment
Parameteriza;on
(mul;-‐
layer
model)
– Building
Energy
Model
(adds
hea;ng/AC
effects
to
BEP)
LSM
Tables
• Proper;es
can
be
changed
in
text
files
(tables)
• VEGPARM.TBL
used
by
Noah
and
RUC
for
vegeta;on
category
proper;es
• SOILPARM.TBL
used
by
Noah
and
RUC
for
soil
proper;es
• LANDUSE.TBL
used
by
5-‐layer
model
• URBPARM.TBL
used
by
urban
models
Ini;alizing
LSMs
• Noah
and
RUC
LSM
require
addi;onal
fields
for
ini;aliza;on
• Soil
temperature
• Soil
moisture
• Snow
liquid
equivalent
• These
are
in
the
Grib
files,
but
are
not
from
observa;ons
• They
come
from
“offline”
models
driven
by
observa;ons
(rainfall,
radia;on,
surface
temperature,
humidity
wind)
Ini;alizing
LSMs
• There
are
consistent
model-‐derived
datasets
for
Noah
and
RUC
LSMs
• Eta/GFS/AGRMET/NNRP
for
Noah
(although
some
have
limited
soil
levels
available)
• RUC
for
RUC
• But,
resolu;on
of
mesoscale
land-‐use
means
there
will
be
inconsistency
in
eleva;on,
soil
type
and
vegeta;on
• This
leads
to
spin-‐up
as
adjustments
occur
in
soil
temperature
and
moisture
• This
spin-‐up
can
only
be
avoided
by
running
offline
model
on
the
same
grid
(e.g.
HRLDAS
for
Noah)
• Cycling
land
state
between
forecasts
also
helps,
but
may
propagate
errors
(e.g
in
rainfall
effect
on
soil
moisture)
ARW only
sst_update=1
Reads
lower
boundary
file
periodically
to
update
the
sea-‐surface
temperature
(otherwise
it
is
fixed
with
;me)
• For
long-‐period
simula;ons
(a
week
or
more)
• wrflowinp_d0n
created
by
real
• Sea-‐ice
can
be
updated
since
Version
3.0
• Vegeta;on
frac;on
update
is
included
– Allows
seasonal
change
in
albedo,
emissivity,
roughness
length
in
Noah
LSM
• usemonalb=.true.
to
use
monthly
albedo
input
Regional
Climate
Op;ons
Provides
Boundary
layer
fluxes
(heat,
moisture,
momentum)
Ver;cal
diffusion
in
whole
column
Planetary
Boundary
Layer
WRF
PBL
Op;ons
(bl_pbl_physics)
• Purpose
is
to
distribute
surface
fluxes
with
boundary
layer
eddy
fluxes
and
allow
for
PBL
growth
by
entrainment
• Classes
of
PBL
scheme
– Turbulent
kine;c
energy
predic;on
(Mellor-‐
Yamada
Janjic,
MYNN,
Bougeault-‐Lacarrere,
TEMF,
QNSE,
CAM
UW)
– Diagnos;c
non-‐local
(YSU,
GFS,
MRF,
ACM2)
• Above
PBL
all
these
schemes
also
do
ver;cal
diffusion
due
to
turbulence
PBL
schemes
in
V3.3
bl_pbl_ Scheme Reference Added
physics
1 YSU Hong, Noh and Dudhia (2006, MWR) 2004
mixing
diff_opt=1
Horizontal
diffusion
acts
along
model
levels
Simpler
numerical
method
with
only
neighboring
points
on
the
same
model
level
Difference
between
diff_opt
1
and
2
diff_opt=2
Horizontal diffusion acts on horizontal gradients
Numerical method includes vertical correction term
using more grid points
ARW only
km_opt
• km_opt
selects
method
for
compu;ng
K
coefficient
– km_opt=1:
constant
(use
khdif
and
kvdif
to
specify
–
idealized)
– km_opt=2:
3d
tke
predic;on
used
to
compute
K
(requires
diff_opt=2)
– km_opt=3:
3d
Smagorisnky
diagnos;c
K
(requires
diff_opt=2)
– km_opt=4:
2d
Smagorinsky
for
horizontal
K
(to
be
used
with
PBL
or
kvdif
for
ver;cal
K
)
ARW only
sfs_opt
• Sub-‐filter-‐scale
stress
model
for
LES
applica;ons
impac;ng
momentum
mixing
(Kosovic,
Mirocha)
– sfs_opt=0
(default)
off
– sfs_opt=1
Nonlinear
Backscafer
and
Anisotropy
(NBA)
op;on
1:
using
diagnos;c
stress
terms
(km_opt=2,3)
– sfs_opt=2
NBA
op;on
2:
using
tke-‐based
stress
terms
(km_opt=2
only)
– Also
m_opt=1
for
added
outputs
of
SGS
stresses
Diffusion
Op;on
Choice
• Real-‐data
case
with
PBL
physics
on
– Best
is
diff_opt=1,
km_opt=4
– This
complements
ver;cal
diffusion
done
by
PBL
scheme
• High-‐resolu;on
real-‐data
cases
(~100
m
grid)
– No
PBL
– diff_opt=2;
km_opt=2,3
(tke
or
Smagorinsky
scheme)
• idealized
cloud-‐resolving
modeling
(smooth
or
no
topography)
– diff_opt=2;
km_opt=2,3
• Complex
topography
with
no
PBL
scheme
– diff_opt=2
is
more
accurate
for
sloped
coordinate
surfaces,
and
prevents
diffusion
up/down
valley
sides
• Note:
WRF
can
run
with
no
diffusion
(diff_opt=0)
ARW only
diff_6th_opt
• 6th
order
op;onal
added
horizontal
diffusion
on
model
levels
– Used
as
a
numerical
filter
for
2*dx
noise
– Suitable
for
idealized
and
real-‐data
cases
– Affects
all
advected
variables
including
scalars
• diff_6th_opt
– 0:
none
(default)
– 1:
on
(can
produce
nega;ve
water)
– 2:
on
and
prohibit
up-‐gradient
diffusion
(befer
for
water
conserva;on)
• diff_6th_factor
– Non-‐dimensional
strength
(typical
value
0.12,
1.0
corresponds
to
complete
removal
of
2*dx
wave
in
a
;me-‐step)
ARW only
Upper
damping
(damp_opt)
Purpose
is
to
prevent
unrealis;c
reflec;ons
of
waves
from
model
top.
Can
be
important
over
high
topography.
Op;ons
• 1:
Upper
level
diffusive
layer
• 2:
Rayleigh
damping
(idealized
only
–
needs
input
sounding)
• 3:
w-‐Rayleigh
damping
(damps
w
only)
All
op;ons
use
– Cosine
func;on
of
height
– Addi;onal
parameters
• zdamp:
depth
of
damping
layer
• dampcoef:
nondimensional
maximum
magnitude
of
damping
Cumulus
Parameteriza;on
Provides
Atmospheric
heat
and
moisture/cloud
tendency
profiles
Surface
sub-‐grid-‐scale
(convec;ve)
rainfall
Cumulus
Schemes
• Use
for
grid
columns
that
completely
contain
convec;ve
clouds
• Re-‐distribute
air
in
column
to
account
for
ver;cal
convec;ve
fluxes
• Updrans
take
boundary
layer
air
upwards
• Downdrans
take
mid-‐level
air
downwards
• Schemes
have
to
determine
• When
to
trigger
a
convec;ve
column
• How
fast
to
make
the
convec;on
act
Parameteriza;ons
of
cumulus
convec;on
WRF
Cumulus
Parameteriza;on
Op;ons
• Cumulus
schemes
fall
into
two
main
classes
– Adjustment
type
(Befs-‐Miller-‐Janjic)
• Relaxes
towards
a
post-‐convec;ve
(mixed)
sounding
– Mass-‐flux
type
(all
others
in
WRF)
• Determines
updran
(and
onen
downdran)
mass
flux
and
other
fluxes
(some;mes
including
momentum
transport)
Cumulus
schemes
in
V3.3
mp_physics Scheme Reference Added
1 Kain-Fritsch Kain (2004, JAM) 2000
99 Old Kain-Fritsch Kain and Fritsch (1990, JAS; 1993 Meteo. 2000
Monogr.)
Deep
Convec;on
• Schemes
work
in
individual
columns
that
are
considered
convec;vely
unstable
• Mass-‐flux
schemes
transport
surface
air
to
top
of
cloud
and
include
subsidence
• Subsidence
around
cloud
warms
and
dries
troposphere
removing
instability
over
;me
• Addi;onally
downdrans
may
cool
PBL
Triggers
• Clouds
only
ac;vate
in
columns
that
meet
certain
criteria
– Presence
of
some
convec;ve
available
poten;al
energy
(CAPE)
in
sounding
– Not
too
much
convec;ve
inhibi;on
(CIN)
in
sounding
(cap
strength)
– Minimum
cloud
depth
from
parcel
ascent
Closures
• Closure
determine
cloud
strength
(mass-‐flux)
based
on
various
methods
– Clouds
remove
CAPE
over
;me
• Specified
CAPE-‐removal
;me
scale
(KF,
Tiedtke,
ZM,
BMJ)
• Quasi-‐equilibrium
(Arakawa-‐Schubert)
with
large-‐scale
destabiliza;on
d(CAPE)/dt
(SAS,
NSAS)
• Column
moisture
convergence
• Low-‐level
large-‐scale
ascent
(mass
convergence)
Ensemble
methods
• G3
and
GD
use
ensemble
of
triggers
and
closures
with
varying
parameters
(effec;vely
144
members)
• Take
mean
of
ensemble
to
feed
back
to
model
• In
principle,
can
be
tuned
to
emphasize
various
members
under
different
condi;ons
Shallow
Convec;on
• Non-‐precipita;ng
shallow
mixing
dries
PBL,
moistens
and
cools
above
• This
can
be
done
by
an
enhanced
mixing
approach
(SAS)
or
mass-‐flux
approach
(KF,
NSAS,
Tiedtke,
G3)
Shallow
Convec;on
• Several
schemes
include
shallow
convec;on
(KF,
SAS
schemes,
G3,
BMJ,
Tiedtke)
– WRF
also
has
UW
Park-‐Bretherton
stand-‐alone
scheme
(shcu_physics=2)
– Note:
TEMF
PBL
op;on
(bl_bl_physics=10)
also
includes
a
mass-‐flux
shallow
convec;on
component
and
some
future
PBL
schemes
may
add
it
• This
development
will
require
making
shallow
convec;onal
independent
of
deep
schemes
in
a
future
WRF
version
Momentum
Transport
• Some
cumulus
parameteriza;ons
also
have
momentum
transport
(SAS,
NSAS,
Tiedtke,
ZM)
• Most
schemes
transport
momentum
as
a
passive
scalar
but
ZM
and
NSAS
include
a
convec;ve
pressure
gradient
term
Cloud
Detrainment
• Most
schemes
detrain
cloud
and
ice
at
cloud
top
(except
BMJ)
• KF
schemes
also
detrain
snow
and
rain
• These
are
then
used
by
the
microphysics
ARW only
cudt
• Time
steps
between
cumulus
scheme
calls
• Typical
value
is
5
minutes
– Note:
for
KF
scheme
this
is
also
used
for
averaging
;me
for
ver;cal
velocity
trigger
– Not
used
by
G3
or
GD
schemes
Cumulus
schemes
in
V3.3
cu_physics Scheme Cores Moisture Momentum Shallow
Tendencies Tendencies Convection
1 Kain-Fritsch Eta ARW NMM Qc Qr Qi Qs no yes
3 Grell-Devenyi ARW Qc Qi no no
Provides
Atmospheric
heat
and
moisture
tendencies
Microphysical
rates
Surface
resolved-‐scale
rainfall
Resolved
clouds
Kessler WSM3
Ferrier
Qv
Qi/Qs/
Qc
Qg
WSM5 Lin et al./WSM6
Qr
WRF
Microphysics
Op;ons
(mp_physics)
• Range
of
levels
of
sophis;ca;on
– Warm
rain
(i.e.
no
ice)
–
Kessler
(idealized)
– Simple
ice
(3
arrays)
–
WSM3
– Mesoscale
(5
arrays,
no
graupel)
–
WSM5
– Cloud-‐scale
single-‐moment
(6
arrays,
graupel)
–
WSM6,
Lin,
Goddard,
SBU,
Eta-‐Ferrier
– Double-‐moment
(8-‐13
arrays)
–
Thompson,
Morrison,
Milbrandt-‐Yau,
WDM5,
WDM6
Microphysics
schemes
in
V3.3
mp_physics Scheme Reference Added
1 Kessler Kessler (1969) 2000
3 WSM3 ARW Qc Qr
7 Goddard ARW Qc Qr Qi Qs Qg
13 SBU-YLin ARW Qc Qr Qi Qs