Modelling the E ﬀ ects of Aerosol on Mei-Yu Frontal Precipitation and Physical Processes

Featured Application: A conceptual model of aerosol-cloud interaction is ﬁrst built to explain the aerosol e ﬀ ect on Mei-Yu frontal precipitation, with possible application to reducing the signiﬁcant uncertainties in weather and climate prediction. Abstract: The Mei-Yu front is a signiﬁcantly important summer precipitation system in eastern Asia. In recent years, anthropogenic air pollution over the Yangtze-Huaihe region of China has been aggravating continuously. A cloud-resolving model coupled with an idealized frontal model is used to investigate the response of aerosols on the Mei-Yu frontal precipitation. The results indicate that increasing droplet concentrations lead to signiﬁcant precipitation enhancement with the current pollution levels in Mei-Yu frontal system. Under the polluted conditions, the enhanced cold-cloud process is of great importance. Moreover, with the “towing” of active cold-cloud process, cold-cloud and warm-cloud processes developed mutually. These account for the complicated and special microphysical mechanism for aerosol impacts on Mei-Yu frontal system. Furthermore, two types of “microphysical-dynamic positive feedback loop” caused by the interactions of various physical processes and e ﬀ ects (direct dynamic e ﬀ ect, frontogenesis e ﬀ ect, and vapor pump e ﬀ ect) can be found in the Mei-Yu precipitation, which in turn reinforce the microphysical processes. The combined e ﬀ ect is to increase Mei-Yu front precipitation. The interaction of microphysical processes and dynamic processes, and the positive feedback loops they create are the main physical mechanisms behind the signiﬁcant impacts of aerosol on Mei-Yu frontal precipitation. This may also be an important feature of climate change in eastern Asia.


Introduction
Air pollution has become a global environmental problem, because of its great adverse impact on human health and the climate [1]. China, especially eastern China, has been facing serious environmental issues like the particulate matter (PM) pollution [2], due to the rapidly expanding economic and industrial developments, causing a prominent increase in aerosol loading over this area [3][4][5].
The rest of this paper proceeds as follows. Section 2 describes the numerical model and the experimental design. Results of the numerical experiments are presented and discussed in Section 3. Section 4 summarizes the physics mechanisms. Finally, conclusions are given in Section 5.

Numerical Framework
In the study, the Advanced Research WRF (ARW) version 3.5 [59] is utilized. To simulate ideal front system, the mode control equations were improved by using Peng et al. [60] method to bring large scale meridional geostrophic wind forces into the momentum equation. We assume that Mei-Yu front is east-west direction and is forced by the meridional geostrophic wind V g , which is only a function of vertical displacement z. To keep balance, geostrophic potential temperature advection term is introduced in thermodynamic equation, and water vapor advection term is introduced in water vapor equation accordingly. This gives where u, v, and w are zonal velocity (along the frontal boundary), meridional velocity (perpendicular to the frontal boundary), and vertical velocity in the z direction, respectively. V g is meridional geostrophic wind velocity, f is Coriolis parameter, θ is the potential temperature of dry air, ρ d is the density of dry air, D * expresses the dissipation term of the physical property, q v is the mixing ratio of water vapor, q m = q c , q r , q i , q s , q g represent the mixing ratio of cloud water, rainwater, cloud ice, snow, and graupel, respectively, ρ = ρ d (1 + q v + q c + q r + q i + q s + q g ) is the total mass density, S θ represents the heating/cooling effect associated with cloud microphysical processes, S q v is the source/sink term for water vapor, and S q m is the source/sink term for hydrometeors. The superscript' represents deviations from the dry, static equilibrium basic-state p(z) and ρ(z). The model domain is L x = 1000 km, L y = 2000 km, L z = 20 km in the f plane, using f = 10 −4 s −1 . The lateral boundaries are periodic in x (west-east) and open in y (south-north). The lower boundary is flat without terrain. To inhibit gravity wave reflection from the upper boundary, Rayleigh damping is applied to the vertical velocity in the uppermost 5 km of the domain [61,62]. Horizontal grid resolution is ∆x = ∆y = 5 km. The number of the vertical levels was 81, with a nearly constant vertical spacing of ∆z ≈ 250 m. Time step is ∆t = 20 s. Morrison double-moment cloud microphysics scheme [63] is used and depicted in Section 2.2. Cumulous cloud parameterization is not considered. The radiative effect of clouds is considered, but the direct radiative effect of aerosol is not. The simulations are integrated for 48 h and model output data is available every 5 min.

Microphysical Scheme
The cloud physics scheme used in this paper follows the two-moment bulk scheme proposed by Morrison et al. [63], but further simplifies the cloud-aerosol interactions by setting cloud droplet number concentrations to reflect changes in aerosol concentration, instead of considering the droplet activation. The scheme takes into account primary production, condensation/evaporation, autoconversion, self-collection, collection between hydrometeor species, melting/freezing, and ice multiplication to describe the warm-cloud and cold-cloud processes in natural rain. However, the aerosol effects on ice nuclei are not included due to the less polluted environment in high levels than near surface. In addition, the saturation adjustment approach is applied to calculate cloud water condensation and evaporation (i.e., excess vapor above water saturation is condensed, or cloud water is evaporated to restore water saturated), with a fixed cloud droplet number concentration. Specifically, if the supersaturation at the advanced time (δ') satisfies δ' > 0, then the condensation rate for cloud droplet (PCD) is computed by Equation (8) as described in the two-moment bulk scheme proposed by Morrison et al. [63]. Instead, if δ' < 0, then the evaporation rate of cloud water is computed by Equation (9), where q sw is the mixing ratio at water saturation, c p is the specific heat of air at constant pressure and L v is the latent heat of evaporation. The output of the scheme contains the number concentrations and mixing ratios of four hydrometeor species (rain, cloud ice, snow, and graupel).
Specific parameter values used are the same as Peng et al. [60]. The distributions are shown in Figure 1.
The basic-state sounding comes from an average measured distribution over the region [110°-125° E, 28°-32° N], taken from a typical Mei-Yu front weather conditions on 12 July 2010 at 00 UTC. Its convective available potential energy (CAPE) was approximately 1000 J kg −1 . Specific parameter values used are the same as Peng et al. [60]. The distributions are shown in Figure 1.  (b) Vertical profiles of the vertical shear of the meridional geostrophic wind (∂V g /∂z, solid line, Units: 10 −3 s −1 ) and the large-scale zonal gradient of the potential temperature (∂θ g /∂x, dotted line, Units: K/100 km). (c) Velocity profile along the frontal boundary (u 0 , isolines, 2.0 m s −1 intervals) and (d) potential temperature profile (θ 0 , isolines, 5 K intervals) plotted with water vapor mixing ratio (q v0 , shaded, 2 g kg −1 intervals).

Experimental Design
With increasing aerosol concentration, the concentration of CCN and hence, cloud droplets, also increase dramatically [64]. Therefore, the concentration of cloud particles can reflect the concentration of aerosol. Given a fixed content of cloud water, the more cloud droplets there are, the smaller the size of droplets could be. For simplicity, cloud-aerosol interactions are not explicitly included in the two-moment bulk scheme described by Morrison et al. [63]; instead, different cloud droplet concentrations are prescribed to imitate the effects of different aerosol concentrations in clean and polluted conditions following some scholars [62,65,66]. For instance, Morrison [61] set the initial cloud droplet concentrations to 250 (750) cm −3 to represent the moderately (highly) polluted conditions in understanding the aerosol effect on deep convections. Similarly, this study establishes nine experiments with the initial cloud droplet number concentrations (N c ) set to 50, 100, 250, 500, 1000, 1500, 2500, 3500, and 5000 (cm −3 ) to reveal the effect of aerosols on Mei-Yu system clouds and precipitation. Note that the highest N c is set to 5000 cm −3 to represent the extreme continental aerosol conditions.

Effects on the Intensity and Distribution of Precipitation
The effect of over-seeding is presented in Weather Modification [67]. Similar results can also be found in our research. Figure 2 shows time series of domain-averaged surface precipitation rate (DPRE). It can be seen in Figure 2a that precipitation begins after t = 14 h in all nine cases. Interestingly, a slight variation of DPRE can be found with the increase of N c before t = 18 h, but after then, such differences reach amplitude peak at t = 21 h. Note that the changes of DPRE are nonmonotonic with the increase of N c . More specifically, DPRE increases with N c ranging from 50 cm −3 to 3500 cm −3 , but decreases slightly when N c ranges from 3500 cm −3 to 5000 cm −3 . In order to understand model behavior over a wide range of conditions, we carried out 9 experiments. However, 5000 cm −3 is more than in-situ observations in recent, which might occur in the future with much higher urban air pollution. For simplicity, here we focus on two typical conditions, one is typical continent condition case (C-case) and the other is polluted continent condition case (P-case) in Mei-Yu front. Previous observations indicated that droplet concentration is about 100~300 cm −3 over continents [68]. Nevertheless, Andreae et al. [69] demonstrated that droplet concentration reaches 2400 cm −3 in polluted clouds over the Amazon by observation. Cloud droplet concentration also reaches 1900 cm −3 in northern China [70]. Therefore, droplet concentration herein is set to 250 cm −3 for C-case, while 2500 cm −3 for P-case due to the severe air pollution in Yangtze River basin [54,71]. A similar set was also adopted to investigate the aerosol effects on squall line in Dong et al. [72]. Figure 2b shows time evolution of DPRE in C-case and P-case. Note that precipitation of both cases starts after t = 14 h, but the onset of surface precipitation in C-case (P-case) is at 14:30 (14:40). Precipitation in C-case starts 10 min earlier than P-case ( Figure 2c). The effect of aerosol varies in different stages of the precipitation. Before t = 18 h, Precipitation rate is slightly smaller in P-case than C-case, but larger after that, and this phenomenon continues to t = 26 h. The greatest difference is about 22% between P-case and C-case in 21-22 h. Domain-mean accumulated precipitation in the two cases is also quite different after 18:00 (see Figure 2d). In order to understand model behavior over a wide range of conditions, we carried out 9 experiments. However, 5000 cm −3 is more than in-situ observations in recent, which might occur in the future with much higher urban air pollution. For simplicity, here we focus on two typical conditions, one is typical continent condition case (C-case) and the other is polluted continent condition case (P-case) in Mei-Yu front. Previous observations indicated that droplet concentration is about 100~300 cm −3 over continents [68]. Nevertheless, Andreae et al. [69] demonstrated that droplet concentration reaches 2400 cm −3 in polluted clouds over the Amazon by observation. Cloud droplet concentration also reaches 1900 cm −3 in northern China [70]. Therefore, droplet concentration herein is set to 250 cm −3 for C-case, while 2500 cm −3 for P-case due to the severe air pollution in Yangtze River basin [54,71]. A similar set was also adopted to investigate the aerosol effects on squall line in Dong et al. [72]. Figure 2b shows time evolution of DPRE in C-case and P-case. Note that precipitation of both cases starts after t = 14 h, but the onset of surface precipitation in C-case (P-case) is at 14:30 (14:40). Precipitation in C-case starts 10 min earlier than P-case ( Figure 2c). The effect of aerosol varies in different stages of the precipitation. Before t = 18 h, Precipitation rate is slightly smaller in P-case than C-case, but larger after that, and this phenomenon continues to t = 26 h. The greatest difference is about 22% between P-case and C-case in 21-22 h. Domain-mean accumulated precipitation in the two cases is also quite different after 18:00 (see Figure 2d). Figure 3 compares the spatial distribution of precipitation in C-case and P-case at three different times. At t = 15 h, the two cases both show a line of precipitation parallel to the front. By t = 18 h, many meso-scale and micro-scale rain cells appear on the frontal rainbands, showing the "core-gape" structure [73] of Mei-Yu front precipitation, but the distribution differences of rainbands are not significant between them. At t = 21 h the differences are clear, the frontal rainband of C-case is more spread-out, while the P-case rainband is closer to the front.
Appl. Sci. 2019, 9, x 8 of 34 Figure 3 compares the spatial distribution of precipitation in C-case and P-case at three different times. At t = 15 h, the two cases both show a line of precipitation parallel to the front. By t = 18 h, many meso-scale and micro-scale rain cells appear on the frontal rainbands, showing the "core-gape" structure [73] of Mei-Yu front precipitation, but the distribution differences of rainbands are not significant between them. At t = 21 h the differences are clear, the frontal rainband of C-case is more spread-out, while the P-case rainband is closer to the front.

Cloud Microphysical Quantities
At t = 14 h, driven by the Mei-Yu front, cloud water begins to form at a height of 2.5-3 km. The content of cloud water in the two cases is the same, while the number of cloud droplets is higher in P-case, resulting in smaller drop size. At t = 14 h 30 min, surface precipitation begins to appear in Ccase. Figure 4 shows the average vertical distribution profile of hydrometeors in the two cases, where Qc, Qr, Qi, Qs, and Qg represent the mixing ratios of cloud water, rain water, ice crystals, snow, and graupel, respectively. It is notable that, at this time, the two cases contain only cloud water and rain water. Qc of P-case is greater than that of C-case, but Qr is less obvious (Figure 4a,b). This indicates that cloud drop-size of P-case is smaller, and less efficient at forming stage of raindrops. Figure 2c shows surface precipitation starts at t = 14 h 40 min in P-case, 10 min later than that of C-case.

Cloud Microphysical Quantities
At t = 14 h, driven by the Mei-Yu front, cloud water begins to form at a height of 2.5-3 km. The content of cloud water in the two cases is the same, while the number of cloud droplets is higher in P-case, resulting in smaller drop size. At t = 14 h 30 min, surface precipitation begins to appear in C-case. Figure 4 shows the average vertical distribution profile of hydrometeors in the two cases, where Q c , Q r , Q i , Q s , and Q g represent the mixing ratios of cloud water, rain water, ice crystals, snow, and graupel, respectively. It is notable that, at this time, the two cases contain only cloud water and rain water. Q c of P-case is greater than that of C-case, but Q r is less obvious (Figure 4a,b). This indicates that cloud drop-size of P-case is smaller, and less efficient at forming stage of raindrops. Figure 2c shows surface precipitation starts at t = 14 h 40 min in P-case, 10 min later than that of C-case.   The symbols Q c , Q r , Q i , Q s , and Q g represent the mixing ratios of cloud water, rain water, ice crystals, snow, and graupel, respectively.
In the early stage of precipitation (t = 15 h), cloud top extends to a height of 8 km (0 • C layer located at 3.75 km height) and ice phase process starts. Even though there is smaller Q r in P-case than in C-case, P-case has larger Q s and Q g (see Figure 5a,b). Overall, P-case has less warm-cloud hydrometeors and more cold-cloud hydrometeors than C-case ( Figure 5c).    The symbols Q c , Q r , Q i , Q s , and Q g represent the mixing ratios of cloud water, rain water, ice crystals, snow, and graupel, respectively.
More cloud droplets with smaller drop-size in P-case suppress the conversion of cloud droplets to raindrops, which inhibits precipitation in the early stage, and, meanwhile, leads to more cloud water detained and foster more snow and graupel formation at upper levels. These two aspects are the main differences in cloud physical characteristics of two cases at the onset of precipitation.
In the developing stage of precipitation (t = 18 h), rain water content exceeds cloud water below a height of 3 km. Snow and graupel content increase dramatically in both cases above a height of 3 km (Figure 6a,b). This shows that cold rain processes are greatly invigorated in the developing stage. Comparing the two cases, the content of rain water is quite similar, but P-case has more hydrometeors of all other species, and larger content of total hydrometeor. The biggest difference is at a height of 6 km (Figure 6c). Note that, in addition to increased ice-phase hydrometeors, the content of cloud water in P-case is much greater now than that at t = 15 h, implying that the condensation process (warm-cloud process) has also been invigorated.
Appl. Sci. 2019, 9, x 10 of 34 In the early stage of precipitation (t = 15 h), cloud top extends to a height of 8 km (0 °C layer located at 3.75 km height) and ice phase process starts. Even though there is smaller Qr in P-case than in C-case, P-case has larger Qs and Qg (see Figure 5a,b). Overall, P-case has less warm-cloud hydrometeors and more cold-cloud hydrometeors than C-case ( Figure 5c).
More cloud droplets with smaller drop-size in P-case suppress the conversion of cloud droplets to raindrops, which inhibits precipitation in the early stage, and, meanwhile, leads to more cloud water detained and foster more snow and graupel formation at upper levels. These two aspects are the main differences in cloud physical characteristics of two cases at the onset of precipitation.
In the developing stage of precipitation (t = 18 h), rain water content exceeds cloud water below a height of 3 km. Snow and graupel content increase dramatically in both cases above a height of 3 km (Figure 6a,b). This shows that cold rain processes are greatly invigorated in the developing stage. Comparing the two cases, the content of rain water is quite similar, but P-case has more hydrometeors of all other species, and larger content of total hydrometeor. The biggest difference is at a height of 6 km ( Figure 6c). Note that, in addition to increased ice-phase hydrometeors, the content of cloud water in P-case is much greater now than that at t = 15 h, implying that the condensation process (warmcloud process) has also been invigorated. In the maximum precipitation period (t = 21 h), ice-phase hydrometeor content of both cases grows dramatically. The two most prominent types of hydrometeors are still snow and graupel In the maximum precipitation period (t = 21 h), ice-phase hydrometeor content of both cases grows dramatically. The two most prominent types of hydrometeors are still snow and graupel ( Figure 7). Figure 8 shows the average zonal vertical distribution of hydrometeors content at t = 21 h. Comparing the two cases, the difference in hydrometeor content is greatest in the mixed-phase layer (0 • C to −20 • C). The difference is mainly in the graupel, snow, and cloud water content. In P-case, maximum hydrometeor mixing ratios are 2 g kg −1 for graupel; 1.5 for snow crystals; and 0.5 for cloud water, which are 2 times, 1.5 times, and 1.7 times of those in C-case, respectively. The content of surface rain water in P-case is roughly double than C-case (Figure 8c,d). Correspondingly, the greatest difference of DPRE (about 22%) appears at this time between the two cases ( Figure 2c).
( Figure 7). Figure 8 shows the average zonal vertical distribution of hydrometeors content at t = 21 h. Comparing the two cases, the difference in hydrometeor content is greatest in the mixed-phase layer (0 °C to −20 °C). The difference is mainly in the graupel, snow, and cloud water content. In P-case, maximum hydrometeor mixing ratios are 2 g kg −1 for graupel; 1.5 for snow crystals; and 0.5 for cloud water, which are 2 times, 1.5 times, and 1.7 times of those in C-case, respectively. The content of surface rain water in P-case is roughly double than C-case (Figure 8c,d). Correspondingly, the greatest difference of DPRE (about 22%) appears at this time between the two cases ( Figure 2c). Overall, the main differences appear in the mixed-phase layer at a height of 6 km, where snow and graupel grow rapidly [74,75]. This suggests that aerosol concentration influences the amounts of cloud water above the 0 °C isotherm, and then influences the production and amounts of graupel and snow in mixed phase layer, thereby affecting the formation and amounts of precipitation in Mei-Yu system. Overall, the main differences appear in the mixed-phase layer at a height of 6 km, where snow and graupel grow rapidly [74,75]. This suggests that aerosol concentration influences the amounts of cloud water above the 0 • C isotherm, and then influences the production and amounts of graupel and snow in mixed phase layer, thereby affecting the formation and amounts of precipitation in Mei-Yu system.

Cloud Microphysical Processes
The differences in precipitation and hydrometeor content caused by aerosols are closely related to microphysical processes. In order to better understand the microphysical differences, we analyze the various sources (sinks) of precipitable hydrometers in both cases. Tables 1-3 give the sources (sinks) of the main precipitable hydrometeors (rain water, graupel, and snow), as well as their

Cloud Microphysical Processes
The differences in precipitation and hydrometeor content caused by aerosols are closely related to microphysical processes. In order to better understand the microphysical differences, we analyze the various sources (sinks) of precipitable hydrometers in both cases. Tables 1-3 give the sources (sinks) of the main precipitable hydrometeors (rain water, graupel, and snow), as well as their associated conversion rates during the onset of precipitation period and the peak precipitation period.

Microphysical Processes of Rain Formation
During the early stage of precipitation (t = 15 h), the formation of rain water mainly comes from collision-coalescence between raindrops and cloud droplets. The second significant sources of rain water at this time include the melting of snow and graupel (see Table 1). Noted that auto-conversion of cloud droplets into raindrops is very important in the initial formation of raindrops, even though it occurs with a very low efficiency. The concentration of cloud droplets in C-case is lower. Accordingly, their size is larger, and they are more efficient at auto-converting into raindrops than in P-case. In fact, the maximum rate of auto-conversion of cloud droplets into raindrops is 25 times in C-case than in P-case (see Table 1). Therefore, raindrops form earlier; then grow larger by collision coalescence and fall faster in C-case, leading to the rainfall to surface earlier (Figure 2c,d). By contrast, the formation of raindrops in P-case is inhibited in the initial stage of the precipitation, due to the initial concentration of cloud droplets is relatively higher and cloud drop-size is smaller. Furthermore, low conversion rate means more retention of cloud water, which will be transported to higher levels by frontal updraft, leading to enhancement of cold-cloud process.
During the maximum precipitation period (t = 21 h), conversion rates of both cases have increased dramatically, and the largest conversion rate is the melting of graupel. At this time, except for the auto-conversion of cloud water into rain, all the other rainfall source items (conversion rates) in P-case is greater than in C-case. The ratios of the maximum conversion rate are as follows: the melting rate of graupel melting is about 2 times; melting rate of snow is 1.3 times; and the collection rate of cloud water by rain is 1.8 times. These results show that cold-rain processes enhance greatly under the polluted conditions.
It is noted that evaporation of raindrops also strongly increases in the maximum precipitation period. Raindrop evaporation in P-case is twice stronger than C-case. By the above analysis, it can be inferred that cold-rain processes play a key role in precipitation enhancement during the maximum precipitation period, especially the melting of graupel and snow. Therefore, it is important to study the formation of graupel and snow, which helps understand the impacts of aerosol on precipitation formation mechanism. Table 1. Source-sink terms of rain (after Morrison and Pinto [76]) and maximum conversion rates in its microphysical processes (10 −3 g kg −1 s −1 ).

Microphysical Processes of Graupel Formation
Melting of graupel is the most important source of rain formation. Table 2 lists the various microphysical process of graupel formation. It can be seen in the table that during the early stage of precipitation, the most significant physical processes of graupel formation are: Accretion of cloud droplets by graupel to form graupel, deposition of graupel, and collection of cloud droplets by snow to form graupel.
During the maximum precipitation period, the production rate of graupel enhanced greatly due to increased upward movement in cold-cloud, as well as the presence of large amounts of supercooled water droplets. The collections of supercooled cloud droplets and raindrops by graupel are two of the biggest source items in the growth process of graupel. The followed important source items are the deposition of graupel, ice-rain collection to graupel, collection of cloud droplets by snow to form graupel, and collection of snow by super-cooled raindrop to form graupel. By contrast, during the early stage of precipitation, P-case has more super-cooled cloud water than C-case (see Figure 4), which contributes to the growth of snow and graupel through riming. In addition, Graupel also grows through the Wagner-Bergeron-Findeisen (WBF) process. These processes cause larger content of graupel in P-case than C-case (see Figure 4). During the maximum precipitation period in P-case, super-cooled cloud droplets are the main participants in forming graupel through accretion growth.

Microphysical Processes of Snow Formation
In addition to graupel's important contribution to the formation of precipitation, snow also plays an important role. Table 3 lists the various microphysical processes that can form snow, and their associated conversion rates (10 −3 g kg −1 s −1 ). During the initial precipitation stage, there are three major processes involved in the growth of snow. They are: the collection of rain-snow to form snow; accretion of super-cooled cloud droplets by snow (both called riming); and deposition of snow crystals. Of these, riming with super-cooled raindrops is of greatest significance. Comparing the conversion rate in two cases, there is more snow riming growth in P-case (0.047) than in C-case (0.036). This is related to the increased amount of cloud water in the mixed-phase region (approximately 4-7 km above the surface).
It is noted that, during the maximum precipitation period, snow growth by deposition has become the largest contributor to snow formation. Riming with super-cooled rain and cloud droplets are second and third respectively. This is perhaps related to two things: increased evaporation of cloud droplets; and increased water vapor transported to higher altitudes by enhanced frontal updrafts in the high humidity environment of the Mei-Yu front. Comparing the two cases, the maximum rate of snow production in all of the above processes is greater in P-case than in C-case. Hence, there is much more snowfall under polluted conditions. This is mainly because the P-case has more supercooled cloud water and rain water to foster the formation of snow. By analyzing the physical processes behind the formation of the main precipitable hydrometeors under different pollution conditions, we can see that raindrops mainly come from the melting behavior of graupel and snow during the maximum precipitation period. The primary growth mechanisms of graupel and snow are riming and deposition. In high-polluted environment, cold-cloud processes are more pronounced, producing more ice particles, which then melt and produce more raindrops. The invigoration of cold-rain processes in a polluted environment is the main microphysical mechanism behind increased Mei-Yu system rainfall production.
Compared with previous studies about frontal systems and severe convective systems, there are some similar results. Increased aerosols generally lead to smaller cloud particle sizes, suppress collision-coalescence process, hinder the warm-rain processes and hence delay rain-drop formation. Meanwhile, there is higher liquid water content (LWC) above the freezing level [30,77,78]. However, there are some differences in these systems. For deep convective system, the retention of higher LWC above the freezing level enhances riming process and suppresses deposition process, leading to precipitation enhancement. For frontal system, the deposition process is enhanced but riming is suppressed, causing precipitation decrease in high aerosol concentration environments [37,79]. Note that both riming and deposition processes are enhanced in Mei-Yu frontal system, especially in developing and vigorous periods. This may be due to two contributors: one is that more cloud water is transported to high levels where it evaporates and is deposited onto ice particles (WBF process); the other is that the invigorated cold-cloud process releases more latent heat, strengthens updrafts to increase the vertical transportation of water vapor. This would increase the degree of water vapor super-saturation to ice at high levels and speed up the growth of ice-phase particles by deposition.

Effects on Latent Heat
Varying aerosol concentrations affect the content and vertical distributions of hydrometeors, which inevitably brings in different latent heat. Latent heating includes three parts: condensation heating; freezing heating; and deposition heating. Likewise, latent cooling includes evaporation cooling, melting cooling, and sublimation cooling. Figure 9 shows an average vertical profile of latent

Effects on Latent Heat
Varying aerosol concentrations affect the content and vertical distributions of hydrometeors, which inevitably brings in different latent heat. Latent heating includes three parts: condensation heating; freezing heating; and deposition heating. Likewise, latent cooling includes evaporation cooling, melting cooling, and sublimation cooling. Figure 9 shows an average vertical profile of latent heating during the early stage of precipitation. It can be seen that the largest source of heating is condensation heating. Secondary sources are deposition heating, and freezing heating. The biggest source of cooling is the evaporation of cloud droplets (Figure 9a,b).
Comparing the two cases, there are only a few differences during the initial stage of precipitation (Figure 9c). At a height of 6 km, the evaporative cooling of cloud water is slightly larger in P-case than in C-case. Correspondingly, deposition heating is slightly larger in P-case than in C-case. This indicates that the higher cloud water content with smaller cloud droplets at high levels in P-case results in more evaporation, leading to the enhanced deposition process to produce ice-phase particles. That is to say, the WBF conversion proceeds faster with the polluted conditions. In addition, at a height of about 6 km in P-case, the freezing heating is a bit higher than that in C-case. This is mainly due to the more cloud droplets with smaller size in P-case being carried to higher levels ( Figures 4 and 5), where they either freeze, or are accreted by ice-phase particles, releasing more latent heat (Figure 9a,b). It is noted that increased latent heating due to freezing and deposition creates a "warmer bed" for cold-cloud development in P-case. ( Figures 4 and 5), where they either freeze, or are accreted by ice-phase particles, releasing more latent heat (Figure 9a,b). It is noted that increased latent heating due to freezing and deposition creates a "warmer bed" for cold-cloud development in P-case. Figure 10. Same as Figure 9, but for t = 21 h. As the cloud system develops and cold-cloud processes strengthen, P-case ice-phase particles release more and more latent heat from freezing and deposition. By t = 16 h, P-case latent heating has already approached C-case latent heating (not shown); At t = 17 h, the warm-cloud latent heat profiles are similar, but cold-cloud heating in P-case is clearly stronger than in C-case (not shown). From here, a shift occurs, and P-case net latent heating begins to exceed C-case. Figure 10 gives the vertical profile of the average latent heating rate in the model domain during the maximum precipitation stage (t = 21 h). Both heating (due to condensation, deposition, and Figure 10. Same as Figure 9, but for t = 21 h. As the cloud system develops and cold-cloud processes strengthen, P-case ice-phase particles release more and more latent heat from freezing and deposition. By t = 16 h, P-case latent heating has already approached C-case latent heating (not shown); At t = 17 h, the warm-cloud latent heat profiles are similar, but cold-cloud heating in P-case is clearly stronger than in C-case (not shown). From here, a shift occurs, and P-case net latent heating begins to exceed C-case. Figure 10 gives the vertical profile of the average latent heating rate in the model domain during the maximum precipitation stage (t = 21 h). Both heating (due to condensation, deposition, and freezing) and cooling (due to sublimation, melting, and evaporation) increase dramatically in the two cases (Figure 10a,b). The height of the net latent heating maximum value has moved up to over 6 km (Figure 10c). This indicates that latent heating in cold-cloud processes increased markedly with the development of the precipitation system. Comparing the two cases, not only freezing and deposition heating, but also condensation heating are larger in P-case. This suggested that P-case water vapor content may already be greater than C-case. Meanwhile, evaporation, melting, and cooling in P-case are all greater than in C-case.

The Dynamic Effects
Tao et al. [78] summarized a great deal of the recent researches regarding the impact of aerosols on convective clouds and precipitation and pointed out that aerosols impact clouds and precipitation by two dynamic mechanisms. The first mechanism is stronger updrafts/downdrafts resulting from enhanced latent heat release. The second mechanism is stronger convergence due to stronger evaporative cooling. These results were mostly derived from researching strong convective systems. Is there a difference between the physical mechanisms pointed out by Tao et al. [78] and those at work in the Mei-Yu frontal system? How do the various physical processes or mechanisms work together? Those questions motivate us to do more research.

Direct Dynamic Effect
The variations of aerosol concentration result in differences in microphysical processes, leading to differences in latent heating and cooling. The vertical cross sections of zonal-mean latent heat and cross-front wind perturbation vectors (v , w) (here, v = v − V g ) are shown in Figure 11. At t = 21 h, the main heating region is located in the ascent region of frontal zones and the main cooling region is located in the descent region of frontal zones. Compared to C-case, P-case heating and cooling are both much greater. In P-case, heating reaches 18 K h −1 and cooling reaches −6 K h −1 , which is 2-3 times larger than C-case. Strong heating invigorates updrafts, while strong cooling invigorates downdrafts, causing both stronger updraft and downdraft in P-case than in C-case (Figure 11). At the same time, P-case convergence in the low-level frontal region and divergence at high altitudes are both greater than C-case ( Figure 12).
Based on the above analyses, under polluted conditions, a large number of small droplets are transported to high levels by frontal updraft, leading to the enhanced riming and WBF processes to produce more ice-phase droplets, which increases latent heat release. As a result, the extra heat reinforces the updrafts and thus enhances cold-cloud processes, leading to the formation of large amount of snow and graupel. On the other hand, with large amounts of snow and graupel falling and melting, rain water content increases below the frontal zone, leading to more evaporation. Stronger melting and evaporation cooling produce a stronger cold pool and strengthens downdrafts, resulting in stronger low-level convergence under the action of frontal wind shear.
Here, the roles of latent heat (strengthens updrafts, downdrafts, and low-level convergence) under polluted conditions involves two mechanisms pointed out by Tao et al. [78], which are the latent heat-dynamic effect and the cool pool effect. We suggest that the two direct effects driven by latent heat release may be collectively referred to as the first dynamic effect/direct dynamic effect (DDE). resulting in stronger low-level convergence under the action of frontal wind shear.
Here, the roles of latent heat (strengthens updrafts, downdrafts, and low-level convergence) under polluted conditions involves two mechanisms pointed out by Tao et al. [78], which are the latent heat-dynamic effect and the cool pool effect. We suggest that the two direct effects driven by latent heat release may be collectively referred to as the first dynamic effect/direct dynamic effect (DDE).

Impacts on Frontogenesis
Significant effects of aerosols on cloud microphysics and dynamics of MCSs have been found [25,30,77]. For the Mei-Yu front, a synoptic-scale system, whether or not, and how, did aerosols associated with thermodynamic and dynamic effects influence the frontal structure and intensity, have not been studied yet.
A significant character of Mei-Yu front is strong humidity gradient across the front. Therefore, it is more suitable to describe the intensity of Mei-Yu fronts by using virtual potential temperature gradient variations [60]. The vertical cross sections of zonal mean virtual potential temperature perturbation (  Figure 13. During the initial stage of precipitation (t = 15 h), the differences of the frontal boundaries and the cross-front circulation between the two cases are negligible. However, there is a slight difference in the intensity of the front (Figure 13a,b). Comparing the tendency of the virtual potential temperature gradients in the two cases, the difference is obvious at high levels around 6 km. While there is less of difference at low levels below 2 km. This is in agreement with the latent heat distribution.

Impacts on Frontogenesis
Significant effects of aerosols on cloud microphysics and dynamics of MCSs have been found [25,30,77]. For the Mei-Yu front, a synoptic-scale system, whether or not, and how, did aerosols associated with thermodynamic and dynamic effects influence the frontal structure and intensity, have not been studied yet.
A significant character of Mei-Yu front is strong humidity gradient across the front. Therefore, it is more suitable to describe the intensity of Mei-Yu fronts by using virtual potential temperature gradient variations [60]. The vertical cross sections of zonal mean virtual potential temperature perturbation , the tendency of the virtual potential temperature gradient ( ∂ ∂t (∇θ v )), and cross-front wind perturbation vectors (v , w) (here, v = v − V g ) are shown in Figure 13.
During the initial stage of precipitation (t = 15 h), the differences of the frontal boundaries and the cross-front circulation between the two cases are negligible. However, there is a slight difference in the intensity of the front (Figure 13a,b). Comparing the tendency of the virtual potential temperature gradients in the two cases, the difference is obvious at high levels around 6 km. While there is less of difference at low levels below 2 km. This is in agreement with the latent heat distribution. Appl. Sci. 2019, 9, x 20 of 34  During the developing stage of precipitation (t = 18 h), the frontal boundaries of both cases move around 20 km to the south. Low-level convergence, high-level divergence, and cross-front ageostrophic circulation all intensify (Figure 13c,d). At this time, the greatest virtual potential temperature perturbations are above 2 K, which is at a height of 6-8 km in both cases. With development of the precipitation system, both high-level heating and low-level cooling increase further, causing the tendency of the virtual potential temperature gradient to increase at both high and low levels. In comparison, the tendencies of virtual potential temperature gradient at a height of 6 km and below 2 km enhanced significantly in P-case. It is noted that the maximum tendencies of virtual potential temperature gradient also exceeds 10 k/100 km · h at low levels. Thus, Frontogenesis associated with Mei-Yu front strengthens significantly at low levels.
At t = 21 h (Figure 13e,f), the frontal boundary in both cases moves slightly to the south. Low-level convergence, high-level divergence, and cross-front vertical circulation are all rather pronounced. By this time, the greatest virtual potential temperature perturbations have exceeded 4 K in both cases, which are still at a height of 6-8 km. The structure of the lager tendency of virtual potential temperature gradient matches with the latent heat structure (Figures 12 and 13e,f).
In comparison, the region of virtual potential temperature perturbations over 4 K was slightly larger in P-case than in C-case. The tendency of virtual potential temperature gradient is still increasing from t = 18 h to 21 h at low levels in P-case. Nevertheless, the tendency of virtual potential temperature gradient has weakened slightly by t = 21h in C-case (Figure 13c,e). This shows a continuous intensification of low-level frontogenesis in P-case. However, in C-case, low-level frontogenesis weakens, and it is much weaker than P-case.
Overall, under polluted conditions, aerosols have little effect on Mei-Yu front structure, but significantly increase the intensity of the low-level Mei-Yu front through its microphysical feedback interactions.
To explore the physical processes behind the influence of aerosols on frontogenesis, we calculate frontogenesis function [80] using the outputs from C-case and P-case numerical experiments. The expression of frontogenesis function is given by: where F represents the frontogenesis function, with F > 0 for frontogenesis and F < 0 for frontolysis; F 1 , F 2 , F 3 , and F 4 represent horizontal divergence term, horizontal deformation term, tilting term associated with vertical motion, and diabatic heating term, respectively. During the maximum precipitation period (21 h), there are three obvious frontogenesis regions at heights of 1-3 km, 4-5 km, and 10-11 km. P-case F-values are larger than C-case (Figure 14a,b). We primarily focus on the low-level effects. Below the height of 3 km, the frontogenesis equation is dominated by the diabatic and vertical velocity terms in both cases. The dispersion term and the horizontal deformation term also have a role. In contrast, all types of frontogenesis are larger in P-case than in C-case. Diabatic frontogenesis is mainly due to temperature gradient increases across the front. These increases are the result of the combined effects of condensation heating ahead of the front (forming a warm pool) and evaporative cooling behind the front (forming a cool pool). This is thermal frontogenesis. The other terms in the equation are all induced by the DDE of latent heat, belonging to dynamic frontogenesis.
Appl. Sci. 2019, 9, x 22 of 34 and vertical velocity terms in both cases. The dispersion term and the horizontal deformation term also have a role. In contrast, all types of frontogenesis are larger in P-case than in C-case. Diabatic frontogenesis is mainly due to temperature gradient increases across the front. These increases are the result of the combined effects of condensation heating ahead of the front (forming a warm pool) and evaporative cooling behind the front (forming a cool pool). This is thermal frontogenesis. The other terms in the equation are all induced by the DDE of latent heat, belonging to dynamic frontogenesis.  ; (e,f) F 3 (contour) and F 4 (shading). F represents the frontogenesis function; F 1 , F 2 , F 3 , and F 4 represent horizontal divergence term, horizontal deformation term, tilting term, and diabatic heating term, respectively. The left panels show C-case, the right panels show P-case.
High aerosol concentrations influence frontogenesis through microphysical effect and DDE. On the one hand, enhanced latent heat release (heating and cooling) produce more intense warm and cool pools, invigorating thermal frontogenesis. On the other hand, the latent heat feedback effect strengthens updrafts ahead of the front and downdrafts behind the front, invigorating convergence near the frontal boundary and, hence, dynamic frontogenesis. Briefly, during precipitation development and maximum precipitation periods, high aerosol concentrations in Mei-Yu systems strengthen both thermal and dynamic low-level frontogenesis through microphysical feedback, leading to an intensified low-level front.

The Frontogenesis Effect
For frontal systems, frontogenesis can trigger ageostrophic circulation and updrafts, lifting vapor to produce clouds and precipitation. At the initial stage of cloud formation, the role of frontal circulation is significant. However, the convective system develops rapidly if it is located in both super-saturated and unstable environments. At this stage, the thermal circulation driving by latent heat dominates and overwhelms frontal circulation [42,81]. Actual circulation is the superimposition of the two (Figure 15). deformation term, tilting term, and diabatic heating term, respectively. The left panels show C-case, the right panels show P-case.
High aerosol concentrations influence frontogenesis through microphysical effect and DDE. On the one hand, enhanced latent heat release (heating and cooling) produce more intense warm and cool pools, invigorating thermal frontogenesis. On the other hand, the latent heat feedback effect strengthens updrafts ahead of the front and downdrafts behind the front, invigorating convergence near the frontal boundary and, hence, dynamic frontogenesis. Briefly, during precipitation development and maximum precipitation periods, high aerosol concentrations in Mei-Yu systems strengthen both thermal and dynamic low-level frontogenesis through microphysical feedback, leading to an intensified low-level front.

The Frontogenesis Effect
For frontal systems, frontogenesis can trigger ageostrophic circulation and updrafts, lifting vapor to produce clouds and precipitation. At the initial stage of cloud formation, the role of frontal circulation is significant. However, the convective system develops rapidly if it is located in both super-saturated and unstable environments. At this stage, the thermal circulation driving by latent heat dominates and overwhelms frontal circulation [42,81]. Actual circulation is the superimposition of the two (Figure 15). Although it is difficult to separate the two types of circulation, ageostrophic circulation generated by frontogenesis is objectively present; and stronger frontogenesis generates stronger ageostrophic circulation. Note that frontal circulation is not the same as typical thermal convective circulation. The interaction of wind shear and the cool pool at the frontal boundary separates updraft and downdraft corridors, forming tilted frontal circulation. This is similar to the slantwise convective circulation generated by wind shear and cool pools that make precipitation systems more "wellorganized", strengthening intensity and prolonging its lifecycle [82]. High pollution aerosol has obvious impacts on both the thermal and dynamic fields of Mei-Yu frontal systems. As clouds and precipitation develop along the Mei-Yu front in the polluted conditions, frontogenesis is markedly Although it is difficult to separate the two types of circulation, ageostrophic circulation generated by frontogenesis is objectively present; and stronger frontogenesis generates stronger ageostrophic circulation. Note that frontal circulation is not the same as typical thermal convective circulation. The interaction of wind shear and the cool pool at the frontal boundary separates updraft and downdraft corridors, forming tilted frontal circulation. This is similar to the slantwise convective circulation generated by wind shear and cool pools that make precipitation systems more "well-organized", strengthening intensity and prolonging its lifecycle [82]. High pollution aerosol has obvious impacts on both the thermal and dynamic fields of Mei-Yu frontal systems. As clouds and precipitation develop along the Mei-Yu front in the polluted conditions, frontogenesis is markedly intensified (Figure 15). Stronger, continuing frontogenesis benefits and maintains cross-front ageostrophic circulation ( Figure 16). It also strengthens and prolongs slantwise updraft, resulting in a more "well-organized" Mei-Yu frontal system. This leads to the more sufficient interactions of the microphysical and dynamic processes influenced by aerosol and enhances precipitation. Hereinafter, "the frontogenesis effect" refers collectively to forced cross-frontal ageostrophic circulation, strengthened and maintained slantwise updraft, and organized precipitation system. The frontogenesis effect is an indirect dynamic effect resulted from microphysical processes, latent heat release, and its DDE. In order to distinguish it from the DDE, we suggest that the indirect dynamic effect may be referred to as the second type of dynamic effect-frontogenesis dynamic effect (FDE).

Third Dynamic Effect: The Water Vapor Pump Effect
Mei-Yu front is an important part of the East Asian monsoon rainfall system, which is associated with plenty of water vapor conditions. From the above analyses, aerosol has significant influence on the dynamic field of Mei-Yu front. In addition, this rainfall system usually maintains long time, which could have important influences on the transportation of water vapor, even though very few studies have investigated the indirect vapor effect of aerosols. Figure 16 shows the vertical cross sections of zonal mean vertical water vapor flux at t = 21 h for both cases. Maximum vertical water vapor flux in P-case is approximately double C-case. Horizontal water vapor flux is also clearly greater in P-case than C-case ( Figure 17). Note that around 200 km south of the rain band, P-case horizontal water vapor flux is still 5% higher than C-case ( Figure 17). This indicates that both vertical and horizontal transportation of water vapor increase significantly of Mei-Yu frontal system under polluted conditions.
From the previous sections, in polluted Mei-Yu frontal environments, both the first dynamic effect (DDE) and second dynamic effect (FDE) of aerosols can enhance low-level convergence and updrafts, which can increase low-level water vapor convergence near the front, leading to the enhanced transportation of water vapor above the frontal zone. On the other hand, in the high humidity conditions associated with the southwest monsoon, strong vertical water vapor transportation and consuming, in turn, increase the horizontal water vapor transportation in the larger-scale domain (Figure 17c). Overall, the increased updrafts and vapor transportation further enhance condensation, deposition, and collision-coalescence processes to produce more precipitable particles and to release more latent heat, which in turn can increase updrafts and further strengthen the positive feedback loop between microphysical and dynamic processes. Thus, the enhanced vapor transportation effect is driven by the first and second dynamic effects, which are rooted to latent heat release. We suggest this effect may be referred to as the third aerosol dynamic effect-vapor pump effect (VPE).
In fact, the VPE has been presented in Section 3.2 and Section 3.3, where it was pointed out that, under polluted conditions, not only freezing and deposition but also condensation is increased ( Figure 10). In other words, as cold-cloud processes strengthen, warm-cloud processes also strengthen. During these processes, enhanced water vapor transportation should play an important role.
This increased water vapor transportation effect is closely related to the Mei-Yu frontal system during monsoon precipitation. Southwest flow brings plenty of water vapor to the Yangtze River Basin, forming precipitation and even long-lasting torrential rainstorms in suitable convergence conditions. Under these environments, the proper convergence near the frontal boundary can drive large-scale domain water vapor to Mei-Yu front. During the developing and mature stage of the Mei-Yu frontal precipitation system, the microphysical-thermal-dynamic coupling effect created by high concentration aerosols causes stronger updrafts, high level divergence and low-level convergence, forming the VPE. This effect enhances water vapor transportation, and further leads to more rainfall production.

Physical Mechanism Discussion
From the analyses of Section 3, it can be noticed that the variation of aerosol concentrations results in the differences of microphysical processes and latent heat release in Mei-Yu frontal system. The differences successively induce dynamic and water vapor variations, which in turn influence microphysical processes. These physical processes create a positive feedback loop, leading to the increased precipitation under polluted conditions in Mei-Yu front. In general, four physical mechanisms are proposed to explain the enhancement of precipitation in Mei-Yu front. The first is the cloud microphysical mechanism, including microphysical processes and the latent heat release. The other three mechanisms are three dynamic effects resulted from the microphysical effects. They are the DDE, FDE, and VPE, respectively.
Previous studies indicated that riming increased and deposition decreased [77] or deposition enhanced and riming suppressed to produce ice phase particles under polluted conditions [37]. However, our finding shows that not only riming but also deposition enhance to produce snow and graupel in Mei-Yu frontal precipitation systems. There are two possible explanations for the enhanced deposition. First, WBF process is enhanced, because more cloud droplets are transported above the 0 • C isotherm where they may evaporate and deposited onto ice-phase particles. Second, the increased cold-cloud processes enhance latent heat release, leading to stronger updrafts and transportation of water vapor to high levels to strengthen deposition. On the other hand, the increased deposition releases more latent heat and increases updrafts that enhance riming process. The mutual promotion of riming and deposition increases the growth of snow and graupel in the strong updraft region and then increases cold-rain process.
During the developing and maximum precipitation stages of Mei-Yu front system, the enhanced melting and evaporation can strengthen dynamic processes at the low-level frontal boundary, and drive and strengthen warm-cloud processes (condensation and collision-coalescence). Likewise, strengthened warm-cloud processes release more condensation heating, strengthening updrafts, which in turn promotes cold-cloud processes. Thus, with the cold-cloud "towing" effect, there is mutual promotion of cold and warm-cloud processes. This could be an even more complicated and peculiar cloud microphysical mechanism of aerosol effects on Mei-Yu front precipitation.
Regarding latent heat release, previous researchers emphasized that increased high-level freezing heating was the major difference under polluted conditions [25,77]. Our research shows that in polluted Mei-Yu frontal cloud systems, high-level heating due to both freezing and deposition processes increases, but deposition heating is much greater than freezing heating. In addition, due to the increased high-level cold-cloud heating, there is increased updraft that not only increases collision-coalescence processes, but also greatly increased vertical transportation of water vapor because of the high humidity environment of Mei-Yu fronts. This leads to deposition and condensation processes all increasing as well ( Figure 10). Furthermore, deposition and condensation processes cause much greater heating than freezing. Because of this, the key mechanism behind increased heating is increased deposition and the condensation processes with the cold-cloud "towing". On the other hand, previous studies pointed that low-level evaporative cooling had a key impact on low-level convergence [30]. Our studies indicated that evaporative cooling is the main driver of low-level convergence, and melting cooling also plays a role (Figures 11 and 12).
Regarding dynamic effects, the similarities between the previous results of convective systems [25,30,34,77,78] and our findings of Mei-Yu front are that increased heating strengthens updrafts, and increased cooling strengthens downdrafts, which in turn strengthens low-level convergence. However, there are also some differences and new insights in our findings. We refer to the stronger updrafts (due to the latent heat-dynamic effect) and convergence (due to the cool pool effect) pointed out by Tao et al. [78] as DDE, which has a very significant effect on the precipitation enhancement in convective clouds and mixed clouds. However, in Mei-Yu system, in addition to DDE, FDE is also an important mechanism for the precipitation enhancement under polluted conditions. Under the role of the frontal vertical wind shear, FDE forces ageostrophic circulation and forms slantwise convective circulation patterns (updrafts and downdrafts are separated) which contributes to strengthening and maintaining updrafts, causing a more "well-organized" Mei-Yu frontal system. And this further strengthens and prolongs precipitation.
Regarding the VPE, Khain et al. [25] mentioned that increased squall line precipitation maybe related to the gathering of water vapor from a larger region into a smaller region, but they did not analyze further. This paper presents a "water vapor pump" effect in Mei-Yu front caused by aerosols in positive feedback loops. It is particularly pointed out that Mei-Yu front is monsoon precipitation system. With the large-scale sufficient water vapor conditions, deep frontal clouds, and long-time precipitation, the long-lasting water vapor convergence and vertical transportation can tow and drive large-scale domain water vapor to Mei-Yu front. The enhanced water vapor transportation effect, in turn, strengthens condensation, deposition and even collision-coalescence processes, leading to more latent heat release, which further enhance the "microphysical-dynamic feedback loop". Note that this effect has not been reported in the previous studies yet.
The four physical mechanisms resulted from high aerosol concentrations in Mei-Yu frontal systems have different influences in different stage of precipitation. The first mechanism, the microphysical mechanism, affects the entire precipitation process. Especially, the microphysical effects on precipitation dominate the early stages of precipitation. With the precipitation development, the other three dynamic effects gradually become active. Among them, DDE is the most direct. FDE are driven by latent heat release and DDE in Mei-Yu frontal environments. The increasing and maintaining low-level frontogenesis can strengthen the intensity of the original front and play an important role in strengthening and prolonging the precipitation. That is to say, FDE only in the frontal precipitation system is more significant. Furthermore, VPE becomes obvious only in Mei-Yu front associated with monsoon systems. It may not be obvious in short-lived isolated convection systems. The aerosol impacts on Mei-Yu frontal clouds and precipitation are summarized in a schematic diagram ( Figure 18). The figure describes the main physical processes through which high aerosol concentration affects Mei-Yu frontal clouds and precipitation via four mechanisms. High aerosol concentrations first influence cloud microphysical processes and their accompanying latent heat release. Increased updrafts, downdrafts, and low-level convergence create the direct dynamic effect. The direct dynamic effect re-invigorates the cloud microphysical processes that create it, forming the first positive feedback loop. In Mei-Yu frontal system, the cool pool effect and direct dynamic effect caused by microphysical processes result in frontogenesis which increases frontal dynamic processes and strengthens the slantwise updrafts, thereby forming the VPE. Increased water vapor The aerosol impacts on Mei-Yu frontal clouds and precipitation are summarized in a schematic diagram ( Figure 18). The figure describes the main physical processes through which high aerosol concentration affects Mei-Yu frontal clouds and precipitation via four mechanisms. High aerosol concentrations first influence cloud microphysical processes and their accompanying latent heat release. Increased updrafts, downdrafts, and low-level convergence create the direct dynamic effect. The direct dynamic effect re-invigorates the cloud microphysical processes that create it, forming the first positive feedback loop. In Mei-Yu frontal system, the cool pool effect and direct dynamic effect caused by microphysical processes result in frontogenesis which increases frontal dynamic processes and strengthens the slantwise updrafts, thereby forming the VPE. Increased water vapor transportation in turn strengthens microphysical processes. Thus, these processes form the second positive feedback loop.
The interaction of several physical processes form and strengthen the two types of "microphysical-dynamic positive feedback loop", leading to the more "well-organized" Mei-Yu front system and enhanced precipitation. The interaction of microphysical processes and frontal dynamic processes, and the two positive feedback loops they create, are the major mechanisms by which aerosols affect Mei-Yu frontal precipitation.

Summary and Conclusions
In this study, the impact of polluted (P-case) conditions versus typical (C-case) continent conditions on an ideal Mei-Yu front was simulated via WRF coupled with a bulk two-moment microphysics scheme [63]. The primary focus was on exploring physical mechanisms of simulated aerosol effects in the face of complex process interactions and feedbacks between the cloud microphysics and dynamics in Mei-Yu frontal system. We summarized three dynamic effects resulted from the microphysical mechanisms and revealed two microphysical/dynamic positive feedback loops in polluted Mei-Yu frontal environments. The main conclusions are as follows: (1) Aerosols play an important role in affecting local precipitation around Mei-Yu fronts. Polluted conditions delay the onset of precipitation but strengthen precipitation during the intense precipitation period and increase the amount of total precipitation. This may be one of the reasons for the increased frequency of the intense Mei-Yu precipitation in recent years. Note that anthropogenic pollution has a significant effect on Mei-Yu front precipitation. This may be an important feature influencing climate change in East Asia. (2) Aerosols affect Mei-Yu frontal precipitation by first changing microphysical processes.
Under polluted conditions, these more numerous and smaller cloud droplets suppressed auto-conversion and collision-coalescence processes, and hence resulted in late raindrop formation, which initially delays the onset of Mei-Yu precipitation. However, a larger amount of cloud water is transported to above the freezing level, increasing ice-phase particle growth by riming and the WBF processes, and producing large amounts of snow and graupel. The melting of snow and graupel is the main reason why rain increases. Furthermore, during the developing and mature stages, cold rain processes are more vigorous under polluted conditions. With cold-cloud processes "towing", warm-cloud progresses are encouraged. Thus, cold-cloud and warm-cloud interactions mutually boost each other. This is perhaps the more complicated and particular microphysical mechanism by which aerosols affects Mei-Yu frontal precipitation. (3) High aerosol concentrations first influence cloud microphysical processes, enhance latent heat release, and then strengthen updrafts, downdrafts, and low-level convergence, forming DDE. In turn, DDE further strengthens microphysical processes and creates the first positive feedback loop. In Mei-Yu frontal environments, the enhanced DDE and the cool pool effect result in stronger frontogenesis, which strengthens frontal dynamic processes. This further increases both vertical and horizontal transportation of water vapor and form VPE, which in turn strengthens microphysical processes and creates the second positive feedback loop. Many physical processes and effects interact with each other, forming and strengthening both microphysical-dynamic feedback loops, leading to the well-organized Mei-Yu front precipitation system. Overall, the combined effect is to increase Mei-Yu front precipitation. The interaction of microphysical processes and dynamic processes, and the positive feedback loops they create, are the main physical mechanisms behind the impacts of aerosol on Mei-Yu frontal precipitation.
Notably, a series of ideal frontal experiments are carried out herein by using a domain-averaged initial sounding of a typical Mei-Yu system. The results could represent the typical Mei-Yu frontal systems, but not all of Mei-Yu frontal temperatures and humidity. On the other hand, the cloud physics scheme used in this paper simplifies the cloud-aerosol interactions by setting cloud droplet concentrations to reflect changes in aerosol concentration. In addition, it does not consider the impacts of aerosols on ice nuclei. The above simplifications may affect the accuracy of our conclusions about the impacts of aerosols on the magnitude of precipitation. This paper investigated the overall impacts of increased aerosols on Mei-Yu front precipitation in highly polluted regions, such as the Yangtze River Basin in China. In reality, Mei-Yu fronts can extend from interior of southern China to southern Japan over varying geography and pollution conditions. The impacts of local atmospheric pollution on the entire Mei-Yu rainband will be investigated in future work.

Conflicts of Interest:
The authors declare no conflict of interest.