Dust aerosol has a widespread impact on air quality and visibility (Huang et al., 2008;Moulin et al., 1998) and modulates the Earth's radiation budget (Pachauri et al., 2014) via scattering and absorbing both shortwave and thermal infrared radiation (Xu et al., 2017), leading to large uncertainties in climate projection. The radiative forcing and spectral fingerprints of dust particles depend on their bulk optical properties, including the extinction coefficient, single-scattering albedo, and phase matrix. While these properties can be computed exactly based on the Lorenz-Mie theory for spherical particles (Mishchenko et al., 2002), most mineral dust particles have highly irregular shapes with a wide range of particle sizes, posing significant challenges to both remote sensing and climate modeling (Wang et al., 2003(Wang et al., , 2004.Various techniques have been developed for estimating non-spherical particle properties, such as the discrete dipole approximation (DDA;Draine & Flatau, 1994;Yurkin et al., 2007), the extended boundary condition method (EBCM) as an implementation of the T-matrix (Mishchenko et al., 1997;Mishchenko & Travis, 1994),