scholarly journals Validation of a Multilag Estimator on NJU-CPOL and a Hybrid Approach for Improving Polarimetric Radar Data Quality

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.

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.


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.


2012 ◽  
Vol 29 (6) ◽  
pp. 772-795 ◽  
Author(s):  
Lei Lei ◽  
Guifu Zhang ◽  
Richard J. Doviak ◽  
Robert Palmer ◽  
Boon Leng Cheong ◽  
...  

Abstract The quality of polarimetric radar data degrades as the signal-to-noise ratio (SNR) decreases. This substantially limits the usage of collected polarimetric radar data to high SNR regions. To improve data quality at low SNRs, multilag correlation estimators are introduced. The performance of the multilag estimators for spectral moments and polarimetric parameters is examined through a theoretical analysis and by the use of simulated data. The biases and standard deviations of the estimates are calculated and compared with those estimates obtained using the conventional method.


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.


2017 ◽  
Vol 56 (4) ◽  
pp. 877-896 ◽  
Author(s):  
Merhala Thurai ◽  
Patrick Gatlin ◽  
V. N. Bringi ◽  
Walter Petersen ◽  
Patrick Kennedy ◽  
...  

AbstractAnalysis of drop size distributions (DSD) measured by collocated Meteorological Particle Spectrometer (MPS) and a third-generation, low-profile, 2D-video disdrometer (2DVD) are presented. Two events from two different regions (Greeley, Colorado, and Huntsville, Alabama) are analyzed. While the MPS, with its 50-μm resolution, enabled measurements of small drops, typically for drop diameters below about 1.1 mm, the 2DVD provided accurate measurements for drop diameters above 0.7 mm. Drop concentrations in the 0.7–1.1-mm overlap region were found to be in excellent agreement between the two instruments. Examination of the combined spectra clearly reveals a drizzle mode and a precipitation mode. The combined spectra were analyzed in terms of the DSD parameters, namely, the normalized intercept parameter NW, the mass-weighted mean diameter Dm, and the standard deviation of mass spectrum σM. The inclusion of small drops significantly affected the NW and the ratio σM/Dm toward higher values relative to using the 2DVD-based spectra alone. For each of the two events, polarimetric radar data were used to characterize the variation of radar-measured reflectivity Zh and differential reflectivity Zdr with Dm from the combined spectra. In the Greeley event, this variation at S band was well captured for small values of Dm (<0.5 mm) where measured Zdr tended to 0 dB but Zh showed a noticeable decrease with decreasing Dm. For the Huntsville event, an overpass of the Global Precipitation Measurement mission Core Observatory satellite enabled comparison of satellite-based dual-frequency radar retrievals of Dm with ground-based DSD measurements. Small differences were found between the satellite-based radar retrievals and disdrometers.


2016 ◽  
Vol 55 (4) ◽  
pp. 829-848 ◽  
Author(s):  
Kiel L. Ortega ◽  
John M. Krause ◽  
Alexander V. Ryzhkov

AbstractThis study is the third part of a paper series investigating the polarimetric radar properties of melting hail and application of those properties for operational polarimetric hail detection and determination of its size. The results of theoretical simulations in Part I were used to develop a hail size discrimination algorithm (HSDA) described in Part II. The HSDA uses radar reflectivity Z, differential reflectivity ZDR, and cross-correlation coefficient ρhv along with melting-level height within a fuzzy-logic scheme to distinguish among three hail size classes: small hail (with diameter D < 2.5 cm), large hail (2.5 < D < 5.0 cm), and giant hail (D > 5.0 cm). The HSDA validation is performed using radar data collected by numerous WSR-88D sites and more than 3000 surface hail reports obtained from the Severe Hazards Analysis and Verification Experiment (SHAVE). The original HSDA version was modified in the process of validation, and the modified algorithm demonstrates probability of detection of 0.594, false-alarm ratio of 0.136, and resulting critical success index (CSI) equal to 0.543. The HSDA outperformed the current operational single-polarization hail detection algorithm, which only provides a single hail size estimate per storm and is characterized by CSI equal to 0.324. It is shown that HSDA is particularly sensitive to the quality of ZDR measurements, which might be affected by possible radar miscalibration and anomalously high differential attenuation.


Sign in / Sign up

Export Citation Format

Share Document