Articles | Volume 21, issue 9
https://doi.org/10.5194/acp-21-6919-2021
https://doi.org/10.5194/acp-21-6919-2021
Research article
 | 
06 May 2021
Research article |  | 06 May 2021

Microphysical processes producing high ice water contents (HIWCs) in tropical convective clouds during the HAIC-HIWC field campaign: evaluation of simulations using bulk microphysical schemes

Yongjie Huang, Wei Wu, Greg M. McFarquhar, Xuguang Wang, Hugh Morrison, Alexander Ryzhkov, Yachao Hu, Mengistu Wolde, Cuong Nguyen, Alfons Schwarzenboeck, Jason Milbrandt, Alexei V. Korolev, and Ivan Heckman
Abstract

Regions with high ice water content (HIWC), composed of mainly small ice crystals, frequently occur over convective clouds in the tropics. Such regions can have median mass diameters (MMDs) <300µm and equivalent radar reflectivities <20 dBZ. To explore formation mechanisms for these HIWCs, high-resolution simulations of tropical convective clouds observed on 26 May 2015 during the High Altitude Ice Crystals – High Ice Water Content (HAIC-HIWC) international field campaign based out of Cayenne, French Guiana, are conducted using the Weather Research and Forecasting (WRF) model with four different bulk microphysics schemes: the WRF single‐moment 6‐class microphysics scheme (WSM6), the Morrison scheme, and the Predicted Particle Properties (P3) scheme with one- and two-ice options. The simulations are evaluated against data from airborne radar and multiple cloud microphysics probes installed on the French Falcon 20 and Canadian National Research Council (NRC) Convair 580 sampling clouds at different heights. WRF simulations with different microphysics schemes generally reproduce the vertical profiles of temperature, dew-point temperature, and winds during this event compared with radiosonde data, and the coverage and evolution of this tropical convective system compared to satellite retrievals. All of the simulations overestimate the intensity and spatial extent of radar reflectivity by over 30 % above the melting layer compared to the airborne X-band radar reflectivity data. They also miss the peak of the observed ice number distribution function for 0.1<Dmax<1 mm. Even though the P3 scheme has a very different approach representing ice, it does not produce greatly different total condensed water content or better comparison to other observations in this tropical convective system. Mixed-phase microphysical processes at −10C are associated with the overprediction of liquid water content in the simulations with the Morrison and P3 schemes. The ice water content at −10C increases mainly due to the collection of liquid water by ice particles, which does not increase ice particle number but increases the mass/size of ice particles and contributes to greater simulated radar reflectivity.

1 Introduction

High concentrations of small ice particles ingested into jet engines can cause power loss and damaging events (Lawson et al.1998; Mason et al.2006). They can also cause air data probe failures (Duviver2010). Regions with high ice water content (HIWC), mainly composed of small ice crystals with median mass diameters (MMDs) as low as 300 µm, frequently occur over oceanic convective systems (Mason and Grzych2011; Ackerman et al.2015; Leroy et al.2016b). Such HIWC regions, with relatively low equivalent radar reflectivities (Ze) (often less than 20 dBZ; Mason et al.2006; Fridlind et al.2015; Protat et al.2016; Wolde et al.2016; Leroy et al.2017), are hard to detect with pilot radars onboard commercial aircraft and are thus potentially hazardous.

In order to explore the processes responsible for the occurrence of HIWC regions and the associated ice crystal properties within tropical convection, the High Altitude Ice Crystals – High Ice Water Content (HAIC-HIWC) international field campaigns (Dezitter et al.2013; Strapp et al.2016a) and the HIWC RADAR campaign (Ratvasky et al.2019) were conducted. The first HAIC-HIWC field campaign took place near Darwin, Australia, during the monsoon season of 2014, with the second out of Cayenne, French Guiana, in May 2015. The HIWC RADAR campaign was conducted out of Florida in August 2015. Data collected during these campaigns have been being analyzed to understand HIWC conditions (Leroy et al.2015, 2016a, 2017; Protat et al.2016; Wolde et al.2016; Korolev et al.2020), to develop warning products that can identify HIWC regions (Yost et al.2018; Bedka et al.2019; Harrah et al.2019; Haggerty et al.2020), and to characterize the high-altitude HIWC environment to assess a new ice crystal aircraft certification envelope.

In general, total condensed water content (TWC) values in these campaigns reached as high as 4.1 g m−3 averaged over 0.93 km (0.5 nautical mile) distance scales, and even up to about 2 g m−3 over 185 km distance scales. Average MMDs in HIWC zones greater than or equal to 1 g m−3 increased with temperature, from ∼326µm at −50C to ∼708µm at −10C (Strapp et al.2021). Leroy et al. (2015, 2016a, 2017) showed that MMDs decrease with increasing TWC and decreasing temperature, indicating small ice crystals are responsible for HIWC regions at high altitudes for both the Darwin and Cayenne datasets. Wolde et al. (2016) found the relationship between the ice water content (IWC) and radar equivalent reflectivity factor followed a power-law fit with coefficients dependent on temperature. However, they found the pilot X-band weather radar on the Canadian National Research Council (NRC) Convair-580 aircraft did not have adequate sensitivity to detect HIWC regions when calibrated using the NRC X-band research radar. Nguyen et al. (2019) proposed a retrieval method for IWC using the specific differential phase (Kdp) and differential reflectivity ratio (Zdr) data from X-band dual-polarization airborne radar. This method was demonstrated to be superior to the power-law fits between IWC and reflectivity as accounting for Zdr reduced the dependency of IWC on the variation in the shapes and orientation of ice particles. Ladino et al. (2017) concluded that secondary ice production (SIP) plays a dominant role in the formation of the observed high concentration ice crystals with ice nucleating particles making only a minor contribution. Korolev et al. (2020) proposed that a new “freezing-drop-shattering” mechanism generated small SIP particles above the melting layer at temperatures between 0 and −15C in both oceanic tropical mesoscale convective systems (MCSs) and midlatitude frontal clouds. In this SIP mechanism, large liquid drops are transported through the melting layer to a supercooled environment by convective turbulent updrafts and then collide with aged ice, freeze, and shatter. Keinert et al. (2020) used a laboratory study to indicate that bubble bursts dominate the SIP in sea salt drizzle droplets while droplet shattering controls the SIP in pure water droplets.

Several numerical simulation studies on tropical MCSs sampled during the HAIC-HIWC projects have been conducted using different numerical models and different microphysics schemes. Franklin et al. (2016) showed that the Met Office Unified Model (UM) with a single-moment microphysics scheme overestimated the radar reflectivity above the freezing level due to the errors of simulated updraft dynamics and particle sizes, hypothesizing that a double-moment microphysics scheme would improve the model's representation of the observed variability of the ice particle size distribution (PSD). The Stanford et al. (2017) WRF simulations of four tropical deep convection events sampled during the HAIC-HIWC Darwin campaign showed that three microphysics schemes (one bin and two double-moment bulk schemes; Lynn et al.2005; Thompson et al.2008; Morrison et al.2009) produced larger MMDs for TWC >1 g m−3 at temperatures between −10 and −40C and a high bias in convective radar reflectivity compared to observations. They hypothesized these differences resulted from errors in parameterized hydrometeor PSDs, single-ice-particle properties (e.g., shape and density), and parameterized microphysical processes. The Qu et al. (2018) simulation of a tropical convective system on 16 May 2015 with the Environment and Climate Change Canada's Global Environmental Multiscale (GEM) model and the Milbrandt–Yau double-moment cloud microphysical scheme (MY2) (Milbrandt and Yau2005) also produced IWC and ice particle number concentration differing from observations during the HAIC-HIWC Cayenne project, which they hypothesized was due to the poor representation of SIP processes in the microphysics scheme.

Other tropical convective clouds have also been observed and simulated. For example, McFarquhar and Heymsfield (1996) found the numbers of smaller particles (D<100µm) close to the convection were 1 order of magnitude higher than the numbers found further away according to data obtained during the Central Equatorial Pacific Experiment (CEPEX). Meanwhile, McFarquhar and Heymsfield (1997) indicated the shapes of the ice PSDs in the tropics substantially differ from those in the midlatitudes, especially at temperatures <-40C. Lohmann et al. (1995) showed that the average IWC simulated by a coarser resolution (125km×125km) general circulation model (GCM) agreed well with the observed IWC during CEPEX, especially with respect to the relationship between IWC and temperature, whereas the model underestimated the variability of simulated IWC within each temperature bin. Chen et al. (1997) indicated the main sources of ice particles are frozen cloud droplets and interstitial aerosol particles, and the number concentration of ice particles is influenced strongly by the amount of condensation nuclei in convective inflows according to the simulations and a sensitivity experiment for cases during CEPEX by using a one-dimensional microphysical model. Ackerman et al. (2015) conducted three 3D cloud-resolving model (CRM) simulations of MCSs observed on 23 January 2006 during the Tropical Warm Pool International Cloud Experiment (TWP-ICE) with bulk and bin microphysics schemes but were not able to produce HIWC regions (IWC >2 g m−3 and Ze<30 dBZ). Lang et al. (2011) greatly reduced model bias of excessively large reflectivity values (e.g., 40 dBZ) in the middle and upper troposphere in the simulation of a continental convective case observed during the Tropical Rainfall Measuring Mission (TRMM) Large-Scale Biosphere–Atmosphere Experiment in Amazonia (LBA) through modifying a single-moment bulk microphysics scheme in the Goddard Cumulus Ensemble model; however, there was much less improvement for an oceanic MCS observed during the TRMM Kwajalein Experiment (KWAJEX).

As indicated by the review above, the numerical studies on the HIWC phenomenon to date have not been able to capture the HIWC phenomenon well. This has been attributed to biases in particle properties, parameterized PSDs, and microphysical processes. The lack of knowledge about processes generating HIWC regions suggests that further numerical simulations are needed to explore the microphysical pathways producing HIWCs. Qu et al. (2018) indicated MY2 greatly overestimated the graupel content and hypothesized that HIWC will be better estimated by the next generation of microphysics schemes (e.g., the Predicted Particle Properties (P3) microphysics scheme; Morrison and Milbrandt2015). Therefore, this study will be the first test of the P3 scheme in simulating the HIWC phenomenon and in comparing P3 to other bulk microphysics schemes to determine whether P3 is better able to predict HIWCs in a high-resolution numerical weather prediction (NWP) context. In this study, a tropical oceanic MCS on 26 May 2015, which was well sampled during the Cayenne field campaign, is simulated at a high resolution with 1 km horizontal grid spacing. The numerical simulation experiments and their evaluation are described in this paper. In an upcoming companion paper, attention will be focused on sensitivity experiments varying some parameters within the microphysics scheme to enhance understanding of processes leading to the formation of small crystals in HIWC regions.

The next section describes the tropical oceanic MCS sampled on 26 May 2015. Section 3 introduces the collected data and how they were processed. Section 4 shows the simulated fields and their evaluation against observations. A summary and conclusions are presented in Sect. 5.

2 Case description

The tropical MCS observed on 26 May 2015 initiated and developed over the tropical Atlantic Ocean north of Cayenne, French Guiana. The MCS was not associated with obvious synoptic-scale flow features, such as clearly identifiable highs or lows. From the soundings at Cayenne, a deep moist absolutely unstable layer (MAUL) existed over the area where the MCS occurred, and the MAUL was maintained as the MCS developed, consistent with the conceptual model of a MCS proposed by Bryan and Fritsch (2000). There were mainly easterly (westerly) winds below (above) 350 hPa. The first convection initiated and developed over the ocean near the coast in the early morning. The convection moved eastward over the course of the day due to the upper westerly winds. Subsequently, new convective cells continually initiated and developed in a similar location and gradually moved eastward to merge with old convective cells that were present over the ocean to form a large and long-lived MCS.

The convective system was sampled by two research aircraft, the French SAFIRE Falcon 20 and Canadian NRC Convair 580, during the HAIC-HIWC field campaign. Figure 1 shows the observed brightness temperature from GOES-13 geostationary satellite channel 4 (10.8 µm) at 10:45 UTC on 26 May 2015, tracks of the two flights (Fig. 1a), and the height of the aircraft above mean sea level and air temperature at the flight levels (Fig. 1b). Both aircraft sampled close to the convective core of the storm as shown by the tracks of the aircraft through the lowest cloud-top brightness temperature (Fig. 1a). The SAFIRE Falcon 20 sampled at three height levels, i.e.,  7, ∼10, and ∼11.5 km, corresponding to temperatures of about −10, −30, and −45C, respectively. The NRC Convair 580 sampled mainly at ∼7 km and a temperature of around −10C (Fig. 1b).

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f01

Figure 1(a) Observed brightness temperature (K, shaded) from GOES-13 geostationary satellite channel 4 (10.8 µm) at 10:45 UTC on 26 May 2015 and flight tracks of SAFIRE Falcon 20 (black thick line) and NRC Convair 580 (blue thick line). (b) Height (km, thick solid curves) of aircraft above mean sea level and air temperature (C, thin dashed curves) at the flight level for SAFIRE Falcon 20 (black) and NRC Convair 580 (blue). The white mark “” in (a) represents the location of the release of a radiosonde at Cayenne.

3 Data and method

3.1 Data

The SAFIRE Falcon 20 was equipped with cloud microphysics instrumentation, including a Cloud Droplet Probe (CDP2), Two Dimensional Stereo Imaging Probe (2D-S), Precipitation Imaging Probe (PIP), and Isokinetic Evaporator Probe (IKP-2, Strapp et al.2016b). The NRC Convair 580 was equipped with an X-band (9.41 GHz) cloud airborne research radar (NAX, Wolde et al.2016) including three antennae (nadir, zenith, and side-looking) and similar cloud microphysics instrumentation. The two optical array probes, 2D-S and PIP, recorded 2D images of ice crystals nominally in the size range of 10–1280 and 100–6400 µm, respectively. The diode resolutions of 2D-S and PIP are 10 and 100 µm, respectively. The size distribution data with uncertainty of 10 %–100 % (Baumgardner et al.2017) are processed following the general approach described in McFarquhar et al. (2017), with only center-in particles accepted, and corrections for out-of-focus particles (Korolev2007), shattered particles (Field and Heymsfield2003; Field et al.2006; Korolev and Field2015), and particles partially within the photodiode array applied (Heymsfield and Parrish1978). Due to a poorly defined depth of field for small particles (Baumgardner et al.2012) and the potential of shattered artifacts only n(D) for Dmax>50µm are considered here. Composite PSDs ranging from 0.05 to 12.845 mm merged from the 2D-S and the PIP were derived at a 5 s time resolution using the 2D-S data below 800 µm, the PIP data above 1200 µm and complementary linear weights for both probes over the range of 800–1200 µm (Fontaine et al.2017; Leroy et al.2017). The IKP-2 bulk TWC probe was designed specifically for these campaigns to measure the high-speed, high-TWC environment, up to at least 10 g m−3 at 200 m s−1, with a target accuracy of 20 % (Strapp et al.2016b; Leroy et al.2017).

Radar reflectivity data from the X-band airborne research radar installed on the NRC Convair 580 (Wolde et al.2016), TWC measured by the IKP-2, and PSDs measured by the 2D-S and PIP installed on the SAFIRE Falcon 20 and on the NRC Convair 580 are used to statistically evaluate the model simulations. It should be noted that observed radar reflectivity profiles were sampled along the flight track of NRC Convair 580, whose horizontal locations differ from those of the flight track of the SAFIRE Falcon 20 (Fig. 1a). Part of observed TWC/PSD samples at −10C and all observed TWC/PSD samples at −45 and −30C are from SAFIRE Falcon 20 (Fig. 1b), implying that observed TWC/PSD at −45 and −30C may be not fully consistent with observed radar reflectivity in the upper levels in this study. Thermodynamic and wind profiles observed by a radiosonde released at Cayenne and the Cayenne GOES-13 Satellite Cloud Products Data (https://doi.org/10.5065/D6NC5ZX6, Langley Research Center2016) are also used.

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f02

Figure 2The model domain configuration (color shaded fields represent terrain elevation, in m). The horizontal grid spacings of d01 and d02 are 3 and 1 km, respectively. The white mark “” represents the location of the release of a radiosonde at Cayenne.

3.2 Model setup

The WRF model version 4.1.3 (Skamarock et al.2019) is used to simulate the tropical oceanic MCS event on 26 May 2015. Two one‐way nested domains with 3 and 1 km horizontal grid spacing and 51 vertical levels are adopted (Fig. 2). The ERA5 reanalysis data available every 1 h with 0.25× 0.25 horizontal grid spacing (https://doi.org/10.5065/BH6N-5N20, last access: 13 October 2019) are used for initial and boundary conditions. The model is run from 00:00 to 18:00 UTC on 26 May 2017 for 18 h with a spin-up time of the first 6 h. Physical parameterization schemes include the revised Rapid Radiative Transfer Model (RRTMG) longwave and shortwave radiation scheme (Iacono et al.2008), the Yonsei University (YSU) planetary boundary layer (PBL) scheme (Hong et al.2006), the MM5 similarity surface layer scheme (Beljaars1995), and the unified Noah land-surface scheme (Tewari et al.2004). The cumulus parameterization scheme is not activated in this study.

Four bulk microphysics schemes, namely the WRF single-moment 6-class (WSM6) microphysics scheme (Hong and Lim2006), the Morrison double-moment scheme (Morrison et al.2009) and the Predicted Particle Properties (P3) microphysics scheme with one- and two-ice options (Morrison and Milbrandt2015; Milbrandt and Morrison2016) are used for separate simulations. The simulations using the WSM6 and Morrison microphysics schemes are referred to as the WSM6 and MORR runs hereafter, respectively. These microphysics schemes are used because the PSDs of ice species are parameterized differently in these schemes. Both the WSM6 and Morrison schemes predict the mixing ratios of five cloud hydrometeor species, including cloud water, rainwater, cloud ice, snow, and graupel1, while the number mixing ratios for all species except cloud water are also predicted in the Morrison scheme. The P3 scheme predicts bulk ice properties (e.g., mean particle density) rather than predicting separate species of ice with fixed properties (e.g., cloud ice, snow, and graupel). P3 uses one or more “free” ice categories to represent all ice-phase hydrometeors, which can eliminate the unphysical “conversion” processes between different traditional ice categories (Morrison and Milbrandt2015; Milbrandt and Morrison2016). In this study, the options of one- and two-ice categories in the P3 scheme are used, referred to as P3-1ICE and P3-2ICE hereafter, respectively. Technically P3-1ICE and P3-2ICE are two configurations of the same scheme, but they have notably different treatments of ice which is the basis on which all the microphysics schemes were chosen, so these are referred to as different schemes. Output data in the model domain d02 with 1 km horizontal grid spacing are analyzed in this study. It should be noted that cloud ice, snow, and graupel in WSM6 and MORR and both categories of ice in P3-2ICE are treated as a single category of ice particles to compare with the observed ice particles, because the observed ice particles are not separated into different categories.

In WSM6, the main microphysical processes associated with ice particles include ice nucleation, deposition, sublimation, homogeneous and heterogeneous freezing, collection by liquid and other ice categories, autoconversion to other ice categories, melting, and sedimentation (Hong and Lim2006). MORR also considers the rime-splintering process in addition to the same microphysical processes as WSM6 (Morrison et al.2009). Both P3-1ICE and P3-2ICE consider the main microphysical processes associated with ice particles including ice nucleation, deposition, sublimation, homogeneous and heterogeneous freezing, collection by liquid, self-collection, melting, and sedimentation, while collision and collection between ice categories and rime-splintering processes are considered only in P3-2ICE (Morrison and Milbrandt2015; Milbrandt and Morrison2016).

3.3 Estimation of X-band radar reflectivity

The computations of simulated radar reflectivity are performed using the Rayleigh approximation which is applicable at the X-band given the size of typical ice particles (Ryzhkov et al.2020). The relations for reflectivity from rain (Zr), graupel (Zg), snow (Zs), and ice (Zi) are derived in detail in Appendix A. The total equivalent radar reflectivity factor (Ze) in units of dBZ can thus be attained using

(1) Z e = 10 × log 10 Z r + Z g + Z s + Z i 1 mm 6 m - 3 .
https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f03

Figure 3Skew-T plot of observed and simulated Cayenne radiosonde data at 12:00 UTC on 26 May 2015. One full wind barb represents 10 knots (∼5.14 m s−1). The location where Cayenne radiosonde was released is shown in Fig. 2.

Download

4 Results

4.1 Evaluation of simulated sounding, brightness temperature, and radar reflectivity

Figure 3 shows a Skew-T plot of the observed and simulated thermodynamic and wind profiles over Cayenne at 12:00 UTC on 26 May 2015. The observed profiles of air temperature and dew-point temperature show a very moist environment especially between 800 and 350 hPa. Regardless of the choice of microphysical scheme, the moist environment is well simulated although the layer between 500 and 350 hPa is drier with a maximum dew point depression (TTd) of ∼5C in P3-2ICE and almost zero in the observations. This upper-level drier layer is associated with the initial condition from ERA5 reanalysis data (not shown). The MORR scheme is moister between 500 and 350 hPa, which is more consistent with observations as the maximum TTd is less than 2 C (Fig. 3). As for the wind profiles, all the simulations predict the observed easterly and westerly winds in the lower and upper troposphere respectively, though there are some biases around 300 hPa and up to 250 hPa, where simulated and observed winds are in the opposite directions and wind speeds differ by ∼12 m s−1.

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f04

Figure 4Simulated and observed brightness temperature (BT, K, shaded) from GOES-13 channel 4 (10.8 µm) at 10:45 UTC on 26 May 2015.

Figure 4 shows simulated and observed brightness temperatures (BT) at 10:45 UTC on 26 May 2015. The simulated BT is calculated using the Community Radiative Transfer Model (CRTM, https://www.jcsda.org/jcsda-project-community-radiative-transfer-model, last access: 1 November 2019) using the assumptions consistent with those in the different microphysics schemes, including the characteristics of the cloud species and PSDs. The storm coverage (BT <232 K, yellow to deep red areas) in MORR (∼59 % of domain) is larger than that in the observations (∼47 % of domain), while WSM6 produces a smaller storm coverage (∼33 % of domain) compared to observations. The storm coverages (BT <232 K) in P3-1ICE (∼51 % of domain) and P3-2ICE (∼46 % of domain) are closer to that of the observations (Fig. 4). The lower brightness temperature areas representing deep convection (BT <212 K, red areas in Fig. 4) are larger in MORR (∼25 % of domain), P3-1ICE (∼20 % of domain), and P3-2ICE (∼13 % of domain) and smaller in WSM6 (∼5 % of domain) than in the observations (∼9 % of domain).

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f05

Figure 5Frequency (%, shaded) of simulated and observed BT from GOES-13 channel 4 from 06:15 to 17:45 UTC on 26 May 2015. Contours represent cumulative frequencies 5 %, 25 %, 50 %, 75 %, and 95 %, respectively.

Download

To examine the storm evolution, the frequency distributions of simulated and observed BT from 06:15 to 17:45 UTC on 26 May 2015 are displayed in Fig. 5. By 10:15 UTC, the dominant frequency (>4 %) of simulated and observed BT is around 280 K, indicating that clear regions dominate the domain in the early stage of MCS. The subdominant frequency (2 %–4 %) around 220 K illustrates the deep convection. This subdominant frequency in all of the simulations is consistent with observations, even though the BT ranges are all slightly different (i.e., 216–232 K in WSM6, 202–226 K in MORR, 200–228 K in P3-1ICE, 206–228 K in P3-2ICE, and 208–230 K in observations). There are 25 % of brightness temperatures <226 K in WSM6, <212 K in MORR, <214 K in P3-1ICE, <216 K in P3-2ICE, and <220 K in the observations at 10:15 UTC, indicating MORR, P3-1ICE, and P3-2ICE overpredict strong convective areas, and WSM6 underpredicts strong convective areas compared to the observations. After 10:15 UTC, the frequency around 220 K becomes dominant in the observations and simulations indicating the MCS develops and deep convective areas enlarge, as the maximum frequency exceeds 10 % in the observations and in MORR. The frequency ranges (>4 %) cover 222–246 K in WSM6, 206–234 K in MORR, 208–238 K in P3-1ICE, 210–236 K in P3-2ICE, and 212–232 K in the observations. There are 50 % of BT <238 K in WSM6, <220 K in MORR, <225 K in P3-1ICE, <228 K in P3-2ICE, and <227 K in the observation at 14:15 UTC. Overall this indicates that the storm coverage in MORR, P3-1ICE, and P3-2ICE is more consistent with the observations than that in WSM6. The frequency of brightness temperatures 214–224 K is over 10 % after 12:30 UTC in MORR, and during 15:00–16:00 UTC in the observations. This means deep convective areas as defined by BT metric are larger in MORR than the observations at an early stage in the system. Overall, all the simulations generally reproduce the storm coverage and evolution of this tropical MCS with the average bias (difference between simulation and observation) in storm coverage (BT <232 K) of -34.3 % in WSM6, ∼30.0 % in MORR, ∼12.9 % in P3-1ICE, and ∼2.3 % in P3-2ICE, indicating there is a relatively larger bias in WSM6 in which the storm coverage is smaller. WSM6 underestimates the observed deep convection areas (BT <212 K) by ∼55.5 % 2 and the overestimates in MORR, P3-1ICE and P3-2ICE are ∼175.4 %, ∼178.2 %, and ∼76.5 %, respectively (Fig. 5).

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f06

Figure 6Diagram of sampling method for radar reflectivity profiles. R=100 km here. The two magenta points at the flight track and model grid point represent the horizontal locations of observed and simulated radar reflectivity profiles. See text for details.

Download

To obtain a statistical comparison between simulated and observed radar reflectivity, contoured frequency by altitude diagrams (CFADs) (Yuter and Houze1995) are used. Because the observed reflectivity is only available along the flight track, equivalent locations for sampling the reflectivities from the modeled fields must be determined. The sampling method here is based on the flight track and observed and simulated BTs. Because there exists a bias between the simulated and observed BTs, they are first spatially normalized respectively within the same region as the model domain d02. Using normalized BTs, here paired samples of observations and simulations are found whose locations are about the same distance from the convective cores. Only the vertical profiles of observed reflectivity in which the observed TWC at the flight locations are larger than 0.1 g m−3 are selected as observational samples for the statistical comparison. Across an area with the horizontal location of observed reflectivity profiles as the center and 100 km as the range (Fig. 6), the simulated vertical profiles of reflectivity at model grid points where simulated TWC at the flight level is larger than 0.1 g m−3 and the normalized simulated brightness temperature is closest to the normalized observed brightness temperature are selected as simulation samples for statistical comparison. It should be noted that radii from 20 to 200 km in 20 km intervals were tested, and the results were similar. A radius of 100 km was adopted because the standard deviation between the observed and simulated reflectively was the lowest when using this radius threshold. The sampling method used here is similar to the one used by Borderies et al. (2018).

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f07

Figure 7Frequency (%, shaded) of simulated (a: WSM6, b: MORR, c: P3-1ICE, d: P3-2ICE) and observed (e) X-band radar reflectivity. Contours represent cumulative frequencies 25 %, 50 %, 75 %, and 95 %, respectively. The gap at 6–7 km in (e) is due to observed radar reflectivity data not available at the first few range gates (within ∼500 m) from the aircraft.

Download

In general, the simulations overestimate the X-band radar reflectivity above the melting layer (∼4.7 km). Figure 7 shows CFADs and cumulative CFADs of simulated and observed X-band radar reflectivity above 5 km. The CFADs are shown only above 5 km, because the formulae for calculating the simulated radar reflectivity in Sect. 3.3 do not consider the effect of melting, and this study mainly focuses on HIWC regions. From Fig. 7e, 95 % of the cumulative observed reflectivities are <30 dBZ above 6 km. The most frequently (frequency >10 %) observed radar reflectivity is around 25 dBZ at heights of 5–7 km and ∼15 and ∼20 dBZ at a height of ∼8 km (Fig. 7e). The simulated radar reflectivity shows broader distributions with larger values than the distributions of the 95 % cumulative observed reflectivities <30 dBZ above 6 km), and maxima in radar reflectivity can reach 50 dBZ, especially for WSM6, P3-1ICE, and P3-2ICE (Fig. 7a, c, and d). There are 95 % of the cumulative simulated reflectivities <44 dBZ in WSM6, <41 dBZ in MORR, <45 dBZ in P3-1ICE, and <47 dBZ in P3-2ICE above 6 km. The simulated radar reflectivity in MORR has a narrower distribution with 70 % of reflectivities between 34 and 42 dBZ at 5 km (Fig. 7b), which better resembles the observation with 70 % of reflectivities between 24 and 36 dBZ at 5 km (Fig. 7e). The other simulations have broader distributions with 70 % of the reflectivities between 30 and 44 dBZ in WSM6 (Fig. 7a), between 17 and 46 dBZ in P3-1ICE (Fig. 7c), and between 25 and 48 dBZ in P3-2ICE (Fig. 7d) at 5 km. The radar reflectivity in all simulations extends above 14 km, whereas the observed radar reflectivity is mainly below 14 km. Examining the reflectivity and Doppler velocity from zenith-viewing Doppler airborne radar shows a peak-to-peak correlation between them. This suggests there may be stronger updrafts in the simulations associated with the higher extended simulated reflectivity. Overall, all the simulations overestimate the intensity and spatial extent of radar reflectivity. By examining each component of reflectivity, the overestimation of radar reflectivity above the melting layer in WSM6 and MORR mainly results from the overprediction of graupel (not shown), which is similar to the tropical MCS simulations of Lang et al. (2011) and Qu et al. (2018). The P3 scheme, which was expected to yield better estimates of HIWCs, does not reduce the biases in simulated radar reflectivities. It should be noted that the NRC Convair 580 operations avoided the cloud regions with high reflectivity due to safety regulations, and thus it did not approach high reflectivity regions (red zones on the pilot’s radar) within 30 nautical miles (∼55.56 km). However, the BT sampling method has been used to minimize these aircraft sampling biases.

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f08

Figure 8Observed (red) and simulated (blue, a: WSM6, b: MORR, c: P3-1ICE, d: P3-2ICE) ice particle size distributions at the level of −45C. The red and blue dashed lines indicate the 25th and 75th percentiles, and the red and blue solid lines represent the median.

Download

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f09

Figure 9Observed (red) and simulated (blue, a: WSM6, b: MORR, c: P3-1ICE, d: P3-2ICE) ice particle size distributions at the level of −30C. The red and blue dashed lines indicate the 25th and 75th percentiles, and the red and blue solid lines represent the median.

Download

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f10

Figure 10Observed (red) and simulated (blue, a: WSM6, b: MORR, c: P3-1ICE, d: P3-2ICE) ice particle size distributions at the level of −10C. The red and blue dashed lines indicate the 25th and 75th percentiles, and the red and blue solid lines represent the median.

Download

4.2 Cloud microphysical properties

Samples are selected to examine the observed and simulated cloud microphysical properties using the same method as used for sampling radar reflectivity profiles. Compared to the observed ice PSDs, the simulated PSDs have different shapes and variability depending on the microphysics scheme. Figures 810 show observed and simulated ice PSDs at the levels of −45, −30, and −10C, respectively. These are the three temperature levels at which the in situ observations were focused. The simulated ice PSDs are the composite PSDs of all the ice-phase hydrometeors predicted in each microphysics scheme.

Table 1Number concentration of median n(D) for 0.1 mm <Dmax<12.845 mm (N0.1−12.845 mm) and contributions of small (0.1 mm <Dmax< 0.3 mm, C0.1−0.3 mm), medium (0.3 mm <Dmax<1 mm, C0.3−1 mm), and large (1 mm <Dmax<12.845 mm, C1−12.845 mm) particles at −45, −30, and −10C. The bold text indicates the dominant contribution.

Download Print Version | Download XLSX

At the −45C level (Fig. 8), compared to the observations, WSM6 underestimates the median number distribution function n(D) by ∼50 % (3×105 m−3 mm−1) near the maximum dimension (Dmax) of 0.1 mm (Fig. 8a), while P3-1ICE and P3-2ICE overestimate the median n(D) by ∼282 % (17×105 m−3 mm−1) and ∼199 % (12×105 m−3 mm−1) respectively at this Dmax (Fig. 8c and d). MORR has similar median n(D) magnitude to the observations with both 6×105 m−3 mm−1 near the Dmax of 0.1 mm (Fig. 8b). WSM6, MORR, and P3-2ICE underpredict the observed number concentration for 0.1 mm <Dmax<12.845 mm (N0.1−12.845 mm) by ∼17 %, ∼50 %, and ∼16 % respectively, while P3-1ICE overpredicts it by 5 times (Table 1). Small particles (0.1 mm <Dmax<0.3 mm) consistently make dominant contributions to the total number concentration in observations and simulations, while P3-1ICE produces too many small particles (∼86.6 %) with the overprediction of total number concentration (Table 1). Compared to the observed PSDs, PSDs in WSM6 and MORR have a smaller spread (Fig. 8a and b), and PSDs in P3-1ICE and P3-2ICE have a larger spread (Fig. 8c and d). Based on the shapes, magnitudes, and spreads of PSDs at −45C, PSDs in MORR are most consistent with the observations among the simulations.

At the −30C level (Fig. 9), all the shapes of simulated PSDs are similar to those at −45C, while the magnitudes are smaller by about an order of magnitude (Fig. 8). WSM6 and MORR have similar PSD characteristics in terms of spread and magnitude. The median n(D) in WSM6 and MORR have similar magnitude (∼105 m−3 mm−1) at Dmax of 0.1 mm compared to the observations (Fig. 9a and b). P3-1ICE and P3-2ICE overestimate the median n(D) by ∼6 times (5.7×105 m−3 mm−1) and ∼264 % (2.5×105 m−3 mm−1) respectively at Dmax of 0.1 mm (Fig. 9c and d). None of the simulations capture the peak of the observed PSD with the median of 2.6×105 m−3 mm−1 near Dmax of 0.3 mm (Fig. 9). There are no obvious peaks in PSDs for 0.1 mm <Dmax<1 mm in WSM6, MORR, and P3-2ICE (Fig. 9a, b, and d). There is a PSD peak with a median n(D) of 9.6×105 m−3 mm−1 near Dmax of 0.17 mm in P3-1ICE (Fig. 9c). Medium particles (0.3 mm <Dmax<1 mm) are dominant in the observations and WSM6, while small particles make the main contributions in MORR, P3-1ICE, and P3-2ICE (Table 1).

At the −10C level (Fig. 10), all simulations underestimate the median PSD for Dmax<1 mm, especially P3-2ICE, with an underestimate of ∼94 % (4.0×104 m−3 mm−1) at Dmax of 0.1 mm. Large particles contribute 58.6 % of N0.1−12.845 mm in P3-2ICE, which is very different from the observations (Table 1). Compared to the observations, MORR has almost the same PSD spread (range between the 25th and 75th percentiles) and median for Dmax of 0.1 mm (Fig. 10b). All of the simulations miss the peak of the observed PSD with a median of 1.5×105 m−3 mm−1 near Dmax of 0.3 mm (Fig. 10). There are enough samples with TWC >1 g m−3 at −10C to examine PSDs for 0.1 g m-3< TWC <1 g m−3 and TWC >1 g m−3 separately. Both results indicate that all of the simulations miss the peak of the observed PSD near Dmax of 0.3 mm (Figs. S1 and S2), which is consistent with that for TWC >0.1 g m−3 (Fig. 10).

Overall, the simulated PSDs at the three temperature levels shown in Figs. 810 have biases in various degrees with respect to their shapes, magnitudes, and spreads compared to observations. The observations are concentrated around the smaller crystal sizes than are the simulations for most of temperature levels, except P3-1ICE at −45 and −30C. These qualitative results do not change through examining the PSDs normalized by TWC at −45, −30, and −10C (not shown). It should be noted that the PSD comparison strongly reflects different assumptions built into the schemes regarding PSD shapes, namely inverse exponential PSDs are assumed in WSM6 and MORR, whereas a gamma PSD with a shape parameter (μ) that varies with the slope (λ) following the observations of Heymsfield (2003) that are used to formulate the P3 scheme.

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f11

Figure 11Violin plots of observed (red) and simulated (blue: WSM6, green: MORR, orange: P3-1ICE, magenta: P3-2ICE) (a) total water content (TWC, g m−3), (b) number concentration (# m−3) within 0.1 mm <Dmax<3 mm (N0.1−3 mm), (c) effective diameter (De, mm), and (d) vertical velocity (m s−1) at temperatures of −10, −30, and −45C. As for the black box-and-whisker plots, the extremes of the whiskers indicate the 5th and 95th percentiles, the lower and upper limits of the boxes correspond to the 25th and 75th percentiles, the dividing line represents the median value, and the white points represent the average value of samples between the 5th and 95th percentiles. The width of the shaded area represents the proportion of the data located there, and only areas between the 5th and 95th percentiles are shown. The numbers below the upper horizontal axis represent the total number of samples. The red triangles above the bottom horizontal axis indicate that the difference between the experiment mean and the observation mean is not statistically significant (not passing the significant test for p<0.05).

Download

Figure 11 shows statistical distributions of observed and simulated TWC, number concentration for particle sizes of 0.1 mm <Dmax<3 mm (N0.1−3 mm), effective diameter (De), and vertical velocity using violin plots (or box-percentile plots, Esty and Banfield2003) at temperatures of −10, −30, and −45C, respectively. The shaded areas of violin plots represent the proportion of the samples outlining the kernel probability densities. It should be noted that the number concentrations for 0.2 mm <Dmax<3 mm were also examined, and the conclusions are consistent. Thus, only the results using N0.1−3 mm are shown here. There are several different definitions of De (McFarquhar and Heymsfield1998). Given the main purpose here is to compare the particle sizes between observation and simulation for simplicity, the definition of De given by

(2) D e = 0.1 mm 12.845 mm k = 1 K D 3 n k ( D ) d D 0.1 mm 12.845 mm k = 1 K D 2 n k ( D ) d D

is used, where K is the number of ice species. Since D is maximum dimension, only one number distribution function that includes all data is used for the observations, and therefore K=1 for observations.

Generally all the simulations, especially MORR, reproduce the TWC reasonably well within the same order of magnitude as observations at the three temperature levels, with biases within 38 % at −10C. In particular at the −10C level the 25th and 75th percentiles of TWC in all the simulations cover the same order of magnitude as the observations (Fig. 11a). The differences in N0.1−3 mm among the simulations are quite large (Fig. 11b). At the three temperature levels, WSM6 and especially MORR underestimate the number concentration (Fig. 11b). This is associated with the underpredicted small particles and overpredicted large particles in WSM6 and MORR compared to the observations (Figs. 810). Thus, WSM6 and MORR produce larger De compared to the observations at the three temperature levels (Fig. 11c). At −30C, the median N0.1−3 mm values in WSM6 (0.4×105 m−3) and MORR (0.2×105 m−3) are underestimated by ∼50 % and ∼75 %, respectively, consistent with the underestimate of particle number near the peak of the observed PSD at Dmax of ∼0.3 mm (Fig. 9a and b). The N0.1−3 mm at −45C in P3-1ICE is about 1 order of magnitude larger than observed (Fig. 11b) mainly due to many more small particles for 0.1 mm <Dmax<0.4 mm (Fig. 8c). Similarly, compared to the observations, P3-1ICE overestimates N0.1−3 mm at −30C, with the median overestimated by ∼129 % (∼105 m−3) (Fig. 11b). This is explained by an overestimate of particle number at Dmax of ∼0.1 mm (Fig. 9c). Accordingly, P3-1ICE produces smaller De than the observations at −45 and −30C with an underestimate of median De by ∼38 % and ∼46 % respectively (Fig. 11c). Due to the larger PSD spreads in P3-2ICE at −45 and −30C (Figs. 8d and 9d), N0.1−3 mm and De in P3-2ICE accordingly have a larger spread than the observations (Fig. 11b and c). This occurs even though the spread of TWC between P3-2ICE and observations is similar (Fig. 11a). At −10C, values of N0.1−3 mm from all of the simulations, in particular P3-2ICE, are about 1 order of magnitude smaller than observed (Fig. 11b), implying larger mean particle size than observed (Fig. 11c). This is mainly attributed to the underestimate of small particle number for Dmax<1 mm, especially near the peak of the observed PSD at Dmax of ∼0.3 mm, and the overestimate of large particles in all of the simulations (Fig. 10). The simulated vertical velocity is in general stronger than in the observations, especially at −45 and −10C (Fig. 11c), corresponding to the higher extent of simulated radar reflectivity (Fig. 7).

4.3 Cloud microphysical processes

As discussed above, all the four microphysics schemes underpredict the number concentration by about 1 order of magnitude at −10C compared to the observations (Fig. 11b), although they predict similar TWC to the observed TWC (Fig. 11a). In this section, the processes producing HIWCs are determined. WSM6 is a single-moment scheme in which the number concentration of ice particles is not predicted directly. Through examining the number concentration derived diagnostically from the water content of each hydrometeor species in WSM6, it is found that their distributions are similar to those in MORR, especially the distributions of ice particles (not shown). Thus, only the double-moment MORR and P3 schemes are examined in detail here. Regions with IWC >1 g m−3 are defined as HIWC regions in this study. There are very limited observed and simulated HIWC samples at temperatures of −45 and −30C (Fig. 11a). Through comparing between simulated sampling profiles with TWC >0.1 g m−3 at −45 and −30C and simulated sampling profiles with TWC >1 g m−3 at −10C, it is found that their main microphysical processes are the same at the same vertical levels. Hence, the profiles of water content, number concentration, and microphysical processes with HIWC regions at −10C are used as examples to discuss here. The subsamples whose observed and corresponding simulated TWCs at −10C are larger than 1 g m−3 are selected from the total samples (1778) at −10C (Fig. 10) to conduct a composite analysis. There are 509, 488, and 427 paired samples selected for MORR, P3-1ICE, and P3-2ICE, respectively (Table 2).

Table 2Averaged TWC, N0.1−3 mm and air vertical velocity (w) of HIWC samples (TWC >1 g m−3) at −10C. The differences in all averages between the simulations and observations pass the significant test for p<0.05.

a Observation. b Model.

Download Print Version | Download XLSX

From Table 2, for the observation–MORR paired HIWC samples the average observed and simulated TWCs at −10C are similar, ∼1.6 and ∼1.8 g m−3, respectively, while the simulated N0.1−3 mm (1.78×104 m−3) is about 1 order of magnitude less than the observed N0.1−3 mm (1.72×105 m−3). The average simulated air vertical velocity for the HIWC points (∼0.32 m s−1) is about twice as large as the observations (∼0.16 m s−1). For the observation–P3-1ICE paired samples, the average observed and simulated TWCs at −10C are ∼1.7 and ∼1.8 g m−3 respectively, while the simulated N0.1−3mm (1.25×104 m−3) is underpredicted compared to the observed N0.1−3 mm (1.98×105 m−3) by ∼94 %. The simulated air vertical velocity (∼1.12 m s−1) is 3.67 times greater than observed (∼0.24 m s−1). Similarly, for the observation–P3-2ICE paired samples, the averaged observed and simulated TWCs at −10C are ∼1.7 and ∼1.9 g m−3 respectively, while the simulated N0.1−3 mm (4.63×103 m−3) is about 2 orders of magnitude less than the observed N0.1−3 mm (2.02×105 m−3), and the simulated air vertical velocity (∼2.05 m s−1) is about 15 times larger than the observation (∼0.13 m s−1). Therefore, both MORR and P3, in particular P3-2ICE, substantially underpredict the ice particle number for 0.1 mm <Dmax<3 mm and overpredict the vertical motion in the HIWC regions, which is associated with stronger and higher-extended simulated radar reflectivity (Fig. 7).

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f12

Figure 12Vertical profiles of averaged water content (a, c, e, solid curves, g m−3) and total number concentration (b, d, f, solid curves, # m−3) of each hydrometeor category in MORR (a, b – cloud water: blue; rain water: red; cloud ice: magenta; snow: orange; graupel: green; and sum of ice species: black), P3-1ICE (c, d – cloud water: blue; rain water: red; and cloud ice of first category: black), and P3-2ICE (e, f – cloud water: blue; rain water: red; cloud ice of first category: orange; cloud ice of second category: green; and sum of two ice categories: black). The thick black dashed lines represent vertical profiles of vertical velocity (m s−1). The vertical thin black dashed lines indicate vertical velocity of 0 m s−1. The horizontal thin black and gray dashed lines indicate the heights of temperatures at 0, −10, −30, and −45C from bottom to top, respectively. The numbers of samples used to calculate the averages are shown in Table 2.

Download

From the vertical profiles of water content for each hydrometeor class in MORR (Fig. 12a), graupel is dominant above the 0 C layer up to ∼8 km height, with snow being the second most important class. Snow is the dominant hydrometeor category above 8 km in terms of mass content. The aforementioned underestimate of ice particle number concentration at −10C is associated with the total number concentration (<4×104 m−3) of cloud ice, snow, and graupel (Fig. 12b). According to Eq. (A17), under the condition of similar simulated and observed TWCs, an underestimate of ice particle number concentration, especially graupel, leads to mass concentrating at the large end of PSD resulting in large reflectivities. Moreover, not only is the vertical motion overpredicted in the HIWC regions as mentioned above, but also the height with air vertical velocity >0 m s−1 in MORR can be up to 15 km with a maximum of ∼1 m s−1 at 13 km (Fig. 12a and b). These characteristics are consistent with the stronger intensity and higher distribution of simulated radar reflectivity shown in Fig. 7b.

The total ice number concentration in P3-1ICE is ∼104 m−3 at −10C, and the maximum is 4×105 m−3 at −45C (Fig. 12d). The air vertical velocity >0 m s−1 in P3-1ICE extends from heights of 1.5 to 15.5 km with a maximum of ∼1.4 m s−1 at 12.5 km (Fig. 12c and d), which is stronger than MORR. Compared to P3-1ICE, there is one more “free” ice category considered in P3-2ICE. The difference in mean mass-weighted diameters between the two ice categories is over 500 µm. The vertical motion in P3-2ICE is stronger with a maximum of ∼2.1 m s−1 at 6 km (Fig. 12e). Although the IWCs are similar between P3-1ICE and P3-2ICE, the total number concentration of the two ice categories in P3-2ICE is less than 5×103 m−3 at −10C, which is less than that in P3-1ICE (Fig. 12d and f). The average liquid water contents (LWCs) in MORR, P3-1ICE, and P3-2ICE at −10C are ∼0.029, ∼0.094, and ∼0.187 g m−3, respectively, while the observed LWC from the Cloud Droplet Probe (CDP) is less than 0.008 g m−3. Thus, all of the simulations, especially P3-2ICE, overpredict LWC at −10C and air vertical velocity above the 0 C layer. The distributions of vertical velocity (Fig. 11d) and Doppler velocity from zenith-viewing Doppler airborne radar (not shown) confirm that the model produces stronger updrafts than observations. It should be noted that both the NRC Convair 580 and SAFIRE Falcon 20 avoided cloud regions with strong updrafts where the presence of liquid phase is expected. However, the BT sampling method has been used to minimize these sampling biases.

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f13

Figure 13Vertical profiles of (a–c) averaged mass conversion rates (solid lines, 10−3 g m−3 s−1) and (d–f) averaged number concentration tendencies (solid lines, # m−3 s−1) of (a, d) cloud ice, (b, e) snow, and (c, f) graupel due to different microphysical processes in MORR. Only the profiles whose column-maximum conversion rates are larger than 10−6 g m−3 s−1 and column-maximum number concentration tendencies are larger than 1 m−3 s−1 are shown in (a–c) and (d–f), respectively. The total microphysics tendencies of cloud ice (a: Qi_TEND, d: Ni_TEND), snow (b: Qs_TEND, e: Ns_TEND), and graupel (c: Qg_TEND, f: Ng_TEND) are shown by black solid curves. The horizontal thin black and gray dashed lines indicate the heights of 0, −10, −30, and −45C from bottom to top, respectively. The numbers of samples used to calculate the averages are shown in Table 2. C: cloud ice; R: rain water; I: cloud ice; S: snow; and G: graupel.

Download

From vertical profiles of microphysical conversion rates in MORR (Fig. 13), the main source terms of cloud ice content are ice nucleation from homogeneous and heterogeneous freezing on aerosol at −45C and vapor deposition at −30 and −10C, and the main sink terms are collection by snow and autoconversion to snow. The net conversion rate of cloud ice (Qi_TEND, sum of all microphysical conversion rates including sedimentation) at −10C is negative (Fig. 13a). The net number concentration tendency of cloud ice (Ni_TEND) at −10C in MORR is -20 m−3 s−1, mainly due to the accretion of cloud ice, autoconversion to snow and sublimation (Fig. 13d). The IWC at −10C mainly consists of snow and graupel (Fig. 12a), and their main source terms are vapor deposition and collection of cloud water (Fig. 13b and c). The net snow conversion rate (Qs_TEND) is 0.1×10-3 g m−3 s−1 at −10C (Fig. 13b), while there is a negative net conversion rate of graupel (Qg_TEND) (-0.6×10-3 g m−3 s−1) mainly due to its large sedimentation tendency at −10C with a rate of -0.8×10-3 g m−3 s−1 (Fig. 13c). Autoconversion of cloud ice to snow (∼5 m−3 s−1) and the sedimentation term (∼1 m−3 s−1) increase the snow particle number; however, they are offset by the self-collection of snow (aggregation), collection of cloud water by snow to form graupel, and sublimation (Fig. 13e). The collection of cloud water by snow to form graupel is the main production term of graupel particle number, while it is offset by the sedimentation and sublimation terms. Finally, the net number concentration tendencies of both snow and graupel (Ns_TEND and Ng_TEND) are near 0 m−3 s−1 at −10C (Fig. 13e and f). Therefore, the collection of cloud water by graupel is the key source term of total IWC at −10C in MORR, which increases the mean mass/size of graupel and does not directly influence its number leading to mass concentrating at the large end of PSD. This is associated with the strong simulated reflectivity above the melting layer.

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f14

Figure 14Vertical profiles of (a, c, e) averaged mass conversion rates (solid lines, 10−3 g m−3 s−1) and (b, d, f) averaged number concentration tendency (solid lines, # m−3 s−1) of ice categories due to different microphysical processes in (a, b) P3-1ICE and (c–f) P3-2ICE. Only the profiles whose column-maximum conversion rates are larger than 10−6 g m−3 s−1 and column-maximum number concentration tendencies are larger than 1 m−3 s−1 are shown in (a, c, e) and (b, d, f), respectively. The total microphysics tendencies of first ice category (a, c Qi_TEND; b, d: Ni_TEND) and second ice category (e: Qi2_TEND; f: Ni2_TEND) are shown by black solid curves. The horizontal thin black and gray dashed lines indicate the heights of 0, −10, −30, and −45C from bottom to top, respectively. The numbers of samples used to calculate the averages are shown in Table 1. C: ice; R: rain water; I: cloud ice.

Download

From vertical profiles of microphysical conversion rates in P3-1ICE (Fig. 14a), the main production terms of ice content are vapor deposition at −45 and −30C, collection of cloud water by ice, vapor deposition, and collection of rain water by ice at −10C. High ice number concentrations at the upper levels in P3-1ICE (Fig. 11b) should be associated with homogeneous freezing of cloud droplets (Fig. 14b). As for the first ice category in P3-2ICE, in addition to the same main production terms as in P3-1ICE, there is another source term merging with the second ice category due to similar mean mass-weighted diameters between the two ice categories (Fig. 14c). The net tendency of the two ice categories in P3-2ICE (i.e., Qi_TEND + Qi2_TEND) is 1.05×10-3 g m−3 s−1 at −10C, which is much larger than that in P3-1ICE (0.02×10-3 g m−3). It is mainly due to the stronger collection of cloud water and rain water by ice in P3-2ICE, which may be associated with the greater cloud water and rain water content at −10C in P3-2ICE than P3-1ICE (Fig. 12c and e). The aforementioned collection of cloud water and rain water by ice does not increase the ice particle number. Although the deposition nucleation can increase the ice particle number in P3-1ICE and P3-2ICE, it is small (less than 0.5 m−3 s−1) at −10C. Merging ice categories does not increase the total ice particle number in P3-2ICE. The sedimentation of ice number is ∼3.3 and ∼3.4 m−3 s−1 at −10C in P3-1ICE and P3-2ICE respectively, which dominates the net ice number concentration tendencies at −10C in both P3-1ICE (∼2.8 m−3 s−1) and P3-2ICE (∼3.1 m−3 s−1). The much lower number concentration in P3-2ICE than P3-1ICE (Table 1) is likely due to aggregation associated with collection between the two ice categories in P3-2ICE (Fig. 14).

To summarize briefly, due to the overprediction of LWC in MORR, P3-1ICE, and P3-2ICE above the melting layer, there exist obvious mixed-phase processes at −10C. The IWC at −10C increases mainly due to the collection of liquid water by ice particles, which does not increase ice particle number but increases the mass and size of ice particles. The lower ice particle numbers in the simulations could also be associated with excessive aggregation and/or missing SIPs, such as collision-induced breakup and “freezing-drop-shattering” proposed by Korolev et al. (2020). The large ice particles and lower ice particle numbers lead to mass concentrating at the large end of PSD contributing to strong simulated radar reflectivity. Introduction of parameterizations for the missing SIPs may be able to overcome some of these model limitations, as will be examined in the future.

https://acp.copernicus.org/articles/21/6919/2021/acp-21-6919-2021-f15

Figure 15As Fig. 7 but for 3 km (d01) simulations.

Download

5 Summary and conclusions

A tropical oceanic convective system observed on 26 May 2015 during the HAIC-HIWC international field campaign based out of Cayenne, French Guiana, was simulated using the WRF model. Observation data from radiosondes, GOES-13 geostationary satellite, airborne radar, and cloud microphysics instrumentation were used to assess the simulated convective system in terms of the thermodynamic and dynamic environment, storm coverage, evolution and structure, and microphysical properties. The major results are summarized as follows:

  1. By comparing simulated and observed soundings, all of the simulations using different microphysics schemes replicate temperature with average bias within 1.6 %, dew-point temperature with average bias within 6 %, wind speed with average bias within 14 %, and wind direction with average bias within 36, with the MORR scheme giving closest agreement with observations.

  2. WRF basically reproduces the coverage and evolution of this tropical MCS based on a comparison between simulated and observed brightness temperature with the average bias in storm coverage (brightness temperature <232 K) by -34.3 % in WSM6, ∼30.0 % in MORR, ∼12.9 % in P3-1ICE, and ∼2.3 % in P3-2ICE. Thus, WSM6 underestimates the storm coverage, and P3-2ICE produces the closest storm coverage to the observation.

  3. In general, all of the simulations overestimate the intensity and spatial extent of radar reflectivity above the melting layer compared to the observed X-band radar reflectivity. There are 95 % of the simulated reflectivities <44 dBZ in WSM6, <41 dBZ in MORR, <45 dBZ in P3-1ICE, and <47 dBZ in P3-2ICE above 6 km, while 95 % of the observed reflectivities are <30 dBZ above 6 km. The radar reflectivity >0 dBZ in all simulations extends above 14 km, whereas the observed radar reflectivity >0 dBZ is mainly below 14 km.

  4. Different microphysics schemes have different shapes, magnitudes, and spreads in the simulation of ice PSDs at the different temperature levels. All of the simulations miss the peak of the observed ice number distribution function for 0.1<Dmax<1 mm.

  5. Both the WSM6 and MORR schemes underestimate the number concentration of ice particles at the temperature levels of −45, −30, and especially −10C with a maximum bias up to 1 order of magnitude, though the simulated total water contents are similar to the observations. This indicates that WSM6 and MORR simulate fewer small particles and more large particles compared to the observations. P3-1ICE overestimates the number concentration with a maximum bias up to 1 order of magnitude, and P3-2ICE generates a larger spread of number concentrations covering about 2 orders of magnitude at temperatures of −45 and −30C. Both schemes, and especially P3-2ICE, underestimate the number concentration by about 1 order of magnitude at −10C. This indicates the P3 scheme produces more large particles at this level.

  6. Mixed-phase processes play an important role at −10C due to the overprediction of LWC in MORR, P3-1ICE, and P3-2ICE above the melting layer. Stronger simulated radar reflectivity in MORR above the melting layer results from large graupel particles associated with greater graupel water content and fewer graupel particles compared with in situ observations. Rapid growth of graupel mass mainly through collecting cloud water but with a limited increase in graupel number mainly by conversion of snow to graupel through the collection of cloud water above the melting layer leads to large mean graupel sizes in MORR. Similarly, in P3-1ICE and P3-2ICE the IWC at −10C increases mainly due to the collection of cloud water and rain water while the net ice number concentration tendencies are near 0, which generates large mean ice particle sizes. The large ice particles generate strong radar reflectivity, partially explaining the bias of simulated radar reflectivity with P3-1ICE and P3-2ICE.

It should be noted that simulations of deep convection at different model resolutions can be much different. To examine the sensitivity of model resolution, CFADs and cumulative CFADs of radar reflectivity are also calculated using the simulation data from the 3 km domain (Fig. 15). Although there are some differences in specific values of reflectivity, the intensity and distribution of reflectivity from the 3 km simulations (Fig. 15) are basically consistent with those from the 1 km simulations (Fig. 7). Although this does not prove that the conclusions with respect to the differences in behavior among the microphysics schemes will be the same at all resolutions, it does at least indicate that the results in this study have some generality. It provides something useful about the microphysics schemes to numerical forecast guidance for HIWCs in current high-resolution NWP models, which are now routinely run at O(3 km) and now more and more often at O(1 km). However, a caveat is that 1 km is not cloud-resolving O(100 m), and thus horizontal entrainment is still not being resolved at this resolution, which affects the amount of liquid available for riming growth in updrafts. Therefore, there still exist uncertainties in the 1 km simulations.

In conclusion, no one microphysics scheme outperforms the other scheme in simulating this tropical oceanic MCS, as is evident from examining the simulated soundings, brightness temperature, radar reflectivity, ice particle size distributions, total water content, and number concentration. The Morrison scheme underestimates the number concentration at different temperature levels compared to the observations. This indicates that large ice particles, especially graupel, are overpredicted in this scheme, which is similar to the Qu et al. (2018) simulation of a different tropical MCS using a different model and microphysics scheme. Even though the P3 scheme has a much different approach for representing ice, it does not produce greatly different TWC or better comparison to the observations using either one- or two-ice categories. This suggests that other aspects need to be considered, such as microphysical process rate formulations or parameters. To enhance understanding of processes leading to the formation of small crystals in HIWC regions, sensitivity experiments varying parameters within the P3 microphysics scheme (e.g., mass-dimensional relations, size distribution parameters, microphysical conversion rates, or representation of different processes like secondary ice production) will be examined in a future paper.

Appendix A: Formulae for X-band radar reflectivity factor

Most cloud models with bulk microphysics parameterizations predict either mass content (one-moment schemes) or mass content and number concentration (two-moment schemes) for a number of hydrometeor categories. In one-moment schemes, number concentration can be obtained diagnostically. They also commonly assume inverse exponential or gamma size distributions of hydrometeors with respect to particle maximum dimension. This allows the use of simple analytical formulae for converting mass content and number concentration to radar reflectivity factor if the scatterers are small compared to the radar wavelength so that the Rayleigh approximation can be used. A gamma size distribution is represented by

(A1) N ( D ) = N 0 D μ e - Λ D ,

where N0, μ, and λ are the intercept, shape, and slope factors respectively. The total number concentration Nt is thus given by

(A2) N t = 0 N ( D ) d D = N 0 Γ ( μ + 1 ) Λ μ + 1 .

A1 Calculation of radar reflectivity factor for rain

For rain, the radar reflectivity factor Z in the Rayleigh approximation is given by the well known formula

(A3) Z = 0 D eq 6 N ( D eq ) d D eq ,

where Deq is the equivolume diameter of raindrop, and Deq=D for liquid water species in bulk schemes. The liquid water content LWC is defined as

(A4) LWC = π 6 ρ w 0 D 3 N ( D ) d D ,

where ρw is the density of water. Taking into account Eq. (A2), integration of Eqs. (A3) and (A4) yields

(A5) Z = N t Γ ( μ + 7 ) Λ 6 Γ ( μ + 1 ) ,

and

(A6) LWC = π ρ w N t Γ ( μ + 4 ) 6 Λ 3 Γ ( μ + 1 ) ,

which results in the following formula for estimating Z from LWC and Nt with an inverse exponential size distribution assumption (μ=0) in WSM6, MORR, and P3:

(A7) Z = 720 LWC 2 π 2 ρ w 2 N t .

A2 Calculation of radar reflectivity factor for ice particles

The radar reflectivity factor for ice is defined as

(A8) Z = | K i | 2 | K w | 2 ρ i 2 0 ρ s 2 ( D eq ) D eq 6 N ( D eq ) d D eq

(Ryzhkov and Zrnić2019, Eq. 5.14), where ρi is the density of solid ice sphere, assumed here to be 0.917 g cm−3, and

(A9) K w = ε w - 1 ε w + 2 ,  and  K i = ε i - 1 ε i + 2 ,

where εw and εi are the dielectric constants of water and solid ice, and |Kw|2=0.930 and |Ki|2=0.176 in this study. Taking into account that the mass of ice particle ms can be expressed as

(A10) m s ( D eq ) = π 6 ρ s ( D eq ) D eq 3 .

Equation (A8) can be written in a different form often used in cloud models (e.g., Hogan et al.2006):

(A11) Z = | K i | 2 | K w | 2 6 π ρ i 2 0 m s 2 ( D eq ) N ( D eq ) d D eq ,

where Deq is the equivalent volume diameter of an ice particle (Ryzhkov and Zrnić2019).

The msDmax relations are commonly used (e.g., Locatelli and Hobbs1974; Mitchell1996; Finlon et al.2019; Ding et al.2020), where Dmax is the maximal dimension of ice particles. These relations are often represented as power-law dependencies

(A12) m s ( D max ) = a D max b .

Then Eq. (A11) can be rewritten as

(A13) Z = | K i | 2 | K w | 2 6 π ρ i 2 0 a 2 D max 2 b N ( D eq ) d D eq = | K i | 2 | K w | 2 6 π ρ i 2 0 a 0 2 D eq 2 b N ( D e q ) d D e q ,

where a0=aηb. In Eq. (A13), the difference between maximal dimension Dmax and equivolume diameter Deq is taken into account with a scaling factor η=Dmax/Deq assuming a constant aspect ratio of ice particles across the size spectrum.

A2.1 Constant mD relation across the size spectrum

In WSM6 and MORR, a and b in the msDmax relations are assumed to be constant across the size spectrum for snow and graupel particles, implying densities of snow and graupel are constant. This leads to the following expressions for Z and ice water content IWC:

(A14) Z = | K i | 2 | K w | 2 6 π ρ i 2 Γ ( 2 b + μ + 1 ) N t Γ ( μ + 1 ) a 0 2 Λ 2 b

and

(A15) IWC = 0 a 0 D eq b N ( D eq ) d D eq = Γ ( b + μ + 1 ) N t Γ ( μ + 1 ) a 0 Λ b ,

so that

(A16) Z = | K i | 2 | K w | 2 6 π ρ i 2 Γ ( μ + 1 ) Γ ( 2 b + μ + 1 ) [ Γ ( b + μ + 1 ) ] 2 IWC 2 N t .

It is important that Eq. (A16) is not sensitive to the variability of the prefactor a in the ms(Dmax) power-law relation (A12) and is minimally affected by the variability of the exponent b in Eq. (A12). In WSM6 and MORR, b=3 and μ=0, and thus Eq. (A16) can be simplified as

(A17) Z = 720 π 2 ρ i 2 | K i | 2 | K w | 2 IWC 2 N t .

A2.2 Variable mD relation across the size spectrum

P3 represents each ice category much differently than WSM6 and MORR. In P3, a and b in the msDmax relations (Eq. A12) are variable across the size spectrum (Morrison and Milbrandt2015). Therefore, radar reflectivity factor is calculated numerically by using Eq. (A13).

Code and data availability

The WRF code is available at https://github.com/wrf-model/WRF (last access: 4 January 2020). Observation data are available at https://data.eol.ucar.edu/master_lists/generated/haic-hiwc_2015 (last access: 26 May 2020). ERA5 reanalysis data are available at https://doi.org/10.5065/BH6N-5N20 (European Centre for Medium-Range Weather Forecasts2019).

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/acp-21-6919-2021-supplement.

Author contributions

YoH, WW, and GMM designed the study. YoH did the calculations, with support from WW, GMM, XW, and HM. AR developed the method for calculating X-band radar reflectivity. MW, CN, AS, AVK, IH, and YaH processed the original observational datasets. YoH wrote the original draft with contributions from all coauthors, and all coauthors contributed to review and editing.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

This work was supported by the National Science Foundation (award numbers: 1213311 and 1842094). Observational data are provided by NCAR/EOL under the sponsorship of the National Science Foundation (https://data.eol.ucar.edu/, last access: 26 May 2020). The authors are grateful to NCAR's Data Support Section for providing ERA5 reanalysis data (https://doi.org/10.5065/BH6N-5N20, last access: 13 October 2019). The authors acknowledge high-performance computing support from Cheyenne (https://doi.org/10.5065/D6RX99HX) provided by NCAR's Computational and Information Systems Laboratory, sponsored by the National Science Foundation. NCAR is sponsored by the National Science Foundation. Some of the computing for this project was performed at the University of Oklahoma (OU) Supercomputing Center for Education and Research (OSCER). The discussions of radar forward simulators with Jiaxi Hu and Djordje Mirkovic, and the discussions of HIWC conditions and IKP with Walter Strapp are greatly appreciated. Major North American funding for flight campaigns was provided by the FAA William Hughes Technical Center and Aviation Weather Research Program, the NASA Aeronautics Research Mission Directorate Aviation Safety Program, the Boeing Co., Environment and Climate Change Canada, the NRC of Canada, and Transport Canada. Major European campaign and research funding was provided by (i) the European Commission Seventh Framework Program in research, technological development and demonstration under grant agreement no. ACP2-GA-2012-314314 and (ii) the European Aviation Safety Agency (EASA) Research Program under service contract no. EASA.2013.FC27. Further funding was provided by the Ice Crystal Consortium.

Financial support

This research has been supported by the National Science Foundation (grant nos. 1213311 and 1842094).

Review statement

This paper was edited by Martina Krämer and reviewed by three anonymous referees.

References

Ackerman, A. S., Fridlind, A. M., Grandin, A., Dezitter, F., Weber, M., Strapp, J. W., and Korolev, A. V.: High ice water content at low radar reflectivity near deep convection – Part 2: Evaluation of microphysical pathways in updraft parcel simulations, Atmos. Chem. Phys., 15, 11729–11751, https://doi.org/10.5194/acp-15-11729-2015, 2015. a, b

Baumgardner, D., Avallone, L., Bansemer, A., Borrmann, S., Brown, P., Bundke, U., Chuang, P. Y., Cziczo, D., Field, P., Gallagher, M., Gayet, J.-F., Heymsfield, A., Korolev, A., Krämer, M., McFarquhar, G., Mertes, S., Möhler, O., Lance, S., Lawson, P., Petters, M. D., Pratt, K., Roberts, G., Rogers, D., Stetzer, O., Stith, J., Strapp, W., Twohy, C., and Wendisch, M.: In situ, airborne instrumentation: Addressing and solving measurement problems in ice clouds, B. Am. Meteorol. Soc., 93, ES29–ES34, https://doi.org/10.1175/BAMS-D-11-00123.1, 2012. a

Baumgardner, D., Abel, S. J., Axisa, D., Cotton, R., Crosier, J., Field, P., Gurganus, C., Heymsfield, A., Korolev, A., Krämer, M., Lawson, P., McFarquhar, G., Ulanowski, Z., and Um, J.: Cloud Ice Properties: In Situ Measurement Challenges, Meteorol. Monogr., 58, 9.1–9.23, https://doi.org/10.1175/AMSMONOGRAPHS-D-16-0011.1, 2017. a

Bedka, K., Yost, C., Nguyen, L., Strapp, J. W., Ratvasky, T., Khlopenkov, K., Scarino, B., Bhatt, R., Spangenberg, D., and Palikonda, R.: Analysis and Automated Detection of Ice Crystal Icing Conditions Using Geostationary Satellite Datasets and In Situ Ice Water Content Measurements, SAE Int. J. Adv. Curr. Pract. Mobility, 2, 35–57, https://doi.org/10.4271/2019-01-1953, 2019. a

Beljaars, A. C.: The parametrization of surface fluxes in large-scale models under free convection, Q. J. Roy. Meteor. Soc., 121, 255–270, 1995. a

Borderies, M., Caumont, O., Augros, C., Bresson, É., Delanoë, J., Ducrocq, V., Fourrié, N., Bastard, T. L., and Nuret, M.: Simulation of W-band radar reflectivity for model validation and data assimilation, Q. J. Roy. Meteor. Soc., 144, 391–403, 2018. a

Bryan, G. H. and Fritsch, J. M.: Moist absolute instability: The sixth static stability state, B. Am. Meteorol. Soc., 81, 1207–1230, 2000. a

Chen, J.-P., McFarquhar, G. M., Heymsfield, A. J., and Ramanathan, V.: A modeling and observational study of the detailed microphysical structure of tropical cirrus anvils, J. Geophys. Res.-Atmos., 102, 6637–6653, 1997. a

Dezitter, F., Grandin, A., Brenguier, J.-L., Hervy, F., Schlager, H., Villedieu, P., and Zalamansky, G.: HAIC (High Altitude Ice Crystals), in: 5th AIAA Atmospheric and Space Environments Conf., San Diego, CA, AIAA, AIAA-2013-2674, 2674, https://doi.org/10.2514/6.2013-2674, 2013. a

Ding, S., McFarquhar, G. M., Nesbitt, S. W., Chase, R. J., Poellot, M. R., and Wang, H.: Dependence of Mass–Dimensional Relationships on Median Mass Diameter, Atmosphere, 11, 756, https://doi.org/10.3390/atmos11070756, 2020. a

Duviver, E.: High Altitude Icing Environment, in: Intl. Air Safety and Climate Change Conf., 8–9 September 2010, Cologne, DE, available at: https://www.easa.europa.eu/newsroom-and-events/events/international-air-safety-climate-change-conference-iascc-2010 (last access: 28 September 2020), 2010. a

Esty, W. W. and Banfield, J. D.: The box-percentile plot, Journal of Statistical Software, 8, 1–14, 2003. a

European Centre for Medium-Range Weather Forecasts: ERA5 Reanalysis (0.25 Degree Latitude-Longitude Grid). Research Data Archive at the National Center for Atmospheric Research [dataset], Computational and Information Systems Laboratory, https://doi.org/10.5065/BH6N-5N20, 2019, updated monthly. a

Field, P. R. and Heymsfield, A. J.: Aggregation and scaling of ice crystal size distributions, J. Atmos. Sci., 60, 544–560, 2003. a

Field, P., Heymsfield, A., and Bansemer, A.: Shattering and particle interarrival times measured by optical array probes in ice clouds, J. Atmos. Ocean. Tech., 23, 1357–1371, 2006. a

Finlon, J. A., McFarquhar, G. M., Nesbitt, S. W., Rauber, R. M., Morrison, H., Wu, W., and Zhang, P.: A novel approach for characterizing the variability in mass–dimension relationships: results from MC3E, Atmos. Chem. Phys., 19, 3621–3643, https://doi.org/10.5194/acp-19-3621-2019, 2019. a

Fontaine, E., Leroy, D., Schwarzenboeck, A., Delanoë, J., Protat, A., Dezitter, F., Grandin, A., Strapp, J. W., and Lilie, L. E.: Evaluation of radar reflectivity factor simulations of ice crystal populations from in situ observations for the retrieval of condensed water content in tropical mesoscale convective systems, Atmos. Meas. Tech., 10, 2239–2252, https://doi.org/10.5194/amt-10-2239-2017, 2017. a

Franklin, C. N., Protat, A., Leroy, D., and Fontaine, E.: Controls on phase composition and ice water content in a convection-permitting model simulation of a tropical mesoscale convective system, Atmos. Chem. Phys., 16, 8767–8789, https://doi.org/10.5194/acp-16-8767-2016, 2016. a

Fridlind, A. M., Ackerman, A. S., Grandin, A., Dezitter, F., Weber, M., Strapp, J. W., Korolev, A. V., and Williams, C. R.: High ice water content at low radar reflectivity near deep convection – Part 1: Consistency of in situ and remote-sensing observations with stratiform rain column simulations, Atmos. Chem. Phys., 15, 11713–11728, https://doi.org/10.5194/acp-15-11713-2015, 2015. a

Haggerty, J. A., Rugg, A., Potts, R., Protat, A., Strapp, J. W., Ratvasky, T., Bedka, K., and Grandin, A.: Development of a Method to Detect High Ice Water Content Environments Using Machine Learning, J. Atmos. Ocean. Tech., 37, 641–663, https://doi.org/10.1175/JTECH-D-19-0179.1, 2020. a

Harrah, S., Strickland, J., Hunt, P., Proctor, F., Switzer, G., Ratvasky, T., Strapp, J. W., Lilie, L., and Dumont, C.: Radar Detection of High Concentrations of Ice Particles-Methodology and Preliminary Flight Test Results, Tech. rep., SAE Technical Paper, https://doi.org/10.4271/2019-01-2028, 2019. a

Heymsfield, A. J.: Properties of tropical and midlatitude ice cloud particle ensembles. Part II: Applications for mesoscale and climate models, J. Atmos. Sci., 60, 2592–2611, 2003. a

Heymsfield, A. J. and Parrish, J. L.: A computational technique for increasing the effective sampling volume of the PMS two-dimensional particle size spectrometer, J. Appl. Meteorol., 17, 1566–1572, 1978. a

Hogan, R. J., Mittermaier, M. P., and Illingworth, A. J.: The retrieval of ice water content from radar reflectivity factor and temperature and its use in evaluating a mesoscale model, J. Appl. Meteorol. Climatol., 45, 301–317, 2006. a

Hong, S.-Y. and Lim, J.-O. J.: The WRF single-moment 6-class microphysics scheme (WSM6), Asia-Pacific J. Atmos. Sci., 42, 129–151, 2006. a, b

Hong, S.-Y., Noh, Y., and Dudhia, J.: A new vertical diffusion package with an explicit treatment of entrainment processes, Mon. Weather Rev., 134, 2318–2341, 2006. a

Iacono, M. J., Delamere, J. S., Mlawer, E. J., Shephard, M. W., Clough, S. A., and Collins, W. D.: Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models, J. Geophys. Res.-Atmos., 113, D13103, https://doi.org/10.1029/2008JD009944, 2008. a

Keinert, A., Spannagel, D., Leisner, T., and Kiselev, A.: Secondary Ice Production upon Freezing of Freely Falling Drizzle Droplets, J. Atmos. Sci., 77, 2959–2967, https://doi.org/10.1175/JAS-D-20-0081.1, 2020. a

Korolev, A.: Reconstruction of the sizes of spherical particles from their shadow images. Part I: Theoretical considerations, J. Atmos. Ocean. Tech., 24, 376–389, 2007. a

Korolev, A. and Field, P. R.: Assessment of the performance of the inter-arrival time algorithm to identify ice shattering artifacts in cloud particle probe measurements, Atmos. Meas. Tech., 8, 761–777, https://doi.org/10.5194/amt-8-761-2015, 2015. a

Korolev, A., Heckman, I., Wolde, M., Ackerman, A. S., Fridlind, A. M., Ladino, L. A., Lawson, R. P., Milbrandt, J., and Williams, E.: A new look at the environmental conditions favorable to secondary ice production, Atmos. Chem. Phys., 20, 1391–1429, https://doi.org/10.5194/acp-20-1391-2020, 2020. a, b, c

Ladino, L. A., Korolev, A., Heckman, I., Wolde, M., Fridlind, A. M., and Ackerman, A. S.: On the role of ice‐nucleating aerosol in the formation of ice particles in tropical mesoscale convective systems, Geophys. Res. Lett., 44, 1574–1582, https://doi.org/10.1002/2016GL072455, 2017. a

Lang, S. E., Tao, W.-K., Zeng, X., and Li, Y.: Reducing the biases in simulated radar reflectivities from a bulk microphysics scheme: Tropical convective systems, J. Atmos. Sci., 68, 2306–2320, 2011. a, b

Langley Research Center (LaRC), NASA: Cayenne GOES-13 Satellite Cloud Products Data, Version 1.0, UCAR/NCAR – Earth Observing Laboratory. https://doi.org/10.5065/D6NC5ZX6, 2016. a

Lawson, R. P., Angus, L. J., and Heymsfield, A. J.: Cloud particle measurements in thunderstorm anvils and possible weather threat to aviation, J. Aircraft, 35, 113–121, 1998. a

Leroy, D., Fontaine, E., Schwarzenboeck, A., Strapp, J. W., Lilie, L., Delanoe, J., Protat, A., Dezitter, F., and Grandin, A.: HAIC/HIWC Field Campaign - Specific Findings on PSD Microphysics in High IWC Regions from In Situ Measurements: Median Mass Diameters, Particle Size Distribution Characteristics and Ice Crystal Shapes, in: SAE 2015 International Conference on Icing of Aircraft, Engines, and Structures, 22–25 June 2015, Prague, Czech Republic, https://doi.org/10.4271/2015-01-2087, 2015. a, b

Leroy, D., Coutris, P., Emmanuel, F., Schwarzenboeck, A., Strapp, J. W., Lilie, L. E., Korolev, A., McFarquhar, G., Dezitter, F., and Grandin, A.: HAIC/HIWC field campaigns-Specific findings on ice crystals characteristics in high ice water content cloud regions, in: 8th AIAA Atmospheric and Space Environments Conference, 13–17 June 2016, Washington, D.C., USA, 4056, https://doi.org/10.2514/6.2016-4056, 2016a. a, b

Leroy, D., Fontaine, E., Schwarzenboeck, A., and Strapp, J.: Ice crystal sizes in high ice water content clouds. Part I: On the computation of median mass diameter from in situ measurements, J. Atmos. Ocean. Tech., 33, 2461–2476, 2016b. a

Leroy, D., Fontaine, E., Schwarzenboeck, A., Strapp, J. W., Korolev, A., McFarquhar, G., Dupuy, R., Gourbeyre, C., Lilie, L., Protat, A., Delanoe, J., Dezitter, F., and Grandin, A.: Ice crystal sizes in high ice water content clouds. Part II: Statistics of mass diameter percentiles in tropical convection observed during the HAIC/HIWC project, J. Atmos. Ocean. Tech., 34, 117–136, 2017. a, b, c, d, e

Locatelli, J. D. and Hobbs, P. V.: Fall speeds and masses of solid precipitation particles, J. Geophys. Res., 79, 2185–2197, https://doi.org/10.1029/JC079i015p02185, 1974. a

Lohmann, U., Roeckner, E., Collins, W. D., Heymsfield, A. J., McFarquhar, G. M., and Barnett, T. P.: The role of water vapor and convection during the Central Equatorial Pacific Experiment from observations and model simulations, J. Geophys. Res.-Atmos., 100, 26229–26245, 1995. a

Lynn, B. H., Khain, A. P., Dudhia, J., Rosenfeld, D., Pokrovsky, A., and Seifert, A.: Spectral (bin) microphysics coupled with a mesoscale model (MM5). Part I: Model description and first results, Mon. Weather Rev., 133, 44–58, 2005. a

Mason, J., Strapp, W., and Chow, P.: The ice particle threat to engines in flight, in: 44th AIAA Aerospace Sciences Meeting and Exhibit, 9–12 January 2006, Reno, Nevada, USA, 206, https://doi.org/10.2514/6.2006-206, 2006. a, b

Mason, J. G. and Grzych, M.: The challenges identifying weather associated with jet engine ice crystal icing, Tech. rep., SAE Technical Paper, SAE International, Chicago, Illinois, USA, 13–17 June 2011, https://doi.org/10.4271/2011-38-0094, 2011. a

McFarquhar, G. M. and Heymsfield, A. J.: Microphysical characteristics of three anvils sampled during the Central Equatorial Pacific Experiment, J. Atmos. Sci., 53, 2401–2423, 1996. a

McFarquhar, G. M. and Heymsfield, A. J.: Parameterization of tropical cirrus ice crystal size distributions and implications for radiative transfer: Results from CEPEX, J. Atmos. Aci., 54, 2187–2200, 1997. a

McFarquhar, G. M. and Heymsfield, A. J.: The definition and significance of an effective radius for ice clouds, J. Atmos. Sci., 55, 2039–2052, 1998. a

McFarquhar, G. M., Baumgardner, D., Bansemer, A., Abel, S. J., Crosier, J., French, J., Rosenberg, P., Korolev, A., Schwarzoenboeck, A., Leroy, D., et al.: Processing of ice cloud in situ data collected by bulk water, scattering, and imaging probes: fundamentals, uncertainties, and efforts toward consistency, Meteorological Monographs, 58, 11–1, 2017. a

Milbrandt, J. and Yau, M.: A multimoment bulk microphysics parameterization. Part I: Analysis of the role of the spectral shape parameter, J. Atmos. Sci., 62, 3051–3064, 2005. a

Milbrandt, J. A. and Morrison, H.: Parameterization of Cloud Microphysics Based on the Prediction of Bulk Ice Particle Properties. Part III: Introduction of Multiple Free Categories, J. Atmos. Sci., 73, 975–995, https://doi.org/10.1175/JAS-D-15-0204.1, 2016. a, b, c

Mitchell, D. L.: Use of Mass- and Area-Dimensional Power Laws for Determining Precipitation Particle Terminal Velocities, J. Atmos. Sci., 53, 1710–1723, https://doi.org/10.1175/1520-0469(1996)053<1710:UOMAAD>2.0.CO;2, 1996. a

Morrison, H. and Milbrandt, J. A.: Parameterization of Cloud Microphysics Based on the Prediction of Bulk Ice Particle Properties. Part I: Scheme Description and Idealized Tests, J. Atmos. Sci., 72, 287–311, https://doi.org/10.1175/JAS-D-14-0065.1, 2015. a, b, c, d, e

Morrison, H., Thompson, G., and Tatarskii, V.: Impact of cloud microphysics on the development of trailing stratiform precipitation in a simulated squall line: Comparison of one-and two-moment schemes, Mon. Weather Rev., 137, 991–1007, 2009. a, b, c

Nguyen, C. M., Wolde, M., and Korolev, A.: Determination of ice water content (IWC) in tropical convective clouds from X-band dual-polarization airborne radar, Atmos. Meas. Tech., 12, 5897–5911, https://doi.org/10.5194/amt-12-5897-2019, 2019. a

Protat, A., Delanoë, J., Strapp, J. W., Fontaine, E., Leroy, D., Schwarzenboeck, A., Lilie, L., Davison, C., Dezitter, F., Grandin, A., and Weber, M.: The Measured Relationship between Ice Water Content and Cloud Radar Reflectivity in Tropical Convective Clouds, J. Appl. Meteorol. Climatol., 55, 1707–1729, https://doi.org/10.1175/JAMC-D-15-0248.1, 2016. a, b

Qu, Z., Barker, H. W., Korolev, A. V., Milbrandt, J. A., Heckman, I., Bélair, S., Leroyer, S., Vaillancourt, P. A., Wolde, M., Schwarzenböck, A., Leroy, D., Strapp, J. W., Cole, J. N. S., Nguyen, L., and Heidinger, A.: Evaluation of a high-resolution numerical weather prediction model's simulated clouds using observations from CloudSat, GOES-13 and in situ aircraft, Q. J. Roy. Meteor. Soc., 144, 1681–1694, https://doi.org/10.1002/qj.3318, 2018. a, b, c, d

Ratvasky, T., Harrah, S., Strapp, J. W., Lilie, L., Proctor, F., Strickland, J., Hunt, P., Bedka, K., Diskin, G., Nowak, J. B., Bui, T. P., Bansemer, A., and Dumont C.: Summary of the High Ice Water Content (HIWC) RADAR Flight Campaigns, Tech. rep., SAE Technical Paper 2019-01-2025, NASA/TM-2020-220306, https://doi.org/10.4271/2019-01-2027, 2019. a

Ryzhkov, A. V. and Zrnić, D. S.: Radar Polarimetry for Weather Observations, Springer, Switzerland, https://doi.org/10.1007/978-3-030-05093-1, 2019. a, b

Ryzhkov, A. V., Snyder, J., Carlin, J. T., Khain, A., and Pinsky, M.: What Polarimetric Weather Radars Offer to Cloud Modelers: Forward Radar Operators and Microphysical/Thermodynamic Retrievals, Atmosphere, 11, 362, https://doi.org/10.3390/atmos11040362, 2020. a

Skamarock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Liu, Z., Berner, J., Wang, W., Powers, J. G., Duda, M. G., Barker, D. M., and Huang, X.-Y.: A description of the Advanced Research WRF Model version 4, Tech. rep., No. NCAR/TN-556+STR, National Center for Atmospheric Research Boulder, Colorado, USA, https://doi.org/10.5065/1dfh-6p97, 2019. a

Stanford, M. W., Varble, A., Zipser, E., Strapp, J. W., Leroy, D., Schwarzenboeck, A., Potts, R., and Protat, A.: A ubiquitous ice size bias in simulations of tropical deep convection, Atmos. Chem. Phys., 17, 9599–9621, https://doi.org/10.5194/acp-17-9599-2017, 2017. a

Strapp, J., Korolev, A., Ratvasky, T., Potts, R., Protat, A., May, P., Ackerman, A., Fridlind, A., Minnis, P., Haggerty, J., Riley, J., Lilie, L., and Isaac, G.: The High Ice Water Content (HIWC) study of deep convective clouds: Science and technical plan, Tech. rep., FAA Rep, DOT/FAA/TC-14/31, 105 pp., available at: http://www.tc.faa.gov/its/worldpac/techrpt/tc14-31.pdf (last access: 6 July 2020), 2016a. a

Strapp, J. W., Lilie, L. E., Ratvasky, T. P., Davison, C. R., and Dumont, C.: Isokinetic TWC Evaporator Probe: Development of the IKP2 and Performance Testing for the HAIC-HIWC Darwin 2014 and Cayenne Field Campaigns, in: 8th AIAA Atmospheric and Space Environments Conference, available at: https://arc.aiaa.org/doi/pdf/10.2514/6.2016-4059 (last access: 6 July 2020), 2016b. a, b

Strapp, J. W., Schwarzenboeck, A., Bedka, K., Bond, T., Calmels, A., Delanoë, J., Dezitter, F., Grzych, M., Harrah, S., Korolev, A., Leroy, D., Lilie, L., Mason, J., Potts, R., Protat, A., Ratvasky, T., Riley, J., and Wolde, M.: Comparisons of Cloud In Situ Microphysical Properties of Deep Convective Clouds to Appendix D/P Using Data from the High-Altitude Ice Crystals-High Ice Water Content and High Ice Water Content-RADAR I Flight Campaigns, SAE Int. J. Aerosp., 14, https://doi.org/10.4271/01-14-02-0007, 2021. a

Tewari, M., Chen, F., Wang, W., Dudhia, J., LeMone, M., Mitchell, K., Ek, M., Gayno, G., Wegiel, J., and Cuenca, R.: Implementation and verification of the unified NOAH land surface model in the WRF model, in: 20th conference on weather analysis and forecasting/16th conference on numerical weather prediction, American Meteorological Society: Seattle, WA, US., available at: https://ams.confex.com/ams/pdfpapers/69061.pdf (last access: 6 July 2020), 2004.  a

Thompson, G., Field, P. R., Rasmussen, R. M., and Hall, W. D.: Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: Implementation of a new snow parameterization, Mon. Weather Rev., 136, 5095–5115, 2008. a

Wolde, M., Nguyen, C., Korolev, A., and Bastian, M.: Characterization of the Pilot X-band radar responses to the HIWC environment during the Cayenne HAIC-HIWC 2015 Campaign, in: 8th AIAA Atmospheric and Space Environments Conference, American Institute of Aeronautics and Astronautics, Washington, D.C., https://doi.org/10.2514/6.2016-4201, 2016. a, b, c, d, e

Yost, C. R., Bedka, K. M., Minnis, P., Nguyen, L., Strapp, J. W., Palikonda, R., Khlopenkov, K., Spangenberg, D., Smith Jr., W. L., Protat, A., and Delanoe, J.: A prototype method for diagnosing high ice water content probability using satellite imager data, Atmos. Meas. Tech., 11, 1615–1637, https://doi.org/10.5194/amt-11-1615-2018, 2018. a

Yuter, S. E. and Houze Jr., R. A.: Three-dimensional kinematic and microphysical evolution of Florida cumulonimbus. Part II: Frequency distributions of vertical velocity, reflectivity, and differential reflectivity, Mon. Weather Rev., 123, 1941–1963, 1995. a

1

The Morrison scheme has an option to represent rimed ice as either hail or graupel which affects the fall speed and density, and here the option for graupel is chosen.

2

The percentage of underestimate/overestimate/underprediction/overprediction in this paper is the ratio of difference between simulation and observation to observation if not otherwise specified.

Short summary
Numerous small ice crystals in the tropical convective storms are difficult to detect and could be potentially hazardous for commercial aircraft. This study evaluated the numerical models against the airborne observations and investigated the potential cloud processes that could lead to the production of these large numbers of small ice crystals. It is found that key microphysical processes are still lacking or misrepresented in current numerical models to realistically simulate the phenomenon.
Altmetrics
Final-revised paper
Preprint