Research Article  Open Access
Ying Zhang, Yu Wang, Guosheng Liu, Jianping Guo, Yuanjian Yang, Rui Li, Yunfei Fu, Liping Liu, "SatelliteBased Assessment of Various Cloud Microphysics Schemes in Simulating Typhoon Hydrometeors", Advances in Meteorology, vol. 2019, Article ID 3168478, 19 pages, 2019. https://doi.org/10.1155/2019/3168478
SatelliteBased Assessment of Various Cloud Microphysics Schemes in Simulating Typhoon Hydrometeors
Abstract
The accurate simulation of typhoon hydrometeors remains a challenge. This study attempts to evaluate the performances of five microphysics schemes (MPSs) in the Weather Research and Forecasting (WRF) model in simulating the supertyphoon Neoguri in July 2014. The observed microwave brightness temperature, as well as retrieved data from the microwave radiometer imager (MWRI) onboard Chinese FY3B satellite, are used to test hydrometeor simulations. In particular, two MWRI radiance indices, including the emission index (EI) and scattering index (SI), are used to assess the performance of five MPSs in simulating liquid and frozen hydrometeors, respectively. Overall, the WRF model can well reproduce the overall pattern of typhoonproduced precipitation, albeit with slightly overestimated precipitation in the inner rainband and underestimated precipitation in the stratiform rainband. Moreover, ice water paths (IWPs) from all five MPS simulations are higher than those estimated from MWRI retrieval in most areas, and the spatial pattern and values of IWP for the National Severe Storms Laboratory doublemoment MPS (NSSL) are much closer to those for MWRI. The NSSL scheme reproduces a more realistic joint histogram distribution of SI and EI than other MPSs do, relative to the observation. Besides, the nonlinear Lucas–Kanade optical flow approach has been used to reflect the horizontal distribution of hydrometeors in the typhoon. The results show that the simulated EI and SI from the five MPSs show a systematic southwest bias of approximately about 10∼20 km and significant intensity bias in the convection area. Further model sensitivity tests confirm that the NSSL scheme generates more realistic graupel and supercooled water close to the observations among all MPSs. The findings suggest that satellite measurements would be helpful to assess MPSs in numeric weather models, especially for hydrometeor distributions in the whole typhoon system.
1. Introduction
Cyclone hydrometeors (e.g., cloud water, rainwater, cloud ice, snow, and graupel), which are related with the cyclone’s diabatic heating and vertical velocity structure, are effective indicators of storm intensity [1]. For instance, the increases or decreases in specific frozen condensate amounts can dramatically modify the surface precipitation distribution of typhoons [2, 3]. When evaporative cooling of cloud droplets and melting of ice particles are removed in the microphysics scheme (MPS), it produces much weaker simulated downdrafts, resulting in much stronger simulated hurricanes [4]. Besides, the terminal velocity of graupel is another factor in tropical cyclone simulation [5]. Generally, the variations in hydrometeors, which is important for tropical cyclone simulations, are largely modulated by the MPSs in models [4, 6–12].
Currently, the MPSs are increasingly available within numerical models such as the Weather Research and Forecasting (ARWWRF) [13] model, ranging from single moment to double moment [12, 14]. Since the doublemoment schemes have an advantage of greater flexibility in size distribution parameters and density assumptions as compared with singlemoment approaches [15], they are often used in community research models in order to offer more faithful representations of the complex microphysical processes in a system [16, 17]. However, it is well established that, in some cases, doublemoment schemes do not offer improved results over singlemoment MPSs. For example, Van Weverberg et al. [18] evaluated mesoscale convective system simulations over the tropical western Pacific using three different MPSs and found that the performance of complex doublemoment schemes demonstrated no superiority over the simpler singlemoment schemes. In addition, considering that most of the existing cloud MPSs in mesoscale models are deduced from observations for midlatitude clouds, rather than from tropical cyclones [19], it remains unclear whether these MPSs related to hydrometeors are suitable and accurate enough in simulating or predicting cyclone. To improve tropical cyclone simulation, therefore, the performance of various MPSs should be compared and evaluated in a comprehensive way.
Until now, there are two ways for the assessment of MPSs based on satellite measurements. One is to directly evaluate modelsimulated hydrometeors against the satelliteretrieved hydrometeors (e.g., [20]). This approach seems straightforward in an explicit way, albeit probably large uncertainties in satellite hydrometeor retrievals. Another is to compare satelliteobserved radiances against those from model simulation. Particularly, satelliteobserved microwave radiance, which can penetrate and provide insights into the vertical profiles of most cyclone hydrometeors, has great advantages in characterizing cyclone hydrometeors and evaluating various MPSs in simulation of tropical cyclone hydrometeors [21]. In addition, microwave radiance signals are sensitive to the assumed bulk density of graupel and ice cloud particle sizes produced by MPSs [22]. Here the modelsimulated radiances are computed via radiative transfer models based on the simulated cloud, temperature, and humidity from numeric models (e.g., [19, 23]). In general, the radiancebased approach is superior in assessing hydrometeor output from models because the radiative transfer model is relatively more accurate in the simulation of radiances [24–29]. For a more comprehensive and reliable satellitebased assessment, we will combine the retrievalbased and radiancebased approach to assess the simulated typhoon hydrometeors in the present study.
The simulated hydrometeors vary significantly depending on the selected MPS in the models [19, 30]. More importantly, the performance of the current MPSs is still controversial in simulating hydrometeors of tropical cyclones. Therefore, the purpose of the present study is to assess the performance of different MPSs in simulation of hydrometeors in a supertyphoon case of Neoguri in 2014, using the measurements of newgeneration polar orbiting satellite of Fengyun3 (FY3) [31, 32]. As the first supertyphoon of 2014, Neoguri occurred during 3–10 July 2014 with a minimum sealevel pressure of 918 hPa and 1minute maximum wind speed of 72 m/s according to Joint Typhoon Warning Center (JTWC) reports. It was a large and powerful tropical cyclone, which directly struck Caroline Island, Guam, and Japan [33, 34]. The remaining of this article is organized as follows. In Section 2, the microwave brightness temperature (TB) observations from FY3 satellite with relevant radiance indices, the WRF model and experiment design for typhoon simulations, the planeparallel microwave radiative transfer model (MWRT) [35, 36] for calculating modelbased radiances, and the nonlinear Lucas–Kanade (LK) optical flow (OF) technique for comparisons are briefly described. The comparison results between observations and model simulations regarding the tracks, intensity, and surface precipitation rate of typhoon will be elaborated and discussed in Section 3. The retrievalbased and radiancebased evaluation of simulated typhoon hydrometeors are performed in Section 4, followed by discussions and conclusions in Sections 5 and 6, respectively.
2. Data and Methodology
2.1. FY3B Observation
As the second satellite in the FY3 series of the newgeneration polarorbit meteorological satellite of China, the FY3B satellite was launched on November 5, 2010, in an 831 km orbit with an inclination of 98.81° and a period of 101.49 min, and detailed descriptions can be found in previous studies [37]. As a passive microwave sensor onboard FY3B, the microwave radiometer imager (MWRI) takes advantage of coincident data at five frequencies with both vertical and horizontal polarization channels: 10.65, 18.7, 23.8, 36.5, and 89 GHz [37–39]. As shown in Table S1 of the supplementary material (SM), there are eight MWRI orbits in total (hereinafter referred to as 0703D, 0705A, 0705D, 0707A, 0707D, 0708A, 0708D, and 0709A) used to evaluate the hydrometeors of typhoon Neoguri in different evolution stages, including formation, developing, mature, and decaying. Moreover, based on the algorithm by Liu and Curry [40], the surface precipitation rate derived from MWRI was obtained to compare with the simulated surface precipitation rate. In addition, the simulated frozen precipitation hydrometeor (graupel and snow) path was compared to retrievals from the Microwave Humidity Sounder and MWRI [39]. Because raw brightness temperatures represent an integration of electromagnetic radiation information from emission and scattering by all of surface and atmospheric constituents, there is a considerable amount of ambiguity in directly interpreting raw brightness temperatures in terms of hydrometeor properties [17, 25, 41]. In order to establish the relationship between radiances and various hydrometeor more clearly, two particular radiance indices are used through quasilinear transformations of brightness temperatures on dualpolarization channels for specific frequency [41–43]. One is the emission index (EI), and the other is the scattering index (SI). Here EI is formulated aswhere represents polarization difference between vertically and horizontally polarized brightness temperatures at 18.7 GHz (∼19) given byand is the corresponding in the same scene but for clearsky situation. It can be noticed that EI would vary from 0 to 1, where small values indicate almost cloudfree field of view, while values near 1 denote opaque liquid clouds due to the strong emission effect at this frequency. As a result, EI is generally related to the content of liquid hydrometeors. On the contrary, SI can be expressed aswhere is the 89 GHz polarizationcorrected temperature defined by Spencer et al. [44] aswhere and are, respectively, vertically and horizontally polarized brightness temperatures at 89 GHz. Similar to , is that at clear sky. In contrast to EI, SI represents the volume scattering effect associated with frozen precipitation hydrometeors left. Thus, a large value of SI indicates the occurrence of abundant solid hydrometeors [25, 43].
2.2. Model Configuration
In order to reveal the difference of hydrometeor simulations among various MPSs, the WRF (Version 3.6.1) [13] model with twoway interactive nesting is used to perform numerical simulations for typhoon Neoguri.
Three multiplenested domains (Figure 1) were constructed, with horizontal grid spacings of 36, 12, and 4 km, with corresponding numbers of grid points: 243 × 264, 463 × 508, and 1042 × 1156, respectively. The vertical height top was set to 50 hPa, and 36 vertical sigma layers were used in the WRF. The two outermost domains (D01 and D02) used five physics schemes: the Rapid Radiative Transfer Model for General Circulation Models scheme (RRTMG) [45] for longwave radiation, the RRTMG scheme for shortwave radiation, the Revised MM5 Monin–Obukhov scheme [46] for surface physics, the Yonsei University scheme (YSU) [47] for the boundary layer parameterization, and the Tiedtke cumulus parameterization scheme [48, 49]. For the innermost domain (D03), the cumulus parameterization scheme was not used because that type of scheme is not designed for horizontal and vertical grid spacings smaller than about 5–10 km [50] as the typical scale of most convective clouds [8, 50].
The initial and boundary conditions used in the WRF model were interpolated from the National Centers for Environmental Prediction Final Operational Global Analysis data with 1° × 1°spatial resolution and 6 hour temporal resolution. These initial parameters were inserted at 1200 UTC each day from 2 to 8 July 2014 and then integrated forward for 36 hours. Since the cyclone hydrometeor information is not included in the initial conditions, lowlevel convergence in these simulations would take some time to spin up from the largescale circulation, and hence the models are not generally reliable for the first 12 h or so [51]. Therefore, only those outputs running after 12 hours are analyzed here.
2.3. Microphysics Schemes
For mesoscale modeling, most MPSs are bulk schemes that predict hydrometeor mass mixing ratios and/or number concentrations but assume size distributions based on empirical exponential or gamma distributions (e.g., [52–55]). In this paper, five different MPSs ranging from single moment to double moment were used in the WRF model, including the WRF SingleMoment 6Class Scheme (WSM6) [56], the Goddard scheme (GSFC) [53], the Thompson scheme (THOM) [54], the Morrison doublemoment scheme (MORR) [55], and the National Severe Storms Laboratory doublemoment scheme (NSSL) [57].
Generally, there are certain differences in hydrometeor schemes among these five MPSs as follows. Both the WSM6 and GSFC, which mainly considered cloud water and ice as being the same uniform size, assumed spherical particles of constant density with an exponential size distribution as described by Lin et al. [52] and Rutledge and Hobbs [58]. Differently, the intercept parameter (N_{ox}) is fixed for all precipitating classes in the GSFC scheme, while the snow intercept parameter varied with a function of temperature in the WSM6 scheme (Table 1), according to the observations by Houze et al. [59]. In the THOM scheme, it was considered that nonspherical snowflake shape with bulk density varies inversely with diameter based on observations [60]; hydrometeor size distribution is represented by a combination of two gamma functions [61], which depended upon observed relationships between predicted snow mass and air temperature to calculate the size distribution moments. In addition, the intercept parameter for graupel varies with the mass of supercooled liquid [16]. For the NSSL scheme [57], it added a frozen hydrometeorhail and employed double moments for considering a more robust treatment of the particlesize distribution (equations S2–S4 in SM). While hail generation occurs only with wet graupel, the hail category was not dominated by highdensity ice. Therefore, graupel particles in the NSSL scheme can range from frozen drops to lowdensity graupel, while there are only negligible amount of hail in the simulated typhoon of NSSL. Finally, for the MORR scheme, it took account for more complicated rain size distributions with different rates of rain evaporation in stratiform and convective regions. Besides, the MORR scheme has options to optimize simulations by accommodating the selection of ice nucleation methods and CCN spectra [62]. Further details of these MPSs are shown in Table 1 and in SM.
 
is the shape factor for the gamma function; N_{o} is the intercept parameter; NA indicates not available; Ra, Sn, Gr, and Ha stand for rain, snow, graupel, and hail, respectively. 
2.4. Microwave Radiative Transfer Simulations
To examine the capability of the various MPSs in replicating the microphysical cloud properties observed by the FY3B MWRI, the model output fields from the WRF including hydrometeor simulations were incorporated into the MWRT model [15, 35, 36]. As such, the modelequivalent MWRI radiance TB at 18.7 GHz and 89.0 GHz frequencies was obtained on model grids. Besides, the extinction and scattering properties of the hydrometeors were determined based on the Mie calculation. To ensure data consistency, the simulated TBs were further interpolated to the MWRI swath resolution for pixelbypixel comparison.
2.5. Nonlinear Lucas–Kanade Optical Flow Approach
The verification method of optical flow (OF) attempts to address whether two fields (e.g., the observed and simulated fields) agree with each other by inferring the mapping function. The mapping provides a visual summary of the “difference” between the two fields [62], which is a technique developed to assess the quality of simulated spatial fields. The notion of an OF is broad. Many variants take a purely algorithmic or nonparametric approach (e.g., [63]), capable of modeling a wide range of flow patterns [64]. Among others, one variant of the OF verification approach called Lucas–Kanade (LK)OF allows for analytic solutions in an objective way [65, 66]. LKOF has the merit of simultaneously estimating three components of errors for a given simulated variable, including its intensity, the direction, and distance of displacement [64]. The original LKOF formulation was linear; then, in order to provide more accurate estimates of the errors, Marzban and Sandgathe [64] extended it into a nonlinear formulation. Here, this nonlinear LKOF approach is employed to compare the observed and simulated radiance indices of EI and SI. Details of the formulas are provided in SM. On the contrary, when carrying out the comparison in the LKOF approach, both the observed and the simulated TBs were regridded to 0.2° × 0.2°resolution.
3. Typhoon Track, Intensity, and Surface Precipitation Simulation
3.1. Typhoon Tracks
The simulation of the typhoon track plays an important role in typhoon modeling studies, as it is the basis for simulating structures of typhoon surface precipitation and cloud hydrometeors. The track errors might be made by the inaccuracy of balanced initial conditions, absence of favorable environmental variables, and/or inadequacy in the other physical processes. Figure 2(a) shows that all five simulated tracks resembled the best track archived by JTWC in most of the typhoon life cycle, except for the formation stage. As illustrated in Table S2 of SM, the simulated typhoons by all the five schemes moved slightly faster than JTWC, which probably is caused by the different azimuthal organizations of convection and precipitation patterns [1, 67]. In Figure 2(b), though MORR produces a relatively better track simulation and has a minimum average error of 69.5 km and a minimum standard deviation of 42.2 km (more statistics can be shown in Table S3), it could be seen that the track simulation for typhoon Neoguri was little sensitive to the choice of MPS.
(a)
(b)
3.2. Typhoon Intensity
Figure 3 shows the time series of mean sea level pressure (MSLP) and maximum surface wind speed (MWS) for JTWC and simulations, respectively. Also, the peak intensity and corresponding peak time of them are shown in Table 2. For simulated MSLP, in all five schemes, the time reaching peak intensity lagged behind the observation (JTWC), and the typhoon intensity was overestimated after 1200 UTC 7 July when Neoguri began to decay. Relatively, the peak times of NSSL were much closer to the observed values, probably because of the lowsimulated translational speed mentioned above. For simulated MWS, all five schemes obviously overestimated MWS during the formation and early development stages (3–4 July). For rapid development and mature stages, WSM6, GSFC, and MORR slightly underestimated it, while NSSL and THOM perform better.
(a)
(b)
 
T denotes the time (UTC) with peak intensity for the typhoon. 
In addition, MSLP and MWS from the WSM6 and GSFC, respectively, are quite similar to each other. Such notable consistent results might be caused by fixing one or several parameters in the semiempirical description of PSD [68]. The aforementioned results, based on simulated track and intensity, show that the simulations based on doublemoment MPS (NSSL and MORR) were somewhat superior to those of the singlemoment or partial doublemoment MPS (WSM6, GSFC, and THOM). The doublemoment MPS offered more faithful representations of the complex microphysical processes in the given weather system [16]. Specially, it can be noted that NSSL more closely matched the JTWC estimates.
3.3. Surface Precipitation
The spatial distribution of simulated surface precipitation is shown in Figure 4, along with that obtained from the FY3BMWRI overpass time instances during the different stages of the typhoon life cycle. In general, the simulated precipitation amounts could capture the trends represented in the observed precipitation. For the precipitation pattern, however, the observed precipitation was mainly concentrated in domains within a radius of about 400 km from the cyclone center during the rapid intensification to the early decaying stage (05/07–08/07), while the associated spatial extent of simulated surface precipitation was slightly larger, and also more isolated precipitation clusters can be found in the all the simulations, especially for WSM6 and GSFC.
Moreover, it is notable that the simulated typhoon exhibited systematically wetter conditions over the strong convective activity areas of the inner rainband than observed, whereas for the outer spiral rainbands, the simulation was noticeably drier, and the diversities in the outlying spiral rainband regions were significant among different MPSs. Given that the strength of northwest Pacific cyclone intensity was remarkably well correlated with surface precipitation within the inner cores of the cyclone [69], it may have led to a simulated typhoon stronger (lower MSLP) than the observed one (JTWC). By comparison, the results of NSSL, with a relatively weaker and smaller distribution of graupel hydrometeor path at the surface precipitation extension (Figure S1 of SM), were generally closer to the observations than those of the other schemes. As showed in Figure S1 of SM, we also note that the spatial distribution pattern of graupel for individual simulation times closely resembled the intense precipitation location, implying that graupel hydrometeors had higher density and faster fall speeds to reduce their losses induced by evaporation and to shorten the residence time of the hydrometeor and thus modulated the surface precipitation [1, 70]. In general, the NSSL scheme’s propensity for graupel production likely explains the ability of this scheme to produce the spatial extent of the precipitation regions.
4. Evaluation of Simulated Hydrometeors in Typhoon Neoguri
In this section, one typical case in the strong stage of typhoon Neoguri on 05/07/2014, fully captured by FY3BMWRI, as shown in Figure 4(g), was chosen for further evaluating the simulated hydrometeors with different MPSs by combining the retrievalbased and radiancebased approaches.
4.1. RetrievalBased Approach and Comparison of Simulated Hydrometeors
Considering the difficulty in retrieving liquid water content for cyclone due to the complexity of microwave signals about liquid precipitation particles, only simulated frozen precipitation path including graupel and snow as ice water path (IWP) is directly evaluated here (Figure 5). It can be seen that simulated IWPs in all schemes are higher than that by retrieval in most areas, especially in outside areas (e.g., spiral cloud belts). Comparatively, the IWP spatial pattern and values for NSSL and WSM6 are generally closer to those for FY3B MWRI, vice versa, for THOM. Moreover, as indicated by white circles in Figure 5, IWP statistics in 1°radius innercore zone, 1°radius ∼2.5°radius middleannular zone, and 2.5°radius ∼4°radius outsideannular zone are given in Table 3, respectively. Obviously, the retrieved IWPs within innercore and middleannular zones (1554 and 1002 g m^{−2}) are obviously higher than those within the outsideannular zone (405 g m^{−2}). However, for all the MPSs, no such apparent IWP decrease from the inside zone to outside zone can be found. The IWP values for NSSL are closer to those retrieved by FY3B MWRI, especially within innercore and outsideannular zones, probably due to less graupel contents (refer to Figure 5 and Table S4 in SM). A potential explanation is that NSSL uses adaptive sedimentation to allow some size sorting and to prevent spurious large particles, which can arise in the larger precipitation hydrometeors [16, 57]. Over all three areas, THOM produces the largest amount of IWPs, among all the schemes, mainly in the form of snow (Figures 6 and S2 in SM; Tables 3 and S4). Wheatley et al. [16] found that falling speed made a large contribution on the sensitivity to precipitation hydrometeors with respect to the corresponding MPS, and the disparity in the snowfall speeds between the THOM (slower) and NSSL (faster) MPSs may account for the different evolutions of snow fields.

(a)
(b)
(c)
(d)
(e)
In general, compared with retrieved IWP, all simulated IWPs of five MPSs are systematically overestimated over all three areas, except for the NSSL over the innercore zone. These noticeable distinctions in the simulated hydrometeors of various MPSs are highlighted by the vertical cross sections of timedomain ensemble means of cloud water, rainwater, cloud ice, snow, and graupel hydrometeor contents, within a radius of 400 km from the typhoon center, from 0000 to 2400 UTC 5 July 2014 (Figure 6). More detailed comparison of hydrometeors among five MPSs can be found below.
The evident differences in the five MPSs are the distributions of liquid hydrometeors (cloud water and rainwater) at high altitudes, especially above 5 km (i.e., supercooled water). In the NSSL simulations, supercooled water of 0.02 g m^{−3} is at about the −40°C level, i.e., ∼12 km horizontally extended over 300 km radius wide in the convective region (Figure 6(e)), whereas the magnitudes of supercooled water contents were relatively smaller in the other MPSs. Based on the combination of radar, passive microwave, and lightning observations, Cecil and Zipser [71] proposed that supercooled water occurs preferentially in outer rainbands, compared to other tropical oceanic precipitation. Meanwhile, Cecil and Zipser [71] and Reinhart et al. [72] also suggested that supercooled water occurs at a height of about 11 to 12 km (∼−40°C) in the Tropical Cyclone cloud band. These observations indicate that NSSL can reproduce a more realistic vertical structure of supercooled water than do the other MPSs. In addition, the rainwater contents of GSFC, WSM6, THOM, and MORR were all similar, while the content for the NSSL was slightly lower than the others.
For cloud ice, the simulation using NSSL/MORR MPS was characterized by more cloud ice in the upper levels. It is evident that cloud ice from the simulations of NSSL and MORR schemes is clearly characterized at a higher altitude (>12 km).While only negligible amounts of cloud ice (<0.02 g m^{−3}) were present in the mid and upper levels of the THOM simulation, indicating that THOMbased simulation may generate less cloud ice. According to previous studies [73–75], this is because THOM did not implement the cloud ice formulation scheme from Meyers et al. [73] and thus resulted in producing less cloud ice contents than MORR although both MORR and THOM employed the ice nucleation approach of Cooper [76] and freeze cloud liquid water into cloud ice from Bigg [77]. In addition, the assumed ice nuclei number concentrations of WSM6 and GSFC usually decreased with the increase of altitude at the vertical direction; that is, more (less) ice can be produced at lower (higher) altitudes from WSM6 and GSFC schemes than those from NSSL and MORR schemes [75, 78].
For snow, in the THOM simulations, aggressive snow content of >0.4 g m^{−3} extended more than 200 km at 5–14 km altitudes in the convective region (Figure 6(c)). For the THOM MPS, the unique snow PSD assumptions were used, which mainly constrained by in situ measurements under winter storm conditions, leading to smaller particle sizes, and the using massdiameter relationship mainly resulted from the fractallike assumption of snow aggregates. Field et al. [61] indicated that these assumptions in THOM might be particularly representative for wintertime precipitation systems rather than typhoons. Instead, the simulated snow content of WSM6 was relatively lower than other four MPSs, possibly because the WSM6 employed a temperaturedependent form and a sink to graupel based on a tunable critical value, causing the discrepancies in snow microphysical processes.
A less compact region of graupel content of 0.05–0.3 g m^{−3} was presented in the simulation using the NSSL MPS (Figure 6(e)), while overwhelming values of graupel (0.05–0.6 g m^{−3}) are presented from 5 to >15 km in the GSFC simulation. The graupel contents of MORR and WSM6 were also apparently greater than those of the NSSL (Figures 6(a) and 6(d)). A possible explanation is that the NSSL scheme’s propensity for graupel allowed graupel particles to range from frozen drops to lowdensity graupel (300–900 kg m^{−3}), whereas the graupel particle densities in all of the WSM6, GSFC, and MORR MPSs were assumed to be constant (Table 1), leading to a disparity in the graupel falling velocity. In addition, the terminal velocity of graupel plays a large part in this sensitivity within tropical cyclone simulation [5], which may also account for the different evolutions of graupel fields. Moreover, graupel stayed at a conspicuous lowest content of ∼0.05–0.2 g m^{−3} in the THOM simulation. Note that NSSL produced negligible hail at about 6 km (∼0.01 g m^{−3}), and the hail generation occurred only with wet graupel, so that the hail was not merely dominated by highdensity ice.
Furthermore, regional average results (Figure S2) show that, above the melting layer (around 5 km), the total cloud hydrometeor profiles showed more pronounced sensitivity to different MPSs, which are also suggested by Rajeevan et al. [79] and Wheatley et al. [16].
4.2. RadianceBased Approach
To further examine the suitability of MPSs in simulating various typhoon hydrometeors, the nonlinear LKOF approach was applied to radiance indices of SI and EI from the observations and simulations (Figures 7 and 8) as the radiancebased evaluation on liquid and frozen hydrometeors, respectively. Their corresponding statistics are summarized in Table 4.
 
Mode and mean are two kinds of metrics for “averages.” 
In terms of EI (Figure 7), the typhoon cloud rain coverage ratio of 78.5% in NSSL, higher than MORR (75.6%), WSM6 (54.8%), GSFC (54.9%), and THOM (68.4%), is the closest to the observed value (85.5%). Moreover, for all five MPSs, the EI spatial distribution was simulated slightly southwest (Figures 7(a)–7(e), the overlaid OF field). With regard to intensity, the simulated EI was generally lower in the outer spiral rainbands and of greater magnitude in the nearcore region than observed EI.
The frequency distributions of the three components of error further demonstrate the difference between simulations and observations. Intensity errors (Figure 7(l)) clearly illustrate that the simulated EI was commonly lower than the observed values, with the same positive mode in all MPSs (Table 4). This result is mainly associated with the underestimation of the high EI areas over the main part of the typhoon, especially for the western area of the exterior spiral cloud band. The GSFC had the minimum average error of 0.05, while THOM, MORR, and NSSL had the maximum of 0.09. The histogram of the EI displacement errors (Figure 7(m)) demonstrates that the typical displacement error was about 10∼20 km, where the largest displacement errors occurred in the boundary of the exterior spiral cloud band (Figures S3(k)–S3(o) in SM). In terms of angle error, the five schemes also had similar modes in the histogram (Figure 7(n)) where all are in the first quadrant (corresponding to the southwest bias of the spatial pattern of the angle errors) ranging from 30° to 60°. In general, these similar systematic biases of EI indicate that all five MPSs underestimated liquid hydrometeors.
With regard to SI (Figure 8), one of the most obvious features is that the simulated SI intensity over the southwest part of the typhoon was significantly higher than the observed intensity, which resulted in the corresponding negative average intensity error, though a slightly positive bias is on the outer cloud edge of the western part of the typhoon (Figures 8(f)–8(j)). Among five MPSs, the NSSL has the smallest mean intensity error of −0.03, while THOM has the largest value of −0.20. In addition, it can be noted that the displacement and directional errors of the various schemes were similar, about 10∼20 km to the southwest (∼45°) of the observations, which is also consistent with the displacement and direction error of EI. These errors may be caused by a model system issue, such as imperfectly balanced initial and boundary conditions or a deficiency in the other models’ physical processes, with the placement of the typhoon.
In order to understand the assignment between liquid and frozen precipitation hydrometeors in typhoon, the joint histograms of SI versus EI were calculated, as shown in Figure 9. Overall, the EI values increased with increasing SI from observations (Figure 9(a)), which implies the fact that liquid hydrometeor would increase with frozen hydrometeors such as graupel and snow in actual typhoons. However, the results from the five MPSs generally disagreed with the observed trends, showing that a slow increase of EI as the SI increased (Figures 9(b) and 9(c)). It indicates that there is a problem in the simultaneous growth process of the simulated frozen and liquid hydrometeors. On the contrary, the observed mode was bounded with the SI values 0.8∼0.9, while the simulated mode extended to higher SI and EI, especially for WSM6, GSFC, and MORR with SI values 1.5. For a given EI value, moreover, the simulated SI value was generally 0.5 greater than the observations. This feature is likely related to the larger amount of frozen precipitation hydrometeors in modeling, especially for large frozen hydrometeor of graupel as the main contributor to SI (Figures 5 and 6; Tables 3 and S4).
In general, the results from NSSL were closer to the observations than those from the other schemes at the SI extension, consistent with the less graupel in NSSL mentioned in Tables 3 and S4. Moreover, the EI values with high occurrence frequency in all five schemes were appreciably (>∼0.15) lower than the observations of MWRI, suggesting that the simulated liquid hydrometeors were less than the actuals, especially for rainwater having strong emission effect. In addition, the highvalue region in simulated joint probability density distribution is more inclined to the SI axis, whereas the observational results tend to the EI axis, suggesting large uncertainties of MPSs in the matching ratio of the simulated frozenliquid hydrometeors. Even so, it is distinct that the NSSL reproduced a more realistic hydrometeor distribution than did the other four schemes.
5. Discussion
To illustrate the contribution of the specific hydrometeor to the radiance indices, the sensitivity was determined by comparing the difference of calculated SI and EI between all with and without this hydrometeor. Note that a hydrometeor species was removed only in the radiative transfer calculations and not in the actual simulations. This sensitivity experiments will show the impact of each hydrometeor species on the simulationderived radiance indices and will enable us to visually and accurately evaluate the deviation contribution of various hydrometeor species in different MPSs.
Taking the stronger stage of typhoon Neoguri on 05/07/2014 as the same example as previous work, the results of the sensitivity experiment were given in Figure 10. It is clear that cloud ice has negligible impact on SI, while graupel and snow have appreciable impacts on it. When snow is removed, the maximum value of SI in the joint histograms of NSSL and THOM (Figure 10(o)(z1)) decreases below 0.9, which is slightly smaller the observed SI value, but the SI ranges of the other schemes (SI > 1.2) are still obviously greater than the observed values, especially for MORR and GSFC. Note that calculations with no snow separated out the effect of graupel in SI, showing that the graupel contents derived from the WSM6, GSFC, and MORR were greatly exaggerated. In addition, SI seems to be somewhat affected by snow. When graupel is removed, the SI extensions of all the schemes are still larger than the observed SI value, indicating that the simulated snow contents are exaggerated (Figure 10).
For all five MPSs, the simulated joint histograms show that cloud water had a negligible impact on EI, which indicates that EI is mainly affected by the emission of rainwater as mentioned above. In this case with the absence of cloud water, the maximum value range of EI was still slightly higher than the observed value, but the EI value at the maximum probability was lower than the observed value, indicating that the five schemes somewhat overestimated strong rainwater (such as in the nearcore region and in strong convective activity areas), while it is possible to underestimate the content of moderate rainwater (in the less active outer spiral rainbands). In addition, supercooled water and graupel have competing effects on the upwelling brightness temperature. While graupel reduces the brightness temperature by scattering upwelling radiation, emissions from supercooled water droplets increase the brightness temperature [71, 80, 81]. Thus, supercooled water masks the scattering signal that results from graupel at the same vertical level. Hong et al. [82] pointed out that the channels at 89 GHz are strongly positive (increasing the brightness temperature, that is, decreasing SI), sensitive to variations in the liquid water content above 5 km (supercooled water). The calculations with no cloud water had little impact on SI for WSM6, GSFC, THOM, and MORR, indicating that there was too little supercooled water in the above four schemes, so that SI was not significantly decreased. Therefore, too little supercooled water may be another reason why the SI extensions of these four schemes were greater than the observed counterparts. However, the spread in the NSSL joint histograms with no cloud water (Figure 10(z3)) is more expansive (∼0.3) along the SIaxis than that with cloud water (Figure 10(y)). This indicates that NSSL, to some extent, could truly generate larger cloud water content and a deep layer of supercooled water extending up to about 11 km (as shown in Figures 6 and S2 in SM), consistent with the results of Cecil and Zipser [71] and Reinhart et al. [72].
Note that the tropical cyclone first and foremost is generated by the complex interactions of the largescale dynamics and cumulus convection [83–85]. To date, a successful forecast of the tropical cyclones strongly depends on how successfully a model can represent the aspects related to these interaction processes. However, the selection of appropriate MPSs responsible for hydrometeor variations makes a difference as well for the accurate simulation of tropical cyclone. Actually, it has been well established that the released latent heat (convective heating) caused by phase transition of cloud hydrometeors is the main source of energy in tropical cyclone [86], and different types of hydrometeors lead to quite different structures and intensities of tropical cyclones through their interactions with radiation [8, 87, 88]. As shown here, even with the same dynamical framework and other physics schemes, the simulations using NSSL as the representative of the doublemoment MPS are much closer to the observations than other MPSs, in term of both hydrometeor distributions and typhoon intensity variations. This should not be regarded as coincidence since the generated hydrometeors in the cloud microphysics process would exert influences on typhoon intensity that is related to dynamics through changing latent heat and radiation.
6. Conclusions
In this study, based on a case of typhoon Neoguri, the ability of typhoon hydrometeor simulation about five MPSs in the WRF Version 3.6 model, including GSFC, WSM6, THOM, MORR, and NSSL, was directly and indirectly evaluated through comparing with observations from FY3B satellite. The following conclusions can be drawn:(1)Comparison with JTWC reports indicates that the simulated track is not very sensitive to the chosen MPSs. Among these schemes, the MORR produces a minimal average track error of 69.5 km. Moreover, NSSL tends to reproduce more actually intensity variations in the whole typhoon life cycle. On the contrary, the simulated typhoon surface precipitation exhibits systematically wetter conditions over the strong convective activity areas of the inner rainband, whereas it exhibits noticeably drier conditions for the stratiform cloud rainband. In comparison, the surface precipitation pattern of NSSL is closer to the observed pattern.(2)In terms of satellite retrievalbased evaluation, simulated IWPs by all MPSs are higher than retrievals in most areas. Notably, the spatial pattern of IWP for NSSL and WSM6 is generally closer to those from FY3B observations. Specially, within the 1°radius innercore zone and 2.5°∼4° outsideannular zones, the IWPs derived from NSSL with less graupel are closer to the retrieved IWP. On the contrary, over almost all areas, THOM produces the largest amount of IWPs, which are mainly in the form of snow.(3)For the radiancebased evaluation, the horizontal distributions of simulated EI and SI were statistically verified using the nonlinear LKOF approach. The results show that both the simulated EI and SI had significant intensity bias, especially in strong convective areas. This implies that the simulated liquid hydrometeor is generally lower in the outer spiral rainbands than the one as observed from satellite. However, the simulated frozen precipitation hydrometeor is obviously greater over the southwest part of the typhoon. In addition, similar displacement and angular errors (a systematic southwest bias of approximately about 10∼20 km) of EI or SI among the five MPSs might indicate a serious model system issue.(4)The SIEI joint histograms show that NSSL reproduces a more realistic matching proportion of frozen and liquid hydrometeors than any of the other four schemes, especially for supercooled cloud water.
Overall, for this case of Neoguri, the doublemoment MPS of NSSL, allowing for greater flexibility in the size distribution and the density of graupel particles with more complex microphysical processes, could offer more reliable representations for the simulation of hydrometeors, as well as the evolution of typhoon.
In summary, our work advances our understanding of the performance of various MPSs of cloud in simulating the hydrometeors of tropical cyclones. Admittedly, the heightrevolved structures of hydrometeors are not evaluated at all, which merits further investigation in the future by combining passive and active satellite observations, including CALIPSO, Cloudsat, TRMM, and GPM [89, 90].
Data Availability
The MWRI data were provided by FENGYUN Satellite Data Center from http://data.nsmc.org.cn/portalsite/default.aspx. The JTWC data were downloaded from https://www.usno.navy.mil/NOOC/nmfcph/RSS/jtwc/best_tracks/. The National Center for Atmospheric Research’s Mesoscale and Microscale Meteorology Division provided the WRF model available at http://www2.mmm.ucar.edu/wrf/users/. The MWRT model was accessed from http://cirrus.met.fsu.edu/research/mwrt.html. All data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to thank Dr. Ted Mansell at the National Severe Storms Laboratory for his great help on understanding the NSSL scheme in the WRF. FY3B data were obtained from FENGYUN Satellite Data Center. This work was supported by the Ministry of Science and Technology of China (grant nos. 2018YFC1507401 and 2017YFC1501402), the National Natural Science Foundation of China (grant nos. 41875024, 41830104, and 41375030), and the Fundamental Research Funds for the Central Universities. This work was also sponsored by the Jiangsu Provincial 2011 program (Collaborative Innovation Center of Climate Change).
Supplementary Materials
Table S1: eight orbit datasets of FY3B, which were obtained from the FY3BMWRI 83 overpass during the different stages of typhoon Neoguri’s life cycle. Table S2: statistics of typhoons translational speed for the model simulations and the best track (JTWC), during 0000 UTC 3 July to 0000 UTC 10 July 2014 at 6 h intervals. Table S3: validation statistics between simulated tracks and best track (JTWC). Table S4: statistics of the domainaveraged hydrometeor path (g/m2) of cloud water (Cw), cloud ice (Ci), rain water (Ra), snow (Sn), graupel (Gr), and hail (Ha), based on matched instances of 0705A, which are derived from five MPSs. Total is total vertically integrated hydrometeor paths. Figure S1: horizontal distribution of modelsimulated graupel hydrometeor path (g/m2) produced with the super typhoon Neoguri, using the MPSs of WSM6 (first column), GSFC (second column), THOM (third column), MORR (fourth column), and NSSL (fifth column). Note that the first to fourth rows correspond to the time of FY3B overpasses of 0703D, 0705A, 0707A, and 0708D. Figure S2: domainaveraged vertical profile of hydrometeor content (g/m3) for the simulations matched 0705A scene captured by FY3B over domains within a radius of 400 km from the cyclone center (the typhoon’s surface minimum pressure), for (a) cloud ice content, (b) snow water content, (c) graupel water content, (d) cloud water content, and (e) rain water content, respectively. Figure S3: the OF ﬁeld and the three components of error in the five simulated ﬁelds. (a)∼(e) The spatial distribution of the emission index (EI) ﬁeld: (a) WSM6, (b) GSFC, (c) THOM, (d) MORR, and (e) NSSL; the overlaid OF ﬁeld (arrows) mapping the simulated to the observed SI ﬁeld (u). The spatial distribution of intensity error (f) WSM6, (g) GSFC, (h) THOM, (i) MORR, and (j) NSSL; the spatial distribution of displacement error (k) WSM6, (l) GSFC, (m) THOM, (n) MORR, and (o) NSSL; the spatial distribution of angular error (p) WSM6, (q) GSFC, (r) THOM, (s) MORR, and (t) NSSL; and histogram of the three components of the simulation error of (v) intensity error (AA), (w) displacement error (DispE), and (x) angular error (see text for details of Figure S4. As in Figure S3, but for the scattering index (SI)). (Supplementary Materials)
References
 S. Pattnaik, C. Inglish, and T. N. Krishnamurti, “Influence of rainrate initialization, cloud microphysics, and cloud torques on hurricane intensity,” Monthly Weather Review, vol. 139, no. 2, pp. 627–649, 2011. View at: Publisher Site  Google Scholar
 R. Rogers, “Convectivescale structure and evolution during a highresolution simulation of tropical cyclone rapid intensification,” Journal of the Atmospheric Sciences, vol. 67, no. 1, pp. 44–70, 2010. View at: Publisher Site  Google Scholar
 J. L. Vigh and W. H. Schubert, “Rapid development of the tropical cyclone warm core,” Journal of the Atmospheric Sciences, vol. 66, no. 11, pp. 3335–3350, 2009. View at: Publisher Site  Google Scholar
 T. Zhu and D.L. Zhang, “Numerical simulation of Hurricane Bonnie (1998). Part II: sensitivity to varying cloud microphysical processes,” Journal of the Atmospheric Sciences, vol. 63, no. 1, pp. 109–126, 2006. View at: Publisher Site  Google Scholar
 C. N. Franklin, G. J. Holland, and P. T. May, “Sensitivity of tropical cyclone rainbands to icephase microphysics,” Monthly Weather Review, vol. 133, no. 8, pp. 2473–2493, 2005. View at: Publisher Site  Google Scholar
 S. J. Lord, H. E. Willoughby, and J. M. Piotrowicz, “Role of a parameterized icephase microphysics in an axisymmetric, nonhydrostatic tropical cyclone model,” Journal of the Atmospheric Sciences, vol. 41, no. 19, pp. 2836–2848, 1984. View at: Publisher Site  Google Scholar
 H. E. Willoughby, H.L. Jin, S. J. Lord, and J. M. Piotrowicz, “Hurricane structure and evolution as simulated by an axisymmetric, nonhydrostatic numerical model,” Journal of the Atmospheric Sciences, vol. 41, no. 7, pp. 1169–1186, 1984. View at: Publisher Site  Google Scholar
 Y. Piotrowicz, “An explicit simulation of tropical cyclones with a triply nested movable mesh primitive equation model: TCM3. Part II: model refinements and sensitivity to cloud microphysics parameterization,” Monthly Weather Review, vol. 130, no. 12, pp. 3022–3036, 2002. View at: Publisher Site  Google Scholar
 Z. Pu, X. Li, C. S. Velden, S. D. Aberson, and W. T. Liu, “The impact of aircraft dropsonde and satellite wind data on numerical simulations of two landfalling tropical storms during the tropical cloud systems and processes experiment,” Weather and Forecasting, vol. 23, no. 1, pp. 62–79, 2008. View at: Publisher Site  Google Scholar
 J. Li, G. Wang, W. Lin, Q. He, Y. Feng, and J. Mao, “Cloudscale simulation study of Typhoon Hagupit (2008) Part I: microphysical processes of the inner core and threedimensional structure of the latent heat budget,” Atmospheric Research, vol. 120121, pp. 170–180, 2013. View at: Publisher Site  Google Scholar
 W.K. Tao, J. J. Shi, S. S. Chen et al., “The impact of microphysical schemes on hurricane intensity and track,” AsiaPacific Journal of Atmospheric Sciences, vol. 47, no. 1, pp. 1–16, 2011. View at: Publisher Site  Google Scholar
 I. Pytharoulis, S. Kartsios, I. Tegoulias et al., “Sensitivity of a mediterranean tropicallike cyclone to physical parameterizations,” Atmosphere, vol. 9, no. 11, p. 436, 2018. View at: Publisher Site  Google Scholar
 W. C. Skamarock, J. B. Klemp, J. Dudhia et al., “A Description of the Advanced Research WRF Version 3. NCAR Technical Note,” NCAR/TN475+STR, Mesoscale and Microscale Meteorology Division, National Center for Atmospheric Research, Boulder, CO, USA, 2008. View at: Google Scholar
 A. L. Molthan and B. A. Colle, “Comparisons of single and doublemoment microphysics schemes in the simulation of a synopticscale snowfall event,” Monthly Weather Review, vol. 140, no. 9, pp. 2982–3002, 2012. View at: Publisher Site  Google Scholar
 W. Roh, M. Satoh, and T. Nasuno, “Improvement of a cloud microphysics scheme for a global nonhydrostatic model using TRMM and a satellite simulator,” Journal of the Atmospheric Sciences, vol. 74, no. 1, pp. 167–184, 2017. View at: Publisher Site  Google Scholar
 D. M. Wheatley, N. Yussouf, and D. J. Stensrud, “Ensemble Kalman filter analyses and forecasts of a severe mesoscale convective system using different choices of microphysics schemes,” Monthly Weather Review, vol. 142, no. 9, pp. 3243–3263, 2014. View at: Publisher Site  Google Scholar
 E.K. Seo and M. I. Biggerstaff, “Impact of cloud model microphysics on passive microwave retrievals of cloud properties. Part II: uncertainty in rain, hydrometeor structure, and latent heating retrievals,” Journal of Applied Meteorology and Climatology, vol. 45, no. 7, pp. 955–972, 2006. View at: Publisher Site  Google Scholar
 K. Van Weverberg, A. M. Vogelmann, W. Lin et al., “The role of cloud microphysics parameterization in the simulation of mesoscale convective system clouds and precipitation in the tropical western Pacific,” Journal of the Atmospheric Sciences, vol. 70, no. 4, pp. 1104–1128, 2013. View at: Publisher Site  Google Scholar
 G. M. McFarquhar, H. Zhang, G. Heymsfield et al., “Factors affecting the evolution of Hurricane Erin (2001) and the distributions of hydrometeors: role of microphysical processes,” Journal of the Atmospheric Sciences, vol. 63, no. 1, pp. 127–150, 2006. View at: Publisher Site  Google Scholar
 D. E. Waliser, J. L. F. Li, C. P. Woods et al., “Cloud ice: a climate model challenge with signs and expectations of progress,” Journal of Geophysical Research: Atmospheres, vol. 114, no. D8, 2009. View at: Publisher Site  Google Scholar
 V. S. Galligani, D. Wang, M. Alvarez Imaz, P. Salio, and C. Prigent, “Analysis and evaluation of WRF microphysical schemes for deep moist convection over southeastern South America (SESA) using microwave satellite observations and radiative transfer simulations,” Atmospheric Measurement Techniques, vol. 10, no. 10, pp. 3627–3649, 2017. View at: Publisher Site  Google Scholar
 S. B. Sieron, E. E. Clothiaux, F. Zhang, Y. Lu, and J. A. Otkin, “Comparison of using distributionspecific versus effective radius methods for hydrometeor singlescattering properties for allsky microwave satellite radiance simulations with different microphysics parameterization schemes,” Journal of Geophysical Research: Atmospheres, vol. 122, no. 13, pp. 7027–7046, 2017. View at: Publisher Site  Google Scholar
 F. Chevallier and P. Bauer, “Model rain and clouds over oceans: comparison with SSM/I observations,” Monthly Weather Review, vol. 131, no. 7, pp. 1240–1255, 2003. View at: Publisher Site  Google Scholar
 F. YunFei, H. ZhiWei, L. TianYi, S. ZhongPing, and W. Yu, “Evaluation of WRF model hydrometeors based on TMI observations using an indirect method,” Atmospheric and Oceanic Science Letters, vol. 5, no. 1, pp. 32–37, 2012. View at: Publisher Site  Google Scholar
 M. Han, S. A. Braun, W. S. Olson, P. O. G. Persson, and J.W. Bao, “Application of TRMM PR and TMI measurements to assess cloud microphysical schemes in the MM5 for a winter storm,” Journal of Applied Meteorology and Climatology, vol. 49, no. 6, pp. 1129–1148, 2010. View at: Publisher Site  Google Scholar
 M. Han, S. A. Braun, T. Matsui, and T. Iguchi, “Comparisons of bin and bulk microphysics schemes in simulations of topographic winter precipitation with radar and radiometer measurements,” Quarterly Journal of the Royal Meteorological Society, vol. 144, no. 715, pp. 1926–1946, 2018. View at: Publisher Site  Google Scholar
 T. Iguchi, J. Santanello, J. J. Shi et al., “Introducing multisensor satellite radiancebased evaluation for regional Earth System modeling,” Journal of Geophysical Research: Atmospheres, vol. 119, no. 13, pp. 8450–8475, 2014. View at: Publisher Site  Google Scholar
 X. Li, W.K. Tao, T. Matsui et al., “Improving a spectral bin microphysical scheme using TRMM satellite observations,” Quarterly Journal of the Royal Meteorological Society, vol. 136, no. 647, pp. 382–399, 2010. View at: Google Scholar
 T. Matsui, X. Zeng, W.K. Tao, H. Masunaga, W. S. Olson, and S. Lang, “Evaluation of longterm cloudresolving model simulations using satellite radiance observations and multifrequency satellite simulators,” Journal of Atmospheric and Oceanic Technology, vol. 26, no. 7, pp. 1261–1274, 2009. View at: Publisher Site  Google Scholar
 Y. Wang, C. N. Long, L. Leung et al., “Evaluating regional cloud‐permitting simulations of the WRF model for the tropical warm pool international cloud experiment (TWP‐ICE), Darwin, 2006,” Journal of Geophysical Research: Atmospheres, vol. 114, no. D21, 2009. View at: Publisher Site  Google Scholar
 J. Yang, P. Zhang, N. Lu, Z. Yang, J. Shi, and C. Dong, “Improvements on global meteorological observations from the current Fengyun 3 satellites and beyond,” International Journal of Digital Earth, vol. 5, no. 3, pp. 251–265, 2012. View at: Publisher Site  Google Scholar
 Z. Yang, N. Lu, J. Shi, P. Zhang, C. Dong, and J. Yang, “Overview of FY3 payload and ground application system,” IEEE Transactions on Geoscience and Remote Sensing, vol. 50, no. 12, pp. 4846–4853, 2012. View at: Publisher Site  Google Scholar
 Y. Han, X. Zou, and F. Weng, “Cloud and precipitation features of Super Typhoon Neoguri revealed from dual oxygen absorption band sounding instruments on board FengYun3C satellite,” Geophysical Research Letters, vol. 42, no. 3, pp. 916–924, 2015. View at: Publisher Site  Google Scholar
 Y.A. Liou, J.C. Liu, C. P. Liu, and C.C. Liu, “Seasondependent distributions and profiles of seven supertyphoons (2014) in the northwestern pacific ocean from satellite cloud images,” IEEE Transactions on Geoscience and Remote Sensing, vol. 56, no. 5, pp. 2949–2957, 2018. View at: Publisher Site  Google Scholar
 G. Liu, “A fast and accurate model for microwave radiance calculations,” Journal of the Meteorological Society of Japan. Ser. II, vol. 76, no. 2, pp. 335–343, 1998. View at: Publisher Site  Google Scholar
 G. Liu, “Approximation of single scattering properties of ice and snow particles for high microwave frequencies,” Journal of the Atmospheric Sciences, vol. 61, no. 20, pp. 2441–2456, 2004. View at: Publisher Site  Google Scholar
 C. Dong, J. Yang, W. Zhang et al., “An overview of a new Chinese weather satellite FY3A,” Bulletin of the American Meteorological Society, vol. 90, no. 10, pp. 1531–1544, 2009. View at: Publisher Site  Google Scholar
 H. Yang, F. Weng, L. Lv et al., “The FengYun3 microwave radiation imager onorbit verification,” IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 11, pp. 4552–4560, 2011. View at: Publisher Site  Google Scholar
 Y. Wang, Y. Fu, X. Fang et al., “Estimating ice water path in tropical cyclones with multispectral microwave data from the FY3B satellite,” IEEE Transactions on Geoscience and Remote Sensing, vol. 52, no. 9, pp. 5548–5557, 2014. View at: Publisher Site  Google Scholar
 G. Liu and J. A. Curry, “Retrieval of precipitation from satellite microwave measurement using both emission and scattering,” Journal of Geophysical Research: Atmospheres, vol. 97, no. D9, pp. 9959–9974, 1992. View at: Publisher Site  Google Scholar
 G. Petty, “Physical retrievals of overocean rain rate from multichannel microwave imagery. Part I: theoretical characteristics of normalized polarization and scattering indices,” Meteorology and Atmospheric Physics, vol. 54, no. 14, pp. 79–99, 1994. View at: Publisher Site  Google Scholar
 G. W. Petty and K. B. Katsaros, “Precipitation observed over the South China Sea by the Nimbus7 scanning multichannel microwave radiometer during winter MONEX,” Journal of Applied Meteorology, vol. 29, no. 4, pp. 273–287, 1990. View at: Publisher Site  Google Scholar
 G. Liu, J. A. Curry, and R.S. Sheu, “Classification of clouds over the western equatorial Pacific Ocean using combined infrared and microwave satellite data,” Journal of Geophysical Research, vol. 100, no. D7, pp. 13811–13826, 1995. View at: Publisher Site  Google Scholar
 R. W. Spencer, H. M. Goodman, and R. E. Hood, “Precipitation retrieval over land and ocean with the SSM/I: identification and characteristics of the scattering signal,” Journal of Atmospheric and Oceanic Technology, vol. 6, no. 2, pp. 254–273, 1989. View at: Publisher Site  Google Scholar
 M. J. Iacono, J. S. Delamere, E. J. Mlawer et al., “Radiative forcing by long‐lived greenhouse gases: calculations with the AER radiative transfer models,” Journal of Geophysical Research: Atmospheres, vol. 113, no. D13, 2008. View at: Publisher Site  Google Scholar
 P. A. Jiménez, J. Dudhia, J. F. GonzálezRouco, J. Navarro, J. P. Montávez, and E. GarcíaBustamante, “A revised scheme for the WRF surface layer formulation,” Monthly Weather Review, vol. 140, no. 3, pp. 898–918, 2012. View at: Publisher Site  Google Scholar
 S.Y. Hong, Y. Noh, and J. Dudhia, “A new vertical diffusion package with an explicit treatment of entrainment processes,” Monthly Weather Review, vol. 134, no. 9, pp. 2318–2341, 2006. View at: Publisher Site  Google Scholar
 M. Tiedtke, “A comprehensive mass flux scheme for cumulus parameterization in largescale models,” Monthly Weather Review, vol. 117, no. 8, pp. 1779–1800, 1989. View at: Publisher Site  Google Scholar
 C. Zhang, Y. Wang, and K. Hamilton, “Improved representation of boundary layer clouds over the southeast pacific in ARWWRF using a modified Tiedtke cumulus parameterization scheme,” Monthly Weather Review, vol. 139, no. 11, pp. 3489–3513, 2011. View at: Publisher Site  Google Scholar
 J. Molinari and M. Dudek, “Parameterization of convective precipitation in mesoscale numerical models: a critical review,” Monthly Weather Review, vol. 120, no. 2, pp. 326–344, 1992. View at: Publisher Site  Google Scholar
 J. W. Wilson, N. A. Crook, C. K. Mueller, J. Sun, and M. Dixon, “Nowcasting thunderstorms: a status report,” Bulletin of the American Meteorological Society, vol. 79, no. 10, pp. 2079–2099, 1998. View at: Publisher Site  Google Scholar
 Y.L. Lin, R. D. Farley, and H. D. Orville, “Bulk parameterization of the snow field in a cloud model,” Journal of Climate and Applied Meteorology, vol. 22, no. 6, pp. 1065–1092, 1983. View at: Publisher Site  Google Scholar
 W.K. Tao, J. Simpson, and M. McCumber, “An icewater saturation adjustment,” Monthly Weather Review, vol. 117, no. 1, pp. 231–235, 1989. View at: Publisher Site  Google Scholar
 G. Thompson, P. R. Field, R. M. Rasmussen, and W. D. Hall, “Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: implementation of a new snow parameterization,” Monthly Weather Review, vol. 136, no. 12, pp. 5095–5115, 2008. View at: Publisher Site  Google Scholar
 H. Hall, G. Thompson, and V. Tatarskii, “Impact of cloud microphysics on the development of trailing stratiform precipitation in a simulated squall line: comparison of one and twomoment schemes,” Monthly Weather Review, vol. 137, no. 3, pp. 991–1007, 2009. View at: Publisher Site  Google Scholar
 S.Y. Hong and J.O. J. Lim, “The WRF singlemoment 6class microphysics scheme (WSM6),” Journal of the Korean Meteorological Society, vol. 42, no. 2, pp. 129–151, 2006. View at: Google Scholar
 E. R. Mansell, C. L. Ziegler, and E. C. Bruning, “Simulated electrification of a small thunderstorm with twomoment bulk microphysics,” Journal of the Atmospheric Sciences, vol. 67, no. 1, pp. 171–194, 2010. View at: Publisher Site  Google Scholar
 S. A. Rutledge and P. Hobbs, “The mesoscale and Microscale structure and organization of clouds and precipitation in midlatitude cyclones. VIII: a model for the “SeederFeeder” process in warmfrontal rainbands,” Journal of the Atmospheric Sciences, vol. 40, no. 5, pp. 1185–1206, 1983. View at: Publisher Site  Google Scholar
 R. A. Houze Jr., P. V. Hobbs, P. H. Herzegh, and D. B. Parsons, “Size distributions of precipitation particles in frontal clouds,” Journal of the Atmospheric Sciences, vol. 36, no. 1, pp. 156–162, 1979. View at: Publisher Site  Google Scholar
 G. Parsons, R. M. Rasmussen, and K. Manning, “Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part I: description and sensitivity analysis,” Monthly Weather Review, vol. 132, no. 2, pp. 519–542, 2004. View at: Publisher Site  Google Scholar
 P. Field, R. Hogan, P. Brown et al., “Parametrization of ice‐particle size distributions for mid‐latitude stratiform cloud,” Quarterly Journal of the Royal Meteorological Society, vol. 131, no. 609, pp. 1997–2017, 2005. View at: Publisher Site  Google Scholar
 C. Marzban, S. Sandgathe, H. Lyons, and N. Lederer, “Three spatial verification techniques: cluster analysis, variogram, and optical flow,” Weather and Forecasting, vol. 24, no. 6, pp. 1457–1471, 2009. View at: Publisher Site  Google Scholar
 C. Lederer and G. C. Craig, “A displacementbased error measure applied in a regional ensemble forecasting system,” Monthly Weather Review, vol. 135, no. 9, pp. 3248–3259, 2007. View at: Publisher Site  Google Scholar
 C. Marzban and S. Sandgathe, “Optical flow for verification,” Weather and Forecasting, vol. 25, no. 5, pp. 1479–1494, 2010. View at: Publisher Site  Google Scholar
 S. Baker and I. Matthews, “Lucaskanade 20 years on: a unifying framework,” International Journal of Computer Vision, vol. 56, no. 3, pp. 221–255, 2004. View at: Publisher Site  Google Scholar
 B. D. Lucas and T. Kanade, “An iterative image registration technique with an application to stereo vision,” in Proceedings of the 7th International Joint Conference on Artificial Intelligence, pp. 121–130, San Francisco, CA, USA, 1981. View at: Google Scholar
 R. Rogers, S. Chen, J. Tenerelli et al., “A numerical study of the impact of vertical shear on the distribution of rainfall in Hurricane Bonnie (1998),” Monthly Weather Review, vol. 131, no. 8, 2003. View at: Publisher Site  Google Scholar
 A. P. Khain, K. D. Beheng, A. Heymsfield et al., “Representation of microphysical processes in cloudresolving models: spectral (bin) microphysics versus bulk parameterization,” Reviews of Geophysics, vol. 53, no. 2, pp. 247–322, 2015. View at: Publisher Site  Google Scholar
 M. S. Deshpande, S. Pattnaik, and P. S. Salvekar, “Impact of cloud parameterization on the numerical simulation of a super cyclone,” Annales Geophysicae, vol. 30, no. 5, pp. 775–795, 2012. View at: Publisher Site  Google Scholar
 S. A. Braun, “Highresolution simulation of Hurricane Bonnie (1998). Part II: water budget,” Journal of the Atmospheric Sciences, vol. 63, no. 1, pp. 43–64, 2006. View at: Publisher Site  Google Scholar
 D. J. Cecil and E. J. Zipser, “Reflectivity, ice scattering, and lightning characteristics of hurricane eyewalls and rainbands. Part II: intercomparison of observations,” Monthly Weather Review, vol. 130, no. 4, pp. 785–801, 2002. View at: Publisher Site  Google Scholar
 B. Reinhart, H. Fuelberg, R. Blakeslee et al., “Understanding the relationships between lightning, cloud microphysics, and airborne radarderived storm structure during hurricane Karl (2010),” Monthly Weather Review, vol. 142, no. 2, pp. 590–605, 2014. View at: Publisher Site  Google Scholar
 M. P. Meyers, P. J. DeMott, and W. R. Cotton, “New primary icenucleation parameterizations in an explicit cloud model,” Journal of Applied Meteorology, vol. 31, no. 7, pp. 708–721, 1992. View at: Publisher Site  Google Scholar
 C. Liu, K. Ikeda, G. Thompson, R. Rasmussen, and J. Dudhia, “Highresolution simulations of wintertime precipitation in the Colorado Headwaters region: sensitivity to physics parameterizations,” Monthly Weather Review, vol. 139, no. 11, pp. 3533–3553, 2011. View at: Publisher Site  Google Scholar
 J. D. McMillen and W. J. Steenburgh, “Impact of microphysics parameterizations on simulations of the 27 october 2010 great salt lakeeffect snowstorm,” Weather and Forecasting, vol. 30, no. 1, pp. 136–152, 2015. View at: Publisher Site  Google Scholar
 W. A. Cooper, “Ice initiation in natural clouds,” Precipitation EnhancementA Scientific Challenge, vol. 43, no. 12, pp. 29–32, 1986. View at: Publisher Site  Google Scholar
 E. Bigg, “The formation of atmospheric ice crystals by the freezing of droplets,” Quarterly Journal of the Royal Meteorological Society, vol. 79, no. 5, pp. 510–519, 1953. View at: Publisher Site  Google Scholar
 S.Y. Hong, J. Dudhia, and S.H. Chen, “A revised approach to ice microphysical processes for the bulk parameterization of clouds and precipitation,” Monthly Weather Review, vol. 132, no. 1, pp. 103–120, 2004. View at: Publisher Site  Google Scholar
 M. Rajeevan, A. Kesarkar, S. B. Thampi, T. N. Rao, B. Radhakrishna, and M. Rajasekhar, “Sensitivity of WRF cloud microphysics to simulations of a severe thunderstorm event over Southeast India,” Annales Geophysicae, vol. 28, no. 2, pp. 603–619, 2010. View at: Publisher Site  Google Scholar
 R. F. Adler, H.Y. M. Yeh, N. Prasad, W.K. Tao, and J. Simpson, “Microwave simulations of a tropical rainfall system with a threedimensional cloud model,” Journal of Applied Meteorology, vol. 30, no. 7, pp. 924–953, 1991. View at: Publisher Site  Google Scholar
 E. A. Smith, H. J. Cooper, X. Xiang, A. Mugnai, and G. J. Tripoli, “Foundations for statisticalphysical precipitation retrieval from passive microwave satellite measurements. Part I: brightnesstemperature properties of a timedependent cloudradiation model,” Journal of Applied Meteorology, vol. 31, no. 6, pp. 506–531, 1992. View at: Publisher Site  Google Scholar
 G. Hong, G. Heygster, J. Miao et al., “Sensitivity of microwave brightness temperatures to hydrometeors in a tropical deep convective cloud system at 89–190 GHz,” Radio Science, vol. 40, no. 4, pp. 1–13, 2005. View at: Publisher Site  Google Scholar
 M. Yanai, “Formation of tropical cyclones,” Reviews of Geophysics, vol. 2, no. 2, pp. 367–414, 1964. View at: Publisher Site  Google Scholar
 Y. Wang and C.C. Wu, “Current understanding of tropical cyclone structure and intensity changes ? a review,” Meteorology and Atmospheric Physics, vol. 87, no. 4, pp. 257–278, 2004. View at: Publisher Site  Google Scholar
 R. A. Anthes, “The dynamics and energetics of mature tropical cyclones,” Reviews of Geophysics, vol. 12, no. 3, pp. 495–522, 1974. View at: Publisher Site  Google Scholar
 B. K. Mahala, P. K. Mohanty, and B. K. Nayak, “Impact of microphysics schemes in the simulation of cyclone phailinusing WRF model,” Procedia Engineering, vol. 116, no. 1, pp. 655–662, 2015. View at: Publisher Site  Google Scholar
 T. N. Krishnamurti, K. S. Yap, and D. K. Oosterhof, “Sensitivity of tropical storm forecast to radiative destabilization,” Monthly Weather Review, vol. 119, no. 9, pp. 2176–2205, 1991. View at: Publisher Site  Google Scholar
 G. C. Craig, “Numerical experiments on radiation and tropical cyclones,” Quarterly Journal of the Royal Meteorological Society, vol. 122, no. 530, pp. 415–422, 2010. View at: Publisher Site  Google Scholar
 J. Guo, H. Liu, F. Wang et al., “Threedimensional structure of aerosol in China: a perspective from multisatellite observations,” Atmospheric Research, vol. 178179, no. 9, pp. 580–589, 2016. View at: Publisher Site  Google Scholar
 T. Chen, J. Guo, Z. Li et al., “ACloudSatPerspective on the cloud climatology and its association with aerosol perturbations in the vertical over eastern China,” Journal of the Atmospheric Sciences, vol. 73, no. 9, pp. 3599–3616, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Ying Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.