Generalization and application of the flux-conservative thermodynamic equations in the AROME model of the ALADIN system

Abstract. General yet compact equations are presented to express the thermodynamic impact of physical parameterizations in a NWP or climate model. By expressing the equations in a flux-conservative formulation, the conservation of mass and energy by the physics parameterizations is a built-in feature of the system. Moreover, the centralization of all thermodynamic calculations guarantees a consistent thermodynamical treatment of the different processes. The generality of this physics–dynamics interface is illustrated by applying it in the AROME NWP model. The physics–dynamics interface of this model currently makes some approximations, which typically consist of neglecting some terms in the total energy budget, such as the transport of heat by falling precipitation, or the effect of diffusive moisture transport. Although these terms are usually quite small, omitting them from the energy budget breaks the constraint of energy conservation. The presented set of equations provides the opportunity to get rid of these approximations, in order to arrive at a consistent and energy-conservative model. A verification in an operational setting shows that the impact on monthly-averaged, domain-wide meteorological scores is quite neutral. However, under specific circumstances, the supposedly small terms may turn out not to be entirely negligible. A detailed study of a case with heavy precipitation shows that the heat transport by precipitation contributes to the formation of a region of relatively cold air near the surface, the so-called cold pool. Given the importance of this cold pool mechanism in the life cycle of convective events, it is advisable not to neglect phenomena that may enhance it.

[1]  John Thuburn,et al.  Some conservation issues for the dynamical cores of NWP and climate models , 2008, J. Comput. Phys..

[2]  Rick Salmon Poisson-Bracket Approach to the Construction of Energy- and Potential-Enstrophy- Conserving Algorithms for the Shallow-Water Equations , 2004 .

[3]  P. R. Bannon Theoretical Foundations for Models of Moist Convection , 2002 .

[4]  V. Masson,et al.  The AROME-France Convective-Scale Operational Model , 2011 .

[5]  S. Ştefănescu,et al.  An overview of the variational assimilation in the ALADIN/France numerical weather‐prediction system , 2005 .

[6]  René Laprise,et al.  The Euler Equations of Motion with Hydrostatic Pressure as an Independent Variable , 1992 .

[7]  Joël Stein,et al.  Deterministic and fuzzy verification methods for a hierarchy of numerical models , 2009 .

[8]  P. Bénard,et al.  Flux-conservative thermodynamic equations in a mass-weighted framework , 2007 .

[9]  Almut Gassmann,et al.  A global hexagonal C‐grid non‐hydrostatic dynamical core (ICON‐IAP) designed for energetic consistency , 2013 .

[10]  Almut Gassmann,et al.  Towards a consistent numerical compressible non‐hydrostatic model using generalized Hamiltonian tools , 2008 .

[11]  David L. Williamson,et al.  Time-Split versus Process-Split Coupling of Parameterizations and Dynamical Core , 2002 .

[12]  K. Ooyama A Dynamic and Thermodynamic Foundation for Modeling the Moist Atmosphere with Parameterized Microphysics. , 2001 .

[13]  Tetsuya Theodore. Fujita,et al.  PRECIPITATION AND COLD AIR PRODUCTION IN MESOSCALE THUNDERSTORM SYSTEMS , 1959 .

[14]  P. Marquet,et al.  On a general definition of the squared Brunt–Väisälä frequency associated with the specific moist entropy potential temperature , 2014, 1401.2379.

[15]  Gérald Desroziers,et al.  Background‐error covariances for a convective‐scale data‐assimilation system: AROME–France 3D‐Var , 2011 .

[16]  Piet Termonia,et al.  Stability and accuracy of the physics—dynamics coupling in spectral models , 2007 .

[17]  Evelyne Richard,et al.  Numerical simulations of three different MAP IOPs and the associated microphysical processes , 2006 .

[18]  P. Marquet On the computation of moist‐air specific thermal enthalpy , 2014, 1401.3125.

[19]  P. Marquet,et al.  Definition of a moist entropy potential temperature: application to FIRE‐I data flights , 2011, 1401.1097.

[20]  Valerio Lucarini,et al.  ENERGETICS OF CLIMATE MODELS: NET ENERGY BALANCE AND MERIDIONAL ENTHALPY TRANSPORT , 2009, 0911.5689.

[21]  Pierre Bénard,et al.  Dynamical kernel of the Aladin–NH spectral limited‐area model: Revised formulation and sensitivity experiments , 2010 .

[22]  Alina Barbu,et al.  The SURFEXv7.2 land and ocean surface platform for coupled or offline simulation of earth surface variables and fluxes , 2012 .

[23]  U. Wacker,et al.  Continuity equations as expressions for local balances of masses in cloudy air , 2003 .

[24]  George H. Bryan,et al.  A Benchmark Simulation for Moist Nonhydrostatic Numerical Models , 2002 .

[25]  Masaki Satoh Conservative Scheme for a Compressible Nonhydrostatic Model with Moist Processes , 2003 .

[26]  Mariano Hortal,et al.  The development and testing of a new two‐time‐level semi‐Lagrangian scheme (SETTLS) in the ECMWF forecast model , 2002 .

[27]  G. Zängl,et al.  The ICON (ICOsahedral Non‐hydrostatic) modelling framework of DWD and MPI‐M: Description of the non‐hydrostatic dynamical core , 2015 .

[29]  Andrew Staniforth,et al.  Aspects of the dynamical core of a nonhydrostatic, deep-atmosphere, unified weather and climate-prediction model , 2008, J. Comput. Phys..

[30]  Nicholas A. Engerer,et al.  Surface Characteristics of Observed Cold Pools , 2008 .

[31]  Pierre Bénard,et al.  Integration of the fully elastic equations cast in the hydrostatic pressure terrain-following coordinate in the framework of the ARPEGE/Aladin NWP system , 1995 .

[32]  Jimy Dudhia,et al.  Conservative Split-Explicit Time Integration Methods for the Compressible Nonhydrostatic Equations , 2007 .

[33]  Véronique Ducrocq,et al.  The Meso-NH Atmospheric Simulation System. Part I: adiabatic formulation and control simulations , 1997 .

[34]  Katsuyuki V. Ooyama,et al.  A Thermodynamic Foundation for Modeling the Moist Atmosphere , 1990 .

[35]  René Laprise,et al.  Consequences of Using the Splitting Method for Implementing Physical Forcings in a Semi-Implicit Semi-Lagrangian Model , 1998 .

[36]  Jan Polcher,et al.  A Proposed Structure for Coupling Tiled Surfaces with the Planetary Boundary Layer , 2004 .

[37]  Jean-François Geleyn,et al.  An Approach for Convective Parameterization with Memory: Separating Microphysics and Transport in Grid-Scale Equations , 2007 .

[38]  A. Gassmann,et al.  How is local material entropy production represented in a numerical model? , 2015 .

[39]  A. Bott Theoretical considerations on the mass and energy consistent treatment of precipitation in cloudy atmospheres , 2008 .

[40]  Richard H. Johnson,et al.  The Relationship of Surface Pressure Features to the Precipitation and Airflow Structure of an Intense Midlatitude Squall Line , 1988 .

[41]  Nigel Wood,et al.  Analysis of the numerics of physics–dynamics coupling , 2002 .

[42]  Eugenia Kalnay,et al.  Rules for Interchange of Physical Parameterizations , 1989 .

[43]  Piet Termonia,et al.  Predicting Small-Scale, Short-Lived Downbursts: Case Study with the NWP Limited-Area ALARO Model for the Pukkelpop Thunderstorm , 2015 .