Dual-polarization Weather Radar Research Articles

Implementation of a UAV-aided calibration method for a mobile dual-polarization weather radar

Weather radar calibration is a crucial factor to be considered for quantitative applications, such as QPE (Quantitative Precipitation Estimation), which is used as input for weather risks management. The present work proposes a novel approach to the end-to-end radar calibration method through the characterization of the radar weighting functions. These are Gaussian functions that model an additional attenuation factor to the radar received power. This approach, based on the inclusion these parameters, allow the obtainment of a calibrated equivalent reflectivity factor expression for a Doppler dual-polarization weather radar that operates in the X band. To calculate these parameters, a UAS (Unmanned Aircraft System) was implemented for suspending the calibration target with a well-defined cross-section and for measuring its inclination due to wind using an IMU (Inertial Measurement Unit). From its measurements, the position of the target can be estimated, which is essential to the characterization of the weighting functions. Their inclusion within the radar equation, alongside the implementation of the angular measurement system highlight the innovation to the traditional radar calibration methodology that does not contemplate them from the explored state-of-the-art. The reflectivity was compared with the measurements from a disdrometer for a moderate rain event. An average reflectivity difference of 0.75 dBZ and a percent bias of 3.3 % were obtained between the expected and estimated measurements when including these functions compared to the 1.51 dBZ and –62.7 % obtained when disregarding them. These experimental results point out that the proposed method can deliver superior accuracy in the reflectivity estimation.

A Weighted Adaptive Range-Averaging Technique to Improve the Precision Consistency of Polarimetric Variable Fields

Abstract The two main metrics for the performance evaluation of radar-variable estimators are the bias and standard deviation (SD) of estimates. Depending on the estimator properties, the bias may increase as the signal-to-noise ratio (SNR) decreases. The standard deviation, however, always rises as the SNR becomes smaller. For instance, if estimates are computed from 16 samples (typically used for WSR-88D surveillance scans) using a rectangular data window and the maximum unambiguous velocity is ∼9 m s−1, the standard deviation of reflectivity estimates increases 1.6 times as the SNR drops from 20 to 2 dB. But for estimates of differential reflectivity, differential phase, and copolar correlation coefficient, SDs increase ∼6.7, ∼6, and ∼54 times, respectively. Hence, this effect impacts the polarimetric variables substantially more than the spectral moments. Additionally, the polarimetric variable SD is also sensitive to the correlation between signals in horizontal and vertical channels leading to reduced data quality in the regions where the correlation coefficient is low. Such increases in the variability of estimates are observable in the fields of dual polarization variables as an increased spatial inhomogeneity (or noisiness) in the areas where radar echoes exhibit low-to-moderate SNRs and/or decreased correlation coefficient. These effects can obscure the visual identification of weather features as well as adversely impact algorithms. Herein, a novel method that applies variable smoothing in the range where the smoothing intensity depends on the SDs of estimates is presented. It applies little or no range averaging in the regions where data SDs are deemed adequate while using more aggressive smoothing in areas where data appear noisy. Significance Statement The noisiness in the fields of polarimetric variables is an issue that has plagued dual-polarization weather radars since their inception. This is because standard deviations of polarimetric variable estimates increase significantly more with the reduction in SNR than those of spectral moment estimates. A typical mitigation approach that indiscriminately averages a fixed number of estimates in the range may lead to unnecessary loss of range resolution in regions where data appearance is satisfactory. Further, such an approach can be inadequate in regions with a high variability of estimates leading to insufficient enhancement of weather feature visibility. In this study, a method that mitigates these issues is proposed.

An early warning approach for the rapid identification of extreme weather disasters based on phased array dual polarization radar cooperative network data.

In recent years, X-band phase array dual polarization weather radar technology has matured. The cooperative networking data from X-band phase array dual polarization weather radar have many advantages compared with traditional methods, namely, high spatial and temporal resolution (approximately 70 seconds in one scan, 30 m in radial distance resolution), wide coverages that can compensate for the observation blind spots, and data fusion technology that is used in the observation overlap area to ensure that the observed precipitation data have spatial continuity. Based on the above radar systems, this study proposes an improved hail and lightning weather disaster rapid identification and early warning algorithm. The improved thunderstorm identification, tracking, analysis, and nowcasting (TITAN) algorithm is used to quickly identify three-dimensional strong convective storm cells. Large sample observation experiment data are used to invert the localized hail index (HDR) to identify the hail position. The fuzzy logic method is used to comprehensively determine the probability of lightning occurrence. The comparative analysis experiment shows that, compared with the live observation data from the ground-based automatic station, the hail and lightning disaster weather warning algorithm developed by this study can increase warning times by approximately 7 minutes over the traditional algorithm, and its critical success index (CSI), false alarm ratio (FAR) and omission alarm ratio (OAR) scores are better than those of the traditional method. The average root mean square error (ARMSE) for identifying hail and lightning locations by this improved method is also significantly better than that of traditional methods. We show that our method can provide probabilistic predictions that improve hail and lightning identification, improve the precision of early warning and support operational utility at higher resolutions and with greater lead times that traditional methods struggle to achieve.

Doppler Velocity De-Aliasing Based on Lag-1 Cross-Correlation for Dual-Polarization Weather Radar

Doppler Velocity De-Aliasing Based on Lag-1 Cross-Correlation for Dual-Polarization Weather Radar

Simulation of Complex Meteorological Target Echoes for Airborne Dual-Polarization Weather Radar Based on Invariant Imbedding T -Matrix

Simulation of Complex Meteorological Target Echoes for Airborne Dual-Polarization Weather Radar Based on Invariant Imbedding T -Matrix

The Kinematic and Microphysical Characteristics of Extremely Heavy Rainfall in Zhengzhou City on 20 July 2021 Observed with Dual-Polarization Radars and Disdrometers

In this study, we utilized dual-polarization weather radar and disdrometer data to investigate the kinematic and microphysical characteristics of an extreme heavy rainfall event that occurred on 20 July 2021, in Zhengzhou. The results are as follows: FY-2G satellite images showed that extremely heavy rainfall mainly occurred during the merging period of medium- and small-scale convective cloud clusters. The merging of these cloud clusters enhanced the rainfall intensity. The refined three-dimensional wind field, as retrieved by the multi-Doppler radar, revealed a prominent mesoscale vortex and convergence structure at the extreme rainfall stage. This led to echo stagnation, resulting in localized extreme heavy rainfall. We explored the formation mechanism of the notable ZDR arc feature of dual-polarization variables during this phase. It was revealed that during the record-breaking hourly rainfall event in Zhengzhou (20 July 2021, 16:00–17:00 Beijing Time), the warm rain process dominated. Effective collision–coalescence processes, producing a high concentration of medium- to large-sized raindrops, significantly contributed to heavy rainfall at the surface. From an observational perspective, it was revealed that raindrops exhibited significant collision interactions during their descent. Moreover, a conceptual model for the kinematic and microphysical characteristics of this extreme rainfall event was established, aiming to provide technical support for monitoring and early warning of similar extreme rainfall events.

Observations of anomalous propagation over waters near Sweden

Abstract. Radio waves propagating in the atmosphere are affected by the prevailing atmospheric state. The state of the atmosphere can cause radio waves to refract more or less towards the ground. When the refractive index of the atmosphere differs from standard atmospheric conditions, the propagation is considered to be anomalous. Radars which are affected by anomalous propagation can observe ground clutter far beyond the radar horizon. In this work, 4.5 years' worth of data from five operational Swedish C-band dual-polarization weather radars are presented. Analyses of the data reveal a strong seasonal cycle and a weaker diurnal cycle in ground clutter from coastal regions across nearby waters. A comparison was drawn between the impacts of anomalous propagation on ground clutter measured with horizontal polarization and vertical polarization, respectively; however, no clear difference was found.

Study on the Backscatter Differential Phase Characteristics of X-Band Dual-Polarization Radar and its Processing Methods

The differential propagation phase (ΦDP) of X-band dual-polarization weather radar (including X-band dual-polarization phased-array weather radar, X-PAR) is important for estimating precipitation and classifying hydrometeors. However, the measured differential propagation phase contains the backscatter differential phase (δ), which poses difficulties for the application of the differential propagation phase from X-band radars. This paper presents the following: (1) the simulation and characteristics analysis of the backscatter differential phase based on disdrometer DSD (raindrop size distribution) measurement data; (2) an improved method of the specific differential propagation phase (KDP) estimation based on linear programming and backscatter differential phase elimination; (3) the effect of backscatter differential phase elimination on the specific differential propagation phase estimation of X-PAR. The results show the following: (1) For X-band weather radar, the raindrop equivalent diameters D > 2 mm may cause a backscatter differential phase between 0 and 20°; in particular, the backscatter differential phase varies sharply with raindrop size between 3.2 and 4.5 mm. (2) Using linear programming or smoothing filters to process the differential propagation phase could suppress the backscatter differential phase, but it is hard to completely eliminate the effect of the backscatter differential phase. (3) Backscatter differential phase correction may improve the calculation accuracy of the specific differential propagation phase, and the optimization was verified by the improved self-consistency of polarimetric variables, correlation between specific differential propagation phase estimations from S- and X-band radar and the accuracy of quantitative precipitation estimation. The X-PAR deployed in Shenzhen showed good observation performance and the potential to be used in radar mosaics with S-band weather radar.

X-Band Radar Attenuation Correction Method Based on LightGBM Algorithm

X-band weather radar can provide high spatial and temporal resolution data, which is essential to precipitation observation and prediction of mesoscale and microscale weather. However, X-band weather radar is susceptible to precipitation attenuation. This paper presents an X-band attenuation correction method based on the light gradient machine (LightGBM) algorithm (the XACL method), then compares it with the ZH correction method and the ZH-KDP comprehensive correction method. The XACL method was validated using observations from two radars in July 2021, the X-band dual-polarization weather radar at the Shouxian National Climatology Observatory of China (SNCOC), and the S-band dual-polarization weather radar at Hefei. During the rainfall cases on July 2021, the results of the attenuation correction were used for precipitation estimation and verified with the rainfall data from 1204 automatic ground-based meteorological network stations in Anhui Province, China. We found that the XACL method produced a significant correction effect and reduced the anomalous correction phenomenon of the comparison methods. The results show that the average error in precipitation estimation by the XACL method was reduced by 39.88% over 1204 meteorological stations, which is better than the effect of the other two correction methods. Thus, the XACL method proved good local adaptability and provided a new X-band attenuation correction scheme.

Correction: Wang et al. Characteristic Analysis of Dual-Polarization Weather Radar Echoes of Convective Precipitation and Snowfall in the Mount Everest Region. Atmosphere 2021, 12, 1671

In the published manuscript [...]

Signal Simulation of Dual-Polarization Weather Radar and Its Application in Range Ambiguity Mitigation

In this paper, a dual-polarization weather radar echo signal simulation method is proposed for the evaluation of the performance enhancement of dual-polarization weather radar systems, the optimization of signal processing algorithms and the improvement of scanning strategies. The actual weather radar base data are used in the simulation as the reference weather scene, which avoids using a complex algorithm for weather modeling. Moreover, based on radar weather equations, the radar system parameters are added into the radar echo signal modeling to establish the relationship between the simulated echo signal and radar system. As a result, the final simulated echo signal not only shows both the time and frequency domain characteristics of the weather target, but also includes the effects of the important performance of the dual-polarization weather radar system. In addition, to evaluate the performance of range ambiguity mitigation using phase coding and batch working modes, two different simulation methods for the radar signal are established on the method above; one is for batch working mode with long-PRT (pulse repetition time) and short-PRT transmission and receiving, and the other is for phase-coded mode with phase-coded transmission and phase-uncoded receiving. Under the same weather scene, the observation of these two different methods of range ambiguity mitigation are simulated and compared. Results show that the performance of the phase coding mode for mitigating range ambiguity is better than that of the batch mode. Obviously, the simulation method can be used to directly show the observation of different algorithms for mitigation range ambiguity under the same weather process, and quickly compare and evaluate the algorithm’s performance, which is not possible for real radars.

Characteristic Analysis of Dual-Polarization Weather Radar Echoes of Convective Precipitation and Snowfall in the Mount Everest Region

This paper introduces the X-band weather radar dual-polarization parameters of isolated convective cell precipitation and meso/microscale snowfall on Mount Everest and presents the first precipitation observations based on dual-polarization weather radar in this area. Compared with the Chengdu Plain, Mount Everest experienced convective precipitation on smaller horizontal and vertical scales with a narrower Zdr probability density spectrum (uniformly distributed around approximately 0). The Zh profile on Mount Everest displayed two peaks, unlike that over the plains, and the precipitation at the strong convective core was denser. Furthermore, during winter snowfall on the northern slope of Mount Everest, when the boundary layer exhibited sufficient water vapor and dynamic uplift, due to the low boundary layer temperature (<0 °C), water vapor produced stratiform clouds in the middle and lower layers (approximately 1.5 km above ground level (AGL)). Water vapor condensation at 1.5–2.5 km AGL led to latent heat release, which increased the temperature of regional stratiform clouds with increasing height. Consequently, the temperature at the stratiform cloud top height (2.5 km AGL) unexpectedly exceeded 0 °C. Additionally, the −20 °C isotherm was at approximately 4 km AGL, indicating that the middle- and upper-layer atmospheric temperatures remained low. Therefore, thermal instability occurred between the stratiform cloud top height and the middle/upper atmosphere, forming convective motion. These findings confirm the occurrence of elevated winter snowfall convection above Mount Everest and may have certain reference value for retrieving raindrop size distributions, quantitatively estimating precipitation, and parameterizing cloud microphysical processes in numerical prediction models for the Qinghai-Tibetan Plateau.

Characteristic analysis of dual-polarization weather radar echoes of clear air and precipitation in the Mount Everest region

X-band dual-polarization weather radar parameters characterized by clear-air turbulence and precipitation over Mount Everest are introduced. As the X-band radar wavelength is short, data quality control is exercised over the obtained differential phase (Φdp), reflectivity (Zh), differential reflectivity (Zdr), and other polarization physical quantities. Based on the X-band dual-polarization weather radar, FY4 satellite, and EC reanalysis data, the 3D structural characteristics of clear-air turbulence, weak precipitation, and meso-microscale strong convective precipitation on the northern slope of Mount Everest from June 2019 to July 2019 are analyzed. Our results reveal that clear-air turbulence on Mount Everest is mainly wet turbulence, and Zh and Zdr are significantly greater than those in the plain area. The ground clutter in some locations exhibits low Zdr and high ρhv characteristics, opposite to those in the plain area. Weak precipitation on Mount Everest has a typical three-layer structure similar to that in the plain area, and when Zh > 25 dBZ, Zh and Zdr display a positive linear correlation. Alternatively, Mount Everest’s atmospheric environment is conducive to triggering strong convective weather, consisting of common isolated convection cells at the microscale. The echo top height and maximum precipitation intensity of strong convection weather on Mount Everest are greater than those in the plain area. However, the whole layer convective thickness is much smaller than that in the plain area. Our preliminary results reveal the raindrop size distribution, quantitative precipitation estimation, and microphysical parameterization scheme of cloud precipitation in the Mount Everest region.

Real-Time Calibration and Monitoring of Radar Reflectivity on Nationwide Dual-Polarization Weather Radar Network

Monitoring calibration bias in reflectivity (ZH) in an operational S-band dual-polarization weather radar is the primary requisite for monitoring and prediction (nowcasting) of severe weather and routine weather forecasting using a weather radar network. For this purpose, we combined methods based on self-consistency (SC), ground clutter (GC) monitoring, and intercomparison to monitor the ZH in real time by complementing the limitations of each method. The absolute calibration bias can be calculated based on the SC between dual-polarimetric observations. Unfortunately, because SC is valid for rain echoes, it is impossible to monitor reflectivity during the non-precipitation period. GC monitoring is an alternative method for monitoring changes in calibration bias regardless of weather conditions. The statistics of GC ZH near radar depend on the changes in radar system status, such as antenna pointing and calibration bias. The change in GC ZH relative to the baseline was defined as the relative calibration adjustment (RCA). The calibration bias was estimated from the change in RCA, which was similar to that estimated from the SC. The ZH in the overlapping volume of adjacent radars was compared to verify the homogeneity of ZH over the radar network after applying the calibration bias estimated from the SC. The mean bias between two radars was approximately 0.0 dB after correcting calibration bias. We can conclude that the combined method makes it possible to use radar measurements, which are immune to calibration bias, and to diagnose malfunctioning radar systems as soon as possible.

Validation of X-band multi-parameter phased array weather radar by comparing data from Doppler weather radar with a parabolic dish antenna

AbstractThe multi-parameter phased array weather radar (MP-PAWR) was the first dual-polarized phased array weather radar to be commissioned in Japan (2017). When conducting a volume scan, the MP-PAWR respectively uses electronic and mechanical scanning in the elevation and azimuth angles to achieve rapid scanning and high-density observations. Although the effectiveness of the MP-PAWR has been demonstrated in case studies, its observation accuracy is yet to be quantitatively analyzed. Therefore, this study compared data of MP-PAWR with that of an operational dual-polarized weather radar with a parabolic-type antenna (X-MP radar) using 2,347,097 data samples obtained over 14 h. The results showed that the observation accuracy of the MP-PAWR was approximately the same as that of the X-MP radar at low elevations. The correlations of observational parameters (radar reflectivity factor, differential resistivity, specific differential phase, and Doppler velocity) between the MP-PAWR and X-MP radar ranged from 0.77–0.99 when MP-PAWR data were recorded within 15 s of the X-MP radar observations. The correlation between the observational parameters of the two radars decreased as the observation time difference between the X-MP radar and MP-PAWR increased. In particular, the correlation coefficients between the specific differential phase and the differential reflectivity were considerably lower than the single-polarization parameter at observation time difference of 240–300 s. By providing high-frequency and high-density dual-polarization observations, the MP-PAWR can contribute to rainfall prediction in Japan and reduce the damage caused by localized, rapidly developing cumulonimbus clouds.

Development of dual-polarization weather radar products display platform

Dual-polarization weather radar is more advanced than doppler weather radar, with higher observation accuracy and more products. Using PHP and Dreamweaver to develop dual-polarization weather radar products display platform. The dual-polarization weather radar products display platform can display ZDR, CC, KDP, HCl, QPE and other dual polarization radar products in real time, and has the function of playing back 30 times products in the past three hours of a single product. Using this platform, forecasters can compare and analyze doppler radar products with dual-polarization radar products to grasp the real-time conditions of short-term strong convective weather better, to improve the accuracy of short-term forecasts, especially disastrous weather forecasts, so as to further exert the advancement and superiority of dual-polarization weather radar. After applications, it shows that the dual-polarization radar products display platform is reliable and stable operation.

Clutter Mitigation Based on Spectral Depolarization Ratio for Dual-Polarization Weather Radars

In this article, a spectral polarimetric filter named moving spectral depolarization ratio (MsDR) filter is proposed to mitigate the effects of clutter and noise for dual-polarization weather radar without cross-polar measurements. The filter is based on a spectral polarimetric variable termed as spectral depolarization ratio (sDR). sDR can be estimated by any polarimetric weather radar, making the MsDR filter widely applicable. Taking advantage of the spectral polarimetric features and the range-Doppler continuity of precipitation, the MsDR filter is implemented in the range-Doppler spectrogram. The performance of the MsDR filter is assessed using data collected by 1) an X-band full-polarimetric radar interfered by narrow-band clutter (both stationary and moving) and sidelobe wind turbine clutter; 2) a C-band operational dual-polarization radar affected by radio frequency interference. In addition, the implementation of the MsDR filter is straightforward, and the computational complexity is relatively low. Hence, the MsDR filter has great potential to be applied in real-time for operational radar systems at different frequencies.

Analysis of the Tornadic Debris Signatures of the Ichihara Tornado in a Typhoon Environment Using Two Operational C-band Dual-Polarization Weather Radars

Analysis of the Tornadic Debris Signatures of the Ichihara Tornado in a Typhoon Environment Using Two Operational C-band Dual-Polarization Weather Radars

Characteristics of the Bright Band Based on Quasi-Vertical Profiles of Polarimetric Observations from an S-Band Weather Radar Network

Bright band (BB) characteristics obtained via dual-polarization weather radars elucidate thermodynamic and microphysical processes within precipitation systems. This study identified BB using morphological features from quasi-vertical profiles (QVPs) of polarimetric observations, and their geometric, thermodynamic, and polarimetric characteristics were statistically examined using nine operational S-band weather radars in South Korea. For comparable analysis among weather radars in the network, the calibration biases in reflectivity (ZH) and differential reflectivity (ZDR) were corrected based on self-consistency. The cross-correlation coefficient (ρHV) bias in the weak echo regions was corrected using the signal-to-noise ratio (SNR). First, we analyzed the heights of BBPEAK derived from the ZH as a function of season and compared the heights of BBPEAK derived from the ZH, ZDR, and ρHV. The heights of BBPEAK were highest in the summer season when the surface temperature was high. However, they showed distinct differences depending on the location (e.g., latitude) within the radar network, even in the same season. The height where the size of melting particles was at a maximum (BBPEAK from the ZH) was above that where the oblateness of these particles maximized (BBPEAK from ZDR). The height at which the inhomogeneity of hydometeors was at maximum (BBPEAK from the ρHV) was also below that of BBPEAK from the ZH. Second, BB thickness and relative position of BBPEAK were investigated to characterize the geometric structure of the BBs. The BB thickness increased as the ZH at BBBOTTOM increased, which indicated that large snowflakes melt more slowly than small snowflakes. The geometrical structure of the BBs was asymmetric, since the melting particles spent more time forming the thin shell of meltwater around them, and they rapidly collapsed to form a raindrop at the final stage of melting. Third, the heights of BBTOP, BBPEAK, and BBBOTTOM were compared with the zero-isotherm heights. The dry-temperature zero-isotherm heights were between BBTOP and BBBOTTOM, while the wet-bulb temperature zero-isotherm heights were close to the height of BBPEAK. Finally, we examined the polarimetric observations to understand the involved microphysical processes. The correlation among ZH at BBTOP, BBPEAK, and BBBOTTOM was high (>0.94), and the ZDR at BBBOTTOM was high when the BB’s intensity was strong. This proved that the size and concentration of snowflakes above the BB influence the size and concentration of raindrops below the BB. There was no depression in the ρHV for a weak BB. Finally, the mean profile of the ZH and ZDR depended on the ZH at BBBOTTOM. In conclusion, the growth process of snowflakes above the BB controls polarimetric observations of BB.

Echo simulation of dual polarization Doppler weather radar based on the physical model

Using the scattering characteristics of particles to simulate the radar echo can supply the test signals close to the real precipitation echo for the weather radar and save the time and cost of the research and development and maintenance of the weather radar. In this paper, the precipitation echo of weather radar is simulated based on the theoretical basis that the falling raindrops have a shape well approximated by an oblate spheroid in the atmosphere. The Marshal-Palmer distribution is applied to describe the raindrop spectrum distribution of precipitation particles. It is assumed that the raindrop particles of different sizes have the random distribution in the radar resolution volume, and then the spatial distribution of precipitation particles in the resolution volume is modeled. The echo signals of horizontal and vertical polarization channels of dual-polarization weather radar are obtained by vector superposition of backscattering echoes of each particle. The experimental results show that this method can describe the microphysical characteristics of precipitation particles more completely and can be used to test the signal processing module of dual-polarization Doppler weather radar.