Numerical Simulation of the Effect of Cloud Condensation Nuclei Concentration on the Microphysical Processes in Typhoon Usagi

3e Weather Research and Forecasting model version 3.2.1 with the Lin microphysics scheme was used herein to simulate super typhoon Usagi, which occurred in 2013. To investigate the effect of the concentration of cloud condensation nuclei (CCN) on the development of typhoon Usagi, a control simulation was performed with a CCN concentration of 100 cm, together with two sensitivity tests: C10 and C1000, having CCN concentrations of 10 cm and 1000 cm, respectively. 3e path, intensity, precipitation, microphysical processes, and the release of latent heat resulting from the typhoon in all three simulations were analyzed to show that an increase in CCN concentration leads to decreases in intensity and precipitation, an increase of the cloudless area in the eye of the typhoon, a more disordered cloud system, and less latent heat released through microphysical processes, especially the automatic conversion of cloud water into rainwater. Overall, an increase in CCN concentration reduces the total latent heat released during the typhoon suggesting that typhoon modification by aerosol injection may be optimized using numerical simulations to ensure the strongest release of latent heat within the typhoon.


Introduction
As one of the most important disaster weather phenomena, torrential rain resulting from typhoons has been studied from different perspectives [1][2][3][4], but with the consistent goal of understanding the generation mechanisms and development of typhoons, and especially the formation of torrential rain, with the help of numerical simulations.In the 1950s, a numerical procedure for the simulation of the development of tropical cyclones based on a barotropic model was proposed by Kasahara [5,6].Next, the development of a typhoon was simulated using a three-level, three-dimensional model, which obtained a spatial structure consistent with observations in [7].e eye-wall structure of a tropical cyclone was simulated by using a nested mesh model in Kurihara and Bender [8].Numerical simulation in Lin et al. [9] established the degree to which the vertical vapor flux is affected by terrain and its importance for the prediction of typhoon precipitation.Simulations also identified that evaporation of marine spray caused by the typhoon influences the heat exchange between the typhoon and the sea surface, leading to the convective inhibition of eye-wall motions, and thus a weakening of the typhoon [10].Mesoscale atmospheric numerical modeling has developed considerably since the 1980s.Beginning around 1997, scientists from several research institutes in the United States jointly developed a new generation of high-resolution mesoscale forecasting models, including the Advanced Research Weather Research and Forecasting (WRF) model [11], which has since been widely used in practical applications.
e impact of cloud condensation nuclei (CCN) on cloud microphysical processes is handled by the cloud microphysics scheme in numerical models.A squall-line process was simulated over Guangdong, China, using the WRF model with three different concentrations of CCN and found that an increase in CCN concentration delays precipitation, weakens the initial precipitation, and lessens the amount of precipitation for high CCN concentrations [12].Furthermore, an increase in CCN concentration leads to increases of cloud water and snow and the decrease of sleet [13].Lin et al. [14] simulated the typhoon Chanchu using the MM5 model to compare the effects of different CCN concentrations; they found that the cloud water growth increases with the concentration of CCN, leading to a positive relationship between the cloud-water mixing ratio and CCN concentration.
Cloud microphysics account for the transformation of hydrometeors, with their overall effect directly impacting the performance of model simulations.Cloud microphysics schemes include, for example, the Lin et al. [15][16][17], WSM3 [18,19], WSM6 [20], Ferrier [21], ompson et al. [22], and Morrison et al. [23] schemes.Niu and Yan [24] simulated a typical precipitation process in the Jianghuai area of China using the WRF model version 2.1.2and compared the predictions of heavy rainfall based on the Lin, WSM3, WSM6, Ferrier, and ompson cloud microphysical schemes, with the Lin scheme giving the best performance.Li et al. [25] simulated deep convection over the Tiwi Islands in northern Australia and compared simulations using both Lin and WSM6 schemes, finding a better agreement between simulation results and observations for the Lin scheme.Khin and Min [26] showed that the Lin microphysical scheme produced more heating within the 300-100 hPa layer compared to both Ferrier and WSM6 schemes, producing storms with greater intensity.
Typhoon Usagi is simulated herein with the help of several sensitivity experiments performed using the WRF model version 3.2.1 to explore the effects of CCN concentration on the typhoon track, intensity, and cloud microphysical processes; these aspects are evaluated using the development and movement of the typhoon, as well as the resulting precipitation.

Simulation Strategy.
e WRF model version 3.2.1 was operated with a Mercator map projection of twin nested grids, with a horizontal grid spacing of 9 km for domain 1 (D01) and 3 km for domain 2 (D02) (Figure 2).e Lin scheme [15] was selected as the cloud microphysics scheme.e rapid radiation transfer model long-wavelength and Goddard shortwavelength radiation parameterizations were used in both domains.e G3 cumulus parameterization was applied in D01, while no cumulus parameterization was applied in D02 because the 3 km grid resolution could not distinguish cumulus processes at the cloud scale.Uncertainty of the cumulus parameterization scheme was much greater than that of the explicit cloud scheme.Typically, the latter scheme describes cloud microphysical processes better and more reasonably represents cloud structure [27].e 6 h FNL (final) analysis data provided by the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) with 1 °horizontal resolution were used as the initial field of the model.e modeling period in the D01 domain amounted to 84 h from 1200 UTC on 19 September 2013 to 0000 UTC on 23 September 2013.
e modeling period in the D02 domain amounted to 60 h from 1200 UTC on 20 September 2013 to 0000 UTC on 23 September 2013.e integral time steps of domains D01 and D02 were 30 s and 10 s, respectively.
In the Lin scheme, the cloud droplet number concentration (N c ), which represents the CCN concentration, was set to 100 cm −3 in the control test (CTL) but varied in the sensitivity tests (C10, C1000) to reflect clean-air and heavy air pollution scenarios, with all other parameters held constant.

e Effect of Cloud Condensation Nuclei Concentration.
Figure 3 compares the simulated paths of all three experiments with the actual path of the typhoon.While the simulated landfall sites in the three experiments deviate from the observed landfall site, the simulated paths are largely consistent, with each other and with the observed path.Figure 4 compares the simulated sea-level pressure in all three experiments with pressures observed.Initially, the simulated intensity was weaker than that observed.Consequently, simulations gave a reduced typhoon speed compared to that observed.e comparison of the three simulations indicates that the lowest sea-level pressure increases with increasing CCN concentration, with the lowest sea-level pressure simulated in C1000 being consistently higher (by 2 hPa) than that in C10, as well as in the C100 case.
erefore, the typhoon intensity in the simulations decreased with increasing CCN concentration.

2
Advances in Meteorology

e Effect of Cloud Condensation Nuclei Concentration on the Usagi Precipitation Simulation. Figure 1 compares
the 24 h accumulated precipitation on land in the CTL experiment with that observed; it shows that the amount of precipitation reached a maximum in the CTL experiment of over 300 mm, which is larger than that observed.In terms Advances in Meteorology of the spatial distribution of precipitation, the area of torrential rain (24 h accumulated precipitation of more than 100 mm) in the CTL simulation contradicts the observations to the west of the Pearl River Delta region (resulting from the deviation of the typhoon path) but is consistent with observations to the east of the Pearl River Delta region.In addition, the heavy rain (24 h accumulated precipitation less than 50 mm) area in the CTL results approximates observations in North Guangdong, but with lower precipitation levels in North and West Guangdong.Contrasting this with the sensitivity tests, in the case of C10, large values of precipitation of more than 300 mm are centered to the west of Shantou City (116.2 °E, 23.2 °N), but with a reduced area in the CTL experiment; this disappears in the results for C1000.Moreover, the region of approximately 200-300 mm precipitation is significantly reduced in C1000.erefore, the precipitation decreased with increasing CCN concentration.

e Effect of Cloud Condensation Nuclei Concentration on Radar Reflectivity and the Distribution of Cloud-Top
Temperature. Figure 5 compares the simulated radar maximum reflectivity with that observed at 0000 UTC on 22 September.While much larger eye-wall regions are evident in all three experiments compared with that observed, the strong radar reflectivity observed in the northern part of the typhoon close to land was captured correctly.e typhoon eye area increased with increasing CCN concentration, according to the three experiments, with C10 producing a larger area of radar reflectivity (>50 dBZ) compared with the CTL and C1000 runs.
Figure 6 shows the distribution of cloud-top temperature in all three experiments at 0000 UTC on 22 September.Clearly, the cloud-top temperature was highest in C1000, followed by the CTL and C10 runs, indicating a decrease of the cloud-top height within the typhoon with increasing CCN concentration.e cloud cover in C10 is more evenly distributed in contrast to the fragmented clouds in the eyewall region in the CTL run and the disconnected regions of large-scale cloud cover in C1000.Because the magnitude of the CCN concentration affects the distribution of cloud cover, this indicates the importance of CCN to the development of typhoon intensity and precipitation.

e Effect of Cloud Condensation Nuclei Concentration on
the Conversion of Cloud Microphysical Parameters.By concentrating on domain D02 (3 km resolution) and periods when the typhoon was still located offshore (from 1700 UTC to 2200 UTC on 21 September) as well as when the typhoon made landfall (from 1500 UTC to 2000 UTC on 22 September), the effect of CCN concentration on typhoon development was analyzed in detail by considering the vertical integral of a particular conversion parameter [28], where [P x ] is the vertical integral of the microphysical process X, ρ is the density of air (kg•m −3 ), Δz is the distance between the adjacent layers in the vertical direction, and p x is the hourly microphysical conversion on each grid Table 1 shows results of the vertical integration using equation ( 1) of the microphysical conversions in domain D02 for all three experiments in both periods.e cloud microphysical conversions within the period when the typhoon made landfall are smaller than those within the period when the typhoon was offshore, which indicates that the typhoon weakened.is is probably related to the increased friction of the underlying surface, as the typhoon passed through the Bashi Channel.
e microphysical processes correspondingly slowed, along with the weakening of the typhoon.Most cloud microphysical conversions decreased with increasing CCN concentration, except for the CEVP process (the conversion of cloud water to vapor related to evaporation).In this case, cloud water decreased, which is conducive to reducing precipitation.e automatic conversion of cloud water to rain in the Lin scheme [15] is calculated by where N 1 and l CW represent the number concentration of cloud droplets and the cloud-water mixing ratio, respectively; l W0 (2 × 10 −3 g•g −1 in the model) represents the automatic conversion threshold; and D 0 (a constant in the model) represents the dispersion of the cloud droplet distribution.In reducing the microphysical conversion rate, RAUT (the process whereby cloud water automatically converts to rain) is weakened significantly.In the first period, for example, RAUT decreases significantly from 1.47 mm•h −1 to 0.01 mm•h −1 with increasing CCN concentration.Equation (2) shows that the increase in CCN concentration (through the increase in N 1 ) inhibits RAUT, leading to the weakening of cloud microphysical processes relevant to rain.From the microphysical perspective, the  Advances in Meteorology increasing CCN concentration results in enhanced competition between CCN for the limited water vapor, leading to a decrease in the cloud droplet size.As a result, the probability that cloud droplets automatically convert to larger rain droplets decreases, which is the core process, when considering the physics of the influence of CCN concentration on precipitation.

e Effect of Cloud Condensation Nuclei Concentration on the Release of Latent Heat during Usagi.
e release of latent heat by various microphysical processes is an important aspect in the development of typhoons.Of the 33 cloud microphysical processes in the Lin scheme, 19 are able to produce latent heat through phase transitions resulting from condensation (Q con ), evaporation (Q evp ), freezing (Q frz ), melting (Q mlt ), deposition (Q dep ), and sublimation (Q sub ); these may be calculated at each integral step in the model to obtain the total heat released over a certain period (Q total ) using where P X (kg•kg −1 •s −1 ) is the transition rate of a cloud microphysical process; L v , L f , and L s (J•kg −1 ) are the latent heats of condensation, freezing, and deposition per unit mass, respectively; and C pm (J•kg −1 •K −1 ) is the specific heat of moist air at constant pressure.e regional time-averaged vertical profile of total latent heat produced from 1200 UTC to 1800 UTC on 22 September in domain D02 (114 °E-116 °E, 20.5 °N-22.5 °N) is presented in Figure 7. is illustrates the consistent trend of two peaks, one of which is located at the height of 2 km, and the other between 5 km and 6 km, on the total latent heating rate in all three experiments.Typically, the latent heating rate decreased rapidly at 2-4 km and increased above 4 km.Further analysis of the latent heat contributions (diagram not shown) indicates that the first peak depends mainly on the degree of latent heat from condensation, while the second peak depends mainly on the degree of latent heat from condensation, freezing, and melting.Comparing the three experiments, only the second peak of C10 is larger than the first one.In addition, the total latent heating rate corresponding to the two peaks in C10 is obviously larger than that in the CTL test, which is slightly larger than that in C1000.Hence, an increase of CCN is beneficial to the suppression of the release of latent heat in a typhoon.e addition of aerosols into the peripheral spiral rain bands of the typhoon is the common mode of typhoon intensity modification, where many CCN reduce the rate of cloud droplets converting into raindrops, which then rise to the height of freezing, creating more supercooled water droplets.e process by which supercooled water droplets convert to ice produces latent heat, triggering convection and enhancing the vertical motion within the spiral cloud belt on the typhoon's periphery.is motion is beneficial to the development of convective clouds and the reduction of the convergence of warm air and heat at the typhoon center.
is process increases the radius of the typhoon eye and decreases the maximum wind speed to conserve momentum, resulting in a weakened typhoon [29].erefore, typhoon modification by aerosol injection may be optimized using numerical simulations to ensure the strongest release of latent heat within typhoons.

Conclusions
e Weather Research and Forecasting model version 3.2.1 with a single nested domain was used in this study to simulate super typhoon Usagi, with a focus on the Lin 6 Advances in Meteorology  Advances in Meteorology scheme as the cloud microphysical scheme for three different CCN concentrations.ese simulations showed the following: (1) e typhoon intensity, precipitation, and cloud-top height decreased with increasing CCN concentration.e distribution of the cloud system also decreased and became more irregular with increasing CCN concentration.(2) An increasing CCN concentration corresponded to a reduced release of latent heat in most cloud microphysical processes as a result of competition between CCN for the limited water vapor, leading to a decrease in cloud droplet size.Because the probability that cloud droplets are automatically converted to larger-sized rain droplets correspondingly decreases, increasing CCN concentration affects RAUT efficiency (automatic cloud water conversion to rainwater), leading to weakening of the cloud microphysical processes related to precipitation.(3) e increasing number of CCN suppressed the release of latent heat within the typhoon.Typically, two peaks occurred on the total latent heating curve, with one peak located just below 2 km, dependent on the degree of latent heat from condensation, and the other located at a height of 5-6 km, dependent on the degree of condensation, freezing, and melting.e existence of latent heat peaks may facilitate a new way of parameterizing and analyzing the extent to which latent heat release modifies the strength of typhoons.
As only numerical simulations of a single typhoon have been carried out here, these simulations do not necessarily represent all circumstances.Clearly, further simulations are necessary to broaden our understanding of typhoon dynamics.

Data Availability
e National Center for Atmospheric Research's Mesoscale and Microscale Meteorology Division provided the WRF-Chem model available at http://www.mmm.ucar.edu/wrf/users.e National Centers for Environmental Prediction Final Global Tropospheric Analysis (NCEP-FNL) data were downloaded from https://rda.ucar.edu/datasets/ds083.2.All data used to support the findings of this study are available from the corresponding author upon request.

Figure 4 :
Figure 4: Comparison of simulated minimum sea-level pressures with observations for the C10, C100, and C1000 runs.

Figure 5 :
Figure 5: Comparison of radar maximum reflectivity in C10 (a), C1000 (b), and C100 (c) runs with observation (d) at 0000 UTC on 22 September.e radar observation is sourced from the China Weather News.

Figure 6 :
Figure 6: Comparison of the distribution of cloud-top temperature in C10 (a), C1000 (b), and C100 (c) runs at 0000 UTC on 22 September.