Application of the shipborne remote sensing supersite OCEANET for profiling of Arctic aerosols and clouds during Polarstern cruise PS106

From 25 May to 21 July 2017, the research vessel Polarstern performed the cruise PS106 to the high Arctic in the region north and northeast of Svalbard. The mobile remote-sensing platform OCEANET was deployed aboard Polarstern. Within a single container, OCEANET houses state-of-the-art remote-sensing equipment, including a multiwavelength Raman polarization lidar PollyXT and a 14channel microwave radiometer HATPRO (Humidity And Temperature PROfiler). For the cruise PS106, the measurements were supplemented by a motion-stabilized 35 GHz cloud radar Mira-35. This paper describes the treatment of technical challenges which were immanent during the deployment of OCEANET in the high Arctic. This includes the description of the motion stabilization of the cloud radar Mira-35 to ensure vertical-pointing observations aboard the moving Polarstern as well as the applied correction of the vessels heave rate to provide valid Doppler velocities. The correction ensured a leveling accuracy of ±0.5 during transits through the ice and an ice floe camp. The applied heave correction reduced the signal induced by the vertical movement of the cloud radar in the PSD of the Doppler velocity by a factor of 15. Low-level clouds, in addition, frequently prevented a continuous analysis of cloud conditions from synergies of lidar and radar within Cloudnet, because the technically determined lowest detection height of Mira-35 was 165 m above sea level. To overcome this obstacle, an approach for identification of the cloud presence solely based on data from the near-field receiver of PollyXT at heights from 50 m and 165 m above sea level is presented. We found low-level stratus clouds, which were below the lowest detection range of most automatic ground-based remote-sensing instruments during 25 % of the observation time. We present case studies of aerosol and cloud studies to introduce the capabilities of the data set. In addition, new approaches for ice crystal effective radius and eddy dissipation rates from cloud radar measurements and the retrieval of aerosol optical and microphysical properties from the observations of PollyXT are introduced.

Abstract. From 25 May to 21 July 2017, the research vessel Polarstern performed the cruise PS106 to the high Arctic in the region north and northeast of Svalbard. The mobile remote-sensing platform OCEANET was deployed aboard Polarstern. Within a single container, OCEANET houses state-of-the-art remote-sensing equipment, including a multiwavelength Raman polarization lidar Polly XT and a 14channel microwave radiometer HATPRO (Humidity And Temperature PROfiler). For the cruise PS106, the measurements were supplemented by a motion-stabilized 35 GHz cloud radar Mira-35. This paper describes the treatment of technical challenges which were immanent during the deployment of OCEANET in the high Arctic. This includes the description of the motion stabilization of the cloud radar Mira-35 to ensure vertical-pointing observations aboard the moving Polarstern as well as the applied correction of the vessels heave rate to provide valid Doppler velocities. The correction ensured a leveling accuracy of ±0.5 • during transits through the ice and an ice floe camp. The applied heave correction reduced the signal induced by the vertical movement of the cloud radar in the PSD of the Doppler velocity by a factor of 15. Low-level clouds, in addition, frequently prevented a continuous analysis of cloud conditions from synergies of lidar and radar within Cloudnet, because the technically determined lowest detection height of Mira-35 was 165 m above sea level. To overcome this obstacle, an approach for identification of the cloud presence solely based on data from the near-field receiver of Polly XT at heights from 50 m and 165 m above sea level is presented. We found low-level stratus clouds, which were below the lowest detection range of most automatic ground-based remote-sensing instruments during 25 % of the observation time.
We present case studies of aerosol and cloud studies to introduce the capabilities of the data set. In addition, new approaches for ice crystal effective radius and eddy dissipation rates from cloud radar measurements and the retrieval of aerosol optical and microphysical properties from the observations of Polly XT are introduced.

Introduction
The Arctic is one of the hotspots of global climate change. This is observed as a change of several parameters, such as the drastic decline of the Arctic sea ice during all seasons, but especially in summer, in both extent and thickness (Meier et al., 2014). Also, in the past 30 years, the mean Arctic near-surface air-temperature anomaly increased at least by a factor of 2 faster compared to the global mean (Serreze and Barry, 2011). These phenomena, summarized by the term Arctic amplification, are assumed to be due to several feedback mechanisms, e.g., the surface albedo feedback, the lapse rate feedback, a change in the meridional atmospheric and oceanic mass and energy transport pattern, and an alteration in cloud cover, aerosol occurrence, and atmospheric moisture content . However, there is still a lack of understanding in the interplay of these feedback mechanisms as well as in quantifying their relative importance and magnitude (Serreze and Barry, 2011;Pithan and Mauritsen, 2014;Goosse et al., 2018).
The radiative effect of clouds is a major source of uncertainty in this matter. Arctic clouds have a high variability in their radiative effects and in their impact on the surface energy balance (Yeo et al., 2018). The cloud-related radiative impacts have been found to be both positive (i.e., clouds have a warming effect) and negative (Goosse et al., 2018). A higher fraction of low clouds, e.g., caused by a decreased sea ice extent, increases the downwelling longwave radiation during polar night and thus induces a positive feedback. A higher fraction of liquid water in mixed-phase clouds due to a warmer climate, on the other hand, was found to increase the cloud albedo. This in turn enhances the reflection of incoming shortwave radiation at the top of the atmosphere during polar day (Goosse et al., 2018) and thus produces a negative feedback. Yet the underlying processes controlling Arctic cloud phase and occurrence -and hence the connected feedback mechanisms driving Arctic amplificationare not completely understood (e.g., Shupe et al., 2013;Kalesse et al., 2016a).
Though they are a key requirement to study the Arctic energy budget, detailed observations of Arctic clouds are still rare. Measurements from ground-based stations, such as the observatories of the International Arctic Systems for Observing the Atmosphere , are of great value, e.g., due to the possibility of conducting long-term observations. However, they are limited to their location and influenced by their surrounding orography. Drifting buoys on the other hand can enter any place in the Arctic ice and thus are very valuable in this harsh environment. Their equipment gets increasingly sophisticated, and for a few years buoys equipped with autonomous lidar systems (Mariage et al., 2017) have been in use, giving them the opportunity to measure vertical profiles of the atmospheric column. But they still are limited in their payload and cannot yet replace measurements from observatories or campaigns. To study Arctic clouds, different aircraft campaigns have been conducted in recent years (e.g., Curry et al., 2000;Jacob et al., 2010;Mc-Farquhar et al., 2011;Wendisch et al., 2019). While airborne measurements yield an unique, accurate description of the observed cloud, they lack the ability to measure continuously the entire tropospheric column over a long period, a feature active remote sensing observations can offer. Given this capability, ground-based remote sensing observations are suitable to investigate the spatiotemporal distribution of clouds (Bühl et al., 2013), their phase partitioning (de Boer et al., 2009;Zhang et al., 2014;Kalesse et al., 2016a), and their interaction with aerosols (Seifert et al., 2010). These data sets serve, e.g., as a basis for model evaluation (Illingworth et al., 2007;Neggers et al., 2019) and radiative transfer calculations Barrientos Velasco et al., 2020). Hence, in addition to the airborne campaigns, several shipborne cam-paigns equipped with remote-sensing instrumentation have been conducted in the past few years in the Arctic (e.g., Uttal et al., 2002;Tjernström et al., 2004Tjernström et al., , 2014Granskog et al., 2018;Wendisch et al., 2019). Observations of spaceborne cloud radar and lidar, as done aboard Cloudsat (Stephens et al., 2008) and CALIPSO (Winker et al., 2010) provide, in addition, a large-scale overview of the Arctic cloud coverage (Liu et al., 2012). But the respective data sets lack information about the lowest cloud levels. Nevertheless, there are still only a few studies of sea-motion-stabilized cloud radars, whose availability is a necessary requirement to also determine vertical cloud dynamics accurately from a shipborne platform.
In order to study the feedback mechanisms causing Arctic amplification, the initiative ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes and Feedback Mechanisms (AC) 3 conducted two complementary field campaigns in the Arctic summer of 2017: Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD), an airborne campaign performed with the research aircraft Polar 5 and Polar 6, and the Physical feedbacks of Arctic planetary boundary layer, Sea ice, Cloud and AerosoL (PASCAL) expedition, deployed on and around the research ice breaker Polarstern (Macke and Flores, 2018;Wendisch et al., 2019). These campaigns took place in May and June 2017 in the regions north and northeast of Svalbard with the aim to combine remote sensing and in situ observations. During PASCAL, a 2-week ice floe camp was performed in the vicinity of Polarstern, and a large number of auxiliary measurements were conducted on the ice. PASCAL was the first part of the split Polarstern cruise PS106, which took place from 25 May until 21 July 2017. During the whole PS106 cruise, measurements with the multiwavelength polarization lidar Polly XT _OCEANET, a 35 GHz cloud radar Mira-35, and a microwave radiometer HATPRO (Humidity And Temperature PROfiler) of the OCEANET platform were conducted aboard Polarstern. Within (AC) 3 the OCEANET observations have the essential role to describe the temporal evolution of the vertical structure of aerosol and clouds in the central Arctic. They constitute the prerequisite for further studies of aerosol-cloud interaction, model evaluation, or radiative transfer modeling, which are partly covered by other subprojects of (AC) 3 . The scope of this study is thus to introduce the instrumentation and data analysis methods which are used to produce the OCEANET-based cloud and aerosol data sets for the cruise PS106.
In Sect. 2 of this paper, a detailed description of the OCEANET instruments and the auxiliary observations is given. The applied motion-stabilization and heave-rate correction of the cloud radar Doppler velocity, the data processing based on the synergistic Cloudnet algorithm (Illingworth et al., 2007), and the development of auxiliary retrievals for processing within Cloudnet are described in Sect. 3. In Sect. 4 the products derived from the OCEANET measurements using Cloudnet are illustrated by means of different case studies from the time period of the ice flow camp. The potential of Polly XT _OCEANET to characterize the freetropospheric aerosol is also highlighted. In addition, a statistical overview of the observed vertical cloud structure with a special focus on low-level clouds during PS106 is presented. A final summary and conclusions are given in Sect. 5.

Instrumentation
During the complete PS106 campaign in the central Arctic in summer 2017 (see Fig. 1), a comprehensive number of remote-sensing instruments were deployed aboard the research vessel (RV) Polarstern to conduct continuous observations of clouds and aerosols. To a large extent, these instruments comprised the OCEANET-Atmosphere observatory (Kanitz et al., 2013a). Additionally, auxiliary instruments for in situ and remote sensing observations installed aboard Polarstern as well as used during a 2-week ice floe-camp, which was performed in the vicinity of the RV, were utilized for the studies presented in here. The first part of PS106 (PS106.1/PASCAL) was accompanied by the ACLOUD aircraft campaign (Ehrlich et al., 2019) about both of which a brief overview is given by Wendisch et al. (2019).
The locations of the OCEANET equipment and of the auxiliary Polarstern instruments deployed during PS106 that were used within this study are depicted in Fig. 2. Table 1 summarizes the technical details of the key equipment applied in the synergistic Cloudnet processing that is further described in Sect. 3.3. In the following, the different instruments will be briefly introduced.

OCEANET
The OCEANET-Atmosphere observatory was already frequently operated aboard Polarstern (Kanitz et al., 2013b;Bohlmann et al., 2018;Yin et al., 2019) -yet, so far only for the transects from the Northern Hemisphere to the Southern Hemisphere (or vice versa) but never in the Polar regions. Its container is by default equipped with the multiwavelength polarization Raman lidar Polly XT _OCEANET (hereafter referred to as Polly XT ), in order to provide continuous profiles of cloud and aerosol properties . Additionally, a 14-channel microwave radiometer (MWR) HATPRO (Rose et al., 2005) for measurements of columnintegrated liquid water and water vapor content and profiles of atmospheric temperature, a standard meteorological station, a pyranometer and a pyrgeometer for incoming shortand longwave radiation observations, and an all-sky camera for passive visible observations of the full sky were installed. During PS106, OCEANET was complemented for the first time by a vertically pointing motion-stabilized 35 GHz polarimetric Doppler cloud radar of type Mira-35 (Görsdorf et al., 2015) for continuous vertically resolved measurements of Doppler spectra produced by vertical cloud motions. The Polly XT system measures profiles of particle backscatter coefficient at three wavelengths (355, 532, and 1064 nm) and profiles of extinction coefficient as well as the linear depolarization ratio at two wavelengths (355 and 532 nm), respectively -for details see Baars et al. (2016). Four near-field channels for detection of elastically and Ramanscattered light from nitrogen molecules are implemented at 355, 387, 532, and 607 nm to enable observations already at low heights, starting at about 50 m above the instrument. An additional channel for detection of Raman-scattered light from water vapor at 407 nm allows the retrieval of the water vapor mixing ratio (Dai et al., 2018) during low-sunlight con-  ditions. From the Polly XT backscatter and extinction measurements, aerosol classification by their optical properties (Müller et al., 2007;Baars et al., 2017) up to the retrieval of particle size distribution and number concentration (Müller et al., 1999;Baars et al., 2012) can potentially be derived. The polarization-sensitive detection channels allow us to distinguish between spherical and nonspherical aerosol and cloud particles (Kanitz et al., 2013a) and, for instance, to separate dust and non-dust particles in mixed aerosol layers (Baars et al., 2011). By applying the shape-detection capabilities of the polarization channels for the discrimination of spherical liquid droplets from nonspherical ice particles, heterogeneous ice formation in mixed-phase clouds can be studied (Seifert et al., 2015). Another application of depolarization observations in mixed-phase cloud studies is the estimation of cloud condensation nuclei (CCN) and icenucleating particle (INP) concentrations . Due to the relatively short wavelengths of the lidar, e.g., compared to the cloud radar, it follows that the lidar is sensitive to rather small particles such as aerosols or small cloud droplets. Also attenuation, especially due to liquid clouds, has to be considered.
The MWR HATPRO provides estimates of the liquid water path (LWP), integrated water vapor (IWV), and humidity and temperature profiles with a temporal resolution of 1 Hz. The MWR measures the emission of radiation from the atmosphere in two frequency bands ranging from 22.24 to 31.4 GHz and from 51.0 to 58.0 GHz at 14 different channels. The MWR data sets shown in this study are based on a retrieval that was created based on a long-term radiosonde data set from Ny-Ålesund, Svalbard, Norway (78.9 • N, 11.8 • E, 11 m height above sea level, a.s.l.; WMO code 6260) according to Löhnert and Crewell (2003).
During PS106, the Ka-band Doppler radar Mira-35 was set up to emit pulses with a width of 208 ns at a pulse rep-etition frequency of 5000 Hz. This corresponds to a vertical resolution of 31.18 m. The upper limit of the measurement range was set to 15 km. The Doppler spectrum was derived from the backscattered signals of 256 consecutive pulses. To allow the correction of the cloud radar data for the vessel movement, the whole spectrum (including noise) has been stored with a temporal resolution of 4 Hz and a Doppler resolution of 0.08 m s −1 . This correction has been done in a postprocessing procedure which is explained in Sect. 3.1. From the profiles of the Doppler spectra, the different Doppler moments such as radar reflectivity, Doppler velocity, and Doppler spectral width were determined as described in Görsdorf et al. (2015). The linear depolarization ratio (LDR) was obtained from the ratio of the radar reflectivity factor observed in the co-and cross-channels of Mira-35 and provides information about the hydrometeor shape (Bühl et al., 2016). In contrast to the lidar, the longer wavelength of operation of the cloud radar defines its sensitivity to range from cloud hydrometeors to slight precipitation. Especially in the case of shallow stratiform clouds, as they dominated the measurements during PS106, attenuation effects can be neglected. The OCEANET data sets of HATPRO, Polly XT , and Mira-35 are publicly available through the Open Access library PANGAEA (Griesche et al., 2020c(Griesche et al., , d, 2019.

Auxiliary instrumentation
Added value of the OCEANET measurements can be obtained when they are accompanied by additional observations. During the 2-week ice floe camp performed in the frame of PASCAL, a tethered balloon site was set up for turbulence and radiation observations , and a network covering 15 pyranometers to determine the spatial variability in the solar radiation was installed (Barrientos Velasco et al., 2020). In the context of this study, the turbulence as determined from the three-dimensional wind vector measured with high temporal resolution of several tens of hertz by an ultrasonic anemometer attached to the tether of the balloon was used. To obtain mass and number concentration as well as optical properties and filter samples of the aerosol at the surface, a container equipped with instrumentation for aerosol in situ measurements was installed on the deck of Polarstern and was measuring continuously during the whole 2-month cruise (Kecorius et al., 2019). Also aboard Polarstern, measurements of the optical thickness of the cloud-free atmosphere were performed using a hand-held Solar Light MICROTOPS sun photometer. The sun photometer measurements are already available through the AErosol RObotic NETwork (AERONET) project. An optical disdrometer, which is part of the OceanRAIN network (Klepp et al., 2018), mounted on the crows nest of the RV was continuously measuring the precipitation rate for different hydrometeor types and size bins. Additionally, launches of Vaisala RS92-SGP radiosondes (Jensen et al., 2016) were conducted every 6 h (shortly before 05:00, 11:00, 17:00, and 23:00 UTC to reach 100 hPa approximately at 06:00, 12:00, 18:00, and 24:00 UTC) to obtain in situ profiles of temperature ( T = 0.5 • C), relative humidity ( RH = 5 %), pressure ( p = 1 hPa), and horizontal wind speed ( v = 0.15 m s −1 ) and direction ( • = 2 • ).

Data processing and synergistic retrievals
Aim of the OCEANET observations from PS106 was to provide a continuous vertically resolved view on cloud and aerosol macro-and microphysical properties in order to enhance the understanding of the Arctic atmosphere system and to support partner projects with data sets for radiative transfer calculations and turbulence studies. To derive continuous products of cloud and aerosol properties, the shipborne OCEANET remote sensing observations were processed using the synergistic retrieval algorithm Cloudnet (Illingworth et al., 2007). In this section, we describe the extension of the standard Cloudnet algorithms by additional simple but operationally applicable products providing estimates of cloud droplet and ice crystal effective radius and the cloud-turbulence parameter eddy dissipation rate (EDR). The procedure for minimizing the influence of the RV motion on the measurement of vertical velocities with Mira-35, which are required for the EDR retrieval, is also explained below.

Correction of vertical-pointing cloud radar observations for ship motion
A structural requirement to derive valid vertical velocity from a Doppler cloud radar is a vertical-pointing radar without self-vertical-velocity component. When the cloud radar is pointing off zenith, the measured vertical-pointing Doppler velocity will be biased by an additional component intro- duced by the horizontal wind. Based on highly resolved horizontal wind data and the radar beam incident angle, a correction is possible for this bias (Wulfmeyer and Janjić, 2005). For PS106, a different approach was chosen. Similar to the approach described by Achtert et al. (2015), the cloud radar was mounted on an active stabilization platform (Fig. 3a), which was in our case a predecessor of the SOMAG AG Jena -GSM 4000 (SOMAG, 2017). This platform actively leveled out the roll and pitch movement of the RV, ideally in a way that no correction of horizontal-wind effects was necessary. Figure 4 shows a comparison of the pitch and roll angle time series during ice breaking conditions from 1 June 2017 07:00 UTC-3 June 2017 08:00 UTC measured by the vessel's own inertial measurement unit (IMU) and directly at the cloud radar. As the platform itself did not provide any position determination, we made use of a single-board computer (Beaglebone Blue) with integrated IMU. During the ice transit and the ice floe camp periods, the stabilization platform ensured an accuracy of the leveling of ±0.5 • . The 2σ standard deviation during the ice transit (ice floe camp) was found to be 0.32 • (0.34 • ) -thus 95 % of the data points show an accuracy of 89.68 • (89.66 • ). During the open-sea passage of RV Polarstern, the accuracy of the stabilization was reduced to around ±1 • with a 2σ standard deviation of 0.7 • .
An additional bias of the true Doppler velocity can occur if the cloud radar itself moves vertically; the vertical velocity superimposes the measured Doppler velocity. In the case of a moving RV, the vertical movement is induced by the RV's heave rate and rotation. The necessary heave correction was done in a postprocessing procedure which will now be introduced.
To enable the correction, the complete cloud radar Doppler spectra as well as the motion data (rotation and translation) of Polarstern were stored with a resolution of 4 and 20 Hz, respectively, throughout the entire cruise. To determine the heave rate of the cloud radar v C z , first the cross product of the rotation vector of Polarstern v P R (angular velocity) and the position of the radar (X R = −5.4, −2.34, 16.35 m) relative to the mass center of the RV was calculated as follows: (1) After that, v R PS had to be transformed into the earth coordinate system by means of the rotation matrices. For the z component, this results in where pitch and roll represent the respective displacement, and • denotes the dot product. The heave rate of the cloud radar v C z was then calculated by the sum of v R z with the z component of the translation vector of Polarstern v P T,z : (3) Figure 5. Flowchart of the heave rate correction.
In Fig. 3b, the time series of v C z for PS106 is shown. The heave rate correction was done by shifting each individual Doppler spectrum opposite to the cloud radar heave rate. An illustration of this procedure is shown in Fig. 5. In an initial step, the cross-correlation between the timestamps of the two data sets, the cloud radar Doppler spectrum and the cloud radar heave rate, was calculated to check for a possible time shift between both data sets. This was found to be 0.25 s. Subsequently, the two values of v C z from before and after the current Doppler spectrum have been linearly interpolated onto its respective time. Finally, the spectrum was shifted according to the number of spectral bins determined by the Doppler resolution ( v DopplerSpectrum = 0.08 m s −1 ) and the interpolated heave rate.
The effect of the heave correction is illustrated in Fig. 6. In Fig. 6a the uncorrected Doppler velocity measured on 30 May 2017 between 00:00 and 01:00 UTC, together with the respective histogram of the velocities, is shown. The RV's movement is visible in both, in the time-height cross section of the Doppler velocity as stripes of enhanced or reduced velocity throughout the whole column as well as in the broadening of the histogram. The same is presented in Fig. 6b but for the corrected Doppler velocity.
To evaluate the effect of the heave correction, we calculated the Fourier spectrum of the corrected and uncorrected Doppler velocity (Fig. 6e, f). Continuous time series of 1 h of Doppler velocity in the highest possible range gates of the cloud were analyzed. To quantify the impact of the heave correction, the integral of the frequency range which was most affected by the ship's movement (0.1-2 Hz) was calculated both for the corrected and for the uncorrected data. The heave correction reduced the fraction of the ship's movement in the power spectral density of the cloud radar Doppler velocity by a factor of 15.

Retrieval of eddy dissipation rate
The rate at which turbulence kinetic energy is transferred from larger eddies into smaller ones and eventually dis- solves into thermal energy is the EDR. This is used as a quantitative proxy of atmospheric turbulence. Several approaches to retrieve the EDR are common. Methods exist for in situ measurements from aircraft- (Nicholls, 1978;Nucciarone and Young, 1991;Meischner et al., 2001), helicopter- (Siebert et al., 2006a), and balloon-borne (Caughey et al., 1979;Siebert et al., 2006b) instruments and for meteorological tower instruments (Caughey et al., 1979;Kaimal et al., 1976;Zhou et al., 1985). Additional retrievals for remote sensing observation have been developed (Borque et al., 2016;Sathe and Mann, 2013). These methods are based on the Doppler velocity structure function derived from vertically pointed Doppler lidar (Frehlich and Cornman, 2002) or Doppler radar (Lothon et al., 2005) or a combination of the width of the Doppler spectrum and the Doppler velocity measurements (Meischner et al., 2001). Other retrievals use time series analyses of vertical velocities from vertical-pointing Doppler radar Kalesse and Kollias, 2013) or Doppler lidar observations (O'Connor et al., 2010).
The presented range of EDR for different cloud conditions also suggests that Arctic clouds might show characteristic differences for varying atmospheric conditions. The vertical alignment of the cloud radar during PS106 enables the determination of EDR from the vertical air motions observed in cloud layers. Below, we thus present a retrieval technique for EDR that can be applied to the OCEANET data set.

EDR from vertical-pointing Doppler velocity power spectra
Assuming the turbulent energy dissipation is a homogeneous and isotropic process, the turbulent energy spectrum S(k) within its inertial subrange is represented according to Borque et al. (2016) by with A = 0.5 the Kolmogorov constant for 1-D wind spectra (Sreenivasan, 1995). k represents the wavenumber, which is related to a length scale L (k = 2π/L) as well as to frequency f with k = f/V h and V h as the horizontal wind speed and assuming a linear wind field. If in a log-log plot the observed spectra within the inertial subrange follow a −5/3 slope, ε can be estimated by where k 0 is the corresponding intercept of the linearized fit. For this study, power spectra of the Doppler velocity with 4 Hz of continuous time series covering 5 min were calculated. To get the best estimate of the respective inertial subrange, the fit was determined by calculating a linear leastsquares regression of the spectrum in 34 different wavenumber intervals. The corresponding wavenumber intervals k i are depicted in Fig. 7a together with the spectrum of the vertical velocity observed on 7 June 2017, from 10:28 to 10:43 UTC. Following Borque et al. (2016), good fit was defined with a slope from the linear regression of −5/3 ± 20 % (−5/3±1/3). If this criterion was matched within more than one wavenumber interval, the mean of all ε i for one spectrum was calculated. In order to evaluate the EDR estimated by cloud radar measurements, it was compared to EDR derived from the tethered balloon . The time periods used for deriving EDR from the tethered balloon were 15 min, during which it was located at a constant height above the ground. As a measure to quantify the uncertainties in the two retrievals, the standard deviation σ of all good fits was calculated.
Three comparisons had been done for situations where the tethered balloon was parked at a constant height within a cloud. One on 7 June 2017, with the tethered balloon being at 380 m height between 10:28 and 10:43 UTC. In Fig. 7a an intercomparison of the power spectrum derived by the tethered balloon measurements (red) with the spectrum derived from the cloud radar Doppler velocity (blue) according to the techniques described above is shown for this period. The other two comparisons were done for 5 June 2017 13:50-14:05 UTC at 330 m height, and for 9 June 2017 09:00-09:15 UTC at 500 m height. The EDR values for these cases from the tethered balloon measurements were 8.90 × 10 −5 ± 1.07×10 −5 and 6.39×10 −5 ±5.48×10 −6 m 2 s −3 , while the cloud radar measurements gave 5.98 × 10 −5 ± 3.53 × 10 −5 and 2.26 × 10 −5 ± 1.64 × 10 −5 m 2 s −3 , respectively. The comparisons of the two retrievals showed that the values differ by a factor of 2-3. This discrepancy is in the order of magnitude, as one could expect due to the spatial distance between the two measurements alone (about 200 m).

Cloud macro-and microphysical properties from instrument synergies
To acquire a data set suitable for the statistical evaluation of the macro-and microphysical properties of clouds observed during PS106, the instrument synergistic approach Cloudnet (Illingworth et al., 2007) was applied. This data set, in addition, serves to realize model evaluations (Illingworth et al., 2007) and radiative transfer calculations, e.g, with the Rapid Radiative Transfer Model for GCMs (RRTMG; Mlawer et al., 1997;Barker et al., 2003;Clough et al., 2005). RRTMG is currently utilized for single-column radiative transfer calculations. The model considers vertical profiles of relative humidity and temperature, standard atmospheric constituent profiles based on Anderson et al. (1986), and cloud macro and microphysical properties of clouds. These assignments include sets of effective radius and mass concentration of liquid and ice hydrometeors. In the following, the approaches for achieving these data set requirements based on the PS106 remote sensing observations are described.

Cloudnet
The instrument synergy approach Cloudnet (Illingworth et al., 2007), which combines the observations from lidar, cloud radar, microwave radiometer, disdrometer, and radiosondes, was used to determine cloud physical properties during PS106. To illustrate this procedure, the Cloudnet approach will now be briefly introduced. The measurements are first averaged on a common grid with a vertical and temporal resolution of 31.18 m (resulting from the cloud radar resolution) and 30 s. To estimate the temperature at the respective time-height pixel, radiosonde-based profiles of thermodynamic variables are interpolated on the Cloudnet grid. If no radiosonde was launched from the RV, but Polarstern was in the vicinity of Svalbard, soundings from Ny Ålesund (Maturilli, 2017) were utilized instead. As a last fall-back option, data from the Global Data Assimilation System model (GDAS, 2017) with a horizontal and vertical resolution of 1 • and 3 h (GDAS1) were used as meteorological input into Cloudnet.
Based on the observations scaled on the Cloudnet grid, a categorization bit mask is derived, which assigns a series of seven distinct features to the observed targets: clear yes/no; liquid yes/no; falling yes/no; wet bulb temperature below 0 • C yes/no; melting layer yes/no; aerosol yes/no; and insects yes/no. The bitwise categorization ensures that each data point is characterized by a defined combination of these features. The detailed definition of the respective categorization bits is beyond the scope of this paper and has already been given by Hogan and O'Connor (2004).
Based on the individual combination of the categorization bits, an atmospheric target classification is derived as follows: "clear sky" is defined as no bit is set for the respective pixel. "Cloud droplets only" is identified by only the droplet bit being set. The falling bit alone identifies "drizzle or rain". Droplet and falling bit together are interpreted as "drizzle/rain and cloud droplets" and the falling and cold bit together as "ice". Droplet, falling, and cold bit combined give "ice and supercooled droplets". The melting bit being set alone identifies "melting ice", and together with the droplet bit, the pixel is defined as "melting ice and cloud droplets". The aerosol and insects bits then are accordingly interpreted as "aerosol", "insects", or "aerosol and insects". Following previous studies (e.g., Shupe, 2011;Mioche et al., 2015), we have defined mixed-phase clouds when (supercooled) liquid water and ice particles are detected in the same data point and when an ice cloud is observed with a liquid or mixed-phase cloud-top layer.
Besides the phase of the cloud, the respective mass concentrations of ice and liquid water are determined where applicable. The liquid water content (LWC) is derived by scaling the MWR liquid water path adiabatically onto the cloud pixels defined as liquid or mixed phase (Frisch et al., 1998;Merk et al., 2016). For pure-liquid data points, the approach of Frisch et al. (2002) is used to derive the cloud droplet effective radius from the observed radar reflectivity factor and liquid water path and an assumed width of the lognormal cloud droplet size distribution (which was, according to Miles et al., 2000, set to 0.35 in our study). The ice water content (IWC) is calculated using an empirical formula from Hogan et al. (2006) relating cloud radar reflectivity Z and temperature T . This approach for IWC is only applied for clouds Cloudnet classified as ice or ice and supercooled droplets. Also, in this step a correction for potential attenuation of the cloud radar signal due to the presence of liquid water is made.

Ice crystal effective radius
As discussed above, Cloudnet offers a variety of retrievals for ice microphysical parameters. Nevertheless, the continuous application of radiative transfer calculations requires a consistent availability of ice and liquid hydrometeor effective radius and mass concentration. While Cloudnet already contains retrievals for effective radius and mass concentration of liquid droplets, as well for ice water content, so far no operational retrieval for ice effective radius is available. We therefore decided for the implementation of a new approach which is based on the combination of a definition of the effective radius as the ratio of the third to the second moment of the particle size distribution (PSD) and an empirical relationship between the visible extinction coefficient α, cloud radar reflectivity Z, and model temperature T . Similarly as for IWC (and α), r e ice is only calculated for data points where Cloudnet classified ice or ice and supercooled droplets.
Using the ratio of the second to the third moment of the PSD, the effective radius r e ice can be related to IWC and α (Delanoë et al., 2007). This yields for r e ice r e ice = 3 2 IWC ρ i α × 10 6 (µm), with ρ i as density of the solid ice (ρ i = 917 kg m −3 ). Both IWC and α have been calculated using empirical relationships between IWC or α and the cloud radar reflectivity Z of a 35 GHz cloud radar and temperature T (Hogan et al., 2006). Finally, we found for the ice crystal effective radius a Z-T relationship: r e ice = 3 2ρ i 10 C ZT ·ZT +C Z ·Z+C T ·T +C × 10 6 (µm), with C ZT = −2.05×10 −4 , C Z = 1.6×10 −3 , C T = −1.71× 10 −2 , and C = −1.52.
To estimate the error of the identified effective radii of the ice crystals, an error propagation of Eq. (7) had be done using the respective error for IWC and α from Hogan et al. (2006).

Detection of low-level stratus clouds
During PS106, frequently low-level stratus clouds (cloud base < 165 m) have been observed. These situations were often associated with a strong attenuation of the lidar beam within the lowest few hundred meters above Polarstern due to the high optical thickness of these clouds. The cloud radar, in turn, has its technical limitation in detecting the lowest part of the boundary layer below 155 m range (165 m a.s.l.). Due to the instrument synergy approach of Cloudnet, this is also the height of the lowest Cloudnet data pixel. Thus, the low-level clouds which occurred during PS106 introduced on the one hand issues to the Cloudnet retrieval due to misinterpretation of attenuated lidar signal as missing signal. On the other hand, since most current statistics of Arctic clouds do not consider clouds in such a low altitude, these clouds tend to be underrepresented in Arctic cloud statistics. To address these issues, we introduce a new Cloudnet classification category called low-level stratus cloud. These clouds were identified by the Polly XT signal-to-noise ratio (SNR; Heese et al., 2010) in the lowest 165 m above sea level. The near-range channels of the Polly XT system have a complete overlap already at 120 m above the instrument , and thus are suitable for the detection of clouds already well below the lowest cloud radar observation height, even though quantitative parameters such as (attenuated) backscatter coefficient from a single elastic backscatter signal cannot yet be determined at these heights.
From visual inspection of the Cloudnet data set, we defined low-level stratus where the SNR exceeded the threshold value of 40. This value was obtained by evaluating signatures of attenuation in the time series of the Cloudnet attenuated backscatter coefficient, increased values of LWP time series, and correlation with the visibility sensor of Polarstern. Since the SNR is not yet range corrected, this threshold crossing at these low altitudes is very likely only due to the occurrence of low-level clouds. The low-level stratus base and top have been derived by simply using the lowest and highest range gate from Polly XT where the SNR exceeded the threshold.

Retrieval of CCN and INP number concentrations
Arctic clouds and their susceptibility to the presence of aerosol are still in the focus of research (Morrison et al., 2012). Based on the measurements of Polly XT , an estimation of cloud condensation nuclei (CCN) and ice-nucleating particle (INP) properties is possible . To do so, profiles of the aerosol backscatter coefficient and depolarization ratio need to be determined. In a second step, these profiles are converted into profiles of the particle extinction coefficient using an appropriate lidar ratio (extinction-to-backscatter ratio).
The direct retrieval of the CCN conversion parameters from the AERONET data (level 2, version 3, inversion products) revealed C 1 = 18.6 cm −3 and exponent d 1 = 0.83 for the range of extinction coefficients from 15 to 300 Mm −1 (500 nm AOD from 0.015 to 0.3 were measured). During the PS106 observations, the aerosol extinction coefficient was mostly around 1-10 Mm −1 in the lower part of the troposphere. The AERONET data for this low range of extinction coefficients indicate that conversion parameters of C 2 = 10 cm −3 , d 2 = 0.9, C 3 = 3.0 cm −3 , and d 3 = 1 would be appropriate. The aerosol in the Arctic is fine-mode dominated and shows Ångstrom exponents (440-870 nm) typically between 0.9 and 1.8 (with an average of 1.5-1.6).

Case studies
Based on near-surface and radiosondes measurements, model data, and satellite observations, Knudsen et al. (2018) gave a detailed description of the synoptic situation during the PASCAL campaign. They defined three periods (cold, warm, and normal period). Three case studies will be presented in the following which are all within the warm pe- riod (WP; 30 May-12 June 2017). These 2 weeks of the WP are characterized by warm and moist air advection from the south crossing Norway and Greenland (Fig. 8).
The three cases between 7 and 9 June 2017 were chosen to demonstrate the potential OCEANET offers. Within these 3 d the near-surface temperature first dropped from −3.5 • C on 7 June 2017 00:00 UTC to −7.5 • C on 8 June 2017 05:00 UTC with an ensuing increase to 1.0 • C on 9 June 2017 22:00 UTC due to warm air advection. These cases were, on the one hand, selected to represent the capabilities of the standard Cloudnet products and of Polly XT and, on the other hand, to illustrate the new products introduced in this study.
4.1.1 Precipitating layered cloud: 9 June 2017 00:00-18:00 UTC ice floe camp An overview of the capability of the OCEANET measurements and its application to analyze cloud and aerosol structures and their interactions will be presented for a complex case that occurred over Polarstern on 9 June 2017 between 00:00 and 18:00 UTC. The radiosonde profiles for this period are shown in Fig. 9 up to a height of 2000 m. The observations of the cloud radar, the lidar, and the MWR are depicted in Fig. 10. The presented day reveals a rather complex situation. Starting at 00:00 UTC, Cloudnet classified a liquid stratocumulus layer between 600 and 900 m height with a cloudtop temperature of −1.5 • C. This layer slowly descended, At around 05:30 UTC, a transformation of the cloud occurred. The LWP increased up to 160 g m −2 and precipitation started, almost reaching the ground (the disdrometer aboard Polarstern showed no precipitation signal, not shown here). Though the LDR showed no increased values, the presence of ice was identified due to detection of enhanced radar reflectivity and vertical velocity and thus a mixed-phase cloud was classified between 05:30 and 06:30 UTC. During this period, the retrievals of the ice and liquid hydrometeor size as proposed in this paper may be influenced by each other. Both retrievals are based on the same quantity, the radar reflectivity, which is characterized by the largest peak in the Doppler spectrum. To tackle this issue, a peak separation of the Doppler spectrum as it is proposed, e.g., by Shupe et al. (2004), Kalesse et al. (2016b would be necessary. This would offer the opportunity to calculate the effective radius of the different hydrometeors species based on their particular reflectivity but is beyond the scope of this paper. At around 06:30 UTC the interpolated temperature of the surrounding radiosonde profiles reached 0 • C, leading to an immediate transition from a mixed-phase to a pure liquid cloud. Therefore, no IWC and no ice effective radius were determined under these conditions. A second transition of the cloud situation during this day is associated with an altocumulus layer which was located above the stratocumulus. Around 09:00 UTC this midlevel cloud layer with a cloud-top temperature of −1 • C occurred Figure 10. OCEANET observations on 9 June 2017 between 00:00 and 18:00 UTC. Panels (a), (b), and (c) show the radar reflectivity factor, Doppler velocity, and linear depolarization ratio. Panels (d) and (e) depict the lidar attenuated backscatter coefficient at 1064 nm and volume depolarization ratio at 532 nm. In panel (f) the liquid water path derived by the microwave radiometer is shown. The dashed vertical lines mark the time of the radiosonde launches on 9 June 2017 (note: the time of the first launch shown in Fig. 9 was before the plotted profiles of the measurements start). at 1900 m height over Polarstern. As this layer increased in geometrical and optical depth, shading effects reduced the cloud-top radiative cooling of the cloud below as already observed in Shupe et al. (2013). This led to a collapse of the EDR in the lower layer at around 12:00 UTC (Fig. 11d) and finally to a dissipation of the cloud. The values for ε in the altocumulus were about the same order of magnitude as for the stratocumulus, indicating that the upper cloud was able to effectively cool to space. Starting at about 14:00 UTC the altocumulus formed a two-layer structure at 1500 and 1200 m, respectively. Due to the shading of the upper layer, the lower one lost its turbulent moment, and the cloud dissipated shortly after. The altocumulus was classified as pure ice cloud, probably due to the fact that the lidar signal was already fully attenuated in the lower layer, which impedes the classification as liquid at an ambient temperature below Figure 11. (a) Cloudnet target classification, (b) ice water content, and (c) liquid water content for 9 June 2017 between 00:00 and 18:00 UTC. Panel (d) shows the time-height profile of EDR calculated from the cloud radar Doppler velocity. The dashed lines mark the time of the radiosonde launches as shown in Fig. 9. 0 • C. At around 11:00 UTC, the temperature exceeded 0 • C, due to which Cloudnet changed its classification from an ice cloud to a liquid cloud. After persisting for about 4 h with rather low EDR, the stratocumulus started to dissipate at around 16:00 UTC. This offered the lidar the opportunity to observe the aerosol structure above Polarstern in the subsequent hours (Fig. 12d).
4.1.2 Aerosol case: 9 June 18:00 UTC-10 June 2017 11:00 UTC ice floe camp Between 9 June 18:00 UTC and 10 June 2017 11:00 UTC, one of the rare cloud-free events of PS106 occurred, and Polly XT observed aerosol layers in the free troposphere (Fig. 12d). The respective radiosonde profiles for this period are shown in Fig. 12a-c. Based on a trajectory analysis of twenty-seven 10 d HYSPLIT back trajectories (Stein et al., 2015), a southern inflow for air masses above the boundary layer is identified for this period. At the 2000 m height level, the trajectories show that these were long-range-transported aerosol layers that originated over continental Europe (Fig. 8) with a high chance of being within the planetary boundary layer at that time. Below 2000 m height, the trajectories in- dicate that pathways mainly crossed the North Sea and the Atlantic Ocean (not shown). The 1064 nm lidar attenuated backscatter coefficient and the 532 nm volume depolarization ratio are shown in Fig. 12d and e. These measurements reveal the existence of three aerosol layers being present over Polarstern on 9 June 2017 around 18:00 UTC. A shallow layer at 500 m height stayed roughly at the same altitude as long it was observed by Polly XT . A second one was visible between 1000 and 1500 m height ascending to 2500 m at 00:00 UTC on 10 June 2017. At 07:00 UTC on 10 June 2017, a liquid cloud formed at the top of this layer. A third aerosol layer with a liquid cloud embedded at 2300 m height, being rather constant in altitude, was present between 19:00 and 21:00 UTC on 9 June 2017.
In Fig. 13, a detailed analysis of the aerosol optical properties as derived from the lidar measurements from the time period of 00:00-02:20 UTC is presented. During this period, two layers were detected and are visible in the profiles of the backscatter coefficient at all three wavelengths (Fig. 13a). The rather strong wavelength dependence of the backscatter coefficient, as shown by the high Ångstrom exponents ( Fig. 13b) in both layers, indicates the presence of small aerosol particles. A back-trajectory analysis and the values for the Ångstrom exponent indicate that air masses transporting polluted aerosol from continental Europe are most probably the source for the upper aerosol layer. The lower aerosol layer on the other hand is most likely a mixture of down-mixed continental and upward-mixed marine aerosol.
Based on the aerosol optical properties retrieved by Polly XT , an estimation of the CCNC was done for all three combinations of conversion factor and extinction exponent as mentioned in Sect. 3.4. We chose the second combination (C 2 = 10 cm −3 and d 2 = 0.9) to illustrate the results in Fig. 13d. The mean value of CCNC for the upper aerosol layer in this case was found to be n CCN,2 = 180 cm −3 . For the two other combinations, we found n CCN,1 = 260 cm −3 and n CCN,3 = 75 cm −3 . For the lower aerosol layer, we found n CCN,2 = 70 cm −3 . For the two other combinations, we found n CCN,1 = 110 cm −3 and n CCN,3 = 25 cm −3 . The uncertainty in this method is up to 200 % (dashed lines in Fig. 13d). In addition, the INPC was calculated for a fixed temperature of −15 • C. The results are shown in Fig. 13e. The INPC for the lower layer was found to be around 0.6 L −1 for this temperature. In the upper layer, n INP went up to values slightly above 1 L −1 . These calculations have an error of a factor of 3 (dashed lines in Fig. 13e) but nevertheless provide a guideline about the conditions of the cloud-relevant aerosol properties on the discussed day. The corresponding thermodynamic profiles from the radiosondes launched during this time period are shown in Fig. 15. In Fig. 16a, the Cloudnet target classification, together with the IWC (b), ice effective radius (c), and ice water path (d) are shown. The ice water path is derived as the integral of the IWC for each profile. This period is chosen to illustrate the low-level stratus cloud detection algorithm and the retrieval of the ice effective radius. At the beginning of the addressed time period, a lowlevel mixed-phase stratiform cloud reaching up to a height of 500 m is present. This layer is visible in both the cloud radar reflectivity and in the lidar attenuated backscatter data. Due to its high optical thickness, this cloud led almost continuously until 23:30 UTC to an attenuation of the lidar signal already close to cloud base (Fig. 14). Only occasionally, backscattered lidar signals from aerosols above the cloud were detected. During this period, the liquid water path var- ied between values of around 0 g m −2 for moments when the lidar was able to detect signal from above the cloud and values up to 100 g m −2 associated with periods when the lidar signal was attenuated already close to cloud base.
In Fig. 17, the derived low-level stratus classification mask (below 165 m height) combined with a simplified Cloudnet target classification mask (above 165 m height) for this period is shown. Red areas depict detected low-level clouds. Blue and green data points indicate clear sky and aerosols, respectively. Though Polly XT detected low-level stratus almost continuously during the case study, this affected the lidar signal most severely during the abovementioned period. After 23:30 UTC, the LWP showed values of around 10 g m −2 , and the cloud lost most of its optical thickness so that the lidar was able to penetrate through the cloud. Well above the low-level stratus, some cirrus clouds formed around 22:30 UTC above 6000 m height. These transformed into a cirrostratus layer at 00:00 UTC, which was present between 6000 and 10 000 m height. Until around 04:30 UTC, this cloud is classified as a pure ice cloud, characterized by LDR values of up to −15 dB and a constantly downward-directed vertical velocity, with a tiny patch of detected liquid at around 03:30 UTC at 6100 m height.
At 04:30 UTC, the cirrostratus dissipated, and another layer started to pass over Polarstern at around 05:00 UTC. This layer with coexisting liquid droplets and ice crystals extended from 5000 up to 7000 m height. While the cloud radar reflectivity factor was higher in this layer compared to the first one, the cloud radar LDR decreased to values of below −20 dB. On top of this layer a supercooled liquid layer was detected by the lidar between 05:00 and 06:00 UTC, characterized by a high attenuated backscatter coefficient and low values of linear depolarization ratio. Additionally, some regions with a high linear depolarization ratio were detected by Polly XT inside the cirrus after 06:30 UTC, probably associated with a mixture of supercooled droplets and ice crystals.
The IWC of the cirrostratus was in a range from 10 −4 to 10 −6 kg m −3 with lowest values at cloud top and highest values at cloud base. The ice effective radius ranged from 30 to 55 µm, and its distribution follows the same pattern as the one of the IWC, as can be expected because both follow a similar reflectivity-temperature relationship.

Cloud statistics
In Fig. 18, an overview of the statistical distribution of the cloud occurrence during PS106 is given. In Fig. 18a, daily statistics of the vertical distribution of low-level stratus are shown. In addition, the frequency of occurrence of this cloud type for each day is illustrated in Fig. 18b. Low-level stratus was detected during a significant period of time on almost every day. The highest frequency of occurrence was observed while the Polarstern was surrounded by sea ice. Rather low values occurred while Polarstern was in the vicinity of Svalbard. In order to assess our retrieval of the low-level clouds, we plotted in comparison the frequency of occurrence of vertical visibility below 1 km. Statistics of the cloud type occurrence are shown in Fig. 18c. The daily frequency of occurrence as well as the total distribution for the complete campaign of low-level stratus clouds (purple), liquid clouds (orange), ice clouds (light blue), mixed-phase clouds (green), multi-layer clouds (dark blue), and cloud-free situations (yellow) is shown. In addition, an analysis of the co-occurrence of low-level stratus and other cloud types was performed and is shown in the rightmost column of Fig. 18c. The rate of coexistence of the respective cloud type together with low-level stratus is indicated by a slightly varied color code.
In total, during 11 % of the time cloud-free conditions were detected by Cloudnet during PS106. The two most prominent cloud types were multi-layer and mixed-phase clouds with an occurrence frequency of 38.5 % and 36 % of the observational time, respectively. Pure ice clouds were present for about 8 % and pure liquid clouds for about 4.5 % of the time, respectively. Single events of the new Cloudnet class low-level stratus cloud were detected during 2.5 % of the time during the 2 month campaign. In addition, 27 % of the observed liquid clouds and 48 % of the ice clouds occurred simultaneously with low-level stratus. Mixed-phase and multi-layer clouds were detected together with low-level stratus clouds during 24 % and 27 % of their respective observational time.
In contrast to Nomokonova et al. (2019), which provided a statistical analysis of the cloud occurrence over Ny Ålesund, Svalbard, for the period between June 2016 and July 2017, we found a higher frequency of single-layer mixedphase clouds at the expense of cloud-free and single-layer liquid clouds when comparing the period of PS106. This may be due to a difference in turbulence as well as a change in the cloud microphysics at locations surrounded by sea ice or open ocean (Young et al., 2016).

Summary and conclusions
A 2-month campaign of RV Polarstern, including an extensive suite of ground-based remote-sensing instruments of the OCEANET platform, has been conducted north and northeast of Svalbard in the Arctic summer of 2017. This study described in detail the deployed instrumentation and the applied processing schemes. Only a few campaigns with comparable equipment have been performed in recent years at these latitudes, e.g., ASCOS, which took place from 2 August to 9 September 2008 , and ASCE in the Arctic summer and early autumn of 2014 . A new feature of PS106 was the deployment of a motion-stabilized vertically pointing 35 GHz cloud radar during and the correction of the Doppler velocity subsequent to the cruise as specified in Sect. 3.1.
For an automatic, seamless analysis of cloud properties from the measured remote sensing time series, the Cloudnet algorithm was utilized. In doing so, new products were developed and applied to the remote sensing data set from PS106. This was done in order to enable the continuous characterization of cloud turbulence by means of EDR and to provide mass concentration and effective radius of ice crystals and liquid water droplets as future input for radiative transfer simulations. Though being well established, applying the Cloudnet algorithm to data from a remote sensing supersite aboard a research vessel in the Arctic reveals new challenges. The movement of the ship has a significant effect on the measured vertical velocity of the cloud radar. To tackle this issue, the cloud radar was mounted on a stabilization platform to guarantee its vertical pointing. The vertical velocity data set was corrected for the heave rate of the ship in a postprocessing procedure subsequent to the cruise.
The motion stabilization was evaluated by means of a small single-board computer mounted on the cloud radar rack. The IMU of the minicomputer measured the residual of the pitch and roll movement after stabilization. We found good stabilization during ice breaking conditions with a leveling precision of ±0.5 • . During rough sea, however, the displacement from zenith was larger, up to ±1 • . Under the strong wave conditions during these time periods it needs also to be considered that the IMU of the orientation sensor used for the cloud radar is based on so-called MEMS (micro-electro-mechanical systems). Such devices are based on spring-mounted capacitor plates, and thus measured pitch and roll angles are affected by translational motions like engine vibrations, etc. As these effects were not investigated in the frame of our study, we conclude that the actual verticalpointing uncertainty range, especially on the open sea, was lower than the one reported by the MEMS sensors, i.e., better than ±1 • . Figure 18. (a) Daily height distribution of the detected low-level stratus during PS106 up to 165 m. (b) Daily fraction of low-level stratus occurrence derived by Polly XT measurements (purple) in comparison to horizontal visibility (HV) below 1 km (grey). (c) Cloud type statistics including low-level stratus clouds during PS106, determined by Cloudnet. Purple indicates the fraction when low-level stratus was determined, while orange indicates liquid clouds, light blue ice clouds, green mixed-phase clouds, dark blue multi-layer clouds, and yellow cloud-free periods. Each column except the last two represents a day of the campaign. The penultimate column represents the total distribution of the different cloud types. The last column distinguishes between the respective cloud type without low-level stratus detected (same color as in the other column) and with an additional low-level stratus detected below (slightly varied color). At the bottom, a rough localization of Polarstern is annotated -green: North Sea (N.S.); light blue: ice transit (Ice t.); dark blue: ice floe; camp; dark green: Svalbard region.
Using the corrected vertical velocities from the cloud radar, the eddy dissipation rates were calculated and evaluated against in-cloud turbulence measurements done by means of a tethered balloon. This intercomparison revealed good agreement between both approaches where the values were within the estimated uncertainty and the expected difference due to the spacial distance. Nevertheless, the tethered balloon approach seems to be the more reliable one, given the smaller standard deviation.
Based on published retrievals of visible extinction coefficient and ice water content, a new approach to derive the effective radius of the ice crystals was introduced. The associated uncertainties, estimated by error propagation, of the ice crystal effective radii are presented in Fig. 16d. On average the uncertainty is about 50 % of the size of the radii themselves, which reflects the strong influence of uncertainties in the underlying observational data on the retrieval. Given the challenges in estimating the effective radius of ice crystals on a continuous basis on the one hand and the necessity of having such values, e.g., for radiative transfer calculations, on the other hand, we consider this estimate to be still in a reasonable range. In Fig. 19, the histogram of the effective ra- Figure 19. Histogram of the ice particle effective radius for PS106. Integration over x yields 100 %. dius for full PS106 is shown. Values range from 20 to 60 µm, with a peak at around 50 µm. This is consistent with other studies of ice effective radius (e.g., Blanchard et al., 2017).
This study revealed, in addition, the relevance of the lowest detection limit of remote-sensing instruments on the representativity of Arctic-cloud statistics. Cloudnet is configured to have its lowest range gate at the lowest detection altitude of the cloud radar, which was 165 m above the ocean surface for PS106. Lower-level cloud layers are thus not identifiable within Cloudnet. In our study, lower cloud structures were identified using the SNR measured by the lidar Polly XT . This ability has been used to study the occurrence of low-level stratus clouds below the first Cloudnet range gate. Sotiropoulou et al. (2014) used a combination of cloud radar and ceilometer measurements to study stratiform Arctic clouds and found that most stable, surface-coupled clouds have a cloud base below 200 m. Yet, so far such clouds have not been considered in many Arctic cloud climatologies derived by remote-sensing instruments. Liu et al. (2012) for example defined low-level clouds as those between 0 and 2000 m, with 960 m above the ground being the height where surface contamination effects on Cloudsat become insignificant, and using a vertical resolution of 240 m. Shupe (2011) summarized cloud statistics from several multiyear data sets derived from ground-based remote sensing observations for different sites in the Arctic. The author specified a height dependence of cloud occurrence down to 300 m by using a combination of lidar and radar. Below 300 m, however, the author provided information about cloud occurrence but without any further specification of the cloud base. Even airborne remote-sensing instruments suffer from the strong ground clutter and thus struggle to deliver information about cloud occurrence below 150 m height above the surface . The autonomous buoys of the IAOOS network equipped with a microlidar observed clouds in 2014 and 2015 in the Arctic with a base height below 500 m during 60 % of their observational time (Mariage et al., 2017).
Our study shows that a higher vertical resolution and reliable signal from very low altitudes are required to characterize the lowest-level cloud layer, which occurs between approximately 50 and 165 m above the ground. Such clouds stay undetected for ground-based in situ sensors (because they are too high) as well as for most automatized groundbased remote-sensing instruments (because they are too low). Future radiative transfer studies should show what the effect of the lowest-level clouds, which occurred during 25 % of the observation time, is on the radiation budget of the region where PS106 was performed.
Future work will confront the observed cloud macro-and microphysical properties as well as the EDR with highresolution model simulations along the PS106 track that have been carried out in the framework of (AC) 3 . The herein introduced remote sensing techniques will also be applied to the data set of the currently ongoing 1-year polar ice drift of RV Polarstern during the MOSAiC project (Shupe et al., 2018), thus providing an unprecedented data set of Arctic aerosol and mixed-phase clouds. This data set will substantially contribute to our understanding of the role of clouds in the current warming of the Arctic climate system.
Author contributions. PS, RE, and JB designed the OCEANET measurements, and PS, RE, JB, MR, and HJG prepared the instruments. HJG, RE, MR, and CBV operated the OCEANET instruments aboard Polarstern. PS and HJG processed the data with Cloudnet and developed the low-level stratus detection algorithm. HJG performed the heave correction and developed the EDR and ice crystal effective radius retrieval with supervision from PS, RE and JB. HB, YZ, and AA provided lidar-based data processing for retrieval of aerosol optical and microphysical properties. AM provided support in the data analysis and manuscript preparation. HJG drafted the paper with contributions from all coauthors.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Arctic mixed-phase clouds as studied during the ACLOUD/PASCAL campaigns in the framework of (AC) 3 (ACP/AMT/ESSD inter-journal SI)". It is not associated with a conference. 172) and by the European Commission's Horizon 2020 research and innovation programme (grant agreement no. 654109).
Review statement. This paper was edited by Joachim Curtius and reviewed by two anonymous referees.