1. Introduction
The Small Cumulus Microphysics Study (SCMS), focused on the onset of precipitation in cumulus clouds, took place in July and August 1995 in the Cape Kennedy area (Florida). Three instrumented aircraft were involved, the University of Wyoming (UWYO) King-Air, the National Center for Atmospheric Research (NCAR)1 C130, and the Météo-France Centre d'Aviation Météorologique (CAM)2 Merlin-IV. An extensive microphysical dataset was collected during the field experiment with Commonwealth Scientific and Industrial Research Organisation (CSIRO) hot wires (King et al. 1978), a particle volume monitor, PVM-100A (Gerber et al. 1994), a cloud drop spectrometer (CDS) (Lawson and Cormack 1995), forward scattering spectrometer probes, FSSP-1003 and the Fast-FSSP (Brenguier et al. 1998), a modified version of the FSSP-100.
Droplet spectra measurements with FSSPs are affected by diverse sources of uncertainty that have already been identified and extensively discussed (Dye and Baumgardner 1984; Baumgardner et al. 1985; Cooper 1988; Brenguier 1989; Baumgardner and Spowart 1990). They are related to the nonuniformity of light intensity in the laser beam, to variations of the size calibration and of the sampling section of the instrument, and to the effects of coincidence in the sampling volume. The measured droplet spectra are not all similarly affected by these instrumental limitations. For example, variations of the sampling section are more significant for small droplets with a diameter close to the detection limit of the instrument; coincidence effects are more severe at high cloud droplet number concentration (CDNC). It is not straightforward in fact to evaluate the accuracy of the measurements when these effects are combined, when measuring a high concentration of small droplets for example. Numerical models of probe functioning have been developed for simulating these combined effects. Cooper (1988) followed a probabilistic approach with the transfer matrix of the instrument that indicates the probability for a droplet, ideally counted in one size class, to be counted in an another one (spectral broadening). The model was also extended to the coincidence of two particles (Cooper 1988). The inversion of the transfer matrix system allows retrieval of the actual droplet spectrum from the measured one. High-level coincidence of three or more droplets has also been simulated with a stochastic, or Monte Carlo, version of the model (Perrin et al. 1998), but there is presently no operational technique for inversion of that stochastic model. The difficulty to reproduce in the laboratory the same conditions as in flight is a serious obstacle to the experimental analysis of the probe response to droplets of different sizes, which is required for validation of the model of probe functioning.
The analysis of the extensive SCMS dataset presented here offers an alternative to the laboratory experiments. The droplet spectra that were collected are unknown, but one can expect that the different instruments converge in average to the actual response. In addition, since measurements were collected in cumulus clouds, one can use the adiabatic prediction as an absolute upper limit of the measured liquid water content (LWC).
The objectives of this study are threefold. First, the analysis aims at providing an evaluation of the data quality with quantification of the variability in the measurements of CDNC, mean volume diameter (MVD), and LWC for the whole SCMS campaign. Second, the diversity of the dataset is used to focus on specific samples, depending on the droplet sizes or the CDNC value, in order to examine separately the diverse sources of uncertainty. Finally, variations of some key parameters are documented for validation of the model of probe functioning.
The SCMS dataset is particularly suited to these objectives with numerous flights, a large diversity of microphysical conditions, and four FSSP probes: three standard versions with different electronic settings and size calibrations, and the Fast-FSSP. The finer size resolution of the Fast-FSSP, and its additional parameters, pulse duration and interarrival time, are particularly useful for the interpretation of the standard probes functioning. The list of flights and instruments considered in this study are indicated in Table 1. Comparisons are made between the Fast-FSSP, FSSP-100, CSIRO, and PVM-100A data, with statistics over the whole campaign, and for each flight separately. The microphysical data are then summarized for the 10 missions flown by the Merlin-IV during SCMS to illustrate the overall accuracy of the campaign.
2. Description of the FSSP-100 and Fast-FSSP
Very detailed descriptions of the FSSP-100 are already available in the literature (Dye and Baumgardner 1984; Baumgardner et al. 1985; Brenguier 1989). The Fast-FSSP is a modified version of the FSSP-100 with new electronics that measure for each detection, the pulse amplitude with 255 size classes instead of 15 in the standard probe, pulse duration and interarrival time from the previous detection with a resolution of 1/16 μs, and a parameter relative to the location where the detected particle crosses the beam (Brenguier et al. 1998). A few typical features of the FSSPs that are important for the following discussions are summarized here.
The FSSP sampling tube, parallel to the airflow, has an internal circular cross section of 40 mm in diameter. It is traversed, perpendicular to the airflow, by a laser beam focused to a diameter of about 0.2 mm at the tube axis. This incident beam is then masked by a dump spot in front of the detection optics. Droplets crossing the beam are detected by measuring the intensity of light scattered forward around the dump spot. Due to the optical setting of the instrument, droplets are, however, detectable only when they cross the incident beam within about 7.5 mm of the focal point. The corresponding cross section of the beam, perpendicular to the airflow, is referred to as the total sensitive area ST (15 × 0.2 mm2 = 3 mm2). The beam intensity is not uniform for all particle intersections of this area, and droplets of the same size can be counted in different size classes, hence producing spectrum broadening. It is therefore more accurate to select for sizing only those droplets that cross the beam through a narrower central section, hereafter referred to as the efficient sampling section (SEFF), where the incident light intensity is more uniform, and where the collected angles of scattering are similar. Beam intensity mapping (Baumgardner and Spowart 1990) shows a cylindrical symmetry, with a bell-shaped profile of light intensity across the beam. On both sides of the focal plane, the beam intensity also decreases along the beam axis because of the beam divergence. It follows that a particle receives the maximum incident light when it crosses the axis plane, perpendicular to the airflow. This occurrence also corresponds to the peak amplitude of the detected pulse that is recorded by the acquisition system. The objective for accurate sizing is thus to select an area centered at the focal point, limited along the beam axis to the nondivergent section of the beam (about 3 mm), and close to its axis (within 50%–60% of the beam diameter).
In the FSSP-100 this selection involves two steps. The first one is based on optical principle, using two detection diodes. The scattered light is split in two directions. One diode collects the signal for sizing. It is referred to as the signal diode. The second one, referred to as the annulus diode, is used for the selection. It is covered at its center by a small mask so that the image of a droplet crossing at the beam focal plane (1 and 3 in Fig. 1a) is projected on the mask and the diode receives less scattered light. On the contrary, the scattering image of particles crossing far from the focal plane (2 in Fig. 1a) is broader, so that scattered light passes around the mask. The comparison between signal and annulus amplitudes thus allows the rejection of particles crossing the beam too far from the focal plane. The selected area is referred to as the depth of field (DOF) cross section SDOF. The droplet counted rate used for CDNC calculation is the number of particles counted per second in DOF. However, droplet sizing requires a second step in order to reject some of the DOF accepted droplets, that is, those crossing the beam edge where the beam intensity is lower. This is performed by comparing the pulse duration to its running mean value. Droplets crossing the beam edge (3 in Fig. 1a), which is a shorter chord thought the beam, exhibit short pulse durations and they are rejected for sizing. This step is referred to as the pulse duration selection. The FSSP-100 acquisition system only records the droplet size information for detections that have been selected for DOF and pulse duration.
In the Fast-FSSP, both steps are accomplished with a slightly different optical principle. The second diode is entirely covered by a mask except for a narrow slit at the center that is oriented parallel to the airflow. A particle crossing out of the focal plane produces a broad image that is partially masked by the slit (2 in Fig. 1b). When the particle crosses at the focal plane, close to the beam edge, it produces a narrow image on the side of the slit (3 in Fig. 1b). Only droplets close to the beam axis and within a short distance of the focal point (1 in Fig. 1b) produce images centered on and smaller than the slit. This setting, originally developed for the FSSP-300 (Baumgardner et al. 1992), thus combines DOF and pulse duration selections. The corresponding droplet rate is used for CDNC and spectral shape calculations.
Consequently, the FSSP-100 and Fast-FSSP sampling sections for CDNC or spectral shape calculations are not identical. In the FSSP-100 the DOF selection ratio RDOF = SDOF/ST is about 20%. The spectral shape is derived from only DOF counts that are also selected for their pulse duration (a selection ratio of about 50%) that is overall REFF = SEFF/ST = 10% of the total counts.
In the Fast-FSSP, CDNC, and spectral shape are derived from the same set of selected counts. The efficient area depends upon the settings of the detection amplifiers. During SCMS it was set to SEFF = 0.13 mm2 (see Fig. 4 in Brenguier et al. 1998); that is, REFF = 4.5%. The Fast-FSSP acquisition system, in contrast to the standard system, records the four parameters: pulse amplitude, duration, interarrival time, and position flag for each detection individually.
3. Measurements of the droplet size distribution and probe limitations
In the FSSP-100 the droplet size information is only available for detections that have been selected for DOF and pulse duration. Because the corresponding sampling section is more difficult to measure than the DOF sampling section SDOF, it is common to proceed in two steps. First, the spectral shape is calculated based on those selected counts, assuming a constant value for Si. Second, the total droplet concentration is derived from the DOF selected counts λDOF, as N = λDOF/(SDOFυa). The droplet spectrum is finally obtained by scaling the previously derived histogram by this total droplet concentration (Σ Ci = N). Therefore, it is common with the FSSP-100 to distinguish between spectral shape (first step) and CDNC (second step) calculations. In particular, problems related to the pulse duration selection only affect the spectral shape calculation. With the Fast-FSSP both steps are combined using SEFF instead of SDOF.
In fact each of the parameters in (1) contributes to the uncertainty in the measurement of a droplet spectrum, and the main difficulty is to assess their respective contributions.
The precise definition of the size class boundaries, ϕi and ϕi+1, is referred to as the size calibration of the probe that relates the measured pulse amplitude to a scale of droplet diameters (Table 2). It is derived from Mie calculations of the forward scattered light integrated over the solid angle collected by the FSSP optics. The adequacy of the standard calibration function is discussed in Dye and Baumgardner (1984). In practice the probes shall be regularly calibrated using glass beads to account for possible variations of the optical alignment and dust contamination of the lenses. An interesting self-calibration technique has been developed for the Fast-FSSP. Because the relationship between droplet diameter and measured scattered light intensity is not monotonic, some values of intensity have a higher probability to be measured. They are precisely related to specific values of droplet diameter. The identification in the dataset of such frequently measured values provides for each flight the calibration coefficients of the instrument (Brenguier et al. 1998).
Despite the DOF and pulse duration selection procedures, the incident light intensity is not perfectly uniform throughout the efficient sampling section, and droplets of the same diameter can still be counted in different size classes depending on where they cross the beam. If the nominal or expected size is defined as the value that is measured when a droplet crosses the beam center, where the light intensity is maximum, such a missizing due to beam inhomogeneities results in spectrum broadening toward smaller sizes than expected (Baumgardner and Spowart 1990). In practice, the size calibration is performed by spreading glass beads randomly through the detection beam providing an intermediate diameter reference, with broadening occurring on both sides of the reference.
The probe sampling section Si is currently defined as a constant for data processing, though laboratory studies suggest that it could be slightly dependent upon the droplet diameter (Dye and Baumgardner 1984). This effect results in a spectrum distortion and errors in CDNC.
The counted rate is lower than the actual rate because of coincidence losses, when two, or more, droplets cross the beam simultaneously. The actual rate can be derived from the counted one using statistical formulas and either the counted rate, the probe activity (sum of pulse durations and electronic delays), a combination of both (Brenguier 1989), or the measured frequency distribution of interarrival times between detections (Baumgardner 1986). A pulse amplitude of coincident droplets is generally higher than the ones the droplets would produce separately. Therefore, a coincidence event reduces the counted rate in the size classes of the coincident droplets and increases the counted rate in a larger size class. Coincidences are thus a source of spectrum broadening toward larger sizes than expected (Perrin et al. 1998).
Airflow simulations have shown that the particle speed through the FSSP sampling tube is lower than the aircraft speed as it is measured with the onboard reference system (Norment 1988). Depending on how the airflow is slowed down in the vicinity of the probe and within its sampling tube, as well as the possibility of inertial separation of droplets of different sizes in the airflow, this effect can result in spectral distortion and CDNC errors. It will not be discussed here because it affects similarly the four FSSPs that are compared in this paper.
The ideal procedure for characterizing the probe's limitations would be in the laboratory to spread droplets of same sizes in various locations of the detection beam. However, the generation of small droplets and their trajectory control is a challenge, especially at a velocity comparable to the current aircraft velocity. For example, the droplet generator described by Wendisch et al. (1996) produces droplets bigger than 15 μm at a speed slightly larger than 1 m s–1. An alternative approach is considered here, which consists in analyzing a large dataset with very different microphysical characteristics in order to address the various issues separately.
a. Size calibration and beam inhomogeneities
The last Merlin-IV SCMS flight (me9514) was devoted to FSSP intercalibration. The three FSSP-100 and the Fast-FSSP were all mounted on the fuselage nose pods so that they were separated by less than 2 m, hence allowing close comparison of the measurements.
Figure 2 shows three examples of droplet spectra sampled in a cumulus cloud at various altitudes: 980, 660, and 1250 m for Figs. 2a,b,c, respectively (the cloud base was at about 450 m). The main difference between the three FSSP-100s is in the first size class, where the CAM shows a greater density than the NCAR and UWYO instruments. In fact only the CAM-FSSP was operated with the delay mode for pulse duration selection (Dye and Baumgardner 1984), an option that is known for increasing the counts in the first size class. The delay mode consists in measuring pulse duration at midamplitude instead of measuring it at a fixed threshold value. It is recommended by PMS to correct for a possible underestimation of the concentration density of the smallest droplets. However, in Fig. 2c the four instruments agree in showing a spectrum of large droplets, with a very low concentration density at diameters smaller than 10 μm. It is unlikely that the actual concentration density specifically increases in the first FSSP size class, as in the CAM-FSSP spectrum, while it is so low in the second class. Moreover, the CAM-FSSP spectra show a systematically high concentration density in the first size class, even in spectra of very large droplets at the upper limit of the diameter range, with a very low concentration density in the intermediate classes. It is thus concluded that the delay mode overestimates the number density in the first CAM-FSSP size class. The comparison with the Fast-FSSP even suggests that FSSP-100 measurements in the normal mode (NCAR- and UWYO-FSSPs) also overestimate the concentration density of the smallest droplets. Such an assessment though cannot be generalized to all FSSP probes because the response of this instrument to small droplets in the first class is very sensitive to the setting of the noise bias offsets in the detection module.
The droplet spectral width is a crucial parameter for the study of droplet condensational growth and the onset of precipitation (Brenguier and Chaumat 2001; Chaumat and Brenguier 2001). Specifically, a significant proportion of droplets much smaller than the mode is a feature that has often been reported in the literature (Hill and Choularton 1985; Bower and Choularton 1988). Many attempts have been made to explain the presence of small droplets by partial evaporation following entrainment-mixing processes, combined with secondary activation of entrained cloud condensation nuclei. However, it is still not yet clear how much of this feature is due to instrumental spectrum broadening. The three examples shown in Fig. 2 are among the narrowest droplet spectra sampled during the flight. Due to the size dispersion effect of the instrument the measured spectrum probably cannot be narrower than the actual one, and it can be surmised that the actual spectra are at least as narrow as measured with the Fast-FSSP. The difference in spectral broadening between the two instruments is attributed partly to the finer size resolution of the Fast-FSSP and partly to the larger efficient section of the FSSP-100 (0.3 mm2) compared to the Fast-FSSP one (0.13 mm2).
Because the detection threshold of the Fast-FSSP was set to 5.6 μm, a value close to the upper limit of the first FSSP-100 size class (Table 2), the FSSP-100 data are now processed without the first size class. Table 3 summarizes characteristic values of the three spectra shown in Fig. 2, with the mean and mean volume diameters, the spectral mode, and the spectral width, calculated as the standard deviation in diameter. The differences in diameters are of the order of 1 μm, but these differences can be as large as 2 μm, for example, between the Fast-FSSP and the CAM-FSSP mean diameters in the third example. For the spectral width, most of the differences are lower than 0.5 μm, with the exception of the Fast-FSSP versus NCAR-FSSP in the second example, and the Fast-FSSP versus the CAM and NCAR-FSSPs in the third example. It can also be noticed that the spectral width measured with FSSP-100 probes is usually larger than that of the Fast-FSSP, a feature that reflects the finer size resolution of the Fast-FSSP, with 255 size classes instead of 15 in the standard instrument, and the efficiency of the slit versus pulse duration selection procedure.
Most of the spectra sampled during SCMS were in fact broader than the examples presented in Fig. 2. In order to generalize the above analysis, MVD values of the 10-Hz sample measured by the four instruments are thus compared in Fig. 3. Figure 3a, for the Fast-FSSP versus CAM-FSSP, includes data from the 10 SCMS flights of the Merlin-IV from me9505 to me9514 (41 957 samples). Figure 3b, for the Fast-FSSP versus the NCAR-FSSP, includes data from two flights of the NCAR C-130, hc9504 and hc9505, and from the intercalibration Merlin-IV flight, me9514 (8454 samples). Finally Fig. 3c, for the Fast-FSSP versus the UWYO-FSSP, is limited to the Merlin-IV intercalibration flight (1614 samples).
The nonlinear size calibration scale of the Fast-FSSP is derived from modeling of Mie scattering and scattered light collection by the Fast-FSSP optics. The calibration coefficients were calculated for each flight, with the self-calibration technique described in Brenguier et al. (1998). The nonlinear calibration of the NCAR-FSSP follows the procedure described in Dye and Baumgardner (1984). The size scales of the UWYO- and CAM-FSSPs are linear (the class width Δϕi is constant), as recommended by PMS, with the standard scale for the UWYO instrument, while the size offset and class width of the CAM-FSSP were adjusted by intercalibration with the Fast-FSSP (Table 2).
Despite the fact that the CAM-FSSP was calibrated versus the Fast-FSSP, Fig. 3a reveals that there are still noticeable discrepancies between the two instruments. Fast-FSSP MVD values are underestimated with respect to the values measured by the CAM-FSSP in the range from 8 to 11 μm, and they are overestimated in the range from 12 to 18 μm. The two instruments only agree for MVD values larger than 20 μm. The same feature is observed in Fig. 3c for the comparison of the Fast-FSSP with the UWYO-FSSP measurements. In Fig. 3b, for the Fast-FSSP versus NCAR-FSSP the default of linearity is less pronounced, hence suggesting that the nonlinear scale recommended by Dye and Baumgardner (1984) is more suited than the linear one. These observations indicate that the choice of the calibration scale for processing FSSP data has a significant impact on measurements of spectral characteristics such as MVD, particularly in the small diameter range where the Mie calibration curve of the instrument is not monotonic. The above analysis reveals that the accuracy of FSSP MVD measurements in the range from 5 to 30 μm is of the order of ±1 μm in standard deviation, and that the standard linear calibration scale introduces an additional positive or negative bias of the same order, depending on the size range. For example, a 1-μm error in MVD at 15 μm leads to an error of 20% in the derived LWC.
b. Variations of the sampling section
In using (1) for CDNC calculations it is always assumed that the sampling section Si does not depend upon droplet size. However, measurements in the laboratory (Dye and Baumgardner 1984) revealed that the DOF increases with decreasing droplet sizes. This feature can be documented with the SCMS dataset by analyzing the ratio of selected to total counts that is equivalent to the DOF acceptance ratio RDOF = SDOF/ST, though the difficulty then is to determine if changes in RDOF with droplet size are due to SDOF or ST variations. Big droplets are producing a detectable pulse throughout the whole sensitive area that is thus constant and only limited by the optical probe setting. Small droplets in contrast are not detected over the whole sensitive area because of the decrease of the incident light intensity far from the focal point. The size detection limit of the instrument is not sharply defined. One should rather consider a detection range from ϕ1 to ϕ2; ϕ1 is the diameter of a droplet producing a detectable pulse only when it crosses the beam at the center, where the light intensity is maximum. Its sampling section, limited to a point, is null. Here, ϕ2 is the diameter of the smallest particle that produces a detectable pulse when it crosses the beam at the limit of the sensitive area. Between ϕ2 and ϕ1 the sensitive area decreases from ST to 0, thus affecting the acceptance ratio RDOF.
Figure 4 shows statistics of RDOF (REFF for the Fast-FSSP) versus MVD, with scatterplots of 10-Hz samples, as well as the mean value and standard deviation, calculated every 1-μm step. In Fig. 4a, statistics are based on samples from the series of C-130 flights from hc9506 to hc9520 for the NCAR-FSSP (32 361 samples), and the complete series of Merlin-IV flights (me9505 to me9514) for the CAM probe (51 119 samples). Figure 4b for the Fast-FSSP is based on samples from two C-130 flights (hc9504 and hc9505) and the complete series of Merlin-IV flights (49 699 samples). The two figures exhibit the same feature with an increase of RDOF and REFF with decreasing MVD values, which is in agreement with the laboratory measurements of Dye and Baumgardner (1984). Consequently, the CDNC values derived with a constant sampling section slightly overestimate the concentration of the small particles with respect to the large ones. Within the range from 10 to 20 μm, however, the Fast-FSSP REFF increases less than the FSSP-100 RDOF.
Figure 5 shows the ratio of the Fast-FSSP REFF to the CAM-FSSP RDOF values of Fig. 4, together with the ratio of their CDNC measured values, as functions of MVD. The CDNC underestimation of the Fast-FSSP with respect to the CAM-FSSP is noticeable in the range from 9 to 13 μm. At smaller droplet sizes the opposite is observed in agreement with the sharp increase of the Fast-FSSP REFF below 8 μm in Fig. 4b. This feature can be interpreted as a decrease of ST below 8 μm which, combined with the increase of SEFF, results in a sharp REFF increase at the lower limit of the probe's diameter range. This feature suggests that the threshold diameter ϕ2 for a constant sensitive area is close to 8 μm.
In fact it is not feasible to directly derive SDOF (SEFF for the Fast-FSSP) from data collected in flight, only RDOF (REFF) can be calculated. It is thus not feasible, from the comparison of Fast-FSSP versus CAM-FSSP, to determine firmly which one of the two probes is in error. This issue is addressed here by comparing both probe's estimations of LWC to the CSIRO measurements. Figure 6 shows the ratio of the FSSP-derived LWC, CAM-FSSP in Fig. 6a and Fast-FSSP in Fig. 6b, to the CSIRO measured one as a function of the FSSP-measured MVD, with mean value and standard deviation every 1-μm step in MVD. The data are the 10-Hz samples from the 10 Merlin-IV flights, with the condition that CDNC values are smaller than 300 cm–3 (35 105 samples) in order to avoid coincidence effects that are discussed in the next section. Assuming droplet sizing is accurate, the figure suggests that the CAM-FSSP overestimates the droplet concentration at MVD values smaller than 12 μm, while the Fast-FSSP slightly underestimates the concentration at MVD values between 9 and 12 μm, and slightly overestimates the concentration at MVD values smaller than 9 μm. It must be noted though that the reliability of the CSIRO probe at small droplet sizes is also questionable (Biter et al. 1987).
These biases in the measurements of small droplet concentrations with FSSPs should not be considered as an incorrect setting of the probe optics or electronics, as it rather seems to be inherent to the optical principle. This feature is in fact well reproduced with the model of probe operation of Perrin et al. (1998). The difference between FSSP-100 and Fast-FSSP arises from the different sampling section selection techniques: mask for the FSSP-100 and slit for the Fast-FSSP. The slight difference between the CAM and the NCAR-FSSPs in Fig. 4a should be attributed to the different gain settings of the signal and annulus amplifiers, a feature that is also well reproduced with the model of probe operation. In summary, the variation of the efficient sampling section in both the FSSP-100 and the Fast-FSSP is a source of uncertainty that mainly affects the measurements of the small droplets, especially at sizes below 8 μm, where the sensitive area is also dependent on droplet size.
The pulse duration selection in the FSSP-100 is a second source of uncertainty in droplet spectra measurements. The procedure is based on the selection of detections with a pulse duration longer than the running mean. The Fast-FSSP that records the pulse duration for each detection is well suited for the examination of this selection procedure. In Fig. 7 the cumulative frequency distribution of pulse duration measured with the Fast-FSSP is reported for three narrow spectra sampled during the Merlin-IV flights me9507 and me9511. The left-hand side shows the cumulative distribution for all the detections (thin line), and for only the slit accepted ones (thick line). The dashed line illustrates the expected cumulative distribution for particles crossing a uniform cylindrical beam of 230-μm diameter. In order to avoid variations due to the aircraft speed, pulse durations are normalized to an aircraft speed of 100 m s–1. The corresponding droplet spectrum for all the detections (thin line) and the slit-selected ones (thick line) is reported on the right-hand side. The values of mean droplet diameter in Figs. 7a,b,c are 8.4, 17, and 26 μm, respectively. The pulse duration distribution exhibits a tail toward long durations that is typical of the coincidence effects (Brenguier and Amodei 1989; Brenguier 1989). This issue will be discussed further. Note also the much sharper pulse duration distribution of the selected pulses that illustrates the efficiency of the slit technique for rejection of particles crossing the beam edge.
In the FSSP-100 the DOF selection procedure with the annulus only selects particles along the beam axis, that is, it includes particles crossing the beam edge. Their expected pulse duration statistics is thus represented by the dashed line in Fig. 7. Pulses with a duration longer than the mean correspond theoretically to particles crossing the beam within the 62% central section (Dye and Baumgardner 1984). When the droplet spectrum is narrow, this mean pulse duration selection is efficient and FSSP-100 spectra compare well with the Fast-FSSP ones, as in Fig. 2. When the droplet spectrum is broad, however, the procedure is no longer reliable because the pulse duration also depends on the particle size. This is illustrated in Fig. 8 where Fast-FSSP measurements of pulse duration from flight me9511 are represented as a function of the droplet diameter. Pulse durations longer than 3 μs that are due to coincidences are not accounted for in the calculation of the mean. The thin line corresponds to all the detections and the thick line to the slit selected ones. Simplified models of the expected variation of pulse duration with particle diameter are also plotted. The dashed line represents the pulse duration of a Gaussian pulse with a variable amplitude, when sampled at a fixed detection threshold: T = K{2 log[A(ϕ)/A0]}1/2, where K is a constant, A0 is the probe's detection threshold, and A(ϕ) is the pulse amplitude for a droplet of diameter ϕ. The dotted line represents the expected increase for particles entering a uniform beam: T = (D + ϕ)/υa, where D is the beam diameter and υa is the aircraft speed.
This dependency is emphasized in Fig. 7 with the complete distribution of pulse duration that is slightly shifted depending on the mean droplet diameter. When the spectrum is broad, the resulting statistics of pulse duration can be regarded as a combination of such distributions, with an intermediate mean value. Consequently, some small droplets crossing within the 62% central region can be rejected, while large droplets crossing the beam outside of this central area can have a pulse duration longer than the mean and as such can be erroneously selected. Since they cross the beam in a region of reduced intensity, they are counted in a smaller size class than expected thus producing a significant spectrum broadening.
The consequence of this is illustrated in Fig. 9 that shows, as in Fig. 2, the comparison of spectra measured with the Fast-FSSP and the three FSSP-100s during the intercalibration flight. However, in Fig. 9 broad bimodal spectra are displayed that exhibit modes at 8 and 15 μm. The four examples are sorted by increasing proportions of the larger mode. In Fig. 9a, the four probes show the same narrow spectrum with a tail toward large particles. As the proportion of big droplets increases, from Fig. 9b to 9d, it is noticeable that the three standard probes measure a significant percentage of the big droplets in the intermediate size class between the two modes, and thus miss the bimodality of the spectra. The figure also reveals that the CAM-FSSP does not perform better than the two other FSSPs in the identification of bimodality, thus suggesting that the delay mode is not efficient at compensating this broadening effect.
c. Effects of coincidences
The coincidence of particles in the detection beam is a limitation inherent to any single particle counter. First of all, it is important to notice that only detections made inside the efficient sampling section can be used for deriving a droplet spectrum (between 4% and 10% of the total), but that all the detections contribute to the probability of coincidence that is proportional to the sensitive volume of the instrument. The FSSP sensitive volume and efficient area are optimized for measurements of CDNC values up to 200 cm–3. At higher concentrations coincidences affect droplet counting and sizing.
Figure 10 shows the comparison of the 10-Hz CDNC-derived values after coincidence correction. In Fig. 10a, data are from the NCAR C-130 flight hc9505 (2542 samples). The Fast-FSSP is corrected using (2) and the NCAR-FSSP is corrected using (3). The NCAR-FSSP total counted rate to use in (3) (total reset) was not measured during flight hc9505. Therefore, it was initially derived from the DOF counted rate by assuming RDOF = 20%. However, the total reset was measured during the following flights and Fig. 4a shows that the actual RDOF value is larger than 20% and that it depends on the droplet diameter. NCAR-FSSP data were thus reprocessed according to the best linear fit in Fig. 4a; that is, RDOF = 37.35-0.47MVD. In Fig. 10b the data (1198 samples) are from the Merlin-IV flight me9514 and the Fast-FSSP is compared to the CAM-FSSP that is corrected with the technique of Brenguier (1989). Finally, Fig. 10c shows the same data as in Fig. 10b, with the CAM-FSSP data corrected using (3). As in the previous section, the first FSSP-100 class was not accounted for.
The Fast and the CAM-FSSP in Fig. 10b agree up to values as high as 800 cm–3, while the NCAR-FSSP and the CAM-FSSP corrected with (3) show a saturation at about 500 cm–3. Such a saturation of CDNC is not supported by the alternative technique based on interarrival time statistics, demonstrating that (3) underestimates the largest CDNC values. The values derived from the statistics of interarrival times in fact agree with the values derived with (2) within 5%, up to values larger than 1000 cm–3 (Burnet and Brenguier 1999).
Figure 11 shows the comparison of LWC values directly measured with the PVM-100A or the CSIRO probe, with the values calculated by integration over the droplet spectra measured with the Fast-FSSP, corrected using (2), during flight hc9505 (2390 samples) in Fig. 11a; with the NCAR-FSSP, corrected using (3), during the same flight (3134 samples) in Fig. 11b; and with the CAM-FSSP, corrected following Brenguier (1989), during the cloud traverses of flight me9508 with the highest CDNC values (422 samples) in Fig. 11c.
These comparisons of LWC are quite surprising and contradictory with the analysis of the CDNC measurements in the previous figure. The Fast- and the CAM-FSSP obviously overestimate LWC compared to either the PVM-100A in Fig. 11a for the Fast-FSSP or the CSIRO probe in Fig. 11c for the CAM-FSSP. This is attested by the calculated adiabatic LWC that provides a theoretical maximum value (for the data shown in Fig. 11c, e.g., sampled at about 600 m above cloud base, this upper limit is 1.4 g m–3, while CAM-FSSP LWC measurements exceed 2 g m–3). The NCAR-FSSP in Fig. 11b slightly underestimates LWC compared to the PVM-100A, but the discrepancy is small with regard to the sharp limitation of the CDNC measurements in Fig. 10a.
These features result from an additional effect of the coincidences; namely, the fact that a detection of coincident droplets is statistically counted in a larger size class than the ones of those same droplets counted individually (Perrin et al. 1998). The distortion is worse in the FSSP-100 because of the pulse duration selection procedure. The mean pulse duration is indeed lengthened by the longer pulses of coincident particles. This means that smaller droplets with short pulse durations can be biased out. The FSSP-100 droplet spectra are therefore shifted toward larger diameters than the Fast-FSSP ones at high coincidences rates. This effect is illustrated in Fig. 12 for four spectra measured at the same level, with increasing CDNC values. Table 4 summarizes characteristic values of the four spectra. The mean volume diameter measured with the Fast-FSSP is almost constant around 10.5 μm, with a spectrum width of about 2.15 μm. In contrast, the spectra measured with the CAM-FSSP exhibit MVD values increasing from 12.5 to 13.8 μm, with a spectral width increasing from 3.03 to 3.42 μm. Therefore, if only CDNC is corrected for coincidence losses, but the spectrum is not corrected for the shift toward larger diameters, the resulting LWC is overestimated. For the NCAR-FSSP corrected with (3), the underestimation of CDNC compensates the distortion of the spectrum.
One can also notice in Table 4 that CDNC values measured with the Fast-FSSP only increase from 757 to 856 cm–3 over these four samples, while the CAM-FSSP CDNC values increase from 995 to 1458 cm–3. This is in contrast with the general agreement between Fast-FSSP and CAM-FSSP CDNC values as illustrated in Fig. 10b. This additional discrepancy is due to variations of the Fast-FSSP efficient sampling section. Figure 13 shows Fast-FSSP REFF values versus MVD in Fig. 13a, and versus the corrected CDNC value in Figs. 13b,c. The data are 10-Hz samples from the same Merlin-IV flight (me9508). The MVD graph shows two groups of samples, below and above 12 μm. The two populations are plotted separately in Figs. 13b,c versus the CDNC value. Figure 13b reveals that REFF decreases with increasing CDNC values for samples with MVD smaller than 12 μm. The coincidences are thus reducing REFF when the droplets are small, hence leading to an underestimation of CDNC. It is likely that this feature is related to the probability for a coincidence of droplets to be selected, when some of the coincident droplets are crossing the beam outside of the efficient area, as discussed by Cooper (1988) for the standard FSSP.
The effects of the coincidence on the spectral shape and CDNC were demonstrated by Perrin et al. (1998) using the stochastic model of probe functioning that is, however, not usable for the retrieval of the actual spectrum from the measured one. Therefore, the FSSP (standard and Fast) spectra measured during SCMS are only corrected for the coincidence losses, but they are not corrected for spectral distortion. As mentioned above, the coincidence correction was applied similarly to all the size classes. Because of the spectral distortion and the variations of the efficient sampling section, the correction factor should in fact be size dependent. The stochastic model of probe functioning is presently used to better understand these features. They are presented here to show that under some circumstances FSSP measurements, either Fast or standard, can be significantly affected by the coincidences.
4. Summary of the Merlin-IV dataset
The microphysical dataset of the 10 Merlin-IV flights is summarized in Fig. 14, with the comparison of Fast-FSSP and CAM-FSSP measurements for CDNC and the CSIRO probe for LWC. The dataset consists of 10-Hz samples from 350 samples during the shortest flight (me9513) up to 7500 samples during the longest one (me9505). For each comparison, the ratio of the 10-Hz measured values is calculated, and its statistics is represented in Fig. 14 by the mean ratio (black dots) and the 20% and 80% percentiles of its cumulative frequency distribution (vertical bars). In addition, each flight is also characterized by the mean and the maximum values of the considered parameter, as measured by the Fast-FSSP. The maximum value is defined as the 98% percentile of the cumulative frequency distribution of the parameter.
The Fast- and the CAM-FSSPs show a very good agreement for the measurements of CDNC with a bias lower than 10% and a dispersion of ±20%, as expected with these instruments (Brenguier et al. 1994). The LWC-derived values show more discrepancy because LWC calculations accumulate errors on CDNC and errors on droplet sizes. In particular, it must be noted that the Fast-FSSP experienced a progressive loss of sensitivity during the campaign that was accounted for by the self-calibration technique (Brenguier et al. 1998). Such a procedure is not applicable to the FSSP-100 because it requires the very fine size resolution of the Fast-FSSP. Consequently, the size calibration of the FSSP-100 was kept constant over the whole campaign duration. Since the sensitivity loss of the Fast-FSSP was probably due to contamination of the optics by dust and sea salt, it is likely that the FSSP-100 was also affected and that it progressively underestimated the droplet sizes. This is reflected in the comparison by the positive LWC bias increasing with time (up to 30% on flight me9514).
This hypothesis is corroborated by the CSIRO probe. Its comparison with the Fast-FSSP-derived LWC shows no trend (Fig. 14d), while the comparison with the CAM-FSSP-derived LWC shows a progressive underestimation by the CAM-FSSP (Fig. 14e). Overall, the Fast-FSSP slightly overestimates LWC with respect to the CSIRO probe. On me9509 the ratio of Fast-FSSP to CSIRO LWC, however, increases up to 70%. This is due to the high CDNC values encountered during this flight that are responsible for a significant distortion of the measured spectra toward larger sizes than expected, hence resulting in the overestimation of LWC by the Fast-FSSP, while the CSIRO probe is not affected by the coincidences. Except for the coincidence effects, LWC measurements stay within an uncertainty range of ±30%, which is remarkable considering the differences between the two probe's operations.
5. Conclusions
The analysis of the large SCMS dataset, combining FSSP-100, Fast-FSSP, CSIRO, and PVM-100A measurements, provides high statistical significance in the characterization of the FSSP limitations. The conclusions are briefly summarized here.
Size calibration: The standard linear calibration of the FSSP-100 introduces a bias in the measurements of the small droplet sizes, overestimation of the diameter below 11 μm, and underestimation of the diameter in the range from 12 to 18 μm. Otherwise, the data from various FSSP-100 instruments agree within a standard deviation of ±1 μm in mean volume diameter, over the whole diameter range.
Beam inhomogeneities: The FSSP-100, with a sampling section for sizing of 0.3 mm2 (10% selection ratio), shows more spectrum broadening due to beam inhomogeneities than the Fast-FSSP, which has a smaller efficient sampling section of 0.13 mm2 (4.5% selection ratio). This results in a 0.5-μm overestimation of the spectral width by the FSSP-100 with respect to the Fast-FSSP. The sampling section can be electronically adjusted and reduced to an area of more uniform laser intensity, but it must be noted that a small sampling section also implies a reduced sampling rate, hence a poor statistical significance of the measured samples.
Variations of the sampling section: The DOF sampling section of the FSSP probes, both standard and Fast, increases with decreasing droplet size. It is not feasible to precisely document the dependence of SDOF with the droplet diameter from data collected in flight because the measured spectra are not monodispersed (MVD is used here as an indication of the mean droplet size) and because SDOF is not directly measurable; only RDOF = SDOF/ST can be recorded. The analysis of the Fast-FSSP dataset suggests that ST is constant down to about 8 μm in droplet diameter and that it decreases substantially for smaller droplets. In contrast SDOF increases continuously with decreasing droplet sizes. These variations are probe dependent and can be significant (e.g., the increase of the FSSP-CAM SDOF from 31 to 8 μm exceeds a factor of 2). Consequently, the concentration density of the small droplets is overestimated with respect to the density of the big ones.
Pulse duration selection: The procedure is efficient at rejecting particles crossing the beam edge when the spectrum is narrow and CDNC is low. The dependence between pulse duration and droplet size, however, introduces a bias when the spectrum is broad. The largest droplets are preferentially selected with respect to the smallest ones, but they are counted in a smaller size class than ideal. This can prevent the detection of bimodal spectra. This effect is accentuated at high CDNC values by the coincidence of particles that lengthen the mean pulse duration. The pulse duration selection procedure is not used in the Fast-FSSP, which is thus not affected by this limitation.
Coincidences: Coincidences of particles in the beam lead to an underestimation of the droplet concentration. That can be accurately compensated by using the statistical correction procedure of Brenguier (1989). Coincidences also produce a distortion of the droplet spectrum toward larger sizes than ideal, which introduces a significant overestimation of the derived LWC. Finally, coincidences can be responsible for variations of the efficient sampling section, depending on the droplet size. It has been shown, for example, that under specific conditions (high CDNC of small droplets), Fast-FSSP measurements of CDNC can be noticeably underestimated.
The size calibration bias, the broadening by beam inhomogeneities, and the coincidence losses can be corrected using straightforward procedures, such as the transfer matrix proposed by Cooper (1988), which linearly relates the measured spectrum to the actual one. However, spectra distortion due to variations of the sampling section, especially the one related to the pulse duration selection and to the coincidence of particles, cannot be linearized. In fact, the transfer matrix depends upon the actual droplet spectrum. Therefore, no operational procedure exists for correction of these effects. They can be precisely simulated with the stochastic model of probe functioning of Perrin et al. (1998), but, as in any stochastic (nondeterministic) model, it is more difficult to invert than the Cooper (1988) matrix model for the retrieval of actual spectra from the measured ones. Alternative approaches are presently tested.
Except for the above-mentioned limitations, the SCMS dataset collected with the Merlin-IV instrumented aircraft is consistent in terms of droplet concentration and liquid water content. A complete catalogue of the analysis is available upon request to the lead author.
Acknowledgments
The authors are grateful to the SCMS teams of the NCAR C-130, especially to Dr. C. Knight and Chris Webster for their support during the field campaign and the data processing, of the UWYO King-Air and of the Météo-France Merlin-IV for their remarkable performance. They acknowledge the crucial contribution of the Météo-France TRAMM team for data processing. This work was supported by Météo-France, CNRS, and the European Union under Grants ENV4-CT-0117 and EVK2-CT-1999-00054.
REFERENCES
Baumgardner, D., 1986: A new technique for the study of cloud microstructure. J. Atmos. Oceanic Technol., 3 , 340–343.
Baumgardner, D., and Spowart M. , 1990: Evaluation of the Forward Scattering Spectrometer Probe. Part III: Time response and laser inhomogeneity limitations. J. Atmos. Oceanic Technol., 7 , 666–672.
Baumgardner, D., Strapp W. J. , and Dye J. E. , 1985: Evaluation of the Forward Scattering Spectrometer Probe. Part II: Corrections for coincidence and dead-time losses. J. Atmos. Oceanic Technol., 2 , 626–632.
Baumgardner, D., Dye J. E. , Gandrud B. W. , and Knollenberg R. G. , 1992: Interpretation of measurements made by the Forward Scattering Spectrometer Probe (FSSP-300) during the airborne Arctic stratospheric expedition. J. Geophys. Res., 97 , 8035–8046.
Biter, C. J., Dye J. E. , Huffman D. , and King W. D. , 1987: The drop-size response of the CSIRO liquid water probe. J. Atmos. Oceanic Technol., 4 , 359–367.
Bower, K. N., and Choularton T. W. , 1988: The effects of entrainment on the growth of droplets in continental cumulus clouds. Quart. J. Roy. Meteor. Soc., 114 , 1411–1434.
Brenguier, J-L., 1989: Coincidence and dead-time corrections for particle counters. Part II: High concentration measurements with an FSSP. J. Atmos. Oceanic Technol., 6 , 585–598.
Brenguier, J-L., . 1993: Observations of cloud microstructure at the centimeter scale. J. Appl. Meteor., 32 , 783–793.
Brenguier, J-L., and Amodei L. , 1989: Coincidence and dead-time corrections for particle counters. Part I: A general mathematical formalism. J. Atmos. Oceanic Technol., 6 , 575–584.
Brenguier, J-L., and Chaumat L. , 2001: Droplet spectra broadening in cumulus clouds. Part I: Broadening in adiabatic cores. J. Atmos. Sci., 58 , 628–641.
Brenguier, J-L., Baumgardner D. , and Baker B. , 1994: A review and discussion of processing algorithms for FSSP concentration measurements. J. Atmos. Oceanic Technol., 11 , 1409–1414.
Brenguier, J-L., Bourrianne T. , Coelho A. , Isbert J. , Peytavi R. , Trevarin D. , and Weschler P. , 1998: Improvements of droplet size distribution measurements with the Fast-FSSP (Forward Scattering Spectrometer Probe). J. Atmos. Oceanic Technol., 15 , 1077–1090.
Burnet, F., and Brenguier J-L. , 1999: Validation of droplet spectra and liquid water content measurements. Phys. Chem. Earth, 24B , 249–254.
Chaumat, L., and Brenguier J-L. , 2001: Droplet spectra broadening in cumulus clouds. Part II: Microscale droplet concentration heterogeneities. J. Atmos. Sci., 58 , 642–654.
Cooper, W. A., 1988: Effects of coincidence on measurements with a Forward Scattering Spectrometer Probe. J. Atmos. Oceanic Technol., 5 , 823–832.
Dye, J. E., and Baumgardner D. , 1984: Evaluation of the Forward Scattering Spectrometer Probe. Part I: Electronic and optical studies. J. Atmos. Oceanic Technol., 1 , 329–344.
Gerber, H., Arends B. G. , and Ackerman A. S. , 1994: New microphysics sensor for aircraft use. Atmos. Res., 31 , 235–252.
Hill, T. A., and Choularton T. W. , 1985: An airborne study of the microphysical structure of cumulus clouds. Quart. J. Roy. Meteor. Soc., 111 , 517–544.
King, W. D., Parkin D. A. , and Handsworth R. J. , 1978: A hot-wire liquid water device having fully calculable response characteristics. J. Appl. Meteor., 17 , 1809–1813.
Lawson, R. P., and Cormack R. H. , 1995: Theoretical design and preliminary tests of two new particle spectrometers for cloud microphysics research. Atmos. Res., 35 , 315–348.
Norment, H. G., 1988: Three-dimensional trajectory analysis of two drop sizing instruments: PMS OAP and PMS FSSP. J. Atmos. Oceanic Technol., 5 , 743–756.
Perrin, T., Brenguier J-L. , and Bourrianne T. , 1998: Modeling coincidence effects in the Fast-FSSP with a Monte-Carlo model. Preprints, Conf. on Cloud Physics, Everett, WA, Amer. Meteor. Soc., 112–115.
Wendisch, M., Keil A. , and Korolev A. V. , 1996: FSSP characterization with monodisperse water droplets. J. Atmos. Oceanic Technol., 13 , 1152–1165.
Summary of the analyzed dataset; “hc” refers to the NCAR C-130 and “me” refers to the CAM Merlin-IV
Calibration scales of the FSSP-100 and the Fast-FSSP. The diameter value (μm) in each line represents the lower boundary of the indicated size class. For the Fast-FSSP, only 15 classes are indicated, with the corresponding class number, for the first flight of the campaign on the NCAR C-130 (hc9504) and for the last one on the Merlin-IV (me9514). The last line contains the upper boundary of the probe's diameter range
Spectral characteristics of the three examples shown in Fig. 2. The spectral width is derived as the spectrum std dev
Characteristics of the four droplet spectra shown in Fig. 12