Evaluation of radar reflectivity factor simulations of ice crystal populations from in situ observations for the retrieval of condensed water content in tropical mesoscale convective systems

This study presents the evaluation of a technique to estimate cloud condensed water content (CWC) in tropical convection from airborne cloud radar reflectivity factors at 94 GHz and in situ measurements of particle size distributions (PSDs) and aspect ratios of ice crystal populations. The approach is to calculate from each 5 s mean PSD and flight-level reflectivity the variability of all possible solutions of m(D) relationships fulfilling the condition that the simulated radar reflectivity factor (T-matrix method) matches the measured radar reflectivity factor. For the reflectivity simulations, ice crystals were approximated as oblate spheroids, without using a priori assumptions on the mass–size relationship of ice crystals. The CWC calculations demonstrate that individual CWC values are in the range ±32 % of the retrieved average CWC value over all CWC solutions for the chosen 5 s time intervals. In addition, during the airborne field campaign performed out of Darwin in 2014, as part of the international High Altitude Ice Crystals/High Ice Water Content (HAIC/HIWC) projects, CWCs were measured independently with the new IKP-2 (isokinetic evaporator probe) instrument along with simultaneous particle imagery and radar reflectivity. Retrieved CWCs from the T-matrix radar reflectivity simulations are on average 16 % higher than the direct CWCIKP measurements. The differences between the CWCIKP and averaged retrieved CWCs are found to be primarily a function of the total number concentration of ice crystals. Consequently, a correction term is applied (as a function of total number concentration) that significantly improves the retrieved CWC. After correction, the retrieved CWCs have a median relative error with respect to measured values of only −1 %. Uncertainties in the measurements of total concentration of hydrometeors are investigated in order to calculate their contribution to the relative error of calculated CWC with respect to measured CWCIKP. It is shown that an overestimation of the concentration by about +50 % increases the relative errors of retrieved CWCs by only +29 %, while possible shattering, which impacts only the concentration of small hydrometeors, increases the relative error by about +4 %. Moreover, all cloud events with encountered graupel particles were studied and compared to events without observed graupel particles. Overall, graupel particles seem to have the largest impact on high crystal number-concentration conditions and show relative errors in retrieved CWCs that are higher than for events without graupel particles

The variability of tropical ice cloud properties as a function of the large-scale context from ground-based radar-lidar observations over Darwin, Australia

The high complexity of cloud parameterizations now held in models puts more pressure on observational studies to provide useful means to evaluate them. One approach to the problem put forth in the modelling community is to evaluate under what atmospheric conditions the parameterizations fail to simulate the cloud properties and under what conditions they do a good job. It is the ambition of this paper to characterize the variability of the statistical properties of tropical ice clouds in different tropical "regimes" recently identified in the literature to aid the development of better process-oriented parameterizations in models. For this purpose, the statistical properties of non-precipitating tropical ice clouds over Darwin, Australia are characterized using ground-based radar-lidar observations from the Atmospheric Radiation Measurement (ARM) Program. The ice cloud properties analysed are the frequency of ice cloud occurrence, the morphological properties (cloud top height and thickness), and the microphysical and radiative properties (ice water content, visible extinction, effective radius, and total concentration). The variability of these tropical ice cloud properties is then studied as a function of the large-scale cloud regimes derived from the International Satellite Cloud Climatology Project (ISCCP), the amplitude and phase of the Madden-Julian Oscillation (MJO), and the large-scale atmospheric regime as derived from a long-term record of radiosonde observations over Darwin. The vertical variability of ice cloud occurrence and microphysical properties is largest in all regimes (1.5 order of magnitude for ice water content and extinction, a factor 3 in effective radius, and three orders of magnitude in concentration, typically). 98 % of ice clouds in our dataset are characterized by either a small cloud fraction (smaller than 0.3) or a very large cloud fraction (larger than 0.9). In the ice part of the troposphere three distinct layers characterized by different statistically-dominant microphysical processes are identified. The variability of the ice cloud properties as a function of the large-scale atmospheric regime, cloud regime, and MJO phase is large, producing mean differences of up to a factor 8 in the frequency of ice cloud occurrence between large-scale atmospheric regimes and mean differences of a factor 2 typically in all microphysical properties. Finally, the diurnal cycle of the frequency of occurrence of ice clouds is also very different between regimes and MJO phases, with diurnal amplitudes of the vertically-integrated frequency of ice cloud occurrence ranging from as low as 0.2 (weak diurnal amplitude) to values in excess of 2.0 (very large diurnal amplitude). Modellers should now use these results to check if their model cloud parameterizations are capable of translating a given atmospheric forcing into the correct statistical ice cloud properties

On the observation of unusual high concentration of small chain-like aggregate ice crystals and large ice water contents near the top of a deep convective cloud during the CIRCLE-2 experiment

During the CIRCLE-2 experiment carried out over Western Europe in May 2007, combined in situ and remote sensing observations allowed to describe microphysical and optical properties near-top of an overshooting convective cloud (11 080 m/−58 °C). The airborne measurements were performed with the DLR Falcon aircraft specially equipped with a unique set of instruments for the extensive in situ cloud measurements of microphysical and optical properties (Polar Nephelometer, FSSP-300, Cloud Particle Imager and PMS 2-D-C) and nadir looking remote sensing observations (DLR WALES Lidar). Quasi-simultaneous space observations from MSG/SEVIRI, CALIPSO/CALIOP-WFC-IIR and CloudSat/CPR combined with airborne RASTA radar reflectivity from the French Falcon aircraft flying above the DLR Falcon depict very well convective cells which overshoot by up to 600 m the tropopause level. Unusual high values of the concentration of small ice particles, extinction, ice water content (up to 70 cm−3, 30 km−1 and 0.5 g m−3, respectively) are experienced. The mean effective diameter and the maximum particle size are 43 μm and about 300 μm, respectively. This very dense cloud causes a strong attenuation of the WALES and CALIOP lidar returns. The SEVIRI retrieved parameters confirm the occurrence of small ice crystals at the top of the convective cell. Smooth and featureless phase functions with asymmetry factors of 0.776 indicate fairly uniform optical properties. Due to small ice crystals the power-law relationship between ice water content (IWC) and radar reflectivity appears to be very different from those usually found in cirrus and anvil clouds. For a given equivalent reflectivity factor, IWCs are significantly larger for the overshooting cell than for the cirrus. Assuming the same prevalent microphysical properties over the depth of the overshooting cell, RASTA reflectivity profiles scaled into ice water content show that retrieved IWC up to 1 g m−3 may be observed near the cloud top. Extrapolating the relationship for stronger convective clouds with similar ice particles, IWC up to 5 g m−3 could be experienced with reflectivity factors no larger than about 20 dBZ. This means that for similar situations, indication of rather weak radar echo does not necessarily warn the occurrence of high ice water content carried by small ice crystals. All along the cloud penetration the shape of the ice crystals is dominated by chain-like aggregates of frozen droplets. Our results confirm previous observations that the chains of ice crystals are found in a continental deep convective systems which are known generally to generate intense electric fields causing efficient ice particle aggregation processes. Vigorous updrafts could lift supercooled droplets which are frozen extremely rapidly by homogeneous nucleation near the −37 °C level, producing therefore high concentrations of very small ice particles at upper altitudes. They are sufficient to deplete the water vapour and suppress further nucleation as confirmed by humidity measurements. These observations address scientific issues related to the microphysical properties and structure of deep convective clouds and confirm that particles smaller than 50 μm may control the radiative properties in convective-related clouds. These unusual observations may also provide some possible insights regarding engineering issues related to the failure of jet engines commonly used on commercial aircraft during flights through areas of high ice water content. However, large uncertainties of the measured and derived parameters limit our observations

A Comparison of Airborne In Situ Cloud Microphysical Measurement with Ground-Based C-Band Radar Observations in Deep Stratiform Regions of African Squall Lines

This study addresses clouds with significant ice water content (IWC) in the stratiform regions downwind of the convective cores of African squall lines in the framework of the French–Indian satellite Megha-Tropiques project, observed in August 2010 next to Niamey (13.5°N, 2°E) in the southwestern part of Niger. The objectives included comparing the IWC–Z reflectivity relationship for precipitation radars in deep stratiform anvils, collocating reflectivity observed from ground radar with the calculated reflectivity from in situ microphysics for all aircraft locations inside the radar range, and interpreting the role of large ice crystals in the reflectivity of centimeter radars through analysis of their microphysical characteristics as ice crystals larger than 5 mm frequently occurred. It was found that, in the range of 20–30 dBZ, IWC and C-band reflectivity are not really correlated. Cloud regions with high IWC caused by important crystal number concentrations can lead to the same reflectivity factor as cloud regions with low IWC formed by a few millimeter-sized ice crystals

An Integrated Approach to Weather Radar Calibration and Monitoring Using Ground Clutter and Satellite Comparisons

The stability and accuracy of weather radar reflectivity calibration are imperative for quantitative applications, such as rainfall estimation, severe weather monitoring and nowcasting, and assimilation in numerical weather prediction models. Various radar calibration and monitoring techniques have been developed, but only recently have integrated approaches been proposed, that is, using different calibration techniques in combination. In this paper the following three techniques are used: 1) ground clutter monitoring, 2) comparisons with spaceborne radars, and 3) the self-consistency of polarimetric variables. These techniques are applied to a C-band polarimetric radar (CPOL) located in the Australian tropics since 1998. The ground clutter monitoring technique is applied to each radar volumetric scan and provides a means to reliably detect changes in calibration, relative to a baseline. It is remarkably stable to within a standard deviation of 0.1 dB (decibels). To obtain an absolute calibration value, CPOL observations are compared to spaceborne radars on board TRMM (Tropical Rainfall Measuring Mission) and GPM (Global Precipitation Measurement) using a volume-matching technique. Using an iterative procedure and stable calibration periods identified by the ground echoes technique, we improve the accuracy of this technique to about 1 dB. Finally, we review the self-consistency technique and constrain its assumptions using results from the hybrid TRMM-GPM and ground echo technique. Small changes in the self-consistency parameterization can lead to 5 dB of variation in the reflectivity calibration. We find that the drop-shape model of Brandes et al. with a standard deviation of the canting angle of 12 degrees best matches our dataset

Coastal Observations of Weather Features in Senegal during the AMMA SOP-3 Period

During 15 August through 30 September 2006, ground and aircraft measurements were obtained from a multi-national group of students and scientists in Senegal. Key measurements were aimed at investigating and understanding precipitation processes, thermodynamic and dynamic environmental conditions, cloud, aerosol and microphysical processes and spaceborne sensors (TRMM, CloudSat/Calipso) validation. Ground and aircraft instruments include: ground based polarimetric radar, disdrometer measurements, a course and a high-density rain gauge network, surface chemical measurements, a 10 m flux tower, broadband IR, solar and microwave measurements, rawinsonde and radiosonde measurements, FA-20 dropsonde, in situ microphysics and cloud radar measurements. Highlights during SOP3 include ground and aircraft measurements of squall lines, African Easterly Waves (AEWs), Saharan Air Layer advances into Senegal, and aircraft measurements of AEWs -- including the perturbation that became Hurricane Isaac

Toward understanding of differences in current cloud retrievals of ARM ground-based measurements

Accurate observations of cloud microphysical properties are needed for evaluating and improving the representation of cloud processes in climate models and better estimate of the Earth radiative budget. However, large differences are found in current cloud products retrieved from ground-based remote sensing measurements using various retrieval algorithms. Understanding the differences is an important step to address uncertainties in the cloud retrievals. In this study, an in-depth analysis of nine existing ground-based cloud retrievals using ARM remote sensing measurements is carried out. We place emphasis on boundary layer overcast clouds and high level ice clouds, which are the focus of many current retrieval development efforts due to their radiative importance and relatively simple structure. Large systematic discrepancies in cloud microphysical properties are found in these two types of clouds among the nine cloud retrieval products, particularly for the cloud liquid and ice particle effective radius. Note that the differences among some retrieval products are even larger than the prescribed uncertainties reported by the retrieval algorithm developers. It is shown that most of these large differences have their roots in the retrieval theoretical bases, assumptions, as well as input and constraint parameters. This study suggests the need to further validate current retrieval theories and assumptions and even the development of new retrieval algorithms with more observations under different cloud regimes

An Improved Instability–Shear Hail Proxy for Australia

We evaluated the performance in Australia of proxies designed to identify atmospheric conditions prone to hail and severe storms. In a convection-resolving but short-duration simulation, proxies that use instability and wind shear thresholds overestimated the probability of hail occurring when compared to the estimated occurrence of surface graupel in the model, particularly in Australia’s tropical north. We used reanalysis data and the Australian Bureau of Meteorology severe storm archive to examine atmospheric conditions at times and locations when hailstorms, other storms, and no storms were reported between January 1979 and March 2021. In instability–shear space, the best discriminator between hail and no-storm times was found to vary predictably with melting-level height, allowing a new proxy to better represent latitudinal trends in atmospheric conditions. We found extra conditions that can be applied to the new proxy to efficiently reduce the number of false alarms. The new proxy outperforms the tested existing proxies for detection of hail-prone conditions in Australia

Measurements of rainfall rate, drop size distribution, and variability at middle and higher latitudes: application to the combined DPR-GMI algorithm

The Global Precipitation Measurement mission is a major U.S.–Japan joint mission to understand the physics of the Earth’s global precipitation as a key component of its weather, climate, and hydrological systems. The core satellite carries a dual-precipitation radar and an advanced microwave imager which provide measurements to retrieve the drop size distribution (DSD) and rain rates using a Combined Radar-Radiometer Algorithm (CORRA). Our objective is to validate key assumptions and parameterizations in CORRA and enable improved estimation of precipitation products, especially in the middle-to-higher latitudes in both hemispheres. The DSD parameters and statistical relationships between DSD parameters and radar measurements are a central part of the rainfall retrieval algorithm, which is complicated by regimes where DSD measurements are abysmally sparse (over the open ocean). In view of this, we have assembled optical disdrometer datasets gathered by research vessels, ground stations, and aircrafts to simulate radar observables and validate the scattering lookup tables used in CORRA. The joint use of all DSD datasets spans a large range of drop concentrations and characteristic drop diameters. The scaling normalization of DSDs defines an intercept parameter NW, which normalizes the concentrations, and a scaling diameter Dm, which compresses or stretches the diameter coordinate axis. A major finding of this study is that a single relationship between NW and Dm, on average, unifies all datasets included, from stratocumulus to heavier rainfall regimes. A comparison with the NW–Dm relation used as a constraint in versions 6 and 7 of CORRA highlights the scope for improvement of rainfall retrievals for small drops (Dm lt; 1 mm) and large drops (Dm gt; 2 mm). The normalized specific attenuation–reflectivity relationships used in the combined algorithm are also found to match well the equivalent relationships derived using DSDs from the three datasets, suggesting that the currently assumed lookup tables are not a major source of uncertainty in the combined algorithm rainfall estimates