More than one hundred days were simulated over very large domains with fine (0.156 km to 2.5 km) grid spacing for realistic conditions to test the hypothesis that storm (kilometer) and large-eddy (hectometer) resolving simulations would provide an improved representation of clouds and precipitation in atmospheric simulations. At scales that resolve convective storms (storm-resolving for short) scales, the vertical velocity variance becomes resolved and a better physical basis is achieved for representing clouds and precipitation. Similar to past studies we find an improved representation of precipitation at kilometer scales, as compared to models with parameterised convection. The main precipitation features (location, diurnal cycle and spatial propagation) are well captured already at kilometer scales, and refining resolution to hectometer scales does not substantially change the simulations in these respects. It does, however, lead to a reduction in the precipitation on the timescales considered-most notably over the Tropical ocean. Changes in the distribution of precipitation, with less frequent extremes are also found in simulations incorporating hecto-meter scales. Hectometer scales appear more important for the representation of clouds, and make it possible to capture many important aspects of the cloud field, from the vertical distribution of cloud cover, to the distribution of cloud sizes, to the diel (daily) cycle. Qualitative improvements, particularly in the ability to differentiate cumulus from stratiform clouds, are seen when reducing the grid spacing from kilometer to hectometer scales. At the hectometer scale new challenges arise, but the similarity of observed and simulated scales, and the more direct 1 connection between the circulation and the unconstrained degrees of freedom make these challenges less daunting. This quality, combined with an already improved simulation as compared to more parameterised models, underpins our conviction that the use and further development of storm-resolving models offers exciting opportunities for advancing understanding of climate and climate change.

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

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

Abstract. A total of 14 years (September 2002 to September 2016) of Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) monthly mean cloud data are used to quantify possible changes in the cloud vertical distribution over the tropical Atlantic. For the analysis multiple linear regression techniques are used. For the investigated time period significant linear changes were found in the domain-averaged cloud-top height (CTH) (−178 m per decade), the high-cloud fraction (HCF) (−0.0006 per decade), and the low-cloud amount (0.001 per decade). The interannual variability of the time series (especially CTH and HCF) is highly influenced by the El Niño–Southern Oscillation (ENSO). Separating the time series into two phases, we quantified the linear change associated with the transition from more La Niña-like conditions to a phase with El Niño conditions (Phase 2) and vice versa (Phase 1). The transition from negative to positive ENSO conditions was related to a decrease in total cloud fraction (TCF) (−0.018 per decade; not significant) due to a reduction in the high-cloud amount (−0.024 per decade; significant). Observed anomalies in the mean CTH were found to be mainly caused by changes in HCF rather than by anomalies in the height of cloud tops themselves. Using the large-scale vertical motion ω at 500 hPa (from ERA-Interim ECMWF reanalysis data), the observed anomalies were linked to ENSO-induced changes in the atmospheric large-scale dynamics. The most significant and largest changes were found in regions with strong large-scale upward movements near the Equator. Despite the fact that with passive imagers such as MODIS it is not possible to vertically resolve clouds, this study shows the great potential for large-scale analysis of possible changes in the cloud vertical distribution due to the changing climate by using vertically resolved cloud cover and linking those changes to large-scale dynamics using other observations or model data.

Top-of-atmosphere (TOA) shortwave (SW) angular distribution models (ADMs) approximate – per angular direction of an imagined upward hemisphere – the intensity of sunlight scattered back from a specific Earth-atmosphere scene. ADMs are, thus, critical when converting satellite-borne broadband radiometry into estimated radiative fluxes. This paper applies a set of newly developed ADMs with a more refined scene definition and demonstrates tenable changes in estimated fluxes compared to currently operational ADMs. Newly developed ADMs use a semi-physical framework to consider cloud-top effective radius, , and above-cloud water vapor, ACWV, in addition to accounting for surface wind speed and clouds’ phase, fraction, and optical depth. In effect, instantaneous TOA SWfluxes for marine liquid-phase clouds had the largest flux differences (of up to 25 W m−2) for lower solar zenith angles and cloud optical depth greater than 10 due to extremes in or ACWV. In regions where clouds had persistently extreme levels of (here mostly for <7μm and >15μm) or ACWV, instantaneous fluxes estimated from Aqua, Terra, and Meteosat 8 and 9 satellites using the two ADMs differed systematically, resulting in significant deviations in daily mean fluxes (up to ±10 W m−2) and monthly mean fluxes (up to ±5 W m−2). Flux estimates using newly developed, semi-physical ADMs may contribute to a better understanding of solar fluxes over low-level clouds. It remains to be seen whether aerosol indirect effects are impacted by these updates.