Articles | Volume 10, issue 2
Geosci. Model Dev., 10, 977–989, 2017
Geosci. Model Dev., 10, 977–989, 2017

Development and technical paper 27 Feb 2017

Development and technical paper | 27 Feb 2017

On the numerical stability of surface–atmosphere coupling in weather and climate models

Anton Beljaars1, Emanuel Dutra1,2, Gianpaolo Balsamo1, and Florian Lemarié3 Anton Beljaars et al.
  • 1European Centre for Medium-Range Weather Forecasts, Shinfield Park, Reading, RG2 9AX, UK
  • 2Instituto Dom Luiz, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisbon, Portugal
  • 3Inria, Univ. Grenoble Alpes, CNRS, LJK, Grenoble 38000, France

Abstract. Coupling the atmosphere with the underlying surface presents numerical stability challenges in cost-effective model integrations used for operational weather prediction or climate simulations. These are due to the choice of large integration time steps compared to the physical timescale of the problem, aiming at reducing computational burden, and to an explicit flux coupling formulation, often preferred for its simplicity and modularity. Atmospheric models therefore use the surface-layer temperatures (representative of the uppermost soil, snow, ice, water, etc.) at the previous integration time step in all surface–atmosphere heat-flux calculations and prescribe fluxes to be used in the surface model integrations. Although both models may use implicit formulations for the time steps, the explicit flux coupling can still lead to instabilities.

In this study, idealized simulations with a fully coupled implicit system are performed to derive an empirical relation between surface heat flux and surface temperature at the new time level. Such a relation mimics the fully implicit formulation by allowing one to estimate the surface temperature at the new time level without solving the surface heat diffusion problem. It is based on similarity reasoning and applies to any medium with constant heat diffusion and heat capacity parameters. The advantage is that modularity of the code is maintained and that the heat flux can be computed in the atmospheric model in such a way that instabilities in the snow or ice code are avoided. Applicability to snow–ice–soil models with variable density is discussed, and the loss of accuracy turns out to be small. A formal stability analysis confirms that the parametrized implicit-flux coupling is unconditionally stable.

Short summary
Coupling an atmospheric model with snow and sea ice modules presents numerical stability challenges in integrations with long time steps as commonly used for weather prediction and climate simulations. Explicit flux coupling is often applied for simplicity. In this paper a simple method is presented to stabilize the coupling without having to introduce fully implicit coupling. A formal stability analysis confirms that the method is unconditionally stable.