RadxPid fuzzylogic

Overview

RadxPid uses polarimetric radar data to provide a best estimate of the dominant particle type at every radar gate. The categorization is done by using a fuzzy logic algorithm, which describes the likelihood that different values of polarimetric variables are caused by different particles. RadxPid uses the method described by Vivekanadan et al. (1999) and this page will highlight the most important concepts in the algorithm.

Fuzzy Logic Basics

As described in Vivekanadan et al. (1999), fuzzy logic methods are similar to probabilistic methods in that they describe uncertainty on a scale from 0 to 1, except fuzzy logic categories do not have hard boundaries. For example, radar gates that have rain or hail can have similar reflectivity values, although rain is more likely to have certain reflectivity values than hail. The fuzzy logic method assigns a likelihood (i.e., a membership value p) that rain or hail is associated with a particular reflectivity value defined by a membership function as shown in the image below. Membership functions can be created for a single polarimetric variable and combinations of polarimetric variables. Thus, membership functions are applied to the polarimetric data at each gate, providing a membership value (from 0 to 1) for each particle type and polarimetric variable/combination. Each membership value is then multiplied by a predetermined weight, the weighted values from each polarimetric variable/combination are summed, producing an aggregation value for each particle category. The particle type that has the largest aggregation value is then the dominant particle type at that gate and is assigned a particular PID value at that gate.

Fuzzy Logic in RadxPid

The membership functions and weights are all set in the PID thresholds file. Due to variations in the polarimetric responses of particles to radar characteristics, different files are required (e.g., S- vs X-band, simultaneous vs alternating transmission). The following examples will all assume a simultaneously transmitting S-band radar (e.g., WSR-88D).

18 PID categories are possible, including categories like heavy rain to hail to insects (lines 14-31). Weights for each variable are assigned on line 63 (e.g., temperature, [math]\displaystyle{ Z_H }[/math], [math]\displaystyle{ Z_{DR} }[/math], [math]\displaystyle{ K_{DP} }[/math], [math]\displaystyle{ LDR }[/math], [math]\displaystyle{ \rho_{HV} }[/math], [math]\displaystyle{ \sigma(Z_{DR}) }[/math], [math]\displaystyle{ \sigma(\phi_{DP}) }[/math].

References

Vivekanandan, J., Zrnic, D. S., Ellis, S. M., Oye, R., Ryzhkov, A. V., & Straka, J. (1999). Cloud Microphysics Retrieval Using S-Band Dual-Polarization Radar Measurements, Bulletin of the American Meteorological Society, 80(3), 381-388. Link