Case Study of a Heavy Snowstorm Associated with an Extratropical Cyclone Featuring a Back-Bent Warm Front Structure

: An extreme snowstorm event that occurred over Heilongjiang and Jilin Provinces on 24–26 November 2013 was related to a cyclone characterized by a back-bent occluded front structure. This study investigates the structure of the back-bent occluded front and snowfall mechanism using multiple observations and NCEP / NCAR 1 ◦ × 1 ◦ reanalysis data in concert with the HYSPLIT model. The main results show that the extreme event was more synoptically governed by the outbreak of the polar vortex and moisture anomaly of the East Sea. The cyclone occurred just ahead of the 500-hPa merged deep trough, and then developed under the e ﬀ ect of the positive vorticity advection ahead of the 500-hPa trough and intense divergence of the upper-level jet. The south-southwest wind strengthened obviously after the merger of the southern and northern branch troughs, which was the main reason behind the cyclone moving northward. The moisture mainly originated from the Sea of Japan, insofar as that dry and cold air in the lower troposphere over the western mainland moistened obviously as it turned southward and passed over the Bohai Sea and the Sea of Japan, supplying abundant moisture for the snowstorm event. The intensity of moisture transport depended on the location and intensity of the cyclone. When the cyclone developed, the dry air continuously intruded into the cyclone’s center, and made a conveyor belt of warm air wrap around it. The dry air gradually changed from descending to ascending motion as it moved ahead of the westerly trough, while the moist air in the northern part of the cyclone moved to the west and south and incorporated into the south of the cyclone center. Warm and moist air was lifted and arrived in the northwestern part of the cyclone after the occluded front’s formation. Frontogenesis within the comma head was enhanced evidently owing to the rotation and deformation. The convergence between the southeast and northeast winds resulted in intense frontogenesis, leading to the enhancement of the front-scale ascent. Strong ascent formed in the comma head of the cyclone, which resulted in intense snowfall.


Introduction
An occluded front is often formed during the development of a stronger cyclone. Operational forecasts show that some extratropical cyclones producing snowstorms over northeastern China in of surface cyclones is a juxtaposition of frontogenesis and instability, where it is conducive to the occurrence of intense snowfall [29].
However, understanding is limited with respect to the evolution, frontal structure and associated snowfall mechanism in a cyclone characterized by a back-bent warm front and wrap-up process, which leads to inaccurate fine-scale forecasting of snow. In operational forecasts, the prediction of snowfall amount is often underestimated, and the locations of intense snowfall are inaccurate, especially in quantitative snow intensity with a big bias.
An intense snowstorm event hit Heilongjiang and Jilin Provinces on 24-26 November 2013, with accumulated snowfall exceeding 20 mm in the east of these provinces and a maximum snowfall of 60 mm in Shuangyashan. The daily snowfall in nine counties over Heilongjiang Province broke the historic record for November since 1961. This extreme event, induced by a rapidly deepening cyclone that formed in northern Jiangsu Province and moved northeastward, brought strong winds and low visibility and disrupted roads and air traffic, seriously affecting people's lives [33,34]. Notably, the development of the cyclone featured a back-bent warm front and warm-core seclusion [2], thus raising the following questions: How did the fronts evolve? What was the structure of the back-bent occluded front? And what effect did it have on snowfall? This paper seeks to answer these questions and provide a reference for forecasting this type of snowstorm.
The rest of this paper is organized as follows: Section 2 describes the data and methods. An overview of the cyclone, snowfall and cloud systems is provided in Section 3. The synoptic-scale background, front evolution and mesoscale characteristics of the event are presented in Section 4. Section 5 explores the snowfall mechanism. And finally, a summary and discussion are presented in Section 6.

Data and Methods
The data used in this paper include conventional surface, upper-air, Jiansanjiang Doppler radar (indicated by the black open circle in Figure 1a), and FY-2E satellite observations; and the Final Analysis (FNL) gridded data of the National Centers for Environmental Prediction (NCEP) with a 6 h interval and 1 • × 1 • horizontal resolution. The FNL data were obtained from the archive (https://rda.ucar.edu/datasets/ds083.2/). FY-2E is Chinese third operational geostationary meteorological satellite, whose Stretched Visible and Infrared Spin Scan Radiometer (S-VISSR) include four Infrared channels (IR1, IR2, IR3, and IR4) and one Visible channel, with a 5 km and 1.25 km resolution at nadir, respectively. The infrared brightness temperature in spectral range between 6.3 and 7.6 µm (IR3) with a 1 h interval was used in this study. The Jiansanjiang radar is a C-band radar, operating in volume coverage pattern 21 (VCP21), scanning nine elevation angles of 0.5 • , 1.5 • , 2.4 • , 3.4 • , 4.3 • , 6.0 • , 9.9 • , 14.6 • , and 19.5 • with a volumetric update time of 6 min during this event. The elevation of the radar antenna is 165 m above mean sea level (MSL) and 90 m above ground level (AGL). The maximum unambiguous ranges for radar reflectivity and radar radial velocity are 300 and 150 km, respectively. The Nyquist velocities are from −24.86 to 24.86 m·s −1 . The Jiansanjing radar's observation time is from 0100 UTC to 0700 UTC during this event.
We used the Air Resources Laboratory Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model, version 4.9 [35,36], developed by the National Oceanic and Atmospheric Administration, to track the moisture transport trajectory and moisture change of the snowstorm event. In this model, the Lagrangian method is applied to calculate the advection and diffusion, and it is commonly used to track the trajectories of air particles, with good results having been achieved [37][38][39][40]. Thus, the moisture transport features and relationship with the cyclone were examined using this model. The central snowfall area (46.5 • -47 • N, 130 • -131.5 • E, indicated by red rectangle in Figure 1a) was selected as the simulation area, and 500 m, 1500 m, and 3000 m were selected as the initial heights in the vertical direction. There are 24 initial points in the simulated track, with 8 initial points at each level. NCEP/GDAS (Global Data Assimilation System) data with a horizontal resolution of 0.5 • × 0.5 • and an interval of 6 h were used as the initial fields and 60 h air parcel back trajectories were calculated Atmosphere 2020, 11, 1272 4 of 21 at the time of initial snowfall (0000 UTC 25 November) and intense snowfall (0600 UTC 25 November), respectively. All air parcel trajectories were recorded every 6 h from the model forecast.

Overview of the Cyclone
The cyclone initially formed over northern Jiangsu Province at 0000 UTC 24 November, with a central pressure of 1012 hPa (Figure 1a). After that, the cyclone moved northeastward and developed intensively, passed over the Yellow Sea and northern Korean Peninsula, moved on over the Sea of Japan and Stanovoy Range, and entered the Sea of Okhotsk by 0000 UTC 27 November. The cyclone developed at a linear deepening rate in the initial 30 h after formation (Figure 1b). Its central pressure was 986 hPa by 0000 UTC 25 November, and it reached its minimum central pressure of 981 hPa at 0600 UTC 25 November. During the subsequent 36 h, the central pressure rose slightly and maintained at lower values. The central pressure of the cyclone decreased by 26 hPa from 0000 UTC 24 to 0000 UTC 25 November at a cyclone-deepening rate of 1.08 Bergeron, meeting the criteria of an explosive cyclone.

Cloud Systems of the Cyclone
The water vapor cloud images ( Figure 2) show the evolution of the cyclone's cloud systems in its development process. At 0600 UTC 24 November, cloud systems had been transformed into comma cloud from the baroclinic leaf as the cyclone initially formed (not shown), and the core of the cyclone was in the area of cold and moist air with cloud top temperatures <−52 °C and 500-hPa relative humidity >90% (Figure 2a). The dark area was located to the southwest of the cyclone, corresponding to the adiabatic descending dry and cold air with 500-hPa relative humidity <70% and cloud top temperatures >−22 °C. By 1800 UTC 24 November, the cyclone had deepened obviously, with more dry air intruding into the cyclone and the core of the cyclone being occupied by the dark area. Common cyclogenesis features were evident in the cloud patterns, including a comma head over Jilin and Heilongjiang Provinces and a dry slot near northern Korean Peninsula (Figure 2b). Compared with Figure 2a, the warm front cloud belts stretched westward and southward around the cyclone center, with the cold front cloud belts moving eastward. As a result, the northern portion of the cold front weakened and the cold and warm fronts in the center of the cyclone did not connect with each other, showing a frontal fracture. The coverage of cloud top temperatures <−52 °C and 500-hPa relative humidity >70% within the warm front cloud belts was expanded to the north and northeast, and the warm and moist air on the northwest side was wrapped around the cyclone center, showing that a back-bent occluded front was formed at this time ( Figure 2b).

Overview of the Cyclone
The cyclone initially formed over northern Jiangsu Province at 0000 UTC 24 November, with a central pressure of 1012 hPa (Figure 1a). After that, the cyclone moved northeastward and developed intensively, passed over the Yellow Sea and northern Korean Peninsula, moved on over the Sea of Japan and Stanovoy Range, and entered the Sea of Okhotsk by 0000 UTC 27 November. The cyclone developed at a linear deepening rate in the initial 30 h after formation (Figure 1b). Its central pressure was 986 hPa by 0000 UTC 25 November, and it reached its minimum central pressure of 981 hPa at 0600 UTC 25 November. During the subsequent 36 h, the central pressure rose slightly and maintained at lower values. The central pressure of the cyclone decreased by 26 hPa from 0000 UTC 24 to 0000 UTC 25 November at a cyclone-deepening rate of 1.08 Bergeron, meeting the criteria of an explosive cyclone.

Cloud Systems of the Cyclone
The water vapor cloud images ( Figure 2) show the evolution of the cyclone's cloud systems in its development process. At 0600 UTC 24 November, cloud systems had been transformed into comma cloud from the baroclinic leaf as the cyclone initially formed (not shown), and the core of the cyclone was in the area of cold and moist air with cloud top temperatures <−52 • C and 500-hPa relative humidity >90% (Figure 2a). The dark area was located to the southwest of the cyclone, corresponding to the adiabatic descending dry and cold air with 500-hPa relative humidity <70% and cloud top temperatures >−22 • C. By 1800 UTC 24 November, the cyclone had deepened obviously, with more dry air intruding into the cyclone and the core of the cyclone being occupied by the dark area. Common cyclogenesis features were evident in the cloud patterns, including a comma head over Jilin and Heilongjiang Provinces and a dry slot near northern Korean Peninsula (Figure 2b). Compared with Figure 2a, the warm front cloud belts stretched westward and southward around the cyclone center, with the cold front cloud belts moving eastward. As a result, the northern portion of the cold front weakened and the cold and warm fronts in the center of the cyclone did not connect with each other, showing a frontal fracture. The coverage of cloud top temperatures <−52 • C and 500-hPa relative humidity >70% within the warm front cloud belts was expanded to the north and northeast, and the Overlaid on the images are the surface low pressure center (L) and the front, the sea level pressure field (black lines; unit: hPa), and the 500-hPa 70% and 90% relative humidity (green lines; unit: %), calculated from NCEP/NCAR reanalysis data with a resolution of 1 • × 1 • . The letter "J" represents the radar station position at Jiansanjiang; the yellow lines AB and CD in (d) denote the location of the cross sections used in Figure 9.
By 0000 UTC 25 November, the dry air had intruded into the core and southwest quadrant of the cyclone, with the cyclone's dry slot evident in coastal areas of eastern Jilin Province and the cloud systems to the north of the occluded front starting to intensify (Figure 2c). The warm front cloud belts stretched southwestward and southward evidently, and wrapped around the cyclone's center (Figure 2c). At 0600 UTC 25 November, the cloud systems to the north of the back-bent occluded front intensified significantly, and the comma head elongated with an east-west orientation and narrowed with a more obvious warp-up feature (Figure 2d). Six hours later, with the movement of the cyclone, the intense cloud systems to the north of the occluded front moved into Russia (not shown), so the snowfall over Heilongjiang Province weakened. In short, the cyclone development was characterized by zonally elongated warm front cloud belts and underwent a wrapping-up process [4,8,9], which is similar to the Shapiro-Keyser cyclone model. That is, the warm front stretched westward and southward as the cyclone developed, and then wrapped around the center of the cyclone.

Overview of Snowfall
Affected by the cyclone, rain first emerged in Liaoning Province at around 0600 UTC 24 November (Figure 3), and snow fell in Heilongjiang, Jilin and northern Liaoning Provinces at 1200 UTC 24 November, with the snowfall event persisting until 0000 UTC 26 November. The distribution of 24 h snowfall shows that the snowstorm occurred mainly in eastern Jilin and southeastern Heilongjiang Provinces, with a maximum snowfall of 26 mm in Antu from 0000 UTC 24 to 0000 UTC 25 November (Figure 3a). The heavy snowfall moved to central and eastern Heilongjiang Province (maximum snowfall amount reaching 57.4 mm in Shuangyashan) from 0000 UTC 25 to 0000 UTC 26 November (Figure 3b), with much of the snow occurring during the period 0000-1200 UTC 25 November and the maximum 6 h snowfall being larger than 20 mm (Figure 3e,f). The observed snowfall in Shuangyashan shows that the strongest snowfall occurred during 0600-0900 UTC 25 November, with a maximum 3 h snowfall of 12.7 mm.

Overview of Snowfall
Affected by the cyclone, rain first emerged in Liaoning Province at around 0600 UTC 24 November (Figure 3), and snow fell in Heilongjiang, Jilin and northern Liaoning Provinces at 1200 UTC 24 November, with the snowfall event persisting until 0000 UTC 26 November. The distribution of 24 h snowfall shows that the snowstorm occurred mainly in eastern Jilin and southeastern Heilongjiang Provinces, with a maximum snowfall of 26 mm in Antu from 0000 UTC 24 to 0000 UTC 25 November (Figure 3a). The heavy snowfall moved to central and eastern Heilongjiang Province (maximum snowfall amount reaching 57.4 mm in Shuangyashan) from 0000 UTC 25 to 0000 UTC 26 November (Figure 3b), with much of the snow occurring during the period 0000-1200 UTC 25 November and the maximum 6 h snowfall being larger than 20 mm (Figure 3e-f). The observed snowfall in Shuangyashan shows that the strongest snowfall occurred during 0600-0900 UTC 25 November, with a maximum 3 h snowfall of 12.7 mm.

Atmospheric Circulation Patterns
The atmospheric circulation patterns were favorable for the extreme snowfall event. Figure 4a,b show the average geopotential height at 300 hPa, the average vertically integrated moisture flux, and its anomalies for 25 November 2013, respectively. The polar vortex was deflected in the eastern hemisphere, with intense negative anomalies of the geopotential height, and a split low center was located over the Novosibirsk Islands (Figure 4a). A blocking high extended from Lake Balkhash to western Siberia, and another high-pressure ridge existed in the Sea of Okhotsk. A deep trough between them stagnated over eastern China. The larger positive and negative anomalies of the geopotential height corresponded with the high-pressure ridges and the low-pressure system, respectively. The split polar vortex was favorable for the cold air to move into the deep trough and affected northeastern China, which is consistent with polar vortices playing an important role in extreme weather events [41][42][43]. There is also abundant moisture over the East Sea during the development of the cyclone (Figure 4b). Associated with the intensified southwesterly flow ahead of the deep trough and the southeasterly flow in the northeast side of the cyclone, above-normal moisture was transported from the Sea of Japan. Enhanced moisture flux led to moisture convergence over the eastern part of the northeastern China, resulting in the snowfall.
Atmosphere 2020, 11, x FOR PEER REVIEW 7 of 21 evolution at Shuangyashan station from 1200 UTC 24 to 0300 UTC 26 Nov 2013 (unit: mm; the interval is 3 h in the daytime and 10 h in the night (indicated by red stars above the bar). The letters "S", "A" and "J' denote the cities of Shuangyashan, Antu and the location of Jiansanjiang radar, respectively.

Atmospheric Circulation Patterns
The atmospheric circulation patterns were favorable for the extreme snowfall event. Figure 4a,b show the average geopotential height at 300 hPa, the average vertically integrated moisture flux, and its anomalies for 25 November 2013, respectively. The polar vortex was deflected in the eastern hemisphere, with intense negative anomalies of the geopotential height, and a split low center was located over the Novosibirsk Islands (Figure 4a). A blocking high extended from Lake Balkhash to western Siberia, and another high-pressure ridge existed in the Sea of Okhotsk. A deep trough between them stagnated over eastern China. The larger positive and negative anomalies of the geopotential height corresponded with the high-pressure ridges and the low-pressure system, respectively. The split polar vortex was favorable for the cold air to move into the deep trough and affected northeastern China, which is consistent with polar vortices playing an important role in extreme weather events [41][42][43]. There is also abundant moisture over the East Sea during the development of the cyclone (Figure 4b). Associated with the intensified southwesterly flow ahead of the deep trough and the southeasterly flow in the northeast side of the cyclone, above-normal moisture was transported from the Sea of Japan. Enhanced moisture flux led to moisture convergence over the eastern part of the northeastern China, resulting in the snowfall.

Synoptic-Scale Background
The heavy snowfall event occurred in an environment that is favorable for cyclogenesis. At 0000 UTC 23 November, a 500-hPa low was present over the central Siberian Plateau, with a longwave trough extending southward from the low center into eastern Mongolia (Figure 5a). A middle branch trough was centered over central Gansu and eastern Qinghai Provinces, with strong cold advection behind its axis. In addition, a strong frontal zone at the base of the two troughs can be observed ( Figure 5a). More to the south, a southern branch trough extended from eastern Sichuan Province to Yunnan Province. A 35 m·s −1 northern jet and a 45 m·s −1 southern jet at 300 hPa were found at the base of the 500-hPa longwave and the middle branch troughs, respectively ( Figure 5b). By 0000 UTC 24 November, the 500-hPa middle branch and the southern branch troughs merged into a deep trough moving eastward, with strong positive vorticity advection ahead of the trough (Figure 5c). An 850-hPa vortex over Shandong Province (not shown) and a surface cyclone over northern Jiangsu Province ( Figure 1a) had occurred just ahead of this 500-hPa merged deep trough. Furthermore, the

Synoptic-Scale Background
The heavy snowfall event occurred in an environment that is favorable for cyclogenesis. At 0000 UTC 23 November, a 500-hPa low was present over the central Siberian Plateau, with a longwave trough extending southward from the low center into eastern Mongolia (Figure 5a). A middle branch trough was centered over central Gansu and eastern Qinghai Provinces, with strong cold advection behind its axis. In addition, a strong frontal zone at the base of the two troughs can be observed ( Figure 5a). More to the south, a southern branch trough extended from eastern Sichuan Province to Yunnan Province. A 35 m·s −1 northern jet and a 45 m·s −1 southern jet at 300 hPa were found at the base of the 500-hPa longwave and the middle branch troughs, respectively ( Figure 5b). By 0000 UTC 24 November, the 500-hPa middle branch and the southern branch troughs merged into a deep trough moving eastward, with strong positive vorticity advection ahead of the trough (Figure 5c). An 850-hPa vortex over Shandong Province (not shown) and a surface cyclone over northern Jiangsu Province ( Figure 1a) had occurred just ahead of this 500-hPa merged deep trough. Furthermore, the northern longwave trough at 500 hPa, with evident cold advection behind its trough axis, lay slightly behind the merged deep trough (Figure 5c). Twelve hours later, the 850-hPa vortex moved over the Yellow Sea, with a west-east oriented shear line extending from the low center and strong southwesterly airflow reaching 33 m·s −1 on the south side of the shear line (Figure 6a), transporting moisture from the Yellow Sea to northeastern China. Accompanying the cold air sweeping to the south and the transport of warm and moist air, the baroclinicity of the atmosphere was intensified, leading to explosive deepening of the surface cyclone ( Figure 2).
Atmosphere 2020, 11, x FOR PEER REVIEW 8 of 21 northern longwave trough at 500 hPa, with evident cold advection behind its trough axis, lay slightly behind the merged deep trough (Figure 5c). Twelve hours later, the 850-hPa vortex moved over the Yellow Sea, with a west-east oriented shear line extending from the low center and strong southwesterly airflow reaching 33 m·s −1 on the south side of the shear line (Figure 6a), transporting moisture from the Yellow Sea to northeastern China. Accompanying the cold air sweeping to the south and the transport of warm and moist air, the baroclinicity of the atmosphere was intensified, leading to explosive deepening of the surface cyclone ( Figure 2).   By 0000 UTC 25 November, the northern longwave trough at 500-hPa and the merged deep trough to the south had merged to form a deep trough with a span of more than 30° of latitude ( Figure  5e), which enlarged the southerly component ahead of the trough and steered the cyclone to move northward. At 300-hPa, the northern and southern jets had merged to form a core over eastern Heilongjiang Province with a wind speed of 50 m·s −1 , with Jilin and Heilongjiang Provinces on the right side of the entrance of the northern jet where strong divergence existed (not shown). The merger of the northern and southern jets occurred at 1200 UTC 24 November (Figure 5d), and the northern jet core intensified obviously after the merger, with eastern Liaoning Province being on the right side of the entrance of the northern upper-level jet. At 850-hPa, the vortex deepened obviously owing to the substantial intensification of the positive vorticity advection ahead of the 500-hPa trough and the enhancing warm advection (Figure 6b). The southeast winds to the north of the shear line increased evidently, with an enlarged strong-wind-speed region, transporting the moisture from the Sea of Japan to Heilongjiang Province (Figure 6b), which supplied abundant water vapor for the snowstorm. The merger of the upper-level jet streaks made the cyclone move northward obviously, which was similar to the approaching of two upper-level jet streaks in the snowstorm event examined by Tian et al. [43] and Gao et al. [18]. The cyclone deepened rapidly under the joint effect of the right By 0000 UTC 25 November, the northern longwave trough at 500-hPa and the merged deep trough to the south had merged to form a deep trough with a span of more than 30 • of latitude (Figure 5e), which enlarged the southerly component ahead of the trough and steered the cyclone to move northward. At 300-hPa, the northern and southern jets had merged to form a core over eastern Heilongjiang Province with a wind speed of 50 m·s −1 , with Jilin and Heilongjiang Provinces on the right side of the entrance of the northern jet where strong divergence existed (not shown). The merger of the northern and southern jets occurred at 1200 UTC 24 November (Figure 5d), and the northern jet core intensified obviously after the merger, with eastern Liaoning Province being on the right side of the entrance of the northern upper-level jet. At 850-hPa, the vortex deepened obviously owing to the substantial intensification of the positive vorticity advection ahead of the 500-hPa trough and the enhancing warm advection (Figure 6b). The southeast winds to the north of the shear line increased evidently, with an enlarged strong-wind-speed region, transporting the moisture from the Sea of Japan to Heilongjiang Province (Figure 6b), which supplied abundant water vapor for the snowstorm. The merger of the upper-level jet streaks made the cyclone move northward obviously, which was similar to the approaching of two upper-level jet streaks in the snowstorm event examined by Tian et al. [43] and Gao et al. [18]. The cyclone deepened rapidly under the joint effect of the right side of the entrance of the 300-hPa jet and ahead of the 500-hPa merged deep trough. By 0600 UTC 25 November, eastern Heilongjiang Province was located just to the left side of the exit of the 300-hPa southern jet, with a maximum divergence of integrated moisture flux of 14 × 10 −5 s −1 (Figure 5f). The intense pumping of the upper-level divergence facilitated the falling of surface pressure and the intensification of ascent, which was then favorable for the deepening of the cyclone and snowfall.
In summary, the extreme event was more synoptically governed by the outbreak of the polar vortex and moisture anomaly of the East Sea. The outburst of the polar vortex intensified the trough and was favorable for cyclogenesis. The snowstorm event occurred in the divergence area of the upper-level jets. Cold air from the west merged with the southern trough, and then merged with the cold air from the polar vortex, leading to the explosive deepening of the cyclone. The merger of the troughs resulted in the increase of the southerly wind component ahead of the trough, which was beneficial to the northward movement of the cyclone. The abundant moisture supply and flow convergence resulted in the snowstorm over northeastern China.

Mesoscale Features of Snowstorm Cloud Systems
The water vapor cloud images ( Figure 2) show that the snowstorm in northeastern China occurred in the intensified cloud belts to the north of the occluded front. Two mesoscale-banded echoes were observed by the Jiansanjiang radar. Thus, the mesoscale features of snowstorm cloud belts are examined from the 6-min resolution reflectivity data. Figure 7a shows two pieces of scattered echo of reflectivity >20 dBZ were first present to the south of the radar scan. After about one hour, the scattered echo was highly organized and two nearby west-east-oriented narrow-banded echoes of reflectivity >30 dBZ appeared in the central part of the echo. The eastern portions of the two narrow-banded echoes had merged by 0548 UTC 25 November, with expansion of reflectivity >35 dBZ (Figure 7c). By 0612 UTC, the two narrow-banded echoes had evolved into a single intense (>30 dBZ) banded echo, with a length of about 200 km, a width of 20-30 km and a maximum reflectivity exceeding 40 dBZ (Figure 7d), showing the formation of a mesoscale banded echo. About 12 min later, the mesoscale banded echo reached the strongest intensity, with a maximum reflectivity of larger than 45 dBZ, and the broadest coverage of reflectivity >40 dBZ during its life span (Figure 7e). The mesoscale banded echo propagated slowly northward and intense reflectivity approached the radar station, but weakened slightly at 0712 UTC (Figure 7f). The life span of the mesoscale echo was about 2 h, contributing more to the total snowfall. The echo gradually weakened, and no observation of the dissipation stage was obtained owing to missing data. This mesoscale band echo was similar to the single band echo defined by Novak et al. (2004), with a linear reflectivity feature 20-200 km in width and >250 km in length, and a minimum reflectivity intensity >30 dBZ that maintained for at least 2 h.
A cross section of the radar reflectivity through the two narrow-banded echoes before the formation of the mesoscale banded echo exhibits two distinct reflectivity cores near the surface, at around 46.3 • N and 46.6 • N, with reflectivity >30 dBZ below 3 km (Figure 7g), corresponding to the two narrow-banded echoes. The northern narrow-banded echo was a little stronger than that in the southern, with the height of the 20 dBZ contour extending to about 6 km and 5 km, respectively, indicating that the two narrow-banded echoes were in different development stages. At 0624 UTC, the mesoscale banded echo formed, and the cross section (Figure 7h) shows that reflectivity had evolved into a single core from multiple cores (Figure 7g), with a reflectivity core of >40 dBZ near the surface and the height of the 20 dBZ contour extending no more than 5 km. The evolution shows that the formation of the mesoscale band echo experienced merger and intensification of radar reflectivity.
In summary, the intense snowfall to the north of the occluded front was closely related to the mesoscale banded echo of reflectivity >30 dBZ in stratiform clouds, which is different from traditional occluded front cloud systems of uniform stratiform cloud precipitation. The snowfall echo was not high, with reflectivity of larger than 30 dBZ below 3 km and a reflectivity core near the surface.

Mechanism of the Snowstorm
Precipitation in the cool season is closely related to moisture, stability, dynamic lift, snowfall efficiency, and temperature [44]. Among these factors, a rich moisture supply and strong ascent are the necessary conditions for snowfall. This section examines the frontal structure and the two key ingredients for snowfall: moisture and frontogenetical forcing.

Frontal Zone Structure and Its Effect on Snowfall
From the above analysis, it can be seen that the intense snowfall was closely related to the backbent occluded front. How, then, did the frontal zone evolve? What was the structure of the occluded front? And what influence did the occluded front have on snowfall?

Horizontal Structure of the Frontal Zone
As shown in Figure 6a, an 850-hPa cyclonic circulation occurred to the right-hand side of a northeast-southwest-oriented and continuous frontal zone, with obvious cold and warm frontal zones near the cyclone center. The corresponding surface cold and warm fronts in the cyclone center connected with each other before the cyclone explosively deepened (Figure 2a). By 0000 UTC 25 November, the vortex had moved northeastward and deepened rapidly, but the temperature gradient of the cold frontal zone near the cyclone center was weakened (Figure 7b). As a result, the initial continuous frontal zone was fractured, with a contracting and narrowing warm frontal zone. It is clear that a surface cold front catching up to a warm front to form an occluded front did not occur instead of forming a fractural front (Figure 2b-c), where the cold and warm fronts in the center of the cyclone did not connect with each other. This stage of the front evolution is called a T-shaped front structure [8], in that the east-moving cold front was nearly perpendicular to the back-bent warm front. The westward (backward) bent warm front twined around the cyclone center to form an occluded front (Figure 2c).

Mechanism of the Snowstorm
Precipitation in the cool season is closely related to moisture, stability, dynamic lift, snowfall efficiency, and temperature [44]. Among these factors, a rich moisture supply and strong ascent are the necessary conditions for snowfall. This section examines the frontal structure and the two key ingredients for snowfall: moisture and frontogenetical forcing.

Frontal Zone Structure and Its Effect on Snowfall
From the above analysis, it can be seen that the intense snowfall was closely related to the back-bent occluded front. How, then, did the frontal zone evolve? What was the structure of the occluded front? And what influence did the occluded front have on snowfall?

Horizontal Structure of the Frontal Zone
As shown in Figure 6a, an 850-hPa cyclonic circulation occurred to the right-hand side of a northeast-southwest-oriented and continuous frontal zone, with obvious cold and warm frontal zones near the cyclone center. The corresponding surface cold and warm fronts in the cyclone center connected with each other before the cyclone explosively deepened (Figure 2a). By 0000 UTC 25 November, the vortex had moved northeastward and deepened rapidly, but the temperature gradient of the cold frontal zone near the cyclone center was weakened (Figure 7b). As a result, the initial continuous frontal zone was fractured, with a contracting and narrowing warm frontal zone. It is clear that a surface cold front catching up to a warm front to form an occluded front did not occur instead of forming a fractural front (Figure 2b,c), where the cold and warm fronts in the center of the cyclone did not connect with each other. This stage of the front evolution is called a T-shaped front structure [8], in that the east-moving cold front was nearly perpendicular to the back-bent warm front. The westward (backward) bent warm front twined around the cyclone center to form an occluded front (Figure 2c).
At 0600 UTC 25 November, the cold frontal zone continued to move eastward and the warm frontal zone further stretched westward, with the increased frontogenesis appearing in the west (back) portion of the cyclone where the northerly airflows, accompanied by the polar cold air, prevailed (Figure 5c). The warm frontal zone contracted and narrowed in the northeast part of the cyclone, whereas the cold front weakened with its northern end extending to about 38 • N (Figure 6c). At 1200 UTC 25 November, the back-bent warm front and polar cold air turned cyclonically and surrounded the low center. A relatively warm center formed around the low center, with a potential temperature of 286 K (Figure 6d), and the coverage of the warm center had enlarged by 1800 UTC 25 November (not shown). Although the potential temperature was higher than that outside of the warm center, it was lower than that in the warm sector, indicating a separation between the warm sector air and the low center. The separation of the warm core occurred in the baroclinic polar air, where the relatively cold and moist air originated from the northeast side of the cyclone rather than the warm sector air in the initial location of the cyclone. Accompanying the warm-core formation, the surface cyclone center was surrounded by the warm front, and the cold front further weakened and moved eastward (not shown). The central pressure increased to 982-hPa, indicating that the cyclone had begun to weaken.
Next, the relative humidity and vertical velocity are analyzed to further discuss the structure of the front in the upper level. Relative humidity, vertical velocity and winds at 600-hPa can also exhibit the evolution of the dark area in the water vapor cloud images ( Figure 2). An east-west oriented dry tongue originating from Gansu, Qinghai, to Shanxi Provinces generated from 0000 UTC 24 November (not shown), and continued to move eastward. At 1800 UTC 24 November, a wedge-shaped dry-air area with relative humidity <60% was located in the southwest portion of the cyclone (Figure 8a), which originated from the northwesterly airflows with descending motion and corresponded to the dark area that extended to the cyclone center in the water vapor cloud images (Figure 2b). The surface cyclone center was located in the area where the relative humidity isolines were dense. The comma head and the cold front cloud areas were in the area of relative humidity >90% that formed by the ascent of warm and moist air, whereas the air to the north of the head of the dry area descended into the cyclone center. At 0000 UTC 25 November (Figure 8b), the 600-hPa-low circulation formed, and ascent intensified within the comma head that was associated with the back-bent occluded front, with an ascent center over southeastern Heilongjiang Province, corresponding to the snowfall. Moreover, a small coverage of ascent presented in the head of the dry-air area. By 0600 UTC 25 November, the coverage of ascent in the head of the dry area expanded largely (Figure 8c), indicating that the dry air was lifted as it arrived ahead of the westerly trough, and was incorporated into the cyclone center through cyclonic rotation (dry air moved westward) at the end of the cyclone's deepening. At the same time, the moist air in the north portion of the cyclone moved to the west and south and was wrapped into the south of the cyclone center through descending motion. Ascent intensified and moved to eastern Heilongjiang Province, with a maximum ascent of −3.2 Pa·s −1 within the comma head, which corresponded to the mesoscale banded echo and resulted in intense snowfall. At 1200 UTC 25 (Figure 8d), the intense ascent moved out of Heilongjiang Province and weakened, leading to a decrease in snowfall.
It is clear that four stages existed during the rapid deepening of the cyclone, including the incipient frontal cyclone, frontal fracture, frontal T-bone and back-bent warm front, and warm-core seclusion, which is similar to the explosive deepening marine cyclone model. Furthermore, the temperatures at other levels showed that the warm core was able to extend to 600-hPa, which is consistent with a cyclone over Mongolia studied by Xiong et al. (2013).

Vertical Structure of the Frontal Zone
To further understand the structure of the back-bent occluded front and its effect on snowfall, zonal (along 46.6° N) and meridional (along 131.5° E) cross sections passing through the snowfall center and the occluded frontal zone are analyzed.
The zonal section along 46.6° N at 0600 UTC 25 November (Figure 9a) shows that deep moist air extended to 300-hPa around the back-bent occluded front and intensive ascent between 700 and 450-hPa, which resulted in the intense cloud area with lower cloud-top temperature within the comma head and the heavy snowfall over eastern Heilongjiang Province. Relative humidity <60% above 600-hPa between 132° E and 136° E corresponded to the descending dry and cold air that rotated cyclonically. A small area of descending motion can be seen below 700-hPa around 135° E, where the distribution of equivalent potential temperature shows that the air was relatively warm, which in turn shows that the head of the dry air had changed from descending to ascending motion. Another deep and moist air layer with ascent was situated to the east of the dry air layer, corresponding to the warm frontal zone in the eastern part of the occluded front.

Vertical Structure of the Frontal Zone
To further understand the structure of the back-bent occluded front and its effect on snowfall, zonal (along 46.6 • N) and meridional (along 131.5 • E) cross sections passing through the snowfall center and the occluded frontal zone are analyzed.
The zonal section along 46.6 • N at 0600 UTC 25 November (Figure 9a) shows that deep moist air extended to 300-hPa around the back-bent occluded front and intensive ascent between 700 and 450-hPa, which resulted in the intense cloud area with lower cloud-top temperature within the comma head and the heavy snowfall over eastern Heilongjiang Province. Relative humidity <60% above 600-hPa between 132 • E and 136 • E corresponded to the descending dry and cold air that rotated cyclonically. A small area of descending motion can be seen below 700-hPa around 135 • E, where the distribution of equivalent potential temperature shows that the air was relatively warm, which in turn shows that the head of the dry air had changed from descending to ascending motion. Another deep and moist air layer with ascent was situated to the east of the dry air layer, corresponding to the warm frontal zone in the eastern part of the occluded front.
The longitudinal section along 131.5 • E at 0600 UTC 25 November (Figure 9b) shows that deep and moist air with ascent was situated to the south of 30 • N, corresponding to cold-front cloud systems. The lines of equivalent potential temperature are downward concave between 43 • N and 45 • N below 500-hPa, indicating that warm air was isolated and corresponded to the occluded front. A deep and moist layer between 44 • N and 54 • N with intensive ascending motion occurred around 45 • N, corresponding to the occluded front, resulting in heavy snowfall over Heilongjiang Province. The longitudinal section along 131.5° E at 0600 UTC 25 November (Figure 9b) shows that deep and moist air with ascent was situated to the south of 30° N, corresponding to cold-front cloud systems. The lines of equivalent potential temperature are downward concave between 43° N and 45° N below 500-hPa, indicating that warm air was isolated and corresponded to the occluded front. A deep and moist layer between 44° N and 54° N with intensive ascending motion occurred around 45° N, corresponding to the occluded front, resulting in heavy snowfall over Heilongjiang Province.
Although intense ascent and a deep moist layer still existed near the occluded front at 1200UTC 25 November, the occluded front frontal zone moved northeastward, with the maximum ascent also moving northward (not shown), resulting in the location of precipitation being further north and the intensity decreasing at the end of the cyclone's development.
In summary, the strong ascent within the comma head was closely related to the snowfall. With the intensification of the 600-hPa-low circulation, the dry air was constantly wrapped into the cyclone center, while the moist air in the north of the cyclone moved to the west and south. The dry and moist air converged within the comma head on the northwest side of the cyclone center and resulted in the enhancement of the frontal zone, which was conducive to the lifting of the moist air to the north of the cyclone, and then the snowfall over eastern Heilongjiang Province.

Moisture Condition
The occurrence of a snowstorm needs an abundant moisture supply. So, what were the features of the moisture transport and relationship with the cyclone for this snowstorm event? Figure 6 illustrates the geopotential height, temperature, and wind at 850-hPa, as well as the divergence of integrated moisture flux from 1000 to 300-hPa. A moisture channel extending from eastern coastal China to eastern Liaoning Province can be clearly seen as the low was just forming, with weak southerly airflows over northeastern China (Figure 6a). A large amount of moisture was transported northward and converged ahead of cold and warm frontal zones, and the moisture flux convergence areas corresponding to cold and warm fronts were connected to each other. Moisture was transported to Liaoning Province from the Yellow Sea and the East China Sea at this time ( Figure  6a). At 0000UTC 25 November, the southeasterly airflow to the north of the shear line strengthened and reached jet intensity owing to the deepening of the vortex (Figure 6b). Northeastern China was located in the northeast quadrant of the vortex and the southeast airflow transported moisture toward northeastern China from the Sea of Japan. Because the northern portion of the cold front was weakened, the moisture flux convergence areas corresponding to cold and warm fronts were disconnected. The northwest-southeast-oriented moisture flux convergence within the warm frontal Although intense ascent and a deep moist layer still existed near the occluded front at 1200 UTC 25 November, the occluded front frontal zone moved northeastward, with the maximum ascent also moving northward (not shown), resulting in the location of precipitation being further north and the intensity decreasing at the end of the cyclone's development.
In summary, the strong ascent within the comma head was closely related to the snowfall. With the intensification of the 600-hPa-low circulation, the dry air was constantly wrapped into the cyclone center, while the moist air in the north of the cyclone moved to the west and south. The dry and moist air converged within the comma head on the northwest side of the cyclone center and resulted in the enhancement of the frontal zone, which was conducive to the lifting of the moist air to the north of the cyclone, and then the snowfall over eastern Heilongjiang Province.

Moisture Condition
The occurrence of a snowstorm needs an abundant moisture supply. So, what were the features of the moisture transport and relationship with the cyclone for this snowstorm event? Figure 6 illustrates the geopotential height, temperature and wind at 850-hPa, as well as the divergence of integrated moisture flux from 1000 to 300-hPa. A moisture channel extending from eastern coastal China to eastern Liaoning Province can be clearly seen as the low was just forming, with weak southerly airflows over northeastern China (Figure 6a). A large amount of moisture was transported northward and converged ahead of cold and warm frontal zones, and the moisture flux convergence areas corresponding to cold and warm fronts were connected to each other. Moisture was transported to Liaoning Province from the Yellow Sea and the East China Sea at this time (Figure 6a). At 0000 UTC 25 November, the southeasterly airflow to the north of the shear line strengthened and reached jet intensity owing to the deepening of the vortex (Figure 6b). Northeastern China was located in the northeast quadrant of the vortex and the southeast airflow transported moisture toward northeastern China from the Sea of Japan. Because the northern portion of the cold front was weakened, the moisture flux convergence areas corresponding to cold and warm fronts were disconnected. The northwest-southeast-oriented moisture flux convergence within the warm frontal zone stretched to the west, with two relatively intense moisture flux convergence areas over eastern Heilongjiang Province and the Sea of Japan (Figure 6b). These two intense moisture flux convergences resulted from confluences between southeast and northeast winds and southwest and southeast winds, respectively.
The backward trajectories show that air parcels originated over western China with slightly different moving tracks at different heights ( Figure 10). The air parcels at 1500 m and 3000 m first moved east or southeast and turned north as they reached the eastern coastal area of the Korean Peninsula, and then rose to arrive in the snowfall area. The air parcels at this two altitudes were initially descending dry air (Figure 10c,e), with specific humidity <2 g·kg −1 , and they moistened as they passing over the Bohai Sea and the Sea of Japan, with specific humidity increasing (Figure 10e). It was shown that the moisture at 1500 m and 3000 m mainly originated from the Bohai Sea and the Sea of Japan, with the southerly flows transporting moisture. The moisture transport of the event is consistent with the work carried out by Sun et al. [19]. They pointed out that the moisture originated primarily from the Yellow Sea and the Sea of Japan. Unlike air parcels at 1500 and 3000 m, the air parcels at an altitude of 500 m first moved northeast to an area further north, and then turned cyclonically as they reached Russia on the east of Heilongjiang Province. Initially, they were moist air with relatively high specific humidity owing to passing over the Bohai Sea (Figure 10a,e).
Atmosphere 2020, 11, x FOR PEER REVIEW 15 of 21 Heilongjiang Province and the Sea of Japan (Figure 6b). These two intense moisture flux convergences resulted from confluences between southeast and northeast winds and southwest and southeast winds, respectively. The backward trajectories show that air parcels originated over western China with slightly different moving tracks at different heights ( Figure 10). The air parcels at 1500 m and 3000 m first moved east or southeast and turned north as they reached the eastern coastal area of the Korean Peninsula, and then rose to arrive in the snowfall area. The air parcels at this two altitudes were initially descending dry air (Figure 10c,e), with specific humidity <2g·kg −1 , and they moistened as they passing over the Bohai Sea and the Sea of Japan, with specific humidity increasing (Figure 10e). It was shown that the moisture at 1500 m and 3000 m mainly originated from the Bohai Sea and the Sea of Japan, with the southerly flows transporting moisture. The moisture transport of the event is consistent with the work carried out by Sun et al. [19]. They pointed out that the moisture originated primarily from the Yellow Sea and the Sea of Japan. Unlike air parcels at 1500 and 3000 m, the air parcels at an altitude of 500 m first moved northeast to an area further north, and then turned cyclonically as they reached Russia on the east of Heilongjiang Province. Initially, they were moist air with relatively high specific humidity owing to passing over the Bohai Sea (Figure 10a,e).  The vertical distribution of the trajectories (Figure 10c) shows that the air parcels at different heights ascended between 0 h and −12 h, and descended before −12 h. Thus, the air parcels between 0 h and −12 h contributed the most to the snowfall, and those before −12 h contributed little. The air parcels at 3000 m before −12 h comprised dry air originating from the mainland. They were obviously moistened, therefore bringing moisture for snowfall when they turned north and passed over the Sea. The air parcels at 500 m and 1500 m were moist during longer period, providing moisture to the snowfall.
At 0600 UTC 25 November, the southeast winds to the north of the shear line intensified, and the intense moisture flux convergence maintained over eastern Heilongjiang Province, disconnecting with that over the Sea of Japan (Figure 6c). Different from 6 h before, the air parcels at 1500 and 3000 m first moved southeast and then turned northwest as they reached the western coast of Japan, traveling longer distances over the Sea of Japan before they rose to arrive in the snowfall area (Figure 10b). As a result, the moisture transport channel became wider and the moisture transport was significantly enhanced, with the specific humidity increased by about 2 g·kg −1 (Figure 10f). Notably, the cyclonic rotation of the air parcels at 500 m was more obvious and some air parcels passed over the Sea of Japan instead of the mainland. The moisture transport and moisture flux convergence intensified significantly over the snowfall area ( Figure 6c). The back trajectories clearly show that the moisture transport was influenced by the cyclone. The air parcels at 1500 m and 3000 m were on the northeast portion of the cyclone, influencing by southeasterly winds, while those at 500 m were on the northwest portion of the cyclone, dominating by northeasterly winds, showing the tilting of the cyclone with height. This can explain why the air parcels at different heights have different moving tracks. The deep cyclone and its baroclinic structure intensified the snowfall.
In summary, the main moisture source over northeastern China originated from the Bohai Sea and the Sea of Japan, and the southeasterly low-level jet continuously transported abundant moisture to the snowstorm area. The moisture transport intensity varied significantly at different times, which mainly depended on the location and intensity of the cyclone.

Analysis of Frontogenesis
The above analysis shows that a back-bent occluded front formed during the rapid deepening of a cyclone, with strong moisture flux convergence and ascent in the northern quadrant of the cyclone. So, what was the main mechanism of ascent? Temperatures and winds at 850-hPa ( Figure 6) show that obvious warm advection was present over northeastern China, which played an important role in generating synoptic-scale ascent. The frontogenetical forcing was also important, so, the frontogenesis function is analyzed below.
Petterssen [45] defined the atmospheric frontogenesis function as the Lagrangian rate of change of the horizontal potential temperature gradient. Frontogenesis can induce an ageostrophic vertical circulation. So, the frontogenesis function is often used to diagnose the vertical motion and it is regarded as a good predictor for snowfall [46,47]. In this paper, the two-dimensional surface frontogenesis function defined by Miller [48] is applied to analyze frontogenesis qualitatively: where u and v denote the horizontal wind component, and θ represents potential temperature. The frontogenesis function is calculated using formula (1). As shown in Figure 11, divergence and descending motion were present to the south of the low center, and obvious frontogenesis was present to the north of the low center, at 1200 UTC 24 November (Figure 11a). An intense frontogenesis associated with the warm shear line was located over the Korean Peninsula. Another weaker southwest-northeast oriented band of frontogenesis passed through Liaoning, Jilin and Heilongjiang Provinces, with maximum frontogenesis in Liaoning Province, where strong ascent existed. This frontogenesis resulted from airflow convergence between southeast and northeast winds in the north portion of the cyclone, which was the vertical extension of the surface warm front (not shown). At 1800 UTC 24 November, the dry air intruded into the cyclone center and a frontal fracture formed, with the low center in the convergence and frontogenesis region. The convergence between southeast and northeast winds was over central Jilin and Heilongjiang Provinces, which resulted in frontogenesis, with a maximum frontogenesis of >3 × 10 −2 K·km −1 ·h −1 and ascent smaller than −0.2 Pa·s −1 in Jilin Province, consistent with the northward-moving snowfall region (Figure 11b).
Atmosphere 2020, 11, x FOR PEER REVIEW 17 of 21 cyclone center and a frontal fracture formed, with the low center in the convergence and frontogenesis region. The convergence between southeast and northeast winds was over central Jilin and Heilongjiang Provinces, which resulted in frontogenesis, with a maximum frontogenesis of >3 × 10 −2 K·km −1 ·h −1 and ascent smaller than −0.2 Pa·s −1 in Jilin Province, consistent with the northwardmoving snowfall region (Figure 11b). At 0000UTC 25 November, the dry air occupied the cyclone center, with the frontogenesis of the low center increasing to 5 × 10 −2 K·km −1 ·h −1 . Accompanying the warm shear line moving northward to 42°-43° N, the southwest and southeast winds on both sides of the warm shear line intensified, resulting in the enhancement of the transport and convergence of the warm and moist air. A large range of frontogenesis formed over the central-eastern area of northeastern China. The convergence between southeast and northeast winds and the wind speed convergence of the southeast winds resulted in another band of frontogenesis extending northward from the low center to eastern Heilongjiang Province. The frontogenesis indicated that the ageostrophic motion was enhanced, and the central-eastern area of northeastern China was covered by convergence and ascent. The strong At 0000 UTC 25 November, the dry air occupied the cyclone center, with the frontogenesis of the low center increasing to 5 × 10 −2 K·km −1 ·h −1 . Accompanying the warm shear line moving northward to 42 • -43 • N, the southwest and southeast winds on both sides of the warm shear line intensified, resulting in the enhancement of the transport and convergence of the warm and moist air. A large range of frontogenesis formed over the central-eastern area of northeastern China. The convergence between southeast and northeast winds and the wind speed convergence of the southeast winds resulted in another band of frontogenesis extending northward from the low center to eastern Heilongjiang Province. The frontogenesis indicated that the ageostrophic motion was enhanced, and the central-eastern area of northeastern China was covered by convergence and ascent. The strong frontogenesis was located at the border regions between Jilin and Heilongjiang Provinces (Figure 11c), and the corresponding snowfall was heavy. At 0600 UTC 25 November, the frontogenesis over central and eastern Heilongjiang Province further enhanced, with an ascent center of <−1.5 Pa·s −1 in eastern Heilongjiang Province (Figure 11d), which was consistent with the intense snowfall. With the bending and wrapping up of the warm front towards the backward of the cyclone center, the back-bent occluded front was formed, with a new convergence area between the southeast and northeast winds within the comma head. It is clear that the intense ascent corresponded to the larger frontogenesis at different times, indicating that frontogenesis played an important role in the ascent.
In summary, with the dry and cold air intruding into the cyclone center and the intensification of the low-level jet, the warm-front frontogenesis obviously enhanced, with the strongest frontogenesis along a convergence area between the southeast and northeast winds; that is, around the surface inverted trough. The increase in frontogenetic secondary circulation and geostrophic deviation led to the intensification of frontal-scale ascent, which was more favorable for the development of snowstorms. After the occluded front formed, the secondary circulation of frontogenesis was strengthened owing to the warm and moist air being lifted to the northwest of the cyclone, where rotation and deformation increased the secondary circulation of frontogenesis, which led to the formation of mesoscale snow belts in the northwest quadrant of the cyclone, resulting in strong snowfall.

Summary and Discussion
A heavy snowfall event occurred in northeastern China on 25-26 November 2013, which was associated with a back-bent occluded front. The synoptic-scale circulation, satellite images, radar echo, frontal structure, and snowfall mechanism were examined. The main conclusions are as follows: (1) The extreme event was more synoptically governed by the outbreak of the polar vortex and moisture anomaly of the East Sea. The outburst of the polar vortex intensified the trough and contributed significantly to cyclogenesis. The snowstorm event occurred in the divergence area of the northern and southern jets. Cold air from the west merged with the southern trough, and then merged with the cold air from the polar vortex to form a deep westerly trough, leading to the explosive deepening of the cyclone ahead of the trough. The increase of the southerly wind component ahead of the trough due to the merger of the troughs steered the cyclone to move northward.
(2) Heavy snowfall was related to the mesoscale-banded echo within the comma head with reflectivity factor of >30 dBZ in stratiform cloud. The heavy snowfall occurred after the occluded front's formation, which is different to the uniform stratiform cloud precipitation in the traditional occluded front cloud system.
(3) The moisture primarily originated from the Bohai Sea and the Sea of Japan, and the dry air from the mainland rotated cyclonically and moistened to provide a supply moisture for snowstorms after passing over the Sea. The moisture mainly come from the middle and lower troposphere. The intensity of the moisture transport was significantly different at different times, mainly depending on the location and intensity of the cyclone.
(4) The dry and relatively warm air was continuously incorporated into the cyclone center, and the warm front stretched towards the rear of the cyclone center. Thus, the dry air gradually changed from descending to ascending motion as it moved ahead of the westerly trough. The moist air in the north of the cyclone moved westward and southward, and descended into the south of the cyclone center. After the occluded front's formation, the warm and moist airflows rose to arrive in the northwest quadrant of the cyclone and the warm-front frontogenesis intensified obviously owing to the rotation and deformation. The strongest frontogenesis was along the convergence area of the southeast and northeast winds, which led to the enhancement of the front-scale intense ascent within the comma head and resulted in strong snowfall.
(5) Due to the shortage of high spatial and temporal resolution observation data, the analysis of the snowfall mechanism was limited. The snowfall event occurred in the warm frontal zone with convective stability (∂θ/∂p < 0). The preliminary study showed that it probably existed moist symmetric and inertial instability, which can explain the observed banded snowfall echoes. However, the resolution of the FNL data is too coarse to resolve the instability. So, the stability mechanism of the event needs to be further revealed based on numerical simulation, which will be discussed in another paper. In addition, we have found many snowstorm cases were associated with the S-K conceptual configuration's cyclones over northeast China, but we can't get a general result through a case study. We will survey many such type cases and hope to obtain a good result in the future. If the result is positive, it may be helpful in pattern recognition in the operational forecast.