Journal topic
Atmos. Chem. Phys., 20, 3459–3481, 2020
https://doi.org/10.5194/acp-20-3459-2020
Atmos. Chem. Phys., 20, 3459–3481, 2020
https://doi.org/10.5194/acp-20-3459-2020

Research article 24 Mar 2020

Research article | 24 Mar 2020

# Low-level mixed-phase clouds in a complex Arctic environment

Low-level mixed-phase clouds in a complex Arctic environment
Rosa Gierens1, Stefan Kneifel1, Matthew D. Shupe2,3, Kerstin Ebell1, Marion Maturilli4, and Ulrich Löhnert1 Rosa Gierens et al.
• 1Institute for Geophysics and Meteorology, University of Cologne, Cologne, Germany
• 2Cooperative Institute for the Research in Environmental Sciences, University of Colorado Boulder, Boulder, Colorado, USA
• 3Earth System Research Laboratory, National Oceanic and Atmospheric Administration, Boulder, Colorado, USA
• 4Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Potsdam, Germany

Correspondence: Rosa Gierens (rgierens@uni-koeln.de)

Abstract

Low-level mixed-phase clouds (MPCs) are common in the Arctic. Both local and large-scale phenomena influence the properties and lifetime of MPCs. Arctic fjords are characterized by complex terrain and large variations in surface properties. Yet, not many studies have investigated the impact of local boundary layer dynamics and their relative importance on MPCs in the fjord environment. In this work, we used a combination of ground-based remote sensing instruments, surface meteorological observations, radiosoundings, and reanalysis data to study persistent low-level MPCs at Ny-Ålesund, Svalbard, for a 2.5-year period. Methods to identify the cloud regime, surface coupling, and regional and local wind patterns were developed. We found that persistent low-level MPCs were most common with westerly winds, and the westerly clouds had a higher mean liquid (42 g m−2) and ice water path (16 g m−2) compared to those with easterly winds. The increased height and rarity of persistent MPCs with easterly free-tropospheric winds suggest the island and its orography have an influence on the studied clouds. Seasonal variation in the liquid water path was found to be minimal, although the occurrence of persistent MPCs, their height, and their ice water path all showed notable seasonal dependency. Most of the studied MPCs were decoupled from the surface (63 %–82 % of the time). The coupled clouds had 41 % higher liquid water path than the fully decoupled ones. Local winds in the fjord were related to the frequency of surface coupling, and we propose that katabatic winds from the glaciers in the vicinity of the station may cause clouds to decouple. We concluded that while the regional to large-scale wind direction was important for the persistent MPC occurrence and properties, the local-scale phenomena (local wind patterns in the fjord and surface coupling) also had an influence. Moreover, this suggests that local boundary layer processes should be described in models in order to present low-level MPC properties accurately.

1 Introduction

The Arctic is warming more rapidly than any other area on Earth due to climate change . It is well established that clouds strongly impact the surface energy budget in the Arctic , but feedback processes that include clouds are not well characterized . Particularly low-level mixed-phase clouds are important for the warming of near-surface air . The multitude of microphysical and dynamical processes within the cloud and the interactions with local and large-scale processes make these mixed-phase clouds difficult to represent in numerical models . Improvements in the process-level understanding are still required to improve the description of low-level mixed-phase clouds in climate models .

Previous studies have shown the prevalence of mixed-phase clouds (MPCs) across the Arctic . MPCs occur in every season, with the highest occurrence in autumn and in the lowest 1 km above the surface, and they can persist from hours to days . The persistent low-level MPCs have a typical structure that consists of one or more supercooled liquid layers embedded in a deeper layer of ice, where liquid is usually found at cloud top, and the ice precipitating from the cloud may sublimate before reaching the ground (Morrison et al.2012, and references therein). Several studies have shown an increase in cloud ice to coincide with increase in cloud liquid, suggesting that ice production is linked to the liquid water in the cloud . The amount of liquid and ice and the partitioning between the condensed phases (i.e., phase partitioning) are important parameters due to their key role in determining the clouds' radiative effect .

A variety of environmental conditions can effect cloud micro- and macro-physical properties. According to simulations over different surface types, changes in surface properties lead to changes in the thermodynamic structure of the atmospheric boundary layer, the extent of dynamical coupling of the cloud to the surface, and the microphysical properties of the MPC . Also, observational evidence on the connection between changes in surface conditions and MPC occurrence as well as thermodynamic structure and droplet number concentration has been found. discovered in a detailed case study that for the MPC in question phase partitioning was affected by the coupling of the cloud to the surface, large-scale advection of different air masses as well as local-scale dynamics. Conversely, did not find differences in cloud water properties between coupled and decoupled clouds. show that at Ross Island, Antarctica, orographic lifting of marine air is likely causing thick MPCs with high ice water content. Changes in aerosol population, especially ice-nucleating particles (INPs), have been found to modulate the ice formation rate . To complicate matters further, the cloud also modifies the boundary layer where it resides by modifying radiative fluxes, generating turbulence (due to cloud-top cooling), and vertically redistributing moisture .

While being common in the entire Arctic, MPCs are most frequently observed in the area around Svalbard and the Norwegian and Greenland seas . report 1 year of ground-based remote sensing observations of clouds at Ny-Ålesund, Svalbard, and find that 20 % of the time single-layer MPCs (defined as single-layer clouds with ice and liquid occurring at any height of the cloud) were present, with the highest frequency in autumn and in late spring–early summer. Svalbard lies in a region where intrusions of warm and moist air from lower latitudes are common , and differences in air mass properties have been associated with differences in ice and liquid water content and particle number concentration in MPCs . Locally, the archipelago exhibits large variations in surface properties (glaciers, seasonal sea ice, and snow cover) as well as orography. There are fewer MPCs over the islands than over the surrounding sea during winter and spring, while during summer and autumn the differences are small , indicating that the islands modify the local boundary layer and the associated clouds. How the orography influences the low MPCs in more detail is difficult to study using spaceborne radars due to the rather big blind zone and considerably large footprint. Aircraft and ground-based remote sensing, together with modeling studies, are better suited for answering this question.

In this paper we investigate persistent low-level mixed-phase clouds (P-MPCs) observed above Ny-Ålesund on the west coast of Svalbard. Mountainous coastlines are common at Svalbard, Greenland, and elsewhere in the Arctic . Ny-Ålesund is an excellent site to study low-level MPCs in such complex environments. The time period considered is June 2016–October 2018, when a cloud radar of the University of Cologne was operating at the French–German Arctic Research Base AWIPEV as part of the project Arctic Amplification: Climate Relevant Atmospheric and Surface Processes, and Feedback Mechanisms (AC)3. A combination of ground-based remote sensing instruments, surface meteorological observations, radiosoundings, and reanalysis data was used to identify and characterize the P-MPCs, to describe the extent of surface coupling, and to evaluate these in the context of wind direction in the area around the station. In addition to providing a description of micro- and macro-physical properties of P-MPCs and their seasonal variation, we aim to identify some of the impacts the coastal location and the mountains have on the observed P-MPCs as well as determine the relevance of surface coupling for cloud properties at the site. In Sect. 2 the measurement site, the instrumentation, and the data products used are introduced, followed by the description of the methodology developed to identify persistent low-level MPCs (Sect. 3.1), coupling of the cloud to the surface (Sect. 3.2), and the approach to describing regional and local wind conditions (Sect. 3.3 and 3.4). The Results and discussion section describes the occurrence of P-MPCs (Sect. 4.1) and their average properties as well as variation under different seasons and dynamical conditions. The relationship between P-MPC and the regional wind direction (Sect. 4.2), different seasons (Sect. 4.3), surface coupling (Sect. 4.4), and local wind conditions (Sect. 4.5) are considered. The results are discussed in the context of atmospheric temperature and humidity and considering previous studies at Ny-Ålesund and other Arctic sites. In the end the main aspects are summarized followed by conclusions in Sect. 5.

2 Observations

## 2.1 Measurement site

The measurements were carried out at the French–German Arctic Research Base AWIPEV in Ny-Ålesund (78.9 N, 11.9 E), located on the west coast of Svalbard, on the south side of Kongsfjorden (Fig. 1). The area is mountainous, featuring seasonal snow cover, a typical tundra system, and glaciers. In the period investigated the sea remained ice-free throughout the year. The local boundary layer is known to be strongly affected by the orography and is often quite shallow with an average mixing layer height below 700 m . Surface layer temperature inversions are common, especially in winter . The mountains reach up to 800 m and strongly influence the wind around Ny-Ålesund. In the free troposphere westerly winds prevail. The wind conditions are described more in detail in Sect. 3.4. Clouds have been found to occur above Ny-Ålesund 60 %–80 % of the time . Clouds generally occur more frequently in summer and autumn and are less common in spring, although the inter-annual variability is large.

Figure 1Topography map for Svalbard (a) and a detailed illustration of the Kongsfjorden area (b). The black star indicates the location of Ny-Ålesund, where the measurements are taken. The domain covered by the circulation weather type (see Sect. 3.3) is shown by the blue rectangle. Topography data by  (a) and the  (b).

## 2.2 Measurements and data products

Most of the measurements and cloud and thermodynamic parameter retrievals utilized were described in detail by and references therein. Here, the most important aspects are summarized, together with additional data products used. A summary of the instrumentation, their specifications, and derived parameters is given in Table 1.

Table 1The instruments used, the most relevant specifications of each measurement, and overview of derived parameters are shown. If the vertical (Δz) or temporal resolution (Δt) is changed from that measured by the instrument, the resolution used in the analysis is given in the last column.

### 2.2.1 Instrumentation

To complement the remote sensing observations, we make use of soundings and standard meteorological parameters measured at the surface. In Ny-Ålesund radiosondes are launched routinely every day at 11:00 UTC, and more often during campaigns . From the surface measurements we utilized temperature, pressure, and wind speed and direction data (technical details in Table 1). The instruments for surface meteorology are continuously maintained by the AWIPEV staff, and all data are quality controlled .

### 2.2.2 Cloudnet

The Cloudnet algorithm combines radar, radiometer, and ceilometer measurements with thermodynamic profiles from a numerical weather prediction (NWP) model to provide the best estimates of cloud properties . The observational data, described in the previous section, are homogenized to a common resolution of 30 s in time and 20 m in the vertical. In the Ny-Ålesund data set, the Global Data Assimilation System 1 (GDAS1, more info at https://www.ready.noaa.gov/gdas1.php, last access: 28 February 2019) was used as the NWP model until the end of January 2017, after which it was replaced by the operational version of the ICON (ICOsahedral Non-hydrostatic) NWP model .

In our work we rely on the target classification product that classifies objects detected in the atmosphere as aerosols, insects, or different types of hydrometeors (cloud droplets, drizzle, rain, ice, melting ice; see Fig. 2 for an example). Radar reflectivity (Ze) and ceilometer β-profiles are used to detect the presence and boundaries of clouds. Cloud phase is distinguished on the basis of Ze, β, temperature (T), and wet bulb temperature; in addition, Doppler velocity from radar is used to position the melting layer. No differentiation is made between ice in a cloud and precipitating ice. While applying this widely accepted methodology, for our study there are two important limitations. Firstly, the detection of liquid within a MPC is based on β, such that if cloud top is not found within 300 m from the height where the ceilometer signal is extinguished, all cloudy bins above this height are classified as ice. Secondly, no method to distinguish supercooled drizzle from ice is available yet .

Figure 2Example of the Cloudnet target classification product from 29 May, 18:00 UTC to 30 May 2018, 12:00 UTC. For the P-MPC, indicated by the gray dashed box, the time series of coupling is also shown. This case was classified as coupled.

### 2.2.3 Derived properties

The amount of liquid and ice in the cloud, and their ratio, is one of the most important properties of MPCs. In addition, humidity supply is a key requirement for cloud formation and continuation. Liquid water path (LWP) and integrated water vapor (IWV) were retrieved from the zenith-pointing observations of the MWR using statistical multivariate linear regression . Coefficients for the retrieval were based on sounding data; more details about the retrieval and corrections applied are given by . Previous studies have found the accuracy of the method to be 20–25 g m−2 .

Ice water content (IWC) was calculated using the ZeT relationship from , where temperature was taken from the same model as used for Cloudnet. The uncertainty of the retrieval is estimated to be −33 % to +50 % for temperatures above −20C. Ice water path (IWP) for P-MPCs was calculated by integrating IWC from the surface to cloud top. Furthermore, the LWP was averaged to 30 s to match the temporal resolution of IWP.

In order to calculate the potential temperature (θ) profile based on the temperature profile retrieved from the MWR elevation scans, an estimate of the pressure profile is required. For this we took the measured surface pressure and used the barometric height formula to estimate pressure at each height. The resulting θ-profiles were compared with the profiles from radiosondes in the period June 2016–October 2018 (not shown). A slight cold bias is present (<0.4 K). The RMSE increases with altitude, but in the lowest 2.5 km the RMSE is still below 1.8 K. For cloud-top temperature, the temperature retrieved from the MWR elevation scan was linearly interpolated between the retrieval levels to cloud-top height.

3 Methods

## 3.1 Identification of persistent low-level mixed-phase clouds

To identify P-MPCs, each profile was evaluated individually to detect low pure-liquid and liquid-topped mixed-phase cloud layers, after which the persistency of the liquid layer was considered. Using Cloudnet target classification, the first step was to identify different cloud layers in each profile. Here a cloud layer refers to a continuous (gaps of fewer than four height bins, corresponding to 80 m, were omitted) layer of cloud droplets and/or ice. Each layer in the profile was classified as ice only, liquid only, or mixed phase. To distinguish between low stratiform and deep multilayered mixed-phase clouds, only profiles with a single liquid layer and the liquid layer being close to cloud top were considered. In practice, the detected upper boundary of the liquid layer was required to be in the uppermost 20 % of the cloud layer. The requirement for liquid being exactly at cloud top was relaxed since the ceilometer signal cannot necessarily penetrate the entire depth of the liquid layer. These criteria (single liquid layer, liquid close to cloud top) were very effective in selecting the desired low-level mixed-phase cloud regime. However, some mid-level clouds also fulfilled the criteria, and therefore we limited cloud-top height to be below 2.5 km. For the remaining profiles, which all contain either liquid-only or liquid-topped mixed-phase clouds, the persistence of the liquid layer was evaluated. We only included clouds where the liquid layer existed for a minimum of 1 h, with gaps ≤5 min. Since the focus of this study is on mixed-phase clouds, we further excluded clouds where no ice was detected. Note that continuous presence of ice was not required; only the cloud liquid had to persist in time. The result is a data set with clouds below 2.5 km where liquid is located at cloud top and persists at least 1 h, and at some point in time ice is associated with the liquid layer. Note that time periods where another cloud layer is found above the P-MPC are not excluded. Figure 2 shows an example of the identified persistent MPC as well as another mixed-phase cloud. Despite the strict criteria, such clouds were present 23 % of the observational time.

In addition to identifying the time periods with P-MPCs present, the Cloudnet data were used to determine the base of the liquid layer and the cloud-top height. A P-MPC case was defined as the time from the beginning to the end of the identified persistent liquid layer. Furthermore, we consider the layer from liquid base to cloud top as the cloud and everything below liquid base to be precipitation. This definition was chosen because the liquid base is well defined from the ceilometer observations. Considering the focus on a persistent liquid layer identified by vertically pointing measurements, the cases included implicitly require either very low wind speeds or a larger cloud field being advected over the site. When another cloud is detected above the P-MPC, the possibility that it contains undetected liquid cannot be excluded, and in these cases the measured LWP cannot be unambiguously attributed to the liquid layer of the P-MPC. Hence, those time periods were flagged to be removed in any analysis of the cloud’s liquid content. Unfortunately, we cannot make the assumption that upper cloud layers would not impact the liquid content of a P-MPC . The presented LWP distributions are therefore only representative for single-layer cases. Furthermore, all columns with liquid precipitation or drizzle were excluded, leading to a loss of data mainly in the summer months. While this is somewhat unavoidable (e.g., when the MWR measurements suffer from a wet radome), it leads to the exclusion of rather warm precipitating P-MPCs from the analysis.

## 3.2 Detecting surface coupling

### 3.2.1 Defining coupling with radiosonde profiles

The thermodynamic coupling of the P-MPC to the surface was determined based on the θ-profile. A quasi-constant profile was taken to indicate a well-mixed layer, while an inversion denotes decoupling between different layers. For the sounding profiles, we simplified the methodology of . The cumulative mean of θ from the liquid layer base height downward is compared to θ at each level below the cloud. If this difference exceeds 0.5 K, the cloud is considered decoupled. Figure 3a and b show two example cases, one for a coupled and one for a decoupled cloud, respectively. Both profiles demonstrate a structure typical for stratiform Arctic MPC: a temperature inversion at cloud top, below which a well-mixed layer is identifiable. In the case of the coupled P-MPC (Fig. 3a), the well-mixed layer extends to the surface. For the decoupled P-MPC (Fig. 3b) the well-mixed layer extends 200 m below the liquid layer base, below which several weaker temperature inversions and a generally stable stratification can be identified.

Figure 3Examples of θ-profiles from sounding and MWR, as well as surface observations: a coupled cloud on 24 October 2017 11:55–12:05 UTC (a) and a decoupled cloud on 1 February 2018 16:47–16:56 UTC (b). The blue shaded area indicates the cloud layer, where cloud base and top are determined as the median values of the Cloudnet-based cloud base and top for the duration of the sounding. The gray dashed line indicates the decoupling height defined from the sounding θ-profile. Comparison of potential temperature profiles from sounding and retrieved from MWR measurements when P-MPCs were present, with height normalized with respect to the liquid layer (c). Comparison of the diagnosed coupling with the new method based on MWR and surface observations and based on sounding profiles (d).

### 3.2.2 New continuous method

To continuously evaluate the coupling of the P-MPC to the surface, we developed a new method based on surface observations and the potential temperature profiles retrieved from the MWR, which are available more frequently, i.e., every 15–20 min, compared to radiosonde data (Sect. 2.2.3). At each time when a MWR θ-profile was available, the cloud was classified as either coupled or decoupled based on a two-step algorithm. First, the stability of the surface layer was evaluated using the measurements of the meteorological station. The premise of this criteria is that if the surface layer is stably stratified, the cloud must be decoupled from the surface as there exists a stable layer between the surface and cloud base. The θ-profile is used as a proxy for stability. If the gradient of the 30 min mean θ between 2 and 10 m was positive (e.g., an inversion was present between 2 and 10 m), the surface layer was considered stably stratified, and therefore the cloud was decoupled. If this was not the case, the second criteria based on the MWR θ-profile was used. We calculate the difference in potential temperature (Δθ) between the surface level and at the height halfway to the liquid base height. If Δθ is below the threshold of 0.5 K, the cloud at this instance was considered coupled, and otherwise it was considered decoupled. The reason for using the height equaling half of the liquid layer base height can be understood by comparing the θ-profiles from sounding and MWR in Fig. 3a and b. While the general shape of the profile can be retrieved from the MWR measurements, it is not possible to resolve sharp inversions or detailed structures of the profile. Yet temperature inversions are very common at the top of P-MPCs. The comparison of MWR profiles with all available soundings when a P-MPC was present (Fig. 3c) shows that the accuracy of the retrieved potential temperature is reduced in the vicinity of the liquid layer top and that the influence extends to below the liquid layer base. At 0.5× liquid base height the impact of cloud-top inversion is smaller than at liquid base and the RMSE is below 1 K, which is why we chose this height to determine the stability of the subcloud layer. Note that it should not be inferred that the method can only detect decoupling occurring in the lowest half of the subcloud layer. When decoupling occurs above the layer explicitly included, it is common that the lower half of the subcloud layer is at least partly stably stratified, as can also be seen in the example of Fig. 3b, prompting a correct decoupling classification.

As the final step, the individual profiles were considered together to define the degree of coupling of each observed P-MPC case. For each detected cloud event, the number of coupled and decoupled profiles were counted. If all profiles were decoupled, the P-MPC was considered fully decoupled. When more than 50 % of the profiles were found to be decoupled, the P-MPC was defined as predominantly decoupled. The rest were considered coupled.

### 3.2.3 Comparison of methods

The performance of the new method for estimating the coupling for each individual profile was evaluated using the soundings as a reference. We restricted the soundings to cases for which the cloud was present from the launch time until the sonde passed a height of 2.5 km (maximum cloud-top height considered). Those soundings were compared to the MWR profile closest (but not more than 20 min away) to the radiosonde launch time. The sounding-based diagnosis found 31 % of the evaluated P-MPCs to be coupled and 69 % to be decoupled, compared to 18 % and 82 %, respectively, for the newly developed method for the corresponding clouds (Fig. 3d). This suggests a tendency in our method towards decoupling. However, the sounding profiles may miss very shallow surface-based inversions. For 24 % of the profiles considered to be coupled based on the radiosondes, the 2 and 10 m temperatures indicate a surface inversion. Classifying these clouds as decoupled instead changes the ratio of coupled and decoupled P-MPCs from the sounding data set to 23 % and 77 %, which is closer to that found with the new method. The main disadvantage of our method is that the temperature profiles retrieved from the MWR measurements do not provide a detailed profile, but rather the general shape of the profile, and so the developed method occasionally fails. Furthermore, the 10 m layer considered for surface stability is rather shallow, and intermittent coupling could occur regardless of the thermodynamic profile structure.

## 3.3 Circulation weather type

Since the local wind direction in the lower troposphere above Ny-Ålesund is heavily influenced by the orography , the circulation weather type based on was applied in order to evaluate cloud properties in the context of the regional wind field. Using 850 hPa geopotential height and shear vorticity from ERA-Interim, the flow at each time (00:00, 06:00, 12:00, and 18:00 UTC) was classified as either W, NW, N, NE, E, SE, S, SW, cyclonic, or anticyclonic. A total of 16 grid points centered around Ny-Ålesund were used, so that the area covered is approximately 300 km in the meridional and 100 km in the zonal direction (77.5–80.5 N, 9.75–14.25 E; see Fig. 1a). The approach aids in assessing whether the observed clouds were advected to the site from the open sea or over the island as well as the proximity of high- and low-pressure systems.

## 3.4 Local wind conditions

The channeling of the free-tropospheric wind along the fjord axis is a typical feature of an Arctic fjord (Svendsen et al.2002; Esau and Repina2012, and references therein). Previous work has found the feature to also be prominent at Kongsfjorden . It is well documented that despite the dominating westerly free-tropospheric wind direction, in Kongsfjorden the near-surface wind tends to blow southeasterly along the fjord axis . This is usually attributed to katabatic forcing of the Kongsvegen glacier about 15 km east-southeast from Ny-Ålesund (Fig. 1), although argued that for typical synoptic conditions the land–sea breeze circulation would be the dominant driver. The secondary mode in surface wind is from the northwest, from the sea towards the island's interior. According to , the northwesterly surface winds are associated with cold air advection that is related to passing low-pressure systems. find this wind direction to be pronounced in June and July, which they associate with sea breeze and the melting of sea ice. In addition, at Ny-Ålesund weak southwesterly surface winds are observed, caused by katabatic flow from Zeppelin mountain range and the Brøggerbreen glacier south of Ny-Ålesund under specific synoptic conditions . The local wind conditions impact the stratification of the local boundary layer . show that during the ARTIST campaign (15 March–16 April 1998 at Ny-Ålesund) stable conditions were mainly observed with southeast wind and hardly ever with northwest wind. Unstable conditions occurred from 90 to 270 and under light wind conditions. Furthermore, large wind shear was observed to generate turbulence and lead to neutral stratification. This brief summary of previous studies demonstrates the complexities of the local wind conditions present at the AWIPEV station.

We cannot properly describe the circulation in Kongsfjorden from our point measurements at Ny-Ålesund or evaluate the drivers behind the local wind, nor are these processes within the scope of our study. However, it is possible that certain wind patterns are associated with phenomena (shear-induced turbulence, drainage flows from mountains and glaciers transporting cold air into the sub-cloud layer) that modify the P-MPC studied. To evaluate whether the local wind patterns modify the P-MPC, we identified the main modes in the 10 m wind direction and combined them with the circulation weather type to create a proxy for different wind regimes. As expected, three modes can be identified in the surface wind (Fig. 4a). The dominating wind direction (85–165) corresponds to the direction out of the fjord to the open sea. Less pronounced but clearly identifiable are the two other modes that indicate flow from the sea into the fjord (270–345) and the katabatic flow from the glaciers south of the station (200–270). Wind speed above 12 m s−1 was only observed between 90 and 120. Seasonal wind roses are provided in Appendix A2. The frequency with which each surface wind mode was associated with the different weather types during the cloud observation period (June 2016–October 2018) is illustrated in Fig. 4b. For most circulation weather types, the southeasterly surface wind dominated and the northwesterly wind was rare. An exception were the weather types N and NW, for which the northwesterly direction was most common.

To illustrate how the weather type and surface wind direction modes correspond to different wind profiles, the average wind direction profiles based on radiosonde data from June 2016 to October 2018 are shown for weather type W (Fig. 4c). When the surface wind direction was northwesterly, the average direction changed only slightly from 280 in the free troposphere to align with the fjord axis at about 310. The largest variation in the lowest 200 m was exhibited by the southwesterly surface wind direction. The most common regime (surface wind from southeast) had an average profile with free-tropospheric wind from the west, turning south and all the way to the southeast (120) in the lowest 300 m. Figure 4c illustrates why a combination of surface and free-tropospheric wind direction is needed to isolate different patterns. Considering the moderate standard deviation in the wind direction profiles shown in Fig. 4c, it is reasonable to assume that each surface wind direction mode together with the weather type, which describes the mean regional wind direction at 850 hPa, describes a certain wind pattern and thus gives a first estimate of the wind conditions around Ny-Ålesund.

## 3.5 Statistical tools

To test the statistical significance of differences between two or more distributions, the Mood's median test to compare the medians in different populations was used (Sheskin2000). This test was chosen because it does not require normally distributed data and the compared samples can be of different sizes. The median of each population is compared to the median of the distribution including all data, and the Pearson's χ2 test is used to test the null hypothesis that medians from different populations are identical. To reject the null hypothesis thus leads to the conclusion that the different populations have different medians.

The data points in the time series of the variables tested (LWP, IWP, cloud-top temperature, and cloud-base height) are correlated with each other and can not as such be used in the statistical test. We assume that each P-MPC case is independent of the others, and for cloud-top temperature and cloud-base height we use the medians for each case for testing. LWP and IWP were found to vary more within each case, and therefore several data points from every P-MPC case were sampled. For this, we estimated the de-correlation timescale as the time where the autocorrelation function, computed for each P-MPC case individually, reaches zero. For the majority of P-MPC cases there were too many gaps in the data to reliably compute the autocorrelation function, and hence no de-correlation timescale could be estimated. From the values available, the median was calculated and then double the median was used as the de-correlation timescale Δtdcr for all cases. For testing, the data were sampled randomly, with a minimum gap of Δtdcr between the sampled data points.

4 Results and discussion

## 4.1 Occurrence of persistent MPC and other clouds

We first examine the frequency of occurrence of different types of clouds in the observation period of the cloud radar (10 June 2016–8 October 2018) considering the 30 s averaged columns of the Cloudnet product. A cloud was found above Ny-Ålesund 76 % of the time measurements were running. The month-to-month variation was considerable, varying from 40 % to over 90 % (Fig. 5). Averaging for all years, cloudiness was slightly higher in summer (June–August; 80 %) and autumn (September–October; 77 %) and lower in spring (March–May; 69 %) and winter (December–February; 74 %). Intra-annual variation is pronounced in autumn, when cloud occurrence frequency varied from 69 % to 84 %. MPCs (defined here as any profile where co-located cloud liquid and ice are found) were present 41 % of the time, with a somewhat higher frequency in autumn. Liquid-only clouds (profiles with cloud droplets without co-located ice) had an overall occurrence frequency of 14 % and a clear seasonal cycle with most liquid-only clouds occurring in summer and hardly any in winter or spring. Thus, the radiatively important cloud liquid was more often found in mixed-phase clouds, although the contribution of liquid-only clouds was notable in summer. All of the presented figures are given relative to the amount of data available. Figure 5a shows the high data coverage obtained, implying that – with the exception of the first and last months – we can give a reliable estimate of the frequency of cloud occurrence within the detection limits of the instruments.

Figure 4Wind rose for 30 min mean 10 m wind for the cloud observation period (June 2016–October 2018) with the three main modes identified (a) and the relative frequency of occurrence for each weather type (b). Wind direction profiles corresponding to each identified near-surface wind direction mode for weather type W based on radiosoundings from June 2016 to October 2018 (c). The line shows the mean wind direction, and the shaded area shows the mean ± standard deviation at each height, estimated using the method by . Data points with wind speed below 0.5 m s−1 were omitted. N gives the number of soundings available for each mean profile.

Figure 5Monthly and total occurrence frequency of clouds in general and selected specific cloud types (see text for definitions) (b) and coverage of Cloudnet data (a).

The persistent low-level mixed-phase clouds (P-MPCs; see Sect. 3.1 for definition) cover 23 % of the data set, highlighting the relevance of this cloud regime. In total 1412 cases of P-MPC were identified. The “all MPC” and the P-MPC occurrences in Fig. 5 are not directly comparable, since the first one refers to individual profiles and the latter is to a large extent defined by a temporally continuous liquid layer and also includes profiles without a mixed-phase layer detected. P-MPCs were most common in summer (32 %) and occurred less often in winter (15 %) and spring (16 %), with autumn being the intermediate season (24 %). The P-MPC occurrence thus follows the seasonal cycle of cloud liquid occurrence .

For defining the persistence of the liquid layer some thresholds needed to be set, including how long gaps were allowed and the minimum duration required. The choices made (5 min and 1 h) were motivated by the aim for a certain cloud regime, namely a stratiform mixed-phase cloud in the boundary layer. A sensitivity test allowing only 2 min gaps in the liquid layer showed the only major difference to be in the occurrence frequency of P-MPCs, while the properties of the clouds or the seasonal cycle of P-MPC occurrence did not differ substantially.

## 4.2 P-MPC properties and regional wind direction

Figure 6 shows the occurrence of each weather type (used to determine the regional free-tropospheric wind direction; see Sect. 3.3) in our period of study, and the fraction of those times when a P-MPC was identified. In general, NE, SE, and NW were less common than the other wind directions. For a given weather type, the fraction of P-MPC occurrence varied considerably. Almost a third of the time when winds were from the west (W), a P-MPC was found at Ny-Ålesund. Weather types S, SW, NW, and anticyclonic were also favorable for P-MPC. Based on an evaluation of sounding profiles, the most common free-tropospheric wind direction for the weather type anticyclonic was west (not shown). On the other hand, winds from north and east (weather types N, NE, and E) brought P-MPCs to the site less often. The weather types which are most commonly associated with P-MPCs can be determined by combining the occurrence frequency of each weather type and its P-MPC fraction (Fig. 6). Consequently, P-MPCs were most often associated with the weather types W, SW, and anticyclonic, which include almost half (48 %) of all profiles.

Figure 6The frequency of occurrence of each weather type and the fraction with P-MPC presence.

The distributions of liquid layer base height, LWP, and IWP and their dependence on wind direction are presented in Fig. 7. The base of the liquid layer was usually between 540 and 1020 m above the surface, with mean and median liquid base height of 860 and 760 m, respectively. The typical P-MPC thus lies above the fjord at a height fairly close to the mountaintops. Fewer P-MPCs were associated with weather types NE, E, and SE, and with mean liquid base heights well above 1 km these were found at larger altitudes than most of the P-MPCs. The mean LWP for P-MPCs was 35 g m−2 with a standard deviation of 45 g m−2. On average most liquid was found in the P-MPC in weather type SW (49 g m−2), and the least was found in weather type NE (12 g m−2). However, the variability within each weather type was larger than the differences between the weather types. The IWP distributions are strongly skewed (Fig. 7c) towards low values. Zeros were ignored, but all nonzero values were included. For all P-MPCs, the mean and median IWPs were 12 and 2.1 g m−2, respectively. Between the different weather types, the mean (median) varied from the 5.6 (1.1) g m−2 of weather type SE to 17 (6.2) g m−2 of weather type NW. The weather types NW, W, and SW stand out in terms of high IWP and have a mean IWP of 16 g m−2. Overall, the westerly weather types (SW, W, and NW) were associated with lower P-MPCs and with more liquid and ice (mean LWP 42 g m−2), while the easterly weather types (SE, E, and NE) were less common, distinctly higher, and connected to the lowest average LWP and IWP.

Figure 7Height of the P-MPC liquid layer base (a), the LWP (b) and IWP (c) distributions for each weather type, and all P-MPCs. The number of P-MPC cases for each weather type is given in (a). The box shows the 25th, 50th, and 75th percentiles, the dot the mean, and the whiskers the 5th and 95th percentiles. The medians for different weather types were found to differ on a 95 % confident level.

Large-scale advection and air mass properties are known to influence MPC properties (Mioche et al.2017; Qiu et al.2018, amongst others). Previous studies suggest that at Svalbard northerly flow is often associated with cold air masses originating from the central Arctic and that southerly flows bring warmer and more humid air from lower latitudes . Furthermore, the open sea west of the Svalbard archipelago might act as a local source of humidity and heat. Here we use temperature at 1.5 km (corresponding to the 850 hPa level) and integrated water vapor (IWV) from the MWR to represent the atmospheric temperature and humidity conditions under which the P-MPCs were occurring. In agreement with previous studies, Fig. 8 shows that the highest average IWV and warmest temperatures were associated with southerly winds, while the lowest average IWV and coldest temperatures were associated with northerly winds. The domain considered for the weather type (Fig. 1a) is too small to describe large-scale advection or air mass origin, but Fig. 8 suggests the weather type is nonetheless a useful proxy for air mass properties. The average IWV and 1.5 km temperature can explain the first-order variation in P-MPC occurrence and LWP between weather types. The south-southwesterly winds are warm and humid and are associated with frequent occurrence of P-MPC with relatively high amounts of liquid, compared to the north-northeasterly winds, which are drier and colder and are associated with less frequent P-MPC occurrence and lower LWP (Figs. 6, 7b, and 8a). Owing to the complexity of ice microphysical processes, such a direct relationship cannot be found between atmospheric humidity and temperature (Fig. 8) and IWP (Fig. 7c). On the other hand, as already noted above, Fig. 7 shows a clear contrast between the properties of easterly and westerly P-MPCs. These differences cannot be explained by the IWV and 1.5 km temperature distributions, which are rather similar for weather types W and E. Hence, atmospheric temperature and humidity are important but not the only relevant forcing for P-MPC at Ny-Ålesund.

Figure 8IWV (a) and 1.5 km temperature (b) for time periods with P-MPC present for each weather type. Boxes and whiskers as in Fig. 7. The medians were found to differ on a 95 % confidence level.

The influence of the island and its orography clearly affects the height of the liquid layer (Fig. 7a). The median altitude of the P-MPC base with easterly winds (weather types NE, E, and SE) was above the height of the mountaintops, suggesting that the clouds were usually advected to the site above the mountains rather than forming locally in the fjord. The P-MPCs associated with easterly winds were also less frequent (Figs. 6 and 7a). If we assume the majority of observed P-MPC being of an advective nature, the low occurrence frequency with winds from the east would imply less cloud formation over the island compared to over the sea or dissipation of cloud fields while being advected over the island. found less low (below 3 km) MPC over land than over sea in the Svalbard region in spring and winter, while in summer and autumn the differences were small. studied liquid-containing clouds in the Arctic and found fewer low (below 3.36 km) liquid-containing clouds above Svalbard than over the surrounding sea in all seasons. Although direct comparison is not possible due to inconsistencies in the observation techniques, cloud sampling, and the considered area, the mentioned studies all indicate that the influence of the Svalbard archipelago decreases the amount of low liquid-bearing clouds.

The combination of the effects of large-scale advection and air mass properties, as well as the influence of the Svalbard archipelago, can provide an explanation for the dependence of the P-MPC properties on weather type presented in Figs. 6 and 7. Southwesterly and westerly free-tropospheric winds were associated with most P-MPCs and the highest average LWP and IWP, likely due to higher amounts of humidity available from lower latitudes. The southeasterly to northeasterly winds had the least P-MPCs and comprise the lowest average LWP and IWP, related to the drier air masses from the north and less favorable conditions for cloud formation over the island. Other mechanisms can be considered to further explain the observed IWP variation. Ice formation could be enhanced in the cold temperatures for weather types N and NE (Fig. 8), whereas the higher IWP for weather types SW, W, and NW might be related to larger amounts of supercooled liquid available in the P-MPCs (Fig. 7b, c) or higher aerosol concentration in air masses advected from lower latitudes.

## 4.3 Seasonality

The seasonal variation in the studied P-MPC properties and atmospheric conditions at Ny-Ålesund are presented in Fig. 9. In agreement with previous studies , the highest average temperature and humidity are found in summer and the lowest in winter and spring (Fig. 9a, b). The height of the P-MPC shows a clear seasonality, with lower liquid base height in summer and higher in winter (Fig. 9c). evaluated 5 years of low-level clouds (cloud base below 2 km) observed at Utqiaġvik (previously known as Barrow), Alaska, and also found a seasonality in cloud height with a minimum in summer. Furthermore, these results are in agreement with the seasonality in cloud height at Ny-Ålesund reported by . The IWP distributions show a clear seasonality, with low values in summer and autumn and a clear maximum in spring (Fig. 9e). The low IWP in summer and autumn (median 0.2 and 1.0 g m−2, respectively) can be attributed to relatively warm temperatures close to 0 C. The median IWP in spring (7.5 g m−2) is almost 2-fold of the median IWP in the winter (4.0 g m−2), which can hardly be attributed to the different temperature conditions (Fig. 9b). The higher IWV in spring compared to winter (Fig. 9a), however, can play a role. Furthermore, the high IWP in spring could be related to the generally higher aerosol loading in the Arctic atmosphere in the late winter and spring, a time period also known as the Arctic haze season .

Figure 9IWV (a), 1.5 km temperature (b), liquid base height (c), LWP (d), and IWP (e) distributions for each season. Only time periods with P-MPCs present are included. Boxes and whiskers as in Fig. 7; the medians were found to differ on a 95 % confident level.

Conversely, the LWP distributions show a minimal seasonality despite the seasonal variation in IWV and 1.5 km temperature related to the P-MPC (Fig. 9a, b, d). The highest (lowest) median LWP in summer and spring (winter) was 24 g m−2 (18 g m−2), and the seasonal mean values varied from 33 to 36 g m−2. Note that this result does not imply a lack of seasonal variability in overall cloud LWP (see Fig. 5 in Nomokonova et al.2019a), only in the specific cloud regime evaluated. One challenge of the algorithm to identify the P-MPC are thick liquid layers where Cloudnet only identifies the lowest parts of the layer as containing liquid. The problem was partly mitigated by relaxing the criteria for liquid presence at cloud top; nonetheless we find cases with a thick liquid layer that do not fulfill the criteria of a liquid-topped mixed-phase layer and the rest of the cloud gets cut off (see Fig. 2 at 12:00 UTC on 30 May 2018). This artificially limits our data set to clouds where the liquid layer is thin enough, and there might be some clouds with more liquid that are not included in our analysis. Considering the LWP distributions were skewed towards lower values (Fig. 7b), these cases are likely to be a minority for the cloud regime considered. However, it is possible that the average LWP is somewhat underestimated. In addition, it could be that the cloud detection algorithm limits the considered cases to a specific LWP regime, which results in the lack of seasonality in the LWP of the P-MPC.

Since P-MPC properties (excluding LWP) as well as atmospheric temperature and humidity vary seasonally, a seasonal dependency in wind direction could explain the weather-type-dependent variations in P-MPC properties found in Sect. 4.2. To examine this possibility, Fig. 10 shows the proportion of P-MPC observations in each season for every weather type. The observation period of 2.5 years from June 2016 to October 2018 together with the seasonal variation in P-MPC occurrence (Fig. 5) lead to the uneven distribution of data between seasons. Overall, the summer months contribute most to the data set. However, there are no extensive differences found between the weather types. Most noteworthy is the high spring and low autumn occurrence of NW, which might contribute to the high IWP for this weather type (Fig. 7c). Furthermore, N and E were relatively more common in winter, N and SE more common in spring, and N less common in autumn. Given the lack of a distinct signal, we believe the seasonal variation in wind direction plays a minor role in the weather-type-dependent differences in P-MPC occurrence and properties described in the previous section.

Figure 10Distribution of seasons in the studied data set for each weather type. Only time periods when a P-MPC was present are included to evaluate the possible impact of wind direction seasonality on cloud properties and occurrence.

We further compare properties of the P-MPCs at Ny-Ålesund and their seasonal variation to observations of similar cloud regimes at other Arctic sites. Only studies that comprise at least 1 year of observations were considered. evaluated MPCs observed at the 1-year-long Surface Heat Budget of the Arctic Ocean (SHEBA) campaign and found an annual average LWP and IWP of 61 and 42 g m−2, respectively. Both IWP and LWP were found to have a maximum in late summer and autumn. The study did not explicitly focus on low-level clouds but found that 90 % of the observed MPCs had a cloud base below 2 km. focused on single-layer mixed-phase stratus at Eureka, Canada, and reported seasonal mean LWP to vary between 10 and 50 g m−2. show monthly mean values for LWP at Utqiaġvik to vary from 10 to 100 g m−2 and for IWP from 10 to 25 g m−2. Similarly to SHEBA, at both Eureka and Utqiaġvik the maximum LWP was found in autumn. However, at Eureka as well as Utqiaġvik a maximum in the amount of ice in MPCs was found in spring as well as autumn. The differences in seasonal cycles of LWP and IWP at different sites could be due to different forcing conditions, in addition to the choice of the cloud regime that might also play a role. included all single-layer clouds below 3 km and found that most of the LWP distribution was within 0 to 100 g m−2, with slightly higher values in the data set from SHEBA than Utqiaġvik. The average figures are comparable to those observed for P-MPC at Ny-Ålesund, although the mean values in our study are at the lower end of the range reported at Utqiaġvik and SHEBA.

Finally, the seasonal variation in P-MPC occurrence is compared with previous studies in the Svalbard region. as well as report most clouds in summer and early autumn above Ny-Ålesund, agreeing with our findings. Conversely, identified most low-level (below 3 km) MPCs in the Svalbard region in autumn and a minimum in occurrence in summer based on the synergy of the measurements from CloudSat and CALIPSO. P-MPCs commonly contain very low amounts of ice, which might be below the sensitivity limit of the satellite observations, explaining some of the disagreement. Furthermore, were missing clouds below 500 m due to the blind zone of CloudSat, and since clouds generally are lower in summer this would lead to a higher fraction of missed clouds in this season. In any case, considering the large month to month variation in cloud occurrence (also shown by Shupe et al.2011), different results when considering different time periods can be expected. Our time series might still not be long enough to give a precise estimate of the seasonal variation in cloud occurrence frequency.

## 4.4 Surface coupling

Figure 11a shows the fraction of observed P-MPCs classified as coupled, predominantly decoupled and fully decoupled in each season. A total of 63 % of all observed P-MPC cases were found to be fully decoupled, and only 15 % were coupled. The degree of coupling had a clear seasonal cycle, with decoupling being the dominant mode in autumn and winter and most coupled P-MPCs occurring in summer. The observed seasonality in the surface coupling of P-MPC could be related to the overall higher lower-tropospheric stability in winter, which could limit the coupling of the cloud. Previous studies have found that the coupling of low Arctic MPCs depends on the proximity of the cloud to the surface since the cloud-driven mixing layer is more likely to reach the surface if the cloud is low . Also, in our data set the median cloud-base height for decoupled P-MPCs (1010 m) is considerably larger than the median cloud-base height of the coupled P-MPCs (620 m) (Fig. 11b). For P-MPC with liquid base heights of more than 1.5 km, coupling to the surface was not observed. P-MPCs were on average higher in winter and lower in summer (Fig. 9c), which could partly explain the seasonal variation in the frequency of surface coupling.

Figure 11The fraction of P-MPC cases classified as coupled, predominantly decoupled, and fully decoupled in each season and for the entire data set (a), and the distribution of the liquid layer base height in the coupling classes (b). Boxes and whiskers as in Fig. 7; the medians were found to differ on a 95 % confident level.

To evaluate the effect surface coupling has on cloud properties, we only considered P-MPC in weather types SW and W in order to limit the different factors in play. These clouds include the full range of coupling states and cover one-third of the data set (Fig. 7a). The coupled P-MPCs had more liquid than the fully and predominantly decoupled P-MPCs (Fig. 12a). The median LWP did not differ significantly between the predominantly and fully decoupled P-MPCs (25 and 28 g m−2, respectively), while the median LWP for coupled cases was clearly larger (47 g m−2). Differences in IWP between the coupling states were small (Fig. 12b). The medians did not vary significantly (from 11 to 12 g m−2), but the larger IWP values (between 30 and 100 g m−2) were less likely for the coupled P-MPC. From the LWP and IWP distributions it follows that the total amount of condensed water (LWP+IWP) was higher for coupled than predominantly or fully decoupled P-MPC. This suggests either a source of humidity from the surface that is not available for the decoupled P-MPC or a smaller sink.

Figure 12Comparison of LWP (a), IWP (b), and cloud-top temperature (c) distributions between P-MPC in weather types W and SE with different degrees of surface coupling. The dashed line and the numbers on top show the median value of each distribution. The bin size for LWP and IWP is 10 g m−2 and for cloud-top temperature 3 C. The medians were found to differ on a 95 % confident level for LWP and cloud-top temperature.

Many ice microphysical processes have a temperature dependency , and the observed differences in LWP and IWP distributions between coupled and decoupled P-MPCs could be caused by different sampling across the temperature range. Observed cloud-top temperatures ranged from −28 to +5C, with most P-MPCs occurring at the warm end of this range (Fig. 12c). As the persistent liquid layer is the defining feature of the P-MPCs, it is not surprising that they occurred more often at warmer temperatures where liquid is generally more abundant. The coldest P-MPCs (cloud-top temperatures below −18C) were always decoupled and occurred in winter and early spring. The cloud-top temperature distributions were very similar for the coupled and predominantly decoupled P-MPCs, suggesting that the differences in IWP and LWP distributions between these two groups can not be explained by a varying frequency of different temperature regimes. The cloud-top temperature distribution of fully decoupled P-MPC differs from that of the predominantly decoupled and coupled clouds by having a larger number of cold cloud tops and a smaller peak at the warm end of the distribution. Yet the IWP and LWP distributions do not differ substantially between fully and predominantly decoupled P-MPCs. Although the observed differences in LWP and IWP between coupled and fully decoupled P-MPCs could be caused by differences in temperature, we cannot explain the differences between predominantly decoupled and coupled clouds or the similarity of the predominantly and fully decoupled P-MPCs simply from the cloud-top temperature distributions.

The analysis presented only included weather types SW and W. These weather types are amongst the weather types with the largest average LWP and IWP. The variation in LWP and IWP between coupled and decoupled P-MPCs for the other weather types would therefore be smaller in absolute numbers. Including all weather types, the medians for LWP for coupled and predominantly and fully decoupled were 34, 22, and 20 g m−2, and the medians for IWP were 7.5, 9.4, and 9.4 g m−2, respectively. The outcome that coupled P-MPC had more liquid and that differences in IWP were small is the same. One needs to keep in mind that these numbers are dominated by the westerly weather types which cover the bulk of the data. It is possible that different relationships between cloud properties of coupled and decoupled clouds would be found for weather types which have distinctly different mean wind conditions. While we cannot conclude that the presented results hold for all situations occurring at Ny-Ålesund, they describe the most common conditions.

The comparison between the coupling detection from sounding and the new method based on MWR and surface observations implied that the new method is more inclined to consider a profile decoupled (Sect. 3.2.3). Yet, the similarity of the LWP and IWP distributions for predominantly and fully decoupled P-MPCs suggests that these groups were very similar. Considering cloud properties, it does not seem that the predominantly decoupled clouds would be mistakenly considered more decoupled than they are. It is possible that our method erroneously classifies weakly coupled P-MPC as predominantly decoupled and that in these cases the interaction with the surface is limited and does not modify the cloud properties, considerably leading to similar LWP and IWP distributions for these clouds and the actually decoupled P-MPC. Accordingly, we conclude that decoupling might be overestimated, but this does not have serious consequences on the results on cloud properties. Considering the different estimates (Figs. 3d and 11a), we can regard 63 %–82 % of the P-MPCs to be decoupled and 15 %–33 % to be coupled. Moreover, intermittent turbulence and the coupling it may lead to are rather challenging for our approach, as the thermodynamic profile takes time to adjust. However, the turbulent transport of heat can be assumed to be similar to the transport of any other scalar. If the turbulence that occurred was too short-lived to modify the temperature profile distinctly, it would also be unlikely to transport great amounts of water vapor or aerosols to the cloud layer.

, , and have evaluated the coupling of low clouds during the ASCOS campaign (August–September 2008) using different methods and slightly different time periods and observed decoupling from the surface 75 %, 72 %, and 76 % of the time, respectively. Their measurements were mostly of clouds above sea ice, and for a shorter time period. The results are therefore not directly comparable with the multiyear statistic presented here. Moreover, the mechanisms that lead to decoupling at ASCOS were likely different than at Ny-Ålesund. Like , but unlike , we found a difference in LWP between coupled and decoupled clouds (Fig. 12a). If we assume that the (sea) surface can provide a source of moisture for the P-MPC, coupling could add moisture to the cloud layer and lead to a higher total water path. Considering the small differences in IWP (Fig. 12b), it does not seem that the surface would be an important source for INPs, or there are some other mechanisms that limit ice formation in coupled clouds where more liquid water is present. The observed seasonality in the surface coupling of P-MPC (Fig. 11a) could be related to the overall higher lower-tropospheric stability in winter, which could limit the coupling of the cloud as well as to the lower cloud-base height in summer (Fig. 9c) that makes it easier for the cloud to couple to the surface due to its proximity.

## 4.5 Local wind patterns around Ny-Ålesund

The effects local winds have on the P-MPC were evaluated using the weather type together with the surface wind direction as a proxy for the wind conditions at Ny-Ålesund (Sect. 3.4). The most common wind situation for the P-MPC at Ny-Ålesund is a southeasterly surface wind underlying westerly and southwesterly upper winds (Figs. 4, 6). Hence, the wind turns from the surface upwards to the almost opposing direction by 1.5 km height (Fig. 4c). Directional wind shear is therefore commonplace for P-MPC at Ny-Ålesund (Fig. 4b), either in or below the cloud layer. The magnitude of the wind direction change varies with the free-tropospheric wind. The only exception are weather types N and NW, for which the most common surface wind is northwesterly, and the wind does not turn, or only turns slightly, with increasing altitude. A further consideration related to the surface wind direction is the history of the boundary layer. The air has experienced very different surface properties when moving from open sea to land with northwesterly surface wind or from mountainous, often snow- and ice-covered terrain to a flat sea surface with southeasterly surface wind.

The influence of local winds on the P-MPC was found to be limited. Figure 13a shows the fraction of time with P-MPC occurring (similarly to Fig. 6) for each weather type and surface wind direction combination. Weather types cyclonic and anticyclonic are somewhat hard to interpret, as these are associated with varying free-tropospheric wind directions above the site and were therefore not included. The low number of cases with southwest and northwest surface wind limits the possibilities to compare different surface wind regimes for most of the weather types. For weather types SW and W the southwest surface wind was associated with higher frequency of cloud occurrence compared to the southeast surface wind. In contrast, for weather type N northwest surface wind had P-MPCs most often and southwest the least. Based on this analysis no overall tendency for a certain surface wind direction or the amount of directional shear between the surface and the free-tropospheric wind to increase or decrease P-MPC occurrence was found.

Figure 13Fraction of time with P-MPC occurring for each surface wind direction and weather type regime (a). The size of the dots represents the amount of data available to compute the value. The fraction of P-MPC cases classified as coupled, predominantly decoupled, and fully decoupled for each surface wind direction mode (b).

Regarding P-MPC properties, no strong relationships with surface wind direction were identified. Only the main findings are summarized here and further details are provided in Appendix A1. Considering weather types N, W, and SW, which have the most cases across different surface wind directions, no statistically significant differences were found in the median liquid base height or cloud-top temperature. The northwest surface wind was associated with the highest median LWP, possibly due to the higher level of humidity available over the open sea. The southwest surface wind was associated with a significantly higher IWP for weather types W and SW (median IWP 16 and 18 g m−2, respectively). However, these variations in LWP and IWP were not found for all three weather types analyzed.

Local winds in Kongsfjorden were quite apparently connected to the coupling of the P-MPC (Fig. 13b). Coupling was most common with northwest surface wind (from the sea) and least common with the southeast surface wind (towards the sea), and the same behavior was found for every season despite the seasonality of both surface wind direction and cloud coupling (see Appendix A2). For the P-MPC to be thermodynamically decoupled from the surface, a stably stratified layer needs to exist between the surface and the cloud base. Argentini et al. report a dependence of surface layer stratification on wind direction (, Fig. 5). Stable conditions were most often found with southeast surface wind, for which only 8 % of the P-MPCs were considered coupled. On the other hand, stable conditions were rare with northwest surface wind, for which 37 % of the P-MPCs were coupled. The near-surface wind from the southwest and southeast is often related to flows from the glaciers (; Sect. 3.4) that bring cold air down to the valley in a shallow layer close to the surface. Such a cold surface layer is very efficient in decoupling the cloud and acts against the cloud-driven turbulence that could otherwise couple the P-MPC to the surface. This effect might be stronger with southeast than southwest surface wind, since the katabatic winds from the southwest are weaker (Fig. 4a). The differences in the coupling of the P-MPC with varying wind conditions can be explained by the differences in stratification of the lower boundary layer under different surface wind conditions. We conclude that the surface wind has the potential to modify the conditions in the boundary layer, which in turn can act to suppress coupling.

The local influence on coupling makes assessing the connection between coupling and cloud properties more challenging. The cloud might have been coupled to the surface while over the sea, and when it was advected into the Kongsfjorden valley the local wind changed in the sub-cloud layer, leading to decoupling. It is also difficult to evaluate coupling and local winds separately, because most coupled clouds were associated with northwest surface wind (Fig. 13b). Coupled P-MPCs had higher LWP than decoupled P-MPCs (Fig. 12a), and P-MPCs associated with northwest surface wind had higher LWP than those occurring with other surface wind directions (Fig. A1b). Perhaps the higher LWP is related to the combined effect of the two: more humidity is available from the open sea than over land and coupling is required for the water vapor to be transported from the surface to the cloud layer. There is a relationship between surface coupling and the local wind conditions at Ny-Ålesund, but to understand the impact of the combined effects on P-MPC properties would require further studies.

5 Conclusions

We present 2.5 years of vertically resolved cloud observations carried out at the AWIPEV station at Ny-Ålesund. Methods to identify persistent low-level mixed-phase clouds (P-MPCs), their coupling to the surface, and the regional and local wind conditions were developed. We found P-MPC to occur 23 % of the time, most often in summer and least often in winter. The cloud base was typically 0.54–1.0 km high, LWP was 6–52 g m−2, and IWP was 0.2–12 g m−2. P-MPCs were found to occur at higher altitudes in winter and lower altitudes in summer. LWP presented a lack of seasonal variation, possibly due to the selection of the cloud regime in this study. On the other hand, IWP had a clear seasonal dependence. IWP was low in the relatively warm months of summer and autumn and had a clear maximum in spring. The frequency of occurrence was found to depend on free-tropospheric wind direction, and most P-MPCs were associated with westerly winds. The height of the cloud was strongly influenced by orography. Less frequent P-MPCs and with higher cloud-base height were found with easterly winds compared to westerly winds, and these clouds had lower LWP and IWP. The most common surface wind direction in Kongsfjorden is from the southeast, but this is typically underlying synoptic winds from westerly directions. Local winds were not found to impact the occurrence or the height of the P-MPCs, but for some free-tropospheric wind directions the surface wind direction was related to variations in LWP and IWP. P-MPCs were mostly decoupled (63 %–82 % of the time), and coupling occurred most often in summer and for clouds close to the surface. Coupled P-MPC had a higher LWP than decoupled P-MPC, but no differences in IWP were found. Furthermore, the local wind patterns appeared to be related to surface coupling; specifically, the P-MPCs with surface wind directions associated with glacier outflows were more commonly decoupled. The variation in median LWP between different wind directions at 850 hPa was larger than the variation found between different surface wind regimes or coupling states. On the other hand, IWP was found to vary with regional and local wind direction as well as season, but no dependency with coupling to the surface was found. We conclude that while the regional to large-scale wind direction was important for P-MPC occurrence and their properties, the local-scale phenomena such as surface coupling and the local flow in the fjord also had an influence.

Our results suggest that the P-MPC water properties can be influenced by the processes in the local boundary layer. The observed LWP values are in the range where the clouds are not yet fully opaque, and changes in LWP will have an impact on the radiative forcing of clouds at Ny-Ålesund . For numerical models to correctly describe low-level MPCs' ice and liquid water content, and hence the radiative effect, the boundary layer dynamics need to be accurately described. In Ny-Ålesund, and in other Arctic fjords, this requires that local wind in the fjord is represented, and thus a description of the orography and key surface properties (temperature, snow cover, etc.) needs to be accounted for in the model.

Long-term data sets are valuable for evaluating models since the evaluation can be carried out in a statistical manner instead of case-by-case basis. The data set presented in this paper can be used for model comparison, to provide insight into model performance regarding low-level MPCs in the complex Arctic fjord environment. In addition, the results presented here provide background information that aids the interpretation of case studies underway from recent measurement campaigns . In this study, the effects of aerosols acting as ice-nucleating particles or cloud condensation nuclei have not been evaluated. Also, the cloud microphysical processes taking place should be considered in more detail. Further work is thus needed to understand the relationships between various processes controlling the properties and development of low-level MPCs at Ny-Ålesund.

Appendix A: Details on the relationship between local wind conditions and P-MPC

## A1 P-MPC properties

The results of the analysis on P-MPC properties for different wind regimes is provided here, and some possible mechanisms are contemplated. The cloud properties associated with different surface wind directions were compared separately for each weather type. Only weather types N, W, and SW were considered (Fig. A1) in order to have a sufficient amount of data (at least 30 cases) in each group being compared (Fig. 13a). The median liquid base height did not differ significantly (on a 95 % confidence level) for any of the three weather types evaluated. The northwest surface wind was associated with the highest median LWP; however, for weather type SW the differences were not statistically significant. For weather type N the median LWP for northwest surface wind was 22 g m−2 compared to 12 and 7.8 g m−2 of southeast and southwest surface winds, respectively. Also, for weather type W the northwest surface wind was associated with the highest median LWP (39 g m−2), however the lowest median LWP was with southeast surface wind (18 g m−2). The median IWP varied insignificantly (from 7.8 to 9.7 g m−2) for weather type N. For weather type SW, the southwest surface wind had the highest median IWP at 18 g m−2, almost double that of the median of southeast (10 g m−2) and northwest (9.1 g m−2) surface winds. Similarly, for weather type W the median IWP for the southwest was 16 g m−2 and only 11 and 9.6 g m−2 for southeast and northwest surface winds. Because of the temperature dependence of many microphysical processes, it would be possible that the observed differences were a result of different temperature regimes dominating in the compared groups. However, no statistically significant difference in the cloud-top temperature distributions was found (Fig. A1d).

Figure A1P-MPC liquid layer base height (a), LWP (b), IWP (c), and cloud-top temperature (d) distributions for selected weather types and surface wind directions. Boxes and whiskers as in Fig. 7. The medians were found to differ (on a 95 % confidence level) in LWP for N and W and IWP for SW and W.

Local conditions evidently modify the wind field in the fjord (Sect. 3.4), but whether this affects the P-MPC is not as easily determined. Although we find some differences in the P-MPC occurrence and properties with different local low-level wind patterns (Figs. 13 and A1), these could also be due to the large-scale conditions related to different local circulation patterns. We here consider some phenomena that might be taking place. The near-surface wind from the southeast could hinder the low P-MPC residing over the sea from advecting into the fjord where the observations were taking place. This would lead to higher cloud-base height for the southeast surface wind regime, or a lower frequency of occurrence, as the lowest P-MPC would be limited. For both weather types N and W the northwest surface wind had the lowest 25 percentiles of the liquid layer base height (lower edge of the boxes in Fig. A1a). Figure 13a gives no indication that the southeast surface wind would have been related to an overall lower frequency of occurrence. Although the lowest P-MPCs were more often associated with northwest surface wind, liquid base height below 400 m was also not that common for this wind regime. Hence, it seems that the southeast surface wind did not substantially prevent the P-MPC on the sea from advecting into the fjord. Considering Figs. 4c and A1a together, the depth of the layer where wind is found to deviate the strongest from the free-tropospheric wind direction is below the median P-MPC base height, and the 25th percentile is above the depth of the layer where on average the wind is in alignment with the surface wind direction. Hence, many of the P-MPCs reside in a layer where the wind direction is changing with altitude, or just above it. The wind shear could induce turbulence, which in turn could affect the properties of P-MPCs, and it might be influencing vertical fluxes of heat, moisture, and aerosols. These kind of processes could explain the differences found in IWP and LWP between different wind regimes. However, to examine these processes would require a more sophisticated description of the local circulation and turbulence in the boundary layer than was used here.

## A2 Seasonality of surface wind direction and P-MPC coupling

Seasonality in near-surface wind and the degree of P-MPC coupling with different surface wind directions are presented. Figure A2 shows the wind rose for each season for 10 m wind in the studied period, June 2016–October 2018. Differences between the seasons are present in the relative importance of the three surface wind modes, in agreement with and . The summer months stand out with more common northwesterly winds, which has previously been attributed to sea breeze . Subsequently, the other directions are less frequent. In autumn and winter the northwesterly winds almost completely disappear. The seasonal variation is likely due to the different degree at which the drivers (e.g., sea breeze circulation, katabatic flow, channeling of free-tropospheric wind along the fjord) act in different seasons.

Figure A2Wind rose for 30 min mean 10 m wind for each season (a–d) in the cloud observation period (June 2016–October 2018). The fraction of P-MPC cases classified as coupled, predominantly decoupled, and fully decoupled for each surface wind direction mode in each season (e–f).

The relationship between surface wind and P-MPC coupling is similar in all seasons except summer (Fig. A2e–h). In winter, autumn, and spring coupling with southeast surface wind was rare or nonexistent. Coupling mostly occurred with northwest surface wind. The reasons follow those given in Sect. 4.5: the more (less) stable stratification of the lower boundary layer associated with the southwest (northwest) wind, probably related to the cold outflow from the glaciers that increase the stability of the sub-cloud layer promoting decoupling. In summer the situation is somewhat different from the other seasons. Southeast wind was still related to the fewest coupled P-MPCs, but the differences between different wind directions were smaller. Furthermore, the coupling frequency with southwest wind was very similar to that of northwest wind. The wind roses for each season (Fig. A2a–d) suggest a variation in boundary layer dynamics in summer, which could be contributing to the altered relationship between surface wind direction and the frequency of P-MPC coupling. Moreover, as discussed in Sect. 4.4, the overall lower stability in the boundary layer as well as the lower cloud-base height in summer enhance surface coupling compared to the other seasons. Hence, local wind conditions seem to have less importance in summer, although the interaction with the local boundary layer is present in all seasons.

Code and data availability
Code and data availability.

The Cloudnet data are available at the Cloudnet website (http://devcloudnet.fmi.fi/, last access: 25 June 2019). The radiosonde data are available in PANGAEA (https://doi.org/10.1594/PANGAEA.845373, for 1993–2014; https://doi.org/10.1594/PANGAEA.875196, for 2015–2016; search term “project:label:AC3 ny-alesund radiosonde” afterwards). The meteorological surface observations are available in PANGAEA under the search term “Continuous meteorological observations at station Ny-Ålesund”. The MWR data are also available in PANGAEA (https://doi.org/10.1594/PANGAEA.902183, ). The software used for the median test was available courtesy of . The cloud microphysical data set is currently under review for PANGAEA (https://doi.org/10.1594/PANGAEA.898556, ). Topography data in Fig. 1 are provided by  (panel a) and the  (panel b). Color maps used in Figs. 1a, 4a, 7a, and A2a–d are provided by .

Author contributions
Author contributions.

RG carried out the method development, statistical analysis, and visualization of the results and prepared the manuscript with contributions from all co-authors. UL, SK, and MS contributed with conceptualization, research supervision, and discussions of the results. KE oversaw data management at the University of Cologne and advised in the analysis and selection of data sets. MM provided the long-term radiosonde data set and insights into the local conditions at Ny-Ålesund.

Competing interests
Competing interests.

Matthew D. Shupe is an editor for the special issue “Arctic mixed-phase clouds as studied during the ACLOUD/PASCAL campaigns in the framework of (AC)3”.

Special issue statement
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.

Acknowledgements
Acknowledgements.

We gratefully acknowledge the funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – project no. 268020496 – TRR 172, within the Transregional Collaborative Research Center “ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes, and Feedback Mechanisms (AC)3”. Contributions by Stefan Kneifel were funded by the German Research Foundation (DFG) under grant KN 1112/2-1 as part of the Emmy-Noether Group “Optimal combination of Polarimetric and Triple Frequency radar techniques for Improving Microphysical process understanding of cold clouds (OPTIMIce)”. Contributions by Matthew Shupe were funded by the U.S. National Oceanic and Atmospheric Administration's Earth System Research Laboratory. We wish to thank the Finnish Meteorological Institute, the Aerosol, Clouds, and Trace Gases Research Infrastructure (ACTRIS), and especially Ewan O'Conner for the Cloudnet algorithm. Furthermore, this work would have not been possible without Tobias Marke's assistance regarding Cloudnet and his providing the circulation weather type, as well as Tatiana Nomokonova, who compiled the cloud microphysical data set and assisted in the use of the MWR data. We further wish to thank Christoph Ritter and Maximilian Mahn for sharing ideas and insightful discussions and Birte Kulla for providing Fig. 1b. Last but not least, we are grateful for the AWIPEV station staff for technical support, maintenance, and operation of the instruments.

Financial support
Financial support.

This research has been supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation, project no. 268020496 and grant no. KN 1112/2-1) and the U.S. National Oceanic and Atmospheric Administration's Earth System Research Laboratory.

Review statement
Review statement.

This paper was edited by Radovan Krejci and reviewed by three anonymous referees.

References

Amante, C. and Eakins, B. W.: ETOPO1 1 Arc-Minute Global Relief Model: Procedures, Data Sources and Analysis, NOAA Technical Memorandum NESDIS NGDC-24, National Geophysical Data Center, https://doi.org/10.7289/V5C8276M, 2009. a, b

Argentini, S., Viola, A. P., Mastrantonio, G., Maurizi, A., Georgiadis, T., and Nardino, M.: Characteristics of the boundary layerat Ny-Ålesund in the Arctic during the ARTIST field experiment, Ann. Geophys., 46, 2, https://doi.org/10.4401/ag-3414, 2003. a, b, c, d, e

Beine, H., Argentini, S., Maurizi, A., Mastrantonio, G., and Viola, A.: The local wind field at Ny-Ålesund and the Zeppelin mountain at Svalbard, Meteorol. Atmos. Phys., 78, 107–113, 2001. a, b, c, d, e, f, g

Brooks, I. M., Tjernström, M., Persson, P. O. G., Shupe, M. D., Atkinson, R. A., Canut, G., Birch, C. E., Mauritsen, T., Sedlar, J., and Brooks, B. J.: The Turbulent Structure of the Arctic Summer Boundary Layer During The Arctic Summer Cloud-Ocean Study, J. Geophys. Res.-Atmos., 122, 9685–9704, 2017. a, b, c

Cesana, G., Kay, J., Chepfer, H., English, J., and De Boer, G.: Ubiquitous low-level liquid-containing Arctic clouds: New observations and climate model constraints from CALIPSO-GOCCP, Geophys. Res. Lett., 39, L20804, https://doi.org/10.1029/2012GL053385, 2012. a

Chang, L., Song, S., Feng, G., Zhang, Y., and Gao, G.: Assessment of the uncertainties in Arctic low-level temperature inversion characteristics in radio occultation observations, IEEE T. Geosci. Remote, 55, 1793–1803, 2017. a

Choi, Y.-S., Ho, C.-H., Park, C.-E., Storelvmo, T., and Tan, I.: Influence of cloud phase composition on climate feedbacks, J. Geophys. Res.-Atmos., 119, 3687–3700, 2014. a

Crameri, F.: Scientific colour-maps, Zenodo, https://doi.org/10.5281/zenodo.1243862, 2018. a

Crewell, S. and Löhnert, U.: Accuracy of Boundary Layer Temperature Profiles Retrieved With Multifrequency Multiangle Microwave Radiometry, IEEE T. Geosci. Remote, 45, 2195–2201, https://doi.org/10.1109/TGRS.2006.888434, 2007. a

Dahlke, S. and Maturilli, M.: Contribution of Atmospheric Advection to the Amplified Winter Warming in the Arctic North Atlantic Region, Adv. Meteorol., 2017, 4928620, https://doi.org/10.1155/2017/4928620, 2017. a, b

De Boer, G., Eloranta, E. W., and Shupe, M. D.: Arctic mixed-phase stratiform cloud properties from multiple years of surface-based measurements at two high-latitude locations, J. Atmos. Sci., 66, 2874–2887, 2009. a, b

Dekhtyareva, A., Holmén, K., Maturilli, M., Hermansen, O., and Graversen, R.: Effect of seasonal mesoscale and microscale meteorological conditions in Ny-Ålesund on results of monitoring of long-range transported pollution, Polar. Res., 37, 1508196, https://doi.org/10.1080/17518369.2018.1508196, 2018. a

Dong, X., Xi, B., Crosby, K., Long, C. N., Stone, R. S., and Shupe, M. D.: A 10 year climatology of Arctic cloud fraction and radiative forcing at Barrow, Alaska, J. Geophys. Res., 115, D17212, https://doi.org/10.1029/2009JD013489, 2010. a

Ebell, K., Nomokonova, T., Maturilli, M., and Ritter, C.: Radiative effect of clouds at Ny-Ålesund, Svalbard, as inferred from ground-based remote sensing observations, J. Appl. Meteorol. Clim., 59, 3–22, https://doi.org/10.1175/JAMC-D-19-0080.1, 2019. a

Eirund, G. K., Possner, A., and Lohmann, U.: Response of Arctic mixed-phase clouds to aerosol perturbations under different surface forcings, Atmos. Chem. Phys., 19, 9847–9864, https://doi.org/10.5194/acp-19-9847-2019, 2019. a

Esau, I. and Repina, I.: Wind climate in Kongsfjorden, Svalbard, and attribution of leading wind driving mechanisms through turbulence-resolving simulations, Adv. Meteorol., 2012, 568454, https://doi.org/10.1155/2012/568454, 2012. a, b, c

Hogan, R. J. and O'Connor, E. J.: Facilitating cloud radar and lidar algorithms: the Cloudnet Instrument Synergy/Target Categorization product, Cloudnet documentation, 1–14, 2004. a, b

Hogan, R. J., Jakob, C., and Illingworth, A. J.: Comparison of ECMWF winter-season cloud fraction with radar-derived values, J. Appl. Meteorol., 40, 513–525, 2001. a

Hogan, R. J., Mittermaier, M. P., and Illingworth, A. J.: The retrieval of ice water content from radar reflectivity factor and temperature and its use in evaluating a mesoscale model, J. Appl. Meteorol. Clim., 45, 301–317, 2006. a

Illingworth, A., Hogan, R., O'connor, E., Bouniol, D., Brooks, M., Delanoë, J., Donovan, D., Eastment, J., Gaussiat, N., Goddard, J., Haeffelin, M., Baltink, H. K., Krasnov, O. A., Pelon, J., Piriou, J., Protat, A., Russchenberg, H. W., Seifert, A., Tompkins, A. M., van Zadelhoff, G., Vinit, F., Willén, U., Wilson, D. R., and Wrench, C. L.: Cloudnet: Continuous evaluation of cloud profiles in seven operational models using ground-based observations, B. Am. Meteorol. Soc., 88, 883–898, 2007. a

Intrieri, J., Fairall, C., Shupe, M., Persson, P., Andreas, E., Guest, P., and Moritz, R.: An annual cycle of Arctic surface cloud forcing at SHEBA, J. Geophys. Res., 107, 8039, https://doi.org/10.1029/2000JC000439, 2002. a

Jackson, R. C., McFarquhar, G. M., Korolev, A. V., Earle, M. E., Liu, P. S., Lawson, R. P., Brooks, S., Wolde, M., Laskin, A., and Freer, M.: The dependence of ice microphysics on aerosol concentration in Arctic mixed-phase stratus clouds during ISDAC and M-PACE, J. Geophys. Res., 117, D15207, https://doi.org/10.1029/2012JD017668, 2012. a

Jenkinson, A. and Collison, F.: An initial climatology of gales over the North Sea, Synoptic Climatology Branch Memorandum, 62, 1–18, 1977. a

Jocher, G., Karner, F., Ritter, C., Neuber, R., Dethloff, K., Obleitner, F., Reuder, J., and Foken, T.: The near-surface small-scale spatial and temporal variability of sensible and latent heat exchange in the Svalbard region: a case study, ISRN Meteorology, 2012. a, b, c, d, e

Kalesse, H., de Boer, G., Solomon, A., Oue, M., Ahlgrimm, M., Zhang, D., Shupe, M. D., Luke, E., and Protat, A.: Understanding Rapid Changes in Phase Partitioning between Cloud Liquid and Ice in Stratiform Mixed-Phase Clouds: An Arctic Case Study, Mon. Weather Rev., 144, 4805–4826, 2016. a

Kay, J. E. and Gettelman, A.: Cloud influence on and response to seasonal Arctic sea ice loss, J. Geophys. Res., 114, D18204, https://doi.org/10.1029/2009JD011773, 2009. a

Kay, J. E., L'Ecuyer, T., Chepfer, H., Loeb, N., Morrison, A., and Cesana, G.: Recent Advances in Arctic Cloud and Climate Research, Current Climate Change Reports, 2, 159–169, 2016. a

Kayser, M., Maturilli, M., Graham, R. M., Hudson, S. R., Rinke, A., Cohen, L., Kim, J.-H., Park, S.-J., Moon, W., and Granskog, M. A.: Vertical thermodynamic structure of the troposphere during the Norwegian young sea ICE expedition (N-ICE2015), J. Geophys. Res.-Atmos., 122, 10855–10872, 2017. a

Keine, C.: Moods Median Test, GitHub, available at: https://www.github.com/ChristianKeine/Moods-Mediantest, last access: 1 February 2019. a

Kim, B.-M., Hong, J.-Y., Jun, S.-Y., Zhang, X., Kwon, H., Kim, S.-J., Kim, J.-H., Kim, S.-W., and Kim, H.-K.: Major cause of unprecedented Arctic warming in January 2016: Critical role of an Atlantic windstorm, Scientific Reports, 7, 40051, https://doi.org/10.1038/srep40051, 2017. a

Klein, S. A., McCoy, R. B., Morrison, H., Ackerman, A. S., Avramov, A., de Boer, G., Chen, M., Cole, J. N., Del Genio, A. D., Falk, M., Foster, M. J., Fridlind, A., Golaz, J.‐C., Hashino, T., Harrington, J. Y., Hoose, C., Khairoutdinov, M. F., Larson, V. E., Liu, X., Luo, Y., McFarquhar, G. M., Menon, S., Neggers, R. A. J., Park, S., Poellot, M. R., Schmidt, J. M., Sednev, I., Shipway, B. J., Shupe, M. D., Spangenberg, D.A., Sud, Y. C., Turner, D. D., Veron, D. E., von Salzen, K., Walker, G. K., Wang, Z., Wolf, A. B., Xie, S., Xu, K.‐M., Yang, F., and Zhang, G.: Intercomparison of model simulations of mixed-phase clouds observed during the ARM Mixed-Phase Arctic Cloud Experiment. I: Single-layer cloud, Q. J. Roy. Meteor. Soc., 135, 979–1002, 2009. a

Knudsen, E. M., Heinold, B., Dahlke, S., Bozem, H., Crewell, S., Gorodetskaya, I. V., Heygster, G., Kunkel, D., Maturilli, M., Mech, M., Viceto, C., Rinke, A., Schmithüsen, H., Ehrlich, A., Macke, A., Lÿpkes, C., and Wendisch, M.: Meteorological conditions during the ACLOUD/PASCAL field campaign near Svalbard in early summer 2017, Atmos. Chem. Phys., 18, 17995–18022, https://doi.org/10.5194/acp-18-17995-2018, 2018. a

Komurcu, M., Storelvmo, T., Tan, I., Lohmann, U., Yun, Y., Penner, J. E., Wang, Y., Liu, X., and Takemura, T.: Intercomparison of the cloud water phase among global climate models, J. Geophys. Res.-Atmos., 119, 3372–3400, 2014. a

Korolev, A. and Isaac, G.: Phase transformation of mixed-phase clouds, Q. J. Roy. Meteor. Soc., 129, 19–38, 2003. a

Küchler, N., Kneifel, S., Löhnert, U., Kollias, P., Czekala, H., and Rose, T.: A W-Band Radar–Radiometer System for Accurate and Continuous Monitoring of Clouds and Precipitation, J. Atmos. Ocean. Tech., 34, 2375–2392, 2017. a

Lamb, D. and Verlinde, J.: Physics and chemistry of clouds, Cambridge University Press, UK, 2011. a

Li, Z., Xu, K.-M., and Cheng, A.: The Response of Simulated Arctic Mixed-Phase Stratocumulus to Sea Ice Cover Variability in the Absence of Large-Scale Advection, J. Geophys. Res.-Atmos., 122, 12335–12352, 2017. a

Löhnert, U. and Crewell, S.: Accuracy of cloud liquid water path from ground-based microwave radiometry 1. Dependency on cloud model statistics, Radio Sci., 38, 8041, https://doi.org/10.1029/2002RS002654, 2003. a, b

Maturilli, M. and Ebell, K.: Twenty-five years of cloud base height measurements by ceilometer in Ny-Ålesund, Svalbard, Earth Syst. Sci. Data, 10, 1451–1456, https://doi.org/10.5194/essd-10-1451-2018, 2018. a, b, c

Maturilli, M. and Kayser, M.: Homogenized radiosonde record at station Ny-Ålesund, Spitsbergen, 1993–2014, PANGAEA, https://doi.org/10.1594/PANGAEA.845373, 2016. a

Maturilli, M. and Kayser, M.: Arctic warming, moisture increase and circulation changes observed in the Ny-Ålesund homogenized radiosonde record, Theor. Appl. Climatol., 130, 1–17, 2017a. a, b, c, d, e, f, g

Maturilli, M. and Kayser, M.: Homogenized radiosonde record at station Ny-Ålesund, Spitsbergen, 2015–2016, PANGAEA, https://doi.org/10.1594/PANGAEA.875196, 2017b. a

Maturilli, M., Herber, A., and König-Langlo, G.: Climatology and time series of surface meteorology in Ny-Ålesund, Svalbard, Earth Syst. Sci. Data, 5, 155–163, https://doi.org/10.5194/essd-5-155-2013, 2013. a

McCoy, D. T., Tan, I., Hartmann, D. L., Zelinka, M. D., and Storelvmo, T.: On the relationships among cloud cover, mixed-phase partitioning, and planetary albedo in GCMs, J. Adv. Model. Earth Sy., 8, 650–668, 2016. a

Mech, M., Kliesch, L.-L., Anhäuser, A., Rose, T., Kollias, P., and Crewell, S.: Microwave Radar/radiometer for Arctic Clouds (MiRAC): first insights from the ACLOUD campaign, Atmos. Meas. Tech., 12, 5019–5037, https://doi.org/10.5194/amt-12-5019-2019, 2019. a

Mioche, G., Jourdan, O., Ceccaldi, M., and Delanoë, J.: Variability of mixed-phase clouds in the Arctic with a focus on the Svalbard region: a study based on spaceborne active remote sensing, Atmos. Chem. Phys., 15, 2445–2461, https://doi.org/10.5194/acp-15-2445-2015, 2015. a, b, c, d, e, f

Mioche, G., Jourdan, O., Delanoë, J., Gourbeyre, C., Febvre, G., Dupuy, R., Monier, M., Szczap, F., Schwarzenboeck, A., and Gayet, J.-F.: Vertical distribution of microphysical properties of Arctic springtime low-level mixed-phase clouds over the Greenland and Norwegian seas, Atmos. Chem. Phys., 17, 12845–12869, https://doi.org/10.5194/acp-17-12845-2017, 2017. a, b, c

Morrison, A., Kay, J., Chepfer, H., Guzman, R., and Yettella, V.: Isolating the liquid cloud response to recent Arctic sea ice variability using spaceborne lidar observations, J. Geophys. Res.-Atmos., 123, 473–490, 2018. a

Morrison, H., Shupe, M. D., Pinto, J. O., and Curry, J. A.: Possible roles of ice nucleation mode and ice nuclei depletion in the extended lifetime of Arctic mixed-phase clouds, Geophys. Res. Lett., 32, L18801, https://doi.org/10.1029/2005GL023614, 2005. a

Morrison, H., Pinto, J. O., Curry, J. A., and McFarquhar, G. M.: Sensitivity of modeled Arctic mixed-phase stratocumulus to cloud condensation and ice nuclei over regionally varying surface conditions, J. Geophys. Res., 113, D05203, https://doi.org/10.1029/2007JD008729, 2008. a, b, c

Morrison, H., De Boer, G., Feingold, G., Harrington, J., Shupe, M. D., and Sulia, K.: Resilience of persistent Arctic mixed-phase clouds, Nat. Geosci., 5, 11–17, 2012. a, b, c

Nomokonova, T. and Ebell, K.: Cloud microphysical properties retrieved from ground-based remote sensing at Ny-Ålesund (10 June 2016–8 October 2018), University of Cologne, PANGAEA, https://doi.org/10.1594/PANGAEA.898556, 2019. a

Nomokonova, T., Ebell, K., Löhnert, U., Maturilli, M., and Ritter, C.: The influence of anomalous atmospheric conditions at Ny-Ålesund on clouds and their radiative effect, Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2019-985, in review, 2019a. a

Nomokonova, T., Ebell, K., Löhnert, U., Maturilli, M., Ritter, C., and O'Connor, E.: Statistics on clouds and their relation to thermodynamic conditions at Ny-Ålesund using ground-based sensor synergy, Atmos. Chem. Phys., 19, 4105–4126, https://doi.org/10.5194/acp-19-4105-2019, 2019b. a, b, c, d, e, f

Nomokonova, T., Ritter, C., and Ebell, K.: HATPRO microwave radiometer measurements at AWIPEV, Ny-Ålesund (2016–2018), PANGAEA, https://doi.org/10.1594/PANGAEA.902183, 2019c. a

Norgren, M. S., de Boer, G., and Shupe, M. D.: Observed aerosol suppression of cloud ice in low-level Arctic mixed-phase clouds, Atmos. Chem. Phys., 18, 13345–13361, https://doi.org/10.5194/acp-18-13345-2018, 2018. a

Norwegian Polar Institute: Kartdata Svalbard 1:100 000 (S100 Kartdata)/Map Data, Norwegian Polar Institute, Data set, https://doi.org/10.21334/npolar.2014.645336c7, 2014. a, b

Pithan, F., Svensson, G., Caballero, R., Chechin, D., Cronin, T. W., Ekman, A. M., Neggers, R., Shupe, M. D., Solomon, A., Tjernström, M., and Wendisch, M.: Role of air-mass transformations in exchange between the Arctic and mid-latitudes, Nat. Geosci., 11, 805–812, 2018. a

Qiu, S., Xi, B., and Dong, X.: Influence of Wind Direction on Thermodynamic Properties and Arctic Mixed-Phase Clouds in Autumn at Utqiaġvik, Alaska, J. Geophys. Res.-Atmos., 123, 9589–9603, 2018. a

Quinn, P., Shaw, G., Andrews, E., Dutton, E., Ruoho-Airola, T., and Gong, S.: Arctic haze: current trends and knowledge gaps, Tellus B, 59, 99–114, 2007. a

Rose, T., Crewell, S., Löhnert, U., and Simmer, C.: A network suitable microwave radiometer for operational monitoring of the cloudy atmosphere, Atmos. Res., 75, 183–200, 2005. a

Savre, J., Ekman, A. M., Svensson, G., and Tjernström, M.: Large-eddy simulations of an Arctic mixed-phase stratiform cloud observed during ISDAC: sensitivity to moisture aloft, surface fluxes and large-scale forcing, Q. J. Roy. Meteor. Soc., 141, 1177–1190, 2015. a

Scott, R. C. and Lubin, D.: Unique manifestations of mixed-phase cloud microphysics over Ross Island and the Ross Ice Shelf, Antarctica, Geophys. Res. Lett., 43, 2936–2945, 2016. a

Sedlar, J., Shupe, M. D., and Tjernström, M.: On the relationship between thermodynamic structure and cloud top, and its climate significance in the Arctic, J. Climate, 25, 2374–2393, 2012. a

Serreze, M. C., Barrett, A. P., Stroeve, J. C., Kindig, D. N., and Holland, M. M.: The emergence of surface-based Arctic amplification, The Cryosphere, 3, 11–19, https://doi.org/10.5194/tc-3-11-2009, 2009. a

Serreze, M. C. and Barry, R. G.: Processes and impacts of Arctic amplification: A research synthesis, Global Planet. Change, 77, 85–96, 2011. a

Sheskin, D.: Handbook of Parametric and Nonparametric Statistical Procedures, 2nd edn., CRC Press, USA, 2000. a

Shupe, M. D.: Clouds at Arctic atmospheric observatories. Part II: Thermodynamic phase characteristics, J. Appl. Meteorol. Clim, 50, 645–661, 2011. a, b

Shupe, M. D. and Intrieri, J. M.: Cloud radiative forcing of the Arctic surface: The influence of cloud properties, surface albedo, and solar zenith angle, J. Climate, 17, 616–628, 2004. a, b, c

Shupe, M. D., Kollias, P., Matrosov, S. Y., and Schneider, T. L.: Deriving mixed-phase cloud properties from Doppler radar spectra, J. Atmos. Ocean. Tech., 21, 660–670, 2004. a

Shupe, M. D., Matrosov, S. Y., and Uttal, T.: Arctic mixed-phase cloud properties derived from surface-based sensors at SHEBA, J. Atmos. Sci., 63, 697–711, 2006. a, b, c

Shupe, M. D., Kollias, P., Persson, P. O. G., and McFarquhar, G. M.: Vertical motions in Arctic mixed-phase stratiform clouds, J. Atmos. Sci., 65, 1304–1322, 2008. a

Shupe, M. D., Walden, V. P., Eloranta, E., Uttal, T., Campbell, J. R., Starkweather, S. M., and Shiobara, M.: Clouds at Arctic atmospheric observatories. Part I: Occurrence and macrophysical properties, J. Appl. Meteorol. Clim, 50, 626–644, 2011. a, b, c, d

Shupe, M. D., Persson, P. O. G., Brooks, I. M., Tjernström, M., Sedlar, J., Mauritsen, T., Sjogren, S., and Leck, C.: Cloud and boundary layer interactions over the Arctic sea ice in late summer, Atmos. Chem. Phys., 13, 9379–9399, https://doi.org/10.5194/acp-13-9379-2013, 2013. a, b, c, d

Solomon, A., Shupe, M. D., Persson, O., Morrison, H., Yamaguchi, T., Caldwell, P. M., and de Boer, G.: The sensitivity of springtime Arctic mixed-phase stratocumulus clouds to surface-layer and cloud-top inversion-layer moisture sources, J. Atmos. Sci., 71, 574–595, 2014. a

Solomon, A., de Boer, G., Creamean, J. M., McComiskey, A., Shupe, M. D., Maahn, M., and Cox, C.: The relative impact of cloud condensation nuclei and ice nucleating particle concentrations on phase partitioning in Arctic mixed-phase stratocumulus clouds, Atmos. Chem. Phys., 18, 17047–17059, https://doi.org/10.5194/acp-18-17047-2018, 2018. a

Solomon, S., Qin, D., Manning, M., Averyt, K., and Marquis, M.: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, vol. 4, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge University Press, Cambridge, UK and New York, NY, USA, 2007. a

Sotiropoulou, G., Sedlar, J., Tjernström, M., Shupe, M. D., Brooks, I. M., and Persson, P. O. G.: The thermodynamic structure of summer Arctic stratocumulus and the dynamic coupling to the surface, Atmos. Chem. Phys., 14, 12573–12592, https://doi.org/10.5194/acp-14-12573-2014, 2014. a, b, c, d

Svendsen, H., Beszczynska-Møller, A., Hagen, J. O., Lefauconnier, B., Tverberg, V., Gerland, S., Børre Ørbæk, J., Bischof, K., Papucci, C., Zajaczkowski, M., Azzolini, R., Bruland, O., and Wiencke, C.: The physical environment of Kongsfjorden–Krossfjorden, an Arctic fjord system in Svalbard, Polar Res., 21, 133–166, 2002. a, b

Wendisch, M., Macke, A., Ehrlich, A., Lüpkes, C., Mech, M., Chechin, D., Dethloff, K., Barientos, C., Bozem, H., Brückner, M., Clemen, H., Crewell, S., Donth, T., Dupuy, R., Ebell, K., Egerer, U., Engelmann, R., Engler, C., Eppers, O., Gehrmann, M., Gong, X., Gottschalk, M., Gourbeyre, C., Griesche, H., Hartmann, J., Hartmann, M., Heinold, B., Herber, A., Herrmann, H., Heygster, G., Hoor, P., Jafariserajehlou, S., Jäkel, E. Jërvinen, E., Jourdan, O., Kästner, U., Kecorius, S., Knudsen, E. M., Köllner, F., Kretzschmar, J., Lelli, L., Leroy, D., Maturilli, M., Mei, L., Mertes, S., Mioche, G., Neuber, R., Nicolaus, M., Nomokonova, T., Notholt, J., Palm, M., van Pinxteren, M., Quaas, J., Richter, P., Ruiz-Donoso, E., Schäfer, M., Schmieder, K., Schnaiter, M., Schneider, J., Schwarzenböck, A., Seifert, P., Shupe, M. D., Siebert, H., Spreen, G., Stapf, J., Stratmann, F., Vogl, T., Welti, A., Wex, H., Wiedensohler, A., Zanatta, M., and Zeppenfeld, S.: The Arctic Cloud Puzzle: Using ACLOUD/PASCAL Multiplatform Observations to Unravel the Role of Clouds and Aerosol Particles in Arctic Amplification, B. Am. Meteorol. Soc., 100, 841–871, 2019. a

Wendish, M., Brückner, M., Burrows, J. P., Crewell, S., Dethloff, K., Ebell, K., Lüpkes, C., Macke, A., Notholt, J., Quaas, J., Rinke, A., and Tegen, I.: Understanding Causes and Effects of Rapid Warming in the Arctic, EOS, 98, https://doi.org/10.1029/2017EO064803, 2017. a

Westbrook, C. D. and Illingworth, A. J.: Evidence that ice forms primarily in supercooled liquid clouds at temperatures $>-\mathrm{27}$C, Geophys. Res. Lett., 38, L14808, https://doi.org/10.1029/2011GL048021, 2011. a

Woods, C., Caballero, R., and Svensson, G.: Large-scale circulation associated with moisture intrusions into the Arctic during winter, Geophys. Res. Lett., 40, 4717–4721, 2013. a

Yamartino, R. J.: A comparison of several “single-pass” estimators of the standard deviation of wind direction, J. Clim. Appl. Meteorol., 23, 1362–1366, 1984. a

Young, G., Jones, H. M., Choularton, T. W., Crosier, J., Bower, K. N., Gallagher, M. W., Davies, R. S., Renfrew, I. A., Elvidge, A. D., Darbyshire, E., Marenco, F., Brown, P. R. A., Ricketts, H. M. A., Connolly, P. J., Lloyd, G., Williams, P. I., Allan, J. D., Taylor, J. W., Liu, D., and Flynn, M. J.: Observed microphysical changes in Arctic mixed-phase clouds when transitioning from sea ice to open ocean, Atmos. Chem. Phys., 16, 13945–13967, https://doi.org/10.5194/acp-16-13945-2016, 2016. a

Zängl, G., Reinert, D., Rípodas, P., and Baldauf, M.: The ICON (ICOsahedral Non-hydrostatic) modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core, Q. J. Roy. Meteor. Soc., 141, 563–579, 2015. a

Zhao, M. and Wang, Z.: Comparison of Arctic clouds between European Center for Medium-Range Weather Forecasts simulations and Atmospheric Radiation Measurement Climate Research Facility long-term observations at the North Slope of Alaska Barrow site, J. Geophys. Res., 115, D23202, https://doi.org/10.1029/2010JD014285, 2010. a, b

Zuidema, P., Baker, B., Han, Y., Intrieri, J., Key, J., Lawson, P., Matrosov, S., Shupe, M., Stone, R., and Uttal, T.: An Arctic springtime mixed-phase cloudy boundary layer observed during SHEBA, J. Atmos. Sci., 62, 160–176, 2005. a