Analysis of the MODIS above-cloud aerosol retrieval algorithm using MCARS

The Multi-sensor Cloud and Aerosol Retrieval Simulator (MCARS) presently produces synthetic radiance data from Goddard Earth Observing System version 5 (GEOS-5) model output as if the Moderate Resolution Imaging Spectroradiometer (MODIS) were viewing a combination of atmospheric column inclusive of clouds, aerosols, and a variety of gases and land–ocean surface at a specific location. In this paper we use MCARS to study the MODIS Above-Cloud AEROsol retrieval algorithm (MOD06ACAERO). MOD06ACAERO is presently a regional research algorithm able to retrieve aerosol optical thickness over clouds, in particular absorbing biomass-burning aerosols overlying marine boundary layer clouds in the southeastern Atlantic Ocean. The algorithm's ability to provide aerosol information in cloudy conditions makes it a valuable source of information for modeling and climate studies in an area where current clear-sky-only operational MODIS aerosol retrievals effectively have a data gap between the months of June and October. We use MCARS for a verification and closure study of the MOD06ACAERO algorithm. The purpose of this study is to develop a set of constraints a model developer might use during assimilation of MOD06ACAERO data. Our simulations indicate that the MOD06ACAERO algorithm performs well for marine boundary layer clouds in the SE Atlantic provided some specific screening rules are observed. For the present study, a combination of five simulated MODIS data granules were used for a dataset of 13.5 million samples with known input conditions. When pixel retrieval uncertainty was less than 30 %, optical thickness of the underlying cloud layer was greater than 4, and scattering angle range within the cloud bow was excluded, MOD06ACAERO retrievals agreed with the underlying ground truth (GEOS-5 cloud and aerosol profiles used to generate the synthetic radiances) with a slope of 0.913, offset of 0.06, and RMSE=0.107. When only near-nadir pixels were considered (view zenith angle within ±20∘) the agreement with source data further improved (0.977, 0.051, and 0.096 respectively). Algorithm closure was examined using a single case out of the five used for verification. For closure, the MOD06ACAERO code was modified to use GEOS-5 temperature and moisture profiles as an ancillary. Agreement of MOD06ACAERO retrievals with source data for the closure study had a slope of 0.996 with an offset of −0.007 and RMSE of 0.097 at a pixel uncertainty level of less than 40 %, illustrating the benefits of high-quality ancillary atmospheric data for such retrievals.

Impact of Giant Sea Salt Aerosol Particles on Precipitation in Marine Cumuli and Stratocumuli: Lagrangian Cloud Model Simulations

The impact of giant sea salt aerosols released from breaking waves on rain formation in marine boundary layer clouds is studied using large eddy simulations (LES). We perform simulations of marine cumuli and stratocumuli for various concentrations of cloud condensation nuclei (CCN) and giant CCN (GCCN). Cloud microphysics are modeled with a Lagrangian method that provides key improvements in comparison to previous LES of GCCN that used Eulerian bin microphysics. We find that GCCN significantly increase precipitation in stratocumuli. This effect is strongest for low and moderate CCN concentrations. GCCN are found to have a smaller impact on precipitation formation in cumuli. These conclusions are in agreement with field measurements. We develop a simple parameterization of the effect of GCCN on precipitation, accretion, and autoconversion rates in marine stratocumuli.

A climatology of open and closed mesoscale cellular convection over the Southern Ocean derived from Himawari-8 observations

Marine atmospheric boundary layer clouds cover vast areas of the Southern Ocean (SO), where they are commonly organized into mesoscale cellular convection (MCC). Using three years of Himawari-8 geostationary satellite observations, open and closed MCC structures are identified using a hybrid convolutional neural network. The results of the climatology show that open MCC clouds are roughly uniformly distributed over the SO storm track across mid-latitudes, while closed MCC clouds are most predominant in the southeast Indian Ocean with a second maximum along the storm track. The ocean polar front, derived from ECMWF-ERA5 sea surface temperature gradients, is found to be aligned with the southern boundaries for both MCC types. Along the storm track, both closed and open MCCs are commonly located in post-frontal, cold air masses. The hourly classification of closed MCC reveals a pronounced daily cycle, with a peak occurring late night/early morning. Seasonally, the diurnal cycle of closed MCC is most intense during the summer months (DJF). Conversely, almost no diurnal cycle is evident for open MCC.

The effect of marine ice nucleating particles on mixed-phase clouds

Shallow marine mixed-phase clouds are important for the radiative balance, but modelling their formation and dynamics is challenging. These clouds depend on boundary layer turbulence and cloud top radiative cooling, which is related to the cloud phase. The fraction of frozen droplets depends on the availability of suitable ice nucleating particles (INPs), which initiate droplet freezing. While desert dust is the dominating INP type in most regions, remote boundary layer clouds are dependent on local marine INP emissions, which are often related to biogenic sources including phytoplankton. Here we use high resolution large eddy simulations to examine the potential effects of marine emissions on boundary layer INP concentrations and their effects on clouds. Surface emissions have a direct effect on INP concentration in a typical well-mixed boundary layer whereas a steep inversion can block the import of background INPs from the free troposphere. The importance of the marine source depends on the background INP concentration, so that marine emissions become dominant with low background concentrations. For the INP budget it is also important to account for INP recycling. Finally, with the high-resolution model we show how ice nucleation hotspots and high INPs concentrations are focused on updraught regions. Our results show that marine INP emissions contribute directly to the boundary layer INP budget and therefore have an influence on mixed-phase clouds.

Analysis of the MODIS Above-Cloud Aerosol Retrieval Algorithm Using MCARS

The Multi-sensor Cloud and Aerosol Retrieval Simulator (MCARS) presently produces synthetic radiance data from Goddard Earth Observing System version 5 (GEOS-5) model output as if the Moderate Resolution Imaging Spectroradiometer (MODIS) was viewing a combination of atmospheric column inclusive of clouds, aerosols and a variety of gases and land/ocean surface at a specific location. In this paper we use MCARS to study the MODIS Above-Cloud AEROsol retrieval algorithm (MOD06ACAERO). MOD06ACAERO is presently a regional research algorithm able to retrieve aerosol optical thickness over clouds, in particular absorbing biomass burning aerosols overlying marine boundary layer clouds in the Southeastern Atlantic Ocean. The algorithm's ability to provide aerosol information in cloudy conditions makes it a valuable source of information for modeling and climate studies in an area where current clear sky-only operational MODIS aerosol retrievals effectively have a data gap between the months of June and October. We use MCARS for a verification and closure study of the MOD06ACAERO algorithm. Our simulations indicate that the MOD06ACAERO algorithm performs well for marine boundary layer clouds in the SE Atlantic provided some specific screening rules are observed. For the present study, a combination of five simulated MODIS data granules was used for a dataset of 13.5 million samples with known input conditions. When pixel retrieval uncertainty was less than 30 %, optical thickness of the underlying cloud layer was greater than 4 and scattering angle range within the cloud bow was excluded, MOD06ACAERO retrievals agreed with the underlying ground truth (GEOS-5 cloud and aerosol profiles used to generate the synthetic radiances) with a slope of 0.913, offset of 0.06, and RMSE = 0.107. When only near-nadir pixels were considered (view zenith angle within +/−20 degrees) the agreement with source data further improved (0.977, 0.051 and 0.096 respectively). Algorithm closure was examined using a single case out of the five used for verification. For closure, the MOD06ACAERO code was modified to use GEOS-5 temperature and moisture profiles as ancillary. Agreement of MOD06ACAERO retrievals with source data for the closure study had a slope of 0.996 with offset −0.007 and RMSE of 0.097 at pixel uncertainty level of less than 40 %, illustrating the benefits of high-quality ancillary atmospheric data for such retrievals.

Exploring microphysical properties of marine boundary layer clouds through network analysis

<p>Aerosol-cloud interactions remain the largest uncertainty in assessments of anthropogenic climate forcing, while the complexity of these interactions require methods that enable abstractions and simplifications that allow their improved treatment in climate models. Marine boundary layer clouds are an important component of the climate system as their large albedo and spatial coverage strongly affect the planetary radiative balance. High resolution simulations of clouds provide an unprecedented understanding of the structure and behavior of these clouds in the marine atmosphere, but the amount of data is often too large and complex to be useful in climate simulations. Data reduction and inference methods provide a way that to reduce the complexity and dimensionality of datasets generated from high-resolution Large Eddy Simulations.</p><p>In this study we use network analysis, (the δ-Maps method) to study the complex interaction between liquid water, droplet number and vertical velocity in Large Eddy Simulations of Marine Boundary Layer clouds. δ-Maps identifies domains that are spatially contiguous and possibly overlapping and characterizes their connections and temporal interactions. The objective is to better understand microphysical properties of marine boundary layer clouds, and how they are impacted by the variability in aerosols. Here we will capture the dynamical structure of the cloud fields predicted by the MIMICA Large Eddy Simulation (LES) model. The networks inferred from the different simulation fields are compared between them (intra-comparisons) using perturbations in initial conditions and aerosol, using a set of four metrics. The networks are then evaluated for their differences, quantifying how much variability is inherent in the LES simulations versus the robust changes induced by the aerosol fields. </p>

Impact of 3D cloud structures on tropospheric NO2 column measurements from UV-VIS sounders

<p>Operational retrievals of tropospheric trace gases from space-borne spectrometers are made using 1D radiative transfer models. To minimize cloud effects generally only partially cloudy pixels are analysed using simplified cloud contamination treatments based on radiometric cloud fraction estimates and photon path length corrections based on oxygen collision pair (O2-O2) or O2A-absorption band measurements. In reality, however, the impact of clouds can be much more complex, involving unresolved sub-pixel clouds, scattering of clouds in neighbouring pixels, and cloud shadow effects, such that 3D radiation scattering from unresolved boundary layer clouds may give significant biases in the trace gas retrievals. In order to quantify this impact, we use the MYSTIC 3D radiative transfer model to generate synthetic data. The realistic 3D cloud fields, needed for MYSTIC input, are generated by the ICOsahedral Non-hydrostatic (ICON) atmosphere model for a region including Germany, the Netherlands and parts of other surrounding countries. The retrieval algorithm is applied to the synthetic data and comparison to the known input trace gas concentrations yields the retrieval error due to 3D cloud effects. <br>In this study, we study NO2, which is a key tropospheric trace gas measured by TROPOMI and the future atmospheric Sentinels (S4 and S5). The work starts with a sensitivity study for the simulations with a simple 2D box cloud. The influence of cloud parameters (e.g., cloud top height, cloud optical thickness), observation geometry, and spatial resolution are studied, and the most significant dependences of retrieval biases are identified and investigated. Several approaches to correct the NO2 retrieval in the cloud shadow are explored and ultimately applied to both synthetic data with realistic 3D clouds and real observations.</p>

Quantification of aerosol effects on marine boundary layer clouds with machine learning

<p>Satellite observations are used in regional machine learning models to quantify sensitivities of marine boundary-layer clouds (MBLC) to aerosol changes.</p><p>MBLCs make up a large part of the global cloud coverage as they are persistently present over more than 20% of the Earth’s oceans in the annual mean.They play an important role in Earth’s energy budget by reflecting solar radiation and interacting with thermal radiation from the surface, leading to a net cooling effect. Cloud properties and their radiative characteristics such as cloud albedo, horizontal and vertical extent, lifetime and precipitation susceptibility are dependent on environmental conditions. Aerosols in their role as condensation nuclei affect these cloud radiative properties through changes in the cloud droplet number concentration and subsequent cloud adjustments to this perturbation. However, the magnitude and sign of these effects remain among the largest uncertainties in future climate predictions.</p><p>In an effort to help improve these predictions a machine learning approach in combination with observational data is pursued:</p><p>Satellite observations from the collocated A-Train dataset (C3M) for 2006-2011 are used in combination with ECMWF atmospheric reanalysis data (ERA5) to train regional Gradient Boosting Regression Tree (GBRT) models to predict changes in key physical and radiative properties of MBLCs. The cloud droplet number concentration (N<sub>d</sub>) and the liquid water path (LWP) are simulated for the eastern subtropical oceans, which are characterised by a high annual coverage of MBLC due to the occurrence of semi-permanent stratocumulus sheets. Relative humidity above cloud, cloud top height and temperature below the cloud base and at the surface are identified as important predictors for both N<sub>d</sub> and LWP.  The impact of each predictor variable on the GBRT model's output is analysed using Shapley values as a method of explainable machine learning, providing novel sensitivity estimates that will improve process understanding and help constrain the parameterization of MBLC processes in Global Climate Models.</p>