scholarly journals Rain microstructure retrievals using 2-D video disdrometer and C-band polarimetric radar

2009 ◽  
Vol 20 ◽  
pp. 13-18 ◽  
Author(s):  
M. Thurai ◽  
V. N. Bringi ◽  
W. A. Petersen

Abstract. Measurements using the 2-D video disdrometer (2DVD) taken during a heavy rainfall event in Huntsville, Alabama, are analysed. The 2DVD images were processed to derive the rain microstructure parameters for each individual drop, which in turn were used as input to the T-matrix method to compute the forward and back scatter amplitudes of each drop at C-band. The polarimetric radar variables were then calculated from the individual drop contribution over a finite time period, e.g., 1 min. The calculated co-polar reflectivity, differential reflectivity, specific differential propagation phase and the co-polar correlation coefficient were compared with measurements from a C-band polarimetric radar located 15 km away. An attenuation-correction method based on the specific differential propagation phase was applied to the co-polar and differential reflectivity data from the C-band radar, after ensuring accurate radar calibration. Time series comparisons of the parameters derived from the 2DVD and C-band radar data show very good agreement for all four quantities, the agreement being sometimes better than the computations using the 1-min drop size distribution and bulk assumptions on rain microstructure (such as mean shapes and model-based assumptions for drop orientation). The agreement is particularly improved in the case of co-polar correlation coefficient since this parameter is very sensitive to variation of shapes as well as orientation angles. The calculations mark the first attempt at utilizing experimentally derived "drop- by-drop" rain microstructure information to compute the radar polarimetric parameters and to demonstrate the value of utilizing the 2-D video disdrometer for studying rain microstructure under various precipitation conditions. Histograms of drop orientation angles as well as the most probable drop shapes and the corresponding variations were also derived and compared with prior results from the 80 m fall "artificial rain" experiment.


2020 ◽  
Vol 12 (1) ◽  
pp. 180
Author(s):  
Shiqing Shao ◽  
Kun Zhao ◽  
Haonan Chen ◽  
Jianjun Chen ◽  
Hao Huang

For the estimation of weak echo with low signal-to-noise ratio (SNR), a multilag estimator is developed, which has better performance than the conventional method. The performance of the multilag estimator is examined by theoretical analysis, simulated radar data and some specific observed data collected by a C-band polarimetric radar in previous research. In this paper, the multilag estimator is implemented and verified for Nanjing University C-band polarimetric Doppler weather radar (NJU-CPOL) during the Observation, Prediction and Analysis of Severe Convection of China (OPACC) field campaign in 2014. The implementation results are also compared with theoretical analysis, including the estimation of signal power, spectrum width, differential reflectivity, and copolar correlation coefficient. The results show that the improvement of the multilag estimator is little for signal power and differential reflectivity, but significant for spectrum width and copolar correlation coefficient when spectrum width is less than 2 ms−1, which implies a large correlation time scale. However, there are obvious biases from the multilag estimator in the regions with large spectrum width. Based on the performance analysis, a hybrid method is thus introduced and examined through NJU-CPOL observations. All lags including lag 0 of autocorrelation function (ACF) are used for moment estimation in this algorithm according to the maximum usable lag number. A case study shows that this hybrid method can improve moment estimation compared to both conventional estimator and multilag estimator, especially for weak weather echoes. The improvement will be significant if SNR decreases or the biases of noise power in the conventional estimator increase. In addition, this hybrid method is easy to implement on both operational and non-operational radars. It is also expected that the proposed hybrid method will have a better performance if applied to S-band polarimetric radars which have twice the maximum useable lags in the same conditions with C-band radars.



2018 ◽  
Vol 2018 ◽  
pp. 1-14 ◽  
Author(s):  
Hongli Li ◽  
Xiangde Xu ◽  
Yang Hu ◽  
Yanjiao Xiao ◽  
Zhibin Wang

Operational Doppler radar observations have potential advantages over other above-surface observations when it comes to assimilation for mesoscale model simulations with high spatial and temporal resolution. To improve the forecast of a heavy frontal rainfall event that occurred in the Yangtze-Huaihe River Basin from 4 July to 5 July 2014 in China, operational radar observations are assimilated by the Local Analysis and Prediction System (LAPS). Radar reflectivity data are used primarily in the LAPS cloud analysis procedure, which retrieves the number of hydrometeors and adjusts the moisture and cloud fields. Radial velocity data are analyzed through the LAPS wind analysis-based successive correction method. A new correction method is developed to correct three-dimensional radar reflectivity data based on hourly surface rain gauge observations. The performance of the correction method is demonstrated by assimilating radar reflectivity observations into LAPS. Experiments with different radar data assimilation are examined. Results show that the assimilation of radar data can effectively correct the background errors and improve the heavy rainfall forecast. The simulated intensity, pattern, and temporal evolution of the heavy rainfall event are better improved with radar reflectivity assimilation, especially when the correction method is implemented to correct radar observations.



2021 ◽  
Vol 13 (16) ◽  
pp. 3060
Author(s):  
Muyun Du ◽  
Jidong Gao ◽  
Guifu Zhang ◽  
Yunheng Wang ◽  
Pamela L. Heiselman ◽  
...  

Polarimetric radar data (PRD) have potential to be used in numerical weather prediction (NWP) models to improve convective-scale weather forecasts. However, thus far only a few studies have been undertaken in this research direction. To assimilate PRD in NWP models, a forward operator, also called a PRD simulator, is needed to establish the relation between model physics parameters and polarimetric radar variables. Such a forward operator needs to be accurate enough to make quantitative comparisons between radar observations and model output feasible, and to be computationally efficient so that these observations can be easily incorporated into a data assimilation (DA) scheme. To address this concern, a set of parameterized PRD simulators for the horizontal reflectivity, differential reflectivity, specific differential phase, and cross-correlation coefficient were developed. In this study, we have tested the performance of these new operators in a variational DA system. Firstly, the tangent linear and adjoint (TL/AD) models for these PRD simulators have been developed and checked for the validity. Then, both the forward operator and its adjoint model have been built into the three-dimensional variational (3DVAR) system. Finally, some preliminary DA experiments have been performed with an idealized supercell storm. It is found that the assimilation of PRD, including differential reflectivity and specific differential phase, in addition to radar radial velocity and horizontal reflectivity, can enhance the accuracy of both initial conditions for model hydrometer state variables and ensuing model forecasts. The usefulness of the cross-correlation coefficient is very limited in terms of improving convective-scale data analysis and NWP.



Author(s):  
Dana M. Tobin ◽  
Matthew R. Kumjian

AbstractA unique polarimetric radar signature indicative of hydrometeor refreezing during ice pellet events has been documented in several recent studies, yet the underlying microphysical causes remain unknown. The signature is characterized by enhancements in differential reflectivity (ZDR), specific differential phase (KDP), and linear depolarization ratio (LDR), and a reduction in co-polar correlation coefficient (ρhv) within a layer of decreasing radar reflectivity factor at horizontal polarization (ZH). In previous studies, the leading hypothesis for the observed radar signature is the preferential refreezing of small drops. Here, a simplified, one-dimensional, explicit bin microphysics model is developed to simulate the refreezing of fully melted hydrometeors, and coupled with a polarimetric radar forward operator to quantify the impact of preferential refreezing on simulated radar signatures. The modeling results demonstrate that preferential refreezing is insufficient by itself to produce the observed signatures. In contrast, simulations considering an ice shell growing asymmetrically around a freezing particle (i.e., emulating a thicker ice shell on the bottom of a falling particle) produce realistic ZDR enhancements, and also closely replicate observed features in ZH, KDP, LDR, and ρhv. Simulations that assume no increase in particle wobbling with freezing produce an even greater ZDR enhancement, but this comes at the expense of reducing the LDR enhancement. It is suggested that the polarimetric refreezing signature is instead strongly related to both the distribution of the unfrozen liquid portion within a freezing particle, and the orientation of this liquid with respect to the horizontal.



Atmosphere ◽  
2019 ◽  
Vol 10 (12) ◽  
pp. 784 ◽  
Author(s):  
Anil Kumar Khanal ◽  
Guy Delrieu ◽  
Frédéric Cazenave ◽  
Brice Boudevillain

The RadAlp experiment aims at developing advanced methods for rain and snow estimation using weather radar remote sensing techniques in high mountain regions for improved water resource assessment and hydrological risk mitigation. A unique observation system has been deployed in the French Alps, Grenoble region. It is composed of a Météo-France operated X-band MOUC radar (volumetric, Doppler and polarimetric) on top of the Mt Moucherotte (1920 m ASL), the X-band XPORT research radar (volumetric, Doppler, polarimetric), a K-band micro rain radar (MRR, Doppler, vertically pointing) and in situ sensors (rain gauges, disdrometers), latter three operated on the Grenoble campus (220 m ASL). Based on the observation that the precipitation phase changes at/below the elevation of mountain-top MOUC radar for more than 60% of the significant events, an algorithm for ML identification has been developed using valley-based radar systems: it uses the quasi vertical profiles of XPORT polarimetric measurements (horizontal and vertical reflectivity, differential reflectivity, cross-polar correlation coefficient) and the MRR vertical profiles of apparent falling velocity spectra. The algorithm produces time series of the altitudes and values of peaks and inflection points of the different radar observables. A literature review allows us to link the micro-physical processes at play during the melting process with the available polarimetric and Doppler signatures, e.g., (i) regarding the altitude differences between the peaks of reflectivity, cross-polar correlation coefficient and differential reflectivity, as well as (ii) regarding the co-variation of the profiles of Doppler velocity spectra and cross-polar correlation coefficient. A statistical analysis of the ML based on 42 rain events (98 h of XPORT data) is then proposed. Among other results, this study indicates that (i) the mean value of the ML width in Grenoble is 610 m with a standard deviation of 160 m; (ii) the mean altitude difference between the horizontal reflectivity and the ρ H V peaks is 90 m and the mean altitude difference between the ρ H V and Zdr peaks is 30 m; (iii) even for the limited rainrate range in the dataset (0–8.5 mm h − 1 ), the “intensity effect” is clear on the reflectivity profile and the ML width, as well as on polarimetric variables such as ρ H V peak value and the Zdr enhancement in the upper part of the profile. On the contrary, the study of both the “density effect” and the influence of the 0   ° C isotherm altitude did not yield significant results with the considered dataset; (iv) a principal component analysis on one hand shows the richness of the dataset since the first 2 PCs explain only 50% of the total variance and on the other hand the added-value of the polarimetric variables since they rank high in a ranking of the total variance explained by individual variables.



2020 ◽  
Vol 12 (22) ◽  
pp. 3711
Author(s):  
Chih-Chien Tsai ◽  
Kao-Shen Chung

Based on the preciousness and uniqueness of polarimetric radar observations collected near the landfall of Typhoon Soudelor (2015), this study investigates the sensitivities of very short-range quantitative precipitation forecasts (QPFs) for this typhoon to polarimetric radar data assimilation. A series of experiments assimilating various combinations of radar variables are carried out for the purpose of improving a 6 h deterministic forecast for the most intense period. The results of the control simulation expose three sources of the observation operator errors, including the raindrop shape-size relation, the limitations for ice-phase hydrometeors, and the melting ice model. Nevertheless, polarimetric radar data assimilation with the unadjusted observation operator can still improve the analyses, especially rainwater, and consequent QPFs for this typhoon case. The different impacts of assimilating reflectivity, differential reflectivity, and specific differential phase are only distinguishable at the lower levels of convective precipitation areas where specific differential phase is found most helpful. The positive effect of radar data assimilation on QPFs can last three hours in this study, and further improvement can be expected by optimizing the observation operator in the future



2017 ◽  
Vol 145 (12) ◽  
pp. 4899-4910 ◽  
Author(s):  
Kendell T. LaRoche ◽  
Timothy J. Lang

A pyrocumulus is a convective cloud that can develop over a wildfire. Under certain conditions, pyrocumulus clouds become vertically developed enough to produce lightning. NEXRAD dual-polarization weather radar and upgraded National Lightning Detection Network (NLDN) data were used to analyze 10 case studies of ash plumes and pyrocumulus clouds from 2013 that either did or did not produce detected lightning. Past research has shown that pyrocumulus cases are most likely to produce lightning when there is a decrease in differential reflectivity (toward 0 dB) and an increase in the correlation coefficient (to >0.8), as measured by polarimetric radar, due to the transition from pure smoke/ash to frozen hydrometeors. All pyrocumulus cases that produced lightning in this study displayed the polarimetric characteristics of rimed ice within their respective clouds. Time series analysis of radar-inferred ash and rimed ice volumes within ash plumes and pyrocumulus clouds showed that NLDN-detected lightning occurred only after the cloud contained significant amounts of precipitation-sized rimed ice. The results suggest that the recently dual-pol-enabled NEXRADs and the more sensitive NLDN network can be used to explore ash plume and pyrocumulus microphysical structure and lightning production.



2013 ◽  
Vol 52 (1) ◽  
pp. 169-185 ◽  
Author(s):  
Qing Cao ◽  
Guifu Zhang ◽  
Ming Xue

AbstractThis study presents a two-dimensional variational approach to retrieving raindrop size distributions (DSDs) from polarimetric radar data in the presence of attenuation. A two-parameter DSD model, the constrained-gamma model, is used to represent rain DSDs. Three polarimetric radar measurements—reflectivity ZH, differential reflectivity ZDR, and specific differential phase KDP—are optimally used to correct for the attenuation and retrieve DSDs by taking into account measurement error effects. Retrieval results with simulated data demonstrate that the proposed algorithm performs well. Applications to real data collected by the X-band Center for Collaborative Adaptive Sensing of the Atmosphere (CASA) radars and the C-band University of Oklahoma–Polarimetric Radar for Innovations in Meteorology and Engineering (OU-PRIME) also demonstrate the efficacy of this approach.



2009 ◽  
Vol 26 (9) ◽  
pp. 1829-1842 ◽  
Author(s):  
Eugenio Gorgucci ◽  
V. Chandrasekar ◽  
Luca Baldini

Abstract A method is proposed to retrieve raindrop shape–size relations from the radar measurements of reflectivity factor Zh, differential reflectivity Zdr, and specific differential phase Kdp at S band. This procedure is obtained using a domain defined by the two variables Kdp/Zh and Zdr where the drop size distribution (DSD) variability is collapsed onto a line and any variation is essentially due to the drop shape variability. To obtain information on the raindrop shape–size relation underlying a set of radar observations, this domain is studied in conjunction with another domain describing the relation between the drop axial ratio (or shape) and its equivolumetric diameter. Using an initial drop shape and choosing a set of DSDs described by a normalized gamma model, polarimetric radar measurements are produced by simulation. An averaged curve of Kdp/Zh versus Zdr is obtained and compared with the same curve obtained from the radar data. By changing the initial axial ratio relation, a procedure of minimization between the two curves is developed to derive the underlying drop shape–size relation governing the radar measurements under consideration. Three sets of radar data collected in different climatic regions are analyzed to evaluate whether there is a unique shape–size relation.



Atmosphere ◽  
2020 ◽  
Vol 11 (7) ◽  
pp. 770
Author(s):  
Matthew Van Den Broeke

Disdrometer and condensation nuclei (CN) data are compared with operational polarimetric radar data for one multicell and one supercell storm in eastern Nebraska on 11 June 2018. The radar was located ~14.3 km from the instrumentation location and provided excellent observation time series with new low-level samples every 1–2 min. Reflectivity derived by the disdrometer and radar compared well, especially in regions with high number concentration of drops and reflectivity <45 dBZ. Differential reflectivity also compared well between the datasets, though it was most similar in the supercell storm. Rain rate calculated by the disdrometer closely matched values estimated by the radar when reflectivity and differential reflectivity were used to produce the estimate. Concentration of CN generally followed precipitation intensity for the leading convective cell, with evidence for higher particle concentration on the edges of the convective cell associated with outflow. The distribution of CN in the supercell was more complex and generally did not follow precipitation intensity.



Sign in / Sign up

Export Citation Format

Share Document