merging the clubb and edmf frameworks in cam6: …
TRANSCRIPT
Merging CLUBB and EDMF in CAM6: Approach, SCM Results and Next
StepsMikael Witte1,2, Marcin Kurowski2, Kay Suselj2, Julio Bacmeister3,
George Matheou4 and Joao Teixeira1,2
1 JIFRESSE, UCLA2 Jet Propulsion Laboratory, Caltech
3 NCAR CGD4 UConn
Why CLUBB and MF?
• ED component should account for subgrid scale (SGS) variability in the environment, e.g. by integrating over an assumed PDF
• Compared with our CAM5 EDMF using Holtslag-Boville PBL scheme for ED (Kurowski et al., 2019), coupling to CLUBB is a more statistically consistent approach
• MF is designed to sample the right tail of SGS distributions of ql, qt(although not a double Gaussian by default)
CAM6 implementation
• MF plumes affect CAM qt, ql via CLUBB’s 5-diagonal prognostic solver for mean fields and turbulent fluxes (solved simultaneously):
• Plume macrophysics (i.e. plume qc and moist area) are diagnosed but not yet coupled to other processes (e.g. radiation, aerosol)• Plume microphysics will be implemented next (stay tuned!)
!𝜑!"#!
Δ𝑡 +1𝜌$
𝜕𝜕𝑧 𝜌$𝑤
%𝜑′&'())!"#! =!𝜑!
Δ𝑡 −1𝜌$
𝜕𝜕𝑧 𝜌$.𝑎* 𝑤*𝜑*%
+,
!
+ 0𝜕 !𝜑𝜕𝑡 -./0*12
SCM results: BOMEX
• Compare SCAM vs. LES, both averaged over hours 4-6• Deep scheme and large-scale precip (MG2) off• Radiation on*• Plume qc and moist area diagnosed & added to CLUBB rcm, cloudfrac
(and therefore to CLDLIQ and CLOUD)• LES results from George Matheou (UConn) with standard GCSS forcing
300 305 310 315 320THETAL [K]
0
500
1000
1500
2000
2500
3000alt
itude
[m]
0 0.005 0.01 0.015 0.02QT [kg/kg]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]0 0.02 0.04 0.06
CLOUD [fraction]0
500
1000
1500
2000
2500
3000
altitu
de [m
]
0 0.5 1CLDLIQ [kg/kg] 10-5
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
-5 0 5 10DCQ [kg/kg/s] 10-8
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
-1 -0.5 0 0.5 1DTCOND [K/s] 10-4
0
500
1000
1500
2000
2500
3000alt
itude
[m]
-20 -10 0 10 20WPTHLP_CLUBB [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]0 50 100 150 200
WPRTP_CLUBB [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
LESnlev=32nlev=60nlev=120
LESnlev=32nlev=60nlev=120
Higher resolution => better agreement with LES (no surprise)
300 305 310 315THETAL [K]
0
500
1000
1500
2000
2500
3000alt
itude
[m]
0 0.005 0.01 0.015 0.02QT [kg/kg]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]0 0.02 0.04 0.06
CLOUD [fraction]0
500
1000
1500
2000
2500
3000
altitu
de [m
]
0 0.5 1CLDLIQ [kg/kg] 10-5
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
0 5 10 15 20DCQ [kg/kg/s] 10-8
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
-1 0 1DTCOND [K/s] 10-4
0
500
1000
1500
2000
2500
3000alt
itude
[m]
-20 -10 0 10 20THETAL flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]0 50 100 150 200
QT flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
LESCLUBB, totalCLUBB+MF, totalCLUBB+MF, CLUBBCLUBB+MF, MF
295 300 305 310 315THETAL [K]
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
0 0.005 0.01 0.015 0.02QT [kg/kg]
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
0 0.02 0.04 0.06CLOUD [fraction]
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
0 0.5 1CLDLIQ [kg/kg] 10-5
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
-40 -20 0 20THETAL flux [W m-2]
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
0 50 100 150 200QT flux [W m-2]
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
LESCLUBBCLUBB+MFCLUBUCLUBU+MF
0 0.5 1 1.5 2MF area-weighted <w> [m/s]
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
0 0.02 0.04 0.06 0.08MF total plume area [fraction]
0
500
1000
1500
2000
2500
3000
altit
ude
[m]
SCM results: ARM shallow Cu
• ARM SGP diurnal cycle of shallow cumulus (Brown et al., 2002)• Non-precipitating• No longer standard SCAM case
• Similar configuration to BOMEX, but uses time-varying IOP surface fluxes and ground temperature
• Transient case vs. BOMEX quasi-steady state
300 305 310 315 320THETAL [K]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
1430 UTC
300 305 310 315 320THETAL [K]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
1730 UTC
300 305 310 315 320THETAL [K]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2030 UTC
300 305 310 315 320THETAL [K]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2330 UTC
0 0.005 0.01 0.015 0.02QT [kg/kg]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
1430 UTC
0 0.005 0.01 0.015 0.02QT [kg/kg]
0
500
1000
1500
2000
2500
3000alt
itude
[m]
1730 UTC
0 0.005 0.01 0.015 0.02QT [kg/kg]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2030 UTC
0 0.005 0.01 0.015 0.02QT [kg/kg]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2330 UTCLESSCAM CLUBBSCAM CLUBB+MF
• ARM
-50 0 50 100THETAL flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
1430 UTC
-50 0 50 100THETAL flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
1730 UTC
-50 0 50 100THETAL flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2030 UTC
-50 0 50 100THETAL flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2330 UTC
0 100 200 300 400QT flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
1430 UTC
0 100 200 300 400QT flux [W m-2]
0
500
1000
1500
2000
2500
3000alt
itude
[m]
1730 UTC
0 100 200 300 400QT flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2030 UTC
0 100 200 300 400QT flux [W m-2]
0
500
1000
1500
2000
2500
3000
altitu
de [m
]
2330 UTCLESCLUBB, totalCLUBB+MF, totalCLUBB+MF, CLUBBCLUBB+MF, MF
2 4 6 8 10 12 140
500
1000
1500
2000
2500
3000al
titud
e [m
]
0
0.2
0.4
0.6
0.8
edm
f_qt
flx_f
rac
[-]
time [hr]
Summary
• Successful initial implementation of CLUBB+MF in SCAM6• MF coupled to CLUBB via prognostic solver for mean fields and turbulent fluxes• Turbulent fluxes increase compared to CLUBB alone• Cloud water mixing ratio and fraction also increase
• MF plumes compensate when CLUBB is reduced to a single Gaussian• ARM diurnal cycle case demonstrates more efficient mixing with MF
• Next steps: plumbing MF to other CAM schemes, adding warm precipto plumes, begin 3D testing
Thanks!