Bulk Microphysical Properties Research Articles

Analysis of Pre-Monsoon Convective Systems over a Tropical Coastal Region Using C-Band Polarimetric Radar, Satellite and Numerical Simulation

Analysis of pre-monsoon convective systems over the southern peninsular India has been performed using C-band radar and numerical simulation. Statistics on the radar polarimetric measurements show that the distribution of differential reflectivity (Zdr) and specific differential phase (Kdp) have much higher spread over convective regions. The distribution of Kdp is almost uniform across the vertical over the stratiform regions. The mean profile of Zdr over stratiform regions shows a distinct local maxima near melting level. A comprehensive analysis has been done on an isolated deep convective system on 13 May 2018. Plan position indicator (PPI) diagrams and satellite measured cloud top temperature demonstrate that pre-monsoon deep convective systems can develop very rapidly within a very short span of time over the region. Heavy precipitation near the surface is reflected in the high value of Kdp (>5° km−1). High values of Zdr (>3 dB) were measured at lower levels indicating the oblate shape of bigger raindrops. A fuzzy logic-based hydrometeor identification algorithm has been applied with five variables (Zh, Zdr, ρhv, Kdp, and T) to understand the bulk microphysical properties at different heights within the storm. The presence of bigger graupel particles near the melting layer indicates strong updrafts within the convective core regions. The vertical ice hydrometeor signifies the existence of a strong electric field causing them to align vertically. Numerical simulation with the spectral bin microphysics (SBM) scheme could produce most of the features of the storm reasonably well. In particular, the simulated reflectivity, graupel mixing ratio and rainfall were in good agreement with the observed values.

Snow microphysical retrieval from the NASA D3R radar during ICE-POP 2018

Abstract. A method is developed to use both polarimetric and dual-frequency radar measurements to retrieve microphysical properties of falling snow. It is applied to the Ku- and Ka-band measurements of the NASA dual-polarization, dual-frequency Doppler radar (D3R) obtained during the International Collaborative Experiments for PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign and incorporates the Atmospheric Radiative Transfer Simulator (ARTS) microwave single-scattering property database for oriented particles. The retrieval uses optimal estimation to solve for several parameters that describe the particle size distribution (PSD), relative contribution of pristine, aggregate, and rimed ice species, and the orientation distribution along an entire radial simultaneously. Examination of Jacobian matrices and averaging kernels shows that the dual-wavelength ratio (DWR) measurements provide information regarding the characteristic particle size, and to a lesser extent, the rime fraction and shape parameter of the size distribution, whereas the polarimetric measurements provide information regarding the mass fraction of pristine particles and their characteristic size and orientation distribution. Thus, by combining the dual-frequency and polarimetric measurements, some ambiguities can be resolved that should allow a better determination of the PSD and bulk microphysical properties (e.g., snowfall rate) than can be retrieved from single-frequency polarimetric measurements or dual-frequency, single-polarization measurements. The D3R ICE-POP retrievals were validated using Precipitation Imaging Package (PIP) and Pluvio weighing gauge measurements taken nearby at the May Hills ground site. The PIP measures the snow PSD directly, and its measurements can be used to derived the snowfall rate (volumetric and water equivalent), mean volume-weighted particle size, and effective density, as well as particle aspect ratio and orientation. Four retrieval experiments were performed to evaluate the utility of different measurement combinations: Ku-only, DWR-only, Ku-pol, and All-obs. In terms of correlation, the volumetric snowfall rate (r=0.95) and snow water equivalent rate (r=0.92) were best retrieved by the Ku-pol method, while the DWR-only method had the lowest magnitude bias for these parameters (−31 % and −8 %, respectively). The methods that incorporated DWR also had the best correlation to particle size (r=0.74 and r=0.71 for DWR-only and All-obs, respectively), although none of the methods retrieved density particularly well (r=0.43 for All-obs). The ability of the measurements to retrieve mean aspect ratio was also inconclusive, although the polarimetric methods (Ku-pol and All-obs) had reduced biases and mean absolute error (MAE) relative to the Ku-only and DWR-only methods. The significant biases in particle size and snowfall rate appeared to be related to biases in the measured DWR, emphasizing the need for accurate DWR measurements and frequent calibration in future D3R deployments.

Regional variation of electrical and lightning properties of thunderclouds during the pre-monsoon season over the north-eastern and eastern part of India

Regional variation of electrical and lightning properties of thunderclouds is investigated by using simultaneous observations of electric field mill (EFM-100) and lightning detector (LD-350) in synergy with their radiance and bulk microphysical properties. The study is carried out over Kohima in the north-eastern and over Rampurhat in the eastern part of India (20o-30o N, 85o-96o E), during the pre-monsoon season (March–May) of 2017–2018. The occurrence of EFM-100 detected thunderclouds is found to be more over Kohima (111) compared to Rampurhat (88). Whereas, the average number of lightning flash is much higher over Rampurhat (~542 flash/thundercloud) compared to Kohima (~32 flash/thundercloud). Over both the regions, the maximum occurrence of thunderclouds is found during afternoon to evening hours. The upper whisker value in the box plot shows that thunderclouds over Rampurhat are associated with higher CAPE value (Kohima: 1100 J Kg−1; Rampurhat: 5000 J Kg−1). The diurnal variation of CAPE supports relatively early occurrence of thunderclouds over Kohima. Thunderclouds over Rampurhat persists for relatively longer duration. It is observed that overall percentage occurrence of CG lightning is more over Kohima (Kohima: −CG: 28.09% & +CG: 09.74%; Rampurhat (–CG: 12.33% & +CG: 08.60%), with relatively higher −CG lightning compared to +CG, over both the stations. On the contrary, the occurrence of IC lightning is more over Rampurhat (Rampurhat: −IC: 31.75% & +IC: 47.30%; Kohima: −IC: 36.44% & +IC: 25.72%). Over Kohima occurrence of −IC is more than the +IC, whereas over Rampurhat, the occurrence of +IC is more than −IC. The overall higher occurrence of total negative and positive discharges over Kohima and Rampurhat respectively is consistence with EFM observations. The prevalence of positive/negative discharge suggests the presence of a tripole charge structure within thunderclouds over both the regions. The upper whisker value in the box plots revels higher CG flash density, IC:CG and larger log(area) over Rampurhat (Kohima: 3.0 × 10−3 flash km−2 h−1 and Rampurhat: 4.25 × 10−3 flash km−2 h−1; Kohima: 7 and Rampurhat: 15; Kohima: 10.5 km2 and Rampurhat: 11.5 km2). The lower whisker value of cloud top temperature reveals prevalence of colder thunderclouds over Rampurhat (Kohima: 220 K; Rampurhat: 184 K). The thunderclouds over Rampurhat are found to be associated with stronger updrafts and mixed-phase processes. The study reveals that, over these two regions, the electrical and lightning properties have distinct characteristics.

A physically based raindrop–cloud droplet accretion parametrization for use in bulk microphysics schemes

AbstractThe accretion process in bulk cloud microphysics schemes can be parametrized using the stochastic collection equation (SCE). In this study, the collection efficiency of each raindrop–cloud droplet pair is applied to the SCE to derive a new accretion parametrization that considers a strong variability of accretion rate depending on the cloud droplet and raindrop size distributions. To evaluate the new accretion parametrization (NP), it is implemented into a cloud‐resolving model, replacing the original accretion parametrization (OP) based on the continuous collection equation. In the idealized simulations of deep convective clouds, NP predicts overall larger accretion rates and smaller autoconversion rates than OP. The resultant high accretion/autoconversion rate ratio in NP increases the mean raindrop size. This induces faster sedimentation of raindrops that is associated with the earlier onset of surface precipitation and also weakens the evaporation cooling of raindrops that can affect the thermodynamics and dynamics of clouds. Meanwhile, for a given pair of rainwater and cloud water mass contents, the accretion rates in NP have a broad distribution while those in OP are less variant, suggesting that the dependence of the accretion rates on bulk microphysical properties other than the mass contents also needs to be properly considered in accretion parametrizations. Most of the aforementioned differences between the two accretion parametrizations found in the idealized simulations are also found in the real‐case simulations of a precipitation event over Bangladesh but the differences are smaller, and the spatial distribution of the accumulated precipitation amount is relatively well predicted by NP compared to OP.

Microphysical Properties of Generating Cells Over the Southern Ocean: Results From SOCRATES

AbstractThe bulk microphysical properties and number distribution functions (N(D)) of supercooled liquid water (SLW) and ice inside and between ubiquitous generating cells (GCs) observed over the Southern Ocean (SO) during the Southern Ocean Clouds Radiation Aerosol Transport Experimental Study (SOCRATES) measured by in situ cloud probes onboard the NCAR/NSF G‐V aircraft are compared. SLW was detected inside all GCs with an average liquid water content of 0.31 ± 0.19 g m−3, 11% larger than values between GCs. The N(D) of droplets (maximum dimension D < 50 μm) inside and between GCs had only slight differences. For ice particles, on the other hand, the mean concentration (median mass diameter) with D > 200 μm inside GCs was 2.0 ± 3.3 L−1 (323 ± 263 μm), 65% (37%) larger than values outside GCs. As D increases, the percentage differences became larger (up to ~500%). The more and larger ice particles inside GCs suggest the GC updrafts provide a favorable environment for particle growth by deposition and riming and that mixing processes are less efficient at redistributing larger particles. The horizontal scale of observed GCs ranged from 200 to 600 m with a mean of 395 ± 162 m, smaller than GC widths observed in previous studies. This study expands knowledge of the microphysical properties and processes acting in GCs over a wider range of conditions than previously available.

Analytical Solutions for Precipitation Size Distributions at Steady State

Abstract Analytical solutions are derived for the steady-state size distributions of precipitating rain and snow particles assuming growth via collection of suspended cloud particles. Application of the Liouville equation to the transfer of precipitating mass through size bins in a “cascade” yields a characteristic gamma distribution with a Marshall–Palmer exponential tail with respect to particle diameter. For rain, the principle parameters controlling size distribution shape are cloud droplet liquid water path and cloud updraft speed. Stronger updrafts lead to greater concentrations of large precipitating drops and a peak in the size distribution. The solutions provide a means for relating size distributions measured in the air or on the ground to cloud bulk microphysical and dynamic properties.

A novel approach for characterizing the variability in mass–dimension relationships: results from MC3E

Abstract. Mass–dimension (m–D) relationships determining bulk microphysical properties such as total water content (TWC) and radar reflectivity factor (Z) from particle size distributions are used in both numerical models and remote sensing retrievals. The a and b coefficients representing m=aDb relationships, however, can vary significantly depending on meteorological conditions, particle habits, the definition of particle maximum dimension, the probes used to obtain the data, techniques used to process the cloud probe data, and other unknown reasons. Thus, considering a range of a,b coefficients may be more applicable for use in numerical models and remote sensing retrievals. Microphysical data collected by two-dimensional optical array probes (OAPs) installed on the University of North Dakota (UND) Citation aircraft during the Mid-latitude Continental Convective Clouds Experiment (MC3E) were used in conjunction with TWC data from a Nevzorov probe and ground-based S-band radar data to determine a and b using a technique that minimizes the chi-square difference between the TWC and Z derived from the OAPs and those directly measured by a TWC probe and radar. All a and b values within a specified tolerance were regarded as equally plausible solutions. Of the 16 near-constant-temperature flight legs analyzed during the 25 April, 20 May, and 23 May 2011 events, the derived surfaces of solutions on the first 2 days where the aircraft-sampled stratiform cloud had a larger range in a and b for lower temperature environments that correspond to less variability in N(D), TWC, and Z for a flight leg. Because different regions of the storm were sampled on 23 May, differences in the variability in N(D), TWC, and Z influenced the distribution of chi-square values in the (a,b) phase space and the specified tolerance in a way that yielded 2.8 times fewer plausible solutions compared to the flight legs on the other dates. These findings show the importance of representing the variability in a,b coefficients for numerical modeling and remote sensing studies, rather than assuming fixed values, as well as the need to further explore how these surfaces depend on environmental conditions in clouds containing ice hydrometeors.

Quantifying uncertainties in radar forward models through a comparison between CloudSat and SPartICus reflectivity factors

AbstractInterpretations of remote sensing measurements collected in sample volumes containing ice‐phase hydrometeors are very sensitive to assumptions regarding the distributions of mass with ice crystal dimension, otherwise known as mass‐dimensional or m‐D relationships. How these microphysical characteristics vary in nature is highly uncertain, resulting in significant uncertainty in algorithms that attempt to derive bulk microphysical properties from remote sensing measurements. This uncertainty extends to radar reflectivity factors forward calculated from model output because the statistics of the actual m‐D in nature is not known. To investigate the variability in m‐D relationships in cirrus clouds, reflectivity factors measured by CloudSat are combined with particle size distributions (PSDs) collected by coincident in situ aircraft by using an optimal estimation‐based (OE) retrieval of the m‐D power law. The PSDs were collected by 12 flights of the Stratton Park Engineering Company Learjet during the Small Particles in Cirrus campaign. We find that no specific habit emerges as preferred, and instead, we find that the microphysical characteristics of ice crystal populations tend to be distributed over a continuum‐defying simple categorization. With the uncertainties derived from the OE algorithm, the uncertainties in forward‐modeled backscatter cross section and, in turn, radar reflectivity is calculated by using a bootstrapping technique, allowing us to infer the uncertainties in forward‐modeled radar reflectivity that would be appropriately applied to remote sensing simulator algorithms.

A study of thunderstorm microphysical properties and lightning flash counts associated with terrestrial gamma‐ray flashes

AbstractThe terrestrial gamma ray flash (TGF) is an emission of highly energetic radiation produced by or at least in close association with lightning. Previous investigations attempted to isolate the production mechanisms and production altitude(s) of TGFs as well as macrophysical characteristics, while thunderstorm microphysical characteristics were largely ignored. This investigation into thunderstorms and their hydrometeor and flash characteristics utilize temporal and spatial coincident satellite passes between the Reuven Ramaty High Energy Solar Spectroscopic Imager and the Tropical Rainfall Measuring Mission to determine the bulk (or footprint) microphysical properties of two types of study events, the thunderstorm complexes which are associated with TGFs (TGF case) and the thunderstorm complexes which did not produce a TGF detected by Reuven Ramaty High Energy Solar Spectroscopic Imager during the pass (non‐TGF case). Results are presented for two different comparison methods. The first case utilizes geographic region weighted by TGF distribution, and the second is based on TGF percentage of occurrence when compared to total flash count of data set. Results show that the associated storms around the TGF location possess differences in the hydrometeor concentrations: cloud liquid water, cloud ice, precipitation water, and precipitation ice. These results take place at different levels of the atmosphere, including the mixed phase region. Additionally, results will show that TGFs are a consistent percentage of observed flashes as the rate of TGFs as a function of Lightning Imaging Sensor flash count is relatively constant.

Properties of small cirrus ice crystals from commercial aircraft measurements and implications for flight operations

Measurements of cloud ice crystal size distributions have been made by a backscatter cloud probe (BCP) mounted on five commercial airliners flying international routes that cross five continents. Bulk cloud parameters were also derived from the size distributions. As of 31 December 2014, a total of 4399 flights had accumulated data from 665 hours in more than 19 000 cirrus clouds larger than 5 km in length. The BCP measures the equivalent optical diameter (EOD) of individual crystals in the 5–90 µm range from which size distributions are derived and recorded every 4 seconds. The cirrus cloud property database, an ongoing development stemming from these measurements, registers the total crystal number and mass concentration, effective and median volume diameters and extinction coefficients derived from the size distribution. This information is accompanied by the environmental temperature, pressure, aircraft position, date and time of each sample. The seasonal variations of the cirrus cloud properties measured from 2012 to 2014 are determined for six geographic regions in the tropics and extratropics. Number concentrations range from a few per litre for thin cirrus to several hundreds of thousands for heavy cirrus. Temperatures range from 205 to 250 K and effective radii from 12 to 20 µm. A comparison of the regional and seasonal number and mass size distributions, and the bulk microphysical properties derived from them, demonstrates that cirrus properties cannot be easily parameterised by temperature or by latitude. The seasonal changes in the size distributions from the extratropical Atlantic and Eurasian air routes are distinctly different, showing shifts from mono-modal to bi-modal spectra out of phase with one another. This phase difference may be linked to the timing of deep convection and cold fronts that lead to the cirrus formation. Likewise, the size spectra of cirrus over the tropical Atlantic and Eastern Brazil differ from each other although they were measured in adjoining regions. The cirrus crystals in the maritime continental tropical region over Malaysia form tri-modal spectra that are not found in any of the other regions measured by the IAGOS aircraft so far, a feature that is possibly linked to biomass burning or dust. Frequent measurements of ice crystal concentrations greater than 1×105 L−1, often accompanied by anomalously warm temperature and erratic airspeed readings, suggest that aircraft often experience conditions that affect their sensors. This new instrument, if used operationally, has the potential of providing real-time and valuable information to assist in flight operations as well as providing real-time information for along-track nowcasting.

Assessment of Global Cloud Datasets from Satellites: Project and Database Initiated by the GEWEX Radiation Panel

Clouds cover about 70% of Earth's surface and play a dominant role in the energy and water cycle of our planet. Only satellite observations provide a continuous survey of the state of the atmosphere over the entire globe and across the wide range of spatial and temporal scales that compose weather and climate variability. Satellite cloud data records now exceed more than 25 years; however, climate data records must be compiled from different satellite datasets and can exhibit systematic biases. Questions therefore arise as to the accuracy and limitations of the various sensors and retrieval methods. The Global Energy and Water Cycle Experiment (GEWEX) Cloud Assessment, initiated in 2005 by the GEWEX Radiation Panel (GEWEX Data and Assessment Panel since 2011), provides the first coordinated intercomparison of publicly available, standard global cloud products (gridded monthly statistics) retrieved from measurements of multispectral imagers (some with multiangle view and polarization capabilities), IR sounders, and lidar. Cloud properties under study include cloud amount, cloud height (in terms of pressure, temperature, or altitude), cloud thermodynamic phase, and cloud radiative and bulk microphysical properties (optical depth or emissivity, effective particle radius, and water path). Differences in average cloud properties, especially in the amount of high-level clouds, are mostly explained by the inherent instrument measurement capability for detecting and/or identifying optically thin cirrus, especially when overlying low-level clouds. The study of long-term variations with these datasets requires consideration of many factors. The monthly gridded database presented here facilitates further assessments, climate studies, and the evaluation of climate models.

Bulk microphysical properties of semi-transparent cirrus from AIRS: a six year global climatology and statistical analysis in synergy with geometrical profiling data from CloudSat-CALIPSO

Abstract. This article presents a retrieval method and a statistical analysis of the bulk microphysical properties of semi-transparent ice clouds using the Atmospheric Infrared Sounder (AIRS). The method relies on spectral differences of cirrus emissivities in the 8–12 μm range and is sensitive to the effective ice crystal diameter (De) and ice water path (IWP) of up to 85 μm and 120 g m−2, respectively. An indication of the most frequent ice crystal habit in the cirrus has been obtained by using separately single scattering properties of column-like and aggregate-like ice crystals in the simulations. Uncertainties due to hypotheses on atmospheric parameters and ice crystal single scattering properties are discussed and the cirrus emissivity and temperature range for the applicability of the method are determined. To be sure that the cirrus only includes ice crystals, one has to restrict the cloud temperature range to Tcld<230 K. On a global scale, these semi-transparent ice clouds (cirrus) represent about 25% of all high clouds and are mainly encountered in the midlatitudes during winter and in the tropics, with an average De and IWP of 52 μm and 27 g m−2, respectively. A comparison with bulk microphysical properties from the TIROS-N Operational Vertical Sounder (TOVS) shows an agreement on global mean values. The addition of spectral information revealed improvements at the limits of the cirrus emissivity range. Collocated Radar-Lidar Geometrical Profiling (GEOPROF) data have been used to study the vertical structure of these clouds and to infer average ice water content (IWC) for cirrus with a small vertical extent. This allowed us to compare and contrast parameterizations of De as functions of IWC and IWP, respectively.

Obtaining Best Estimates for the Microphysical and Radiative Properties of Tropical Ice Clouds from TWP-ICE In Situ Microphysical Observations

AbstractBest estimates of the bulk microphysical and radiative properties (ice water content, visible extinction, effective radius, and total concentration) are derived for three case studies of tropical ice clouds sampled during the Tropical Warm Pool International Cloud Experiment (TWP-ICE). Two case studies are aged cirrus clouds produced by deep convection (the so-called 27/01 and 29/01 cases), and the third (“02/02”) is a fresh anvil produced by deep convective activity over the Tiwi Islands. Using crystal images obtained by a Cloud Particle Imager (CPI), it is observed that small ice particles (with maximum dimensionD< 50–100μm) were predominantly quasi spherical, with the degree of nonsphericity increasing rapidly in the 50 <D< 100-μm range. ForD> 100μm, the aged cirrus clouds were predominantly characterized by bullet rosettes and aggregates of bullet rosettes, plates, and columns. In contrast, the fresh anvil had more frequent occurrences of plates, columns, aggregates of plates, and occasionally capped columns. The impact of shattering of large ice crystals on probe tips and the overall quality of the TWP-ICE in situ microphysical measurements are assessed. It is suggested that shattering has a relatively small impact on the CPI and cloud droplet probe (CDP) TWP-ICE data and a large impact on the Cloud Aerosol Spectrometer data, as already documented by others. It is also shown that the CPI size distributions must be multiplied by a factor of 4 to match those of the cloud imaging probe (CIP) for maximum dimension larger than 100μm (taken as a reference). A technique [named Best Estimate of Area and Density (BEAD)] to minimize errors associated with the density (ρ)–Dand projected area (A)–Dassumptions in bulk microphysics calculation is introduced and applied to the TWP-ICE data. The method makes direct use of the frequency of occurrence of each particle habit as classified from the CPI data and prescribedρ–DandA–Drelationships from the literature. This approach produces ice water content (IWC) estimates that are virtually unbiased relative to bulk measures obtained from a counterflow spectrometer and impactor (CSI) probe. In contrast, the use ofρ–DandA–Drelationships for single habits does produce large biases relative to the CSI observations: from −50% for bullet rosettes to +70%–80% for aggregates. The so-called width, length, area, and perimeter (WLAP) technique, which also makes use of individual CPI images, is found to produce positively biased IWCs (by 40% or so), and has a standard deviation of the errors similar to the BEAD technique. The impact of the large variability of the size distributions measured by different probe combinations on the bulk microphysical properties is characterized. The mean fractional differences with respect to the CSI measurements are small for the CPI + CIP, CPI, and CDP + CIP combinations (2.2%, −0.8%, and −1.1%, respectively), with standard deviations of the fractional differences ranging from 7% to 9%. This result provides an independent validation of the CPI scaling factor. The fractional differences produced between the CPI + CIP, CPI, and CDP + CIP combinations for extinction, effective radius, and total concentration are 33%, 10%–20%, and 90%, respectively, with relatively small standard deviations of 5%–8%. The fractional difference on total concentration varies greatly over the concentration range though, with values being larger than a factor of 2 for total concentrations smaller than 40 L−1, but reducing to 10%–20% for concentrations larger than 100 L−1. Therefore, caution should be exercised when using total concentrations smaller than 60–80 L−1as references for radar–lidar retrieval evaluation.

A method for comparing properties of cirrus clouds in global climate models with those retrieved from IR sounder satellite observations

A methodology to compare cloud properties simulated by global climate models with those retrieved from observations by satellite-based infrared (IR) sounders has been developed. The relatively high spectral resolution in the CO2 absorption band of these instruments leads to especially reliable cirrus properties, day and night. Additionally, bulk microphysical properties can be retrieved for semi-transparent cirrus, based on the observed spectral emissivity differences between 8 and 11 µm. The particular intention of this study is to compare macro-and microphysical properties of high cloudiness as represented by the model simulations and the satellite data. For this purpose, a method has been developed to process the model output to be comparable to the satellite measurements, as in other observational simulator packages (for example the ISCCP-simulator). This simulator method takes into account i) the differences in horizontal resolution of the model and the observations, ii) the specific observation time windows, iii) the determination of the pressure of a cloud system, identified with the pressure at the middle of the uppermost cloud, and iv) the selection of high clouds with specific cloud optical thickness ranges for the microphysical property retrieval using IR sounder data. Applying this method to simulations by the global climate model ECHAM and TOVS satellite observations has important effects. The frequency of high clouds selected from the model output by using the method is significantly smaller than the total frequency of high cloudiness in the model. Largest differences occur around the equator where the zonal mean frequency of high cloudiness is reduced by about 30 % (relative change). The selection method is essential for the comparison of modelled and observed microphysical properties of high clouds. The selection of high clouds from the ECHAM simulation according to the optical thickness range of the TOVS data results in a reduction of the mean water path of high clouds by factors of more than 3 compared to the case where also high clouds of other optical thicknesses are considered. Furthermore, the selection by optical thickness causes a significant increase in the mean effective cloud particle diameter. These changes significantly reduce the differences between the simulation and the observations. The method can also be applied for comparisons with other IR sounder climatologies such as from AIRS and IASI. © Gebruder Borntraeger, Stuttgart 2010.

Rapid-Scan Super-Resolution Observations of a Cyclic Supercell with a Dual-Polarization WSR-88D

Abstract In recent years, there has been widespread interest in collecting and analyzing rapid updates of radar data in severe convective storms. To this end, conventional single-polarization rapid-scan radars and phased array radar systems have been employed in numerous studies. However, rapid updates of dual-polarization radar data in storms are not widely available. For this study, a rapid scanning strategy is developed for the polarimetric prototype research Weather Surveillance Radar-1988 Doppler (WSR-88D) radar in Norman, Oklahoma (KOUN), which emulates the future capabilities of a polarimetric multifunction phased array radar (MPAR). With this strategy, data are collected over an 80° sector with 0.5° azimuthal spacing and 250-m radial resolution (“super resolution”), with 12 elevation angles. Thus, full volume scans over a limited area are collected every 71–73 s. The scanning strategy was employed on a cyclic nontornadic supercell storm in western Oklahoma on 1 June 2008. The evolution of the polarimetric signatures in the supercell is analyzed. The repetitive pattern of evolution of these polarimetric features is found to be directly tied to the cyclic occlusion process of the low-level mesocyclone. The cycle for each of the polarimetric signatures is presented and described in detail, complete with a microphysical interpretation. In doing so, for the first time the bulk microphysical properties of the storm on small time scales (inferred from polarimetric data) are analyzed. The documented evolution of the polarimetric signatures could be used operationally to aid in the detection and determination of various stages of the low-level mesocyclone occlusion.

Microphysical characterisation of West African MCS anvils

AbstractDeep convection in the Tropics is the source of large tropospheric extended clouds usually called anvils. These anvils may produce precipitation (stratiform region of deep convective storms), and also cirrus shields persisting from several hours to several days. Anvils impact the radiation budget, they induce a storage term in the whole water budget which is still relatively poorly quantified, and dynamical feedbacks may be induced in the upper troposphere. The AMMA field campaign which was held over West Africa during the 2006 wet season provided a unique opportunity to document the microphysics of these anvils from unprecedented airborne observations.Airborne in situ measurements of ice crystals and observations from a 95 GHz Doppler radar are used to characterize the microphysical properties of tropical anvils. The dataset is binned into stratiform and cirriform regions. Some data in the stratiform regions were likely obtained close enough to convective cores that the particles may have grown primarily within those cores. The data obtained over the continent and over the ocean are also characterized separately. Particle habit and growth processes are inferred from an examination of the collected particle images, from quantitative comparisons of 95 GHz reflectivities calculated from the in situ microphysical observations with the measured radar reflectivities, and from a statistical analysis of the two‐dimensional particle images. The predominant precipitation particles above the 0 °C isotherm in the stratiform anvil region are rimed aggregates. These rimed aggregates seem to get less dense and of smaller diameter when moving rearward of the system towards the cirriform region. The retrieved density laws (assumed to be power laws) lie close to the relationship for rimed particles of Locatelli and Hobbs.Particle size distributions in tropical anvils are also studied. The exponential shape seems to be a good approximation for these particle size distributions overall. The decrease in concentration with diameter is also found to be faster for cirriform regions than for stratiform regions. Normalising these particle size distributions produces a relatively invariant shape (in agreement with earlier studies), with however an increased variability for the smallest and largest values of the normalised diameter.The characterisation of the bulk microphysical properties using these in situ microphysical observations shows that the ice water content, the effective radius and the reflectivity‐weighted fall velocity generally increase with air temperature, in agreement with earlier studies. These parameters are found to be systematically smaller on average in cirriform regions than in stratiform regions, and this is true at all temperatures. These values are then compared with statistical relationships used in cloud‐resolving models and general circulation models, since a realistic representation of microphysics in models is very important to understand not only the processes at work, but the dynamical feedbacks and effects on climate. Large differences are found, the current parametrizations being unable to reproduce the large values of the considered microphysical parameters. Copyright © 2010 Royal Meteorological Society

Optical properties of a vertically inhomogeneous mid-latitude mid-level mixed-phase altocumulus in the infrared region

A vertically inhomogeneous mid-latitude mixed-phase altocumulus cloud was observed around 17:26 UTC on Oct. 14, 2001 during the 9th Cloud Layer Experiment (CLEX9). In this study the microphysical and optical properties of this cloud are investigated on the basis of in-situ observed vertical profiles of particle size and habit distributions. Two cloud models, assuming that the cloud properties were vertically homogeneous and inhomogeneous, are adopted to derive the bulk optical and radiative properties of this cloud. The observed microphysical properties are combined with the theoretical solutions to the scattering and absorption properties of individual cloud particles to determine the bulk optical properties at various heights within the cloud layer. The single-scattering properties of spherical liquid water droplets and nonspherical ice crystals are obtained from the Lorenz–Mie theory and an existing database, respectively. The bulk microphysical and optical properties associated with the inhomogeneous model depend strongly on the height above the cloud-base whereas the dependence is smoothed out in the case of the homogeneous model. Furthermore, the transfer of infrared radiation is simulated in conjunction with the two cloud models. It is shown that the brightness temperatures at the top of the atmosphere in the case of the homogeneous model can be 1.5% (3.8 K) higher than their counterpart associated with the inhomogeneous cloud model. This result demonstrates that the effect of the vertical inhomogeneity of a mixed-phase cloud on its radiative properties is not negligible.

Retrieval of Cirrus Microphysical Properties with a Suite of Algorithms for Airborne and Spaceborne Lidar, Radar, and Radiometer Data

AbstractAlgorithms are developed to convert data streams from multiple airborne and spaceborne remote sensors into layer-averaged cirrus bulk microphysical properties. Radiometers such as the Moderate-Resolution Imaging Spectroradiometer (MODIS) observe narrowband spectral radiances, and active remote sensors such as the lidar on the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite and the millimeter radar on CloudSat will provide vertical profiles of attenuated optical backscatter and radar reflectivity. Equivalent airborne remote sensors are also routinely flown on the NASA WB-57F and ER-2 aircraft. Algorithms designed to retrieve cirrus microphysical properties from remote sensor data must be able to handle the natural variability of cirrus that can range from optically thick layers that cause lidar attenuation to tenuous layers that are not detected by the cloud radar. An approach that is adopted here is to develop an algorithm suite that has internal consistency in its formulation and assumptions. The algorithm suite is developed around a forward model of the observations and is inverted for layer-mean cloud properties using a variational technique. The theoretical uncertainty in the retrieved ice water path retrieval is 40%–50%, and the uncertainty in the layer-mean particle size retrieval ranges from 50% to 90%. Two case studies from the Cirrus Regional Study of Tropical Anvils and Cirrus Layers (CRYSTAL) Florida Area Cirrus Experiment (FACE) field campaign as well as ground-based cases from the Atmospheric Radiation Measurement Program (ARM) are used to show the efficacy and error characteristics of the algorithms.

Cirrus Clouds and the Large-Scale Atmospheric State: Relationships Revealed by Six Years of Ground-Based Data

AbstractThe properties of cirrus clouds observed at the Atmospheric Radiation Measurement (ARM) Climate Research Facility (ACRF) in Oklahoma are documented from a nearly continuous 6-yr record of 35-GHz cloud radar data. Cirrus frequency over the ACRF is 23% and 28% of the time in the warm (May–September) and cold seasons (November–March), respectively, with maxima and minima during the period studied of 30% and 16% in the warm season and 34% and 24% in the cold seasons. Cirrus, as defined here, reveal a seasonal oscillation in their macroscale properties that can be traced to the seasonal deepening of the troposphere in the Southern Plains region. While the average bulk microphysical properties do not change significantly from season to season, the variability of certain parameters demonstrates seasonal change. It is shown that the properties of cirrus clouds vary perceptively with the large-scale vertical motion. Using NCEP–NCAR reanalysis data to define the large-scale meteorological state when cirrus are observed at the ACRF, the authors find that cirrus tend to exist within a maximum in upper-tropospheric humidity and downstream of the peak upper-tropospheric vertical motion. Cirrus that exist in large-scale ascent upstream of the synoptic-scale middle-tropospheric ridge axis are shown to have higher water contents than cirrus that exist in large-scale subsidence downstream of the ridge axis, although the overall nature of the statistical distributions of water contents do not change greatly, suggesting that it may be difficult to parameterize the properties of cirrus based solely on large-scale vertical motion. The layer-mean particle size, on the other hand, shows no such sensitivity to the large-scale vertical motion.

Bulk Scattering Properties for the Remote Sensing of Ice Clouds. Part I: Microphysical Data and Models

Abstract This study reports on the use of in situ data obtained in midlatitude and tropical ice clouds from airborne sampling probes and balloon-borne replicators as the basis for the development of bulk scattering models for use in satellite remote sensing applications. Airborne sampling instrumentation includes the two-dimensional cloud (2D-C), two-dimensional precipitation (2D-P), high-volume precipitation spectrometer (HVPS), cloud particle imager (CPI), and NCAR video ice particle sampler (VIPS) probes. Herein the development of a comprehensive set of microphysical models based on in situ measurements of particle size distributions (PSDs) is discussed. Two parameters are developed and examined: ice water content (IWC) and median mass diameter Dm. Comparisons are provided between the IWC and Dm values derived from in situ measurements obtained during a series of field campaigns held in the midlatitude and tropical regions and those calculated from a set of modeled ice particles used for light-scattering calculations. The ice particle types considered in this study include droxtals, hexagonal plates, solid columns, hollow columns, aggregates, and 3D bullet rosettes. It is shown that no single habit accurately replicates the derived IWC and Dm values, but a mixture of habits can significantly improve the comparison of these bulk microphysical properties. In addition, the relationship between Dm and the effective particle size Deff, defined as 1.5 times the ratio of ice particle volume to projected area for a given PSD, is investigated. Based on these results, a subset of microphysical models is chosen as the basis for the development of ice cloud bulk scattering models in Part II of this study.