Climate models exhibit major radiative biases over the Southern Ocean owing to a poor representation of mixed‐phase clouds. This study uses the remote‐sensing dataset from the Measurements of Aerosols, Radiation and Clouds over the Southern Ocean (MARCUS) campaign to assess the ability of the Weather Research and Forecasting (WRF) model to reproduce frontal clouds off Antarctica. It focuses on the modeling of thin mid‐level supercooled liquid water layers which precipitate ice. The standard version of WRF produces almost fully glaciated clouds and cannot reproduce cloud top turbulence. Our work demonstrates the importance of adapting the ice nucleation parameterization to the pristine austral atmosphere to reproduce the supercooled liquid layers. Once simulated, droplets significantly impact the cloud radiative effect by increasing downwelling longwave fluxes and decreasing downwelling shortwave fluxes at the surface. The net radiative effect is a warming of snow and ice covered surfaces and a cooling of the ocean. Despite improvements in our simulations, the local turbulent circulation related to cloud‐top radiative cooling is not properly reproduced, advocating for the need to develop a parameterization for top‐down convection to capture the turbulence‐microphysics interplay at cloud top.

Snow distribution in complex alpine terrain and its evolution in the future climate is important in a variety of applications including hydropower, avalanche forecasting and freshwater resources. However, it is still challenging to quantitatively forecast precipitation, especially over complex terrain where the interaction between local wind and precipitation fields strongly affects snow distribution at the mountain ridge scale. Therefore, it is essential to retrieve high-resolution information about precipitation processes over complex terrain. Here, we present very-high-resolution Weather Research and Forecasting model (WRF) simulations (COSMO–WRF), which are initialized by 2.2 km resolution Consortium for Small-scale Modeling (COSMO) analysis. To assess the ability of COSMO–WRF to represent spatial snow precipitation patterns, they are validated against operational weather radar measurements. Estimated COSMO–WRF precipitation is generally higher than estimated radar precipitation, most likely due to an overestimation of orographic precipitation enhancement in the model. The high precipitation amounts also lead to a higher spatial variability in the model compared to radar estimates. Overall, an autocorrelation and scale analysis of radar and COSMO–WRF precipitation patterns at a horizontal grid spacing of 450 m show that COSMO–WRF captures the spatial variability normalized by the domain-wide variability in precipitation patterns down to the scale of a few kilometers. However, simulated precipitation patterns systematically show a lower variability on the smallest scales of a few hundred meters compared to radar estimates. A comparison of spatial variability for different model resolutions gives evidence for an improved representation of local precipitation processes at a horizontal resolution of 50 m compared to 450 m. Additionally, differences of precipitation between 2830 m above sea level and the ground indicate that near-surface processes are active in the model.

Accurately simulating snow cover dynamics and the snow–atmosphere coupling is of major importance for topics as wide-ranging as water resources, natural hazards, and climate change impacts with consequences for sea level rise. We present a new modelling framework for atmospheric flow simulations for cryospheric regions called CRYOWRF. CRYOWRF couples the state-of-the-art and widely used atmospheric model WRF (the Weather Research and Forecasting model) with the detailed snow cover model SNOWPACK. CRYOWRF makes it feasible to simulate the dynamics of a large number of snow layers governed by grain-scale prognostic variables with online coupling to the atmosphere for multiscale simulations from the synoptic to the turbulent scales. Additionally, a new blowing snow scheme is introduced in CRYOWRF and is discussed in detail. CRYOWRF's technical design goals and model capabilities are described, and the performance costs are shown to compare favourably with existing land surface schemes. Three case studies showcasing envisaged use cases for CRYOWRF for polar ice sheets and alpine snowpacks are provided to equip potential users with templates for their research. Finally, the future roadmap for CRYOWRF's development and usage is discussed.

High-resolution (< 1 km) atmospheric modeling is increasingly used to study precipitation distributions in complex terrain and cryosphere–atmospheric processes. While this approach has yielded insightful results, studies over annual timescales or at the spatial extents of watersheds remain unrealistic due to the computational costs of running most atmospheric models. In this paper we introduce a high-resolution variant of the Intermediate Complexity Atmospheric Research (ICAR) model, HICAR. We detail the model development that enabled HICAR simulations at the hectometer scale, including changes to the advection scheme and the wind solver. The latter uses near-surface terrain parameters which allow HICAR to simulate complex topographic flow features. These model improvements clearly influence precipitation distributions at the ridge scale (50 m), suggesting that HICAR can approximate processes dependent on particle–flow interactions such as preferential deposition. A 250 m HICAR simulation over most of the Swiss Alps also shows monthly precipitation patterns similar to two different gridded precipitation products which assimilate available observations. Benchmarking runs show that HICAR uses 594 times fewer computational resources than the Weather Research and Forecasting (WRF) atmospheric model. This gain in efficiency makes dynamic downscaling accessible to ecohydrological research, where downscaled data are often required at hectometer resolution for whole basins at seasonal timescales. These results motivate further development of HICAR, including refinement of parameterizations used in the wind solver and coupling of the model with an intermediate-complexity snow model.

In this study, near-surface snow and graupel dynamics from formation to deposition are analyzed using WRF in a large-eddy configuration. The results reveal that a horizontal grid spacing of ≤50 m is required to resolve local orographic precipitation enhancement, leeside flow separation, and thereby preferential deposition. At this resolution, precipitation patterns across mountain ridges show a high temporal and spatial variability. Simulated and observed event-mean snow precipitation across three mountain ridges in the upper Dischma valley (Davos, Switzerland) for two precipitation events show distinct patterns, which are in agreement with theoretical concepts, such as small-scale orographic precipitation enhancement or preferential deposition. We found for our case study that overall terrain–flow–precipitation interactions increase snow accumulation on the leeward side of mountain ridges by approximately 26%–28% with respect to snow accumulation on the windward side of the ridge. Cloud dynamics and mean advection may locally increase precipitation on the leeward side of the ridge by up to about 20% with respect to event-mean precipitation across a mountain ridge. Analogously, near-surface particle–flow interactions, that is, preferential deposition, may locally enhance leeward snow precipitation on the order of 10%. We further found that overall effect and relative importance of terrain–flow–precipitation interactions are strongly dependent on atmospheric humidity and stability. Weak dynamic stability is important for graupel production, which is an essential component of solid winter precipitation. A comparison to smoothed measurements of snow depth change reveals a certain agreement with simulated precipitation across mountain ridges.