Microphysical Retrievals from Dual-Polarization Radar Measurements at X Band

Abstract The recent advances in attenuation correction methodology are based on the use of a constraint represented by the total amount of the attenuation encountered along the path shared over each range bin in the path. This technique is improved by using the inner self-consistency of radar measurements. The full self-consistency methodology provides an optimization procedure for obtaining the best estimate of specific and cumulative attenuation and specific and cumulative differential attenuation. The main goal of the study is to examine drop size distribution (DSD) retrieval from X-band radar measurements after attenuation correction. A new technique for estimating the slope of a linear axis ratio model from polarimetric radar measurements at attenuated frequencies is envisioned. A new set of improved algorithms immune to variability in the raindrop shape–size relation are presented for the estimation of the governing parameters characterizing a gamma raindrop size distribution. Simulations based on the use of profiles of gamma drop size distribution parameters obtained from S-band observations are used for quantitative analysis. Radar data collected by the NOAA/Earth System Research Laboratory (ESRL) X-band polarimetric radar are used to provide examples of the DSD parameter retrievals using attenuation-corrected radar measurements. Retrievals agree fairly well with disdrometer data. The radar data are also used to observe the prevailing shape of raindrops directly from the radar measurements. A significant result is that oblateness of drops is bounded between the two shape models of Pruppacher and Beard, and Beard and Chuang, the former representing the upper boundary and the latter the lower boundary

Rainfall rate retrieval in presence of path attenuation using C-band polarimetric weather radars

Weather radar systems are very suitable tools for the monitoring of extreme rainfall events providing measurements with high spatial and temporal resolution over a wide geographical area. Nevertheless, radar rainfall retrieval at C-band is prone to several error sources, such as rain path attenuation which affects the accuracy of inversion algorithms. In this paper, the so-called rain profiling techniques (namely the surface reference method FV and the polarimetric method ZPHI) are applied to correct rain path attenuation and a new neural network algorithm is proposed to estimate the rain rate from the corrected measurements of reflectivity and differential reflectivity. A stochastic model, based on disdrometer measurements, is used to generate realistic range profiles of raindrop size distribution parameters while a T-matrix solution technique is adopted to compute the corresponding polarimetric variables. A sensitivity analysis is performed in order to evaluate the expected errors of these methods. It has been found that the ZPHI method is more reliable than FV, being less sensitive to calibration errors. Moreover, the proposed neural network algorithm has shown more accurate rain rate estimates than the corresponding parametric algorithm, especially in presence of calibration errors

Retrieval of lower-order moments of the drop size distribution using CSU-CHILL X-band polarimetric radar: a case study

The lower-order moments of the drop size distribution (DSD) have generally been considered difficult to retrieve accurately from polarimetric radar data because these data are related to higher-order moments. For example, the 4.6th moment is associated with a specific differential phase and the 6th moment with reflectivity and ratio of high-order moments with differential reflectivity. Thus, conventionally, the emphasis has been to estimate rain rate (3.67th moment) or parameters of the exponential or gamma distribution for the DSD. Many double-moment “bulk” microphysical schemes predict the total number concentration (the 0th moment of the DSD, or M0) and the mixing ratio (or equivalently, the 3rd moment M3). Thus, it is difficult to compare the model outputs directly with polarimetric radar observations or, given the model outputs, forward model the radar observables. This article describes the use of double-moment normalization of DSDs and the resulting stable intrinsic shape that can be fitted by the generalized gamma (G-G) distribution. The two reference moments are M3 and M6, which are shown to be retrievable using the X-band radar reflectivity, differential reflectivity, and specific attenuation (from the iterative correction of measured reflectivity Zh using the total Φdp constraint, i.e., the iterative ZPHI method). Along with the climatological shape parameters of the G-G fit to the scaled/normalized DSDs, the lower-order moments are then retrieved more accurately than possible hitherto. The importance of measuring the complete DSD from 0.1 mm onwards is emphasized using, in our case, an optical array probe with 50 µm resolution collocated with a two-dimensional video disdrometer with about 170 µm resolution. This avoids small drop truncation and hence the accurate calculation of lower-order moments. A case study of a complex multi-cell storm which traversed an instrumented site near the CSU-CHILL radar is described for which the moments were retrieved from radar and compared with directly computed moments from the complete spectrum measurements using the aforementioned two disdrometers. Our detailed validation analysis of the radar-retrieved moments showed relative bias of the moments M0 through M2 was 0.9. Both radar measurement and parameterization errors were estimated rigorously. We show that the temporal variation of the radar-retrieved mass-weighted mean diameter with M0 resulted in coherent “time tracks” that can potentially lead to studies of precipitation evolution that have not been possible so far

Use of Dual Polarization Radar in Validation of Satellite Precipitation Measurements: Rationale and Opportunities

Dual-polarization weather radars have evolved significantly in the last three decades culminating in the operational deployment by the National Weather Service. In addition to operational applications in the weather service, dual-polarization radars have shown significant potential in contributing to the research fields of ground based remote sensing of rainfall microphysics, study of precipitation evolution and hydrometeor classification. Furthermore the dual-polarization radars have also raised the awareness of radar system aspects such as calibration. Microphysical characterization of precipitation and quantitative precipitation estimation are important applications that are critical in the validation of satellite borne precipitation measurements and also serves as a valuable tool in algorithm development. This paper presents the important role played by dual-polarization radar in validating space borne precipitation measurements. Starting from a historical evolution, the various configurations of dual-polarization radar are presented. Examples of raindrop size distribution retrievals and hydrometeor type classification are discussed. The quantitative precipitation estimation is a product of direct relevance to space borne observations. During the TRMM program substantial advancement was made with ground based polarization radars specially collecting unique observations in the tropics which are noted. The scientific accomplishments of relevance to space borne measurements of precipitation are summarized. The potential of dual-polarization radars and opportunities in the era of global precipitation measurement mission is also discussed

Method to combine spaceborne radar and radiometric observations of precipitation, A

2010 Fall.Includes bibliographical references.This dissertation describes the development and application of a combined radar-radiometer rainfall retrieval algorithm for the Tropical Rainfall Measuring Mission (TRMM) satellite. A retrieval framework based upon optimal estimation theory is proposed wherein three parameters describing the raindrop size distribution (DSD), ice particle size distribution (PSD), and cloud water path (cLWP) are retrieved for each radar profile. The retrieved rainfall rate is found to be strongly sensitive to the a priori constraints in DSD and cLWP; thus, these parameters are tuned to match polarimetric radar estimates of rainfall near Kwajalein, Republic of Marshall Islands. An independent validation against gauge-tuned radar rainfall estimates at Melbourne, FL shows agreement within 2% which exceeds previous algorithms' ability to match rainfall at these two sites. The algorithm is then applied to two years of TRMM data over oceans to determine the sources of DSD variability. Three correlated sets of variables representing storm dynamics, background environment, and cloud microphysics are found to account for approximately 50% of the variability in the absolute and reflectivity-normalized median drop size. Structures of radar reflectivity are also identified and related to drop size, with these relationships being confirmed by ground-based polarimetric radar data from the North American Monsoon Experiment (NAME). Regional patterns of DSD and the sources of variability identified herein are also shown to be consistent with previous work documenting regional DSD properties. In particular, mid-latitude regions and tropical regions near land tend to have larger drops for a given reflectivity, whereas the smallest drops are found in the eastern Pacific Intertropical Convergence Zone. Due to properties of the DSD and rain water/cloud water partitioning that change with column water vapor, it is shown that increases in water vapor in a global warming scenario could lead to slight (1%) underestimates of a rainfall trends by radar but larger overestimates (5%) by radiometer algorithms. Further analyses are performed to compare tropical oceanic mean rainfall rates between the combined algorithm and other sources. The combined algorithm is 15% higher than the version 6 of the 2A25 radar-only algorithm and 6.6% higher than the Global Precipitation Climatology Project (GPCP) estimate for the same time-space domain. Despite being higher than these two sources, the combined total is not inconsistent with estimates of the other components of the energy budget given their uncertainties

High-precision measurements of the co-polar correlation coefficient: non-Gaussian errors and retrieval of the dispersion parameter µ in rainfall

The co-polar correlation coefficient (ρhv) has many applications, including hydrometeor classification, ground clutter and melting layer identification, interpretation of ice microphysics and the retrieval of rain drop size distributions (DSDs). However, we currently lack the quantitative error estimates that are necessary if these applications are to be fully exploited. Previous error estimates of ρhv rely on knowledge of the unknown "true" ρhv and implicitly assume a Gaussian probability distribution function of ρhv samples. We show that frequency distributions of ρhv estimates are in fact highly negatively skewed. A new variable: L = -log10(1 - ρhv) is defined, which does have Gaussian error statistics, and a standard deviation depending only on the number of independent radar pulses. This is verified using observations of spherical drizzle drops, allowing, for the first time, the construction of rigorous confidence intervals in estimates of ρhv. In addition, we demonstrate how the imperfect co-location of the horizontal and vertical polarisation sample volumes may be accounted for. The possibility of using L to estimate the dispersion parameter (µ) in the gamma drop size distribution is investigated. We find that including drop oscillations is essential for this application, otherwise there could be biases in retrieved µ of up to ~8. Preliminary results in rainfall are presented. In a convective rain case study, our estimates show µ to be substantially larger than 0 (an exponential DSD). In this particular rain event, rain rate would be overestimated by up to 50% if a simple exponential DSD is assumed

Confronting the Challenge of Modeling Cloud and Precipitation Microphysics

In the atmosphere, microphysics refers to the microscale processes that affect cloud and precipitation particles and is a key linkage among the various components of Earth\u27s atmospheric water and energy cycles. The representation of microphysical processes in models continues to pose a major challenge leading to uncertainty in numerical weather forecasts and climate simulations. In this paper, the problem of treating microphysics in models is divided into two parts: (i) how to represent the population of cloud and precipitation particles, given the impossibility of simulating all particles individually within a cloud, and (ii) uncertainties in the microphysical process rates owing to fundamental gaps in knowledge of cloud physics. The recently developed Lagrangian particle‐based method is advocated as a way to address several conceptual and practical challenges of representing particle populations using traditional bulk and bin microphysics parameterization schemes. For addressing critical gaps in cloud physics knowledge, sustained investment for observational advances from laboratory experiments, new probe development, and next‐generation instruments in space is needed. Greater emphasis on laboratory work, which has apparently declined over the past several decades relative to other areas of cloud physics research, is argued to be an essential ingredient for improving process‐level understanding. More systematic use of natural cloud and precipitation observations to constrain microphysics schemes is also advocated. Because it is generally difficult to quantify individual microphysical process rates from these observations directly, this presents an inverse problem that can be viewed from the standpoint of Bayesian statistics. Following this idea, a probabilistic framework is proposed that combines elements from statistical and physical modeling. Besides providing rigorous constraint of schemes, there is an added benefit of quantifying uncertainty systematically. Finally, a broader hierarchical approach is proposed to accelerate improvements in microphysics schemes, leveraging the advances described in this paper related to process modeling (using Lagrangian particle‐based schemes), laboratory experimentation, cloud and precipitation observations, and statistical methods

Inferring microstructure and turbulence properties in rain through observations and simulations of signal spectra measured with Doppler–polarimetric radars

F. J. Yanovsky, Inferring microstructure and turbulence properties in rain through observations and simulations of signal spectra measured with Doppler–polarimetric radars (pp. 501-542). Book Chapter in: Polarimetric Detection, Characterization, and Remote Sensing, Springer, 2011.Doppler radars are able to measure important parameters of the target velocity. In contrast, polarimetric radars are very sensitive to features of the target shape and orientation relative to the radar beam direction. This chapter describes a novel Doppler–polarimetric approach to radar remote sensing. The combination of the Doppler ability and polarization diversity in the radar technology enables more comprehensive investigations of objects and phenomena in radar coverage. The discussion is adapted to the case of atmospheric remote sensing. A special case of cloud and precipitation observations is considered in greater detail. Mathematical models of signals and spectra of Doppler–polarimetric returns are discussed. It is demonstrated (theoretically, by simulation, and by real data processing) that important parameters of dynamic characteristics and microstructure of meteorological objects can be retrieved from Doppler–polarimetric observations. These results lead to new interesting and important applications like turbulence intensity measurement, drop size distribution estimation, recognition of type of scatterers, detection of hail zones, etc