Numerical Simulation of a Heavy Rainstorm in Northeast China Caused by the Residual Vortex of Typhoon 1909 (Lekima)

: From 14 to 17 August 2019, a heavy rainstorm occurred in Northeast China due to the combined inﬂuence of the residual vortex of typhoon 1909 (Lekima) and cold air intrusion. Based on the precipitation data of China Meteorological observation stations, surface and upper charts, HMW-8 satellite images, NCEP/NCAR 0.25 ◦ × 0.25 ◦ reanalysis data and WRF4.0 numerical prediction model are used to carry out numerical simulations. According to the weather situation and numerical simulation results, the cause of 1 h severe precipitation is thoroughly studied. Results show that: (1) The strong precipitation area is located on the lee in the northwest of Changbai Mountain. There is a convergence area in the middle of the troposphere, and a strong divergence area in the upper troposphere, with remarkable topographic effect, and the west divergence column inclines on the east convergence column.


Introduction
Northeast China is located in the northernmost part of China. In summer, it is affected by the interaction of westerly zone weather systems, high latitude, subtropical and tropical These NCEP FNL (Final) operational global analysis and forecast data are on 0.25-degree by 0.25-degree grids prepared operationally every 6 h. This product is from the Global Data Assimilation System (GDAS), which continuously collects observational data from the Global Telecommunications System (GTS), and other sources, for many analyses. The FNLs are made with the same model which NCEP uses in the Global Forecast System (GFS), but the FNLs are prepared about an hour or so after the GFS is initialized. The FNLs are delayed so that more observational data can be used. The GFS is run earlier in support of time-critical forecast needs, and uses the FNL from the previous 6 h cycle as part of its initialization. The analyses are available on the surface, at 26 mandatory (and other pressure) levels from 1000 millibars to 10 millibars, in the surface boundary layer and at some sigma layers, the tropopause and a few others. Parameters include surface pressure, sea level pressure, geopotential height, temperature, sea surface temperature, soil values, ice cover, relative humidity, u-and v-winds, vertical motion, vorticity and ozone.

WRF Simulations
The Advanced Research WRF (ARW) modeling system has been in development for the past eighteen years. The current release is Version 4.0, available since June 2018. The ARW is designed to be a flexible, state-of-the-art atmospheric simulation system that is portable and efficient on available parallel computing platforms. The ARW is suitable for use in a broad range of applications across scales ranging from meters to thousands of kilometers. Used in this study is performed by NCAR, which is publicly available at 0.25 • × 0.25 • reanalysis data. One-way nested-grids are employed ( Figure 1). The centers of Domain1, Domain2 are at (45.0 • N, 130.0 • E), with their grid spacings of 9 and 3 km respectively. The grid sizes are 480 × 480 and 721 × 661 correspondingly. There are 33 layers in the vertical direction. The simple ice phase scheme and Kain-Fritsch cumulus parameterization scheme are used [13]. The integrations are performed for 84 h from 0000 UTC 14 to 1200 UTC 17 August 2019.

Typhoon Track and Surface Precipitation
According to the typhoon information provided by the tropical cyclone data center of the China Meteorological Administration (http://www.typhoon.org.cn), and based on the moving path of typhoon 1909 (Lekima) (Figure 2a), which occurred in the Pacific Ocean east of the Philippines on 4 August 2019, at 1745 UTC 9, it landed on Wenling City, Zhejiang Province, and moved northward gradually, then passing through Zhejiang and Jiangsu Provinces. At about 0400 UTC 11, "Lekima" set sail from Lianyungang City, Jiangsu Province. At about 1250 UTC 11, it landed in Qingdao City, Shandong Province, and at that time, the maximum wind force near the center of the typhoon 1909 was 23 m/s. It moved into the Bohai Sea at 2100 UTC 11, and weakened to the tropical depression at 0000 UTC 13, then stopped numbering at 0600 UTC that day.

Typhoon Track and Surface Precipitation
According to the typhoon information provided by the tropical cyclone data center of the China Meteorological Administration (http://www.typhoon.org.cn), and based on the moving path of typhoon 1909 (Lekima) (Figure 2a), which occurred in the Pacific Ocean east of the Philippines on 4 August 2019, at 1745 UTC 9, it landed on Wenling City, Zhejiang Province, and moved northward gradually, then passing through Zhejiang and Jiangsu Provinces. At about 0400 UTC 11, "Lekima" set sail from Lianyungang City, Jiangsu Province. At about 1250 UTC 11, it landed in Qingdao City, Shandong Province, and at that time, the maximum wind force near the center of the typhoon 1909 was 23 m/s. It moved into the Bohai Sea at 2100 UTC 11, and weakened to the tropical depression at 0000 UTC 13, then stopped numbering at 0600 UTC that day.

Typhoon Track and Surface Precipitation
According to the typhoon information provided by the tropical cyclone data center of the China Meteorological Administration (http://www.typhoon.org.cn), and based on the moving path of typhoon 1909 (Lekima) (Figure 2a), which occurred in the Pacific Ocean east of the Philippines on 4 August 2019, at 1745 UTC 9, it landed on Wenling City, Zhejiang Province, and moved northward gradually, then passing through Zhejiang and Jiangsu Provinces. At about 0400 UTC 11, "Lekima" set sail from Lianyungang City, Jiangsu Province. At about 1250 UTC 11, it landed in Qingdao City, Shandong Province, and at that time, the maximum wind force near the center of the typhoon 1909 was 23 m/s. It moved into the Bohai Sea at 2100 UTC 11, and weakened to the tropical depression at 0000 UTC 13, then stopped numbering at 0600 UTC that day. From 14 to 17 August 2019, affected by the northward movement of the residual circulation of typhoon 1909 (Lekima) and the cold air intrusion, a large range of rainstorms was suffered in Northeast China.
According to the synoptic situation, the precipitation process can be divided into two stages. The first stage is from 12 UTC 14 to 12 UTC 16, which is caused by the interaction between the residual vortex and the inverted trough of typhoon 1909, together with the cold air intrusion, and the rain belt extends from southwest to northeast (Figure 2b). The horizontal width of the rainband composed by the rainfall larger than 50 mm is about 100 km, and the length is about 600 km. The rainfall centers are embedded in many blocks, which are typical meso-β-scale characteristics. There are eleven stations whose rainfalls are greater than 100 mm, six stations in Heilongjiang Province, four stations in Jilin Province and one station in Liaoning Province separately. The rainstorm center 154.4 mm is located in Huadian, Jilin Province. The maximum 1 h rainfall is 53.2 mm, which occurs at the period of 19-20 UTC on 14 August. At the second stage, from 12 UTC 16 to 12 UTC 17, the residual vortex center of typhoon 1909 entered Heilongjiang Province. The precipitation was caused by the addition of residual circulation of typhoon 1910(Krosa) and the southeast movement of the upper cold vortex into Liaoning Province (Figure 2c). There are two rainfall centers, one located at the central junction of Liaoning and Jilin Provinces, and the other at the central junction of Jilin and Heilongjiang Provinces, distributed in blocks. The rainstorm center 71.0 mm is located at Liaoyuan, Jilin Province.
From the 1 h intensive precipitation of the automatic weather stations, it can be seen that there are three rainfall peaks during the period from 00 UTC 14 to 12 UTC 17, the first According to the synoptic situation, the precipitation process can be divided into two stages. The first stage is from 12 UTC 14 to 12 UTC 16, which is caused by the interaction between the residual vortex and the inverted trough of typhoon 1909, together with the cold air intrusion, and the rain belt extends from southwest to northeast (Figure 2b). The horizontal width of the rainband composed by the rainfall larger than 50 mm is about 100 km, and the length is about 600 km. The rainfall centers are embedded in many blocks, which are typical meso-β-scale characteristics. There are eleven stations whose rainfalls are greater than 100 mm, six stations in Heilongjiang Province, four stations in Jilin Province and one station in Liaoning Province separately. The rainstorm center 154.4 mm is located in Huadian, Jilin Province. The maximum 1 h rainfall is 53.2 mm, which occurs at the period of 19-20 UTC on 14 August. At the second stage, from 12 UTC 16 to 12 UTC 17, the residual vortex center of typhoon 1909 entered Heilongjiang Province.   Figure 3a,b show the synoptic situation at 500 hPa at 1900 UTC on 14 August 2019, there are two cold vortices and one warm ridge locate in the middle and high latitudes of 110-150 • E, the east (west) cold vortex is 5640 gpm, −18 • C (5600 gpm, −20 • C), which is in the Bering Sea (the Lake Baikal). The west cold vortex was split into two 5640 gpm centers, another small 5640 gpm center is separated to the southwest. As the transverse trough moves to the south and merges into the North China trough, the cold air behind the shallow trough continuously diffuses and moves toward the southeast. For the strong subtropical high, its 5880 gpm isoline reaches the sea of Japan and Sakhalin, which is a huge warm zone. From Figure 4, the HMW-8 satellite image at 2000 UTC on 14 August 2019 shows that double arc-shaped rain belts appeared in Northeast China, and the south branch was over Changbai mountain, and the north one was between Wanda mountain and Sanjiang plain, separately. From the location of the weather systems, the south one is related to the residual vortex of typhoon 1909, and the north one is related to the airflow at the edge of the subtropical high.  Figure 3a,b show the synoptic situation at 500 hPa at 1900 UTC on 14 August 2019, there are two cold vortices and one warm ridge locate in the middle and high latitudes of 110-150° E, the east (west) cold vortex is 5640 gpm, −18 °C (5600 gpm, −20 °C), which is in the Bering Sea (the Lake Baikal). The west cold vortex was split into two 5640 gpm centers, another small 5640 gpm center is separated to the southwest. As the transverse trough moves to the south and merges into the North China trough, the cold air behind the shallow trough continuously diffuses and moves toward the southeast. For the strong subtropical high, its 5880 gpm isoline reaches the sea of Japan and Sakhalin, which is a huge warm zone.   From Figure 4, the HMW-8 satellite image at 2000 UTC on 14 August 2019 shows tha double arc-shaped rain belts appeared in Northeast China, and the south branch was ove Changbai mountain, and the north one was between Wanda mountain and Sanjiang plain separately. From the location of the weather systems, the south one is related to the resid ual vortex of typhoon 1909, and the north one is related to the airflow at the edge of th subtropical high.

Analysis of Water Vapor Condition
At 1412, the residual vortex of typhoon 1909 reaches the Yellow Sea about 150 km southeast of Liaodong Peninsula, and its inverted trough extends northward, and converges with the southwest airflow at the edge of the subtropical high, which transports water vapor and energy to the northeast, and the center of the severe tropical storm 1910 (Krosa) is close to the Kyushu Island of Japan, its outer cloud system is transported from Kyushu, South Korea and the North Yellow Sea to the southeast of Liaoning Province (where the residual vortex located), and water vapor is conveyed to the northeast. After 1412, there are many convergence centers along the inverted trough, and it is located in central Liaoning, Jilin and Heilongjiang Provinces. At 1419, at 850 hPa, Figure 5a shows that the water vapor starts from the center of typhoon 1910, and transports northward along the southerly airflow on the east side of the inverted trough, as does the convergence center of water vapor (Figure 5b). At 1419, the center was -21 × 10 −8 g·hPa −1 ·cm −2 ·s −1 , at the junction of Liaoning and Jilin Provinces, at 1420, the central part −18 × 10 −8 g·hPa −1 ·cm −2 ·s −1 changed from massive to banded. At 1421, the −19 × 10 −8 g·hPa −1 ·cm −2 ·s −1 , which is at the center of Jilin Province, is located at the center of the rainstorm. Therefore, the center of typhoon 1910 transports water vapor to the northwest along the inverted trough, and under the organization of the residual vortex of typhoon 1909, it provides abundant water vapor conditions for the generation of heavy rainfall. Consequently, the first time 1 h intensive precipitation peak of 53.2 mm occurred in Huachuan, central and Eastern Jilin Province.

Analysis of Water Vapor Condition
At 1412, the residual vortex of typhoon 1909 reaches the Yellow Sea about 150 km southeast of Liaodong Peninsula, and its inverted trough extends northward, and converges with the southwest airflow at the edge of the subtropical high, which transports water vapor and energy to the northeast, and the center of the severe tropical storm 1910 (Krosa) is close to the Kyushu Island of Japan, its outer cloud system is transported from Kyushu, South Korea and the North Yellow Sea to the southeast of Liaoning Province (where the residual vortex located), and water vapor is conveyed to the northeast. After 1412, there are many convergence centers along the inverted trough, and it is located in central Liaoning, Jilin and Heilongjiang Provinces. At 1419, at 850 hPa, Figure 5a shows that the water vapor starts from the center of typhoon 1910, and transports northward along the southerly airflow on the east side of the inverted trough, as does the convergence center of water vapor (Figure 5b). At 1419, the center was -21 × 10 −8 g•hPa −1 •cm −2 •s −1 , at the junction of Liaoning and Jilin Provinces, at 1420, the central part −18 × 10 −8 g•hPa −1 •cm −2 •s −1 changed from massive to banded. At 1421, the −19 × 10 −8 g•hPa −1 •cm −2 •s −1 , which is at the center of Jilin Province, is located at the center of the rainstorm. Therefore, the center of typhoon 1910 transports water vapor to the northwest along the inverted trough, and under the organization of the residual vortex of typhoon 1909, it provides abundant water vapor conditions for the generation of heavy rainfall. Consequently, the first time 1 h intensive precipitation peak of 53.2 mm occurred in Huachuan, central and Eastern Jilin Province. At 1600, at 500 hPa, the split cold vortex (5640 gpm) reaches upper Hebei and Liaoning Provinces. At 1610, at 850 hPa, the residual vortex of typhoon 1909 arrives at the Middle East of Jilin Province, the 500 hPa cold vortex(5640 gpm) moves toward southeast, convective systems transport to the east of Liaoning Province, the transportation of water vapor is not sufficient, the west dry air intrudes to eastward, and converges with the wet air delivered from the east. At 1600, at 500 hPa, the split cold vortex (5640 gpm) reaches upper Hebei and Liaoning Provinces. At 1610, at 850 hPa, the residual vortex of typhoon 1909 arrives at the Middle East of Jilin Province, the 500 hPa cold vortex(5640 gpm) moves toward southeast, convective systems transport to the east of Liaoning Province, the transportation of water vapor is not sufficient, the west dry air intrudes to eastward, and converges with the wet air delivered from the east. Figure 6a shows that on the vertical section of 42.5 • N latitudinal relative humidity, there is an 80% vertical quasi-saturated water vapor column around 124 • E that reaches 400 hPa. Meanwhile, because there is a dry area above 700 hPa to the west of 121 • E, and a long and narrow relative dry column around 125 • E is above 500 hPa, the humidity gradient on both sides of the vertical saturated water column is very large. As the relative dry column is placed on the lower saturated water vapor block around 125 • E, it is favorable for the formation of convective instability. The 70% water vapor column between 126-129 • E reaches 300 hPa. Additionally, there is a long and narrow east-west dry zone between 700-850 hPa to the east of 130 • E. Therefore, the convection on the west side of the typhoon-inverted trough is more vigorous than that on the east side (i.e., near the west edge of the subtropical high). As a result, the precipitation on the west side of the typhoon-inverted trough is larger than that on the east side.
Atmosphere 2021, 12, x FOR PEER REVIEW 10 of 20 Figure 6a shows that on the vertical section of 42.5° N latitudinal relative humidity, there is an 80% vertical quasi-saturated water vapor column around 124° E that reaches 400 hPa. Meanwhile, because there is a dry area above 700 hPa to the west of 121° E, and a long and narrow relative dry column around 125° E is above 500 hPa, the humidity gradient on both sides of the vertical saturated water column is very large. As the relative dry column is placed on the lower saturated water vapor block around 125° E, it is favorable for the formation of convective instability. The 70% water vapor column between 126-129° E reaches 300 hPa. Additionally, there is a long and narrow east-west dry zone between 700-850 hPa to the east of 130° E. Therefore, the convection on the west side of the typhoon-inverted trough is more vigorous than that on the east side (i.e., near the west edge of the subtropical high). As a result, the precipitation on the west side of the typhooninverted trough is larger than that on the east side.  Figure 6b shows that, on the vertical profile of the meridional relative humidity along 126.5° E, the 70% water vapor columns near 39-42° N and 46-47° N both reach 300 hPa, the dry region of 300-850 hPa to the north around 46° N aloft tilted from north to south. The 90% saturated water vapor column between 39° N and 42° N reaches 520 hPa, the wet center relative between 42° N and 45° N with its humidity more than 90% is below 700 hPa, and there is a dry column less than 50% that reaches 500 hPa. Therefore, the upper layer is dry and the lower layer is wet, which is conducive to the formation of convective instability in the lower layer there. As a result, it is beneficial to increase precipitation in this area. Figure 7a shows that, along the 42.5° N zonal vertical profile of water vapor flux, water vapor is mainly concentrated in the lower troposphere, below 850 hPa, and it is transported from east to west. Figure 7b shows that, on the vertical profile of water vapor flux along the meridional direction of 126.5° E, the water vapor is transported from south to north, concentrated in the south of 45° N, 35-40° N is larger, and it reaches 500 hPa among 40-45° N. Figure 7c shows that, along 42.5° N the zonal vertical profile of water vapor flux divergence, water vapor is concentrated in the lower troposphere below 600 hPa around 124° E, and below 700 hPa around 127° E. Figure 7d shows that, along the meridional direction of 126.5° E, several water vapor convergence centers were transported from south to north, there are two convergence centers of water vapor at 39-40° N and 42-43° N to the south of 45° N below 700 hPa respectively.  Figure 6b shows that, on the vertical profile of the meridional relative humidity along 126.5 • E, the 70% water vapor columns near 39-42 • N and 46-47 • N both reach 300 hPa, the dry region of 300-850 hPa to the north around 46 • N aloft tilted from north to south. The 90% saturated water vapor column between 39 • N and 42 • N reaches 520 hPa, the wet center relative between 42 • N and 45 • N with its humidity more than 90% is below 700 hPa, and there is a dry column less than 50% that reaches 500 hPa. Therefore, the upper layer is dry and the lower layer is wet, which is conducive to the formation of convective instability in the lower layer there. As a result, it is beneficial to increase precipitation in this area. Figure 7a shows that, along the 42.5 • N zonal vertical profile of water vapor flux, water vapor is mainly concentrated in the lower troposphere, below 850 hPa, and it is transported from east to west. Figure 7b shows that, on the vertical profile of water vapor flux along the meridional direction of 126.5 • E, the water vapor is transported from south to north, concentrated in the south of 45 • N, 35-40 • N is larger, and it reaches 500 hPa among 40-45 • N. Figure 7c shows that, along 42.5 • N the zonal vertical profile of water vapor flux divergence, water vapor is concentrated in the lower troposphere below 600 hPa around 124 • E, and below 700 hPa around 127 • E. Figure 7d shows that, along the meridional direction of 126.5 • E, several water vapor convergence centers were transported from south to north, there are two convergence centers of water vapor at 39-40 • N and 42-43 • N to the south of 45 • N below 700 hPa respectively. Atmosphere 2021, 12, x FOR PEER REVIEW 11 of 20

Dynamic Mechanism
After 1415, there is a belt, consisting of serval vortex convergence centers from Bohai Bay to North Korea, which extends northeast along the southwest jet at the east of the inverted trough on the west edge of the subtropical high. The typhoon-inverted trough and the upper trough converge in the middle of Jilin Province, which is the place of the rainstorm center. The convergence centers and the divergence centers match with each other, which is favorable for the formation of the vertical secondary circulations. Figure 8a shows the maximum divergence center at 850 hPa is −6 × 10 −5 s −1 , it is in the middle of Liaoning and Jilin Province. At 1610, typhoon 1910 reaches the north of Japan, the cold vortex is west of Heilongjiang and Jilin Provinces. There is a convergence center in the middle of the saddle field and near the circulation of the upper cold vortex, the

Dynamic Mechanism
After 1415, there is a belt, consisting of serval vortex convergence centers from Bohai Bay to North Korea, which extends northeast along the southwest jet at the east of the inverted trough on the west edge of the subtropical high. The typhoon-inverted trough and the upper trough converge in the middle of Jilin Province, which is the place of the rainstorm center. The convergence centers and the divergence centers match with each other, which is favorable for the formation of the vertical secondary circulations. Figure 8a shows the maximum divergence center at 850 hPa is −6 × 10 −5 s −1 , it is in the middle of Liaoning and Jilin Province. At 1610, typhoon 1910 reaches the north of Japan, the cold vortex is west of Heilongjiang and Jilin Provinces. There is a convergence center in the middle of the saddle field and near the circulation of the upper cold vortex, the center is in the east of Jilin Province. There is a convergence center at the bottom of the cold vortex, which is −6 × 10 −5 s −1 , with a westerly flow, at the south of Liaoning and the Bohai Sea. From 1706 to 1708, the 500 hPa cold vortex moves to the Heilongjiang and Jilin, the southwest jet in front of the trough begins to strengthen, and there appears to be a strong convergence center −6 × 10 −5 s −1 in Jilin and the east of Heilongjiang Province.  Figure 8d shows, at 1419, between the height of 1500-6000 m, there is a positive vertical wind shear belt located at the west side of the inverted trough, at the junction of eastern Hebei and western Liaoning. The high-value center 28 m/s is located at the junction of Inner Mongolia, Liaoning and Jilin Provinces, and a small center of vertical wind shear 2 m/s at the north edge of the subtropical high is located at the south of Sihote mountain, and the rainstorm center is located on the line between the maximum center and the minimum center of vertical wind shear, which is close to the right side of the large value center. The results show that the right side of the downshear center is favorable for the occurrence of heavy rainfall. Figure 9a shows at 126.5 • E, there is a rising column with the center value of 0.1 m/s at the height of 450-600 hPa, and at the west side, there is an ascending column with the center value of 0.15 m/s at the height of 400-500 hPa, it shows that there is an eastward-moving convection cell at the west of the rainstorm center, and the intensity is strengthened. The results show that the convective system in front of the westerly trough moves eastward and replenishes new energy and water vapor. Figure 9b shows, at 300-600 hPa, 39-43 • N, there is a mass updraft center with the center value of 0.15 m/s, whose horizontal width is around 400 km and vertical height is around 4 km. The middle-level vortex may be the function of the typhoon residual vortex. Figure 9c shows there is a zonal westerly belt with several gale cores embedded in the upper 100-300 hPa, a 20 m/s gale core to the west of 125 • E is caused by the northward cold vortex descending southward, and the gale core between 130-135 • E is 25 m/s, which is caused by the southwest airflow in front of the westerly trough under the blocking of the subtropical high, indicating that there is a westerly jet eastward continuously in the latitudinal direction. Near the ground, east of 126 • E, due to the east side of the typhoon-inverted trough, there is a maximum easterly airflow core between 131-135 • E. The vertical wind shear between the upper and lower layers is between 132-134 • E, which is bigger than 40 m/s. From 120 to 126 • E, the easterly flow gradually decreases, and from 300 hPa to the ground, the westerly flow sinks. Figure 9d shows the jet core with a southerly wind of 30 m/s at 100-350 hPa is caused by the southwest airflow in front of the upper trough. From 39-43 • N, the upper trough is deep, and the surface in front of the westerly trough is a southerly wind. From 43 • N to the north, the lower layer is located at the west of the typhoon-inverted trough, which is the northerly wind. The upper-level in front of the upper trough is the southerly wind. There is a vertical wind shear in the upper and lower layers, and it inclines northward with the height.

Thermal and Instability Characteristics
At 1419, at 850 hPa, the high energy tongue is transported from the center of the typhoon 1910 along the southwest airflow of the typhoon-inverted trough and the edge of the subtropical high to the northeast. Figure 10a shows, at 1419, along 42.5 • N the zonal vertical profile of pseudo equivalent potential temperature and specific humidity, between 120-124 • E, at the height of 400-850 hPa θ se is the low-value area, less than 320 K, which indicates that dry and cold air intrudes from the west. The specific humidity above 700 hPa is less than 1 g/kg. From the west to the east of 120 • E, the dry and cold air is obliquely intruded into the high humidity center between 124-135 • E from the middle troposphere. The maximum value of θ se at the center is greater than 352 k, and the specific humidity at 850 hPa is more than 14 g/kg. At 1421, the maximum vertical velocity between 42-44 • N 124-126 • E is more than 1 m/s. The center area of the low value of θ se is located at 700 hPa, which indicates that the area below is a convective instability area. It can be seen that this area is in the convective instability area before and when the rainstorm occurs. As long as the upward movement exists, the convective systems will be triggered and the convective upward motion will be maintained and strengthened. Figure 10b shows, at 1419, the meridional vertical profile of pseudo equivalent potential temperature and specific humidity along 126.5 • E, and it can be seen that to the south of 44 • N under 700 hPa, θ se is greater than 348 K and the warm moist air with specific humidity greater than 16 g/kg extends 47 • N northward with the southerly wind. There is a high-value area at 41-43 • N. The low value of θ se is between the height of 500-700 hPa to the north of 45 • N, the θ se is 316 K at and the altitude is about 650 hPa. The intensified dry and cold air intrudes obliquely from the middle and upper levels. The front area at the junction of warm and humid air is stable at 41-47 • N. The center of the low value of θ se is located at 3-4 km height. The convective instability is strengthened in the middle and low troposphere, which leads to more intense precipitation. As the unstable state continues, there is still strong precipitation during 1420-1422.

Thermal and Instability Characteristics
At 1419, at 850 hPa, the high energy tongue is transported from the center of the typhoon 1910 along the southwest airflow of the typhoon-inverted trough and the edge of the subtropical high to the northeast. Figure 10a shows, at 1419, along 42.5° N the zonal vertical profile of pseudo equivalent potential temperature and specific humidity, between 120-124° E, at the height of 400-850 hPa θse is the low-value area, less than 320 K, which indicates that dry and cold air intrudes from the west. The specific humidity above 700 hPa is less than 1 g/kg. From the west to the east of 120° E, the dry and cold air is obliquely intruded into the high humidity center between 124-135° E from the middle troposphere. The maximum value of θse at the center is greater than 352 k, and the specific warm and humid air is stable at 41-47° N. The center of the low value of θse is located at 3-4 km height. The convective instability is strengthened in the middle and low troposphere, which leads to more intense precipitation. As the unstable state continues, there is still strong precipitation during 1420-1422.  Figure 11 shows the topographic distribution map of Northeast China, the first stage precipitation of the residual vortex of typhoon 1909 (Lekima) occurred in Huadian of Jilin Province, which is located at the leeward slope of the northwest part of Changbai Mountain, and the topographic effect is significant. However, the second and third stages of precipitation occurred in Gaizhou of Liaoning Province and Huachuan of Heilongjiang Province respectively, which are plain areas with very flat terrain, and the role of topography is relatively small [12].  Figure 11 shows the topographic distribution map of Northeast China, the first stage precipitation of the residual vortex of typhoon 1909 (Lekima) occurred in Huadian of Jilin Province, which is located at the leeward slope of the northwest part of Changbai Mountain, and the topographic effect is significant. However, the second and third stages of precipitation occurred in Gaizhou of Liaoning Province and Huachuan of Heilongjiang Province respectively, which are plain areas with very flat terrain, and the role of topography is relatively small [12]. In order to further study the role of the topography of Changbai Mountain at the first heavy rainfall stage, a meridional section along 126.5° E was made ( Figure 12) [14]. Figure  12a shows that the divergence distribution at 1419 the whole layer above the high terrain is almost a southerly flow before the typhoon landing, and the updraft is not obvious.

Topographic Effect
There is a small convergence area −10 × 10 −5 s −1 in the lower layer of the windward slope on the south side of the high terrain. On the leeward slope of Changbai Mountain, there −5 −1 Figure 11. Topographic height map of Northeast China.
In order to further study the role of the topography of Changbai Mountain at the first heavy rainfall stage, a meridional section along 126.5 • E was made ( Figure 12) [14]. Figure 12a shows that the divergence distribution at 1419 the whole layer above the high terrain is almost a southerly flow before the typhoon landing, and the updraft is not obvious. There is a small convergence area −10 × 10 −5 s −1 in the lower layer of the windward slope on the south side of the high terrain. On the leeward slope of Changbai Mountain, there is a convergence area −10 × 10 −5 s −1 in the middle troposphere, and a large strong divergence area 15 × 10 −5 s −1 in the upper troposphere. Figure 12b shows, at 1420, as typhoon 1909 moves closer to the land, affected by the terrain uplift, the convective systems were uplifted continuously on the windward slope of the south side of the high terrain. The range of convergence area in the middle troposphere of the leeward slope is enlarged. Figure 12c shows, at 1421, the convective systems are still excited continuously on the windward slope of the south side of the high terrain. However, driven by the northward movement of the upper southerly weather system, the upper troposphere divergence column moves northward, and the height of the convergence column below it decreases. As a result, the convective intensity in the rainstorm area is weakened.  Figure 13a shows, at 1419, that there are two convection systems along latitudinal 42.5 • N in the vertical direction between 120-130 • E. The convective system in the 123-125 • E is caused by the westerly trough, with the convergence column in the west and the divergence column in the east, the divergence column tilts above the convergence column to form an inclined convective system. The convective systems are excited by the typhoon-inverted trough in the east of 126-128 • E, and the convective system delivers from the west, when they pass through the windward slope, because of the terrain effect and the vertical upward motion of the western edge of the subtropical high is very strong, the vertical motion is enhanced due to the strengthened convection and the pumping effect of the upper troposphere divergence center. Figure 13b shows, at 1420, with the westerly upper trough moving closer to the mountainous area, due to the topographic uplift, the convective system is uplifted on the windward slope of the west side of the highland. Figure 13c shows, at 1421, due to the topographic uplift, that the convective system is uplifted in the windward slope on the west side of the highland. Due to the weakening of the convective convergence column in the west, and the vertical upward motion of the subtropical high edge decreases, which leads to the weakening of convective activity and the decrease in precipitation in the rainstorm area.

Conclusions
A heavy rainstorm occurred in Northeast China from 14 to 17 August 2019. In this paper, a synoptic method is used to study the weather situation, and the numerical simulation results are studied. Main conclusions are drawn as follows: 1. The precipitation in the first stage is caused by the interaction between the residual vortex of typhoon 1909 (Lekima) and its inverted trough and the upper trough. The rain belt extends from southwest to northeast, and the rainfall center has typical meso-β-scale characteristics. At the second stage, the upper cold vortex falls to the south, the precipitation is caused by the interaction between the cold vortex and the residual circulation of typhoon 1909 and 1910. 2. Before the occurrence of the first 1 h precipitation peak, the upper cold vortex was split into two parts, then the cold air flows southward along with the transverse trough, thus forming an extrusion between the subtropical high. Due to the organization of the residual vorticity of typhoon 1909, a vortex-like belt composed of con-

Conclusions
A heavy rainstorm occurred in Northeast China from 14 to 17 August 2019. In this paper, a synoptic method is used to study the weather situation, and the numerical simulation results are studied. Main conclusions are drawn as follows: 1.
The precipitation in the first stage is caused by the interaction between the residual vortex of typhoon 1909 (Lekima) and its inverted trough and the upper trough. The rain belt extends from southwest to northeast, and the rainfall center has typical meso-β-scale characteristics. At the second stage, the upper cold vortex falls to the south, the precipitation is caused by the interaction between the cold vortex and the residual circulation of typhoon 1909 and 1910.

2.
Before the occurrence of the first 1 h precipitation peak, the upper cold vortex was split into two parts, then the cold air flows southward along with the transverse trough, thus forming an extrusion between the subtropical high. Due to the organization of the residual vorticity of typhoon 1909, a vortex-like belt composed of convergence centers is formed near the rainstorm center. There is a positive vorticity belt in the typhoon-inverted trough. The center of heavy rain is located in front of the positive vorticity center. The secondary circulation forms between the ascending area of the inverted trough, and the sinking center at the inner edge of the subtropical high and the west side of the inverted trough. The water vapor is concentrated in the middle and low troposphere. Starting from the center of typhoon 1910, the water vapor and energy are transported to the northeast along the southwest airflow of the typhoon-inverted trough and the edge of the subtropical high. The upper trough intersects with the residual vortex trough of typhoon 1909, and convective systems are triggered, thus resulting in three times consecutive 1 h heavy rainfalls in the central and eastern part of Jilin Province.

3.
There is a zone tilted downward from northwest to southeast that contains dry and cold air, and it is in a convectively unstable state. Between the vertical direction of 1500-6000 m, the heavy rainfall center is located on the line between the maximum and minimum centers of vertical wind shear, which is close to the right of the maximum center. Upper westerly momentum propagates downward, the center of the rainstorm is located at the west of the downshear, and the upper southerly momentum propagates downward. The system tilts northward with the height, and the rainstorm area is located at the south side of the downshear. The divergent column tilts above the convergence column thus form the inclined convective systems.

4.
The heavy rainfall area is located in the leeward slope of the northwest of Changbai Mountain. The strong pumping effect is formed between the upper strong divergence area and the middle convergence area, which strengthens the vertical velocity, and the terrain effect is significant. The upper cold vortex fell to the south, and the circulation of typhoon 1910 merges into the residual vortex of typhoon 1909, which results in the second 1 h precipitation peak. The fusion of typhoon 1909, 1910 residual vortices and the cold vortices at the lower layer, which results in the third 1 h precipitation peak.
Author Contributions: Y.W. initiated the study and did the analysis of meteorological conditions, numerical simulation, dynamic diagnosis and analysis; wrote the original draft. T.W. contributes to full text discussion, checking and proofreading. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. Due to the amount of data is very large, the data are not publicly available.