Condensation clouds in substellar atmospheres have been widely inferred from spectra and photometric variability. Up until now, their horizontally averaged vertical distribution and mean particle ...size have been largely characterized using models, one of which is the eddy diffusion-sedimentation model from Ackerman and Marley that relies on a sedimentation efficiency parameter, fsed, to determine the vertical extent of clouds in the atmosphere. However, the physical processes controlling the vertical structure of clouds in substellar atmospheres are not well understood. In this work, we derive trends in fsed across a large range of eddy diffusivities (Kzz), gravities, material properties, and cloud formation pathways by fitting cloud distributions calculated by a more detailed cloud microphysics model. We find that fsed is dependent on Kzz, but not gravity, when Kzz is held constant. fsed is most sensitive to the nucleation rate of cloud particles, as determined by material properties like surface energy and molecular weight. High surface energy materials form fewer, larger cloud particles, leading to large fsed (>1), and vice versa for materials with low surface energy. For cloud formation via heterogeneous nucleation, fsed is sensitive to the condensation nuclei flux and radius, connecting cloud formation in substellar atmospheres to the objects' formation environments and other atmospheric aerosols. These insights could lead to improved cloud models that help us better understand substellar atmospheres. For example, we demonstrate that fsed could increase with increasing cloud base depth in an atmosphere, shedding light on the nature of the brown dwarf L/T transition.
This study attempts a new identification of mechanisms of secondary ice production (SIP) based on the observation of small faceted ice crystals (hexagonal plates or columns) with typical sizes ...smaller than 100µm. Due to their young age, such small ice crystals can be used as tracers for identifying the conditions for SIP. Observations reported here were conducted in oceanic tropical mesoscale convective systems (MCSs) and midlatitude frontal clouds in the temperature range from 0 to −15◦C and heavily seeded by aged ice particles. It was found that in both MCSs and frontal clouds, SIP was observed right above the melting layer and extended to higher altitudes with colder temperatures. The roles of six possible mechanisms to generate the SIP particles are assessed using additional observations. In most observed SIP cases, small secondary ice particles spatially correlated with liquid-phase, vertical updrafts and aged rimed ice particles. However, in many cases, neither graupel nor liquid drops were observed in the SIP regions, and therefore, the conditions for an active Hallett–Mossop process were not met. In many cases, large concentrations of small pristine ice particles were observed right above the melting layer, starting at temperatures as warm as −0.5◦C. It is proposed that the initiation of SIP above the melting layer is stimulated by the recirculation of large liquid drops through the melting layer with convective turbulent updrafts. After re-entering a super cooled environment above the melting layer, they impact with aged ice, freeze, and shatter. The size of the splinters generated during SIP was estimated as 10µm or less. A principal conclusion of this work is that only the freezing drop-shattering mechanism could be clearly supported by the airborne in situ observations.
Recent studies have shown that, in response to a surface warming, the marine tropical low-cloud cover (LCC) as observed by passive-sensor satellites substantially decreases, therefore generating a ...smaller negative value of the top-of-the-atmosphere (TOA) cloud radiative effect (CRE). Here we study the LCC and CRE interannual changes in response to sea surface temperature (SST) forcings in the GISS model E2 climate model, a developmental version of the GISS model E3 climate model, and in 12 other climate models, as a function of their ability to represent the vertical structure of the cloud response to SST change against 10 years of CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) observations. The more realistic models (those that satisfy the observational constraint) capture the observed interannual LCC change quite well (deltaLCC/deltaSST= -3.49±1.01%K negative 1 superscript vs. deltaLCC/deltaSSTsubscript obs= -3.59±0.28%Knegative 1 superscript) while the others largely underestimate it (deltaLCC/deltaSST= -1.32±1.28%Knegative 1 superscript). Consequently, the more realistic models simulate more positive shortwave (SW) feedback (deltaCRE/deltaSST=2.60±1.13Wmnegative 2 superscript Knegative 1 superscript) than the less realistic models (CRE/SST=0.87±2.63Wm2K1), in better agreement with the observations (deltaCRE/deltaSSTsubscript obs=3±0.26Wmnegative 2 superscript Knegative 1 superscript ), although slightly underestimated. The ability of the models to represent moist processes within the planetary boundary layer (PBL) and produce persistent stratocumulus (Sc) decks appears crucial to replicating the observed relationship between clouds, radiation and surface temperature. This relationship is different depending on the type of low clouds in the observations. Over stratocumulus regions, cloud-top height increases slightly with SST, accompanied by a large decrease in cloud fraction, whereas over trade cumulus (Cu) regions, cloud fraction decreases everywhere, to a smaller extent.
This study investigates effects of drizzle and cloud horizontal inhomogeneity on cloud effective radius (re) retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS). In order to ...identify the relative importance of various factors, we developed a MODIS cloud property retrieval simulator based on the combination of large‐eddy simulations (LES) and radiative transfer computations. The case studies based on synthetic LES cloud fields indicate that at high spatial resolution (∼100 m) 3‐D radiative transfer effects, such as illumination and shadowing, can induce significant differences between retrievals ofre based on reflectance at 2.1 μm (re,2.1) and 3.7 μm (re,3.7). It is also found that 3‐D effects tend to have stronger impact onre,2.1 than re,3.7, leading to positive difference between the two (Δre,3.7−2.1) from illumination and negative Δre,3.7−2.1from shadowing. The cancellation of opposing 3‐D effects leads to overall reasonable agreement betweenre,2.1 and re,3.7 at high spatial resolution as far as domain averages are concerned. At resolutions similar to MODIS, however, re,2.1 is systematically larger than re,3.7when averaged over the LES domain, with the difference exhibiting a threshold‐like dependence on bothre,2.1and an index of the sub‐pixel variability in reflectance (Hσ), consistent with MODIS observations. In the LES cases studied, drizzle does not strongly impact reretrievals at either wavelength. It is also found that opposing 3‐D radiative transfer effects partly cancel each other when cloud reflectance is aggregated from high spatial resolution to MODIS resolution, resulting in a weaker net impact of 3‐D radiative effects onre retrievals. The large difference at MODIS resolution between re,3.7 and re,2.1 for highly inhomogeneous pixels with Hσ> 0.4 can be largely attributed to what we refer to as the “plane‐parallelrebias,” which is attributable to the impact of sub‐pixel level horizontal variability of cloud optical thickness onre retrievals and is greater for re,2.1 than re,3.7. These results suggest that there are substantial uncertainties attributable to 3‐D radiative effects and plane‐parallelre bias in the MODIS re,2.1retrievals for pixels with strong sub‐pixel scale variability, and theHσ index can be used to identify these uncertainties.
Key Points
Sub‐pixel inhomogeneity can cause strong uncertainty in MODIS cloud re retrieval
There is likely significant bias in the MODIS re,2.1 retrievals.
Sub‐pixel cloud inhomogeneity index can be used to identify this bias
Marine cold air outbreaks (CAOs) commonly form overcast cloud decks that transition into broken cloud fields downwind, dramatically altering the local radiation budget. In this study, we investigate ...the impact of frozen hydrometeors on these transitions. We focus on a CAO case in the NW Atlantic, the location of the multi-year flight campaign ACTIVATE (Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment). We use MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, version 2) reanalysis fields to drive large eddy simulations with mixed-phase two-moment microphysics in a Lagrangian framework. We find that transitions are triggered by substantial rain (rainwater paths >25 g m−2), and only simulations that allow for aerosol depletion result in sustained breakups, as observed. Using a range of diagnostic ice nucleating particle concentrations, Ninp, we find that increasing ice progressively accelerates transitions, thus abbreviating the overcast state. Ice particles affect the cloud-topped boundary layer evolution, primarily through riming-related processes prior to substantial rain, leading to (1) a reduction in cloud liquid water, (2) early consumption of cloud condensation nuclei, and (3) early and light precipitation cooling and moistening below cloud. We refer to these three effects collectively as “preconditioning by riming”. Greater boundary layer aerosol concentrations available as cloud condensation nuclei (CCN) delay the onset of substantial rain. However, cloud breakup and low CCN concentration final stages are found to be inevitable in this case, due, primarily, to liquid water path buildup. An ice-modulated cloud transition speed suggests the possibility of a negative cloud–climate feedback. To address prevailing uncertainties in the model representation of mixed-phase processes, the magnitude of ice formation and riming impacts and, thereby, the strength of an associated negative cloud–climate feedback process, requires further observational evaluation by targeting riming hot spots with in situ imaging probes that allow for both the characterization of ice particles and abundance of supercooled droplets.
Here, we compare the available data to the predictions of our own set of atmospheric and evolution models that have been extensively tested against observations of field L and T dwarfs, including the ...reddest L dwarfs. Unlike some previous studies, we require mutually consistent choices for effective temperature, gravity, cloud properties, and planetary radius. We find that the cloud properties of the HR 8799 planets are not unusual but rather follow previously recognized trends, including a gravity dependence on the temperature of the L to T spectral transition-some reasons for which we discuss. We find that the inferred mass of planet b is highly sensitive to whether or not we include the H- and the K-band spectrum in our analysis. Solutions for planets c and d are consistent with the generally accepted constraints on the age of the primary star and orbital dynamics.
Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size ...distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part 1 addressed collision-coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision-coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.
We present an algorithm for the retrieval of cloud droplet size distribution parameters (effective radius and variance) from the Research Scanning Polarimeter (RSP) measurements. The RSP is an ...airborne prototype for the Aerosol Polarimetery Sensor (APS), which was on-board of the NASA Glory satellite. This instrument measures both polarized and total reflectance in 9 spectral channels with central wavelengths ranging from 410 to 2260nm. The cloud droplet size retrievals use the polarized reflectance in the scattering angle range between 135° and 165°, where they exhibit the sharply defined structure known as the rain- or cloud-bow. The shape of the rainbow is determined mainly by the single scattering properties of cloud particles. This significantly simplifies both forward modeling and inversions, while also substantially reducing uncertainties caused by the aerosol loading and possible presence of undetected clouds nearby. In this study we present the accuracy evaluation of our algorithm based on the results of sensitivity tests performed using realistic simulated cloud radiation fields.
► We present an algorithm for retrieval of cloud droplet size from polarized rainbow. ► Polarized rainbow shape is determined by single scattering of cloud droplets. ► The algorithm is designed for the airborne Research Scanning Polarimeter (RSP). ► RSP measures polarized and total reflectance at 9 wavelengths between 410 to 2260nm. ► We evaluate the algorithm using tests on simulated cloud radiation fields.
Some of the global warming from anthropogenic greenhouse gases is offset by increased reflection of solar radiation by clouds with smaller droplets that form in air polluted with aerosol particles ...that serve as cloud condensation nuclei. The resulting cooling tendency, termed the indirect aerosol forcing, is thought to be comparable in magnitude to the forcing by anthropogenic CO2, but it is difficult to estimate because the physical processes that determine global aerosol and cloud populations are poorly understood. Smaller cloud droplets not only reflect sunlight more effectively, but also inhibit precipitation, which is expected to result in increased cloud water. Such an increase in cloud water would result in even more reflective clouds, further increasing the indirect forcing. Marine boundary-layer clouds polluted by aerosol particles, however, are not generally observed to hold more water. Here we simulate stratocumulus clouds with a fluid dynamics model that includes detailed treatments of cloud microphysics and radiative transfer. Our simulations show that the response of cloud water to suppression of precipitation from increased droplet concentrations is determined by a competition between moistening from decreased surface precipitation and drying from increased entrainment of overlying air. Only when the overlying air is humid or droplet concentrations are very low does sufficient precipitation reach the surface to allow cloud water to increase with droplet concentrations. Otherwise, the response of cloud water to aerosol-induced suppression of precipitation is dominated by enhanced entrainment of overlying dry air. In this scenario, cloud water is reduced as droplet concentrations increase, which diminishes the indirect climate forcing.
We present a generalization of the binary-value Markovian model previously used for statistical characterization of cloud masks to a continuous-value model describing 1D fields of cloud optical ...thickness (COT). This model has simple functional expressions and is specified by four parameters: the cloud fraction, the autocorrelation (scale) length, and the two parameters of the normalized probability density function of (non-zero) COT values (this PDF is assumed to have gamma-distribution form). Cloud masks derived from this model by separation between the values above and below some threshold in COT appear to have the same statistical properties as in binary-value model described in our previous publications. We demonstrate the ability of our model to generate examples of various cloud-field types by using it to statistically imitate actual cloud observations made by the Research Scanning Polarimeter (RSP) during two field experiments.