1 of 27

LES and Climate Modeling with COMBLE

GISS ModelE3-devel

GISS ModelE2.1

source: Andrew Ackerman / GISS

Ann Fridlind • NASA GISS

with contributions from many ...

subtropical �Sc decks

2 of 27

Background: Clouds in Earth system models (ESMs)

  • across ESMs, differences in cloud-related physics contribute substantially to large differences in simulated climate change
    • during ESM development, cloud physics parameters are commonly slightly tuned to bring simulations into line with targets such as preindustrial TOA radiative balance (e.g., Schmidt et al. GMD 2017)
    • as Earth warms in a 2X CO2 experiment, ESMs predict differing cloud changes, which feed back on warming (cloud-climate feedbacks)
    • e.g., CMIP6 ESM-predicted range of 1.8–5.6 K surface warming was found to be inconsistent with other estimates of 1.5–4.5 K, "tied to the physical representation of clouds" (Zelinka et al. GRL 2020)
    • clouds participate in a complex coupled Earth system environment

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

3 of 27

ModelE3 development approach

Field campaigns 🡪 LES 🡪 SCM

CALIPSO

Global data 🡪 ESM tuning

GMAO/cubed-sphere

ACTIVATE Flight RF13

1 March 2020

mixed-phase cold-air outbreak

Elsaesser et al., in prep.

Tornow et al. (ACP 2021)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

4 of 27

NASA GISS large-eddy simulation (LES) model

  • Distributed Hydrodynamic Aerosol-Radiation Model for Atmospheres (DHARMA)
    • second-order forward-in-time process-splitting with advection in semi-Lagrangian spirit (Stevens and Bretherton, 1997)
    • dynamic Smagorinsky subgrid-scale mixing (Kirkpatrick et al., 2006)
    • two-stream radiative transfer (Toon et al., 1989)
    • microphysics options
      • condensational adjustment (liquid only)
      • two-moment mixed-phase cloud microphysics after Morrison et al., 2009 (up to five species) with prognostic supersaturation (after Morrison and Grabowski, 2008) and optional substitution of alternative schemes for individual processes
      • fixed droplet number or multi-modal one-moment prognostic aerosol
      • size-resolved microphysics (Ackerman et al., 1995; Fridlind et al., 2004) with flexible treatment of multiple non-spherical ice classes (Fridlind et al., 2017)
    • can be run as a column model (Fridlind et al. 2015) or a parcel model (Ackerman et al., 2015)
    • used in scientific studies and as an engineering tool for ESM development

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

5 of 27

Field campaigns —> LES —> Single-column model (SCM)

5

Conditions

Case study

Aerosol aware?

dry convective boundary layer

idealized [Bretherton and Park 2009]

dry stable boundary layer

GABLS1 [Cuxart et al. 2006]

marine stratocumulus

DYCOMS-II RF02 [Ackerman et al. 2009]

observed (2 modes)

marine trade cumulus (shallow)

BOMEX [Siebesma et al. 2003]

marine trade cumulus (deep, raining)

RICO [van Zanten et al. 2011]

marine stratocumulus-to-cumulus *

SCT [Sandu and Stevens 2011]

continental cumulus ^

RACORO [Vogelmann et al. 2015]

observed profile (3 modes)

Arctic mixed-phase stratus

M-PACE [Klein et al. 2009]

observed (2 modes)

Antarctic mixed-phase stratus *

AWARE [Silber et al. 2019, 2021, 2022]

estimated (1 mode)

tropical deep convection

TWP-ICE [Fridlind et al. 2012]

observed profile (3 modes)

mid-latitude synoptic cirrus *

SPARTICUS [cf. Mühlbauer et al. 2014]

mid-latitude cold-air outbreak *^

ACTIVATE [Tornow et al., 2021, 2022, in prep.]

observed profile (3 modes)

high-latitude cold-air outbreak *^

COMBLE [Tornow et al., in prep.]

observed/estimated profiles (3 modes, 1 INP)

marine cumulus and congestus *^

CAMP2Ex [Stanford et al., in prep.]

observed profiles (3 modes)

subtropical marine deep convection *^

SEAC4RS [Stanford et al., in prep.]

observed profiles (TBD)

continental sea breeze convection *^

TRACER [Matsui et al., in prep.]

observed profiles (TBD)

*Lagrangian (cf. Neggers JAMES 2015, Pithan et al. NatGeo 2019)

^ensemble (cf. Neggers et al. JAMES 2019)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

6 of 27

Moving into the Lagrangian reference frame

  • Lagrangian and aerosol-aware go hand-in-hand

• SHEBA [Morrison et al. 2011]

Vertically pointing mm-wavelength radar

Fridlind and Ackerman [Ch. 7 in Mixed-Phase Clouds, Ed. C. Andronache, 2018]

• ISDAC [Ovchinnikov et al. 2014]

• M-PACE [Klein et al. 2009]

M-PACE [Klein et al. 2009]

Eulerian �(stationary)

Lagrangian �(moving with cloudy air)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

7 of 27

Background: Ice formation in supercooled clouds

  • esp. large uncertainties remain in supercooled cloud physics
    • ice spontaneously melts at 0°C but does not spontaneously form until near –38°C owing to energetic barriers to forming a more organized state
    • within that temperature range, supercooled liquid appears to be a gateway to persistent weak ice formation, attributable to ice-nucleating particle activation (Silber et al. 2020 based on NSA and AWARE data)
    • however, ice crystals are commonly orders of magnitude more abundant than can be explained by that weak pathway, especially in clouds with riming and/or large drops (Rangno and Hobbs, 2002; Korolev et al., 2021)
    • exactly how much more abundant is often unknown owing to lack of reliable measurements (Korolev et al. 2021; Morrison et al. 2020)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

8 of 27

Ice formation approach in ModelE3

  • Only physically-based mechanisms and parameterizations
  • Avoid unnecessary complexity
  • Each mechanism should be demonstrably active in observed case studies
  • Heterogeneous freezing mechanisms should be linkable to aerosol properties
  • But start with diagnostic INP and tune with machine learning
    • DeMott et al. 2010 * fscale_iifn
    • fscale_iifn < 1 can crudely account for efficient precipitation scavenging (Fridlind et al. JAS 2012)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

9 of 27

Mechanism

Include?

Comments

Primary

homogeneous freezing

Y

aerosol, cloud droplets, rain

immersion freezing

Y

aerosol, cloud droplets

deposition freezing

Y

aerosol

contact freezing

N

lab/field support currently lacking

Secondary

rime-splintering

Y

lab/field support poorly constrained

drop fragmentation*

N

lab/field support currently lacking

ice-ice collisions*

N

lab/field support currently lacking

Other common elements

Bigg [PPSB 1953]

N

no link to aerosol properties

Bergeron enhancement

N

ice vapor growth already included

* additive to Gettelman and Morrison (2015)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

10 of 27

Highly supercooled drizzle over Antarctica

AWARE campaign case study (Silber et al. JGR 2019)

CTT ≈ –25°C

Lubin et al.�(BAMS 2019)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

11 of 27

AWARE case study

  • SCM performs quite well
  • stable conditions common�(Silber et al. ACP 2020; GRL 2021)

see Silber et al. [GMD 2022]

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

12 of 27

ACTIVATE LES case study selection

  • choose 2020-2022 flights with greatest fetch offshore

Tornow et al. �(in prep.)

MBL �aerosol

FT �aerosol

DHARMA LES

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

13 of 27

COMBLE LES case study selection

  • choose CAO back trajectories passing near Zeppelin

source: Christian Lackner

source: �Paul DeMott

Tornow et al. (in prep.; see AMS 12B.5)

Ranking

Intensity

Data Availability

1st

Mar 12-13

May 7

2nd

Mar 28-29

Apr 25

3rd

Feb 2-6

May 11-12

4th

Jan 4

Apr 9-10

5th

Dec 1-2

Mar 12-13

6th

Apr 9-10

Feb 2-6

7th

Feb 23-24

Dec 1-2

8th

Dec 31

Feb 23-24

9th

Jan 21-22

Dec 31

10th

Dec 9

Mar 28-29

11th

May 11-12

Dec 9

12th

May 7

Jan 4

13th

Apr 25

Jan 21-22

INP at Andenes

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

14 of 27

COMBLE observational constraint

Silber et al. �(in prep.)

  • use EMC2 (Silber et al. GMD 2022) to evaluate LES vs ground-based radar + lidar

LWP confined to cell cores

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

15 of 27

COMBLE observational constraint

Silber et al. �(in prep.)

  • use EMC2 (Silber et al. GMD 2022) to evaluate LES vs CALIPSO satellite
  • LES clouds too deep + dense

DHARMA LES

CALIPSO

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

16 of 27

Now onto machine learning tuning...

  • LES modeler: We don't even understand how most atmospheric ice crystals are formed!
  • Climate modeler: We need to deliver a climate model now...

GISS ModelE3-devel

GISS ModelE2.1

source: Andrew Ackerman / GISS

subtropical �Sc decks

17 of 27

Tuning Protocol

17

  • scale_iifn is one of 45 parameters taken to be poorly constrained
  • LES/SCM used to estimate parameter ranges
  • new satellite datasets used in tuning
  • satellite dataset uncertainties are specified

Metrics (36 in total)

Data Source

Radiation (Longwave [LW], Shortwave [SW])

CERES-EBAF-Ed4.1

Cloud Radiative Forcing (LWcrf, SWcrf)

CERES-EBAF-Ed4.1

Column Water Vapor (CWV)

*Obs4MIPS RSS, G-VAP

Specific Humidity profiles (qv)

*Obs4MIPS AIRS, MLS

Temperature profiles (T)

*Obs4MIPS AIRS, MLS, GNSS-RO

Total Liquid Water Path (TLWP)

*MAC-LWP, GPM/TRMM

Total Ice Water Path (TIWP)

*CloudSat, MODIS

Total Precipitation (Pr)

*GPCP, GPM/TRMM

Convective Precipitation (Prc)

GPM/TRMM

Total Cloud Cover (TCC)

CloudSat/CALIPSO, ISCCP

Low (Shallow Cu, StratoCu) Cloud Cover

CloudSat/CALIPSO

Cloudtop Droplet Number Concentration

*MODIS (Bennartz, Grosvenor)

Surface Wind (W)

*WindSat, QuikSCAT

Liquid-to-ice transition Temperature/Height

CALIPSO

source: Greg Elsaesser

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

18 of 27

correlation�with mean of�obs metrics

*

*

*

source: Greg Elsaesser

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

19 of 27

19

ModelE3 emulator based on 450 1-year atmosphere runs

Latin Hypercube sampling in a 45-dimensional parameter state space. Lots of empty state space; emulator (neural network) fills in the gaps.

P1

P2

P3

Example Penalty State Space

Transect for any given model metric

source: Marcus van Lier-Walqui

20 of 27

After The Machine

  • photo of white board at GISS

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

21 of 27

Aerosol indirect effect and ECS from E3 candidates

  • AIE from 2000-2010 AMIP runs, PD minus PI offline aerosol for droplet activation only
  • ECS from 30-year Q-flux PI runs

21

source: Andy Ackerman

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

22 of 27

ModelE3 supercooled cloud fraction vs CALIPSO

  • COSP simulator modified to see ”precipitation”
  • note: cloud ice is continuous with precipitating ice (e.g., Fridlind et al. JAS 2012)
  • ”precipitation” also affects cloud feedbacks across ModelE3s

Cesana et al. (GRL 2021, Fig. S6)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

23 of 27

AWARE case study

  • SCM performs quite well
  • stable conditions common�(Silber et al. ACP 2020; GRL 2021)

see Silber et al. [GMD 2022]

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

24 of 27

A new ground-based lidar/radar simulator: EMC2

  • Earth Model Column Collaboratory
  • Python open source, community code base
  • tool to evaluate supercooled cloud fraction, cloud base and surface precipitation, ...

Silber, Jackson, Collis et al. (GMD, 2021)

from surface:

sounding approach

from surface:

lidar approach

from space

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

25 of 27

COSP simulator tested on LES/SCM AWARE case

Cesana et al. (GRL 2021, Fig. S1)

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

26 of 27

Ground-based with EMC2

Stanford et al. (ACP, 2023)

McFarquahar

et al. 2021

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov

27 of 27

SCM

LES

ESM

PARCEL + 1D

satellite + COSP

long-term ground + EMC2

AEROICESTUDY

KIT laboratory SIP

calibration

aircraft field campaigns

long-term ground

development

development

• ESM development workhorse

• pre-calibration tool

• simulator testbed

• cloud feedback analysis tool

• primary and secondary ice formation

+ rain formation and mesoscale structure

+ gravity waves, surface fluxes, ice properties,

aerosol-cloud interactions, ...

DOE ARM Summer School • 22 May 2024 • ann.fridlind@nasa.gov