A process-based evaluation of the Intermediate Complexity Atmospheric Research Model (ICAR) 1.0.1
Authors: Horak Johannes , Hofer Marlis , Gutmann Ethan , Gohm Alexander , Rotach Mathias W. 
Published: October 2020

Abstract. The verification of models in general is a non-trivial task and can, due to epistemological and practical reasons, never be considered as complete. As a consequence, a model may yield correct results for the wrong reasons, i.e. by a different chain of processes than found in observations. While in the atmospheric sciences guidelines and strategies exist to maximize the chances that models are correct for the right reasons, these are mostly applicable to full-physics models, such as numerical weather prediction models. The Intermediate Complexity Atmospheric Research (ICAR) model is an atmospheric model employing linear mountain wave theory to represent the wind field. In this wind field atmospheric quantities, such as temperature and moisture are advected and a microphysics scheme is applied to represent the formation of clouds and precipitation. This study conducts an in-depth process-based evaluation of ICAR, employing idealized simulations to increase the understanding of the model and develop recommendations to maximize the probability that its results are correct for the right reasons. To contrast the obtained results from the linear-theory-based ICAR model to a full-physics model, idealized simulations with the Weather Research and Forecasting (WRF) model are conducted. The impact of the developed recommendations is then demonstrated with a case study for the South Island of New Zealand. The results of this investigation suggest three modifications to improve different aspects of ICAR simulations. The representation of the wind field within the domain improves when the dry and the moist Brunt-Väisälä frequencies are calculated in accordance to linear mountain wave theory from the unperturbed base state rather than from the time-dependent perturbed atmosphere. Imposing boundary conditions at the upper boundary different to the standard zero gradient boundary condition is shown to reduce errors in the potential temperature and water vapor fields. Furthermore, the results show that there is a lowest possible model top elevation that should not be undercut to avoid influences of the model top on cloud and precipitation processes within the domain. The method to determine the lowest model top elevation is applied to both the idealized simulations as well as the real terrain case study. Notable differences between the ICAR and WRF simulations are observed across all investigated quantities such as the wind field, water vapor and hydrometeor distributions, and the distribution of precipitation. The case study indicates a large shift in the precipitation maximum for the ICAR simulation employing the developed recommendations in contrast to an unmodified version of ICAR. The cause for the shift is found in influences of the model top on cloud formation and precipitation processes in the ICAR simulations. Furthermore, the results show that when model skill is evaluated from statistical metrics based on comparisons to surface observations only, such analysis may not reflect the skill of the model in capturing atmospheric processes such as gravity waves and cloud formation.