HICAR: An intermediate-complexity atmospheric model capable of resolving ridge-scale snow deposition processes over large mountain ranges

Snow models rely on accurate meteorological input data at the spatial scales at which they operate. However, even the highest resolution operational atmospheric models often run at horizontal resolutions at least an order of magnitude coarser than most snow models. Different downscaling techniques can be employed to bridge this scale gap, typically being sorted into either statistical or dynamical techniques. Statistical techniques often rely on temporally invariant spatial patterns or simplistic conceptual relationships, making them computationally cheap but prone to errors. Dynamical downscaling generally offers a counterpoint to this tradeoff: stronger physical basis but more computational demand. Efforts have been made to optimize this tradeoff of dynamic downscaling, reducing computational demand while maintaining physical accuracy of predicted variables as well as the interdependency of downscaled variables such as winds and precipitation. The Intermediate Complexity Atmospheric Research (ICAR) model recently demonstrated an ability to match precipitation patterns from WRF, but with computational costs at least two orders of magnitude lower. While promising, these results from a 4km comparison did not translate to finer spatial resolutions often needed as input to snow models.

Thus, we introduce the High-resolution Intermediate Complexity Atmospheric Research Model (HICAR), a new variant of the ICAR model developed for spatial resolutions as high as 50m. Relative to a traditional atmospheric model like WRF, HICAR maintains the orders-of-magnitude reduction in computational demand which ICAR displayed, while resolving terrain-induced effects on the wind field not seen in ICAR. This is achieved through a novel combination of adjustments to a background wind field based on terrain descriptors with a wind solver. The solver enforces a mass-conservation constraint on the 3D wind field. These modifications successfully mimic dynamic effects such as flow blocking, ridge-crest speed up, and lee-side recirculation to be captured in the resulting wind field. These features are of particular importance for resolving snow deposition patterns, where the snow particles are particularly susceptible to advection by the near-surface flow field. We validate the accuracy of HICAR’s flow features using a wind LiDAR deployed in complex terrain and show a comparison between flow fields from HICAR and WRF at a horizontal resolution of 50 m. These comparisons demonstrate HICAR’s ability to resolve terrain-induced modifications to the flow field which result in increased heterogeneity of ridge-scale snowfall patterns. To this point, preliminary comparisons of snow deposition patterns in complex terrain between the HICAR and WRF models are offered. With this new model, physically-based downscaling of precipitation and other atmospheric variables which preserves their interdependencies is made available for high-resolutions (100m) and large-spatial extents (10,000 km²) which are often demanded by operational land-surface models.