Representing the number and mass of cloud and aerosol particles independently in a climate, weather prediction or air quality model is important in order to simulate aerosol direct and indirect effects on radiation balance. Here we introduce the first configuration of the UK Met Office Unified Model in which both cloud and aerosol particles have “double-moment” representations with prognostic number and mass. The GLObal Model of Aerosol Processes (GLOMAP) aerosol microphysics scheme, already used in the Hadley Centre Global Environmental Model version 3 (HadGEM3) climate configuration, is coupled to the Cloud AeroSol Interacting Microphysics (CASIM) cloud microphysics scheme. We demonstrate the performance of the new configuration in high-resolution simulations of a case study defined from the CLARIFY aircraft campaign in 2017 near Ascension Island in the tropical southern Atlantic. We improve the physical basis of the activation scheme by representing the effect of
existing cloud droplets on the activation of new aerosol, and we also discuss the effect of unresolved vertical velocities. We show that neglect of these two competing effects in previous studies led to compensating errors but realistic droplet concentrations. While these changes lead only to a modest improvement in model performance, they reinforce our confidence in the ability of the model microphysics code to simulate the aerosol–cloud microphysical interactions it was designed to represent. Capturing these interactions accurately is critical to simulating aerosol effects on climate.
Introduction
Shallow marine clouds are an important source of uncertainties in climate forcing and sensitivity. Representing aerosol effects on these
clouds is a priority for climate modelling efforts worldwide. In this paper, we describe model simulations in a 2∘×2∘ region of the tropical southern
Atlantic Ocean near Ascension Island. The simulations are performed with
the UK Met Office Unified Model (UM), with double-moment aerosol
microphysics driving double-moment bulk cloud microphysics at
500 m horizontal resolution. The cloud and aerosol microphysics parameterisations form, or are intended to form in future, part of the atmosphere model code used for climate simulations in Coupled Model Intercomparison Project (CMIP) experiments and for operational weather forecasts across the Unified Model partnership. The rest of the model code is also used in the HadGEM3-GC3.1 configuration for CMIP6 and in current operational numerical weather prediction configurations. As well as testing the aerosol and chemistry component of the model at a higher resolution than has been attempted before, we study and suggest improvements to the
performance of the aerosol activation scheme. We evaluate the second day of a 2 d long simulation against CLARIFY (CLouds and Aerosol Radiative Impacts and Forcing) aircraft measurements on 19 August 2017.
A series of recent field campaigns have focused on the tropical
south-east Atlantic ocean, which hosts one of the planet's largest
stratocumulus decks, and is the destination for much of the biomass
burning aerosol that originates from central and southern
Africa. The prevailing winds, which are south-easterly in the boundary
layer and easterly in the free troposphere, advect smoke over the
ocean where slow subsidence causes the smoke to mix with the clouds. The
smoke can have large direct, semi-direct and indirect radiative effects on the regional climate . Here we focus on the indirect effect of the aerosols.
In a previous study, evaluated aerosol transport
and microphysics in global and convection-permitting simulations of the south-east Atlantic. They built on earlier work with the same aerosol microphysics scheme employed at high spatial resolution by . In this paper, we start from an updated
version of the same atmospheric model (and increase the resolution
further) to approach the cloud-resolving scale. We also increase the sophistication
of the cloud microphysics scheme from single-moment to a double-moment scheme
in order to study aerosol–cloud interactions in more
detail. The resulting model is the first configuration of the Unified
Model with fully double-moment aerosol and cloud
microphysics. Double-moment microphysics, including prognostic cloud
droplet number concentration, is important to enable good
representations of processes such as aerosol activation and droplet
settling at high spatial and temporal resolutions. The coupling to the double-moment interactive aerosol microphysics scheme enables aerosol-induced
variability in the droplet number concentration. We evaluate the aerosol and cloud microphysics in the new model configuration in this paper, paying particular attention to shortcomings of the simulations that are specific to aerosol–cloud interactions or to simulating aerosols at high spatial resolution. We also highlight some underlying issues with aerosols in the coarse-resolution global climate model that drives our regional simulations, which we will address in future work.
Compared to single-moment cloud microphysics schemes, double-moment schemes have been shown previously to improve the representation of stratiform rain in NWP (numerical weather prediction) simulations and to reduce a range of
biases in high-resolution climate models . The CASIM (Cloud AeroSol Interacting Microphysics) double-moment microphysics scheme we use here is that published previously . In the last of these articles, evaluate CASIM for deep convective
clouds and compare it to a reduced single-moment version of the same scheme and to the different single-moment cloud microphysics scheme used in the
operational version of the model. For the case they study, the
CASIM double-moment microphysics initially performed better
than the single-moment schemes, although these gave comparable
performance to the double-moment scheme when tuned.
In climate simulations with resolutions coarser than around 0.5∘, updraught velocities in shallow clouds are almost entirely unresolved, and convection is parameterised. The activation of aerosols to form cloud droplets requires the supersaturation of water vapour relative to aerosols and hydrometeors to be diagnosed or parameterised. Typically, supersaturation is calculated by imposing an updraught speed (or a distribution of updraught speeds) derived from diagnostics of the sub-grid turbulence rather than the grid-box mean
updraught speed, which is close to zero. A single cloud droplet number concentration per grid
box is thus produced and used in the prediction of rain rates (via an
autoconversion parameterisation) and cloud albedo.
As the resolution of the simulation is increased into the “terra
incognita” or “gray zone” , a higher fraction of
the turbulence in the boundary layer is resolved; this happens until, at the large
eddy simulation (LES) scale of the order of 10 m, we assume for the purposes of this paper that the spatial variability of prognosed updraughts would be a good representation of reality. The turbulence starts to be resolved when the effective grid resolution is below about 4 times the height of the boundary layer , which is typically 4–8 km. There is, therefore, a point at which it is no longer necessary to use an updraught speed diagnosed from sub-grid turbulence in the activation scheme, and the grid-box
mean can be used instead. Prior CASIM simulations have also used grid-box mean updraught speeds to activate aerosols at horizontal resolutions ranging from 250 m to 1 km, and a similar approach has been taken in the Regional Atmospheric Modeling System (RAMS) and some Weather Research and Forecasting (WRF) simulations . The scale invariance of activation
schemes has been tested before down to horizontal resolutions of
1 km. However, these resolutions are still much coarser than typical LES resolutions; therefore, the full variability in updraught speeds will not be resolved. By comparing
near-cloud-resolving and LES simulations,
developed a bootstrapping parameterisation which enables the estimation of the fraction of the variance in updraught that is resolved.
In existing clouds, activation of new droplets will often be negligible but not always, for example if the updraught strengthens towards the top of the cloud. In a detailed model, suggested that in-cloud activation leads to a bimodal cloud droplet size spectrum and is an important factor in accelerating rain formation by broadening the droplet
spectrum . In a case study of deep clouds, the large eddy simulations of suggested that aerosol concentrations in the boundary layer were too low to explain observed droplet number concentrations; therefore, aerosols must be entrained from
the mid and upper troposphere and activated inside the clouds. More recent detailed modelling studies have continued to investigate secondary activation . To fully capture effects of in-cloud activation on the droplet size distribution, a detailed
microphysics scheme is needed, such as a size-bin-resolving scheme or
a super-droplet model. To accurately simulate the full
dynamical response of deep clouds to aerosols, very high spatial
resolution may be needed, which in turn requires relatively short time steps. If the model time step is shorter than the “relaxation timescale” (the timescale for supersaturation production
to be balanced by condensation onto existing droplets), a bias will occur unless supersaturation is represented as a prognostic variable . This timescale is typically a few seconds, depending on the vertical velocity and droplet spectrum. Prognostic supersaturation is desirable in simulations requiring accurate calculations of the latent heat released by condensation , but the short time steps needed mean it is very expensive to treat supersaturation prognostically , although this can be mitigated to some extent with semi-analytic approaches . The relaxation timescale is short (less than 10 s) in the case of thick, polluted clouds with many or large droplets and/or low updraught speeds, and it is long in the case of very clean clouds or strong updraughts . If, however, the model time step is longer than the
relaxation timescale (and we verify that this is the case in our study), one can then assume that supersaturation produced during a time step leads to condensation immediately, so the relative humidity is 100 %
at the end of the time step. This assumption, termed “saturation
adjustment”, may lead to biases in the latent heat released by condensation and in the evaporation of clouds. However, it is much simpler and computationally cheaper than treating supersaturation prognostically, and it works well when updraught velocities are not fully resolved, as in weather prediction and climate models.
In our CASIM microphysics scheme, saturation adjustment is applied. It has sometimes been assumed that prognostic supersaturation is generally part of double-moment microphysics schemes e.g.. However, compared several single- and double-moment bulk microphysics schemes including CASIM with a bin microphysics scheme in a single-column framework. The bin scheme treated supersaturation prognostically, while most of the bulk schemes did not. The double-moment bulk schemes with saturation adjustment that they tested were in closer agreement with the bin microphysics scheme than the single-moment schemes. Their conclusions were substantiated further by . While useful, prognostic supersaturation is not essential for a double-moment microphysics scheme to improve on a single-moment scheme.
In our model, in each time step the activation scheme is rerun assuming there are no existing droplets. If the new droplet concentration is greater than the existing droplet concentration, the old droplet concentration is overwritten by the new one. Changes to the cloud fraction in the grid box may also change the droplet concentration, as discussed later. A similar procedure is followed in the widely used scheme. In this study, we follow the suggestion by , and others to improve on this procedure by accounting for existing cloud droplets when new droplets are activated, assuming a supersaturation that results from a balance between production (updraught) and loss (condensation on existing droplets). This assumption can be seen as a natural extension of the saturation adjustment assumption and should apply in the same conditions. When it was tested in WRF-Chem with cloud microphysics, found that it improved simulated wet scavenging. Our implementation is aware of and consistent with our simulated sub-grid cloud fraction. We emphasise that we seek here to improve the accuracy of
our existing model, but our improvements will not enable it to compete
with detailed studies of aerosol–cloud interactions that employ spectral bin microphysics or prognostic supersaturation.
We attempt to improve the activation scheme further by considering unresolved sub-grid-scale updraught velocities. We examine the suitability of the parameterisation of for our model, but we are not yet able to implement it explicitly and instead derive an ad hoc correction factor suitable for our case study. We show, unsurprisingly, that accounting for existing cloud droplets in the activation scheme reduces the cloud droplet number, while our correction factor for sub-grid-scale updraughts increases it. This may explain why previous studies with CASIM have successfully
produced realistic droplet concentrations with realistic CCN (cloud condensation nuclei), despite
underestimating the updraught and ignoring existing droplets. We examine
the implications of our improvements
for the cloud droplet spatial distribution and size distribution as well as for
rain formation. We evaluate our simulations
against CLARIFY aircraft measurements in 2017 to confirm that they are realistic
and to identify directions for future developments.
Case study and aircraft measurements
The CLARIFY aircraft campaign took place from 16 August to 7 September 2017. The BAE-146 aircraft of the Facility for Airborne Atmospheric Measurements (FAAM) was based at Ascension Island during this period. Extensive sampling of biomass burning aerosol interacting with clouds was achieved during 24 flights usually of around 3.5–4 h duration. The aircraft flew in all directions around Ascension Island, according to where aerosol and cloud transitions could be identified or in order to pass under satellite tracks.
The research aircraft was fitted with a comprehensive suite of
thermodynamic, radiometric, cloud physics and aerosol instrumentation
for the CLARIFY field campaign. Ambient air temperature was measured
using a non-deiced Rosemount/Goodrich type-102 total air temperature
sensor. Atmospheric water vapour was measured using a WVSS-II
near-infrared tunable diode laser absorption spectrometer fed from a
standard flush-mounted inlet . The temperature
and humidity measurements are used to calculate relative humidity in
cloud-free air. The temperature measurement in cloud was subject to
significant wetting effects that led to a cold bias, as illustrated by . The derived in-cloud relative
humidity is therefore set to 100 % in this work. Zonal, meridional
and vertical wind components were derived from the five-port turbulence probe located on the aircraft radome . On the transit at 5180 m altitude from Ascension Island to the cloud, the mean vertical velocity was -0.010ms-1. The small magnitude of this mean velocity compared to the in-cloud vertical velocities we discuss later suggests the probe is sufficiently well calibrated for our analysis, so we do not subtract any baseline offset from the observed updraughts in the evaluation we present in Sect. .
The size distribution of aerosol particles was measured using a passive cavity aerosol spectrometer probe (PCASP, Droplet Measurement Technologies) for nominal diameters between about 0.1 and 3 µm. We applied a complex refractive index of ri=1.54-0.027i, which is appropriate for biomass burning aerosol during the CLARIFY time period , and recomputed the bin boundaries for the PCASP instrument. This resulted in changes to the locations of bin centres, compared to the nominal values from the manufacturer, of usually around 5 % but sometimes up to 20 % in the diameter range below 1 µm. For smaller particles with diameters between 0.03 and 0.3 µm, we also used a scanning mobility particle sizer (SMPS). The total number concentrations of aerosols with diameters larger than about 2.5 nm were measured with a TSI
3786 condensation particle counter (CPC). Aerosol data were only used when in cloud- and precipitation-free air. We determined this using the standard deviation of raw power on the Nevzorov total water content probe , where a power greater than 3.0 mW (∼1.5×10-4gm-3) indicates cloud conditions, following . An additional safety window of 5 s (∼500m) either side of positively identified cloud was applied to
account for diffuse cloud edges and imperfect temporal and spatial
synchronisation between probes and data recording systems. The cloud droplet size distribution (between 2 and 52 µm diameter) was measured with a cloud droplet probe (CDP) that was calibrated with a 10-point bead
calibration . Precipitation-sized particles were measured using a 2D stereo (2DS) probe (10 to 1280 µm diameter) and a cloud imaging probe (CIP-100) (100 µm to 6.4 mm diameter). Data from the CDP, 2DS and CIP-100 were combined to produce a composite PSD (particle size distribution) at a 1 Hz sampling frequency, following the method of . Elsewhere in this paper, the cloud droplet number concentration (CDNC) and liquid water content (LWC) are calculated using the CDP data only. For cloud measurements, in-cloud conditions were determined using a liquid water content threshold of 0.01 gkg-1 calculated
by integrating the cloud droplet number size distribution from the CDP.
The most comprehensive sampling of a cloud feature took place on 19 August. The aircraft flew south of the island to sample a large precipitating cloud structure around 1.5∘ in size. Biomass burning aerosol was detected both within and just above the boundary layer. Cloud-top height peaked at around 2.5 km altitude. The cloud deck was sampled, as shown in Fig. , in a series of five straight-and-level aircraft trajectories along a line of strong radar echoes seen on the aircraft weather radar.
Flight pattern of the CLARIFY flight on 19 August 2017 used in the evaluations presented here (represented by the green line). The left plot is superimposed on imagery from MODIS on the TERRA satellite that corresponds to the morning of the flight, showing the cloud deck that was sampled. The
red box corresponds to the domain of the 500 m resolution simulation. The right panel shows the path of the aircraft as a function of height and latitude.
Taken together, the aircraft observations described later in the paper suggest that the boundary layer is decoupled or cumulus-like, with a stratocumulus cap above large precipitating cumulus clouds beneath. The observed stratocumulus clouds could also be detrained remnants of cumulus. In the simulations, the cumulus sometimes seems to extend up to the top of the boundary layer. However, the size of the cloud feature is larger than average
for stratocumulus-to-cumulus transition clouds, and the deepest cumulus appeared to be organised linearly. Based on this and on geostationary
satellite imagery (not shown), we would describe the cloud pattern as the beginnings
of Flowers rather than Sugar, Fish or Gravel as labelled by.
In order to validate our simulations, we are also able to draw on surface
measurements from the Atmospheric Radiation Measurement (ARM) site at Ascension Island and satellite
observations of cloud droplet concentration and liquid water path. To obtain these, the same procedure as used by was followed.
Model setup
To establish a regional simulation at 500 m spatial resolution, we set up three configurations of the Unified Model (UM): a global model, a 7 km resolution regional model and the 500 m resolution model. The code for all three models is almost identical, except for the cloud microphysics, as described later. The global model is used to produce lateral boundary conditions for the 7 km resolution model, and this in turn provides boundary conditions to the 500 m model. The higher-resolution models do not feed back to the lower-resolution models. (Our setup is usually described as “one-way nesting”.) To illustrate the setup, simulated accumulation-mode aerosol number concentrations are shown in the 7 km model with the global model (top, Fig. a, b) and 500 m model (bottom, Fig. c, d). The configuration of the three models is summarised in Table and described in this section.
Simulated number concentration of accumulation-mode aerosol (abbreviated as “Acc. mode” in the legend labels) in the boundary layer (averaged from the surface to 750 m altitude) on the left and in the lower free troposphere (averaged from 2900 to 3400 m altitude) on the right (in cm-3) at midday UTC on 19 August. The 500 m regional domain is shown as a red square on all four plots; the plots (a, b) show the global and 7 km models with the domain of the plots (c, d) as a red dotted square, while panels (c, d) show the 7 km model and the 500 m model.
Summary of differences between model configurations used. Unless explicitly mentioned, affected by recent minor updates or dictated by the name of the model configuration as described in the appropriate documentation paper, the code for global and regional simulations is identical. In the “Cloud microphysics” row, 1M signifies a single-moment scheme and 2M a double-moment scheme. In the “Aerosol activation” row, ARG stands for the parameterisation of . This parameterisation has two implementations, labelled DIAG when cloud droplet number concentration is diagnostic and PROG when it is prognostic. The microphysics sub-step is marked with an asterisk, because we wish to emphasise that the sub-step is not applied to condensation or evaporation, which are treated by the cloud
parameterisation and not by the microphysics scheme. However, it does apply to precipitation-related processes such as autoconversion and accretion, for example.
Configuration orGlobal7 km model500 m modelparameterisationConfiguration nameGA7.1RA1RA1Number of lat × long grid cells324×432450×670450×450Time step20 min2 min20 sMicrophysics sub-step*2 minnone10 sNumber of vertical levels85 (to 85 km)140 (to 40 km)140 (to 40 km)Code version11.210.811.3Cloud microphysics1M 1M 2M (CASIM; )Aerosol microphysicsGLOMAP GLOMAP GLOMAP Aerosol activationARG (DIAG)ARG (DIAG)ARG (PROG)Sub-grid cloudPC2
Our global model setup is similar to that used by and identical to that used in the intercomparison study of . It follows the GA7.1 configuration , which is the global climate configuration submitted to CMIP6, labelled HadGEM3-GC3.1. The horizontal resolution is
0.8∘×0.55∘ (N216), and there are 70 vertical levels from the surface to 85 km altitude. The horizontal winds are nudged to ERA-Interim starting at 1700 m altitude and ramping up to full strength over the next 453 m of altitude. The relaxation timescale is set to the frequency of the ERA reanalysis files (6 h). The nudging method follows that of . The global model is initialised from the
model for August 2016 used by and run through to 2017.
This global model drives a 7 km resolution regional simulation centred at -16∘, 0∘ (latitude, longitude), with 670 longitude grid boxes, 450 latitude grid boxes and 140 vertical levels from the surface to 40 km altitude. This simulation is initialised on 17 August at 00:00 UTC from the global model and run until the end of 19 August. We use the RA1 configuration of the UM with some settings borrowed from the GA7.1 configuration and UM version 10.8. During the 3 d simulation, air masses in the boundary layer advect from approximately the south-eastern corner to the north-western corner of the model domain. Both the global and 7 km resolution models use the single-moment cloud microphysics scheme of and the double-moment GLOMAP aerosol microphysics scheme . GLOMAP stands for GLObal Model of Aerosol Processes. Because the aerosol microphysics is a double-moment scheme, cloud droplet number concentration can be represented to some extent in the cloud microphysics and radiation schemes but as a diagnostic variable rather than a prognostic. We use “diagnostic” here to indicate that the diagnostic droplet number concentration is calculated from the simulated aerosols, updraught speed and temperature every time step without reference to the droplet concentration in the previous time step, while the prognostic droplet concentration is retained in memory from one time step to the next and advected by the simulated wind fields, though it may also be updated if the simulated aerosol concentration or updraught speed changes. The details are explained in the next section.
The 500 m resolution simulation uses CASIM double-moment cloud microphysics as well as GLOMAP aerosol microphysics. This simulation is driven by the 7 km resolution model and has 450 grid boxes by 450 grid boxes in latitude and longitude, with the same 140 vertical levels as the
7 km resolution simulation. It is centred on -11∘, -14.5∘ (latitude, longitude). Like the 7 km resolution simulation, we use the RA1 configuration, but a more recent UM version, 11.3, is used,
as this has the latest iteration of the CASIM microphysics. Other differences compared to version 10.8 are expected
to have only minor effects. This simulation is initialised from
the global model on 18 August at 00:00 UTC. By 19 August, all of the air
masses that advect into the domain from the boundaries will have been
simulated by the 7 km resolution model rather than the global model for
at least 2 d. Therefore, the resolved wet scavenging processes
evident in Fig. will have affected the aerosol concentrations, and the higher resolution will have had time to affect the winds. In the domain averages we present, we exclude the 20 grid boxes nearest to the domain boundaries to remove the
transition region between the 7 km and 500 m resolution simulations. The
number 20 is arbitrary and chosen by eye; it corresponds to around 30 min of advection time for a wind speed of 5 ms-1: enough time
to produce some more resolved turbulence but not enough time for full
mixing of the boundary layer. We are able to run 500 m resolution simulations driven by the 65 km resolution global model without the intermediate-resolution nest, but then we would need to exclude more grid boxes at the boundaries of the innermost simulation to allow the high-resolution structure to spin up.
The boundary layer scheme we use in all three models is based on that of , which is blended with the Smagorinsky-type scheme from the Met Office Large Eddy Model as grid spacing decreases, as described by . The scheme is expected to be dominated by the scheme rather than the Smagorinksy-type scheme, even at the 500 m spatial resolution, while their contributions would be approximately equal at 250 m resolution.
Following the GA7.1 and RA1 configurations of the UM, the global model uses the PC2 subgrid cloud scheme of , while the regional models use the subgrid cloud scheme of . All three models employ the
area adjustment approach of , so the cloud
fraction seen by the microphysics is the mean of that in three sub-layers of each vertical level in each grid cell, while that seen by the radiation code is the maximum.
In the global model and in the higher-resolution regional models, anthropogenic and natural aerosol emissions are taken from the CMIP5 database , except for biomass burning emissions, which are from the Fire Energetics and Emissions Research (FEER) inventory for August 2017. In addition, sea spray and dust emissions are represented interactively using the parameterisations of and , respectively. The offline-oxidants configuration of the United Kingdom Chemistry and Aerosols (UKCA) chemistry is used together with dust from the CLASSIC aerosol scheme, as in the HadGEM3-GC3.1 climate model .
Simulations at all resolutions are fully coupled to the standard radiative transfer scheme in the UM via the RADAER module for the aerosols . Direct and semi-direct effects of the absorbing aerosols on the cloud are therefore included in the simulations but will
be fairly small in this period due to the relatively low aerosol concentrations and are not the focus of this work.
To illustrate how updraught speeds in clouds are resolved by simulations with different grid sizes, we additionally perform sensitivity simulations at 200 m, 1.5 km and 3 km horizontal resolutions, which are all driven by the 7 km resolution model. These simulations are centred on the same location as the 500 m resolution simulation and have the same configuration, with three exceptions. In the case of the simulation at 200 m resolution, we switch off the sub-grid cloud fraction scheme and assume clouds are fully resolved by the model. In the case of the simulations at 1.5 and 3 km resolutions, we use 224 grid cells in the horizontal directions, instead of the 450 we use for the 200 and 500 m simulations, to save CPU time. Lastly, the time steps in the simulations were adjusted from the 20 s used in the 500 m model to 15 s for the 200 m model, 60 s for the 1500 m model and 120 s for the 3000 m model.
Our setup can be viewed as an update of that documented by . However, as the 500 m resolution model configuration is different to that published by (e.g. most aerosol-related settings are upgraded to GA7.1 from GA6.1), many of the tunings used in our previous simulations are no
longer required. In the wet scavenging code,
changed a parameter designed to represent the fraction of the area of
a grid box over which rain occurs from 30 % to 100 %. We
reverse this change, because while it is still more likely that entire 500 m grid boxes are raining than entire global model grid boxes, we do not account for the evaporation of rain returning scavenged aerosols to the atmosphere. The 30 % parameter can therefore be thought of as the fraction of rain which does not evaporate. We have not verified the accuracy of this assumption, which will clearly depend on the regime studied and should be revisited in future. We no longer tune the dry deposition velocity. The biomass burning emissions diameter (specifically the number geometric mean
diameter) is still 120 nm instead of the default of
150 nm for GLOMAP. This diameter is shown by to give aerosol dry diameters in reasonable agreement (within 40 %) with measurements from the parallel NASA ORACLES campaign (ObseRvations of Aerosols above CLouds and their intEractionS), although the diameter is still slightly overestimated compared to ORACLES. For example, in the lower free troposphere most affected by smoke, the simulated dry diameters are biased 37 % high. The mass of organic carbon in the same location, which dominates the overall aerosol mass, is biased high by 8 %.
Activation, microphysics and coupling
In the global and 7 km models, a diagnostic activation scheme based on the parameterisation of is used to calculate the droplet number concentration . We refer to this as the ARG
(DIAG) scheme later in the paper. Once per time step, the scheme calculates the droplet concentration at cloud base and imposes it on grid boxes that are
above cloud base and still in the same cloud. The droplet concentration does not depend on the concentration in previous time steps. These diagnostic droplet concentrations are used to calculate autoconversion rates in the single-moment cloud microphysics scheme of .
In our 500 m resolution simulations, aerosols activate in the CASIM code to form cloud droplets also using the “ARG” parameterisation of . Other examples of the use of this parameterisation in models at or near the cloud-resolving scale are documented by (their Table 3). The treatment, subsequently referred to as the “ARG (PROG) activation scheme” where PROG refers to the prognostic droplet concentration, is called once per time step in grid boxes when a non-zero mass of water is condensing from the vapour phase. The updraught speed used is set equal to the grid-box mean updraught speed or 0.001 ms-1 (whichever is higher). We reduce this threshold from the original 0.1 ms-1 in this paper, which avoids an unphysical spike in the distribution of droplets (Fig. S1 in the Supplement) and an underestimation of the frequency of low droplet concentrations. Instead of the tuned values used by , the hygroscopicities used are now those recommended by , including a kappa value for organic carbon of 0.2. If the number of droplets activated in a time step exceeds the
number of droplets already existing in that grid cell, the droplet
concentration used in the microphysics and radiation schemes is updated
to the new value. If, on the other hand, the cloud fraction in the
grid box goes down, the droplet concentration is altered in proportion, as discussed later. Cloud droplet number concentrations in our 500 m resolution simulation are also calculated diagnostically in each time step by the ARG (DIAG) scheme , using the same procedure as in the global model described above. These number concentrations are not used by the model's microphysics or radiation schemes, so in future simulations this parameterisation could be switched off. However, for this study we
leave it switched on in order to examine its scale invariance and to compare the predicted droplet concentrations with those from ARG (PROG).
In the sub-grid cloud parameterisation, described by , condensation of water vapour onto cloud droplets is treated with “saturation
adjustment”, so supersaturation is not prognostic and droplets are assumed to be in equilibrium at the end of each model time step. We note that because we use a sub-grid cloud scheme, the grid-box mean equilibrium relative humidity in clouds may be below 100 %. In the CASIM microphysics scheme, autoconversion and accretion are handled by the parameterisation of . In our simulations the clouds are entirely warm phase; for a description of the representation of cold clouds in CASIM, see . In addition to our reduction of the minimum updraught speed used in the activation scheme, we also change the cloud droplet size
distribution assumed by the bulk scheme from an exponential distribution to a gamma distribution. This change is explained further in the context of the evaluation in Sect. .
The coupling from the GLOMAP aerosol microphysics code to the CASIM cloud
microphysics proceeds simply by passing the aerosol mass and number in the
soluble Aitken, accumulation and coarse modes (and the volume-weighted
kappa values of these modes) to the activation scheme in CASIM. The kappa values are parameters which describe the hygroscopicity of an aerosol chemical component . The CASIM microphysical process rates are coupled to GLOMAP aerosols following the procedure used in the default configuration of the UM with single-moment microphysics from . The autoconversion and accretion rates are summed and passed back to the aerosol microphysics code to determine the rate of removal of aerosols in droplets by rain, while the rain and snow rates are used to determine the rate of impaction scavenging of aerosol by precipitation. Autoconversion and accretion also reduce the prognostic droplet number concentration. The liquid water content is
used in the calculation of the rate of conversion of sulfur dioxide to
aerosol-phase sulfate inside cloud droplets in the GLOMAP module of the code. The CASIM microphysics code has the capability to simulate cloud microphysical processing of aerosol ,
e.g. the reduction of aerosol number concentration when cloud
droplets collide and coalesce, or the increase in aerosol number
concentration when rain evaporates. However, there is no capability to
track the composition of the aerosol inside hydrometeors during processing nor to perform aqueous chemistry in the CASIM module. For simplicity and to save on computational expense, therefore, we do not keep track of aerosols in
hydrometeors separately to aerosols outside them, and we keep the aqueous
chemistry in the GLOMAP code. Aerosols that activate to form cloud droplets are only removed if they are wet scavenged (i.e. they form rain), and if this
happens it is irreversible: the evaporation of rain does not return
aerosols to the atmosphere. A coupled GLOMAP-CASIM double-moment model
that includes cloud microphysical and chemical processing of aerosol is deferred to future work.
Developments to the activation scheme
The procedure for aerosol activation adopted in both ARG (DIAG) and ARG (PROG) activation schemes is to activate aerosols once per time step as if no cloud were present whenever there is a tendency for water mass to condense. In ARG (DIAG), the number of droplets that exist in the box before the activation scheme is run is not stored, so the new value of the droplet concentration is used regardless of the previous concentrations. In
ARG (PROG), the number of droplets already in the grid box is
stored for the double-moment CASIM microphysics, and if the new number exceeds the old, the number of droplets is increased to the new value. The overwriting of old droplet concentrations by new droplet concentrations if they exceed the old is the procedure of , and others, but only activated aerosols at cloud base and assumed cloud droplet concentrations were uniform in columns within a cloud, which is a procedure since followed by the ARG (DIAG) activation scheme . The procedure in CASIM, as in models such as CAM5.0 that employ the activation scheme or its successors, is to do activation via ARG (PROG) at all vertical levels within the cloud. In CASIM, the maximum supersaturation is diagnosed assuming no hydrometeors are present in each time step, while above cloud base in CAM5.0 a fixed supersaturation of 0.3 % is assumed . The use of the grid-scale mean updraught in the ARG (PROG) activation scheme, as in RAMS , differs from the scheme of , which was written with climate model resolution in mind and therefore uses a turbulent sub-grid-scale updraught instead of the grid-scale mean, which is mostly unresolved at low resolution.
There are two possible mechanisms by which the ARG (PROG) scheme for
the double-moment CASIM microphysics may overestimate the impact of in-cloud activation, leading to overestimated cloud droplet concentration. First, the effect of existing cloud droplets on supersaturation is neglected by default, but activation is still repeated every time step at all levels in the cloud. Second, provided a cloud does not evaporate, the cloud droplet number produced depends on the maximum updraught speed over the cloud's lifetime rather than the mean updraught speed, though it is not clear if this mechanism leads to an overestimate or is correct.
Existing cloud droplets certainly affect the supersaturation of water vapour in the cloud. In models like ours with saturation adjustment, there is an assumption that the concentrations of water vapour and liquid water reach equilibrium over one model time step, which in our case is 20 s. In pre-existing clouds with sufficiently high liquid water content and droplet number, we may assume the sink of water vapour to unactivated aerosols will be
negligible; therefore, we may write, after and
others, e.g. and , an equation for the time evolution of the supersaturation s using the notation of :
dsdt=α(T)w-γ*GNr‾s.
Here w is updraught velocity, α(T) is a thermodynamic term that relates the updraught to the tendency for water vapour to condense as it cools,
N is the droplet number, r‾ is the number mean droplet radius,
γ* is another term which follows from the thermodynamics of
rising moist air with assumptions detailed in Chap. 12 of , and G is the growth coefficient, which depends on the diffusivity of water vapour in air and on the thermal conductivity of the air. The prescription in Eq. () is valid only for warm-phase clouds; describe a more general approach for mixed-phase clouds. We correct the diffusivity following the
size-independent formulation of except with an accommodation coefficient of 1 as recommended by . If (hypothetically) the system were not in equilibrium and w, T, N and r‾ were constant in time, the supersaturation s(t) could then be approximated by e.g.s(t)=seq+s0-seqexp-γ*GNr‾(t-t0),
where s=s0 at t=t0. In this equation, γ*GNr‾-1 may be interpreted as τ,
the relaxation timescale for the supersaturation, and in liquid clouds with
sufficiently high water content, the supersaturation relevant to aerosol activation is given by the equilibrium or “quasi-steady” value :
seq=αwτ.
Therefore, the concentration of newly activated aerosol in each
aerosol mode i, Nd,i,new, is related to the concentration of aerosol in that mode, Na,i, in the same way as in the activation parameterisations by
Nd,i,new=12Na,i1+Erf2lnseqsc,i32ln(σi),
where σi is the mode width and
sc,i=ra,i-1.54A327Bi
is the critical supersaturation. The coefficient A is a function of temperature, and Bi is a function of the particle hygroscopicity; ra,i is the geometric mean radius of the ith aerosol mode, and the equation assumes that the aerosols are internally mixed . The effect on the supersaturation of the condensation or evaporation of rain is currently neglected in these simulations. While the rain water mass is
non-negligible compared to the cloud liquid water mass, above the cloud base (where it matters) the product of the rain number and the radius is less than 1 % of the product of the cloud number and cloud radius at least in the clouds we study here, so the effect on the relaxation time is negligible. At cloud base, the rain mass concentration can exceed the cloud droplet mass, but the small number of rain drops still means that the effect of rain on the relaxation time is negligible.
In well-established clouds with high liquid water content, , suggested that using Eq. () should be a better approximation for the supersaturation than the maximum supersaturation, smax, generated by activation
parameterisations such as ARG. tested this approximation in large eddy simulations and found the maximum supersaturation, seq, calculated assuming equilibrium with existing droplets to be a much better approximation than that derived using a precursor to the ARG parameterisation , which is valid in the approximately same
conditions as ARG: at cloud base. In WRF-Chem, found that implementing the suggestion improved simulated wet scavenging. We use the quasi-steady-state equation only when it produces a lower supersaturation than the ARG parameterisation. The more detailed microphysics scheme of also uses a diagnostic parameterisation similar to that of ARG at cloud base and a different approach above. However, inside clouds represent supersaturation prognostically (without saturation adjustment) in contrast to our cruder quasi-steady approximation.
In a general circulation or mesoscale model with a sub-grid cloud fraction scheme, partially cloudy grid boxes must be accounted for. By default in
ARG (PROG) within CASIM, the grid-box mean change in cloud droplet number concentration in each time step is
ΔN‾d,t+1Δt=Ft+1×maxÑd,t+1-Ñd,tFt/Ft+1,0Δt,
where Δt is the model time step (20 s in our tests)
Ñd,t+1 denotes the newly calculated in-cloud (denoted by tilde) cloud droplet number concentration, Ñd,t represents that calculated in the previous time step, N‾d,t+1 denotes the newly calculated grid-box mean (denoted by overbar) cloud droplet number concentration and Ft+1 is the current cloud fraction. The Ñd is calculated by running the ARG parameterisation over the whole grid box on the assumption that it is completely cloudy, so the real
grid-box mean cloud droplet number concentration is this value multiplied by the cloud fraction. In some models, Δt is set to an activation timescale rather than the model time step (20 min in the case of ). The multiplication by cloud fraction is done after taking the maximum instead of before taking the maximum to handle the case where the cloud fraction increases but the cloud droplet number decreases. The activation scheme is not run when the cloud fraction decreases. When
a cloud evaporates, there is a homogeneous mixing assumption: the
cloud is assumed to evaporate uniformly across the grid box
so that all of the cloud droplets get smaller as the liquid water
content decreases, and they are only removed when
the cloud fraction or mass of liquid water reaches thresholds close to
zero (10-10kgkg-1 for liquid water content and 10-12
for cloud fraction).
In our improved activation scheme, we assume that in each partially
cloudy grid box the total number of droplets that will either be activated or remain activated in a new time step (denoted by t+1) is the sum of those activated inside the old cloud and those activated in any new cloud that
forms. Therefore, we replace Ñd,t+1 in
Eq. () by
Ñd,cloud,t+1Ft+Ft+1-FtÑd,ARG,t+1Ft+1.
The first term in Eq. () represents aerosols that will be activated at equilibrium in the existing cloud. This term comprises both aerosols that are already contained in large cloud droplets and aerosols that may be newly activated in the existing cloud due to increasing updraught speeds and supersaturations (“secondary activation”). The second term in Eq. () represents aerosols that will be activated in any additional new cloud that forms in the grid box. We thus obtain
ΔN‾d,t+1Δt=Ft+1×maxÑd,cloud,t+1Ft+Ft+1-FtÑd,ARG,t+1Ft+1-Ñd,tFtFt+1,0Δt,
where Ft is the cloud fraction calculated at the end of the previous time step and then advected, Ñd,cloud,t+1 is the in-cloud droplet number calculated using Eq. () for the supersaturation and Ñd,ARG is the result of the ARG parameterisation (or a similar treatment such as that of or ). In our model, Ft is available even when the
sub-grid cloud scheme is diagnostic rather than prognostic, and the equation is not used if Ft+1-Ft<0: if this is the case, we use instead
ΔN‾d,t+1Δt=Ft+1-FtÑd,tΔt.
If the cloud fraction is below 5 % or the in-cloud supersaturation calculated from the cloud droplets via Eq. () is higher than that calculated with the ARG parameterisation, we revert to Eq. (), using only the ARG scheme. The equations for the various conditions are summarised in Table .
Equations for activation to apply in different possible
situations as described in the text.
ConditionEquationApplicationFt+1-Ft>0, Ft+1>0.05 and Scloud<SargEq. ()Activation inside existing cloudsFt+1-Ft>0, Ft+1<0.05 or Scloud>SargEq. ()Activation in newly forming cloudsFt+1-Ft<0Eq. ()Evaporating clouds
With all of the double-moment approaches we are aware of, over the
lifetime of a cloud, the in-cloud droplet number is more likely to increase than it is to decrease, because the droplet number in each grid box is overwritten if the new droplet number exceeds the old. Therefore, the fraction of activated aerosol will end up corresponding to the highest updraught speed seen during that lifetime. divided the
number of additional new droplets activated by two to help compensate for this. This “ratcheting” mechanism is not necessarily unrealistic, since
turbulence-induced upward fluctuations in updraught speed in clouds may well
activate more droplets, while downward fluctuations around a positive
mean updraught are unlikely to lead to large existing droplets completely evaporating. However, the ratcheting means the double-moment scheme
should produce more droplets on average than a single-moment scheme
provided the number concentration of droplets is diagnosed from the aerosol concentration in the single-moment scheme. For example, if the ARG (DIAG) activation scheme were fed by the same updraught speeds and the same aerosol as a double-moment scheme, fewer droplets on average would be produced. The droplet concentrations shown later (in Fig. ) suggest that indeed the ARG (PROG) scheme does produce more droplets, but the comparison is complicated by the ARG (DIAG) procedure of setting vertically constant
droplet concentrations above cloud base.
In the version of CASIM we are using, aerosol processing is not
included, so aerosols are not removed from the gas phase or tracked
separately in the cloud phase. Therefore, in our improved scheme, once a cloud forms and the supersaturation used in the activation scheme is reduced due to the existing droplets, the number of aerosols activated in subsequent time steps within the cloud will be smaller than the number activated the first time the cloud formed. Because the droplet number is only updated if it is higher than the previous droplet number, this will not generally cause the droplet number to decrease, but it will artificially hinder secondary activation from increasing the droplet number. In other words, while our improvements are designed to reduce activation above cloud base, we may have reduced it too much because of the lack of processing. This potential bias could be avoided in future work when aerosol processing is reintroduced to the model.
Aerosol activation is affected by the degree to which updraught speeds in clouds are resolved by a model. If the grid-box mean updraught speed is used in the activation scheme and not all updraughts are resolved, one might expect an underestimate in the concentrations of activated aerosols to result. This bias may counteract the effect of overestimating in-cloud activation discussed earlier. The result of the study of is that the ratio of the variance of simulated updraughts to the variance of real updraughts is given by
R=1-ΔxZmlE1+aΔxZmlE2ΔxZmlE1+bΔxZmlE2+c
for grid resolution Δx; a parameter Zml proportional to the boundary layer height as described in Table S4 in the Supplement; and
fitted constants E1=2.59, E2=1.34, a=7.95, b=8.00, and c=1.05. An additional correction factor f is required to account for the difference between effective and actual grid resolution. This factor varies from model to model, but for the UM it is a factor of 4 correction to the variance . The constant of proportionality between Zml and the boundary layer height depends on the cloud regime. The updraught velocity in the activation scheme can then be set to
wact=wfR
to correct for unresolved fluctuations in the velocity distribution, following Eq. (21) in .
In the cloudy part of the model domain, the mixed layer height is
around 2000 m, so the formulation, accounting for effective grid resolution, suggests a scaling factor of approximately 3 (see Table S4). The
comparison of the PDF of updraught speeds between the model and observations suggests a factor that varies from 1.5 to 4.2, depending on altitude and liquid water content. We give more weight to results at cloud base and to the data binned by liquid water content in Fig. S6, so we use a factor of 2. We do not implement the scheme in full as it appears to exaggerate the scaling factor required for the model resolutions and boundary layer we are studying. Moreover, we note that scaling all our updraughts up by a factor of 2 doubles the mean updraught speed used in the activation
scheme as well as the standard deviation. For the cloud we simulate, at cloud base and at cloud
top this leads
to poorer agreement of the mean updraught speed with observations, while
in the middle of the cloud the agreement improves. The potential
for such a scaling to introduce bias requires further study.
Evaluation of the double-moment microphysics
In this section we evaluate the configuration of the model with double-moment microphysics without the additional developments described in the previous section (to establish a baseline for testing these additional developments). All of our evaluation is performed with a snapshot of the model at 12:00 UTC on 19 August. The aircraft took off at 10:01 UTC, entered the domain of the 500 m model at 11:03 UTC, left it at 13:10 UTC and landed at 13:44 UTC. We assume that changes in meteorological or aerosol variables in the region between 11:03 and 13:10 UTC are small.
Satellite liquid water path and cloud droplet number
Figure shows the liquid water path and cloud
droplet number in the simulation domain compared to the MODIS
instrument on the TERRA satellite . The calculation of droplet number concentration from MODIS effective radius, cloud-top height and optical depth is described by and follows and . The cloud feature is reproduced by the model, but it is substantially smaller than that observed by the
satellite and is displaced to the south-east. The south-easterly
displacement suggests the wind speed in the boundary layer in the 7 km model is likely a little slower than in reality. The liquid water path
seems to match the satellite well in the 7 km simulation. In the 500 m
simulation, the in-cloud liquid water path is also realistic, but the relative frequency of intermediate liquid water paths between 100 and
300 gm-2 is underestimated compared to the satellite. We were able to increase the cloud cover in the 500 m resolution model by tuning the critical minimum relative humidity in the sub-grid cloud scheme, as shown in Fig. S2. However, to maintain consistency with the official RA1 UM configuration , we retain the default values of these parameters in the simulations that follow.
Liquid water path (LWP, a, b, c) and cloud droplet number concentration (CDNC, d, e, f) compared to MODIS on 19 August at midday. Instantaneous output from the 500 m resolution simulation is shown in (a, d), the 7 km resolution simulation in (b, e), and MODIS TERRA in (c, f). The sizes of isolated MODIS pixels are increased during regridding for ease of viewing. There is not always a valid MODIS CDNC retrieval where there is a valid liquid water path, so the cloud extent can only be inferred from (c). The green trace in (d) shows the path of the aircraft.
The retrieval of cloud droplet number does not yield data in
all cloudy pixels due to confusion with overlying cirrus cloud, which
prevents the retrieval of cloud-top temperature. However, from the
valid pixels, the cloud droplet concentration in the 7 km simulation is
realistic. In the 500 m
simulation, a large number of isolated high
peaks in the droplet concentration are seen, while there are no such
high peaks in the
measurements. In the simulation, these peaks correspond to high
updraught speeds at some point in the history of the cloud. It is possible that these peaks correspond to the few grid cells where our saturation adjustment assumption has led to excessive release of latent heat, temporarily strengthening the updraughts. In practice, this variability may not bias the area average: Fig. S3 shows the same figure, except with the central subfigures replaced by the droplet concentration and water path in the 500 m simulation regridded to 5 km resolution, which is not dissimilar to the apparent resolution of the relatively sparsely sampled MODIS cloud droplet number concentration. The simulations in these subfigures show better
agreement with the observations.
Temperature and humidity
Figure shows the aircraft measurements of temperature
and relative humidity compared to the global and 7 km resolution simulations in Fig. a and b. These measurements are supported by Fig. S4, which shows the corresponding data from radiosondes at Ascension Island compared to the regional and global simulations. The simulations clearly underestimate the boundary layer height (by around 400 m), with the 7 km simulation performing slightly worse than the 500 m simulation. The grid spacing at the top of the boundary layer is approximately 80 m. This
underestimate will affect our comparison of simulated cloud properties, especially close to cloud top, to observations, as we discuss later. The underestimate is partly due to an underestimate of the boundary layer height by the ERA-Interim meteorology used to nudge the global simulation, and it is exacerbated by the 7 km resolution simulation. We verified that the
bias persists in the more up-to-date ERA5 reanalysis. We were able to reduce the underestimate by about 100 m by tuning entrainment parameters in the global model, but as this gain is small relative to the overall 400 m bias, we retain the default configuration in the following analysis. We also do not try to correct for this discrepancy in the evaluation by comparing lower altitudes in the model to higher altitudes in observations, as making such a correction would lead to further complications from the different temperatures at the different altitudes. The radiosondes from Ascension Island, which lies just outside the domain of the 500 m resolution simulation, indicate that the boundary layer height
there is lower than the boundary layer height inside the 500 m simulation domain (by approximately 500 m in the observations and 300 m in the 7 km resolution model).
FAAM aircraft observations of temperature (a), relative humidity (b), and out-of-cloud Aitken (AIT) and accumulation-mode (ACC) aerosol number concentration (c) with the corresponding model means over the 500 m domain. The aircraft data are used whenever they are in the domain of the 500 m resolution simulation. Variability across the model domain is assumed here to be small, but for accumulation-mode aerosol, it may be inferred from Fig. . Dotted lines indicate 1 standard deviation. Aerosol data inside clouds are removed as described in Sect. .
As well as the error in the boundary layer height, the temperature is
generally underestimated by the simulations by around 1.5 ∘C
above the boundary layer, although the simulations do reproduce the hint of a secondary inversion at around 3.5 km altitude reasonably well. The relative
humidity is slightly underestimated in most of the boundary layer compared
to the aircraft data. Unlike the humidity derived from the aircraft, the radiosondes record lower humidities than the model. This discrepancy may be due to imperfect matching between the clouds in the model grid boxes and the observations, as we also saw in Fig. , or due to effects from Ascension Island, which is not included in the model. Above cloud, all
simulations produce an elevated relative humidity in the moderately polluted aerosol layer, which is in good agreement with the aircraft measurements in Fig. .
Aerosol number concentration and size distribution
The simulated vertical profiles of aerosol number concentration are compared to observations in Fig. c. The observed accumulation-mode aerosol number concentration in the boundary layer varies from 350 to 700 cm-3 below 1000 m altitude and from 200 to 400 cm-3 between 1000 and 2000 m. The domain mean in the simulation (around 700 cm-3 below 1000 m and 600 cm-3 above) is biased high compared to the observations. The Aitken-mode number concentration above the boundary layer is overpredicted by the model by a larger amount, i.e. more than a factor of 2. The overestimate is due to excessive new particle formation in the upper troposphere in the global simulation. The full vertical profile of Aitken-mode particle concentrations in the global model is shown in Fig. S5. The bias is also present in the evaluation of . The excessive new particle formation is itself due at least in part to an overestimate of sulfur dioxide concentrations in the model (not shown). The overestimation of the Aitken-mode number concentration is likely responsible in part for overestimates in the cloud droplet concentration at the top of the simulated clouds (around 2215 m), as discussed in Sect. . It will be
important to study the scale invariance of nucleation-mode microphysics at high resolution and the biases in its parent model, which also affect the UK contributions to the CMIP6 experiments, in future work.
The observed aerosol size distribution is shown in Fig. . As suggested by the number concentrations in Fig. , the Aitken-mode concentration is consistently overestimated, but the accumulation mode dry diameter is simulated well (generally within 30 % of observations), though it appears to be overestimated by a larger amount at the 1900 m level, which is the altitude of the cloud layer. (Note that cloudy grid boxes are excluded from the average.) The
limited amount of aerosol data available at the level of the clouds
precludes a more quantitative comparison at this altitude. This
overestimation may indicate too much aqueous sulfate production or
too little rainout of the larger particles.
The dry size distributions (out of cloud) observed by the aircraft from a combination of PCASP (black) and SMPS (grey) instruments are shown at different altitudes compared to the model Aitken (orange) and accumulation mode (red) number concentrations, with the total shown as a red dotted line. The simulated coarse mode is not shown. All data shown are out-of-cloud data at standard temperature and pressure. The simulated data are a mean over the domain at 12:00 UTC. Only aircraft data from straight-and-level legs at the specified altitude are shown.
Liquid water content
The aircraft targeted areas of strong radar reflectivity, and so its sampling of liquid water content is biased towards thicker clouds. This bias must be addressed as we compare the observations with the whole range of model grid cells in our domain, as discussed by . The distribution of liquid water content at the five altitudes sampled by the aircraft is compared to the distribution at the same altitudes of the cloudy grid boxes of the 500 m resolution model in Fig. . Unsurprisingly, the observed
liquid water contents are skewed towards high values, while the
simulated liquid water contents are mostly much lower. To compare vertical velocities and cloud droplet number concentrations more fairly between
simulations and observations, for some of our analysis we split the samples at a liquid water content of 0.15 gkg-1. The difference in the distribution of liquid water content between the model and the observations
within each of these bins is then much smaller than the difference over the full range of liquid water content (Fig. ). In some of
the subsequent evaluation, we focus on observations and simulation
data made at liquid water contents above 0.15 gkg-1, as this is where most of the observations lie. The histograms of vertical velocities in two bins of liquid water content in Fig. S7 show that, as expected, high liquid water content is associated with updraughts, while low liquid water
content is more likely in downdraughts.
Frequencies of observed in-cloud liquid water contents from the cloud droplet probe compared to the liquid water contents in cloudy grid cells in the 500 m resolution simulation at the five altitudes marked. We include all model grid cells at the specified altitude within the simulation domain with a liquid water content greater than 0.01 gkg-1 at the instant of 12:00 UTC on 19 August in the histogram, except for the grid cells within 10 km of the domain boundary. The histograms of the simulated liquid water contents are scaled to contain the same number of entries as the histograms of the observations (so, for example, there are 11 model grid boxes plotted at 2550 m altitude).
Updraughts
We compare simulated grid-box mean updraught speeds in our 7 km
and 500 m resolution models (and in the 200 m
sensitivity simulation) to observations in Fig. . We include all
grid boxes where the overall cloud liquid water content is greater than
0.01 gkg-1 irrespective of the cloud fraction. Clearly, the
width of the distribution is underestimated substantially, even in the
200 m simulation. An underestimate is expected, because the 200 and 500 m model grid boxes do not resolve all the updraughts. On the other hand, the aircraft records vertical velocity at 32 Hz, corresponding to approximately one
measurement for every 4 m it traverses; therefore, it
should represent all of the turbulence that would be captured by a
large eddy simulation. In our model, we therefore rely on parameterised boundary
layer mixing. At 200 m resolution, the width of the spatial distribution of updraught
speeds is slightly wider than at 500 m (as shown later in Table ), and the frequency of cloudy grid cells at 2550 m altitude is increased. The increased cloudiness at high altitude is presumably because the in-cloud updraught speeds are better resolved. However, the domain-averaged boundary layer height in these simulations is almost the same as that at 500 m resolution. The low bias in the boundary layer height is mainly due to the driving 7 km resolution and global simulations. In the 7 km simulations, almost no variability in updraught speed in these clouds is resolved.
Updraught speed in clouds in 32 Hz CLARIFY observations and in the model (before changes to the activation scheme) at altitudes of 800, 1450, 1900, 2215 and 2550 m from left to right. Model data are from instantaneous simulation output of the 200 m, 500 m and the 7 km resolution simulations at 12:00 UTC on 19 August 2017. Only the area of the 7 km resolution simulation that overlaps with the 500 m resolution simulation is included. Data are selected as being in-cloud data if the liquid water content exceeds 0.01 gkg-1 in both the model and in the observations. The histograms of simulated vertical velocity are scaled so their maxima match the maximum of the histogram of the observed vertical velocity.
Model resolution explains only part of the underestimated variability
in vertical velocity. The other reason the updraught widths
are too narrow is sampling biases: the sampling of clouds with high radar
reflectivity discussed in the previous section. We compare the
moments of the simulated and observed updraught distribution in bins of cloud liquid water content (in Fig. S6). The distributions and their moments when split into two bins are shown in Fig. S7, which also includes the 200 m resolution simulation, and Table S2. When split by liquid water content, the
width of the vertical velocity distributions in the model matches
observations better, although it is still underestimated. Before
splitting by liquid water content, the mean vertical velocity at each
altitude was generally within 0.1 ms-1 of zero, but after
the split this is no longer the case, as indicated in Table S2. We note that the two variables are not independent: large positive liquid water
contents are correlated to large positive vertical velocities, because high
vertical velocities imply high supersaturations. The reader is also reminded that the 2215 m altitude level is probably more representative of cloud top in the model than the 2550 m level, due to the underestimated height of the boundary layer.
If we approximately correct for the sampling bias associated with the
path of the aircraft by only considering vertical velocities where liquid water content exceeds 0.15 gkg-1 and also smooth the observed vertical velocity distribution (to show the mean observed vertical velocity
over approximately 500 m), we would expect good agreement of
the model with observations. When we do this, as shown in Fig. S8, we do see
substantially improved agreement of the model with observations at cloud base and cloud top but not in the middle of the cloud: the model still
underestimates the variability in updraught speed. The residual bias may
be a sampling artefact we did not successfully remove, or it may
indicate that there is not enough simulated convection.
Droplet concentrations
The simulated in-cloud droplet number concentration is calculated by
first dividing the grid-box mean droplet number concentration by the
cloud fraction and then removing any grid cells where the in-cloud mean
liquid water content is below 0.01 gkg-1 or where the cloud fraction is below 0.05. With the ARG (PROG) activation scheme, Table shows that the mean droplet concentration is overestimated by a factor of 2 to 3, depending on altitude in the clouds with higher liquid water content,
and by a larger factor in clouds with low liquid water content. The
overestimate is likely to be due in part to the substantial
overestimate in aerosol concentrations in the accumulation mode and
in part due to the activation scheme and warm rain representation (discussed later). Figure shows the spatial distributions of the
simulated and observed droplet concentrations (using the whole sample and not split by liquid water content). We do not plot droplet concentrations at 2550 m altitude due to the low sample size. The figure shows the variability (the width of the distribution) is also overestimated, although the frequency of very low droplet concentrations is underestimated. Table shows that the
overestimate in mean droplet concentration is more severe than
the figure suggests, because the aircraft targeted areas of strong
radar reflectivity and therefore preferentially sampled high liquid
water contents and high droplet concentrations. At 2215 m altitude, it seems likely that the overestimated mean and width of the droplet concentration spatial distribution is at least partly due to activation of the Aitken mode, given that the number concentration of the Aitken mode is substantially overestimated as discussed in Sect. .
Cloud droplet number concentrations in all clouds and in two bins of liquid water content. The 500 m resolution, 7 km resolution and global models are shown, with the ARG (DIAG) and ARG (PROG) activation schemes where they are run. In both model and observations, the threshold liquid water content to define a cloud is 0.01 gkg-1. All data are shown above the first horizontal line; the model grid cells and observations with in-cloud liquid water content above 0.15 gkg-1 are shown above the second horizontal line, and those with liquid water content below 0.15 gkg-1 are shown below. The global model, labelled “Global DIAG”, is not separated by liquid water content as the number of grid boxes in the domain is too small.
AltitudeObs.500 m PROG500 m DIAG7 km DIAGGlobal DIAG(cm-3) 800200253144118180145020622988921801900144227799214522151383011922501802550101370364––800249512280180–1450247416164139–1900173421131127–2215158383234250–2550126388338––800117220127116–1450891836986–1900451746578–221545156118––255014182647––
Cloud droplet number concentration in CLARIFY observations at 1 Hz; in the 500 m resolution model as predicted by ARG (PROG), labelled “PROG”; in CASIM and by ARG (DIAG) with its updraught PDF, labelled “DIAG”; and in the 7 km resolution model as predicted by ARG (DIAG), labelled “DIAG 7” at altitudes of 800, 1450, 1900, and 2215 m from left to right. The histograms of simulated cloud droplet number concentration are scaled to contain the same number of entries as the 1 Hz observations. The histograms show droplet numbers over the full range of in-cloud liquid water contents. The observations are compared to grid boxes sampled from the whole model domain, except for the 20 grid boxes nearest the boundaries. Values are grid-box means rather than in-cloud means. Strictly, the aircraft traverses 500 m in about 4 s, so 1 Hz is too high a sampling frequency for comparison to the model. However, down-sampling to 0.25 Hz reduces the data sample without changing the shape of the distribution.
The ARG (DIAG) activation scheme in the same simulation also overestimates the variability in droplet concentration, but the mean is closer to observations. Most likely, ARG (DIAG) produces fewer droplets than ARG (PROG),
because the ARG (PROG) droplet number depends on the highest updraught in the
history of the cloud, while the ARG (DIAG) droplet number depends on the
updraught at the particular instant the droplet number is diagnosed, as
discussed in Sect. . In ARG (DIAG) at 7 km resolution, by contrast, the variability is underestimated,
and the mean droplet concentration is also underestimated by a factor
often around 2. This ARG (DIAG) activation scheme (so far only used at
convection-permitting resolution by )
also produces a spike in the very lowest bin. This spike is the result of a
lack of scale awareness in the scheme as it is currently coded. It corresponds to the minimum droplet concentration being assigned 5 cm-3, because the characteristic updraught speed is out of the range allowed in the code. (A more detailed discussion is given in the Supplement.) Grid boxes with out-of-range updraughts
are seen in the 7 km resolution simulations as well as the
500 m simulations. We describe this effect in more detail (and propose a fix to this specific problem) in the Supplement.
The lack of scale invariance in the ARG (DIAG) diagnostic activation scheme also means the 7 km simulation does not yield the same distribution of cloud droplet concentrations as the 500 m simulation. At 1900 m
altitude the simulations agree on the domain-mean droplet
concentration: it is 86 cm-3 in the 500 m simulation and 92 cm-3 in
the 7 km simulation. However, at 800 m altitude they do
not agree: the droplet concentration is 181 cm-3 in the 500 m
simulation and 117 cm-3 in the 7 km
simulation. The reduced variability in the 7 km simulation compared to the 500 m
simulation is expected, but further work is needed to ensure the means
are consistent. However, because the diagnostic droplet concentrations in the
7 km model do not feed through the lateral boundaries of the
500 m simulations, biased droplet concentrations in our
7 km model should not substantially affect the
500 m model.
The spike at 5 cm-3 is not observed in the
global climate model, and partly because of this, there is no substantial low bias in the mean droplet
concentration. Other reasons for the lack of low bias in the climate
model could be the lack of resolved wet scavenging of aerosols. A
global evaluation of droplet number in the climate model is presented
by .
We considered the possible effect of biases in the underlying ARG
algorithm on the results of the ARG (PROG) activation scheme by
comparing it to a cloud parcel model with an explicit activation scheme following Köhler theory . Figure S9 shows the fraction of accumulation-mode
aerosols activated at cloud base as a function of updraught speed from
simulations with our model, which are run at 200 m resolution with no sub-grid
cloud fraction (for ease of interpretation). We ignored possible
contributions from the Aitken and coarse modes. Superposed on the
figure, we show simulations with the parcel model in red and the
predictions of the ARG algorithm run offline in orange. We used the
Pyrcel parcel model of to perform the
parcel model and offline ARG calculations. We set the accumulation-mode
aerosol and thermodynamic parameters to be consistent with those included
in, or simulated by, the Unified Model. We assumed the aerosol was composed of ammonium sulfate (kappa =0.61) and additionally ran the parcel model again using a kappa value of 0.2 instead of 0.61 to see how the results would change if the aerosol was instead organic carbon. The activated fractions from the parcel model and the ARG parameterisation agree to within 20 % for both hygroscopicities, confirming that the ARG parameterisation is appropriate. The fractions in the UM are more scattered, because cloud droplet number is prognostic and can undergo autoconversion, accretion and sedimentation.
Cloud and rain particle size distributions
We compare the cloud and rain size distributions between the
model and the observations in Fig. . We find the model represents large cloud droplets well but overestimates the number of small droplets and underestimates the number of large rain drops. The cloud droplet size distributions are broadly similar for cloud liquid water contents
below and above 0.15 gkg-1, as shown as Fig. S11. The simulated surface rainfall amount is compared to observations from the Global Precipitation Measurement (GPM) satellite
dataset in Fig. S12. Despite the poor representation of the rain size distribution observed by the aircraft, the total surface rainfall rate simulated by the model is in generally good agreement with the satellite
data. Simulated vertical profiles of rain mass and number concentration are
evaluated in Fig. , which is shown later as it includes model
developments discussed in Sect. . This evaluation confirms that, compared to the aircraft data, the rain mass is generally well
simulated, but the number concentration of small drops is overestimated
and that of very large drops underestimated. We note that the
autoconversion and accretion parameterisations from (or the parameters of the gamma size distribution) could in principle be
tuned to attempt to address this bias, but we did not attempt to do
so. Such tuning would likely only be applicable to the clouds that were sampled by the aircraft, which, as we have discussed, are not necessarily even a fair sample of the clouds in our model domain as the aircraft was targeting areas of high radar reflectivity .
Evaluation of cloud and rain droplet size distributions (in micrometres) at the altitudes sampled by the aircraft. The simulations are from the 500 m resolution model. The observed size distributions are shown by solid lines and the simulated distributions by dotted lines. The observation data are combined from the cloud droplet probe (CDP), two-dimensional stereo probe (2DS) and cloud imaging probe (CIP). The CDP data and simulated cloud droplet concentrations are filtered to be in cloud, i.e. to include only samples with LWC>0.01gkg-1, while the CIP data and simulated rain are not.
In the simulations shown in Fig. , we have used a gamma size
distribution instead of the default exponential cloud droplet size distribution assumed by CASIM and published by , and others. We find (not
shown) that once both the moments for droplet mass and number concentration are predicted by the model, adjusting the prescribed size
distribution shape (e.g. changing from an exponential to a gamma distribution) without changing these moments does not strongly affect the important results of the microphysics
parameterisation, such as rain formation rates, in the case we study, though this may not be true in other clouds. However, in Sect. we calculate the number mean droplet radius in order to determine supersaturation, and this is sensitive to the shape. For example, an exponential size distribution has a number mean
radius equal to one-third of the effective radius and 0.55 times the
volume mean radius, while the more realistic gamma size distribution
with μ=5 has a number mean radius equal to 75 % of the effective
radius, and the effective radius is larger than the volume mean radius
by 15 %. This 15 % is still larger than the 8 % found in observations but in much better agreement than the 66 % factor by which the effective radius of the exponential distribution is greater than its volume mean radius. The gamma size distribution is given by
N(D)=Ndλμ+1Dμe-λDΓ(μ+1)
for droplet diameter D, shape parameter λ and droplet concentration Nd. The size distributions shown in Fig. and used in our subsequent simulations have μ=5. The use of μ=5 is inspired by the relationship of
used by , which yields μ=5 for a
cloud droplet number concentration of 240 cm-3, which is close to the
concentrations observed here. The exponential size distribution, the size distribution used by and the gamma distribution with μ=5
are compared in Fig. S10, and the number mean droplet radii that result
are tabulated in Table S3. The gamma size distribution we use
overestimates the number of small droplets compared to observations,
but it captures the number of larger cloud droplets very well.
Simulations by and others showed that in principle a bimodal droplet spectrum could arise as a result of secondary droplet nucleation. We have some indications, from Fig. , that updraughts may accelerate in our clouds, leading to secondary activation, as strong updraughts, above 4 ms-1, are more frequently observed at higher altitudes. It is then tempting to claim that the droplet spectrum
in Fig. is a bimodal distribution with peaks at 7 and about 20 µm. When we filter the observed size distribution to exclude low liquid water contents in which secondary activation is not expected, we find (on the right of Fig. S11) that at high altitudes in the cloud the slightly bimodal distribution remains. However, it could easily be an instrument artefact, and it could probably also arise from different aerosol types, entrainment or collision–coalescence processes. Thus it is very speculative to suggest secondary activation is actually taking place. In our
simulations, by contrast, there is no mechanism for secondary activation to
lead to a bimodal size spectrum: the new droplets must fit into the existing size spectrum definition.
Results of improvements to the ARG (PROG) activation scheme
For our 500 m simulations, Fig. and Table show that when we correct the ARG (PROG) activation scheme in our model to account for the effect on supersaturation of existing droplets and the effect of unresolved vertical velocities, we obtain relatively similar cloud droplet numbers to those we started with, with a slight improvement to the agreement of the model with observations. Overall, it looks like the factor 2 correction to the updraught (as suggested by observations at cloud base) produces cloud droplet number in good agreement with the default model for this resolution. However, as well as being resolution dependent, this finding may be specific to the type of cloud we study and will likely be influenced by other
sources of bias – for example, the assumption that the aerosols are internally mixed, the assumptions about hygroscopicity or errors in
simulating the aerosol size distribution.
Liquid water content, maximum supersaturation “Smax” calculated via the ARG parameterisation and the in-cloud equilibrium “EQM”, supersaturation relaxation time “tau”, instantaneous concentration of new droplets activated “CDNC (inst)” and prognostic droplet number concentration “CDNC (prog)” at 12:00 UTC on 19 August. Each plot shows the old version of the 500 m resolution model, labelled “Default”; the case where the reduction of supersaturation in clouds is accounted for, “in-cloud SS”, and the cases where this reduction is accounted for and the updraught speed in the activation scheme is increased by a factor 2, “in-cl SS, w×2”. In-cloud medians are shown (so the liquid water content is above 0.01 gkg-1, and we divide by the cloud fraction), and except in the case of liquid water content and prognostic cloud droplet number (a, f), the plots also only show means over grid cells where a positive water mass is condensing. The observed median liquid water content and droplet concentrations are shown as black dots in panels (a, f). Where shown, error bars or shading indicates the interquartile range of the sampled observations of the model grid cells. The plot of the relaxation time, tau, only includes grid cells where the relaxation time is used – i.e. where the cloud fraction exceeds 0.05 and the equilibrium supersaturation in the cloudy part of the grid cell is lower than the ARG supersaturation. We also show shaded in grey the vertical band containing 95 % of the total liquid water content, because otherwise the eye is drawn to the higher altitudes. Above the shaded region, there are only a few cloudy grid cells, which do not have a strong impact on the domain-mean
radiative properties, because the tops of most of the clouds in the
domain are at or slightly below the top of the shaded region and not above it.
The mean maximum supersaturation from the ARG parameterisation is around
0.25 % at cloud base in the default simulation and around 0.35 % when the updraught speed is multiplied by two, which is in line with expectations . The maximum diagnosed supersaturations in new clouds (calculated by the ARG parameterisation) and in existing clouds calculated with the quasi-steady-state assumption are shown as in-cloud domain-mean vertical profiles in Fig. and as in-cloud histograms in Fig. S13. They demonstrate that in-cloud activation is clearly not negligible even in this relatively shallow cloud: the mean in-cloud maximum supersaturation can exceed 0.1 %. However, in the cloud we study, it is lower than the 0.3 % prescribed for convective clouds above cloud base by . The standard deviation of the maximum supersaturations (not shown for clarity) is also around 0.1 % for the case when the updraught speed is multiplied by two and 0.05 % when the updraught is not corrected. The mean in-cloud value is similar to the 0.1 % found in a smaller cumulus cloud by . Very occasionally, the relaxation time can exceed the model time step of 20 s, as shown in the histograms in Fig. S13 and in this case a bias will result, but the number of grid boxes in which this happens is not significant.
When we account for the effect of existing cloud droplets on supersaturation in the activation scheme, the instantaneous concentration of new droplets activated is always lower above cloud base than when we do not. Increasing the updraught speed increases the activation at cloud base, leading to more activation at cloud base, but for a factor of 2 increase in updraught, the number of droplets activated inside the cloud remains lower than in the original model. This difference in vertical profile propagates to the
prognostic concentration of cloud droplets, and the consequence is a lower cloud droplet concentration where it is most important for radiative
transfer – where the domain-mean liquid water content is highest close
to cloud top – compared to the droplet concentration at cloud base.
The spatial distributions of cloud droplet number concentration, shown in Fig. , match the aircraft slightly better after our modifications to the activation scheme, as the PDF of the cloud droplet number concentration assembled from the model grid cells narrows slightly. However, the width of the PDF is still biased to be too wide and sometimes too concentrated at low droplet number concentrations.
Comparison of in-cloud cloud droplet number distributions in the five straight-and-level runs at altitudes of 800, 1450, 1900, 2215 and 2550 m from left to right. Each plot shows the old version of the model, labelled “Default”; the case where the reduction of supersaturation in clouds is accounted for, “in-cloud SS”; and the cases where this reduction is accounted for and the updraught speed in the activation scheme is increased by a factor 2, “in-cl SS, w×2”. Only model grid boxes and CDP observations where the liquid water content exceeds 0.15gkg-1 are shown.
Table shows a comparison of the simulated and observed width of the spatial distribution of the in-cloud updraught speed before scaling. For the 500 m resolution simulation, the factor of 2 correction is clearly appropriate, except at 1900 m altitude where it is too small. By contrast the correction of a factor of 3 suggested by the parameterisation of would be too large. We also ran simulations with 200 m, 1.5 km and 3 km resolutions. For the 3 km simulation, the scheme would suggest that a scaling factor of about 8 is needed, and we find similarly that for most altitudes this would
increase the width of the updraught PDF substantially beyond that which
is observed (Table and Fig. S13). Similarly, at 200 m resolution the factor is between 2 and 3, which is also too high. Further work is needed to understand the reasons for this overestimated correction factor fully, but the underlying premise that the updraught may be scaled up to account for its unresolved fraction may still be appropriate.
Standard deviation of in-cloud updraught speed in observations and in
simulations at four different resolutions on 19 August 2017 at
12:00 UTC.
Our improved activation scheme gives a modest improvement to model performance in the cloud we chose, but we have not demonstrated whether it would work in any cloud. In a thin cloud with high cloud fraction but low liquid water content, the relaxation time might reach values comparable to the model time step, which will most likely lead to biases. To avoid unphysical results, the quasi-steady-state scheme could be switched off if the relaxation time is equal to or exceeds the time step, and the ARG scheme could be used by itself instead. In our simulation, the relaxation time is either well below the time step or the quasi-steady-state supersaturation is higher than the ARG
supersaturation, and in this case the quasi-steady-state supersaturation is not used. However, there is no guarantee that these conditions would always be satisfied with different time step lengths or in different cloud types.
The effect on cloud microphysics of the changes to the activation scheme is
detailed in Fig. . Reducing the cloud droplet number
concentration by accounting for the effect of existing clouds on
supersaturation increases rain number and mass concentration, unsurprisingly; increasing the updraught speed in the activation scheme reduces rain formation in the very highest clouds. However, the impact of the results on microphysics is small: there is, in particular, no substantial impact on cloud liquid water content in these simulations, and there is no effect on the altitude at which rain forms.
Mean in-cloud cloud droplet number concentration in the three simulations in Figs. and (original, supersaturation corrected in clouds, and supersaturation corrected and updraught speed increased) and in the observations from the FAAM aircraft cloud droplet probe. Model grid boxes and CDP observations where the liquid water content exceeds 0.15 gkg-1 are shown above the horizontal line and those with lower liquid water contents below.
Rain and the rain formation process for observations (black dots) and model grid cells (coloured lines) where the liquid water content exceeds 0.15 gkg-1 in the three simulations shown in Figs. and (original, supersaturation corrected in clouds, and supersaturation corrected and updraught speed increased). In CASIM as in many microphysics schemes, rain formation begins with an autoconversion rate, labelled “Autoconv”. Rain droplets then accumulate liquid water by accretion. All plots except (f) show in-cloud medians where we only include grid boxes with cloud fraction greater than 0.05, and the grid-box value is divided by the cloud fraction. Where shown, error bars and blue shading represent the interquartile range, while the grey shading shows the region where 95 % of the domain-mean liquid water content resides. The last plot shows the domain-mean rain mass mixing ratio with no liquid water content threshold applied.
Discussion and conclusions
We have coupled the GLOMAP and CASIM two-moment aerosol and two-moment cloud microphysics components within the Unified Model and tested the resulting model at 500 m horizontal resolution against CLARIFY campaign data. The new configuration is intended to lead to improved simulations of aerosol–cloud interactions. We made some additional developments to the model, which is summarised in Table .
In the case study we simulated, smoke and marine aerosols are emitted and propagated to the neighbourhood of Ascension
Island, where they interact with clouds and are scavenged by
rain. The CLARIFY flight on 19 August 2017 has proved a useful
test bed for model evaluation. The five level in-cloud aircraft flight
segments allowed large data samples to be obtained with
state-of-the-art instrumentation. Additional profiles and saw-tooth
segments are also available for analysis. The cloud deck, which
is just under 200 km across, is substantial but self-contained, allowing plenty of data to be obtained
both in and out of cloud. The case study can be simulated
reasonably well with our numerical weather prediction model. While other cloud cases have
been measured in more detail, this is one of the best cases from the
south-east Atlantic stratocumulus-to-cumulus transition zone and would be
suitable for future detailed studies of aerosol–cloud interactions.
Our evaluation of aerosol and cloud microphysics is ambitious as it relies on our hierarchy of simulations accurately representing the complete cycle of aerosol emissions, transport and
deposition at the synoptic scale, as well as cloud microphysical
processes. We highlighted some
areas where the models perform well but also some
shortcomings. The 500 m resolution model simulates the boundary layer clouds
realistically, although the top of the boundary layer is low and the
cloud we focus on is not quite simulated in the right
place. Unfortunately, the low boundary layer top means the matching of
the simulated cloud to observations at the same altitudes leads to an
imperfect comparison of other cloud properties such as droplet
concentration and updraught speed. However, accumulation-mode aerosol, cloud and rain number concentrations are realistic: simulated mean accumulation-mode aerosol and cloud droplet number concentrations are generally within a factor of 2 of observations, although larger discrepancies exist in certain cloud regimes; for example, droplet number concentration is more severely overestimated in clouds with low liquid water contents. Rain mass is also simulated realistically (within a factor of 2), but rain number concentration is overestimated by around a factor of 3. Further work on the global model aerosol code is needed to address the substantial overprediction of Aitken-mode aerosol number concentrations. Cloud liquid water content and rain mass compare reasonably well to satellite
measurements. We found the ARG (DIAG) aerosol activation
requires a fix to work once updraught speeds are partially or fully resolved.
We examined two possible improvements to the ARG (PROG) aerosol
activation scheme in CASIM: the correction to the supersaturation for existing cloud and the correction to the updraught speeds. A gamma distribution with μ=5 was found to be needed for the first of these improvements to work instead of the exponential distribution usually used for the CASIM cloud droplet size distribution, and so all of our results are presented with the gamma distribution. Further work is needed to ensure the correction scales with resolution correctly before it can be implemented online. Conversely, the ARG (PROG) activation scheme will break at coarser model resolution, as it is not possible to apply the correction unless some fraction of the updraughts are resolved. To get
around this, a sub-grid vertical velocity of the form suggested
by (or a PDF as in ARG (DIAG)) could be implemented into
CASIM and switched on for low grid resolutions (probably coarser than
around 3 km, depending on the cloud type) so that aerosols
will still activate when the grid-box mean updraught speed is zero. We also recommend reducing
the minimum updraught speed in the ARG (PROG) activation scheme from 0.1 to 0.001 ms-1 to avoid unphysical spikes in the distribution of cloud droplet number.
Summary of model developments and adjustments documented in this article.
SectionDevelopmentSection GLOMAP-CASIM coupling via activationSection GLOMAP-CASIM coupling via scavenging/chemistrySection Adjustment of μ in cloud droplet size distributionSection Reduction of CASIM threshold updraught for ARG (PROG) activationSection In-cloud droplet activation in ARG (PROG) activation schemeSection Sub-grid updraught correction in ARG (PROG) activation scheme only implemented as a case-specific tuning factor
We have now made the model more physically sound by reducing
incorrect in-cloud activation and using vertical velocities to
activate aerosol that is closer to the real vertical velocities
observed and without introducing any new tuning parameters or
computational expense. For high-resolution simulations, we believe the approach to secondary activation we adopted following , , and is relatively generalisable and could be readily applied to other models that represent sub-grid cloud fraction, e.g. WRF-Chem with
or microphysics and other activation
schemes such as that of . We note that the microphysics uses an updraught speed that is already adjusted
for sub-grid turbulence and so has no need of the
correction. Implementing only the correction to the supersaturation
(and no correction to the updraughts) would lead to lower simulated droplet concentrations.
There remain arguments for either switching off droplet nucleation
above cloud base altogether, assuming the cloud droplet concentration is uniform in vertical columns above cloud base, or assuming the supersaturation is constant above cloud base as done in convective clouds in CAM5 by . In-cloud activation leads to a broadening of the cloud droplet size distribution above cloud base, but with our approach, any new droplets activated above cloud base must follow the same size distribution as existing droplets . However, we have still improved on the original method where the importance of activation above cloud base was substantially exaggerated.
The improvements to the CASIM microphysics and its ARG (PROG)
activation scheme may lead to more reliable simulations of aerosol
activation and aerosol–cloud interactions that depend more on the
concentration of fresh aerosols at cloud base and less on the more
processed aerosols found higher in the cloud than previous
simulations using similar models, while still allowing activation
above cloud base. However, we emphasise that we have not tested
the improved model in deep or cold clouds, and we do not
suggest it could replace more sophisticated activation parameterisations, as some of which already calculate in-cloud supersaturation prognostically e.g.. In particular, the explicit prognostic calculation is likely to be needed for calculations of convective invigoration since it is likely that latent heat release and cloud evaporation rates will be erroneous in bulk models with saturation adjustment e.g.. However, more approximate calculations, used in many if not most convection-permitting studies, may benefit from the improvements to the procedure for aerosol activation we tested, with no additional computational cost.
Data availability
Data from the FAAM aircraft are available on the CEDA (Centre for Environmental Data Analysis) repository http://archive.ceda.ac.uk/ (last access: August 2019, ). Simulation data are stored at the Met Office and available from the corresponding author on reasonable request. The Met Office Unified Model is available for use under licence. A number of research organisations and national meteorological services use the UM in collaboration with the Met Office to undertake basic atmospheric process research, produce forecasts, develop the UM code, and build and evaluate Earth system models. For further information on how to apply for a licence, see
https://www.metoffice.gov.uk/research/approach/modelling-systems/unified-model/index (last access: 1 June 2020, ).
The supplement related to this article is available online at: https://doi.org/10.5194/acp-20-10997-2020-supplement.
Author contributions
HG, KSC, PRF and AAH developed the concepts and ideas for the
direction of the paper. HG, PRF, DPG, AAH and JW coupled the GLOMAP and
CASIM aerosol and cloud microphysics codes. SJA, PB, KB, IC, ZC, JT,
and HW collected and calibrated the measurement data. HG carried out and
analysed the model runs, performed the model evaluation with input
from PRF, SJA and PB, and wrote the paper with input and comments
from KSC, PRF, SJA, PB, ZC, DPG, and AAH.
Competing interests
The authors declare that they have no conflict of interest.
Special issue statement
This article is part of the special issue “New observations and related modelling studies of the aerosol–cloud–climate system in the Southeast Atlantic and southern Africa regions (ACP/AMT inter-journal SI)”. It is not associated with a conference.
Acknowledgements
We thank the Facility for Airborne Atmospheric
Measurements, who operate the aircraft and some of the instrumentation
used for the CLARIFY campaign, and everyone involved with the campaign.
Radiosonde data were obtained from the Atmospheric Radiation Measurement (ARM) Climate Research Facility, a U.S. Department of Energy Office of Science user facility sponsored by the Office of Biological and Environmental Research. The Terra/MODIS cloud optical depth and effective radius datasets were
acquired from the Level-1 and Atmosphere Archive & Distribution System (LAADS) Distributed Active Archive Center (DAAC), located in the Goddard Space Flight Center in Greenbelt, Maryland, USA (https://ladsweb.nascom.nasa.gov/, last access: 4 July 2018). We acknowledge use of the Monsoon2 system, a collaborative facility supplied under the Joint Weather and Climate Research Programme, a strategic partnership between the UK Met Office and NERC. We thank our two anonymous reviewers for their helpful suggestions, which improved the paper.
Financial support
This research has been supported by the Natural Environment Research Council (grant no. NE/L013479/1).
Review statement
This paper was edited by Joshua Schwarz and reviewed by two anonymous referees.
ReferencesAbdul-Razzak, H. and Ghan, S. J.: A parameterization of aerosol activation: 2. Multiple aerosol types, J. Geophys. Res., 105, 6837–6844, 10.1029/1999JD901161, 2000.Abel, S. J. and Boutle, I. A.: An improved representation of the raindrop size distribution for single-moment microphysics schemes, Q. J. Roy. Meteor. Soc., 138, 2151–2162, 10.1002/qj.1949, 2012.Árnason, G. and Brown, P. S.: Growth of cloud droplets by condensation: A problem in computational stability, J. Atmos. Sci., 28, 72–77, 10.1175/1520-0469(1971)028<0072:GOCDBC>2.0.CO;2, 1971.Barrett, P. A., Blyth, A., Brown, P. R. A., and Abel, S. J.: The structure of turbulence and mixed-phase cloud microphysics in a highly supercooled altocumulus cloud, Atmos. Chem. Phys., 20, 1921–1939, 10.5194/acp-20-1921-2020, 2020.Bellouin, N., Mann, G. W., Woodhouse, M. T., Johnson, C., Carslaw, K. S., and Dalvi, M.: Impact of the modal aerosol scheme GLOMAP-mode on aerosol forcing in the Hadley Centre Global Environmental Model, Atmos. Chem. Phys., 13, 3027–3044, 10.5194/acp-13-3027-2013, 2013.Boers, R., Acarreta, J. R., and Gras, J. L.: Satellite monitoring of the first indirect aerosol effect: Retrieval of the droplet concentration of water clouds, J. Geophys. Res., 111, D22208, 10.1029/2005JD006838, 2006.Boutle, I. A. and Morcrette, C. J.: Parametrization of area cloud fraction, Atmos. Sci. Lett., 11, 283–289, 10.1002/asl.293, 2010.Boutle, I. A., Eyre, J. E. J., and Lock, A. P.: Seamless stratocumulus simulation across the turbulent gray zone, Mon. Weather Rev., 142, 1655–1668, 10.1175/MWR-D-13-00229.1, 2014.Brown, A. R.: The sensitivity of large-eddy simulations of shallow cumulus convection to resolution and subgrid model, Q. J. Roy. Meteor. Soc., 125, 469–482, 10.1002/qj.49712555405, 1999.Bush, M., Allen, T., Bain, C., Boutle, I., Edwards, J., Finnenkoetter, A., Franklin, C., Hanley, K., Lean, H., Lock, A., Manners, J., Mittermaier, M., Morcrette, C., North, R., Petch, J., Short, C., Vosper, S., Walters, D., Webster, S., Weeks, M., Wilkinson, J., Wood, N., and Zerroukat, M.: The first Met Office Unified Model–JULES Regional Atmosphere and Land configuration, RAL1, Geosci. Model Dev., 13, 1999–2029, 10.5194/gmd-13-1999-2020, 2019.CEDA (Centre for Environmental Data Analysis): The Natural Environment Research Council's Data Repository for Atmospheric Science and Earth Observation, available at: http://archive.ceda.ac.uk/, last access: August 2019.Clark, T. L.: Numerical Modeling of the Dynamics and Microphysics of Warm Cumulus Convection, J. Atmos. Sci., 30, 857–878, 10.1175/1520-0469(1973)030<0857:NMOTDA>2.0.CO;2, 1973.Costantino, L. and Bréon, F.-M.: Aerosol indirect effect on warm clouds over South-East Atlantic, from co-located MODIS and CALIPSO observations, Atmos. Chem. Phys., 13, 69–88, 10.5194/acp-13-69-2013, 2013.
Dearden, C.: Investigating the simulation of cloud microphysical processes in numerical models using a one-dimensional dynamical framework, Atmos. Sci. Lett., 10, 207–214, 10.1002/asl.239, 2009.Fan, J., Leung, L. R., Li, Z., Morrison, H., Chen, H., Zhou, Y., Qian, Y., and Wang, Y.: Aerosol impacts on clouds and precipitation in eastern China: Results from bin and bulk microphysics, J. Geophys. Res., 117, D00K36, 10.1029/2011JD016537, 2012.Fan, J., Rosenfeld, D., Zhang, Y., Giangrande, S. E., Li, Z., Machado, L. A. T., Martin, S. T., Yang, Y., Wang, J., Artaxo, P., Barbosa, H. M. J., Braga, R. C., Comstock, J. M., Feng, Z., Gao, W., Gomes, H. B., Mei, F., Pöhlker, C., Pöhlker, M. L., Pöschl, U., and de Souza, R. A. F.: Substantial convection and precipitation enhancements by ultrafine aerosol
particles, Science, 359, 411–418, 10.1126/science.aan8461, 2018.Field, P. R. and Furtado, K.: How biased Is aircraft cloud sampling?, J. Atmos. Ocean. Tech., 33, 185–189, 10.1175/JTECH-D-15-0148.1, 2016.Fountoukis, C. and Nenes, A.: Continued development of a cloud droplet formation parameterization for global climate models, J. Geophys. Res., 110, D11212, 10.1029/2004JD005591, 2005.Freud, E., Rosenfeld, D., and Kulkarni, J. R.: Resolving both entrainment-mixing and number of activated CCN in deep convective clouds, Atmos. Chem. Phys., 11, 12887–12900, 10.5194/acp-11-12887-2011, 2011.Fridlind, A. M., Ackerman, A. S., Jensen, E. J., Heymsfield, A. J., Poellot, M. R., Stevens, D. E., Wang, D., Miloshevich, L. M., Baumgardner, D., Lawson, R. P., Wilson, J. C., Flagan, R. C., Seinfeld, J. H., Jonsson, H. H., VanReken, T. M., Varutbangkul, V., and Rissman, T. A.: Evidence for the predominance of mid-tropospheric aerosols as subtropical anvil cloud nuclei, Science, 304, 718–722, 10.1126/science.1094947, 2004.Furtado, K., Field, P. R., Luo, Y., Liu, X., Guo, Z., Zhou, T., Shipway, B. J., Hill, A. A., and Wilkinson, J. M.: Cloud microphysical factors affecting simulations of deep convection during the presummer rainy season in Southern
China, J. Geophys. Res.-Atmos., 123, 10477–10505, 10.1029/2017JD028192, 2018.Gettelman, A. and Morrison, H.: Advanced two-Moment bulk microphysics for global models. Part I: off-Line tests and comparison with other schemes, J. Climate, 28, 1268–1287, 10.1175/JCLI-D-14-00102.1, 2015.Ghan, S. J., Abdul-Razzak, H., Nenes, A., Ming, Y., Liu, X., Ovchinnikov, M., Shipway, B., Meskhidze, N., Xu, J., and Shi, X.: Droplet nucleation: physically-based parameterizations and comparative evaluation, J. Adv. Model. Earth Sy., 3, M10001, 10.1029/2011MS000074, 2011.Gong, S. L.: A parameterization of sea-salt aerosol source function for sub- and super-micron particles, Global Biogeochem. Cy., 17, 1097, 10.1029/2003GB002079, 2003.Gordon, H., Field, P. R., Abel, S. J., Dalvi, M., Grosvenor, D. P., Hill, A. A., Johnson, B. T., Miltenberger, A. K., Yoshioka, M., and Carslaw, K. S.: Large simulated radiative effects of smoke in the south-east Atlantic, Atmos. Chem. Phys., 18, 15261–15289, 10.5194/acp-18-15261-2018, 2018.Grabowski, W. W.: Representation of turbulent mixing and buoyancy reversal in bulk cloud models, J. Atmos. Sci., 64, 3666–3680, 10.1175/JAS4047.1, 2007.Grabowski, W. W. and Jarecka, D.: Modeling condensation in shallow nonprecipitating convection, J. Atmos. Sci., 72, 4661–4679, 10.1175/JAS-D-15-0091.1, 2015.Grabowski, W. W. and Morrison, H.: Modeling condensation in deep convection, J. Atmos. Sci., 74, 2247–2267, 10.1175/JAS-D-16-0255.1, 2017.Grabowski, W. W. and Wang, L.-P.: Growth of cloud droplets in a turbulent environment, Annu. Rev. Fluid Mech., 45, 293–324, 10.1146/annurev-fluid-011212-140750, 2013.Grosvenor, D. P. and Wood, R.: The effect of solar zenith angle on MODIS cloud optical and microphysical retrievals within marine liquid water clouds, Atmos. Chem. Phys., 14, 7291–7321, 10.5194/acp-14-7291-2014, 2014.Grosvenor, D. P., Field, P. R., Hill, A. A., and Shipway, B. J.: The relative importance of macrophysical and cloud albedo changes for aerosol-induced radiative effects in closed-cell stratocumulus: insight from the modelling of a case study, Atmos. Chem. Phys., 17, 5155–5183, 10.5194/acp-17-5155-2017, 2017.Guichard, F. and Couvreux, F.: A short review of numerical cloud-resolving models, Tellus A, 69, 1373578, 10.1080/16000870.2017.1373578, 2017.Hall, W. D.: A detailed microphysical model within a two-dimensional dynamic framework: model description and preliminary results, J. Atmos. Sci., 37, 2486–2507, 10.1175/1520-0469(1980)037<2486:ADMMWA>2.0.CO;2, 1980.Heymsfield, A. J., Dye, J. E., and Biter, C. J.: Overestimates of entrainment from wetting of aircraft temperature sensors in cloud, J. Appl. Meteorol., 18, 92–95, 10.1175/1520-0450(1979)018<0092:OOEFWO>2.0.CO;2, 1979.Heymsfield, A. J., Bansemer, A., Heymsfield, G., and Fierro, A. O.: Microphysics of maritime tropical convective updrafts at temperatures from -20∘ to -60∘, J. Atmos. Sci., 66, 3530–3562, 10.1175/2009JAS3107.1, 2009.
Hill, A. A., Dobbie, S., and Yin, Y.: The impact of aerosols on non-precipitating marine stratocumulus. I: Model description and prediction of the indirect effect, Q. J. Roy. Meteor. Soc., 134, 1143–1154, 2008.Hill, A. A., Shipway, B. J., and Boutle, I. A.: How sensitive are aerosol-precipitation interactions to the warm rain representation?, J. Adv. Model. Earth Sy., 7, 987–1004, 10.1002/2014MS000422, 2015.Honnert, R., Masson, V., and Couvreux, F.: A diagnostic for evaluating the representation of turbulence in atmospheric models at the kilometric scale, J. Atmos. Sci., 68, 3112–3131, 10.1175/JAS-D-11-061.1, 2011.Huffman, G., Bolvin, D., Braithwaite, D., Hsu, K., Joyce, R., and Xie, P.: Integrated Multi-satellitE Retrievals for GPM (IMERG), version 4.4, NASA's Precipitation Processing Center, 10.5067/GPM/IMERG/3B-HH/05, 2014.Ichoku, C. and Ellison, L.: Global top-down smoke-aerosol emissions estimation using satellite fire radiative power measurements, Atmos. Chem. Phys., 14, 6643–6667, 10.5194/acp-14-6643-2014, 2014.Khain, A. and Lynn, B.: Simulation of a supercell storm in clean and dirty atmosphere using weather research and forecast model with spectral bin microphysics, J. Geophys. Res., 114, D19209, 10.1029/2009JD011827, 2009.Khain, A. P., Phillips, V., Benmoshe, N., and Pokrovsky, A.: The Role of Small Soluble Aerosols in the Microphysics of Deep Maritime Clouds, J. Atmos. Sci., 69, 2787–2807, 10.1175/2011JAS3649.1, 2012.Khain, A. P., Beheng, K. D., Heymsfield, A., Korolev, A., Krichak, S. O., Levin, Z., Pinsky, M., Phillips, V., Prabhakaran, T., Teller, A., van den Heever, S. C., and Yano, J.-I.: Representation of microphysical processes in cloud-resolving models: Spectral (bin) microphysics versus bulk parameterization, Rev. Geophys., 53, 247–322, 10.1002/2014RG000468, 2015.Khairoutdinov, M. and Kogan, Y.: A New cloud physics parameterization in a Large-Eddy Simulation model of marine stratocumulus, Mon. Weather Rev., 128, 229–243, 10.1175/1520-0493(2000)128<0229:ANCPPI>2.0.CO;2, 2000.Kogan, Y. L. and Martin, W. J.: Parameterization of bulk condensation in numerical cloud models, J. Atmos. Sci., 51, 1728–1739, 10.1175/1520-0469(1994)051<1728:POBCIN>2.0.CO;2, 1994.Korolev, A., Strapp, J. W., Isaac, G. A., and Emery, E.: Improved airborne
hot-wire measurements of ice water content in clouds, J. Atmos. Ocean. Tech., 30, 2121–2131, 10.1175/JTECH-D-13-00007.1, 2013.Korolev, A. V.: The influence of supersaturation fluctuations on droplet size spectra formation, J. Atmos. Sci., 52, 3620–3634, 10.1175/1520-0469(1995)052<3620:TIOSFO>2.0.CO;2, 1995.Korolev, A. V. and Mazin, I. P.: Supersaturation of water vapor in clouds, J. Atmos. Sci., 60, 2957–2974, 10.1175/1520-0469(2003)060<2957:SOWVIC>2.0.CO;2, 2003.Köhler, H.: The nucleus in and the growth of hygroscopic droplets, T. Faraday Soc., 32, 1152–1161, 10.1039/TF9363201152, 1936.Kuhlbrodt, T., Jones, C. G., Sellar, A., Storkey, D., Blockley, E., Stringer, M., Hill, R., Graham, T., Ridley, J., Blaker, A., Calvert, D., Copsey, D., Ellis, R., Hewitt, H., Hyder, P., Ineson, S., Mulcahy, J., Siahaan, A., and Walton, J.: The low‐resolution version of HadGEM3 GC3.1: Development and evaluation for global climate, J. Adv. Model. Earth Sy., 10, 2865–2888, 10.1029/2018MS001370, 2020.Laaksonen, A., Vesala, T., Kulmala, M., Winkler, P. M., and Wagner, P. E.: Commentary on cloud modelling and the mass accommodation coefficient of water, Atmos. Chem. Phys., 5, 461–464, 10.5194/acp-5-461-2005, 2005.Lamarque, J.-F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z., Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D., Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M., Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and van Vuuren, D. P.: Historical (1850–2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols: methodology and application, Atmos. Chem. Phys., 10, 7017–7039, 10.5194/acp-10-7017-2010, 2010.Lebo, Z. J., Morrison, H., and Seinfeld, J. H.: Are simulated aerosol-induced effects on deep convective clouds strongly dependent on saturation adjustment?, Atmos. Chem. Phys., 12, 9941–9964, 10.5194/acp-12-9941-2012, 2012.Lock, A. P., Brown, A. R., Bush, M. R., Martin, G. M., and Smith, R. N. B.: A New Boundary Layer Mixing Scheme. Part I: Scheme Description and Single-Column Model Tests, Mon. Weather Rev., 128, 3187–3199, 10.1175/1520-0493(2000)128<3187:ANBLMS>2.0.CO;2, 2000.Lohmann, U.: Possible aerosol effects on ice clouds via contact nucleation, J. Atmos. Sci., 59, 647–656, 10.1175/1520-0469(2001)059<0647:PAEOIC>2.0.CO;2, 2002.Lu, Z., Liu, X., Zhang, Z., Zhao, C., Meyer, K., Rajapakshe, C., Wu, C., Yang, Z., and Penner, J. E.: Biomass smoke from southern Africa can significantly enhance the brightness of stratocumulus over the southeastern Atlantic Ocean, P. Natl. Acad. Sci. USA, 115, 2924–2929, 10.1073/pnas.1713703115, 2018.Malavelle, F. F., Haywood, J. M., Field, P. R., Hill, A. A., Abel, S. J., Lock, A. P., Shipway, B. J., and McBeath, K.: A method to represent subgrid-scale updraft velocity in kilometer-scale models: Implication for aerosol activation, J. Geophys. Res.-Atmos., 119, 4149–417, 10.1002/2013JD021218, 2014.Mann, G. W., Carslaw, K. S., Spracklen, D. V., Ridley, D. A., Manktelow, P. T., Chipperfield, M. P., Pickering, S. J., and Johnson, C. E.: Description and evaluation of GLOMAP-mode: a modal global aerosol microphysics model for the UKCA composition-climate model, Geosci. Model Dev., 3, 519–551, 10.5194/gmd-3-519-2010, 2010.Martin, G. M., Johnson, D. W., and Spice, A.: The measurement and parameterization of effective radius of droplets in warm stratocumulus clouds, J. Atmos. Sci., 51, 1823–1842, 10.1175/1520-0469(1994)051<1823:TMAPOE>2.0.CO;2, 1994.Met Office: CLARIFY: in-situ airborne observations by the FAAM BAE-146 aircraft, Met Office, Natural Environment Research Council, Facility for Airborne Atmospheric Measurements, Centre for Environmental Data Analysis, available at: https://catalogue.ceda.ac.uk/uuid/38ab7089781a4560b067dd6c20af3769 (last access: August 2019), 2017.Met Office: Unified Model, available at: http://www.metoffice.gov.uk/research/modelling-systems/unified-model, last access: 1 June 2020.Miltenberger, A. K., Field, P. R., Hill, A. A., Rosenberg, P., Shipway, B. J., Wilkinson, J. M., Scovell, R., and Blyth, A. M.: Aerosol–cloud interactions in mixed-phase convective clouds – Part 1: Aerosol perturbations, Atmos. Chem. Phys., 18, 3119–3145, 10.5194/acp-18-3119-2018, 2018.Ming, Y., Ramaswamy, V., Donner, L. J., Phillips, V. T. J., Klein, S. A., Ginoux, P. A., and Horowitz, L. W.: Modeling the Interactions between Aerosols and Liquid Water Clouds with a Self-Consistent Cloud Scheme in a General Circulation Model, J. Atmos. Sci., 64, 1189–1209, 10.1175/JAS3874.1, 2007.Morrison, H. and Gettelman, A.: A new two-moment bulk stratiform cloud microphysics scheme in the Community Atmosphere Model, Version 3 (CAM3). Part I: description and numerical tests, J. Climate, 21, 3642–3659, 10.1175/2008JCLI2105.1, 2008.Morrison, H. and Grabowski, W. W.: Modeling supersaturation and subgrid-scale mixing with two-moment bulk warm microphysics, J. Atmos. Sci., 65, 792–812, 10.1175/2007JAS2374.1, 2008.Morrison, H., Curry, J. A., and Khvorostyanov, V. I.: A New Double-Moment Microphysics Parameterization for Application in Cloud and Climate Models. Part I: Description, J. Atmos. Sci., 62, 1665–1677,
10.1175/JAS3446.1, 2005.Morrison, H., Thompson, G., and Tatarskii, V.: Impact of cloud microphysics on the development of trailing stratiform precipitation in a simulated squall line: comparison of one- and two-moment schemes, Mon. Weather Rev., 137, 991–1007, 10.1175/2008MWR2556.1, 2009.Mulcahy, J. P., Jones, C., Sellar, A., Johnson, B., Boutle, I. A., Jones, A., Andrews, T., Rumbold, S. T., Mollard, J., Bellouin, N., Johnson, C. E., Williams, K. D., Grosvenor, D. P., and McCoy, D. T.: Improved aerosol processes and effective radiative forcing in HadGEM3 and UKESM1, J. Adv. Model. Earth Sy., 10, 2786–2805, 10.1029/2018MS001464, 2018.Mulcahy, J. P., Johnson, C., Jones, C. G., Povey, A. C., Scott, C. E., Sellar, A., Turnock, S. T., Woodhouse, M. T., Abraham, N. L., Andrews, M. B., Bellouin, N., Browse, J., Carslaw, K. S., Dalvi, M., Folberth, G. A., Glover, M., Grosvenor, D., Hardacre, C., Hill, R., Johnson, B., Jones, A., Kipling, Z., Mann, G., Mollard, J., O'Connor, F. M., Palmieri, J., Reddington, C., Rumbold, S. T., Richardson, M., Schutgens, N. A. J., Stier, P., Stringer, M., Tang, Y., Walton, J., Woodward, S., and Yool, A.: Description and evaluation of aerosol in UKESM1 and HadGEM3-GC3.1 CMIP6 historical simulations, Geosci. Model Dev. Discuss., 10.5194/gmd-2019-357, in review, 2020.Nenes, A. and Seinfeld, J. H.: Parameterization of cloud droplet formation in global climate models, J. Geophys. Res., 108, 4415, 10.1029/2002JD002911, 2003.Peers, F., Francis, P., Fox, C., Abel, S. J., Szpek, K., Cotterell, M. I., Davies, N. W., Langridge, J. M., Meyer, K. G., Platnick, S. E., and Haywood, J. M.: Observation of absorbing aerosols above clouds over the south-east Atlantic Ocean from the geostationary satellite SEVIRI – Part 1: Method description and sensitivity, Atmos. Chem. Phys., 19, 9595–9611, 10.5194/acp-19-9595-2019, 2019.Petersen, G. N. and Renfrew, I. A.: Aircraft-based observations of air–sea fluxes over Denmark Strait and the Irminger Sea during high wind speed conditions, Q. J. Roy. Meteor. Soc., 135, 2030–2045, 10.1002/qj.355, 2009.Petters, M. D. and Kreidenweis, S. M.: A single parameter representation of hygroscopic growth and cloud condensation nucleus activity, Atmos. Chem. Phys., 7, 1961–1971, 10.5194/acp-7-1961-2007, 2007.Phillips, V. T. J., Donner, L. J., and Garner, S. T.: Nucleation processes in deep convection simulated by a cloud-system-resolving model with double-moment bulk microphysics, J. Atmos. Sci., 64, 738–761, 10.1175/JAS3869.1, 2007.Pinsky, M. B. and Khain, A. P.: Effects of in-cloud nucleation and turbulence on droplet spectrum formation in cumulus clouds, Q. J. Roy. Meteor. Soc., 128, 501–533, 10.1256/003590002321042072, 2002.Planche, C., Mann, G. W., Carslaw, K. S., Dalvi, M., Marsham, J. H., and Field, P. R.: Spatial and temporal CCN variations in convection-permitting aerosol microphysics simulations in an idealised marine tropical domain, Atmos. Chem. Phys., 17, 3371–3384, 10.5194/acp-17-3371-2017, 2017.Platnick, S., King, M., Wind, G., Ackerman, S., Menzel, P., and Frey, R.: MODIS Atmosphere L2 Cloud Product (06_L2), NASA MODIS Adaptive Processing System, Goddard Space Flight Center, USA, 10.5067/MODIS/MYD06_L2.006, 2015.Politovich, M. K. and Cooper, W. A.: Variability of the supersaturation in cumulus clouds, J. Atmos. Sci., 45, 1651–1664, 10.1175/1520-0469(1988)045<1651:VOTSIC>2.0.CO;2, 1988.Possner, A., Zubler, E., Lohmann, U., and Schär, C.: The resolution dependence of cloud effects and ship-induced aerosol-cloud interactions in marine stratocumulus, J. Geophys. Res.-Atmos., 121, 4810–4829, 10.1002/2015JD024685, 2016.
Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and Precipitation, Springer, Dordrecht, 1997.Rosenberg, P. D., Dean, A. R., Williams, P. I., Dorsey, J. R., Minikin, A., Pickering, M. A., and Petzold, A.: Particle sizing calibration with refractive index correction for light scattering optical particle counters and impacts upon PCASP and CDP data collected during the Fennec campaign, Atmos. Meas. Tech., 5, 1147–1163, 10.5194/amt-5-1147-2012, 2012.Rothenberg, D. and Wang, C.: Metamodeling of Droplet Activation for Global Climate Models, J. Atmos. Sci., 73, 1255–1272, 10.1175/JAS-D-15-0223.1, 2016.Saleeby, S. M. and Cotton, W. R.: A large-droplet mode and prognostic number concentration of cloud droplets in the Colorado State University Regional Atmospheric Modeling System (RAMS). Part I: module descriptions and supercell test simulations, J. Appl. Meteorol., 43, 182–195, 10.1175/1520-0450(2004)043<0182:ALMAPN>2.0.CO;2, 2004.Segal, Y., Pinsky, M., Khain, A., and Erlick, C.: Thermodynamic factors influencing bimodal spectrum formation in cumulus clouds, Atmos. Res., 66, 43–64, 10.1016/S0169-8095(02)00172-2, 2003.Seiki, T., Kodama, C., Noda, A. T., and Satoh, M.: Improvement in Global Cloud-System-Resolving Simulations by Using a Double-Moment Bulk Cloud Microphysics Scheme, J. Climate, 28, 2405–2419, 10.1175/JCLI-D-14-00241.1, 2015.Shinozuka, Y., Saide, P. E., Ferrada, G. A., Burton, S. P., Ferrare, R., Doherty, S. J., Gordon, H., Longo, K., Mallet, M., Feng, Y., Wang, Q., Cheng, Y., Dobracki, A., Freitag, S., Howell, S. G., LeBlanc, S., Flynn, C., Segal-Rosenhaimer, M., Pistone, K., Podolske, J. R., Stith, E. J., Bennett, J. R., Carmichael, G. R., da Silva, A., Govindaraju, R., Leung, R., Zhang, Y., Pfister, L., Ryoo, J.-M., Redemann, J., Wood, R., and Zuidema, P.: Modeling the smoky troposphere of the southeast Atlantic: a comparison to ORACLES airborne observations from September of 2016, Atmos. Chem. Phys. Discuss., 10.5194/acp-2019-678, in review, 2019.Shipway, B. J. and Hill, A. A.: Diagnosis of systematic differences between multiple parametrizations of warm rain microphysics using a kinematic framework, Q. J. Roy. Meteor. Soc., 138, 2196–2211, 10.1002/qj.1913, 2012.Smith, R. N. B.: A scheme for predicting layer clouds and their water content in a general circulation model, Q. J. Roy. Meteor. Soc., 116, 435–460, 10.1002/qj.49711649210, 1990.Squires, P.: The Growth of Cloud Drops by Condensation. I. General Characteristics, Aust. J. Sci. Res. Ser. A, 5, 59–86, 10.1071/CH9520059, 1952.Stevens, B., Feingold, G., Cotton, W. R., and Walko, R. L.: Elements of the Microphysical Structure of Numerically Simulated Nonprecipitating Stratocumulus, J. Atmos. Sci., 53, 980–1006, 10.1175/1520-0469(1996)053<0980:EOTMSO>2.0.CO;2, 1996.Stevens, B., Bony, S., Brogniez, H., Hentgen, L., Hohenegger, C., Kiemle, C., L'Ecuyer, T. S., Naumann, A. K., Schulz, H., Siebesma, P. A., Vial, J., Winker, D. M., and Zuidema, P.: Sugar, gravel, fish and flowers: Mesoscale cloud patterns in the trade winds, Q. J. Roy. Meteor. Soc., 146: 141–152, 10.1002/qj.3662, 2019.Stevens, R. G., Loewe, K., Dearden, C., Dimitrelos, A., Possner, A., Eirund, G. K., Raatikainen, T., Hill, A. A., Shipway, B. J., Wilkinson, J., Romakkaniemi, S., Tonttila, J., Laaksonen, A., Korhonen, H., Connolly, P., Lohmann, U., Hoose, C., Ekman, A. M. L., Carslaw, K. S., and Field, P. R.: A model intercomparison of CCN-limited tenuous clouds in the high Arctic, Atmos. Chem. Phys., 18, 11041–11071, 10.5194/acp-18-11041-2018, 2018.Telford, P. J., Braesicke, P., Morgenstern, O., and Pyle, J. A.: Technical Note: Description and assessment of a nudged version of the new dynamics Unified Model, Atmos. Chem. Phys., 8, 1701–1712, 10.5194/acp-8-1701-2008, 2008.
Thompson, G.: Advances in a microphysics parameterization to predict
supercooled liquid water and application to aircraft icing, PhD thesis, University of Colorado, Boulder, United States, 2016.Twomey, S.: The nuclei of natural cloud formation part II: The supersaturation in natural clouds and the variation of cloud droplet concentration, Geofisica pura e applicata, 43, 243–249, 10.1007/BF01993560, 1959.Vance, A. K., Abel, S. J., Cotton, R. J., and Woolley, A. M.: Performance of WVSS-II hygrometers on the FAAM research aircraft, Atmos. Meas. Tech., 8, 1617–1625, 10.5194/amt-8-1617-2015, 2015.Walters, D., Baran, A. J., Boutle, I., Brooks, M., Earnshaw, P., Edwards, J., Furtado, K., Hill, P., Lock, A., Manners, J., Morcrette, C., Mulcahy, J., Sanchez, C., Smith, C., Stratton, R., Tennant, W., Tomassini, L., Van Weverberg, K., Vosper, S., Willett, M., Browse, J., Bushell, A., Carslaw, K., Dalvi, M., Essery, R., Gedney, N., Hardiman, S., Johnson, B., Johnson, C., Jones, A., Jones, C., Mann, G., Milton, S., Rumbold, H., Sellar, A., Ujiie, M., Whitall, M., Williams, K., and Zerroukat, M.: The Met Office Unified Model Global Atmosphere 7.0/7.1 and JULES Global Land 7.0 configurations, Geosci. Model Dev., 12, 1909–1963, 10.5194/gmd-12-1909-2019, 2019.Wang, H., Easter, R. C., Rasch, P. J., Wang, M., Liu, X., Ghan, S. J., Qian, Y., Yoon, J.-H., Ma, P.-L., and Vinoj, V.: Sensitivity of remote aerosol distributions to representation of cloud–aerosol interactions in a global climate model, Geosci. Model Dev., 6, 765–782, 10.5194/gmd-6-765-2013, 2013.West, R. E. L., Stier, P., Jones, A., Johnson, C. E., Mann, G. W., Bellouin, N., Partridge, D. G., and Kipling, Z.: The importance of vertical velocity variability for estimates of the indirect aerosol effects, Atmos. Chem. Phys., 14, 6369–6393, 10.5194/acp-14-6369-2014, 2014.Wilson, D. R. and Ballard, S. P.: A microphysically based precipitation scheme for the UK Meteorological Office Unified Model, Q. J. Roy. Meteor. Soc., 125, 1607–1636, 10.1002/qj.49712555707, 1999.Wilson, D. R., Bushell, A. C., Kerr-Munslow, A. M., Price, J. D., and Morcrette, C. J.: PC2: A prognostic cloud fraction and condensation scheme. I: Scheme description, Q. J. Roy. Meteor. Soc., 134, 2093–2107, 10.1002/qj.333, 2008.Woodward, S.: Modeling the atmospheric life cycle and radiative impact of mineral dust in the Hadley Centre climate model, J. Geophys. Res., 106, 18155–18166, 10.1029/2000JD900795, 2001.
Wyngaard, J. C.: Toward Numerical Modeling in the “Terra Incognita”, J. Atmos. Sci., 61, 1816–1826, 10.1175/1520-0469(2004)061<1816:TNMITT>2.0.CO;2, 2004.Yang, Q., Easter, R. C., Campuzano-Jost, P., Jimenez, J. L., Fast, J. D., Ghan, S. J., Wang, H., Berg, L. K., Barth, M. C., Liu, Y., Shrivastava, M. B., Singh, B., Morrison, H., Fan, J., Ziegler, C. L., Bela, M., Apel, E., Diskin, G. S., Mikoviny, T., and Wisthaler, A.: Aerosol transport and wet scavenging in deep convective clouds: A case study and model evaluation using a multiple passive tracer analysis approach, J. Geophys. Res.-Atmos., 120, 8448–8468, 10.1002/2015JD023647, 2015.Zuidema, P., Redemann, J., Haywood, J., Wood, R., Piketh, S., Hipondoka, M., and Formenti, P.: Smoke and Clouds above the Southeast Atlantic: Upcoming Field Campaigns Probe Absorbing Aerosol’s Impact on Climate, B. Am. Meteorol. Soc., 97, 1131–1135, 10.1175/BAMS-D-15-00082.1, 2016.