Introduction
Aerosol particles play an important indirect role in the climate system by
modifying cloud micro- and macrophysical properties, which is referred to as
aerosol–cloud interactions . An increase in
aerosols supplies more numerous cloud condensation nuclei, resulting in
numerous and smaller cloud droplets leading to brighter clouds, which is
known as the “albedo effect” . Smaller cloud droplets
suppress the onset of precipitation in warm clouds due to the less efficient
collision–coalescence process, resulting in a longer cloud lifetime, which
is known as the “lifetime effect” . There has been much
discussion about the climatic impacts of aerosol-induced modulation of water
clouds, which are particularly sensitive to aerosol perturbations
e.g.,.
However, quantitative estimates of radiative forcing with regard to
aerosol–cloud–precipitation–climate interactions remain uncertain, as
reported in the Fifth Assessment Report of the Intergovernmental Panel on
Climate Change (2013).
One of the most important factors that quantify the magnitude of
aerosol–cloud interactions is the response of the cloud liquid water path
(LWP) to aerosol perturbations. This factor also characterizes aerosol
impacts on the global hydrological cycle through its representation of the
aerosol effect on precipitation efficiency. This effect is represented in
general circulation models (GCMs) as aerosol-induced changes in rainwater
production from cloud water, which are parameterized with a bulk microphysics
as the so-called autoconversion process. The water conversion rate by this
process (Paut) is generally given as a function of the liquid water
content (Lc) and cloud droplet number concentration (Nc) as
Paut∼Lcα×Nc-β,
where α and β are prescribed constants
e.g.,. The Nc is then
somehow related to the aerosol number concentration (Na). In GCMs, Eq. () provides the only pathway through which aerosols
modulate precipitation formation and, thus, the cloud lifetime. Note that
GCMs also partly include the opposing processes (decreasing LWP due to
enhancement in evaporation) via the so-called direct and semi-direct effect
e.g.,. Given that rainwater production is
always suppressed with increasing Nc according to Eq. (), GCMs tend to increase the LWP uniformly with increasing
Nc for stratiform clouds.
On the other hand, some observational studies have shown two pathways of LWP
responses to perturbed aerosols, i.e., both increasing and decreasing
tendencies of LWP with increasing aerosols
; the mechanisms for these opposing
responses cannot be understood by a simple microphysical argument alone, but
are likely to relate to macrophysical and meteorological factors as well
e.g.,.
found that there are processes that modify the cloud
geometric thickness to aerosol perturbations in such a way that cancels the
aerosol indirect effect at sufficiently long timescales. Such a compensation
mechanism is currently considered one of the “buffering effects”
, which generate the opposite result to the original
hypotheses of cloud albedo and lifetime effects, for the cloud system as a
whole . Despite its critical
importance to accurate climate simulations, the operation of this mechanism
at the global scale remains poorly understood.
To determine the mechanisms involved in the competition between the
“lifetime effect” and “buffering effect”, the complexity in aerosol
effects on clouds needs to be untangled at a fundamental process-level. For
this purpose, GCMs should be evaluated extensively against observations in
the context of their process representations, which are key to the
aerosol–cloud–precipitation interaction.
In this study, we analyze results from both GCM and A-Train data, with a
particular focus on their discrepancies in the key indices of aerosol–cloud
interactions relating to fundamental processes. The factors examined are the
susceptibilities of cloud optical thickness (τc), droplet effective
radius (re), and LWP to Nc. To focus directly on the cloud physical
parameters, we use Nc as an aerosol proxy rather than Na
e.g.. A satellite-based study by reported
that cloud susceptibilities show similar results whether aerosol index,
aerosol optical depth (AOD), or Nc are applied as an aerosol proxy (see
their supplementary information). Given the fundamental relationship of
τc ∝ LWP / re, the susceptibilities are related as follows
:
dlnτcdlnNc=-dlnredlnNc+dlnLWPdlnNc,
where the first and second terms on the right side of Eq. ()
represent the “albedo effect” and the “lifetime effect”, respectively.
Equation () has the advantage that it can quantify the
contributions from the two effects that determine the aerosol impact on cloud
radiative properties. As discussed here and also in recent studies
, the two terms in Eq. () are
related to representations of different processes. This approach makes it
easier to understand the mechanisms that determine the resultant magnitude of
aerosol indirect forcing in the context of relevant processes
.
The aim of the study is to clarify the fundamental source of uncertainty in
process representations of aerosol–cloud–precipitation interactions in GCMs
for stratiform and shallow cumulative warm clouds (excluding deep convective
thick clouds or ice clouds; see Sect. 2). Given that the aerosol–cloud
interaction processes are also influenced by both macrophysics (e.g.,
environmental conditions, dynamical regime, cloud type) and microphysics, we
also place an emphasis on the importance of macrophysics
e.g.,.
Data
MIROC–SPRINTARS
A global climate model, Model for Interdisciplinary Research On Climate
(MIROC) version 5.2 , was used in this study. The
interactions of the main tropospheric aerosols (i.e., black carbon, organic
matter, soil dust, sea salt, sulfate, and the precursor gases of sulfate)
with cloud–precipitation microphysics and radiation–climate effects are
incorporated in the aerosol module, Spectral Radiation-Transport Model for
Aerosol Species (SPRINTARS; ),
which is coupled with MIROC (MIROC–SPRINTARS).
The cloud macro- and microphysics framework in MIROC–SPRINTARS is based on a
prognostic large-scale condensation scheme, which explicitly considers
subgrid-scale variability of clouds . This probability distribution function (PDF)-based
prognostic cloud scheme couples with the ice microphysics scheme proposed by
. MIROC–SPRINTARS treats both cloud droplets and ice
crystals as a two-moment bulk microphysics scheme . The
nucleation of cloud droplets is parameterized by the scheme of
, and the process of cloud-to-rain water conversion
is diagnosed based on the autoconversion scheme. Rainwater
is not a prognostic variable in the current version of MIROC–SPRINTARS.
We extracted warm-phase low clouds (> 273.15 K in whole cloud
layers) from every 6 h instantaneous output for 5 full years; as a
result, 1 595 753 warm cloud samples were obtained. The horizontal and
vertical resolutions were T42 (approximately 2.8∘ × 2.8∘
in latitude and longitude) and 20 layers, respectively. A more detailed
description of the model and its settings are documented in
.
CloudSat and MODIS
We used the synergistic satellite data sets of the CloudSat and MODIS, which
are both part of the A-Train constellation .
The data products, 2B-TAU , 2B-GEOPROF
, and ECMWF-AUX were used for the
period June 2006 to April 2011, i.e., a total of 5 full years. This
facilitated the construction of stable statistics with a horizontal
resolution (2.5∘ grid boxes in this study) close to the GCM output. We
defined the cloud layer as where the cloud mask value is greater than 30 from the
2B-GEOPROF product, which means a good, or strong, echo with high-confidence
detection . The analysis was restricted to single-layer
water clouds; in total, 7 872 426 cloud samples were obtained.
The LWP was derived from the MODIS-retrieved optical thickness and effective
radius using the following equation for an adiabatically stratified cloud
:
LWP=59τcre.
Nc was also calculated based on an adiabatic assumption
as
Nc=2B3Γeff1/2LWP1/2re3,
where B =(3/4πρw)1/3= 0.0620 kg-1/3m,
ρw is the density of liquid water, and Γeff is the
adiabatic rate (gm-3km-1) of increase in the liquid water
content with height (see also for more details of
derivation). We note that satellite data inherently include uncertainties
stemming from retrieval assumptions, which are not replicated in the model
output. Although it could be a part of the reason for discrepancies between the
model and observations, this would mostly be canceled when susceptibilities
of cloud and precipitation to aerosol loading are evaluated by a logarithmic
form. This study applies the uncertainty thresholds of < 5 and < 1 µm for τc and re from the 2B-TAU product,
respectively, which contributes to the reduction of the retrieval uncertainty of
Nc described above as much as possible .
To examine the cloud-to-rain conversion process, the conversion rate
(Pconv) contributed by both autoconversion (collision–coalescence of
cloud droplets) and accretion (collision of cloud droplets by raindrops) was
derived from the approximation suggested by . This method
is established by the continuous collection equation
using observed drop size distributions. Pconv was estimated from MODIS
LWP and CloudSat mean cloud-layer radar reflectivity Z‾ as
Pconv=c1LWPZ‾H[Z-Zc],
where c1=κ2/26 is a coefficient from the collection kernel
with κ2=1.9×1011 cm-3 s-1 and sixth moment factor with radar reflectivity. H[Z-Zc] is
the Heaviside step function to exclude the cases that are less than the critical
radar threshold Zc of -15 dBZe for which the conversion
process is negligible . Although this formulation is
based on marine stratocumulus cases from DYCOMS-II measurements
, it is applicable for global analysis to study
aerosol–cloud interactions in drizzling
light-rain cases (Z‾< 0dBZe). The
parameterization and assumptions used in this method (Eq. )
are also valid for comparison between observations and model simulation
. This brings valuable understanding for microphysical
conversion processes and its timescales, which matches the scope of our
study.
Lower-tropospheric stability (LTS) was derived from the ECMWF-AUX product as the
difference in potential temperatures between 700 hPa and the surface
, and is used for a metric of macroscopic thermodynamic
conditions.
Results
Precipitation microphysics
A commonly known problem in GCMs associated with low cloud precipitation
microphysics is the timing of precipitation and its frequency and/or intensity
. This issue is also related to the magnitudes
of the aerosol indirect effect, i.e., dependency of precipitation on Nc.
As a proxy for this, we use the precipitation susceptibility (Sp)
metric, defined as
Sp=-dlnRdlnNc,
where R is the rain rate. Observed values of R are derived from CloudSat
radar reflectivity through the Z-R relationship ,
while model values of R are obtained as the large-scale precipitation rate.
We apply a threshold of R > 0.14 mm day-1, which is
equivalent to a radar reflectivity of -15 dBZe
, for precipitation in the model to enable a fair comparison
with satellite observations. The Sp metric is useful for examining the
aerosol impact on precipitation .
Figure a shows the behavior of Sp as a function of the LWP
obtained from MIROC and A-Train satellite observations. The satellite Sp
increases with increasing LWP up to around LWP ∼ 450 g m-2; further increases in LWP result in a decrease in
Sp. This behavior can be interpreted as follows .
Sp is low for a low LWP because clouds cannot generate much rainwater,
regardless of the aerosol loading. At a high LWP, Sp is also low because
precipitation is dominant, regardless of the aerosol loading, due to the
abundant LWP. In other words, the LWP value at which the Sp peaks
corresponds to the turning point where the water conversion process shifts
from the autoconversion regime to the accretion regime, as suggested by
previous studies .
On the other hand, the Sp amplitude predicted by MIROC is smaller than
the satellite results for a wide range of LWPs, and Sp remains high even
after its values peak near LWP ∼ 450 g m-2. This is
mainly because the autoconversion parameterization assumes a constant
dependency on Nc (i.e., β) regardless of the LWP, as is clear from
Eq. (). This also leads to a significant overestimation of
Sp for LWP < 100 g m-2, which means that the model
readily generates rainwater even when only a small amount of cloud water is
present.
To understand the uncertainty in the conversion process from cloud water to
rainwater in more detail, we define a new metric, the “susceptibility of
microphysical conversion (Sconv)” as
Sconv=-dlnPconvdlnNc.
This metric represents how aerosol burdens suppress rainwater production. In
satellite analysis, Pconv can be estimated from the method proposed by
, as shown in Eq. (). In the model,
Pconv is obtained as a native output of the process rate. The method of
was compared with a native model output of the process
rate in a global cloud-resolving model (CRM) by . The study
showed that the radar reflectivity is a gross measure of the water conversion
timescale, supporting the underlying assumption of .
This implies that the Sconv, which represents the timescale dependency
on Nc, can be compared between satellite observations and model
simulations although absolute values of Pconv can be different between
them. Although the use of satellite simulators would be helpful for more
direct comparison between model and observations, it is left as the subject
of future work.
Susceptibilities of (a) precipitation Sp and (b) microphysical
conversion Sconv as a function of the liquid water path (LWP) for the
MIROC–SPRINTARS results and A-Train observations. The left axis shows the
value of the susceptibility (refer to the line graph), and the right (red) axis
shows the probability distribution function for each LWP bin (refer to the
bar chart).
As shown in Fig. b, MIROC overestimates Sconv,
particularly in the lower LWP range. This means that the model generates
precipitation at a higher frequency, even at low LWPs, compared to
observations, which is mainly because the autoconversion in the model is too
rapid , as described above. Consequently, the
PDF of the LWP is biased toward lower
values because cloud water is depleted quickly by the rapid surface
precipitation. Alternatively, it is also possible that the model has biases
in the condensation processes, which lead to lower LWP and, thus, result in a
lower autoconversion rate. These tendencies in the model are strongly related
to unrealistically light rain that is too frequent, which is a common problem
in GCMs , including MIROC.
Besides this, Sconv can also be biased from the error of cloud geometric
thickness due to insufficient vertical resolution in GCMs. In addition to the
microphysical aspects mentioned above, biases in macrophysical structure are
also related to model performances, which will be discussed later (see Sect. 3.3).
Global distribution of dlnLWP/dlnNc from
(a, c) MIROC–SPRINTARS and (b, d) A-Train satellite estimations for
non-precipitating and precipitating clouds, respectively. The threshold of
the large-scale precipitation rate of 0.14 mm day-1 is used to
distinguish between non-precipitating or precipitating events (see text for
details).
Cloud susceptibilities
The response of cloud liquid water to aerosol perturbations determines the
cloud lifetime via the modification of cloud fraction ,
and is thus related to global hydrological cycles as well as radiation budget
e.g.,. As such, it is of great importance
to global climate studies to understand why there are two competing
mechanisms reported in the literature regarding the pathways of LWP, which
cause the LWP to either increase or decrease in response to an increase in
aerosols.
Figure shows the geographical distributions of
LWP susceptibility (i.e., the second term on the right side in Eq. ) obtained from the MIROC simulation and satellite retrievals.
The model produces positive values of dlnLWP/dlnNc
almost everywhere across the globe in both non-precipitating and precipitating
cases, which indicates that the LWP systematically increases with an
increasing aerosol burden. Even when the model applies AOD or hygroscopic
Na burden as an aerosol proxy instead of Nc, we obtain the similar
results (i.e., globally enhanced LWP). This result is expected from the model
parameterization of the cloud lifetime effect (Eq. ), which
monotonically delays the onset of precipitation in polluted conditions. This
is also a characteristic common to other GCMs, as reported in a recent study
. In contrast, the satellite-derived LWP susceptibility has a
coherent geographical pattern that includes both increasing and decreasing
responses, which is quite different from the model results. The decreasing
response occurs over the tropics and subtropics where more convective cloud
is dominant. The increasing response is apparent mainly over the midlatitudes
and regions where low clouds are dominant .
Global mean susceptibility (60∘ S–60∘ N) of τc,
re, and the LWP to Nc. The MIROC result is shown in orange and the
A-Train observation is in blue.
Nevertheless, it should be noted that Fig. c captures the
horizontal distribution of LWP susceptibility, whose pattern is very similar
to observations. That is, the relationship becomes weaker towards the
tropics, although the sign is still different. One of the possible mechanisms
is the dominance of cloud dynamical processes with high natural variability
over tropical and/or subtropical oceans rather than microphysical modifications by
aerosols . The same processes observed from
satellites could be at work in the model, and hence it might be related to
the parameterization of subgrid-scale variability. However, this is not
always true, particularly in non-precipitating cases (Fig. a),
so we must interpret the mechanisms carefully with further analysis in
future.
Figure b and d show that the geographical
patterns are qualitatively similar between the non-precipitating and
precipitating conditions, whereas the contrast between the two is slightly
different. More specifically, the value of dlnLWP/dlnNc is smaller in the precipitating condition than in the
non-precipitating case, which implies a smaller effect of aerosols when
precipitation occurs. However, it is noteworthy that the positive response of
dlnLWP/dlnNc in the non-precipitating condition has
a negative value in the precipitating conditions over East Asia, the eastern
United States, and Europe, where the anthropogenic aerosol burden is severe.
This suggests that aerosols act to prolong the cloud lifetime in
non-precipitating conditions, while they enhance cloud evaporation or can be
a precipitation driver in precipitating conditions which ultimately result in
less cloud water. These two competing mechanisms are reasonably consistent
with theories suggested by recent studies
, which propose the existence of a
buffering effect in the cloud system that results in smaller-magnitude
aerosol–cloud interactions. These comparisons suggest that the model does
not appropriately represent the buffering effect which compensates for the
positive responses of the LWP to aerosol perturbations. Current GCMs which
use similar parameterization frameworks (e.g., autoconversion) therefore
inherently overestimate the aerosol indirect effect as reported by previous
studies .
Figure summarizes the relationship among the
“process-oriented metrics” corresponding to each term in Eq. (). The global mean susceptibilities (averaged from
60∘ S to 60∘ N) for each term were calculated from the individual
susceptibility in each grid box in which more than 10 warm cloud samples were
obtained, which contributed to the reduction of statistical noise.
The cloud susceptibility of τc to Nc in MIROC is approximately
twice as large as in the A-Train results. The significant bias is decomposed
into contributions due to the “albedo effect” -dlnre/dlnNc and the “lifetime effect” dlnLWP/dlnNc.
Although the former is underestimated in the model compared with A-Train
estimations, its sign is positive and consistent with observations. The
overestimation of τc-susceptibility in the model is therefore
attributed to a positive response of the LWP, which is in stark contrast to
the slight negative responses in satellite observations.
reported a wide diversity in the relationship between the LWP and Nc
among nine AeroCom GCMs, with all models showing an enhanced response of LWP
to increased Nc. This further causes uncertainties in the estimation of
radiative forcing .
As Figs. and indicate, the discrepancy in the
LWP response between the model and observations can be a critical source of
model uncertainty, which causes a bias in climate responses via the
aerosol–cloud–precipitation–climate interaction.
Dependency of LWP responses on meteorology
Another key question regarding the LWP response for aerosol perturbations is
how and to what extent it depends on macrophysics such as cloud regimes and
thermodynamic conditions in the real atmosphere. Recent studies have
suggested that the source of uncertainty in the LWP response could be
attributed to differences in meteorology, different cloud types and regimes,
or more theoretical reasons, based on satellite observations
, large-eddy simulation (LES;
), and GCM intercomparison
.
To address this question, we examine the dependency of the LWP susceptibility
on both column maximum radar reflectivity (Zmax) and LTS as shown in
Fig. . Given that the horizontal axis characterizes the rain
regime (i.e., non-precipitating, drizzling, or precipitating) and the
vertical axis represents the thermodynamical stability conditions (i.e.,
unstable, intermediate, or stable), the diagram illustrates how the LWP
response to perturbed aerosols varies as a function of both rain
characteristics and stability conditions, thus providing a way to classify
cloud susceptibility according to macroscopic and meteorological conditions.
Figure clearly shows a systematic variation of the cloud
susceptibility under the two conditions. Positive responses of the LWP to
Nc are dominant in the non-precipitating and stable environments, while
negative responses can be seen in precipitating and unstable conditions. The
top-left region in the diagram corresponds to a stratocumulus regime in the
marine boundary layer. Because this type of cloud typically produces light
precipitation (i.e., drizzle; ) rather than heavy
precipitation it depletes a large amount of cloud water, and the aerosols
ingested into this type of cloud effectively act to enhance cloud water
storage, resulting in a positive response of LWP to an increased aerosol
loading. In contrast, the right-bottom region in the diagram corresponds to a
more convective cumulus regime, which is present mainly over the tropics.
This type of cloud is characterized by a relatively fast precipitation
timescale , and is favorable for cloud water
evaporation due to the larger extent of entrainment mixing ,
which results in negative responses of the LWP.
Susceptibility matrix of the LWP response to Nc as a function
of column maximum radar reflectivity (Zmax) and lower-tropospheric
stability (LTS) based on A-Train satellite data.
It is interesting that there is a positive correlation in the top-right
region even though precipitation occurs. One possible interpretation of this
is that the water-vapor supply is dominant over the loss of cloud water by
precipitation, and this type of cloud may correspond to the sustained frontal
precipitation (precipitating nimbostratus) systems found mainly over
midlatitude oceanic regions, where water vapor is abundant. In pristine and clean
environments, referred to as “aerosol-limited” conditions
, aerosols ingested into clouds will tend to store the cloud
water but also simultaneously produce more rain due to abundant water
mass. We note that it is just a speculation at this stage, and it might be
related to background aerosol number and environmental conditions (cf. Sect. 4 for more discussion). It is also noteworthy that the bottom-left region
displays negative susceptibilities even though precipitation does not occur.
Non-precipitating clouds in a significantly unstable environment would
correspond to inland and/or daytime cumulus. This tendency agrees with the results
of a previous study that focused on the non-precipitating
cumulus regime, and suggests a mechanism whereby the LWP decreases with
increased aerosol loading via evaporation–entrainment feedback. This results
in a loss of cloud water without precipitation.
Although the model version of the LWP-susceptibility diagram is not shown, it
will indicate a positive value in the matrix overall, as is obvious from Fig. . The mechanisms proposed above must be confirmed by more
detailed examinations using GCM and CRM with satellite simulators, or using
high-resolution process modeling, such as LES, in future studies. However,
the observation-based findings described above strongly suggest that rigorous
studies focusing on macrophysical conditions, including regional
characteristics of meteorological factors, in addition to microphysical
conditions, are indispensable for better understanding the response of the
aerosol–cloud–precipitation interaction.
Summary and discussion
We explored the source of discrepancy in the aerosol–cloud–precipitation
interaction for warm clouds between an aerosol–climate model and A-Train
satellite retrieval. The instantaneous model output was analyzed using as
many samples as possible to provide reliable statistics and fair comparisons
with satellite observations.
We found critical biases in the model in the response of the LWP to aerosol
perturbations. The model predicted a monotonic increase in the LWP across the
globe, in contrast to the observations that clearly showed a regional
variation of the LWP response that either increased or decreased with an
increasing aerosol loading. This variability in cloud susceptibility observed
by the A-Train was closely related to differences in meteorological factors,
such as cloud regimes and thermodynamic conditions. For example, stratiform
clouds under stable conditions had a tendency to increase the LWP given
aerosol perturbations, while cumulus clouds over an unstable environment
tended to decrease the LWP as the aerosol loading increased. The
bidirectional responses of LWP (both positive and negative) found in
satellite observations in different aerosol concentrations might be related
to the concept of “optimal aerosol concentration (Nop)”, recently
suggested by . More specifically, in case of Na<Nop, clouds tend to be deeper with larger liquid mass, referred
to as cloud invigoration e.g.,, for increased aerosol
loading, whereas the case of Na>Nop would be favorable for
cloud suppression due to enhanced entrainment and evaporation. This could
lead the bidirectional LWP susceptibilities, although we cannot mention the
exact mechanisms at this stage because Nop also depends on both cloud
geometric scale and environmental conditions
.
This can explain why previous studies have reported conflicting results for
the LWP response, with either an increase or decrease with increasing aerosol
loading e.g.,. Previous
studies have focused on different study regions and/or targets, which has
resulted in different cloud responses due to the different mechanisms of
aerosol–cloud–precipitation interaction. This means that global-mean cloud
susceptibility is not very meaningful in constraining the
aerosol–cloud–precipitation interaction. Future studies should consider the
regional dependence of the susceptibility metrics .
The monotonic increase in the LWP with increasing aerosol loading in the
model is attributed to the autoconversion scheme, which assumes only
suppression of rainwater generation to account for the traditional cloud
lifetime effect without its compensation, and does not take meteorological
conditions into account. Mechanisms that can decrease the LWP in polluted
conditions, such as the enhancement of evaporation due to entrainment mixing
, are not incorporated in our model. This means
that the model fails to represent the buffering effect for the
aerosol–cloud–precipitation interaction . Moreover, the
model overestimates Sconv around low LWPs compared with A-Train
satellite retrievals due to uncertainties in process rate parameterization
. This is evidence that the autoconversion in the model is
too fast, which results in the LWP having a high dependency on the Nc.
This bias in the model is consistent with a previous study that reported a
higher LWP susceptibility in GCMs due to their diagnostic treatment of
rainwater .
In future studies, the aerosol–cloud–precipitation framework must be
expanded to represent the effect of environmental conditions in a flexible
manner, in addition to the microphysics. Current microphysical frameworks
without such macrophysical aspects bring highly sensitive aerosol–cloud
interactions, although they vary to some extent depending on the
autoconversion scheme .
also reported a wide diversity in the LWP response to
Nc among various GCMs, and concluded that their inconsistency could
mainly be attributed to their different representations of the autoconversion
process. However, it is also true that different choices of microphysical
scheme alone do not significantly improve the model biases in both cloud
physics and cloud radiative effects , and these two
requirements sometimes contradict each other . This is due
to the arbitrary nature of tuning and assumptions (e.g., artificial threshold
parameters, diagnostic treatment of rain), which is a bottleneck in GCMs
. Recently, a fundamental model
improvement was achieved by introducing a prognostic precipitation framework
, which represents important progress in
process representations for more realistic cloud and precipitation
microphysics. This improvement is expected to overcome some of the common
problems in GCMs, such as the overestimation of the aerosol indirect effect
and spurious light rain .
Furthermore, a representation of subgrid-scale fluctuations has also been
critical problem in GCMs. Although the magnitude as well as sign of
LWP susceptibility differs between the model and observations, the horizontal
pattern is similar in precipitating conditions. The parameterization of
subgrid-scale variability may partly contribute to the weakening of the aerosol roles
by capturing the large natural variability of clouds, especially over
tropical and/or subtropical oceans , which would lead
to a more realistic representation of cloud dynamical processes. For example,
showed that both positive and negative LWP responses
can be represented in even a GCM framework, by the PDF-based macrophysics
parameterization, called “Cloud Layers Unified By Binormals
CLUBB;”. estimated a weighting
factor of process rate equations to consider the subgrid effects based on
A-Train retrievals unless the accretion process is significantly underestimated.
The interaction between microphysics and subgrid-scale dynamics
(microphysics–dynamics interactions) in GCMs is therefore one of the
indispensable processes for incorporating buffering effects and for improving
model physics as a whole.
Although this study focused only on warm-phase clouds, our findings regarding
different cloud responses to aerosol perturbations between GCMs and satellite
observations will assist future model development for more accurate climate
simulations. Further studies should also contain an extension of the research
target from liquid to mixed and/or iced clouds, and from a process-level to a cloud
system to understand the whole cloud system response to aerosol
perturbations, taking into account the buffered system morphology.