ALBERT

Description: Based on results of 11 yr of heterogeneous ice nucleation experiments at the Aerosol Interaction and Dynamics in the Atmosphere (AIDA) chamber in Karlsruhe, Germany, a new empirical parameterization framework for heterogeneous ice nucleation was developed. The framework currently includes desert dust and soot aerosol and quantifies the ice nucleation efficiency in terms of the ice nucleation active surface site (INAS) approach. The immersion freezing INAS densities nS of all desert dust experiments follow an exponential fit as a function of temperature, well in agreement with an earlier analysis of AIDA experiments. The deposition nucleation nS isolines for desert dust follow u-shaped curves in the ice saturation ratio–temperature (Si–T) diagram at temperatures below about 240 K. The negative slope of these isolines toward lower temperatures may be explained by classical nucleation theory (CNT), whereas the behavior toward higher temperatures may be caused by a pore condensation and freezing mechanism. The deposition nucleation measured for soot at temperatures below about 240 K also follows u-shaped isolines with a shift toward higher Si for soot with higher organic carbon content. For immersion freezing of soot aerosol, only upper limits for nS were determined and used to rescale an existing parameterization line. The new parameterization framework is compared to a CNT-based parameterization and an empirical framework as used in models. The comparison shows large differences in shape and magnitude of the nS isolines especially for deposition nucleation. For the application in models, implementation of this new framework is simple compared to that of other expressions.

Description: The response of clouds to changes in the aerosol concentration is complex and may differ depending on the cloud type, the aerosol regime, and environmental conditions. In this study, a novel technique is used to systematically modify the environmental conditions in realistic convection-resolving simulations for cases with weak and strong large-scale forcing over central Europe with the Consortium for Small-Scale Modeling (COSMO) model. Besides control runs with quasi-operational settings, initial and boundary temperature profiles are modified with linear increasing temperature increments from 0 to 5 K between 3 and 12 km AGL to represent different amounts of convective available potential energy (CAPE) and relative humidity. The results show a systematic decrease of total precipitation with increasing cloud condensation nuclei (CCN) concentrations for the cases with strong synoptic forcing caused by a suppressed warm-rain process, whereas no systematic aerosol effect is simulated for weak synoptic forcing. The effect of increasing CCN tends to be stronger in the simulations with increased temperatures and lower CAPE. While the large-scale domain-averaged responses to increased CCN are weak, the precipitation forming over mountainous terrain reveals a stronger sensitivity for most of the analyzed cases. Our findings also demonstrate that the role of the warm-rain process is more important for strong than for weak synoptic forcing. The aerosol effect is largest for weakly forced conditions but more predictable for the strongly forced cases. However, more accurate environmental conditions are much more important than accurate aerosol assumptions, especially for weak large-scale forcing.

Description: Clouds and aerosols contribute the largest uncertainty to current estimates and interpretations of the Earth’s changing energy budget. Here we use a new-generation large-domain large-eddy model, ICON-LEM (ICOsahedral Non-hydrostatic Large Eddy Model), to simulate the response of clouds to realistic anthropogenic perturbations in aerosols serving as cloud condensation nuclei (CCN). The novelty compared to previous studies is that (i) the LEM is run in weather prediction mode and with fully interactive land surface over a large domain and (ii) a large range of data from various sources are used for the detection and attribution. The aerosol perturbation was chosen as peak-aerosol conditions over Europe in 1985, with more than fivefold more sulfate than in 2013. Observational data from various satellite and ground-based remote sensing instruments are used, aiming at the detection and attribution of this response. The simulation was run for a selected day (2 May 2013) in which a large variety of cloud regimes was present over the selected domain of central Europe. It is first demonstrated that the aerosol fields used in the model are consistent with corresponding satellite aerosol optical depth retrievals for both 1985 (perturbed) and 2013 (reference) conditions. In comparison to retrievals from ground-based lidar for 2013, CCN profiles for the reference conditions were consistent with the observations, while the ones for the 1985 conditions were not. Similarly, the detection and attribution process was successful for droplet number concentrations: the ones simulated for the 2013 conditions were consistent with satellite as well as new ground-based lidar retrievals, while the ones for the 1985 conditions were outside the observational range. For other cloud quantities, including cloud fraction, liquid water path, cloud base altitude and cloud lifetime, the aerosol response was small compared to their natural variability. Also, large uncertainties in satellite and ground-based observations make the detection and attribution difficult for these quantities. An exception to this is the fact that at a large liquid water path value (LWP 〉 200 g m−2), the control simulation matches the observations, while the perturbed one shows an LWP which is too large. The model simulations allowed for quantifying the radiative forcing due to aerosol–cloud interactions, as well as the adjustments to this forcing. The latter were small compared to the variability and showed overall a small positive radiative effect. The overall effective radiative forcing (ERF) due to aerosol–cloud interactions (ERFaci) in the simulation was dominated thus by the Twomey effect and yielded for this day, region and aerosol perturbation −2.6 W m−2. Using general circulation models to scale this to a global-mean present-day vs. pre-industrial ERFaci yields a global ERFaci of −0.8 W m−2.

Description: Secondary ice production via processes like rime splintering, frozen droplet shattering, and breakup upon ice hydrometeor collision have been proposed to explain discrepancies between in-cloud ice crystal and ice-nucleating particle numbers. To understand the impact of this additional ice crystal generation on surface precipitation, we present one of the first studies to implement frozen droplet shattering and ice–ice collisional breakup parameterizations in a mesoscale model. We simulate a cold frontal rainband from the Aerosol Properties, PRocesses, And InfluenceS on the Earth's Climate campaign and investigate the impact of the new parameterizations on the simulated ice crystal number concentrations (ICNC) and precipitation. Near the convective regions of the rainband, contributions to ICNC can be as large from secondary production as from primary nucleation, but ICNCs greater than 50 L−1 remain underestimated by the model. The addition of the secondary production parameterizations also clearly intensifies the differences in both accumulated precipitation and precipitation rate between the convective towers and non-convective gap regions. We suggest, then, that secondary ice production parameterizations be included in large-scale models on the basis of large hydrometeor concentration and convective activity criteria.

Description: Multilayer clouds (MLCs) occur more often in the Arctic than globally. In this study we present the results of a detection algorithm applied to radiosonde and radar data from an 1-year time period in Ny-Ålesund, Svalbard. Multilayer cloud occurrence is found on 29 % of the investigated days. These multilayer cloud cases are further analysed regarding the possibility of ice crystal seeding, meaning that an ice crystal can survive sublimation in a subsaturated layer between two cloud layers when falling through this layer. For this we analyse profiles of relative humidity with respect to ice to identify super- and subsaturated air layers. Then the sublimation of an ice crystal of an assumed initial size of r=400 µm on its way through the subsaturated layer is calculated. If the ice crystal still exists when reaching a lower supersaturated layer, ice crystal seeding can potentially take place. Seeding cases are found often, in 23 % of the investigated days (100 % includes all days, as well as non-cloudy days). The identification of seeding cases is limited by the radar signal inside the subsaturated layer. Clearly separated multilayer clouds, defined by a clear interstice in the radar image, do not interact through seeding (9 % of the investigated days). There are various deviations between the relative humidity profiles and the radar images, e.g. due to the lack of ice-nucleating particles (INPs) and cloud condensation nuclei (CCN). Additionally, horizontal wind drift of the radiosonde and time restriction when comparing radiosonde and radar data cause further deviations. In order to account for some of these deviations, an evaluation by manual visual inspection is done for the non-seeding cases.

Description: The Arctic climate is changing; temperature changes in the Arctic are greater than at midlatitudes, and changing atmospheric conditions influence Arctic mixed-phase clouds, which are important for the Arctic surface energy budget. These low-level clouds are frequently observed across the Arctic. They impact the turbulent and radiative heating of the open water, snow, and sea-ice-covered surfaces and influence the boundary layer structure. Therefore the processes that affect mixed-phase cloud life cycles are extremely important, yet relatively poorly understood. In this study, we present sensitivity studies using semi-idealized large eddy simulations (LESs) to identify processes contributing to the dissipation of Arctic mixed-phase clouds. We found that one potential main contributor to the dissipation of an observed Arctic mixed-phase cloud, during the Arctic Summer Cloud Ocean Study (ASCOS) field campaign, was a low cloud droplet number concentration (CDNC) of about 2 cm−3. Introducing a high ice crystal concentration of 10 L−1 also resulted in cloud dissipation, but such high ice crystal concentrations were deemed unlikely for the present case. Sensitivity studies simulating the advection of dry air above the boundary layer inversion, as well as a modest increase in ice crystal concentration of 1 L−1, did not lead to cloud dissipation. As a requirement for small droplet numbers, pristine aerosol conditions in the Arctic environment are therefore considered an important factor determining the lifetime of Arctic mixed-phase clouds.

Description: A regional modeling study on the impact of desert dust on cloud formation is presented for a major Saharan dust outbreak over Europe from 2 to 5 April 2014. The dust event coincided with an extensive and dense cirrus cloud layer, suggesting an influence of dust on atmospheric ice nucleation. Using interactive simulation with the regional dust model COSMO-MUSCAT, we investigate cloud and precipitation representation in the model and test the sensitivity of cloud parameters to dust–cloud and dust–radiation interactions of the simulated dust plume. We evaluate model results with ground-based and spaceborne remote sensing measurements of aerosol and cloud properties, as well as the in situ measurements obtained during the ML-CIRRUS aircraft campaign. A run of the model with single-moment bulk microphysics without online dust feedback considerably underestimated cirrus cloud cover over Germany in the comparison with infrared satellite imagery. This was also reflected in simulated upper-tropospheric ice water content (IWC), which accounted for only 20 % of the observed values. The interactive dust simulation with COSMO-MUSCAT, including a two-moment bulk microphysics scheme and dust–cloud as well as dust–radiation feedback, in contrast, led to significant improvements. The modeled cirrus cloud cover and IWC were by at least a factor of 2 higher in the relevant altitudes compared to the noninteractive model run. We attributed these improvements mainly to enhanced deposition freezing in response to the high mineral dust concentrations. This was corroborated further in a significant decrease in ice particle radii towards more realistic values, compared to in situ measurements from the ML-CIRRUS aircraft campaign. By testing different empirical ice nucleation parameterizations, we further demonstrate that remaining uncertainties in the ice-nucleating properties of mineral dust affect the model performance at least as significantly as including the online representation of the mineral dust distribution. Dust–radiation interactions played a secondary role for cirrus cloud formation, but contributed to a more realistic representation of precipitation by suppressing moist convection in southern Germany. In addition, a too-low specific humidity in the 7 to 10 km altitude range in the boundary conditions was identified as one of the main reasons for misrepresentation of cirrus clouds in this model study.

Description: We perform a model intercomparison of summertime high Arctic (〉 80∘ N) clouds observed during the 2008 Arctic Summer Cloud Ocean Study (ASCOS) campaign, when observed cloud condensation nuclei (CCN) concentrations fell below 1 cm−3. Previous analyses have suggested that at these low CCN concentrations the liquid water content (LWC) and radiative properties of the clouds are determined primarily by the CCN concentrations, conditions that have previously been referred to as the tenuous cloud regime. The intercomparison includes results from three large eddy simulation models (UCLALES-SALSA, COSMO-LES, and MIMICA) and three numerical weather prediction models (COSMO-NWP, WRF, and UM-CASIM). We test the sensitivities of the model results to different treatments of cloud droplet activation, including prescribed cloud droplet number concentrations (CDNCs) and diagnostic CCN activation based on either fixed aerosol concentrations or prognostic aerosol with in-cloud processing. There remains considerable diversity even in experiments with prescribed CDNCs and prescribed ice crystal number concentrations (ICNC). The sensitivity of mixed-phase Arctic cloud properties to changes in CDNC depends on the representation of the cloud droplet size distribution within each model, which impacts autoconversion rates. Our results therefore suggest that properly estimating aerosol–cloud interactions requires an appropriate treatment of the cloud droplet size distribution within models, as well as in situ observations of hydrometeor size distributions to constrain them. The results strongly support the hypothesis that the liquid water content of these clouds is CCN limited. For the observed meteorological conditions, the cloud generally did not collapse when the CCN concentration was held constant at the relatively high CCN concentrations measured during the cloudy period, but the cloud thins or collapses as the CCN concentration is reduced. The CCN concentration at which collapse occurs varies substantially between models. Only one model predicts complete dissipation of the cloud due to glaciation, and this occurs only for the largest prescribed ICNC tested in this study. Global and regional models with either prescribed CDNCs or prescribed aerosol concentrations would not reproduce these dissipation events. Additionally, future increases in Arctic aerosol concentrations would be expected to decrease the frequency of occurrence of such cloud dissipation events, with implications for the radiative balance at the surface. Our results also show that cooling of the sea-ice surface following cloud dissipation increases atmospheric stability near the surface, further suppressing cloud formation. Therefore, this suggests that linkages between aerosol and clouds, as well as linkages between clouds, surface temperatures, and atmospheric stability need to be considered for weather and climate predictions in this region.

Description: Two-dimensional horizontal fields of cloud optical thickness τ derived from airborne measurements of solar spectral, cloud-reflected radiance are compared with semi-idealized large eddy simulations (LESs) of Arctic stratus performed with the Consortium for Small-scale Modeling (COSMO) atmospheric model. The measurements were collected during the Vertical Distribution of Ice in Arctic Clouds (VERDI) campaign carried out in Inuvik, Canada, in April/May 2012. The input for the LESs is obtained from collocated airborne dropsonde observations of a persistent Arctic stratus above the sea-ice-free Beaufort Sea. Simulations are performed for spatial resolutions of 50 m (1.6 km × 1.6 km domain) and 100 m (6.4 km × 6.4 km domain). Macrophysical cloud properties, such as cloud top altitude and vertical extent, are well captured by the COSMO simulations. However, COSMO produces rather homogeneous clouds compared to the measurements, in particular for the simulations with coarser spatial resolution. For both spatial resolutions, the directional structure of the cloud inhomogeneity is well represented by the model. Differences between the individual cases are mainly associated with the wind shear near cloud top and the vertical structure of the atmospheric boundary layer. A sensitivity study changing the wind velocity in COSMO by a vertically constant scaling factor shows that the directional, small-scale cloud inhomogeneity structures can range from 250 to 800 m, depending on the mean wind speed, if the simulated domain is large enough to capture also large-scale structures, which then influence the small-scale structures. For those cases, a threshold wind velocity is identified, which determines when the cloud inhomogeneity stops increasing with increasing wind velocity.