Stratocumulus Regions Research Articles

Investigating the sign of stratocumulus adjustments to aerosols in the ICON global storm-resolving model

Abstract. Aerosols can cause brightening of stratocumulus clouds, thereby cooling the climate. Observations and models disagree on the magnitude of this cooling, partly because of the aerosol-induced liquid water path (LWP) adjustment, with climate models predicting an increase in the LWP and satellites observing a weak decrease in response to increasing aerosols. With higher-resolution global climate models, which allow the simulation of mesoscale circulations in which stratocumulus clouds are embedded, there is hope to start bridging this gap. In this study, we present boreal summertime simulations conducted with the ICOsahedral Non-hydrostatic (ICON) global storm-resolving model (GSRM). Compared to geostationary satellite data, ICON produces realistic cloud coverage in the stratocumulus regions; however, these clouds look cumuliform, and the sign of LWP adjustments disagrees with observations. We investigate this disagreement with a causal approach, which combines time series with knowledge of cloud processes, allowing us to diagnose the sources of observation–model discrepancies. The positive ICON LWP adjustment results from a superposition of processes, with an overestimated positive response due to (1) precipitation suppression, (2) a lack of wet scavenging, and (3) cloud deepening under a weak inversion, despite (4) small negative influences from cloud-top entrainment enhancement. We also find that precipitation suppression and entrainment enhancement occur at different intensities during the day and the night, implying that daytime satellite studies suffer from selection bias. This causal methodology can guide modelers on how to modify model parameterizations and setups to reconcile conflicting studies concerning the sign and magnitude of LWP adjustments across different data sources.

Satellite-based analysis of top of atmosphere shortwave radiative forcing trend induced by biomass burning aerosols over South-Eastern Atlantic

This study investigates long-term changes in the shortwave direct aerosol radiative effect (DARE) at the top of the atmosphere (TOA) induced by biomass burning aerosol (BBA) transported from southern Africa to the south-eastern Atlantic (SEA) stratocumulus region during extended fire seasons. The evolution since 2002 of aerosol, cloud properties, and TOA shortwave outgoing radiation from advanced passive satellite sensors are presented, as well as the observational trend in clear-sky DAREclr and the retrieval trend in all-sky DAREall. Supplemented by chemical transport model simulations, we estimate that DAREclr has become more negative (−0.09 ± 0.06 W m−2 yr−1) due to increased aerosol presence in SEA. Meanwhile, DAREall has become more positive ( + 0.04 ± 0.15 W m−2 yr−1) due to aerosols in cloudy sky regions. This study reveals satellite capabilities in capturing complex BBA-cloud-solar radiation interactions for accurate radiative forcing estimates and projections.

How Cloud Droplet Number Concentration Impacts Liquid Water Path and Precipitation in Marine Stratocumulus Clouds—A Satellite-Based Analysis Using Explainable Machine Learning

Aerosol–cloud–precipitation interactions (ACI) are a known major cause of uncertainties in simulations of the future climate. An improved understanding of the in-cloud processes accompanying ACI could help in advancing their implementation in global climate models. This is especially the case for marine stratocumulus clouds, which constitute the most common cloud type globally. In this work, a dataset composed of satellite observations and reanalysis data is used in explainable machine learning models to analyze the relationship between the cloud droplet number concentration (Nd), cloud liquid water path (LWP), and the fraction of precipitating clouds (PF) in five distinct marine stratocumulus regions. This framework makes use of Shapley additive explanation (SHAP) values, allowing to isolate the impact of Nd from other confounding factors, which proved to be very difficult in previous satellite-based studies. All regions display a decrease of PF and an increase in LWP with increasing Nd, despite marked inter-regional differences in the distribution of Nd. Polluted (high Nd) conditions are characterized by an increase of 12 gm−2 in LWP and a decrease of 0.13 in PF on average when compared to pristine (low Nd) conditions. The negative Nd–PF relationship is stronger in high LWP conditions, while the positive Nd–LWP relationship is amplified in precipitating clouds. These findings indicate that precipitation suppression plays an important role in MSC adjusting to aerosol-driven perturbations in Nd.

Cloud properties and their projected changes in CMIP models with low to high climate sensitivity

Abstract. Since the release of the first Coupled Model Intercomparison Project version 6 (CMIP6) simulations, one of the most discussed topics is the higher effective climate sensitivity (ECS) of some of the models, resulting in an increased range of ECS values in CMIP6 compared to previous CMIP phases. An important contribution to ECS is the cloud climate feedback. Although climate models have continuously been developed and improved over the last few decades, a realistic representation of clouds remains challenging. Clouds contribute to the large uncertainties in modeled ECS, as projected changes in cloud properties and cloud feedbacks also depend on the simulated present-day fields. In this study, we investigate the representation of both cloud physical and radiative properties from a total of 51 CMIP5 and CMIP6 models. ECS is used as a simple metric to group the models, as the sensitivity of the physical cloud properties to warming is closely related to cloud feedbacks, which in turn are known to have a large contribution to ECS. Projected changes in the cloud properties in future scenario simulations are analyzed by the ECS group. In order to help with interpreting the projected changes, model results from historical simulations are also analyzed. The results show that differences in the net cloud radiative effect as a reaction to warming among the three model groups are driven by changes in a range of cloud regimes rather than individual regions. In polar regions, high-ECS models show a weaker increase in the net cooling effect of clouds, due to warming, than the low-ECS models. At the same time, high-ECS models show a decrease in the net cooling effect of clouds over the tropical ocean and the subtropical stratocumulus regions, whereas low-ECS models show either little change or even an increase in the cooling effect. Over the Southern Ocean, the low-ECS models show a higher sensitivity of the net cloud radiative effect to warming than the high-ECS models.

Global observations of aerosol indirect effects from marine liquid clouds

Abstract. Interactions between aerosols and liquid clouds are one of the largest sources of uncertainty in the historical radiative forcing of climate. One widely shared goal to reduce this uncertainty is to decompose radiative anomalies arising from aerosol–cloud interactions into components associated with changes in cloud-droplet number concentration (Twomey effect), liquid-water-path adjustments, and cloud-fraction adjustments. However, there has not been a quantitative foundation for simultaneously estimating these components with global satellite observations. Here we present a method for assessing shortwave radiative flux anomalies from the Twomey effect and cloud adjustments over ocean between 55∘ S and 55∘ N. We find that larger aerosol concentrations are associated with widespread cloud brightening from the Twomey effect, a positive radiative adjustment from decreasing liquid water path in subtropical stratocumulus regions, and a negative radiative adjustment from increasing cloud fraction in the subtropics and midlatitudes. The Twomey effect and total cloud adjustment have contributed −0.77 ± 0.25 and −1.02 ± 0.43 W m−2, respectively, to the effective radiative forcing since 1850 over the domain (95 % confidence). Our findings reduce uncertainty in these components of aerosol forcing and suggest that cloud adjustments make a larger contribution to the forcing than is commonly believed.

Dominant Cloud Controlling Factors for Low‐Level Cloud Fraction: Subtropical Versus Extratropical Oceans

AbstractTo improve cloud feedback understanding and simulation, observations have been used to quantify the rate of change of cloud radiative properties as a function of specific environmental metrics (or cloud controlling factors; CCFs). The study focuses on low‐level cloud dominated regions during 2006–2010. For each ocean gridpoint, Spearman's rank correlation coefficients of daily mean observed cloud fraction versus (a) 10‐m wind, (b) sensible heat flux (SHF), (c) sea surface temperature, (d) estimated inversion strength (EIS), (e) 850 hPa vertical velocity, and (f) the M parameter () are sorted to identify the dominant CCFs in both extratropics (30°–60°N/S) and subtropics (30°S–30°N). A novel map for visualizing dominant CCFs for low‐level cloud fraction reveals that: dominates in the subtropical stratocumulus regions while 10‐m winds dominate in shallow cumulus regions but in the extratropics, a different inversion structure diagnostic (EIS) dominates, while SHF dominates in western boundary current areas.

Sensitivities of cloud radiative effects to large-scale meteorology and aerosols from global observations

Abstract. The radiative effects of clouds make a large contribution to the Earth's energy balance, and changes in clouds constitute the dominant source of uncertainty in the global warming response to carbon dioxide forcing. To characterize and constrain this uncertainty, cloud-controlling factor (CCF) analyses have been suggested that estimate sensitivities of clouds to large-scale environmental changes, typically in cloud-regime-specific multiple linear regression frameworks. Here, local sensitivities of cloud radiative effects to a large number of controlling factors are estimated in a regime-independent framework from 20 years (2001–2020) of near-global (60∘ N–60∘ S) satellite observations and reanalysis data using statistical learning. A regularized linear regression (ridge regression) is shown to skillfully predict anomalies of shortwave (R2=0.63) and longwave cloud radiative effects (CREs) (R2=0.72) in independent test data on the basis of 28 CCFs, including aerosol proxies. The sensitivity of CREs to selected CCFs is quantified and analyzed. CRE sensitivities to sea surface temperature and estimated inversion strength are particularly pronounced in low-cloud regions and generally in agreement with previous studies. The analysis of CRE sensitivities to three-dimensional wind field anomalies reflects the fact that CREs in tropical ascent regions are mainly driven by variability of large-scale vertical velocity in the upper troposphere. In the subtropics, CRE is sensitive to free-tropospheric zonal and meridional wind anomalies, which are likely to encapsulate information on synoptic variability that influences subtropical cloud systems by modifying wind shear and thus turbulence and dry-air entrainment in stratocumulus clouds, as well as variability related to midlatitude cyclones. Different proxies for aerosols are analyzed as CCFs, with satellite-derived aerosol proxies showing a larger CRE sensitivity than a proxy from an aerosol reanalysis, likely pointing to satellite aerosol retrieval biases close to clouds, leading to overestimated aerosol sensitivities. Sensitivities of shortwave CREs to all aerosol proxies indicate a pronounced cooling effect from aerosols in stratocumulus regions that is counteracted to a varying degree by a longwave warming effect. The analysis may guide the selection of CCFs in future sensitivity analyses aimed at constraining cloud feedback and climate forcings from aerosol–cloud interactions using data from both observations and global climate models.

Distinct regional meteorological influences on low-cloud albedo susceptibility over global marine stratocumulus regions

Abstract. Marine stratocumuli cool the Earth effectively due to their high reflectance of incoming solar radiation and persistent occurrence. The susceptibility of cloud albedo to droplet number concentration perturbations depends strongly on large-scale meteorological conditions. Studies focused on the meteorological dependence of cloud adjustments often overlook the covariability among meteorological factors and their geographical and temporal variability. We use 8 years of satellite observations sorted by day and geographical location to show the global distribution of marine low-cloud albedo susceptibility. We find an overall cloud brightening potential for most marine warm clouds, which is more pronounced over subtropical coastal regions. A weak cloud darkening potential in the annual mean is evident over the remote SE Pacific and SE Atlantic. We show that large-scale meteorological fields from the ERA5 reanalysis data, including lower-tropospheric stability, free-tropospheric relative humidity, sea surface temperature, and boundary layer depth, have distinct covariabilities over each of the eastern subtropical ocean basins where marine stratocumuli prevail. This leads to a markedly different annual cycle in albedo susceptibility over each basin. Moreover, we find that basin-specific regional relationships between key meteorological factors and albedo susceptibilities are absent in a global analysis. Our results stress the importance of considering the geographical distinctiveness of temporal meteorological covariability when scaling up the local-to-global response of cloud albedo to aerosol perturbations.

Cloud Climatologies from Global Climate Models—A Comparison of CMIP5 and CMIP6 Models with Satellite Data

Abstract Simulating clouds with global climate models is challenging as the relevant physics involves many nonlinear processes covering a wide range of spatial and temporal scales. As key components of the hydrological cycle and the climate system, an evaluation of clouds from models used for climate projections is an important prerequisite for assessing the confidence in the results from these models. Here, we compare output from models contributing to phase 6 of the Coupled Model Intercomparison Project (CMIP6) with satellite data and with results from their predecessors (CMIP5). We use multiproduct reference datasets to estimate the observational uncertainties associated with different sensors and with internal variability on a per-pixel basis. Selected cloud properties are also analyzed by region and by dynamical regime and thermodynamic conditions. Our results show that for parameters such as total cloud cover, cloud water path, and cloud radiative effect, the CMIP6 multimodel mean performs slightly better than the CMIP5 ensemble mean in terms of mean bias, pattern correlation, and relative root-mean square deviation. The intermodel spread in CMIP6, however, is not reduced compared to CMIP5. Compared with CALIPSO-ICECLOUD data, the CMIP5/6 models overestimate cloud ice, particularly in the lower and middle troposphere, partly due to too high ice fractions for given temperatures. This bias is reduced in the CMIP6 multimodel mean. While many known biases such as an underestimation in cloud cover in stratocumulus regions remain in CMIP6, we find that the CMIP5 problem of too few but too reflective clouds over the Southern Ocean is significantly improved.

Linking Large-Scale Double-ITCZ Bias to Local-Scale Drizzling Bias in Climate Models

Abstract Tropical precipitation in climate models presents significant biases in both the large-scale pattern (i.e., double intertropical convergence zone bias) and local-scale characteristics (i.e., drizzling bias with too frequent drizzle/convection and reduced occurrences of no and heavy precipitation). By untangling the coupled system and analyzing the biases in precipitation, cloud, and radiation, this study shows that local-scale drizzling bias in atmospheric models can lead to large-scale double-ITCZ bias in coupled models by inducing convective-regime-dependent biases in precipitation and cloud radiative effects (CRE). The double-ITCZ bias consists of a hemispherically asymmetric component that arises from the asymmetric SST bias and a nearly symmetric component that exists in atmospheric models without the SST bias. By increasing light rain but reducing heavy rain, local-scale drizzling bias induces positive (negative) precipitation bias in the moderate (strong) convective regime, leading to the nearly symmetric wet bias in atmospheric models. By affecting the cloud profile, local-scale drizzling bias induces positive (negative) CRE bias in the stratocumulus (convective) regime in atmospheric models. Because the stratocumulus (convective) region is climatologically more pronounced in the southern (northern) tropics, the CRE bias is deemed to be hemispherically asymmetric and drives warm and wet (cold and dry) biases in the southern (northern) tropics when coupled to ocean. Our results suggest that correcting local-scale drizzling bias is critical for fixing large-scale double-ITCZ bias. The drizzling and double-ITCZ biases are not alleviated in models with mesoscale (0.25°–0.5°) or even storm-resolving (∼3 km) resolution, implying that either large-eddy simulation or fundamental improvement in small-scale subgrid parameterizations is needed.

The impact of sampling strategy on the cloud droplet number concentration estimated from satellite data

Abstract. Cloud droplet number concentration (Nd) is of central importance to observation-based estimates of aerosol indirect effects, being used to quantify both the cloud sensitivity to aerosol and the base state of the cloud. However, the derivation of Nd from satellite data depends on a number of assumptions about the cloud and the accuracy of the retrievals of the cloud properties from which it is derived, making it prone to systematic biases. A number of sampling strategies have been proposed to address these biases by selecting the most accurate Nd retrievals in the satellite data. This work compares the impact of these strategies on the accuracy of the satellite retrieved Nd, using a selection of in situ measurements. In stratocumulus regions, the MODIS Nd retrieval is able to achieve a high precision (r2 of 0.5–0.8). This is lower in other cloud regimes but can be increased by appropriate sampling choices. Although the Nd sampling can have significant effects on the Nd climatology, it produces only a 20 % variation in the implied radiative forcing from aerosol–cloud interactions, with the choice of aerosol proxy driving the overall uncertainty. The results are summarised into recommendations for using MODIS Nd products and appropriate sampling.

Environmental Response in Coupled Energy and Water Cloud Impact Parameters Derived from A-Train Satellites, ERA-Interim, and MERRA-2

Abstract Understanding the connections between latent heating from precipitation and cloud radiative effects is essential for accurately parameterizing cross-scale links between cloud microphysics and global energy and water cycles in climate models. Although commonly examined separately, this study adopts two cloud impact parameters (CIPs), the surface radiative cooling efficiency Rc and atmospheric radiative heating efficiency Rh, that explicitly couple cloud radiative effects and precipitation to characterize how efficiently precipitating cloud systems influence the energy budget and water cycle using A-Train observations and two reanalyses. These CIPs exhibit distinct global distributions that suggest cloud energy and water cycle coupling are highly dependent on cloud regime. The dynamic regime ω500 controls the sign of Rh, whereas column water vapor (CWV) appears to be the larger control on the magnitude. The magnitude of Rc is highly coupled to the dynamic regime. Observations show that clouds cool the surface very efficiently per unit rainfall at both low and high sea surface temperature (SST) and CWV, but reanalyses only capture the former. Reanalyses fail to simulate strong Rh and moderate Rc in deep convection environments but produce stronger Rc and Rh than observations in shallow, warm rain systems in marine stratocumulus regions. Although reanalyses generate fairly similar climatologies in the frequency of environmental states, the response of Rc and Rh to SST and CWV results in systematic differences in zonal and meridional gradients of cloud atmospheric heating and surface cooling relative to A-Train observations that may have significant implications for global circulations and cloud feedbacks. Significance Statement Studying climate change requires understanding coupled interactions between clouds, precipitation, and their environment. Here we calculate two parameters to reveal how efficiently clouds can heat the atmosphere or cool the surface per unit rain. The satellite observations and reanalyses show similar global patterns, but there are some differences in areas of deep convection and low cloud regions. Examination of these parameters as a function of their environment shows that reanalyses cool the atmosphere too much per unit rain in environments with low sea surface temperatures and water vapor. Vertical velocity determines whether clouds heat or cool the atmosphere. Both observations and reanalyses suggest that water vapor is the stronger control on how much clouds heat the atmosphere per unit rain.

Albedo susceptibility of northeastern Pacific stratocumulus: the role of covarying meteorological conditions

Abstract. Quantification of the radiative adjustment of marine low clouds to aerosol perturbations, regionally and globally, remains the largest source of uncertainty in assessing current and future climate. One of the important steps towards quantifying the role of aerosol in modifying cloud radiative properties is to quantify the susceptibility of cloud albedo and liquid water path (LWP) to perturbations in cloud droplet number concentration (Nd). We use 10 years of spaceborne observations from the polar-orbiting Aqua satellite to quantify the albedo susceptibility of marine low clouds to Nd perturbations over the northeast (NE) Pacific stratocumulus (Sc) region. Mutual information analysis reveals a dominating control of cloud state (e.g., LWP and Nd) on low-cloud albedo susceptibility, relative to the meteorological states that drive these cloud states. Through a LWP–Nd space decomposition of albedo susceptibilities, we show clear separation among susceptibility regimes (brightening or darkening), consistent with previously established mechanisms through which aerosol modulates cloud properties. These regimes include (i) thin non-precipitating clouds (LWP < 55 g m−2) that exhibit brightening (occurring 37 % of the time), corresponding to the Twomey effect; (ii) thicker non-precipitating clouds, corresponding to entrainment-driven negative LWP adjustments that manifest as a darkening regime (36 % of the time); and (iii) another brightening regime (22 % of the time) consisting of mostly precipitating clouds, corresponding to precipitation-suppression LWP positive adjustments. Overall, we find an annual-mean regional low-cloud brightening potential of 20.8±2.68 W m−2 ln(Nd)−1, despite an overall negative LWP adjustment for non-precipitating marine stratocumulus, owing to the high occurrence of the Twomey–brightening regime. Over the NE Pacific, clear seasonal covariabilities among meteorological factors related to the large-scale circulation are found to play an important role in grouping conditions favorable for each susceptibility regime. When considering the covarying meteorological conditions, our results indicate that for the northeastern Pacific stratocumulus, clouds that exhibit the strongest brightening potential occur most frequently within shallow marine boundary layers over a cool ocean surface with a stable atmosphere and a dry free troposphere above. Clouds that exhibit a darkening potential associated with negative LWP adjustments occur most frequently within deep marine boundary layers in which the atmospheric instability and the ocean surface are not strong and warm enough to produce frequent precipitation. Cloud brightening associated with warm-rain suppression is found to preferably occur either under unstable atmospheric conditions or humid free-tropospheric conditions that co-occur with a warm ocean surface.

Modeled and observed properties related to the direct aerosol radiative effect of biomass burning aerosol over the southeastern Atlantic

Abstract. Biomass burning smoke is advected over the southeastern Atlantic Ocean between July and October of each year. This smoke plume overlies and mixes into a region of persistent low marine clouds. Model calculations of climate forcing by this plume vary significantly in both magnitude and sign. NASA EVS-2 (Earth Venture Suborbital-2) ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) had deployments for field campaigns off the west coast of Africa in 3 consecutive years (September 2016, August 2017, and October 2018) with the goal of better characterizing this plume as a function of the monthly evolution by measuring the parameters necessary to calculate the direct aerosol radiative effect. Here, this dataset and satellite retrievals of cloud properties are used to test the representation of the smoke plume and the underlying cloud layer in two regional models (WRF-CAM5 and CNRM-ALADIN) and two global models (GEOS and UM-UKCA). The focus is on the comparisons of those aerosol and cloud properties that are the primary determinants of the direct aerosol radiative effect and on the vertical distribution of the plume and its properties. The representativeness of the observations to monthly averages are tested for each field campaign, with the sampled mean aerosol light extinction generally found to be within 20 % of the monthly mean at plume altitudes. When compared to the observations, in all models, the simulated plume is too vertically diffuse and has smaller vertical gradients, and in two of the models (GEOS and UM-UKCA), the plume core is displaced lower than in the observations. Plume carbon monoxide, black carbon, and organic aerosol masses indicate underestimates in modeled plume concentrations, leading, in general, to underestimates in mid-visible aerosol extinction and optical depth. Biases in mid-visible single scatter albedo are both positive and negative across the models. Observed vertical gradients in single scatter albedo are not captured by the models, but the models do capture the coarse temporal evolution, correctly simulating higher values in October (2018) than in August (2017) and September (2016). Uncertainties in the measured absorption Ångstrom exponent were large but propagate into a negligible (<4 %) uncertainty in integrated solar absorption by the aerosol and, therefore, in the aerosol direct radiative effect. Model biases in cloud fraction, and, therefore, the scene albedo below the plume, vary significantly across the four models. The optical thickness of clouds is, on average, well simulated in the WRF-CAM5 and ALADIN models in the stratocumulus region and is underestimated in the GEOS model; UM-UKCA simulates cloud optical thickness that is significantly too high. Overall, the study demonstrates the utility of repeated, semi-random sampling across multiple years that can give insights into model biases and how these biases affect modeled climate forcing. The combined impact of these aerosol and cloud biases on the direct aerosol radiative effect (DARE) is estimated using a first-order approximation for a subset of five comparison grid boxes. A significant finding is that the observed grid box average aerosol and cloud properties yield a positive (warming) aerosol direct radiative effect for all five grid boxes, whereas DARE using the grid-box-averaged modeled properties ranges from much larger positive values to small, negative values. It is shown quantitatively how model biases can offset each other, so that model improvements that reduce biases in only one property (e.g., single scatter albedo but not cloud fraction) would lead to even greater biases in DARE. Across the models, biases in aerosol extinction and in cloud fraction and optical depth contribute the largest biases in DARE, with aerosol single scatter albedo also making a significant contribution.

Vertical grid refinement for stratocumulus clouds in the radiation scheme of the global climate model ECHAM6.3-HAM2.3-P3

Abstract. In this study, we implement a vertical grid refinement scheme in the radiation routine of the global aerosol–climate model ECHAM-HAM, aiming to improve the representation of stratocumulus clouds and address the underestimation of their cloud cover. The scheme is based on a reconstruction of the temperature inversion as a physical constraint for the cloud top. On the refined grid, the boundary layer and the free troposphere are separated and the cloud's layer is made thinner. The cloud cover is recalculated either by conserving the cloud volume (SC-VOLUME) or by using the Sundqvist cloud cover routine on the new grid representation (SC-SUND). In global climate simulations, we find that the SC-VOLUME approach is inadequate, as there is a mismatch, in most cases, between the layer of the inversion and the layer of the stratocumulus cloud, which prevents its application and is itself likely caused by an overly low vertical resolution. Additionally, we find that the occurrence frequency of stratocumulus clouds is underestimated in ECHAM-HAM, limiting a priori the potential benefits of a scheme like SC-VOLUME targeting only cloud amount when present. With the SC-SUND approach, the possibility for new clouds to be formed on the refined grid results in a large increase in mean total cloud cover in stratocumulus regions. In both cases, however, the changes exerted in the radiation routine are too weak to produce a significant improvement in the simulated stratocumulus cloud cover. We investigate and discuss the reasons behind this. The grid refinement scheme could be used more effectively for this purpose if implemented directly in the model's cloud microphysics and cloud cover routines, but other possible ways forward are also discussed.

Vertical structure of cloud radiative heating in the tropics: confronting the EC-Earth v3.3.1/3P model with satellite observations

Abstract. Understanding the coupling of clouds to large-scale circulation is one of the grand challenges for the global climate research community. In this context, realistically modelling the vertical structure of cloud radiative heating (CRH) and/or cooling in Earth system models is a key premise to understand this coupling. Here, we evaluate CRH in two versions of the European Community Earth System Model (EC-Earth) using retrievals derived from the combined radar and lidar data from the CloudSat and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellites. One model version is also used with two different horizontal resolutions. Our study evaluates large-scale intraseasonal variability in the vertical structure of CRH and cloud properties and investigates the changes in CRH during different phases of the El Niño–Southern Oscillation (ENSO), a process that dominates the interannual climate variability in the tropics. EC-Earth generally captures both the intraseasonal and meridional pattern of variability in CRH over the convectively active and stratocumulus regions and the CRH during the positive and negative phases of ENSO. However, two key differences between model simulations and satellite retrievals emerge. First, the magnitude of CRH, in the upper troposphere, over the convectively active zones is up to twice as large in the models compared to the satellite data. Further dissection of net CRH into its shortwave and longwave components reveals noticeable differences in their vertical structure. The shortwave component of the radiative heating is overestimated by all model versions in the lowermost troposphere and underestimated in the middle troposphere. These over- and underestimates of shortwave heating are partly compensated by an overestimate of longwave cooling in the lowermost troposphere and heating in the middle troposphere. The biases in CRH can be traced back to disagreement in cloud amount and cloud water content. There is no noticeable improvement of CRH by increasing the horizontal resolution in the model alone. Our findings highlight the importance of evaluating models with satellite observations that resolve the vertical structure of clouds and cloud properties.

Processes controlling the vertical aerosol distribution in marine stratocumulus regions – a sensitivity study using the climate model NorESM1-M

Abstract. The vertical distribution of aerosols plays an important role in determining the effective radiative forcing from aerosol–radiation and aerosol–cloud interactions. Here, a number of processes controlling the vertical distribution of aerosol in five subtropical marine stratocumulus regions in the climate model NorESM1-M are investigated, with a focus on the total aerosol extinction. A comparison with satellite lidar data (CALIOP, Cloud–Aerosol Lidar with Orthogonal Polarization) shows that the model underestimates aerosol extinction throughout the troposphere, especially elevated aerosol layers in the two regions where they are seen in observations. It is found that the shape of the vertical aerosol distribution is largely determined by the aerosol emission and removal processes in the model, primarily through the injection height, emitted particle size, and wet scavenging. In addition, the representation of vertical transport related to shallow convection and entrainment is found to be important, whereas alterations in aerosol optical properties and cloud microphysics parameterizations have smaller effects on the vertical aerosol extinction distribution. However, none of the alterations made are sufficient for reproducing the observed vertical distribution of aerosol extinction, neither in magnitude nor in shape. Interpolating the vertical levels of CALIOP to the corresponding model levels leads to better agreement in the boundary layer and highlights the importance of the vertical resolution.

Direct and semi-direct radiative forcing of biomass-burning aerosols over the southeast Atlantic (SEA) and its sensitivity to absorbing properties: a regional climate modeling study

Abstract. Simulations are performed for the period 2000–2015 by two different regional climate models, ALADIN and RegCM, to quantify the direct and semi-direct radiative effects of biomass-burning aerosols (BBAs) in the southeast Atlantic (SEA) region. Different simulations have been performed using strongly absorbing BBAs in accordance with recent in situ observations over the SEA. For the July–August–September (JAS) season, the single scattering albedo (SSA) and total aerosol optical depth (AOD) simulated by the ALADIN and RegCM models are consistent with the MACv2 climatology and MERRA-2 and CAMS-RA reanalyses near the biomass-burning emission sources. However, the above-cloud AOD is slightly underestimated compared to satellite (MODIS and POLDER) data during the transport over the SEA. The direct radiative effect exerted at the continental and oceanic surfaces by BBAs is significant in both models and the radiative effects at the top of the atmosphere indicate a remarkable regional contrast over SEA (in all-sky conditions), with a cooling (warming) north (south) of 10 ∘S, which is in agreement with the recent MACv2 climatology. In addition, the two models indicate that BBAs are responsible for an important shortwave radiative heating of ∼0.5–1 K per day over SEA during JAS with maxima between 2 and 4 km a.m.s.l. (above mean sea level). At these altitudes, BBAs increase air temperature by ∼0.2–0.5 K, with the highest values being co-located with low stratocumulus clouds. Vertical changes in air temperature limit the subsidence of air mass over SEA, creating a cyclonic anomaly. The opposite effect is simulated over the continent due to the increase in lower troposphere stability. The BBA semi-direct effect on the lower troposphere circulation is found to be consistent between the two models. Changes in the cloud fraction are moderate in response to the presence of smoke, and the models differ over the Gulf of Guinea. Finally, the results indicate an important sensitivity of the direct and semi-direct effects to the absorbing properties of BBAs. Over the stratocumulus (Sc) region, DRE varies from +0.94 W m−2 (scattering BBAs) to +3.93 W m−2 (most absorbing BBAs).

Diurnal Variation in Cloud Liquid Water Path Derived from Five Cross-Track Microwave Radiometers Onboard Polar-Orbiting Satellites

The Advanced Microwave Sounding Unit (AMSU)-A/Advanced Technology Microwave Sounder (ATMS) onboard the National Oceanic Atmospheric Administration (NOAA)-18/-19, MetOp-A/-B, and Suomi National Polar-orbiting Partnership satellites provide global observations of the cloud Liquid Water Path (LWP) almost 10 times a day. This study explores the possibility of capturing the diurnal cycle of the LWP. An inter-satellite cross-calibration is first carried out using a double-difference method. A remapping is then used to obtain the AMSU-A-like LWP to account for beam shape discrepancies between the ATMS and AMSU-A. We finally examine the diurnal cycle of the LWP over the Southeast Pacific Ocean using the ATMS and AMSU-A data from the five satellites mentioned above. Results show that the remapped ATMS results agree well with the AMSU-A results at the same local time over a stratocumulus region. LWP retrievals from multiple satellite cross-track microwave radiometers can well reproduce the diurnal variation characteristics of LWP in 2015 over the East Pacific Ocean, including the seasonal variation of the diurnal variation. This study presents the first step toward merging LWP data from all ATMS and AMSU-A radiometers and will be of interest to many researchers studying LWP-related weather and climate changes, especially considering the possible loss of higher-resolution microwave-frequency conical-scanning sensors in the coming years.

Investigating physical constraints on climate feedbacks using a perturbed parameter ensemble

A large parameter-perturbed ensemble (PPE) of the Met Office climate model is used to explore the relationship between radiative feedbacks and the present-day simulation of the associated physical processes. We highlight three tropical regimes (deep convection over ocean and land, and marine stratocumulus) in which the same set of processes drives the present-day simulation of clouds and their feedbacks. In each case, the amount of the dominant cloud types reduces in response to warming and the reduction is approximately proportional to the amount simulated in the present day. In deep convective regions, convective process parameters lead the spread among multiple contributing processes, with vegetation processes contributing as well for the land regions. Multiple parameters, such as boundary layer processes, drive stratocumulus regions. However, the low-thick clouds are systematically overestimated, suggesting a structural error in their process representations which would limit the efficacy of the constraint. The influence of convection is largely confined to the tropical deep convective regions in the present day but extends to mid-latitudes under warming. Because of this, contributing processes to the spread in the present-day and the response are different in the extra-tropics, making it much more difficult to establish links between the present-day and the feedback within the region. This suggests that identifying a constraint on convective processes in the tropics for the present-day simulations could constrain both the tropical feedbacks and feedbacks in the extra-tropics. A parameter representing deep-convective entrainment links the present-day tropical mean high cloud and clear-sky longwave flux to their feedbacks in our model, suggesting a potential process constraint from the observations. Understanding and improving the detailed processes controlling feedbacks is ultimately only possible in individual models. Different process-based constraints might be inferred for different models. The approach described here could usefully extended to other single model ensembles and the collective understandings could be valuable for improving model process and feedbacks more generally.