representing in-cloud oxidation of sulfur in a particle-based cloud-microphysics scheme

Transcription

1 representing incloud oxidation of sulfur in a particlebased cloudmicrophysics scheme Anna Jaruga, Hanna Paw lowska, Sylwester Arabas Institute of Geophysics Faculty of Physics, University of Warsaw 8 th Conference on Atmospheric Chemistry AMS, New Orleans, January 206

2 Eulerian microhysics available liquid water is divided into histogram bins N r

3 Eulerian microhysics N available liquid water is divided into histogram bins for aerosol cloud interaction studies 2D histogram is created with wet and dry radius r r d

4 Eulerian microhysics N available liquid water is divided into histogram bins for aerosol cloud interaction studies 2D histogram is created with wet and dry radius parametrisations of source and sink terms are applied to each bin r r d

5 Eulerian microhysics r d N r available liquid water is divided into histogram bins for aerosol cloud interaction studies 2D histogram is created with wet and dry radius parametrisations of source and sink terms are applied to each bin new histogram dimensions needed for chemical compounds...

6 Lagrangian microphysics superdroplets in the domain

7 Lagrangian microphysics superdroplets in the domain attributes:

8 Lagrangian microphysics superdroplets in the domain attributes: location

9 Lagrangian microphysics superdroplets in the domain attributes: location wet radius

10 Lagrangian microphysics superdroplets in the domain attributes: location wet radius dry radius

11 Lagrangian microphysics 3 * 28 * 9 * 7 * * 52 * 25 * 57 * 504 * 5 * 88 * superdroplets in the domain attributes: location wet radius dry radius multiplicity 5 * 75 *

12 Lagrangian microphysics 3 * 28 * 9 * 7 * * 52 * 25 * 57 * 504 * 5 * 88 * superdroplets in the domain attributes: location wet radius dry radius multiplicity... 5 * 75 *

13 Cloud microphysics MaxwellMason equation of condensational growth for each superdroplet using κkoehler parametrisation of higroscopicity (Petters & Kreidenweis, 2007) CCN activation condensational growth

14 Cloud microphysics MaxwellMason equation of condensational growth for each superdroplet using κkoehler parametrisation of higroscopicity (Petters & Kreidenweis, 2007) collisions for each superdroplet by a MonteCarlo coalescence scheme (Shima et al. 2009) CCN activation condensational growth collision coalescence

15 Cloud microphysics MaxwellMason equation of condensational growth for each superdroplet using κkoehler parametrisation of higroscopicity (Petters & Kreidenweis, 2007) collisions for each superdroplet by a MonteCarlo coalescence scheme (Shima et al. 2009) sedimentation of each superdroplet (Khvorostyanov & Curry, 2002) CCN activation condensational growth collision coalescence precipitation wet deposition droplet deactivation

16 Example results with collisions (2D kinematic setup) aerosol concentration [mg ] singleeddy velocity field setup: Grabowski & Lebo (ICMW 202) 2D prescribed flow advection: libmpdata++ (2pass FCT) µphysics: libcloudph++

17 Example results with collisions (2D kinematic setup) x cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ]

18 Example results with collisions (2D kinematic setup) x o cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ]

19 Example results with collisions (2D kinematic setup) x o o cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ]

20 Example results with collisions (2D kinematic setup) x o o o cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ]

21 Example results with collisions (2D kinematic setup) x o o o o cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ]

22 2 2 cell particlederived spectra x o o(a) o o o [mg m ] [mg m ] 0 wet radius dry radius cloud water mixing ratio [g/kg] (c) particle radius [ m] wet radius dry radius [mg m ] [mg m ] (b) 0 wet radius dry radius rain drop spec. conc. [mg ] (d) particle radius [ m] 0. wet radius dry radius aerosol concentration [mg ] a b e c g f d h i j [mg m ] (e) 0 0. (g) 0. wet radius dry radius 0. particle radius [ m] [mg m ] (f) 0 0. particle radius [ m] (h) 0. wet radius dry radius 0. particle radius [ m] particle radius [ m]

23 representing incloud oxidation of sulfur in a particlebased cloudmicrophysics scheme CCN activation condensational growth collisional growth precipitation wet deposition droplet deactivation

24 representing incloud oxidation of sulfur in a particlebased cloudmicrophysics scheme CCN activation condensational growth collisional growth aqueous chemistry precipitation wet deposition droplet deactivation

25 Lagrangian microphysics + aqueous chemistry 3 * 28 * 9 * 7 * * 52 * 25 * 57 * 504 * 5 * 88 * superdroplets in the domain with attributes: location wet radius dry radius multiplicity 5 * 75 *

26 Lagrangian microphysics + aqueous chemistry 3 * 28 * 5 * 9 * 7 * * 52 * 75 * 25 * 57 * 504 * 5 * 88 * superdroplets in the domain with attributes: location wet radius dry radius multiplicity mass of chemical compounds within droplets: H 2 O SO 2, O 3, H 2 O 2, H 2 O CO 2, H 2 O NH 3, HNO 3, HSO 3, SO 3 2, HCO 3, CO 3 2 NO 3, NH 4 +, H +, OH, HSO 4, SO 4 2

27 Aqueous chemistry dissolving of trace gases SO 2 O 3 HNO CO 3 2 NH 3 H 2 O 2

28 Aqueous chemistry dissolving of trace gases dissociation SO 2 O 3 HNO CO 3 2 NH 3 H 2 O 2 SO 3 2 HSO 3 HCO 3 NO NH 4 CO 3

29 Aqueous chemistry dissolving of trace gases dissociation oxidation SO 2 O 3 HNO CO 3 2 NH 3 H 2 O 2 SO 3 2 HSO 3 HCO 3 NO NH 4 CO 3 3 H 2 O 2 HSO 4 SO 4 2

30 Aqueous chemistry dissolving of trace gases dissociation oxidation gas uptake is treated as nonequilibrum process SO 2 O 3 HNO CO 3 2 NH 3 H 2 O 2 SO 3 2 HSO 3 HCO 3 NO NH 4 CO 3 3 H 2 O 2 HSO 4 SO 4 2

31 Aqueous chemistry dissolving of trace gases dissociation oxidation O SO 3 2 HNO CO 3 2 gas uptake is treated as nonequilibrum process gasliquid and liquidliquid equilibra are computed for all superdroplets NH 3 H 2 O 2 SO 3 2 HSO 3 HCO 3 NO NH 4 CO 3 3 H 2 O 2 HSO 4 SO 4 2

32 Aqueous chemistry dissolving of trace gases dissociation oxidation O SO 3 2 HNO CO 3 2 NH 3 H 2 O 2 HSO 3 SO 3 2 HCO 3 NO NH 4 CO 3 gas uptake is treated as nonequilibrum process gasliquid and liquidliquid equilibra are computed for all superdroplets oxidation of sulfur is computed only for cloud droplets 3 H 2 O 2 HSO 4 SO 4 2

33 Aqueous chemistry dissolving of trace gases dissociation oxidation O SO 3 2 HNO CO 3 2 NH 3 H 2 O 2 HSO 3 SO 3 2 HCO 3 NO NH 4 CO 3 gas uptake is treated as nonequilibrum process gasliquid and liquidliquid equilibra are computed for all superdroplets oxidation of sulfur is computed only for cloud droplets no adjustment for highionic strength of some droplets 3 H 2 O 2 HSO 4 SO 4 2

34 Aqueous chemistry 3 dissolving of trace gases dissociation oxidation O SO 3 2 HNO CO 3 2 NH 3 H 2 O 2 HSO 3 SO 3 2 HCO 3 NO NH 4 CO 3 H 2 O 2 HSO 4 SO 4 2 gas uptake is treated as nonequilibrum process gasliquid and liquidliquid equilibra are computed for all superdroplets oxidation of sulfur is computed only for cloud droplets no adjustment for highionic strength of some droplets impacts condensation via dry radii, but no impact on κ (yet)

35 Validation of the chemistry module adiabatic parcel simulations from model comparison study by Kreidenweis et al 2003

36 Validation of the chemistry module adiabatic parcel simulations from model comparison study by Kreidenweis et al 2003 initial lognormal monomodal dry aerosol size distribution

37 Validation of the chemistry module adiabatic parcel simulations from model comparison study by Kreidenweis et al 2003 initial lognormal monomodal dry aerosol size distribution ds(vi)/dlog (D) [μg /m 3 log (size interval)] particle diameter [μm]

38 Example results with chemistry (2D kinematic setup) aerosol concentration [mg ] singleeddy velocity field e f g h i j setup: Grabowski & Lebo (ICMW 202) 2D prescribed flow advection: libmpdata++ (2pass FCT) µphysics: libcloudph++

39 Example results with chemistry (2D kinematic setup) x cloud droplet effective radius [μm] dry radius [μm] aerosol concentration [mg ] e i g f j h e i g f j h e i g f j h

40 Example results with chemistry (2D kinematic setup) x o cloud droplet effective radius [μm] dry radius [μm] aerosol concentration [mg ] e i g f j h e i g f j h e i g f j h

41 Example results with chemistry (2D kinematic setup) x o o cloud droplet effective radius [μm] dry radius [μm] aerosol concentration [mg ] e i g f j h e i g f j h e i g f j h

42 Example results with chemistry (2D kinematic setup) x o o o cloud droplet effective radius [μm] dry radius [μm] aerosol concentration [mg ] e i g f j h e i g f j h e i g f j h

43 Example results with chemistry (2D kinematic setup) x o o o o cloud droplet effective radius [μm] dry radius [μm] aerosol concentration [mg ] e i g f j h e i g f j h e i g f j h

44 2 2 cell particlederived spectra x o o(e) o o o [mg μm ] [mg μm ] wet radius 0 dry radius cloud droplet effective radius [μm] e 0. i (g) 0 0. g 0. f h particle radius [μm] j wet radius dry radius [mg μm ] [mg μm ] (f) 0 (h) e i dry radius [μm] g 0. particle radius [μm] j wet radius dry radius f wet radius dry radius h aerosol concentration [mg ] e f g h i j [mg μm ] (i) 0 wet radius dry radius 0. particle radius [μm] 0. [mg μm ] (j) 0 wet radius dry radius 0. particle radius [μm] 0. particle radius [μm] 0. particle radius [μm]

45 Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes

46 Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes reusable: design: no assumptions on dimensionality or dyncore type documentation: API described in the paper/manual legal/practical matters: open source, GPL, hosted on github

47 Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes reusable: design: no assumptions on dimensionality or dyncore type documentation: API described in the paper/manual legal/practical matters: open source, GPL, hosted on github Lagrangian scheme optionally GPUresident (via Thrust)

48 Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes reusable: design: no assumptions on dimensionality or dyncore type documentation: API described in the paper/manual legal/practical matters: open source, GPL, hosted on github Lagrangian scheme optionally GPUresident (via Thrust) compact code (5000 lines of code)

49 Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes reusable: design: no assumptions on dimensionality or dyncore type documentation: API described in the paper/manual legal/practical matters: open source, GPL, hosted on github Lagrangian scheme optionally GPUresident (via Thrust) compact code (5000 lines of code) equipped with Python bindings

50 Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes reusable: design: no assumptions on dimensionality or dyncore type documentation: API described in the paper/manual legal/practical matters: open source, GPL, hosted on github Lagrangian scheme optionally GPUresident (via Thrust) compact code (5000 lines of code) equipped with Python bindings Thank you for your attention!