Diabatic evolution of clouds in a Lagrangian framework: turbulence, vorticity dynamics and precipitation effects

Lead Research Organisation: University of St Andrews
Department Name: Mathematics and Statistics

Abstract

The importance of clouds in weather and climate has long been recognised, yet accurate cloud modelling is immensely difficult due to the wide range of energetic, spatial and temporal scales involved in cloud formation, growth and decay, and the highly nonlinear nature of these turbulent cloud processes. Indeed, the turbulent behaviour of clouds is responsible for many of the uncertainties in atmospheric models, affecting cloud formation and thereby disrupting the global energy balance in climate simulations, as well as the timing and intensity of precipitation in weather forecasts. Weather and climate models fail to resolve the details of the interactions between clouds and their environment and suffer from a crude representation of important microphysical processes, such as rain and snow formation.

These processes can be studied in detail using `large eddy simulation' (LES), a widely-used computational approach employing a fixed discrete grid - a so called `Eulerian' approach. Typically, LES uses grid spacings below 100 metres to resolve such cloud-environment interactions. However, LES still suffers from substantial sensitivity to numerical mixing. This is due to the highly nonlinear nature of the dynamics and thermodynamics of clouds. For the formation of precipitation, regions of high liquid water content are crucial, and numerical mixing can prevent the occurrence of such regions in Eulerian models.

We have recently developed a promising new computational model, MPIC (Moist Parcel-In-Cell), which largely overcomes these numerical errors. MPIC deals with the dynamics of clouds in an essentially `Lagrangian' framework, i.e. by explicitly following parcels of fluid rather than approximating this motion on a fixed grid as in LES. The MPIC model represents both dynamics and processes explicitly using Lagrangian parcels that carry a volume, circulation and thermodynamic properties (e.g. potential temperature and moisture content). This approach accurately preserves key parcel properties and avoids problems due to numerical mixing at small scales that are inherent to Eulerian models.

Our research aims to adapt the MPIC model for use in realistic atmospheric case studies as well as advance our understanding of the fundamental dynamics of cloud turbulence. In particular, we will investigate how discontinuous diabatic forcing associated with condensation and evaporation modifies turbulence relative to adiabatic conditions. This research presents a novel opportunity to improve our theoretical understanding of generic features of cloud processes. It is especially timely due to the development of a massively parallel version of MPIC, which can resolve detailed processes far beyond the reach of existing numerical models.

The potential impacts of this research are vast. Since entrainment in convective clouds strongly influences heavy rainfall, new discoveries about entrainment should enable improved flood forecasting, potentially saving lives, property and business assets. Beyond this project, aspects of the MPIC model may be adopted by numerical weather models to improve the representation of cloud processes and rainfall events. The pathway to this is first to engage scientists at the Met Office, who have agreed to be a project partner. Aspects of the MPIC model will need to be integrated into the latest Met Office model, and the proposed research will make significant progress towards this goal.

The proposed research will make substantial advancements in our understanding of fundamental turbulent processes, ultimately leading to improved weather and global climate models. Extensions to the MPIC model will make it well-suited to a range of geophysical applications including cloud simulations and atmospheric chemistry experiments, making it attractive to a wide scientific and industrial audience.

Planned Impact

The proposed research will exploit a recent computational modelling breakthrough, developed in an EPSRC-funded feasibility study and eCSE project, to advance our understanding of clouds, a critically important component of weather, climate and environmental modelling.

The work will benefit the community of scientists studying a variety of problems related to turbulent clouds and their transport and mixing properties, and using models to make environmental predictions. These include challenges of very high socio-economic importance, notably prediction of extreme rain, global climate modelling, and pollution dispersion. This community will benefit through the provision of the model, and associated research demonstrating its capabilities.

We will engage scientists in universities and forecasting centres by holding a series of workshops, which will train participants to use the modelling approach and develop further applications. We will make all software available through a public online repository, together with documentation and examples of use. We will also promote MPIC using detailed visualisations developed by a trained software developer.

Research conducted as part of the proposal, such as studying entrainment occurring at the edges of clouds, is directly relevant to the representation of clouds in Numerical Weather Prediction models. A particular application is improving forecasts for heavy precipitation, which is strongly influenced by entrainment. Our work will contribute to improved forecasts: earlier warnings save lives, property and business assets. The annual economic benefit to the UK of the Met Office and its forecasts has been calculated to be £3 billion.

While our research focuses on mathematical aspects of cloud evolution, the changes to the MPIC model that we will implement will make it possible in principle to embed MPIC in a larger-scale model for fine-scale forecasting. There is further potential for developing the model for use in weather and climate prediction, as a subgrid model or as part of the dynamical core. These ideas will be developed with our project partner, the Met Office. While development of such a model in the timescale of this project is unrealistic, the approach offers new directions for the years to come, when pressures on computing resources will demand new ideas.

The proposed framework is particularly suitable for problems involving aerosols or chemistry, as it models advection very efficiently. Aerosol-induced radiative forcing and cloud feedbacks are the largest cause of uncertainty in current estimates of climate sensitivity (IPCC report, 2014). A more accurate representation of cloud processes in climate models will allow for better quantitative arguments to be made to influence mitigation strategies to be adopted globally.

Worldwide, air pollution is known to cause a significant number of deaths. Being able to model the transport and dispersion of pollutants out of the near-surface atmospheric layer is an important part of the problem. Using the modelling approach proposed, it will be possible to efficiently model the motion and behaviour of the pollutants in clouds. The model could be further developed for real-time dispersion modelling. The pathway to this impact will be through scientists at the Met Office and the National Centre for Atmospheric Sciences in Leeds.

The applicability of our approach may extend much further and have impacts in areas far from cloud turbulence. One topical area concerns the mixing of bio-geochemical tracers by turbulent fluid motions in the oceanic sub-mesoscales (100m-10km). Our approach may enable a much more accurate representation of this flow, which would have significant impact. The pathway to impact is first through engagement with experts, already under way, and subsequently with interested stakeholders.

Publications

10 25 50

publication icon
Frey M (2023) The 3D Elliptical Parcel-In-Cell (EPIC) method in Journal of Computational Physics: X

publication icon
Frey M (2022) EPIC: The Elliptical Parcel-In-Cell method in Journal of Computational Physics: X

 
Description We developed a major extension to a computational method for simulating the complex dynamics of clouds in the atmosphere. The novelty of the method is it use of deformable ellipsoidal parcels which move with the flow and contain all of the essential information needed (temperature, humidity, etc) to represent the full three-dimensional, time-evolving flow. The method offers major advantages over existing methods for studying clouds, including the critically important ability to faithfully represent very small scale turbulent motion. A paper is expected to be submitted shortly to the Journal of Computational Physics X, where we published a previous prototype method developed for two-dimensional flows (which are assumed independent of one coordinate direction).
Exploitation Route The associated software will be made available to the wider atmospheric dynamics community upon publication of the paper. Moreover, the method can be further extended into domains outside the original scope of this project, e.g. to magnetised flows in the sun. With little modification, the method could be adapted to study the upper ocean, whose complex mixing processes are vital to understanding climate.
Sectors Digital/Communication/Information Technologies (including Software),Environment

 
Description We have been involved in detailed discussions with the UK Met Office (UKMO), specifically the Convection Division led by Dr Alison Stirling, most recently in early December 2022. The UKMO is interested in applying our new computational method, and we have jointly worked out the best way forward to extend our method for realistic applications in weather and climate. This is an evolving impact in its early stages.
First Year Of Impact 2022
Sector Digital/Communication/Information Technologies (including Software),Environment
Impact Types Societal,Policy & public services

 
Title EPIC - the Elliptical Parcel In Cell method 
Description A new numerical method for simulating 2D stratified flows under the Boussinesq approximation. Uses parcels to represent all quantities and an underlying grid for efficiency. 
Type Of Technology Software 
Year Produced 2022 
Open Source License? Yes  
Impact This method has been shown to greatly improve the modelling of density stratified flows, reaching higher accuracy and efficiency than standard methods. 
 
Title Pseudospectral code repository 
Description This is a general set of numerical codes for simulating two-dimensional stratified flows under the Boussinesq approximation. The domain is periodic in x, but bounded by free slip walls in y. 
Type Of Technology Software 
Year Produced 2022 
Open Source License? Yes  
Impact This software was used to test the new EPIC method for simulating two-dimensional stratified flows. 
 
Title matt-frey/epic: EPIC version 0.12.3 
Description Updated version of EPIC (Elliptical Parcel-in-Cell) method. 
Type Of Technology Software 
Year Produced 2022 
Open Source License? Yes  
Impact The 2D version of the model as used in a publication (https://doi.org/10.1016/j.jcpx.2022.100109). We will soon submit an article with the 3D version to a journal. 
URL https://zenodo.org/record/6997439#.ZAG67YDP2V4
 
Title matt-frey/ps3d: ps3d version 0.0.11 
Description Pseudo-spectral code for turbulent flows in 3D. 
Type Of Technology Software 
Year Produced 2023 
Open Source License? Yes  
Impact The code was used in the publication "The stability of inviscid Beltrami flow between parallel free-slip impermeable boundaries" (https://doi.org/10.1017/jfm.2022.1007). We also use this software in a recent work which will be submitted to JCP-X soon. 
URL https://doi.org/10.5281/zenodo.7694990
 
Description A fully Lagrangian approach to modelling convection (5th Convective scale workshop) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Lighning talk on EPIC at the 5th Convective scale workshop (Convection processes session), organised by the Met Office Unified Model partnership
Year(s) Of Engagement Activity 2022
 
Description Clouds and coherent structures (Coherent Structures in Astro-Geo-Turbulence Workshop, Flatiron Institute) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact A group of experts attended this talk, which was part of a workshop in preparation for a grant submission.
Year(s) Of Engagement Activity 2022
 
Description Essentially Lagrangian simulation of turbulence, mixing and phase transitions with the Ellipsoidal Parcel-in-Cell model (Leipzig Meteorological Colloquium) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Seminar talk at the Leipzig Meteorological Colloquium (held jointly between the Leibniz Institute for Tropospheric Research and University of Leipzig). About 15 participants (partially online).
Year(s) Of Engagement Activity 2022
 
Description Seminar given by David Dritschel entitled "A new Lagrangian approach to modelling turbulent flows", given as part of the Isaac Newton Institute programme "Advances in geophysical and astrophysical turbulence (TURW05)" on 18 May 2022 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact This seminar was given at the Isaac Newton Institute where I was invited to participate in a programme on turbulence attended by a wide range of international experts in the field. There was much interest in the talk, in particular people were surprised that our new approach could be so effective in modelling cloud dynamics.
Year(s) Of Engagement Activity 2022
 
Description Seminar given by David Dritschel entitled "What is a cloud?" given at the School Research Day on 1 December 2022 at the University of St Andrews 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Professional Practitioners
Results and Impact The entire School of Mathematics and Statistics attended a School Research Day where I was able to showcase our grant activities and research to date.
Year(s) Of Engagement Activity 2022
 
Description Seminar given by M Frey at Paul Scherrer Institute, "Recent advances in the modelling of incompressible turbulent two and three-dimensional flows", 12 October 2022 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Seminar given at Paul Scherrer Institute (Switzerland) in the Laboratory for Simulation and Modelling.
Year(s) Of Engagement Activity 2022
 
Description Seminar given by M Frey at the University of Exeter, "EPIC: the Elliptical Parcel-In-Cell Method", 25 January 2022 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Seminar given at the University of Exeter, which is a leading centre for atmospheric research, and which is strongly engaged in analogous modelling activities.
Year(s) Of Engagement Activity 2022
 
Description Talk by M Frey at ParaCon Plenary Meeting, "EPIC: the Elliptical Parcel-In-Cell Method", 13 December 2021 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Talk at the plenary meeting of another major UK project investigating cloud convection, attended by experts from many universities and the Met Office.
Year(s) Of Engagement Activity 2021
 
Description Talk given by D Dritschel at the Plenary Meeting of the European Project ASTuS (Paris), "EPIC: the Elliptical-Parcel-In-Cell Method", 14 January 2022 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact This talk reached a new audience engaged in a European Project "ASTuS" studying smoke clouds associated with major fires and their influence on the climate. Dritschel is associated with this European project, and the EPIC method developed in the present EPSRC project is likely to be useful for simulating these smoke clouds, whose dynamics are very much unknown.
Year(s) Of Engagement Activity 2022
 
Description Talk given by M Frey at the 34th Scottish Fluid Mechanics Meeting, "EPIC - Elliptic Parcel-in-Cell", 27 May 2021 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Talk at a major venue for fluid dynamics research. Significant interest was shown by those attending.
Year(s) Of Engagement Activity 2021
 
Description Talk given by M Frey at the project plenary meeting, "EPIC - Elliptic Parcel-in-Cell", 25 May 2021 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Talk presenting results of the grant to date to the plenary meeting. This reached the Met Office who is a project partner.
Year(s) Of Engagement Activity 2021
 
Description Talk given by M Frey at the project plenary meeting, "EPIC - Elliptic Parcel-in-Cell", 27 November 2021 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Talk at a project plenary meeting discussing progress so far. Attended by the Met Office.
Year(s) Of Engagement Activity 2021
 
Description Talk given by M Frey at the project plenary meeting, "EPIC: the Elliptical Parcel-In-Cell Method", 28 November 2022 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Talk at a project plenary meeting discussing progress so far. Attended by the Met Office.
Year(s) Of Engagement Activity 2022
 
Description Talk given by M Frey, D Dritschel and S Boeing at the project plenary meeting, "EPIC: the Elliptical Parcel-In-Cell Method", 28 November 2022 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Talk at a project plenary meeting discussing progress so far. Attended by the Met Office.
Year(s) Of Engagement Activity 2022
 
Description Vortex Group Seminar given by M Frey in St Andrews, "EPIC: the Elliptical Parcel-In-Cell Method", 24 November 2021 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Professional Practitioners
Results and Impact Talk given to my local research group discussing the latest developments in our project.
Year(s) Of Engagement Activity 2021
 
Description Vortex Group Seminar given by M Frey in St Andrews, "EPIC: the Elliptical Parcel-In-Cell Method", 6 April 2022 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Professional Practitioners
Results and Impact Talk given to Prof Dr David Dritschel's local research group discussing the latest developments in our project
Year(s) Of Engagement Activity 2022
 
Description Vortex Group Seminar given by M Frey in St Andrews: "EPIC - Elliptic Parcel-in-Cell", 28 April 2021 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Professional Practitioners
Results and Impact Description of new numerical method to local research group (around 10 people)
Year(s) Of Engagement Activity 2021