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
Eulerian microhysics available liquid water is divided into histogram bins N r
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
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
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...
Lagrangian microphysics superdroplets in the domain
Lagrangian microphysics superdroplets in the domain attributes:
Lagrangian microphysics superdroplets in the domain attributes: location
Lagrangian microphysics superdroplets in the domain attributes: location wet radius
Lagrangian microphysics superdroplets in the domain attributes: location wet radius dry radius
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 *
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 *
Cloud microphysics MaxwellMason equation of condensational growth for each superdroplet using κkoehler parametrisation of higroscopicity (Petters & Kreidenweis, 2007) CCN activation condensational growth
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
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
Example results with collisions (2D kinematic setup) aerosol concentration [mg ] singleeddy velocity field 40 20 80 60 40 20 0 setup: Grabowski & Lebo (ICMW 202) 2D prescribed flow advection: libmpdata++ (2pass FCT) µphysics: libcloudph++
Example results with collisions (2D kinematic setup) x cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ].4 40 20 0.8 0.4 0. 80 60 40 0.2 20 0 0
Example results with collisions (2D kinematic setup) x o cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ].4 40 20 0.8 0.4 0. 80 60 40 0.2 20 0 0
Example results with collisions (2D kinematic setup) x o o cloud water mixing ratio [g/kg] rain drop spec. conc. [mg ] aerosol concentration [mg ].4 40 20 0.8 0.4 0. 80 60 40 0.2 20 0 0
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 ].4 40 20 0.8 0.4 0. 80 60 40 0.2 20 0 0
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 ].4 40 20 0.8 0.4 0. 80 60 40 0.2 20 0 0
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) 0. 0 0. 0. particle radius [ m] wet radius dry radius.4 0.8 0.4 0.2 0 [mg m ] [mg m ] (b) 0 wet radius dry radius rain drop spec. conc. [mg ] (d) 0. 0 0. 0. particle radius [ m] 0. wet radius dry radius aerosol concentration [mg ] a b e c g f d h i j 40 20 80 60 40 20 0 [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]
representing incloud oxidation of sulfur in a particlebased cloudmicrophysics scheme CCN activation condensational growth collisional growth precipitation wet deposition droplet deactivation
representing incloud oxidation of sulfur in a particlebased cloudmicrophysics scheme CCN activation condensational growth collisional growth aqueous chemistry precipitation wet deposition droplet deactivation
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 *
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
Aqueous chemistry dissolving of trace gases SO 2 O 3 HNO CO 3 2 NH 3 H 2 O 2
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 3 + 2 NH 4 CO 3
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 3 + 2 NH 4 CO 3 3 H 2 O 2 HSO 4 SO 4 2
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 3 + 2 NH 4 CO 3 3 H 2 O 2 HSO 4 SO 4 2
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 3 + 2 NH 4 CO 3 3 H 2 O 2 HSO 4 SO 4 2
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 3 + 2 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
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 3 + 2 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
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 3 + 2 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)
Validation of the chemistry module adiabatic parcel simulations from model comparison study by Kreidenweis et al 2003
Validation of the chemistry module adiabatic parcel simulations from model comparison study by Kreidenweis et al 2003 initial lognormal monomodal dry aerosol size distribution
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)] 0. 0. particle diameter [μm]
Example results with chemistry (2D kinematic setup) aerosol concentration [mg ] singleeddy velocity field e f 40 20 g h i j 80 60 40 20 0 setup: Grabowski & Lebo (ICMW 202) 2D prescribed flow advection: libmpdata++ (2pass FCT) µphysics: libcloudph++
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 20 8 6 4 2 8 6 e i g f j h 0. 9 8 7 6 5 e i g f j h 40 20 80 60 40 4 2 4 3 20 0
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 20 8 6 4 2 8 6 e i g f j h 0. 9 8 7 6 5 e i g f j h 40 20 80 60 40 4 2 4 3 20 0
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 20 8 6 4 2 8 6 e i g f j h 0. 9 8 7 6 5 e i g f j h 40 20 80 60 40 4 2 4 3 20 0
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 20 8 6 4 2 8 6 e i g f j h 0. 9 8 7 6 5 e i g f j h 40 20 80 60 40 4 2 4 3 20 0
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 20 8 6 4 2 8 6 e i g f j h 0. 9 8 7 6 5 e i g f j h 40 20 80 60 40 4 2 4 3 20 0
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 20 8 6 4 2 8 6 4 2 [mg μm ] [mg μm ] (f) 0 (h) 0. 0 0. e i dry radius [μm] g 0. particle radius [μm] j wet radius dry radius f wet radius dry radius h 0. 9 8 7 6 5 4 3 aerosol concentration [mg ] e f g h i j 40 20 80 60 40 20 0 0. [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]
Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes http://libcloudphxx.igf.fuw.edu.pl/
Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes http://libcloudphxx.igf.fuw.edu.pl/ 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
Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes http://libcloudphxx.igf.fuw.edu.pl/ 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)
Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes http://libcloudphxx.igf.fuw.edu.pl/ 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)
Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes http://libcloudphxx.igf.fuw.edu.pl/ 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
Final remarks available online as a part of libcloudph++ library of cloud microphysics schemes http://libcloudphxx.igf.fuw.edu.pl/ 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!