Microphysical and optical properties of precipitating drizzle and ice particles obtained from alternated lidar and in situ measurements

During the international ASTAR experiment (Arctic Study of Aerosols, Clouds and Radiation) carried out from Longyearbyen (Spitsbergen) from 10 May to 11 June 2004, the AWI (Alfred Wegener Institute) Polar 2 aircraft was equipped with a unique combination of remote and in situ instruments. The airborne AMALi lidar provided downward backscatter and Depolarisation ratio profiles at 532 nm wavelength. The in situ instrumental setup comprised a Polar Nephelometer, a Cloud Particle Imager (CPI) as well as a Nevzorov and standard PMS probes to measure cloud particle properties in terms of scattering characteristics, particle morphology and size, and in-cloud partitioning of ice/water content. The objective of the paper is to present the results of a case study related to observations with ice crystals precipitating down to supercooled boundary-layer stratocumulus. The flight pattern was predefined in a way that firstly the AMALi lidar probed the cloud tops to guide the in situ measurements into a particular cloud formation. Three kinds of clouds with different microphysical and optical properties have therefore been quasi-simultaneously observed: (i) water droplets stratiform-layer, (ii) drizzle-drops fallstreak and (iii) precipitating ice-crystals from a cirrus cloud above. The signatures of these clouds are clearly evidenced from the in situ measurements and from the lidar profiles in term of backscatter and Depolarisation ratio. Accordingly, typical lidar ratios, i.e., extinction-to-backscatter ratios, are derived from the measured scattering phase function combined with subsequent particle shapes and size distributions. The backscatter profiles can therefore be retrieved under favourable conditions of low optical density. From these profiles extinction values in different cloud types can be obtained and compared with the direct in situ measurements.


Introduction
The impact of clouds on the radiation balance of our planet has been identified as a major unsolved problem in climate research (see among others Norris, 2005).It is generally recognized that inadequate parameterisation of the ice-phase, mixed-phase, and polluted clouds is one of the greatest sources of uncertainty in the modeling of climate processes (Sun and Shine, 1995;Gregory and Morris, 1996).Accurate determination of the microphysical and optical cloud properties, including the identification of particle phase and subsequent water partition (liquid and/or ice), is crucial for better understanding of the cloud formation and evolution particularly in Arctic regions (Curry et al., 1996) as well as the effects of anthropogenic emissions.Detailed observations of cloud microphysical and optical characteristics are also required for the validation of retrieved cloud parameters obtained from remote sensing devices (Korolev and Isaac, 1999).Indeed, remote sensing of cloud properties has become an increasingly determinant field, because retrieval algorithms used from ground (see among others Shupe et al., 2005) and from airborne platforms (Brogniez et al., 2004;Tinel et al., 2005) are now being adapted and implemented for widespread use from space like the A-Train project which includes the CALIPSO satellite (Stephens et al., 2002).The strategy of validation of retrieved products yet has to be developed for mixed-phase and ice clouds.In the case of water clouds, the spherical shape of the droplets considerably simplifies the validation of the retrieved effective diameter with the exception for highly contaminated particles by atmospheric pollutants.In the latter case the non-absorption hypothesis is no longer valid.Because of the occurrence of both spherical droplets and ice crystals in mixed-phase clouds and the complexity of the ice crystal shapes in such clouds as well as in glaciated (cirrus) clouds, the determination of an "effective dimension" needs to be reviewed with detailed in situ information on the shape of these particles and on the related effects on scattering properties.
J.-F.Gayet et al.: Microphysical and optical properties of precipitating drizzle and ice particles In this paper we present the results of a case study (5 June 2004) which has been documented during the ASTAR experiment (Herber et al., 2004) from observations carried out onboard the Polar 2 aircraft operated by the Alfred Wegener Institute (AWI).These observations relate alternated lidar and in situ measurements of ice crystals precipitating down to supercooled boundary-layer stratocumulus.The flight pattern consisted firstly to probe the cloud tops by lidar then to proceed to in situ measurements at lower levels in the cloud.The instruments used for the remote and in situ observations are presented first, together with the evaluation of the measurement errors.

Instrumentation
The Airborne Mobile Aerosol Lidar (AMALi) has been used onboard the Polar-2 aircraft to probe the backscatter properties of the atmosphere below the aircraft.This instrument, the performances and the retrieving technique for the endproducts determination have been thoroughly described by Stachlewska (2006a).
The instruments used for the determination of microphysical and optical properties of Arctic clouds during ASTAR include four independent techniques: (1) the Polar Nephelometer, (2) the PMS 2D-C, (3) the Cloud Particle Imager (CPI) and ( 4) the Nevzorov probe.The combination of these techniques provides a description of particles within a diameter range varying from a few micrometers (typically 3 µm) to several millimeters.The method of data processing, the reliability of the Polar Nephelometer and PMS 2D-C instruments and the uncertainties of the derived microphysical and optical parameters have been described in detail by Gayet et al. (2004).

Remote measurements
The Airborne Mobile Aerosol Lidar (AMALi) was developed by the Lidar Group of the Alfred Wegener Institute (AWI) in Potsdam (Stachlewska et al., 2004).The system is designed to provide remote sensing measurements at 1064 nm and 532 nm wavelengths with Depolarisation information at 532 nm.Onboard the Polar-2 aircraft, the AMAli was operated in the nearly nadir-looking configuration.The data acquisition system recorded single-shot lidar signals providing high temporal and hence horizontal resolution.Measured lidar profiles deliver qualitative information about properties of the atmosphere between the flight altitude and the ground or sea surface.To assure eye-safety, the laser beam divergence was large (2.6 mrad).The overlap range of the field of view of the telescope (3.1 mrad) and the laser beam determines a distance between the aircraft and the first reliable point in lidar profiles, which was at 235 m.The effects of the multiple scattering are neglected.The resolutions of lidar profiles used in this paper were chosen 7.5 m for the vertical resolution and 15 s for the temporal resolution.The latter one corresponds to a horizontal resolution of 1.2 km for the aircrafts cruising speed of about 80 m/s.
The method of data processing of the AMALi measurements to retrieve the backscatter and extinction coefficients and the Depolarisation ratio are discussed with details in the Appendix along with the evaluation of the errors.

Microphysical and optical measurements
The Polar Nephelometer (Gayet et al., 1997) measures the scattering phase function of an ensemble of cloud particles (i.e., water droplets or ice crystals or a mixture of these particles ranging from a few micrometers to about 1 mm diameter), which intersect a collimated laser beam near the focal point of a parabolic mirror.The light source is a high-power (1.0 W) multimode laser diode operating at λ=804 nm.The data acquisition system of the airborne version of the Polar Nephelometer is designed to provide a continuous sampling volume by integrating the measured signals of each of the detectors over a selected period.For instance, the sampling volume (v) is determined by the sampling surface (10-mm long and 5-mm diameter beam) multiplied by the Polar-2 cruise speed of approximately 80 m s −1 , i.e.: 500 cm 3 for an acquisition frequency of 10 Hz.This means that the detection threshold is 0.5 particle per litre at this frequency.Direct measurement of the scattering phase function allows particle types (water droplets or ice crystals) to be distinguished and calculations of the optical parameters to be performed (extinction coefficient and asymmetry parameter, see Gayet et al., 2002).The accuracies on extinction coefficient and asymmetry parameters derived from the Polar Nephelometer are estimated to be within 25% and 4%, respectively (Gayet et al., 2004).
The PMS 2D-C instrument provides information on crystal size and shape for the size range 25-800 µm.Because large ice crystals as larger as several millimeters have been observed during the ASTAR experiment, the 'reconstructed' method of 2-D data processing (Heymsfield and Parrish, 1978) has been used in this study.The 'reconstructed' method allows an extrapolation of the data to particles larger than 800 µm.The sampling surfaces have been derived accordingly for the calculation of microphyical parameters.In order to improve the statistical significance of low particle concentrations, a 5-s running mean was applied.Irregular ice particles were the most predominant crystals sampled in arctic clouds during ASTAR.Therefore, the bulk parameters were calculated assuming the surface-equivalent diameter relationships given by Heymsfield (1972) and Locatelli and Hobbs (1974).As the sensitivity of the probe to small particles decreases with airspeed (i.e.∼70 to 90 m s −1 with the Polar-2 aircraft), the five-first channels (up to 125 µm) were corrected according to the results of Baumgardner and Korolev (1997) assuming that these corrections take into account the miss and/or the under-sizing of the particles evidenced by Strapp et al. (2001) and by Lawson et al. (2006).The use of the first six 2D-C channels was justified as they were found reliable in terms of statistical significance by applying the Principal Component Analysis procedure (Shcherbakov et al., 2005).The accuracy of derived microphysical parameters is considerably affected by inherent shortcomings of probes and data processing.Considering the sampling time of 5 s, the total random uncertainties were estimated as ranging from 75% to 100% for PMS probe data (see details in Gayet et al., 2002).These errors may be considerably reduced by taking averages over longer periods.

Other measurements
In this study the data from the Cloud Particle Imager (CPI, Lawson et al., 1998) will only be used to give indication on the type and/or the shape of particles, which have been observed.Of course, microphysical parameters can also be derived from the CPI data (Lawson et al., 1998).As for the Nevzorov instrument (Korolev et al., 1998), the liquid water content (LWC) and the ice water content (IWC) can be deduced according to calibration relationships.The accuracies on these quantities have been evaluated to 30%.

Cloud situation and flight procedure
The observations discussed in this paper were obtained on 5 June 2004 between 07:40-11:35 UT over the Storfjorden area which is located in the South-East of the Svalbard Archipelago (see Fig. 1).From the AVHRR/NOAA15 satellite images (09:16 UT), the weather situation over the observation area was characterized by a multi-layers cloudy system, i.e., scattered cirrus clouds at upper levels with a low optical depth (rather transparent in the visible channel) and a broken stratiform layer at lower levels.The vertical sounding derived from the airborne measurements revealed that the boundary layer cloud extended from about 1200 m/−9 • C to1400 m/−11 • C with a wind rather homogeneous both in the magnitude (6±1 m/s) and the direction (5-30 • ) in the sampled altitude range (up to 3000 m: the aircraft ceiling).
The flight consisted, first, to perform nadir remote measurements with the AMALi lidar at the aircraft ceiling along a flight distance of 200 km while heading strait into southern part of Storfjorden (see Fig. 1).Afterwards, the aircraft descended to 1400 m/−11 • C and in-situ measurements were carried out at this level in order to document the cloud microphysical and optical properties along a horizontal distance of about 200 km and with a heading North oriented (see Fig. 1).
The 1400 m flight level has been chosen according to onboard real-time lidar observations which revealed interesting cloud features, i.e. stratiform layer and ice crystals precipitating (Stachlewska et al., 2006b).
The lidar measurements will be first presented to overview the backscatter cloud properties, then the cloud microphysi- cal and optical properties obtained at a single flight level will be discussed in order to give some insights about the validation of retrievals results in terms of particle phase (water droplets and/or ice crystals) and extinction coefficient.

Lidar observations
Figures 2a and b display respectively the vertical profiles (between 3000 m and the sea surface) of the lidar signal (background and range corrected) and the Depolarisation ratio at 532 nm during a flight sequence of about 200 km long, i.e. between 08:57 and 09:42 UT.The line superimposed on each Figs. 2 indicates the flight level of the in situ measurements (see Sect. 4) performed about one hour after the sequence of remote sensing observations.The analysis of the results on Fig. 2 highlights a broken cloud layer vertically extended from about 1200 m and 1600 m with strong backscattering signals and low Depolarisation ratio (1-5%).This feature reveals the presence of a stratiform cloud with water (spherical) supercooled droplets.In these cloud regions, the laser beam is rapidly attenuated due to very high concentration of scatters.Additionally, multiple scattering must be expected under these conditions.Hence, the lidar signal and Depolarisation ratio below these clouds and down to the sea level cannot be exploited (i.e.blank areas on Fig. 2b).An interesting feature, evidenced between 09:15 and 09:20, is the fallstreak signature attributed to precipitating freezing-drizzle because of a Depolarisation ratio as low as a few percent, typical for spherical particles.Subsequently, the cloud layer  was optically thinner compared to the other cloud parts, thus the lidar information can be processed down to the sea level.A similar situation was described by Gayet et al. (1994) where lidar data were inverted in drizzle precipitation situation which occurred in a boundary-layer stratiform cloud.At uppermost levels (i.e. from about 1600 m to 3000 m) and from the beginning of the flight sequence to 09:20, the high values of the Depolarisation ratio (25% to 50%) give indication that different types of particles were detected that is attributed to irregularly shaped ice crystals falling from the cirrus cloud above observed from satellite (AVHRR images).
The final part of the measurement between 09:36 and 09:42 UT was taken in a cloud-free background atmosphere.The noisy character of the Depolarisation ratio profiles observed in the clear atmosphere areas is due to lower target concentrations compared to aerosol-rich atmosphere or in clouds cases, resulting in a weaker backscatter signal at perpendicular 532 nm wavelength.The representative mean value of Depolarisation ratio of 5% in the clear parts of atmosphere above the water cloud and in cloud-free atmosphere between 09:25 and 09:42 UT was found.C level.C25: particle concentration of particles with D>25 µm, Deff: effective diameter, Ext: extinction coefficient, g: asymmetry parameter and LWC : liquid water content.The shape classification of particles larger than 50 µm from CPI images is also reported with three categories, i.e. water spherical droplets, regular shape (columns, plates, . . . ) and irregular shape particles.

In situ measurements
that may be related to the cloud properties inferred from the lidar observations depicted in Figs. 2 at the corresponding level.As a matter of fact, the examination of Fig. 3 shows that the aircraft samples first (10:27-10:31) a cloud which exhibits typical (supercooled) water droplets signature since the asymmetry parameter value is 0.845.The liquid water content (LWC from the Nevzorov instrument) and the concentration of drops larger than 25 µm (C25) reach 0.20 g/m 3 and 1000 l −1 , respectively.It should be noticed that very similar cloud properties are observed during the last part of the flight sequence (10:40-10:54).These cloud characteristics may be related to the broken stratiform water droplet layer detected by the AMALi in which the laser beam is rapidly attenuated due to very high concentration of particles.Indeed the cloud microphysical parameters retrieved from the measured scattering phase function (Gayet et al., 2002) are: 130 cm −3 , 0.25 g/m 3 and 17.3 µm for the droplet concentration, liquid water content and effective diameter respectively.We note in passing the good correspondence between the two liquid water content measurements from the Nevzorov probe and the Polar Nephelometer, i.e., 0.20 g/m 3 and 0.25 g/m 3 , respectively.Furthermore, at the fringe of the first sampled cloud patch (i.e. between 10:30 and 10:31:30), the CPI revealed the occurrence of drizzle droplets with diameter up to 500 µm (see examples of CPI images in Fig. 5).Therefore, the www.ann-geophys.net/25/1487/2007/Ann.Geophys., 25, 1487-1497, 2007 Coming back to Fig. 3, the results show that the middlepart of the sequence (10:35-10:42) exhibits quite different microphysical and optical cloud properties.Indeed, an optically thin layer is evidenced with low values of the asymmetry parameter, extinction and concentration of particles with large effective diameter.These characteristics reveal the occurrence of large ice crystals as confirmed by the CPI analysis (see examples in Fig. 5) and could also be related to ice particles precipitating from the levels above as inferred from the lidar data.Figure 4b, with a similar representation as in Fig. 4a, displays the results observed in the ice crystals region between 10:35 and 10:42 (see Fig. 3) which indicates the properties of the optically thin ice layer with mean values of the ice water content, extinction and effective diameter of 10 l −1 , 9.0 mg m −3 , 0.2 km −1 and 270 µm, respectively.The scattering phase function exhibits a smoothed feature typical of irregular ice-crystal signature (see Gayet et al., 1998) with an asymmetry parameter value of 0.771.
Particle Backscatter Coefficient x 10 -5 (m -1 srad -1 ) Height a.s.l.(m) Figure 5 Fig. 5. Vertical profiles of particle backscatter coefficient retrieved for the selected regions in Fig. 3 of: (i) drizzle fallstreak (dashed line), (ii), ice crystals (dotted line) and (iii) cloud-free atmosphere (solid line) from the AMALi measurements at 532 nm.Some CPI images exemplify water drops and ice crystals to related-profiles.

An attempt to compare remote and in situ observations
We recall that the particle backscatter coefficient profiles were obtained using the iterative method (Stachlewska, 2006a, see Appendix) with lidar ratios evaluated by using the method of Jourdan et al. (2003) applied to the scattering phase functions measured by the Polar Nephelometer and displayed on Figs.4a and b.At a wavelength of 532 nm the lidar ratios (i.e., extinction-to-backscatter ratio) are 14.5 sr and 13.8 sr for the drizzle drops and the ice crystals, respectively.We note in passing that similar values of lidar ratios are found for both water droplets and ice crystals and agree well to values expected in water and cirrus clouds (Sassen et al., 1989;Ansmann et al., 1992;Chen et al., 2002).Three vertical profiles of particle backscatter coefficient obtained from the lidar signals in selected regions of Fig. 2 are displayed on Fig. 5.In order to get statistically representative observations on large volumes, these profiles represent mean values calculated over 2 mn for the drizzle-dropfallstreak sample at 09:19 UT (see Fig. 2) and over 7 mn in the ice crystal layer at 09:13 UT (see Fig. 2).The profile at 09:37 related to a cloud-free atmosphere has been calculated as a background profile which allows to better identify the cloud heights and the transition between the cloud particles and background aerosols.
The results in Fig. 5 show that in the drizzle region the particle backscatter profile can be retrieved down to the sea level despite the large peak of the backscatter coefficient (4.5×10 −5 m −1 sr −1 ) evidenced between 1300 m and 1400 m.This layer is likely a remnant part of the strat-iform water cloud in which the drizzle precipitation was formed by coalescence process and then has washed-out the cloud.At lower levels, much lower values of the backscatter coefficient characterize the drizzle precipitation, i.e. from 0.4×10 −5 m −1 sr −1 to 0.9×10 −5 m −1 sr −1 .Table 1 summarizes the retrieved particle backscatter coefficients related to the drizzle cloud parts and the corresponding particle extinction coefficient calculated by multiplication with indicated lidar ratios.Compared to the extinction coefficient derived from the Polar Nephelometer during the 10:30-10:31:30 sequence (see mean value and standard deviation on Fig. 4a and in Table 1) the remote values are smaller (0.65 km −1 versus 3.5 km −1 ) but are overlapping within the range of horizontal variations of the in-situ measurements and within the given uncertainties.
The profile of the retrieved backscattering coefficient related to the ice crystal layer between 1400 m and 3000 m (at 09:13 UT, see Fig. 5) highlights vertical variation of the layer structure between 0.4×10 −5 m −1 sr −1 and 0.6×10 −5 m −1 sr −1 or between 0.06 and 0.08 km −1 in term of extinction coefficient.Comparisons in Table 1 also show smaller remote values versus direct measurements (0.20 km −1 ) with variation ranged within the in situ variability of the horizontal layer structure.
In conclusion, the extinction values significantly disagree between the two instruments for both cases (drizzle particles and ice crystals).The main reason to explain these discrepancies could be the non-accurate co-location of the two measurements.Indeed, a time separation of about 80 mn between lidar and in situ observations can lead to large differences in the results due both to the time-evolution of the cloud and to the airmass advection.A lower sensitivity of the in situ measurements for the detection of small particle concentration may also explain parts of the differences in the compared extinction values.Nevertheless, the signatures of these kinds of clouds are clearly evidenced from the in situ measurements and from the lidar profiles in terms of backscatter and Depolarisation ratio.With the prospect of the validation of the retrieved products (i.e.extinction coefficient, particle phase, effective diameter, . . . ) of the new generation of the A-Train satellites (including CALIPSO), an airborne plate-form for in situ observations, like the one presented in this paper, could be an pertinent tool.The experiment strategy should consider co-located observations with the satellite traces in order to get relevant results.

Conclusions
During the international ASTAR experiment (Arctic Study of Aerosols, Clouds and Radiation) carried out from Longyearbyen (Spitsbergen) a case study related to observations with ice crystals precipitating down to supercooled boundary layer stratocumulus has been performed from airborne measurements provided by an unique combination of remote and in situ instruments.Quasi-simultaneous alternated remote observations of clouds vertical and horizontal structures and Depolarisation effects due to presence of ice particles were combined with in-situ microphysics and optical observations of cloud particles.Three kinds of clouds with different microphysical and optical properties have therefore been alternately observed: (i) water droplets stratiform layer, (ii) drizzle-drops fallstreak and (iii) precipitating ice crystals from cirrus above.The signatures of these clouds are clearly evidenced from the in situ measurements and from the lidar profiles in term of backscatter and Depolarisation ratio.Accordingly, typical lidar ratios are derived from the measured scattering phase function combined with subsequent particle shapes and size distributions.The remote extinction profiles can therefore be retrieved under favourable conditions of low optical density and compared with the direct in situ measurements, though a further improvement of the methodology needs more coordinated observations in space and time.With the prospect of the validation of the retrieved products (i.e.extinction coefficient, particle phase, effective diameter, . . . ) of the new generation of the A-Train satellites (including CALIPSO), an airborne plate-form for in situ observations, like the one presented in this paper, could be an pertinent tool.calculation of the particle backscatter coefficient profile using the Klett's inversion β part Klett (h).Then, the calculated value of the particle backscatter coefficient at the height close to the aircraft β part Klett (h f ) is compared to the estimatedβ part (h 0 ).If they disagreed by more than 1%, the value of β part Klett (h ng ) is modified according to the Newton method.The β part Klett (h) profile is recalculated and both values were compared again.These iterations with changed β part Klett (h ng ) are performed until β part Klett (h f ) and β part (h 0 ) matched.For these calculations a lidar ratio has to be assumed.For most of the AMALi data under clean clear-sky conditions, the constant lidar ratio of 20-25 sr was assumed (Stachlewska, 2006b).Because the lidar ratio depends on the particle phase and shape, typical values were evaluated using the approach reported by Jourdan et al. (2003) from measured phase functions by the Polar Nephelometer (see Sect. 5).At the wavelength of 532 nm, the lidar ratios are of 14.5 sr and of 13.8 sr for the water and ice clouds, respectively.
The particle extinction coefficients profiles α part (h) were calculated by multiplication of the β part Klett (h) with a lidar ratio B(h).It was assumed that the lidar ratio profile depends on the type of particles in a cloud.In the areas where the Depolarisation ratio was large and hence indicated an existence of the ice particles, the lidar ratio corresponding to the value obtained from in situ measurements for the ice crystals was applied.In the areas of low Depolarisation ratio, the lidar ratio corresponding to the in situ value obtained for the water particles was used.
When a lidar probes in optically thick aerosol layers or clouds, measurements are affected by presence of higher order scattering.In such cases the Depolarisation ratio can be as large as 50%.For nadir-aiming short-range AMALi we assume that the amount and order of the multiple scattering depend mainly on the range from the lidar, the particle size and the optical depth of the particle layer.The optical depth threshold for a layer from the aircraft altitude downwards was assumed in this case study.The value of the threshold was of 1 for the case study presented here.In the areas were cloud system was thicker, the retrieval was not possible (for example, below the water cloud).

A2 Derivation of the Depolarisation ratio
Depolarisation ratio DR(h) profiles were calculated as the ratio of the background-corrected signals in perpendicular S ⊥ (h) and parallel S II (h) channels.The DR(h) profiles were normalised just below the overlap range of 235 m to a value of 1.44% assumed for the Depolarisation that is due to the air anisotropy (including the Cabannes line and the rotational Raman lines).Identical detectors with exactly the same settings (PMT high voltage 850 V, no neutral density filters) were used.For instrumental calibration, the gain ratio G 1 =P 1 /P 2 was obtained for the raw signals P 1 and P 2 measured on each detector.Then, the detectors were physically exchanged and the raw signals P 1 and P 2 were measured to obtain G 2 =P 1 /P 2 .The correction gain G=(G 1 +G 2 )/2 was not dependent on the polarization and small changes of source light intensity during the measurement.The dependence on changes in surrounding temperature was negligible (less than 1%/C • ).High quality of optics used for the polarization separation and accurate adjustment of optical elements assured very low cross-talk between the perpendicular and parallel channels, i.e., less than 10 −6 of "wrong" polarization was transmitted on each component.In order to avoid the effect of horizontally oriented ice plates on the Depolarisation ratio, the AMALi was always pointed in nearly the nadir-looking direction (aircraft's pitch angle ∼7 • ).Linear Depolarisation ratio measurements provided a way to distinguish well between spherical DR≈0 and non-spherical 0<DR<1 particles if the assumption of single scattering holds.

A3 Discussion on the errors
In case of the AMALi lidar, for the height chosen just below the overlap range of 235 m, the error associated with neglecting the transmittance term varies from 0.7% for the low particle extinction of 1.5×10 −5 m −1 , which is typical for a clear Arctic atmosphere, up to 2.8% for the particle extinction of 6.×10 −5 m −1 .In the cloud case discussed in this paper, the latter uncertainty was considered.Uncertainties in the particle backscatter coefficient were estimated accordingly to the error propagation.In the case of the nadir-looking shortrange lidar, the largest term of the error propagation of the backscatter coefficient profiles is due to the wrong assumption of the backscatter calibration value β part (h 0 ), while the wrong assumption of the lidar ratio contributes much less to the uncertainty (small values of dβ/dB).Due to the short range for the AMALi measurements (the flight altitude < 3 km, and the signal-to-nose ratio >35) the assumption of the constant lidar ratio does not introduce a significant error.The accuracy of the backscatter coefficient calculated with the iterative approach for the different arctic aerosol types was discussed in details by Stachlewska (2006a).According to that study, the highest error value was found to be less than the particle backscatter error of 2.0×10 −7 m −1 sr −1 .
The molecular backscatter coefficient profiles used for the calibration of the particle backscatter coefficient profiles are obtained as a climatology monthly mean of air density and temperature profiles for the Arctic region from the daily radiosonding data base of the AWI Koldewey Station in Ny Ålesund, Spitsbergen.The molecular backscatter profile accuracy is of ±4.2×10 −8 m −1 sr −1 .For the cloudy areas the uncertainty of this calibration is negligible.
As for the extinction coefficient profiles, the deviation of the estimated lidar ratio value from the real one is transformed directly into an offset of the extinction coefficient.

Fig. 1 .
Fig. 1.Map of the Svalbard Archipelago and location of the Polar 2 aircraft flight (grey line) on 5 June 2004 during the ASTAR campaign.The ASTAR flights were made from Longyearbyen (LYR).

Fig. 2 .
Fig. 2. (a) Vertical cross-section of the range corrected AMALi lidar signal at 532 nm for measurements carried out between 08:57 and 09:42 UT on 5 June 2004.(b) Corresponding vertical cross-section of the Depolarisation ratio at 532nm.The blank areas represent nonexploited data due to strongly attenuated signal by water droplet clouds (see text for explanation).

Figure 3
Figure 3 displays time-series of cloud parameters measured between 10:25 and 10:55 UT at a single flight level (1400 m/−11• C).The parameters are: the particle concentration (C25) and the effective diameter (Deff) measured by the PMS 2D-C probe, the extinction coefficient (Ext) and the asymmetry parameter (g) derived from the Polar Nephelometer measurements and the liquid water content (LWC) deduced from the Nevzorov instrument.Shape classification of particles larger than 50 µm from CPI images is also shown in Fig.3with three categories, i.e. water spherical droplets, regular shape (columns, plates, . . . ) and irregular shape particles.Of course, due to the time differences in sampling times (about 1 h 20 min) the in situ measurements cannot be directly related to the remote observations in Fig.2because of both the natural time-evolution of the cloud and the airmass advection.Nevertheless at the flight level, the wind being roughly parallel-oriented to the flight trajectory (5 • , see Fig.1) the airmass advection is about 30 km according to the sampling time difference and a wind speed of 6 m/s.This leads us to the conclusion that the aircraft sampled the same cloud system which may retain similar features between the two flight sequences.Indeed, Fig.3shows that the aircraft sampled distinct cloud patches with different features Figure 4

Table 1 .
Values of the particle backscatter and extinction coefficient with uncertainties (error discussion in text) retrieved from the AMALi observations in three different clouds.The values of the particle extinction coefficient and the lidar ratio derived from the Polar Nephelometer measurements (mean values and standard deviations) are also reported.