Aerosol processing by clouds - development of a multifaceted object-oriented numerical simulation framework
|prof. dr hab. Hanna Pawłowska||Project leader|
|dr Sylwester Arabas||Co-investigator|
|dr Piotr Dziekan||Co-investigator|
|dr Anna Jaruga||Co-investigator|
|dr Maciej Waruszewski||Co-investigator|
Research project objectives
The project aims at developing a numerical model for studying aerosol processing by clouds. It responds to the pressing need for the development of effective tools for research of aerosol processing by clouds, for which even the most sophisticated traditional cloud microphysics approaches are either inadequate or impractical.
Aerosols and clouds have been identified to contribute significantly to the Earth's radiation budget and therefore play an important, yet still poorly understood, role in the radiation imbalance of projected climate variations (IPCC, 2007). Clouds process atmospheric aerosol by taking part in its wet deposition and cloud condensation nuclei regeneration through evaporation of cloud droplets and drizzle. Understanding the way clouds process aerosol is needed to comprehensively describe the interactions between aerosol, clouds, precipitation and radiation. The complexity of these interactions makes numerical simulation the only viable tool for investigating how natural and anthropogenic air pollution impacts the environment.
A Large Eddy Simulation (LES) framework capable of representing comprehensively processes leading to aerosol processing by clouds will be constructed by extending the existing cloud-simulation software package (icicle) being developed at the Institute of Geophysics, Faculty of Physics at the University of Warsaw (IGFUW). Up to now, the icicle framework allows performing only two- dimensional kinematic simulations with a prescribed flow field using a novel particle-based aerosol-cloud microphysics. The newly-developed icicle modules will constitute new implementations of the key components of the NCAR LES-type anelastic cloud model EULAG. The pioneering nature of the project is embodied in the coupling of an efficient description of cloud dynamics (EULAG’s MPDATA- based anelastic framework) with a novel treatment of aerosol-cloud-precipitation microphysics (the Super-Droplet method) within a robust numerical modelling tool. Neither MPDATA-based anelastic solvers nor particle-based microphysics modules are commonly available in presently developed models for cloud research, partially owing to the complexity of these algorithms.
It is planned also to extend the icicle implementation of the Super-Droplet scheme with representation of in-droplet aqueous chemical reactions.
The proposed development will be carried out in compliance with free and open-source software engineering practices employing modern coding techniques.
Research project impact
Coupling detailed particle-based microphysics with a state-of-the-art fluid dynamics solver within an open-source tool engineered using modern programming techniques will bring to the community a novel and effective tool for studying aerosol-cloud interactions. The openness and clarity of the source code, together with the completeness of the documentation will make the code auditable, and reusable by the community making the developed software a tool for genuinely reproducible research.
The proposed development aims at advancing both research and education. It will generate novel and important results readily publishable in peer-reviewed high-impact journals, and provide documented computer programs suitable for future Masters's and PhD projects.
The new icicle modules will be developed in collaboration with Drs Piotr Smolarkiewicz and Wojciech Grabowski from the National Center for Atmospheric Research (NCAR) in Boulder, Colorado, USA. NCAR partners are internationally recognized experts in applications of computational fluid dynamics, with a wide-ranging portfolio of applications to all scale atmospheric problems. They both have a record of close collaboration with the Institute of Geophysics (IGFUW), Faculty of Physics, and University of Warsaw (their joint Alma Mater).
The current proposal is a consequent step in realizing a synergistic program, mutually beneficial to NCAR and IGFUW and begins a new chapter in the NCAR-IGFUW collaboration by starting a first-join software-development endeavour.
- Hill, A.A., Lebo Z.J., Andrejczuk M., Arabas S., Dziekan P., Field P., Gettelman A., Hoffmann F., Pawlowska H., Onishi R. and Vie B., 2023, Toward a numerical benchmark for warm rain processes, Journal of the Atmospheric Sciences, 10.1175/JASD- 21-0275.1.
- Grabowski W. W., Dziekan P., and Pawlowska H., 2018, Lagrangian condensation microphysics with Twomey CCN activation, Geoscientific Model Development, vol. 11, pp. 103-120, 10.5194/gmd-11-1-2018
- Jaruga A., and Pawlowska H., 2018, libcloudph++ 2.0: aqueous-phase chemistry extension of particle-based cloud microphysics scheme, Geoscientific Model Development, vol. 11, pp. 3623-3645, 10.5194/gmd-11-3623-2018
- Waruszewski, M., Kühnlein, C., Pawlowska, H., Smolarkiewicz, P.K., 2018, MPDATA: Third-order accuracy for variable flows, Journal of Computational Physics, vol. 359, pp. 361-379, 10.1016/j.jcp.2018.01.005
- Arabas S. and Shima S., 2017, On the CCN (de)activation nonlinearities, Nonlinear Processes in Geophysics, vol. 24, pp. 535-542, 10.5194/npg-24-535-2017
- Dziekan P. and Pawlowska H., 2017, Stochastic coalescence in Lagrangian cloud microphysics, Atmospheric Chemistry and Physics, vol. 17 (2), pp. 13509-13520, 10.5194/acp-17-13509-2017
- Arabas S., Jaruga A., Pawlowska H., and Grabowski W. W., 2015, libcloudph++ 1.0: a single-moment bulk, double-moment bulk, and particle-based warm-rain microphysics library in C++, Geoscientific Model Development, vol. 8, pp. 1677-1707, doi: 10.5194/gmd-8-1677-2015
- Jarecka D., Jaruga A., Smolarkiewicz P.K., 2015, A spreading drop of shallow water, Journal of Computational Physics, vol. 289, pp. 53–61, 10.1016/j.jcp.2015.02.003
- Jaruga A., Arabas S., Jarecka D., Pawlowska H., Smolarkiewicz P. K., Waruszewski M., 2015, libmpdata++ 1.0: a library of parallel MPDATA solvers for systems of generalised transport equations, Geoscientific Model Development, vol. 8, pp. 1005-1032, doi:10.5194/gmd-8-1005-2015