When light (a laser beam for OAPs) illuminates a cloud particle, a shadow image can be observed on a screen at the rear of the particle. The formed image depends on the diffraction, refraction, and transmission of light by the particle. As a first approximation, it is convenient to neglect the refraction and transmission of light by the particle, i.e., considering the cloud particle as an opaque particle. It should be noted that ice particles which allow significant light transmission will have additional sources of error that are not captured in this experiment. A further assumption is that the diffraction pattern produced by an opaque particle is accurately described by the diffraction pattern produced by an opaque planar object representing the cross section of the particle. Laboratory studies showed that these approximations work well for out-of-focus transparent spherical particles e.g.,. For this study, we assume that it is also true for ice cloud particles.

The diffraction pattern of an opaque planar shape can be computed with different theoretical methods. computed the diffraction pattern of an opaque disc using the Maggi–Rubinowicz representation of the Helmholtz–Kirchhoff diffraction integral . Also, the Maggi–Rubinowicz method (MRM) can be adapted to other planar shapes. At the same time, the method is quite time-consuming and the analytical parameterization of the non-circular particle contour needs to be developed for each shape. In this study, we employed the method proposed in , which is based on angular spectrum theory (AST; see Appendix ). Despite the apparent differences, “the angular spectrum approach and the first Rayleigh–Sommerfeld solution yield identical predictions of diffracted fields” p. 61. Recall that the Helmholtz–Kirchhoff diffraction integral and the Rayleigh–Sommerfeld formula can be considered different formulations of the general scalar diffraction theory e.g.,. Note also that the validity region of Rayleigh–Sommerfeld solutions includes the Fresnel approximation and the Fraunhofer approximation regions of validity e.g.,p. 362.

Fast Fourier transform (FFT) can be utilized for the numerical implementation of the AST. Also, FFT-based algorithms are easy to implement and effective. That is why the AST is extensively employed in different domains, including simulations of diffraction patterns e.g.,. According to our computations, the AST-FFT modeling is several orders of magnitude faster than MRM-based calculations. AST-FFT modeling needs as input information a binary matrix representing the opaque shape of the particle and its location within the optical path. The latter is defined by the distance Z from the particle to the position at which we want to compute the diffraction pattern or the distance from the object plane to the particle when using an optical system. We performed numerical simulations of diffraction patterns recorded by a binary OAP probe, which in this study corresponds to a 2D-S probe described in Sect. .

As noted by , the diffraction pattern by an opaque disc can be presented as a function of only one dimensionless variable Zd=λZR2, with λ the wavelength and R the radius of the disc. Note that Zd=1NF is the inverse of the well-known Fresnel number NF. Figure  shows the diffraction pattern by an opaque disc at a specific distance Zd=1.0. Figure axes ξx and ξy are coordinates normalized by R.

We notice that a bright spot, called Poisson's spot, appears at the center of the diffraction pattern shadow image. The orange dashed line represents the 50 % intensity threshold generally applied in binary (monoscale) OAP probes. For this specific case, it can be seen that an OAP operating with this threshold will produce a “donut” image with an external diameter that exceeds by 30 % the true disc diameter.

Figure , calculated in analogy to Fig. 2 of , shows the light intensity along a radius of the diffraction pattern produced by an opaque disc with a normalized radial coordinate ξ=rR as a function of Zd.

With increasing distance Zd from the object plane (where Zd=0), the opaque disc is more and more out-of-focus and the diffraction pattern shows a Poisson spot which continuously increases and an external diameter which generally increases with Zd, albeit with small oscillations for Zd<2. From Zd=0.15 the disc shadow can appear with a magnification greater than 110 %. This arbitrary limit defines the separation of in-focus from out-of-focus images. As Zd continues to increase, the diffraction pattern becomes more and more blurry, and at some point the light intensity no longer falls below a specific triggering level at any point of the pattern. This means for a 50 % intensity triggering level of a binary OAP probe that the disc shadow totally disappears beyond Zdmax=8.17 . This distance delimits the region of the depth of field (DoF) of the particle . Beyond this limit, a particle is out-of-DoF and is no longer recorded on a 50 % intensity triggering level instrument. As Zmax=R2Zdmax/λ, the DoF increases with the square of the particle size and the inverse of the wavelength.

The two-dimensional stereo (2D-S) probe is an OAP that records the diffraction pattern of particles illuminated by a laser beam with wavelength λ=783 nm on a 128-photodiode array. As a particle crosses the laser beam of the instrument, slices of 128 pixels are recorded one after another at a specific frequency which is chosen such that the pixel size in the direction of the moving particle is 10 µm. The size of the pixels in the array direction can slightly vary from one instrument to another. Based on our own calibration using a spinning glass disc with imprinted opaque disc shapes (described further below), we found a mean value of 11.4 µm for the pixel size parallel to the array for the 2D-S used in this study. Note that the 2D-S has been sent to the SPEC Inc. company for a complete check prior to the tests for this study. Images recorded are monochromatic images based on a 50 % intensity triggering level. The transmitting optics consist of a single-mode fiber-coupled diode laser and beam shaping optics . Thus, we can assume that a particle is illuminated by a monochromatic plane wave. The receiving system consists of imaging optics and a linear photodiode array. The imaging optical system is based on a Keplerian telescope design. The photodiode array is positioned in the focal plane of the back lens (the eyepiece) in the image space. The object plane is the conjugate plane, that is, the focal plane of the front lens (the objective) in the object space. Based on the instrument optics, the object plane is located in the middle of the laser beam between the two arms of the probe (Fig. a).

Therefore, distance Z discussed in the previous section corresponds to the distance of the particle to the object plane (located at Z=0). The diffraction image of a particle crossing the laser beam at Z=x is identical to the diffraction image produced by a particle crossing the laser beam at Z=-x.

Considering opaque discs crossing the 2D-S laser beam, Fig.  shows the evolution of the DoF limit Zmax and the in-focus limit Z110% as a function of the variable opaque disc diameter D.

A particle with D=500 µm is seen in-focus, i.e., Dedge<110%D, as long as the distance of the particle to the object plane is Z<1.2 cm. From Z=1.2 to Z=3.1 cm (arm limit), the particle is out-of-focus; i.e., the particle is still detected (at least one photodiode triggered), but its image is progressively deformed due to diffraction. For a particle with D=50 µm, the in-focus zone is very small (Z110%=0.01 cm) and the particle is no longer detected beyond Zmax=0.65 cm. Particles larger than 109 µm should be always seen by the 2D-S since the DoF is starting to exceed the distance between the probe arms (6.2 cm). However, particles of that size are potentially observed with more or less important distortion in the out-of-focus domain. Particles larger than 806 µm should be imaged by the 2D-S without important distortion, since Z110%>3.1 cm for sizes larger than 806 µm.

In the following third section, we compare theoretical diffraction patterns of different opaque shapes with experimental measurements of the 2D-S probe. Therefore, several spinning glass discs with various chrome opaque particle shapes imprinted on the glass disc surfaces were used (Fig. b). This has been performed for opaque disc shapes in the past e.g., and most recently also for plate and rosette-type particles . Here we present results for two opaque planar particle shapes: rectangle-type shapes which represent the cross section of a columnar particle and “H”-type shapes which represent the cross section of a capped columnar particle. Particle shapes were imprinted with three different orientations: 0, 45, and 90∘. Note that the chosen shapes represent only one particular projection of a 3-D particle (columnar or capped columnar particles) on a 2-D plane.

Figure  shows results for a 1:2 (width : height) short columnar particle with four different sizes (the well-known results for opaque discs are reported in Appendix ).

In each column, on the left-hand side are shown the theoretical 2D-S records and on the right-hand side are shown the images recorded by the 2D-S experimentally. The theoretical 2D-S records are obtained by AST-FFT simulations (Sect. ) with a modeling resolution of 1 µm, 50 % intensity threshold, and subsequently 10×11.4 µm pixelization. Note that images are framed according to the particle size; i.e., smaller particle images are up-scaled with respect to larger ones. Each line shows results for particles at a specific distance Z from the object plane (Z=0 cm) to a distance very close to the arm of the probe (Z=3 cm). Results obtained for negative values of Z look very similar (not shown). Note that different particle orientations (0, 90∘; see Sect. ) were chosen for the measurements in order to avoid splitting of a particle into two or more images by the probe's image separator, at least in one orientation. A striking result is that measurements obtained with the 2D-S probe are in really good agreement with the theoretical diffraction simulation results, not only in terms of diffraction pattern, but also with respect to the DoF limit which is illustrated by the disappearing image at approximately the same Z value in theory and the corresponding measurement. We notice that several Poisson spots can appear in the same particle shadow, as observed for example for the 100×200 µm particle passing at Z=2 cm in both the simulation and the measurement. In contrast to opaque disc particles which preserve their contour shape, columns recorded by the 2D-S can end up in 2-D images that are very different from a column, as seen for the 75×150 µm particle at Z=1.5 cm. The image resembles more a disc with Poisson's spot than a column. Also, particle images can present small patterns detached from the main particle, as seen at Z=1.5 and Z=2.5 cm for the 100×200 µm particle. Approaching the DoF limit, the columnar particle is split into two symmetrical particles with the shape of a “crescent moon”. The evolution of the particle shape with Z with a 0.01 µm step is shown in the Supplement videos (available here: columnar particle with sizes of 25×50 µm (10.5446/40656), 50×100 µm (10.5446/40658), 75×150 µm (10.5446/40660), and 100×200 µm (10.5446/40662)).

Figure  then shows results for a 1:4 (width : height) elongated column with three different sizes.

Comparison between theory and measurements again shows very good agreement. Also for this crystal geometry, diffraction can produce patterns that are very different compared to the initial shape. In addition, we notice that these patterns look very different than those found in Fig.  for the shorter columns. Here, the diffraction pattern turns into a capped columnar shape as Z increases (see, e.g., 75×300 µm particle at Z=3 cm) and ends by splitting into two symmetrical particles when Z approaches the DoF limit (see, e.g., 50×200 µm particle with Z>2 cm). This is also illustrated in the Supplement videos (available here: columnar particle with sizes of 25×100 µm (10.5446/40657), 50×200 µm (10.5446/40659), and 75×300 µm (10.5446/40661)).

Figure  shows results for four capped columnar particles with different ratios and sizes.

Once again, the very good agreement between simulated and recorded images is striking. The diffraction pattern for one individual capped columnar particle can adopt many different image shapes, as can be seen for example for the 150×250 µm capped column. We notice that several small patterns may be detached from the main particle as a function of distance Z from the object plane. We also notice that capped columnar particles can fall apart into three distinct “large” parts in the produced binary diffraction image as seen for the 200×150 µm particle at Z=2 cm. The Supplement videos (available here: capped columnar particle with sizes of 75×125 µm with 25×25 µm bar (10.5446/40652), 100×75 µm with 50×25 µm bar (10.5446/40653), 150×250 µm with 50×50 µm bar (10.5446/40654), and 200×150 µm with 100×50 µm bar (10.5446/40655)) are very helpful for visualizing the evolution of the diffraction pattern.

As a particle moves away from the object plane, we notice that its image becomes more and more roundish regardless of its initial shape. The information of the real shape of the particle ends up being lost as the diffraction patterns progressively adopt a more circular form. This is particularly striking on videos. This means that any particle shape far from the object plane produces a more and more circular diffraction pattern which no longer allows us to identify the original shape of the respective particle. See for example the theoretical 75×125 µm capped columnar particle diffraction pattern at Z=1 cm (Fig. ) or the theoretical and measured 75×150 µm short column diffraction pattern at Z=1.5 cm (Fig. ). Depending on the initial shape, the binary diffraction pattern is generally broken into two, sometimes three image parts of similar sizes when approaching the DoF limit. Depending on the orientation of the particle, these two or three particles can be interpreted as different particles by the 2D-S probe. Indeed, the 2D-S probe stores a new particle as soon as a slice with 128 white pixels (image separator set by probe) is found. Because these two or three “particles” are registered as very close particles in time, algorithms will generally remove these particles from the data set, thereby considering them to be particles generated from shattering . This can create significant differences when particles are not properly recorded or deleted by the shattering algorithm before reaching the theoretical DoF limit. For example, the 75×150 µm short columnar particle (Fig. ; second column) might be removed by the shattering algorithm beyond Z=1.59 cm (see video in the Supplement available here: 10.5446/40660), whereas its theoretical DoF exceeds the distance between the arms. Therefore, it seriously affects the sample volume which is based on the DoF limit and consequently affects the retrieved concentrations.

Another interesting remark based on these results is that an out-of-focus image of a distinct particle shape can closely resemble another particle of a very different shape. As an example, note that the diffraction pattern of a 75×300 µm columnar particle at Z=2.5 cm (Fig. ) looks a lot like an in-focus capped columnar particle in the object plane. Also, the diffraction pattern of a 75×125 µm capped columnar particle at Z=1 cm (Fig. ) looks like a disc with Poisson's spot, meaning that an out-of-focus capped columnar ice particle can be interpreted as a droplet faintly out-of-focus.

In this section, we present some diffraction simulation results for four chosen particle shapes to illustrate the evolution of the particle diameter, including uncertainty evaluation. Our purpose here is neither to present an exhaustive list of results related to each shape nor to quantify the uncertainty of the probe in an absolute manner.

The size of the particle from a 2-D image has no absolute definition. Several definitions are used in the literature with different pros and cons depending on the objective of the study. In the study presented here, we illustrate results with two commonly used particle size definitions: the surface-equivalent diameter Deq and the maximum diameter Dmax. Deq is defined as the diameter of a disc with the same surface as the analyzed particle image. Dmax is defined as the diameter of the smallest circle encompassing the particle image e.g.,. Note that all triggered pixels of the image are considered here for the computation of Deq and Dmax, even when the pattern is split into several fragments (see Sect. ).

Figure  shows the theoretical evolution of Deq (solid blue line) and Dmax (solid green line), applying the 50 % threshold to a simulated 1×1 µm pixel image pattern, as a function of Z for four different particles: two rectangles with very similar true Deq (119.7 and 112.8 µm) but different aspect ratios (1:2 and 1:4, or short and elongate columns, respectively) and two capped column-type particles with the same shape but different sizes.

Dashed lines show the true Deq and Dmax particle diameters. The produced binary images from diffraction simulations are presented on top of each sub-figure of Fig.  for a few distinct distances of Z (Z=0, 0.5, 1.0, …, and 3.0 cm). Blue and green shadow areas in this figure show Deq and Dmax for the theoretical records by the 2D-S (pixel of 10×11.4 µm). Indeed, a spread for each size would appear due to the discrete pixel effect. As a photodiode needs to be shadowed from at least 50 % light intensity to be triggered, the number of triggered pixels will depend on the position of the particle shadow on the photodiode array of the probe. To account for this discrete pixel effect, particles are systematically shifted over one 10×11.4 µm 2D-S pixel by increments of 1 µm in both directions. This discrete pixel effect is illustrated in the Supplement video (available here: 10.5446/40663), showing the 150×250 µm capped columnar particle at Z=0.5 cm. We note that the pixels at the edge can be triggered or not as the particle is shifted, which then affects the particle size.

At first, we discuss diffraction simulation results of short and elongated columns presented in Figs. a and b. Both figures illustrate the evolution of Deq and Dmax as both particles pass the laser beam between Z=0 and Z=3.1 cm. The evolution of Deq is rather smooth with varying Z, compared to Dmax evolution, which is much more oscillating. Still, both columnar particles have a comparable true Deq (119.7 µm versus 112.8 µm), and both Dmax and Deq evolution show some nice similarities. Nevertheless, between Z=0 and Z=3.1 cm, for the short column (Fig. a) Deq and Dmax from theoretical binary images are always greater than the corresponding true Deq (119.7 µm) and Dmax (168.6 µm) values, whereas the Deq and Dmax retrieved from binary diffraction images of the elongated column (Fig. b) are lying on either side of the theoretical value lines, Deq=112.8 and Dmax=207.4 µm. Repeated abrupt decreases in measured Dmax are related to changes in the outer binary pixel ensembles of the diffraction pattern which are drifting away from the particle center as Z increases. This then leads to sudden loss of outer pixel and related decrease in Dmax. Detaching pixel ensembles are frequently observed, with two to three produced sub-images stemming from the diffraction pattern of one particle. This can lead to virtually large Dmax (e.g., Dmax=240 µm for the short column at Z=3.0 cm corresponding to the true Dmax of 168.6 µm).

Secondly, the two capped column-type particles are discussed. With increasing Z, the size of the larger capped column (Fig. c) increases for both diameter definitions Deq and Dmax. Deq of this particle when observed at a distance of 3 cm is about 210 µm, which exceeds by 50 µm the true Deq of 159.6 µm of the particle. Also, Dmax generally increases with increasing Z, however with few transient smaller diameter decreases at some distances depending on the evolution of the diffraction pattern details at the edges of the binary particle image. The smaller capped columnar particle (Fig. d) disappears before reaching Z=3.1 cm (arm limit). The apparent particle size first grows with Z, then shrinks continuously in terms of Deq (more abruptly in terms of Dmax), followed by another phase of slight increase in Dmax, before the particle then completely disappears in both diameter definitions at a distance Z of roughly 1.8 cm. Again, it can be noticed that in general Deq changes more gently with Z, as compared to Dmax. Close to the DoF limit, which is estimated from diffraction simulations for the particle in Fig. d as the position Z where the binary image disappears (roughly at Z=1.8 cm), Deq is underestimated, whereas Dmax continues to be overestimated as the particle image is formed by a few distant pixels. We note also that DoF limits from calculations of Zmax=D2Zdmax/(4λ) produces different DoF limit values for Deq and Dmax size definitions for the same particle . The DoF limit values are Zmax=1.66 (using Deq) and Zmax=4.13 cm (using Dmax). These different DoF limit estimations (factor of 2.5 between both calculations) using the above equation for this particle means different sample volumes, and finally different concentrations. For this particle, using the DoF limit estimations with Deq would be closer to the DoF limit (Z=1.8 cm) found from diffraction simulation. Moreover, as the uncertainty in Deq for out-of-focus particles is relatively small (Table 2) compared to Dmax (Fig. 8), it should be a relatively good option to estimate the DoF limit. However, our arbitrary diameter definition used in the classical DoF limit calculation remains questionable, since ice particles are primarily non-spherical.

Furthermore, for both columnar and capped columnar particles, it is evident that the discrete pixel effect (shadow areas) is almost negligible with respect to the diameter variability along the Z distance according to the diffraction simulations.

Finally, Table  summarizes Fig.  in terms of maximum, minimum, and average Deq and Dmax diameters over the whole distance between the two arms of the 2D-S probe compared to the true Deq and Dmax. For these four particles, we note that uncertainty for particles with a DoF limit beyond the arm limit spans from -23 % to +35 % in Deq and from -15 % to +48 % in Dmax. For the smallest particle, uncertainty spans from -85 % to +39 % in Deq and from -88 % to +69 % in Dmax. For this small particle, with a DoF limit smaller than the arm limit, the lower bound of the uncertainty has been calculated for a 1-pixel particle.

In this section, we compare particle size distributions retrieved theoretically from diffraction pattern simulations and experimentally measured by the 2D-S probe. For the measurements, a spinning disc (Fig. b) has been utilized with imprinted short columnar particles shown in Fig. . The spinning disc contains four different particle sizes, all of them imprinted in three different orientations (0, 45, and 90∘) and repeated six times. Therefore, each particle size should be seen 18 times at each revolution of the spinning disc. We simulate and measure 10 s of disc spinning at 108.6±0.2 rps (≈9.5 m s-1 in equivalent particle speed), which should result for each of the four short columns in about 19 500 images; 9.5 m s-1 already represents the maximum equivalent speed of particles on the rotating disc, which is small compared to aircraft speeds and therefore does not allow us to study possible effects of electronic response time related to disc speed.

Figure a shows theoretical results obtained from diffraction patterns by simulating each of the four particles 19 500 times.

Each of the three orientations (0, 45, and 90∘) accounts for one-third of the contribution of each of the four particles to the particle size distribution. The true equivalent diameter Deq of a 100×200 µm particle is 159.6 µm. At Z=0 cm, the particle image projected onto the photodiode array is undistorted. However, we see in Fig. a that the diameter is not retrieved perfectly and is also slightly varying. Depending on the position of the shadow particle projected onto the photodiode array, Deq can vary from 152 to 169 µm, which is due to taking into account the discrete pixel effect. Therefore, even in the best case when particles cross the laser beam of the 2D-S in the object plane and without any noise, this particle can be recorded with an uncertainty in Deq of up to 6 %, only due to the discrete pixel effect. By moving the 100×200 µm particle away from the object plane to Z=1 cm, then 2 cm, and finally 3 cm, we notice that Deq first starts to increase, and then decreases. The described behavior depends on particle size and shape, but has a common feature: particle Deq generally starts to increase and then decreases until it disappears when reaching the DoF limit, as shown for the small capped column in Fig. d.

Figure b shows the size distribution of the four short columnar particles measured by the 2D-S. At Z=0 cm, the number of counted particles (yellow curve with four distinct modes attributed to four particle sizes) for the four particle sizes (from larger to smaller ones) is 19 488, 19 484, 19 464, and 29 490. Note that most of the very small particles (<4 pixels at Z=0 cm) are due to dust on the spinning disc (the first four bars of the orange histogram). The number of counted particles smaller than 4 pixels at Z=0 cm is 12 029. At Z=0 cm, the four observed particle modes can be clearly attributed to the four different short columnar particle sizes. Distributions show reasonably narrow peaks. The position of the individual peaks is on average 5 µm below the expected respective true Deq, which is half of the pixel size. Experimentally, the discrete photodiodes in the array and the 50 % occultation criterion introduce a digitization uncertainty of roughly 1 size resolution (=10 µm for the 2D-S) depending upon where the particle passes across the array . We recall also that we found a mean value of 11.4 µm for the pixel size along the photodiode array based on our own calibration using a spinning glass disc with a printed opaque disc shape of 800 µm in diameter. Smaller size effects could be that different particles cross different photodiodes at different positions of the linear photodiode array, which is due to the fact that the rotation axis and the spinning glass disc center are not perfectly coaxial. Indeed, all photodiodes may not have an identical response. With increasing values of Z>0 cm, distributions of individual modes are getting broader. Whereas at Z=1 cm the three larger particles are observed, at Z=2 and Z=3 cm solely the two larger particles have their DoF beyond these distances. For the larger short columnar particle (100×200 µm), the four peaks corresponding to the four distances Z are well located compared to theoretical results, and each peak contains 19500±20 records. The second larger short columnar particle (75×150 µm) shows quite good agreement with simulations for peak position at Z=1 and Z=2 cm, and contains 19 514 and 18 548 records, respectively. The missing particles at Z=2 cm appear in the very small mode centered around 90 µm. In this small mode, we found 1924 particles which correspond to 962 particles split into two parts (Fig. ). Results at Z=3 cm show two distinct modes centered on 72 and 102 µm, well below the theoretical peak at 128 µm. The 102 µm mode is due to split particles with 45 and 90∘ orientations which were not separated by the probe into two images. This mode is not centered as the theoretical mode (128 µm) because the split particles show fewer triggered pixels in measurements compared to theoretical simulations (see Fig. ). This finding may be due to the fact that the probe's 50 % threshold does not perfectly coincide with the theoretical 50 % threshold. This effect may particularly impact the resulting probe images when approaching the DoF limit. The 72 µm mode is due to split particles with 0∘ orientation such that a separator can be set by the probe which generates two smaller particles. The sum of the record number in the larger mode corresponding to 45 and 90∘ orientations (13 090) and half of the smaller mode corresponding to 0∘ orientation (6403) yields in total 19 483. The 50×100 µm short column also shows two modes at Z=1 cm for the same reason. Particles smaller than 40 µm in Deq at Z≥ 1 cm are not dust particles on the spinning disc since their DoF limit is smaller than Z=1 cm. Actually, it can be seen from visual inspection of consecutive images that these small particles result from the separation of very small patterns from the main particles at the edge of the diffraction image. This effect is particularly striking at Z=3 cm and is consistent with findings published in the literature. For example, fragmented diffraction patterns of spherical droplets traversing the sample area near the edges of the DoF were shown in the work by ; diffraction fringes around out-of-focus images measured by CIP were underscored by . “Reacceptance” algorithms e.g., should address rigorously the problem of intact, i.e., not shattered, but fragmented particles.

We stated in Sect.  that the uncertainty in Deq due to diffraction in the out-of-focus region is far more important than the discrete pixel effect uncertainty. Table  shows the theoretical and measured uncertainty ΔDeq=ε±σ for the four opaque short columnar particles shown in Fig. , where ε represents the relative mean error and σ its relative standard deviation, both with respect to the true Deq. We notice that standard deviation σ, both in theory and measurements, is generally small compared to the mean error ε for out-of-focus particles, except when measurements present a two-modal distribution as seen for the 75×150 µm short columnar particle at Z=2 and Z=3 cm, and for the 50×100 µm short columnar particle at Z=1 cm. For these specific cases, we notice that particles with true Deq larger than 100 µm can have uncertainties up to 27%±1% theoretically and up to 28%±13% in measurements when split particles appear. Particles with true Deq smaller than 100 µm (the two smaller short columns of Fig. ) do not show uncertainties larger than derived for the two larger particles in this table because uncertainties are solely quantified for a few discrete Z distances. Uncertainties for particles smaller than 100 µm would dramatically increase as Z approaches the DoF limit as shown in Fig. d.