Dual-polarization Weather Research Articles

Contextualizing Polarimetric Retrievals of Boundary Layer Height Using State-of-the-Art Boundary Layer Profiling

Abstract Knowledge about the depth of the planetary boundary layer (PBL) is crucial for a variety of applications, but direct observations of PBL depth are spatiotemporally sparse. Recent studies have proposed using operational dual-polarization weather radars to observe the evolution of PBL depth by capitalizing on unique differential reflectivity (ZDR) signatures of Bragg scatter at the top of the PBL. While this approach appears promising and cost-effective, uncertainties remain about the representativeness of these estimates and how its efficacy may vary by geography and climatology. To address these outstanding uncertainties, this study compares collocated observations collected from two WSR-88D radars and two state-of-the-art mobile boundary layer profiling systems and evaluates the proposed methodology over the full diurnal cycle. Results indicate good overall correspondence between the profiling- and radar-based PBL depth estimates, with an abrupt divergence during the early evening transition and large discrepancies overnight. Relatively large root-mean-square-deviations (RMSDs) coupled with small biases match expectations when comparing spatially averaged data with point observations during PBL growth, which capture frequent fluctuations. A qualitative examination of the radar data reveals signatures of elevated residual layers, clouds, and ground clutter, all of which can obfuscate the desired surface-based PBL signal but which may have their own utility. The prominence of the Bragg scatter signal is found to be correlated with the observed moisture gradient at the top of the PBL, reflecting climatological variability that should be considered. These findings motivate further work to improve the automated detection of Bragg scatter layers from polarimetric radar data. Significance Statement Knowledge of the height of the planetary boundary layer matters for weather forecasting, air quality, and renewable energy production. Currently, boundary layer height measurements are taken at select locations twice a day. However, a method to use the existing national network of polarimetric weather radars for this purpose has been proposed. This work evaluates this method against specialized boundary layer measurements. The results show that the method is generally reliable during the daytime and could be used for a variety of applications including climatologies and model evaluation. There remain a number of situational caveats, including residual turbulence, clouds/precipitation, ground clutter, and certain meteorological environments, that may require modification of the approach and need to be considered in future work.

Quantitative Precipitation Estimation Using Weather Radar Data and Machine Learning Algorithms for the Southern Region of Brazil

In addressing the challenges of quantitative precipitation estimation (QPE) using weather radar, the importance of enhancing the rainfall estimates for applications such as flash flood forecasting and hydropower generation management is recognized. This study employed dual-polarization weather radar data to refine the traditional Z–R relationship, which often needs higher accuracy in areas with complex meteorological phenomena. Utilizing tree-based machine learning algorithms, such as random forest and gradient boosting, this research analyzed polarimetric variables to capture the intricate patterns within the Z–R relationship. The results highlight machine learning’s potential to improve the precision of precipitation estimation, especially under challenging weather conditions. Integrating meteorological insights with advanced machine learning techniques is a remarkable achievement toward a more precise and adaptable precipitation estimation method.

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

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

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.

Estimation of extreme precipitation events in Estonia and Italy using dual-polarization weather radar quantitative precipitation estimations

Abstract. Evaluating extreme rainfall for a certain location is commonly considered when designing stormwater management systems. Rain gauge data are widely used to estimate rainfall intensities for a given return period. However, the poor spatial and temporal resolution of operational gauges is the main limiting factor. Several studies have used rainfall estimates based on weather radar horizontal reflectivity (Zh), but they come with a great caveat: while proven reliable for low or moderate rainfall rates, they are subject to major errors in extreme rainfall and convective cases. It is widely known that C-band weather radar can underestimate precipitation intensity due to signal attenuation or overestimate it due to hail and clutter contamination. From the late 1990s, dual-polarization weather radar started to become operational in the national surveillance radar network in Europe, providing innovative quantitative precipitation estimation (QPE) based on polarimetric variables. This study circumvents Zh shortcomings by using specific differential-phase (Kdp) data from operational dual-polarization C-band weather radars. The rain intensity estimates based on a specific differential-phase data are immune to attenuation and less affected by hail contamination. In this study, for the first time, QPEs based on polarimetric observations by operational C-band weather radars and without any rain gauge adjustments are analyzed. The purpose is to estimate return periods for 1 h rainfall total computed from polarimetric weather radar data using non-adjusted QPEs based on R(Zh,Kdp) data and to compare the results with those derived using R(Zh) and rain gauge data. Only the warm period during the year is considered here, as most of the extreme precipitation events for such a duration occur for both places studied (Italy and Estonia) at this time. Limiting the dataset to warm periods also allows us to use the radar-based rainfall quantitative precipitation estimations, which are more reliable than the snowfall ones. Data from operational dual polarimetric C-band weather radar sites are used from both Italy and Estonia. Given climatologically homogeneous regions, this study demonstrates that polarimetric weather radar observations can provide reliable QPEs compared to single-polarization estimates with respect to rain gauges and that they can provide a reliable estimation of return periods of 1 h rainfall total, even for relatively short time series.

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.

A Uniformity Index for Precipitation Particle Axis Ratios Derived from Radar Polarimetric Parameters for the Identification and Analysis of Raindrop Areas

A uniformity index for the axis ratios (Uar) derived from dual polarization weather radar data is proposed for raindrop area identification and analysis. The derivation of this new parameter is based on radar scattering simulations and assumptions. Uar is between 0 and 1 and can be calculated from the differential reflectivity (ZDR) and the copolar correlation coefficient (ρhv), which reflects the uniformity of the axis ratio (r) of the particle group. For raindrops, Uar is close to 1 under ideal conditions, but is clearly different from that of ice particles whose value is close to 0. Studies conducted during two convective weather events observed by X-band and S-band radar are presented to show the Uar features. In convective areas, high Uar presents a U-shaped vertical structure. One branch corresponds to the ZDR column, while the other branch is located at the rear of the convective cloud zone and is lower in altitude, representing the process of ice particles melting into raindrops and then being transported upward by a strong updraft. In stratiform cloud areas, a more than 95% overall identification ratio is obtained when the threshold of Uar is set to 0.2~0.3 for discriminating rain layers.

Distribution and evolution of hydrometeors in the stratiform cloud with embedded convection in the autumn ITCZ precipitation in Xisha: Case study

In autumn, the clouds over the South China Sea contain more cloud water and cloud ice. Intertropical Convergence Zone sometimes can strengthen and move north, causing heavy precipitation in the northern South China Sea. To reveal the distribution and evolution of hydrometeors in the Intertropical Convergence Zone precipitation clouds, a rainfall process occurred in Xisha and surrounding regions on 16 October 2021 was analyzed by utilizing S-band dual-polarization weather radar data and fuzzy logic algorithm. The classified hydrometeors showed that drizzle, rain, and dry snow were the three most abundant types, while dry crystal was less, indicating deposition and aggregation were more active in the marine environment with sufficient water vapor. The relative content of drizzle and dry snow changed oppositely to that of rain particles, suggesting the coalescence of drizzle and the transformation of dry snow were important processes affecting the formation of rain particles. The precipitation clouds were characterized by stratiform clouds with embedded convections. The strong updraft in the convective clouds transported liquid water upward, thus dry snow could collide with supercooled water and rime to form graupel, then graupel melted to form large raindrops below the 0°C layer. However, compared with continental convective clouds, the riming was weaker. In the stratiform clouds, the ascending motion was weak, no graupel was generated, and the 0°C-layer bright band indicated that dry snow could directly contribute to the small raindrops by melting. This study revealed the evolution of hydrometeors in the Intertropical Convergence Zone precipitation clouds and found that the increase in raindrop size in convective clouds was caused by the combined effects of stronger coalescence and riming.

Attenuation Correction in Weather Radars for Snow

Weather radars play a prominent role in remote sensing of the atmosphere. Various fields, such as meteorology and hydrology, rely on accurate weather radar data as input for their models. Different hydrometeors present during a weather event influence the amount of attenuation encountered by the radar signal. Attenuation correction for dual-polarization weather radar data is necessary to improve the radar products and get accurate measurements. Most of the existing attenuation correction research is associated with rain hydrometeors. Currently, research that addresses the attenuation correction of snow in weather radars is limited. Although it is known that attenuation of radar signals when it encounters rain is much greater than that for snow, attenuation for all hydrometeors needs to be addressed for accurate radar estimates. In this research work, the attenuation of different hydrometeors is studied using signal simulations. Various factors which influence attenuation, such as the elevation angle and particle size distribution, are considered, and the results are presented. An attenuation correction algorithm that uses the hydrometeor classification and specific differential phase products from the DROPS2.0 algorithm is introduced. Signal simulations are employed to obtain the relationship between specific attenuation and specific differential phase for different hydrometeors used in the proposed algorithm. The attenuation correction method is applied to X-band and Ku-band radar data. Path integrated attenuation of about 8 dB was observed in the snow case discussed from Ku band radar data. The method proposed for attenuation shows promising results at both frequency bands.

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.

Effects of Rough Hail Scattering on Polarimetric Variables

We use a commercially available full-wave electromagnetic (EM) tool to model scattering off rough hail. Through modeling various axis ratios and surface roughness, we systematically evaluate their effect on the polarimetric variables used for hail size discrimination. Our results produce the differential reflectivity and the copolar correlation coefficient typically not achieved using forward operators. We compare our results with available information in the literature. The inclusion of the additional shape parameter brings new insights into problems associated with the polarimetric variable for hailstone size gauging. Finally, using dual-polarization weather radars’ (WSR-88Ds’) observations of a hail storm on two opposite sides, we hypothesize that the source of negative differential reflectivity is wet oblate hail in the resonant size range.

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.