Adaptive Mesh Modelling of the Global Atmosphere

Lead Research Organisation: University of Reading
Department Name: Meteorology

Abstract

Current weather and climate forecasting models use fixed, uniform resolution, so that all regions are represented using the same grid spacing. This has enabled efficient, valuable forecasts, but it may be possible to improve the representation of the atmosphere by varying the grid spacing depending on the local weather conditions, that is, adaptive mesh modelling. This could improve the representation of, for example, heavy precipitation or mountain ranges. The overarching themes of this research are to try to make adaptive mesh modelling sufficiently accurate and efficient for operational use and give a realistic assessment of whether this is possible. It is first proposed to study vertical discretisation; the splitting of the atmosphere vertical layers around mountains. Climate models currently use terrain following layers which are draped over mountains like blankets, becoming squashed at summits and spaced out in valleys. This is accurate for the coarse resolution used in climate models since the layers are never very steep. Once resolution increases, steeper slopes are simulated and accuracy can be compromised. Therefore cut cells have been used in some small scale models, whereby the mountains break through flat layers. Cut cells have different accuracy compromises. Novel hybrids between terrain following layers and cut cells will be created in order to minimise both sources of error and compared with discretisations which do not use layers. If the computational mesh can be aligned with the flow, with wide mesh spacing in directions in which conditions change little and fine spacing in directions in which conditions change rapidly, then accuracy can be maximised for the cost. The creation of these anisotropic meshes first requires anisotropic refinement criteria and then an algorithm to generate a mesh with the desired properties. Both of these will be addressed. There has been little work done on where to increase resolution in order to improve atmospheric simulation using the minimum extra computing resources, ie mesh refinement criteria. This is a challenging topic due to the complex inter-dependencies between the weather in different regions and the impossibility of resolving everything in a single global model, even with adaptive meshing. A new anticipative mesh refinement technique will be developed, involving a coarse resolution prediction of where resolution will be needed followed by the adapted mesh simulation. The refined mesh will therefore be in place before small scale features appear, improving accuracy. Cost will be reduced since the expensive mesh adaptation needs only to be done infrequently. Clouds and precipitation are crucial aspects of atmospheric circulation and can be influential on the weather thousands of kilometres away. However they are usually predicted poorly in climate models, with some improvement in higher resolution weather forecasting models. If adaptive meshes can be used to resolve the deepest clouds associated with deep convection, then atmospheric simulation may be improved. However criteria to identify in advance where convection will need high resolution have not before been created. In climate models, deep convection cannot be resolved so schemes have been developed which estimate in which grid boxes sub-grid scale convection may occur. In this project it is planned to use criteria from existing convection schemes in order to enable refinement of deep convection before it breaks out. The work on vertical discretisation is necessary for higher horizontal resolution around mountains with either adaptive meshes or static grids. The other work involves creating optimal meshes for adaptive modelling, with optimal cell shapes and orientations, mesh refinement frequency and criteria and novel strategies for resolving deep convection before it breaks out. These achievements could greatly improve the fidelity and competitiveness of adaptive mesh modelling of the global atmosphere.

Publications

10 25 50
publication icon
Chen Y (2017) Comparison of dimensionally split and multi-dimensional atmospheric transport schemes for long time steps in Quarterly Journal of the Royal Meteorological Society

publication icon
Weller H (2014) Non-orthogonal version of the arbitrary polygonal C-grid and a new diamond grid in Geoscientific Model Development

 
Description I have made some major advances to numerical methods for models of the global atmosphere relating to the efficiency of model on modern computers, time-stepping techniques, the accurate representation of mountains and techniques for including adaptive meshing.

I have developed new numerical methods suitable for quasi-uniform grids of the sphere (ie not latitude-longitude grids) which are necessary for massively parallel computers. Some of the promising developments were found to be unsuitable for weather forecasting due to low accuracy. I identified the reasons for low accuracy of some numerical methods which was important for the Met Office when designing their next model.

As well as needing quasi-uniform grids, next generation weather and climate forecasting models need time-stepping schemes that do not entail too much communication between processor of massively parallel computers. I have created a new technique for combining different time-stepping techniques in order to reduce communication without reducing accuracy, adapting some techniques in the mathematical literature to be directly applicable to atmospheric models.

Operational weather and climate models use terrain following grids, with layers of the grid rising and falling with the terrain. Traditionally this has lead to spurious circulations around mountains in calm conditions. I have developed a new numerical method which guarantees that pressure gradients will not produce these spurious circulations. With a PhD student, I have also created a new kind of cut-cell grid in which grid layers intersect mountains rather than following them. This is not a new idea. Our new contribution is to create a type of cut cell grid that does not need a very small time-step for stability.

I have also started to work on moving meshes for adapting to changing meteorological conditions. Moving meshes rather than h-adaptivity (adding and removing points) should lead to better load balancing on parallel computers and will not lead to jumps in the solution when the mesh is adapted. I have created a new technique for generating adaptive moving meshes on the sphere, solving the Monge-Ampere equation, which has not been solved before on the surface of a sphere.
Exploitation Route My work on accurate forecasts over mountains would be useful for other centres that write weather and climate forecasting models
Sectors Aerospace, Defence and Marine,Environment

 
Description The findings of my research have been used to inform the design of the new Met Office weather and climate forecasting model.
First Year Of Impact 2012
Sector Aerospace, Defence and Marine
Impact Types Policy & public services

 
Description Faculty half studentships
Amount £1 (GBP)
Organisation University of Reading 
Sector Academic/University
Country United Kingdom
Start 10/2012 
End 10/2015
 
Description Isaac Newton Institute Programme on Multiscale Numerics for the Atmosphere and Ocean
Amount £1 (GBP)
Organisation Isaac Newton Institute for Mathematical Sciences 
Sector Academic/University
Country United Kingdom
Start 08/2012 
End 12/2013
 
Title AtmosFOAM 
Description A suite of numerical models and methods for atmospheric modelling using emerging numerical techniques 
Type Of Material Computer model/algorithm 
Year Produced 2014 
Provided To Others? Yes  
Impact Papers and collaborative code development with students 
URL https://github.com/AtmosFOAM/
 
Title Atmospheric Model using a moving mesh 
Description Model to solve the shallow water equations and linear advection equation using a moving mesh. 
Type Of Material Computer model/algorithm 
Year Produced 2017 
Provided To Others? Yes  
Impact None yet 
URL https://github.com/AtmosFOAM/AMMM
 
Description Paracon: Joint project between the UK Met Office and the Universities of Reading, Exeter, Leeds and Cambridge 
Organisation Meteorological Office UK
Country United Kingdom 
Sector Academic/University 
PI Contribution Developing an atmospheric convection scheme using conditional averaging
Collaborator Contribution Understanding convection
Impact None yet
Start Year 2015
 
Description Test cases for multi-scale numerics 
Organisation Naval Postgraduate School, Monterrey CA
Country United States 
Sector Academic/University 
PI Contribution The PI (Hilary Weller) was the lead organiser of the Newton Institute Programme on "Multiscale numerics for the Atmosphere and Ocean" from September to December 2012, attended by over 100 scientists and mathematicians. During this programme, Hilary lead an effort to evaluate local mesh refinement strategies. Mesh refinement by itself will often not improve the simulated climate nor reduce global errors, but instead is used to allow the model to resolve additional processes in select regions. But it remains critically important that the mesh refinement strategy does not harm the larger scales such as geostrophic balance and Rossby waves that are well resolved by global uniform resolution models. In this effort, first a suite of idealized test cases were developed to make it easy to evaluate the impact of mesh refinement on these large scales. The tests were then undertaken by numerous groups and work is ongoing in completing the tests and writing up the work.
Collaborator Contribution Other partners undertook the test cases
Impact gmd-2014-99 Submitted on 31 May 2014 The spectral element method on variable resolution grids: Evaluating grid sensitivity and resolution-aware numerical viscosity O. Guba, M. Taylor, P. Ullrich, J. Overfelt, and M. Levy
Start Year 2012
 
Title AMMM - Adaptive Moving Mesh Methods for Atmospheric Modelling 
Description Code for numerical solution of the equations governing mesh movements and numerical methods for solution of equations of atmospheric motion on moving meshes. 
Type Of Technology Software 
Year Produced 2016 
Open Source License? Yes  
Impact Scientific publications 
URL https://github.com/AtmosFOAM/AMMM
 
Title AtmosFOAM 
Description A set of library routines and applications for simulating the atmosphere using arbitrary meshes 
Type Of Technology Software 
Year Produced 2015 
Open Source License? Yes  
Impact This is a research tool and has enabled the research for all of my publications. 
URL https://github.com/AtmosFOAM
 
Description Climate Modelling Summer School 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Postgraduate students
Results and Impact Lead lecturer at NCAS Climate Modelling Summer school

don't know
Year(s) Of Engagement Activity 2011,2013
URL https://www.ncas.ac.uk/index.php/en/climate-modelling-summer-school
 
Description Climate Modelling summer school 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Postgraduate students
Results and Impact Summer school for post-graduate students and early career scientists to learn about climate modelling. My part in the summer school was to teach programming and the mathematics of climate modelling.
Year(s) Of Engagement Activity 2015
URL https://www.ncas.ac.uk/index.php/en/climate-modelling-summer-school
 
Description Dynamical Core Model Intercomparison Project (DCMIP) 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Series of lectures at the Dynamical Core Model Intercomparison Project workshop in order to train post-graduates and weather forecast model developers in the latest modelling techniques.
Year(s) Of Engagement Activity 2016
URL https://www.earthsystemcog.org/projects/dcmip-2016/
 
Description London Mathematical Society Popular Lecture 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Schools
Results and Impact Numerous questions afterwards from interested sixth formers about mathematical modelling of cliamte change. A UG student also asked about PhD prospects

none
Year(s) Of Engagement Activity 2011
URL http://www.rmets.org/events/london-mathematical-society-popular-lectures-2011-0
 
Description Mathematics Residential Summer school 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Schools
Results and Impact This summer school is run by The Prince's Teaching Institute, a charity set up to re-invigorate secondary school teachers by
getting them out of school and re-inspiring them in their subjects. I gave a talk at this summer school describing climate modelling to school teachers and providing material for climate modelling exercises in classrooms.
Year(s) Of Engagement Activity 2015
URL http://www.princes-ti.org.uk/about/news/residential-summer-school-begins-in-cambridge
 
Description School visit to Caversham primary school 
Form Of Engagement Activity Participation in an open day or visit at my research institution
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Schools
Results and Impact 60 pupils attended a careers morning

none
Year(s) Of Engagement Activity 2014
 
Description Seminar at the National Oceanography Centre 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Professional Practitioners
Results and Impact Seminar at the National Oceanography Centre, Liverpool
on
Atmospheric Modelling on Arbitrary Grids

The advantages of using arbitrary grids for atmospheric modelling were described and the latest outcomes from the Gung Ho Project to design the dynamical core of the next Met Office weather and climate model were summarised

none
Year(s) Of Engagement Activity 2013
URL http://noc.ac.uk/seminars/liverpool/atmospheric-modelling-arbitrary-grids
 
Description Talk at Mathematics for Planet Earth CDT Jamboree 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Postgraduate students
Results and Impact Workshop of the Mathematics for Planet Earth Centre for Doctoral Training
Year(s) Of Engagement Activity 2018
URL https://mpecdt.org/mathematics-of-planet-earth-jamboree-2018/
 
Description Talk at Prospects in Mathematics meeting 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Undergraduate students
Results and Impact Annual meeting for undergraduate students informing them about PhD opportunities in mathematics
Year(s) Of Engagement Activity 2015
URL http://homepages.lboro.ac.uk/~maap/LMS/