Cloud-base vertical velocity statistics: a comparison between an atmospheric mesoscale model and remote sensing observations

The statistics of cloud-base vertical velocity simulated by the non-hydrostatic mesoscale model AROME are compared with Cloudnet remote sensing observations at two locations: the ARM SGP site in Central Oklahoma, and the DWD observatory at Lindenberg, Germany. The results show that, as expected, AROME significantly underestimates the variability of vertical velocity at cloud-base compared to observations at their nominal resolution; the standard deviation of vertical velocity in the model is typically 4-6 times smaller than observed, and even more during the winter at Lindenberg. Averaging the observations to the horizontal scale corresponding to the physical grid spacing of AROME (2.5 km) explains 70-80% of the underestimation by the model. Further averaging of the observations in the horizontal is required to match the model values for the standard deviation in vertical velocity. This indicates an effective horizontal resolution for the AROME model of at least 4 times the physically-defined grid spacing. The results illustrate the need for special treatment of sub-grid scale variability of vertical velocities in kilometer-scale atmospheric models, if processes such as aerosol-cloud interactions are to be included in the future

Explicit representation of subgrid variability in cloud microphysics yields weaker aerosol indirect effect in the ECHAM5-HAM2 climate model

Peer reviewe

Aerosol-landscape-cloud interaction : Signatures of topography effect on cloud droplet formation

Long-term in situ measurements of aerosol–cloud interactions are usually performed in measurement stations residing on hills, mountains, or high towers. In such conditions, the surface topography of the surrounding area can affect the measured cloud droplet distributions by increasing turbulence or causing orographic flows and thus the observations might not be representative for a larger scale. The objective of this work is to analyse, how the local topography affects the observations at Puijo measurement station, which is located in the 75 m high Puijo tower, which itself stands on a 150 m high hill. The analysis of the measurement data shows that the observed cloud droplet number concentration mainly depends on the cloud condensation nuclei (CCN) concentration. However, when the wind direction aligns with the direction of the steepest slope of the hill, a clear topography effect is observed. This finding was further analysed by simulating 3-D flow fields around the station and by performing trajectory ensemble modelling of aerosol- and wind-dependent cloud droplet formation. The results showed that in typical conditions, with geostrophic winds of about 10 m s−1, the hill can cause updrafts of up to 1 m s−1 in the air parcels arriving at the station. This is enough to produce in-cloud supersaturations (SSs) higher than typically found at the cloud base of  ∼  0.2 %), and thus additional cloud droplets may form inside the cloud. In the observations, this is seen in the form of a bimodal cloud droplet size distribution. The effect is strongest with high winds across the steepest slope of the hill and with low liquid water contents, and its relative importance quickly decreases as these conditions are relaxed. We therefore conclude that, after careful screening for wind speed and liquid water content, the observations at Puijo measurement station can be considered representative for clouds in a boreal environment

Modelling mixed-phase clouds with the large-eddy model UCLALES-SALSA

The large-eddy model UCLALES-SALSA, with an exceptionally detailed aerosol description for both aerosol number and chemical composition, has been extended for ice and mixed-phase clouds. Comparison to a previous mixed-phase cloud model intercomparison study confirmed the accuracy of newly implemented ice microphysics. A further simulation with a heterogeneous ice nucleation scheme, in which ice-nucleating particles (INPs) are also a prognostic variable, captured the typical layered structure of Arctic mid-altitude mixed-phase cloud: a liquid layer near cloud top and ice within and below the liquid layer. In addition, the simulation showed a realistic freezing rate of droplets within the vertical cloud structure. The represented detailed sectional ice microphysics with prognostic aerosols is crucially important in reproducing mixed-phase clouds

UCLALES–SALSA v1.0: a large-eddy model with interactive sectional microphysics for aerosol, clouds and precipitation

Challenges in understanding the aerosol–cloud interactions and their impacts on global climate highlight the need for improved knowledge of the underlying physical processes and feedbacks as well as their interactions with cloud and boundary layer dynamics. To pursue this goal, increasingly sophisticated cloud-scale models are needed to complement the limited supply of observations of the interactions between aerosols and clouds. For this purpose, a new large-eddy simulation (LES) model, coupled with an interactive sectional description for aerosols and clouds, is introduced. The new model builds and extends upon the well-characterized UCLA Large-Eddy Simulation Code (UCLALES) and the Sectional Aerosol module for Large-Scale Applications (SALSA), hereafter denoted as UCLALES-SALSA. Novel strategies for the aerosol, cloud and precipitation bin discretisation are presented. These enable tracking the effects of cloud processing and wet scavenging on the aerosol size distribution as accurately as possible, while keeping the computational cost of the model as low as possible. The model is tested with two different simulation set-ups: a marine stratocumulus case in the DYCOMS-II campaign and another case focusing on the formation and evolution of a nocturnal radiation fog. It is shown that, in both cases, the size-resolved interactions between aerosols and clouds have a critical influence on the dynamics of the boundary layer. The results demonstrate the importance of accurately representing the wet scavenging of aerosol in the model. Specifically, in a case with marine stratocumulus, precipitation and the subsequent removal of cloud activating particles lead to thinning of the cloud deck and the formation of a decoupled boundary layer structure. In radiation fog, the growth and sedimentation of droplets strongly affect their radiative properties, which in turn drive new droplet formation. The size-resolved diagnostics provided by the model enable investigations of these issues with high detail. It is also shown that the results remain consistent with UCLALES (without SALSA) in cases where the dominating physical processes remain well represented by both models

Turbulent structure and scaling of the inertial subrange in a stratocumulus-topped boundary layer observed by a Doppler lidar

The turbulent structure of a stratocumulus-topped marine boundary layer over a 2-day period is observed with a Doppler lidar at Mace Head in Ireland. Using profiles of vertical velocity statistics, the bulk of the mixing is identified as cloud driven. This is supported by the pertinent feature of negative vertical velocity skewness in the sub-cloud layer which extends, on occasion, almost to the surface. Both coupled and decoupled turbulence characteristics are observed. The length and timescales related to the cloud-driven mixing are investigated and shown to provide additional information about the structure and the source of the mixing inside the boundary layer. They are also shown to place constraints on the length of the sampling periods used to derive products, such as the turbulent dissipation rate, from lidar measurements. For this, the maximum wavelengths that belong to the inertial subrange are studied through spectral analysis of the vertical velocity. The maximum wavelength of the inertial subrange in the cloud-driven layer scales relatively well with the corresponding layer depth during pronounced decoupled structure identified from the vertical velocity skewness. However, on many occasions, combining the analysis of the inertial subrange and vertical velocity statistics suggests higher decoupling height than expected from the skewness profiles. Our results show that investigation of the length scales related to the inertial subrange significantly complements the analysis of the vertical velocity statistics and enables a more confident interpretation of complex boundary layer structures using measurements from a Doppler lidar.Peer reviewe

Spatial distributions and seasonal cycles of aerosol climate effects in India seen in a global climate-aerosol model

Climate-aerosol interactions in India are studied by employing the global climate-aerosol model ECHAM5-HAM and the GAINS inventory for anthropogenic aerosol emissions. Model validation is done for black carbon surface concentrations in Mukteshwar and for features of the monsoon circulation. Seasonal cycles and spatial distributions of radiative forcing and the temperature and rainfall responses are presented for different model setups. While total aerosol radiative forcing is strongest in the summer, anthropogenic forcing is considerably stronger in winter than in summer. Local seasonal temperature anomalies caused by aerosols are mostly negative with some exceptions, e.g., parts of northern India in March-May. Rainfall increases due to the elevated heat pump (EHP) mechanism and decreases due to solar dimming mechanisms (SDMs) and the relative strengths of these effects during different seasons and for different model setups are studied. Aerosol light absorption does increase rainfall in northern India, but effects due to solar dimming and circulation work to cancel the increase. The total aerosol effect on rainfall is negative for northern India in the months of June-August, but during March-May the effect is positive for most model setups. These differences between responses in different seasons might help converge the ongoing debate on the EHPs and SDMs. Due to the complexity of the problem and known or potential sources for error and bias, the results should be interpreted cautiously as they are completely dependent on how realistic the model is. Aerosol-rainfall correlations and anticorrelations are shown not to be a reliable sole argument for deducing causality

Technical note: Emulation of a large-eddy simulator for stratocumulus clouds in a general circulation model

Here we present for the first time a proof of concept for an emulation-based method that uses a large-eddy simulations (LESs) to present sub-grid cloud processes in a general circulation model (GCM). We focus on two key variables affecting the properties of shallow marine clouds: updraft velocity and precipitation formation. The LES is able to describe these processes with high resolution accounting for the realistic variability in cloud properties. We show that the selected emulation method is able to represent the LES outcome with relatively good accuracy and that the updraft velocity and precipitation emulators can be coupled with the GCM practically without increasing the computational costs. We also show that the emulators influence the climate simulated by the GCM but do not consistently improve or worsen the agreement with observations on cloud-related properties, although especially the updraft velocity at cloud base is better captured. A more quantitative evaluation of the emulator impacts against observations would, however, have required model re-tuning, which is a significant task and thus could not be included in this proof-of-concept study. All in all, the approach introduced here is a promising candidate for representing detailed cloud- and aerosol-related sub-grid processes in GCMs. Further development work together with increasing computing capacity can be expected to improve the accuracy and the applicability of the approach in climate simulations.</p

Implementation of the sectional aerosol module SALSA2.0 into the PALM model system 6.0: model development and first evaluation

Urban pedestrian-level air quality is a result of an interplay between turbulent dispersion conditions, background concentrations, and heterogeneous local emissions of air pollutants and their transformation processes. Still, the complexity of these interactions cannot be resolved by the commonly used air quality models. By embedding the sectional aerosol module SALSA2.0 into the large-eddy simulation model PALM, a novel, high-resolution, urban aerosol modelling framework has been developed. The first model evaluation study on the vertical variation of aerosol number concentration and size distribution in a simple street canyon without vegetation in Cambridge, UK, shows good agreement with measurements, with simulated values mainly within a factor of 2 of observations. Dispersion conditions and local emissions govern the pedestrian-level aerosol number concentrations. Out of different aerosol processes, dry deposition is shown to decrease the total number concentration by over 20&thinsp;%, while condensation and dissolutional increase the total mass by over 10&thinsp;%. Following the model development, the application of PALM can be extended to local- and neighbourhood-scale air pollution and aerosol studies that require a detailed solution of the ambient flow field.</p

A model intercomparison of CCN-limited tenuous clouds in the high Arctic

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 1cm‾³. 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