Increasing the spatial resolution of cloud property retrievals from Meteosat SEVIRI by use of its high-resolution visible channel: Evaluation of candidate approaches with MODIS observations

This study presents and evaluates several candidate approaches for downscaling observations from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) in order to increase the horizontal resolution of subsequent cloud optical thickness (τ) and effective droplet radius (reff) retrievals from the native ≈ 3km×3km spatial resolution of the narrowband channels to ≈ 1km×1km. These methods make use of SEVIRI's coincident broadband high-resolution visible (HRV) channel. For four example cloud fields, the reliability of each downscaling algorithm is evaluated by means of collocated 1km×1km MODIS radiances, which are reprojected to the horizontal grid of the HRV channel and serve as reference for the evaluation. By using these radiances, smoothed with the modulation transfer function of the native SEVIRI channels, as retrieval input, the accuracy at the SEVIRI standard resolution can be evaluated and an objective comparison of the accuracy of the different downscaling algorithms can be made. For the example scenes considered in this study, it is shown that neglecting high-frequency variations below the SEVIRI standard resolution results in significant random absolute deviations of the retrieved τ and reff of up to ≈ 14 and ≈ 6μm, respectively, as well as biases. By error propagation, this also negatively impacts the reliability of the subsequent calculation of liquid water path (WL) and cloud droplet number concentration (ND), which exhibit deviations of up to ≈ 89gm-2 and ≈ 177cm-3, respectively. For τ , these deviations can be almost completely mitigated by the use of the HRV channel as a physical constraint and by applying most of the presented downscaling schemes. Uncertainties in retrieved reff at the native SEVIRI resolution are smaller, and the improvements from downscaling the observations are less obvious than for τ. Nonetheless, the right choice of downscaling scheme yields noticeable improvements in the retrieved reff. Furthermore, the improved reliability in retrieved cloud products results in significantly reduced uncertainties in derived WL and ND. In particular, one downscaling approach provides clear improvements for all cloud products compared to those obtained from SEVIRI's standard resolution and is recommended for future downscaling endeavors. This work advances efforts to mitigate impacts of scale mismatches among channels of multiresolution instruments on cloud retrievals. © Author(s) 2020

Evaluation of satellite-based aerosol datasets and the CAMS reanalysis over the ocean utilizing shipborne reference observations

Reliable reference measurements over the ocean are essential for the evaluation and improvement of satelliteand model-based aerosol datasets. Within the framework of the Maritime Aerosol Network, shipborne reference datasets have been collected over the Atlantic Ocean since 2004 with Microtops Sun photometers. These were recently complemented by measurements with the multi-spectral GUVis- 3511 shadowband radiometer during five cruises with the research vessel Polarstern. The aerosol optical depth (AOD) uncertainty estimate of both shipborne instruments of ±0:02 can be confirmed if the GUVis instrument is cross calibrated to the Microtops instrument to account for differences in calibration, and if an empirical correction to account for the broad shadowband as well as the effects of forward scattering is introduced. Based on these two datasets, a comprehensive evaluation of aerosol products from the Moderate Resolution Imaging Spectroradiometer (MODIS) flown on NASA's Earth Observing System satellites, the Spinning Enhanced Visible and Infrared Imager (SEVIRI) aboard the geostationary Meteosat satellite, and the Copernicus Atmosphere Monitoring Service reanalysis (CAMS RA) is presented. For this purpose, focus is given to the accuracy of the AOD at 630 nm in combination with the Ångström exponent (AE), discussed in the context of the ambient aerosol type. In general, the evaluation of MODIS AOD from the official level-2 aerosol products of C6.1 against the Microtops AOD product confirms that 76% of data points fall into the expected error limits given by previous validation studies. The SEVIRI-based AOD product exhibits a 25% larger scatter than the MODIS AOD products at the instrument's native spectral channels. Further, the comparison of CAMS RA and MODIS AOD versus the shipborne reference shows similar performance for both datasets, with some differences arising from the assimilation and model assumptions. When considering aerosol conditions, an overestimation of AE is found for scenes dominated by desert dust for MODIS and SEVIRI products versus the shipborne reference dataset. As the composition of the mixture of aerosol in satellite products is constrained by model assumptions, this highlights the importance of considering the aerosol type in evaluation studies for identifying problematic aspects. © Author(s) 2020

View angle dependence of MODIS liquid water path retrievals in warm oceanic clouds

We investigated the view angle dependence of domain mean Moderate Resolution Imaging Spectroradiometer (MODIS) liquid water path (LWP) and that of corresponding cloud optical thickness, effective radius, and liquid cloud fraction as proxy for plane-parallel retrieval biases. Independent Advanced Microwave Scanning Radiometer–EOS LWP was used to corroborate that the observed variations with sun-view geometry were not severely affected by seasonal/latitudinal changes in cloud properties. Microwave retrievals showed generally small (<10%) cross-swath variations. The view angle (cross-swath) dependence of MODIS optical thickness was weaker in backscatter than forward scatter directions and transitioned from mild ∩ shape to stronger ∪ shape as heterogeneity, sun angle, or latitude increased. The 2.2 µm effective radius variations always had a ∪ shape, which became pronounced and asymmetric toward forward scatter in the most heterogeneous clouds and/or at the lowest sun. Cloud fraction had the strongest and always ∪-shaped view angle dependence. As a result, in-cloud MODIS cloud liquid water path (CLWP) showed surprisingly good view angle (cross-swath) consistency, usually comparable to that of microwave retrievals, due to cancelation between optical thickness and effective radius biases. Larger (20–40%) nadir-relative increases were observed in the most extreme heterogeneity and sun angle bins, that is, typically in the polar regions, which, however, constituted only 3–8% of retrievals. The good consistency of MODIS in-cloud CLWP was lost for gridbox mean LWP, which was dominated by the strong cloud fraction increase with view angle. More worryingly, MODIS LWP exhibited significant and systematic absolute increases with heterogeneity and sun angle that is not present in microwave LWP

The sub-adiabatic model as a concept for evaluating the representation and radiative effects of low-level clouds in a high-resolution atmospheric model

The realistic representation of low-level clouds, including their radiative effects, in atmospheric models remains challenging. A sensitivity study is presented to establish a conceptual approach for the evaluation of low-level clouds and their radiative impact in a highly resolved atmospheric model. Considering simulations for six case days, the analysis supports the notion that the properties of clouds more closely match the assumptions of the sub-adiabatic rather than the vertically homogeneous cloud model, suggesting its use as the basis for evaluation. For the considered cases, 95.7% of the variance in cloud optical thickness is explained by the variance in the liquid water path, while the droplet number concentration and the sub-adiabatic fraction contribute only 3.5% and 0.2% to the total variance, respectively. A mean sub-adiabatic fraction of 0.45 is found, which exhibits strong inter-day variability. Applying a principal component analysis and subsequent varimax rotation to the considered set of nine properties, four dominating modes of variability are identified, which explain 97.7% of the total variance. The first and second components correspond to the cloud base and top height, and to liquid water path, optical thickness, and cloud geometrical extent, respectively, while the cloud droplet number concentration and the sub-adiabatic fraction are the strongest contributors to the third and fourth components. Using idealized offline radiative transfer calculations, it is confirmed that the shortwave and longwave cloud radiative effects exhibit little sensitivity to the vertical structure of clouds. This reconfirms, based on an unprecedented large set of highly resolved vertical cloud profiles, that the cloud optical thickness and the cloud top and bottom heights are the main factors dominating the shortwave and longwave radiative effect of clouds and should be evaluated together with radiative fluxes using observations to attribute model deficiencies in the radiative fluxes to deficiencies in the representation of clouds. Considering the different representations of cloud microphysical processes in atmospheric models, the analysis has been further extended and the deviations between the radiative impact of the single- and double-moment schemes are assessed. Contrasting the shortwave cloud radiative effect obtained from the double-moment scheme to that of a single-moment scheme, differences of about similar to 40 Wm(-2) and significant scatter are observed. The differences are attributable to a higher cloud albedo resulting from the high values of droplet number concentration in particular in the boundary layer predicted by the double-moment scheme, which reach median values of around similar to 600 cm(-3)

Metrics for the evaluation of warm convective cloud fields in a large-eddy simulation with Meteosat images

The representation of warm convective clouds in atmospheric models and satellite observations can considerably deviate from each other partly due to different spatial resolutions. This study aims to establish appropriate metrics to evaluate high-resolution simulations of convective clouds by the ICON Large-Eddy Model (ICON-LEM) with observations from Meteosat SEVIRI over Germany. The time series and frequency distributions of convective cloud fraction and liquid water path (LWP) are analyzed. Furthermore, the study focuses on size distributions and decorrelation scales of warm convective cloud fields. The investigated metrics possess a pronounced sensitivity to the apparent spatial resolution. At the fine spatial scale, the simulations show higher occurrence frequencies of large LWP values and a factor of two to four smaller convective cloud fractions. Coarse-graining of simulated fields to the optical resolution of Meteosat essentially removes the differences between the observed and simulated metrics. The distribution of simulated cloud sizes compares well with the observations and can be represented by a power law, with a moderate resolution sensitivity. A lower limit of cloud sizes is identified, which is 8–10 times the native grid resolution of the model. This likely marks the effective model resolution beyond which the scaling behaviour of considered metrics is not reliable, implying that a further increase in spatial resolution would be desirable to better resolve cloud processes below 1 km. It is finally shown that ICON-LEM is consistent with spatio-temporal decorrelation scales observed with Meteosat having values of 30 min and 7 km, if transferred to the true optical satellite resolution. However, the simulated Lagrangian decorrelation times drop to 10 min at 1 km resolution, a scale covered by the upcoming generation of geostationary satellites

Algorithms and uncertainties for the determination of multispectral irradiance components and aerosol optical depth from a shipborne rotating shadowband radiometer

The 19-channel rotating shadowband radiometer GUVis-3511 built by Biospherical Instruments provides automated shipborne measurements of the direct, diffuse and global spectral irradiance components without a requirement for platform stabilization. Several direct sun products, including spectral direct beam transmittance, aerosol optical depth, Ångström exponent and precipitable water, can be derived from these observations. The individual steps of the data analysis are described, and the different sources of uncertainty are discussed. The total uncertainty of the observed direct beam transmittances is estimated to be about 4% for most channels within a 95% confidence interval for shipborne operation. The calibration is identified as the dominating contribution to the total uncertainty. A comparison of direct beam transmittance with those obtained from a Cimel sunphotometer at a land site and a manually operated Microtops II sunphotometer on a ship is presented. Measurements deviate by less than 3 and 4% on land and on ship, respectively, for most channels and in agreement with our previous uncertainty estimate. These numbers demonstrate that the instrument is well suited for shipborne operation, and the applied methods for motion correction work accurately. Based on spectral direct beam transmittance, aerosol optical depth can be retrieved with an uncertainty of 0.02 for all channels within a 95% confidence interval. The different methods to account for Rayleigh scattering and gas absorption in our scheme and in the Aerosol Robotic Network processing for Cimel sunphotometers lead to minor deviations. Relying on the cross calibration of the 940 nm water vapor channel with the Cimel sunphotometer, the column amount of precipitable water can be estimated with an uncertainty of ±0.034 cm

Increasing Resolution and Resolving Convection Improve the Simulation of Cloud-Radiative Effects Over the North Atlantic

Clouds interact with atmospheric radiation and substantially modify the Earth's energy budget. Cloud formation processes occur over a vast range of spatial and temporal scales, which make their thorough numerical representation challenging. Therefore, the impact of parameter choices for simulations of cloud-radiative effects is assessed in the current study. Numerical experiments are carried out using the ICOsahedral Nonhydrostatic (ICON) model with varying grid spacings between 2.5 and 80 km and with different subgrid-scale parameterization approaches. Simulations are performed over the North Atlantic with either one-moment or two-moment microphysics and with convection being parameterized or explicitly resolved by grid-scale dynamics. Simulated cloud-radiative effects are compared to products derived from Meteosat measurements. Furthermore, a sophisticated cloud classification algorithm is applied to understand the differences and dependencies of simulated and observed cloud-radiative effects. The cloud classification algorithm developed for the satellite observations is also applied to the simulation output based on synthetic infrared brightness temperatures, a novel approach that is not impacted by changing insolation and guarantees a consistent and fair comparison. It is found that flux biases originate equally from clear-sky and cloudy parts of the radiation field. Simulated cloud amounts and cloud-radiative effects are dominated by marine, shallow clouds, and their behavior is highly resolution dependent. Bias compensation between shortwave and longwave flux biases, seen in the coarser simulations, is significantly diminished for higher resolutions. Based on the analysis results, it is argued that cloud-microphysical and cloud-radiative properties have to be adjusted to further improve agreement with observed cloud-radiative effects. © 2020. The Authors

The day-to-day co-variability between mineral dust and cloud glaciation: a proxy for heterogeneous freezing

To estimate the global co-variability between mineral dust aerosol and cloud glaciation, we combined an aerosol model reanalysis with satellite retrievals of cloud thermodynamic phase. We used the CALIPSO-GOCCP product from the A-Train satellite constellation to assess whether clouds are composed of liquid or ice and the MACC reanalysis to estimate the dust mixing ratio in the atmosphere. Night-time retrievals within a temperature range from +3 to −42 ∘C for the period 2007–2010 were included. The results confirm that the cloud thermodynamic phase is highly dependent on temperature and latitude. However, at middle and high latitudes, at equal temperature and within narrow constraints for humidity and static stability, the average frequency of fully glaciated clouds increases by +5 to +10 % for higher mineral dust mixing ratios. The discrimination between humidity and stability regimes reduced the confounding influence of meteorology on the observed relationship between dust and cloud ice. Furthermore, for days with similar mixing ratios of mineral dust, the cloud ice occurrence frequency in the Northern Hemisphere was found to be higher than in the Southern Hemisphere at −30 ∘C but lower at −15 ∘C. This contrast may suggest a difference in the susceptibility of cloud glaciation to the presence of dust. Based on previous studies, the differences at −15 ∘C could be explained by higher feldspar fractions in the Southern Hemisphere, while the higher freezing efficiency of clay minerals in the Northern Hemisphere may explain the differences at −30 ∘C