Snowfall in the northeastern part of South Korea is the result of complex snowfall mechanisms due to a highly contrasting terrain combined with nearby warm waters and three synoptic pressure patterns. All these factors together create unique combinations, whose disentangling can provide new insights into the microphysics of snow on the planet. This study focuses on the impact of wind flow and topography on the microphysics drawing of 20 snowfall events during the ICE-POP 2018 (International Collaborative Experiment for PyeongChang 2018 Olympic and Paralympic winter games) field campaign in the Gangwon region. The vertical structure of precipitation and size distribution characteristics are investigated with collocated MRR (micro rain radar) and PARSIVEL (particle size velocity) disdrometers installed across the mountain range. The results indicate that wind shear and embedded turbulence were the cause of the riming process dominating the mountainous region. As the strength of these processes weakens from the mountainous region to the coastal region, riming became less significant and gave way to aggregation. This study specifically analyzes the microphysical characteristics under three major synoptic patterns: air–sea interaction, cold low, and warm low. Air–sea interaction pattern is characterized by more frequent snowfall and vertically deeper precipitation systems on the windward side, resulting in significant aggregation in the coastal region, with riming featuring as a primary growth mechanism in both mountainous and coastal regions. The cold-low pattern is characterized by a higher snowfall rate and vertically deep systems in the mountainous region, with the precipitation system becoming shallower in the coastal region and strong turbulence being found in the layer below 2 km in the mountainous upstream region (linked with dominant aggregation). The warm-low pattern features the deepest system: precipitation here is enhanced by the seeder–feeder mechanism with two different precipitation systems divided by the transition zone (easterly below and westerly above). Overall, it is found that strong shear and turbulence in the transition zone is a likely reason for the dominant riming process in the mountainous region, with aggregation being important in both mountainous and coastal regions.

This study uses high-resolution radar and surface observations to investigate the finescale structural evolution of airflow and precipitation over complex terrain in the Pyeongchang area, South Korea. The Taebaek Mountain range (TMR) runs parallel to the northeastern coast of South Korea, with a perpendicular ridge known as the Pyeongchang branch (PCB). The objective of this study was to identify the mechanisms of wintertime precipitation over these topographic features during the passage of a low pressure system (LPS) through the southern Korean Peninsula. The analysis indicates that intense precipitation occurred over the southwestern and northeastern sides of the TMR during stage I but only over the northeastern side during stage II. The prevailing southwesterly winds were dominated by warm advection associated with the LPS over the PCB during stage I. These prevailing southwesterly winds locally enhanced precipitation on the southwestern end of the PCB; multiple influences of mountain waves, airflow convergence, and drifted particles are possible factors for causing precipitation on the northeastern side of the TMR. During stage II, the prevailing winds changed from easterlies to northeasterlies offshore from Gangneung. The easterly winds decelerated and were deflected locally along the mountainous coast, and this blocked zone interacted with the oncoming flow to trigger a precipitation band. Consequently, the northeasterly winds helped push the precipitation band toward the coast, causing heavy precipitation in Gangneung. The observational evidence presented shows that the interaction of temporally changing winds accompanying the movement of an LPS over topography is a critical factor for determining the distribution and intensity of precipitation.

This study evaluates the performance of several cloud microphysics parameterizations in simulating surface precipitation for two snowstorm cases during the International Collaborative Experiment held at the PyeongChang 2018 Olympics and Winter Paralympic Games (ICE-POP 2018) field campaign. We compared four different schemes in the Weather Research and Forecasting (WRF) model, namely the Double-Moment 6-class (WDM6), the WRF Single-Moment 6-class (WSM6), and Thompson and Morrison parameterizations. Both WSM6 and WDM6 overestimated the precipitation amount for the shallow precipitation system because of the substantial amount of cloud ice, mostly generated by the deposition process. The simulated precipitation amount and distribution for the deep precipitation system showed no noticeable differences in the different cloud microphysics parameterizations. However, the simulated hydrometeor type at the surface using WSM6 and WDM6 showed good agreement with observations for all cases. The accuracy of the mean mass-weighted terminal velocity of cloud ice (⁠ $VI¯$⁠) applied in WSM6 and WDM6 is ±20%. The number concentration of cloud ice and the ice microphysics processes are newly retrieved with 1.2 times increased $VI¯$⁠. For the shallow snowstorm, the precipitation amount was reduced by approximately 8% because of the inefficient deposition and its effects on the subsequent ice microphysical processes, such as the accretion of cloud ice by snow and the conversion from cloud ice to snow.

Ground-based radar and radiometer data observed during the 2017–2018 winter season over the Pyeongchang area on the east coast of the Korean Peninsula were used to simultaneously estimate both the cloud liquid water path and snowfall rate for three types of snow clouds: near-surface, shallow, and deep. Surveying all the observed data, it is found that near-surface clouds are the most frequently observed cloud type with an area fraction of over 60 %, while deep clouds contribute the most in snowfall volume with about 50 % of the total. The probability distributions of snowfall rates are clearly different among the three types of clouds, with the vast majority hardly reaching 0.3 mm h⁻¹ (liquid water equivalent snowfall rate) for near-surface, 0.5 mm h⁻¹ for shallow, and 1 mm h⁻¹ for deep clouds. However, the liquid water paths in the three types of clouds all have the substantial probability to reach 500 g m⁻². There is no clear correlation found between snowfall rate and the liquid water path for any of the cloud types. Based on all observed snow profiles, brightness temperatures at Global Precipitation Measurement Microwave Imager (GPM/GMI) channels are simulated, and the ability of a Bayesian algorithm to retrieve snowfall rate is examined using half the profiles as observations and the other half as an a priori database. Under an idealized scenario, i.e., without considering the uncertainties caused by surface emissivity, ice particle size distribution, and particle shape, the study found that the correlation as expressed by R² between the “retrieved” and “observed” snowfall rates is about 0.32, 0.41, and 0.62, respectively, for near-surface, shallow, and deep snow clouds over land surfaces; these numbers basically indicate the upper limits capped by cloud natural variability, to which the retrieval skill of a Bayesian retrieval algorithm can reach. A hypothetical retrieval for the same clouds but over ocean is also studied, and a major improvement in skills is found for near-surface clouds with R² increasing from 0.32 to 0.52, while a smaller improvement is found for shallow and deep clouds. This study provides a general picture of the microphysical characteristics of the different types of snow clouds and points out the associated challenges in retrieving their snowfall rate from passive microwave observations.

Evaluating satellite ability in capturing sudden natural disasters such as heavy snowstorms is a topic of societal interest. This paper presents a rapid qualitative analysis of an intense snowfall in Madrid using data from the Global Precipitation Measurement (GPM) mission, specifically the GPM IMERG (Integrated Multi-satellitE Retrievals for GPM) Late Precipitation L3 Half Hourly 0.1° × 0.1° V06 estimates of precipitation (IMERG-Late), and Sentinel-2 imagery. The main research question addressed is the consistency of ground observations, model outputs and satellite data, a topic of major interest for an appropriate and timely societal response to severe weather episodes. Indeed, the choice of the ‘Late’ product over the IMERG ‘Final’ or other GPM datasets was motivated by the availability of data for near real-time response to the storm. Additionally, the 30-min temporal resolution of the product would in principle allow for a detailed analysis of the dynamic processes involved in the snowstorm. Using several complementary data sources, it is shown that optical remote sensing sensors (Sentinel) add value to existing ground data and that is invaluable for rapid response to severe meteorological events such as Filomena. Regarding the GPM precipitation radar, the sampling of the GPM-core satellite was insufficient to provide the IMERG algorithm with enough quality data to correctly represent the actual sequence of precipitation. Without corrections, the total precipitation differs from observations by a factor of two. The difficulties of retrieving precipitation with radiometers over snow-covered surfaces was a major factor for the mismatch. Thus, the calibrated precipitation product did not fully capture the historic storm, and neither did the IR-based element of the IMERG-Late product, which is a neural network merging of microwave and infrared data. It follows that increased temporal resolution of spaceborne microwave sensors and improved retrieval of precipitation from radiometers are critical in order to provide a complete account of these sorts of extreme, significant, short-duration cases. Otherwise, the high-quality, radar and radiometer data feeding the high temporal resolution algorithms simply slip through the grasp of the ascending and descending orbits, leaving little quality data to be interpolated into successive overpasses.

An accurate classification of the precipitation type is important for forecasters, particularly in the winter season. We explored the capability of three supervised machine learning (ML) methods (decision tree, random forest, and support vector machine) to determine ground precipitation types (no precipitation, rain, mixed, and snow) for winter precipitation. We provided information on the particle characteristics within a radar sampling volume and the environmental condition to the ML model with the simultaneous use of polarimetric radar variables and thermodynamic variables. The ML algorithms were optimized using predictor selection and hyperparameter tuning in order to maximize the computational efficiency and accuracy. The random forest (RF) had the highest skill scores in all precipitation types and outperformed the operational scheme. The spatial distribution of the precipitation type from the RF model showed a good agreement with the surface observation. As a result, RF is recommended for the real-time precipitation type classification due to its easy implementation, computational efficiency, and satisfactory accuracy. In addition to the validation, this study confirmed the strong dependence of precipitation type on wet-bulb temperature and a 1000–850 hPa layer thickness. The results also suggested that the base heights of the radar echo are useful in discriminating non-precipitating area.