The Use of Superspheroids as Surrogates for Modeling Electromagnetic Wave Scattering by Ice Crystals

: Electromagnetic wave scattering by ice particles is commonly modeled by deﬁning representative habits, including droxtals, columns, plates, and aggregates, although actual particles in the atmosphere can be even much more complex. In this study, we examined a superspheroidal approximation method for modeling electromagnetic wave scattering by ice crystals. Superspheroid can be associated with a shape index ( SI ) deﬁned by the particle volume and average projected area. Corresponding to realistic ice crystals, suitable superspheroid models with the same SI (that means, identical volume and average projected area) and aspect ratio can be identiﬁed as surrogates for optical property calculations. We systematically compared the optical properties of ice crystals and superspheroids at 33 microwave bands in the range of 3–640 GHz and at three representative visible or infrared wavelengths (0.66, 2.13, and 11 µ m). It was found that the single-scattering properties of compact ice crystal habits and their superspheroidal model particles were quite close. For an aggregate with sparse distribution of elements, a superspheroid model produces relatively large errors because the aspect ratio may not be sufﬁcient to describe a particle shape. However, the optical similarity of a superspheroid and an aggregate is still encouraging.


Introduction
Cirrus clouds cover approximately 20% of the earth's surface, and their radiation characteristics depend on their microphysical properties, including the shape, refractive index, and size distribution, of ice crystal particles [1][2][3][4][5]. In addition to cirrus clouds, precipitation clouds or ice fogs also contain various ice crystals and snowflakes depending on the ambient temperature and moisture conditions [6,7]. Therefore, substantial research has been devoted to studying the microphysics of ice crystals and their optical properties to explore the impact of ice particles on the earth's radiation budget, weather, and climate system [8][9][10][11][12][13].
Canonical ice crystal habits include hexagonal columns, plates, and aggregates [3,14]. However, ice crystals in the atmosphere can be much more complicated than canonical models [8,[15][16][17]. Moreover, ice crystals have a wide range of particle sizes from a few microns to thousands of microns. Because of the diverse particle shapes and a wide range of particle sizes, it has been a great challenge to accurately compute the optical properties of ice particles for radiative transfer simulations, data assimilation, and remote sensing applications.
Early attempts at modeling electromagnetic wave scattering by ice crystals were based on an assumption of spherical particle shapes [1,12], followed by nonspherical approximations including simple geometries, such as spheroids and cylinders [18][19][20][21][22]. These explorations were greatly motivated by the robustness of the Lorenz-Mie theory [23] and the extended boundary condition method [24,25], along with partial motivation by possible optical similarities between model particles and realistic ice crystal habits. For example, Lee et al. [26] proposed a method to simulate the light scattering characteristics of ice crystals using a circular cylinder, considering that the sharp edges of hexagonal ice particles may not be a decisive factor in light scattering calculation in the near-infrared range. The results showed that hexagonal ice particles and circular cylinders have optical similarity in the near-infrared band. Mishchenko and Macke [27] exploited the optical similarity between hexagonal column and circular cylinder to determine how big the particle size should be to have 46 • halo.
However, these approximations have gradually disappeared because of the drawback of shape approximations being hard to evaluate, and a variety of computational methods have been developed that can deal efficiently with more realistic ice crystal habits. For small to moderate sizes of ice particles, numerous numerically exact methods have been developed, including the discrete dipole approximation [28,29], the finite-difference timedomain method [30,31], the pseudo-spectral time-domain method [32,33], the boundary element method [34], and the invariant imbedding T-matrix method (IITM) [35][36][37][38]. For large ice crystals, the geometric optics or physical optics methods are now common choices in research [39][40][41][42][43][44][45][46].
A few databases of the optical properties of ice crystals have been made available using the abovementioned computational methods. For example, Hess [47] used hexagonal plates and columns to establish a database of optical properties of ice crystals, including phase functions, extinction cross sections, single-scattering albedos, and asymmetry factors. Yang et al. [14] developed a comprehensive database of scattering and absorption characteristics of ice crystals in the spectrum range of 0.2-100 µm, and 11 randomly oriented ice crystal habits with maximum dimensions of 2-10,000 µm were considered. Bi and Yang [48] further improved the accuracy of this database. Hong [49] evaluated the scattering and absorption characteristics of nonspherical ice crystals in ice clouds using the discrete dipole approximation method in the microwave band of 89-340 GHz. It was found that the singlescattering characteristics of ice crystals were sensitive to their shapes. Kim [50] reported single-scattering properties of randomly oriented snow particles at microwave frequencies.
Liu [51] established the single-scattering characteristics of 11 ice crystals with different shapes at 15-340 GHz. Ding et al. [52] developed a database of 12 kinds of scattering characteristics of ice crystals in the microwave spectral region and compared the conditions at 160, 200, 230, and 270 K. It was confirmed that the single-scattering characteristics of ice crystals in the microwave range depend on the temperature. More recently, Eriksson et al. [53] developed a single-scattering property database at microwave and submillimeter wavelengths with 34 ice habits. It is expected that more precise results of ice crystal properties can be obtained using more accurate ice crystal shape characterizations, although trial and tuning are needed in most applications. This paper presents a new approach to modelling electromagnetic scattering by ice crystals. Specifically, we explore the use of superspheroids as surrogates of realistic ice habits for computing optical properties, including the extinction efficiency, single-scattering albedo, and scattering matrix. The present approach is similar to earlier optical modeling methods, such as spherical, spheroidal, and cylindrical approximations. However, superspheroids have one more degree of freedom than spheroids, which has been found to be useful and convenient for parameterizing the optical properties of nonspherical aerosol particles [54][55][56][57], even though the geometry of superspheroids can be quite different from realistic particle shapes. Physically, superspheroids can be understood as simplified models that capture optically important parameters of ice crystals (the volume, area, and aspect ratio) for optical modeling. The present research is also influenced by previous studies in the use of irregular particle models (e.g., Voronoi particles [58,59] and Gaussian random spheres [60]) to represent realistic ice crystals.
Since comprehensive light-scattering databases are available from the ultraviolet to microwave bands, a question naturally arises as to whether it is still necessary to study superspheroidal approximation. As a matter of fact, databases only include a set of welldefined ice crystal geometries. Realistic ice particles can be much more complex, and it is challenging to manage databases with increasing number of ice crystal habits. This may also induce additional difficulty in habit selection for specific applications. In the case of ice crystals, the advantage of superspheroids is that the variation of particle shape is continuous. The model may be useful in estimating random errors associated with particle complexities and offer a complementary choice for optical modeling of ice crystals.
In view of this, this paper proposes identifying suitable superspheroidal models as surrogates of well-defined ice crystal models for computing their optical properties. We only considered randomly oriented particles in this study. The modeling capability was tested and evaluated by comparing the extinction efficiency factor, scattering efficiency factor, and scattering matrix against available databases. In particular, we highlight a new definition of the shape index (SI) that is valuable for finding an appropriate superspheroid with close optical similarity to a predefined ice crystal.
The remainder of this paper is organized as follows. Section 2 discusses the superspheroidal model and definition of related parameters. Sections 3 and 4 compare the computational results of superspheroids and ice crystals. The comparisons were systematically carried out in the microwave band and at a few visible and infrared wavelengths. A discussion is given in Section 5, and Section 6 concludes this study.

Superspheroidal Models and Shape Index
The superspheroidal equation is defined as follows [61,62]: where a, b, and c are semiaxes aligned with the x, y, and z directions in the Cartesian coordinate system, respectively, and n is referred to as the roundness parameter. Following previous studies for modeling the optical properties of dust and sea salt aerosols [55,56], we mainly studied the superspheroidal particles by assuming a = b. The particle shape can be changed by controlling the two parameters of aspect ratio (ξ) and roundness parameter (n). The aspect ratio (ξ) is defined as follows: ( For the roundness parameter n < 1, a factor 2 (1−n)/2 is introduced because the largest horizontal dimension is larger than a.
From a physical perspective of electromagnetic wave scattering by particles, the particle can be thought of as being composed of a number of electrical dipoles [23]. The total number of dipoles is determined by the wavelength and the particle volume, and the spatial distribution of dipoles is characterized by the particle shape. As the particle size increases, the extinction cross sections are fundamentally determined by the average projected area [23], which also affects the Fraunhofer diffraction. Therefore, microphysical properties, such as the aspect ratio (as a quantity to characterize the dipole distribution), volume, and average projected area of particles, have an important impact on the optical properties of particles. Thus, to identify a surrogate for a given ice crystal particle, the most suitable candidate is the superspheroid that has the same aspect ratio, volume, and projected area as the ice crystal. For this purpose, a shape index SI is defined as follows: where V is the volume and S is the average projected area of the particle. According to Equation (3), if SI is the same for two different particles, then their volume and average projected area are identical. Figure 1 shows the variation of SI with respect to the aspect ratio and the roundness of particles. The SI of spherical particles is 1 and that of superspheroidal particles is less than 1. In addition, for the same SI, superspheroidal models also allow additional freedom to change the aspect ratio. Therefore, with a predefined ice crystal model, by defining SI, superspheroidal models can be changed to approach the aspect ratio, volume, and average projected area of ice crystal particles concurrently. Figure 2 shows representative ice crystals and their superspheroidal counterparts. For the shapes of droxtal, 8-column aggregate, 5-plate aggregate, and 10-plate aggregate with the same aspect ratio and SI parameters, the superspheroids simulated the ice crystals. For the column and plate, their shapes changed with the size of the particles, so the superspheroids corresponding to three ice crystal particles with different sizes are listed here. The SIs of ice crystals and the aspect ratios and roundness parameters of superspheroid counterparts are given in Table 1.    , and 10-plate (f) aggregates are simulated by superspheroids of a single shape. The shape of superspheroids corresponding to column (b) and plate (d) changes with the particle size D. D is defined as the maximum dimension and can be found in Yang et al. [14]. Here, D is taken as 500, 3500, and 8000 μm for example.

Optical Properties
Ice clouds in the atmosphere present as an ensemble of ice crystals of different sizes. Therefore, it is necessary to compute the bulk scattering characteristics of ice crystals per unit volume. The scattering coefficient bsca and extinction coefficient bext are given as follows [23]: where Csca is the scattering cross section, which is the product of the average projected area S and the scattering efficiency factor Qsca. Cext and Qext are the extinction cross section and the extinction efficiency factor, respectively. D is the maximum dimension of ice crystals, and the size distribution n(D)dD represents the particle number concentration within (D,D+dD) per volume. In computations, D1 = 2 μm and D2 = 10,000 μm. For randomly oriented superspheroids, the relationship between the Stokes vectors of an incident wave , , , and a scattering wave , , , is given by the scattering matrix [23]: , and 10-plate (f) aggregates are simulated by superspheroids of a single shape. The shape of superspheroids corresponding to column (b) and plate (d) changes with the particle size D. D is defined as the maximum dimension and can be found in Yang et al. [14]. Here, D is taken as 500, 3500, and 8000 µm for example.

Optical Properties
Ice clouds in the atmosphere present as an ensemble of ice crystals of different sizes. Therefore, it is necessary to compute the bulk scattering characteristics of ice crystals per unit volume. The scattering coefficient b sca and extinction coefficient b ext are given as follows [23]: where C sca is the scattering cross section, which is the product of the average projected area S and the scattering efficiency factor Q sca . C ext and Q ext are the extinction cross section and the extinction efficiency factor, respectively. D is the maximum dimension of ice crystals, and the size distribution n(D)dD represents the particle number concentration within (D, D+dD) per volume. In computations, D 1 = 2 µm and D 2 = 10,000 µm. For randomly oriented superspheroids, the relationship between the Stokes vectors of an incident wave [I, Q, U, V] inc and a scattering wave [I, Q, U, V] sca is given by the scattering matrix [23]: where θ is the scattering angle in the range of 0 to 180 • . Because of the symmetry of the superspheroid, eight off-diagonal elements of the scattering matrix are all 0. The incident and scattering Stokes vectors are defined relative to the scattering plane, where I is the total intensity, Q and U describe the linear polarization intensity, and V describes the circular polarization intensity. The bulk scattering matrix is given by [63]: where P ij is an element in the scattering matrix, C sca is the scattering cross section, and b sca is the scattering coefficient. In this paper, we used the standard gamma distribution [63]: where D is the maximum dimension, D m is the effective diameter, and v e is the effective variance. The superspheroidal approximation was examined by referencing the optical properties of droxtals, columns, plates, 8-column aggregates, 5-plate aggregates, and 10-plate aggregates from optical property databases [14,48,52]. The optical properties of the superspheroidal particles were computed from the IITM and the improved geometric optics method (IGOM) [44]. The applicability regimes of superspheroidal approximations can be examined by comparing the extinction efficiency, scattering efficiency, asymmetry factor, and phase matrix of well-defined ice crystal habits and their superspheroid counterparts.

Microwave Bands
The experimental temperature was assumed to be 160 K, and the frequency range was 3-640 GHz with 33 bands. The temperature-dependent complex refractive indices were taken from Iwabuchi and Yang [64]. For each ice crystal model, a total of 24 particle sizes in the range of 2-10,000 µm were selected. The IITM was used for all the computations in the microwave bands. The optical properties examined included the extinction efficiency factor Q ext , scattering efficiency factor Q sca , phase matrix element P ij , and asymmetry factor g.

Superspheroids for Droxtal
According to the SI and the aspect ratio values of a droxtal (see Table 1), a superspheroidal model with an aspect ratio of 0.9 and roundness parameter of 1.6 was identified. We note that superspheroids with roundness parameters of 0.5 and 1.6 both satisfy the condition. The results of two superspheroids were found to be close. The roundness of 1.6 was selected because of a slightly better comparison to the results of droxtal.
Although the droxtal particle and the superspheroid particle have obviously different geometries, they have identical volume, average projected area, and aspect ratios. Figure 3 shows the simulation results of the extinction efficiency, scattering efficiency, and asymmetry factor at 33 bands in the frequency range of 3-640 GHz. The maximum dimensions were 3000 and 9000 µm (the sizes were selected for the method examination; note that realistic droxtal ice crystal may not have such large particle sizes). We found that the trends of these results were the same as those of the superspheroid, and the values were very close. Because the particle size was fixed, the size parameter increased with the frequency. Both the droxtals and superspheroids showed similar resonance effect as the size parameter varied. At the lower end of the frequency, the size parameter was quite small and the shape effect was negligible, so the agreement was better. With the increase of size parameters, the particle shape was expected to exert a greater influence on the optical properties of particles. However, in the comparisons, the differences between droxtal and superspheroids were small. This means that the volume, projected area, and aspect ratio are optically important parameters that determine the optical properties. Therefore, it is reasonable to use superspheroids as surrogates of droxtal ice particles. Next, we examined the scattering matrix. Figure 4 compares six nonzero scatteringmatrix elements of droxtals and superspheroids. The frequency was 220 GHz, and the maximum dimensions were 500, 4000, and 9000 μm. For the first column, where the maximum dimension was 500 μm, the nonsphericity effect was almost invisible, as shown by At the lower end of the frequency, the size parameter was quite small and the shape effect was negligible, so the agreement was better. With the increase of size parameters, the particle shape was expected to exert a greater influence on the optical properties of particles. However, in the comparisons, the differences between droxtal and superspheroids were small. This means that the volume, projected area, and aspect ratio are optically important parameters that determine the optical properties. Therefore, it is reasonable to use superspheroids as surrogates of droxtal ice particles.
Next, we examined the scattering matrix. Figure 4 compares six nonzero scatteringmatrix elements of droxtals and superspheroids. The frequency was 220 GHz, and the maximum dimensions were 500, 4000, and 9000 µm. For the first column, where the maximum dimension was 500 µm, the nonsphericity effect was almost invisible, as shown by element P 22 /P 11 (=1), which is an indicator of nonsphericity. The scattering matrix elements from the two different models were almost identical and showed Rayleigh scattering behavior.
Remote Sens. 2021, 13, 1733 9 of 34 element P22/P11 (=1), which is an indicator of nonsphericity. The scattering matrix elements from the two different models were almost identical and showed Rayleigh scattering behavior.
For the second and third columns, as the size parameter increased, all scattering matrix elements showed strong oscillations, and the influence of particle shape on the scattering characteristics became visible. The scattering matrix elements of superspheroids and droxtals showed some differences. However, the trends of the scattering matrix elements with respect to the scattering angle for the two geometries were obviously consistent. This means that the superspheroid model offers a good approximation of droxtal geometries for computing the scattering phase matrix. Note that the spherical approximation failed to produce most features (see Yang et al. [65]). To compute the bulk scattering optical properties, a standard gamma distribution with an effective maximum dimension of 4000 μm and an effective variance of 0.1 was For the second and third columns, as the size parameter increased, all scattering matrix elements showed strong oscillations, and the influence of particle shape on the scattering characteristics became visible. The scattering matrix elements of superspheroids and droxtals showed some differences. However, the trends of the scattering matrix elements with respect to the scattering angle for the two geometries were obviously consistent. This means that the superspheroid model offers a good approximation of droxtal geometries for computing the scattering phase matrix. Note that the spherical approximation failed to produce most features (see Yang et al. [65]).
To compute the bulk scattering optical properties, a standard gamma distribution with an effective maximum dimension of 4000 µm and an effective variance of 0.1 was used to integrate the single-scattering properties. The results in the 220 GHz band are shown in Figure 5. As expected, most oscillation features in Figure 4 were cancelled out, resulting in a closer agreement between the results computed for the droxtal and superspheroid models. Note that P 22 /P 11 values deviated more from unity for the superspheroid than the droxtal. The scattering and extinction coefficients for the droxtal and the superspheroid were close (see Table 2).
Remote Sens. 2021, 13, 1733 10 of 34 used to integrate the single-scattering properties. The results in the 220 GHz band are shown in Figure 5. As expected, most oscillation features in Figure 4 were cancelled out, resulting in a closer agreement between the results computed for the droxtal and superspheroid models. Note that P22/P11 values deviated more from unity for the superspheroid than the droxtal. The scattering and extinction coefficients for the droxtal and the superspheroid were close (see Table 2).

Superspheroids for Columns
The aspect ratio of column particles varies with the particle size. The relationship between the aspect ratio and size is given as follows [66,67]: where a is the semiwidth, and L is the length. As such, superspheroids for different particle sizes are also different. Figure 6 shows a comparison of the extinction efficiency factor, scattering efficiency factor, and asymmetry factor computed for columns and superspheroids in the frequency range of 3-640 GHz. Maximum dimensions of 3000 and 7000 µm were chosen for comparison. Close agreement was again obtained between the results from the two different models.

Superspheroids for Columns
The aspect ratio of column particles varies with the particle size. The relationship between the aspect ratio and size is given as follows [66,67]: where a is the semiwidth, and L is the length. As such, superspheroids for different particle sizes are also different. Figure 6 shows a comparison of the extinction efficiency factor, scattering efficiency factor, and asymmetry factor computed for columns and superspheroids in the frequency range of 3-640 GHz. Maximum dimensions of 3000 and 7000 μm were chosen for comparison. Close agreement was again obtained between the results from the two different models.  The scattering matrix was also compared for ice columns and superspheroids with multiple sizes, as shown in Figure 7. With the increase of particle size, the error increased gradually, but the trends were consistent. The simulation error of phase function P 11 was visibly small for most particle sizes. Figure 8 shows the bulk scattering matrix elements at 220 GHz. Similar to the droxtal case, after the size integration, the differences in the scattering matrix elements were reduced to some extent, and better agreement was reached. The scattering and extinction coefficients for the column and the superspheroid were close (see Table 2).
Remote Sens. 2021, 13, 1733 12 of 34 The scattering matrix was also compared for ice columns and superspheroids with multiple sizes, as shown in Figure 7. With the increase of particle size, the error increased gradually, but the trends were consistent. The simulation error of phase function P11 was visibly small for most particle sizes. Figure 8 shows the bulk scattering matrix elements at 220 GHz. Similar to the droxtal case, after the size integration, the differences in the scattering matrix elements were reduced to some extent, and better agreement was reached. The scattering and extinction coefficients for the column and the superspheroid were close (see Table 2).

Superspheroids for Plates
The relationship between the aspect ratio of a plate and its maximum dimension is given as follows [68]: Figure 9 compares the efficiency factors and asymmetry factor of superspheroids and plates with different sizes at 3-640 GHz. As can be seen, for D of 3000 μm, the efficiency factor of superspheroids was a little bit smaller in the frequency band of 500-640 GHz, and the asymmetry factor was slightly larger in the frequency band of 3-100 GHz. For D of 7000 μm, the efficiency factor of the superspheroids was lower at 300-400 GHz, but the asymmetry factor was fairly close. In general, the two results showed good agreement. Figure 10 compares the results of scattering matrix elements for small, medium, and large sizes in the 220 GHz band. It was found that the simulation of P11 and the other elements from superspheroidal models also fit the ice plate results well. Taking Dm = 4000 μm and ve = 0.1 as an example, the bulk scattering results at 220 GHz are shown in Figure 11. The scattering and extinction coefficients for the plate and superspheroid are given in Table 2.

Superspheroids for Plates
The relationship between the aspect ratio of a plate and its maximum dimension is given as follows [68]: Figure 9 compares the efficiency factors and asymmetry factor of superspheroids and plates with different sizes at 3-640 GHz. As can be seen, for D of 3000 µm, the efficiency factor of superspheroids was a little bit smaller in the frequency band of 500-640 GHz, and the asymmetry factor was slightly larger in the frequency band of 3-100 GHz. For D of 7000 µm, the efficiency factor of the superspheroids was lower at 300-400 GHz, but the asymmetry factor was fairly close. In general, the two results showed good agreement. Figure 10 compares the results of scattering matrix elements for small, medium, and large sizes in the 220 GHz band. It was found that the simulation of P 11 and the other elements from superspheroidal models also fit the ice plate results well. Taking D m = 4000 µm and v e = 0.1 as an example, the bulk scattering results at 220 GHz are shown in Figure 11. The scattering and extinction coefficients for the plate and superspheroid are given in Table 2.

Superspheroids for an 8-Column Aggregate
The 8-column aggregate is a more complex particle shape than droxtals and hexagonal columns/plates. It is currently used as the ice crystal model in the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 product [69]. Considering the aspect ratio, average projected area, and volume, the selected superspheroid model had an aspect ratio of 2 and a roundness parameter of 2.7. The scattering characteristics of superspheroids and the 8-column aggregate at each frequency band were compared. Figure 12 compares the extinction, scattering, and asymmetry factor of the aggregate and a single superspheroid. In contrast to the previous comparisons, the differences were relatively large. The fundamental reason may be that the aspect ratio is not sufficient to describe the particle shape. The mutual interactions between particle elements may also be important. However, from comparison of the scattering matrix elements (Figure 13), the optical similarity is still encouraging. Figure 14 shows the bulk scattering characteristics of the superspheroids and the aggregate particle at 220 GHz. In this figure, D m is 7000 µm and v e is 0.1. It can be seen that the overall trend of the two results was consistent, and the error was smaller. The scattering and extinction coefficients for the 8-column aggregate and the superspheroid counterpart are given in Table 2.

Superspheroids for an 8-Column Aggregate
The 8-column aggregate is a more complex particle shape than droxtals and hexagonal columns/plates. It is currently used as the ice crystal model in the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 product [69]. Considering the aspect ratio, average projected area, and volume, the selected superspheroid model had an aspect ratio of 2 and a roundness parameter of 2.7. The scattering characteristics of superspheroids and the 8-column aggregate at each frequency band were compared. Figure 12 compares the extinction, scattering, and asymmetry factor of the aggregate and a single superspheroid. In contrast to the previous comparisons, the differences were relatively large. The fundamental reason may be that the aspect ratio is not sufficient to describe the particle shape. The mutual interactions between particle elements may also be important. However, from comparison of the scattering matrix elements (Figure 13), the optical similarity is still encouraging. Figure 14 shows the bulk scattering characteristics of the superspheroids and the aggregate particle at 220 GHz. In this figure, Dm is 7000 μm and ve is 0.1. It can be seen that the overall trend of the two results was consistent, and the error was smaller. The scattering and extinction coefficients for the 8-column aggregate and the superspheroid counterpart are given in Table 2.

Superspheroids for Sparse Aggregate Shapes: 5-Plate and 10-Plate Aggregates
The optical properties of 5-plate and 10-plate aggregates are shown in Figure 15. Different from the 8-column aggregate, the two plate aggregates were sparsely distributed. According to Fox [70], the mass-dimension relationship of the 8-column aggregate is quite different from other models and observations, which emphasizes the necessity of studying sparse ice aggregates for optical simulations because their mass-dimension relationship could be more reasonable. The optical similarities of these sparse shapes and their superspheroidal counterparts might not be as good as those of compact ice crystals because the use of the aspect ratio to describe the particle shape could be oversimplified. Although each element can be modeled as a superspheroid, the aggregation process could complicate the calculation of optical properties, which is against the original motivation of using a simple modeling approach.
According to the SI and the aspect ratio of the 5-plate aggregate, a superspheroid model with an aspect ratio of 3.4 and roundness of 2.2 was selected. For the 10-plate aggregate, the superspheroid model's aspect ratio was 4.0, and the roundness was 2.7. Figure 15 compares Qsca, Qext, and g between the superspheroids and the two aggregate models in 33 bands from 3 to 640 GHz. The extinction and scattering coefficients for the two different shapes had large differences, particularly for larger frequencies. However, the asymmetry parameters of the aggregates and their superspheroidal counterparts were close. Similar to the 8-column aggregate case, the scattering matrix elements of the superspheroid model and the plate aggregates resembled each other (see Figure 16). The phase function P 11 of the superspheroid was smoother than that of the ice crystal model, and the backward scattering was lower.

Superspheroids for Sparse Aggregate Shapes: 5-Plate and 10-Plate Aggregates
The optical properties of 5-plate and 10-plate aggregates are shown in Figure 15. Different from the 8-column aggregate, the two plate aggregates were sparsely distributed. According to Fox [70], the mass-dimension relationship of the 8-column aggregate is quite different from other models and observations, which emphasizes the necessity of studying sparse ice aggregates for optical simulations because their mass-dimension relationship could be more reasonable. The optical similarities of these sparse shapes and their superspheroidal counterparts might not be as good as those of compact ice crystals because the use of the aspect ratio to describe the particle shape could be oversimplified. Although each element can be modeled as a superspheroid, the aggregation process could complicate the calculation of optical properties, which is against the original motivation of using a simple modeling approach. According to the SI and the aspect ratio of the 5-plate aggregate, a superspheroid model with an aspect ratio of 3.4 and roundness of 2.2 was selected. For the 10-plate aggregate, the superspheroid model's aspect ratio was 4.0, and the roundness was 2.7. Figure  15 compares Qsca, Qext, and g between the superspheroids and the two aggregate models Figure 15. Comparison of the extinction efficiency, scattering efficiency, and asymmetry factor of a plate aggregate ice crystal and its corresponding superspheroid in the frequency range of 3-640 GHz. From left to right: 5-plate aggregate with maximum dimension of 3000 µm, 5-plate aggregate with maximum dimension of 7000 µm, 10-plate aggregate with maximum dimension of 3000 µm, and 10-plate aggregate with maximum dimension of 7000 µm.  The bulk scattering matrix was computed using the gamma size distribution with D m of 7000 µm and v e of 0.1 for the superspheroids and the two aggregate shapes, as shown in Figure 17. Obviously, the oscillation seen in Figure 16 was smoothed after integration. The −P 12 /P 11 also had a large difference in the case of 10-plate aggregates. The scattering and extinction coefficients for plate aggregates and their superspheroid counterparts are given in Table 2. Note that all the above simulations were in the microwave frequency range. The largest size parameter of all ice crystals was located in the resonant scattering region. In the shortwave regions, the size parameter could be much larger. The shape effect could be pronounced in geometric optics regimes. We will examine the superspheroidal approximation in a few short wavelengths in the next section.

Visible and Infrared Wavelengths
We considered three wavelengths (0.66, 2.13, and 11 μm), at which the refractive indices were, 1.3078 + 1.66 × 10 , 1.2677 + 5.255 × 10 , and 1.0086 + 2.478 × 10 , respectively. For each ice crystal model, a total of 45 particle sizes in the range of 2-100 μm were selected. According to the data library from Yang et al. [14], the IITM was used for small size parameters (<~20) and the IGOM (>~20) was used for large size parameters. For the IGOM calculations, three roughness conditions were considered. We followed the surface roughness definition given in [14]; the roughness parameters (σ) were 0, 0.03, and 0.5, associated with smooth, weakly roughened, and strongly roughened ice crystals, Note that all the above simulations were in the microwave frequency range. The largest size parameter of all ice crystals was located in the resonant scattering region. In the shortwave regions, the size parameter could be much larger. The shape effect could be pronounced in geometric optics regimes. We will examine the superspheroidal approximation in a few short wavelengths in the next section.

Visible and Infrared Wavelengths
We considered three wavelengths (0.66, 2.13, and 11 µm), at which the refractive indices were, 1.3078 + i1.66 × 10 −8 , 1.2677 + i5.255 × 10 −4 , and 1.0086 + i2.478 × 10 −1 , respectively. For each ice crystal model, a total of 45 particle sizes in the range of 2-100 µm were selected. According to the data library from Yang et al. [14], the IITM was used for small size parameters (<~20) and the IGOM (>~20) was used for large size parameters. For the IGOM calculations, three roughness conditions were considered. We followed the surface roughness definition given in [14]; the roughness parameters (σ) were 0, 0.03, and 0.5, associated with smooth, weakly roughened, and strongly roughened ice crystals, respectively. Note that the calculation of IITM did not consider the surface roughness because the size parameter was relatively smaller. Figure 18 shows the comparison of phase matrix elements of superspheroid and 8-column aggregate at 0.66 µm. The maximum dimension of the ice particle was chosen as 100 µm. The first column shows the results of the smooth particle (σ = 0). Obviously, the phase function of 8-column aggregate had ice halos at 22 and 46 degrees, which were invisible in the case of superspheroidal particle. The difference is understandable because the superspheroidal particle does not have similar abjoint facets of ice crystals causing the geometric optics halo features. However, as the degree of surface roughness increased, the results of superspheroid and 8-column aggregate were closer. Because the roughened ice crystals have been proven to be more reasonable in developing satellite cloud product, this comparison gives us more confidence in applying the superspheroidal approximation in practical applications. In addition, strong halo peaks in the phase functions stemming from light scattering by smooth ice particles do not appear frequently in reality [27].   Figure 19 shows the phase matrix elements for a droxtal (the first column), 8-column aggregate (the second column), and 5-plate aggregate (the third column) at wavelengths of 0.66 µm. We only considered the roughness parameter of 0.5. The results with a maximum dimension of 100 µm were again taken as an example. Figures 20 and 21 are similar to Figure 19 but for the wavelength of 2.13 and 11 µm. It was found that the superspheroids had much better performance in the droxtal case, and the two curves almost coincided. For the 8-column aggregate, the results of P 11 and −P 12 /P 11 at 0.66 and 2.13 µm were better, and the P 22 /P 11 at 60-150 • were slightly higher than those of the ice crystal model. When the wavelength was 11 µm, the fitting result of the ice crystal model was smoother, but the trend of the two fitting results was similar. For the 5-plate aggregate, the fitting results were better at 0.66 and 2.13 µm, but the P 44 /P 11 of superspheroids were slightly lower than that of the ice crystal model, especially for the backscattering (120 to 180 • ). Note that particles were highly absorptive at 11 µm, but the surface roughness had almost no effect.    We believe that SI and aspect ratio capture essential factors that determine the optical properties. The shape difference of ice crystals and superspheroids could cause additional differences of the optical properties of large particle sizes. However, the condition of surface roughness can inhibit the special morphology effect. Thus, it is expected that better comparison can be made when the surface roughness condition is applied. When the particle is highly absorbing, the diffraction and external reflection dominate the scattering pattern. An observable difference between an ice crystal aggregate and superspheroids is understandable because of distinct multiple external reflections.
The extinction and scattering efficiency factors and asymmetry factor of three shapes are compared in Figures 22-24. Figure 22 shows the IITM results for small-sized smooth particles. This comparison is similar to those in the microwave regions because size parameters are close. The size parameter is defined as kD, where k is the wavenumber. Note that the oscillation of the extinction curve fundamentally stems from the interference between the diffraction and the forward transmission. For a single and compact particle, it is easier to see the oscillation features (see the left droxtal panel). Figure 23 shows the results of large ice crystals computed from the IGOM at the wavelength of 0.66 µm. Two roughness conditions (0 and 0.5) were considered for comparison. As expected from the previous phase matrix comparisons, the optical similarity between ice crystals and superspheroids was much better for roughened surfaces (the lower three rows) than smooth surfaces (the upper three rows). Figure 24 is similar to Figure 23  The extinction and scattering efficiency factors and asymmetry factor of three shapes are compared in Figures 22-24. Figure 22 shows the IITM results for small-sized smooth particles. This comparison is similar to those in the microwave regions because size parameters are close. The size parameter is defined as kD, where k is the wavenumber. Note that the oscillation of the extinction curve fundamentally stems from the interference between the diffraction and the forward transmission. For a single and compact particle, it is easier to see the oscillation features (see the left droxtal panel). Figure 23 shows the results of large ice crystals computed from the IGOM at the wavelength of 0.66 μm. Two roughness conditions (0 and 0.5) were considered for comparison. As expected from the previous phase matrix comparisons, the optical similarity between ice crystals and superspheroids was much better for roughened surfaces (the lower three rows) than smooth surfaces (the upper three rows). Figure 24 is similar to Figure 23 but for the wavelength 2.13 μm. Similar conclusions can be obtained from Figures 23 and 24.  Figure 22. Comparison of the extinction efficiency, scattering efficiency, and asymmetry factor of three ice crystal models (droxtal (left), 8-column aggregate (middle), and 5-plate aggregate (right)) and their superspheroidal counterparts at the wavelength of 0.66 μm. The particle surface is smooth (σ = 0). Figure 22. Comparison of the extinction efficiency, scattering efficiency, and asymmetry factor of three ice crystal models (droxtal (left), 8-column aggregate (middle), and 5-plate aggregate (right)) and their superspheroidal counterparts at the wavelength of 0.66 µm. The particle surface is smooth (σ = 0). Figure 23. Comparison of the extinction efficiency, scattering efficiency, and asymmetry factor of three ice crystal models (droxtal (left), 8-column aggregate (middle), and 5-plate aggregate (right)) and their superspheroidal counterparts at the wavelength of 0.66 μm. The size parameter is from 55 to 620. For the first three rows, σ = 0, and for the last three rows, σ = 0.5.

Discussions
From the comparisons presented, it is encouraging to find that the single-scattering optical properties of ice crystals can be reasonably reproduced using superspheroidal approximations. Different from the conventional method with the definition of realistic ice crystals shapes, the superspheroid model simplifies the complex nonspherical model by using two parameters: the aspect ratio and the roundness parameter. A continuous change of the shape parameters can be used to construct a shape space with a large range of SIs to represent the complexity of realistic ice crystals, which might not be resolved by a few well-defined ice particle shapes. By comparing the scattering efficiency factor, extinction efficiency factor, asymmetry factor, and scattering matrix elements, we confirmed the promising applicability of using superspheroids as surrogates for compact ice crystals. Although for sparse aggregates, the extinction and scattering efficiencies factors in a resonance region could be underestimated, the asymmetry factor and scattering matrix can be approximated with reasonable accuracy.
The aggregation shapes of ice crystals could be much more complex. There may be no aggregate model that can adequately represent realistic ice crystals shapes. In this re-

Discussions
From the comparisons presented, it is encouraging to find that the single-scattering optical properties of ice crystals can be reasonably reproduced using superspheroidal approximations. Different from the conventional method with the definition of realistic ice crystals shapes, the superspheroid model simplifies the complex nonspherical model by using two parameters: the aspect ratio and the roundness parameter. A continuous change of the shape parameters can be used to construct a shape space with a large range of SIs to represent the complexity of realistic ice crystals, which might not be resolved by a few well-defined ice particle shapes. By comparing the scattering efficiency factor, extinction efficiency factor, asymmetry factor, and scattering matrix elements, we confirmed the promising applicability of using superspheroids as surrogates for compact ice crystals. Although for sparse aggregates, the extinction and scattering efficiencies factors in a resonance region could be underestimated, the asymmetry factor and scattering matrix can be approximated with reasonable accuracy.
The aggregation shapes of ice crystals could be much more complex. There may be no aggregate model that can adequately represent realistic ice crystals shapes. In this regard, superspheroidal particles with different SIs and aspect rations can be used for sensitivity studies. However, the current approach does not diminish the importance of existing databases based on well-defined ice habits. Databases for canonical ice crystal habits represent the true values for specific habits. Furthermore, it is impossible to use superspheroids for producing particular optical phenomena for smooth ice crystals, such as ice halos, although a higher degree of surface roughness generally results in a smoother phase function.

Conclusions
In this paper, superspheroid models are proposed as possible surrogates for representative ice crystal habits to compute the optical properties. The single-scattering properties of ice crystals and superspheroids, including the extinction/scattering efficiency factor, asymmetry factor, and the scattering matrix elements, were systematically compared for representative ice habits, including droxtals, columns, plates, 8-column aggregates, 5-plate aggregates, and 10-plate aggregates. The shape index (SI) was found to be quite useful for finding a superspheroid that has high optical similarity to realistic ice crystal shapes.
In both the microwave bands and the visible and the infrared wavelengths, a suitable superspheroid model was identified for each ice crystal, and the single-scattering properties of the superspheroid were close or similar to those of a realistic ice crystal habit.
Given the large range of SIs and aspect ratios of superspheroids, the proposed approach is highly promising for modeling the optical properties of ice crystals for atmospheric radiative transfer, remote sensing, and radiative forcing studies. In particular, for realistic ice crystals that cannot be represented by regular ice crystal habits, superspheroids could be used for sensitivity studies and random error estimations. It is believed that a database of superspheroids for ice crystal modeling could provide another valuable alternative for atmospheric scattering and radiative transfer studies.
Superspheroids have already been shown to be powerful for modeling aerosols like dust and sea salt [55,56]. This study indicates that the superspheroidal model may provide a more generalized approach for modeling most atmospheric particles. However, they should not be considered as models that can always perfectly model ice crystals. The performance is expected to worsen as the particle shape becomes sparsely distributed. Thus, some amount of trial and error is unavoidable.